abaqus 固化變形的案例
碳纖維復合材料平板固化翹曲變形 ￥25
碳纖維復合材料平板固化翹曲變形，內附inp文件，ODB文件及操作視頻
NX_samcef mecano復合材料固化變形回彈_制造工藝分析
NX_samcef mecano復合材料固化變形回彈_制造工藝分析
視頻材料為復合材料的固化變形回彈建模與分析，視頻可作為操作例程。
百度網盤：http://pan.baidu.com/s/1qXSgCEw
優酷：
基于Abaqus的復合材料固化成型過程中的熱-固化數值模擬 ￥99
復合材料固化成型過程中，許多材料參數都是與溫度場及固化度相關的，因此模擬復合材料固化成型時首先需要知道溫度場和固化度的變化情況。
溫度和固化度這兩部分是相互耦合的，復合材料固化過程的熱傳導需要考慮固化放熱的影響
式中，ρc為復合材料密度；Cc為復合材料比熱容，λ為導熱系數，T為溫度，t為時間；Q為熱生成率
式中，ρr為樹脂密度；Vf為纖維體積分數；Hr為樹脂放熱；α為固化度；固化反應速率
其中
式中，K為自催化模型反應速率常數；A為頻率因子；ΔE為活化能；R為理想氣體常數。
數值模擬過程中主要用到SDVINI、FILM、DISP、HETVAL及USDFLD子程序。
1) SDVINI和USDFLD子程序主要用來定義初始狀態變量，并且兩者可以互相替代。
2) FILM子程序用來定義熱傳導第三類邊界條件中的對流換熱系數和環境溫度。
3) DISP用來定義熱傳導第一類邊界條件，當熱交換系數非常大時，DISP和FILM定義的邊界效果相近。
4) HETVAL用來定義材料內部產生的熱量，該程序是連接熱傳導和固化動力學方程的關鍵。
使用的材料屬性見下表
仿真得到的固化度和溫度變化結果見下圖
[1]丁安心. 熱固性樹脂基復合材料固化變形數值模擬和理論研究[D].
[2]喬巍,姚衛星,馬銘澤.復合材料殘余應力和固化變形數值模擬及本構模型評價[J].材料導報,2019,33(24):4193-4198.
考慮粘彈性本構的固化仿真http://www.yqgqt.org.cn/content/post/1283755
大家有問題可以私信或者聯系QQ1653004885
附件中為子程序和inp文件
展開 基于粘彈性本構模型的熱固性樹脂基復合材料固化變形數值仿真模型
背景介紹
熱固性樹脂基復合材料在制件成型過程中會產生殘余應力，引起固化變形，從而增加裝配和制造的難度，因此，合理預測預制件固化過程中的殘余應力的發展具有重要意義。
早期的研究主要集中于彈性理論來研究復材的固化成型，現今，越來越多的文獻考慮了樹脂的固化放熱以及材料的各向異性等因素的影響，發展了基于粘彈性模型的數值仿真計算方法，證明了粘彈性的結果固化變形量小于線彈性的結果，且樹脂含量越高的復材，其粘彈性效果越明顯。
RTM成型工藝示意圖
二。粘彈性模型在Abaqus中的實現
本文作者在參考文獻【1】的基礎上，使用廣義Maxwell粘彈性本構模型，聯合編寫了HETVAL、USDFLD、DISP、UMAT及UEXPAN子程序，在abaqus軟件平臺中實現了復材固化成型的仿真模擬，其基本編程思路如下圖所示：
其中，最關鍵的粘彈性本構公式為：
參考上述公式和子程序的編寫流程，可以完成上述模型。最后得到仿真Mises應力云圖和S33云圖如下：
得到的S33關于時間的曲線趨勢如下所示：
該曲線結果和文獻有出入，但是榮的文獻中關于底數的取值有錯誤，亦即下列公式的底數應以e為底數，而不是10
【1】
基于黏彈性本構模型的熱固性樹脂基復合材料固化變形數值仿真模型.pdf
最后，歡迎大家關注“320科技工作室”微信公眾號，有相關需求可以添加管理員聯系方式~
展開 
基于Abaqus的光固化3D打印結構仿真
數字化光處理技術(DLP)是一種面成型的光固化打印方式，也是一種廣泛使用的3D打印技術。DLP的技術原理是在光源作用下使液態樹脂發生聚合反應固化成型。掃描完一層后，平臺上升或下降一個切片層厚度，樹脂補充完成后進行下一層的固化，新固化層與上一層緊密結合在一起，如此層層疊加即可完成三維結構的構建。
我們在之前的文章里介紹了復合材料固化變形的相關內容。樹脂固化過程中會出現化學-熱-變形多場之間的相互耦合，固化產生的化學收縮以及熱應變會導致復合材料結構內部產生較大的內應力，并導致結構形狀發生改變。光固化3D打印結構變形與復合材料固化變形本質上是類似的，都是由樹脂的固化收縮和熱應變導致內部產生殘余應力，釋放邊界約束后結構發生回彈變形。與復合材料固化變形相比，光固化需要額外考慮光照對固化速率的影響，一定程度上增加了分析的復雜性。
本文將從固化動力學模型、材料本構以及建模方法三方面展開介紹。
固化動力學模型
光固化過程中，樹脂的固化速率與溫度和光照強度相關。本文采用同時考慮溫度和光照強度影響的自催化固化模型：
式中k為反應速率常數，
固化過程會同時釋放出大量熱量，熱量與固化速率的關系如下
用到的子程序：sdvini(設置初始固化度),usdfld(更新固化度和固化速率)，hetval(內部生熱)
材料本構
我們在之前的文章中介紹了復合材料固化變形的幾種本構，包括線彈性本構、路徑依賴本構以及粘彈性本構。這里我們采用指數形式的本構來描述材料剛度與固化度之間的關系。
展開 基于Abaqus的復合材料固化仿真模擬
利用abaqus進行結構力學仿真已經十分普遍，但有關預測復合材料固化過程中變形及殘余應力的內容相對較少。復合材料的固化行為可分解為熱傳導、固化交聯反應、樹脂流動-壓實、固化變形四個重要環節，涉及Hetval、Uexpan及Umat等子程序，內容繁復且不易理解，下面將簡要介紹各個環節及所使用的子程序。
溫度場研究一直是復合材料構件制造中的一個主要研究熱點，溫度場通過影響樹脂的固化度和固化速率對材料熱物理性能、化學收縮、殘余應力等產生影響，因此針對熱-化學耦合的仿真研究比較多。此過程中主要使用的子程序有Hetval、Film及USDFLD等，在莊茁先生的著作中有部分源代碼，在此不再贅述。
Film子程序簡介
該子程序在熱交換分析中用來定義非均勻的對流換熱系數和環境溫度，可以定義的變量有H（1）、H（2）及sink。
H（1）用于定義節點上的對流換熱系數，如果沒有定義，那么將被初始化為0。
H（2）用于定義節點上對流換熱系數相對于表面溫度的變化率，通過定義這個值，可以提高非線性分析中的收斂速度。
sink用于定義環節溫度，如果沒有定義，那么也將被初始化為0。
Hetval子程序簡介
該子程序用來在傳熱分析中定義內部樹脂固化生熱，可以定義和更新的變量有FLUX（1）、FLUX（2）及STATEV。
FLUX（1）用于定義當前節點上的熱流密度。
FLUX（2）用于定義單位溫度變化導致的熱流密度的改變速度。
展開 abaqus纖維復合材料V型件熱固化模型 ￥300
image_process=/format,webp" data-initial-src="https://img.jishulink.com/202504/attachment/3abd0e07f98244bba27d1fa72f4d8fc6.jpg" data-extentions-extra-ocr-id="09c40cd3c8cb6996aa550e0087e7d638">
</figure>
</figure><p>abaqus纖維復合材料V型件熱固化模型，chile子程序，內附CAE，inp及ODB文件，chile模型，操作步驟詳解</p>
展開 abaqus纖維復合材料平板件熱固化仿真，chile模型! ￥200
abaqus纖維復合材料平板件熱固化仿真，chile模型!
內附模型，子程序，學習資料，ODB文件!
Abaqus基于粘彈性本構的復合材料固化成型仿真
復合材料制件成型過程中，由于材料自身的各向異性、樹脂基體的化學收縮反應以及模具作用等因素的影響，導致制件成型過程中產生殘余應力，引起固化變形，從而增加制造成本和裝配難度。因此，合理預測制件固化過程中殘余應力的發展，計算制件的固化變形量，成為降低制造成本、提高生產效率的重要手段。
復合材料固化成型仿真主要包括三個部分：熱-化學模型，固化動力學方程和固化本構。http://www.yqgqt.org.cn/content/post/1261705中介紹了固化成型過程中的熱化學模型和固化動力學方程。為了進一步研究復合材料的固化變形過程，本文又引入了粘彈性本構模型，采用完全熱力耦合的分析方法，預測了復合材料的固化變形。
目前常用的固化本構模型包括：線彈性模型，路徑依賴模型和粘彈性本構模型。
Zocher等提出的粘彈性本構模型其本構關系和應力增量方程為：
其中
式中St_im是歷史狀態變量
其中，增量步內的折算時間
式中，Cu_ij和Cf_ij分別為完全松弛剛度和未松弛剛度；aT、Wm和τm分別為轉換因子、權重系數和松弛時間。松弛時間和權重因子如下
通過Umat子程序編寫粘彈性本構模型，結合Hetval、Disp等子程序進行固化成型過程分析。有限元模型如下圖所示，包括復合材料及模具。在回彈分析時，通過Model Change 移除模具。
固化過程中的溫度和固化度關系的關系如圖所示
計算得到的溫度和應力的關系如圖所示
固化過程中的應力場如下圖所示
移除模具后，可以得到復合材料的回彈變形如圖所示
有相關需求歡迎通過微信公眾號聯系我們。
展開 abaqus碳纖維復合材料熱固化模擬-球形件模型 ￥400
abaqus碳纖維復合材料熱固化模擬，球形件模型，chile模型，內附inp，CAE，ODB模型
Abaqus基于粘彈性本構的復合材料固化成型仿真
Abaqus基于粘彈性本構的復合材料固化變形分析
復合材料制件成型過程中，由于材料自身的各向異性、樹脂基體的化學收縮反應以及模具作用等因素的影響，導致制件成型過程中產生殘余應力，引起固化變形，從而增加制造成本和裝配難度。因此，合理預測制件固化過程中殘余應力的發展，計算制件的固化變形量，成為降低制造成本、提高生產效率的重要手段。
復合材料固化成型仿真主要包括三個部分：熱-化學模型，固化動力學方程和固化本構。http://www.yqgqt.org.cn/content/post/1261705中介紹了固化成型過程中的熱化學模型和固化動力學方程。為了進一步研究復合材料的固化變形過程，本文又引入了粘彈性本構模型，采用完全熱力耦合的分析方法，預測了復合材料的固化變形。
目前常用的固化本構模型包括：線彈性模型，路徑依賴模型和粘彈性本構模型。
Zocher等提出的粘彈性本構模型其本構關系和應力增量方程為：
其中
式中St_im是歷史狀態變量
其中，增量步內的折算時間
式中，Cu_ij和Cf_ij分別為完全松弛剛度和未松弛剛度；aT、Wm和τm分別為轉換因子、權重系數和松弛時間。松弛時間和權重因子如下
通過Umat子程序編寫粘彈性本構模型，結合Hetval、Disp等子程序進行固化成型過程分析。有限元模型如下圖所示，包括復合材料及模具。在回彈分析時，通過Model Change 移除模具。
固化過程中的溫度和固化度關系的關系如圖所示
計算得到的溫度和應力的關系如圖所示
固化過程中的應力場如下圖所示
移除模具后，可以得到復合材料的回彈變形如圖所示
有關于子程序二次開發或者復材仿真的問題可以聯系QQ1653004885或者關注CAE320公眾號
展開 
批量提取Abaqus的節點坐標(初始坐標、指定Step下的變形量、變形后節點坐標) ￥40
<h2>摘要</h2><p>本文介紹如何使用Python腳本二次開發來批量提取ABAQUS輸出數據庫（ODB）文件中指定Step下的Set節點集變形量。通過詳細的步驟說明、代碼示例和圖片展示，您將學會如何使用該腳本，自動化輸出CSV文件包含（Node Label；Step Name、Increment、Step Time，U1，U2）。</p><p>如果還需要按Increment提取每個增量下的變形后的節點坐標的話，在提取變形量的基礎上，與初始坐標進行簡單的計算就可以求得坐標。 （備注：該代碼只提取了x，y方向的變形量）</p><h2>1. 問題描述</h2><p>在工程仿真和分析領域，提取ABAQUS輸出數據庫（ODB）文件中的節點集變形量是一項常見任務。然而，手動提取這些數據是一項繁瑣且容易出錯的工作。因此，需要一種自動化的方法來批量提取指定步驟下按節點集組織的變形量數據。</p><h2>2. 實例展示</h2><p>假設我們有一個名為`example.odb`的ODB文件，其中包含名為`Step-x`的步驟和名為`Set-x`的節點集。運行以上代碼后，腳本會自動將該步驟下節點集的變形量提取出來，并保存為`NodalDisplacement.csv`文件。
展開 Abaqus+PyQt+Python平面變形歐拉角計算
下面以簡單例子介紹平面變形、指向歐拉角的計算，包括絕對歐拉角、相對歐拉角。
1 簡化模型
下面的六面體為表面殼模型，下面由三段梁支持，三段梁分別沿X、Y、Z軸向。六個面的厚度不同，在上側3個面施加不同的壓力，如下左圖所示。位移云圖如下右圖所示。
2 計算要求
計算六面體上面3個面的變形歐拉角，包括3個面的絕對歐拉角，平面2、3相對與平面1的相對歐拉角。平面1、2、3如下圖所示。
3 數據處理
使用平面節點坐標、位移數據計算平面變形歐拉角。可以使用Python腳本輸出平面節點編號、節點坐標（X、Y、Z）、節點位移（U1、U2、U3），如下圖所示。下圖為平面1的10個工況的數據文件，打開的文本文件中7列數據為節點編號、坐標、位移。
三個平面10個工況的節點數據文件如下圖所示。每個文件中包含一個工況一個平面的節點編號、坐標、位移數據。
4 絕對歐拉角計算
使用PyQt+Python開發了一個簡單的小軟件，計算絕對歐拉角、相對歐拉角。
首先計算各平面的絕對歐拉角。
計算平面1的10個工況的絕對歐拉角。
平面1變形的絕對歐拉角計算結果如下圖所示。
伴隨絕對歐拉角計算結果，軟件同時寫出了平面變形前后的坐標系數據，如下圖。每行18個數據，每3個數據為一個坐標軸向量，變形前后2個坐標系，6個坐標軸，18個數據。
5 相對歐拉角計算
利用計算絕對歐拉角時得到的坐標系文件，計算平面變形相對歐拉角，如下圖所示，計算平面2相當于平面1、平面3相對與平面1的相對歐拉角。
計算結果如下圖所示。
6 小結
上述軟件用的算法申請了發明專利，軟件申請了軟著。
展開 Abaqus模擬橡膠大變形
Abaqus為用戶提供了多種本構關系來模擬超彈性材料，這種材料具有高度非線性，當Abaqus進行模擬時假設這種材料是具有彈性、各向同性，并且同時考慮幾何非線性效應。與材料的剪切柔度相比，對于大多數類似橡膠的固體材料，其可壓縮性非常小，當分析對象為平面應力問題、殼、薄膜、梁、桁架、或者鋼筋等，這個問題不值得關注。但是對于固體、平面應變或者軸對稱問題卻不能忽略。對此，Abaqus/Standard提供了雜交單元來模擬超彈性材料中完全的不可壓縮行為。
橡膠材料力學性能的描述方法主要為兩類：一類是認為橡膠為連續介質的現象學描述；另一類是基于熱力學統計的方法。基于連續介質力學的本構模型主要有Polynomial、Reduce Polynomial、Ogden模型等，其中Mooney-Rivlin模型是 Polynomial的特殊形式，Neo-Hookean 模型是Reduce Polynomial的特殊形式。基于熱力學統計主要有Arruda-Boyce和Van der Waals等本構模型。本文利用Abaqus模擬大變形的橡膠，具體步驟如下。
1、在Abaqus/CAE Sketch模塊中作出模型草圖，如圖1所示，然后在Part模塊中分別建立Push、Rubber、Base三個部件。其中Push為解析剛體，Base為離散剛體。
圖1 草圖
2、在Property模塊中定義橡膠的屬性，采用Mooney-Rivlin模型，參數如圖2所示，然后賦給Rubber部件。
圖2 橡膠參數設置
3、裝配，定義分析步，采用默認的場輸出和歷史輸出。為了保證剛開始能夠較容易收斂，設置分析步初始增量步為0.01，打開幾何非線性。
圖3 分析步定義
4、定義接觸對：Push下表面和橡膠表面，Base上表面和橡膠表面。
展開 ABAQUS 單向拉伸大變形模擬
靜態模擬一種軟材料POE的單向拉伸，拉伸應變希望到300%，但是總是在100%就失敗了。不知道哪里出了問題，有沒有高手幫幫忙。
