Difference between revisions of "Common Material Properties"

From OSUPDOCS
Jump to navigation Jump to search
 
(154 intermediate revisions by the same user not shown)
Line 2: Line 2:


== Basic Properties ==
== Basic Properties ==
These are basic material properties.
{| class="wikitable"
|-
! Property !! Description !! Units !! Default
|-
| rho || The material's initial density. || [[ConsistentUnits Command#Legacy and Consistent Units|density units]] || 0.001
|-
| Cv ||  The constant-volume heat capacity. It is used when doing [[Thermal Calculations#Conduction|conduction calculations and/or coupled mechanical energy]] and by some material constitutitive laws. You do not need to enter a constant pressure heat capacity (Cp). Instead, when [[Thermal Calculations|MPM thermodynamics]] are done correctly, the simulations will automatically account for all boundary condition effects (constant volume, constant pressure, or mixed) on effective heat capacity. || [[ConsistentUnits Command#Legacy and Consistent Units|heat capacity units]] || 1
|-
| kCond || Thermal conductivity for isotropic materials ([[Material Models|anisotropic materials]] will have alternate properties for setting thermal conductivity tensor). || [[ConsistentUnits Command#Legacy and Consistent Units|conductivity units]] || 0
|-
| csat || The saturation concentration potential as a weight friction from 0 to 1. it is only used when doing [[Diffusion Calculations|diffusion]] calculations. || none || 1
|-
| beta || Moisture expansion coefficient for isotropic materials ([[Material Models|anisotropic materials]] will have alternate properties for setting moisture expansion tensor). || strain/(wt fraction) || 0
|-
| D || Solvent diffusion constant for isotropic materials ([[Material Models|anisotropic materials]] will have alternate properties for setting diffusion tensor). || [[ConsistentUnits Command#Legacy and Consistent Units|diffusion units]] || 0
|-
| largeRotation || To use (1) or not use (0) polar decomposition when calculating rotations in small-stain materials. This option applies only to [[Material Models#Linear Elastic Small Strain Materials|linear elastic small strain materials]] and to [[Material Models#Elastic-Plastic Small Strain Materials|elastic-plastic small strain materials]]. Option 1 finds [[MPM Methods and Simulation Timing#Incremental Deformation Gradient|incremental deformation gradient]] use the selected number of terms and evaluates the rotational part of that increment using polar decomposition. The stress update is then rotated according to the decomposed rotation. In contrast, option 0 finds [[MPM Methods and Simulation Timing#Incremental Deformation Gradient|incremental deformation gradient]] from a linear expansion (<tt>k<sub>max</sub></tt>=1) and evaluates the rotation component using second order (in 2D) or first order (in 3D) approximate polar decomposition. The stress update is rotated by standard hypoelastic methods. Both methods can handle large total rotations (provided they are incrementally small). Option 1 may be more accurate, but it is less efficient. || none || 0
|-
| matDamping || Sets custom particle damping to apply only to particles of this material. Its value replaces global [[Damping Options#Grid and Particle Damping|particle damping setting]] and must be a constant (function of time not allowed). In XML files, set matDamping using <tt><PDamping>matDamping</PDamping></tt> || 1/[[ConsistentUnits Command#Legacy and Consistent Units|time units]] || none
|-
| color || Sets the color of the material. The color is used in material point method plots material type in [[NairnFEAMPM]] and in [[NairnFEAMPMViz]]. If no color is provided, a color will be picked from the current spectrum using the material number. In scripted files, this property takes four arguments being red, green, blue, and alpha values between 0.0 and 1.0. A single argument means to set gray level between 0.0 and 1.0 (with alpha=1.0). Three arguments means set red, green, and blue with alpha=1.0. In <tt>XML</tt> files, the color is set with "red", "green", "blue" and "alpha" attributes (and the element's content is ignored). || none || none
|}


== Fracture Toughness Properties ==
== Fracture Toughness Properties ==
Line 11: Line 36:
! Property !! Description !! Units !! Default
! Property !! Description !! Units !! Default
|-
|-
| JIc || Critical energy release rate fracture toughness for mode I. It is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criteria 2, 3, or 7]]. For criterion 2, it is only used if <tt>initTime</tt> is not specified. It is also used to set toughness of traction law materials.|| J/m<sup>2</sup> || none
| JIc || Critical energy release rate fracture toughness for mode I. It is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criteria 2, 3, or 7]]. For criterion 2, it is only used if <tt>initTime</tt> is not specified. It is also used to set toughness of traction law materials.|| [[ConsistentUnits Command#Legacy and Consistent Units|energy release units]] || none
|-
|-
| JIIc || Critical energy release rate fracture toughness for mode II. It is currently only used to set toughness of traction law materials. || J/m<sup>2</sup> || none
| JIIc || Critical energy release rate fracture toughness for mode II. It is currently only used to set toughness of traction law materials. || [[ConsistentUnits Command#Legacy and Consistent Units|energy release units]] || none
|-
|-
| KIc || Critical mode I stress intensity factor. It is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criteria 1, 4, or 5]]. || MPa m<sup>1/2</sup> || none
| KIc || Critical mode I stress intensity factor. It is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criteria 1, 4, or 5]]. || [[ConsistentUnits Command#Legacy and Consistent Units|stress intensity units]] || none
|-
|-
| KIIc || Critical mode II stress intensity factor. It is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criteria 1, 4, or 5]]. || MPa m<sup>1/2</sup> || none
| KIIc || Critical mode II stress intensity factor. It is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criteria 1, 4, or 5]]. || [[ConsistentUnits Command#Legacy and Consistent Units|stress intensity units]] || none
|-
|-
| KIexp || Exponent p in the elliptical criteria for crack growth. It is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 5]]. || none || 2
| KIexp || Exponent p in the elliptical criteria for crack growth. It is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 5]]. || none || 2
Line 23: Line 48:
| KIIexp || Exponent q in the elliptical criteria for crack growth. It is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 5]]. || none || 2
| KIIexp || Exponent q in the elliptical criteria for crack growth. It is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 5]]. || none || 2
|-
|-
| delIc || Critical crack opening displacement for mode I. Only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 6]]. It is also used by traction-law materials. || mm || none
| delIc || Critical crack opening displacement for mode I. Only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 6]]. It is also used by traction-law materials. || [[ConsistentUnits Command#Legacy and Consistent Units|length units]] || none
|-
| delIIc || Critical crack opening displacement for mode II. Only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 6]]. It is also used by traction-law materials. || [[ConsistentUnits Command#Legacy and Consistent Units|length units]] || none
|-
|-
| delIIc || Critical crack opening displacement for mode II. Only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 6]]. It is also used by traction-law materials. || mm || none
| initTime || The time when crack propagation starts. It is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 2]]. For criterion 2, when <tt>initTime</tt> is specified, takes precedence over the <tt>JIc</tt> property. || [[ConsistentUnits Command#Legacy and Consistent Units|alt time units]] || none
|-
|-
| initTime || The time when crack propagation starts. It is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 2]]. For criterion 2, when <tt>initTime</tt> is specified, takes precedence over the <tt>JIc</tt> property. || ms || none
| speed || The crack speed in steady state crack propagation. This speed, however, is only active for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 2]]. (also used in [[Crack Propagation Commands#Crack Propagation Criterion|criterion 3]] as an initial crack speed, but that criterion is not meant for general use) || [[ConsistentUnits Command#Legacy and Consistent Units|alt velocity units]] || 1
|-
|-
| speed || The crack speed in steady state crack propagation. This speed, however, is only active for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 2]]. (also used in [[Crack Propagation Commands#Crack Propagation Criterion|criterion 3]] as an initial crack speed, but that criterion is not meant for general use) || m/sec || 1
| maxLength || The maximum crack length for steady state crack propagation. The simulation will stop soon after crack reaches the input length. This length, however, is only active for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 2]]. || [[ConsistentUnits Command#Legacy and Consistent Units|length units]] || none
|-
|-
| nmix || An exponent used in mixed-modes failure of some traction laws. || none || 1
| nmix || An exponent used in mixed-modes failure of some traction laws. || none || 1
|-
|-
| gamma || Crack surface energy. It is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 3]], which is in development and not meant for general use. || J/m<sup>2</sup> || none
| allowsCracks || Set to 1 to allow the material to have cracks or 0 to not allow cracks. The default is 1 for non-rigid materials and 0 for rigid materials. Setting a non-rigid material to 0 lets that material enter a crack in another material such as to simulate cutting, to wedge open a crack, or model hydraulic fracturing. The non-cracking material can open the crack in the other material by contact on its crack surfaces. || none || varies
|-
| p || Irreversible energy term in energy balance crack growth and probably should be between 0 and 1. It is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 3]], which is in development and not meant for general use. || none || 1
|-
| gain || Gain used in feedback loop to maintain energy balance during crack growth. It is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 3]], which is in development and not meant for general use. || 1e5 || 1
|}
|}


