(新)平面弹性力学有限元源程序(FORTRAN)——精品.docVIP

平面弹性力学有限元源程序（FORTRAN） 说明： 1 基本控制参数信息： NG,NE,MC,NX,NB,EO,VO,DENSITY ,T(共计5个整形数，4个实型数)? ?NG：结构的结点总数；? ?NE：结构的单元总数；? ?MC：平面问题的类型，MC=0，为平面应力，MC=1，为平面应变；? ?NX：荷载工况数；? ?NB：支承位移数；? ?EO：材料弹性模量(Pa)；? ?VO：材料泊松比；??DENSITY ：容重(N/m3)??T ：材料厚度(m)； 2 打印输出控制参数： NWA,NEW,NWK,NWP(4个整形数)??等于1时，输出，否则不输出。 3 单元结点信息： (K,(IJK(I,K),I=1,3),K=1,NE) (每行4个整形数，共计NE行)??K：单元号；??IJK(1,K)：K单元I结点编号；??IJK(2,K)：K单元J结点编号；??IJK(3,K)：K单元K结点编号； 4??结点坐标信息： ((K,XY(1,K),XY(2,K)),K=1,NG)(每行3个整形数，共计NG行)??K：结点号??XY(1,K)：K结点X坐标；??XY(2,K)：K结点Y坐标； 5 支承信息： ((K,MB(1,K),MB(2,K),ZB(K)),K=1,NB)(每行3个整形数，1个实型数，共计NB行)??K：支承位移序号；??MB(1,K)：第K个支承位移所在的结点号；??MB(2,K)：第K个支承位移的坐标方向；??ZB(K)：??第K个支承位移的数值； 6 按NX荷载工况数输入荷载信息： 每一荷载工况如下?: ?(1) NF,NP,NM(3个整型数)? ???NF：集中荷载个数；? ???NP：分布荷载个数；? ???NM：计自重单元数；?? (2) 若NF≠0，则输入下面数据? ?? ?K,MF(1,K),MF(2,K),ZF(K)(每行3个整形数，1个实型数，共计NF行)? ?? ? K：集中荷载序号；? ?? ? MF(1,K)：第K个集中荷载作用的结点号；? ?? ? MF(2,K)：第K个集中荷载的坐标方向；? ?? ? ZF(K)：??第K个集中荷载的数值； (3) 若NP≠0，则输入下面数据? ?? ?K,MP(1,K),MP(2,K),ZP(K)(每行3个整形数，1个实型数，共计NP行)? ?? ? K：分布荷载序号；? ?? ? MP(1,K)：第K个分布荷载作用的结点号；? ?? ? MP(2,K)：第K个分布荷载的坐标方向；? ?? ? ZP(K)：??第K个分布荷载的数值；? ?(4) 若NM≠0，则输入下面数据? ?? ?若NM≥NE，则表示计所有单元的自重，不需输入计自重的单元号；? ?? ?若NMNE，则需要输入计自重的单元号； $DEBUG ? ? ? ? ?PROGRAM PLANE ? ? ? ? ?IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) ? ? ? ? ?ALLOCATABLE::IJK(:,:),XY(:,:),BCA(:,:),SK(:,:),STR(:,:),MB(:,:),ZB(:),B(:) ? ? ? ? ?ALLOCATABLE::DELD(:,:,:),TOD(:,:),DELST(:,:,:),TOST(:,:),DELSUP(:,:),TOTSUP(:) ? ? ? ? ? ? ? ? DIMENSION EK(6,6) ? ?? ???CHARACTER PN*40,FN*12 ? ?? ???WRITE(*,(A)) 本程序为计算平面问题的有限元程序 ? ?? ???WRITE(*,(A)) 特点：(1)采用三结点三角形单元； ? ?? ???WRITE(*,(A)) ? (2)采用等带宽存贮技术； ? ?? ???WRITE(*,(A)) ? ? (3)采用高斯消元法解线性方程组。 ? ?? ???WRITE(*,(/A)) 输入计算问题名(PN): ? ?? ???READ(*,(A)) PN ? ?? ???CALL FNAME(PN,.DAT,FN) ? ?? ???WRITE(*,(2A)) ??输入数据文件名为：,FN ? ?? ???OPEN(5,FILE=FN,STATUS=OLD) ? ?? ???CALL FNAME(PN,.OUT,FN) ? ?? ???WRITE(*,(/2A)) 结果输出数据文件名为： ,FN ? ?? ???OPEN(6,FILE=FN,STATUS=UNKNOWN) ? ?? ???CALL FNAME(PN,.OU1,FN) ? ?? ???WRITE(*,(/2A)) 参数输出数据文件名为： ,FN ? ?? ???OPEN