/MAT/LAW108 (SPR_GENE)

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  1. LOCAL_UNIT_SYSTEM:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/2
units for material and property
                  kg                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/PROP/TYPE23/6680004/2
SPOTWELD_NO_RUPTURE
#    Imass                          Area             Inertia   skew_ID   sens_ID    Isflag     
         2                          8E+3                .002         0         0         0 
/MAT/LAW108/6680004/2
TYPE8
#            Density
                1E-6
#    Ifail    Iequil    Ifail2
                   0             
#              KTens               CTens               ATens               BTens               DTens
                 1.8                   0                   0                   0                   0
#  fct_ID1     HTens   fct_ID2   fct_ID3   fct_ID4                 delta_minTens       delta_maxTens
                   0         0         0         0                             0                   0
#                  F                   E              Ascale             Hscalex
                   0                   0                   0                   0
#              KTens               CTens               ATens               BTens               DTens
                  .3                   0                   0                   0                   0
#  fct_ID1     HTens   fct_ID2   fct_ID3   fct_ID4                 delta_minTens       delta_maxTens
                   0         0         0         0                             0                   0
#                  F                   E              Ascale             Hscalex
                   0                   0                   0                   0
#              KTens               CTens               ATens               BTens               DTens
                 1.8                   0                   0                   0                   0
#  fct_ID1     HTens   fct_ID2   fct_ID3   fct_ID4                 delta_minTens       delta_maxTens
                   0         0         0         0                             0                   0
#                  F                   E              Ascale             Hscalex
                   0                   0                   0                   0
#                  K                   C                   A                   B                   D
    114.649681528662                   0                   0                   0                   0
#       N1         H        N2        N3        N4                     theta_min           theta_max
                   0         0         0         0                             0                   0
#                  F                   E              Ascale             Hscalex
                   0                   0                   0                   0
#                  K                   C                   A                   B                   D
    114.649681528662                   0                   0                   0                   0
#       N1         H        N2        N3        N4                     theta_min           theta_max
                   0         0         0         0                             0                   0
#                  F                   E              Ascale             Hscalex
                   0                   0                   0                   0
#                  K                   C                   A                   B                   D
    114.649681528662                   0                   0                   0                   0
#       N1         H        N2        N3        N4                     theta_min           theta_max
                   0         0         0         0                             0                   0
#                  F                   E              Ascale             Hscalex

#   ISTRAT             ASTRAT 
         0                   0
#enddata

1. When used with /PROP/TYPE23 (SPR_MAT), this material law has the same behavior as spring property /PROP/TYPE8 (SPR_GENE); except that in this material the mass is calculated from the density and volume, or density, area and length.
2. Inputs repeated for each degree of freedom (DOF) $i$ are defined with the following directions:
   • $i$ =1: tension/compression
   • $i$ =2: shear xy
   • $i$ =3: shear xz
   • $i$ =4: torsion
   • $i$ =5: bending y
   • $i$ =6: bending z
3. This spring should only be used when the spring has zero initial length or when one of the spring nodes is constrained in all directions. For other situations, the equilibrium is insured for forces, but not for moments. Thus, the moment calculation can be incorrect. If the initial spring length is not zero, use /PROP/TYPE13 (SPR_BEAM) or /MAT/LAW113 (SPR_BEAM).
4. The spring has six DOF ${\delta }^1,{\delta }^2,{\delta }^3,{\theta }^4,{\theta }^5,{\theta }^6$ computed in a local coordinate system. The local system can be defined for each element using /SPRING Skew_ID. If Skew_ID is not defined for each element, the local system defined in the /PROP/TYPE23 (SPR_MAT) Skew_ID is used. If a system is not defined in the element or property, then the global system is used.
5. The sign (tension or compression) of the spring force depends on the relative displacement of Node N2 with respect to Node N1. If the movement of node N2 relative to node N1 is in the positive direction of the spring system, then the spring is in tension. If the movement of node N2 relative to node N1 is in the negative direction of the spring system, then the spring is in compression.
6. If I_equil = 0 (no equilibrium), then:(1)
   $$\text{f}\left(\theta \right)=M_{2y}=-M_{1y}$$

   $M_{2y}$ is moment in $Y$ by N₂

   $M_{1y}$ is moment in $Y$ by N₁

   If I_equil = 1 (force and moment equilibrium), then:(2)
   $$-M_{1y}\ne M_{2y}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}-M_{1z}\ne M_{2z}$$
   (3)
   $$\text{f}\left(\theta \right)=\frac{M_{2y}-M_{1y}}{2}$$

   $M_{2y}$ is moment in $Y$ by N₂

   $M_{1y}$ is moment in $Y$ by N₁

   $M_{2z}$ is moment in $Z$ by N₂

   $M_{1z}$ is moment in $Z$ by N₁

