有限元分析程序设计

oformbandedglobalstiffnessmatrix SUBROUTINEBAND_K(N_DOF,N_BAND,N_ELE,IE,N_NODE,IJK_ELE,X,Y,PE,&?PR,PT,AK)DIMENSIONIJK_ELE(500,3),X(N_NODE),Y(N_NODE),AKE(6,6),AK(500,100)N_DOF=2*N_NODE?DO40I=1,N_DOF?DO40J=1,N_BAND40?AK(I,J)=0DO50IE=1,N_ELE?CALLFORM_KE(IE,N_NODE,N_ELE,IJK_ELE,X,Y,PE,PR,PT,AKE)?DO50I=1,3?DO50II=1,2?IH=2*(I-1)+II?IDH=2*(IJK_ELE(IE,I)-1)+II?DO50J=1,3?DO50JJ=1,2?IL=2*(J-1)+JJ?IZL=2*(IJK_ELE(IE,J)-1)+JJ?IDL=IZL-IDH+1?IF(IDL.LE.0)THEN?ELSE?AK(IDH,IDL)=AK(IDH,IDL)+AKE(IH,IL) ENDIF50?CONTINUE?RETURN?tocalculatetheareaofelementSUBROUTINECAL_AREA(IE,N_NODE,IJK_ELE,X,Y,AE)?DIMENSIONIJK_ELE(500,3),X(N_NODE),Y(N_NODE)?I=IJK_ELE(IE,1)?J=IJK_ELE(IE,2)?K=IJK_ELE(IE,3)?XIJ=X(J)-X(I)?YIJ=Y(J)-Y(I)?XIK=X(K)-X(I)?YIK=Y(K)-Y(I)?AE=(XIJ*YIK-XIK*YIJ)/2.0?RETURN?tocalculatetheelasticmatrixofelement