Chapt.11 (Finite Element Analysis)
Chapt.11 (Finite Element Analysis)
Bar element
Model with multiple bar elements
Consistent and lumped-mass matrices
Analysis of trusses
Torsional elements
Beam elements
MATLAB applications
Chapter review
Outline
This chapter introduces the finite element method that provides a more
accurate system description than that used to develop a lumped-parameter
model but that for complex systems is easier to solve than partial differential
equations and more accurate than their finite-difference approximations.
f1 k u1 k u2 , f 2 k u1 k u2
f1 k k u1 EA 1 1 u1
f k
k u2 L 1 1 u2
2
Mass Matrix
d 2u ( x )
Static displacement of bar : EA 0 (0 x L)
dx 2
Solution : u ( x) a bx
Shape functions :
x x
S1 ( x) 1 , S2 ( x)
L L
An approximate solution of EOM of bar
1 L L 1 L
Au12 (t ) S12 ( x)dx Au1 (t )u2 (t ) S1 ( x)S 2 ( x)dx Au22 (t ) S 22 ( x)dx
2 0 0 2 0
Evaluation of integral :
1 1
KE AL[u12 (t ) u1 (t )u2 (t ) u22 (t )] uT (t )Mu(t )
6 2
Velocity vector :
u (t )
u(t ) 1
u2 (t )
This result is dependent on the assumed mode shape Full Ship Model
EA 1 1
K
L 1 1
Shape function used are equivalent to linear force-deflection relation used to derive
stiffness matrix K
So, derivation of K based on potential energy should give same result
Natural frequency :
1 3E 1.7321 E
n
L L
x x
u ( x, t ) 1 u1 (t ) u2 (t )
L L
x
u2 (t )
L
This means 1st row and 1st column of M, K can strike out
Striking out row and column corresponding to the displacement of a fixed node
is a procedure often used as a quick way to reduce the number of equations
If bar has an axial force f(t) applied at its free end, EOM :
AL EA
u2 u2 f (t )
3 L
One element
Two elements
Three elements
Three equations, two unknown u2, u3 , reduce to just two equations by adding Eqs. (1), (2)
m1 m2 m2
u2 u3 (k1 k2 )u2 k2u3 0 (4)
3 3 6
Solve frequencies and mode shapes for case with constant cross section and length
A1 A2 A, L1 L2 L / 2, m1 m2 m / 2, k1 k2 k 2EA / L
FEM predicts 2.6%, 19% higher for 1st and 2nd modes than distributed-parameter model
Both models predict same mode shapes for first two modes
Ei Ai
ki , mi Ai Li
Li
Banded structure due to that a given element is affected only by its adjacent elements
A main diagonal, one upper sub-diagonal, one lower sub-diagonal
1 0 0 1 0 0
AL
0 (1 1) 0 0 2 0
AL
ML
2 2
0 0 1 0 0 2
Boundary conditions : u1 u3 0
2EA
Stiffness matrix : K
L
Lumped-mass matrix : M L AL
2 EA
Equation of motion with lumped-mass matrix : ALu2 u2 0
L
1 2 E 1.41 E
Frequency : (1.41 – 1.57)/1.57
L L = 10.19%
1.57 E
Lowest frequency for distributed-parameter model :
L
Lumped-mass approximation :
Leads to a diagonal mass matrix
More convenient in models having a large number of elements
More easily inverted without excessive numerical error
Introduce modeling error
Choice between use of consistent versus lumped-mass matrices is usually not clear
Stiffness Matrix
GJ
Force-deflection relation of a uniform rod having a net applied torque T : T k
L
GJ EA
Torsional stiffness : k k (Longitudinal vibration)
L L
Mass Matrix
d 2 ( x)
Static displacement of rod : GJ 0 (0 x L)
dx 2
Solution : ( x) a bx
0 2 t 2 0
Evaluation of integral :
1 1 12 (t ) 1 (t )2 (t ) 22 (t ) 1 T
KE JL[1 (t ) 1 (t ) 2 (t ) 2 (t )] JL
2 2
θ (t )Mθ(t )
6 2 3 3 3 2
Velocity vector :
(t )
θ(t ) 1
2 (t )
Mass matrix of bar element :
JL 2 1
M
6 1 2
FEM predicts 10% higher than that by 1st mode of distributed-parameter model
Boundary conditions :
v(0, t ) v1 (t )
v(0, t )
1 (t )
x
v( L, t ) v2 (t )
v( L, t )
2 (t )
x
Shape Functions
Solution of distributed-parameter model approximately :
v( x, t ) a(t ) b(t ) x c(t ) x 2 d (t ) x3
Shape function : Si ( x)
2 3
x x
S1 ( x) 1 3 2
L L
2 3
x x
S2 ( x) x 2 L L
L L
2 3
x x
S3 ( x ) 3 2
L L
2 3
x x
S4 ( x) L L
L L
AL
Evaluation of integral : KE vT Mv
420
Velocity vector : v1
v 1
v2
2
Stiffness Matrix
Similar approach to derive stiffness matrix K from potential energy expression :
2
L v
2
1
PE EI 2 dx vT Kv
2 0
x
I : Area moment of inertia of the section
12 6 L 12 6 L
2
EL 6 L 4 L 6 L 2 L
2
K 3
L 12 6 L 12 6 L
2
6 L 2 L 6 L 4 L
2
Matrix Representation
For future reference, when treat beams having multiple elements
M Mb K Kb
M Ta K Ta
M b M c K b K c
EL 12 6 L EL 12 6 L EL 12 6 L
Ka
L3 6 L 4 L2 L3 6 L 2 L2 L3 6 L 4 L2
, K b , K c
Equation of motion :
AL 4L2 3L2 1 EI 4L2 2L2 1 0
2
420 3L2 4 L2 2 L3 2L 4 L2 2 0
840 EI
, j j eit (2 4 2 )1 (3 2 )2 0
AL 4
(3 2 )1 (2 4 2 )2 0
Frequency equation :
7 4 22 2 3 2 0
2 0.142857 , 2 3
10.9544 EI 50.9544 EI
Natural frequency : 1 ,
L2 A 2
L2 A
For distributed-parameter model :
n
2
EI
n
L A 1st mode :
(10.9544 9.8696)/ 9.8696
9.8696 EI 39.4784 EI
1 , = 10.99%
L2 A 2
L2 A
2nd mode :
(50.9544 39.4784)/ 39.4784
= 29.07%
2
Applying boundary conditions :
Fixed condition : v1 1 0
Equation of motion :
v v 0
Mc1 2 K c1 2
2 2 0
K 3
L 12 3Li 12 3L L3 12 3Li 12 3L
3L L2 / 2 3L L2 3L L / 2 3L
2
L2
Equations of motion :
Forcing terms : f 2 kv2 , f3 P, T2 T3 0
v2 v2 f 2 0
T 0
M K 2 2
2
v3 v3 f3 P
3 3 T3 0
AL 156 22 L EL 12 6 L
M
420 22 L 4 L2 L3 6 L 4 L2
, K
Force : f ( x, t ) w
L x 2 x
3
wL
f 2 w 3 2 dx
0
L L 2
L x 2 x
3
wL2
T2 w L L dx
0
L L 12
Equations of motion :
wL
AL 156 22 L v2 EI 12 6 L v2 2
420 22 L 4 L2 2 L3 6 L 4 L2 2 wL2
12
Need not evaluate integrals for f1, T1 since they do not appear in equations of
motion because of particular boundary conditions