dp模型的案例
abaqus二維切削DP本構模型 ￥10
abaqus二維切削DP本構模型
solid186與solid185單元結果對比下載
為克服經典DP模型的一些缺點，ANSYS開發了擴展DP模型，簡稱EDP,
打開ANSYS HELP，查看Plane182/183、Solid185/186支持的材料模型列表，發現已經沒有經典DP材料模型，而與DP材料有關主要有兩個：
1、Drucker-Prager concrete
2、Extended Drucker-Prage
Plane182/183單元支持材料本構模型
Solid185/186單元支持材料本構模型
顯然，若要使用Plane182/183、Solid185/186來模擬巖土，則必須使用Extended Drucker-Prage，也即擴展的DP材料模型，翻開ANSYS對EDP材料模型的介紹，發現其參數與傳統DP材料模型輸入不同，傳統DP模型需要輸入內摩擦角、粘聚力、膨脹角三個參數，而對于EDP材料模型，使用時除了定義屈服函數外還需要定義流動準則，根據不同的屈服函數與流動準則，輸入的參數有一定差異性，如下表所示：
所以若要使用EDP材料模型，就涉及到經典DP模型到EDP模型的參數轉換問題，也即如何通過已知的內摩擦角、粘聚力計算得到EDP所需要的參數數值。
通過對比經典DP模型和EDP模型的屈服函數，發現EDP模型中的線性屈服函數與DP模型的屈服函數形式上相似且屈服面形狀也相同，通過參數等效替換，可得到EDP參數的計算公式如下：
C1=6sin(a)/[3-sin(a)];
C2=6Ccos(a)/[3-sin(a)];
上式中，C1為EDP模型第一個參數，根據ANSYS HELP中的英文名稱，可解釋為應力敏感度，C2為EDP模型的第二個參數，可解釋為屈服強度，a代表已知的摩擦角，C代表已知的粘聚力。
展開 基于ANSYS的樁土分析模型 ￥15
對于土體的單元類型采用ANSYS中提供的SOLID45實體單元類型，它是一種三維六面體單元，可用于建立各向同性固體力學問題的模型。SOLID45實體單元有8個節點，每個節點有沿X、Y、Z三個方向的平移自由度，在單元的各個側面可施加分布式載荷。在求解分析大位移、大應變、塑性和屈服等方面的問題時，SOLID45單元求解的輸出結果包括節點位移，各個方向的主應力、正應力、剪應力及總應變等。
土體的本構模型采用ANSYS中提供的Drucker-Prager模型，簡稱DP模型，該模型對MC模型的屈服面函數作了適當的修改并且考慮了體積力對屈服的影響，易于程序的編制和進行數值計算，可用于顆粒狀的材料，例如：土壤、巖石、混凝土等[34][41-43]。除了DP模型以外，土體的本構模型還有線彈性模型、DC模型、MC模型等。線彈性模型遵從胡克定律，只有兩個參數，只是簡單的應力應變關系，無法描述土的很多特征；DC模型是一種非線性彈性模型，只是單純的采用了彈性理論，而未曾涉及到塑性理論，著重于對應力-應變簡單的描述，因而沒有反映出土體的很多重要性質，例如土體的剪脹性、球應力對剪應變的影響等[47,48]；MC模型是一種彈-理想塑性模型，采用了彈塑性理論，涉及到了土體的五個參數，能夠較好的描述土體的破壞狀態，但沒有考慮到應力歷史的影響及區分加荷與卸荷[45,46]。
混凝土單元類型采用ANSYS中的SOLID65實體單元類型，它是在SOLID45的基礎上專門開發出來用于建立鋼筋混凝土或混凝土材料問題的有限元模型。同樣，它也有8個節點，每個節點同樣有沿X、Y、Z方向的三個平移自由度，主要用于單元受壓破碎、受拉開裂等問題方面的模擬分析 [38]。混凝土的本構模型同樣采用DP模型，定義其參數。
附件包括一個分析文檔，另有兩個a樁和C樁的建模分析流程。
展開 通過UMAT實現基于DP屈服準則的改進西原模型的三維粘彈塑性（蠕變）本構模型
傳統西原模型是目前可以比較好地描述巖石蠕變過程曲線的元件模型，但是，西原模型使用的元件為黏彈、黏塑性元件（如圖1），難以描述巖石屈服破壞后進入加速階段的蠕變變形。滑坡預報，特別是臨滑預報在地質災害防治領域具有重要意義。
通過編寫abaqus UMAT子程序，可得到如下結果：
（1）應力狀態較小時，僅發生彈性應變和粘彈性應變，最后隨時間趨于穩定值。
（2）單元屈服時，發生粘彈塑性應變，應變隨加載時長逐漸增加，但尚未達到觸發應變，曲線呈現兩階段特性。
（3）隨著加載時長的增加，應變進一步增加，超越觸發應變后，進入快速蠕變階段，應變快速增加，曲線呈現三階段蠕變特性。
參考文獻：
[1] 齊亞靜, 姜清輝, 王志儉, 等. 改進西原模型的三維蠕變本構方程及其參數辨識[J]. 巖石力學與工程學報, 2012, 31(2): 347-355.
[2] 沈才華, 張兵, 王媛, 等. 基于DP屈服準則的西原本構模型及其運用[J]. 地下空間與工程學報, 2016, 12(2): 402-407.
最后，有需要歡迎通過微信公眾號聯系我們。
微信公眾號：320科技工作室。
展開 
ABAQUS不收斂的原因一
因此,相對來說,建議考慮使用DP模型。尤其是ABAQUS中有非線性的DP模型,可以在一定程度上克服線性DP在剪拉區面積過大的情況,因此可以在一定程度上減小開挖卸載時土體的回彈。
一般材料相關的不收斂提示基本就是類似以下:
The plasticity/creep/connector friction algorithm did not converge at 364 points
碰到這個錯誤,一般直接原因就是材料應變太大了,塑性計算迭代不收斂。但這并不意味著材料參數給得不合適或材料強度太弱,很有可能是你的接觸、約束、荷載或邊界出了問題,導致計算中出現的特別大的位移。或者是初始條件(如初始地應力)出了問題,程序沒有計算,直接就報了這個錯誤。可以在下面位置看到計算不收斂的單元位置:
后處理 -> Tools->Job Diagnostics
不收斂需要單一調每個參數確定哪里有問題。
另外，在使用Soils分析步進行孔壓—應力分析時:
a. 墻裂建議各位注意單位,應力單位最好選擇kPa或MPa,不要用Pa,否則可能會遇到各種無腦錯誤提示。
b. 墻裂建議給定孔壓邊界,否則會遇到DOF. 8極大的情況(孔壓就是第8自由度)。
轉自公眾號——ABAQUS大世界
旨在分享，若侵即刪.
展開 ansys workbench鋼筋混凝土建模方法
首先說明下，比較少接觸鋼筋混凝土的理論分析或試驗，本文主要是一個學習的過程，可能很多說法存在問題，但是本文所提及的模型都是一步一步做過的，數據也是盡可能的準確，如有錯誤，歡迎指正。如果某個模型較多人感興趣，再出一期詳細的。
參考文獻：1、周炬《Ansys Workbench有限元分析實例詳解》2、公眾號：搬磚2號叉會腰3、公眾號：ansys結構院4、ansys官方、YouTube等資料。
本文小結:
1、 Mw或DPC+HSD模型，可以說是官方首推的方法，workbench最適用的方法，其solid185和solid186(混凝土)和reinf單元（鋼筋）完美適合用（workbench 2020r2以后版本推出，鋼筋采用此單元，鋼筋與混凝土節點自動耦合），和《混規》GB50010的本構模型相比，DP模型區分了彈性段，強化段，軟化段，殘余應力段。未屈服前按照彈性材料處理，屈服后根據用戶選擇的HSD模型進行計算。中國規范中在峰值拉壓應變前后本構模型為冪函數，HSD模型中的Expotential HSD和中國規范為接近，實際中既可以采用指數函數的HSD也可以采用線性的HSD來進行計算。方法1是王新敏老師推薦的方法。
2、損傷-塑性微平面模型（CPT215單元）在模擬混凝土軟化、下降段方面，優于solid65（壓根就沒有），Mw或DPC（通常采用solid185、186），之前看到一個消息，說官方不建議在wb中使用，但是我用WB2024R1測試，沒啥問題，可以與renif單元聯合使用，相比方法1，需要在WB中插入命令流。其最大的缺點是涉及12個參數，很容易調對一個行為，另一個行為出錯。
展開 ansys模擬鋼管混凝土
然后再來一個
MP,PRXY,1,0.20
MP,DENS,1,25E3
TB,MELA,1,0,9
TBTEMP,0
TBPT,,0.000100,3.3e6
TBPT,,0.000500,8.82e6
TBPT,,0.001000,15.11e6
TBPT,,0.001500,18.89e6
TBPT,,0.002000,20.15e6
TBPT,,0.002500,19.31e6
TBPT,,0.003000,18.47e6
TBPT,,0.003500,17.63e6
TBPT,,0.004000,16.79e6
我看，覺得沒有理由這樣做啊，希望大家給我答疑啊
===========
我這幾天做了帶約束拉桿的方鋼管在水平荷載和垂直荷載下分析，
當solid65只定義concr破壞準則時，打開壓脆，可以收斂，峰值和試驗曲線基本相同，力－位移關系與試驗曲線相差較大，
當同時定義本構關系，
我嘗試用了2種，一種是隨動強化模型，一種是dp模型
采用mkin模型時，本構關系采用了honstad本構關系，也試了韓林海老師的方鋼管混凝土的約束混凝土本構，都試了考慮下降段和不考慮下降段的情況，
對dp材料模型，利用國內文獻上出現的兩種不同的計算c1,c2,c3的方法（計算結果相差較大）
這些嘗試都沒有達到最終的收斂效果，我參考的試驗算例的變形位移角比較大，達到6%。我基本上沒辦法得到收斂的結果，我采用的是位移控制加載
前面的回貼好像已經有人做出來，而且效果不錯，希望更詳細介紹一下，謝謝
==========
1、solid45單元本身不帶考慮開裂與壓碎功能的，所以您打開和關閉壓碎、拉裂，對他來說，對沒有任何意義，不信，您可以檢查一下您的結果，壓根就看不到任何的裂縫！
展開 最新項目需求丨在線派單，先到先得！
點擊鏈接搶單：https://www.yqgqt.org.cn/requirement/detail/8273
項目九：
預算范圍：3000搶單中
發布時間：2023-02-24
工期：一周左右
使用的軟件： matlab
需求描述：
使用matlab編程 一個表面隨機生成微凸體 另一個表面也隨機生成微凸體 使兩個表面進行接觸 期間微凸體之間會發生正接觸和側接觸 需要判斷哪個先接觸 接觸的樣式是什么
點擊鏈接搶單：https://www.yqgqt.org.cn/requirement/detail/8273
項目十：
預算范圍：400搶單中
發布時間：2023-02-23
工期：15天
使用的軟件： FLAC
需求描述：
FLAC模擬飽和土體強夯，提供FLAC命令流就行，可以采用二維模型。土體半徑50m，第一層土厚度5m，第二層土厚度15m，第一層土體滲透系數10-3m/s，第二層土體滲透系數10-7m/s，土體都是飽和的，土體本構模型可以采用DP模型或者MC，或者MCC。落錘速度6m/s先錘擊一次，然后提起來再以6m/s錘擊同一位置，要能計算出土體內部孔壓增長和土體累積變形。
點擊鏈接搶單：https://www.yqgqt.org.cn/requirement/detail/8270
如何開通項目接單？
第一步：點擊個人主頁>修改>完成職場信息和教育信息
第二步：返回主頁>項目接單>立即開通
大家在接單過程中遇到的一些問題或者想了解更多項目需求，掃碼添加客服，客服將會為您詳細解答~
展開 免費視頻 | ABAQUS土木行業高級培訓
課程內容
第一講：巖土綜述
課程內容介紹
介紹巖土工程系列課程內容
巖土綜述
介紹巖土工程概況、BIM與CAE的關系
第二講：緒論
介紹巖土工程經典與現代設計方法、物理實驗與數值分析、現實本構理論要求
第三講：物理實驗
介紹物理試驗測試、試驗要求與本構模型參數的標定
第四講：本構模型（I）
MC模型
介紹應力空間和應力不變量、彈性理論、土體的塑性行為、Mohr-Coulomb模型本構
DP模型
介紹擴展的Drucker-Prager模型
第五講：本構模型（II）
DPC模型
介紹修正的Drucker-Prager/Cap模型
Clay模型
介紹臨界面的塑性本構模型
Jointed模型
介紹節理材料本構模型
第六講：多孔介質分析
理論
介紹基本假定和有效應力、應力平衡和流動連續性、分析的類型與用法
例子
介紹飽和問題與部分飽和問題，講解具體示例
第七講：模擬技術
介紹單元技術、地應力狀態、無限域理論與用法
介紹孔隙流體域表面的相互作用、單元添加與移動、材料磨損與網格自適應技術、水力壓裂等應用
擴展學習：混凝土和鋼筋混凝土
介紹無鋼筋混凝土力學行為、Abaqus中的混凝土本構模型
介紹加強筋（鋼筋）在Abaqus中的應用、應用實例演示
如何領取
展開 【STKO/OpenSEES】STKO在邊坡工程中的應用
Drucker-Prager模型（有時也稱作廣義von Mises（米賽斯）模型）修正了Mohr-Coulomb屈服函數，消除了由尖角造成的奇異點。不同于Mohr-Coulomb模型，Drucker-Prager屈服面的拐角是光滑的，在主應力空間中成圓錐形。和MC模型類似，DP屈服面也取決于有效平均應力σm。在該有限元分析計算軟件中使用的DP模型只適用于三軸拉伸實驗，即其投影在偏應力平面上的屈服軌跡外接于Mohr-Coulomb六邊形的內角點（θ = -300），其中θ為Lode角。
2.3 單元剖分
本文使用了Quad結構化單元。在有限元分析中，實際劃分出來的網格很多是不規則形狀的單元且形狀各不相同，造成位移插值函數也不相同，求解計算效率較低。通過將形狀不規整的實際單元和形狀規整的標準單元之間建立坐標映射關系，在標準單元上去研究它的性質，卻不改變原來的問題，映射稱為參數映射，單元稱為參數單元，如果坐標變換和試函數里的形函數、插值節點完全相同，則變換就稱為等參數變換，參數單元稱等參數單元。
2.4 邊界條件
邊坡底部固定，兩側鉸支，ConstraintHandler對象確定如何在分析中強制執行約束方程。約束方程強制執行 DOF 的指定值，或 DOF 之間的關系。邊坡受到重力作用，在quad單元中進行設置。
2.5分析
在 OpenSees 中，分析是一個由組件對象聚合而成的對象。定義在模型上執行的分析類型的是組件對象。靜力分析。輸出任意節點的位置。
展開 勘察設計領域的工程仿真方案
ANSYS中可以考慮不同的土的本構模型，包括Drucker－Prager模型，Mohr－Coulomb模型，Extended DP模型，CAP模型等。
地下工程
地下工程最大的特點就是它的賦存環境和工程本身的不確定因素太多，這包括工程區域的地質構造、圍巖產狀和巖性、巖體力學參數、地應力場構成、地下水分布以及開挖支護等施工因素的不確定性等。正是由于這些不確定因素的存在，使得每一個地下工程都具有自己的特殊性，尤其是像大型水電站地下廠房洞室群這樣的地下工程，它的洞室規模巨大，尺度和斷面形狀各異、布置復雜、縱橫交錯，地質環境千差萬別，因而每一個新的地下水電站洞室群，對設計和施工都是新的挑戰。
砼壩溫度場與應力場全過程仿真模型
隧道工程設計
隧道結構的動靜力學計算是一項比較困難的課題。地層巖土介質和隧道結構相互作用過程相當復雜。只有具有規則幾何形狀和理想的材料特性，且載荷形式與邊界條件是簡單的線彈性體系，才能得到較為精確的解答。但是，對于非線性巖土體內的連續或不連續介質和任意幾何外形的隧道結構，其力學計算必須借助有限元法，采用ANSYS可模擬隧道的開挖、支護及連續施工過程模擬。
水工結構工程
隨著近年來國家基本建設力度的加大,工程建設中遇到的邊坡穩定問題相應的增多。事實上，在邊坡的施工及運營過程中，邊坡失穩現象時有發生，且邊坡一旦失穩，其往往帶來嚴重后果。因此，對邊坡進行模擬和分析判斷其穩定性，并制定出有效的防護與治理措施，具有十分重要的理論與工程實踐意義。
展開 
abaqus顯式動力學在煤層頂板垮落中的應用
abaqus中的剪切損傷模型+DP塑性準則可以很好的模擬煤層頂板的垮落破壞，但計算時容易產生不收斂問題。
abaqus中的顯式動力學不存在不存在收斂性問題。
所以，剪切損傷模型+DP塑性準則+顯式動力學=完美的煤層頂板仿真模擬。
模擬效果如下：
視頻連接：
仿真任務說明書：
計算任務描述：模擬在煤層巷道開挖后，煤層頂板的垮落情況。
仿真計算所采用的設備基本情況：華碩四核筆記本電腦，酷睿i5處理器，2G內存。
計算模型的處理技術：將實際的三維地質情況，假定為平面應變彈性力學問題；對煤層頂板的力學行為假定為彈性-DP塑性-剪切損傷。
計算機的耗時情況：5min。
仿真結果：在煤層巷道開挖后，煤層頂板可以完全垮落，放頂采煤效果100
展開 ANSYS常見運行錯誤及其解決辦法
主要是在一些邊線的節點會有載荷，邊界條件的約束同時施加在上面，一個會覆蓋掉另外一個 一般忽略該警告
NO.0037
solid model data is contaminated
實體模型被污染了（布爾操作中經常出現）
原因就是布爾操作中出現運算錯誤，實體模型被污染。
解決辦法：１、修補模型。
２、最好的還是重新建立模型。
關于這個，還有一種方法是將模型用write命令寫出來，然后用read命令讀進去就好了
NO.0038
建了一個三維結構的模型模型中土體材料使用ansys中的dp模型，可是每次一使用dp模型時，計算怎么也不收斂，我用mcheck命令檢查了一下單元，出現了類似以下的好多warning：
*** WARNING ***
SUPPRESSED MESSAGE
CP =
4.828
TIME= 20:58:40
Element 3526 node 5851 is part of at least 2 distinct sets of exterior
element faces.
This may indicate that the attached elements are
connected in an unusual manner.
Element 3526 node 5851同時是兩個（及以上）單元外表面的一部分。這意味著，相互連接的單元式以一種不尋常的方式連接起來的。
也就是說，網格質量過差，導致不收斂。
展開 關于ANSYS中初始地應力施加方法的介紹
答：為了使模型更加的符合實際。
ANSYS中對于初始地應力的平衡沒有類似設計軟件那般（例如Midas NX）方便，也即是我們在做類似基坑開挖分析之前，首先要進行地應力平衡操作，ANSYS中主要分為兩步進行：
第一步：原始應力場計算，導出地應力文件；
第二步：新建模型，導入地應力文件，施加重力，平衡地應力。
關于導入與導出的命令流，ANSYS以前老版本是采用Iswrite與isfile命令，新版本可采用Inistate命令，但是Iswrite和isfile依然可以用，只是幫助文件已經沒了這兩個命令的解釋，其用法同inistate，可具體查詢Help。
通過上兩步操作，能達到位移清零，真實應力不清零的效果，下面以一個小例子來進行說明操作過程。
某二維地塊，長度50m，高度20m，需進行基坑開挖操作，操作之前，需進行初始地應力的平衡，材料為中風化砂質泥巖，彈性模量取1200Mpa,粘聚力取450Kpa，摩擦角取30度，采用DP材料模型。
為詳細對比每一步操作過程，本次分了三個模型進行操作演示：
模型一：初始應力場計算，導出初始應力文件；
模型二：導入初始應力文件，不施加重力荷載；
模型三：導入初始應力文件，施加重力荷載。
一、模型一
只施加重力荷載，求解命令流如下，
!============
!建模
!
展開 ANSYS常用命令
(14) Tb, lab, mat, ntemp,npts,tbopt,eosopt 定義非線性材料特性表
Lab: 材料特性表之種類
Bkin: 雙線性隨動強化
Biso: 雙線性等向強化
Mkin: 多線性隨動強化(最多5個點)
Miso: 多線性等向強化（最多100個點）
Dp: dp模型
Mat: 材料號
Ntemp: 數據的溫度數
對于bkin: ntemp缺省為6
miso: ntemp缺省為1，最多20
biso: ntemp缺省為6，最多為6
dp: ntemp, npts, tbopt 全用不上
Npts: 對某一給定溫度數據的點數
?(15) TBTEMP,temp,kmod 為材料表定義溫度值
temp: 溫度值
kmod: 缺省為定義一個新溫度值
如果是某一整數，則重新定義材料表中的溫度值
注意：此命令一發生，則后面的TBDATA和TBPT均指此溫度，應該按升序
若Kmod為crit, 且temp為空，則其后的tbdata數據為solid46,shell99,solid191中所述破壞準則
如果kmod為strain,且temp為空，則其后tbdata數據為mkin中特性。
? TBDATA, stloc, c1,c2,c3,c4,c5,c6
給當前數據表定義數據（配合tbtemp,及tb使用）
stloc: 所要輸入數據在數據表中的初始位置，缺省為上一次的位置加1
每重新發生一次tb或tbtemp命令上一次位置重設為1，
（發生tb后第一次用空閑此項，則c1賦給第一個常數）
?
展開 