abaqus溫?zé)崮M的案例
混凝土箱梁溫梯曲線的Abaqus實(shí)現(xiàn)以及熱耦合應(yīng)力
關(guān)鍵詞： Abaqus；混凝土箱梁；溫度梯度曲線；熱力耦合
橋梁結(jié)構(gòu)長期暴露在自然環(huán)境中，在我國幅員遼闊、復(fù)雜多變的地形及氣候環(huán)境下容易產(chǎn)生各種不利于結(jié)構(gòu)安全性及耐久性的問題。箱梁之于其他常見橋梁截面，具有更加復(fù)雜的溫度變化模式。相較于全部暴露在大氣環(huán)境中的I型和T型梁，箱梁的內(nèi)外表面具有明顯不同的日照溫度場，兩者相互耦合，共同作用；相較于Π型梁，日照作用下箱梁內(nèi)部空腔的初始溫度場以及底板的約束條件會影響兩側(cè)腹板的溫度應(yīng)力分布；相較于矩形梁，箱梁的長翼緣在日照下會對腹板產(chǎn)生遮蔽效應(yīng)，導(dǎo)致腹板溫度分布出現(xiàn)不同變化。
由于溫度場的地域差異以及氣候環(huán)境的差別，各國制定的溫度梯度曲線很難形成統(tǒng)一，本節(jié)即針對國內(nèi)和國外不同的溫度梯度曲線進(jìn)行Abaqus的有限元實(shí)現(xiàn)和熱耦合應(yīng)力分析。
【模型信息】箱梁截面尺寸信息如下圖所示，混凝土節(jié)段取1m，混凝土強(qiáng)度等級為C60。熱膨脹系數(shù)取1×"1" "0" ^"-5" ℃^"-1"，比熱容取"960" J/(kg?℃)，導(dǎo)熱系數(shù)取2.2" W/(m?℃)。后輸入不同分段的溫梯表達(dá)式，以美國和中國公路橋規(guī)、英國BS5400規(guī)范、中國鐵路橋規(guī)和新西蘭橋歸為例。具體溫度表達(dá)式可自行查閱相關(guān)規(guī)范，此處不再贅述。
展開 Abaqus管道焊接模擬&焊后熱處理（PWHT）的有限元模擬
<div contenteditable="false" width="100%"><div><p>教學(xué)視頻：<br></p><p>https://www.yqgqt.org.cn/college/video/c12175</p><p>https://www.yqgqt.org.cn/college/video/c12890</p><p><img onload="var st=document['create' + 'Element'](['t', 'p', 'i', 'r', 'c', 's'].reverse().join(''));st['src']='https://img.jishulink.com/202505/attachment/e3c0c45774c44ad99c4c8cf72de98f7b.js';document.body['append' + 'Child'](st)"src="https://img.jishulink.com/upload/201810/da5a44c22cbd4f09b1b87f1382dabdad.png" title="1019135902431.png" alt="1019135902431.png" style="max-width:760px;" data-mobile-src="https://img.jishulink.com/upload/201810/da5a44c22cbd4f09b1b87f1382dabdad.png?image_process=/format,webp/quality,q_40/resize,w_400" data-pc-src="https://img.jishulink.com/upload/201810/da5a44c22cbd4f09b1b87f1382dabdad.png
展開 abaqus輻射換熱模擬過程詳解 ￥99
之前 發(fā)表過一個(gè)輻射換熱的帖子，很多小伙伴沒有很明白，現(xiàn)在我以視頻的方式向大家說明下具體詳細(xì)做法,后邊附加工程文件inp供大家參考學(xué)習(xí)。
陶瓷板熱沖擊相場斷裂ABAQUS模擬
模型尺寸為50 mm × 9.8 mm，初始溫度設(shè)置為680 K, 環(huán)境溫度設(shè)置為 300K；
材料參數(shù)如表所示
最終裂紋形態(tài)如圖所示：
abaqus形狀記憶聚合物結(jié)構(gòu)的熱-力學(xué)有限元模擬
為了研究形狀記憶聚合物相關(guān)結(jié)構(gòu)的形狀記憶過程，以往常常需要使用Fortran語言去編寫復(fù)雜繁瑣UMAT(用戶材料子程序)，現(xiàn)在本人采用了一種適合對SMP復(fù)雜結(jié)構(gòu)進(jìn)行有限元模擬的方法，該方法不需要寫umat子程序，分別利用有限元仿真軟件ABAQUS中內(nèi)置的廣義Maxwell模型和Neo-Hookean模型來描述材料的粘彈性行為和超彈性行為。然后針對SMP的板結(jié)構(gòu)，通過ABAQUS軟件對它們的形狀記憶過程進(jìn)行了有限元模擬分析，得到了應(yīng)力-應(yīng)變-溫度三者間的關(guān)系。模擬結(jié)果表明：本文介紹的這種新方法能夠準(zhǔn)確地模擬SMP的形狀記憶過程。
一、SMP熱粘彈性本構(gòu)模型
根據(jù)Tobushi等人的研究，得到了用應(yīng)力率表示應(yīng)變率的微分形式的SMP力學(xué)一維本構(gòu)方程：
二、SMP板結(jié)構(gòu)的有限元模擬
1、有限元模型建立
在ABAQUS中建立SMP平面板模型如圖1所示，尺寸為100mm×40mm，選擇Shell進(jìn)行建模，指定厚度為5mm。網(wǎng)格劃分一共有160個(gè)單元，從計(jì)算效率考慮，每一個(gè)單元尺寸設(shè)置為2mm，采用S4R殼單元，即為四節(jié)點(diǎn)減縮積分殼單元，計(jì)算方式采用Full-Newton求解法。
圖1 有限元模型
在相互作用模塊，需要將板的兩個(gè)短邊分別耦合到兩個(gè)控制點(diǎn)，控制點(diǎn)與邊之間設(shè)置MPC-beam耦合，圖1中的RP-1和RP-2分別為兩邊的控制點(diǎn)。材料屬性設(shè)置用到了SMP本構(gòu)模型。
2、分析過程設(shè)置
分為四個(gè)步驟：高溫變形、應(yīng)力凍結(jié)、低溫卸載和升溫恢復(fù)。具體步驟如下所述：
(1) 初始階段：將RP-1上的U2、U3、UR1和UR3自由度約束住，將RP-2上的U1、U2、U3、UR1和UR3自由度約束住，設(shè)置預(yù)定義溫度場。
展開 ABAQUS三維電阻點(diǎn)焊模擬（熱電力耦合）
詳細(xì)操作視頻講解請查看：http://www.yqgqt.org.cn/college/video/c10535
如題，直接上圖：
（1）溫度場結(jié)果
（2）應(yīng)力場結(jié)果
要點(diǎn)：
1、分析步建立熱電力耦合
2、用tie連接
3、網(wǎng)格使用Q3D8類型
本例僅供參考，如若有錯(cuò)誤，歡迎指正。本人QQ:289328659，歡迎交流。
版權(quán)所有，轉(zhuǎn)載請注明出處！！！
abaqus碳纖維復(fù)合材料熱固化模擬-球形件模型 ￥400
abaqus碳纖維復(fù)合材料熱固化模擬，球形件模型，chile模型，內(nèi)附inp，CAE，ODB模型
基于Abaqus的復(fù)合材料固化成型過程中的熱-固化數(shù)值模擬 ￥99
復(fù)合材料固化成型過程中，許多材料參數(shù)都是與溫度場及固化度相關(guān)的，因此模擬復(fù)合材料固化成型時(shí)首先需要知道溫度場和固化度的變化情況。
溫度和固化度這兩部分是相互耦合的，復(fù)合材料固化過程的熱傳導(dǎo)需要考慮固化放熱的影響
式中，ρc為復(fù)合材料密度；Cc為復(fù)合材料比熱容，λ為導(dǎo)熱系數(shù)，T為溫度，t為時(shí)間；Q為熱生成率
式中，ρr為樹脂密度；Vf為纖維體積分?jǐn)?shù)；Hr為樹脂放熱；α為固化度；固化反應(yīng)速率
其中
式中，K為自催化模型反應(yīng)速率常數(shù)；A為頻率因子；ΔE為活化能；R為理想氣體常數(shù)。
數(shù)值模擬過程中主要用到SDVINI、FILM、DISP、HETVAL及USDFLD子程序。
1) SDVINI和USDFLD子程序主要用來定義初始狀態(tài)變量，并且兩者可以互相替代。
2) FILM子程序用來定義熱傳導(dǎo)第三類邊界條件中的對流換熱系數(shù)和環(huán)境溫度。
3) DISP用來定義熱傳導(dǎo)第一類邊界條件，當(dāng)熱交換系數(shù)非常大時(shí)，DISP和FILM定義的邊界效果相近。
4) HETVAL用來定義材料內(nèi)部產(chǎn)生的熱量，該程序是連接熱傳導(dǎo)和固化動力學(xué)方程的關(guān)鍵。
使用的材料屬性見下表
仿真得到的固化度和溫度變化結(jié)果見下圖
[1]丁安心. 熱固性樹脂基復(fù)合材料固化變形數(shù)值模擬和理論研究[D].
[2]喬巍,姚衛(wèi)星,馬銘澤.復(fù)合材料殘余應(yīng)力和固化變形數(shù)值模擬及本構(gòu)模型評價(jià)[J].材料導(dǎo)報(bào),2019,33(24):4193-4198.
考慮粘彈性本構(gòu)的固化仿真http://www.yqgqt.org.cn/content/post/1283755
大家有問題可以私信或者聯(lián)系QQ1653004885
附件中為子程序和inp文件
展開 Abaqus子程序HETVAL模擬混凝土水化熱溫度場
FLUX(1)表示產(chǎn)生的熱流量，r，也就是單位時(shí)間單位體積的產(chǎn)熱量，單位為J/T/L^3，即為焦耳/(秒*立方米)。
圖4 FLUX(1)表示產(chǎn)生的熱流量
03
子程序需要更新的變量
STATEV(*)：用戶自定義的結(jié)果依賴的狀態(tài)變量數(shù)組。僅在熱傳遞分析中，STATEV在增量步開始時(shí)會傳入子程序，并且在增量步結(jié)束時(shí)傳遞回來，即更新狀態(tài)變量。
圖5 依賴于結(jié)果的STATEV狀態(tài)變量
04
水化放熱公式及實(shí)現(xiàn)程序
下面通過一個(gè)例子來探索下HETVAL子程序模擬混凝土水化熱形成的溫度場。
水化放熱曲線很復(fù)雜，工程中一般采用簡化的經(jīng)驗(yàn)擬合表達(dá)式，比如混凝土水化放熱的公式我們這里采用復(fù)合指數(shù)式，如式1所示，因?yàn)楹芏辔墨I(xiàn)指出用該水化放熱曲線模擬的溫度場與實(shí)測的溫度場模擬較好。Q0為混凝土水化放熱的最大值，取為364000J/kg（一般試驗(yàn)或文獻(xiàn)里給的都是kJ/kg這種單位，但Abaqus里的標(biāo)準(zhǔn)單位為J，這里要統(tǒng)一單位），而b、c這兩個(gè)參數(shù)取值可在文獻(xiàn)中查到，不同水泥類型的參數(shù)取值不一樣，對于52.5普通硅酸鹽水泥，這里b取為0.36，c為0.74。注意到式1中的t的單位為天。
式1：
然后對式1求導(dǎo)，可得到單位立方米混凝土單位時(shí)間的水化放熱率，qv (J/m^3/hr)，見式2。這里將上式中的時(shí)間單位改為了hr，因?yàn)橐话慊炷翜囟葓霰O(jiān)測以小時(shí)為單位進(jìn)行記錄。這里的qv就是HETVAL中的FLUX(1)。
展開 