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

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

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

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

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