第3講-簡(jiǎn)單桿系w03ban

有限單元法及軟件應(yīng)用

山東建筑大學(xué)2011進(jìn)度安排1有限元方法概述2數(shù)理力學(xué)基礎(chǔ)3簡(jiǎn)單桿系結(jié)構(gòu)有限元法4彈性力學(xué)平面有限元方法5等參元和高斯積分6空間問(wèn)題有限元法7梁結(jié)構(gòu)單元8板殼問(wèn)題有限元法9結(jié)構(gòu)動(dòng)力問(wèn)題有限元法進(jìn)度安排10材料非線性問(wèn)題11幾何非線性問(wèn)題12熱傳導(dǎo)問(wèn)題13有限元Fortran程序設(shè)計(jì)14ANSYS有限元軟件期末考試桿系結(jié)構(gòu)的有限元法3.1概述3.2彈簧系統(tǒng)的剛度矩陣3.3拉壓直桿的有限元分析3.4梁的有限元分析3.5剛架的有限元分析3.1概述3.1.1桿系結(jié)構(gòu)定義由有限根桿件在它們的端點(diǎn)處相互連接而成的結(jié)構(gòu)分類平面桿系：各桿軸線和外力作用線在一個(gè)平面內(nèi)空間桿系：各桿軸線和外力作用線不在一個(gè)平面內(nèi)工程中常見(jiàn)類型拉壓直桿，桁架（平面和空間），梁（簡(jiǎn)支懸臂梁等），剛架（平面和空間）當(dāng)結(jié)構(gòu)長(zhǎng)度尺寸比兩個(gè)截面方向的尺寸大得很多時(shí)，這類結(jié)構(gòu)稱為桿件。由桿件組成的結(jié)構(gòu)體系稱為桿系。由桁桿組成的桿系稱為桁架；由梁組成的桿系稱為剛架。若桿系和作用力均位于同一平面內(nèi)，則稱為平面桁架或平面剛架，否則稱為空間桁架或空間剛架。3.1概述3.1.2桿系單元定義桿系結(jié)構(gòu)中的桿件、梁、柱等稱為桿系單元。連接的點(diǎn)稱為節(jié)點(diǎn)。桿系單元為一維單元。結(jié)構(gòu)離散一般原則：桿系的交叉點(diǎn)、邊界點(diǎn)、集中力作用點(diǎn)、桿件截面尺寸突變處等都應(yīng)該設(shè)置節(jié)點(diǎn)，節(jié)點(diǎn)之間的桿件即構(gòu)成單元。F節(jié)點(diǎn)1節(jié)點(diǎn)2單元①節(jié)點(diǎn)3節(jié)點(diǎn)2單元②桿系結(jié)構(gòu)：梁、拱、框架、桁架等，它們?？呻x散成桿元和梁元。○○○○○○○○○梁拱框架○○○○○桁架3.1概述3.1.2桿系單元分類桁架單元：桁架中的桿件剛架單元：剛架中的桿件區(qū)別：桁架節(jié)點(diǎn)：鉸節(jié)點(diǎn)傳遞力！剛架節(jié)點(diǎn)：剛節(jié)點(diǎn)

