资源描述：
《平面弹性力学有限元源程序(FORTRAN).doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、平面弹性力学有限元源程序（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单元结点
2、信息：(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：
3、支承位移序号； 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个集中荷载作用的结点号； M
4、F(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，则表示计所有单元的自重，不需输入计自重的单元号； 若NM5、DEBUG PROGRAMPLANE IMPLICITREAL*8(A-H,O-Z),INTEGER(I-N) ALLOCATABLE::IJK(:,:),XY(:,:),BCA(:,:),SK(:,:),STR(:,:),MB(:,:),ZB(:),B(:) ALLOCATABLE::DELD(:,:,:),TOD(:,:),DELST(:,:,:),TOST(:,:),DELSUP(:,:),TOTSUP(:) DIMENSIONEK(6,6) CHARACTERPN*40,FN
6、*12 WRITE(*,'(A)')'本程序为计算平面问题的有限元程序' WRITE(*,'(A)')'特点：(1)采用三结点三角形单元；' WRITE(*,'(A)')' (2)采用等带宽存贮技术；' WRITE(*,'(A)')' (3)采用高斯消元法解线性方程组。' WRITE(*,'(/A)')'输入计算问题名(PN):' READ(*,'(A)')PN CALLFNAME(PN,'.DAT',FN) WRITE(*,'(2A)')' 输入数据文件
7、名为：',FN OPEN(5,FILE=FN,STATUS='OLD') CALLFNAME(PN,'.OUT',FN) WRITE(*,'(/2A)')'结果输出数据文件名为：',FN OPEN(6,FILE=FN,STATUS='UNKNOWN') CALLFNAME(PN,'.OU1',FN) WRITE(*,'(/2A)')'参数输出数据文件名为：',FN OPEN(7,FILE=FN,STATUS='UNKNOWN') READ(5,*)NG,N
8、E,MC,NX,NB,EO,VO,DENSITY,T WRITE(6,120)NG,NE,MC,NX,NB WRITE(6,130)EO,VO,DENSITY,T
