Abaqus課程設計論文

1、學校代碼：10128學 號：201220907039 課程論文題 目：球體接觸問題的有限元分析學生姓名：李勁波 學 院：理學院班 級：工程力學12-1班 指導教師：韋廣梅、周承恩 2015年9月26日內蒙古工業(yè)大學課程設計（論文）任務書課程名稱： 有限單元法課程設計 學院：理學院 班級： 力學12-1班 _ 學生姓名： 李勁波 學號：201220907039 指導教師： 韋廣梅、周承恩 一、題目球體接觸問題的有限元分析二、目的與意義有限單元法課程設計是“有限單元法”課程的課程設計，本課程設計的目的是使學生應用和拓展所學“有限單元法”課程的理論，使學生學會用大型商用有限元軟件進行有限元分析，為畢

2、業(yè)論文的完成和畢業(yè)后從事專業(yè)工作做好計算分析的準備。 三、要求（包括原始數(shù)據(jù)、技術參數(shù)、設計要求、圖紙量、工作量要求等） 1、內容：橢球與平面接觸2、原始數(shù)據(jù)：1.半無限大體，橢球三半軸為50，50，70mm。50mm半軸與半無限大體表面垂直。接觸壓力為30kN。 材料參數(shù)：E=210 Gpa，=0.33、要求： 1.建立有限元模型，研究接觸區(qū)域形狀，應力分布。求二者內部最大的應力及位置。2.如果材料為理想彈塑性，屈服極限350MPa，分析接觸區(qū)域形狀及應力應變分布情況。四、工作內容、進度安排 1、熟悉有限元商用軟件（1周）。2、在有限元商用軟件平臺上對題目進行設計（1周）。3、結果分析、報告

3、形成（1周）。五、主要參考文獻1、莊茁，基于ABAQUS的有限元分析和應用，清華大學出版社，20092、王勖成。有限單元法（FINITE ELEMENT METHOD）北京：清華大學出版社，2005；3、梁力，宋瑞強等。新型掛梁結構恒載力學特性的有限元數(shù)值分析，遼寧：自然科學出版社，2002第23卷。4、劉齊茂，李微?；谟邢迒卧ǖ匿?混凝土組合梁截面優(yōu)化設計 西安建筑科技大學學報，2005年12月，37審核意見系（教研室）主任（簽字） 指導教師下達時間 年 月 日指導教師簽字：_摘 要許多工程問題都涉及兩個或多個部件之間的接觸，在這些問題中，當兩個物體彼此接觸時，垂直于接觸面上的力作用在兩

4、個物體上。如果在接觸面之間存在摩擦，可能產生剪力以阻止物體的切向運動。本文主要是運用ABAQUS軟件對橢球體與剛性半無限大體接觸問題進行探討，因為橢球體是前后左右都對稱的三維實體，所以對它取四分之一分析，這樣不僅計算方便而且結果直觀。先將材料定義為彈性材料通過前處理建模、計算、后處理來分析橢球體與大體之間的接觸，運用多種方法討論接觸問題。結果得到的數(shù)據(jù)遠超彈性屈服極限，后將材料重新賦值為彈塑性，仿照彈性的分析方法對彈塑性進行分析，得到的結果完全正確。關鍵詞：接觸；ABAQUS；彈塑性；屈服極限AbstractMany engineering problems involve contact b

5、etween two or more components, in these problems, when two objects come into contact with each other, on the surface of the perpendicular to the contact force between two objects. If there is friction between contact surface, shear may be produced to prevent the tangential motion of the object. This

6、 article mainly using ABAQUS software to ellipsoid and rigid semi-infinite discusses the general contact problem, because the ellipsoid is before and after the left and right sides is symmetrical three-dimensional entity, so take a quarter analysis to it, that not only the calculation results of con

7、venient and intuitive. Before the adoption of the first material is defined as elastic material modeling, calculation, after processing to analyze between ellipsoid and general contact, using a variety of methods to discuss contact problem. Than elastic yield limit, the resulting data will reopen as

