abaqus 顯卡 計(jì)算的案例
顯卡隨機(jī)振動(dòng)疲勞仿真計(jì)算
本文將以顯卡模型為例，闡述如何使用ANSYS Mechanical聯(lián)合ANSYS nCode DesignLife進(jìn)行隨機(jī)振動(dòng)疲勞仿真。
大咖慧網(wǎng)絡(luò)培訓(xùn)
2022年5月24日-26日，安世亞太大咖慧推出電子行業(yè)疲勞壽命專題線上培訓(xùn)，專題講座包含：隨機(jī)振動(dòng)載荷下支撐構(gòu)件疲勞壽命評(píng)估、PCB電路板中的焊點(diǎn)可靠性分析、PCB電路板疲勞壽命分析內(nèi)容，不容錯(cuò)過。
報(bào)名方式
案例背景描述
計(jì)算幾何為簡(jiǎn)化的顯卡模型，見下圖。PCB板、金屬支架材料為結(jié)構(gòu)鋼，其余構(gòu)件簡(jiǎn)化為鋁合金。金屬支架左側(cè)3端面固定支撐，隨機(jī)振動(dòng)載荷類型為G加速度譜，方向?yàn)閅向，具體數(shù)值見圖，計(jì)算該工況下的疲勞壽命。
圖 1模型
圖2 G加速度譜
1、仿真流程搭建
為提升計(jì)算效率，本例采用MSUP諧響應(yīng)分析聯(lián)合nCode進(jìn)行隨機(jī)振動(dòng)疲勞仿真。具體模塊搭建如下：
圖 3仿真流程
注：使用該方法進(jìn)行隨機(jī)振動(dòng)疲勞計(jì)算時(shí)，需先將nCode模塊拖拽至“modal”模塊處，然后再將“Harmonic Response”的solution與“nCode”的solution相連，完成流程搭建。若直接將“Harmonic Response”與“nCode”相連，在后續(xù)提交計(jì)算時(shí)，軟件會(huì)提示沒有材料數(shù)據(jù)，無法進(jìn)行求解計(jì)算。
2、仿真關(guān)鍵步驟
1) 模態(tài)計(jì)算
模態(tài)計(jì)算的結(jié)果-固有頻率為結(jié)構(gòu)的固有屬性，所以在仿真中只需設(shè)置結(jié)構(gòu)的約束條件和提取的階數(shù)即可。
展開 技術(shù)干貨丨如何選用顯卡 (GPU) 開展EDEM的計(jì)算？
但對(duì)于消費(fèi)級(jí)GPU顯卡，64位雙精度運(yùn)算由于使用場(chǎng)景較少 (大多數(shù)為一般通用計(jì)算任務(wù))，廠商為了降低售價(jià)提高銷量，削弱了消費(fèi)級(jí)GPU的FP64能力 (但Nvidia RTX 30系以上消費(fèi)級(jí)顯卡基本保留部分FP64能力)，導(dǎo)致這些顯卡雙精度計(jì)算能力受限。</p><p><br></p><p>對(duì)于這些顯卡，在使用時(shí)需要考慮更換精度模式保證計(jì)算速度。對(duì)于Nvidia頂尖顯卡，如A100 (800)、H100 (800) 等，這些顯卡專供HPC科學(xué)計(jì)算，F(xiàn)P64能力不會(huì)削減，因此直接采用雙精度模式計(jì)算即可。而對(duì)于桌面級(jí)顯卡例如 Nvidia RTX 30/40 系顯卡 (如4090等)，雙精度模式受限，需適當(dāng)降低精度保證計(jì)算效率。尤其是在試算或者粗略計(jì)算大概結(jié)果的場(chǎng)景，快速計(jì)算十分必要。此時(shí)，降低精度快速出結(jié)果，對(duì)于業(yè)務(wù)而言至關(guān)重要。</p><p><br></p><p>混合精度模式 (hybrid precision) 為Nvidia特色技術(shù)，可根據(jù)GPU狀況，智能對(duì)浮點(diǎn)數(shù)據(jù)指定單雙精度策略。可以在效率和精度之間取得最佳平衡。</p><p><br></p><p><strong>多GPU使用策略</strong></p><p><br></p><p>最新的EDEM版本已支持調(diào)度多個(gè)顯卡 (GPU) 進(jìn)行并行計(jì)算。依舊以小學(xué)作比喻，這就仿佛將一個(gè)口算任務(wù)分配給兩個(gè)小學(xué)完成。</p><p><br></p><p>一般而言，多GPU僅在問題規(guī)模特別大時(shí)使用。因?yàn)?em>計(jì)算過程中，GPU之間會(huì)通信，通信成本較大。因此，盡可能使用單個(gè)GPU進(jìn)行計(jì)算。</p><p><br></p><p>但是如果問題規(guī)模特別大，例如數(shù)百萬顆粒的時(shí)候，單GPU可能面臨存儲(chǔ)空間不足的情況 (離散單元法不僅計(jì)算量大，存儲(chǔ)空間占用也很大)。
展開 ABAQUS混凝土損傷塑性模型-C30EXCEL計(jì)算表格（含計(jì)算結(jié)果）
這是根據(jù)GB50010-2010中混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范中的混凝土本構(gòu)模型，結(jié)合文獻(xiàn)所述的損傷因子定義，編制的計(jì)算C30混凝土非彈性應(yīng)變和損傷因子的EXCEL表格。也是邊學(xué)變做，希望能和大家多交流。
C30砼本構(gòu)（損傷塑性模型）.rar
ABAQUS混凝土損傷塑性損傷因子計(jì)算依據(jù).rar
ABAQUS 在計(jì)算到step3后計(jì)算中斷，變形過大
我在原本可以計(jì)算的模型的基礎(chǔ)上修改的，只是刪了幾個(gè)樓板與梁連接的栓釘
Mesh Free-眼鏡剛強(qiáng)度校核計(jì)算，附Abaqus計(jì)算結(jié)果對(duì)比
采用Mesh Free對(duì)某品牌眼鏡整體剛度、強(qiáng)度進(jìn)行校核，如下圖所示，Mesh Free支持在不用做幾何清理的前提下進(jìn)行計(jì)算分析，導(dǎo)入模型部件可以包含細(xì)節(jié)特征，比如螺釘上的倒角。
眼鏡定義了5種線彈性材料：鏡架主體采用鈦合金；眼鏡片采用樹脂；螺釘?shù)冗B接件采用鋼；鼻托和鏡片扎線采用兩種不同的尼龍材料。
眼鏡腿一只固定，另一只向上掰，加力1N，模擬分析此種工況下眼鏡整體結(jié)構(gòu)的剛度、強(qiáng)度。
Mesh Free所有接觸面定義為完全剛性連接，Abaqus作同樣處理，不考慮非線性因素，對(duì)比二者的線性計(jì)算結(jié)果。
Mesh Free給出的眼鏡最大變形為23.92mm，Abaqus的結(jié)果為23.46mm。
Mesh Free給出的眼鏡最大應(yīng)力為303.4MPa，Abaqus的結(jié)果為308.3MPa。
談?wù)凪esh Free使用感受：
雖然我常用ABQ，但是不得不說，對(duì)于包含細(xì)節(jié)幾何特征的復(fù)雜裝配結(jié)構(gòu)建模分析，Mesh Free真的要比Abaqus高效的多。
據(jù)我了解Mesh Free的非線性也在大力的開發(fā)之中，目前已經(jīng)支持經(jīng)典塑性材料非線性、邊界條件非線性也可以設(shè)置滑動(dòng)和一般的摩擦接觸。
對(duì)不熟悉常規(guī)有限元操作的結(jié)構(gòu)設(shè)計(jì)人員來說，不用幾何清理、不用劃網(wǎng)格是極好的體驗(yàn)。
關(guān)鍵是Mesh Free的結(jié)果也確實(shí)很準(zhǔn)，目前的CAE無非是追求更準(zhǔn)的基礎(chǔ)上算的更快，這兩點(diǎn)Mesh Free無疑是滿足的。
Mesh Free
Abaqus
展開 ABAQUS混凝土損傷塑性模型-C30EXCEL計(jì)算表格（含計(jì)算結(jié)果） ￥3.7
計(jì)算表格（如下）中標(biāo)黃部分的參數(shù)可自行設(shè)定后，EXCEL程序會(huì)自動(dòng)計(jì)算“抗拉強(qiáng)度、非彈性應(yīng)變、受拉損傷因子”。
隨機(jī)振動(dòng)分析-abaqus（附一個(gè)電池包計(jì)算案例） ￥20
四、如何將時(shí)域隨機(jī)振動(dòng)曲線轉(zhuǎn)換得到功率譜密度曲線
五、 隨機(jī)振動(dòng)分析理論
附.常見功率譜密度曲線給出形式
附.以dB/oct形式給出的功率譜密度曲線如何計(jì)算
附.國(guó)標(biāo)中定義的PSD譜總均方根加速度值是如何計(jì)算的？
六. 隨機(jī)振動(dòng)分析案例-abaqus
第一步：計(jì)算結(jié)構(gòu)模態(tài)，輸出位移和應(yīng)力。
第二步：隨機(jī)振動(dòng)分析
2.1 定義輸出頻率上下限和模態(tài)阻尼
2.2 定義PSD載荷及加載
2.3 定義輸出
2.4 隨機(jī)振動(dòng)計(jì)算頭文件設(shè)置
2.5 隨機(jī)振動(dòng)分析結(jié)果
2.6 隨機(jī)振動(dòng)σ應(yīng)力結(jié)果評(píng)價(jià)
Abaqus接觸非線性在有限元計(jì)算分析中的應(yīng)用 附莊茁ABAQUS非線性有限元分析與實(shí)例下載
來源：有限元在線
ABAQUS的非線性主要在有三種：幾何非線性，材料非線性以及接觸非線性。接觸非線性在ABAQUS的有限元計(jì)算分析中應(yīng)用非常廣泛，特別是動(dòng)態(tài)顯式的求解，只要模型中包含兩個(gè)以上相互接觸的部件，就要用到接觸非線性。
ABAQUS接觸非線性的設(shè)置主要在Interation模塊中完成，設(shè)置接觸的屬性時(shí)，可以設(shè)置摩擦系數(shù)，阻尼系數(shù)，損壞，失效準(zhǔn)則等非線性參數(shù)，如圖1所示。
如圖2所示，在接觸定義界面，可以選擇通用接觸、面-面接觸、自接觸等各種非線性接觸方式。
在接觸編輯界面，可以選擇機(jī)械約束方式為運(yùn)動(dòng)學(xué)接觸算法，或是懲罰接觸方式，還可選擇滑移方式為有限滑移或小滑移，如圖3所示。
這是對(duì)模型定義非線性接觸后得到的分析結(jié)果，以供參考。
下載地址：莊茁ABAQUS非線性有限元分析與實(shí)例
展開 有限元理論基礎(chǔ)及Abaqus內(nèi)部實(shí)現(xiàn)方式研究系列27： Abaqus內(nèi)部計(jì)算和顯示的應(yīng)變
（1）顯示應(yīng)變：Abaqus計(jì)算完畢后得到導(dǎo)入結(jié)果，在后處理中查看，應(yīng)變E11=8.528e-1,E22=-5.173e-1如下：
（2）計(jì)算應(yīng)變：Abaqus中采用UMAT子程序，利用我們的子程序調(diào)試插件DUS調(diào)試UMAT，在Visual Studio中查看dStran的值，發(fā)現(xiàn)在計(jì)算完應(yīng)變后，進(jìn)入U(xiǎn)MAT時(shí)，E11=8.528e-1,E22=-5.173e-1，調(diào)試如下：
可以發(fā)現(xiàn)殼單元Abaqus的計(jì)算應(yīng)變和顯示應(yīng)變一樣，猜測(cè)都是對(duì)數(shù)應(yīng)變。
1.5.3 iSolver的應(yīng)變
iSolver中采用自帶材料進(jìn)行計(jì)算，材料參數(shù)和UMAT的輸入完全一致。
為了計(jì)算和Abaqus完全一致，iSolver也采用對(duì)數(shù)應(yīng)變計(jì)算方式，得到的應(yīng)變顯示如下，可發(fā)現(xiàn)和Abaqus完全一致。
==總結(jié)==
由上可以看到，在實(shí)際計(jì)算中，對(duì)體單元，Abaqus和iSolver都采用變形率積分方式來計(jì)算應(yīng)變，對(duì)殼單元，Abaqus和iSolver都采用對(duì)數(shù)應(yīng)變。一般理論書都認(rèn)為Abaqus是因?yàn)閷?duì)數(shù)應(yīng)變計(jì)算復(fù)雜才采用別的應(yīng)變，但個(gè)人認(rèn)為應(yīng)該不是這個(gè)原因，因?yàn)?em>Abaqus對(duì)體單元為了顯示對(duì)數(shù)應(yīng)變，依然重新計(jì)算了一遍，說明Abaqus體單元采用變形率是有其它原因的，具體什么原因我也沒研究清楚，歡迎探討。
如果有任何其它疑問或者項(xiàng)目合作意向，也歡迎聯(lián)系我們：
snowwave02 From www.yqgqt.org.cn
email: snowwave02@qq.com
以往的系列文章：
1.7.1 ========第一階段========
第一篇：S4殼單元?jiǎng)偠染仃囇芯俊?http://www.yqgqt.org.cn/content/post/338859
第二篇：S4殼單元質(zhì)量矩陣研究。
展開 Abaqus+PyQt+Python平面變形歐拉角計(jì)算
下面以簡(jiǎn)單例子介紹平面變形、指向歐拉角的計(jì)算，包括絕對(duì)歐拉角、相對(duì)歐拉角。
1 簡(jiǎn)化模型
下面的六面體為表面殼模型，下面由三段梁支持，三段梁分別沿X、Y、Z軸向。六個(gè)面的厚度不同，在上側(cè)3個(gè)面施加不同的壓力，如下左圖所示。位移云圖如下右圖所示。
2 計(jì)算要求
計(jì)算六面體上面3個(gè)面的變形歐拉角，包括3個(gè)面的絕對(duì)歐拉角，平面2、3相對(duì)與平面1的相對(duì)歐拉角。平面1、2、3如下圖所示。
3 數(shù)據(jù)處理
使用平面節(jié)點(diǎn)坐標(biāo)、位移數(shù)據(jù)計(jì)算平面變形歐拉角。可以使用Python腳本輸出平面節(jié)點(diǎn)編號(hào)、節(jié)點(diǎn)坐標(biāo)（X、Y、Z）、節(jié)點(diǎn)位移（U1、U2、U3），如下圖所示。下圖為平面1的10個(gè)工況的數(shù)據(jù)文件，打開的文本文件中7列數(shù)據(jù)為節(jié)點(diǎn)編號(hào)、坐標(biāo)、位移。
三個(gè)平面10個(gè)工況的節(jié)點(diǎn)數(shù)據(jù)文件如下圖所示。每個(gè)文件中包含一個(gè)工況一個(gè)平面的節(jié)點(diǎn)編號(hào)、坐標(biāo)、位移數(shù)據(jù)。
4 絕對(duì)歐拉角計(jì)算
使用PyQt+Python開發(fā)了一個(gè)簡(jiǎn)單的小軟件，計(jì)算絕對(duì)歐拉角、相對(duì)歐拉角。
首先計(jì)算各平面的絕對(duì)歐拉角。
計(jì)算平面1的10個(gè)工況的絕對(duì)歐拉角。
平面1變形的絕對(duì)歐拉角計(jì)算結(jié)果如下圖所示。
伴隨絕對(duì)歐拉角計(jì)算結(jié)果，軟件同時(shí)寫出了平面變形前后的坐標(biāo)系數(shù)據(jù)，如下圖。每行18個(gè)數(shù)據(jù)，每3個(gè)數(shù)據(jù)為一個(gè)坐標(biāo)軸向量，變形前后2個(gè)坐標(biāo)系，6個(gè)坐標(biāo)軸，18個(gè)數(shù)據(jù)。
5 相對(duì)歐拉角計(jì)算
利用計(jì)算絕對(duì)歐拉角時(shí)得到的坐標(biāo)系文件，計(jì)算平面變形相對(duì)歐拉角，如下圖所示，計(jì)算平面2相當(dāng)于平面1、平面3相對(duì)與平面1的相對(duì)歐拉角。
計(jì)算結(jié)果如下圖所示。
6 小結(jié)
上述軟件用的算法申請(qǐng)了發(fā)明專利，軟件申請(qǐng)了軟著。
展開 Abaqus-原來顯式計(jì)算也可以這么快
本文詳細(xì)介紹了 Abaqus/explicit 中的準(zhǔn)靜態(tài)分析及其提高求解速度的方法。
我如何知道我的模擬是否是準(zhǔn)靜態(tài)的？
準(zhǔn)靜態(tài)問題是通常可以使用 Abaqus/Standard 解決的問題之一，但由于觸或材料的復(fù)雜性，可能難以收斂，從而導(dǎo)致大量迭代。具有挑戰(zhàn)性的非線性準(zhǔn)靜態(tài)問題通常涉及：
接觸條件非常復(fù)雜，Abaqus/Standard 可能會(huì)因接觸問題而無法收斂；
非常大的變形可能導(dǎo)致嚴(yán)重的網(wǎng)格變形；
例如，通常在金屬成形分析中，我們面臨這樣的困難：用 Abaqus/Standard 來模擬這樣的問題確實(shí)很困難。
示例：模擬拉深過程中的撕裂
Abaqus/Explicit顯式準(zhǔn)靜態(tài)分析問題
Abaqus/Explicit 對(duì)于高度非線性靜態(tài)（準(zhǔn)靜態(tài)分析）問題的建模更加有效。對(duì)于涉及接觸和金屬成型等非常大變形的三維問題尤其如此。
應(yīng)用 Abaqus/Explicit 模擬準(zhǔn)靜態(tài)事件需要特別考慮。在自然時(shí)間段內(nèi)對(duì)過程進(jìn)行建模在計(jì)算上是不切實(shí)際的。從字面上看，需要數(shù)百萬次的時(shí)間增量。因此，我們?cè)谀M中人為地提高加載過程的速度以獲得經(jīng)濟(jì)的解決方案。
使用 Abaqus/Explicit 獲得經(jīng)濟(jì)的準(zhǔn)靜態(tài)分析解決方案的兩種方法是：
1、人為提高加載速率
我們可以通過提高加載率來人為地減少該過程的時(shí)間尺度。增加的加載速率會(huì)減少模擬的時(shí)間。將加載速率提高 f 倍，分析速度提高 f 倍。
2、采用質(zhì)量縮放
它增加了穩(wěn)定時(shí)間增量的大小，因此完成作業(yè)所需的增量更少。人為地將材料密度（質(zhì)量縮放）增加 f*f 倍，可以將分析速度提高 f 倍。
在本文中，我們的重點(diǎn)是提高加載速率的方法。
展開 
Abaqus批量計(jì)算工具
功能：
每個(gè)inp任務(wù)放到單獨(dú)的文件夾中，不會(huì)導(dǎo)致計(jì)算結(jié)果文件都混在一起
支持 .for 文件的子程序
輸出計(jì)算成功或失敗信息
4個(gè)文件的求解結(jié)果示例
剛開發(fā)完成，沒怎么測(cè)試，可能有bug，歡迎大家反饋。
Abaqus批量計(jì)算.zip
abaqus并行計(jì)算中CPU超限的解決辦法
abaqus并行計(jì)算中CPU超限的解決辦法-Abaqus Error: The number of cpus (16) exceeds the number of cpus available(8）
# 1.說明，
CPU配置為i9-9900K，8核心16線程
# 2. 出現(xiàn)問題：
abaqus job=XX user=XX.for input=XX.inp cpus=16 int
出現(xiàn)報(bào)錯(cuò)：Abaqus Error: The number of cpus (16) exceeds the number of cpus available (8).
# 3. 問題解析
這種情況主要是配置文件關(guān)于并行計(jì)算設(shè)置導(dǎo)致的，尤其是ABAQUS2023，
# 4.解決辦法：
使用everything找到路徑下：X:\SIMULIA\EstProducts\2023\win_b64\SMA\site\custom_v6.env的
custom_v6.env文件，其中X代表盤符，
修改前
在文件最后加一句話：
import socket
mp_host_list=[[socket.gethostname(),16],]
其中16代表自己電腦的線程數(shù)，一定要與自己電腦的線程一樣
修改后
展開 基于ABAQUS的橡膠懸置膠合件剛度仿真計(jì)算
本文將以一個(gè)懸置膠合件仿真的實(shí)例講解一下如何利用ABAQUS來獲取其三個(gè)方向的靜態(tài)特性。所用膠合件的數(shù)模圖如圖1所示。其設(shè)計(jì)圖紙上標(biāo)注的三向剛度如表1所示，膠料硬度是邵氏50±5度。
圖1 膠合件結(jié)構(gòu)
表1設(shè)計(jì)要求
1、 網(wǎng)格劃分
采用HYPERMESH對(duì)圖一懸置進(jìn)行網(wǎng)格劃分到的有限元模型如圖2所示。
2、材料設(shè)置
把劃分好的網(wǎng)格導(dǎo)入ABAQUS中，設(shè)置其材料參數(shù)，由于不同本構(gòu)模型對(duì)橡膠懸置膠合件剛度計(jì)算結(jié)果有一定的影響。結(jié)合何小靜，上官文斌發(fā)表的《橡膠隔振器靜態(tài)力－ 位移關(guān)系計(jì)算方法》一文的研究結(jié)果表明，Mooney-Rivlin 模型的計(jì)算精度最高，其相對(duì)誤差均小于10%，所以本文采用M-R模型進(jìn)行計(jì)算。50度膠料的M-R材料常數(shù)C10=0.2969，C01=0.0584。
3、剛度求解
3.1求解X方向剛度
按表 1要求，做如下設(shè)置：在Z方向先預(yù)載8mm，再在X向加載500N。取值0~5.6mm,對(duì)X向靜剛度進(jìn)行求解。
求得的力和位移關(guān)系見表2所示，用表中數(shù)據(jù)進(jìn)行畫圖差值可得到圖3所示的X向靜剛度為38N/mm,與設(shè)計(jì)值非常接近，其變形云圖見圖4所示
表2 X向力和位移關(guān)系表
圖3 X向剛度差值結(jié)果
圖4 X向云變形圖
3.2求解Y方向剛度
按表 1要求，做如下設(shè)置：在Z方向先預(yù)載8mm，再在X向加載1000N。取值2~4mm,對(duì)Y向靜剛度進(jìn)行求解。
求得的力和位移關(guān)系見表3所示，用表中數(shù)據(jù)進(jìn)行畫圖差值可得如圖5到Y(jié)向靜剛度為98N/mm,與設(shè)計(jì)值80N/mm有一定差異，見圖4。其變形云圖見圖6.
展開 Abaqus中陶瓷本構(gòu)模型及其數(shù)值計(jì)算應(yīng)用
2 數(shù)值計(jì)算軟件中本構(gòu)模型
陶瓷由于其波速高、模量大，具有良好的抗侵徹性能，在各類型裝甲設(shè)計(jì)中被廣泛應(yīng)用。而JH本構(gòu)形式簡(jiǎn)單，易于理解，已成為Abaqus、LS-DYNA和Autodyn等商用軟件的內(nèi)嵌本構(gòu)模型，可一定程度上滿足日常使用及工程計(jì)算要求。
對(duì)于陶瓷材料Abaqus幫助中給出了3種本構(gòu)模型，Extended Drucker-Prager本構(gòu)（以下簡(jiǎn)稱DP本構(gòu)）、JH-2和JHB本構(gòu)模型。DP本構(gòu)多用來模擬巖土材料（粒狀土壤和巖石），擴(kuò)展DP本構(gòu)給出的應(yīng)力與壓力的關(guān)系也與JH本構(gòu)中未損傷時(shí)應(yīng)力與壓力的關(guān)系類似，其損傷段定義采用等效塑性應(yīng)變與應(yīng)力三軸度的對(duì)應(yīng)關(guān)系進(jìn)行定義，狀態(tài)方程采用Mie-Grüneisen形式（詳見Abaqus相應(yīng)部分幫助）。
Abaqus官方幫助中給出的JHB本構(gòu)模型參數(shù)如表1所示。其中標(biāo)紅部分與Abaqus幫助（2021版本）不同，應(yīng)為幫助原文疏漏。
表1 JHB本構(gòu)模型參數(shù)
JHB本構(gòu)模型的應(yīng)力與壓力關(guān)系主要分為完整（Intact）和損傷（Failed）兩部分，表1中下標(biāo)帶有 i 的即為完整部分相應(yīng)參數(shù)，下標(biāo) f 即代表?yè)p傷部分參數(shù)；雖然JHB本構(gòu)模型公式中考慮了脆性材料的相變特性，表1標(biāo)藍(lán)部分參數(shù)應(yīng)為對(duì)應(yīng)的相變參數(shù)，但幫助中全部設(shè)置為0，推知官方幫助中給出的這組參數(shù)不能考慮陶瓷相變的影響。
Abaqus官方幫助中給出的JH-2本構(gòu)模型參數(shù)如表2所示。
表2 JH-2本構(gòu)模型參數(shù)
JH-2本構(gòu)模型以無量綱形式描述了應(yīng)力和壓力的關(guān)系，以Hugoniot極限下的壓力對(duì)壓力變量進(jìn)行了無量綱化。
JHB本構(gòu)模型狀態(tài)變量如表3所示。
展開 