资源描述:
《平面弹性力学有限元源程序(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