资源描述：
《桁架结构计算》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、subroutineestifi(MM,IA1,XY,EK1)!求单刚的子程序REALLIJDIMENSIONEK1(4,4),IA1(5,2),XY(4,2)E1=1.0F1=1.0I=IA1(MM,1)X1=XY(I,1)Y1=XY(I,2)I=IA1(MM,2)X2=XY(I,1)Y2=XY(I,2)LIJ=SQRT((X2-X1)*(X2-X1)+(Y2-Y1)*(Y2-Y1))C=(X2-X1)/LIJS=(Y2-Y1)/LIJF=F1*E1/LIJEK1(1,1)=F*C*CEK1(2,1)=F*C*
2、SEK1(2,2)=F*S*SEK1(3,1)=-F*C*CEK1(3,2)=-F*C*SEK1(3,3)=F*C*CEK1(4,1)=-F*C*SEK1(4,2)=-F*S*SEK1(4,3)=F*S*CEK1(4,4)=F*S*SDO1I=1,4DO1J=I+1,41EK1(I,J)=EK1(J,I)ENDSUBROUTINEGS(N,A,X)!高斯消元法DIMENSIONX(N),A(N,N+1)DO10K=1,N-1DO10I=K+1,NDO10J=K+1,N+110A(I,J)=A(I,J)-A(I,K
3、)*A(K,J)/A(K,K)X(N)=A(N,N+1)/A(N,N)DO30K=N-1,1,-1X(K)=0.0DO20J=K+1,N20X(K)=X(K)+A(K,J)*X(J)30X(K)=(A(K,N+1)-X(K))/A(K,K)ENDsubroutineXL(MM,IA1,XY,XX,FF)!求内力的子程序REALLIJ,FFDIMENSIONIA1(5,2),XY(4,2),XX(8)E1=1.0F1=1.0I=IA1(MM,1)X1=XY(I,1)Y1=XY(I,2)J=IA1(MM,2)X2=X
4、Y(J,1)Y2=XY(J,2)LIJ=SQRT((X2-X1)*(X2-X1)+(Y2-Y1)*(Y2-Y1))C=(X2-X1)/LIJS=(Y2-Y1)/LIJF=F1*E1/LIJU1=XX(2*I-1)V1=XX(2*I)U2=XX(2*J-1)V2=XX(2*J)FF=F*(C*(U2-U1)+S*(V2-V1))ENDPROGRAMX!主程序INTEGERM,N,IA(5,2)REALA(4,4),XY1(4,2),K(8,9),U(8),F(8),X1(8),F1(8),K1(8,8),F2(5)
5、,F0DATAIA/1,2,3,1,2,2,3,4,3,4/!单元的构成节点（按列输入结果）DATAXY1/0.0,0.0,5.0,5.0,0.0,5.0,5.0,0.0/!节点坐标DATAU/1,0,1,1,1,1,0,0/!节点位移的初始值DATAF/0,0,10,0,0,-20,0,0/!节点的初始外载N1=5N2=8DOM=1,5M1=MCALLestifi(M,IA,XY1,A)!调用单刚子程序II=IA(M1,1)JJ=IA(M1,2)K(2*II-1,2*II-1)=K(2*II-1,2*II-1)
6、+A(1,1)K(2*II-1,2*II)=K(2*II-1,2*II)+A(1,2)K(2*II-1,2*JJ-1)=K(2*II-1,2*JJ-1)+A(1,3)K(2*II-1,2*JJ)=K(2*II-1,2*JJ)+A(1,4)K(2*II,2*II-1)=K(2*II,2*II-1)+A(2,1)K(2*II,2*II)=K(2*II,2*II)+A(2,2)K(2*II,2*JJ-1)=K(2*II,2*JJ-1)+A(2,3)K(2*II,2*JJ)=K(2*II,2*JJ)+A(2,4)K(2*
7、JJ-1,2*II-1)=K(2*JJ-1,2*II-1)+A(3,1)K(2*JJ-1,2*II)=K(2*JJ-1,2*II)+A(3,2)K(2*JJ-1,2*JJ-1)=K(2*JJ-1,2*JJ-1)+A(3,3)K(2*JJ-1,2*JJ)=K(2*JJ-1,2*JJ)+A(3,4)K(2*JJ,2*II-1)=K(2*JJ,2*II-1)+A(4,1)K(2*JJ,2*II)=K(2*JJ,2*II)+A(4,2)K(2*JJ,2*JJ-1)=K(2*JJ,2*JJ-1)+A(4,3)K(2*JJ,
8、2*JJ)=K(2*JJ,2*JJ)+A(4,4)ENDDODOM=1,N2DON=1,N2K1(M,N)=K(M,N)ENDDOENDDODON=1,N2!给最后一列赋值K(N,N2+1)=F(N)ENDDODON=1,N2!实现1-0的处理过程IF(U(N)<10E-2)THEN!判断U(N)的大小DOM=1,N2K(N,M)=0K(M,N)=0ENDDOK(N,N)