8、signment for elastic-plastic materials, imitates the elastic analysis method of elastic-plastic analysis, get the right result.Key words: Contact; ABAQUS; Elastic-plastic; Yield limit 目 錄第一章 緒論1.1 有限單元法課程設計目的與任務11.2 接觸問題概述11.3 有限元法概述1第二章 有限元工具ABAQUS介紹及應用方法32.1 ABAQUS簡介32.2 ABAQUS功能簡介3第三章 橢球平面接觸問題有限元

9、分析7 3.1 平面球接觸問題描述7 3.2 用ABAQUS建立有限元模型7 3.2.1 建立實體7 3.2.2 賦值材料與裝配部件8 3.2.3 創(chuàng)建分析步8 3.2.4 定義接觸和邊界條件9 3.2.5 劃分網(wǎng)格93.2.6 計算及后處理103.2.7 理想彈塑性分析11第四章 結論與收獲14參考文獻15第1章 緒論1.1有限單元法課程設計目的與任務目的：有限單元法課程設計是“有限單元法”課程的拓展部分。本課程設計的目的是使學生結合所學的力學知識和基礎理論，應用和拓展“有限單元法”課程，初步學會用大型商用有限元軟件進行力學分析，為以后從事專業(yè)工作做好計算分析的準備。任務：使用有限元商用軟件

10、進行力學分析。1.2 接觸問題概述彈性力學中研究的兩接觸物體受壓力后產生局部應力和變形的問題的統(tǒng)稱。軸承、凸輪機構、齒輪、硬度計、軋鋼機的軋輥、橋梁支座和剛性壓頭等在使用中都有接觸問題。接觸問題曾是應用數(shù)學力學家面臨的一大難題。在過去一百多年的研究中，產生了一系列有效的數(shù)學方法。 材料在接觸區(qū)的變形受到各個方面的限制，因而處于三向應力狀態(tài)。接觸應力具有明顯的局部性質，而且總是隨著離開接觸區(qū)距離的增大而迅速衰減。一般在接觸表面中心的壓應力最大。影響接觸應力的因素很多，其中主要的有：殘余應力（由于接觸點附近應力很大，有時會進入塑性狀態(tài)而引起殘余應力）；熱應力（兩接觸面相對滑動時摩擦會引起熱應力）；

