第二章平面桁架有限元及程序設(shè)計(jì)

1、有限單元法及程序設(shè)計(jì)有限單元法及程序設(shè)計(jì)第一章 緒論第二章 平面桁架有限元分析第三章 平面剛架有限元分析第四章 平面問(wèn)題有限元分析有限單元法及程序設(shè)計(jì)有限單元法及程序設(shè)計(jì)第一章第一章 緒論緒論1.1 有限單元法的概念1.2 有限單元法的基本步驟1.3 常用有限單元分析軟件第一章第一章 緒論緒論數(shù)值方法（模擬） 有限元方法 邊界元方法 有限差分方法 有限體積法 無(wú)網(wǎng)格方法1.1 有限單元法的概念有限單元法的概念基本思路：借助數(shù)學(xué)和力學(xué)知識(shí)基本思路：借助數(shù)學(xué)和力學(xué)知識(shí), ,利用計(jì)算機(jī)技術(shù)解決工利用計(jì)算機(jī)技術(shù)解決工程技術(shù)問(wèn)題程技術(shù)問(wèn)題有限元分析是利用數(shù)學(xué)近似分析方法對(duì)真實(shí)物理系統(tǒng)（幾何、載荷工況）進(jìn)

2、行模擬，利用簡(jiǎn)單而又相互作用的元素，即單元，用有限數(shù)量的未知量去逼近無(wú)限未知量的真實(shí)系統(tǒng)。有限單元法(FEM)是20世紀(jì)50年代以來(lái)隨著計(jì)算機(jī)的廣泛應(yīng)用而發(fā)展起來(lái)的一種現(xiàn)代數(shù)值解法。該方法首先應(yīng)用在連續(xù)力學(xué)領(lǐng)域飛機(jī)結(jié)構(gòu)靜、動(dòng)態(tài)特性分析中。隨后很快就廣泛應(yīng)用于求解傳導(dǎo)、電磁場(chǎng)、流體力學(xué)等連續(xù)性問(wèn)題。Rxyo真實(shí)半圓周長(zhǎng)：真實(shí)半圓周長(zhǎng)：RL近似半圓周長(zhǎng)：近似半圓周長(zhǎng)：NilL數(shù)學(xué)：數(shù)學(xué)：離散化的手段，將偏微分方程、變分方程變成代離散化的手段，將偏微分方程、變分方程變成代數(shù)方程；數(shù)方程；計(jì)算曲線長(zhǎng)度：計(jì)算曲線長(zhǎng)度：將曲線離散成將曲線離散成N段，每段用直線長(zhǎng)度代替段，每段用直線長(zhǎng)度代替曲線長(zhǎng)度，直線長(zhǎng)

3、度總和即可曲線長(zhǎng)度，當(dāng)曲線長(zhǎng)度，直線長(zhǎng)度總和即可曲線長(zhǎng)度，當(dāng)N為無(wú)窮大時(shí)，為無(wú)窮大時(shí)，即為精確解。即為精確解。1.1 有限單元法的概念有限單元法的概念01234567891001234Arc LengthNumber1.1 有限單元法的概念有限單元法的概念圓周率的計(jì)算方法：圓周率的計(jì)算方法：RlRLNi小數(shù)點(diǎn)后面小數(shù)點(diǎn)后面20位：位：3.141592653589793238461.1 有限單元法的概念有限單元法的概念力學(xué)：力學(xué)：離散化的手段，利用單元力學(xué)特性，將連續(xù)體離離散化的手段，利用單元力學(xué)特性，將連續(xù)體離散成有限個(gè)單元的組合結(jié)構(gòu)；散成有限個(gè)單元的組合結(jié)構(gòu)；齒輪的應(yīng)力分布：齒輪的應(yīng)力分布：

4、將齒輪分割成很多單元，由結(jié)點(diǎn)相連將齒輪分割成很多單元，由結(jié)點(diǎn)相連并傳遞求解信息，用有限個(gè)連續(xù)體的組合代替原來(lái)的齒并傳遞求解信息，用有限個(gè)連續(xù)體的組合代替原來(lái)的齒輪結(jié)構(gòu)。輪結(jié)構(gòu)。1.1 有限單元法的概念有限單元法的概念電動(dòng)鏟運(yùn)機(jī)電動(dòng)鏟運(yùn)機(jī)1.1 有限單元法的概念有限單元法的概念液壓挖掘機(jī)動(dòng)臂液壓挖掘機(jī)動(dòng)臂1.1 有限單元法的概念有限單元法的概念在建的大連國(guó)際貿(mào)易中心大廈在建的大連國(guó)際貿(mào)易中心大廈(78層，層，342米米)63863個(gè)梁柱單元；個(gè)梁柱單元；34180個(gè)結(jié)點(diǎn)個(gè)結(jié)點(diǎn)；1.1 有限單元法的概念有限單元法的概念在建的大連市體育中心在建的大連市體育中心1.1 有限單元法的概念有限單元法的概念

5、1.1 有限單元法的概念有限單元法的概念溪洛渡拱壩溪洛渡拱壩 (高高297米米)456 個(gè)壩體單元，個(gè)壩體單元，1040個(gè)個(gè)地基單元，地基單元，2163 個(gè)結(jié)點(diǎn)；個(gè)結(jié)點(diǎn)；1.1 有限單元法的概念有限單元法的概念1.11.1 有限單元法的概念有限單元法的概念導(dǎo)管架碼頭結(jié)構(gòu)1.11.1 有限單元法的概念有限單元法的概念高樁碼頭高樁碼頭下部結(jié)構(gòu)高樁碼頭上部結(jié)構(gòu)1.11.1 有限單元法的概念有限單元法的概念游艇碼頭結(jié)構(gòu)1.11.1 有限單元法的概念有限單元法的概念重力式碼頭結(jié)構(gòu)有限元模型有限元模型 是真實(shí)系統(tǒng)理想化的數(shù)學(xué)抽象是真實(shí)系統(tǒng)理想化的數(shù)學(xué)抽象。真實(shí)系統(tǒng)真實(shí)系統(tǒng)有限元模型有限元模型有限元模型有限

