狀態(tài)方程

在某些情況下，為了準確模擬材料行為，需要EOS。EOS通過將壓力計算為密度的函數(shù)，也許還有能量和/或溫度，來確定材料的靜水壓或體積行為。需要EOS時，通常是具有這些特點：應(yīng)變率非常高，材料壓力遠遠超過屈服應(yīng)力，以及沖擊波的傳播。當然，這些現(xiàn)象通常是相互關(guān)聯(lián)的。

*EOS_LINEAR_POLYNOMIAL或*EOS_GRUNEISEN可能是非氣態(tài)材料最常用的EOS狀態(tài)方程。Gruneisen參數(shù)可用于包括金屬在內(nèi)的許多材料。

總應(yīng)力是偏應(yīng)力和壓力的總和。平均應(yīng)力（sig1+sig2+sig3）/3等于壓力。不采用EOS的本構(gòu)模型直接計算總應(yīng)力。在這些模型中，總應(yīng)力的壓力分量僅基于體積應(yīng)變。例如，對于彈性材料，p＝K*mu，其中K是體積模量，mu = rho/rho0 - 1。

需要附帶EOS的材料模型只計算應(yīng)力的偏分量，即強度行為，而EOS計算總應(yīng)力的壓力分量，即靜水壓力行為。

請注意，EOS只能用于連續(xù)介質(zhì)單元（殼類型為13、14或15的*ELEMENT_SHELL或*ELEMENT_SOLID）。這些單元使用在960版本用戶手冊第II卷第6頁開始的表格中為EOS指定的*MAT材料之一。

如果您使用的材料模型需要EOS，則可以通過使用*EOS_LINEAR_POLYNOMIAL并將C1設(shè)置為體積模量和彈性模量，將所有其他C項設(shè)置為0來實現(xiàn)簡單的體積行為。只有當應(yīng)變率為低到中等時，才會推薦這種方法。車輛碰撞中的應(yīng)變率屬于中等。

Zukas（1990，John Wiley and Sons）編輯的《高速沖擊動力學》一書是關(guān)于高應(yīng)變率下材料行為的一本很好的參考書。

大約50種材料的EOS參數(shù)在“Equation of State and Strength Properites of Selected Materials”書中給出，Danial J.Steinberg，Lawrence Livermore國家實驗室，1991年（1996年發(fā)布的變更1），UCRL-MA-106439。（LLNL不提供副本。）

關(guān)于*EOS_TABULATED_COMPACTION和*EOS_TABULATED：手冊不是很具體。說明如下：

• eVi項（曲線的橫坐標）表示ln（相對體積），因此在壓縮時為負。

• eVi=ln（相對體積）值應(yīng)按降序給出，即首先是拉伸（正）值，最后是壓縮（負）值。

• 壓縮時壓力為正。如果gamma=0，Ci等于加載曲線上的壓力。因此Ci應(yīng)該具有與eVi相反的代數(shù)符號。

  
  

額外的歷史變量

當讀入d3plot數(shù)據(jù)庫時，額外的歷史變量（參見*DATABASE_EXTENT_BINARY中的NEIPH和NEIPS）可以通過LS-POST中選擇?Fcomp>Misc>history var#n，繪制顯示出該變量隨時間的變化曲線。通過選擇 History > Scalar 來創(chuàng)建額外歷史變量的時間歷程圖。

在用戶材料子程序應(yīng)用中，采用額外的歷史變量優(yōu)勢明顯，因為用戶可以知道額外的歷史變量是什么。在LS-DYNA的970以前版本，umat子程序中的第一個歷史變量不會作為歷史變量#1存儲在d3plot數(shù)據(jù)庫中。存儲位置取決于許多因素，例如子程序是矢量化的還是非矢量化的，單元是殼還是實體類型等。

//?來自Lee Bindeman關(guān)于該主題的更多信息

