平面桁架静力计算程序
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
E=2.1e11;
A=1E-4;
L1=1;L2=PlaneTrussElementLength(1,1,2,0);L3=1;L4=1;
L5=PlaneTrussElementLength(0,1,1,0);L6=1;L7=1;
theta1=0;theta2=-45;theta3=180;theta4=-90;theta5=-45;theta6=-90;theta7=180;
k1=PlaneTrussElementStiffness(E,A,L1,theta1);
k2=PlaneTrussElementStiffness(E,A,L2,theta2);
k3=PlaneTrussElementStiffness(E,A,L3,theta3);
k4=PlaneTrussElementStiffness(E,A,L4,theta4);
k5=PlaneTrussElementStiffness(E,A,L5,theta5);
k6=PlaneTrussElementStiffness(E,A,L6,theta6);
k7=PlaneTrussElementStiffness(E,A,L7,theta7);
K=zeros(10,10);
K=PlaneTrussAssemble(K,k1,1,2);
K=PlaneTrussAssemble(K,k2,2,3);
K=PlaneTrussAssemble(K,k3,3,4);
K=PlaneTrussAssemble(K,k4,2,4);
K=PlaneTrussAssemble(K,k5,1,4);
K=PlaneTrussAssemble(K,k6,1,5);
K=PlaneTrussAssemble(K,k7,4,5);
K1 = [K(1:5,1:5) K(1:5,7:8);K(7:8,1:5) K(7:8,7:8)];
f=[0;0;0;-200000;0;0;0];
U1=K1\f;
U=[U1(1:5);0;U1(6:7);0;0];
F=K*U;
u1=[U(1);U(2);U(3);U(4)];
u2=[U(3);U(4);U(5);U(6)];
u3=[U(5);U(6);U(7);U(8)];
u4=[U(3);U(4);U(7);U(8)];
u5=[U(1);U(2);U(7);U(8)];
u6=[U(1);U(2);U(9);U(10)];
u7=[U(7);U(8);U(9);U(10)];
sigma1=PlaneTrussElementStress(E,L1,theta1,u1);
sigma2=PlaneTrussElementStress(E,L2,theta2,u2);
sigma3=PlaneTrussElementStress(E,L3,theta3,u3);
sigma4=PlaneTrussElementStress(E,L4,theta4,u4);
sigma5=PlaneTrussElementStress(E,L5,theta5,u5);
sigma6=PlaneTrussElementStress(E,L6,theta6,u6);
sigma7=PlaneTrussElementStress(E,L7,theta7,u7);