6、元模型1.1 有限單元法的概念有限單元法的概念單元 網(wǎng)格劃分中每一個(gè)小的塊體結(jié)點(diǎn) 確定單元形狀、單元之間相互聯(lián)結(jié)的點(diǎn)結(jié)點(diǎn)力 單元上結(jié)點(diǎn)處的結(jié)構(gòu)內(nèi)力載荷 作用在單元結(jié)點(diǎn)上的外力（集中力、分布力）約束 限制某些結(jié)點(diǎn)的某些自由度 單元單元單元單元載荷載荷結(jié)結(jié)點(diǎn)點(diǎn)結(jié)點(diǎn)力結(jié)點(diǎn)力約束約束有限元單元模型中幾個(gè)重要概念有限元單元模型中幾個(gè)重要概念1.1 有限單元法的概念有限單元法的概念1.2 有限單元法的基本步驟有限單元法的基本步驟一、結(jié)構(gòu)離散化一、結(jié)構(gòu)離散化將求解域離散成單元表示的組合體，單元以結(jié)點(diǎn)相連；二、選擇插值函數(shù)（位移模式）二、選擇插值函數(shù)（位移模式）選擇插值函數(shù)來(lái)表達(dá)單元內(nèi)場(chǎng)變量的變化規(guī)律三、形成

7、單元性質(zhì)的矩陣方程三、形成單元性質(zhì)的矩陣方程場(chǎng)變量：標(biāo)量 向量 力、位移 張量 應(yīng)力、應(yīng)變單元?jiǎng)偠确匠蹋簡(jiǎn)卧Y(jié)點(diǎn)力與單元結(jié)點(diǎn)位移之間的關(guān)系；1.2 有限單元法的基本步驟有限單元法的基本步驟四、形成整體系統(tǒng)的矩陣方程四、形成整體系統(tǒng)的矩陣方程綜合所有單元性質(zhì)的矩陣，形成整體系統(tǒng)的矩陣方程；五、約束處理，求解系統(tǒng)方程五、約束處理，求解系統(tǒng)方程引入邊界條件，求解結(jié)點(diǎn)上的未知場(chǎng)變量；六、其他參數(shù)的計(jì)算六、其他參數(shù)的計(jì)算利用已經(jīng)求解的場(chǎng)變量，計(jì)算其他場(chǎng)變量；1.3 常用的有限元分析軟件常用的有限元分析軟件(1) ANSYS功能強(qiáng)大、模塊多、比較通用；土木工程：CivilFEM是商業(yè)化比較早的一個(gè)軟件（收

8、購(gòu)了一些其他軟件公司）(2) MSC產(chǎn)品系列多、通用軟件、二次開(kāi)發(fā)功能強(qiáng)；土木工程：MSC.MARC系列，Patran, MSC Nastran 及AdamsAdams -多體動(dòng)力學(xué)；多體動(dòng)力學(xué)；Actran-聲學(xué)仿真；聲學(xué)仿真；Easy5-控制仿真工具；控制仿真工具；Marc-非線性；非線性；SimXpert-多學(xué)科仿真；多學(xué)科仿真；MSC Nastran-結(jié)構(gòu)化結(jié)構(gòu)化與多學(xué)科與多學(xué)科FEA；Dytran-顯式動(dòng)力學(xué)與流固耦合；顯式動(dòng)力學(xué)與流固耦合；MSC Fatigue-基于基于FE的耐久性的耐久性仿真工具；仿真工具；Sinda-高級(jí)熱分析解決方案；高級(jí)熱分析解決方案；Digimat-非線

9、性，多尺非線性，多尺度的材料與結(jié)構(gòu)建模平臺(tái)；度的材料與結(jié)構(gòu)建模平臺(tái)；SimDesigner-CAD嵌入式多學(xué)科仿真；嵌入式多學(xué)科仿真；Patran-有限元分析解決方案；有限元分析解決方案；SimManager-仿真數(shù)據(jù)和流程管理仿真數(shù)據(jù)和流程管理1.31.3 常用的有限元分析軟件常用的有限元分析軟件(3) ADINA復(fù)雜非線性問(wèn)題與復(fù)雜動(dòng)力問(wèn)題；具有直接求解和迭代求解兩種流固耦合分析方法Automatic Dynamic Incremental Nonlinear Analysis (4) ABAQUS （專注結(jié)構(gòu)分析）（專注結(jié)構(gòu)分析）主要應(yīng)用于結(jié)構(gòu)動(dòng)力彈塑性分析；二次開(kāi)發(fā)功能；HKS公司產(chǎn)品

10、1.31.3 常用的有限元分析軟件常用的有限元分析軟件(5) SAP2000Systems Applications and Products in Data Processing通用結(jié)構(gòu)分析通用結(jié)構(gòu)分析 第二章第二章 平面桁架有限元分析及程序設(shè)計(jì)平面桁架有限元分析及程序設(shè)計(jì)2.1 平面桁架單元的離散2.2 平面桁架單元分析2.3 結(jié)點(diǎn)平衡與整體剛度矩陣的集成2.4 邊界條件的處理2.5 單元內(nèi)力與支座反力的計(jì)算2.6 平面桁架有限元程序設(shè)計(jì)有限單元法及程序設(shè)計(jì)有限單元法及程序設(shè)計(jì)解題方法方法方法1：節(jié)點(diǎn)法節(jié)點(diǎn)法BaFP1ACDaaaFP方法方法2：截面法截面法靜定桁架回顧回顧第二章第二章 平

