本文選題：地下工程 + Mohr-Coulomb模型�。� 參考：《華南理工大學(xué)》2015年碩士論文

【摘要】：地下工程中涉及的材料主要是巖土類材料,巖土的性質(zhì)非常復(fù)雜,屬于三相體摩擦類材料。由于地質(zhì)體具有不連續(xù)性、非均質(zhì)性、各向異性、存在地下水和初應(yīng)力等特征,因此很多的巖土工程問題很難通過解析法求解,為此可以采用有限元等數(shù)值方法來解決;巖土具有明顯的彈塑性特征,須進(jìn)行彈塑性分析才能得到切合實(shí)際的計(jì)算結(jié)果�？紤]以上因素并為了拓寬有限元軟件ANSYS在地下工程的運(yùn)用范圍,本文利用ANSYS的二次開發(fā)功能—用戶可編程特性UPFs,編制Mohr-Coulomb模型子程序usermat.F,對巖土結(jié)構(gòu)進(jìn)行材料非線性有限元分析。本文首先,在理解彈塑性理論中平面應(yīng)變問題、屈服條件、硬化模型、流動法則、加卸載準(zhǔn)則等基本概念的基礎(chǔ)上。通過Mohr-Coulomb定律,導(dǎo)出用應(yīng)力不變量表示的屈服準(zhǔn)則,考慮各向同性硬化模型和相關(guān)聯(lián)流動法則,采用簡單的處理方法消去無法確定塑性流動方向的奇異點(diǎn),推導(dǎo)了便于數(shù)值計(jì)算的彈塑性剛度矩陣。然后,根據(jù)usermat子程序有應(yīng)力更新和輸出一致切線算子矩陣的任務(wù)。詳細(xì)闡述了求解應(yīng)力更新的回退算法中彈性預(yù)測步和塑性修正步這兩個(gè)步驟的數(shù)學(xué)和力學(xué)原理;通過增量型有限元說明一致切線算子矩陣即彈塑性剛度矩陣。據(jù)此,采用結(jié)構(gòu)化程序設(shè)計(jì)方法,編制了子程序usermat.F。通過編譯鏈接程序,將Mohr-Coulomb模型子程序?qū)階NSYS材料庫中。最后,通過兩個(gè)算例來驗(yàn)證本文模型的正確性和可靠性。通過算例一證明了子程序usermat.F能完成應(yīng)力更新的任務(wù),并在宏觀上能反映土體從開始受力到剪切坡壞的整個(gè)過程。通過算例二,說明了自編的Mohr-Coulomb模型能達(dá)到ANSYS自帶DP模型相近的效果;另外從安全因素和土所具有的本質(zhì)屬性來考慮,對于地下工程,采用彈塑性分析才能真實(shí)反映巖土的力學(xué)特性。
[Abstract]:The materials involved in underground engineering are mainly geotechnical materials. The properties of rock and soil are very complex and belong to three-phase friction materials. Because geological bodies have the characteristics of discontinuity, heterogeneity, anisotropy, groundwater and initial stress, many geotechnical engineering problems are difficult to be solved by analytical method, which can be solved by finite element method. The rock and soil have obvious elastic-plastic characteristics, so it is necessary to carry out elastic-plastic analysis to obtain practical results. Considering the above factors and in order to widen the application scope of finite element software ANSYS in underground engineering, this paper makes use of the secondary development function of ANSYS-user programmable characteristic UPFs, compiles the Mohr-Coulomb model subprogram usermat.F, carries on the material nonlinear finite element analysis to the geotechnical structure. In this paper, firstly, the basic concepts of plane strain problem, yield condition, hardening model, flow rule, loading and unloading criterion in elastic-plastic theory are understood. By means of Mohr-Coulomb 's law, the yield criterion expressed by stress invariants is derived. The isotropic hardening model and the associated flow rule are considered, and the singularity points in which the direction of plastic flow can not be determined are eliminated by a simple treatment method. The elastoplastic stiffness matrix which is easy to calculate numerically is derived. Then, according to the usermat subroutine, there are the tasks of updating the stress and outputting the uniform tangent operator matrix. The mathematical and mechanical principles of elastic prediction step and plastic correction step in the regression algorithm for stress updating are described in detail, and the uniform tangent operator matrix, i.e. the elastic-plastic stiffness matrix, is explained by the incremental finite element method. Based on this, the subprogram usermat.F. By compiling the link program, the Mohr-Coulomb model subroutine is imported into the ANSYS material library. Finally, two examples are given to verify the correctness and reliability of the model. An example is given to prove that the subroutine usermat.F can complete the task of stress renewal and can reflect the whole process from the beginning to the failure of the shear slope. The second example shows that the self-made Mohr-Coulomb model can achieve the similar effect of ANSYS self-contained DP model, in addition, considering the safety factors and the nature of soil, for underground engineering, The elastoplastic analysis can truly reflect the mechanical properties of rock and soil.
【學(xué)位授予單位】：華南理工大學(xué)
【學(xué)位級別】：碩士
【學(xué)位授予年份】：2015
【分類號】：TU43

