Lab – exercise problem 1
v(x=0) = 0 a=0
v(x=L) =0
49047 Finite Element Analysis
1
Solution:
49047 Finite Element Analysis
2
Substitute and perform second derivative, we obtain:
and seek the minimum of the total potential energy (set the first derivative = 0).
Use MATLAB to perform integration and derivative
4*EI*L*c - (L^3*q0)/12 + 6*EI*L^2*d + 8*EI*L^3*e + 10*EI*L^4*f
6*EI*L^2*c - (2*L^4*q0)/15 + 12*EI*L^3*d + 18*EI*L^4*e + 24*EI*L^5*f
8*EI*L^3*c - (L^5*q0)/6 + 18*EI*L^4*d + (144*EI*L^5*e)/5 + 40*EI*L^6*f
10*EI*L^4*c - (4*L^6*q0)/21 + 24*EI*L^5*d + 40*EI*L^6*e + (400*EI*L^7*f)/7
360EI
49047 Finite Element Analysis
4
Lab – exercise problem 2
Note: Ignore friction between the carts and surface
49047 Finite Element Analysis
5
Total internal energy Total external energy
Spring elongation
Displacement
Does not contribute to Spring elongation 2
Applied load
Spring elongation
Solution: Total potential energy
P U
1 1 1
P = k1 D12 k2 D2 D1 k3 D3 D2 P1 D1 P2 D2 P3 D3
2 2
2 2 2
and seek the minimum of the total potential energy.
P P
K1 D1 K 2 D2 D1 P1 =0
D1 D1
P P
K 2 D2 D1 K 3 D3 D2 P2 =0
D2 D2
P P
K 3 D3 D2 P3 =0
D3 D3 What if P3 = 0 ????
49047 Finite Element Analysis
7
Rearrange equations in matrix form: k Di P
2k k 0 D1 P
k
2k k D2 P K1=K2=K3 = k
0 k k D3 P
Solve Di:
3P
D1 =
k
2P
D2 =
k
6P Try to solve by Matlab
D3 =
k
49047 Finite Element Analysis
8