傳遞力和力矩！3.1概述3.1.3桿系單元的有限元分析與平面問(wèn)題和空間問(wèn)題比較，基本流程完全相同；具體計(jì)算細(xì)節(jié)需要按照桿系單元的特性來(lái)進(jìn)行。3.2彈簧系統(tǒng)的剛度矩陣一、單個(gè)彈簧的剛度矩陣彈簧的作用力向量位移向量為為了求出剛度矩陣，將彈簧系統(tǒng)看成兩個(gè)簡(jiǎn)單的系統(tǒng)，(1)只有節(jié)點(diǎn)1變形，節(jié)點(diǎn)2固定。(2)只有接點(diǎn)2變形，接點(diǎn)1固定，如圖3-3(b)所示。(3)根據(jù)線彈簧系統(tǒng)的疊加原理，疊加上述兩種情形作用在節(jié)點(diǎn)1上的合力作用在節(jié)點(diǎn)2上的合力從上式可以看出，這一剛度矩陣是對(duì)稱的另外，這一矩陣是奇異的，即它的行列式的值等于零。二、組合彈簧的剛度矩陣(1)先令只允許節(jié)點(diǎn)1有位移考慮彈簧1-2，由靜力平衡條件有節(jié)點(diǎn)3沒(méi)有力作用，因此只允許節(jié)點(diǎn)2有位移(2)令即彈簧1-2的伸長(zhǎng)量與彈簧2-3的縮短量相等，對(duì)彈簧1-2有拉力對(duì)彈簧2-3有壓力分別對(duì)兩彈簧求靜力平衡有(3)最后令只允許節(jié)點(diǎn)3有位移類似情況(1)(4)合成。對(duì)整個(gè)系統(tǒng)來(lái)說(shuō)有3個(gè)節(jié)點(diǎn)，每個(gè)節(jié)點(diǎn)只有一個(gè)方向的位移。因此方程式為在知道了單個(gè)彈簧單元的剛度矩陣后，是否可以利用它來(lái)直接疊加出多個(gè)彈簧串聯(lián)系統(tǒng)得總剛度矩陣呢？答案是肯定的。仍以上面兩彈簧系統(tǒng)為例，說(shuō)明如何疊加剛度矩陣。疊加剛度矩陣1、首先，寫(xiě)出每個(gè)彈簧單元的受力方程和單元?jiǎng)偠染仃噯卧?：?jiǎn)卧?：2、由于整個(gè)系統(tǒng)有3個(gè)節(jié)點(diǎn)(位移)，將上述方程擴(kuò)大成3階方程，有3、按矩陣相加原則，將兩式相加，有

在上面的疊加單元?jiǎng)偠染仃嚨倪^(guò)程中，我們用到了矩陣擴(kuò)大的辦法。這對(duì)于只有幾個(gè)單元的簡(jiǎn)單系統(tǒng)比較方便，但是隨著系統(tǒng)中單元數(shù)目的增多，相應(yīng)擴(kuò)大后的矩陣就相當(dāng)大，擴(kuò)大后的非零元素在矩陣中的位置就不很清楚。

為了清楚起見(jiàn)，一般采用按節(jié)點(diǎn)號(hào)將相應(yīng)單元的剛度矩陣中元素寫(xiě)到總剛度矩陣中去的辦法疊加。

對(duì)于具有n個(gè)節(jié)點(diǎn)的彈簧系統(tǒng)，由于每個(gè)節(jié)點(diǎn)只有一個(gè)可能的位移方向，(或只有一個(gè)自由度)，因此整個(gè)系統(tǒng)有n

個(gè)自由度，相應(yīng)的總剛度矩陣應(yīng)該是我們先寫(xiě)成一個(gè)空的矩陣,將單元?jiǎng)偠染仃嚢磫卧墓?jié)點(diǎn)號(hào)寫(xiě)到空矩陣中去。第1個(gè)單元的節(jié)點(diǎn)號(hào)為1和2，則單元?jiǎng)偠染仃囋賹⒌诙€(gè)單元(節(jié)點(diǎn)號(hào)為2和3)的剛度矩陣疊加到總剛度矩陣的第2、3行的第2、3列上去，即以上面兩彈簧系統(tǒng)為例，總剛度矩陣應(yīng)該是方程求解如前所述，剛度矩陣是一個(gè)奇異陣，即，它的行列式的值為零，矩陣的逆不存在。對(duì)應(yīng)線性代數(shù)方程組無(wú)定解。物理解釋：對(duì)整個(gè)系統(tǒng)的位移都沒(méi)有加以限制，這樣在任何外力的作用下，系統(tǒng)都會(huì)像剛體一樣產(chǎn)生整體運(yùn)動(dòng)