11、面桁架有限元分析及程序設(shè)計(jì)平面桁架有限元分析及程序設(shè)計(jì)解題方法：力法和位移法超靜定桁架P1234 1l1l2lPNN21cos2cos142411113414AElN222224AElN第二章第二章 平面桁架有限元分析及程序設(shè)計(jì)平面桁架有限元分析及程序設(shè)計(jì)如圖所示桁架，求各桿軸力。力的平衡條件：位移的協(xié)調(diào)方程：PN1N1N21桿和3桿位移：2桿位移：超靜定桁架41412111411111coscoscosvkvlAElAENNy4242222422222vkvlAElAENNy12111coslAEk2222lAEk 第二章第二章 平面桁架有限元分析及程序設(shè)計(jì)平面桁架有限元分析及程序設(shè)計(jì)1桿軸

12、力豎向分量：2桿軸力：式中： 和 為桿件的剛度系數(shù)；1k2k物理意義：物理意義：4點(diǎn)產(chǎn)生單位位移，桿端產(chǎn)生的豎向桿端力；由桿件的物理性質(zhì)和幾何性質(zhì)決定；V4為第4節(jié)點(diǎn)豎向位移P1234 1l1l2lPNNyy212cos)2(coscos2114111kkPkvkNNyPvkk421)2(2142kkPv2124222kkPkvkN第二章第二章 平面桁架有限元分析及程序設(shè)計(jì)平面桁架有限元分析及程序設(shè)計(jì)超靜定桁架代入平衡方程：結(jié)構(gòu)的整體剛度系數(shù)位移法求解超靜定結(jié)構(gòu)。位移法求解超靜定結(jié)構(gòu)。離散原則：每個(gè)結(jié)點(diǎn)離散后還是一個(gè)結(jié)點(diǎn)，每個(gè)桿件離散后變成一個(gè)單元1結(jié)構(gòu)的離散化：盡量將結(jié)構(gòu)離散成數(shù)量最少的等截

13、面直結(jié)構(gòu)的離散化：盡量將結(jié)構(gòu)離散成數(shù)量最少的等截面直桿單元桿單元234562.1 平面桁架單元的離散平面桁架單元的離散9個(gè)單元，6個(gè)結(jié)點(diǎn)1234567815121413111616個(gè)單元，8個(gè)結(jié)點(diǎn)2.2.1 局部坐標(biāo)系下的單元?jiǎng)偠染仃嚲植孔鴺?biāo)系的建立局部坐標(biāo)系的建立iE，A，ljyxe 軸：沿單元的桿軸方向；軸：沿單元的桿軸方向； x2.2 平面桁架的單元分析平面桁架的單元分析 軸：從軸：從 軸軸逆逆時(shí)針旋轉(zhuǎn)時(shí)針旋轉(zhuǎn)90。yx 原點(diǎn)：以第一個(gè)結(jié)點(diǎn)為坐標(biāo)原點(diǎn)；原點(diǎn)：以第一個(gè)結(jié)點(diǎn)為坐標(biāo)原點(diǎn)； 桿端位移：桿端位移：ijeijeijyxeiuju桿端力：桿端力：iUjU符號(hào)：符號(hào)：與坐標(biāo)系的方向一致為

14、正，反之為負(fù)與坐標(biāo)系的方向一致為正，反之為負(fù)。單元右端桿端力：?jiǎn)卧叶藯U端力：?jiǎn)卧蠖藯U端力：?jiǎn)卧蠖藯U端力：?jiǎn)卧獞?yīng)力：?jiǎn)卧獞?yīng)力：?jiǎn)卧獞?yīng)變：?jiǎn)卧獞?yīng)變：右結(jié)點(diǎn)固定右結(jié)點(diǎn)固定結(jié)點(diǎn)位移：結(jié)點(diǎn)位移：luj左結(jié)點(diǎn)固定左結(jié)點(diǎn)固定桿的受力分為兩種情況：桿的受力分為兩種情況：2.2 平面桁架的單元分析平面桁架的單元分析luiiiuu 0ju0iujjuu iulEEluEEjiiulAEAUjiulAEAUijulAEAUjjulAEAU任意情況（左右結(jié)點(diǎn)均有變形）即為以上兩種狀態(tài)的疊加：任意情況（左右結(jié)點(diǎn)均有變形）即為以上兩種狀態(tài)的疊加：jijjiiulAEulAEUulAEulAEU jiejijjjii

15、jiijijiuukuukkkkuulAEUU1111 jijiuuSuulAEN11 11lAES桿端力為：桿端力為： ek式中式中 為單元?jiǎng)偠染仃嚍閱卧獎(jiǎng)偠染仃?局部坐標(biāo)系局部坐標(biāo)系)桿單元軸力為：桿單元軸力為： S式中式中 為單元應(yīng)力為單元應(yīng)力(廣義廣義)矩陣；矩陣；2.2 平面桁架的單元分析平面桁架的單元分析1,1ijijuuuAENAAEAEull單元桿端力方程：?jiǎn)卧獥U端力方程： jijjiiulAEulAEUulAEulAEU jiejijiuukuulAEUU1111 jijiuuSuulAEN1100 jiVVjiVV000000000000 jivvjivvjivvjivv0

16、0桿端位移：桿端位移：ijeijeiuju桿端力：桿端力：iUjU0jivv0jiVV單元軸力：?jiǎn)卧S力： 2.2 平面桁架的單元分析平面桁架的單元分析桿端位移和桿端力桿端位移和桿端力符號(hào)：符號(hào)：與坐標(biāo)系的方向一致為正，反之為負(fù)與坐標(biāo)系的方向一致為正，反之為負(fù)。桿端力：桿端力：iujujUjVjviViviUxyNNiuiviUiVjujUjVjv2.2 平面桁架的單元分析平面桁架的單元分析2.2.2 整體坐標(biāo)系下的單元?jiǎng)偠染仃嚾艟植孔鴺?biāo)系與整體坐標(biāo)系重合，則整體坐標(biāo)系下的單元?jiǎng)偠染仃嚺c局部坐標(biāo)下的單元?jiǎng)偠染仃囅嗤?。若局部坐?biāo)系與整體坐標(biāo)系不重合，如下圖所示：桿端位移：桿端位移：桿端力：桿端力