【參考文獻(xiàn)】

相關(guān)期刊論文 前10條

1 張柔雷;;巖體層狀材料非線性問題的無迭代全量解法[J];工程力學(xué);1990年03期

2 張吉萍,陳虬,王文星;材料非線性問題的遺傳算法與神經(jīng)網(wǎng)絡(luò)耦合分析[J];貴州工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版);2005年01期

3 劉勇,GlaucioH.Paulino,梁利華;基于彈粘塑性一致切線算子概念的隱式邊界元法[J];中國科學(xué)E輯:技術(shù)科學(xué);1999年06期

4 曲圣年,殷有泉;塑性力學(xué)的Drucker公設(shè)和Ильюшин公設(shè)[J];力學(xué)學(xué)報(bào);1981年05期

5 段云嶺;非線性方程組的解法:局部弧長法[J];力學(xué)學(xué)報(bào);1997年01期

6 茹忠亮;趙洪波;;彈塑性計(jì)算中Mohr-Coulomb準(zhǔn)則屈服面奇異點(diǎn)處理方法[J];河南理工大學(xué)學(xué)報(bào)(自然科學(xué)版);2013年04期

7 洪開榮;;我國隧道及地下工程發(fā)展現(xiàn)狀與展望[J];隧道建設(shè);2015年02期

8 范秀昌，，王秀琴，李明澤;等效應(yīng)力─等效應(yīng)變曲線的實(shí)驗(yàn)測量[J];天津理工學(xué)院學(xué)報(bào);1996年02期

9 馬志圣;關(guān)于正則曲面的定義[J];四川師范大學(xué)學(xué)報(bào)(自然科學(xué)版);1994年04期

10 鄭鑫元;二維彈塑性有限元程序設(shè)計(jì)[J];上海汽輪機(jī);1994年03期

相關(guān)博士學(xué)位論文 前3條

1 由廣明;巖土彈塑性通用數(shù)值平臺及其應(yīng)用研究[D];同濟(jì)大學(xué);2007年

2 嚴(yán)祖文;基于ANSYS二次開發(fā)的樁土相互作用的三維非線性有限元分析[D];天津大學(xué);2007年

3 黃林沖;多孔介質(zhì)隧道圍巖穩(wěn)定性分析的程序研發(fā)與應(yīng)用[D];中南大學(xué);2009年

相關(guān)碩士學(xué)位論文 前9條

1 張文完;基于三剪統(tǒng)一強(qiáng)度理論的ANSYS二次開發(fā)及其應(yīng)用[D];南昌大學(xué);2010年

2 鄧建霞;有自由面滲流分析的加密高斯點(diǎn)單元傳導(dǎo)矩陣調(diào)整法[D];四川大學(xué);2004年

3 王偉;二維彈塑性土層的波動數(shù)值模擬[D];中國地震局工程力學(xué)研究所;2005年

4 劉連民;大型地下洞室抗震分析及ANSYS二次開發(fā)[D];天津大學(xué);2006年

5 凡建偉;柔性大直徑水平受荷樁承載特性及計(jì)算方法研究[D];重慶交通大學(xué);2009年

6 楊永俊;不同地區(qū)黃土動強(qiáng)度特性試驗(yàn)研究[D];西北農(nóng)林科技大學(xué);2009年

7 趙清華;油層中孔隙介質(zhì)本構(gòu)方程的建立[D];哈爾濱工業(yè)大學(xué);2009年

8 羅建飛;基于非線性本構(gòu)模型的地下工程高效有限元算法研究及應(yīng)用[D];華南理工大學(xué);2012年

9 張朋;ANSYS二次開發(fā)及其在地下洞室分析中的應(yīng)用[D];蘭州交通大學(xué);2013年

本文編號：1914116