7. Force and moment computation.
   • $\delta$ is a translational DOF, the force in direction $\delta$ is computed as:

     $\mathrm{F}\left(\delta \right)=\mathrm{f}\left(\frac{{\delta }^i}{Ascal{e}_i}\right)\left[A_i+B_i\ln \left|\frac{{\dot{\delta }}^i}{D_i}\right|+E_i\mathrm{g}\left(\frac{{\dot{\delta }}^i}{F_i}\right)\right]+C_i{\dot{\delta }}^i+Hscal{e}_i\mathrm{h}\left(\frac{{\dot{\delta }}^i}{F_i}\right)$ with $i$ =1,2,3

   • $\theta$ is a rotational DOF, the moment is computed as:

     $\mathrm{M}\left(\theta \right)=\mathrm{f}\left(\frac{{\theta }^i}{Ascal{e}_i}\right)\left[A_i+B_i\ln \left|\frac{{\dot{\theta }}^i}{D_i}\right|+E_i\mathrm{g}\left(\frac{{\dot{\theta }}^i}{F_i}\right)\right]+C_i{\dot{\theta }}^i+Hscal{e}_i\mathrm{h}\left(\frac{{\dot{\theta }}^i}{F_i}\right)$ with $i$ =4,5,6

     Where,

     • ${\delta }^i$ (with $-l_0<{\delta }^i<+\infty$ ) is the difference between the current length $l$ and the initial length $l_0$ of the spring element for corresponding translational DOF.
     • ${\theta }^i$ is the relative angle for corresponding rotational DOF in radians.
     • For linear springs, $\mathrm{f}\left(\delta \right),\text{}\mathrm{g}\left(\dot{\delta }\right)$ and $\mathrm{h}\left(\dot{\delta }\right),\left(\mathrm{f}\left(\theta \right),\mathrm{g}\left(\dot{\theta }\right)\text{and}\mathrm{h}\left(\dot{\theta }\right)\right)$ are zero functions and $A_i$ , $B_i$ , $E_i$ and $Hscal{e}_i$ are not taken into account.
     • If stiffness function $\mathrm{f}\left(\delta \right)$ or $\mathrm{f}\left(\theta \right)$ is requested, then $K$ is used as a slope for unloading only.
     • If $K$ is lower than the maximum slop of function $\mathrm{f}\left(\delta \right)$ or $\mathrm{f}\left(\theta \right)$ ( $K$ is not consistent with the maximum slope of the curve), $K$ is set to the maximum slope of the curve.
8. Note that material density and inertia in /PROP/TYPE23 must not be null to ensure that the nodes connected by such a spring get a non-zero mass and inertia, unless they are secondary nodes of a rigid body, in which case unrealistically high stiffness and damping values are avoided.

   If the nodes connected by such a spring get a zero mass and inertia, they can still be secondary nodes of a rigid body, except for rigid bodies activated by a sensor. The rigid body cannot be set to OFF in the Radioss Engine.

   It is not possible to use elementary time step for springs with null density or null inertia in the model. Nodal Time Step or Advanced Mass Scaling must be used. If control time step is not set to Nodal Time Step, neither Advanced Mass Scaling in the Radioss Engine nor Nodal Time Step will be turned ON automatically.