Finite Element Method18．图示一结构网络,共有 10 个结点,9 个单元和四个支杆,在结点10 作用竖向荷载 10KN,按平面应力问题计算各结点位移和应力.设E=1,μ=0,厚度 t=1*******************************************************************************程序如下：c PLAN STRESS/STRAIN PROBLAM ANALYSIS PROGRAM DIMENSION LND(500,3),X(200),Y(200), &JZ(50,3),PJ(50,3),P(500), &TL(20),AK(500,100),AKE(6,6) OPEN(6,FILE='1.txt',STATUS='OLD') OPEN(8,FILE='2.txt',STATUS='OLD') READ(6,10)TL10 FORMAT(20A2) READ(6,*)NJ,NE,NZ,NPJ,IPS WRITE(8,10)TL IF(IPS.EQ.0)WRITE(8,20) IF(IPS.EQ.1)WRITE(8,30)20 FORMAT(/1X,'PLATE STRESS PROBLEM')30 FORMAT(/1X,'PLATE STRAIN PROBLEM') NPJ0=NPJ IF(NPJ.EQ.0)NPJ=1 CALL INPUT(NJ,NE,NZ,ND,NPJ,NPJ0,IPS,E,PR, &T,V,LND,X,Y,JZ,PJ) N=2*NJ DO 40 I=1,N DO 40 J=1,ND40 AK(I,J)=0.0 DO 50 IE=1,NE CALL MKE(IE,NJ,LND,X,Y,E,PR,T,AKE) DO 50 I=1,3 DO 50 II=1,2 IH=2*(I-1)+II IDH=2*(LND(IE,I)-1)+II DO 50 J=1,3 DO 50 JJ=1,2 L=2*(J-1)+JJ IL=2*(LND(IE,J)-1)+JJ IDL=IL-IDH+1 IF(IDL.LE.0)GOTO 50 AK(IDH,IDL)=AK(IDH,IDL)+AKE(IH,L)50 CONTINUE CALL MF(NJ,NE,NPJ,NPJ0,N,T,V,LND,X,Y,PJ,P) CALL RKR(NZ,ND,N,JZ,AK,P) CALL SLOV(NJ,N,ND,AK,P) CALL MADE(NE,NJ,N,E,PR,LND,X,Y,P) CLOSE(6) CLOSE(8) STOP END SUBROUTINE INPUT(NJ,NE,NZ,ND,NPJ,NPJ0,IPS,E,PR,T,V, &LND,X,Y,JZ,PJ) DIMENSION LND(500,3),X(200),Y(200),JZ(50,3),PJ(50,3) READ(6,*)E,PR,T,V READ(6,*)((LND(I,J),J=1,3),I=1,NE) READ(6,*)(X(I),Y(I),I=1,NJ) READ(6,*)((JZ(I,J),J=1,3),I=1,NZ) IF(NPJ0.EQ.0) GOTO 10 READ(6,*)((PJ(I,J),J=1,3),I=1,NPJ)10 ND=0 DO 20 IE=1,NE DO 20 I=1,3 DO 20 J=1,3 IW=IABS(LND(IE,I)-LND(IE,J)) IF(ND.GE.IW) GOTO 20 ND=IW20 CONTINUE ND=(ND+1)*2 IF(IPS.EQ.0)GOTO 30 E=E/(1.0-PR) PR=PR/(1.0-PR)30 RETURN END SUBROUTINE MKE(IE,NJ,LND,X,Y,E,PR,T,AKE) DIMENSION LND(500,3),X(200),Y(200),B(3,6),D(3,3),S(3,6), &AKE(6,6) CALL MA(IE,NJ,LND,X,Y,AE) CALL MD(E,PR,D) CALL MB(IE,NJ,LND,X,Y,AE,B) DO 10 I=1,3 DO 10 J=1,6 S(I...
