平面有限元程序輸入說明

1、PFEM,EXE平面問題有限元程序的數(shù)據(jù)文件輸入格式：基本輸入數(shù)據(jù)文件名：*.DAT計(jì)算結(jié)果輸出數(shù)據(jù)文件名：*.OUT中間參數(shù)輸出數(shù)據(jù)文件名：*.OU1下面按在文件中的前后順序進(jìn)行說明：1 基本控制參數(shù)信息：NG,NE,MC,NX,NB,EO,VO,DENSITY ,T(共計(jì)5個(gè)整形數(shù)，4個(gè)實(shí)型數(shù)) NG：結(jié)構(gòu)的結(jié)點(diǎn)總數(shù)； NE：結(jié)構(gòu)的單元總數(shù)； MC：平面問題的類型，MC=0，為平面應(yīng)力，MC=1，為平面應(yīng)變； NX：荷載工況數(shù)； NB：支承位移數(shù)； EO：材料彈性模量(Pa)； VO：材料泊松比； DENSITY ：容重(N/m3) T ：材料厚度(m)；2 打印輸出控制參數(shù)：NWA,NE

2、W,NWK,NWP(4個(gè)整形數(shù)) 等于1時(shí)，輸出，否則不輸出。3 單元結(jié)點(diǎn)信息：(K,(IJK(I,K),I=1,3),K=1,NE) (每行4個(gè)整形數(shù)，共計(jì)NE行) K：?jiǎn)卧?hào)； IJK(1,K)：K單元I結(jié)點(diǎn)編號(hào)； IJK(2,K)：K單元J結(jié)點(diǎn)編號(hào)； IJK(3,K)：K單元K結(jié)點(diǎn)編號(hào)；4 結(jié)點(diǎn)坐標(biāo)信息：(K,XY(1,K),XY(2,K),K=1,NG)(每行3個(gè)整形數(shù)，共計(jì)NG行) K：結(jié)點(diǎn)號(hào) XY(1,K)：K結(jié)點(diǎn)X坐標(biāo)； XY(2,K)：K結(jié)點(diǎn)Y坐標(biāo)；5 支承信息：(K,MB(1,K),MB(2,K),ZB(K),K=1,NB)(每行3個(gè)整形數(shù)，1個(gè)實(shí)型數(shù)，共計(jì)NB行) K：支承

3、位移序號(hào)； MB(1,K)：第K個(gè)支承位移所在的結(jié)點(diǎn)號(hào)； MB(2,K)：第K個(gè)支承位移的坐標(biāo)方向；（x向輸入1；y向輸入0）； ZB(K)： 第K個(gè)支承位移的數(shù)值；6 按NX荷載工況數(shù)輸入荷載信息：每一荷載工況如下 (1) NF,NP,NM(3個(gè)整型數(shù)) NF：集中荷載個(gè)數(shù)； NP：分布荷載個(gè)數(shù)； NM：計(jì)自重單元數(shù)； (2) 若NF0，則輸入下面數(shù)據(jù) K,MF(1,K),MF(2,K),ZF(K)(每行3個(gè)整形數(shù)，1個(gè)實(shí)型數(shù)，共計(jì)NF行) K：集中荷載序號(hào)； MF(1,K)：第K個(gè)集中荷載作用的結(jié)點(diǎn)號(hào)； MF(2,K)：第K個(gè)集中荷載的坐標(biāo)方向；（x向輸入1；y向輸入2）； ZF(K)：

4、第K個(gè)集中荷載的數(shù)值； (3) 若NP0，則輸入下面數(shù)據(jù) K,MP(1,K),MP(2,K),ZP(K)(每行3個(gè)整形數(shù)，1個(gè)實(shí)型數(shù)，共計(jì)NP行) K：分布荷載序號(hào)； MP(1,K)：第K個(gè)分布荷載作用的i結(jié)點(diǎn)號(hào)；（x坐標(biāo)較小的結(jié)點(diǎn)號(hào)） MP(2,K)：第K個(gè)分布荷載作用的j結(jié)點(diǎn)號(hào)；（x坐標(biāo)較大的結(jié)點(diǎn)號(hào)）x坐標(biāo)相同時(shí)，y坐標(biāo)較小的在前！ ZP(K)： 第K個(gè)分布荷載的數(shù)值；（斜荷載的正負(fù)與y向分量的正負(fù)同？） ZP(K)： 第K個(gè)分布荷載的數(shù)值；（正荷載的正負(fù)與坐標(biāo)方向同？） (4) 若NM0，則輸入下面數(shù)據(jù) 若NMNE，則表示計(jì)所有單元的自重，不需輸入計(jì)自重的單元號(hào)； 若NM<NE，則

