Difference between revisions of "Transversely Isotropic Viscoelastic Material"
| Line 15: | Line 15: | ||
<math>\mathbf{C}(t) = \left[\begin{array}{cccccc} | <math>\mathbf{C}(t) = \left[\begin{array}{cccccc} | ||
K_T(t)+G_T(t) & K_T(t)-G_T(t) & \ell(t) & 0 & 0 & 0 \\ | K_T(t)+G_T(t) & K_T(t)-G_T(t) & \ell(t) & 0 & 0 & 0 \\ | ||
K_T(t)-G_T(t) & K_T+G_T & \ell & 0 & 0 & 0 \\ | K_T(t)-G_T(t) & K_T+G_T & \ell(t) & 0 & 0 & 0 \\ | ||
\ell(t) & \ell(t) & n(t) & 0 & 0 & 0 \\ | \ell(t) & \ell(t) & n(t) & 0 & 0 & 0 \\ | ||
0 & 0 & 0 & G_A(t) & 0 & 0 \\ | 0 & 0 & 0 & G_A(t) & 0 & 0 \\ | ||
| Line 22: | Line 22: | ||
\end{array}\right]</math> | \end{array}\right]</math> | ||
Here <math>K_T(t)</math> is the plane strain, bulk modulus, <math>G_T(t)</math> is the transverse shear modulus, <math>G_A(t)</math> is the axial shear modulus, and <math>n(t)</math> and <math>\ell(t)</math> give time-dependence of the ''C<sub>11</sub>'' and ''C<sub>12</sub>=C<sub>13</sub>'' elements of the stiffness matrix (as [[ | Here <math>K_T(t)</math> is the plane strain, bulk modulus, <math>G_T(t)</math> is the transverse shear modulus, <math>G_A(t)</math> is the axial shear modulus, and <math>n(t)</math> and <math>\ell(t)</math> give time-dependence of the ''C<sub>11</sub>'' and ''C<sub>12</sub>=C<sub>13</sub>'' elements of the stiffness matrix (as [[Transversely Isotropic Material#Transverse 1|defined here]]). | ||
Revision as of 16:38, 7 January 2021
Constitutive Law
(This material is available only in OSParticulas because it is still in development)
This anisotropic MPM material is a small strain, linear viscoelastic material that extends the Viscoelastic Material to model anisotropic viscoelasticity.
The stress (σ) and strain (ε) are related by:
[math]\displaystyle{ \sigma(t) = \mathbf{C}(t) * \varepsilon(t) }[/math]
Here [math]\displaystyle{ * }[/math] indicates convolution (or Boltzman's superposition) between time-dependent stiffness tensor ([math]\displaystyle{ \mathbf{C}(t) }[/math]) and strain tensor. In Voight-notation with unique axis in the z direction, the time-dependent stiffness tensor is
[math]\displaystyle{ \mathbf{C}(t) = \left[\begin{array}{cccccc} K_T(t)+G_T(t) & K_T(t)-G_T(t) & \ell(t) & 0 & 0 & 0 \\ K_T(t)-G_T(t) & K_T+G_T & \ell(t) & 0 & 0 & 0 \\ \ell(t) & \ell(t) & n(t) & 0 & 0 & 0 \\ 0 & 0 & 0 & G_A(t) & 0 & 0 \\ 0 & 0 & 0 & 0 & G_A(t) & 0 \\ 0 & 0 & 0 & 0 & 0 & G_T(t) \end{array}\right] }[/math]
Here [math]\displaystyle{ K_T(t) }[/math] is the plane strain, bulk modulus, [math]\displaystyle{ G_T(t) }[/math] is the transverse shear modulus, [math]\displaystyle{ G_A(t) }[/math] is the axial shear modulus, and [math]\displaystyle{ n(t) }[/math] and [math]\displaystyle{ \ell(t) }[/math] give time-dependence of the C11 and C12=C13 elements of the stiffness matrix (as defined here).