/PROP/TYPE6 (SOL_ORTH)

This property is only available for 8-node linear solid elements (/BRICK), tetrahedron elements (/TETRA4 and /TETRA10), and 2D solid elements (/QUAD). Quadratic bricks (/BRIC20 and /SHEL16) and pentahedron elements (/PENTA6) are not compatible with this property.

Format

Definitions

Example 1

Uses skew

Figure 1.

#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
unit for prop
#              MUNIT               LUNIT               TUNIT
                  kg                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/SKEW/FIX/1
New SKEW 1
#                 OX                  OY                  OZ
                   0                 100                   0
#                 X1                  Y1                  Z1
                   1                   0                  -1
#                 X2                  Y2                  Z2
                   0                   1                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  2. GEOMETRICAL SETS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/PROP/SOL_ORTH/1/2
SOL_ORTH example
#   Isolid    Ismstr               Icpre  Itetra10     Inpts   Itetra4    Iframe                  dn
        14         0                   1         0         0         0         0                   0
#                q_a                 q_b                   h
                   0                   0                   0
#                 Vx                  Vy                  Vz   skew_ID        Ip     Iorth
                   0                   0                   0         1         0         0
#                phi
                   0
#             dt_min
                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

Example 2

Uses I_p =1 and angle $\varphi$ to get same material direction (fiber direction) m1.

Figure 2.

#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
unit for prop
#              MUNIT               LUNIT               TUNIT
                  kg                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/PROP/SOL_ORTH/1/2
SOL_ORTH example
#   Isolid    Ismstr               Icpre  Itetra10     Inpts   Itetra4    Iframe                  dn
        14         0                   1                   0                   0                   0
#                q_a                 q_b                   h
                   0                   0                   0
#                 Vx                  Vy                  Vz   skew_ID        Ip     Iorth
                   0                   0                   0         0         1         0
#                phi
                  45
#             dt_min
                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

Comments

1. I_solid - Solid elements formulation
   • I_solid =17, brick deviatoric behavior is computed using 8 Gauss points, but the bulk behavior can be chosen with I_cpre, and compatible with all solid type material laws.
   • I_solid =24 (HEPH) solid elements use a physical hourglass formulation that is similar the hourglass formulation used by I_shell =24 (QEPH) shell elements. This hourglass formulation gives better results than the viscous hourglass formulation used by I_solid = 1 or 2.
   • I_solid =14 (HA8) is locking-free general solid formulation. Example: I_npts =222 is an 8 Gauss integration points solid. HA8 formulation is compatible with all orthotropic and isotropic material laws.
   • I_solid =18, the I_cpre and I_smstr default values depend on the material and are recommended values:

     Default	Material Laws
     Icpre = 2	2, 21, 22, 23, 24, 27, 36, 52, 79, 81, 84
     Icpre = 3	12, 14, 15, 25, 28, 50, 53, 68, and If $\nu \le 0.49$ , then 1, 13, 16, 33, 34, 35, 38, 40, 41, 70 and 77
     Icpre =1	All other laws and If $\nu \ge 0.49$ , then 1, 13, 16, 33, 34, 35, 38, 40, 41, 70, and 77
     Ismstr = 10	38, 42, 62, 69, 82, 88, 92, 94, 95
     Ismstr = 11	70
     Ismstr = 1	28
     Ismstr = 2	All other laws

2. When using the automatic setting option I_smstr = I_cpre = I_frame=-1, the values for these options are defined using the best options based on the element formulation, element type, and material. Alternatively, defining I_smstr = I_cpre = I_frame=-2 will overwrite the values for these options defined in this property with the best value based on element type and material law. To see the values defined by Radioss, review the “PART ELEMENT/MATERIAL PARAMETER REVIEW” section of the Starter output file.
3. Small strain:
   • If the small strain option is set (I_smstr=1 or 3), the strains and stresses used in material laws and output to time history and animation result files are engineering strains and stresses. Otherwise, they are true strains and stresses.
   • The Radioss Engine option /DT/BRICK/CST will only work for brick property sets with I_smstr =2 and 12.
   • The flag I_smstr =10 and 12 are only compatible with material LAW28 which uses total strain formulation.
   • I_smstr=12 is compatible with /DT/BRICK/CST, and total strain will be switched to small total strain, but not like the case of I_smstr=2, there is slight discontinuity of stresses during the passage.
   • Starting with version 2017, Lagrangian elements whose volume becomes negative during a simulation will automatically switch strain formulations to allow the simulation to continue. When this occurs, a WARNING message will be printed in the Engine output file. The following options are supported.

     Element Type and Formulation	Strain Formulation	Negative Volume Handling Method
     /BRICK Isolid=1, 2, 14, 17, 24 /TETRA4 /TETRA10	Full geometric nonlinearities. Ismstr = 2, 4.	Switch to small strain using element shape from cycle before negative volume.
     	Lagrange type total strain . Ismstr = 10, 12.	Lagrange type total strain with element shape at time=0.0.

4. I_cpre - Constant pressure formulation flag
   • I_cpre =1 is used to prevent volumetric locking in incompressible or quasi-incompressible material. For this case, the stress tensor is decomposed into a spherical and deviatoric part. Reduced integration is then used for the spherical part so that the pressure remains constant.
   • I_cpre =2 is only available for elasto-plastic laws. To prevent volume locking, additional terms with Poisson’s coefficient are added to the strain. When in the material is still elastic and thus compressible, the Poisson’s coefficient terms are small. As the material becomes plastic and thus incompressible, the Poisson’s coefficient terms increase to prevent volume locking. Refer to the Radioss Theory Manual for additional explanation.
