Marc > Building A Model > Material Library
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX''">XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX''">   
Material Library
The Materials application defines Marc materials which are later associated to the elements of the model in the Element Properties application described in the next section, Element Properties.
The following tables outlines the available options that can be created for Structural, Thermal, and Coupled analyses.
 
Isotropic/Orthotropic/Anisotropic
Constitutive Model
2D Conditions
Method
Plane Stress / Thin Shell
Plane Strain / Axisymmetric
Thick Shell
Axisymmetric with Twist
Axisymmetric Shell
None (Isotropic and 3D cases)
Entered Values
User Subs. ANELAS ANEXP (Anisotropic Only)
Constitutive Model
Failure Criterion
Failure Option
Hill
Hoffman
Tsai-Wu
Maximum Strain
Maximum Stress
User Sub. UFAIL
Default
Progressive Failure
Constitutive Model
Model
Domain Type
Number of Terms
Neo-Hookean
Mooney-Rivlin
Full 3rd Order
Time
Frequency
1
Ogden
Foam
Time
1 - 6
Arruda-Boyce
Gent
Time
1
User Sub. (UELASTOMER)
Ogden
Foam-Invariants
Foam-Principals
Foam-Invariants (Deviatoric Split)
Foam-Principals (Deviatoric Split)
Constitutive Model
Thermal Expansion
Stress-Strain Law
Entered Values
User Sub. ANEXP
User Sub. HYPELA
User Sub. HYPELA2 (Grad/Rot)
User Sub. HYPELA2 (Grad/Str)
User Sub. HYPELA2 (All Input)
User Sub. UBEAM
Constitutive Model
Shift Function
No Function
Williams-Landel-Ferry
Power Series Expansion
Narayanaswamy Model
User Sub. TRSFAC
Constitutive Model
Method
Creep
Power Law - Piecewise
User Sub.CRPLAW
Constitutive Model
Dmping
Constitutive Model
Method
Entered Values
User Subs. ANKOND ORIENT
Constitutive Model
Memory Model
Mechanical (Auricchio)
Thermal Mechanical
Constitutive Model
Damage Type
Damage Model
Damage
Elastic/Plastic
No Nucleation
Plastic Strain Control Nucleation
Stress Control Nucleation
User Sub. UVOIDN
Elastomer (Rubber)
(Isotropic Only)
Additive Decomposition
Multiplicative Decompostion
User Sub. UELDAM
Simple
(Isotropic Only)
Yield- User Sub. UDAMAG
Yield/Youngs Mod. (UDAMAG)
Constitutive Model
Method
Entered Values
User Subs. UCRACK...
Constitutive Model
Method
Fitted
Predicted
Table
Constitutive Model
Method
Yada
User Sub. UGRAIN
Constitutive Model
Model
Linear
Cam Clay
User Sub.HYPELA
Constitutive Model
Method
Entered Values
User Sub. UPOWDR
Constitutive Model
Model
Entered Values
Entered Values
Magnetostatic (p. 109)
Entered Values
User Sub UMU
 
 
 
Piezoelectric (p. 109)
Stress Based
Strain Based
 
 
 
 
Isotropic/Orthotopic/Anisotropic
Constitutive Model
Type
Hardening Rule
Yield Criteria
Strain Rate Method
Elastic-Plastic
Isotropic
Kinematic
Combined
von Mises
Hill Yield
Barlat
Linear Mohr-Coulomb (Isotropic Only)
Parabolic Mohr-Coulomb (Isotropic Only)
Buyukozturk Concrete (Isotropic Only)
Oak Ridge National Lab
2-1/4 Cr-Mo ORNL
Reversed Plasticity ORNL
Full Alpha Reset ORNL
Generalized Plasticity
Piecewise Linear
Cowper-Symonds
Power Law (Isotropic only)
Rate Power Law (Isotropic only)
Johnson-Cook (Isotropic only)
Kumar (Isotropic only)
Chaboche (Isotropic only)
Viscoplastic (UVSCPL) (Isotropic, Orthotropic only)
Constitutive Model
Type
Hardening Rule
Yield Criteria
Strain Rate Method
Perfectly Plastic
None
von Mises
Linear Mohr-Coulomb
Hill Yield
Barlat
Linear Mohr-Coulomb (Isotropic Only)
Parabolic Mohr-Coulomb (Isotropic Only)
Buyukozturk Concrete (Isotropic Only)
Oak Ridge National Lab
2-1/4 Cr-Mo ORNL
Reversed Plasticity ORNL
Full Alpha Reset ORNL
Generalized Plasticity
Piecewise Linear
Cowper-Symonds
 
Rigid-Plastic (Isotropic only)
Power Law
Rate Power Law
Johnson-Cook
Kumar
Piecewise Linear
None
Piecewise Linear
Cowper-Symonds
Material Input Properties
This is an example of one of many Input Properties forms that can appear when defining material properties. There is a Constitutive Model plus other optional selections followed by places for input of specific property parameters.
For each material type, see the following pages: Isotropic (p. 81), 2D Orthotropic (p. 101), 3D Orthotropic (p. 101), 2D Anisotropic (p. 81), 3D Anisotropic (p. 81), or Composite (p. 110). For thermal material property definitions see (p. 94).
 
Note:  
For Coupled analysis, the thermal properties are also presented along with the structural. The thermal properties are listed in Thermal - Isotropic / Orthotropic / Anisotropic.
Elastic - Isotropic / Orthotropic / Anisotropic
This input data creates the ISOTROPIC and INITIAL STATE keyword options.
 
Elastic - Isotropic
Description
Method (Coupled only)
User Subs. ANKOND ORIENT - writes a 1 to the 4th field of the 3rd datablock of the ISOTROPIC option. Entered Values allows for the properties in this table to be entered.
Elastic Modulus
Defines the elastic modulus. It is entered in the first data field on the fourth card of the ISOTROPIC option. This property is generally required. May vary with temperature via a defined material field and placed on 4b data block of the TEMPERATURE EFFECTS option.
Poisson’s Ratio
Defines the Poisson’s ratio. It is entered in the second data field on the fourth card of the ISOTROPIC option. This property is generally required. May vary with temperature via a defined material field and placed on 5b data block of the TEMPERATURE EFFECTS option.
Density
Defines the mass density. It is entered in the third data field on the fourth card of the ISOTROPIC option. This property is optional.
Coefficient of Thermal
Expansion
Defines the instantaneous coefficient of thermal expansion. This is entered in the fourth data field on the fourth card of the ISOTROPIC option. This property is optional. May vary with temperature via a defined material field and placed on 6b data block of the TEMPERATURE EFFECTS option.
Reference Temperature
Defines the reference temperature for the thermal expansion coefficient. It is entered in the first data field on the fourth card of the INITIAL STATE option. This property is optional. When defining temperature dependent properties, this is the reference temperature from which values will be extracted or interpolated for the WORK HARD and STRAIN RATE options. See note below.
Cost per Unit Volume
For design optimization, entered on the 7th field of the 4th data block of the ISOTROPIC option.
Cost per Unit Mass
For design optimization, entered on the 8th field of the 4th data block of the ISOTROPIC option.
Latent Heat vs Solidus Temp.
Latent Heat vs Liquidus Temp.
(Coupled only)
Both of these should be present. If one is missing you must treat all the temperature values as zero for the missing one. When both are present, they must reference Temperature material fields and they must all have exactly the same number of latent heats in them (with the same values). For Coupled analysis, the TEMPERATURE EFFECTS option is written with the values in block 11b and the number of latent heats in field 9 of block 2b.
This input data creates the ORTHOTROPIC and INITIAL STATE keyword options. The required properties vary based on dimension and element type which for a 2D Orthotropic option can be set to either Plane Stress/Thin Shell, Plane Strain/Axisymmetric, Thick Shell, Axisymmetric with Twist, or Axisymmetric Shell.
Elastic - Orthotropic
Description
Method (Coupled only)
User Subs. ANKOND ORIENT - writes a 1 to the 4th field of the 3rd datablock of the ORTHOTROPIC option. Entered Values allows for the properties in this table to be entered.
Elastic Modulus 11/22/33
Defines the elastic moduli in the element’s coordinate system. They are entered in the first through third data fields on the fourth card of the ORTHOTROPIC option. This is required data. May vary with temperature via a defined material field and placed on 5b, 6b, and 7b data blocks of the ORTHO TEMP option.
Poisson’s Ratio 12/23/31
Defines the Poisson’s ratios relative to the element’s coordinate system. They are entered in the fourth through sixth data fields on the fourth card of the ORTHOTROPIC option. This is required data. May vary with temperature via a defined material field and placed on 8b, 9b, and 10b data blocks of the ORTHO TEMP option.
Shear Modulus 12/23/31
Defines the shear moduli relative to the element’s coordinate system. They are entered in the first through third data fields on the fifth card of the ORTHOTROPIC option. This is required data. May vary with temperature via a defined material field and placed on 11b, 12b, and 13b data blocks of the ORTHO TEMP option.
Coefficient of Thermal
Expansion 11/22/33
Defines the instantaneous coefficients of thermal expansion relative to the element’s coordinate system. They are entered in the fourth through sixth data fields on the fifth card of the ISOTROPIC option. These properties are optional. This is required data. May vary with temperature via a defined material field and placed on 14b, 15b, and 16b data block of the ORTHO TEMP option.
Reference Temperature
Defines the reference temperature for the thermal expansion coefficient. It is entered in the first data field on the fourth card of the INITIAL STATE option. When defining temperature dependent properties, this is the reference temperature from which values will be extracted or interpolated for the WORK HARD and STRAIN RATE options. See note below.
Density
Defines the mass density which is an optional property. It is entered in the seventh data field on the fourth card of the ORTHOTROPIC option.
Cost per Unit Volume
For design optimization, entered on the 7th field of the 5th data block of the ORTHOTROPIC option.
Cost per Unit Mass
For design optimization, entered on the 8th field of the 5th data block of the ORTHOTROPIC option.
Latent Heat vs Solidus Temp.
Latent Heat vs Liquidus Temp.
(Coupled only)
Both of these should be present or none. If one is missing the temperature values are treated as zero for the missing one. When both are present, they must reference Temperature material fields and they must all have exactly the same number of latent heats in them (with the same values). For Coupled analysis, the TEMPERATURE EFFECTS option is written with the values in block 11b and the number of latent heats in field 9 of block 2b.
This input data creates the ANISOTROPIC and INITIAL STATE keyword options. The required properties vary based on dimension and element type which for a 2D Anisotropic option can be set to either Plane Stress/Thin Shell, Plane Strain/Axisymmetric, Thick Shell, Axisymmetric with Twist, or Axisymmetric Shell.
 
