假設(shè)一均勻直桿，彈性模量E,，截面積A,，長(zhǎng)度2L兩端固定,，在中間位置施加一向右的載荷P,，求直桿內(nèi)應(yīng)力,。

材料力學(xué)解法：假設(shè)P的作用點(diǎn)發(fā)生了一微小位移X

求解完畢,。

有限元解法思路：

將直桿劃分為2個(gè)單元,，用這連個(gè)離散單元的組合來(lái)逼近連續(xù)體直桿，設(shè)單元編號(hào),，結(jié)點(diǎn)編號(hào)1,，2，3,，單元上軸向位移與軸向內(nèi)力如上圖 u表示結(jié)點(diǎn)位移 N表示結(jié)點(diǎn)內(nèi)力,。

假設(shè)x為單元 內(nèi)任意一點(diǎn)則位移表達(dá)式為（以下的位移表達(dá)式為假設(shè)的，其實(shí)就是有限元中的單元形函數(shù)）

根據(jù)胡克定律,，對(duì)于單元,，由軸向力引發(fā)的位移

同樣對(duì)于單元                    

對(duì)單元，受力分析得出結(jié)點(diǎn)的受力情況

將上述單元的結(jié)點(diǎn)受力表達(dá)式寫成矩陣結(jié)構(gòu)

同樣的方式針對(duì)單元,，那么單元2 的結(jié)點(diǎn)受力表達(dá)式的矩陣結(jié)構(gòu)是

單元的剛度可表示為

在結(jié)點(diǎn)2處的平衡條件      

各結(jié)點(diǎn)的位移連續(xù)條件   

單元的剛度矩陣方程可以組合到如下形式（整體剛度矩陣方程）

可寫成如下形式  P = Kd         K ———— 結(jié)構(gòu)總剛度矩陣    d ———— 結(jié)構(gòu)結(jié)點(diǎn)位移向量   P ———— 結(jié)構(gòu)結(jié)點(diǎn)載荷向量

邊界條件   d1=d3 =0   代入求得

帶入各單元內(nèi)力方程

到現(xiàn)在位置我們并沒(méi)有使用黃色背景部分的公式,，它起什么作用呢？目前我們能得到的是各結(jié)點(diǎn)的位移,，各點(diǎn)的受力,，以及桿的內(nèi)力及內(nèi)應(yīng)力（此時(shí)因?yàn)閭�(gè)單元內(nèi)部的應(yīng)力處處相等），那如何求得單元內(nèi)部各處的位移呢,？那就得使用黃色部分的內(nèi)容了,。比如單元中間位置發(fā)生的位移是多少

，它不僅可以用作單元的內(nèi)插值函數(shù),，即把單元內(nèi)任意一點(diǎn)的位移用結(jié)點(diǎn)位移表示,，而且可作為加權(quán)余量法中的加權(quán)函數(shù),，可以處理外載荷,，將分布力等效為結(jié)點(diǎn)上的集中力和力矩，此外可用于等參單元的坐標(biāo)變換,。

路過(guò)

雷人

握手

2 鮮花

雞蛋

剛表態(tài)過(guò)的朋友 (2 人)