梁單元-有限元分析
2017-04-08 by:CAE仿真在線 來源:互聯(lián)網(wǎng)
有限元法的基本思想是將結構離散化,用有限個容易分析的單元來表示復雜的對象,單元之間通過有限個節(jié)點相互連接,然后根據(jù)變形協(xié)調條件綜合求解。由于單元的數(shù)目是有限的,節(jié)點的數(shù)目也是有限的,所以稱為有限元法(FEM,Finite Element Method)。是隨著電子計算機的發(fā)展而迅速發(fā)展起來的一種彈性力學問題的數(shù)值求解方法。
有限元法是最重要的工程分析技術之一。它廣泛應用于彈塑性力學、斷裂力學、流體力學、熱傳導等領域。有限元法是60年代以來發(fā)展起來的新的數(shù)值計算方法,是計算機時代的產(chǎn)物。雖然有限元的概念早在40年代就有人提出,但由于當時計算機尚未出現(xiàn),它并未受到人們的重視。
隨著計算機技術的發(fā)展,有限元法在各個工程領域中不斷得到深入應用,現(xiàn)已遍及宇航工業(yè)、核工業(yè)、機電、化工、建筑、海洋等工業(yè),是機械產(chǎn)品動、靜、熱特性分析的重要手段。早在70年代初期就有人給出結論:有限元法在產(chǎn)品結構設計中的應用,使機電產(chǎn)品設計產(chǎn)生革命性的變化,理論設計代替了經(jīng)驗類比設計。目前,有限元法仍在不斷發(fā)展,理論上不斷完善,各種有限元分析程序包的功能越來越強大,使用越來越方便。
所謂梁桿結構是指其長度比橫截面尺寸大很多的梁和桿件、以及由它們組成的系統(tǒng),這一類結構的應力、應變和位移都是一個坐標的函數(shù),所以屬于一維單元問題。
特點:桿件位于一個平面內,桿件間用鉸節(jié)點連接,作用力也在該平面內。
單元特性:只承受拉力或壓力。
單元劃分:常采用自然單元劃分。即以兩個鉸接點之間的桿件作為一個單元。為使桁架桿件只產(chǎn)生軸力,桁架的計算常作以下假定:
①桁架中每根桿件的兩端由理想鉸聯(lián)結;
②每根桿件的軸線必須是直線;
③所有桿件的軸線都只交于所聯(lián)理想鉸的幾何中心。
④荷載均只作用于理想鉸的幾何中心。
在此條件下所算得的各種應力稱為主應力。實際上各種桁架結構不可能完全滿足上述各假定,因而桿件將產(chǎn)生彎曲,由這種彎曲而在桿件中所引起的軸向應力稱為次應力。鋼桁架如設計得較合理,次應力一般不太大;鋼筋混凝土或預應力混凝土桁架,由于其結點剛度大、桿件粗短,次應力的影響不能忽視。
特點:與平面桁架的不同在于桿件之間以剛性鉸接點作為連接。
單元特性:除承受拉、壓軸向載荷外,還要承受剪力和彎矩。
單元劃分:除采用自然單元劃分外。即以兩個剛接點之間的桿件作為一個單元。同時也要考慮在集中載荷作用位置、截面突變位置和分布載荷突變的位置增設節(jié)點。
多用鋼筋混凝土或鋼材建造,剛架結構多為空間剛架,但許多空間剛架可分解為平面剛架進行計算。剛架可分為靜定剛架和超靜定剛架,但工程中所應用的主要是超靜定剛架。
特點:簡化了的剛架。
特性:不考慮軸向變形,只考慮垂直于軸線的線位移和軸線平面內的轉角位移。
單元劃分:同剛架。
A.粱在變形前垂直于粱軸線的橫截面,變形后仍然為平面(剛性橫截面假定);
B.橫截面變形后的平面仍與變形后的軸線相互垂直。
也就是說歐拉梁忽略了剪切變形和轉動慣量,認為初始垂直于中性軸的截平面在變形時仍保持為平面垂直于中性軸(Kirchhoff假設),即認為截面的轉動等于撓度曲線切線的斜率。適用于梁的高度遠小于跨度情況下。
A.粱在變形前垂直于粱軸線的橫截面,變形后仍然為平面(剛性橫截面假定);
B.由于歐拉-伯努利梁的第二個假設忽略了梁的剪切變形,對于有效長度較短或復合材料梁板橋時,忽略剪切變形是不妥的,鐵木辛柯提出讓梁的應力應變關系得到滿足,即考慮剪切變形與轉動慣量。
在鐵木辛柯梁中,需要考慮橫向剪切變形影響的情況,如高度相對跨度不太小的高梁。此時梁內的橫向剪切力Q所產(chǎn)生的剪切變形將引起梁的附加撓度,并使原來垂直于中性面的截面變形后不再與中性面垂直,且發(fā)生翹曲。
節(jié)點位移向量
![梁單元-有限元分析ansys仿真分析圖片1](http://m.09bigdata.com/i/g/08/08987e7b3dea87026e0566c557e2ad24r.jpg)
(1-1)
節(jié)點力向量
![梁單元-有限元分析ansys workbanch圖片2](http://m.09bigdata.com/i/g/c1/c1135efa4a7295bb9484a957b386755fr.jpg)
(1-2)
每個單元四個自由度,所以可以有4個待定系數(shù)
![梁單元-有限元分析ansys workbanch圖片3](http://m.09bigdata.com/i/g/df/df2d077e73bec1924b8ae55dabeb85ear.jpg)
(1-3)
有材料力學知,撓度與轉角之間的關系是:
![梁單元-有限元分析ansys workbanch圖片4](http://m.09bigdata.com/i/g/bf/bf0e87a32795be3747f83969a5fc9a7er.jpg)
(1-4)
將i、j節(jié)點的位移代入,求得待定系數(shù):
![梁單元-有限元分析ansys分析圖片5](http://m.09bigdata.com/i/g/1d/1d3f7ca7a743cbda9f210a92255b7675r.jpg)
(1-5)
所以
![梁單元-有限元分析ansys分析圖片6](http://m.09bigdata.com/i/g/0b/0b495a1211b4388c898063ce8ac969d8r.jpg)
(1-6)
其中[N]為形函數(shù)矩陣
![梁單元-有限元分析ansys分析圖片7](http://m.09bigdata.com/i/g/30/30c90a1acd306e94f4fcedf4a5cae680r.jpg)
![梁單元-有限元分析ansys分析圖片8](http://m.09bigdata.com/i/g/44/44b4b8e2fe9a419be4a81178998c6a04r.jpg)
![梁單元-有限元分析ansys結構分析圖片9](http://m.09bigdata.com/i/g/5c/5cd6cea348250dd2b595f92534fb4523r.jpg)
(1-7)
![梁單元-有限元分析ansys結構分析圖片10](http://m.09bigdata.com/i/g/34/347c8bc0e0ecad640c652f5963eee4d6r.jpg)
它也滿足前面形函數(shù)的性質。
梁彎曲變形時,若忽略剪切的影響,由材料力學知道:
![梁單元-有限元分析ansys結構分析圖片11](http://m.09bigdata.com/i/g/7c/7cd89f516ab47ba3bd2c77d5570be975r.jpg)
(1-8)
將它代入到幾何方程,得:
![梁單元-有限元分析ansys結果圖片12](http://m.09bigdata.com/i/g/25/25d21c79d17c90e202276e7f8d60bd6er.jpg)
(1-9)
換入位移插值公式
![梁單元-有限元分析ansys結果圖片13](http://m.09bigdata.com/i/g/fc/fcde6515389af0659039195a16e08943r.jpg)
(1-10)
單向應力狀態(tài),所以
![梁單元-有限元分析ansys結果圖片14](http://m.09bigdata.com/i/g/32/321579c452ca903ab0dc9a8cf6b33f69r.jpg)
(1-11)
![梁單元-有限元分析ansys workbanch圖片15](http://m.09bigdata.com/i/g/0f/0fdee73362d7bb6720968172e2d3c575r.jpg)
(1-12)
代入[B]矩陣和[D]=E,并考慮
![梁單元-有限元分析ansys workbanch圖片16](http://m.09bigdata.com/i/g/14/1400a6203dfc1d913464d5931880e79fr.jpg)
![梁單元-有限元分析ansys workbanch圖片17](http://m.09bigdata.com/i/g/d6/d6c239c9369498dc07bf63cbedeaecacr.jpg)
(1-13)
其中:
![梁單元-有限元分析ansys培訓課程圖片18](http://m.09bigdata.com/i/g/c2/c2498f1e9ccc722e273c2244cb215ef6r.jpg)
——剪切影響系數(shù),As為有效的抗剪切面積。
有一外伸梁及其承載情況如圖所示,其中,
![梁單元-有限元分析ansys培訓課程圖片19](http://m.09bigdata.com/i/g/98/981bb3b486d40c6851fdb65655dc6ff9r.jpg)
,
![梁單元-有限元分析ansys培訓課程圖片20](http://m.09bigdata.com/i/g/20/200c9ac402d4fcbb7acbc8d437113849r.jpg)
,
![梁單元-有限元分析ansys圖片21](http://m.09bigdata.com/i/g/4d/4da451a49051b04aa12123f88dad413ar.jpg)
。對該梁進行分析,畫出彎矩圖和剪力圖。
![梁單元-有限元分析ansys圖片22](http://m.09bigdata.com/i/g/50/502e0693ca9e85b5d26527d42d05e99br.jpg)
圖1. 外伸梁結構示意圖
用材料力學計算所得剪力和彎矩圖如下(以供對照):
![梁單元-有限元分析ansys圖片23](http://m.09bigdata.com/i/g/69/6995cc2ad158d7e6847ceb17119605c8r.jpg)
圖2. 外伸梁的剪力圖和彎矩圖
利用有限元軟件ANSYS對此外伸梁及其承載情況進行分析。
步驟概括:
1)創(chuàng)建梁的各個節(jié)點;
2)顯示各個節(jié)點。
2.定義單元類型和材料
1)定義單元類型-2D elastic 3;
2)定義材料特性- Isotropic, 在EX后的文本框內輸入數(shù)值207e5作為彈性模量;
3)定義幾何參數(shù)。
1)創(chuàng)建單元;
2)顯示單元資料。
1)施加集中載荷F;
2)施加彎矩M;
3)施加分布載荷q。
1)改變分析類型;
2)求解。
1)顯示梁單元變形結果;
2)建立元素結果表:
①創(chuàng)建單元表,計算節(jié)點彎矩;
②創(chuàng)建單元表,計算節(jié)點剪力。
3)列出所有表格資料:
①列出資料;
②畫剪力圖;
![梁單元-有限元分析ansys結果圖片24](http://m.09bigdata.com/i/g/63/6338186c286a75bc6a7d759ee29b0306r.jpg)
圖3. 外伸梁的剪力圖
③畫彎矩圖。
![梁單元-有限元分析ansys結果圖片25](http://m.09bigdata.com/i/g/91/91745b127705f3b2f63b49cebea180edr.jpg)
圖4. 外伸梁的彎矩圖
將理論的彎矩圖和剪力圖與ansys軟件分析出的結果是一致的,通過對比,我們可以看出有限元的優(yōu)越性。
相關標簽搜索:梁單元-有限元分析 Ansys有限元培訓 Ansys workbench培訓 ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術咨詢 ansys基礎知識 ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