17、：桿端位移：桿端位移：i 結(jié)點(diǎn)：結(jié)點(diǎn)：j 結(jié)點(diǎn)：結(jié)點(diǎn)：iidxu sincos222)()(ijijyyxxl)()(ijijijijdydyyydxdxxxldl)()()()(ijijijijdydylyydxdxlxxdliidyv jjdxu jjdyv )()(jijivvluulldl)()(jijivvuulAEEAN2.2 平面桁架的單元分析平面桁架的單元分析設(shè)桿件的長(zhǎng)度為 l ，則：兩邊微分：由于桿件的變形產(chǎn)生位移：因此，桿件應(yīng)變?yōu)椋簵U件軸力為：符號(hào)：符號(hào)：桿件軸力以拉為正，壓為負(fù)桿件軸力以拉為正，壓為負(fù)。)()(ijijdydydxdxNNUicos eejjiijjiie

18、kvuvulAEVUVUF22222222NNVisinNNUjcosNNVjsin桿件的結(jié)點(diǎn)力為：因此，桿件結(jié)點(diǎn)力向量為： ek式中 是整體坐標(biāo)系下的單元?jiǎng)偠染仃嚕?.2 平面桁架的單元分析平面桁架的單元分析1,11,1iijjiijjuvAENuvluvAEulv jijjjiijiijiekkkkFFF iiiVUF jjjVUF iiivu jjjvu 22lAEkii 22lAEkjj jjiijjiievuvulAEVUVUF22222222寫成分塊矩陣形式：式中：2.2 平面桁架的單元分析平面桁架的單元分析 lAEScossin 22lAEkkjiij（1 1）單元?jiǎng)偠认禂?shù)）單元

19、剛度系數(shù)kij的意義的意義j自由度自由度(結(jié)點(diǎn)結(jié)點(diǎn))產(chǎn)生的單位桿端位移引起的產(chǎn)生的單位桿端位移引起的i自由度自由度(結(jié)點(diǎn)結(jié)點(diǎn))的桿端力的桿端力（2 2）單元?jiǎng)偠染仃囀菍?duì)稱矩陣）單元?jiǎng)偠染仃囀菍?duì)稱矩陣反力互等定理反力互等定理式中：桿件單元的應(yīng)力矩陣為：桿件單元的應(yīng)力矩陣為：?jiǎn)卧獎(jiǎng)偠染仃嚨男再|(zhì)單元?jiǎng)偠染仃嚨男再|(zhì)2.2 平面桁架的單元分析平面桁架的單元分析（3 3）單元?jiǎng)偠染仃囈话闶遣豢赡娴模﹩卧獎(jiǎng)偠染仃囈话闶遣豢赡娴?.2.3 單元坐標(biāo)轉(zhuǎn)換矩陣單元坐標(biāo)轉(zhuǎn)換矩陣iujujUjVjviViviUxy jjiijjiievuvulAEVUVUF00000101000001012.2 平面桁架的單元分析平

20、面桁架的單元分析取任意桿件取任意桿件, ,建立如圖所示的局部坐標(biāo)系：建立如圖所示的局部坐標(biāo)系：桿端力：桿端力：iuiviUiVjujUjVjv桿端位移：桿端位移：2.2.3 單元坐標(biāo)轉(zhuǎn)換矩陣單元坐標(biāo)轉(zhuǎn)換矩陣iujujUjVjviViviUxyxyiViViUxyxyiUsincosiiiVUUcossiniiiVUVcossinsincosijjjjjVUVVUU2.2 平面桁架的單元分析平面桁架的單元分析桿端力：桿端力：iuiviUiVjujUjVjv桿端位移：桿端位移：在上圖中在上圖中, ,建立如圖所示的整體坐標(biāo)系：建立如圖所示的整體坐標(biāo)系：以以i i結(jié)點(diǎn)為例：結(jié)點(diǎn)為例：cosiUsini

21、VsiniUcosiV同理，對(duì)于同理，對(duì)于j j 結(jié)點(diǎn)：結(jié)點(diǎn)：jjiijjiiVUVUVUVU00000000 eeFTF 00000000Tcossin TTT I 1 TTT2.2 平面桁架的單元分析平面桁架的單元分析寫成矩陣形式：寫成矩陣形式：因此：因此：其中，其中，TT為轉(zhuǎn)換矩陣：為轉(zhuǎn)換矩陣：轉(zhuǎn)換矩陣的性質(zhì)轉(zhuǎn)換矩陣的性質(zhì)轉(zhuǎn)換矩陣是正交矩陣；轉(zhuǎn)換矩陣是正交矩陣； 00000000000001010000010100000000lAEke eeT eeeeeeFTFkkT eeeeTeeekTkTTkTF1 TkTkeTe22222222lAE同理，位移也存在轉(zhuǎn)換關(guān)系：同理，位移也存在轉(zhuǎn)換

22、關(guān)系：代入局部坐標(biāo)系下代入局部坐標(biāo)系下的剛度方程：的剛度方程：2.2 平面桁架的單元分析平面桁架的單元分析與利用微分得到的單元在總體坐標(biāo)下的剛度方程相同xyF1、對(duì)總體結(jié)點(diǎn)位移和單元進(jìn)行編碼；2、單元局部坐標(biāo)系下的剛度矩陣;123例：例：如圖所示平面桁架，桿長(zhǎng)為l，截面積為A，求三個(gè)單元在整體坐標(biāo)系下的剛度矩陣。 kkk 212300232100002123002321T2.2 平面桁架的單元分析平面桁架的單元分析0000010100000101lAE3、單元整體坐標(biāo)系分析:o6021cos23sin解：解：?jiǎn)卧w坐標(biāo)系下的剛度矩陣為: TkTkT 212300232100002123002

