TAGGED: ansysautodyn, autodyn, equationofstate, explicitautodyn


April 29, 2024 at 8:55 pmShuanTai YehSubscriber
Hi everyone,
I'm currently working on implementing an user defined EOS (Equation of State) using a lookup table similar to SESAME EOS in Autodyn. I've managed to modify the userdefined EOS in theÂ
mdeos_user_1.f90
file, but I'm unsure about how to include a table in the EOS. If anyone has any advice or suggestions, I would greatly appreciate your assistance.Thanks in advance!

April 30, 2024 at 9:50 pmChris QuanAnsys Employee
You can refer to the Section 7.1. Subroutine MDSTR_USER_1 from Autodyn User's Subroutines Guide:Â
7.1. Subroutine MDSTR_USER_1 (ansys.com)
It has an example on how to use onedimentional array in Autodyn User Subroutines. Below is the definition of array in this example:
! FOR EACH REAL INPUT PARAMETER, ASSIGN DATA ! ('name ' L, T,M,H, val, min,max,default,0,required) EQ%PAR(1)=PRMT (1,'Shear Modulus',1,2,1,0,ZERO,ZERO,BIG,ZERO ,0,1) ! THIS LINE MUST ALWAYS EXIST EQ%PAR(2)=PRMT (2,'EPS #1' , 0, 0,0,0,ZERO,ZERO,BIG,ZERO ,0,0) EQ%PAR(3)=PRMT (3,'EPS #2' , 0, 0,0,0,ZERO,ZERO,BIG,ZERO ,0,0) EQ%PAR(4)=PRMT (4,'EPS #3' , 0, 0,0,0,ZERO,ZERO,BIG,ZERO ,0,0) EQ%PAR(5)=PRMT (5,'YIELD #1' ,1,2,1,0,ZERO,ZERO,BIG,ZERO ,0,0) EQ%PAR(6)=PRMT (6,'YIELD #2' ,1,2,1,0,ZERO,ZERO,BIG,ZERO ,0,0) EQ%PAR(7)=PRMT (7,'YIELD #3' ,1,2,1,0,ZERO,ZERO,BIG,ZERO ,0,0)
Be aware that L, T, M, H in the Fortran statements above represent the exponent order of the unit Length, Time, Mass, and Temprature where a parameter has.
For example, shear modulus has the unit of Force/Area = Mass*Acceleration/Area.
If the system of unit is SI units of mks, the unit of shear modulus = Kg*m/s^2/m^2 = Kg*s^2*m^1. Thus, L=1, T=2, M=1, and H=0 in the parameter definition of shear modulus.Â In the user subroutine "mdeos_user_1.f90", you can
 Define the size of a table in the MODULE EOS_USER_1.Â
 Initialize the table in the SUBROUTINE INIT_EOS_USER_1.
 Assign values to the table in the SUBROUTINE SET_EOS_USER_1.
Â

 You must be logged in to reply to this topic.
 Error messages ‘The demanded solver was not found in the specified path’
 CESE Chemistry solver fails in 2D Axisymmetric case
 Problems with fatigue calculations
 How start learning LSDyna
 Setting up contacts between an ALE body and a Lagrange body
 Fatigue calculation scaling solution
 Successful explicit analysis
 LSDyna Restart
 Expected Fatigue cycles shows 0
 Add materials from Workbench Lsdyan using engineering data

341

160

130

117

99
Â© 2024 Copyright ANSYS, Inc. All rights reserved.