Aerostructures Chapter 2
Aerostructures Chapter 2
i.e.
1.29 + 8.14 =
2
72 + 4xy
3Q
2 150 300
from which
Q = 95 100 N = 95.1 kN.
4p
P (2,3)
3p
3p
x
O
4p
4p
2p
Fig. S.2.1
=
E
E
E
2p
3p
2.75p
y =
+
=
E
E
E
x =
(i)
(ii)
so that u =
3.5p
x + f1 (y)
E
(iii)
15
16
Solutions Manual
so that v =
2.75p
y + f2 (x)
E
(iv)
4p
v u
f2 (x) f1 (y)
=
+
=
+
G
x
y
x
y
Suppose
f1 (y)
=A
y
then
f1 ( y) = Ay + B
(v)
(vi)
3.5p
x + Ay + B
E
(vii)
and
2.75p
y + Cx + D
(viii)
E
Since the origin of the axes is xed in space it follows that when x = y = 0, u = v = 0.
Hence, from Eqs (vii) and (viii), B = D = 0. Further, the direction of Ox is xed in space
so that, when y = 0, v/x = 0. Therefore, from Eq. (viii), C = 0. Thus, from Eqs (1.28)
and (vii), when x = 0.
v=
4p
u
=
=A
y
G
Eqs (vii) and (viii) now become
u=
3.5p
4p
x+ y
E
G
(ix)
2.75p
y
E
From Eq. (1.50), G = E/2(1 + ) = E/2.5 and Eq. (ix) becomes
p
u = (3.5x + 10y)
E
v=
(x)
(xi)
23p
E
and
8.25p
(from Eq. (x))
E
The point P therefore moves at an angle to the x axis given by
v=
= tan1
8.25
= 19.73
23
S.2.2
An Airy stress function, , is dened by the equations (Eqs (2.8)):
x =
2
y2
y =
2
x 2
xy =
2
x y
and has a nal form which is determined by the boundary conditions relating to a
particular problem.
Since
= Ay3 + By3 x + Cyx
4
=0
x 4
4
=0
y4
(i)
4
=0
x 2 y2
2
= 6Ay + 6Byx
y2
2
=0
x 2
2
= 3By2 C
=
x y
(ii)
y =
(iii)
xy
(iv)
The distribution of shear stress in a rectangular section beam is parabolic and is zero
at the upper and lower surfaces. Hence, when y = d/2, xy = 0. Thus, from Eq. (iv)
B=
4C
3d 2
(v)
17
18
Solutions Manual
d/2
(3By2 C)t dy = P
d/2
which gives
2t
Bd 3
Cd
+
8
2
=P
3P
2td
(vi)
2P
td 3
(vii)
At the free end of the beam where x = l the bending moment is zero and thus x = 0
for any value of y. Therefore, from Eq. (ii)
6A + 6Bl = 0
whence
A=
2Pl
td 3
(viii)
12Pl
12P
y 3 xy
td 3
td
or
12P(l x)
y
(ix)
td 3
Equation (ix) is the direct stress distribution at any section of the beam given by
simple bending theory, i.e.
My
x =
I
3
where M = P(l x) and I = td /12.
The shear stress distribution given by Eq. (iv) is
x =
xy =
6P 2 3P
y
td 3
2td
or
6P 2 d 2
(x)
xy = 3 y
td
4
Equation (x) is identical to that derived from simple bending theory and may be found
in standard texts on stress analysis, strength of materials, etc.
S.2.3
The stress function is
=
w
(15h2 x 2 y 5x 2 y3 2h2 y3 + y5 )
20h3
Then
2
w
=
(30h2 y 10y3 ) = y
x 2
20h3
w
2
=
(30x 2 y 12h2 y + 20y3 ) = x
2
y
20h3
w
2
=
(30h2 x 30xy2 ) = xy
x y
20h3
4
=0
x 4
4
w
=
(120y)
4
y
20h3
4
w
=
(60y)
2
2
x y
20h3
Substituting in Eq. (2.9)
4 = 0
so that the stress function satises the biharmonic equation.
The boundary conditions are as follows:
At y = h, y = w and xy = 0 which are satised.
At y = h, y = w and xy = 0 which are satised.
At x = 0, x = w/20h3 (12h2 y + 20y3 ) = 0.
Also
h
h
h
w
(12h2 y + 20y3 )dy
20h3 h
w
=
[6h2 y2 + 5y4 ]hh
20h3
=0
x dy =
19
20
Solutions Manual
Finally
h
h
h
w
(12h2 y2 + 20y4 )dy
20h3 h
w
=
[4h2 y3 + 4y5 ]hh
20h3
=0
x y dy =
S.2.4
The Airy stress function is
=
p
[5(x 3 l 2 x)(y + d)2 (y 2d) 3yx(y2 d 2 )2 ]
120d 3
Then
4
=0
x 4
4
3pxy
= 3
4
y
d
4
3pxy
=
2
2
x y
2d 3
3pxy 3pxy
3 =0
2d 3
d
2
px
=
[5y(x 2 l 2 ) 10y3 + 6d 2 y]
2
y
20d 3
pl
(10y3 + 6d 2 y)
20d 3
x 1 dy
pl
=
20d 3
d
d
(10y3 + 6d 2 y)dy = 0
Thus the stress function satises the boundary conditions for axial load in the x direction.
Also, the direct stress, y , is given by (see Eqs (2.8))
y =
2
px
= 3 ( y3 3yd 2 2d 3 )
x 2
4d
px
(d 3 + 3d 3 2d 3 ) = 0
4d 3
y =
px 3
(d 3d 3 2d 3 ) = px
4d 3
and when y = +d
Thus, the stress function satises the boundary conditions for load in the y direction.
The shear stress, xy , is given by (see Eqs (2.8))
xy =
p
2
=
[5(3x 2 l 2 )( y2 d 2 ) 5y4 + 6y2 d 2 d 4 ]
x y
40d 3
When x = 0
xy =
p
[5l 2 ( y2 d 2 ) 5y4 + 6y2 d 2 d 4 ]
40d 3
xy 1 dy =
p
40d 3
d
d
pl2
6
1 2
lpl l
2 3
i.e.
xy (x = 0) =
pl2
6d
S.2.5
The stress function is
=
w
(10c3 x 2 15c2 x 2 y + 2c2 y3 + 5x 2 y3 y5 )
40bc3
21
22
Solutions Manual
Then
w
2
=
(12c2 y + 30x 2 y 20y3 ) = x
y2
40bc3
w
2
=
(20c3 30c2 y + 10y3 ) = y
x 2
40bc3
w
2
=
(30c2 x + 30xy2 ) = xy
x y
40bc3
4
=0
x 4
4
w
=
(120y)
y4
40bc3
4
w
=
(60y)
x 2 y2
40bc3
Substituting in Eq. (2.9)
4 = 0
so that the stress function satises the biharmonic equation.
On the boundary, y = +c
w
xy = 0
y =
b
At y = c
y = 0
At x = 0
x =
Then
xy = 0
w
(12c2 y 20y3 )
40bc3
c
w
x dy =
(12c2 y 20y3 )dy
40bc3 c
c
w
[6c2 y2 5y4 ]cc
=
40bc3
=0
c
i.e. the direct stress distribution at the end of the cantilever is self-equilibrating.
The axial force at any section is
c
c
w
x dy =
(12c2 y + 30x 2 y 20y3 )dy
40bc3 c
c
w
=
[6c2 y2 + 15x 2 y2 5y4 ]cc
40bc3
=0
i.e. no axial force at any section of the beam.
c
w
(12c2 y2 20y4 )dy
40bc3 c
w
[4c2 y3 4y5 ]cc = 0
=
40bc3
i.e. the beam is a cantilever beam under a uniformly distributed load of w/unit area with
a self-equilibrating stress application at x = 0.
S.2.6
From physics, the strain due to a temperature rise T in a bar of original length L0 and
nal length L is given by
=
L0 (1 + T ) L0
L L0
=
= T
L0
L0
x =
2(1 + )
xy
E
E
x y
E x 2
x 2
x 2
E y2
y2
y2
or
2(1 + )
2 y
2 y
2 xy
2 x
2 x
+
+ E 2 T
=
x 2
y2
x y
y2
x 2
2 xy
2 y
= 2
x y
y
Hence
2
2 y
2 xy
2 x
= 2
x y
x
y2
(i)
23
24
Solutions Manual
and
2
2 y
2 xy
2 x
= 2 2
x y
x
y
Thus
2 y
2 y
2 x
2 x
=
+
+ E 2 T
x 2
y2
x 2
y2
2
2
+
x 2
y2
(x + y ) + E 2 T = 0
and since
x =
2
y2
y =
2
2
+ 2
2
x
y
2
x 2
2 2
+ 2
y2
x
+ E 2 T = 0
or
2 ( 2 + ET ) = 0
S.2.7
The stress function is
=
3Qxy Qxy3
4a
4a3
Then
2
= 0 = y
x 2
3Qxy
2
= 3 = x
y2
2a
2
3Q 3Qy2
=
= xy
x y
4a
4a3
Also
4
=0
x 4
4
=0
y4
4
=0
x 2 y2
y=0
y = a
x =
y = +a
3Q
2a
x = 0
3Q
x =
2a
y = a
y = +a
x =
3Q
2a
y2
3Q
1 2 , i.e. parabolic
=
4a
a
y2
3Q
1 2 dy
4a
a
i.e.
SF = Q.
The resultant bending moment at x = a is
=
x y dy
3Qay2
dy
2a3
i.e.
BM = Qa.
25