23、321T kk TkTkT33333131333331314lAE4、單元整體坐標(biāo)系分析:o0 IT 5、單元整體坐標(biāo)系分析:33333131333331314lAEo12021cos23sinxyF123例：如圖所示平面桁架，桿長(zhǎng)為例：如圖所示平面桁架，桿長(zhǎng)為l，截面積為，截面積為A，求結(jié)構(gòu)的剛度矩陣。，求結(jié)構(gòu)的剛度矩陣。2.3 結(jié)點(diǎn)平衡與整體剛度矩陣的集成結(jié)點(diǎn)平衡與整體剛度矩陣的集成 33333131333331314lAEk 0000010100000101lAEk 33333131333331314lAEkjjjiijiikkkkjjjiijiikkkkjjjiijiikkkk1、單元

24、整體坐標(biāo)系下剛度矩陣分塊解：解：2.3.1 結(jié)點(diǎn)的平衡方程結(jié)點(diǎn)的平衡方程2.3 結(jié)點(diǎn)平衡與整體剛度矩陣的集成結(jié)點(diǎn)平衡與整體剛度矩陣的集成 jjjijijkkF jjjijijkkF 01PFFjjxyF1ij jF iF jF jF iF jFij2 2、結(jié)點(diǎn)、結(jié)點(diǎn)1 1的平衡方程：的平衡方程： 1PFFjj結(jié)點(diǎn)結(jié)點(diǎn)1 1的受力狀態(tài)為（如右圖）的受力狀態(tài)為（如右圖）: :結(jié)點(diǎn)結(jié)點(diǎn)1 1的平衡條件為的平衡條件為: :由單元的剛度方程：由單元的剛度方程：由單元的剛度方程：由單元的剛度方程： 12jjjikk 13jjjikk2.3 結(jié)點(diǎn)平衡與整體剛度矩陣的集成結(jié)點(diǎn)平衡與整體剛度矩陣的集成 1131

25、2Pkkkkjjjijjji 1321Pkkkkjijijjjj代入結(jié)點(diǎn)代入結(jié)點(diǎn)1 1的平衡條件的平衡條件: :3 3、結(jié)點(diǎn)、結(jié)點(diǎn)2 2的平衡方程：的平衡方程： 12ijiijijiiiikkkkF 2PFFii 32ijiijijiiiikkkkF同理，結(jié)點(diǎn)同理，結(jié)點(diǎn)2 2的平衡條件為的平衡條件為: :由單元的剛度方程：由單元的剛度方程：由單元的剛度方程：由單元的剛度方程：2.3 結(jié)點(diǎn)平衡與整體剛度矩陣的集成結(jié)點(diǎn)平衡與整體剛度矩陣的集成代入結(jié)點(diǎn)代入結(jié)點(diǎn)2 2的平衡條件的平衡條件: :4 4、結(jié)點(diǎn)、結(jié)點(diǎn)3 3的平衡方程：的平衡方程：同理，結(jié)點(diǎn)同理，結(jié)點(diǎn)3 3的平衡條件為的平衡條件為: :由單元

26、的剛度方程：由單元的剛度方程：由單元的剛度方程：由單元的剛度方程： 23212Pkkkkijiiijii 2321Pkkkkijiiiiij 3PFFij 32jjjijjjijijkkkkF 13ijiijijiiiikkkkF2.3 結(jié)點(diǎn)平衡與整體剛度矩陣的集成結(jié)點(diǎn)平衡與整體剛度矩陣的集成代入結(jié)點(diǎn)代入結(jié)點(diǎn)3 3的平衡條件的平衡條件: :5 5、系統(tǒng)的平衡方程：、系統(tǒng)的平衡方程： 31332Pkkkkijiijjji 3321Pkkkkiijjjiij 3321 Pkkkkiijjjiij 1321 Pkkkkjijijjjj 2321 Pkkkkijiiiiij2.3 結(jié)點(diǎn)平衡與整體剛度矩

27、陣的集成結(jié)點(diǎn)平衡與整體剛度矩陣的集成321321 PPPkkkkkkkkkkkkiijjjiijijiiiiijjijijjjj PK 6 6、結(jié)構(gòu)整體剛度方程、結(jié)構(gòu)整體剛度方程寫成矩陣形式，即可得到結(jié)構(gòu)的整體剛度方程寫成矩陣形式，即可得到結(jié)構(gòu)的整體剛度方程 333231232221131211112221121211111221212222 K KKKKKKKKKkkkkkkkkkkkk其中，其中，KK為結(jié)構(gòu)的為結(jié)構(gòu)的整體剛度矩陣整體剛度矩陣；2.3.2 整體剛度矩陣的集成步驟1 1、定位、定位單元單元結(jié)點(diǎn)結(jié)點(diǎn)編號(hào)編號(hào) 2 2、累加、累加整體整體結(jié)點(diǎn)結(jié)點(diǎn)編號(hào)編號(hào) ijmn單元?jiǎng)偠认禂?shù)單元?jiǎng)偠?/p>

28、系數(shù) 整體剛度系數(shù)整體剛度系數(shù)2.3 結(jié)點(diǎn)平衡與整體剛度矩陣的集成結(jié)點(diǎn)平衡與整體剛度矩陣的集成eiikmmKeijkmnKejiknmKejjknnKTenm單元定位向量單元定位向量)()(ji)2() 1 (xyF1232.3 結(jié)點(diǎn)平衡與整體剛度矩陣的集成結(jié)點(diǎn)平衡與整體剛度矩陣的集成 jjjiijiikkkkk jjjiijiikkkkk jjjiijiikkkkk例：求上例平面桁架結(jié)構(gòu)的整體剛度矩陣；例：求上例平面桁架結(jié)構(gòu)的整體剛度矩陣；1 1、定位、定位單元：?jiǎn)卧?K1 2 3123(2)(1)(2) (1) 2 2、累加、累加12k11k21k22kT122.3 結(jié)點(diǎn)平衡與整體剛度矩

