桁架结构的计算

桁架结构的计算subroutineestifi(MM,IA1,XY,EKI)!求单刚的子程序REALLIJDIMENSIONEKI(4,4),IA1(5,2),XY(4,2)El=1.0Fl=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-Yl)*(Y2-Y1))0(X2-XI)/LIJS=(Y2-Y1)/LIJF=F1*E1/LIJEK1(1,1)=F*C*CEKI(2,1)=F*C*SEKI(2,2)=F*S*SEK1(3,1)=-F*C*CEK1(3,2)=-F*C*SEKI(3,3)=F*C*CEK1(4,1)=-F*C*SEKI(4,2)=-F*S*SEKI(4,3)=F*S*CEKI(4,4)=F*S*SDO11=1,4DO1J=I+1,41EK1(I,J)=EK1(J,I)ENDSUBROUTINEGS(N,A,X)!高斯消元法DIMENSIONX(N),A(N,N+l)DO10K=1,N-1DO10I=K+1,NDO10J=K+1,N+110A(I,J)=A(I,J)-A(I,K)*A(K,J)/A(K,K)X(N)=A(N,N+1)/A(N,N)DO30K=N-1,1,-1X(K)=0.0DO20尸K+ON20X(K)=X(K)+A(K,J)*X(J)30X(K)=(A(K,N+l)-X(K))/A(K,K)ENDsubroutineXL(MM,IA1,XY,XX,FF)!求内力的子程序REALLIJ,FFDIMENSIONIA1(5,2),XY(4,2),XX(8)El=1.0Fl=1.0X1=XY(I,1)Y1=XY(I,2)J=IA1(MM,2)X2=XY(J,1)Y2二XY(J,2)LIJ=SQRT((X2-X1)*(X2-X1)+(Y2-Y1)*(Y2-Y1))0(X2-XI)/LIJS=(Y2-Y1)/LIJF=F1*E1/LIJU1=XX(2*1-1)V1=XX(2*I)U2=XX(2*JT)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),XI(8),Fl(8),KI(8,8),F2⑸,FODATAIA/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/!节点的初始外载Nl=5N2=8DOM=l,5Ml二MCALLestifi(M,IA,XY1,A)!调用单刚子程序II=IA(M1,1)JJ=IA(Ml,2)K(2*11-1,2*11-1)=K(2*11-1,2*H-1)+A(1,1)K(2*11-1,2*II)=K(2*11-1,2*II)+A(1,2)K(2*11-1,2*JJ-1)=K(2*II-1,2*JJ-1)+A(1,3)K(2*11-1,2*JJ)=K(2*H-1,2*JJ)+A(1,4)K(2*H,2*H-1)=K(2*H,2*II-1)+A(2,1)K(2*II,2*H)=K(2*n,2*11)+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*H,2*JJ)+A(2,4)K(2*JJ-1,2*11-1)=K(2*JJ-1,2*II-1)+A(3,1)K(2*JJ-1,2*11)=K(2*JJ-1,2*H)+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*H)=K(2*JJ,2*11)+A(4,2)K(2*JJ,2*JJ-1)=K(2*JJ,2*JJ-1)+A(4,3)K(2*JJ,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)=1.0ENDIFENDDO!实现1-0的处理过程CALLGS(N2,K,X1)!调用高斯消元的子程序print*，”节点位移是:“DOM=l,8WRITE(*,200)XI(M)ENDDODOM=1,N2DON=1,N2F1(M)=F1(M)+K1(M,N)*X1(N)ENDDOENDDOprint*，"支座反力为:“DOM=1,N2PRINT100,Fl(M)ENDDODOM=1,N1CALLXL(