大作业matlab编程-ypsank truss第一题本程序可以解决任意二元杆

1、K(2*i-1,2*i)=K(2*i-1,2*i)+K(2*i-1,2*j-1)=K(2*i-1,2*i)=K(2*i-1,2*i)+K(2*i-1,2*j-1)=K(2*i-1,2*j-1)+K(2*i-1,2*j)=K(2*i-1,2*j)+ 杨交通学件的受力问题（力E=input(E:); A=input(A:); for i=1:1:nstr=xnum2str(i的值 str=ynum2str(i的值 K(2*i,2*i-1)=K(2*i,2*i-1)+K(2*i,2*i)=K(2*i,2*i)+K(2*i,2*j-1)=K(2*i,2*j-1)+ K(2*i,2*j)=K(2*i,2

2、*j)+K(2*j-1,2*i-1)=K(2*j-1,2*i-1)+K(2*j-1,2*i)=K(2*j-1,2*i)+K(2*j-1,2*j-1)=K(2*j-1,2*j-1)+K(2*j-1,2*j)=K(2*j-1,2*j)+ K(2*j,2*i-1)=K(2*j,2*i-1)+ K(2*j,2*i)=K(2*j,2*i)+K(2*j,2*j-1)=K(2*j,2*j-1)+ K(2*j,2*j)=K(2*j,2*j)+k1(4,4); K(2*m-1,2*m-1)=K(2*m-1,2*m-1)+ K(2*m-1,2*m)=K(2*m-1,2*m)+ K(2*m-1,2*n-1)=K(2

3、*m-1,2*n-1)+ K(2*m-1,2*n)=K(2*m-1,2*n)+K(2*m,2*m-1)=K(2*m,2*m-1)+ K(2*m,2*m)=K(2*m,2*m)+K(2*m,2*n-1)=K(2*m,2*n-1)+K(2*m,2*n)=K(2*m,2*n)+K(2*n-1,2*m-1)=K(2*n-1,2*m-1)+ K(2*n-1,2*m)=K(2*n-1,2*m)+K(2*n-1,2*n-1)=K(2*n-1,2*n-1)+ K(2*n-1,2*n)=K(2*n-1,2*n)+K(2*n,2*m-1)= K(2*n,2*m-1)+K(2*n,2*m)=K(2*n,2*m)+K

4、(2*n,2*n-1)=K(2*n,2*n-1)+K(2*n,2*n)=K(2*n,2*n)+ C1=cos(a1); S1= sin(a1); C2=cos(a2); S2= sin(a2);k1=E*A/l1*C1*C1C1*S1 -C1*S1;C1*S1S1*S1-S1*S1;-C1*C1-C1*S1C1*C1C1*S1-C1*S1-S1*S1C1*S1k2=E*A/l2*C2*C2C2*S2 -C2*S2;C2*S2S2*S2-S2*S2;-C2*C2-C2*S2C2*C2C2*S2-C2*S2-S2*S2C2*S2K(2*i-1,2*i-1)=K(2*i-1,2*i-1)+E=inp

5、ut(E:); A=input(A:); n=input(n:); ly=input(ly:);p=input(输入上桁架节点处的竖向荷载P:);for j=1:1:n:4*n)K(4*n+3:4*n+4,4*n+3:4*n+4); forb=input输入单元右端节点编号ifif p=input( ifp=0; p=input( ifp=0; PCG if nargin = 4eps=1.0e-r1=b-A*x0; iM=inv(M); z1 = iM*r1; = z1;=tol=while=x=x0+ r2=r1-alpha*A*p; z2 = iM*r2;=dot(r2,z2)/dot(

6、r1,z1); = z2+belta*p;=n+ tol=norm(x-x0); x0 = x;%r1=r2; z1=1.0e+003第一题（算例验证输入弹性模量输入截面积输入节点总数n:3输入x1 的值: 0 输入y1 的值: 0 输入x2 的值: 5 输入y2 的值: 0 输入x3 的值: 5 输入y3 的值: -5输入荷载 f:10KN中间节点位移单元力00第二题（算例验证第二题（算例验证输入弹性模量桁架下部编号依次为 1、3、5、7桁架上部编号依次为 2、4、6、8输入截面积A:0.16m2输入桁架全长ls:20m输入分段个数n:20m输入桁架高ly:1输入上桁架节点处的竖向荷载输入单

7、元右端节点编号附录function y = dg(E,A,L,a) C = cos(a);S=y = E*A/L*C*C C*S -C*C -C*S ; C*S -C*S-S*S-C*C -C*S C*C C*S ; -C*S -S*S C*S 集成总刚的函数functiony=K(2*i-1,2*i-1)=K(2*i-1,2*i-1)+K(2*i-1,2*i)=K(2*i-1,2*i)+K(2*i-1,2*j-1)=K(2*i-1,2*j-1)+K(2*i-1,2*j)=K(2*i-1,2*j)+ K(2*i,2*i-1)=K(2*i,2*i-1)+K(2*i,2*i)=K(2*i,2*i)+K(2*i,2*j-1)=K(2*i,2*j-1)+ K(2*j,2*j)=K(2*j,2*j)+K(2*i,2*j)=K(2*i,2*j)+K(2*j-1,2*i-1)=K(2*j-1,2*i-1)+K(2*j-1,2*i)=K(2*j-1,2*i)+K(2*j-1,2*j-1)=K(2*j-1,2*j-1)+K(2*j-1,2*j)=K(2*j-1