為使方程組有定解，只需要給系統(tǒng)加大上一定的約束(約束條件、邊界條件)。例如，假設(shè)節(jié)點(diǎn)1固定不動(dòng)從而得到定解。解上述方程可得到各個(gè)節(jié)點(diǎn)的位移，利用所求的位移就可計(jì)算出每個(gè)彈簧所受力的大小。從上面的推導(dǎo)計(jì)算中，可以總結(jié)出用有限元求解彈簧系統(tǒng)受力問(wèn)題的基本步驟:1)形成每個(gè)單元的剛度矩陣2)由每個(gè)單元的剛度矩陣按節(jié)點(diǎn)號(hào)疊加成整個(gè)系統(tǒng)的剛度矩陣3)引入約束條件；4)以節(jié)點(diǎn)位移為未知量求解線性代數(shù)方程組5)用每個(gè)單元的力—位移關(guān)系求單元力。例題3-1如圖3-11所示為三彈簧系統(tǒng)，，節(jié)點(diǎn)1和4固定，在節(jié)點(diǎn)2和節(jié)點(diǎn)3處施加軸向力，求節(jié)點(diǎn)2、節(jié)點(diǎn)3的位移和節(jié)點(diǎn)1、4處的作用力。解：1)求出單個(gè)彈簧的剛度矩陣對(duì)于彈簧1-2對(duì)于彈簧2-3對(duì)于彈簧3-4將這些單個(gè)彈簧的剛度矩陣擴(kuò)展成2)列出邊界條件和結(jié)構(gòu)矩陣方程：為了驗(yàn)證結(jié)果的正確性，我們進(jìn)行受力平衡驗(yàn)證：3.3拉壓直桿的有限元分析3.3.1拉壓直桿（單元描述）幾何形狀：等截面A,長(zhǎng)度為l載荷：沿軸線方向分布節(jié)點(diǎn)：2個(gè)局部坐標(biāo)系：沿軸線定義的一維坐標(biāo)系ox節(jié)點(diǎn)坐標(biāo)

在x軸的坐標(biāo)：xi,xj節(jié)點(diǎn)位移（自由度）

沿x軸的位移：ui

,uj單元節(jié)點(diǎn)位移列陣ijlxuiuj單元分析單元分析的目的是為整體分析做準(zhǔn)備，單元分析就是建立單元桿端力和桿端位移之間的關(guān)系，即單元?jiǎng)偠确匠獭卧治龅囊话悴襟E如下：1)用廣義坐標(biāo)法或試湊法建立形函數(shù)，從而建立滿足變形協(xié)調(diào)的單元位移場(chǎng)，即單元內(nèi)任意一點(diǎn)的位移用單元節(jié)點(diǎn)位移(對(duì)桿件單元為桿端位移)來(lái)表示；2)由幾何方程建立單元應(yīng)變場(chǎng)，即單元內(nèi)任意一點(diǎn)的應(yīng)變由節(jié)點(diǎn)位移來(lái)表示；4)用變形體虛位移原理或最小勢(shì)能原理建立單元?jiǎng)偠确匠?，獲得單元?jiǎng)偠染仃嚭蛦卧刃Ч?jié)點(diǎn)載荷矩陣。3)由物理方程建立單元應(yīng)力場(chǎng)，也即單元內(nèi)任意一點(diǎn)的應(yīng)力由節(jié)點(diǎn)位移來(lái)表示；

拉壓桿單元，已知等直桿件桿長(zhǎng)為橫截面面積為材料彈性模量為所受軸向分布載荷集度為桿端位移分別為桿端力分別記為1、建立位移場(chǎng)設(shè)局部坐標(biāo)系下桿中任意點(diǎn)a的坐標(biāo)為因?yàn)橹挥袃蓚€(gè)邊界條件因此桿軸任意一點(diǎn)的位移可假設(shè)為由結(jié)點(diǎn)位移條件

x=0

時(shí)，u=u1

x=l