29、陣的集成結(jié)點(diǎn)平衡與整體剛度矩陣的集成1 1、定位、定位單元單元 ：(2)(1)(2) (1) 2 2、累加、累加1 1、定位、定位單元單元 ：2 2、累加、累加 K1 2 312312k11k12k22k (1) (2) (1)(2)11k12k21k22k K1 2 312312k11k12k22k11k12k21k22k11k12k21k22kT32T132.3 結(jié)點(diǎn)平衡與整體剛度矩陣的集成結(jié)點(diǎn)平衡與整體剛度矩陣的集成 666564636261565554535251464544434241363534333231262524232221161514131211K KKKKKKKKKKKK

30、KKKKKKKKKKKKKKKKKKKKKKKK 333231232221131211112221121211111221212222 K KKKKKKKKKkkkkkkkkkkkk結(jié)構(gòu)的整體剛度矩陣為結(jié)構(gòu)的整體剛度矩陣為: :將每個(gè)字塊展開(kāi)，結(jié)構(gòu)的整體剛度矩陣為將每個(gè)字塊展開(kāi)，結(jié)構(gòu)的整體剛度矩陣為: :下標(biāo)表示自下標(biāo)表示自由 度 編 號(hào)由 度 編 號(hào)下標(biāo)表示下標(biāo)表示結(jié)點(diǎn)編號(hào)結(jié)點(diǎn)編號(hào)2.3 結(jié)點(diǎn)平衡與整體剛度矩陣的集成結(jié)點(diǎn)平衡與整體剛度矩陣的集成1 1、整體剛度矩陣的性質(zhì)、整體剛度矩陣的性質(zhì)2.3.2 整體剛度矩陣的集成方法整體剛度矩陣的集成方法 333231232221131211K KKK

31、KKKKKK（1 1）剛度系數(shù)）剛度系數(shù)Kij的意義的意義分塊矩陣：分塊矩陣：（2 2）單元?jiǎng)偠染仃囀菍?duì)稱矩陣）單元?jiǎng)偠染仃囀菍?duì)稱矩陣反力互等定理反力互等定理j結(jié)點(diǎn)結(jié)點(diǎn)產(chǎn)生的單位桿端位移引起的產(chǎn)生的單位桿端位移引起的i結(jié)點(diǎn)結(jié)點(diǎn)的桿端力；的桿端力；j自由度自由度產(chǎn)生的單位桿端位移引起的產(chǎn)生的單位桿端位移引起的i自由度自由度的桿端力；的桿端力；不分塊矩陣：不分塊矩陣：666564636261565554535251464544434241363534333231262524232221161514131211 KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK2.3 結(jié)點(diǎn)平衡

32、與整體剛度矩陣的集成結(jié)點(diǎn)平衡與整體剛度矩陣的集成思考題：思考題：如何由剛度矩陣元素的意義確定整體剛度矩陣每個(gè)元素的如何由剛度矩陣元素的意義確定整體剛度矩陣每個(gè)元素的組成？組成？22K33 K11 K12k12 K13 K 1111kk2222 kk1122kk21K31K32K21k21k23 K2.3 結(jié)點(diǎn)平衡與整體剛度矩陣的集成結(jié)點(diǎn)平衡與整體剛度矩陣的集成1234練習(xí)題：練習(xí)題：利用整體剛度矩陣的意義確定以下桁架剛度矩陣元素（分利用整體剛度矩陣的意義確定以下桁架剛度矩陣元素（分塊）的組成。塊）的組成。11 K23 K13 K24 K44 K34 K2.3.2 整體剛度矩陣的集成步驟1 1、

33、定位、定位單元自由度編號(hào)單元自由度編號(hào) 2 2、累加、累加整體自由度編號(hào)整體自由度編號(hào) ijmn2.3 結(jié)點(diǎn)平衡與整體剛度矩陣的集成結(jié)點(diǎn)平衡與整體剛度矩陣的集成klrsTesrnm 自由度定位向量自由度定位向量)()()()(lkji)4()3()2() 1 (單元?jiǎng)偠认禂?shù)單元?jiǎng)偠认禂?shù) 整體剛度系數(shù)整體剛度系數(shù)eiikmmKeijkmnKejiknmKejjknnK練習(xí)題：利用整體剛度矩陣的意義確定以下桁架剛練習(xí)題：利用整體剛度矩陣的意義確定以下桁架剛度矩陣元素（自由度）的組成。度矩陣元素（自由度）的組成。11 2 3 4T21 2 5 6T35 6 7 8T43 4 7 8T57 8 1 2

34、T63 4 5 6T2.4.1 結(jié)點(diǎn)邊界條件2.4 邊界條件的處理邊界條件的處理結(jié)點(diǎn)可以自由變形，整體結(jié)點(diǎn)力等于對(duì)應(yīng)的外荷載。結(jié)點(diǎn)可以自由變形，整體結(jié)點(diǎn)力等于對(duì)應(yīng)的外荷載。1 1、自由變形的結(jié)點(diǎn)、自由變形的結(jié)點(diǎn)2 2、約束結(jié)點(diǎn)或給定了結(jié)點(diǎn)位移的數(shù)值、約束結(jié)點(diǎn)或給定了結(jié)點(diǎn)位移的數(shù)值約束結(jié)點(diǎn)：約束結(jié)點(diǎn)：給定結(jié)點(diǎn)位移：給定結(jié)點(diǎn)位移：0iu0iviiubiivc332211332211666564636261565554535251464544434241363534333231262524232221161514131211 VUVUVUvuvuvuKKKKKKKKKKKKKKKKKKKKKKKKK

