AE 5310-003,004
Finite Element method Prof. B.P. Wang
Final Examination
DUE: Tuesday, December 8, 2016 at the beginning of class
Name:____________,CID:_______________
[Link]
Rules:
(1) Open books and notes
(2) Work independently
(3) Summarize your answers in the answer sheet
(4) Put all work in aSolution Book with page numbers
(5) Include all long hand calculation in the solution book
(6) Do not submit excessive Matlab outputs. For MASTRAN solutions,
just hand in (a) input data (b) output displacement and element
force, reactions (c) a plot with model and deformation
1
Name:____________,CID:_______________
Answer Sheet (Report any value less than 10-12 as zero)
Problem 1(20%) (pages ___ to ____)
(1) (6%)Largest A1=
U3x= ___________, U3y=______________
(2) (12%) U3x= ___________, U3y=______________
(3) (2%) U3x= ___________, U3y=______________
Problem 2(16%) (pages ___ to ____)
V(free end) = ______________
Problem 3 (24%)(pages ___ to ____)
(a) Nodes N4 N7 N3
4
7
3
(b) Along 4-7-3
N4(s,t) =
N7(s,t) =
N3(s,t) =
(c) J=Jacobian matrix = [ ]
[ ]
det(J)=
2
Problem 4(16) (pages ___ to ____)
(a) P: (X,Y)=( , )
(b) Nbar=[N1 n2 N3 N4 ]
=[ , , , ]
Sum =
(c) J=Jacobian Matrix =[ ]
[ ]
(d) J(at point P) = [ ]
[ ]
Problem 5(24%) (pages ___ to ____)
(a) Displacement at P
Use 1 Q4 element = [ , ]
Use 2 triangular element (let these two elements be element a
and b)
Element a solution =[ , ]
Element b solution =[ , ]
(b) Stress at P
Use 1 Q4 element = [ , . ]
Use 2 triangular element (let these two elements be element a
and b)
Element a solution =[ , , ]
Element b solution =[ , , ]
(c) Z-strain at P
Use 1 Q4 element = [ ]
3
Use 2 triangular element (let these two elements be element a
and b)
Element a solution =[ ]
Element b solution =[ ]
Problem 6(20%) (pages ___ to ____)
(a) Displacement responses at A
Ux Uy theta
Long Hand Solution
MASTRAN Solution
Use _______ elements.
(b) Displacement responses at P
Ux Uy theta
Long Hand Solution
MASTRAN Solution
(c) Reactions at B
Rx Ry Mz
Long Hand Solution
MASTRAN Solution
Reactions at C
Rx Ry Mz
Long Hand Solution
MASTRAN Solution
4
Problem 1 (20%) Two bar truss: long hand solutions
For the 2-bar truss shown below,
(a) Find the displacement at node 3 by finite element method use
long hand solution. Try several large values of A1 and report
the largest A1 you can get a solution.
(b) Solve the problem by treating element 1 as a rigid element use
MPC.
(c) Validate the nodal displacement by the formula given in the
Appendix.
E = 107 psi
Y
2 L = CID*10 in.
A1 ->Inf, A2=X(1)
(areas in in^2)
H H = X(2) in.
Py Px = 10000 lb.
Py = 2000 lb.
A1, E
X
3 Px
1 L
Use X=[ 1 CID*5]
Py = 2000 lb.
5
Appendix: Truss Deflection Formula
()
=
Where
Uj=displacement of dof j
Le=length of truss element
Fe=force in element e due to applied load
fe(j) = force in element e due to a unit load applied in dof j, no other
loads
Ae = cross sectional area of element e
E = Youngs modulus
6
Problem 2 (16%) Two span beam one element solution
Find the deflection at the free end of the beam shown in Figure
2 use finite element method. Model the system by 1 beam element.
Use the following data:
a=5;b=2*CID+10;c=CID+20 all in inch; EI=100000 lb-in^2 ;M=2000
in-lb. (Long hand Solution)
a b c
A B M C
Figure 2
7
Problem 3 (24%) Q8 element
For the following 8-noded Q8 element
(a) Calculate the values of N4, N7 and N3 at nodes 4, 7 and
3.
(b) Find analytical form of the shape function along the
edge 4-7-3. Use these equations to calculate N4, N7 and
N3 at nodes 4, 7 and 3.
(c) Find the Jacobian matrix at node 3 and find its
determinant.
Data for the Q8 element
Node Coordinates
Number x y s t
1 0 0 -1 -1
2 20 0 1 -1
3 40 40 1 1
4 X4 35 -1 1
5 10 0 0 -1
6 35 30 1 0
5 X7 Y7 0 1
8 X8 Y8 -1 0];
Where X4=-.1*CID-4;
Nodes 7, 8 are mid-side nodes (at the
center of the sides)
8
Problem 4 (16%): Q4 Element (Long hand solution)
For the Q4 element shown in Figure 4, the nodal coordinates are:
0 10
100 0
110 150
-1.5*CID 100
(a) Find the (x,y) coordinates of point P which is the
intersection point of line 1-3 and 2-4.
(b) Find the shape function matrix Nbar at point P
What is the sum of elements in Nbar(i.e, N1+N2+N3+N4 = ?)
(c) Find the JAC (Jacobian matrix) for this element (in terms of
s,t)
(d) Evaluate JAC at point P.
9
Problem 5 (24%)
For the quadrilateral plate shown in Figure 5, the coordinates
of A to D are:
0 10
100 0
110 150
-1.5*CID 100
Find the displacement and stresses and z-component of the
strain vector at point P (the intersection point of line 1-3 and
2-4.) by finite element method use the following models
(a) A Q-4 element
(b) Two triangular elements.
The plate has a thickness of 0.5 in, E=30e7 psi with loading
Px=2000 lb and Py=5000 lb at point C.
The plate has the following BC's: at point A, Ux=Uy=0, at point
B, Ux=0.
Solve the problem by longhand calculation use ISO-Tools.
10
Problem 6(20%): 2-member frame
For the 2-member frame shown below, find the displacement
response at points A and P (mid-point between B and A) and
reactions at the clamped points at B and C. Solve the problem
use finite element method. (a) By long hand calculation, model
the structure by 2 elements. (b) Solve the problem by MASTRAN.
= 37 ()
Circular cross section with d=3 (in) for both members.
A Z0
a2
a4
a1
C
a3
Figure 6 2-member frame
(Note: the figure is not to scale)
11
Dimension data:
b= 5 + CID
a1 = 4b
a2 = 3b
a3 = 3b
a4 = 4b
Loading data:
Uniform downward load along BA with 5 (lb/in)
Horizontal load (Linearly varying) along AC with 20 (lb/in) to 5
(lb/in)
12
13