5、需要輸入計(jì)自重的單元號(hào)；例題驗(yàn)證1 (選自徐芝倫編彈性力學(xué)簡(jiǎn)明教程(第二版)設(shè)有對(duì)角受壓的正方形薄板，荷載沿厚度均勻分布，為2N/m，由于xz和yz面均為該薄板的對(duì)稱面，所以只需取1/4部分作為計(jì)算對(duì)象，圖1(a,b)，將該對(duì)象劃分為4個(gè)單元，共有6個(gè)結(jié)點(diǎn)。對(duì)稱面上的結(jié)點(diǎn)沒有垂直于對(duì)稱面的位移分量，因此，在1、2、4三個(gè)結(jié)點(diǎn)上設(shè)置了水平鏈桿支座，在4、5、6三個(gè)結(jié)點(diǎn)設(shè)置了鉛直鏈桿支座，如圖1c所示。SUBROUTINE ABC(NG,NE,IJK,XY,BCA,NWA)IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N)DIMENSION X(2,5),XY(2,NG),

6、IJK(3,NE),B(7),BCA(7,NE) SUBROUTINE SUMK(K,NE,NT,ND,IJK,EK,SK) IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N) DIMENSIONIJ(3),EK(6,6),IJK(3,NE),SK(NT,ND)SUBROUTINE KE(K,NE,T,V,A1,A2,IJK,BCA,EK,NWE)IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N)DIMENSION IJ(3),IJK(3,NE),BCA(7,NE),EK(6,6)SUBROUTINE PF(NF,NT,B)IMPLICIT RE

7、AL*8(A-H,O-Z),INTEGER(I-N)DIMENSION MF(2,NF),ZF(NF),B(NT) SUBROUTINE FNAME(PN,FN2,FN) CHARACTER PN*40,FN2*4,FN*12! 去掉PN中前面的空格SUBROUTINE PP(NP,NG,NT,T,XY,B)IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N)DIMENSION MP(2,NP),ZP(NP),XY(2,NG),B(NT)SUBROUTINE PG(NM,NE,NT,DENSITY,T,IJK,BCA,B)IMPLICIT REAL*8(A-H,O-Z),

8、INTEGER(I-N)DIMENSION IJK(3,NE),BCA(7,NE),MG(NM),B(NT)SUBROUTINE DBC(NB,ND,NT,NX,NX1,MB,ZB,SK,B)IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N)DIMENSION MB(2,NB),ZB(NB),SK(NT,ND),B(NT)SUBROUTINE CALND(NG,NE,IJK,ND)IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N)DIMENSION IJK(3,NE)SUBROUTINE TRES(IP,NG,NE,NX,NT,B,STR,DEL

9、D,TOD,DELST,TOST)IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N)DIMENSION B(NT),STR(3,NE),DELD(2,NG,NX),TOD(2,NG),DELST(3,NE,NX),TOST(3,NE)SUBROUTINE SUPT(IP,NG,NE,NB,NX,T,V,A1,A2,IJK,MB,BCA,DELD,DELSUP,TOTSUP) IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N)DIMENSION IJK(3,NE),MB(2,NB),BCA(7,NE),DELD(2,NG,NX),EK(6,6)DI

10、MENSION D1(6),DELSUP(NB,NX),TOTSUP(NB)SUBROUTINE CLEAR1(N,B) IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N)DIMENSION B(N)SUBROUTINE INPUT(NG,NE,NB,IJK,XY,MB,ZB)IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N)DIMENSION IJK(3,NE),XY(2,NG),MB(2,NB),ZB(NB)SUBROUTINE GAUSS(NT,ND,NX1,NX,SK,B,NWK,NWP)IMPLICIT REAL*8(A-H,O-Z),

11、INTEGER(I-N)DIMENSION SK(NT,ND),A(NT,NT),B(NT)SUBROUTINE CHECK(NT,ND,SK)IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N)DIMENSION SK(NT,ND)SUBROUTINE STRESS(NE,NT,V,A1,A2,IJK,BCA,B,STR)IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N)DIMENSION D1(6),D2(3),B(NT),S(3,6),IJK(3,NE),BCA(7,NE),STR(3,NE)SUBROUTINE OUTPUT(IP,NG,NE,NB,NX,MB,DELD,TOD,DELST,TOST,DELSUP,TOTSUP)IMPLICIT REAL*8(A-H,O-Z),INTEGER(I-N)DIMENSION MB(2,NB),DELD(2,NG,NX),TOD(2,NG)DIMENSION DELST(3,NE,NX),TOST(3,NE),DELSUP(NB,NX),TO