Elastic - Anisotropic
Description
Method
User Subs. ANELAS ANEXP ...- writes a 1 to 4th field of 3rd datablock of the ANISOTROPIC option - datablocks 4a-f not written.  Entered Values allows for the properties in this table to be entered.
Stress-Strain Matrix, Cij
Defines the upper right portion of the symmetric stress-strain matrix relative to the element’s coordinate system. They are entered on the 4a, 4b and 4c card of the ANISOTROPIC option.
Coefficient of Thermal
Expansion 11/22/33/12/23/31
Defines the instantaneous coefficients of thermal expansion relative to the element’s coordinate system. They are entered on the 4d card of the ANISOTORPIC option, and are optional properties.
Reference Temperature
Defines the reference temperature for the thermal expansion coefficient. It is entered in the first data field on the fourth card of the INITIAL STATE option. When defining temperature dependent properties, this is the reference temperature from which values will be extracted or interpolated for the WORK HARD and STRAIN RATE options. See note below.
Density
Defines the mass density which is an optional property. It is entered in the fourth data field on the fourth card of the ANISOTROPIC option.
Cost per Unit Volume
For design optimization, entered on the 7th field of the 4th data block of the ANISOTROPIC option.
Cost per Unit Mass
For design optimization, entered on the 8th field of the 4th data block of the ANISOTROPIC option.
Latent Heat vs Solidus Temp.
Latent Heat vs Liquidus Temp.
(Coupled only)
Both of these should be present. If one is missing you must treat all the temperature values as zero for the missing one. When both are present, they must reference Temperature material fields and they must all have exactly the same number of latent heats in them (with the same values). For Coupled analysis, the TEMPERATURE EFFECTS option is written with the values in block 11b and the number of latent heats in field 9 of block 2b.
Note:  
Note on reference temperature. If the reference temperature is left blank, zero is assumed. If the reference temperature does not fall between temperature values defined for work hardening or strain rate, the highest or lowest values will be used depending on whether the reference temperature is greater or lower than the given temperature range. If it falls inbetween, then values are interpolated. For Structural analysis, if Nodal LBC Temperatures (POINT TEMP) also exist then the INITIAL STATE will not be written since this is incompatible.
Failure - Isotropic / Orthotropic / Anisotropic
This input data creates the FAIL DATA option. The first data field of the fourth card is set to either HILL, HOFFMAN, TSAI-WU, MX STRAIN (maximum strain), MX STRESS (maximum stress) or User Sub. UFAIL. A number of the following input properties will appear depending on the material type and options set. Note that there are three Failure constitutive models: Failure, Failure 2, and Failure 3. This means that you can have up to three failure criteria per material model
.
Failure Criteria -
Hill, Hoffman, Tsai-Wu, Maximum Stress/Strain
Description
Failure Option
Progressive Failure - writes a one (1) in the 3rd field of the 3rd data block of the FAIL DATA option for each criterion defined with this option set.
Max Tensile Stress X, Y & Z
Defines the tension stress (or strain) limits in the element’s coordinate system. 2nd, 4th and 6th fields of 4th datablock of FAIL DATA option, respectively.
Max Compressive
Stress X, Y & Z
Defines the compression stress (or strain) limits in the element’s coordinate system. 3rd, 5th, and 7th field of 4th datablock of FAIL DATA option. Absolute values are used.
Max Shear Stress XY, YZ, ZX
Defines the shear stress (or strain) limits. 1st, 2nd and 3rd fields of 5th datablock of FAIL DATA option, respectively.
Failure Index
4th field of 5th datablock of FAIL DATA option.
Interactive Term XY, YZ, & ZX
Defines the stress interaction parameters. 5th, 6th, and 7th fields of 5th datablock of FAIL DATA option.
Note:  
When User Sub. UFAIL is used, no input data is necessary and the word UFAIL is written in the 4th data block of the FAIL DATA option.
Hyperelastic - Isotropic
The following Hyperelastic models can be created.
 
Caution:  
If one of these constitutive models exists and is active, the Elastic or Plastic constitutive models must be turned off (made inactive) otherwise ISOTROPIC, WORK HARD and MOONEY or some other hyperelastic option will be written to the input file which will cause an incompatibility in the analysis.
 
Neo-Hookean,
Mooney-Rivlin,
Full 3rd Order Invariant
Time Domain
Description
Strain Energy Function, C10, C01, C11, C20, C30
Strain energy densities as a function of the strain invariants in the material. Creates MOONEY option; 1st, 2nd, 5th, 6th, and 7th fields of 4th data block, respectively. May vary with temperature via a defined material field and placed on 4b data block of the TEMPERATURE EFFECTS option.
Density
Defines the mass density which is an optional property. It is entered in the third data field on the fourth card of the MOONEY option.
Coefficient of Thermal
Expansion
Defines the instantaneous coefficient of thermal expansion. This is entered in the fourth data field on the fourth card of the MOONEY option. This property is optional. May vary with temperature via a defined material field and placed on 6b data block of the TEMPERATURE EFFECTS option.
Bulk Modulus
8th field of 4th data block of MOONEY option.
Reference Temperature
Defines the reference temperature for the thermal expansion coefficient. It is entered in the first data field on the fourth card of the INITIAL STATE option.
For Neo-Hookean, Mooney-Rivlin and Full 3rd Order in the Frequency Domain the additional
inputs are:
 
Neo-Hookean
Frequency Domain
Description
, Real and Imaginary
Creates PHI-COEFFICIENTS option. One PHI-COEFFICIENTS option is created for each pair of real and imaginary PHIs that has input. Input is a material field of frequency versus value. This frequency, real and imaginary phi coefficients are entered into the 1st, 2nd, and 3rd fields of the 3rd data block respectively.
 
