/PROP/TYPE51

Description

The input can be either "substack" or by ply, but not both.

Format

Definition

Comments

1. The stack property (TYPE51) is used in combination with /PROP/TYPE19 (PLY) to create composites properties through the ply-based definition.
   • Each ply in the stack is defined using Pply_ID_i or Pply_ID_j and can only be referenced once in the stack. If ply ids are used more than once an error message will output.
   • Each ply of a shell element can contain a different (≥ 1) number of integration points through its ply thickness. That is, Npt_ply =3 in /PROP/TYPE19 (PLY)
   • Different integration scheme (uniformly distribution or Gauss integration scheme) could be defined with option I_int inside this property.
2. I_shell, I_sh3n – 4-node and 3-node shell formulation flag
   • I_shell =1,2,3,4 (Q4): original 4 nodes Radioss shell with hourglass perturbation stabilization.
   • I_shell =24 (QEPH): formulation with hourglass physical stabilization for general use.
   • I_shell =12 (QBAT): modified BATOZ Q4γ24 shell with four Gauss integration points and reduced integration for in-plane shear. No hourglass control is needed for this shell.
   • I_sh3n =30 (DKT18): BATOZ DKT18 thin shell with three Hammer integration points.
   • The hourglass formulation is visco-elastic for Q4 shells.
3. I_smstr - Small strain formulation
   • Small strain formulation is activated from time t = 0, if I_smstr =1. It may be used for a faster preliminary analysis, but the accuracy of results is not ensured. Any shell for which $\text{Δ}t<\text{Δ}{t}_{\min }$ can be switched to a small strain formulation by Radioss Engine option /DT/SHELL/CST, except if I_smstr=4.
   • If I_smstr =1, the strains and stresses which are given in material laws are engineering strains and stresses; otherwise they are true strains and stresses.
4. h_m, h_f, and h_r- Hourglass coefficients
   • h_m, h_f, and h_r are only used for Q4 shells. They must have a value between 0 and 0.05.
   • For I_shell=3, default values of h_m and h_r are 0.1 with larger values possible.
5. d_n - Shell numerical damping coefficient
   It is only used for I_shell =12 or 24:

   • for I_shell =24 (QEPH),d_n is used for hourglass stress calculation.
   • for I_shell =12 (QBAT), d_n is used for all stress terms, except transverse shear.
   • for I_sh3n =30 (DKT18), d_n is only used for membrane.
6. I_thick - Shell resultant stresses calculation flag
   • If I_thick= 1, the small strain option is automatically deactivated in the corresponding type of element.
7. I_drill - Drilling degree of freedom stiffness flag
   • Drilling DOF stiffness is recommended for implicit solutions especially for Riks method and bending dominated problems.
   • I_drill is available for QEPH, QBAT (I_shell =12, 24), and standard triangle (C0) shell elements (I_sh3n = 1, 2).
8. Anisotropy direction definition.
   The reference vector $V$ is defined as following according flag I_P:

   • If I_P=0 and skew_ID = 0, the reference vector $V$ is defined with VX, VY and VZ.
   • If I_P=0 and skew_ID ≠ 0, the reference vector $V$ is the first direction (local X) of the local coordinate system skew_ID.
   • If I_P = 20, the reference vector $V$ is defined with the node N1 and N2 of the shell elements.
   • If I_P = 22, the reference vector $V$ is the first direction (local X) of the local coordinate system skew_ID. Vector components VX, VY and VZ are ignored.
   • If I_P = 23, the reference vector $V$ is defined with VX, VY and VZ. Local coordinate system skew_ID is ignored.
   The reference vector $V$ is projected on the shell element plane and becomes the vector $V^{\prime }$ . Then for each layer, the 1^st material direction (m1) is vector $V^{\prime }$ turned ${\varphi }_i$ degrees (turns positive direction around shell normal $n$ ).

   
   
   Figure 3.

   The hierarchy order to define the reference vector $V$ is:

   • shell element ( ${\varphi }_s$ > 0 in either /SHELL or /SH3N)
   • initial state card (/INISHE/ORTHO)
   • shell property

   In case of reference metrics, the orientation for directions of anisotropy must be defined with the reference geometry, not the initial one. The 2^nd material direction m2 is derived from direction m1 with the angle ${\alpha }_i$ (defined in /PROP/TYPE19). If ${\alpha }_i=90°$ , the layer is orthotropic.