時(shí)，u=u2從而可得到a截面位移為那么a截面位移形函數(shù)形函數(shù)矩陣桿端位移向量或者節(jié)點(diǎn)位移向量形函數(shù)具有如下性質(zhì)：1)本端為1，它端為02)單元內(nèi)任意一點(diǎn)總和為12、應(yīng)變分析B為應(yīng)變矩陣或者幾何矩陣。3、應(yīng)力分析單元的應(yīng)變能4、建立單元?jiǎng)偠确匠谈鶕?jù)單元上的載荷，可得單元的外力勢(shì)能為由此可得到單元總勢(shì)能為表示單元等效節(jié)點(diǎn)載荷表示局部坐標(biāo)單元?jiǎng)偠染仃囉煽倓?shì)能一階變分等于零對(duì)于均質(zhì)等截面較支桿來(lái)說(shuō)，它的剛度值也可直接由材料力學(xué)中力與變形的關(guān)系獲得。當(dāng)單元上受有均布滿跨軸向載荷p時(shí)節(jié)點(diǎn)載荷等效當(dāng)單元距離1節(jié)點(diǎn)處有一集中軸向載荷圖示所示桁架試求1-2桿和1-4桿單元的局部坐標(biāo)單元?jiǎng)偠染仃?-2桿：抗拉剛度1-4桿：抗拉剛度3.3拉壓直桿的有限元分析3.3.2位移模式單元位移模式的推導(dǎo)位移模式形函數(shù)ijlxuiuja1a23.3拉壓直桿的有限元分析3.3.3應(yīng)變應(yīng)變分量

拉壓直桿只有軸向應(yīng)變：幾何方程的推導(dǎo)3.3拉壓直桿的有限元分析3.3.4應(yīng)力應(yīng)力分量

拉壓直桿只有軸向應(yīng)力：物理方程的推導(dǎo)3.3拉壓直桿的有限元分析3.3.5單元?jiǎng)偠染仃?.3拉壓直桿的有限元分析3.3.6單元節(jié)點(diǎn)等效載荷（軸向載荷）集中力根據(jù)離散的要求，集中力直接施加在所處節(jié)點(diǎn)上體力軸向分布載荷q(x)推導(dǎo)依據(jù)：

面力按照集中載荷施加在面所在的節(jié)點(diǎn)上3.3拉壓直桿的有限元分析3.3.7平面桁架的有限元分析網(wǎng)格離散單元分析整體分析400mm300mmXY20kN25kN①②③④1234取桿件與桿件交點(diǎn)、集中力作用點(diǎn)、桿件與支承的交點(diǎn)為節(jié)點(diǎn)。相鄰兩節(jié)點(diǎn)間的桿件段是單元。節(jié)點(diǎn)編號(hào)時(shí)力求單元兩端點(diǎn)號(hào)差最小。

ij3.3拉壓直桿的有限元分析3.3.7平面桁架的有限元分析網(wǎng)格離散單元分析：在局部坐標(biāo)系下建立單元平衡方程整體分析：在整體坐標(biāo)系下組裝整體平衡方程因此，組裝過(guò)程中需要兩個(gè)坐標(biāo)系之間的轉(zhuǎn)換：整體坐標(biāo)系:OXY局部坐標(biāo)系:OxyxyXYOα有限元中的坐標(biāo)系有整體坐標(biāo)系和局部坐標(biāo)系。對(duì)于一個(gè)結(jié)構(gòu)，整體坐標(biāo)系一般只有一個(gè)；而局部坐標(biāo)系有很多個(gè)，一個(gè)單元就有一個(gè)局部坐標(biāo)。局部坐標(biāo)系每一個(gè)單元的規(guī)定都是相同的，這樣，同類型單元?jiǎng)偠染仃囅嗤?。XY○○○○○Pxy桿系結(jié)構(gòu)單元主要有鉸接桿單元和梁?jiǎn)卧獌煞N類型。它們都只有2個(gè)節(jié)點(diǎn)i、j。約定：?jiǎn)卧鴺?biāo)系的原點(diǎn)置于節(jié)點(diǎn)i；節(jié)點(diǎn)i到j(luò)的桿軸（形心軸）方向?yàn)閱卧鴺?biāo)系中x軸的正向。y軸、z軸都與x軸垂直，并符合右手螺旋法則。對(duì)于梁?jiǎn)卧?，y軸和z軸分別為橫截面上的兩個(gè)慣性主軸。xyzij··3.3.7平面桁架的有限元分析整體坐標(biāo)系OXY：節(jié)點(diǎn)位移為Ui