35、KKKKKKKKKKK2.4.2 邊界條件的處理方法2.4 邊界條件的處理邊界條件的處理1 1）劃行劃列法）劃行劃列法處理方法：處理方法：若第若第i i個(gè)自由度位移為零，則將總剛第個(gè)自由度位移為零，則將總剛第i i行和第行和第i i列劃列劃掉，剛度矩陣相應(yīng)降低一階。掉，剛度矩陣相應(yīng)降低一階。0 00 00 0優(yōu)點(diǎn)：優(yōu)點(diǎn)：簡(jiǎn)單易行，矩陣降階，減小計(jì)算工作量；簡(jiǎn)單易行，矩陣降階，減小計(jì)算工作量；缺點(diǎn)：缺點(diǎn)：矩陣行列、位移、荷載向量需重新編號(hào)，程序?qū)崿F(xiàn)比較復(fù)雜；矩陣行列、位移、荷載向量需重新編號(hào)，程序?qū)崿F(xiàn)比較復(fù)雜；只適用于約束結(jié)點(diǎn)情況；只適用于約束結(jié)點(diǎn)情況；2.4 邊界條件的處理邊界條件的處理332

36、211332211 666564636261565554535251464544434241363534333231262524232221161514131211 VUVUVUvuvuvuKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK2 2）0 0、1 1置換法（填置換法（填0 0置置1 1法）法）處理方法：處理方法：將與約束自由度對(duì)應(yīng)整體剛度矩陣對(duì)角線元素全部置換將與約束自由度對(duì)應(yīng)整體剛度矩陣對(duì)角線元素全部置換成成1 1，相應(yīng)行和列其他元素置換成，相應(yīng)行和列其他元素置換成0 0，將同一行荷載分量置換成，將同一行荷載分量置換成0 0。0 00 00 01 11 11

37、 10 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 0適用條件：適用條件：只適用于約束結(jié)點(diǎn)，不適用給定位移邊界條件；只適用于約束結(jié)點(diǎn)，不適用給定位移邊界條件；2.4 邊界條件的處理邊界條件的處理332211332211666564636261565554535251464544434241363534333231262524232221161514131211 VUVUVUvuvuvuKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK3 3）乘大數(shù)法）乘大數(shù)法處

38、理方法：處理方法：將與約束自由度對(duì)應(yīng)整體剛度矩陣對(duì)角線元素乘以一個(gè)將與約束自由度對(duì)應(yīng)整體剛度矩陣對(duì)角線元素乘以一個(gè)大數(shù)大數(shù)N N (10(101010-10-101515) )，將同一行荷載分量置換成，將同一行荷載分量置換成N N與對(duì)角線元素的乘積與對(duì)角線元素的乘積與給定位移之積。與給定位移之積。NNNNK33b1NK44b2NK66b3優(yōu)點(diǎn)：優(yōu)點(diǎn)：處理工作量小，適用于給定位移情況；處理工作量小，適用于給定位移情況；0002.5 單元軸力及支座反力的計(jì)算單元軸力及支座反力的計(jì)算單元軸力：?jiǎn)卧S力： (1)整體坐標(biāo)系下的單元桿端位移整體坐標(biāo)系下的單元桿端位移(2)整體坐標(biāo)系下的單元軸力整體坐標(biāo)系

39、下的單元軸力 e e eeSN2.5.1 單元軸力的計(jì)算位移和荷載向量分解為自由結(jié)點(diǎn)和約束結(jié)點(diǎn)兩部分，剛度矩陣相應(yīng)位移和荷載向量分解為自由結(jié)點(diǎn)和約束結(jié)點(diǎn)兩部分，剛度矩陣相應(yīng)分塊，如下所示：分塊，如下所示：bababbbaabaaPPKKKK bbbababKKP bbbRFP bbbbababFKKR2.5.2 支座反力的計(jì)算其中：其中： a為自由位移為自由位移; ; b為約束位移為約束位移; ; aP為外荷載為外荷載; ; bP為約束自由度結(jié)點(diǎn)力為約束自由度結(jié)點(diǎn)力; ;因此：因此： bR為支座反力為支座反力; ; bF為約束自由度結(jié)點(diǎn)荷載為約束自由度結(jié)點(diǎn)荷載; ;其中：其中：2.5 單元軸力

40、及支座反力的計(jì)算單元軸力及支座反力的計(jì)算1、計(jì)算分析題、計(jì)算分析題平面桁架結(jié)構(gòu)網(wǎng)格如圖所示，已知平面桁架結(jié)構(gòu)網(wǎng)格如圖所示，已知EA=1500 kN，采用乘大數(shù)方法引，采用乘大數(shù)方法引入支撐條件，試求后處理法引入支撐條件后的滿陣存貯的整體剛度入支撐條件，試求后處理法引入支撐條件后的滿陣存貯的整體剛度矩陣矩陣 K。結(jié)構(gòu)坐標(biāo)系下單元?jiǎng)偠染仃囉?jì)算公式如下：。結(jié)構(gòu)坐標(biāo)系下單元?jiǎng)偠染仃囉?jì)算公式如下： 例題例題 22222222SCSSCSCSCCSCSCSSCSCSCCSCLEAke式中：式中：sin,cosSC3213myx4m2、計(jì)算分析題、計(jì)算分析題按照有限元法的計(jì)算步驟，求圖示桁架結(jié)構(gòu)各桿軸力。已

