usdfld子程序開發的案例
ABAQUS用戶子程序USDFLD應用實例
一、引言
為滿足用戶多樣化的分析需求，ABAQUS提供了眾多子程序接口。其中USDFLD可用于復雜材料模型的建模分析，例如粉末冶金、損傷分析等。本文以某藥物粉末的模壓為例，介紹了ABAQUS用戶子程序USDFLD的使用方法。
二、用戶子程序USDFLD簡介
1.USDFLD即Use Defined Field，通過將材料屬性定義為場變量fi的函數，來實現復雜材料行為的建模
①ABAQUS/Standard中的大多數材料屬性都可定義為場變量的函數。
②USDFLD允許用戶在單元的每個積分點處定義fi。
展開 Abaqus子程序USDFLD使用案例講解
ABAQUS中使用USDFLD子程序的應用例子 ￥10
今天來解讀一個在ABAQUS中使用USDFLD子程序實現損傷計算的例子。
USDFLD子程序的界面如下：
需要定義的變量
數組FIELD（NFIELD）包含當前材料（積分）點的場變量fi。
——在當前增量步的末尾，這些傳入了由節點插值得到的值，由*INITIAL CONDITIONS選項或*FIELD選項定義。
——更新的fi用于計算材料屬性值，其作為場變量的函數。更新的fi傳遞到其他用戶子程序中（CREEP、HETVAL、UEXPAN，UHARD、UHYPEL、UMAT、UMATHT和UTRS），在該材料點可以被調用。
可能需要定義的變量
• 包含依賴于解的狀態變量的數組STATEV（NSTATV），能夠在USDFLD中定義。
——這些在增量步的開始作為值傳入。
——在該子程序中STATEV在所有情況下都能夠更新，所有更新的值被傳遞到其他用戶子程序中（CREEP、HETVAL、UEXPAN，UHARD、UHYPEL、UMAT、UMATHT和UTRS），在該材料點可以被調用。
——與該材料點相關的狀態變量的數量由*DEPVAR選項定義。
• 建議的新時間步長與現在正在使用的時間步長（DTIME）的比例PNEWDT，可以給定。
——該變量允許用戶在ABAQUS中的自動時間增量步長算法中提供輸入。
展開 Abaqus通過USDFLD子程序進行泥巖的應變軟化模擬
本文基于泥巖的三軸壓縮試驗曲線，建立考慮應變軟化特性的泥巖彈塑性本構模型，使用Abaqus及其子程序對泥巖的三軸壓縮試驗進行了數值模擬。泥巖在受壓過程中主要經過了5個階段，即壓密階段、彈性變 形、應變硬化、應變軟化、殘余階段。泥巖應變軟化模型如下所示。
式中，ξ為強度參數，ξp為峰值強度參數，ξr為殘余階段強度參數，η為應變軟化參數，η*為殘余階段的應變軟化參數初始值。對于三軸壓縮試驗，η用塑性剪切應變來表示
塑性屈服準則采用Mohr-Coulomb準則，則粘聚力和內摩擦角的參數演化可以用下式表示
在巖石的塑性變形過程中會產生比較明顯的剪脹現象，而用來描述這一現象的較常用的力學參數就是剪脹角 Ψ，
上述模型可以通過USDFLD子程序進行實現，流程圖如下
有限元模型如下圖所示
計算得到的應力云圖及不同圍壓下的載荷位移響應如下所示
參考文獻：張力偉,賈善坡,鄒江濤,舒婧曦.泥巖的峰后軟化力學模型.中國科技論文,2016,11(21):2456-2461
有關于abaqus子程序開發的相關問題可以通過公眾號聯系我們.
公眾號: 320科技工作室
展開 
USDFLD子程序復合材料強度分析 ￥20
USDFLD子程序復合材料強度分析
abaqus基于usdfld子程序的內聚力疲勞模型
基于usdfld的內聚力疲勞模型
之前在Abaqus Cohesive單元的疲勞UMAT － 技術鄰 (jishulink.com) 上介紹了通過umat子程序來編寫內聚力疲勞本構，實現裂紋疲勞擴展的方法。在實際計算中經常出現不收斂的情況，因此重新編寫了雙線性本構下的usdfld內聚力疲勞子程序。
在不考慮疲勞損傷的情況下，單一裂紋模式雙線性內聚力本構如圖所示
混合模式下的斷裂準則采用BK準則
損傷萌生和失效對應的等效張開位移分別為
和
疲勞損傷采用roe提出的損傷演化方程
考慮疲勞損傷后的內聚力本構如圖所示
這里同樣假設卸載以及法向壓縮不會累積疲勞損傷。
建立三點彎曲模型對疲勞裂紋擴展進行了模擬，計算結果如圖所示
斷裂過程
跨中載荷位移曲線
損傷演化過程
跨中底部單元的應力應變關系
更新預告：早期混凝土熱-濕-力多場耦合分析，編寫了基于水化度理論和考慮熱學參數變化的溫度場計算子程序（umatht和film）; 考慮溫度對濕度擴散系數影響的濕度場計算子程序（umatht和film）; 基于成熟度理論和雙冪徐變函數的應力場子程序（umat）。
展開 Abaqus用戶子程序USDFLD調用GETVRM 返回值為0
.0)THEN WRITE(6,*) 'REQUEST ERROR IN USDFLD FOR ELEMENT NUMBER ', 1 NOEL,'INTEGRATION POINT NUMBER ',NPT ENDIF C RETURN END
[原創]Abaqus UFIELD 和 USDFLD子程序詳解
UFIELD和USDFLD子程序詳解
———公眾號‘CAE仿真實驗室’出品
UFIELD和USDFLD是Abaqus極具迷惑性的兩個子程序，這兩個家伙長得就比較像兄弟，而且都是用來自定義場變量的，同時還有一個狀態變量state variable摻和在中間，再加上幫助里面的解釋也基本上不是人話，很容易就把人搞二了，今天帖主就來理一理這兩個子程序。
1、自定義場變量
Abaqus中場變量可以先簡單的理解為每個單元或者每個節點都有一個值的變量，像結果輸出中的溫度、應力和應變等。有時候我們也需要定義Abaqus中不存在的一些稀奇古怪，有意義或無意義的場變量，比如濃度場，固化度場，損傷度。場變量最場用于定義變化的材料參數，即讓材料屬性如密度，彈性模量和一個場變量相關，通過更新場變量的值來改變材料參數（這一點和溫度相關的材料參數本質上是一致的），而場變量的更新在這兩個子程序中均可以完成，UFIELD是用來指定預定義場變量的，USDFLD對積分點的場變量重新定義。
展開 Abaqus通過USDFLD子程序進行泥巖的應變軟化模擬
本文基于泥巖的三軸壓縮試驗曲線，建立考慮應變軟化特性的泥巖彈塑性本構模型，使用Abaqus及其子程序對泥巖的三軸壓縮試驗進行了數值模擬。泥巖在受壓過程中主要經過了5個階段，即壓密階段、彈性變 形、應變硬化、應變軟化、殘余階段。泥巖應變軟化模型如下所示。
式中，ξ為強度參數，ξp為峰值強度參數，ξr為殘余階段強度參數，η為應變軟化參數，η*為殘余階段的應變軟化參數初始值。對于三軸壓縮試驗，η用塑性剪切應變來表示
塑性屈服準則采用Mohr-Coulomb準則，則粘聚力和內摩擦角的參數演化可以用下式表示
在巖石的塑性變形過程中會產生比較明顯的剪脹現象，而用來描述這一現象的較常用的力學參數就是剪脹角 Ψ，
上述模型可以通過USDFLD子程序進行實現，流程圖如下
有限元模型如下圖所示
計算得到的應力云圖及不同圍壓下的載荷位移響應如下所示
參考文獻：張力偉,賈善坡,鄒江濤,舒婧曦.泥巖的峰后軟化力學模型.中國科技論文,2016,11(21):2456-2461
有關于abaqus子程序開發的相關問題可以聯系扣扣1653004885或者關注cae320公眾號
展開 ABAQUS案例—帶USDFLD子程序的復合材料強度分析 ￥3
本案例（附件中的inp文件及for文件）講述了采用ABAQUS進行含USDFLD子程序的復合材料的強度分析。
在讀入inp后在JOB模塊輸入for文件即可運行。
Abaqus聯用USDFLD+HETVAL子程序分析慢速烤燃 ￥99
采用HETVAL子程序的一步分解慢速烤燃模型請參考：https://www.yqgqt.org.cn/content/post/1300451
然而有時反應并不是一步完成的，需要多步反應模型，此時HETVAL并不夠用。為此，調用強大的USDFLD子程序跟HETVAL聯用解決問題。當然，USDFLD非常強大，本例僅演示其與HETVAL的聯用并分析慢速烤燃過程。
本例提供HETVAL、USDFLD子程序、cae及inp源文件，咨詢請聯系我 QQ：180280578。
附件包括詳細建模過程、材料參數、注意事項及子程序。
展開 
abaqus利用usdfld子程序在焊接降溫階段改變材料參數 ￥29.9
abaqus利用usdfld子程序在焊接降溫階段改變材料參數，在tig焊接仿真中，焊接后，焊縫組織強度變高，為了提高仿真的精度，在焊接時，在升溫結束開始降溫時提高焊縫的材料性能。子程序包括了采用高斯熱源dflux和usdfld聯合仿真。（該usdfld也可用在混凝土在升溫降溫不同的材料參數）
USDFLD子程序實現循環荷載作用下土體的剛度弱化
哪位大佬用USDFLD子程序編寫過循環荷載土體剛度弱化
[原創]Abaqus傳熱分析HETVAL和USDFLD子程序聯合
在Thermal 選項下添加Heat Generation
5、邊界和分析步
給定初始溫度30 、通過hetval子程序添加熱生成率 1e7W/m3，下端面全部約束。選擇*Coupled Temperature-displacement 分析步，進行熱力學分析
6、建立作業
此時已經定義了USDFlD和HETVAL子程序，但是USDFLD子程序還沒有激活，需要通過* *INITIAL CONDITIONS, TYPE=FIELD來激活，導出iob inp文件 ，在STEP前面加上如下命令*INITIAL CONDITIONS, TYPE=FIELD, VARIABLE=1，如截圖所示
最后通過cmd命令提交job
7、結果展示
初始溫度30度
通過Hetval 加載熱生成率1e7W/m3，計算完溫度32.84
熱應力
理論計算
根據比熱容 即可得到
7800*0.5*452*=1e7*0.5*1
可見和計算結果完全吻合，說明了熱源加載的正確性
8、總結
1、兩個子程序寫在一個for文件中，保證程序的正確性
2、兩個子程序都要在material模塊設置，并且USDFLD要通過initial condition關鍵字激活，否則不起作用
3、Initial condition 不支持CAE操作 所以通過cmd 提交job
4、USDFLD中按實際需求更改STATE(1)的值
歡迎小伙伴們掃碼關注我的公眾號‘CAE仿真實驗室’哦，我們將定期給大家分享有限元仿真干貨,謝謝~
展開 Abaqus用戶子程序USDFLD實例詳解-復合材料層壓板漸進損傷強度
Abaqus用戶子程序USDFLD實例詳解--復合材料層壓板漸進損傷強度分析
1、用USDFLD子程序漸進損傷分析方法分析層壓板的強度。
2、層壓板由25層單層板組成，單層板厚度為0.15mm。層壓板的鋪層順序、單層板的材料屬性見下表。
3、層壓板幾何尺寸：
厚度 T = 3.75mm;
寬度 W = 12mm；
長度 L >1.5W。
L可以任取，主要研究材料性能，為研究大部件性能做參考。
4、加載方式：沿長度方向施加壓縮載荷，在本例中施加位移載荷，讀取相應的力。
詳細步驟：
code.zip
Abaqus用戶子程序USDFLD實例詳解.rar
展開 