MATLAB 计算空间桁架

应用MATLAB计算空间桁架结构，有多个子函数，点击第一个即可进行计算

clear
E=200e9;
A1=0.001;
A2=0.002;
A3=0.001;
L1=SpaceTELength（00-4050）;
L2=SpaceTELength（-300050）;
L3=SpaceTELength（004050）;
%1单元刚度阵
thetax1=acos（0/L1）*180/pi
thetay1=atan（4/5）*180/pi
thetaz1=atan（5/4）*180/pi
k1=SpaceTEStiffness（EA1L1thetax1thetay1thetaz1）
%2单元刚度阵
thetax2=acos（3/L2）*180/pi;
thetay2=atan（3/5）*180/pi;
thetaz2=90;
k2=SpaceTEStiffness（EA2L2thetax2thetay2thetaz2）
%3单元刚度阵
thetax3=90;
thetay3=atan（4/5）*180/pi;
thetaz3=180-atan（5/4）*180/pi;
k3=SpaceTEStiffness（EA3L3thetax3thetay3thetaz3）
%组装刚度阵
K=zeros（1212）;
K=SpaceTAssemble（Kk114）
K=SpaceTAssemble（Kk224）
K=SpaceTAssemble（Kk334）

f=[12;0;0]
k=K（10:1210:12）
%计算位移（高斯消去法）
u=k\f;
U=[0;0;0;0;0;0;0;0;0;u];
%节点力
F=K*U
%单元内力
u1=[U（1）;U（2）;U（3）;U（10）;U（11）;U（12）];
u2=[U（4）;U（5）;U（6）;U（10）;U（11）;U（12）];
u3=[U（7）;U（8）;U（9）;U（10）;U（11）;U（12）];
sigma1=SpaceTEStress（EL1thetax1thetay1thetaz1u1）
sigma2=SpaceTEStress（EL2thetax2thetay2thetaz2u2）
sigma3=SpaceTEStress（EL3thetax3thetay3thetaz3u3）


line（[00][04][-40]）;
line（[-30][04][00]）;
line（[00][04][40]）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-03-02 12:42  空间桁架\
     文件         118  2010-05-14 09:59  空间桁架\SpacTELength.m
     文件        1364  2010-05-14 11:02  空间桁架\SpaceTAssemble.asv
     文件        1364  2010-05-14 11:03  空间桁架\SpaceTAssemble.m
     文件         178  2010-05-14 09:50  空间桁架\SpaceTEForce.m
     文件         118  2010-05-14 10:01  空间桁架\SpaceTELength.m
     文件         229  2010-05-14 10:24  空间桁架\SpaceTEStiffness.m
     文件         175  2010-05-14 09:53  空间桁架\SpaceTEStress.m
     文件         459  2010-05-12 19:22  空间桁架\SpingElementStiffness.asv
     文件        1169  2010-05-31 14:24  空间桁架\spacehangjia.m