-
Notifications
You must be signed in to change notification settings - Fork 188
/
problem10.1.m
26 lines (26 loc) · 868 Bytes
/
problem10.1.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
E=210e6;
G=84e6;
A=2e-2;
Iy=10e-5;
Iz=20e-5;
J=5e-5;
k1=SpaceFrameElementStiffness(E,G,A,Iy,Iz,J,0,0,0,0,5,0);
k2=SpaceFrameElementStiffness(E,G,A,Iy,Iz,J,0,0,4,0,5,4);
k3=SpaceFrameElementStiffness(E,G,A,Iy,Iz,J,4,0,4,4,5,4);
k4=SpaceFrameElementStiffness(E,G,A,Iy,Iz,J,4,0,0,4,5,0);
k5=SpaceFrameElementStiffness(E,G,A,Iy,Iz,J,0,5,0,0,5,4);
k6=SpaceFrameElementStiffness(E,G,A,Iy,Iz,J,0,5,4,4,5,4);
k7=SpaceFrameElementStiffness(E,G,A,Iy,Iz,J,4,5,4,4,5,0);
k8=SpaceFrameElementStiffness(E,G,A,Iy,Iz,J,0,5,0,4,5,0);
K=zeros(48,48);
K=SpaceFrameAssemble(K,k1,1,5);
K=SpaceFrameAssemble(K,k2,2,6);
K=SpaceFrameAssemble(K,k3,3,7);
K=SpaceFrameAssemble(K,k4,4,8);
K=SpaceFrameAssemble(K,k5,5,6);
K=SpaceFrameAssemble(K,k6,6,7);
K=SpaceFrameAssemble(K,k7,7,8);
K=SpaceFrameAssemble(K,k8,5,8);
k=K(25:48,25:48);
f=[0;0;0;0;0;0;0;0;0;0;0;0;-15;0;0;0;0;0;0;0;0;0;0;0];
u=k\f