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