41、知：按照有限元法的計(jì)算步驟，求圖示桁架結(jié)構(gòu)各桿軸力。已知：EA=10 kn。結(jié)構(gòu)坐標(biāo)系下的單元?jiǎng)偠染仃囉?jì)算公式如下：。結(jié)構(gòu)坐標(biāo)系下的單元?jiǎng)偠染仃囉?jì)算公式如下： ： 例題例題 22222222SCSSCSCSCCSCSCSSCSCSCCSCLEAke式中：式中：sin,cosSC2.6 程序設(shè)計(jì)程序設(shè)計(jì)程序設(shè)計(jì)原則完整性擴(kuò)充性兼容性邏輯性可讀性可維護(hù)性1.模塊化2.6 程序設(shè)計(jì)程序設(shè)計(jì)程序設(shè)計(jì)流程圖（程序框圖）什么是流程圖？“程序流程圖”常簡(jiǎn)稱為“流程圖”，是一種傳統(tǒng)的算法表示法，程序流程圖是人們對(duì)解決問(wèn)題的方法、思路或算法的一種描述。它利用圖形化的符號(hào)框來(lái)代表各種不同性質(zhì)的操作，并用流程線來(lái)連

42、接這些操作。 2、如何畫流程圖編碼和單元測(cè)試 這個(gè)階段的任務(wù)是程序員根據(jù)目標(biāo)系統(tǒng)的性質(zhì)和實(shí)際環(huán)境，選取一種適當(dāng)?shù)母呒?jí)程序設(shè)計(jì)語(yǔ)言(必要時(shí)用匯編語(yǔ)言)，把詳細(xì)設(shè)計(jì)的結(jié)果翻譯成用選定的語(yǔ)言書寫的程序，并且仔細(xì)測(cè)試編寫出的每一個(gè)模塊。 程序員在書寫程序模塊時(shí)，應(yīng)使它的可讀性、可理解性和可維護(hù)性良好。 綜合測(cè)試 這個(gè)階段的任務(wù)是通過(guò)各種類型的測(cè)試，使軟件達(dá)到預(yù)這個(gè)階段的任務(wù)是通過(guò)各種類型的測(cè)試，使軟件達(dá)到預(yù)定的要求。定的要求。 最基本的測(cè)試是集成測(cè)試和驗(yàn)收測(cè)試。集成測(cè)試是根據(jù)最基本的測(cè)試是集成測(cè)試和驗(yàn)收測(cè)試。集成測(cè)試是根據(jù)設(shè)計(jì)的軟件結(jié)構(gòu)，把經(jīng)單元測(cè)試的模塊按某種選定的策略設(shè)計(jì)的軟件結(jié)構(gòu)，把經(jīng)單元測(cè)試的

43、模塊按某種選定的策略裝配起來(lái)，在裝配過(guò)程中對(duì)程序進(jìn)行必要的測(cè)試。驗(yàn)收測(cè)裝配起來(lái)，在裝配過(guò)程中對(duì)程序進(jìn)行必要的測(cè)試。驗(yàn)收測(cè)試是按照需求規(guī)格說(shuō)明書的規(guī)定，由用戶對(duì)目標(biāo)系統(tǒng)進(jìn)行試是按照需求規(guī)格說(shuō)明書的規(guī)定，由用戶對(duì)目標(biāo)系統(tǒng)進(jìn)行驗(yàn)收。驗(yàn)收。 通過(guò)對(duì)軟件測(cè)試結(jié)果的分析可以預(yù)測(cè)軟件的可靠性；反通過(guò)對(duì)軟件測(cè)試結(jié)果的分析可以預(yù)測(cè)軟件的可靠性；反之，根據(jù)對(duì)軟件可靠性的要求也可以決定測(cè)試和調(diào)試過(guò)程之，根據(jù)對(duì)軟件可靠性的要求也可以決定測(cè)試和調(diào)試過(guò)程什么時(shí)候可以結(jié)束。什么時(shí)候可以結(jié)束。 在進(jìn)行測(cè)試的過(guò)程中，應(yīng)該用正式的文檔把測(cè)試計(jì)劃、在進(jìn)行測(cè)試的過(guò)程中，應(yīng)該用正式的文檔把測(cè)試計(jì)劃、詳細(xì)測(cè)試方案以及實(shí)際測(cè)試結(jié)果保存下來(lái)

44、，作為軟件配置詳細(xì)測(cè)試方案以及實(shí)際測(cè)試結(jié)果保存下來(lái)，作為軟件配置的一部分。的一部分。 2.6 程序設(shè)計(jì)程序設(shè)計(jì)2.6.1 程序框圖程序框圖輸入數(shù)據(jù)輸入數(shù)據(jù)單元局部剛度單元局部剛度坐標(biāo)轉(zhuǎn)換矩陣坐標(biāo)轉(zhuǎn)換矩陣單元整體剛度單元整體剛度集成整體剛度矩陣元素集成整體剛度矩陣元素約束條件處理、解方程約束條件處理、解方程計(jì)算單元軸力、約束反力計(jì)算單元軸力、約束反力單元單元循環(huán)循環(huán)包括單元、結(jié)點(diǎn)、材包括單元、結(jié)點(diǎn)、材料、荷載、約束數(shù)據(jù)料、荷載、約束數(shù)據(jù)2.6 程序設(shè)計(jì)程序設(shè)計(jì)2.6.2 程序說(shuō)明程序說(shuō)明1 1、總體剛度矩陣的半帶寬存儲(chǔ)、總體剛度矩陣的半帶寬存儲(chǔ)總體剛度矩陣：總體剛度矩陣： 對(duì)稱稀疏矩陣；對(duì)稱稀疏矩陣； nm1nnnmmnmmKKKK00TenmeiikmmKeijkmnKejiknmKejjknnKwb) 1(max*2nmbw半帶寬：半帶寬：總體剛度矩陣集成：總體剛度矩陣集成：2.6 程序設(shè)計(jì)程序設(shè)計(jì)行號(hào)：行號(hào)：nm1 nnnmmnmmKKKK00wbwbji1jji wbwb01i1jjii1ij1jii