平面三角形3節(jié)點有限元程序

1、平面三角形３結點有限元程序平面三角形３結點有限元程序1、程序名、程序名：FEMT3.F，FEMT3.EXE2、程序功能、程序功能本程序能計算彈性力學的平面應力問題和平面應變問題；能考慮自重和結點集中力兩種荷載的作用，在計算自重時y軸取垂直向上為正；能處理非零已知位移，如支座沉降的作用。主要輸出的內容包括：結點位移、單元應力、主應力、第一主應力與x軸的夾角以及約束結點的支座反力。程序采用VisualFtran5.0編制而成，輸入數據全部采

2、用自由格式。3、程序流程及框圖、程序流程及框圖圖１1程序流程圖輸入數據信息形成K形成R計算，輸出應力解Kδ＝R，輸出δ啟動停機NI——問題類型標識，0為平面應力問題，1為平面應變問題；NL——受荷載作用的結點的數目；NG——考慮自重作用為1，不計自重為0；ND——非零已知位移結點的數目；NC——要計算支座約束反力的結點數目。⑵材料信息，共NM條，每條依次輸入EO，VO，W，tEO——彈性模量（kNm2）；VO——泊松比；W——材料容重（

3、kNm3）；t——單元厚度（m）。這些信息都存放在數組AE（4，NM）中。⑶坐標信息，共NP條，每條依次輸入IP，X，YIP——結點號；X，Y——分別為結點的x坐標和y坐標。坐標信息存放在數組X(2，NP)中。⑷單元信息，共NE條，每條依次輸入JE，L，Io，Jo，MoJE——單元號；L——為該單元的材料類型號。Io，Jo，Mo——分別為該單元i，j，m的整體編碼。單元信息存放在數組MEO(4NE)中。⑸約束信息共NR條，每條依次輸入一

4、個數??????IPx?y?IP——結點號；Ix，Iy——分別為該結點的約束情況，如果方向受約束時填0，如果自由則填1。⑹荷載信息，共NL條，每條依次輸入yIP，Fx，FyIP——結點號；Fx，Fy——分別為該結點的x，y方向的荷載分量（kN）。結點號存放在數組NF（NL）中，結點荷載分量存放在數組FV（2，NL）中。⑺若ND＞0，輸入非零已知位移信息，共ND條，每條依次輸入IP，ux，uyIP——結點號；ux，uy——分別為該結點x，