11、潤滑，它影響熱應力及切向載荷，還可能產生動壓油膜而影響接觸應力的大小及分布；接觸面的幾何形狀偏差等。目前，在解決彈性體接觸問題方面已采用有限元法，通過它可計算各種形狀、材料和載荷下的接觸問題，并且所得的結果精度較高。1.3 有限元法概述有限元分析（FEA，F(xiàn)inite Element Analysis）利用數(shù)學近似的方法對真實物理系統(tǒng)（幾何和載荷工況）進行模擬。還利用簡單而又相互作用的元素，即單元，就可以用有限數(shù)量的未知量去逼近無限未知量的真實系統(tǒng)。有限元分析是用較簡單的問題代替復雜問題后再求解。它將求解域看成是由許多稱為有限元的小的互連子域組成，對每一單元假定一個合適的(較簡單的）近似解，然

12、后推導求解這個域總的滿足條件(如結構的平衡條件），從而得到問題的解。這個解不是準確解，而是近似解，因為實際問題被較簡單的問題所代替。由于大多數(shù)實際問題難以得到準確解，而有限元不僅計算精度高，而且能適應各種復雜形狀，因而成為行之有效的工程分析手段。14第二章 有限元工具ABAQUS介紹及應用方法2.1 ABAQUS簡介ABAQUS是一種有限元素法軟件，用于機械、土木、電子等行業(yè)的結構和場分析。ABAQUS早年屬于美國HKS公司的產品，于2000年代中期賣給了達索公司，該軟件又被稱為達索SIMULIA。ABAQUS非常適合用作科學研究。作為一種功能強大的有限元分析軟件，ABAQUS在商業(yè)有限元軟件

13、中占有極其重要的位置。從簡單的線彈性問題到復雜的幾何非線性和材料非線性問題均獲得了廣泛應用。其有效性不論在工程應用還是在科學研究方面均得到了驗證。ABAQUS包含了豐富的單元庫和材料庫，能夠模擬各種材料受力和變形行為，特別是其提供了UMAT和UEL，為研究者提供了開發(fā)平臺。2.2 ABAQUS功能簡介 ABAQUS有兩個主求解器模塊 ABAQUS/Standard 和 ABAQUS/Explicit。ABAQUS 還包含一個全面支持求解器的圖形用戶界面，即人機交互前后處理模塊 ABAQUS/CAE 。 ABAQUS 對某些特殊問題還提供了專用模塊來加以解決。ABAQUS 被廣泛地認為

14、是功能最強的有限元軟件，可以分析復雜的固體力學結構力學系統(tǒng)，特別是能夠駕馭非常龐大復雜的問題和模擬高度非線性問題。 ABAQUS 不但可以做單一零件的力學和多物理場的分析，同時還可以做系統(tǒng)級的分析和研究。 ABAQUS 的系統(tǒng)級分析的特點相對于其他的分析軟件來說是獨一無二的。由于 ABAQUS 優(yōu)秀的分析能力和模擬復雜系統(tǒng)的可靠性使得 ABAQUS 被各國的工業(yè)和研究中所廣泛的采用。 ABAQUS 產品在大量的高科技產品研究中都發(fā)揮著巨大的作用。1. Part模塊 ABAQUS/CAE:Part modulePartCreate，可以創(chuàng)建Planar和Axisym

15、-metric類型幾何模型，幾何模型的類型可以是變形體和剛體，幾何屬性可以為實體、殼、線和點。這里值得注意的是軸對稱問題被單獨提出來進行分析，而平面應力和平面應變問題則被包含在2D Planar中。 生成Part后，為了對Part不同部分賦予不同材料參數(shù)、接觸的定義、規(guī)則網(wǎng)格劃分等原因，往往需要對Part進行切割，切割一般需要進行兩步操作才能完成，首先是定義切割平面，定義方式:Part moduleToolsDatum進行定義，對于切割平面很多時候是采用Offsetfrom principal plane(xy plane, xz plane和

16、yz plane),Offset from plane和3 points。工作平面定義完成后就可以進行切割了，切割的操作方式為:Part moduleToolsPartition. 2. Property模塊 該模塊主要用于定義材料的本構模型，ABAQUS中封裝了大量可用于土木工程的本構模型，通過PropertyMaterialCreateMechanical進行定義。完成力學定義行為并給定材料名稱后，接著通過PropertySectionCreate定義斷面(section)，建立與材料屬性相關聯(lián)的斷面，然后通過Prope

17、rtyAssignSection對不同Part或Part的不同部分(可定義為set)進行材料賦值。這里需要強調的是對干梁單元需要通過PropertyProfileCreate建立梁的幾何斷面，同時還要通過PropertyAssignBeam Section Orientation對梁單元斷面方向進行定義。 3. Assembly模塊 任何一個結構都可以視為一個Instance，它由很多個Part構成，Part間連續(xù)、接觸或通過螺栓等形式連接。通過AssemblyInstanceCreate建立由不同Part構成的In-stance，可以對Ins

18、tance中的Part進行移動、旋轉等操作，更主要的是可以對Instance中的Part進行布爾運算，如AssemblyInstanceMerge/Cut，在合并操作中對相交邊界有刪除(Remove)和保留(Retain)兩個選項，后一選項在很多時候用得更多，因為進行Merge/Cut操作后，先前的多個Part生成了一個新的Part，選擇“保留”選項后可以進人到Property模塊對同一Part不同部分賦予不同材料參數(shù)。4. Step模塊 進人StepCreate，分析類型(Procedure type)有兩類，其一為Linear perturbation,包括B

19、uckle,Frequency,Static,Linear pertubation和Steadystate dynamics,Direct四部分，應用更多的是頻率的計算和振型的提??；其二為General，幾乎包含了所有的分析類型，具體包含Coupled temp-displacement，用干溫度位移場的禍合分析，如熱輻射使結構溫度升高，因變溫而使受約束結構產生溫度應力；Coupled thermal-electric，主要用于壓電材料的分析;Dynamic, Implicit，隱式動力分析；Dynamic, Explicit，顯式動力分析；D

20、ynamic, Temp-disp, Explicit，顯式溫度位移分析;Geostatic，地應力場計算，在巖土工程中地應力場的準確確定決定了后繼計算結果的有效性；Heat transfer，熱傳導計算分析；Mass diffusion；Soils，土的固結分析；Static, General，一般的靜力分析；Static,Riks，弧長法靜力分析；Visco，粘彈性及蠕變分析。荷載步中增量步分為兩種:自動增量法和固定時間增量法，同時在荷載步中要給出最大迭代次數(shù)等參量。5. Interaction模塊 結構中各個構件間的連

21、接形式千差萬別，連接形式對結構在外荷載作用下響應的影響至關重要。在該模塊中，主要包含四部分，具體為Interaction, Constraint, Connector, Special oInteractionInteractionProperty，定義接觸面間的接觸屬性，包括定義摩擦系數(shù)等參量;由InteractionInteractionCreate創(chuàng)建由接觸面和目標面構成的接觸對，應用最多的是Surface-to-surface contact接觸問題，如結構物與土之間的接觸問題。由InteractionConstraintCreate,

22、Tie,Coupling,Shell-to-solid coupling,Em-bedded region, Equation等約束條件，其中Tie約束適用于約束面間網(wǎng)格劃分不一致但變形又連續(xù)的情況。而Embedded region功能可能模擬加筋對基體的增強功能，如鋼筋混凝土實體結構(不是采用殼或梁進行模擬)，分別創(chuàng)建鋼筋和混凝土模型，然后將鋼筋埋置到混凝土結構中去。InteractionConnectorCreate，主要用于構件間的U型連接、焊接、鉸接等連接的模擬，這一點在機械工程中具有重要的應用價值。InteractionSpecialInert

23、ia，用于定義結構中的集中質量和慣性矩，這主要應用在懸索橋梁及高層建筑的地震分析中。InteractionSpecialCrack，用于定義結構中的裂紋及裂尖，為后繼網(wǎng)格劃分生成奇異單元奠定基礎。InteractionSpecialSpring/Dash-pot，定義彈簧和阻尼器。6. Mesh模塊網(wǎng)格劃分前需要對Part定義種子(Seeds)，通過MeshSeedEdge By Number (Edge By Size, Edge Biased)完成種子的定義，對幾何模型的邊可以均勻劃分，也可按等比數(shù)列進行劃分。

24、具體劃分的數(shù)目和方法要通過數(shù)值模擬試驗確定，最終與理論解或試驗結果進行對比以確定經濟合理的網(wǎng)格劃分規(guī)模。7. Job模塊 通過JobJobCreate創(chuàng)建Job的名字，分析完成后存儲的相關文件的主文件名均為該Job名字。在Submission中包含Job type, Run Mode和Submit Time；在General下主要是User subroutine file的選擇;在Memory下主要是內存的劃分;Parallelization設置并行計算的一些內容；Precision定義問題分析的精度。8.

25、0;Visualization模塊 后處理部分主要包括變形圖顯示(VisualizationPlotDeformed shape)，輸出變量的云狀圖(VisualizationPlotContours).輸出結果包括場變量輸出(Visualization-Result-Field output)和與時間有關的變量的輸出(VisualizationResultHistory output).也可以將后處理結果導出到外部文件中去(Visualization ReportXY, Field output)，以便進一步編輯和處理。&

26、#160;9. Sketch模塊 通過該模塊可以生成輪廓線或由外部文件導人生成輪廓線，對于土木工程來說，用得比較多的是Autocad繪圖軟件，在Autocad中生成的二維線轉存為擴展名為DXF的文件，然后在ABAQUS/CAE下選擇FileImportSketch就可導人到Sketch模塊中。在創(chuàng)建Part時進人到Sketch界面下，由SketchAddSketch就可以將先前由Autocad導人的Sketch顯示于當前界面下，進而創(chuàng)建所需要的Part。第三章 橢球平面接觸問題有限元分析3.1 橢球平面接觸問題描述橢球三半軸為50，50，70mm。50mm半軸與半無限大體表

27、面垂直。接觸壓力為30kN。 材料參數(shù)：E=210 GPa，=0.3。用ABAQUS建一個橢球與相對較大的半無限解析剛體如圖3-1，因為橢球是前后左右對稱的圖形，所以可以取四分之一進行接觸分析。圖3-1 模型圖3.2 用ABAQUS建立有限元模型3.2.1 建立實體在part模塊繪制四分之一橢球（圖3-2）和相對大的半無限大體（圖3-3） 圖3-2 橢球 圖3-3 大體3.2.2 賦值材料與裝配部件在Property模塊對橢球進行材料賦值。先是創(chuàng)建一個材料命名為steel，材料參數(shù)：E=210 GPa，=0.3的彈性材料(圖3-4），然后賦值給橢球體。因為大體為彈性剛體所以不能賦值材料默認即可

28、。進入Assembly 模塊，創(chuàng)建實體橢球和大體，通過平移使橢球與大體接觸（圖3-5）3.2.3 創(chuàng)建分析步 圖3-4 創(chuàng)建材料 進入step模塊，創(chuàng)建兩個分析步step-1、step-2。對step-1選擇默認，在step-2增量步中選擇最小迭代為0.1(圖3-6）。然后對Field output Requests Manager(場變量輸出管理器)和 History output Requests Manager(歷史變量輸出管理器)進行設置輸出所需要的數(shù)據(jù)。 圖3-5 裝配實體 圖3-6 創(chuàng)建分析步3.2.4 定義接觸和邊界條件 在Interaction模塊中定義大體上表面和橢球曲面小滑

29、移接觸，然后在load模塊中對配件進行邊界定義，創(chuàng)建四個邊界條件如圖3-7所示，對橢球的兩個截面和大體的參考點進行約束，然后在上頂點施加一個向下0.1的位移（BC-4），最后對橢球施加一個向下的集中力得到如圖3-8所示的圖形 圖3-7 定義邊界條件 圖3-8 邊界條件圖3.2.5 劃分網(wǎng)格進入Mesh模塊，對橢球體進行單元網(wǎng)格劃分。本題主要是研究橢球與大體接觸部分的應力、應變狀態(tài)，所以應該對接觸部分進行細分，而其他不重要的部分粗分，這樣可以使要研究對象能夠得到最大限度的精度，也節(jié)約了計算時間和降低了復雜度。首先用seed part 對橢球進行粗糙劃分，然后再用seed Edges對其采用single的細劃分，最后通過Assign Mesh、Mesh part的處理得到劃分好的橢球體（圖3-9）。圖3-9 劃分網(wǎng)格3.2.6 計算及后處理上面的步驟完成之后進入Job模塊，創(chuàng)建一個job然后對它進行Submit及Monitor得到橢球與大體接觸的有限元計算結果如圖3-10，然后點擊Results進入Visualization模塊，得到彈性材料Tresca云圖（圖3-11），通過放大下頂點得到接觸處的彈性Mises云圖（圖3-12）圖3-10