狀態方程的案例
LS-DYNA的狀態方程模型
狀態方程模型
17.1 狀態方程形式1：Linear Polynomial
這個多項式狀態方程，單位初始體積的內能呈線性，E由
（17.1.1）
其中C0，C1，C2，C3，C4，C5和C6是用戶定義常數。
（17.1.2）
V是相對體積，在膨脹單元中，的系數設為零，即：
線性多項式狀態方程可用伽馬定律狀態方程來模擬氣體。這可以通過設置來實現：
和
其中是比熱的比率。壓力則由下式給出：
請注意，E的單位是壓力的單位
17.2 狀態方程形式2：JWL High Explosive
JWL狀態方程將壓力定義為相對體積，V，以及單位初始體積的內能，E，的函數：
（17.2.1）
其中，ω、A、B、R1和R2為用戶定義的輸入參數。這個狀態方程通常用于在涉及金屬加速度的應用中確定烈性炸藥的爆轟產物的壓力。該方程的輸入參數由Dobratz [1981]給出了各種高爆炸材料的輸入參數，該狀態方程與爆炸燃燒（材料模型8）材料模型一起使用，該模型決定了爆炸單元的點火時間。
17.3 狀態公式形式3：Sack “Tuesday” High Explosives
爆轟產物的壓力根據相對體積V和單位初始體積的內能E給出，如[Woodruff 1973]：
（17.3.1）
其中A1、A2、A3、B1、B2為用戶定義的輸入參數
該狀態方程與爆炸燃燒（材料模型8）材料模型一起使用，該模型決定了爆炸單元的點火時間。
展開 LS-DYNA中的操作及設置（二）（有效塑性應變，環境變量，狀態方程）
setenv LSTC_MEMORY auto
setenv LSTC_MEMORY heap
三、狀態方程（Equation of state）
在某些情況下，需要使用狀態方程來精確模擬材料的變形行為。狀態方程可以通過計算材料所受壓力與密度（有時還有能量和溫度）之間的關系來確定材料的變形行為。需要使用狀態方程的情形主要有應變率非常高、材料所受壓力遠高于屈服應力以及沖擊波的傳播等。實際上，這些情況一般都是同時出現的。
對于非氣態材料來說，*EOS_LINEAR_POLYNOMIAL 和 *EOS_GRUNEISEN是最常用的兩種狀態方程。Gruneisen的參數對于包括金屬在內的許多材料都是適用的。
在物體受力時，總應力是偏應力和壓力的總和，平均應力(sig1 + sig2 + sig3)/3等于壓力。對于不考慮狀態方程的本構模型，程序會直接計算主應力，主應力的壓力分量只與體積應變有關。例如，對于彈性材料來說，p = K * mu，其中K為體積模量，mu = rho/rho0 - 1。
對于考慮狀態方程的模型來說，材料本身的本構模型會計算總應力的偏應力分量，而狀態方程則會計算壓力分量。
注意，狀態方程只適用于連續介質單元（*ELEMENT_SHELL with shell type 13, 14, or 15 or *ELEMENT_SOLID），并且材料模型為需要EOS的*MAT_。
如果你在使用需要EOS的本構模型，可以利用*EOS_LINEAR_POLYNOMIAL來實現簡單的體積行為（bulk behavior），此時需要設置C1為體積模量，其他參數均為0。
展開 復雜隱式極限狀態方程可靠性敏度分析新方法及其工程應用
隱式極限狀態方程
復雜隱式極限狀態方程可靠性敏度分析新方法及其工程應用.pdf
ABAQUS空氣材料狀態方程？
空氣材料的狀態方程數值US-UP和單位制有關系嗎？具體怎么換算？
高壓固體狀態方程
YTRAN狀態方程庫不是很全
只有多項式,TAIT,JWL和理想氣體.這些都是用于流體的.對于做固體材料在高壓下的瞬態模擬的兄弟們來說這些是遠遠不夠的.
好在DYTRAN提供了用戶子程序這條路徑.
兄弟作了一個Gruneison狀態方程子程序,供大家參考.
170318-.doc
170317-EOS_Gruneison.rar
LSDYNA的grun狀態方程
狀態方程的解釋，不是原創
想學習更多的知識，請聯系我們！
微信公眾號：名稱：“DR有限元”
號碼：“hello_cae”
炸藥爆轟產物JWL狀態方程
關于炸藥爆轟產物狀態方程，很多文獻都有給出，但不同文獻之間存在明顯區別，尤其是關于內能密度和相對體積，有的采用初始內能密度，有的采用內能密度，有的采用初始相對體積，有的采用相對體積，莫衷一是。
在此拋磚引玉，請各位發表高見！
一些常用的參數和狀態方程，注意單位
材料模型與狀態方程.doc
25種材料狀態仿真、Johnson-cook本構方程、Johnson-cook失效模型參數 ￥49.99
25中金屬材料的狀態方程和Johnson-cook本構和Johnson-cook斷裂失效參數，囊括了鋁，銅，鋼，鈦，鉛，鎢等常見的材料，完整的D1-D5參數，稀缺資源，具有較高的參考價值。
轉，關于ABAQUS在流固耦合方面的應用
針對ABAQUS有限元分析軟件中的流固耦合功能，簡述了其中理想氣體狀態方程的各參數含義以及流固耦合的分析要點。文章通過ABAQUS仿真分析和理論數值解的對比，證明了ABAQUS軟件計算理想氣體狀態方程的可信性，最后介紹其在某型號彈上的分析應用。
隨著有限元技術的發展和用戶要求的提高，各大有限元軟件都含有流固耦合模塊，其主要用于液體、理想氣體和JWL的模擬，本文著重介紹ABAQUS中理想氣體狀態方程的功能和應用。為了驗證ABAQUS理想氣體狀態方程模擬氣體壓縮的正確性，首先利用其模擬簡單的氣體壓縮過程，并獲得該過程中氣體的狀態變化曲線（仿真曲線）；然后通過matlab求解該模型理論上的氣體狀態方程，并在Matlab上獲得數值解（理論曲線）。將仿真曲線和理論曲線進行對比，發現二者非常吻合，證明了ABAQUS模擬理想氣體狀態的可信性。在此基礎上，將其用在某型號彈上的流固耦合分析。
在設置模型過程中，活塞和氣體之間的接觸通過inp文件的關鍵字實現，經過實踐證明，這樣的定義方式可以有效避免氣體的泄露。其定義過程為：
*******************為氣體分配材料屬性**********************
*Initial Conditions, type=VOLUME FRACTION
Set-inner-gas, gas-1.gas-inner, 1.
展開 Simulink中4種電機建模方式
同理，對電機運動微分方程進行拉式變化有：
Te(S)=(J*S+B)*w(S)
其傳遞函數為：
w(S)/Te(S)=1/(J*S+B)
使用Transfer Fcn模塊，完成以上的傳遞函數搭建，模型如下圖。
03
狀態方程建模
狀態方程建模，就是根據現代控制理論中的狀態方程和輸出方程，使用Simulink中State-Space模塊完成建模。
狀態方程，描述系統狀態變量與系統輸入之間關系的一階微分方程組。
輸出方程，描述系統輸出向量與系統狀態向量和系統輸入向量之間的函數關系式。
狀態方程和輸出方程可表達為：
其中，A是系統的狀態矩陣，B是系統的輸入矩陣，C是系統的輸出矩陣，D是系統的直接轉移矩陣。
一般，可以通過對微分方程進行變換，得到狀態方程和輸出方程。
展開 
詳解LS-DYNA爆炸仿真計算的模型與算法
0
3
狀態方程
空氣
空氣采用理想氣體狀態方程描述：
其中P 為空氣壓力；γ 為多方指數；ρ 為空氣現時密度，ρ0為初始密度；E為空氣能量密度。理想氣體狀態方程可以通過設置LS-DYNA 程序中預定義的線性多項式狀態方程的相關常數得到。
爆炸產物
目前存在的爆炸產物狀態方程，按照是否顯含產物組分可以分為兩類，一類是顯含產物組分的狀態方程，如BKW 方程、KHT 方程等，另一類是不顯含產物組分的狀態方程，如多方方程、JWL 方程等，其中JWL 方程是數值計算中應用最為廣泛的狀態方程。
JWL 方程是1968 年由美國LLNL 的E. L. Lee 在H. Jones 和M. L. Wilkins 的工作基礎上提出的半經驗狀態方程。JWL 等熵方程形式如下：
其中為相對比容，即現時比容與初始比容之比，A 、B 、C 、1 R 、2 R 、ω 為六個參數。
根據熱力學第一定律和等熵條件，比內能增量，將(0.8)式代入并積分可得比內能。
根據式(0.8)和(0.9)，可消去參數C ，得到
其中為以初始體積表示的能量密度，初始能量密度。
展開 詳解LS-DYNA爆炸仿真計算的模型與算法
0
3
狀態方程
空氣
空氣采用理想氣體狀態方程描述：
其中P 為空氣壓力；γ 為多方指數；ρ 為空氣現時密度，ρ0為初始密度；E為空氣能量密度。理想氣體狀態方程可以通過設置LS-DYNA 程序中預定義的線性多項式狀態方程的相關常數得到。
爆炸產物
目前存在的爆炸產物狀態方程，按照是否顯含產物組分可以分為兩類，一類是顯含產物組分的狀態方程，如BKW 方程、KHT 方程等，另一類是不顯含產物組分的狀態方程，如多方方程、JWL 方程等，其中JWL 方程是數值計算中應用最為廣泛的狀態方程。
JWL 方程是1968 年由美國LLNL 的E. L. Lee 在H. Jones 和M. L. Wilkins 的工作基礎上提出的半經驗狀態方程。JWL 等熵方程形式如下：
其中為相對比容，即現時比容與初始比容之比，A 、B 、C 、1 R 、2 R 、ω 為六個參數。
根據熱力學第一定律和等熵條件，比內能增量，將(0.8)式代入并積分可得比內能。
根據式(0.8)和(0.9)，可消去參數C ，得到
其中為以初始體積表示的能量密度，初始能量密度。
展開 Dyna求解的工程爆破模擬—鋼筋混凝土結構抗爆
13.介質及炸藥匹配材料模型、狀態方程、及ALE屬性
（1）清空HyperMesh模型；
（2）導入MAIN.k文件；
（3）介質匹配材料模型；
（4）介質匹配ALE單元屬性；
（5）介質匹配狀態方程；
（6）炸藥匹配材料模型；
（7）炸藥匹配ALE單元屬性；
（8）炸藥匹配狀態方程；
（10）導出設置結束模型。
14.提交計算
我們應用TCE工具提供的調用Dyna求解器功能，提交計算文件。求解器功能如何生成到HyperMesh界面上，接下來的文章會進行介紹，我們這里只使用。TCE插件請進群：551922835獲取。
TNT炸*藥_多點延時起爆【JWL狀態方程】 ￥89.9
<p>在歐拉域內有三個TNT炸*藥，設定依次起爆</p><p><img src="https://img.jishulink.com/upload/202102/18e45d82dd0a451588d687b3acc6111c.png" alt="微信截圖_20210203215011.png"></p><p>時間間隔如下所示：</p><div contenteditable="false" width="100%"><img src="https://img.jishulink.com/upload/202102/8f8af335ee82455cab7fb131b9c38b39.png" title="微信截圖_20210203215114.png" alt="微信截圖_20210203215114.png" style="max-width:760px;" data-mobile-src="https://img.jishulink.com/upload/202102/8f8af335ee82455cab7fb131b9c38b39.png?image_process=/format,webp/quality,q_40/resize,w_400" data-pc-src="https://img.jishulink.com/upload/202102/8f8af335ee82455cab7fb131b9c38b39.png?image_process=/format,webp/quality,q_40/resize,w_760" data-initial-src="https://img.jishulink.com/upload/202102/8f8af335ee82455cab7fb131b9c38b39.png">
</div><p>正常情況下會按照設定時間，依次起爆
展開 