Line 48: Line 71:
! Property !! Description !! Units !! Default
! Property !! Description !! Units !! Default
|-
|-
| criterion || To set a custom crack propagation criterion. You can use any [[Crack Propagation Commands#Crack Propagation Criterion|propagation option]], but must enter using the numeric setting. Setting this property in a material overrides the [[Crack Propagation Commands#Propagate Command|default propagation criterion setting]]. || none || none
| criterion || To set a custom crack propagation criterion. You can use any [[Crack Propagation Commands#Crack Propagation Criterion|propagation option]]. Setting this property in a material overrides the [[Crack Propagation Commands#Propagate Command|default propagation criterion setting]]. This command is for scripted files only; see [[#Propagation Properties in XML Files|below]] to set <tt>criterion</tt> in <tt>XML</tt> files. || none || none
|-
|-
| direction || To set a custom crack propagation direction. You can use any [[Crack Propagation Commands#Crack Propagation Direction|direction option]], but must enter using the numeric setting. Setting this property in a material overrides the [[Crack Propagation Commands#Propagate Command|default propagation direction setting]]. || none || none
| direction || To set a custom crack propagation direction. You can use any [[Crack Propagation Commands#Crack Propagation Direction|direction option]]. Setting this property in a material overrides the [[Crack Propagation Commands#Propagate Command|default propagation direction setting]]. This command is for scripted files only; see [[#Propagation Properties in XML Files|below]] to set <tt>direction</tt> in <tt>XML</tt> files.|| none || none
|-
|-
| traction || To set a custom [[Crack Propagation Commands#Traction Law in Wake of Propagation|traction law to create for crack propagation]] in this material. A traction law set in a material overrides the [[Crack Propagation Commands#Propagate Command|default traction law]]. The traction law must be set by number (1, 2, <i>etc.</i>) where the number corresponds to the material number, in order, as defined by <tt>Material</tt> commands. You cannot use a material ID. Numerical values are needed because you may be using this option before the traction law material has been created. || none || none
| traction || To set a custom [[Crack Propagation Commands#Traction Law in Wake of Propagation|traction law to create for crack propagation]] in this material. A traction law set in a material overrides the [[Crack Propagation Commands#Propagate Command|default traction law]]. The traction law can be set by material ID (if the traction law has already been defined) or by number (if it is not defined yet). This command is for scripted files only; see [[#Propagation Properties in XML Files|below]] to set <tt>traction</tt> in <tt>XML</tt> files. || none || none
|-
| altcriterion || Same as "criterion" property above except that it applies to the [[Crack Propagation Commands#Alternate Propagation Criterion|alternate propagation criterion]] for the material || none || none
|-
| altdirection || Same as "direction" property above except that it applies to the [[Crack Propagation Commands#Alternate Propagation Criterion|alternate propagation criterion]] for the material || none || none
|-
| alttraction || Same as "traction" property above except that it applies to the [[Crack Propagation Commands#Alternate Propagation Criterion|alternate propagation criterion]] for the material || none || none
|-
|-
| xGrow || This property along with <tt>yGrow</tt> (if on one given the other is set to 0) specify a unit vector for a constant crack growth direction. t is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 2]] and then only if that criterion is using its [[Crack Propagation Commands#Crack Propagation Direction|default propagation direction]]. The result is a constant crack growth direction regardless of stress state or crack tip orientation. Any input vector will be normalized to a unit vector. If a constant crack growth direction with a fixed crack is located precisely on grid lines, it is possible the crack algorithm will not recognize the crack plane. Is it better to move such a crack slightly off grid lines. || none || none
| xGrow || This property along with <tt>yGrow</tt> (if on one given the other is set to 0) specify a unit vector for a constant crack growth direction. t is only used for crack propagation by [[Crack Propagation Commands#Crack Propagation Criterion|criterion 2]] and then only if that criterion is using its [[Crack Propagation Commands#Crack Propagation Direction|default propagation direction]]. The result is a constant crack growth direction regardless of stress state or crack tip orientation. Any input vector will be normalized to a unit vector. If a constant crack growth direction with a fixed crack is located precisely on grid lines, it is possible the crack algorithm will not recognize the crack plane. Is it better to move such a crack slightly off grid lines. || none || none
|-
|-
| yGrow || Crack growth direction - see <tt>xGrow</tt> above. || none || none
| yGrow || Crack growth direction - see <tt>xGrow</tt> above. || none || none
|-
| constantTip || Set to 0 or 1.  The default of 0 means the crack tip will track the [[Crack Propagation Commands#Crack Tip Materials|material around the crack tip]]. Changing it to 1 means crack tips with this material will always use this material even if the crack propagates into another material. The default 0 allows modeling crack growth in composites with fracture properties changing as cracks move between materials. Using 1 allows modeling multiple cracks in the same material having different fracture propertie by using the following steps:
<ol>
<li>Define multiple materials that are identical except for their fracture and/or propagation properties and all with constantTip=1.
<li>Create a model and uses one of the materials  for all material points.
<li>For each crack tip, assign its material to be the material in step 1 with the appropriate fracture properties.
</ol>
|| none || 0
|}
|}


== Script Notes ==
=== Propagation Properties in XML Files ===


<p>The following material properties are common to all types of materials, although not used by all materials. These properties are only used in MPM analyses.
In <tt>XML</tt> files, the criterion, direction, and traction properties (and the analogous alternate propagation properties) are set differently. To set crack propagation criteria, you use instead
</p>


<ul>
<Propagate criterion='(critNum)' direction='(dirNum)' traction='(traction)'/>
<AltPropagate criterion='(critNum)' direction='(dirNum)' traction='(traction)'/>


<li><code>rho</code> for density in g/cm^3.</li>
where the settings are the same as defined in the [[Crack Propagation Commands#Propagate Command|default crack propagation commands]] (or the [[Crack Propagation Commands#Alternate Propagation Criterion|alternate propagation command]]), but the <tt>XML</tt> element is now used within a <tt><Material></tt> definition instead of within  the [[Crack Settings|<tt><Cracks></tt> element in the <tt><MPMHeader></tt>]].


<li><code>Cp</code> or <code>Cv</code> for heat capacity (in J/(kg-K)). Heat capacity is used when doing <a href="conduction.html">conduction calculations and/or coupled mechancial energy</a> adn by some material models. Both <code>Cp</code> and <code>Cv</code> set the same (and only) heat capacity used in the MPM calculations. MPM is a strain-drive method, which implies the heat capacity should be interpreted as the constant-volume heat capacity.</li>
<tt>XML</tt> files set xGrow, yGrow, and constantTip as ordinary properties and they function as described [[#Crack Propagation Properties|above]].


<li><a name="csat"></a><code>csat</code> The saturation concentration (in weight fraction) for diffusion calculations.</li>
== Contact Properties ==


<li><code>maxLength</code> for maximum crack length (in mm) for steady state crack propagation. The simulation will stop soon after crack reaches the input length. This limit, however, is only active when there is crack propagation with <a href="propagate.html">criterion 2</a>.</li>
These properties can set custom the contact mechanics between two specific materials. If these properties are not used, material-to-material contact will be determined by the global setting made using a <tt>ContactMM</tt> command (for scripted files) or a <tt>&lt;Friction&gt;</tt> command within the [[Multimaterial MPM#Multimaterial Mode Input Commands|<tt><MultimaterialMode></tt> element]] in the <tt><MPMHeader></tt> (for XML files). For more details see these commands under [[Friction#Friction in Multimaterial MPM|friction settings]] and [[Imperfect Interfaces#Imperfect Interfaces in Multimaterial MPM|interface setting]].


<li><code>altcriterion</code> is identical to <code>criterion</code>, but applies instead to the alternate crack growth criterion in the <a href="altpropagate.html">AltPropagate Command</a>.</li>
{| class="wikitable"
|-
! Property !! Description !! Units !! Default
|-
| Contact<br><Friction> || In scripted files, use the Contact property with two arguments — <tt>(lawID)</tt> and <tt>(matID)</tt>. It sets a custom [[Contact Laws|contact law]] (by its material ID in <tt>(lawID)</tt>) to use for contact between this material and the material defined by <tt>(matID)</tt>. You only need to enter a Contact property for one of the two materials in each pair. Since you have to provide the other material's ID, you should use the Contact command in the second material after the first one has already been defined (and after the [[Contact Laws|contact law]] has been defined as well).<br>
In XML files use the command:<br> &nbsp;&nbsp;&nbsp;<tt>&lt;Friction law='(lawnum)' lawname='(lawID)' mat='(matnum)' matname='(matID)'/&gt;</tt><br>to set the custom [[Contact Laws|contact law]] by number or ID for both the law and the other material. If both number and ID are used, the ID will take precedence.  || none || none
|-
| shareMatField || The value selects another material to share the same velocity field such that the two materials move together and interact by perfect contact (see below for how to have more than two materials share one field). Although materials in the same velocity field interact by perfect contact, they can interact with materials in different velocity fields by [[Multimaterial MPM|multimaterial mode]] contact or interface laws. Finally, materials sharing velocity fields must use "compatible" fields  (''e.g.'', all rigid or all nonrigid); an error  will occur (with an explanation) if you try to share velocity fields of incompatible materials.  || none || none
|}


<li><code>altdirection</code> is identical to <code>direction</code>, but applies instead to the alternate crack growth direction in the <a href="altpropagate.html">AltPropagate Command</a>.</li>
In scripted files, the other material is specified by its material ID, which means the Friction and Interface commands must be used in the secondly-defined materials (such that material ID for the first material is available). In <tt>XML</tt> files, the second material in Friction and Interface  commands is defined by [[Material Command Block|number (or by name)]] using a <tt>mat</tt> or <tt>matname</tt> attribute. You only need a Friction or an Interface command in one material for each pair of materials with custom contact properties.


<li><code>alttraction</code> is identical to <code>traction</code>, but applies instead to the alternate crack growth traction law material in the <a href="altpropagate.html">AltPropagate Command</a>.</li>
To share velocity fields, you first create one "base material" that is not shared and then any number that share that field with a shareMatField property. The base material is specified in the shareMatField property by material ID in scripted files, but must be by number (as the value of the property command)  in <tt>XML</tt> files. To specify custom friction or interfaces between shared materials and other materials, you use any material in the shared block; if you use custom commands for more then one material in a block, only the last one will be used.


<li><code>xGrow</code> and <code>yGrow</code>: These properties (one or both with default 0) specify a unit vector for a constant crack growth direction. They only apply for crack propagation with <a href="propagate.html">criterion 2</a> and then only when that criterion selects the default propagation direction. The result is a constant crack growth direction regardless of stress state or crack tip orientation. Any input vector will be normalized to a unit vector. If constant crack growth directions with a fixed crack are located precisely on grid lines, it is possible for the crack algorithm to not recognize the crack plane. Is it better to move such a crack slightly off grid lines.
=== Deprecated Contact Properties ===
</li>


<li><a name="artificialvisc"></a><code>ArtificialVisc</code> can be <code>on</code> or <code>off</code> to turn artificial viscosity on. When it is on, it adds a pressure, Q, related to velocity gradient on the particle, but only when it is compressing. The equation is
The following material properties for setting contact properties are deprecated. They should be replaced by the <tt>Contact</tt> property above. Note that in XML files, contact properties are still with with a <tt>&lt;Friction&gt;</tt> command but prior attributes are deprecated to be replaced by attributes to specify the desired [[Contact Laws|contact law]].
<blockquote>
Q = &Delta;x*|D<sub>kk</sub>|*(A1*C + A2*&Delta;x*|D<sub>kk</sub>|)
</blockquote>
where &Delta;x is the cell size of the mesh, |D<sub>kk</sub>| is the relative volume change rate, C is the bulk wave speed in the material, and A1 and A2 are adjustable constants, which are entered using the above <code>avA1</code> and <code>avA2</code> properties. The defaults are artificial viscosity off, A1=0.2, and A2=2.0. It is common, but not required, to choose A2 = 10*A1.<br>
&nbsp;&nbsp;&nbsp;&nbsp;The artificial viscosity property is supported in some materials (as documented under each material). If you use these commands in a material that does not support, an error will result.
</li>


<li><code>avA1</code> and <code>avA2</code> to set the coefficients when artficial viscosity is on.</li>
{| class="wikitable"
|-
! Property !! Description !! Units !! Default
|-
| Friction || A Friction property within a material definition can define custom frictional properties for [[Multimaterial MPM|multimaterial mode MPM]] contact between the current material and another material. This property takes two parameters; the first is the same as for the standard [[Friction|Friction command]] and the second gives the other material. || none || none
|-
| Interface || An Interface property within a material definition can define custom imperfect interface parameters properties for [[Multimaterial MPM|multimaterial mode MPM]] contact between the current material and another material. This property takes four parameters; the first three are the same as for a standard [[Imperfect Interfaces#Imperfect Interfaces in MPM|ImperfectInterface command]] (which is actually a <tt><Friction></tt> element in <tt>XML</tt> files) and the fourth gives the other material.|| none || none
|}


<li><code>color 0.1,0.5,0.2,1.0</code>:
== Artificial Viscosity ==
Sets the color of the material with the four arguments being red, green, blue, and alpha values between 0.0 and 1.0. A single argument can be used instead to set gray level between 0.0 and 1.0 (with alpha=1.0). Three arguments set red, green, and blue with alpha=1.0. The color is used in <a href="../movieplots.html">material point method plots</a> of <a href="../quantities.html">material type</a>. If no color is provided, a color will be picked from the current spectrum using the material number.
</li>


<li><code>Friction</code>: A <a href="friction.html">Friction command</a> within a material definition can define custom frictional properties for <a href="multimaterialmode.html">multimaterial mode MPM</a> contact between the current material and another material.</li>
Some materials support artificial viscosity to dampen pressure waves. When it is on, it adds a pressure, <i>Q</i>, related to velocity gradient on the particle, but only when it is compressing. The equation is


<li><code>Interface</code>: An <a href="imperfectinterface.html">Interface command</a> within a material definition can define custom imperfect interface parameters properties for <a href="multimaterialmode.html">multimaterial mode MPM</a> contact between the current material and another material.</li>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<math>{Q\over\rho} = \Delta x|D_{kk}\bigl|(A_1C + A_2\Delta x|D_{kk}|\bigr)</math>


</ul>
where &Delta;<i>x</i> is the cell size of the mesh, |<i>D<sub>kk</sub></i>| is the relative volume change rate (''i.e.'' trace of the velocity gradient), <i>C</i> is the bulk wave speed in the material, and ''A''<sub>1</sub> and ''A''<sub>2</sub> are adjustable constants.


== XML Notes ==
{| class="wikitable"
|-
! Property !! Description !! Units !! Default
|-
| ArtificialVisc || Set to "on" or "off" to activate artificial viscosity. In <tt>XML</tt> files, an <tt><Artificial/></tt> command turns it on and its absence keeps the default setting of "off". || none || off
|-
| avA1 || The ''A''<sub>1</sub> constant in the artificial viscosity law || none || 0.2
|-
| avA2 || The ''A''<sub>2</sub> constant in the artificial viscosity law || none || 2.0
|}


<p>There are some material properties that are common to all material types. These properties are only used in NairnMPM and are not allowed in NairnFEA. Some of these properties only work for selected materials (refer to each material for possible details). The common properties are
The artificial viscosity property is supported in some isotropic materials (because the theory assumes isotropy). If you use these commands in a material that does not support it, an error will result. The following materials currently support artificial viscosity:
</p>


<dl>
* [[Isotropic, Hyperelastic-Plastic Material|HEIsotropic]] and its [[Material Models#Material Class Hierarchy|subclasses]]
** [[Isotropic, Hyperelastic-Plastic Mie-Grüneisen Material|HEMGEOSMaterial]]
* [[Mooney Material|Mooney]]
* [[Neo-Hookean Material|Neohookean]]
* [[Tait Liquid Material|Tait Liquid]]
* [[Ideal Gas Material|Ideal Gas]]
* [[JWLPlusPlus Material|JWLPlusPlus]]
* [[Viscoelastic Material|Viscoelastic]]
* [[Isotropic, Elastic-Plastic Material|IsoPlasticity]]


<dt>&lt;rho&gt;1.5&lt;/rho&gt;</dt>
== Poroelasticity Properties ==
<dd>The material density (in g/cm<sup>3</sup>) (can be used in NairnFEA files, but it is ignored).</dd>
   
<dt>&lt;Cp&gt;1.5&lt;/Cp&gt; or &lt;Cv&gt;1.5&lt;/Cv&gt;</dt>
<dd>The material's heat capacity (in J/(kg-K)). It is used when doing <a href="thermal.html">conduction calculations and/or coupled mechanical energy</a> and by some material constitutitive laws. Both <code>Cp</code> and <code>Cv</code> set the same (and only) heat capacity used in the MPM calculations. MPM is a strain-drive method, which implies the heat capacity should be interpreted as the constant-volume heat capacity.</dd>
   
<dt><a name="csat"></a>&lt;csat&gt;0.5&lt;/csat&gt;</dt>
<dd>The material's saturation concentration entered as a weight fraction. Only used when doing diffusion calculations.</dd>


<dt>&lt;Propagate criterion='1' direction='0' traction='1'/&gt;<br>
Some materials support [[Poroelasticity Calculations|poroelasticity calculations]] and the properties in this section control pore pressure flow between particles and coupling between stress and strain and pore pressure. The properties to use depend on symmetry of the parent material.
&lt;AltPropagate criterion='7' direction='4' traction='1'/&gt;</dt>
<dd>These properties let each material have a different primary and secondary propagation criterion and/or directions than the ones defined in the <a href="mpmheader.html">MPM Header</a>. Set these criterion as described for the <a href="mpmheader.html#propagate"><code>Propagate</code></a> and <a href="mpmheader.html#altpropagate">AltPropagate</a> commands in the <a href="mpmheader.html">MPM Header</a>. To activate propagation:
<ol>
<li>The <a href="mpmheader.html">MPM Header</a> must have a <a href="mpmheader.html#propagate"><code>Propagate</code></a> command.</li>
<li>The <a href="mpmheader.html">MPM Header</a> can optionally have an <a href="mpmheader.html#altpropagate"><code>AltPropagate</code></a> command.</li>
<li>Any material can override the default primary and/or secondary criteria by using a <code>Propagate</code> and/or <code>AltPropagate</code> command within its material defintion. Even if the <a href="mpmheader.html">MPM Header</a> does not have an <a href="mpmheader.html#altpropagate"><code>AltPropagate</code></a> command, a material can add a secondary criterion with its own <code>AltPropagate</code> command.</li>
</ol>
</dd>
   
<dt>&lt;maxLength&gt;20&lt;/maxLength&gt;</dt>
<dd>Define the maximum length of the crack (in mm). The crack growth will stop when this length is reached. This limit, however, is only active when there is crack propagation with <a href="mpmheader.html#cracks">criterion</a> 2.</dd>


<dt>&lt;xGrow&gt;20&lt;/xGrow&gt; and &lt;yGrow&gt;0&lt;/yGrow&gt;</dt>
=== Isotropic Poroelasticity Properties ===
<dd>These properties (one or both with default 0) specify a unit vector for a constant crack growth direction. They only apply for crack propagation with <a href="mpmheader.html#cracks">criterion</a> 2 and then only when that criterion selects the default propagation direction. The result is constant crack growth direction regardless of stress state or crack tip orientation. Any input vector will be normalized to a unit vector. If constant crack growth directions with a <a href="cracks.html">fixed crack</a> are located precisely on grid lines, it is possible for the crack algorithm to not recognize the crack plane. It is better to move such a crack sligthly off grid lines.
</dd>


<dt><a name="artificialvisc"></a>&lt;ArtificialVisc/&gt;, &lt;avA1&gt;.2&lt;/avA1&gt;, and &lt;avA2&gt;2.0&lt;/avA2&gt;</dt>
{| class="wikitable"
<dd>This first empty tag turns artificial visocsity on, which adds a pressure, Q, related to velocity gradient on the particle, but only when it is compressing. The equation is
|-
<blockquote>
! Property !! Description !! Units !! Default
Q = &Delta;x*|D<sub>kk</sub>|*(A1*C + A2*&Delta;x*|D<sub>kk</sub>|)
|-
</blockquote>
| Ku || Undrained bulk modulus. It must be greater than the material bulk modulus. || [[ConsistentUnits Command#Legacy and Consistent Units|pressure units]] || none
where &Delta;x is the cell size of the mesh, |D<sub>kk</sub>| is the relative volume change rate, C is the bulk wave speed in the material, and A1 and A2 are adjustable constants, which are entered using the above &lt;avA1&gt; and &lt;avA2&gt; commands. The defaults are off, A1=0.2, and A2=2.0. It is common, but not required, to choose A2 = 10*A1.<br>
|-
&nbsp;&nbsp;&nbsp;&nbsp;The artificial viscosity property is supported in some materials (as documented under each material). If you use these commands in a material that does not support, an error will result.
| alphaPE || The poroelasticity Biot coefficient that translates fraction of volume change the results in pore pressure change. It must be between 0 and 1 || none || 0
</dd>
|-
| Darcy || Darcy law permittivity for the material || [[ConsistentUnits Command#Legacy and Consistent Units|length units<sup>2</sup>]] || 0
|}


<dt><a name="friction"></a>&lt;Friction mat='2'&gt;0.3&lt;/Friction&gt;
=== Transversely Isotropic Poroelasticity Properties ===  
</dt>
<dd>This command within a material definition allows multimaterial mode to have different contact laws between different pairs of materials. Set the contact law as explained for <a href="mpmheader.html#friction">Friction command</a> in the <a href="mpmheader.html#cracks">Cracks section of the MPM header</a> except that a setting less than -10 means to assume materials are always in contact and revert to single material mode. The <code>mat</code> attribute defines contact mechanics between this material and the one listed (self contact is ignored). Alternatively, the other material can be specified by name using a <code>matname</code> attribute
as explained <a href="definematl.html#usemats">here</a>. If a node has more than two materials, the contact law used will be one between each material and the other material with the most volume.<br>
&nbsp;&nbsp;&nbsp;&nbsp;You only need to set friction on one of the materials in a pair. If you set both, the setting for the second material will override the one from the first material.
</dd>


<dt>&lt;color red='0.1' green='0.5' blue='0.3' alpha='1.0/&gt;
{| class="wikitable"
</dt>
|-
<dd>Sets the color of the material with red, green, blue, and alpha values between 0.0 and 1.0. The color may be used by post-analysis software, but it has no affect on the calculations.
! Property !! Description !! Units !! Default
</dd>
|-
| Ku || Undrained bulk modulus. It must be greater than the material bulk modulus. ||  [[ConsistentUnits Command#Legacy and Consistent Units|pressure units]] || none
|-
| alphaAPE<br>alphaTPE || The axial and transverse poroelasticity Biot coefficients that translate how strains results in pore pressure change. They must be between 0 and 1 || none || 0
|-
| DarcyA<br>DarcyT || The axial and transverse Darcy law permittivities for the material || [[ConsistentUnits Command#Legacy and Consistent Units|length units<sup>2</sup>]] || 0
|}
 
=== Orthotropic Poroelasticity Properties ===


</dl>
{| class="wikitable"
|-
! Property !! Description !! Units !! Default
|-
| Ku || Undrained bulk modulus. It must be greater than the material bulk modulus. ||  [[ConsistentUnits Command#Legacy and Consistent Units|pressure units]] || none
|-
| alphaxPE<br>alphayPE<br>alphazPE || The poroelasticity Biot coefficients for the three orthotropic directions that translate how strains results in pore pressure change. They must be between 0 and 1 || none || 0
|-
| alphaRPE<br>alphaZPE<br>alphaTPE || For cylindrical, orthotropic materials, these are radial, axial, and hoop Biot coefficients. They are synonyms for orthotropic coefficients (with R=x, Z=y, and T=z). || none || 0
|-
| Darcyx<br>Darcyy<br>Darcyz || The Darcy law permittivities for the three orthotropic directions material || [[ConsistentUnits Command#Legacy and Consistent Units|length units<sup>2</sup>]] || 0
|-
| DarcyR<br>DarcyZ<br>DarcyT || For cylindrical, orthotropic materials, these are radial, axial, and hoop Darcy law permittivities for the material. They are synonyms for orthotropic coefficients (with R=x, Z=y, and T=z). || [[ConsistentUnits Command#Legacy and Consistent Units|length units<sup>2</sup>]] || 0
|}

Latest revision as of 13:56, 4 January 2021

These material properties are common to all types of materials used in MPM simulations.

Basic Properties

These are basic material properties.

Property Description Units Default
rho The material's initial density. density units 0.001
Cv The constant-volume heat capacity. It is used when doing conduction calculations and/or coupled mechanical energy and by some material constitutitive laws. You do not need to enter a constant pressure heat capacity (Cp). Instead, when MPM thermodynamics are done correctly, the simulations will automatically account for all boundary condition effects (constant volume, constant pressure, or mixed) on effective heat capacity. heat capacity units 1
kCond Thermal conductivity for isotropic materials (anisotropic materials will have alternate properties for setting thermal conductivity tensor). conductivity units 0
csat The saturation concentration potential as a weight friction from 0 to 1. it is only used when doing diffusion calculations. none 1
beta Moisture expansion coefficient for isotropic materials (anisotropic materials will have alternate properties for setting moisture expansion tensor). strain/(wt fraction) 0
D Solvent diffusion constant for isotropic materials (anisotropic materials will have alternate properties for setting diffusion tensor). diffusion units 0
largeRotation To use (1) or not use (0) polar decomposition when calculating rotations in small-stain materials. This option applies only to linear elastic small strain materials and to elastic-plastic small strain materials. Option 1 finds incremental deformation gradient use the selected number of terms and evaluates the rotational part of that increment using polar decomposition. The stress update is then rotated according to the decomposed rotation. In contrast, option 0 finds incremental deformation gradient from a linear expansion (kmax=1) and evaluates the rotation component using second order (in 2D) or first order (in 3D) approximate polar decomposition. The stress update is rotated by standard hypoelastic methods. Both methods can handle large total rotations (provided they are incrementally small). Option 1 may be more accurate, but it is less efficient. none 0
matDamping Sets custom particle damping to apply only to particles of this material. Its value replaces global particle damping setting and must be a constant (function of time not allowed). In XML files, set matDamping using <PDamping>matDamping</PDamping> 1/time units none
color Sets the color of the material. The color is used in material point method plots material type in NairnFEAMPM and in NairnFEAMPMViz. If no color is provided, a color will be picked from the current spectrum using the material number. In scripted files, this property takes four arguments being red, green, blue, and alpha values between 0.0 and 1.0. A single argument means to set gray level between 0.0 and 1.0 (with alpha=1.0). Three arguments means set red, green, and blue with alpha=1.0. In XML files, the color is set with "red", "green", "blue" and "alpha" attributes (and the element's content is ignored). none none

Fracture Toughness Properties

These properties set material properties that determine the fracture toughness of the material and control various aspects of crack propagation.

Property Description Units Default
JIc Critical energy release rate fracture toughness for mode I. It is only used for crack propagation by criteria 2, 3, or 7. For criterion 2, it is only used if initTime is not specified. It is also used to set toughness of traction law materials. energy release units none
JIIc Critical energy release rate fracture toughness for mode II. It is currently only used to set toughness of traction law materials. energy release units none
KIc Critical mode I stress intensity factor. It is only used for crack propagation by criteria 1, 4, or 5. stress intensity units none
KIIc Critical mode II stress intensity factor. It is only used for crack propagation by criteria 1, 4, or 5. stress intensity units none
KIexp Exponent p in the elliptical criteria for crack growth. It is only used for crack propagation by criterion 5. none 2
KIIexp Exponent q in the elliptical criteria for crack growth. It is only used for crack propagation by criterion 5. none 2
delIc Critical crack opening displacement for mode I. Only used for crack propagation by criterion 6. It is also used by traction-law materials. length units none
delIIc Critical crack opening displacement for mode II. Only used for crack propagation by criterion 6. It is also used by traction-law materials. length units none
initTime The time when crack propagation starts. It is only used for crack propagation by criterion 2. For criterion 2, when initTime is specified, takes precedence over the JIc property. alt time units none
speed The crack speed in steady state crack propagation. This speed, however, is only active for crack propagation by criterion 2. (also used in criterion 3 as an initial crack speed, but that criterion is not meant for general use) alt velocity units 1
maxLength The maximum crack length for steady state crack propagation. The simulation will stop soon after crack reaches the input length. This length, however, is only active for crack propagation by criterion 2. length units none
nmix An exponent used in mixed-modes failure of some traction laws. none 1
allowsCracks Set to 1 to allow the material to have cracks or 0 to not allow cracks. The default is 1 for non-rigid materials and 0 for rigid materials. Setting a non-rigid material to 0 lets that material enter a crack in another material such as to simulate cutting, to wedge open a crack, or model hydraulic fracturing. The non-cracking material can open the crack in the other material by contact on its crack surfaces. none varies

Crack Propagation Properties

The setting of crack propagation properties are done differently for scripted and XML files. For scripted commands, you can set the following material properties:

Property Description Units Default
criterion To set a custom crack propagation criterion. You can use any propagation option. Setting this property in a material overrides the default propagation criterion setting. This command is for scripted files only; see below to set criterion in XML files. none none
direction To set a custom crack propagation direction. You can use any direction option. Setting this property in a material overrides the default propagation direction setting. This command is for scripted files only; see below to set direction in XML files. none none
traction To set a custom traction law to create for crack propagation in this material. A traction law set in a material overrides the default traction law. The traction law can be set by material ID (if the traction law has already been defined) or by number (if it is not defined yet). This command is for scripted files only; see below to set traction in XML files. none none
altcriterion Same as "criterion" property above except that it applies to the alternate propagation criterion for the material none none
altdirection Same as "direction" property above except that it applies to the alternate propagation criterion for the material none none
alttraction Same as "traction" property above except that it applies to the alternate propagation criterion for the material none none
xGrow This property along with yGrow (if on one given the other is set to 0) specify a unit vector for a constant crack growth direction. t is only used for crack propagation by criterion 2 and then only if that criterion is using its default propagation direction. The result is a constant crack growth direction regardless of stress state or crack tip orientation. Any input vector will be normalized to a unit vector. If a constant crack growth direction with a fixed crack is located precisely on grid lines, it is possible the crack algorithm will not recognize the crack plane. Is it better to move such a crack slightly off grid lines. none none
yGrow Crack growth direction - see xGrow above. none none
constantTip Set to 0 or 1. The default of 0 means the crack tip will track the material around the crack tip. Changing it to 1 means crack tips with this material will always use this material even if the crack propagates into another material. The default 0 allows modeling crack growth in composites with fracture properties changing as cracks move between materials. Using 1 allows modeling multiple cracks in the same material having different fracture propertie by using the following steps:
  1. Define multiple materials that are identical except for their fracture and/or propagation properties and all with constantTip=1.
  2. Create a model and uses one of the materials for all material points.
  3. For each crack tip, assign its material to be the material in step 1 with the appropriate fracture properties.
none 0

Propagation Properties in XML Files

In XML files, the criterion, direction, and traction properties (and the analogous alternate propagation properties) are set differently. To set crack propagation criteria, you use instead 

<Propagate criterion='(critNum)' direction='(dirNum)' traction='(traction)'/>
<AltPropagate criterion='(critNum)' direction='(dirNum)' traction='(traction)'/>

where the settings are the same as defined in the default crack propagation commands (or the alternate propagation command), but the XML element is now used within a <Material> definition instead of within the <Cracks> element in the <MPMHeader>.

XML files set xGrow, yGrow, and constantTip as ordinary properties and they function as described above.

Contact Properties

These properties can set custom the contact mechanics between two specific materials. If these properties are not used, material-to-material contact will be determined by the global setting made using a ContactMM command (for scripted files) or a <Friction> command within the <MultimaterialMode> element in the <MPMHeader> (for XML files). For more details see these commands under friction settings and interface setting.

Property Description Units Default
Contact
<Friction>		 In scripted files, use the Contact property with two arguments — (lawID) and (matID). It sets a custom contact law (by its material ID in (lawID)) to use for contact between this material and the material defined by (matID). You only need to enter a Contact property for one of the two materials in each pair. Since you have to provide the other material's ID, you should use the Contact command in the second material after the first one has already been defined (and after the contact law has been defined as well).

In XML files use the command:
   <Friction law='(lawnum)' lawname='(lawID)' mat='(matnum)' matname='(matID)'/>
to set the custom contact law by number or ID for both the law and the other material. If both number and ID are used, the ID will take precedence. || none || none

shareMatField The value selects another material to share the same velocity field such that the two materials move together and interact by perfect contact (see below for how to have more than two materials share one field). Although materials in the same velocity field interact by perfect contact, they can interact with materials in different velocity fields by multimaterial mode contact or interface laws. Finally, materials sharing velocity fields must use "compatible" fields (e.g., all rigid or all nonrigid); an error will occur (with an explanation) if you try to share velocity fields of incompatible materials. none none

In scripted files, the other material is specified by its material ID, which means the Friction and Interface commands must be used in the secondly-defined materials (such that material ID for the first material is available). In XML files, the second material in Friction and Interface commands is defined by number (or by name) using a mat or matname attribute. You only need a Friction or an Interface command in one material for each pair of materials with custom contact properties.

To share velocity fields, you first create one "base material" that is not shared and then any number that share that field with a shareMatField property. The base material is specified in the shareMatField property by material ID in scripted files, but must be by number (as the value of the property command) in XML files. To specify custom friction or interfaces between shared materials and other materials, you use any material in the shared block; if you use custom commands for more then one material in a block, only the last one will be used.

Deprecated Contact Properties

The following material properties for setting contact properties are deprecated. They should be replaced by the Contact property above. Note that in XML files, contact properties are still with with a <Friction> command but prior attributes are deprecated to be replaced by attributes to specify the desired contact law.

Property Description Units Default
Friction A Friction property within a material definition can define custom frictional properties for multimaterial mode MPM contact between the current material and another material. This property takes two parameters; the first is the same as for the standard Friction command and the second gives the other material. none none
Interface An Interface property within a material definition can define custom imperfect interface parameters properties for multimaterial mode MPM contact between the current material and another material. This property takes four parameters; the first three are the same as for a standard ImperfectInterface command (which is actually a <Friction> element in XML files) and the fourth gives the other material. none none

Artificial Viscosity

Some materials support artificial viscosity to dampen pressure waves. When it is on, it adds a pressure, Q, related to velocity gradient on the particle, but only when it is compressing. The equation is

      [math]\displaystyle{ {Q\over\rho} = \Delta x|D_{kk}\bigl|(A_1C + A_2\Delta x|D_{kk}|\bigr) }[/math]

where Δx is the cell size of the mesh, |Dkk| is the relative volume change rate (i.e. trace of the velocity gradient), C is the bulk wave speed in the material, and A1 and A2 are adjustable constants.

Property Description Units Default
ArtificialVisc Set to "on" or "off" to activate artificial viscosity. In XML files, an <Artificial/> command turns it on and its absence keeps the default setting of "off". none off
avA1 The A1 constant in the artificial viscosity law none 0.2
avA2 The A2 constant in the artificial viscosity law none 2.0

The artificial viscosity property is supported in some isotropic materials (because the theory assumes isotropy). If you use these commands in a material that does not support it, an error will result. The following materials currently support artificial viscosity:

Poroelasticity Properties

Some materials support poroelasticity calculations and the properties in this section control pore pressure flow between particles and coupling between stress and strain and pore pressure. The properties to use depend on symmetry of the parent material.

Isotropic Poroelasticity Properties

Property Description Units Default
Ku Undrained bulk modulus. It must be greater than the material bulk modulus. pressure units none
alphaPE The poroelasticity Biot coefficient that translates fraction of volume change the results in pore pressure change. It must be between 0 and 1 none 0
Darcy Darcy law permittivity for the material length units2 0

Transversely Isotropic Poroelasticity Properties

Property Description Units Default
Ku Undrained bulk modulus. It must be greater than the material bulk modulus. pressure units none
alphaAPE
alphaTPE		 The axial and transverse poroelasticity Biot coefficients that translate how strains results in pore pressure change. They must be between 0 and 1 none 0
DarcyA
DarcyT		 The axial and transverse Darcy law permittivities for the material length units2 0

Orthotropic Poroelasticity Properties

Property Description Units Default
Ku Undrained bulk modulus. It must be greater than the material bulk modulus. pressure units none
alphaxPE
alphayPE
alphazPE		 The poroelasticity Biot coefficients for the three orthotropic directions that translate how strains results in pore pressure change. They must be between 0 and 1 none 0
alphaRPE
alphaZPE
alphaTPE		 For cylindrical, orthotropic materials, these are radial, axial, and hoop Biot coefficients. They are synonyms for orthotropic coefficients (with R=x, Z=y, and T=z). none 0
Darcyx
Darcyy
Darcyz		 The Darcy law permittivities for the three orthotropic directions material length units2 0
DarcyR
DarcyZ
DarcyT		 For cylindrical, orthotropic materials, these are radial, axial, and hoop Darcy law permittivities for the material. They are synonyms for orthotropic coefficients (with R=x, Z=y, and T=z). length units2 0
Retrieved from "http://osupdocs.forestry.oregonstate.edu/index.php?title=Common_Material_Properties&oldid=7688"