如圖1所示的單位厚度變截面懸臂梁，左端固定，右端自由，上邊界承受20N/m的均布荷載，右下端承受60N的集中荷載，材料彈性模量為E=30E7pa，泊松比u=0.3，現(xiàn)用4節(jié)點(diǎn)平面四邊形等參元單元計(jì)算程序分析該平面應(yīng)力問(wèn)題。

圖1

步驟1：劃分網(wǎng)格，進(jìn)行節(jié)點(diǎn)和單元編號(hào)，如圖2。

圖2

步驟2：在各個(gè)EXCEL表格中輸入相關(guān)數(shù)據(jù)。

1.在文件名位Material的excel中輸入彈性模量和泊松比，注意先后順序，如圖3所示。

圖3

2.在文件名為Node的excel中輸入節(jié)點(diǎn)坐標(biāo)，輸入順序與網(wǎng)格劃分時(shí)的節(jié)點(diǎn)編號(hào)順序相同，且不必再輸入節(jié)編號(hào)，如圖4所示。

圖4

3.在文件名為Element的excel中輸入各個(gè)單元的節(jié)點(diǎn)編號(hào)，節(jié)點(diǎn)編號(hào)按逆時(shí)針順序輸入，單元輸入順序與網(wǎng)格劃分的單元編號(hào)順序相同，且不必再輸入單元的編號(hào)，如圖5所示。

圖5

4.在文件名為BoundaryCondition的excel中輸入節(jié)點(diǎn)的約束狀態(tài)，0表示約束，1表示自由，B列為X方向，C列表示Y方向，如圖6所示。

圖6

5.在文件名為ConcentratedLoad的excel中輸入施加在節(jié)點(diǎn)上的集中荷載， B列為X方向，C列表示Y方向，如圖7所示。

圖7

6.在文件名為UniformLoad的excel中輸入均布荷載，????A列單元編號(hào)，B、C列為承受均布荷載的邊兩端節(jié)點(diǎn)編號(hào)，D列表示為X方向施加的均布力，E列表示Y方向施加的均布力，如圖8所示。

圖8

步驟3：運(yùn)行計(jì)算程序，將MATLAB工作目錄設(shè)置FEM，并運(yùn)行main便可以得到計(jì)算結(jié)果，如圖9所示。

圖9

比較圖10—圖13，可以看出該程序的計(jì)算結(jié)果與ABAQUS的計(jì)算結(jié)果是一致的，驗(yàn)證了該程序的準(zhǔn)確性。

圖10 MATLAB計(jì)算結(jié)果(部分節(jié)點(diǎn)位移)

圖11 MATLAB計(jì)算結(jié)果(部分節(jié)點(diǎn)應(yīng)力)

圖12 ABAQUS計(jì)算結(jié)果(部分節(jié)點(diǎn)位移與部分節(jié)點(diǎn)應(yīng)力)

總結(jié)

以上是生活随笔為你收集整理的MATLAB怎么表示均布荷载,MATLAB平面四节点等参元分析程序说明的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。