MCE 526 - Theory of Elasticity Dr. M.

Sahin
2016-2017 Spring Semester Ankara Yildirim Beyazit University

MIDTERM EXAM III SOLUTIONS

Name :

17
Student No : Department :

1. [35] Consider a plate with sides 2a and b in x and y directions respectively. The stress
distribtion is given as

3 1

20
σx = px3 y − 2c1 xy + c2 y, σy = pxy 3 − 2px3 y, τxy = − px2 y 2 + c1 y 2 + px4 + c3
2 2
a) Obtain the corresponding stress function

b) Verify that this field is a solution for a thin plate of thickness t

c) Find the resultant normal and boundary forces (Py and Vx ) along the edges y = 0 and y = b
of the plate

Solution
6
52
1 a) The stress function
Integrate τxy wrt y first

∂ 2 φ (x, y) ∂ 2 φ (x, y) ∂φ (x, y)

Z Z
τxy (x, y) = − =⇒ τxy (x, y) dy = − dy = −
∂x∂y ∂x∂y ∂x

c1 px2 px4
Z      
∂φ (x, y) 3 2 2 2 1 4 3
=− − px y + c1 y + px + c3 dy = −y − −y c3 + +f (x)
∂x 2 2 3 2 2
CE

Take derivative of the previous equation wrt to x to get σy i.e df (x)

dx
= 0 =⇒ f (x) = C5

∂ 2 φ (x, y) df (x) df (x)

σy = 2
= −2px3 y + pxy 3 + = pxy 3 − 2px3 y =⇒ =0
∂x dx dx
find f (x) after ordinary integration wrt to x

∂φ (x, y) 1 2 3 y3 1 4
= px y − c1 − px y − c3 y + C5
M

∂x 2 3 2
Integrate wrt to x to get φ (x, y)
∂φ(x,y)
∂x

c1 px2 px4
Z     
∂φ (x, y)
Z
3
φ (x, y) = dx = C5 − y − − y c3 + dx + g (y)
∂x 3 2 2
py 5 px3 3 c1 y 3
 
φ (x, y) = − x + y + x C5 − − c3 y + g (y)
10 6 3

1
MCE 526 - Theory of Elasticity Dr. M. Sahin
2016-2017 Spring Semester Ankara Yildirim Beyazit University

Take the derivative wrt y twice

py 5 px3 3 c1 y 3
 
φ (x, y) = − x + y + x C5 − − c3 y + g (y)
10 6 3

17
∂φ (x, y) px5 px3 2
y + x −c1 y 2 − c3 + g 0 (y)


=− +
∂y 10 2
∂ 2 φ (x, y)
= xy −2c1 + px2 + g 00 (y) = σx = px3 y − 2c1 xy + c2 y


∂y 2

1 1
=⇒ g 00 (y) = c2 y =⇒ g 0 (y) = c2 y 2 + C6 =⇒ g (y) = c2 y 3 + C6 y + C7

20
2 6
py 5 px3 3 c1 y 3
 
1
φ (x, y) = − x + y + x C5 − − c3 y + c2 y 3 + C 6 y + C 7
10 6 3 6

1b) Use equilibrium equation σij,j = 0 or bihormanic equation ∇4 φ = 0 to show that the Airy
stress function is solution for the plate
∂σx ∂τxy
∂x
+

+
∂y
∂τxy ∂σy
=0

=0
6 =⇒ 3pyx2 − 2c1 y − 3pyx2 + 2c1 y = 0 =⇒ OK

=⇒ −3py 2 x + 2px3 + 3pxy 2 − 2px3 = 0 =⇒ OK

52
∂x ∂y

1c) The stresses at these lines

σx = px3 y − 2c1 xy + c2 y σy = pxy 3 − 2px3 y σy (x, y = 0) = 0

Z
Py (y = 0) = σy (x, y = 0) dy = 0
CE

Z Z a
pxb3 − 2px3 b dx = 0


Py (y = b) = σy (x, y = b) dy =
−a
a
1
Z

a
b2 x − 2x3 xdx = pb b2 x2 − x4 −a = 0


Py (y = b) = pb
−a 2
Z Z
px3 y − 2c1 xy + c2 y dx = 0


Vx = σx (x, y = 0) dx =
M

Z a Z a
px3 b − 2c1 xb + c2 b dx


Vx = σx (x, y = b) dx =
−a −a
a a
x4
Z 
3
dx = b p − c1 x2 + c2 x


Vx = b px − 2c1 x + c2 = 2bc2 a
−a 4 −a

2
MCE 526 - Theory of Elasticity Dr. M. Sahin
2016-2017 Spring Semester Ankara Yildirim Beyazit University

2. [20] A stress distribution is given as

1
σx = kxy σy = kx σz = νkx (1 + ν) τxy = − ky 2 τyz = τxz = 0 k = constant

17
2
show that this distribution satisfy the plane stress formulation relations

Solution:

20
6
52
3. [20] Show that the Airy stress function
 
q0 1 3 2
Φ(x, y) = 3 x (y −3b y + 2b )− y y −2b
2 3 2 3 2


8b 5

satisfies the compatibility condition. Determine the stress field and find what problem it
corresponds to when applied to the region −b ≤ y ≤ b and x = 0, L
CE

Solution: ∇4 φ = 0 OK

6b2
 
q0 y 2 2 q0 3 2 3

3q0 x 2 2


σxx =− 3 − − 3x + 2y σyy = 2b − 3b y + y τxy = b − y
4b 5 4b3 4b3
M

3
MCE 526 - Theory of Elasticity Dr. M. Sahin
2016-2017 Spring Semester Ankara Yildirim Beyazit University

4. [25] The Airy stress function for cantilever beam subjected free end shear force P is given as
φ (x, y) = Axy + Bxy 3
Determine stress componenents
y

17
P
c
x
c

t=1
L

20
Solution
a) B.C.s. Traction BC.s,

τxy (x, ±c) = 0 (0.1)

Z c
P = τxy (0, y) dy (0.2)
−c

The displacement BC.s u (L, y) = 0 v (L, y) = 0

The stresses
6
52
∂ 2φ ∂ 2φ ∂ 2φ
= σxx = 6Bxy = σyy = 0 − = τxy = −A − 3By 2
∂y 2 ∂x2 ∂y∂x
Apply BCs to the stress equations
• The 1st BC τxy (x, −c) = 0 − A − 3Bc2 = 0

• The 2nd BC τxy (x, c) = 0 − A − 3Bc2 = 0 =⇒ A = −3Bc2

• The 3rd BC
CE

Z c Z c
c
A + 3By 2 dy = Ay + By 3 −c = 2c A + Bc2



P =− τxy (0, y) dy =
−c −c

2Bc3 + 2Ac = +B2c3 − 6Bc3 = P

The solution of the system
P 3P
B=− A=
M

4c3 4c
6P 3P
σxx = 6Bxy = − 3
xy = − 3 xy
4c 2c
1 y2
 
2 3
τxy = −A − 3By = − P − 3
4 c c