WRITE(6,10) E,PR,T,V

10 FORMAT(/6X,'E',10X,'PR',9X,'T',9X,'V'/,4F10.2) READ(5,*)((LND(I,J),J=1,3),I=1,NE)（结点编码） WRITE(6,20)

20 FORMAT(/1X,'ELEMENT INFORMATION'/3X,'ELEM',3X, * 'I J K'/)

WRITE(6,30)(I,(LND(I,J),J=1,3),I=1,NE)

30 FORMAT(1X,4I5)

READ(5,*)(X(I),Y(I),I=1,NJ)（结点坐标）

WRITE(6,40)

40 FORMAT(/1X,'COORDINATES OF NODES'/3X,'NODES', * 8X,'X',13X,'Y')

WRITE(6,50)(I,X(I),Y(I),I=1,NJ)

50 FORMAT(1X,I5,2E15.6)

READ(5,*)((JR(I,J),J=1,3),I=1,NS)（约束信息） WRITE(6,60)

60 FORMAT(/1X,'CONSTRAINED NODES'/3X,'NODE',3X,'X',4X,'Y') WRITE(6,70)((JR(I,J),J=1,3),I=1,NS)

70 FORMAT(1X,3I5)

READ(5,*)((PJ(I,J),J=1,3),I=1,NPJ)（荷载信息） WRITE(6,80)

80 FORMAT(/1X,'LOAD CASES'/3X,'NODE',8X,'X',13X,'Y') WRITE(6,90)((PJ(I,J),J=1,3),I=1,NPJ)

90 FORMAT(1X,F5.0,2E15.6)

100 NW=0（半带宽）

DO 110 IE=1,NE

DO 110 I=1,3

DO 110 J=1,3

IW=IABS(LND(IE,I)-LND(IE,J))

IF(NW.LT.IW) THEN

NW=IW

ENDIF

110 CONTINUE

NW=(NW+1)*2

IF(IPS.NE.0) THEN

E=E/(1.0-PR*PR)

PR=PR/(1.0-PR)

ENDIF

END

*---------------------------------------------------------

C SUBPROGRAM-2

C CALCULATE ELEMENT STIFFNESS MATRIX

SUBROUTINE STE(IE,NJ,NE,LND,X,Y,E,PR,T,KE)

DIMENSION LND(NE,3),X(NJ),Y(NJ),B(3,6),D(3,3) REAL KE(6,6)