Solutions Manual
To Accompany the Book
MATLAB Guide to Finite Elements
By
Peter I. Kattan
Problem 2.1:
» k1=SpringElementStiffness(200)
k1 =
200 -200
-200 200
» k2=SpringElementStiffness(250)
k2 =
250 -250
-250 250
» K=zeros(3,3)
K =
0 0 0
0 0 0
0 0 0
» K=SpringAssemble(K,k1,1,2)
K =
200 -200 0
-200 200 0
0 0 0
» K=SpringAssemble(K,k2,2,3)
K =
200 -200 0
-200 450 -250
0 -250 250
» k=K(2,2)
k =
450
» f=[10]
f =
10
» u=k\f
u =
0.0222
» U=[0 ; u ; 0]
U =
0
0.0222
0
» F=K*U
F =
-4.4444
10.0000
-5.5556
» u1=[0;u]
u1 =
0
0.0222
» f1=SpringElementForces(k1,u1)
f1 =
-4.4444
4.4444
» u2=[u ; 0]
u2 =
0.0222
0
» f2=SpringElementForces(k2,u2)
f2 =
5.5556
-5.5556
Problem 2.2:
» k1=SpringElementStiffness(170)
k1 =
170 -170
-170 170
» k2=SpringElementStiffness(170)
k2 =
170 -170
-170 170
» k3=SpringElementStiffness(170)
k3 =
170 -170
-170 170
» k4=SpringElementStiffness(170)
k4 =
170 -170
-170 170
» K=zeros(4,4)
K =
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
» K=SpringAssemble(K,k1,1,2)
K =
170 -170 0 0
-170 170 0 0
0 0 0 0
0 0 0 0
» K=SpringAssemble(K,k2,2,3)
K =
170 -170 0 0
-170 340 -170 0
0 -170 170 0
0 0 0 0
» K=SpringAssemble(K,k3,2,3)
K =
170 -170 0 0
-170 510 -340 0
0 -340 340 0
0 0 0 0
» K=SpringAssemble(K,k4,3,4)
K =
170 -170 0 0
-170 510 -340 0
0 -340 510 -170
0 0 -170 170
» k=K(2:4,2:4)
k =
510 -340 0
-340 510 -170
0 -170 170
» f=[0 ; 0 ; 25]
f =
0
0
25
» u=k\f
u =
0.1471
0.2206
0.3676
» U=[0;u]
U =
0
0.1471
0.2206
0.3676
» F=K*U
F =
-25.0000
0.0000
0.0000
25.0000
» u1=[0;U(2)]
u1 =
0
0.1471
» f2=SpringElementForces(k1,u1)
f2 =
-25.0000
25.0000
» u2=[U(2);U(3)]
u2 =
0.1471
0.2206
» f2=SpringElementForces(k2,u2)
- 1
- 2
前往页