Ogden
Description
Bulk Modulus K
Creates OGDEN option; 1st field of 4th data block.
Density
2nd field of 4th data block of OGDEN option.
Coefficient of Thermal
Expansion
3rd field of 4th data block of OGDEN option.
Reference Temperature
Creates INITIAL STATE option. Defines the reference temperature for the thermal expansion coefficient.
Modulus 1
1st field of 6th data block of OGDEN option.
Exponent 1
2nd field of 6th data block of OGDEN option.
Note:  
Modulus 1 and Exponent 1 will repeat for the Number of Terms and will increment as such, e.g., Modulus 2, Exponent 2 - Modulus 3, Exponent 3, etc. Same comment applies to FOAM option for repeating terms.
Foam
Description
Density
Creates FOAM option; 2nd field of 4th data.
Coefficient of Thermal
Expansion
3rd field of 4th data block of FOAM option.
Reference Temperature
Creates INITIAL STATE option. Defines the reference temperature for the thermal expansion coefficient.
Modulus 1
1st field of 6th data block of FOAM option.
Deviatoric Exponent 1
2nd field of 6th data block of FOAM option.
Volumetric Exponent 1
3rd field of 6th data block of FOAM option.
Arruda-Boyce
Description
NKT
Creates the ARRUDBOYCE option: 1st field of 4th data block. May vary with temperature via a defined material field and placed on 4b data block of the TEMPERATURE EFFECTS option.
Chain Length
2nd field of 4th data block of ARRUDBOYCE option. May vary with temperature via a defined material field and placed on 5b data block of the TEMPERATURE EFFECTS option.
Bulk Modulus
5th field of 4th data block of ARRUDBOYCE option.
Density
3rd field of 4th data block of ARRUDBOYCE option.
Coefficient of Thermal
Expansion
4th field of 4th data block of ARRUDBOYCE option.
Reference Temperature
Creates INITIAL STATE option. Defines the reference temperature for the thermal expansion coefficient.
Gent
Description
Tensile Modulus
Creates the GENT option: 3rd field of 4th data block. May vary with temperature via a defined material field and placed on 4b data block of the TEMPERATURE EFFECTS option.
Maximum 1st Invariant
4th field of 4th data block of GENT option. May vary with temperature via a defined material field and placed on 5b data block of the TEMPERATURE EFFECTS option.
Bulk Modulus
5th field of 4th data block of GENT option.
Density
1st field of 4th data block of GENT option.
Coefficient of Thermal
Expansion
2nd field of 4th data block of GENT option.
Reference Temperature
Creates INITIAL STATE option. Defines the reference temperature for the thermal expansion coefficient.
User Sub. UELASTOMER
Description
Domain Type
The User Sub. UELASTOMER can be used with the Ogden or Foam model. If Ogden is selected, this places a 3 in the 3rd field of the 3rd datablock of the OGDEN option. If a Foam model is selected, it places a 1, 2, 3, or 4, respectively, in the 4th field of the 3rd datablock of the FOAM option. No terms are required if this user subroutine is selected for either Ogden or Foam.
Bulk Modulus K
Creates OGDEN option; 1st field of 4th data block.
Density
2nd field of 4th data block of OGDEN option. OR Creates FOAM option; 2nd field of 4th data.
Coefficient of Thermal
Expansion
3rd field of 4th data block of OGDEN option. OR 3rd field of 4th data block of FOAM option.
Reference Temperature
Creates INITIAL STATE option. Defines the reference temperature for the thermal expansion coefficient.
Note:  
Marc may force you to use a Herrmann formulated element when using some Hyperelastic constitutive models.
Hypoelastic - Isotropic
The following Hypoelastic models can be created. The HYPOELASTIC option is written to the input file. This constitutive model requires the use of user subroutines as explained below.
Hypoelastic
Description
Thermal Expansion
 
User Sub. ANEXP: This places a 1 in 2nd field of the 3rd data block of the HYPOELASTIC option. Otherwise it is zero (default).
Stress-Strain Law
User Sub. HYPELA or UBEAM flags use of the HYPELA or UBEAM user subroutines which is default and a zero is placed in the 3rd field of the 3rd data block of the HYPOELASTIC option. If HYPELA2 is selected, the 3rd field is set according to Rotation (Grad/Rot), Stretch Ratio (Grad/Str) or Both (All Input) which puts a 1, 2, or 3, respectively in the 3rd field of the 3rd data block.
Density
Defines the mass density which is an optional property. It is entered in the 1st data field on the fourth card of the HYPOELASTIC option and in the 6th field for Coupled or Thermal analysis.
Coefficient of Thermal
Expansion
Defines the instantaneous thermal expansion coefficient which is an optional property. It is entered in the 2nd data field on the fourth card of the HYPOELASTIC option.
Conductivity
Defines the thermal conductivity which is an optional property. It is entered in the 3rd data field on the fourth card of the HYPOELASTIC option.
Specific Heat
Defines the specific heat which is an optional property. It is entered in the 4th data field on the fourth card of the HYPOELASTIC option.
Reference Temperature
Defines the reference temperature for the thermal expansion coefficient. It is entered in the first data field on the fourth card of the INITIAL STATE option.
Emissivity
Defines the emissivity which is an optional property. It is entered in the 7th data field on the fourth card of the HYPOELASTIC option.
A TEMPERATURE EFFECTS option is written for items above that accept temperature dependent field references.
Viscoelastic - Isotropic / Orthotropic
This input data creates the VISCELPROP, VISCELMOON, VISCELOGDEN, or VISCELORTH options. The Prony series are defined in Fields - Tables as material properties with time (relaxation time) as their independent variable and then selected here as input properties. All inputs must have the same number of time points (at the same times) in the referenced fields. The following equations may be useful when creating the Prony series for the bulk and shear moduli: This also supports the SHIFT FUNCTION option for Thermo-Rheologically simple viscoelastic materials. The SHIFT FUNCTION is written for ISOTROPIC, ORTHOTROPIC, MOONEY, OGDEN, ARRUDA-BOYCE, & GENT models if present in the defined material.
Viscoelastic - Isotropic
Description
Shift Function
Enters a 1, 2, 3, or -1 in the 2nd field of the 3rd data block of SHIFT FUNCTION to specify the type of function: Williams-Landel-Ferry, Power Serires, Narayanaswamy, User Sub. TRSFAC. If the latter, no other data blocks are required. Input properties for the different shift functions are listed in this table.
Shear Constant
If a material field of time vs. value is supplied, will create a VISCELPROP option. This is valid when an Elastic and/or Plastic constitutive model is present. Fills out 1st and 2nd fields of 4th data block for the number of terms present in the field.
Bulk Constant
Same as above. Fills out 3rd and 4th fields of 4th data block for the number of terms present in the field. (Field code 5)
Energy Function Multiplier
Defines the duration effect on the hyperelastic model as a multiplier to the strain energy density function. If a material field of time vs. value is supplied, will create a VISCELMOON option. This is valid when a Hyperelastic constitutive model for Neo-Hookean, Mooney-Rivlin, Full 3rd Order, Arruda-Boyce, or Gent is present. Fills out the 4th data block for the number of terms present in the field. (Field code 5)
Deviatoric Multiplier
If a material field of time vs. value is supplied, will create a VISCELOGDEN option. This is valid when a Hyperelastic constitutive model of Ogden is present. Fills out 1st and 2nd fields of 4th data block for the number of terms present in the field. (Field code 5)
Dilatational Multiplier
Same as above. Fills out 3rd and 4th fields of 4th data block for the number of terms present in the field. (Field code 5)
Solid Coeff of Thermal Exp
If input is supplied, will create a VISCEL EXP option; 2nd field of 3rd data block.
Liquid Coeff of Thermal Exp
3rd field of 3rd data block of VISCEL EXP option.
Reference Temperature
For all Shift Functions except None, 4th field of 3rd data block of SHIFT FUNCTION option.
Constant C1
For Shift Function 1 only - Field 1, 4th data block
Constant C2
For Shift Function 1 only - Field 2, 4th data block
Constant Coefficients Co-Cm
For Shift Function 2 only - data block 4 - must be defined by a 1D material field where the independent value is arbitrary. The first value is Co and the number of field entries is placed in 3rd field of 3rd data block.
Activation Energ/ Gas Const.
For Shift Function 3 only - field 5, data block 3
Structural Relax. Ref. Temp.
For Shift Function 3 only - field 8, data block 3
Fraction Parameter
For Shift Function 3 only - field 6, data block 3
Abs Temperature Shift
For Shift Function 3 only - field 7, data block 3
Weighting Factors
For Shift Function 3 only - data blocks 4 & 5 where this is defined by a material time field. Weighing factor values are written to data block 4, and time values are written to datablock 5.
Note:  
Instantaneous values are entered for the elastic model, and the difference between the instantaneous value and the summation of the values in the series is the long-term property value.
Viscoelastic - Orthotropic
Description
Shift Function
Enters a 1, 2, 3, or -1 in the 2nd field of the 3rd data block of SHIFT FUNCTION to specify the type of function: Williams-Landel-Ferry, Power Serires, Narayanaswamy, User Sub. TRSFAC. If the latter, no other data blocks are required. Input properties for the different shift functions are listed in the table above for Isotropic.
Youngs Modulus, E11/E22/E33
Defines the duration effects on the elastic moduli. This information is entered on the 2nd, 3rd, and 4th fields of the 4th datablock of the VISCELORTH option, and is optional. This is only valid when an Elastic and/or Plastic constitutive model is present.
Poissons Ratio 12/23/31
Defines the duration effects on the Poisson’s ratios. This information is entered on the 5th, 6th, and 7th fields of the 4th datablock of the VISCELORTH option, and is optional.
Shear Modulus G12/G23/G31
Defines the duration effects on the shear moduli. This information is entered on the fifth card of the VISCELORTH option, and is optional.
Solid Coeff of Thermal Exp
Same as for Isotropic
Liquid Coeff of Thermal Exp
Same as for Isotropic
Creep - Isotropic / Orthotropic / Anisotropic
The following input is for the Creep constitutive model. This places a CREEP option in the input file.
Creep
Description
Method
User Sub. CRPLAW - writes a zero in the 5th field of the 2nd data block of the CREEP option. No other data blocks beyond are written. User subroutine UCRPLW will automatically get called if it exists if Implicit creep is set.
Power Law - Piecewise allows for input of the material properties in the table below.
Coefficient
Creates the CREEP option. It is compatible with all other constitutive models except Viscoelastic and Hyperelastic. This is 5th field in 2nd data block.
Exponent of Temperature
1st field of 3rd data block.
Temperature vs. Creep Strain
References a material field of temperature vs. value. Overrides Exponent of Temperature if present. Fills out 3rd data block.
Exponent of Stress
1st field of 4th data block.
Creep Strain vs. Stress
References a material field of stress vs. value. Overrides Exponent of Stress if present. Fills out 4th data block.
Exponent of Creep Strain
1st field of 5th data block.
Strain Rate vs. Creep Strain
References a material field of strain rate vs. value. Overrides Exponent of Creep Strain if present. Fills out 5th data block.
Exponent of Time
1st field of 6th data block.
Time vs. Creep Strain
References a material field of time vs. value. Overrides Exponent of Time if present. Fills out 6th data block.
Back Stress
For implicit creep - goes on 5th field of 4th data block of ISOTROPIC option and can vary with strain and/or temperature via a field definition in which case the WORK HARD and/or TEMPERATURE EFFECTS options may be written also.
Damping - Isotropic / Orthotropic / Anisotropic
The following input is for Damping constitutive model. If any one of these values is present, they are placed on a DAMPING option and the element to which the material is associated are referenced. This option is used for harmonic analysis and direct transient dynamic integration only.
Damping
Description
Raleigh Mass Matrix Multiplier
1st field of 4th data block of DAMPING option.
Raleigh Stiff Matrix Multiplier
2nd field of 4th data block of DAMPING option.
Numerical Damping Multiplier
3rd field of 4th data block of DAMPING option.
Thermal - Isotropic / Orthotropic / Anisotropic
This input data creates the ISOTROPIC keyword option for heat transfer analysis.
Thermal - Isotropic
Description
Method
 
User Subs. ANKOND ORIENT - writes a 1 to 2nd field of 3rd datablock of the ISOTROPIC option.
Conductivity
Defines the thermal conductivity. It is entered in the first data field on the fourth card of the ISOTROPIC option. This property is required. May vary with temperature via a defined material field and placed on 9b data block of the TEMPERATURE EFFECTS option.
Specific Heat
Defines the specific heat per unit mass which is an optional property. It is entered in the second data field on the fourth card of the ISOTROPIC option. May vary with temperature via a defined material field and placed on 10b data block of the TEMPERATURE EFFECTS option.
Density
Defines the mass density which is an optional property. It is entered in the third data field on the fourth card of the ISOTROPIC option.
Emissivity
Defines the emmisivity property (5th field of the 5a data block of the ISOTROPIC option). May vary with temperature via a defined material field and placed on 12b data block of the TEMPERATURE EFFECTS option.
Latent Heat vs Solidus Temp.
Latent Heat vs Liquidus Temp.
Both of these should be present or none. If one is missing the temperature values are treated as zero for the missing one. When both are present, they must reference Temperature material fields and they must all have exactly the same number of latent heats in them (with the same values). For Heat Transfer, the TEMPERATURE EFFECTS option is written with the values in the 5b data block. Field 3 of the 2b data block contains the number of latent heats in the fields.
This input data creates the ORTHOTROPIC keyword option for heat transfer analysis.
 
Thermal - Orthotropic
Description
Method
User Subs. ANKOND ORIENT - writes a 1 to 2nd field of 3rd datablock of ORTHOTROPIC option.
Conductivity 11/22/33
Defines the thermal conductivity in the element’s coordinate system. These are entered in the 1st through 3rd data fields on the 4th datablock of the ORTHOTROPIC option, and are required properties.
Specific Heat
Defines the specific heat per unit mass which is an optional property. It is entered in the fifth data field on the fourth card of the ORTHOTROPIC option.
Density
Defines the mass density. It is entered in the fourth data field on the fourth card of the ORTHOTROPIC option. This property is optional.
Emissivity
Defines the emmisivity property (1st field of the 5th data block of the ORTHOTROPIC option). May vary with temperature via a defined material field and placed on 11b data block of the ORTHO TEMP option.
Latent Heat vs Solidus Temp.
Latent Heat vs Liquidus Temp.
Both of these should be present. If one is missing you must treat all the temperature values as zero for the missing one. When both are present, they must reference Temperature material fields and they must all have exactly the same number of latent heats in them (with the same values). For Heat Transfer, the TEMPERATURE EFFECTS option is written with the values in the 5b data block. Field 3 of the 2b data block contains the number of latent heats in the fields.
This input data creates the ANISOTROPIC keyword option for heat transfer analysis.
 
Thermal - Anisotropic
Description
Method
User Subs. ANKOND ORIENT - writes a 1 to 2nd field of 3rd datablock of the ANISOTROPIC option - datablock 4a not written.
Conductivity 11/22/33
Defines the thermal conductivity in the element’s coordinate system. These are entered on the 4a datablock of the ANISOTROPIC option, and are required properties.
Specific Heat
Defines the specific heat per unit mass which is an optional property. It is entered in the 2nd data field on the 4th datablock of the ANISOTROPIC option.
Density
Defines the mass density which is an optional property. It is entered in the 1st data field on the 4th datablock of the ANISOTROPIC option.
Emissivity
Defines the emmisivity property (3rd field of the 4th data block of the ANISOTROPIC option). May vary with temperature via a defined material field and placed on 11b data block of the ORTHO TEMP option.
Latent Heat vs Solidus Temp.
Latent Heat vs Liquidus Temp.
Both of these should be present. If one is missing you must treat all the temperature values as zero for the missing one. When both are present, they must reference Temperature material fields and they must all have exactly the same number of latent heats in them (with the same values). For Heat Transfer, the TEMPERATURE EFFECTS option is written with the values in the 5b data block. Field 3 of the 2b data block contains the number of latent heats in the fields.
Plastic - Isotropic
This input data can create the WORK HARD, TEMPERATURE EFFECTS, STRAIN RATE and the ISOTROPIC keyword options, with the 2nd data field of the 3rd data block of the latter set to VON MISES, LIN MOHRC, PBL MOHRC, BUY MOHRC, NORM ORNL, CRMO ORNL, REVP ORNL, ARST ORNL, GEN-PLAST, RIGID, or VISCO PLAS depending on the Yield Criteria set. One or more of the following input properties will appear depending on the options set:
For Hardening Rules = Isotropic, Kinematic, and Combined, properties for each combination are:
 
Von Mises
Linear Mohr-Coulomb
Parabolic Mohr-Coulomb
Buyukozturk Concrete
ORNL Models
General Plasticity
Description
Stress vs. Plastic Strain

or
Yield Stress
Defines the uniaxial tensile stress versus plastic strain by reference to a tabular field. The field is selected from the Field Definition list. The field is created using the Fields application. See Fields - Tables. It is entered on the third card of the WORK HARD option. For Perfectly Plastic models, only a Yield Stress needs to be entered. See Caution on page 100 below.
Extracts yield stress from first data point from field (zero plastic stain at the reference temperature) for the 5th field of 4th data block of ISOTROPIC option. Can also be temperature dependent which creates TEMPERATURE EFFECTS option.
Can also be strain rate dependent if Strain Rate Method is Piecewise Linear. Accepts field of yield stress vs. strain rate and creates STRAIN RATE option with DATA in 2nd field. Data is input in data block 3 for Option B.
10th Cycle Yield Stress vs.
Plastic Strain

or
10th Cycle Yield Stress
Accepts field of 10th cycle yield stress vs. plastic strain and creates WORK HARD option. Goes on same WORK HARD option as Stress vs. Plastic Strain. 7th field of 4th data block of ISOTROPIC option also extracted from first value of field. Can be temperature dependent also and reference temperature field which creates TEMPERATURE EFFECTS option (data block 7b). For Perfectly Plastic models, only a 10th Cycle Yield Stress needs to be entered.
or 10th Cycle Slope Data
Same as or Break Point Slope Data except for 10th Cycle Yield vs. Strain.
Coefficient C
Visible if Strain Rate Method is Cowper-Symonds. Creates STRAIN RATE option with COWPER in 2nd field. Data is placed in data block 3 for Option C.
Inverse Exponent P
Visible if Strain Rate Method is Cowper-Symonds. Creates STRAIN RATE option with COWPER in 2nd field. Data is placed in data block 3 for Option C.
Alpha
When set to Linear Mohr-Coulomb, defines the slope of the yield surface in square root J2 versus J1 space. It is entered in the sixth data field, on the fourth card of the ISOTROPIC option. This property is required.
Beta
When set to Parabolic Mohr-Coulomb, defines the beta parameter in the equation that defines the parabolic yield surface in square root J2 versus J1 space. It is entered in the sixth data field on the fourth card of the ISOTROPIC option. This property is required.
Note:  
2 1/4 Cr-Mo ORNL, Reversed Plasticity ORNL, Full Alpha Reset ORNL are the same as Oak Ridge National Labs. Generalized Plasticity is the same as Von Mises.
Hill Yield
Barlat
Description
Stress vs. Plastic Strain

or
Yield Stress
Same as table above.
Kinematic Ratio
This is only writen if the Hardening Rule is set to Combined and is written to the 6th field of the 4th data block for ISOTROPIC, the 2nd field of the 6th data block for ORTHOTROPIC, and 3rd field of the 4th data block for ANISOTROPIC.
Stress 11, 22, 33 Yield Ratio
Stress 12, 23, 13 Yield Ratio
These are property words for Hill Yield criterion and are writen to fields 1-6 of the 5th datablock for ISOTROPIC, fields 3-8 of the 6th data block for ORTHOTROPIC, and fields 1-6 or the 4e data block for ANISOTROPIC.
M, C1, C2, C3, C6
These are property words for Barlat criterion and are writen to fields 1-5 of the 5th datablock for ISOTROPIC, fields 3-7 of the 6th data block for ORTHOTROPIC, and and fields 1-5 or the 4e data block for ANISOTROPIC.
For the rest of the Hardening Rules, the input properties are as shown. No WORK HARD or STRAIN RATE options are created with these.
 
Power Law &
Rate Power Law
Description
Coefficient A
1st field of 6th data block of ISOTROPIC option.
Coefficient B
3rd field of 6th data block of ISOTROPIC option.
Exponent M
2nd field of 6th data block of ISOTROPIC option.
Exponent N
4th field of 6th data block of ISOTROPIC option.
Initial Equivalent Strain
5th field of 6th data block of ISOTROPIC option (Power Law). Not used in pre Marc 2005.
Minimum Yield Stress
5th field of 6th data block of ISOTROPIC option (Rate Power Law). Not used in pre Marc 2005.
All the above properties can be temperature dependent if Use Tables is ON and Marc 2005 or later.
Johnson-Cook
Description
Coefficient A
1st field of 8th data block of ISOTROPIC option.
Coefficient B
2nd field of 8th data block of ISOTROPIC option.
Coefficient C
4th field of 8th data block of ISOTROPIC option.
Exponent M
5th field of 8th data block of ISOTROPIC option.
Exponent N
3rd field of 8th data block of ISOTROPIC option.
Initial Strain Rate
8th field of 8th data block of ISOTROPIC option.
Room Temperature
6th field of 8th data block of ISOTROPIC option.
Melt Temperature
7th field of 8th data block of ISOTROPIC option.
Kumar
Description
Coefficient B0
1st field of the 7a data block of ISOTROPIC option.
Coefficient A
2nd field of the 7a data block of ISOTROPIC option. Not necessary if B1-B3 is supplied.
Coefficient B1 - B3
3rd - 5th fields of the 7a data block of ISOTROPIC option. Not necessary if A is supplied.
Coefficient N
1st field of the 7b data block of ISOTROPIC option. Not necessary if B4-B6 is supplied.
Coefficient B4 - B5
2nd - 4th fields of the 7b data block of ISOTROPIC option. Not necessary if N is supplied.
Note:  
Perfectly Plastic is identical to Elastic-Plastic except that no hardening rules apply. Thus no WORK HARD options are created; only ISOTROPIC and STRAIN RATE options with TEMPERATURE EFFECTS, if requested. Stress vs Plastic Strain is replaced with Yield Stress data only as is 10th Cycle Yield vs. Strain replaced with 10th Cycle Yield Stress data. Thus no tabular data is necessary.
Note:  
Rigid-Plastic is identical to Elastic Plastic for Hardening Rules: Power Law, Rate Power Law, Johnson-Cook, and Kumar. Piecewise Linear is identical to Von Mises. The difference here is that the ISOTROPIC option is written and does not contain E
or nu. If an Elastic constitutive model has been created it is ignored, or that is, those values are ignored (elasticity is ignored). A RIGID identifier is placed in the ISOTROPIC option.
Caution:  
In general, you should use true stress vs natural log of plastic strain when defining plasticity curves.
The first value of plastic strain in a stress-strain field must be zero. The corresponding yield stress for this zero plastic strain is placed in the ISOTROPIC option as the Tensile Yield Stress. If yield stress can vary with temperature, the first data point in the field must be the temperature at this yield stress, which will be placed in the TEMPERATURE EFFECTS option, unless you are using the TABLE format, in which the fully defined fields will be converted to equivalent TABLES.
The stress-strain field causes the WORK HARD, DATA option to be written if the first pair of data points of the given field is: (zero, nonzero) This indicates that true stress vs natural log plastic strain data has been supplied. This is consistent with default functionality of Marc. However, if the first data point pair is detected to be (nonzero, nonzero), then this indicates that the engineering stress/strain curve has been given, where the strain is the total strain. Thus the data is converted from engineering stress/strain to true stress/strain before writing the data to the input file. In any case, stress/strain data must begin at the yield stress. In other words, the first pair of data points cannot both be zero. If conversion is necessary, the following formulation is used:
s = Engineering Stress, e = Engineering Strain,
s = True Stress, et = True Total Strain, ee = True Elastic Strain,
ep = True Plastic Strain, E = Young’s Modulus
Plastic - Orthotropic / Anisotropic
This input data can create the ORTHOTROPIC, or ANISOTROPIC, plus WORK HARD, ORTHO TEMP, and STRAIN RATE options. The second data field on the third card of the ORTHOTROPIC or ANISOTROPIC options is set to the corresponding yield criteria.
 
Note:  
All of the Yield Criteria / Hardening Rules have identical inputs as for Isotropic - Plastic materials. The input property values are placed in the equivalent location on the ORTHOTROPIC or ANISOTROPIC options. The only difference is noted here for von Mises yield criteria.
Plastic - von Mises
Description
Stress vs. Plastic Strain
or
Tensile Yield Stress
Same as description for Isotropic Elastic-Plastic - creates WORK HARD, ORTHO TEMP and STRAIN RATE options. Yield Stress is extracted from 1st data point - 1st field of 6th data block of ORTHOTROPIC option or 2nd field on the 4th data block of the ANISOTROPIC option. Temperature field reference creates ORTHO TEMP option. If Strain Rate Method is Piecewise Linear, accepts field of yield stress vs. strain rate and creates STRAIN RATE option with DATA in 2nd field. Data is input in data block 3 for Option B.
Or defines an isotropic yield stress. It is entered in the first data field on the sixth card of the ORTHOTROPIC option and is a required property when the plasticity type is Perfectly Plastic.
Note:  
Perfectly Plastic is identical to Elastic-Plastic except that no hardening rules apply. Thus no WORK HARD options are created. Stress vs Plastic Strain is replaced with Yield Stress data only as is 10th Cycle Yield vs. Strain replaced with 10th Cycle Yield Stress data. Thus no tabular data is necessary.
 
Shape Memory - Isotropic
This input data creates the SHAPE MEMORY keyword option.
Shape Memory
Description
Memory Model
Either a Mechanical (Auricchio’s) model or a Thermal-Mechanical model is written. These are options to the constitutive model. Datablock 3, field 2. Note: Reference temperature values taken from the Elastic constitutive model.
Property Word
Description (Mechanical - Auricchio’s)
Young’s Modulus &
Poisson’s Ratio
These must be defined in an Elastic constitutive model. Thus an Elastic constitutive model must exist in order to write a SHAPE MEMORY option for the Mechanical option. Block 4b, 1st and 2nd fields, respectively.
Sigma AS_s
Block 4b, field 3.
Sigma AS_f
Block 4b, field 4.
Sigma SA_s
Block 4b, field 5.
Sigma SA_f
Block 4b, field 6.
Epsilon L (0.0 ~ 1.0)
Block 5b, field 1.
Alpha (0.0 ~ 0.10)
Block 5b, field 2.
Martensite Slope
Block 5b, field 4.
Austenite Slope
Block 5b, field 5.
Property Word
Description (Thermal-Mechanical)
Young’s Modulus
Poisson’s Ratio
Coefficient of Thermal Expansion
Initial Yield Stress
Mass Density
(Austenite)
Block 4a, fields 1-5, respectively
Young’s Modulus
Poisson’s Ratio
Coefficient of Thermal Expansion
Initial Yield Stress
Mass Density
(Martensite)
Block 5a, fields 1-5, respectively
Martensite Start Temperature
Block 6a, field 1.
Martensite Finish Temperature
Block 6a, field 2.
Martensite Slope
Block 6a, field 3.
Austenite Start Temperature
Block 6a, field 4.
Austenite Finish Temperature
Block 6a, field 5.
Austenite Slope
Block 6a, field 6.
Deviatoric Trans. Strain
Block 7a, field 1.
Volumetric Trans. Strain
Block 7a, field 2.
Twinning Stress
Block 7a, field 3.
Stress Dependency Coefficient g-A
Block 8a, field 1.
Exponent g-B
Block 8a, field 2.
Coefficient g-C
Block 8a, field 3.
Exponent g-D
Block 8a, field 4.
Coefficient g-E
Block 8a, field 5.
Exponent g-F
Block 8a, field 6.
Nondimensionalizign Stress g-O
Block 9a, field 1.
Cut Off Value g-max
Block 9a, field 2.
Stress at g-max
Block 9a, field 3.
Damage - Isotropic / Orthotropic / Anisotropic
Below is the Damage constitutive model and writes the DAMAGE option. This is a constitutive model valid for the types listed above and can reference ISOTROPIC, ORTHOTROPIC, ANISOTROPIC options or one of the Hyperelastic options: MOONEY, OGDEN, GENT, ARRUDA-BOYCE, but not both. So if a Hyperelastic model is active, and the Damage model below is 4,5, or 6, it should reference the Hyperelastic model; if it is 0-3, 9 or 10 it should reference the Isotropic, Orthotropic, or Anisotropic materials.
Damage
Description
Damage Type
Damage Model
For Isotropic, all models are valid. For Orthotropic and Anisotropic only models 0-3 and 9/10 are valid. The given model number is written to the 2nd datablock of the DAMAGE option (the valid property words are indicated):
0 - No Nucleation (1-5)
1 - Strain Controlled Nucleation (1-6,8,9)
2 - Stress Controlled Nucleation (1-5, 7-9)
3 - User Sub UVOIDN (1-5)
4 - Rubber - additive decomposition (10-17, 24)
5 - Rubber - multiplicative decomp. (18-24)
6 - User Sub UELDAM (none)
9 - Simplified Yield - User Sub UDAMAG (none)
10 - Simplified Yield/E - User Sub UDAMAG (none)
1st Yield Surface Multiplier
(1) 1st field, 4a data block of DAMAGE option.
2nd Yield Surface Multiplier
(2) 2nd field, 4a data block
Initial Void Volume Fraction
(3 3rd field, 4a data block)
Critical Void Volume Fraction
(4) 4th field, 4a data block
Failure Void Volume Fraction
(5) 5th field, 4a data block
Mean Strain for Nucleation
(6) 7th field, 4a data block
Mean Stress for Nucleation
(7) 7th field, 4a data block
Standard Deviation
(8) 8th field, 4a data block
Volume Fraction of Void
Nucleation
(9) 9th field, 4a data block
1st Scale Factor - Cont. Damage
(10) 1st field, 4b data block
1st Relax Factor - Cont. Damage
(11) 2nd field, 4b data block
2nd Scale Factor - Cont. Damage
(12) 3rd field, 4b data block
2nd Relax Factor - Cont. Damage
(13) 4th field, 4b data block
1st Scale Factor - Discont.
Damage
(14) 5th field, 4b data block
1st Relax Factor - Discont.
Damage
(15) 6th field, 4b data block
2nd Scale Factor - Discont.
Damage
(16) 7th field, 4a data block
2nd Relax Factor - Discont.
Damage
(17) 8th field, 4a data block
1st Scale Factor
(18) 1st field, 4c data block
1st Proportional Term
(19) 2nd field, 4c data block
1st Relax Rate Constant
(20) 3rd field, 4c data block
2nd Scale Factor
(21) 4thfield, 4c data block
2nd Proportinal Term
(22) 5th field, 4c data block
2nd Relax Rate Constant
(23) 6th field, 4c data block
Scale Factor @ Infinity
(24) 3rd field, 3rd data block
Cracking - Isotropic
Below is the Cracking constitutive model for concrete cracking and writes the CRACK DATA option.
Cracking
Description
Method
Either Entered Values or User Sub. UCRACK... If user subroutine is specified, CRACK DATA may not have to be written - needs investigation.
Critical Stress
1st field, 3rd data block of CRACK DATA
Softening Modulus
2nd field, 3rd data block
Crushing Strain
3rd field, 3rd data block
Shear Retention
4th field, 3rd data block
Forming Limit - Isotropic / Orthotropic / Anisotropic
Below is the Forming Limit constitutive model addition for Isotropic, Orthotropic, and Anisotropic material categories. This writes the FORMING LIMIT option.
Forming Limit
Description
Method
Either Fitted, Predicted, or Table. A zero, 1, or 2 is written to the 1st field of the 2nd data block, respectively.
C0-C1 and D1-D4
Data block 3a and 4a for Option 0 (Method - Fitted)
Strain Hardening Exponent
Thickness Coefficient
Data block 3b for Option 1 (Method - Predicted)
Forming Limit Diagram
Data block 3c of Option 2 (Method - Table). Reference value always 1.0. Must use a TABLE option for this as it must reference a Strain field.
Grain Size - Isotropic
Below is the Grain Size constitutive model for Isotropic model only. This writes the GRAIN SIZE and MATERIAL DATA options.
Grain Size
Description
Method
Either Yada or User Sub UGRAIN. A 1 or -1, respectively, in 2nd field of 3rd data block of GRAIN SIZE option.
Initial Grain Size
Data block 4, 1st field
C1-C5
Data block 4, fields 2-6.
Activation Energy (Q)
This is written to the MATERIAL DATA option (1st field, 4th data block) where the GRAIN SIZE material ID is referenced in the MATERIAL DATA option.
Soil - Isotropic
Below is the Soil constitutive model addition for Isotropic and Orthotropic models only. This writes the SOIL option and if necessary, the INITIAL POROSITY, INITIAL VOID RATIO, INITIAL PC and SPECIFIC WEIGHT options.
Soil
Description
Model
Either Linear, Cam Clay, or User Sub. HYPELA. This is indicated in the 2nd field of the 3rd data block by entering LINEAR, NON LINEAR (user sub. HYPELA) or CAMCLAY. If a Plastic model is also defined, this overrides this option and the Plastic model setting will write either VON MISES, LIN MOHRC, or PLB MOHRC for von Mises, Linear Mohr-Coulomb or Parabolic Mohr-Coulomb yield models. For orthotropic models, the ORTHOTROPIC keyword is written.
Dynamic Viscosity
Data block 4, 8th field
Fluid Density
Data block 4, 7th field
Fluid Bulk Modulus
Data block 4, 7th field
Permeability
Data block 5, 1st field
Compression Ratio
Data block 5, 2nd field
Recompression Ratio
Data block 5, 3rd field
Critical State Curve Slope
Data block 5, 4th field
Young’s Modulus
Poisson’s Ratio
Mass Density
Coefficient of Thermal Expansion
These values get placed in the 1st-4th fields of datablock 4. If any of these values reference a temperature field, the TEMPERATURE EFFECTS is written (or TABLES if Use Tables is ON). Or for Orthotropic properties, they are placed in the 4th, 5th, and 6th datablocks.
Yield Stress
This value comed from a Plastic constitutive model. If this model is not available, then zero is written for the Yield Stress. If a Perfectly Plastic model is available, the Yield Stress is placed in the 5th field of the 4th datablock. If a stress-strain field is available, then the WORK HARD option is written (or TABLE) with this value being the reference value at zero plastic strain.
Initial Porosity
Initial Void Ratio
Initial Preconsolidation Pressure
Gravity Constants in 1st-3rd coordinate directions
These properties are written to the INITIAL POROSITY, INITIAL VOID RATIO, INITIAL PC, and SPECIFIC WEIGHT options, respectively and are assigned to the same elements as this material.
Powder - Isotropic
Below is the Powder constitutive model for Isotropic model only. This writes the POWDER, RELATIVE DENSITY, and DENSITY EFFECTS options.
Powder
Description
Method
Either Entered Values or User Sub. UPOWDR. If the latter is seletect, then no POWDER option (or RELATIVE DENSITY, DENSITY EFFECTS) options are written. Everything is taken care of in the UPOWDR routine supposedly.
Material Prop. Gama
Data block 4, 6th field
Material Prop. Beta
Data block 4, 7th field
Powder Viscosity
Data block 4, 8th field
Gamma Coef. 1-4
Data block 6
Beta Coef. 1-4
Data block 7
Initial Relative Density
This goes on the RELATIVE DENSITY option. Note that for shell elements, the integration points have to be written also.
Young’s Modulus
Poisson’s Ratio
Mass Density
Coefficient of Thermal Expansion
These come from an Elastic constitutive model, which must be defined also in addition to the Powder model. These values get placed in the 1st-4th fields of datablock 4. If any of these values reference a temperature field, the TEMPERATURE EFFECTS is written (or TABLES if Use Tables is ON). If the first two (or last two for Coupled analysis) reference a Strain field, then the DENSITY EFFECTS, DATA option is written with the density effects field written to the appropriate block of the option. This is written in an identical way to the TEMPERATURE EFFECTS, DATA option. We are using the Strain field to indicate a Density field in this case since Density fields are not yet supported in Patran Fields application. Of course if Use Tables is ON, then TABLES are used and not TEMP/DENSITY EFFECTS.
Yield Stress
This value comed from a Plastic constitutive model. If this model is not available, then zero is written for the Yield Stress. If a Perfectly Plastic model is available, the Yield Stress is placed in the 5th field of the 4th datablock. If a stress-strain field is available, then the WORK HARD option is written (or TABLE) with this value being the reference value at zero plastic strain.
Electrostatic - Isotropic/Orthotropic
Below is the Electrostatic constitutive model for Isotropic and Orthotropic models only. This writes the ISOTROPIC, ELECTROSTA or ORTHOTROPIC, ELECTROSTA options, respectively.
 
Powder
Description
Permittivity, Permittivity 11/22/33
Values written to the above mention options.
Electrodynamic - Isotropic/Orthotropic/Anisotropic
Below is the Electrodynamic constitutive model for Isotropic, Orthotropic, and Anisotropic models. This writes the ISOTROPIC, THERMAL or ORTHOTROPIC, THERMAL options, respectively.
 
Powder
Description
Resistivity,
Resistivity 11/12/13/22/23/33
Values written to the above mention options.
Magnetostatic - Isotropic/Orthotropic
Below is the Magnetostic constitutive model for Isotropic and Orthotropic models. This writes the ISOTROPIC or ORTHOTROPIC options, respectively for magnetostatics.
 
Powder
Description
Permeability,
Permeability 11/22/33
Inverse Permeability,
Inverse Permeability 11/22/33
Values written to the above mention options.
Hn-Bn / Bn-Hn Curve
These curves are defined under the Field application using a Magnetic material field.
Piezoelectric - Isotropic/Orthotropic/Anisotropic
Below is the Piezoelectric constitutive model for Isotropic, Orthotropic, and Anisotropic models. This writes the ISOTROPIC or ORTHOTROPIC or ANISOTROPIC options, respectively for piezoelectic
 
Powder
Description
Piezoelectric Constants
Electric Permitivity 11/22/33
Values written to the above mention options.
 
Composite - Homogeneous
The following composite material types may also be defined as shown in this table.
The Composite forms are used to create new materials by combining existing materials. All of the composite materials, with the exception of the laminated composites, can be assigned to elements, as
any homogeneous material, through the element property forms. For the laminated composites, the section thickness is entered indirectly through the definition of the stack, and the Homogeneous option, on the Element Properties for shells, plates and beam, must be changed to Laminate to avoid reentry of this information.
For details on entering data on the Composite forms, refer to the Composite Materials Construction (p. 110) in the Patran Reference Manual.
For all composite types except Composite - Laminate, an equivalent set of properties are entered in the ANISOTROPIC keyword option when an Marc input file is created. For Composite - Laminate the COMPOSITE option is used.
 
Caution:  
It is extremely important that when you define a layup (in the form on the next page), that it be done from top to bottom. Think of the top layer of the layup as being the top row of the spreadsheet and you should have no problems. As an example of how important this is, consider a cantilevered flat plate subject to an axial load with two layers. The top layer is extremely flexible compared to the bottom layer, which is relatively much stiffer than the top. Due to the shear forces created between the layers, the vertical deflection should tend to favor the side of the stiffer layer, thus the plate should bend down. If the layer is defined from bottom to top instead of top to bottom, you will get what appears to be the opposite answer where the deflection bends up. The answers are correct in both cases. The problem is how you defined the layup.
Composite - Laminate
This form appears when Composite is the selected Object and Laminate is the selected Method in the Materials application. Use this form to create the COMPOSITE keyword option.
The difference between the "Total" option and the "Total - %thicknesses" option is that the former requires that the user give actual thickness values of each ply and the latter requires each ply thickness to be given as a percentage of the total layup thickness. This is the prefered method when applying the composite material to solid (CHEXA) elements or 2D solid element (axisymmetric, plane strain).
The 3rd field of the 3rd datablock of the COMPOSITE option is either 0 or 1, respectively for the same options.
 
Caution:  
See the caution on the previous page. Layers must be defined from top to bottom.
Constitutive Model Status
A single material may contain multiple constitutive models. The constitutive model used is determined by the Constitutive Model Status. Patran will use all constitutive models active when the analysis is submitted. Redundant or unneeded constitutive models should be rendered inactive.
 
Note:  
The modifications are not saved until Apply button is pressed.
Experimental Data Fitting
This is a very useful tool available under the Tools pull-down menu from the main Patran form and is only available if the Analysis Preference is set to Marc.
The tool is used to curve fit experimentally derived raw elastomeric material data and fit a number of material models to the data. This data can then be saved as constitutive hyperelastic and/or viscoelastic models for use in an Marc analysis. The operation of curve fitting is done in three basic steps corresponding to the actions in the Action pull-down menu.
1. Import Raw Data - data is read from standard ASCII files and stored in Patran in the form of a field (table).
2. Select Test Data - the fields from the raw data are associated to a test type.
3. Calculate Properties - the curve fit is done to the selected test data; coefficients are calculated based on the selected material model; curve fit is graphically displayed and the properties can be saved as a constitutive model for a later analysis.
 
Note:  
Strain input should be engineering strain to give reasonable results.
The Ogden Formulation was first given in the paper "Large Deformation Isotropic Elasticity - on the Correlation of Theory and Experiment for Incompressible Rubberlike Solids", R.W. Ogden, Proc.R.Soc.Lond.A., Vol. 326, 526-584 (1972). The curve fitting determines ( mu_n, alpha_n ) pairs. These constants are material constants and may not represent physical values for rubbers since during the curve fitting process, certain calculations are made with the assumption of imcompressibility. The most important issue during data fitting is to make sure that the data fit is sufficiently close.
The Foam Model (see - Storåkers, B., On Material Representation and Constitutive Branching in Finite Compressible Elasticity, Journal of the Mechanics and Physics of Solids, vol.34, no.2, pp. 125-145, 1986.) is a compressible Ogden formulation and should be used for materials going through large volumetric deformations. The curve fitting calculates sets of ( mu_n, alpha_n, beta_n ) coefficients where the Beta coefficients represent to some extent a measure of foam compressibility. The Planar (Pure) Shear and Simple Shear responses are identical to the Ogden Formulation since the motion is isochoric; therefore, use of either Pure or Simple shear experiments to determine the Beta coefficients is pointless. The model works well in compression (densification).
When using the foam model, note that like the Ogden formulation, it is acceptable to get different parameters for the fit as long as the fit is correct and the also yields a positive definite strain energy function for the range of the fit. (A positive definite strain energy function means that the material matrix derived from it will not have a negative Jacobian through the range of deformation). If a negative Jacobian occurs during the analyis, this may cause an exit 1005 or 1009 which signifies "inside-out elements".
The beta coefficients (which represent some measure of compressibility) may vary since there are more than one way to handle the strain energy attributed to the volumetric deformation. For the foam model, compressibility (in the form of fictive poisson's ratio) is included and in the test data, the independent stretch and volume ratios would need to be considered.
Finally, it is highly recommended that mathematical checks be used for all data fitting, especially for the Ogden and Foam formulations.
 
Import Raw Data
You can import the raw materials data by following these general steps:
Keep in mind the following points and considerations when importing raw data:
1. You can skip any number of header lines in the raw data file by setting the Header Lines to Skip data box.
2. You may edit the raw data file after selecting it by using the Edit File... button. The editor is Notepad on Windows platforms and vi on UNIX platforms unless you change the environment variable P3_EDITOR to reference a different editor. The editor must be in the user’s path or the entire pathname must be referenced.
3. Raw data files may have up to three columns of data. By default the first column of data is the independent variable value. The second column is the measured data, and the last column can be the area reduction or volumetric data. More than three columns is not accepted. If the third column is blank, the material is considered incompressible.
4. If you have cross-sectional area reduction data in the third column, you can give it an optional field name also by turning ON the Area Data toggle and supplying an Area Field Name. If you have three columns of data and this toggle is OFF, the third column is still detected and read and two fields are created. This results in a _C1 and _C2 being appended to the New Field name.
5. The data may be space, tab, or comma delimited.
6. If for some reason the independent and dependent columns need to be interchanged, you can turn the Switch Ind./Dep. Columns toggle ON. Check your imported fields before proceeding to ensure they are correct. This is done in the Fields application.
7. When you press the Apply button, you will be taken to the second step. If you need to import more than one file, you will have to reset the Action pull-down.
Select Test Data
Once raw test data is imported, you must associate them with particular test types or modes by following these steps:
Keep in mind the following points and considerations when selecting test data:
1. Typical stress-strain data for Deformation Mode tests are referenced in the Primary column. If you have volumetric data, these are entered in the Secondary column databoxes and are optional.
2. For Viscoelastic (time relaxation data), you must turn ON the ViscoElastic toggle. Only viscoelastic curve fitting will be done in this case. To return to Deformation Mode, turn this toggle OFF.
3. Damage models are not yet supported.
4. When you press the Apply button, you will be taken to the third step.
Calculate Properties
Once test data has been associated to a test type or mode the curve fit function is performed by following these steps:
Keep in mind the following points and considerations when calculating properties:
1. The plots are appended to the existing XY Window until you press the Unpost Plot button. You can turn the Append function ON/OFF under the Plot Parameters... form.
2. By default, all the deformation modes are plotted along with the raw data even if raw data has not been supplied for those modes. This is very important. These additional modes are predicted for you. You should always know your model’s response to each mode of deformation due to the different types of stress states. For example, a rule of thumb for natural rubber and some other elastomers is that the tensile tension biaxial response should be about 1.5 to 2.5 times the uniaxial tension response.
3. You can turn ON/OFF these additional modes or any of the curves under the Plot Parameters button as well as change the appearance of plot. More control and formatting of the plot can be done under the XY Plot application on the Patran application switch on the main form.
4. Viscoelastic constitutive models are useless without a Hyperelastic constitutive model also. Be sure your model has both defined under the same material name if you use viscoelastic properties.
5. You may actually change the coefficient values in the Coefficients spread sheet if you wish to see the effect they have on the curve fit. Select one of the cells with the coefficient you wish to change, then type in a new coefficient value in the Coefficient Value data box and press the Return or Enter key. Then press the Plot button again. If you press the Apply button, the new values will be saved in the supplied material name.
6. For viscoelastic relaxation data, the Number of Terms used in the data fit should, as a rule of thumb, be as many as there are decades of data.
7. A number of Optional and Plot Parameters are available to message the data and control the curve fitting. See the table below for more detailed descriptions.
 
Optional Parameters
Description
Uniaxial Test
Biaxial Test
Planar Shear Test
Only available for Ogden and Foam models. Defines whether area or volumetric data was measured.
Mathematical Checks
OFF by default. Only available for Ogden and Foam models.
Positive Coefficients
OFF by default. Will force positive coefficients to be determined if ON. Available for all Model types.
Extrapolate
Left/Right Bounds
OFF by default. If ON, the Left and Right Bounds databoxes will become available to enter data to extrapolate results to. Available for all Model types.
Error
Can be set to Relative (default) or Absolute. Good for all Model types.
Error Limit
Only available for Ogden, Foam, Arruda-Boyce, and
Gent Models.
# of Iterations
Only available for Ogden, Foam, Arruda-Boyce, and
Gent Models.
Convergence Tolerance
Only available for Ogden, Foam, Arruda-Boyce, and Gent Models. This can have a significant difference in the calculated coefficients and the plots.
Use Fictive Coefficient
Fictive Coeff.
Only valid for Foam. Allows you to enter a fictive Poison’s ratio for use in the data fit.
Append Curves
Curves will be appended to existing plot. If OFF, plot will be cleared each time.
X/Y Axis Options
Plot data in linear or logarithmic fashion.
Modes
Turns ON/OFF each respective mode including the raw
data plot.