有限元法ppt課件

《有限元法ppt課件》由會員分享，可在線閱讀，更多相關《有限元法ppt課件（109頁珍藏版）》請在裝配圖網上搜索。

工程有限單元法 1 課程介紹 一 課程內容 1 有限元法理論基礎 2 有限元軟件ANSYS應用 二 學習方法 理論與實踐相結合 即通過應用有限元分析實際問題來掌握有限元理論 三 學時數 36學時 理論學時 上機學時 四 考核方式 平時成績 報告成績 工程有限單元法 2 第一章概述 1 1有限元法概述有限元法誕生于20世紀中葉 隨著計算機技術和計算方法的發(fā)展 已成為計算力學和計算工程科學領域里最為有效的方法 它幾乎適用于求解所有連續(xù)介質和場的問題 工程有限單元法 3 一 什么是有限元法 有限元法是將連續(xù)體理想化為有限個單元集合而成 這些單元僅在有限個節(jié)點上相連接 即用有限個單元的集合來代替原來具有無限個自由度的連續(xù)體 工程有限單元法 4 有限元方法是分析連續(xù)體的一種很有效的近似計算方法 是計算機問世以后迅速發(fā)展起來的一種廣泛用于工程結構建模與分析的方法 說明工程實際問題與計算方法息息相關 自然現象的背后都對應有相關的物理本質與事物規(guī)律 用數學方法對物理本質與事物規(guī)律進行描述可以得到普適性定律和特定性定理 以及各種形式的 如代數 微分或積分 數學方程 即數學模型 工程有限單元法 5 對于一個實際的工程問題 建立數學模型時 不僅需要根據實際物理背景采用有效的數學方法 還要考慮求解的效率 結果的精度以及方法的適用性等因素 即分析方法 常用的分析方法有 1 對線性的 邊界規(guī)則的簡單問題 一般可以利用解析法 得到精確解 2 對于許多實際工程問題 由于研究系統的龐大 使得微分方程 邊界和初始條件的復雜性大大增加 一般難以得到它的精確解 對非線性的 邊界不規(guī)則等問題 一般不存在精確的解析解 只能利用數值法 如 有限差分法FDM 有限元方法FEM等 得到近似解 工程有限單元法 6 有限元方法的發(fā)展 首先 有限元方法在航空結構分析中取得了明顯的成效1941年 Hrenikoff利用框架分析法 frameworkmethod 分析平面彈性體 將平面彈性體描述為桿和梁的組合體 1943年 Courant在采用三角形單元及最小勢能原理研究扭轉問題時 利用分片連續(xù)函數在子域中近似描述未知函數此后 有限元方法在固體力學 溫度場和溫升應力 流體力學 流固耦合 水彈性 問題 均有發(fā)展 工程有限單元法 7 現如今 有限元法廣泛應用于航空航天 汽車工業(yè) 橋梁 建筑 電子產品 重型機械 微機電系統 生物醫(yī)學等設計過程中的結構與力學分析 實例1 EMA 火箭發(fā)動機 衛(wèi)星 雷達 工程有限單元法 8 實例2 汽車 工程機械 工程有限單元法 9 工程有限單元法 10 工程有限單元法 11 二 有限元法的基本思想 有限元法的基本思想是 分與合 分 是為了劃分單元 進行單元分析 合 則是為了集合單元 對整體結構進行綜合分析 結構離散 單元分析 整體求解 工程有限單元法 12 2 1有限元法的實現過程 工程有限單元法 13 1 對象離散化當研究對象為連續(xù)介質問題時 首先需要將所研究的對象進行合理的離散化分割 即根據精度預期或經驗將連續(xù)問題進行有限元分割 2 單元分析有限元方法的核心工作是單元分析 通過分析各單元的結點力與結點位移之間的關系和邊界條件 以便建立單元剛度矩陣 3 構造總體方程將單元剛度矩陣組成總體方程剛度矩陣 且總體方程應滿足相鄰單元在公共結點上的位移協調條件 即整個結構的所有結點載荷與結點位移之間應存在相互的變量關系 工程有限單元法 14 4 解總體方程在求解有限元模型時 應考慮總體剛度方程中引入的邊界條件 以便得到符合實際情況的唯一解 5 輸出結果有限元模型求解結束后 可通過數值解序列或由其構成的圖形顯示研究對象的物理結構變形情況以及各種物理量間的變化關系 如通過列表顯示各種數據信息 用等值線分布圖顯示等受力點 或動畫顯示各種量的變化過程 工程有限單元法 15 1 直接方法直接方法是指直接從結構力學引伸得到 直接方法具有簡單 物理意義明確 易于理解等特點 2 變分方法變分方法是一種最常用的方法之一 主要用于線性問題的模型建立 3 加權殘值法對于線性自共軛形式方程 加權殘值法可得到和變分法相同的結果 如得到一個對稱的剛度矩陣 對于那些 能量泛函 不存在的問題 主要是一些非線性問題和依賴于時間的問題 加權殘值法是一種很有效的方法 2 2建立有限元方程的常用方法 工程有限單元法 16 通常 實際工程問題可分為線性問題和非線性問題 邊界規(guī)則與不規(guī)則問題 有限元法其實是非線性問題 如圖右所示 2 3有限元法與工程求解問題的關系 工程有限單元法 17 三 有限元法的基本步驟 無論對于什么樣的結構 有限元分析過程都是類似的 其基本步驟為 1 研究分析結構的特點 包括結構形狀與邊界 載荷工況等 2 將連續(xù)體劃分成有限單元 形成計算模型 包括確定單元類型與邊界條件 材料特性等 工程有限單元法 18 3 以單元節(jié)點位移作為未知量 選擇適當的位移函數來表示單元中的位移 再用位移函數求單元中的應變 根據材料的物理關系 把單元中的應力也用位移函數表示出來 最后將作用在單元上的載荷轉化成作用在單元上的等效節(jié)點力 建立單元等效節(jié)點力和節(jié)點位移的關系 這一過程就是單元特性分析 工程有限單元法 19 4 利用結構力的平衡條件和邊界條件把各個單元按原來的結構重新連接起來 集合成整體的有限元方程 求解出節(jié)點位移 重點 對于不同的結構 要采用不同的單元 但各種單元的分析方法又是一致的 工程有限單元法 20 四 有限元法的學習路線 從最簡單的平面結構入手 由淺入深 介紹有限元理論及其相關應用 工程有限單元法 21 五 有限元法的發(fā)展與應用 有限元法不僅能應用于結構分析 還能解決歸結為場問題的工程問題 從二十世紀六十年代中期以來 有限元法得到了巨大的發(fā)展 為工程設計和優(yōu)化提供了有力的工具 工程有限單元法 22 一 算法與有限元軟件 從二十世紀60年代中期以來 進行了大量的理論研究 不但拓展了有限元法的應用領域 還開發(fā)了許多通用或專用的有限元分析軟件 理論研究的一個重要領域是計算方法的研究 主要有 大型線性方程組的解法 非線性問題的解法 工程有限單元法 23 目前應用較多的通用有限元軟件如下表 另外還有許多針對某類問題的專用有限元軟件 例如金屬成形分析軟件Deform Autoform 焊接與熱處理分析軟件SysWeld等 工程有限單元法 24 二 應用實例 有限元法已經成功地應用在以下一些領域 固體力學 包括強度 穩(wěn)定性 震動和瞬態(tài)問題的分析 傳熱學 電磁場 流體力學 工程有限單元法 25 轉向機構支架的強度分析 劉道勇 東風汽車工程研究院動 用MSC Nastran完成 工程有限單元法 26 基于ANSYS的齒輪嚙合仿真 工程有限單元法 27 第2章彈性力學基本方程及平面問題的有限元法 工程有限單元法 28 2 1彈性力學簡介 本課程中的有限單元法理論要用到彈性力學的某些基本概念和基本方程 將簡單介紹這些概念和方程 作為彈性力學有限單元法的預備知識 工程有限單元法 29 彈性力學 區(qū)別與聯系 材料力學 1 研究的內容 基本上沒有什么區(qū)別 彈性力學也是研究彈性體在外力作用下的平衡和運動 以及由此產生的應力和變形 2 研究的對象 有相同也有區(qū)別 材料力學基本上只研究桿 梁 柱 軸等桿狀構件 即長度遠大于寬度和厚度的構件 彈性力學雖然也研究桿狀構件 但還研究材料力學無法研究的板與殼及其它實體結構 即兩個尺寸遠大于第三個尺寸 或三個尺寸相當的構件 工程有限單元法 30 彈性力學 區(qū)別與聯系 材料力學 3 研究的方法 有較大的區(qū)別 雖然都從靜力學 幾何學與物理學三方面進行研究 但是在建立這三方面條件時 采用了不同的分析方法 材料力學是對構件的整個截面來建立這些條件的 因而要常常引用一些截面的變形狀況或應力情況的假設 這樣雖然大大簡化了數學推演 但是得出的結果往往是近似的 而不是精確的 而彈性力學是對構件的無限小單元體來建立這些條件的 因而無須引用那些假設 分析的方法比較嚴密 得出的結論也比較精確 所以 我們可以用彈性力學的解答來估計材料力學解答的精確程度 并確定它們的適用范圍 工程有限單元法 31 彈性力學 區(qū)別與聯系 材料力學 例如 材料力學在研究有孔的拉伸構件通常就假定拉應力在凈截斷面均勻分布 工程有限單元法 32 彈性力學 區(qū)別與聯系 材料力學 總之 彈性力學與材料力學既有聯系又有區(qū)別 它們都同屬于固體力學領域 但彈性力學比材料力學 研究的對象更普遍 分析的方法更嚴密 研究的結果更精確 因而應用的范圍更廣泛 但是 彈性力學也有其固有的弱點 由于研究對象的變形狀態(tài)較復雜 處理的方法又較嚴謹 因而解算問題時 往往需要冗長的數學運算 但為了簡化計算 便于數學處理 它仍然保留了材料力學中關于材料性質的假定 工程有限單元法 33 彈性力學基本方程 一 彈性力學中的幾個基本概念 1 體力 是分布于物體體積內的外力 如重力 磁力 慣性力等 單位體積內的體力亦可分解為三個成分 用記號X Y Z表示 2 面力 是分布于物體表面的力 如靜水壓力 一物體與另一物體之間的接觸壓力等 單位面積上的表面力通常分解為平行于座標軸的三個成分 用記號來表示 工程有限單元法 34 3 內力 平均應力和應力 1 內力 Internalforces 是物體本身不同部分之間相互作用的力 2 平均應力 theaveragestress 設作用在包含P點某一個截面mn上的單元面積 elementaryarea A上的力為 F 則 F A稱為 A上的平均應力 3 應力 如果假設內力分布連續(xù) 命 A無限減小并趨向P點 則 F A將趨向一個極限p 這個極限P就叫做物體在截面mn上 在P點的應力 彈性體受外力以后 其內部將產生應力 工程有限單元法 35 內力 平均應力和應力的概念 工程有限單元法 36 4 正應力和切應力的概念正應力 應力在作用截面法線方向的分量 切應力 應力在作用截面切線方向的分量 正平行六面體應力 從物體中取出一個微小的正平行六面體 它的棱邊分別平行于三個坐標軸 長度分別為dx dy dz 正平行六面體應力如圖所示 工程有限單元法 37 1 應力的表示正應力用 表示 它的下標表示作用方向 如 x表示正應力沿著x方向 剪應力用 表示 它有兩個下標 例如 xy表示剪應力作用在垂直x軸的平面上 但沿著y方向 2 應力的符號如果一個截面的外法線沿著坐標軸的正方向 這個面就稱為正面 這個面上的應力就以沿著坐標軸的正方向為正 沿著坐標軸的負方向為負 工程有限單元法 38 這個應力符號的規(guī)定與材料力學的不同 在材料力學中 正應力的符號為拉為正 壓為負 而剪應力為正面向下的為正 負面向上為正 或用右手法則確定 右手姆指沿面的外法線時 其余四個手指反時針為正 順時針為負 材料力學中正的剪應力 彈性力學中正的剪應力 工程有限單元法 39 剪應力互等定律作用在兩個互相垂直的面上并且垂直于該兩面交線的剪應力是互等的 大小相等 正負號也相同 因此剪應力記號的兩個角碼可以對調 工程有限單元法 40 可以證明 如果這六個量在P點是已知的 就可以求得經過該點的任何面上的正應力和剪應力 因此 這六個量可以完全確定該點的應力狀態(tài) 它們就稱為在該點的應力分量 一般說來 彈性體內各點的應力狀態(tài)都不相同 因此 描述彈性體內應力狀態(tài)的上述六個應力分量并不是常量 而是坐標x y z的函數 六個應力分量的總體 可以用一個列矩陣來表示 工程有限單元法 41 5 形變和正應變 剪應變的概念 1 形變 形狀的改變 它包含長度和角度的改變 2 正應變 各線段單位長度的伸縮 以伸長為正 縮短為負 3 剪應變 各線段之間的直角的改變 6 位移是指位置的移動 它在x y和z軸上的投影用u v和w 來表示 它的符號是沿坐標軸正向為正 沿坐標軸負向為負 工程有限單元法 42 二 彈性力學中關于材料性質的基本假定 1 連續(xù)性 假定物體是連續(xù) 即整個物體的體積都被組成這個物體的介質所填滿 不留任何空隙 這樣 物體內的物理量 例如應力形變和應變 才可能是連續(xù)的 才可以用連續(xù)函數來表示 2 完全彈性 假定物體是完全彈性的 所謂彈性 是指物體在引起形變的外力被除去以后能恢復原形的性質 而完全彈性是指物體能完全恢復原形而沒有任何剩余變形 3 均勻性 假定物體是均勻的 整個物體由同一材料組成 4 各向同性 假定物體是各向同性的 即物體的彈性性質在所有各個方向都相同 符合以上四個假定的物體 稱為理想彈性體 工程有限單元法 43 5 小變形假定 假定物體的位移和形變是微小的 即物體的位移遠小于物體原來的尺寸 而且應變和轉角都遠小于1 因此 本課程所討論的問題 都是理想彈性體的小變形問題 工程有限單元法 44 三 彈性力學的研究方法 位移邊界條件 邊界條件 應力邊界條件 工程有限單元法 45 彈性力學的基本變量 工程有限單元法 46 彈性力學的基本方程 平衡方程 由物體的受力平衡條件建立的方程 工程有限單元法 47 彈性力學的基本方程 幾何方程 由物體的受力變形后 各應變分量和位移分量的關系建立的方程 工程有限單元法 48 彈性力學的基本方程 物理方程 由物體材料本身的物理特性建立的方程 其中E 彈性模量 泊松比 G 剪切彈性模量 且對各向同性材料 工程有限單元法 49 在限元法中 物理方程可表示為 工程有限單元法 50 彈性力學的基本方程 邊界條件 工程有限單元法 51 四 彈性力學問題的解法 空間彈性力學問題共有15個方程 3個平衡方程 6個幾何方程 6個物理方程 其中包括6個應力分量 6個應變分量 3個位移分量 共有15個未知函數 在給定邊界條件時 問題是可解的 彈性力學問題的提法是 給定作用在物理全部邊界或內部的作用 求解物理由此產生的應力場和位移場 工程有限單元法 52 按照三種不同的邊界條件 彈性力學問題可分為應力邊界條件問題 位移邊界問題和混合邊界 由于有限元模型是對實際結構的反映 對有限元模型施加合適的載荷條件和邊界條件 是正確求解有限元解的關鍵 工程有限單元法 53 根據先求出的基本未知量的不同 彈性力學問題有三種方法 1 應力法 以應力分量作為基本未知量 此時將一切未知量和基本方程都轉換為用應力表示 求得應力分量后 由物理方程求應變分量 再由幾何方程求出位移分量 2 位移法 以位移分量作為基本未知量 此時將一切未知量和基本方程都轉換為用位移表示 求得位移分量后 用幾何方程求應變分量 再由物理方程求應力分量 目前 有限元法中多采用位移法的思想 3 混合法 采用各點的一部分位移分量和一部分應力分量作為基本未知量 混合求解 工程有限單元法 54 五 虛功原理及虛功方程 圖1 8a示一平衡的杠桿 對C點寫力矩平衡方程 圖1 8b表示杠桿繞支點C轉動時的剛體位移圖 綜合可得 即 上式是以功的形式表述的 表明 圖a的平衡力系在圖b的位移上作功時 功的總和必須等于零 這就叫做虛功原理 55 虛功原理 進一步分析 當杠桿處于平衡狀態(tài)時 和這兩個位移是不存在的 但是如果某種原因 例如人為地振一下讓它傾斜 一定滿足上式的關系 將這個客觀存在的關系抽象成一個普遍的原理 去指導分析和計算結構 對于在力的作用下處于平衡狀態(tài)的任何物體 不用考慮它是否真正發(fā)生了位移 而假想它發(fā)生了位移 由于是假想 故稱為虛位移 那么 物體上所有的力在這個虛位移上的總功必定等于零 這就叫做虛位移原理 也稱虛功原理 在圖1 8a中的和所作的功就不是發(fā)生在它本身 狀態(tài)a 的位移上 因為它本身是平衡的 不存在位移 而是在狀態(tài) b 的位移上作的功 可見 這個位移對于狀態(tài) a 來說就是虛位移 亦即是狀態(tài) a 假象的位移 工程有限單元法 56 虛功原理 必須指出 虛功原理的應用范圍是有條件的 它所涉及到的兩個方面 力和位移并不是隨意的 對于力來講 它必須是在位移過程中處于平衡的力系 對于位移來講 雖然是虛位移 但并不是可以任意發(fā)生的 它必須是和約束條件相符合的微小的剛體位移 還要注意 當位移是在某個約束條件下發(fā)生時 則在該約束力方向的位移應為零 因而該約束力所作的虛功也應為零 這時該約束力叫做被動力 如圖1 8中的反力 由于支點C沒有位移 故所作的虛功對于零 反之 如圖1 8中的和是在位移過程中作功的力 稱為主動力 因此 在平衡力系中應當分清楚哪些是主動力 哪些是被動力 而在寫虛功方程時 只有主動力作虛功 而被動力是不作虛功的 工程有限單元法 57 虛功原理與虛功方程 虛功原理表述如下 在力的作用下處于平衡狀態(tài)的體系 當發(fā)生與約束條件相符合的任意微小的剛體位移時 體系上所有的主動力在位移上所作的總功 各力所作的功的代數和 恒對于零 虛功原理用公式表示為 這就是虛功方程 其中P和相應的代表力和虛位移 工程有限單元法 58 虛功原理 用于彈性體的情況 虛功方程是按剛體的情況得出的 即假設圖1 8的杠桿是絕對剛性 沒有任何的變形 因而在方程中沒有內功項出現 而只有外功項 將虛功原理用于彈性變形時 總功W要包括外力功 T 和內力功 U 兩部分 即 W T U 內力功 U 前面有一負號 是由于彈性體在變形過程中 內力是克服變形而產生的 所有內力的方向總是與變形的方向相反 所以內力功取負值 根據虛功原理 總功等于零得 T U 0外力虛功T 內力虛功U彈性力學中的虛功原理可表達為 在外力作用下處于平衡狀態(tài)的彈性體 如果發(fā)生了虛位移 那么所有的外力在虛位移上的虛功 外力功 等于整個彈性體內應力在虛應變上的虛功 內力功 工程有限單元法 59 六 兩種平面問題 彈性力學可分為空間問題和平面問題 嚴格地說 任何一個彈性體都是空間物體 一般的外力都是空間力系 因而任何實際問題都是空間問題 都必須考慮所有的位移分量 應變分量和應力分量 但是 如果所考慮的彈性體具有特殊的形狀 并且承受的是特殊外力 就有可能把空間問題簡化為近似的平面問題 只考慮部分的位移分量 應變分量和應力分量即可 平面應力問題平面應變問題 工程有限單元法 60 平面應力問題 厚度為t的很薄的均勻木板 只在邊緣上受到平行于板面且不沿厚度變化的面力 同時 體力也平行于板面且不沿厚度變化 以薄板的中面為xy面 以垂直于中面的任一直線為Z軸 由于薄板兩表面上沒有垂直和平行于板面的外力 所以板面上各點均有 另外由于平板很薄 外力又不沿厚度變化 可認為在整個薄板內各點均有 于是 在六個應力分量中 只需要研究剩下的平行于XOY平面的三個應力分量 即 所以稱為平面應力問題 工程有限單元法 61 平面應力問題 應力矩陣 1 2 可以簡化為 工程有限單元法 62 物理方程 1 10 中后兩式可見 這時的剪應變 由物理方程 1 10 中的第三式可見 一般 并不一定等于零 但可由及求得 在分析問題時不必考慮 于是只需要考慮三個應變分量即可 于是應變矩陣 1 3 2 簡化為 工程有限單元法 63 平面應力問題 物理方程 1 10 簡化為 轉化成應力分量用應變分量表示的形式 工程有限單元法 64 平面應力問題 將 1 21 式用矩陣方程表示 它仍然可以簡寫為 彈性矩陣 D 則簡化為 工程有限單元法 65 平面應力問題 只有三個應變分量需要考慮 所以幾何方程 1 3 簡化為 工程有限單元法 66 平面應力問題 彈性體的虛功方程 1 17 簡化為 工程有限單元法 67 平面應變問題 一縱向 即Z向 很長 且沿橫截面不變的物體 受有平行于橫截面而且不沿長度變化的面力和體力 如圖1 11所示 由于物體的縱向很長 在力學上可近似地作為無限長考慮 截面尺寸與外力又不沿長度變化 當以任一橫截面為xy面 任一縱線為Z軸時 則所有一切應力分量 應變分量和位移分量都不沿Z方向變化 它們都只是x和y的函數 此外 在這一情況下 由于對稱 任一橫截面都可以看作對稱面 所有各點都只會有x和y方向的位移而不會有Z方向的位移 即w 0因此 這種問題稱為平面位移問題 但習慣上常稱為平面應變問題 工程有限單元法 68 平面應變問題 既然w 0 而且u及v又只是x和y的函數 由幾何方程 1 3 1 可見 于是只剩下三個應變分量 幾何方程仍然簡化為方程 1 24 工程有限單元法 69 平面應變問題 因為由物理方程 1 11 中后兩式可見又由物理方程 1 11 中的第三式可見 在平面應變問題中 雖然 但一般并不等于零 不過它可以由及求得 在分析問題時不必考慮 于是也就只有三個應力分量需要考慮 工程有限單元法 70 平面應變問題 物理方程 1 11 簡化為 工程有限單元法 71 平面應變問題 將 1 25 式用矩陣方程表示 它仍然可以簡寫為 彈性矩陣 D 則為 工程有限單元法 72 平面應變問題 平面應變問題 由于在Z方向沒有外力 應力和應變也不沿Z方向變化 所以虛功方程 1 25 仍然適用 其中的t可以取為任意數值 但必須是這個t范圍內的外力 需要說明一下 工程中有許多問題很接近于平面應變問題 如受內壓力的圓管 滾柱軸承中的滾柱等等 但它們的沿Z向長度都不是無限長的 故在靠近兩端的部分 其應力應變狀態(tài)比較復雜 并不符合平面應變問題的條件 因此將這類問題當作平面應變問題來考慮時 對于離開兩端有一定距離的地方 得出的結果還是相當滿意的 但對靠近兩端的部位 卻有較大的出入 往往需要加以處理 工程有限單元法 73 平面應力問題與平面應變問題 對于兩種平面問題 幾何方程都是 1 24 虛功方程都是 1 25 物理方程都是 工程有限單元法 74 平面應力問題與平面應變問題 對于平面應力情況下的彈性矩陣 應該采用 1 23 式 而對于平面應變則采用 1 28 式 還可注意 在 1 23 式中 若將E改換為 將改換為 就得出公式 1 28 工程有限單元法 75 平面應力問題與平面應變問題 在兩種平面問題中 如果 則和1 3中 1 4 式相似 由幾何方程的積分得出 其中及分別代表彈性體沿x及y方向的剛體移動 而代表彈性體繞Z軸的剛體轉動 工程有限單元法 76 2 2平面問題的有限元法 工程有限單元法 77 有限單元法的基本思路 1 把物體分成有限大小的單元 單元間用節(jié)點相連接 2 把單元節(jié)點的位移作為基本未知量 在單元內的位移 設成線性函數 或其它函數 保證在單元內和單元間位移連接 3 將節(jié)點的位移與節(jié)點的力聯系起來 4 列出節(jié)點的平衡方程 得出以節(jié)點位移表達的平衡方程組 5 求解代數方程組 得出各節(jié)點的位移 根據節(jié)點位移求出各單元中的應力 有限單元法的基本未知量是節(jié)點位移 用節(jié)點的平衡方程來求解 工程有限單元法 78 彈性力學平面問題的有限單元法包括三個主要步驟 1 離散化2 單元分析3 單元綜合 1 離散化有限單元法的基礎是用所謂有限個單元的集合體來代替原來的連續(xù)體 因而必須將連續(xù)體簡化為由有限個單元組成的離散體 對于平面問題 最簡單 因而最常用的單元是三角形單元 這些單元在節(jié)點處用鉸相連 荷載也移置到節(jié)點上 成為節(jié)點荷載 在節(jié)點位移或其某一分量可以不計之處 就在節(jié)點上安置一個鉸支座或相應的連桿支座 工程有限單元法 79 2 單元分析對三角形單元 建立節(jié)點位移與節(jié)點力之間的轉換關系 節(jié)點位移 節(jié)點力 80 2 單元分析 單元剛度矩陣取節(jié)點位移作基本未知量 由節(jié)點位移求節(jié)點力 其中 轉換矩陣稱為單元剛度矩陣 單元分析的主要目的就是要求出單元剛度矩陣 單元分析的步驟可表示如下 工程有限單元法 81 3 單元綜合將離散化了的各個單元合成整體結構 利用節(jié)點平衡方程求出節(jié)點位移 在位移法中 主要的任務是求出基本未知量 節(jié)點位移 為此需要建立節(jié)點的平衡方程 工程有限單元法 82 i點總的節(jié)點力應為 根據節(jié)點的平衡條件 得單元e的節(jié)點力 可按式 2 2 用節(jié)點位移表示 代入得到用節(jié)點位移表示的平衡方程 每個可動節(jié)點有兩個未知位移 有兩個平衡方程 所以方程總數與未知位移總數相等 可以求出所有的節(jié)點位移 單元綜合的目的就是要求出節(jié)點位移 節(jié)點位移求出后 可進一步求出各單元的應力 工程有限單元法 83 2 2 1平面問題的離散化 對任何工程平面構件進行有限元分析 首先都是從簡化其幾何形狀 繪出其平面簡圖入手 連續(xù)體的離散化就是單元網格劃分 平面問題中最常用的單元是三角形和矩形單元 總之 通過單元劃分 載荷移置以及約束簡化 就形成了有限元模型 工程有限單元法 84 在劃分單元時 應注意以下幾點 1 單元類型的選擇 主要取決于結構的幾何形狀 施加的載荷類型和要求的計算精度 2 單元的大小 即網格的疏密 從有限元的理論上講 單元劃分越細 節(jié)點布置越多 計算結果精度越高 但相應要求計算機容量也增大 計算時間也增加 3 單元有疏有密 對結構的不同部位可采用不同大小的單元 4 不同厚度或不同材料處 應取作為單元的邊界線 而且在該處附近的單元還應劃分的小一些 以盡可能反映出邊線兩側應力的突變情況 5 預留載荷位置 在分布載荷集度變化處和集中力作用處 應布置節(jié)點 以利加載 并且其附近的單元也應劃分的小些 以反映此處的應力變化 工程有限單元法 85 2 2 2單元位移函數 如果彈性體的位移分量是坐標的已知函數 則可用幾何方程求應變分量 再從物理方程求應力分量 但對一個連續(xù)體 內部各點的位移變化情況很難用一個簡單函數來描繪 有限單元法的基本原理是分塊近似 即將彈性體劃分成若干細小網格 在每一個單元范圍內 內部各點的位移變化情況可近似地用簡單函數來描繪 對每個單元 可以假定一個簡單函數 用它近似表示該單元的位移 這個函數稱為位移函數 或稱為位移模式 位移模型 位移場 對于平面問題 單元位移函數可以用多項式表示 多項式中包含的項數越多 就越接近實際的位移分布 越精確 但選取多少項數 要受單元型式的限制 工程有限單元法 86 三節(jié)點三角形單元 六個節(jié)點位移只能確定六個多項式的系數 所以平面問題的3結點三角形單元的位移函數如下 所選用的這個位移函數 將單元內部任一點的位移定為座標的線性函數 位移模式很簡單 位移函數寫成矩陣形式為 工程有限單元法 87 最終確定六個待定系數 工程有限單元法 88 令 下標i j m輪換 簡寫為 I 是單位矩陣 N 稱為形態(tài)矩陣 Ni稱為位移的形態(tài)函數 工程有限單元法 89 選擇單元位移函數時 應當保證有限元法解答的收斂性 即當網格逐漸加密時 有限元法的解答應當收斂于問題的正確解答 因此 選用的位移模式應當滿足下列兩方面的條件 1 必須能反映單元的剛體位移和常量應變 6個參數到反映了三個剛體位移和三個常量應變 2 必須保證相鄰單元在公共邊界處的位移連續(xù)性 線性函數的特性 工程有限單元法 90 例題 圖示等腰三角形單元 求其形態(tài)矩陣 N 工程有限單元法 91 由三角形的面積 工程有限單元法 92 本節(jié)利用幾何方程 物理方程 實現用結點位移表示單元的應變和單元的應力 用結點位移表示單元的應變的表達式為 B 矩陣稱為幾何矩陣 2 2 3單元應變和應力 工程有限單元法 93 對于平面應力問題 工程有限單元法 94 2 2 4單元剛度矩陣 單元節(jié)點力與單元位移的關系式 稱為單元剛度方程組 工程有限單元法 95 單元剛度矩陣的性質 1 單元剛度矩陣中每個元素有明確的物理意義 2 剛度矩陣是對稱矩陣 3 剛度矩陣是奇異矩陣 另外 單元剛度矩陣取決于 1 單元的位移函數 2 單元的幾何參數 3 單元的材料性質 工程有限單元法 96 2 2 5單元等效節(jié)點載荷 連續(xù)彈性體離散為單元組合體時 為簡化受力情況 需把彈性體承受的任意分布的載荷都向節(jié)點移置 分解 而成為結點載荷 如果彈性體受承受的載荷全都是集中力 則將所有集中力的作用點取為節(jié)點 就不存在移置的問題 集中力就是節(jié)點載荷 但實際問題往往受有分布的面力和體力 都不可能只作用在節(jié)點上 因此 必須進行載荷移置 如果集中力的作用點未被取為節(jié)點 該集中力也要向結節(jié)移置 將載荷移置到節(jié)點上 必須遵循靜力等效的原則 靜力等效是指原載荷與節(jié)點載荷在任意虛位移上做的虛功相等 在一定的位移模式下 移置結果是唯一的 且總能符合靜力等效原則 工程有限單元法 97 在線性位移模式下 對于常見的一些載荷 可以通過簡單的虛功計算 得出所需的載荷列矩陣 均質等厚度的三角形單元所受的重力 把1 3的重力移到每個節(jié)點 工程有限單元法 98 例 總載荷的2 3移置到節(jié)點i 1 3移置到節(jié)點j 與原載荷同向 工程有限單元法 99 載荷向節(jié)點的移置 可以用普遍公式來表示 體力的移置分布面力的移置在線性位移模式下 用直接計算法簡單 非線性模式下 要用普遍公式計算 工程有限單元法 100 2 2 6總剛度矩陣 K 為總剛度矩陣 R 為節(jié)點力分量矩陣 為節(jié)點位移分量矩陣 總剛度矩陣性質 1 總剛度矩陣也是對稱矩陣 2 總剛度矩陣呈稀疏帶狀分布 3 總剛度矩陣奇異矩陣 工程有限單元法 101 2 2 7邊界約束條件 有限元法中通常采用兩種方法 劃行劃行法和乘大數法 其中前者適用于簡單的手算練習 后者適合于實際問題的計算機處理 工程有限單元法 102 2 2 8解題步驟與算例 有限元法的一般分析步驟如下 1 首先繪出結構的幾何簡圖 在此基礎上將結構離散 2 其次進行單元分析 3 組集總剛度矩陣 4 最終求單元應力和節(jié)點應力 工程有限單元法 103 算例講解 P27 工程有限單元法 104 2 2 9計算結果處理 有限元中計算結果主要包括位移和應力兩方面 其中位移可根據計算結果中的節(jié)點位移分量畫出結構的位移圖 而對于應力計算結果必須進行整理 方法有 1 形心法 2 繞節(jié)點法 3 二單元法 工程有限單元法 105 2 2 10平面高階單元 為了提高有限元法計算結果的精度 除了增加單元數目外 還常采用具有較高次位移函數的單元 即高階單元 常用的四節(jié)點矩形單元和六節(jié)點三角形單元 工程有限單元法 106 1 四節(jié)點任意四邊形等參數單元 任意四結點四邊形單元 四結點正方形單元 工程有限單元法 107 1 八節(jié)點任意四邊形等參數單元 四邊形八結點單元 八結點基本單元 工程有限單元法 108 3 應用等參單元應注意以下幾點問題 1 各向長度的相對大小 單元長度之比不宜相差太大 接近正方形的單元誤差最小 長寬比很大 誤差也很大 2 棱邊的曲折 應使單元邊上沒有折點 如邊上不可避免有折點 應使棱邊只有凸出的折點 3 棱邊的夾角 盡量接近90度 4 棱邊上節(jié)點的間距 盡量均勻 工程有限單元法 109