平面彈性力學有限元源程序(FORTRAN)（doc X頁）

《平面彈性力學有限元源程序(FORTRAN)（doc X頁）》由會員上傳分享，免費在線閱讀，更多相關內容在工程資料-天天文庫。

1、平面彈性力學有限元源程序（FORTRAN）說明：1基本控制參數(shù)信息：NG,NE,MC,NX,NB,EO,VO,DENSITY,T(共計5個整形數(shù)，4個實型數(shù))??NG：結構的結點總數(shù)；??NE：結構的單元總數(shù)；??MC：平面問題的類型，MC=0，為平面應力，MC=1，為平面應變；??NX：荷載工況數(shù)；??NB：支承位移數(shù)；??EO：材料彈性模量(Pa)；??VO：材料泊松比；??DENSITY：容重(N/m3)??T：材料厚度(m)；2打印輸出控制參數(shù)：NWA,NEW,NWK,NWP(4個整形數(shù))??等于1時，輸出，否則不輸出。3單元結點信息：(K,(IJK(I,K),I

2、=1,3),K=1,NE)(每行4個整形數(shù)，共計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個整形數(shù)，共計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個整形數(shù)，1個實型數(shù)，共計NB行)??K：支承位移序號；??MB(1,K)：第K個支承位移所在的結點號；??M

3、B(2,K)：第K個支承位移的坐標方向；??ZB(K)：??第K個支承位移的數(shù)值；6按NX荷載工況數(shù)輸入荷載信息：每一荷載工況如下?:?(1)NF,NP,NM(3個整型數(shù))????NF：集中荷載個數(shù)；????NP：分布荷載個數(shù)；????NM：計自重單元數(shù)；??(2)若NF≠0，則輸入下面數(shù)據(jù)????K,MF(1,K),MF(2,K),ZF(K)(每行3個整形數(shù)，1個實型數(shù)，共計NF行)????K：集中荷載序號；????MF(1,K)：第K個集中荷載作用的結點號；????MF(2,K)：第K個集中荷載的坐標方向；????ZF(K)：??第K個集中荷載的數(shù)值；(3)若NP≠0

4、，則輸入下面數(shù)據(jù)????K,MP(1,K),MP(2,K),ZP(K)(每行3個整形數(shù)，1個實型數(shù)，共計NP行)????K：分布荷載序號；????MP(1,K)：第K個分布荷載作用的結點號；????MP(2,K)：第K個分布荷載的坐標方向；????ZP(K)：??第K個分布荷載的數(shù)值；??(4)若NM≠0，則輸入下面數(shù)據(jù)????若NM≥NE，則表示計所有單元的自重，不需輸入計自重的單元號；????若NM

5、ALLOCATABLE::IJK(:,:),XY(:,:),BCA(:,:),SK(:,:),STR(:,:),MB(:,:),ZB(:),B(:)?????ALLOCATABLE::DELD(:,:,:),TOD(:,:),DELST(:,:,:),TOST(:,:),DELSUP(:,:),TOTSUP(:)????????DIMENSIONEK(6,6)??????CHARACTERPN*40,FN*12??????WRITE(*,'(A)')'本程序為計算平面問題的有限元程序'??????WRITE(*,'(A)')'特點：(1)采用三結點三角形單元；'?????

6、?WRITE(*,'(A)')'?(2)采用等帶寬存貯技術；'??????WRITE(*,'(A)')'??(3)采用高斯消元法解線性方程組。'??????WRITE(*,'(/A)')'輸入計算問題名(PN):'??????READ(*,'(A)')PN??????CALLFNAME(PN,'.DAT',FN)??????WRITE(*,'(2A)')'??輸入數(shù)據(jù)文件名為：',FN??????OPEN(5,FILE=FN,STATUS='OLD')??????CALLFNAME(PN,'.OUT',FN)??????WRITE(*,'(/2A)')'結果輸出數(shù)據(jù)文件名

7、為：',FN??????OPEN(6,FILE=FN,STATUS='UNKNOWN')??????CALLFNAME(PN,'.OU1',FN)??????WRITE(*,'(/2A)')'參數(shù)輸出數(shù)據(jù)文件名為：',FN??????OPEN(7,FILE=FN,STATUS='UNKNOWN')????????READ(5,*)NG,NE,MC,NX,NB,EO,VO,DENSITY,T????????WRITE(6,120)NG,NE,MC,NX,NB????????WRITE(6,130)EO,VO,DENSITY,T?