Solution Homework4
Solution Homework4
Solutions
Problem 1
Consider a one-element triangular mesh shown in Fig. 1. The boundary
conditions are as follows. The edge BC is constrained in y and traction free in x,
whereas the edge AB is constrained in x and traction free in y. The edge AC is
subject to traction normal to the edge as shown in Fig. 1. Assume Young’s
modulus E= 3 × 10 11 Pa and Poisson’s ratio v=0.3.
a. Construct the stiffness matrix.
b. Calculate the global force matrix.
c. Solve for the unknown displacement matrix and calculate the stress at (1.5,
1.5).
3 (0,3)
(1)
1 (0,0) 2 (3,0)
Solution
We will use the plane strain and plane stress conditions to solve this problem. As
shown in Fig.1, x 1(1 ) = 0 , y1( 1) = 0 , x2(1 ) = 3 , y2(1 ) = 0 , x 3(1 ) = 0 , and y3(1 ) = 3 .
(1 + v)(1 − 2v)
0 0 (1 − 2v) / 2 0 0 0.2
0
0
0.2600
0.7799 (1 ) −10
; ε (1.5,1.5 ) = B d = 10 0.2600
(1 ) −10 (1)
Thus, d = 10
0 0
0
0.7799
15
σ (1 )
(1.5,1.5) = D (1)
ε (1 )
(1.5,1.5) = 15 .
0
Problem 2
For the two-dimensional loaded plate shown in Fig. 2, determine the
displacement of nodes 1 and 2 and the element stresses using plane stress
conditions. Body force may be neglected in comparison with the external forces.
3 2
e=2
2m
e=1
Thickness t=0.5m
E=30*106Pa, ν=0.25 1
4 x
3m
Fig. 2
Solution
As shown in Fig.2, for elment 1 (node 1, 2, 3), x 1(1) = 3 , y1( 1) = 0 , x 2(1) = 3 , y2(1) = 2 ,
x 3(1 ) = 0 , and y3(1) = 0 . For element 2 (node 3, 4, 2), x 1( 2 ) = 0 , y1( 2 ) = 2 , x 2( 2 ) = 0 ,
y2( 2 ) = 0 , x 3( 2 ) = 3 , and y3( 2 ) = 2 .
1 v 0 3.2 0.8 0
E
D (1 )
=D (2)
= 2
v 1 0 = 10 0.8 3.2 0
7
1−v
0 0 (1 − v) / 2 0 0 1.2
0
r
1y
0
− 1000
The force vector is f Γ +r=
r3 x .
r3 y
r
4x
r4 y
u1 x
Thus d E = 0 , d F = u2 x ,
u2 y
1x 2x 2y
0.983 − 0.45 0.2 1 x
K F = 10 − 0.45 0.983 0 2 x ,
7
0.2 0 1.4 2 y
0
and f F = 0 .
− 1000
For element 1,
1.913
0
0.638 - 93.3
−5
0. 875 (1) −5
d = 10 , ε = B d = 10 − 3.718 , σ = D ε = - 1138.7 .
(1 ) (1) (1) (1 ) (1) ( 1 )
For element 2,
0
0
0.292 93.4
0
d ( 2 ) = 10 −5 ( 2) (2) (2) −5
, ε = B d = 10 0 , σ = D ε = 23.4 .
(1) (1) ( 1 )