文档介绍：C《有限单元法基本原理和数值方法》一书的源程序C*********************************************************************C*CPL----*C*C*********************************************************************CC--------输入数据顺序--------=,MC,NX,NB,ND,EO,VO,E结构单元总数CMC计算控制类型参数CMC=0平面应力C=D结构刚度矩阵的半带宽CEO弹性模量CVO泊松比CT单元(结构),NWK,NWP,NWD5整型C输出控制参数C=1输出C=(3,NE)单元结点编码数组3×NE整型CIJM(1,I);IJM(2,I);IJM(3,I)C第I个三角形单元的结点编号,(2,NG)结构结点坐标数组2×(2,N,ZB(N2×NB整型,NB实型CMB(1,I)---第I个给定位移所在的结点号CMB(2,I)=1--给定X方向位移C=0--给定Y方向位移CZB(I)----给定位移值(以坐标正向为正),F-----作用在结点上的集中载荷(坐标方向)P-----作用均布侧压的单元边数C若NF>(2,NF),ZF(NF)2×NF整型,NF实型CMF(1,I)---第I个集中载荷所在的结点号CMF(2,I)=1--给定X方向集中力C=0--给定Y方向集中力CZF(I)-----作用的集中力值C若NP>(2,NP),ZP(NP)2×NP整型,NP实型CMP(1,I)----第I个载荷作用边的起始结点号CMP(2,I)----第I个载荷作用边的起始结点号CZP(I)------若NX>1重复7.-9.(NX-1)次CC最后NG=--------输出数据顺序--------(3,NE)单元结点编码数组3×NE整型CIJM(1,I);IJM(2,I);IJM(3,I)C第I个三角形单元的结点编号,(2,NG)结构结点坐标数组2×若NWA=1,,B(7)单元参数NE行,1×NE整型,7×NE实型C每行结构为:'NE='+单元号+Bi+Bj+Bm+Ci+Cj+Cm+ACC若NWE=1,,EK(6×6)单元刚度阵NE行,1×NE整型,6×6×NE实型C每行结构为:'NE='+单元号+EK(若NWK=1,(NT,ND)结构刚度矩阵NT×ND=2NG×若NWD=1,,B结点位移数据NG行,1×NG整型,2×NG实型C每行结构为:单元号+U+,S2,S3,X1,X2,CTAC单元应力数据6×NE实型C分别代表σx,σy,τxy,σ1,若NX>1重复6.-7.(NX-1)次CC--------可调数组分配--------CC实型数组C(100000)整型数组IA((1)XY(2,NG)IA(1)IJM(3,(N1)ZB(NIA(M1)MB(2,MCC(N2)BCA(7,NE)IA(M2)MF(2,NCC(N3)SK(NT,ND)IA(M3)MP(2,(N4)F(NT)IA(MEND)(N5)ZF((N6)ZP((NEND)--------程序停止代码--------C0正常停止C111数组C越界C222数组C/IA越界C333单元面积非正C444结构刚度矩阵主元非正C----------------------------------------------------------------------CC主程序CDIMENSIONC(500000),IA(50000),EK(36)CHARACTER*12IN,OUTCIN和OUT为输入文件和输出文件的文件名WRITE(*,*)''WRITE(*,*)'PLEASEINPUTTHEINPUT-FILENAME(A<12)'WRITE(*,*)''READ(*,5)INC输入输入文件的文件名WRITE(*,*)''WRITE(*,*)'PLEASEINPUTTHEOUTPUT-FILENAME(A<12)'WRITE(*,*)''READ(*,5)OUTC输入输出文件的文件名5FORMAT(A12)OPEN(5,FILE=IN,STATUS='OLD')OPEN(6,FILE=OUT,STATUS='UNKNOWN')C打开对应的输入和输出文件10READ(5,*)NGIF()STOPC输入