9. I_pos - ply position
   • In this property, ply position defined with Z_i and option I_pos .
   • Z_i values are real ply positions in the local Z axis (negative and positive value are allowed and Z_i defines the position of the middle of the each ply).
   • Z₀ only used with I_pos=2, it indicated distance from the shell element reference middle surface to the bottom surface of the whole ply or stack. More info see below I_pos=2.
   • I_pos = 0: ply positions are automatically calculated with regard to ply thicknesses.
     • The input thickness “Thick” may be equal to zero as the true shell thickness is always calculated from the ply thicknesses and positions of ply defined by Z_i.
     • If the input thickness “Thick” is different from 0.0, a warning will be displayed if the calculated thickness from the plies is different from the input thickness.

       
       
       Figure 4.

   • I_pos = 1: All ply positions and thickness in the element are user-defined (with in /PROP/PLY and Z_i).
     • Global thickness “Thick” is not checked, as it does not need to be equal to the sum of ply thickness.
     • Multiple plies are allowed to have the same space position.
     • This option is not recommended unless the number of ply is constant for all elements.
   • I_pos = 2: the shell element mid-surface is at Z₀ from the bottom of the ply layout.

     
     
     Figure 5.

   • I_pos = 3: the top of the ply layout is coincident with the element mid-surface.

     
     
     Figure 6.

   • I_pos = 4: the bottom of the ply layout is coincident with the element mid-surface.

     
     
     Figure 7.

10. Define by substack
    • If the substack definition is used, each substack definition begins with the keyword “SUB” and the substack identification number.
    • All the plies until the next keyword “SUB” or the keyword “INT” is encountered belong to the same substack.
    • All substacks must be defined before defining the “INT connections
    • If an element contains plies belonging to several substacks, the connection between the substacks is defined with an INT connection: the ply Pply_ID_t of the first substack is connected to the ply Pply_ID_b of the second substack.
    • For more detail, refer to above examples for this property.
11. Material for each ply
    • Material for each ply defined with mat_ID in /PROP/TYPE19 (PLY).
    • Different material laws could be defined for each ply.
    • The material law number given in /PART will be used to define the mass and the sound speed of the composite, as well as the interface stiffness.
    • This property is compatible with XFEM (crack propagation) using /FAIL/JOHNSON and /FAIL/TBUTCHER.
12. P_thick_fail, P_thickl_{fail_i} and F_weight_i parameters are not compatible with failure defined within the material law itself, such as plastic failure strain in LAW36.
13. Element suppression rules used with /FAIL models:
    • A single Gauss point in each layer is deleted when relative thickness of failed integration points is greater than prescribed value P_thickl_{fail_i} for this layer:

      ${\displaystyle {\sum }_{i}Thick{l}_i}\ge P\_thick{l}_{fail\_i}$ with $i=1,N_i$

      Where,

      $N_i$

      The number of in-thickness integration points in the layer, defined in corresponding PLY property.

      $Thick{l}_i$

      The ratio of thickness of integration point versus layer thickness.

    • Each single layer is turned OFF when all in-plane Gauss integration points in the layer are deleted.
    • The shell element is deleted when the following criterion is met:

      ${\left[{{\displaystyle {\sum }_i\left(Thic{k}_i\cdot F\_weigh{t}_i\right)}}^{F_{exp}}\right]}^{\frac{1}{F_{exp}}}\ge P\_thic{k}_{fail}$ with $i=1,N_i$

      Where,

      $Thic{k}_i$

      Relative thickness of failed layer i over total element thickness.

    • By default, the element failure parameters defined in the property are used and local Ifail_sh settings defined at the failure model level are ignored. If the P_thickl_{fail_i} parameters per ply are not defined, the element failure settings defined in the local failure models are applied to the corresponding layer. The P_thick_fail parameter must be defined in the property otherwise the default value is used.