451006
Linear Bar Element
E=210e6;
L=100/10;
A1=(13+12)/2;
A2=(12+11)/2;
A3=(11+10)/2;
A4=(10+9)/2;
A5=(9+8)/2;
A6=(8+7)/2;
A7=(7+6)/2;
A8=(6+5)/2;
A9=(5+4)/2;
A10=(4+3)/2;
k1=LinearBarElementStiffness(E,A1,L)
k2=LinearBarElementStiffness(E,A2,L)
k3=LinearBarElementStiffness(E,A3,L)
k4=LinearBarElementStiffness(E,A4,L)
k5=LinearBarElementStiffness(E,A5,L)
k6=LinearBarElementStiffness(E,A6,L)
k7=LinearBarElementStiffness(E,A7,L)
k8=LinearBarElementStiffness(E,A8,L)
k9=LinearBarElementStiffness(E,A9,L)
k10=LinearBarElementStiffness(E,A10,L)
K=zeros(11,11);
K=LinearBarAssemble(K,k1,1,2);
K=LinearBarAssemble(K,k2,2,3);
K=LinearBarAssemble(K,k3,3,4);
K=LinearBarAssemble(K,k4,4,5);
K=LinearBarAssemble(K,k5,5,6);
K=LinearBarAssemble(K,k6,6,7);
K=LinearBarAssemble(K,k7,7,8);
K=LinearBarAssemble(K,k8,8,9);
K=LinearBarAssemble(K,k9,9,10);
K=LinearBarAssemble(K,k10,10,11);
k=K(2:11,2:11)
f=[0;0;0;0;0;0;0;0;0;50]
u=k\f
U=[0 ;u(1:10)]
u1=[0 ; U(1)]
u2=[U(1) ; U(2) ]
u3=[U(2) ; U(3) ]
u4=[U(3) ; U(4) ]
u5=[U(4) ; U(5) ]
u6=[U(5) ; U(6) ]
u7=[U(6) ; U(7) ]
u8=[U(7) ; U(8) ]
u9=[U(8) ; U(9) ]
u10=[U(9) ; U(10) ]
sigma1=BarElementStress(E,L,u1)
sigma2=BarElementStress(E,L,u2)
451006
sigma3=BarElementStress(E,L,u3)
sigma4=BarElementStress(E,L,u4)
sigma5=BarElementStress(E,L,u5)
sigma6=BarElementStress(E,L,u6)
sigma7=BarElementStress(E,L,u7)
sigma8=BarElementStress(E,L,u8)
sigma9=BarElementStress(E,L,u9)
sigma10=BarElementStress(E,L,u10)
strain1=BarElementStrain(L,u1)
strain2=BarElementStrain(L,u2)
strain3=BarElementStrain(L,u3)
strain4=BarElementStrain(L,u4)
strain5=BarElementStrain(L,u5)
strain6=BarElementStrain(L,u6)
strain7=BarElementStrain(L,u7)
strain8=BarElementStrain(L,u8)
strain9=BarElementStrain(L,u9)
strain10=BarElementStrain(L,u10)
k1 =
262500000 -262500000
-262500000 262500000
k2 =
241500000 -241500000
-241500000 241500000
k3 =
220500000 -220500000
-220500000 220500000
k4 =
199500000 -199500000
-199500000 199500000
k5 =
178500000 -178500000
-178500000 178500000
k6 =
451006
157500000 -157500000
-157500000 157500000
k7 =
136500000 -136500000
-136500000 136500000
k8 =
115500000 -115500000
-115500000 115500000
k9 =
94500000 -94500000
-94500000 94500000
k10 =
73500000 -73500000
-73500000 73500000
k =
Columns 1 through 6
504000000 -241500000 0 0 0 0
-241500000 462000000 -220500000 0 0 0
0 -220500000 420000000 -199500000 0 0
0 0 -199500000 378000000 -178500000 0
0 0 0 -178500000 336000000 -157500000
0 0 0 0 -157500000 294000000
0 0 0 0 0 -136500000
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
Columns 7 through 10
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
-136500000 0 0 0
252000000 -115500000 0 0
-115500000 210000000 -94500000 0
0 -94500000 168000000 -73500000
451006
0 0 -73500000 73500000
f =
0
0
0
0
0
0
0
0
0
50
u =
1.0e-05 *
0.0190
0.0398
0.0624
0.0875
0.1155
0.1472
0.1839
0.2272
0.2801
0.3481
U =
1.0e-05 *
0
0.0190
0.0398
0.0624
0.0875
0.1155
0.1472
0.1839
0.2272
0.2801
0.3481
u1 =
0
0
451006
u2 =
1.0e-06 *
0
0.1905
u3 =
1.0e-06 *
0.1905
0.3975
u4 =
1.0e-06 *
0.3975
0.6243
u5 =
1.0e-06 *
0.6243
0.8749
u6 =
1.0e-05 *
0.0875
0.1155
u7 =
1.0e-05 *
0.1155
0.1472
u8 =
1.0e-05 *
451006
0.1472
0.1839
u9 =
1.0e-05 *
0.1839
0.2272
u10 =
1.0e-05 *
0.2272
0.2801
sigma1 =
sigma2 =
4.0000
sigma3 =
4.3478
sigma4 =
4.7619
sigma5 =
5.2632
sigma6 =
5.8824
sigma7 =
6.6667
451006
sigma8 =
7.6923
sigma9 =
9.0909
sigma10 =
11.1111
strain1 =
strain2 =
1.9048e-08
strain3 =
2.0704e-08
strain4 =
2.2676e-08
strain5 =
2.5063e-08
strain6 =
2.8011e-08
strain7 =
3.1746e-08
strain8 =
3.6630e-08
451006
strain9 =
4.3290e-08
strain10 =
5.2910e-08
Published with MATLAB® R2012b
