abaqus 子程序 熱源模型的案例
abaqus傳統(tǒng)攪拌摩擦焊接熱源Fortran子程序和模型inp文件 ￥19.89
abaqus傳統(tǒng)攪拌摩擦焊接熱源Fortran子程序和模型inp文件
Abaqus模擬焊接（雙橢球熱源）的方法與步驟詳解（內(nèi)含完整雙橢球熱源子程序） ￥1.7
Abaqus模擬焊接（雙橢球熱源）的方法與步驟詳解（內(nèi)含完整雙橢球熱源子程序）
三角形熱源 磨削子程序 abaqus
SUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS,
1 JLTYP,TEMP,PRESS,SNAME)
C
INCLUDE 'ABA_PARAM.INC'
C
DIMENSION FLUX(2), TIME(2), COORDS(3)
CHARACTER*80 SNAME
qm = 3.2321*10000
lc = 2
dx = (4000/60)*time(1)
distance=coords(1)-dx
if(distance>=0.0.and.distance<=2)then
Flux(1) = 2*qm*distance/2
else
Flux(1)=0
end if
RETURN
END
展開(kāi) abaqus_子程序_高斯熱源 ￥15
前期準(zhǔn)備(這個(gè)可以百度，有更詳細(xì)的安裝教程)：
首先需要安裝兩個(gè)軟件（InterFortran和vs studio，注意ABAQUS與軟件版本號(hào)的兼容性，不然會(huì)出現(xiàn)配置不成功的情況）。鑒于題主用的是abauqs2016，此處以abauqs2016為例子簡(jiǎn)述步驟。
Abaqus2016的兼容版本為VS2012，InterFortran2013（需要lic許可證，要不然只能用一年_maybe）.
安裝過(guò)程中需要記錄各自的安裝地址，軟件安裝完成后，在abaqus中的一個(gè)launcher.bat文件中將上述兩個(gè)文件的路徑寫(xiě)入該文件（注意路徑要正確，同時(shí)嚴(yán)格按照格式書(shū)寫(xiě)，禁止使用中文字符，禁止多莫名其妙的空格）。
顯示上圖所示即配置成功。
建模（建模過(guò)程僅敘述關(guān)鍵點(diǎn)）：
材料屬性：必須輸入用于溫度場(chǎng)計(jì)算的熱物性參數(shù)。將絕對(duì)零度設(shè)定為-273.15，此時(shí)計(jì)算出來(lái)的溫度單位為℃。
如果只是單純的溫度場(chǎng)分析,將分析步類型設(shè)定為熱傳遞，分析步時(shí)間按照熱源移動(dòng)速度，模型大小計(jì)算得出。
邊界條件：根據(jù)需求設(shè)定輻射、散熱系數(shù)。荷載選擇表面熱流，作用區(qū)域定位在熱源照射面。預(yù)定義場(chǎng)設(shè)定初始溫度。
網(wǎng)格類型為DC3D8。
提交計(jì)算：
提交作業(yè)，在通用模塊用戶子程序一欄鍵入.for子程序路徑。
嗯，大概就這樣。
付費(fèi)內(nèi)容包括子程序文件，CAE模型，聯(lián)系方式。
展開(kāi) 
abaqus焊接雙熱源同步焊接子程序
熱源子程序
ABAQUS任意路徑移動(dòng)熱源Dflux子程序編寫(xiě) ￥20
ABAQUS復(fù)雜路徑雙橢球體熱源Dflux子程序，直線-圓弧-斜線，平面坐標(biāo)變換
Abaqus子程序系列：UMDFLUX（定義多個(gè)非均勻分布熱源）
子程序UMDFLUX
· 可用于描述是位置、時(shí)間、溫度等函數(shù)的多個(gè)移動(dòng)或靜止熱通量；
· 可以使用和更新單元解相關(guān)的變量；
· 每個(gè)單元可以調(diào)用一次，來(lái)描述單元內(nèi)容開(kāi)始點(diǎn)和結(jié)束點(diǎn)之間的移動(dòng)熱源；
· 施加熱通量為點(diǎn)移動(dòng)熱源，單位是能量單位（JT-1）;
· 在用戶子程序外，不需要體積積分；
1. 與子程序DFLUX的區(qū)別
子程序DFLUX
· 可以用來(lái)在熱傳導(dǎo)或質(zhì)量擴(kuò)散分析中，定義一個(gè)非均勻分布的熱通量，可以是位置，時(shí)間，溫度，單元號(hào)，積分點(diǎn)號(hào)等的函數(shù)；
· 對(duì)于分析中定義了基于單元或基于表面（僅僅熱傳導(dǎo)）的非均勻分布熱通量，每個(gè)熱通量積分點(diǎn)調(diào)用；
· 忽略任何可能出現(xiàn)的與非均勻分布通量定義相關(guān)的幅值參考;
· 將節(jié)點(diǎn)作為一階傳熱單元、一階溫度-位移耦合單元、一階熱-電-結(jié)構(gòu)耦合單元和質(zhì)量擴(kuò)散單元的通量積分點(diǎn)。
SUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS,JLTYP,TEMP,PRESS,SNAME)
C
INCLUDE 'ABA_PARAM.INC'
C
DIMENSION FLUX(2), TIME(2), COORDS(3)
CHARACTER*80 SNAME
user coding to define FLUX(1) and FLUX(2)
RETURN
END
2.
展開(kāi) Abaqus移動(dòng)加載高斯熱源子程序，適于焊接、3D打印領(lǐng)域 ￥2.5
Abaqus移動(dòng)加載高斯熱源子程序，適于焊接、3D打印領(lǐng)域，fortran編制的子程序見(jiàn)“付費(fèi)后”的附件中。
利用Abaqus用戶子程序實(shí)現(xiàn)制動(dòng)盤(pán)熱分析中的熱源移動(dòng)
在進(jìn)行制動(dòng)器熱分析時(shí)，若要分析整個(gè)制動(dòng)過(guò)程或多次制動(dòng)后的結(jié)果，直接采用接觸摩擦生熱的方式受到制動(dòng)時(shí)間以及大位移接觸的影響，計(jì)算較困難，若采用結(jié)構(gòu)靜態(tài)傳熱又無(wú)法實(shí)現(xiàn)隨制動(dòng)盤(pán)的轉(zhuǎn)動(dòng)熱源的移動(dòng)。Abaqus中可通過(guò)對(duì)結(jié)構(gòu)設(shè)置質(zhì)量流率（MASS FLOW），采用熱傳導(dǎo)分析實(shí)現(xiàn)熱源移動(dòng)。需要注意的是，MASS FLOW不支持CAE，需要通過(guò)關(guān)鍵字設(shè)置。下面以汽車制動(dòng)盤(pán)熱分析說(shuō)明整個(gè)過(guò)程及MASS FLOW的使用方法。
首先僅建立制動(dòng)盤(pán)的模型即可，不需要制動(dòng)片模型。如下圖所示：
分別設(shè)置制動(dòng)盤(pán)的材料屬性，設(shè)置熱傳導(dǎo)分析步（Heat transfer）。邊界為初始溫度場(chǎng)，對(duì)流換熱、熱輻射等。載荷為表面熱流密度，施加在初始制動(dòng)片與制動(dòng)盤(pán)接觸的面積內(nèi)。其中需要考慮熱流密度隨制動(dòng)盤(pán)速度的變化，對(duì)流換熱系數(shù)隨速度的變化，以及對(duì)流換熱系數(shù)沿制動(dòng)盤(pán)徑向的變化等。如下圖所示為熱流密度施加的區(qū)域。
接著設(shè)置質(zhì)量流率（MASS FLOW），首先需要指定質(zhì)量流率施加的區(qū)域，此模型中選擇制動(dòng)盤(pán)的所有節(jié)點(diǎn)，將其設(shè)置為一個(gè)set，給定相應(yīng)的name（后續(xù)將用到此set）。確認(rèn)材料、分析步、邊界、載荷等設(shè)置正確后可輸出inp文件添加質(zhì)量流率關(guān)鍵字或直接在CAE中添加關(guān)鍵字。
在step后添加關(guān)鍵字，如下圖所示。添加完成后保存inp文件。
接著是最關(guān)鍵的部分，質(zhì)量流率需要使用用戶子程序實(shí)現(xiàn)，子程序中需要給出由制動(dòng)盤(pán)速度變化引起的質(zhì)量流率的變化，以及沿制動(dòng)盤(pán)徑向的質(zhì)量流率的變化，同時(shí)若為多次制動(dòng)，需要區(qū)分制動(dòng)、停止的過(guò)程，僅在制動(dòng)過(guò)程施加質(zhì)量流率。具體的用戶子程序如下所示:
用戶子程序
求解時(shí)調(diào)用此子程序計(jì)算即可。整個(gè)分析過(guò)程可看到熱源在制動(dòng)盤(pán)上移動(dòng)，計(jì)算后的溫度場(chǎng)結(jié)果如下圖所示。
來(lái)源：有限元在線的博客，版權(quán)歸作者所有。
展開(kāi) ABAQUS焊接模擬-Python編寫(xiě)移動(dòng)高斯熱源子程序-不帶生死單元
ABAQUS平板對(duì)接-不帶生死單元。Python編寫(xiě)移動(dòng)高斯熱源子程序（包括高斯面熱源、雙橢球熱源）
模型作如下假設(shè)：材料為各向同性材料，不考慮熔池流動(dòng)及相變影響。
考慮到過(guò)來(lái)學(xué)習(xí)的大多都是和我一樣的學(xué)生黨，因此設(shè)置了一個(gè)大家都能接受的價(jià)格。
如果視頻中有什么錯(cuò)誤或沒(méi)講清的大家可以留言！！
abaqus剪切修正GTN模型的VUMAT子程序開(kāi)發(fā)
這里對(duì)相應(yīng)的算法進(jìn)行簡(jiǎn)要說(shuō)明：
NH-GTN模型
屈服函數(shù)：
其中等效孔洞體積分?jǐn)?shù)定義為：
孔洞體積分?jǐn)?shù)包含新孔隙形核，原有空隙生長(zhǎng)以及剪切相關(guān)的等效體積分?jǐn)?shù)增加：
形核，生長(zhǎng)，剪切相關(guān)體積分?jǐn)?shù)的演化遵循：
其中：
剪切效應(yīng)的修正，考慮應(yīng)力狀態(tài)的影響
參數(shù)的物理含義如下
通過(guò)將文獻(xiàn)中的數(shù)值算法編程實(shí)現(xiàn)在VUMAT子程序中，可以用來(lái)實(shí)現(xiàn)對(duì)延性金屬材料在不同應(yīng)力狀態(tài)下的損傷演化進(jìn)行合理的數(shù)值預(yù)測(cè)，應(yīng)用于金屬成型領(lǐng)域（沖壓，軋制，擠壓等）
預(yù)測(cè)修正后的模型應(yīng)該在簡(jiǎn)單拉伸情況下于abaqus自帶的GTN模型保持相同的損傷和其他狀態(tài)變量的分布，并在剪切情況中損傷發(fā)展顯著高于abaqus自帶的模型（自帶的模型忽略了剪切效應(yīng)）。（為了進(jìn)行對(duì)比使用于自帶的本構(gòu)相同的硬化方式，模擬中使用了相同的質(zhì)量縮放，但質(zhì)量縮放容易產(chǎn)生數(shù)值振蕩，模擬的拉伸曲線存在波動(dòng)。）
初步模擬結(jié)果：
拉伸情況（abaqus-VUMAT）
應(yīng)力分情況
孔洞體積分?jǐn)?shù)
剪切模型（abaqus-VUMAT）
不同變形時(shí)刻的應(yīng)力分布
T=0.1s
局部放大圖
T=0.5s
局部放大圖
T=0.6s
局部放大圖
可以看到模型在拉伸預(yù)測(cè)中與原始模型保持一致，而在剪切修正后損傷發(fā)展顯著快于原始模型，利用作者提出的方法可以應(yīng)用于復(fù)雜應(yīng)力狀態(tài)下金屬材料的損傷分析，相關(guān)參數(shù)部分參考文獻(xiàn)，其中Kw=3.T1=0.2，T2=0.7.模擬結(jié)果符合文獻(xiàn)所提出方法的基本趨勢(shì)。
最后，如果有相關(guān)需要?dú)g迎通過(guò)公眾號(hào)“320科技工作室”聯(lián)系我們。
展開(kāi) 
Abaqus 修正劍橋模型的vumat子程序開(kāi)發(fā)
巖土工程中廣泛采用的莫爾–庫(kù)侖準(zhǔn)則，能較好地描述巖土材料的強(qiáng)度特性和破壞行為，而涉及變形和固結(jié)分析問(wèn)題時(shí)該模型從理論和實(shí)踐中均有一定的不足，基于等向加載試驗(yàn)建立起來(lái)的劍橋模型則能較好地描述黏土的彈塑性變形特性。
傳統(tǒng)相關(guān)聯(lián)修正劍橋本構(gòu)模型的屈服面方程、硬化準(zhǔn)則和流動(dòng)法則分別為
由于式(1)所確定的屈服軌跡在p平面是一個(gè) 圓，不能反映巖土介質(zhì)拉壓不等(S-D)效應(yīng)，而且劍 橋模型是基于正常固結(jié)狀態(tài)試驗(yàn)推導(dǎo)而來(lái)，子午面上臨界狀態(tài)線通過(guò)應(yīng)力坐標(biāo)原點(diǎn)，表現(xiàn)為不考慮土 體黏聚力的純摩擦型本構(gòu)，而大多數(shù)的巖土介質(zhì)具 有一定的黏聚力，屬摩擦–黏聚型材料。考慮到莫 爾–庫(kù)侖準(zhǔn)則有此特征，故將莫爾–庫(kù)侖準(zhǔn)則與劍橋模型相結(jié)合，得
聯(lián)合式(1)，(4)，(5)，有改進(jìn)的屈服函數(shù)：
根據(jù)上述理論，可以通過(guò)以下流程圖完成VUMAT子程序編寫(xiě)
通過(guò)單胞模型計(jì)算得到得結(jié)果如下圖所示。
[1]袁克闊,陳衛(wèi)忠,于洪丹,譚賢君,趙武勝,李香玲.考慮黏聚特性和拉壓不等效應(yīng)的修正劍橋模型及數(shù)值實(shí)現(xiàn)[J].巖石力學(xué)與工程學(xué)報(bào),2012,31(08):1574-1579.
展開(kāi) 帶內(nèi)熱源的UMATHT子程序入門(mén)小例子 ￥5
自己探索了一段時(shí)間，下面將在abaqus幫助文檔里的例子基礎(chǔ)上。結(jié)合傳熱學(xué)的理論，介紹在UMATHT實(shí)現(xiàn)內(nèi)熱源生熱，并與用abaqus的HETVAL子程序的結(jié)果做對(duì)比，驗(yàn)證正確性。（注：這個(gè)例子很簡(jiǎn)單，大家可以根據(jù)需要進(jìn)一步開(kāi)發(fā)）</p><h2>1.導(dǎo)熱微分方程</h2><p> 為了更好得理解UMATHT需要定義變量，先介紹一下一維導(dǎo)熱微分方程（三維可類似推廣）。
展開(kāi) Abaqus修正GTN模型的VUMAT子程序
Gurson-Tvergaard-Needleman ( GTN) 模型是研究金屬損傷的重要工具。GTN 模型通過(guò)孔洞體積分?jǐn)?shù)的演變來(lái)判斷材料的失效， 但不適用于剪切斷裂為主的韌性斷裂。本文在GTN模型中引入剪應(yīng)力的影響，編寫(xiě)了相關(guān)的VUMAT子程序。
GTN模型的屈服函數(shù)可以用下式表示
其中q1，q2是模型參數(shù)，取q1=1.5，q2=1,σ0為等效應(yīng)力，p為靜水應(yīng)力，q為Mises等效應(yīng)力；f為空洞的體積分?jǐn)?shù)。
p和q可以通過(guò)徑向返回算法得到
應(yīng)變控制的孔洞形核系數(shù)
GTN模型可以通過(guò)以下4個(gè)方程進(jìn)行描述
Nahshon and Hutchinson考慮了剪應(yīng)力對(duì)模型的影響
于是孔隙體積分?jǐn)?shù)的演化可以通過(guò)下式描述
仿真計(jì)算得到的結(jié)果如下圖所示
有問(wèn)題私信或者關(guān)注cae320公眾號(hào)
展開(kāi) ABAQUS三維hill48彈塑性模型VUmat子程序（彈性為正交各向異性） ￥388
1.ABAQUS三維hill48彈塑性模型VUmat子程序
2.彈性階段為正交各項(xiàng)異性材料
3.hill48和正交各項(xiàng)異性材料參數(shù)參考ABAQUS靜力模塊自帶的模型參數(shù)
4.發(fā)貨方式為百度網(wǎng)盤(pán)鏈接，包含子程序及上面跑的兩個(gè)模型相關(guān)文件，包含Cae，inp文件，odb文件等
5.ABAQUS版本為2024，低版本可以利用導(dǎo)入inp文件的方式運(yùn)行及修改
6.可以免費(fèi)答疑三次，后續(xù)添加你自己的模型或者相關(guān)參數(shù)等輔導(dǎo)都可以優(yōu)惠。
展開(kāi) 