,Vi(i,j)局部坐標(biāo)系Oxy：節(jié)點(diǎn)位移為ui

,uj

則有：3.3拉壓直桿的有限元分析ijxyXYOαUjVjUiViujui從整體坐標(biāo)到局部坐標(biāo)的坐標(biāo)變換矩陣[T]寫(xiě)成矩陣形式xyiia+abXYXYxybacdcd-ca取——i節(jié)點(diǎn)在單元坐標(biāo)系中的位移向量——i節(jié)點(diǎn)在結(jié)構(gòu)坐標(biāo)系中的位移向量x對(duì)X、Y的方向余弦y對(duì)X、Y的方向余弦同理可得單元j節(jié)點(diǎn)在單元坐標(biāo)系和結(jié)構(gòu)坐標(biāo)系中的位移向量：有組合上述結(jié)果，得平面桿單元的局部坐標(biāo)單元位移和整體坐標(biāo)單元位移之間關(guān)系：

i、j兩節(jié)點(diǎn)間的位移變換關(guān)系互不耦合。上式可寫(xiě)成坐標(biāo)變換矩陣[T]的計(jì)算式：3.3拉壓直桿的有限元分析3.3.7平面桁架的有限元分析推導(dǎo)：注意：局部坐標(biāo)系下的應(yīng)力和應(yīng)變3.3拉壓直桿的有限元分析3.3.7平面桁架的有限元分析因此，單元?jiǎng)偠染仃囋诰植孔鴺?biāo)系和整體坐標(biāo)系下的變換式：3.3拉壓直桿的有限元分析3.3.7平面桁架的有限元分析在整體坐標(biāo)系下的單元?jiǎng)偠染仃嚍椋盒纬烧w剛度矩陣、引入約束條件、求解位移和單元內(nèi)力例題講解空間桿單元（3DLINK8）ijxylz（1）單元位移向量124536（2）形函數(shù)（3）應(yīng)變矩陣

（4）應(yīng)力矩陣(5)等價(jià)節(jié)點(diǎn)力(6)單元坐標(biāo)單元?jiǎng)偠染仃?/p>

對(duì)于等截面空間桿單元，3.4梁的有限元分析3.4.1純彎梁?jiǎn)卧?單元描述)幾何形狀：長(zhǎng)度l，橫截面為A。材料屬性：彈性模量E，橫截面的慣性矩為I。節(jié)點(diǎn)：i,j共2個(gè)局部坐標(biāo)系：oxy3.4梁的有限元分析材料力學(xué)基礎(chǔ)知識(shí)1xy1y3.4梁的有限元分析3.4.1純彎梁?jiǎn)卧?單元描述)局部坐標(biāo)系節(jié)點(diǎn)坐標(biāo)值：xi=0,xj=l節(jié)點(diǎn)位移值：撓度vi和轉(zhuǎn)角θi節(jié)點(diǎn)力：彎距Mi和剪力Qi因此，單元位移列陣：?jiǎn)卧d荷列陣：3.4梁的有限元分析3.4.2

位移模式

代入單元兩個(gè)節(jié)點(diǎn)的坐標(biāo)和位移條件，即可求解四個(gè)待定常數(shù)a1-a4：梁?jiǎn)卧獌?nèi)一點(diǎn)有2個(gè)位移：v、，因?yàn)椋?dv/dx；僅一個(gè)位移是獨(dú)立的，取v

