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

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

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

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

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