ABAQUS usdfld子程序的案例
abaqus利用usdfld子程序在焊接降溫階段改變材料參數(shù) ￥29.9
abaqus利用usdfld子程序在焊接降溫階段改變材料參數(shù)，在tig焊接仿真中，焊接后，焊縫組織強(qiáng)度變高，為了提高仿真的精度，在焊接時(shí)，在升溫結(jié)束開始降溫時(shí)提高焊縫的材料性能。子程序包括了采用高斯熱源dflux和usdfld聯(lián)合仿真。（該usdfld也可用在混凝土在升溫降溫不同的材料參數(shù)）
ABAQUS案例—帶USDFLD子程序的復(fù)合材料強(qiáng)度分析 ￥3
本案例（附件中的inp文件及for文件）講述了采用ABAQUS進(jìn)行含USDFLD子程序的復(fù)合材料的強(qiáng)度分析。
在讀入inp后在JOB模塊輸入for文件即可運(yùn)行。
ABAQUS中使用USDFLD子程序的應(yīng)用例子 ￥10
今天來(lái)解讀一個(gè)在ABAQUS中使用USDFLD子程序實(shí)現(xiàn)損傷計(jì)算的例子。
USDFLD子程序的界面如下：
需要定義的變量
數(shù)組FIELD（NFIELD）包含當(dāng)前材料（積分）點(diǎn)的場(chǎng)變量fi。
——在當(dāng)前增量步的末尾，這些傳入了由節(jié)點(diǎn)插值得到的值，由*INITIAL CONDITIONS選項(xiàng)或*FIELD選項(xiàng)定義。
——更新的fi用于計(jì)算材料屬性值，其作為場(chǎng)變量的函數(shù)。更新的fi傳遞到其他用戶子程序中（CREEP、HETVAL、UEXPAN，UHARD、UHYPEL、UMAT、UMATHT和UTRS），在該材料點(diǎn)可以被調(diào)用。
可能需要定義的變量
• 包含依賴于解的狀態(tài)變量的數(shù)組STATEV（NSTATV），能夠在USDFLD中定義。
——這些在增量步的開始作為值傳入。
——在該子程序中STATEV在所有情況下都能夠更新，所有更新的值被傳遞到其他用戶子程序中（CREEP、HETVAL、UEXPAN，UHARD、UHYPEL、UMAT、UMATHT和UTRS），在該材料點(diǎn)可以被調(diào)用。
——與該材料點(diǎn)相關(guān)的狀態(tài)變量的數(shù)量由*DEPVAR選項(xiàng)定義。
• 建議的新時(shí)間步長(zhǎng)與現(xiàn)在正在使用的時(shí)間步長(zhǎng)（DTIME）的比例PNEWDT，可以給定。
——該變量允許用戶在ABAQUS中的自動(dòng)時(shí)間增量步長(zhǎng)算法中提供輸入。
展開 Abaqus用戶子程序USDFLD實(shí)例詳解-復(fù)合材料層壓板漸進(jìn)損傷強(qiáng)度
Abaqus用戶子程序USDFLD實(shí)例詳解--復(fù)合材料層壓板漸進(jìn)損傷強(qiáng)度分析
1、用USDFLD子程序漸進(jìn)損傷分析方法分析層壓板的強(qiáng)度。
2、層壓板由25層單層板組成，單層板厚度為0.15mm。層壓板的鋪層順序、單層板的材料屬性見下表。
3、層壓板幾何尺寸：
厚度 T = 3.75mm;
寬度 W = 12mm；
長(zhǎng)度 L >1.5W。
L可以任取，主要研究材料性能，為研究大部件性能做參考。
4、加載方式：沿長(zhǎng)度方向施加壓縮載荷，在本例中施加位移載荷，讀取相應(yīng)的力。
詳細(xì)步驟：
code.zip
Abaqus用戶子程序USDFLD實(shí)例詳解.rar
展開 
[原創(chuàng)]Abaqus UFIELD 和 USDFLD子程序詳解
UFIELD和USDFLD子程序詳解
———公眾號(hào)‘CAE仿真實(shí)驗(yàn)室’出品
UFIELD和USDFLD是Abaqus極具迷惑性的兩個(gè)子程序，這兩個(gè)家伙長(zhǎng)得就比較像兄弟，而且都是用來(lái)自定義場(chǎng)變量的，同時(shí)還有一個(gè)狀態(tài)變量state variable摻和在中間，再加上幫助里面的解釋也基本上不是人話，很容易就把人搞二了，今天帖主就來(lái)理一理這兩個(gè)子程序。
1、自定義場(chǎng)變量
Abaqus中場(chǎng)變量可以先簡(jiǎn)單的理解為每個(gè)單元或者每個(gè)節(jié)點(diǎn)都有一個(gè)值的變量，像結(jié)果輸出中的溫度、應(yīng)力和應(yīng)變等。有時(shí)候我們也需要定義Abaqus中不存在的一些稀奇古怪，有意義或無(wú)意義的場(chǎng)變量，比如濃度場(chǎng)，固化度場(chǎng)，損傷度。場(chǎng)變量最場(chǎng)用于定義變化的材料參數(shù)，即讓材料屬性如密度，彈性模量和一個(gè)場(chǎng)變量相關(guān)，通過(guò)更新場(chǎng)變量的值來(lái)改變材料參數(shù)（這一點(diǎn)和溫度相關(guān)的材料參數(shù)本質(zhì)上是一致的），而場(chǎng)變量的更新在這兩個(gè)子程序中均可以完成，UFIELD是用來(lái)指定預(yù)定義場(chǎng)變量的，USDFLD對(duì)積分點(diǎn)的場(chǎng)變量重新定義。
展開 ABAQUS用戶子程序USDFLD應(yīng)用實(shí)例
一、引言
為滿足用戶多樣化的分析需求，ABAQUS提供了眾多子程序接口。其中USDFLD可用于復(fù)雜材料模型的建模分析，例如粉末冶金、損傷分析等。本文以某藥物粉末的模壓為例，介紹了ABAQUS用戶子程序USDFLD的使用方法。
二、用戶子程序USDFLD簡(jiǎn)介
1.USDFLD即Use Defined Field，通過(guò)將材料屬性定義為場(chǎng)變量fi的函數(shù)，來(lái)實(shí)現(xiàn)復(fù)雜材料行為的建模
①ABAQUS/Standard中的大多數(shù)材料屬性都可定義為場(chǎng)變量的函數(shù)。
②USDFLD允許用戶在單元的每個(gè)積分點(diǎn)處定義fi。
展開 Abaqus子程序USDFLD使用案例講解
Abaqus通過(guò)USDFLD子程序進(jìn)行泥巖的應(yīng)變軟化模擬
本文基于泥巖的三軸壓縮試驗(yàn)曲線，建立考慮應(yīng)變軟化特性的泥巖彈塑性本構(gòu)模型，使用Abaqus及其子程序對(duì)泥巖的三軸壓縮試驗(yàn)進(jìn)行了數(shù)值模擬。泥巖在受壓過(guò)程中主要經(jīng)過(guò)了5個(gè)階段，即壓密階段、彈性變 形、應(yīng)變硬化、應(yīng)變軟化、殘余階段。泥巖應(yīng)變軟化模型如下所示。
式中，ξ為強(qiáng)度參數(shù)，ξp為峰值強(qiáng)度參數(shù)，ξr為殘余階段強(qiáng)度參數(shù)，η為應(yīng)變軟化參數(shù)，η*為殘余階段的應(yīng)變軟化參數(shù)初始值。對(duì)于三軸壓縮試驗(yàn)，η用塑性剪切應(yīng)變來(lái)表示
塑性屈服準(zhǔn)則采用Mohr-Coulomb準(zhǔn)則，則粘聚力和內(nèi)摩擦角的參數(shù)演化可以用下式表示
在巖石的塑性變形過(guò)程中會(huì)產(chǎn)生比較明顯的剪脹現(xiàn)象，而用來(lái)描述這一現(xiàn)象的較常用的力學(xué)參數(shù)就是剪脹角 Ψ，
上述模型可以通過(guò)USDFLD子程序進(jìn)行實(shí)現(xiàn)，流程圖如下
有限元模型如下圖所示
計(jì)算得到的應(yīng)力云圖及不同圍壓下的載荷位移響應(yīng)如下所示
參考文獻(xiàn)：張力偉,賈善坡,鄒江濤,舒婧曦.泥巖的峰后軟化力學(xué)模型.中國(guó)科技論文,2016,11(21):2456-2461
有關(guān)于abaqus子程序開發(fā)的相關(guān)問(wèn)題可以通過(guò)公眾號(hào)聯(lián)系我們.
公眾號(hào): 320科技工作室
展開 abaqus基于usdfld子程序的內(nèi)聚力疲勞模型
基于usdfld的內(nèi)聚力疲勞模型
之前在Abaqus Cohesive單元的疲勞UMAT － 技術(shù)鄰 (jishulink.com) 上介紹了通過(guò)umat子程序來(lái)編寫內(nèi)聚力疲勞本構(gòu)，實(shí)現(xiàn)裂紋疲勞擴(kuò)展的方法。在實(shí)際計(jì)算中經(jīng)常出現(xiàn)不收斂的情況，因此重新編寫了雙線性本構(gòu)下的usdfld內(nèi)聚力疲勞子程序。
在不考慮疲勞損傷的情況下，單一裂紋模式雙線性內(nèi)聚力本構(gòu)如圖所示
混合模式下的斷裂準(zhǔn)則采用BK準(zhǔn)則
損傷萌生和失效對(duì)應(yīng)的等效張開位移分別為
和
疲勞損傷采用roe提出的損傷演化方程
考慮疲勞損傷后的內(nèi)聚力本構(gòu)如圖所示
這里同樣假設(shè)卸載以及法向壓縮不會(huì)累積疲勞損傷。
建立三點(diǎn)彎曲模型對(duì)疲勞裂紋擴(kuò)展進(jìn)行了模擬，計(jì)算結(jié)果如圖所示
斷裂過(guò)程
跨中載荷位移曲線
損傷演化過(guò)程
跨中底部單元的應(yīng)力應(yīng)變關(guān)系
更新預(yù)告：早期混凝土熱-濕-力多場(chǎng)耦合分析，編寫了基于水化度理論和考慮熱學(xué)參數(shù)變化的溫度場(chǎng)計(jì)算子程序（umatht和film）; 考慮溫度對(duì)濕度擴(kuò)散系數(shù)影響的濕度場(chǎng)計(jì)算子程序（umatht和film）; 基于成熟度理論和雙冪徐變函數(shù)的應(yīng)力場(chǎng)子程序（umat）。
展開 Abaqus用戶子程序USDFLD調(diào)用GETVRM 返回值為0
.0)THEN WRITE(6,*) 'REQUEST ERROR IN USDFLD FOR ELEMENT NUMBER ', 1 NOEL,'INTEGRATION POINT NUMBER ',NPT ENDIF C RETURN END
[原創(chuàng)]Abaqus傳熱分析HETVAL和USDFLD子程序聯(lián)合
HETVAL和USDFLD子程序聯(lián)合
---- 公眾號(hào)‘CAE仿真實(shí)驗(yàn)室’出品
1、HETVAL簡(jiǎn)介
Abaqus傳熱分析中，Hetval子程序常用于定義內(nèi)部熱源，該子程序可以定義材料內(nèi)部的生熱量(如相變、高分子結(jié)晶生熱)，并且能夠調(diào)用state variable 狀態(tài)變量，可以與USDFLD子程序聯(lián)合使用。其函數(shù)體如下：
SUBROUTINE HETVAL(CMNAME,TEMP,TIME,DTIME,STATEV,FLUX,
1 PREDEF,DPRED)
INCLUDE 'ABA_PARAM.INC'
CHARACTER*80 CMNAME
DIMENSION TEMP(2),STATEV(*),PREDEF(*),TIME(2),FLUX(2),
1 DPRED(*)
User coding to define FLUX and update statev !此處定義你的Flux() 更新狀態(tài)變量
RETURN
END
2、模型 1mx1mx0.5m的方塊，進(jìn)行熱應(yīng)力耦合分析，選用C3D8T單元
3、Hetval 和USDFLD聯(lián)合
USDFLD和HETVAL寫在一個(gè)for文件中，HETVAL中通過(guò)State(1)來(lái)定義熱源，而state(1)是通過(guò)USDFLD更新，實(shí)現(xiàn)數(shù)據(jù)傳遞
4、材料
設(shè)置材料為鋼，給定力、熱屬性。
展開 
Abaqus通過(guò)USDFLD子程序進(jìn)行泥巖的應(yīng)變軟化模擬
本文基于泥巖的三軸壓縮試驗(yàn)曲線，建立考慮應(yīng)變軟化特性的泥巖彈塑性本構(gòu)模型，使用Abaqus及其子程序對(duì)泥巖的三軸壓縮試驗(yàn)進(jìn)行了數(shù)值模擬。泥巖在受壓過(guò)程中主要經(jīng)過(guò)了5個(gè)階段，即壓密階段、彈性變 形、應(yīng)變硬化、應(yīng)變軟化、殘余階段。泥巖應(yīng)變軟化模型如下所示。
式中，ξ為強(qiáng)度參數(shù)，ξp為峰值強(qiáng)度參數(shù)，ξr為殘余階段強(qiáng)度參數(shù)，η為應(yīng)變軟化參數(shù)，η*為殘余階段的應(yīng)變軟化參數(shù)初始值。對(duì)于三軸壓縮試驗(yàn)，η用塑性剪切應(yīng)變來(lái)表示
塑性屈服準(zhǔn)則采用Mohr-Coulomb準(zhǔn)則，則粘聚力和內(nèi)摩擦角的參數(shù)演化可以用下式表示
在巖石的塑性變形過(guò)程中會(huì)產(chǎn)生比較明顯的剪脹現(xiàn)象，而用來(lái)描述這一現(xiàn)象的較常用的力學(xué)參數(shù)就是剪脹角 Ψ，
上述模型可以通過(guò)USDFLD子程序進(jìn)行實(shí)現(xiàn)，流程圖如下
有限元模型如下圖所示
計(jì)算得到的應(yīng)力云圖及不同圍壓下的載荷位移響應(yīng)如下所示
參考文獻(xiàn)：張力偉,賈善坡,鄒江濤,舒婧曦.泥巖的峰后軟化力學(xué)模型.中國(guó)科技論文,2016,11(21):2456-2461
有關(guān)于abaqus子程序開發(fā)的相關(guān)問(wèn)題可以聯(lián)系扣扣1653004885或者關(guān)注cae320公眾號(hào)
展開 Abaqus子程序代碼分享
Advances in Engineering Software, 105: 9-16 (2017)
13、Abaqus2Matlab_Toolbox2.00.zip
14、ABAQUS USDFLD SUBROUTINE FOR FUNCTIONALLY GRADED ELEMENTS
ABAQUS USDFLD子程序用于功能性漸變?cè)?ABAQUS USDFLD Subroutine for the implementation of a continuous variation of the material elastic properties between integration points.
ABAQUS USDFLD子例程，用于在積分點(diǎn)之間連續(xù)改變材料的彈性。
Paper: E. Martínez-Pa?eda. On the finite element implementation of functionally graded materials. Materials, 12(2): 287 (2019)
14、FGMcode.zip
15、NON-LINEAR X-FEM CODE FOR MSG PLASTICITY
MSG可塑性的非線性X-FEM代碼
Non-linear eXtendend Finite Element code developed in MATLAB including the following material models: linear elasticity, von Mises plasticity and conventional mechanism-based strain gradient (CMSG) plasticity.
展開 Abaqus聯(lián)用USDFLD+HETVAL子程序分析慢速烤燃 ￥99
采用HETVAL子程序的一步分解慢速烤燃模型請(qǐng)參考：https://www.yqgqt.org.cn/content/post/1300451
然而有時(shí)反應(yīng)并不是一步完成的，需要多步反應(yīng)模型，此時(shí)HETVAL并不夠用。為此，調(diào)用強(qiáng)大的USDFLD子程序跟HETVAL聯(lián)用解決問(wèn)題。當(dāng)然，USDFLD非常強(qiáng)大，本例僅演示其與HETVAL的聯(lián)用并分析慢速烤燃過(guò)程。
本例提供HETVAL、USDFLD子程序、cae及inp源文件，咨詢請(qǐng)聯(lián)系我 QQ：180280578。
附件包括詳細(xì)建模過(guò)程、材料參數(shù)、注意事項(xiàng)及子程序。
展開 USDFLD子程序在基于ABAQUS開展多孔介質(zhì)（油氣工程為例）流-固耦合分析中的應(yīng)用
為此，本教程開發(fā)啦USDFLD子程序用于實(shí)現(xiàn)上述目的。
1、隨著油氣資源的開發(fā)，儲(chǔ)層等多孔介質(zhì)應(yīng)力和孔壓等會(huì)發(fā)生變化。同時(shí)，油氣儲(chǔ)層等多孔介質(zhì)的屬性受多種因素影響，例如應(yīng)力和孔隙壓力；
2、ABAQUS默認(rèn)條件下，Property模塊只能設(shè)置恒定材料屬性，沒(méi)法反映儲(chǔ)層特性隨應(yīng)力和孔壓等因素對(duì)材料屬性的影響；
3、USDFLD子程序是最常用的子程序，通過(guò)設(shè)置和編程可以實(shí)現(xiàn)對(duì)油氣開發(fā)過(guò)程中材料屬性的動(dòng)態(tài)控制，進(jìn)而更準(zhǔn)確的模擬工程實(shí)際情況；
4、使用USDFLD子程序時(shí)，常規(guī)屬性（彈性模量等力學(xué)屬性）可以通過(guò)界面直接完成，而滲流參數(shù)（滲透系數(shù)和孔隙比）則需要修改INP文件或Keywords實(shí)現(xiàn)，本教程給出了實(shí)現(xiàn)方法；
5、通常，使用USDFLD開展有限元模擬時(shí)基本全是設(shè)置一個(gè)場(chǎng)變量，本教程給出了應(yīng)力和孔壓同時(shí)（兩個(gè)場(chǎng)變量）影響儲(chǔ)層滲透率時(shí)的USDFLD設(shè)置方法；
6、本教程可以用于油氣開發(fā)過(guò)程，也可以用于涉及到多孔介質(zhì)流固耦合分析的其他領(lǐng)域；
7、第一次錄制視頻教程，有瑕疵和紕漏，請(qǐng)大家提出講得不清楚的地方，或不理解的地方，以便在后續(xù)過(guò)程中更新教程.
感謝大家的支持！！
HTTP：USDFLD子程序在基于ABAQUS開展多孔介質(zhì)（油氣工程為例）流-固耦合分析中的應(yīng)用
展開 