平面三角形3节点有限元程序

=1,2Р 10 JR(J,I)=1Р DO 20 I=1,NRР READ(4,*) IP,IX,IYР JR(1,IP)=IXР JR(2,IP)=IYР 20 CONTINUEР N=0Р DO 30 I=1,NPР DO 30 J=1,2Р IF (JR(J,I)) 30,30,25Р 25 N=N+1Р JR(J,I)=NР 30 CONTINUEР DO 55 J=1,NMР READ (4,*)jj,(AE(I,jj),I=1,4)РC WRITE(*,910)jj,(AE(I,jj),I=1,4)Р?IF(NI.eq.1)then Р AE(1,jj)=AE(1,jj)/(1.0-AE(2,jj)*AE(2,jj)) Р AE(2,jj)=AE(2,jj)/(1.0-AE(2,jj)) Р?endifР55 CONTINUEР910 FORMAT (/20X,'MATERIAL PROPERTIES'/(3X,I5,4(1x,E8.3)))Р RETURNР ENDРC **********************************************Р SUBROUTINE CBAND (MA,JR,MEL)Р DIMENSION MA(*),JR(2,*),MEL(3,*),NN(6)Р COMMON /CA/ NP,NE,NM,NRР COMMON /CC/ N,MX,NHР DO 65 I=1,NР65 MA(I)=0Р DO 90 IE=1,NEР DO 75 K=1,3Р IEK=MEL(K,IE)Р DO 95 M=1,2Р JJ=2*(K-1)+MР NN(JJ)=JR(M,IEK)Р95 CONTINUEР75 CONTINUEР L=NР DO 80 I=1,6Р NNI=NN(I)Р IF(NNI.EQ.0) GO TO 80