有限单元法FORTRAN程序设计作业

**2)/AEРDO 30 I=1,3РDO 30 J=1,3РM=2*(I-1)+1РN=2*(J-1)+1РEK(M,N)=D*(B(I)*B(J)+C(I)*C(J)*(1-UM)/2)РEK(M,N+1)=D*(UM*B(I)*C(J)+C(I)*B(J)*(1-UM)/2)РEK(M+1,N)=D*(UM*C(1)*B(J)+B(I)*C(J)*(1-UM)/2)Р30 EK(M+1,N+1)=D*(C(I)*C(J)+B(I)*B(J)*(1-UM)/2)РRETURNРENDРSUBROUTINE STE (XYE,QE,NE,LS,ME,E,UM,T)РDIMENSION XYE(LS,NE),QE(ME),SG(3),B(3),C(3),S(3,6)РB(1)=XYE(2,2)-XYE(2,3)РB(2)=XYE(2,3)-XYE(2,1)РB(3)=XYE(2,1)-XYE(2,2)РC(1)=XYE(1,3)-XYE(1,2)РC(2)=XYE(1,1)-XYE(1,3)РC(3)=XYE(1,2)-XYE(1,1)РAE=(B(2)*C(3)-B(3)*C(2))/2РD=E/2/(1-UM**2)/AEРDO 220 I=1,3РS(1,2*I-1)=D*B(I)РS(2,2*I-1)=D*UM*B(I)РS(3,2*I-1)=D*(1-UM)*C(I)/2РS(1,2*I)=D*UM*C(I)РS(2,2*I)=D*C(I)Р220 S(3,2*I)=D*(1-UM)*B(I)/2РDO 250 I=1,3РSG(I)=0РDO 250 K=1,MEР250 SG(I)=SG(I)+S(I,K)*QE(K)РWRITE (9,300) SGР300 FORMAT (5X,'SEGMA',3E20.4)РRETURNРENDР输入文本:Р输出文本: