桁架结构的计算.docVIP

subroutine estifi(MM,IA1,XY,EK1) ! 求单刚的子程序 REAL LIJ DIMENSION EK1(4,4),IA1(5,2),XY(4,2) E1=1.0 F1=1.0 I=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)/LIJ S=(Y2-Y1)/LIJ F=F1*E1/LIJ EK1(1,1)=F*C*C EK1(2,1)=F*C*S EK1(2,2)=F*S*S EK1(3,1)=-F*C*C EK1(3,2)=-F*C*S EK1(3,3)=F*C*C EK1(4,1)=-F*C*S EK1(4,2)=-F*S*S EK1(4,3)=F*S*C EK1(4,4)=F*S*S DO 1 I=1,4 DO 1 J=I+1,4 1 EK1(I,J)=EK1(J,I) END SUBROUTINE GS(N,A,X) ! 高斯消元法 DIMENSION X(N),A(N,N+1) DO 10 K=1,N-1 DO 10 I=K+1,N DO 10 J=K+1,N+1 10 A(I,J)=A(I,J)-A(I,K)*A(K,J)/A(K,K) X(N)=A(N,N+1)/A(N,N) DO 30 K=N-1,1,-1 X(K)=0.0 DO 20 J=K+1,N 20 X(K)=X(K)+A(K,J)*X(J) 30 X(K)=(A(K,N+1)-X(K))/A(K,K) END subroutine XL(MM,IA1,XY,XX,FF) !求内力的子程序 REAL LIJ,FF DIMENSION IA1(5,2),XY(4,2),XX(8) E1=1.0 F1=1.0 I=IA1(MM,1) X1=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)) C=(X2-X1)/LIJ S=(Y2-Y1)/LIJ F=F1*E1/LIJ U1=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)) END PROGRAM X !主程序 INTEGER M,N,IA(5,2) REAL A(4,4),XY1(4,2),K(8,9),U(8),F(8),X1(8),F1(8),K1(8,8),F2(5),F0 DATA IA/1,2,3,1,2,2,3,4,3,4/ !单元的构成节点（按列输入结果） DATA XY1/0.0,0.0,5.0,5.0,0.0,5.0,5.0,0.0/ !节点坐标 DATA U/1,0,1,1,1,1,0,0/ !节点位移的初始值 DATA F/0,0,10,0,0,-20,0,0/ !节点的初始外载 N1=5 N2=8 DO M=1,5 M1=M CALL estifi(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)+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*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