abaqus屈服變形的案例
一個有意思的材料本構模型設計方案，拉伸變形采用von Mises屈服，壓縮側 cap屈服本構模型設計。
原始文獻：《Mechanical modelling of indentation-induced densification in amorphous silica》
該文章為了模擬非晶態二氧化硅的壓縮力學性能，把拉伸與壓縮分開處理：拉伸側采用熟悉的 von Mises 屈服，壓縮側則切換到 cap 屈服面。這樣的設計，正好對應了非晶二氧化硅在壓痕加載下“既會發生剪切塑性，又會發生永久致密化”的真實特征。
分享這個代碼的主要原因：一方面，它很適合做玻璃、非晶材料、壓痕問題中的壓力敏感塑性分析；另一方面，它也是學習 cap 模型、致密化硬化和隱式本構積分的一個很好的范例。論文結果表明，這一模型能夠較好復現實驗載荷—位移曲線以及壓痕致密化分布，不過需要明確指出的是，當前模型暫時還沒有考慮剪切硬化，因此更適合用于理解“壓痕致密化”這一核心機制，而不是直接覆蓋所有復雜失效問題。作為一份用于科研復現和二次開發的代碼，我覺得它很有參考價值。
展開 Abaqus 采用YLD2000-2D屈服準則的UMAT子程序
Barlat在2003年提出了專門針對平面應力問題的各向異性屈服準則，該屈服準則對于各向異性材料具有很高的精度，得到了廣泛的應用。
YLD2000-2D屈服面示意圖
Yld2000-2d屈服準則由下式給出
其中
矩陣X′和X″的元素分別由柯西應力的下列線性變換獲得
L′和L″的分量由下式求得
積分算法采用徑向返回算法，該方法是穩健和精確的。
當彈性試算超出屈服面時，則需要進行塑性修正
使其滿足
公式9可以通過牛頓法進行迭代求解。
計算的應力應變曲線如下圖所示
B, F. Barlat A , et al. "Plane stress yield function for aluminum alloy sheets—part 1: theory." International Journal of Plasticity 19. 9(2003):1297-1319.
王海波, 萬敏, 閻昱,等. 屈服準則在有限元軟件中實現的正確性驗證[J]. 固體力學學報, 2010, 031(002):173-180.
展開 批量提取Abaqus的節點坐標(初始坐標、指定Step下的變形量、變形后節點坐標) ￥40
<h2>摘要</h2><p>本文介紹如何使用Python腳本二次開發來批量提取ABAQUS輸出數據庫（ODB）文件中指定Step下的Set節點集變形量。通過詳細的步驟說明、代碼示例和圖片展示，您將學會如何使用該腳本，自動化輸出CSV文件包含（Node Label；Step Name、Increment、Step Time，U1，U2）。</p><p>如果還需要按Increment提取每個增量下的變形后的節點坐標的話，在提取變形量的基礎上，與初始坐標進行簡單的計算就可以求得坐標。 （備注：該代碼只提取了x，y方向的變形量）</p><h2>1. 問題描述</h2><p>在工程仿真和分析領域，提取ABAQUS輸出數據庫（ODB）文件中的節點集變形量是一項常見任務。然而，手動提取這些數據是一項繁瑣且容易出錯的工作。因此，需要一種自動化的方法來批量提取指定步驟下按節點集組織的變形量數據。</p><h2>2. 實例展示</h2><p>假設我們有一個名為`example.odb`的ODB文件，其中包含名為`Step-x`的步驟和名為`Set-x`的節點集。運行以上代碼后，腳本會自動將該步驟下節點集的變形量提取出來，并保存為`NodalDisplacement.csv`文件。
展開 針對平面應力問題的YLD2000-2D屈服準則及其在ABAQUS中UMAT子程序的實現
Barlat在2003年提出了專門針對平面應力問題的各向異性屈服準則，該屈服準則對于各向異性材料具有很高的精度，得到了廣泛的應用。
YLD2000-2D屈服面示意圖
Yld2000-2d屈服準則由下式給出
其中
矩陣X′和X″的元素分別由柯西應力的下列線性變換獲得
L′和L″的分量由下式求得
積分算法采用徑向返回算法，該方法是穩健和精確的。
當彈性試算超出屈服面時，則需要進行塑性修正
使其滿足
公式9可以通過牛頓法進行迭代求解。
計算的應力應變曲線如下圖所示
B, F. Barlat A , et al. "Plane stress yield function for aluminum alloy sheets—part 1: theory." International Journal of Plasticity 19. 9(2003):1297-1319.
王海波, 萬敏, 閻昱,等. 屈服準則在有限元軟件中實現的正確性驗證[J]. 固體力學學報, 2010, 031(002):173-180.
最后，有需要歡迎通過微信公眾號聯系我們。
微信公眾號：320科技工作室。
展開 
ABAQUS中mises應力云圖顯示的最大值還不到屈服應力值為啥還有PEEQ值
ABAQUS中mises應力云圖顯示的最大值還不到屈服應力值為啥還有PEEQ值，PEEQ云圖有變形值
Abaqus+PyQt+Python平面變形歐拉角計算
下面以簡單例子介紹平面變形、指向歐拉角的計算，包括絕對歐拉角、相對歐拉角。
1 簡化模型
下面的六面體為表面殼模型，下面由三段梁支持，三段梁分別沿X、Y、Z軸向。六個面的厚度不同，在上側3個面施加不同的壓力，如下左圖所示。位移云圖如下右圖所示。
2 計算要求
計算六面體上面3個面的變形歐拉角，包括3個面的絕對歐拉角，平面2、3相對與平面1的相對歐拉角。平面1、2、3如下圖所示。
3 數據處理
使用平面節點坐標、位移數據計算平面變形歐拉角。可以使用Python腳本輸出平面節點編號、節點坐標（X、Y、Z）、節點位移（U1、U2、U3），如下圖所示。下圖為平面1的10個工況的數據文件，打開的文本文件中7列數據為節點編號、坐標、位移。
三個平面10個工況的節點數據文件如下圖所示。每個文件中包含一個工況一個平面的節點編號、坐標、位移數據。
4 絕對歐拉角計算
使用PyQt+Python開發了一個簡單的小軟件，計算絕對歐拉角、相對歐拉角。
首先計算各平面的絕對歐拉角。
計算平面1的10個工況的絕對歐拉角。
平面1變形的絕對歐拉角計算結果如下圖所示。
伴隨絕對歐拉角計算結果，軟件同時寫出了平面變形前后的坐標系數據，如下圖。每行18個數據，每3個數據為一個坐標軸向量，變形前后2個坐標系，6個坐標軸，18個數據。
5 相對歐拉角計算
利用計算絕對歐拉角時得到的坐標系文件，計算平面變形相對歐拉角，如下圖所示，計算平面2相當于平面1、平面3相對與平面1的相對歐拉角。
計算結果如下圖所示。
6 小結
上述軟件用的算法申請了發明專利，軟件申請了軟著。
展開 ABAQUS CEL(例11) 地震工況下的邊坡大變形模擬 ￥70
ABAQUS CEL(例11) 地震工況下的邊坡大變形模擬
一、建模技術
地震工況下邊坡可能失穩進而出現滑坡現象，為避免模擬滑坡時網格產生的畸變問題，采用耦合歐拉拉格朗日法(CEL)進行滑坡的大變形模擬；土體本構采用摩爾庫倫模型；采用模型底部小范圍內的周期性荷載模擬地震荷載。
二、模型及部分結果展示
圖1：藍色為邊坡；紅色為空氣層
圖2：網格的劃分
圖3：賦予模型初始應力
圖4：土體達到地應力平衡時的應力分布
圖5：土體底部的地震荷載施加區域
圖6：所施加的周期性荷載（地震荷載）
圖7：邊坡因地震荷載產生的位移
圖8：地震波產生的區域
展開 ABAQUS模擬多道次變形的變量繼承方法
一、引言
使用ABAQUS進行多道次加工時，往往牽扯道次之間變量的繼承（如晶粒尺寸、累積損傷等），這對多道次變形模擬結果的準確性有較大的影響。本文以VUHARD子程序及簡單的熱壓縮模型為例，分享雙道次壓縮之間的晶粒尺寸的繼承方法。
基于ABAQUS的橡膠密封圈大變形仿真分析
橡膠圈的材料選取、形狀的設計及受力大小對其密封性能有很大的影響，然而在實際壓縮試驗過程中很難觀測到其受力變形的瞬態大變形行為。通過ABAQUS有限元分析可以得到橡膠圈的受力變形過程，對產品的設計及優化具有較大的幫助，也有利于縮短研發周期，降低經濟成本。
2模型建立
模型采用常用的橡膠材料與模具裝配模型，如圖1所示。整個建模過程與后續的有限元分析中均采用統一的mm單位制。
圖1 模型基本尺寸
3有限元分析
本案例的有限元分析是在ABAQUS 2017平臺上全程進行的。運用Standard/Explicit分析模塊，之后進入Part模塊創建上述分析模型。建立的有限元模型如圖2所示。模型中主要涉及兩種材料模型，橡膠本構已經很成熟了，選用超彈性Mooney-Rivlin本構，模具使用鋼鐵本構，輸入基本的物理參數即可。橡膠圈及鋼鐵本構參數分別如圖3、4所示。之后定義接觸及邊界條件完成有限元模型的前處理操作。
圖2有限元模型
圖3橡膠圈本構參數
圖4模具本構參數
4結果與討論
模型的后處理操作是在Abaqus/CAE的Visualization模塊，模型求解完成后對云圖只顯示材料填充區域云圖，此時，橡膠材料就從一開始的圓形被壓縮成類似于矩形的形狀，如圖5所示。
圖5應力云圖
5結論
本案例針對橡膠圈進行了一個簡單的大變形分析，從應力云圖來看，仿真結果很好模擬了橡膠圈在壓縮時候的大變形行為，后續可以單獨提取最大變形處的應力應變曲線等，對產品的設計有一定的參考意義。
展開 abaqus彈塑性粒狀材料的有限變形
Finite deformation of an elastic-plastic granular material.rar
ABAQUS網格劇烈變形導致計算終止
模型就是封面那樣，是一個螺紋連接二維對稱模型，我要算多孔金屬塑性，就是GTN模型，然后我在1這邊加了一個位移載荷，2這邊加一個U1的位移轉角，3這里是完全約束。
各相互接觸的螺紋牙之間都設置的面與面接觸，具體參數沒有設置，都是用的系統默認。
但是就算我載荷設置的再小，都會出現這種錯誤
The ratio of deformation speed to wave speed exceeds 1.0000 in at least one element. This usually indicates an error with the model definition. Additional diagnostic information may be found in the message file.
展開 
如何在ABAQUS對可視化后，生成的變形圖云圖，顯示等值線？
<p>1、點擊：視圖切面管理器</p><div contenteditable="false" width="100%">
<figure class="figure-image" data-img="https://img.jishulink.com/202312/attachment/b98ee6c3399c4d5bb4e7ff856e4e1ec9.png" style="text-align: center">
<img src="https://img.jishulink.com/202312/attachment/b98ee6c3399c4d5bb4e7ff856e4e1ec9.png" data-mobile-src="https://img.jishulink.com/202312/attachment/b98ee6c3399c4d5bb4e7ff856e4e1ec9.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/202312/attachment/b98ee6c3399c4d5bb4e7ff856e4e1ec9.png?image_process=/format,webp" data-initial-src="https://img.jishulink.com/202312/attachment/b98ee6c3399c4d5bb4e7ff856e4e1ec9.png">
</figure>
</div><p>2、點擊：create創建</p><div contenteditable="false" width="100%">
<figure class="figure-image" data-img="https://img.jishulink.com
展開 Abaqus模擬橡膠大變形
Abaqus為用戶提供了多種本構關系來模擬超彈性材料，這種材料具有高度非線性，當Abaqus進行模擬時假設這種材料是具有彈性、各向同性，并且同時考慮幾何非線性效應。與材料的剪切柔度相比，對于大多數類似橡膠的固體材料，其可壓縮性非常小，當分析對象為平面應力問題、殼、薄膜、梁、桁架、或者鋼筋等，這個問題不值得關注。但是對于固體、平面應變或者軸對稱問題卻不能忽略。對此，Abaqus/Standard提供了雜交單元來模擬超彈性材料中完全的不可壓縮行為。
橡膠材料力學性能的描述方法主要為兩類：一類是認為橡膠為連續介質的現象學描述；另一類是基于熱力學統計的方法。基于連續介質力學的本構模型主要有Polynomial、Reduce Polynomial、Ogden模型等，其中Mooney-Rivlin模型是 Polynomial的特殊形式，Neo-Hookean 模型是Reduce Polynomial的特殊形式。基于熱力學統計主要有Arruda-Boyce和Van der Waals等本構模型。本文利用Abaqus模擬大變形的橡膠，具體步驟如下。
1、在Abaqus/CAE Sketch模塊中作出模型草圖，如圖1所示，然后在Part模塊中分別建立Push、Rubber、Base三個部件。其中Push為解析剛體，Base為離散剛體。
圖1 草圖
2、在Property模塊中定義橡膠的屬性，采用Mooney-Rivlin模型，參數如圖2所示，然后賦給Rubber部件。
圖2 橡膠參數設置
3、裝配，定義分析步，采用默認的場輸出和歷史輸出。為了保證剛開始能夠較容易收斂，設置分析步初始增量步為0.01，打開幾何非線性。
圖3 分析步定義
4、定義接觸對：Push下表面和橡膠表面，Base上表面和橡膠表面。
展開 ABAQUS 單向拉伸大變形模擬
靜態模擬一種軟材料POE的單向拉伸，拉伸應變希望到300%，但是總是在100%就失敗了。不知道哪里出了問題，有沒有高手幫幫忙。
abaqus凍土路基的溫度-水分-變形多場耦合分析
凍土路基產生的變形與內地路基產生的變形不同。凍土路基隨著季節的交替發生凍結與融化的同時路面會產生相應的變形，并且這樣的變形隨著時間的推移還在持續不斷的變化。在同一路基橫斷面處，由于凍土路基溫度場和水分場分布的不同，路基表面會產生不均勻變形，即在道路橫向發生了變形。在青藏公路的不同路段，由于不同的路基填料、不同的路基高度、不同的多年凍土類型以及不同的路側積水等情況，會使得凍土路基形成縱向的波浪變形。
1 路基溫度場
溫度場的控制方程如下所示
由于凍土路基會存在凍結和融化過程，這就會伴隨著相變熱的產生，因此需要在傳統溫度控制方程中額外考慮相變熱的的影響。
路基的溫度場邊界比較復雜，本文采用第二類和第三類邊界條件，考慮太陽輻射、對流換熱和地面有效輻射的影響。太陽輻射主要影響大氣溫度變化，這里采用下式描述大氣溫度變化
對流換熱則采用下式描述
建立如圖所示的有限元模型
可以計算得到路基的溫度場分布和一年中路基的溫度變化如圖所示
2 水分場分析
凍土路基的變形與水的凍結和融化息息相關。所以分析凍土路基的變形時必須考慮水場分布的影響。
路基中水分場遷移可以通過達西定律來描述
由于凍土路基中，水分凍結后，水分會發生遷移，因此需要考慮相變對水分遷移的影響。
計算得到的飽和度分布如圖所示
3 變形場分析
凍土路基的變形包括融沉變形和車載變形。進行變形場分析時，采用摩爾庫倫準則
路面的車輛載荷采用脈沖載荷來模擬，如下圖所示
同時，水分的凍結時會產生凍脹變形，因此需要考慮凍脹率的影響。這里凍脹率選擇為0.03。
結合溫度場分析和水分場分析可以獲得路基的變形結果。
展開 