abaqus剛度修正的案例
引入位移邊界條件和整體剛度矩陣的修正
消除了整體剛度矩陣的奇異性后，才能從方程組（32）求解結點位移。在一般情況下，所考慮問題的邊界往往已有一定的約束條件，排除了剛體運動的可能性。否則，可適當指定某些結點的位移值，以避免計算機存儲作大的更動。下面就介紹兩種比較簡單的引入已知結點位移的方法。
1、對角元素改l法
這種方法是把結點的指定值置入方程給(32)，保持方程仍是2nX2n階，而將K和P修正。例如，若指定結點i在y方向位移vi的值，則令K中的元素Ki,i為1，而第i行和i列的其余元素都為零。P中的第i個元素則用位移v的已知值代入，P中的其他各行元素都減去結點位移的指定值和原來K中這行的相應行元素的乘積。
為了說明這一引進結點已知位移的過程，我們來考察下面只有四個方程的簡單例子。方程(32)展開成如下的形式
設這個系統中結點位移u1和u2被指定為
當引用上述方法后，方程(50)就變成
然后，就用這組維數不變的方程來求解所有的結點位移。顯然，其解答為u1=β1、u2=β3；
v1、v2仍為原方程的解答。
這種方法最適用于給定零位移，此時除將給定的零值位移修改對應的載荷陣元(如例中令Px1=0，Px2=0)外，其他載荷陣中的元素不必作任何修正。
2、對角元素乘大數法
此法是將K中與指定結點位移有關的主對角元素乘上一個大數，例2xl015，同時將p的
對應元素換上結點位移指定值與同一個大數的乘積。實際上，這種方法就是使得K中相應行的修正項遠大于非修正項。
展開 從C3D8的uel源代碼入門Abaqus的uel編寫, 更新B-Bar修正 ￥99
/202403/b4a813b2b2c5b8764d7fbbf45e643b3b.png"><img src="https://img.jishulink.com/msimage/202403/18d7ba93ba85e4f55fe830ab64d17f31.png"></h1><h1>uel源代碼</h1><p>2024/03/27更新:已加入B-bar算法進行剛度矩陣的修正,可看到uel的計算結果已于Abaqus計算結果一致</p>
展開 ABAQUS用戶定義單元UEL與VUEL從入門到放棄系列2 ￥50
本構關系的張量表達式為：
其中，廣義應變的有限元格式可表示為：
其中，電場強度E為負的電勢梯度：
則廣義應變列陣記為：
那么，單元的剛度矩陣可以表示為：
其中Kuu為C3D8原本的剛度矩陣，Kuv與Kvu為壓電耦合剛度矩陣，Kvv為電場的廣義剛度矩陣。
ABAQUS中的C3D8和C3D8E都是做了一些剛度修正的，比如C3D8為了防止單元自鎖，采用了B-Bar方法，得出的剛度矩陣是介于C3D8和C3D8R之間的值，同樣的C3D8E也有一些類似的修正，以下我將提供一個不包含修正的版本，對ABAQUS剛度修正方法感興趣的朋友可以去拿去跟ABAQUS CAE對比。
展開 Abaqus修正GTN模型的VUMAT子程序
Gurson-Tvergaard-Needleman ( GTN) 模型是研究金屬損傷的重要工具。GTN 模型通過孔洞體積分數的演變來判斷材料的失效， 但不適用于剪切斷裂為主的韌性斷裂。本文在GTN模型中引入剪應力的影響，編寫了相關的VUMAT子程序。
GTN模型的屈服函數可以用下式表示
其中q1，q2是模型參數，取q1=1.5，q2=1,σ0為等效應力，p為靜水應力，q為Mises等效應力；f為空洞的體積分數。
p和q可以通過徑向返回算法得到
應變控制的孔洞形核系數
GTN模型可以通過以下4個方程進行描述
Nahshon and Hutchinson考慮了剪應力對模型的影響
于是孔隙體積分數的演化可以通過下式描述
仿真計算得到的結果如下圖所示
有問題私信或者關注cae320公眾號
展開 
Abaqus 修正劍橋模型的vumat子程序開發
傳統相關聯修正劍橋本構模型的屈服面方程、硬化準則和流動法則分別為
由于式(1)所確定的屈服軌跡在p平面是一個 圓，不能反映巖土介質拉壓不等(S-D)效應，而且劍 橋模型是基于正常固結狀態試驗推導而來，子午面上臨界狀態線通過應力坐標原點，表現為不考慮土 體黏聚力的純摩擦型本構，而大多數的巖土介質具 有一定的黏聚力，屬摩擦–黏聚型材料?？紤]到莫 爾–庫侖準則有此特征，故將莫爾–庫侖準則與劍橋模型相結合，得
聯合式(1)，(4)，(5)，有改進的屈服函數：
根據上述理論，可以通過以下流程圖完成VUMAT子程序編寫
通過單胞模型計算得到得結果如下圖所示。
[1]袁克闊,陳衛忠,于洪丹,譚賢君,趙武勝,李香玲.考慮黏聚特性和拉壓不等效應的修正劍橋模型及數值實現[J].巖石力學與工程學報,2012,31(08):1574-1579.
展開 基于Abaqus的修正劍橋模型的vumat子程序開發
傳統相關聯修正劍橋本構模型的屈服面方程、硬化準則和流動法則分別為
由于式(1)所確定的屈服軌跡在p平面是一個 圓，不能反映巖土介質拉壓不等(S-D)效應，而且劍 橋模型是基于正常固結狀態試驗推導而來，子午面上臨界狀態線通過應力坐標原點，表現為不考慮土 體黏聚力的純摩擦型本構，而大多數的巖土介質具 有一定的黏聚力，屬摩擦–黏聚型材料?？紤]到莫 爾–庫侖準則有此特征，故將莫爾–庫侖準則與劍橋模型相結合，得
聯合式(1)，(4)，(5)，有改進的屈服函數：
根據上述理論，可以通過以下流程圖完成VUMAT子程序編寫
通過單胞模型計算得到得結果如下圖所示。
[1]袁克闊,陳衛忠,于洪丹,譚賢君,趙武勝,李香玲.考慮黏聚特性和拉壓不等效應的修正劍橋模型及數值實現[J].巖石力學與工程學報,2012,31(08):1574-1579.
最后，有需要歡迎通過微信公眾號聯系我們。
微信公眾號：320科技工作室。
展開 abaqus模擬隧道開挖，土體采用修正劍橋本構 ￥50
使用abaqus模擬隧道分階段開挖，土體開挖一層襯砌激活一層，土體本構采用修正劍橋本構，模型部分設置和結果見圖片
abaqus剪切修正GTN模型的VUMAT子程序開發
這里對相應的算法進行簡要說明：
NH-GTN模型
屈服函數：
其中等效孔洞體積分數定義為：
孔洞體積分數包含新孔隙形核，原有空隙生長以及剪切相關的等效體積分數增加：
形核，生長，剪切相關體積分數的演化遵循：
其中：
剪切效應的修正，考慮應力狀態的影響
參數的物理含義如下
通過將文獻中的數值算法編程實現在VUMAT子程序中，可以用來實現對延性金屬材料在不同應力狀態下的損傷演化進行合理的數值預測，應用于金屬成型領域（沖壓，軋制，擠壓等）
預測修正后的模型應該在簡單拉伸情況下于abaqus自帶的GTN模型保持相同的損傷和其他狀態變量的分布，并在剪切情況中損傷發展顯著高于abaqus自帶的模型（自帶的模型忽略了剪切效應）。（為了進行對比使用于自帶的本構相同的硬化方式，模擬中使用了相同的質量縮放，但質量縮放容易產生數值振蕩，模擬的拉伸曲線存在波動。）
初步模擬結果：
拉伸情況（abaqus-VUMAT）
應力分情況
孔洞體積分數
剪切模型（abaqus-VUMAT）
不同變形時刻的應力分布
T=0.1s
局部放大圖
T=0.5s
局部放大圖
T=0.6s
局部放大圖
可以看到模型在拉伸預測中與原始模型保持一致，而在剪切修正后損傷發展顯著快于原始模型，利用作者提出的方法可以應用于復雜應力狀態下金屬材料的損傷分析，相關參數部分參考文獻，其中Kw=3.T1=0.2，T2=0.7.模擬結果符合文獻所提出方法的基本趨勢。
最后，如果有相關需要歡迎通過公眾號“320科技工作室”聯系我們。
展開 ABAQUS路面材料使用修正burgers模型時總是出現編譯錯誤
C 瞬態溫度場下修正Burgers模型UMAT子程序源代碼
C
C 給狀態變量數組賦初值為零，調用ABAQUS子程序SDVINI
C GIVE STATEV THE INITIAL VALUE OF ZERO
C
SUBROUTINE SDVINI(STATEV,COORDS,NSTATV,NCRDS,NOEL,NPT,LAYER,KSPT)
C
INCLUDE 'ABA_PARAM.INC'
C
DIMENSION STATEV(NSTATV),COORDS(NCRDS)
C
DO K=1,NSTATV
STATEV(K)=0.0
END DO
C
RETURN
END
C 瞬態溫度場下修正Burgers模型UMAT子程序
C UMAT FOR MODIFIED BURGERS MODEL
C
SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT,
1 DRPLDE,DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,
2 CMNAME,NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,
3 PNEWDT,CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)
C
INCLUDE 'ABA_PARAM.INC'
C
CHARACTER*80 CMNAME
C
DIMENSION STRESS(NTENS),STATEV(NSTATV),DDSDDE
展開 ABAQUS 小應變分析(例4) 修正劍橋模型(MCC)模擬粘土的比奧固結 ￥67
ABAQUS 小應變分析(例4) 修正劍橋模型(MCC)模擬粘土的比奧固結
修正劍橋模型(MCC)被廣泛應用于粘土的滲流固結當中，能較準確地預測因滲流固結導致的土體沉降，有效應力變化，孔隙水壓力和孔隙比(e)的變化。該模型模擬粘土(Clay)在受荷作用下土體的固結，粘土為均質粘土，其先期固結壓力為200kPa，在實施地應力平衡后，土體頂部施加50kPa的固結壓力。土表面為自由滲流邊界。
建模及結果展示：
模型位移邊界及頂部50kPa荷載
模型網格劃分
土表面在50kPa荷載下隨時間的沉降位移曲線
土的孔隙水壓力分布
土的孔隙比(e)的分布
下加載面修正劍橋模型及其在Abaqus中umat子程序的實現
模型驗證
下面給出一個簡單的應用 UMAT 子程序的算例，將模型劃分為多個單元， 分別采用位移的加載方式和力的加載方式，并分別與試驗數據和 ABAQUS 自帶 的修正劍橋模型計算結果進行對比，以此說明 UMAT 子程序可以應用于多個單 元的模型，對于位移和力的加載方式均適用，并證明 UMAT 子程序適用于超固 結土體的力學特性分析和地基的固結沉降分析。
從而證明了編寫的子程序的有效性和實用性。
最后，大家有關于Abaqus二次開發的任何需求可以添加管理員微信號：CAE320，同時也歡迎大家關注“320科技工作室”的微信公眾號，掃一掃二維碼即可關注~~
展開 
abaqus fe-safe疲勞分析CH03-損傷累積與修正因子
課程目錄
CH01-疲勞理論概述
CH02-fe-safe軟體介面
WORKSHOP01-帶孔洞圓管疲勞分析
CH03-損傷累積與修正因子
WORKSHOP02-轉軸多重形式荷載之疲勞分析
CH04-EN法介紹
WORKSHOP03-托架中低週期疲勞分析
視頻鏈接：https://www.yqgqt.org.cn/college/video/c15144
ABAQUS增材制造AlSi10Mg修正JC本構二次開發（論文復現）
[圖片]
【ABAQUS模態動力學】Composite&abaqus 預應力模態分析&輸出單元剛度矩陣
劃分網格
定義邊界條件
5.2 預應力模態
預應力模態，按我的理解就是，假設t=0時，結構的剛度矩陣、質量矩陣為M0,K0；t=t1時，結構（分析對象）收到外部激勵的作用，使得結構的剛度矩陣，質量矩陣發生改變，結構的剛度矩陣、質量矩陣為M1,K1；t=t2時，開始進行特征值提取，此時求解的是t1狀態的結果。
從上面這個理解出發，ABAQUS預應力模態只要在frequency分析步之前進行General,Static分析步,打開NLGeom選項(分析過程中剛度矩陣會不斷變化)。
提取單元剛度矩陣：
【ABAQUS 二次開發筆記】輸出單元剛度矩陣 - hayden_william - 博客園
以上均為我的一點理解，不一定完全正確，本文僅作為個人學習記錄之用，其他概不負責。
展開 ABAQUS 剛度矩陣 ￥1000
我有個abaqus的問題，你們幫幫我出出主意。 是這樣的：
1，TestElement.py 是編寫的測試代碼，可輸出8結點線性單元的剛度矩陣。
2，abaqus文件Job-testing.inp ，運行可以輸出單元剛度矩陣。
問題是：他們的結點坐標/排序，材料參數都一致， 但得到的剛度矩陣就不一樣。
需要：代碼輸出的剛度矩陣與abaqus得到的剛度矩陣一致或者基本一致。