5. Co-rotational formulation:

   For I_solid =1 or 2, and I_frame =2, the stress tensor is computed in a co-rotational coordinate system. This formulation is more accurate if large rotations are involved, at the expense of higher computation cost. It is recommended in case of elastic or visco-elastic problems with important shear deformations. Co-rotational formulation is compatible with 8 node bricks. Co-rotational formulation is also compatible with bi-dimensional and axisymmetric analysis (/QUAD) element).

6. d_n - Numerical damping and h - hourglass viscosity coefficient
   • Numerical damping d_n is used in the hourglass stress calculation for I_solid=24 (HEPH) solid elements.
   • When comparing results between I_solid=24 and I_solid =1 or 2 where d_n=h, the numerical damping is $\left(2/3\right)\times {10}^{-3}$ times smaller for I_solid =24 than I_solid =1 or 2.
7. Output for post-processing
   • For post-processing solid element stress, refer to /ANIM/BRICK/TENS for animation and /TH/BRICK for plot files.
   • In animation file, if elements are using co-rotational formulation,
     • I_solid =14, 24 (co-rotational frame always used)
     • I_frame =2 with I_solid =1, 2, 17

     Then the stress output is represented in material orthotropic coordinate system defined in /PROP/SOL_ORTH. Otherwise (co-rotational formulation not used), then the stress components (SIGX, SIGY, SIGZ, SIGXY, SIGYZ, and SIGXZ) are expressed in global coordinate system.

   • In plot files, the stress components SX, SY, SZ, SXY, SYZ, and SXZ are expressed in the global frame and the stress tensors components LSX, LSY, LSZ, LSXY, LSYZ, and LSXZ are expressed in the orthotropic frame (refer to /TH/BRIC for post-processing solid element stress in plot files).
8. Isoparametric systems (r-s-t)
   For 8 node bricks (I_solid =1 or 2), 4-node tetrahedron and 10-node tetrahedron, the orthotropic system rotates like the orthogonalized isoparametric system. Attention must be paid to the orientation of the orthotropic system in case of large shear.

   
   
   Figure 3.

   r, s, t: isoparametric frame

   r: center of (1, 2, 6, 5) to center of (4, 3, 7, 8)

   s: center of (1, 2, 3, 4) to center of (5, 6, 7, 8)

   t: center of (1, 4, 8, 5) to center of (2, 3, 7, 6)

   
   
   Figure 4.

9. Orthotropy direction
   For 3D solid elements, there are three different ways to define orthotropy direction:

   • With I_p= 0 and skew_ID ≠ 0, skew is used.

     Then no reference plane is used; skew is taken directly as the orthotropic system (in this case r =x, s =y, and t =z). x-direction is orthotropy direction 1 and y-direction is orthotropy direction 2.

   • With I_p = 1, 2 or 3, orthogonalized isoparametric system (r’-s’-t’) and orthotropic angle are used:
     The orthotropic system initial orientation (1-2-3) is defined with respect to the initial orthogonalized isoparametric system (r’-s’-t’) , as:

     
     
     Figure 5.

     
     
     Figure 6.

     In this case, the orthotropic system initial orientation is defined the same way as for bricks, I_solid = 0, 1 or 2 (that is with respect to the orthogonalized isoparametric system), and knowledge of the co-rotational system orientation is unnecessary to input the orthotropic system initial orientation.

   • With I_p = 11, 12 or 13, orthogonalized isoparametric system (r’-s’-t’) and reference vector $V$ are used:

     In this case, Global vector $V$ may be used to define the orthotropy direction. If this reference vector is orthogonal to plane, first axis of the plane is taken as orthotropy direction.

   For 2D solid elements (/QUAD), there are two different ways to define orthotropy direction:

   • With I_p = 1, isoparametric system (r-s) and orthotropic angle $\varphi$ are used. 10
   • With I_p = 11, Reference vector $V$ is used: Global vector $V$ may be used to define the orthotropy direction. If this reference vector is orthogonal to plane, first axis of the plane (r-s) is taken as orthotropy direction.
10. $\varphi$ - orthotropic angle

    For 2D solid elements (/QUAD):

    Orthotropic angle $\varphi$ is defined with respect to the first direction of the isoparametric frame (r-s).

    
    
    Figure 7.

    For 3D solid elements:

    If I_frame =2 (co-rotational formulation), the orthotropic system rotates like the co-rotational system. A co-rotational system is an orthogonalization of isoparametric systems (r-s-t) that has the same orientation whatever the permutation of r, s, t.

11. Solid to SPH properties (Sol2SPH)
    • When using Sol2SPH, solid elements are converted to SPH particles when a solid is deleted due to contact, a material failure criteria or time step criteria.
    • The number activated of SPH particles depends on parameter Ndir. The particles properties are computed using the sphpart_ID part number.
    • Skew definition is not required in the SPH property as the skew definition and orientation is automatically transmitted from the solid to the particles. It is not advised to use the same SPH part ID for an isotropic and orthotropic Sol2SPH part.

    The option Sol2SPH is only compatible with I_solid = 1, 2 or 24, I_frame = 1 or 2.

12. The I_solid flag is not used with 4-node (/TETRA4) or 10-node (/TETRA10) tetrahedron elements.

    4-node tetrahedron with I_tetra4 = 1 and 10-node tetrahedron I_tetra10 = 2 are compatible with all small strain formulation I_smstr.

    4-node tetrahedron with I_tetra4 = 1000 and all 10-node tetrahedron I_tetra10 = 2, 1000 are compatible with I_smstr = 10, 11 and 12.