。3.4梁的有限元分析3.4.3應(yīng)變3.4梁的有限元分析3.4.4應(yīng)力3.4.5單元?jiǎng)偠染仃噯卧胶夥匠蹋?.4梁的有限元分析3.4.6等效節(jié)點(diǎn)載荷若存在集中力或者集中力矩，將作用點(diǎn)取為節(jié)點(diǎn)若存在分布載荷，按照虛功等效的原則進(jìn)行計(jì)算適用情況：截面高度小于長(zhǎng)度的1/5的桿系結(jié)構(gòu)。原因：?jiǎn)卧奈灰颇Ｊ?，決定了沒(méi)有考慮剪切撓度。xyijlq(x)將形函數(shù)矩陣[N]代入上式，積分可得分布荷載的等效結(jié)點(diǎn)力。表1給出了幾種特殊情況的等價(jià)節(jié)點(diǎn)力。荷載分布QiMiQjMjql/2ql2/12ql/2-ql2/123ql/20ql2/307ql/20-ql2/20ql/45ql2/96ql/4-5ql2/96ijqqijqij幾種橫向分布荷載等價(jià)節(jié)點(diǎn)力表1例：有一變截面梁，一端固定，另一端鉸支。梁長(zhǎng)為2l，固支端的截面盡寸為b×1.6h，鉸支端的截面尺寸為b×h。梁上作用均布載荷p0。求梁端的約束反力。xy離散化將梁劃分成2個(gè)單元，3個(gè)結(jié)點(diǎn)。每個(gè)單元長(zhǎng)度為,截面取平均截面。單元?jiǎng)偠染仃噄j1223對(duì)號(hào)入座，組合整體剛度矩陣123荷載等效結(jié)點(diǎn)力向量約束反力向量123總荷載向量引入邊界條件將整體平衡方程中對(duì)應(yīng)的1、2、5行和總剛中1、2、5列刪去，得解方程組，得結(jié)點(diǎn)位移值將結(jié)點(diǎn)位移值代入整體平衡方程，可得約束反力3.4梁的有限元分析3.4.7應(yīng)用實(shí)例12312kN/m1m1m3.5剛架的有限元分析3.5.1平面剛架

對(duì)于小變形問(wèn)題，可以認(rèn)為軸向變形和彎曲變形互不影響，因此，位移模式和形函數(shù)可以分別按一維拉壓桿單元和彎剪平面梁?jiǎn)卧慕Y(jié)果相互獨(dú)立的兩種變形形式軸向拉壓面內(nèi)彎曲因此：剛架單元＝桿單元＋梁?jiǎn)卧植孔鴺?biāo)系:oxyz3.5剛架的有限元分析3.5.1平面剛架兩個(gè)坐標(biāo)系：局部坐標(biāo)系整體坐標(biāo)系3.5.2平面剛架單元（單元描述：局部坐標(biāo)系下）節(jié)點(diǎn)位移軸向位移橫向位移繞z軸的轉(zhuǎn)角節(jié)點(diǎn)載荷軸向力剪力彎矩3.5剛架的有限元分析ijxyijxy2356l14F2F3F5F6lF1F4（1）單元節(jié)點(diǎn)位移和單元節(jié)點(diǎn)荷載①單元節(jié)點(diǎn)位移②單元節(jié)點(diǎn)荷載剛架的有限元分析因此，局部坐標(biāo)系下：?jiǎn)卧?jié)點(diǎn)位移列陣：?jiǎn)卧?jié)點(diǎn)載荷列陣：3.5.2平面剛架單元（單元描述：局部坐標(biāo)系下）3.5剛架的有限元分析3.5.3單元?jiǎng)偠染仃?局部坐標(biāo)系下)3形函數(shù)式中形函數(shù)[N]為：其中,單元分析：局部坐標(biāo)系整體分析：整體坐標(biāo)系因此，需要進(jìn)行坐標(biāo)轉(zhuǎn)換。3.5剛架的有限元分析

由于1

、2

、4、5的性質(zhì)和平面鉸接桿相同，因而有相