大作業(yè)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)+ 楊交通學(xué)件的受力問題（力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(輸入上桁架節(jié)點處的豎向荷載P:);for j=1:1:n:4*n)K(4*n+3:4*n+4,4*n+3:4*n+4); forb=input輸入單元右端節(jié)點編號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第一題（算例驗證輸入彈性模量輸入截面積輸入節(jié)點總數(shù)n:3輸入x1 的值: 0 輸入y1 的值: 0 輸入x2 的值: 5 輸入y2 的值: 0 輸入x3 的值: 5 輸入y3 的值: -5輸入荷載 f:10KN中間節(jié)點位移單元力00第二題（算例驗證第二題（算例驗證輸入彈性模量桁架下部編號依次為 1、3、5、7桁架上部編號依次為 2、4、6、8輸入截面積A:0.16m2輸入桁架全長ls:20m輸入分段個數(shù)n:20m輸入桁架高ly:1輸入上桁架節(jié)點處的豎向荷載輸入單

7、元右端節(jié)點編號附錄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 集成總剛的函數(shù)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