abaqus修正單元的案例
SACS軟件單元計算長度修正指南：不再困惑Ly與Lz
計算長度取值一般為不同水平撐桿之間的間距，不再是桿件的幾何長度
情況3：叉撐
導管架叉撐的情況比較復雜，很難準確判斷支撐邊界是什么，一般情況下認為叉撐對平面內是有效支撐，叉撐對平面外是無效支撐，
平面內（Ly）：取叉撐交點到主腿的距離，若桿件沒被其他桿件打斷，則不需修改
平面外（Lz）：取主腿之間的距離
情況4：平臺梁單元
對于上部平臺的梁單元：一般不需要修改計算長度
因為平臺梁通常采用工字形截面，平面內是強軸，平面外是弱軸；穩定問題通常由平面外控制，而平面外計算長度就是幾何長度
在SACS中修改計算長度很簡單：
選中需要修改的桿件——>打開桿件屬性對話框——>在對應位置輸入修正后的值
特別注意：一定要確保你修改的是正確的方向！弄反了可能導致嚴重后果。
總 結
計算長度反映的是桿件真實的支撐情況，而不是簡單的幾何長度。當你對一根桿件的支撐條件有清晰認識時，自然就知道該如何修正了。軟件是工具，您的判斷才是設計的核心。
希望這篇文章能幫你理清思路，對Ly和Lz不再迷茫。
展開 從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修正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.
最后，有需要歡迎通過微信公眾號聯系我們。
微信公眾號：320科技工作室。
展開 
Abaqus 修正劍橋模型的vumat子程序開發
傳統相關聯修正劍橋本構模型的屈服面方程、硬化準則和流動法則分別為
由于式(1)所確定的屈服軌跡在p平面是一個 圓，不能反映巖土介質拉壓不等(S-D)效應，而且劍 橋模型是基于正常固結狀態試驗推導而來，子午面上臨界狀態線通過應力坐標原點，表現為不考慮土 體黏聚力的純摩擦型本構，而大多數的巖土介質具 有一定的黏聚力，屬摩擦–黏聚型材料。考慮到莫 爾–庫侖準則有此特征，故將莫爾–庫侖準則與劍橋模型相結合，得
聯合式(1)，(4)，(5)，有改進的屈服函數：
根據上述理論，可以通過以下流程圖完成VUMAT子程序編寫
通過單胞模型計算得到得結果如下圖所示。
[1]袁克闊,陳衛忠,于洪丹,譚賢君,趙武勝,李香玲.考慮黏聚特性和拉壓不等效應的修正劍橋模型及數值實現[J].巖石力學與工程學報,2012,31(08):1574-1579.
展開 abaqus模擬隧道開挖，土體采用修正劍橋本構 ￥50
使用abaqus模擬隧道分階段開挖，土體開挖一層襯砌激活一層，土體本構采用修正劍橋本構，模型部分設置和結果見圖片
下加載面修正劍橋模型及其在Abaqus中umat子程序的實現
模型驗證
下面給出一個簡單的應用 UMAT 子程序的算例，將模型劃分為多個單元， 分別采用位移的加載方式和力的加載方式，并分別與試驗數據和 ABAQUS 自帶 的修正劍橋模型計算結果進行對比，以此說明 UMAT 子程序可以應用于多個單 元的模型，對于位移和力的加載方式均適用，并證明 UMAT 子程序適用于超固 結土體的力學特性分析和地基的固結沉降分析。
從而證明了編寫的子程序的有效性和實用性。
最后，大家有關于Abaqus二次開發的任何需求可以添加管理員微信號：CAE320，同時也歡迎大家關注“320科技工作室”的微信公眾號，掃一掃二維碼即可關注~~
展開 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增材制造AlSi10Mg修正JC本構二次開發（論文復現）
[圖片]
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 fe-safe疲勞分析CH03-損傷累積與修正因子
課程目錄
CH01-疲勞理論概述
CH02-fe-safe軟體介面
WORKSHOP01-帶孔洞圓管疲勞分析
CH03-損傷累積與修正因子
WORKSHOP02-轉軸多重形式荷載之疲勞分析
CH04-EN法介紹
WORKSHOP03-托架中低週期疲勞分析
視頻鏈接：https://www.yqgqt.org.cn/college/video/c15144
abaqus里的非線性薄層單元，零厚度cohesive單元，goodman接觸單元等的基本形式是什么？如何構建與應用？
在使用Abaqus，Comsol等軟件進行薄層區域的力學分析過程中，例如在研究水壓致裂、裂縫擴展，接觸粘結滑移的這類薄層力學性質時，我們經常需要采用應力-相對位移(σ-u)關系，而不是傳統本構描述的應力-應變(σ-ε)關系來描述，例如Abaqus里面的Cohesive單元，Goodman單元，以及Comsol里的彈性薄層（在后面我把這類單元統稱為增量非線性力學薄層）。這類單元厚度非常小甚至為0，薄層兩側的節點（單元）用一組力（應力）與相對位移的關系方程聯系起來，例如給出一個形式最為簡單的典型應力-位移方程
此方程描述了1,2,3方向（通常是法向和兩個切向）上相對位移與應力的關系，應力與相對位移呈線性關系，類似于“線性彈簧”。但是對于土-結構接觸、裂縫的張開閉合這類問題，線性方程已經不足以準確描述這些物理量之間的關系，這時就需要引入增量非線性方程來構建薄層單元。
引入增量非線性薄層的概念之前，首先介紹一下全量非線性薄層以理解非線性的概念，首先給出以下公式
這是一個全量非線性薄層，其非線性的表現可以用下面幾個例子體現，
對比①和②項，可以發現僅存在3方向上的位移變化的情況下，1，2方向上的力也會發生改變，體現了彈簧三個方向力學性質的非獨立性，對比①和③項，可以發現力的大小并不和位移大小成正比，也就是非線性特征。
所以對于增量非線性方程，就是把應力-位移關系方程寫成應力增量-位移增量的關系方程，例如
寫成微分形式的好處是，可以體現出應力路徑對位移結果的影響，也就是類似于“塑性”特征（所以所有的彈塑性本構也都是增量方程）。但是對于此類微分方程的求解，必須給定一個力的初始值。
展開 BCC點陣結構梁單元Abaqus壓縮仿真模擬-顯示動力學質量縮放 ￥19.89
本文通過abaqus顯示動力學的方法對BCC結構進行壓縮仿真模擬，同時為減小計算量，采用梁單元模擬點陣結構，壓頭設置為剛性面，添加質量縮放，加快運算速度，為點陣結構壓縮模擬提供一種便捷方法。
1. 建立BCC點陣模型,以單胞尺寸5X5X5為例。
a.首先建立立方體實體，然后對實體進行處理，得到點陣單胞點陣結構。
b.建立單胞BCC梁單元點陣模型,然后進行刪除面的操作，得到單胞BCC點陣結構，接下來進行陣列操作，得到最大外形尺寸為25x25x25的點陣壓縮模擬試件。
C.建立剛性壓板，設置參考點，模擬萬能試驗機壓頭，剛性單元不參與計算，不影響計算結果，加快運算速度。
2. 裝配，按壓縮試驗進行裝配，從上到下依次為壓板-點陣-壓板。
3.設置材料屬性，本文為鈦合金TC4，密度4.43e-9Tone/mm3,彈性模量為118000MPa，泊松比0.3，應力應變值見下表所示。
設置截面屬性Beam，定義截面半徑0.5mm
指派截面，定義方向[1,2,3]，完成材料屬性設置。
4.設置分析步Dynamic,Explicit，時間設置為5s，以每秒1mm的速度進行壓縮模擬，開啟質量縮放為1e-5，歷程輸出勾選位移和力，以便輸出力-位移曲線，然后計算相應的應力-應變曲線。
5.設置相互作用-切向行為和法向行為，摩擦系數為0.3，設置通用接觸。
以下部分為付費部分
展開 【JY】Abaqus“殼”單元概述與應用（二）——固體殼單元
寫在前文
在有限元分析中，單元類型的選擇對計算結果的精度和效率有著決定性影響，尤其對于復合材料結構和薄壁結構的分析更是如此。
Abaqus 作為主流的有限元分析軟件，提供了多種固體殼單元類型以滿足不同工程需求。連續實體殼單元 (CSS8)、非協調元 (C3D8I) 和連續殼單元 (SC8R) 是 Abaqus 中常用于復合材料和薄壁結構分析的三種單元類型，各自具有獨特的理論基礎和適用場景。
相關閱讀：
【JY】Abaqus殼單元概述與應用（一）
除了上述采用類實體單元的“殼”單元外，還有完全的殼單元，如S4R 單元，是 Abaqus 中最常用的常規殼單元之一，為 4 節點減縮積分殼單元，基于經典殼理論，適用于各類薄壁結構的線性與非線性分析，尤其在大變形和接觸問題中表現穩定，將該單元作為對比基準，對上述實體類“殼”單元進行對比分析。
本文旨在對這三種單元類型進行深入比較研究，從理論基礎、自由度、材料本構、積分方案、閉鎖敏感性、計算成本等多個維度展開分析，為工程實踐中的單元選擇提供參考。特別是針對復合材料分析、金屬薄壁結構模擬以及混合建模等應用場景，探討這三種單元的適用性差異，并分析它們在幾何非線性情況下的計算成本和精度表現。
單元類型基本原理與特點
2.1 連續實體殼單元 (CSS8)
連續實體殼單元 (CSS8) 是一種介于 C3D8I (非協調元) 和 SC8R (連續殼單元) 之間的特殊一階單元，由 Vu-Quoc 和 Tan 于 2003 年提出，后集成于 SIMULIA 2017 及以后的版本。它是一種三維單元，具有以下基本特點：
幾何與自由度：CSS8 為 8 節點六面體單元，僅有位移自由度 (無轉動自由度，與實體單元一致)，與實體單元混合建模時易于處理連接過渡。
展開 