下题是桁架结构：
其中，F p =1，试求各杆的轴力
clear
[ length1,theta1,T1 ] = Geometry( 1,0, 1,1) % 由结点坐标求六个单元几何信息 [ length2,theta2,T2 ] = Geometry( 0,0, 1,0) [ length3,theta3,T3 ] = Geometry( 0,0, 0,1) [ length4,theta4,T4 ] = Geometry( 0,1, 1,1) [ length5,theta5,T5 ] = Geometry( 1,0, 0,1) [ length6,theta6,T6 ] = Geometry( 0,0, 1,1) g1=[1 2 3 4] % 六个单元定位向量
g2=[0 0 1 2] g3=[0 0 0 0] g4=[0 0 3 4] g5=[1 2 0 0] g6=[0 0 3 4]
丢手绢E=1;A=1; % 材料信息 k1=LocalElementStiff(E,A,length1) % 局部单刚 k2=LocalElementStiff(E,A,length2) k3=LocalElementStiff(E,A,length3) k4=LocalElementStiff(E,A,length4) k5=LocalElementStiff(E,A,length5) k6=LocalElementStiff(E,A,length6)
K1=T1'*k1*T1 % 整体单刚
右手戒指的戴法和意义K2=T2'*k2*T2 K3=T3'*k3*T3
（6）
3
2
1
4
K4=T4'*k4*T4
K5=T5'*k5*T5
K6=T6'*k6*T6
GK=zeros(4) % 由六个单元和相应定位向量聚集总刚GK=ElementAsmble1(GK,K1,g1)
GK=ElementAsmble1(GK,K2,g2)
GK=ElementAsmble1(GK,K3,g3)
GK=ElementAsmble1(GK,K4,g4)
GK=ElementAsmble1(GK,K5,g5)
GK=ElementAsmble1(GK,K6,g6)
P=[1 0 0 1]' % 荷载信息
U=GK\P % 求解方程
家常素菜100道素菜U1=ElementDisp(U,g1) % 定位向量提取六个单元的结点位移
U2=ElementDisp(U,g2)
把我变成坏孩子
U3=ElementDisp(U,g3)描写草的诗句
U4=ElementDisp(U,g4)
U5=ElementDisp(U,g5)
U6=ElementDisp(U,g6)
f1=k1*T1*U1 % 求局部坐标下的单元杆端内力
f2=k2*T2*U2
f3=k3*T3*U3
f4=k4*T4*U4
f5=k5*T5*U5
f6=k6*T6*U6
下题是钢架结构：
（4,0,5）
（1,2,3）
（0,0,0）
clear
[ length1,theta1,T1 ] = Geometry( 0,0, 0,10) % 由结点坐标求两个单元几何信息[ length2,theta2,T2 ] = Geometry( 0,0, 10,0)
g1=[1 2 3 0 0 0] % 两个单元定位向量
g2=[1 2 3 4 0 5]
E=128000;A=2;I=7/640 % 材料信息
k1=LocalElementStiff(E,A,I,length1) % 局部单刚
k2=LocalElementStiff(E,A,I,length2)
K1=T1'*k1*T1
K2=T2'*k2*T2
GK=zeros(5) % 由两个单元和相应定位向量聚集总刚GK=ElementAsmble1(GK,K1,g1)
GK=ElementAsmble1(GK,K2,g2)
P=[10 0 0 0 20]' % 荷载信息
U=GK\P % 求解方程
U1=ElementDisp(U,g1) % 定位向量提取两个单元的结点位移
U2=ElementDisp(U,g2)
f1=k1*T1*U1 % 求局部坐标下的单元杆端内力
f2=k2*T2*U2软件测试总结
故M图：
Q图：
N图：
下题是两个连续梁单元的结构
clear
length1 = Geometry( 0,0, 6,0) % 由结点坐标求2个单元几何信息length2 = Geometry( 6,0, 12,0)
g1=[0 1] % 2个单元定位向量
g2=[1 2]
E=1,I=1 % 材料信息
k1=LocalElementStiff(E,I,length1) % 局部单刚
k2=LocalElementStiff(E,I,length2)
K1=k1 % 整体单刚
K2=k2
GK=zeros(2) % 由2个单元和相应定位向量聚集总刚
GK=ElementAsmble1(GK,K1,g1)饱满的意思
GK=ElementAsmble1(GK,K2,g2)
P=[50 0]' % 荷载信息
三十六计故事
U=GK\P % 求解方程
U1=ElementDisp(U,g1) % 定位向量提取2个单元的整体坐标结点位移U2=ElementDisp(U,g2)
u1=U1 % 局部坐标结点位移
u2=U2
f1=k1*u1 % 求局部坐标下的单元杆端内力
f2=k2*u2
故M图为：