“當使用矢量化子程序（即umat46v而非umat46）和一個用于3D實體單元的材料時，無論用戶定義的材料是否正交各向異性，都有6個歷史變量自動用于轉(zhuǎn)換矩陣的6項中。當材料為正交各向異性（IORTHO=1）時，這6個變量是自動分配的；然而當材料不是正交各向異性（IORTHO=0）時，這些變量不是自動分配的，必須由用戶定義的材料輸入來分配。因此，如果材料使用46個歷史變量，則需要設(shè)置NHV=52。為了將46個歷史變量寫入d3plot文件，需要通過在*DATABASE_EXTENT_BINARY上設(shè)置NEIPH=52來請求52個額外的歷史變量。

后處理時，歷史變量1到6將包含轉(zhuǎn)換矩陣項。如果材料是各向同性的，這些都將為零。歷史變量7到52將包含子程序中的歷史變量1到46。

對于二維材料（對于殼單元），上述規(guī)則會發(fā)生更改。在這種情況下，只存儲了2個轉(zhuǎn)換項，因此只需要分配和請求2個額外的歷史變量。

因為上述這種混亂的情況，在版本970修訂版2903及更高版本后做了修正。

有了這個修正，就不再需要為NHV分配額外的歷史變量，也不再需要在d3plot文件中請求6（或2）個額外歷史變量。您只需要分配相同數(shù)量的要使用的歷史變量，并請求寫入d3plot文件的數(shù)量。如果材料是各向同性的（IORTHO=0），則將省略轉(zhuǎn)換項，然后用戶子程序中的歷史變量編號將與d3plot文件中的編號相匹配。然而，如果材料是正交各向異性的（IORTHO=1），那么6個（或2個）變換項將被寫入d3plot文件，因此用戶子程序和d3plot中的歷史變量編號將不匹配。”

力的輸出

//?可以通過以下幾種方式獲得受力情況

• 使用*DATABASE_CROSS_SECTION?和?*DATABASE_SECFORC可以獲得一個橫截面上的內(nèi)力和內(nèi)力矩。注意，在使用set選項設(shè)置橫截面時，必須提供用于定義橫截面的節(jié)點集以及橫截面某一側(cè)的至少一個單元集。

• 使用*DATABASE_NODAL_FORCE_GROUP?和?*DATABASE_NODFOR可以獲得施加在節(jié)點或節(jié)點集上的外部力

• 使用*DATABASE_BNDOUT可以獲得邊界反力

• 使用*DATABASE_SPCFORC可以獲得由*BOUNDARY_SPC定義的邊界上的反力

• *DATABASE_RCFORC可以輸出接觸力合力。在單面接觸中，還需要使用*CONTACT_FORCE_TRANSDUCER_(option)來提取接觸力合力

• 單個節(jié)點所受的力由*DATABASE_NCFORC輸出。必須在*CONTACT中設(shè)置SPR, MPR中至少一個輸出控制選項為1，以便向程序指定輸出側(cè)為接觸面主面?zhèn)然驈拿鎮(zhèn)取?/span>

以上數(shù)據(jù)分別由程序輸出到ASCII文件SECFORC, NODFORC, BNDOUT, SPCFORC,?RCFORC和NCFORC中。

節(jié)點力

下表總結(jié)了各種ASCII輸出文件是否包含由各種載荷或邊界條件引起的節(jié)點力。

NCFORC和RCFORC分別包含節(jié)點接觸力和接觸合力。只有當設(shè)置了*CONTACT（Card 2）中的打印標志（SPR、MPR）時，才會包含在NCFORC中。

相關(guān)閱讀

LS-DYNA電池行業(yè)應(yīng)用介紹（三）

LS-DYNA電池行業(yè)應(yīng)用介紹（二）

LS-DYNA電池行業(yè)應(yīng)用介紹（一）

LS-DYNA | 如何將隱式計算轉(zhuǎn)為顯式計算？

LS-DYNA | 如何查找和消除初始穿透？

摩爾芯創(chuàng)專注于為硅基光電子、電力電子、高科技半導(dǎo)體等行業(yè)提供先進的電子設(shè)計自動化（EDA）和計算機輔助工程（CAE）協(xié)同解決方案；提供從光學、光電子學、電磁場、結(jié)構(gòu)、流體、多物理場耦合等全面的工業(yè)軟件應(yīng)用解決方案和咨詢服務(wù)。

官網(wǎng)：rfteuxon.cn

電話：15521163312(微信同號)

郵箱：wenye@mooreda.com.cn