Difference between revisions of "Anisotropic, Elastic-Plastic Material"
| (38 intermediate revisions by the same user not shown) | |||
| Line 1: | Line 1: | ||
== Constitutive Law == | == Constitutive Law == | ||
This [[Material Models|MPM material]] is identical to an [[Orthotropic Material|orthotropic material]] in the elastic regime, but can plastically deform according to a built-in, anistropic Hill yielding criterion. The Hill plastic yield | This [[Material Models|MPM material]] is identical to an [[Orthotropic Material|orthotropic material]] in the elastic regime, but can plastically deform according to a built-in, anistropic Hill yielding criterion. The two Hill plastic yield criteria available are <tt>HillStyle=1</tt>: | ||
| | ||
| Line 7: | Line 7: | ||
+ 2L\tau_{yz}^2 + 2M\tau_{xz}^2 + 2N\tau_{xy}^2} = H(\alpha)</math> | + 2L\tau_{yz}^2 + 2M\tau_{xz}^2 + 2N\tau_{xy}^2} = H(\alpha)</math> | ||
where σ and τ are normal and shear stresses in the material axis system after rotation from the analysis coordinates, <math>H(\alpha)</math> is a hardening law, and <math>\alpha</math> is a plastic hardening variable. The remaining constants are determined by the yield stresses: | and <tt>HillStyle=2</tt>: | ||
| |||
<math>F(\sigma_{yy}-\sigma_{zz})^2 + G(\sigma_{xx}-\sigma_{zz})^2 + H(\sigma_{yy}-\sigma_{xx})^2 | |||
+ 2L\tau_{yz}^2 + 2M\tau_{xz}^2 + 2N\tau_{xy}^2 = \bigl(H(\alpha)\bigr)^2</math> | |||
where σ and τ are normal and shear stresses in the material axis system after rotation from the analysis coordinates, <math>H(\alpha)</math> is a hardening law, and <math>\alpha</math> is a plastic hardening variable. The reason for two styles is discussed [[#Why Two Hill Criteria|below]]. The remaining constants are determined by the yield stresses: | |||
| | ||
| Line 36: | Line 42: | ||
| | ||
<math>{\sigma_{Y,ii}\over 1+R} \le \sigma_{Y,kk} \le {\sigma_{Y,ii}\over 1-R}</math> | <math>{\sigma_{Y,ii}\over 1+R} \le \sigma_{Y,kk} \le {\sigma_{Y,ii}\over 1-R}</math> | ||
=== Example Yield Strengths === | |||
Some examples in various types of materials follow. If one direction is prevented from yielding by setting its yield strength to ∞, then | Some examples in various types of materials follow. If one direction is prevented from yielding by setting its yield strength to ∞, then | ||
| Line 42: | Line 50: | ||
<math>\sigma_{Y,jj} = \infty \quad\implies\quad R = 0 \quad\implies\quad \sigma_{Y,kk} = \sigma_{Y,ii} </math> | <math>\sigma_{Y,jj} = \infty \quad\implies\quad R = 0 \quad\implies\quad \sigma_{Y,kk} = \sigma_{Y,ii} </math> | ||
In other words, the other two direction must have the same yield stress. If two directions have the same yield strength, then the third direction must greater than half that value: | In other words, the other two direction must have the same yield stress. If two directions have the same yield strength, then the third direction must be greater than half that value: | ||
| | ||
| Line 55: | Line 63: | ||
Because the thickness direction may be much weaker than the in-plane direction, this limitation limits modeling thickness-direction plasticity. Either Hill plasticity needs to be replaced with another anisotropic yielding model or it is realistic that thickness direction of such materials has very little plasticity. They may still have plasticity in shear. | Because the thickness direction may be much weaker than the in-plane direction, this limitation limits modeling thickness-direction plasticity. Either Hill plasticity needs to be replaced with another anisotropic yielding model or it is realistic that thickness direction of such materials has very little plasticity. They may still have plasticity in shear. | ||
=== Hardening Laws === | === Hardening Laws === | ||
There hardening laws | There currently available hardening laws: | ||
| | ||
| Line 68: | Line 75: | ||
| | ||
<math>{\rm Exponential:} \qquad H(\alpha) = 1 + \frac{K}{\beta}(1-\exp(-\beta\alpha))</math> | <math>{\rm Exponential:} \qquad H(\alpha) = \left\{ \begin{array}{ll} | ||
1 + \frac{K}{\beta}(1-\exp(-\beta\alpha)) & \alpha \le \alpha_{max} \\ | |||
1 + K^*\alpha & \alpha \gt \alpha_{max}\end{array} \right.</math> | |||
Note, that because <math>H'(\alpha)</math> in the first law is infinite for <math>\alpha=0</math> when ''n'' is less than one, this law becomes unstable for ''n'' less than about 0.7. | |||
The parameter <math>\alpha_{max}</math> is to prevent exponential hardening from devolving into elastic-plastic deformation. When this parameter is used, the hardening is converted to linear hardening once <math>\alpha</math> exceeds <math>\alpha_{max}</math> and the slope of the linear hardening, <math>K^*</math>, is set equal to the slope of the exponential hardening law at <math>\alpha_{max}</math>: | |||
| |||
<math>K^* = K \exp(-\beta\alpha_{max})</math> | |||
Thus, both the hardening law and its slope are continuous at <math>\alpha_{max}</math>. One way to pick <math>\alpha_{max}</math> is to choose a minimum relative slope compared to initial slope and then set: | |||
| |||
<math>\alpha_{max} = -\frac{\ln s_{rel}}{\beta} \qquad{\rm where}\qquad s_{rel} = \frac{K^*}{K}</math> | |||
== Why Two Hill Criteria == | |||
The two Hill criteria described [[#Constitutive Law|above]] can be written with a symmetric, fourth-rank tensor as | |||
| |||
<math>\sqrt{\mathbf{A}\cdot \sigma \cdot \mathbf{A}} = H(\alpha) | |||
\quad {\rm and} \quad \mathbf{A}\cdot \sigma \cdot \mathbf{A} = \bigl(H(\alpha)\bigr)^2</math> | |||
where elements of <math>\mathbf{A}</math> are derived from material properties ''F'' through ''N''. A discussion with implementation details for <tt>HillStyle=1</tt> (the first one with a square root) can be found in Simo and Hughes.<ref name="SH"/> Unfortunately, that derivation has a step that multiplies by <math>\mathbf{A}^{-1}</math> but the tensor for Hill plasticity is singular. Despite this inconsistency, the overall implementation seems to work well. | |||
To avoid this problem, the criteria can use <tt>HillStyle=2</tt> instead (with squared terms). Implementation by this approach<ref name="dB/> does not need to use <math>\mathbf{A}^{-1}</math>. It is now the default style for this material type. | |||
== Material Properties == | == Material Properties == | ||
| Line 94: | Line 125: | ||
| nhard || Hardening law n parameter. If n>0, the Nonlinear 1 law is used; if n<0, the Nonlinear 2 law is used with abs(n). || dimensionless || 1 | | nhard || Hardening law n parameter. If n>0, the Nonlinear 1 law is used; if n<0, the Nonlinear 2 law is used with abs(n). || dimensionless || 1 | ||
|- | |- | ||
| exphard || Hardening law β parameter in the exponential hardening law || dimensionless || 0 | | exphard || Hardening law β parameter in the exponential hardening law. Use of this parameter also switches to using exponential hardening || dimensionless || 0 | ||
|- | |||
| alphaMax || Optional <math>\alpha_{max}</math> where exponential hardening converts to linear hardening || dimensionless || ∞ | |||
|- | |||
| HillStyle || Use 1 for square-root criterion<ref name="SH"/> or 2 for squared criterion<ref name="dB"/> || dimensionless || 2 | |||
|- | |- | ||
| [[Common Material Properties#Basic Properties|largeRotation]] || If used, this setting is ignored and material always uses 1 || dimensionless || 1 (fixed) | | [[Common Material Properties#Basic Properties|largeRotation]] || If used, this setting is ignored and material always uses 1 || dimensionless || 1 (fixed) | ||
| Line 110: | Line 145: | ||
== Examples == | == Examples == | ||
== References == | |||
<references> | |||
<ref name="SH">J. C. Simo and T. J. R. Hughes, "Computational Inelasticity", Springer-Verlag, New York, 1997 (page 96)</ref> | |||
<ref name="dB">R. de Borst and P. H. Feenstra. "Studies in anisotropic plasticity with reference to the Hill criterion," ''Int. J. Numer. Meth. Engng'', '''29''', 315–336 (1990).</ref> | |||
</references> | |||
Latest revision as of 13:14, 25 August 2023
Constitutive Law
This MPM material is identical to an orthotropic material in the elastic regime, but can plastically deform according to a built-in, anistropic Hill yielding criterion. The two Hill plastic yield criteria available are HillStyle=1:
[math]\displaystyle{ \sqrt{F(\sigma_{yy}-\sigma_{zz})^2 + G(\sigma_{xx}-\sigma_{zz})^2 + H(\sigma_{yy}-\sigma_{xx})^2 + 2L\tau_{yz}^2 + 2M\tau_{xz}^2 + 2N\tau_{xy}^2} = H(\alpha) }[/math]
and HillStyle=2:
[math]\displaystyle{ F(\sigma_{yy}-\sigma_{zz})^2 + G(\sigma_{xx}-\sigma_{zz})^2 + H(\sigma_{yy}-\sigma_{xx})^2 + 2L\tau_{yz}^2 + 2M\tau_{xz}^2 + 2N\tau_{xy}^2 = \bigl(H(\alpha)\bigr)^2 }[/math]
where σ and τ are normal and shear stresses in the material axis system after rotation from the analysis coordinates, [math]\displaystyle{ H(\alpha) }[/math] is a hardening law, and [math]\displaystyle{ \alpha }[/math] is a plastic hardening variable. The reason for two styles is discussed below. The remaining constants are determined by the yield stresses:
[math]\displaystyle{ F = {1\over 2}\left({1\over \sigma_{Y,yy}^2} + {1\over \sigma_{Y,zz}^2} - {1\over \sigma_{Y,xx}^2}\right) \qquad\qquad L = {1\over 2\tau_{Y,yz}^2} }[/math]
[math]\displaystyle{ G= {1\over 2}\left({1\over \sigma_{Y,xx}^2} + {1\over \sigma_{Y,zz}^2} - {1\over \sigma_{Y,yy}^2}\right) \qquad\qquad M = {1\over 2\tau_{Y,xz}^2} }[/math]
[math]\displaystyle{ H = {1\over 2}\left({1\over \sigma_{Y,xx}^2} + {1\over \sigma_{Y,yy}^2} - {1\over \sigma_{Y,zz}^2}\right) \qquad\qquad N = {1\over 2\tau_{Y,xy}^2} }[/math]
where σY and τY are yield stresses for loading in the indicated direction. The yield stresses have to be selected such that the plastic potential is positive semidefinite. Analysis shows that all tensile yield stresses must satisfy:
[math]\displaystyle{ \left|{1\over \sigma_{Y,ii}^2} - {1\over \sigma_{Y,jj}^2} \right| \le {1\over \sigma_{Y,kk}^2} \le {1\over \sigma_{Y,ii}^2} + {1\over \sigma_{Y,jj}^2} }[/math]
where i, j, and k are any combination or x, y, and z. In more practical terms, if two yield stresses are related by some ratio:
[math]\displaystyle{ R = {\sigma_{Y,ii}^2\over \sigma_{Y,jj}^2} \quad\rm{or} \quad {1\over \sigma_{Y,jj}^2} = {R\over \sigma_{Y,ii}^2} }[/math]
and [math]\displaystyle{ R\le1 }[/math] (i.e., [math]\displaystyle{ \sigma_{Y,ii}\le\sigma_{Y,jj} }[/math]) then the third yield stress is bracketed by:
[math]\displaystyle{ {\sigma_{Y,ii}\over 1+R} \le \sigma_{Y,kk} \le {\sigma_{Y,ii}\over 1-R} }[/math]
Example Yield Strengths
Some examples in various types of materials follow. If one direction is prevented from yielding by setting its yield strength to ∞, then
[math]\displaystyle{ \sigma_{Y,jj} = \infty \quad\implies\quad R = 0 \quad\implies\quad \sigma_{Y,kk} = \sigma_{Y,ii} }[/math]
In other words, the other two direction must have the same yield stress. If two directions have the same yield strength, then the third direction must be greater than half that value:
[math]\displaystyle{ \sigma_{Y,jj} = \sigma_{Y,ii}\quad\implies\quad R = 1 \quad\implies\quad \sigma_{Y,kk} \ge {\sigma_{Y,ii}\over 2} }[/math]
If two directions are similar, such as in in-plane, random fiber composites (e.g., medium density fiber board, particle board, fiberglass composites, or paperboard), the yield strength in the thickness direction must still be:
[math]\displaystyle{ \sigma_{Y,jj} \sim \sigma_{Y,ii} \quad\implies\quad R \sim 1 \quad\implies\quad \sigma_{Y,zz} \gtrapprox \frac{1}{2}\min(\sigma_{Y,xx},\sigma_{Y,yy}) }[/math]
Because the thickness direction may be much weaker than the in-plane direction, this limitation limits modeling thickness-direction plasticity. Either Hill plasticity needs to be replaced with another anisotropic yielding model or it is realistic that thickness direction of such materials has very little plasticity. They may still have plasticity in shear.
Hardening Laws
There currently available hardening laws:
[math]\displaystyle{ {\rm Nonlinear\ 1:} \qquad H(\alpha) = 1 + K \alpha^n }[/math]
[math]\displaystyle{ {\rm Nonlinear\ 2:} \qquad H(\alpha) = (1 + K \alpha)^n }[/math]
[math]\displaystyle{ {\rm Exponential:} \qquad H(\alpha) = \left\{ \begin{array}{ll} 1 + \frac{K}{\beta}(1-\exp(-\beta\alpha)) & \alpha \le \alpha_{max} \\ 1 + K^*\alpha & \alpha \gt \alpha_{max}\end{array} \right. }[/math]
Note, that because [math]\displaystyle{ H'(\alpha) }[/math] in the first law is infinite for [math]\displaystyle{ \alpha=0 }[/math] when n is less than one, this law becomes unstable for n less than about 0.7.
The parameter [math]\displaystyle{ \alpha_{max} }[/math] is to prevent exponential hardening from devolving into elastic-plastic deformation. When this parameter is used, the hardening is converted to linear hardening once [math]\displaystyle{ \alpha }[/math] exceeds [math]\displaystyle{ \alpha_{max} }[/math] and the slope of the linear hardening, [math]\displaystyle{ K^* }[/math], is set equal to the slope of the exponential hardening law at [math]\displaystyle{ \alpha_{max} }[/math]:
[math]\displaystyle{ K^* = K \exp(-\beta\alpha_{max}) }[/math]
Thus, both the hardening law and its slope are continuous at [math]\displaystyle{ \alpha_{max} }[/math]. One way to pick [math]\displaystyle{ \alpha_{max} }[/math] is to choose a minimum relative slope compared to initial slope and then set:
[math]\displaystyle{ \alpha_{max} = -\frac{\ln s_{rel}}{\beta} \qquad{\rm where}\qquad s_{rel} = \frac{K^*}{K} }[/math]
Why Two Hill Criteria
The two Hill criteria described above can be written with a symmetric, fourth-rank tensor as
[math]\displaystyle{ \sqrt{\mathbf{A}\cdot \sigma \cdot \mathbf{A}} = H(\alpha) \quad {\rm and} \quad \mathbf{A}\cdot \sigma \cdot \mathbf{A} = \bigl(H(\alpha)\bigr)^2 }[/math]
where elements of [math]\displaystyle{ \mathbf{A} }[/math] are derived from material properties F through N. A discussion with implementation details for HillStyle=1 (the first one with a square root) can be found in Simo and Hughes.[1] Unfortunately, that derivation has a step that multiplies by [math]\displaystyle{ \mathbf{A}^{-1} }[/math] but the tensor for Hill plasticity is singular. Despite this inconsistency, the overall implementation seems to work well.
To avoid this problem, the criteria can use HillStyle=2 instead (with squared terms). Implementation by this approach[2] does not need to use [math]\displaystyle{ \mathbf{A}^{-1} }[/math]. It is now the default style for this material type.
Material Properties
| Property | Description | Units | Default |
|---|---|---|---|
| yldxx | Yield stress for axial loading in the x direction | pressure units | ∞ |
| yldyy | Yield stress for axial loading in the y direction | pressure units | ∞ |
| yldzz | Yield stress for axial loading in the z direction | pressure units | ∞ |
| yldxy | Yield stress for shear loading in the x-y plane | pressure units | ∞ |
| yldxz | Yield stress for shear loading in the x-z plane | pressure units | ∞ |
| yldyz | Yield stress for shear loading in the y-z plane | pressure units | ∞ |
| Khard | Hardening law K parameter | dimensionless | 0 |
| nhard | Hardening law n parameter. If n>0, the Nonlinear 1 law is used; if n<0, the Nonlinear 2 law is used with abs(n). | dimensionless | 1 |
| exphard | Hardening law β parameter in the exponential hardening law. Use of this parameter also switches to using exponential hardening | dimensionless | 0 |
| alphaMax | Optional [math]\displaystyle{ \alpha_{max} }[/math] where exponential hardening converts to linear hardening | dimensionless | ∞ |
| HillStyle | Use 1 for square-root criterion[1] or 2 for squared criterion[2] | dimensionless | 2 |
| largeRotation | If used, this setting is ignored and material always uses 1 | dimensionless | 1 (fixed) |
| (other) | All other properties are identical to the properties for an othotropic material. | varies | varies |
Although default tensile yield stresses are all infinite, one of them must be finite to use this material. Similarly, the combination of properties must satisfy conditions for positive definiteness described above. If not, an error message will appear and simulation will not run.
If only nhard is provided, the modeling will use Nonlinear 1 or Nonlinear 2 hardening law. If only exphard is provided, it will use the exponential hardening law. If both are given, the last one provided determines the hardening law.
History Variables
The one history variable is the value of plastic hardening variable [math]\displaystyle{ \alpha }[/math]. This variable can be archived as history variable 1.