{"id":293199,"date":"2023-07-15T06:57:45","date_gmt":"2023-07-15T06:57:45","guid":{"rendered":"\/forum\/forums\/topic\/how-to-write-ideal-gas-equation-of-state-in-ansys-autodyn-using-fortran\/"},"modified":"2023-07-15T06:57:45","modified_gmt":"2023-07-15T06:57:45","slug":"how-to-write-ideal-gas-equation-of-state-in-ansys-autodyn-using-fortran","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/how-to-write-ideal-gas-equation-of-state-in-ansys-autodyn-using-fortran\/","title":{"rendered":"How to write ideal gas equation of state in ansys autodyn using fortran?"},"content":{"rendered":"
How to write ideal gas equation of state (EOS) in autodyn using fortran?<\/div>\n
 <\/div>\n
I modified the Autodyn user subroutine mdeos_user_1.f90, writing a custom ideal state equation. Ea0 (E0) is treated as a constant, 2.068 × 105Kj\/kg, and the simulation is still more consistent with the no user subroutines results. My code is as follows\uff0c<\/div>\n
IF_SIMULTANEOUS = 0<\/div>\n
ARHO = (1.4_real8-ONE) * E0 * DEN<\/div>\n
BRHO = ZERO<\/div>\n
DADETA = ((1.4_real8 – ONE)* E0)\/ RHOREF<\/div>\n
DBDETA = ZERO<\/div>\n
******<\/div>\n
(1) Properties of Air<\/div>\n
The air and TNT gas product media were assumed to behave as an ideal gas, thus the default Equations of state (1) were specified. The air was specified with an internal energy of 2.068 × 105Kj\/kg in AUTODYN 2D.<\/div>\n
?<\/div>\n
where, ‘ρ’ is the density of gas, ‘cv’ and ‘cp’ are specific heats at constant volume and pressure, respectively; ‘T’ is the gas temperature, ‘P’ is the pressure, ‘γ’ is the ratio of specific heats and ′Ea0′ is the specific energy. The air model was taken from AUTODYN 2D material library as shown in Table 2.<\/div>\n
Reference\uff1aNumerical Simulations of Exposed Circular Surface of Rolled<\/em><\/div>\n
Homogeneous Armor Steel Plates Subjected to Blast Loadings by Using AUTODYN 2D<\/em><\/div>\n","protected":false},"template":"","class_list":["post-293199","topic","type-topic","status-closed","hentry","topic-tag-ansys-autodyn-1","topic-tag-autodyn"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t