/PROP/TYPE22 (TSH_COMP)

Block Format Keyword This property set is used to define the composite thick shell property set.

Format

(1)	(2)	(3)	(5)	(6)	(7)	(8)	(9)
/PROP/TYPE22/`prop_ID`/`unit_ID` or /PROP/TSH_COMP/`prop_ID`/`unit_ID`
`prop_title`
`I`_solid	`I`_smstr		`I`_cstr	`I`_npts	`I`_int		`d`_n
`q`_a		`q`_b
`V`_X		`V`_Y	`V`_Z		`skew_ID`	`I`_orth	`I`_pos
`A`_shear

For each layer (integration point) per line:

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
$ϕ {}_{i}$		`t`_i/`t`		`Z`_i		`mat_ID`_i

Last card:

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
$Δ t_{\min}$

Definitions

Field	Contents	SI Unit Example
`prop_ID`	Property identifier. (Integer, maximum 10 digits)
`unit_ID`	Unit Identifier. (Integer, maximum 10 digits)
`prop_title`	Property title. (Character, maximum 100 characters)
`I`_solid	Solid elements formulation flag. = 0 Use value /DEF/SOLID. = 14 HA8 locking-free 8-node thick shell, co-rotational, full integration, variable number of Gauss points in all directions. = 15 HSEPH/PA6 thick shell (8-node and 6-node respectively), co-rotational, under integrated (1-point in-plan quadrature) with physical stabilization, variable number of integration points in thickness direction. (Integer)
`I`_smstr	Small strain formulation flag. 4 = -1 Automatically define the best value based on element type and material law. = 0 (Default) Used value in /DEF_SOLID. = 1 Small strain from time =0. = 2 Full geometric nonlinearities with possible small strain formulation in Radioss Engine (/DT/Eltyp/Keyword3/Iflag). = 3 Simplified small strain formulation from time =0 (non-objective formulation). = 4 Default, if /DEF_SOLID is not defined Full geometric nonlinearities (/DT/BRICK/CST has no effect). (Integer)
`I`_cstr	Constant stress formulation flag. Only valid for `I`_solid = 14. = 001 Reduced stress integration in t direction. = 010 Reduced stress integration in s direction. = 100 Reduced stress integration in r direction. (Integer)
`I`_npts	Number of integration points. 2 = j 1 ≤ j ≤ 200 for `I`_solid =15 = ijk 2 ≤ i,j,k ≤ 9 for `I`_solid =14 (Integer) Where, i Number of integration points in r direction. j Number of integration points in s direction. k Number of integration points in t direction.
`I`_int	Number of layers when 9 < number of layers ≤ 200. Only valid for `I`_solid = 14. 4 (Integer)
`d`_n	Numerical damping for stabilization. Only valid for `I`_solid = 15. Default = 0.1 (Real)
`q`_a	Quadratic bulk viscosity. Default = 1.10 (Real) Default = 0.0 for /MAT/LAW70
`q`_b	Linear bulk viscosity. Default = 0.05 (Real) Default = 0.0 for /MAT/LAW70
`A`_shear	Shear factor. Default = 1.0 (Real)
`V`_X	X component for reference vector. Default = 1.0 (Real)
`V`_Y	Y component for reference vector. Default = 0.0 (Real)
`V`_Z	Z component for reference vector. Default = 0.0 (Real)
`skew_ID`	Skew identifier. If the local skew has been defined, its X-axis replaces the reference vector (`V`_X, `V`_Y, and `V`_Z will be ignored). (Integer)
`I`_orth	Orthotropic system formulation flag for reference vector. = 0 (Default) The first axis of orthotropy is maintained at constant angle with respect to the orthonormal co-rotational element coordinate system. = 1 The first orthotropy direction is constant with respect to a non-orthonormal isoparametric coordinates. (Integer)
`I`_pos	Layer positioning flag for reference vector. = 0 (Default) Layer positions are automatically calculated with regard to layer thicknesses partition. The coherence of global thickness with the sum of layer thicknesses is automatically checked. = 1 All layer positions in the element thickness are user-defined. Multiple layers may have the same special position. (Integer)
$Δ t_{\min}$	Minimum time step. Default = 10⁶ (Real)	$[s]$
$ϕ {}_{i}$	Angle for layer I. (Real)	$[\deg]$
`t`_i/`t`	Relative thickness of layer i. `t`_i The thickness of `i_th` layer. `t` The total thickness. (Real)
`Z`_i	Z position (normalized by the thickness) of layer i (-0.5 ≤ `Z`_i ≤ 0.5). Default = 0.0 (Real)
`mat_ID`_i	Material identifier for layer I. (Integer)

Example

3 layers (I_npts=333), each layer defining different material direction (fiber direction) m1 with vector

V

and angle

ϕ

. Plane (1', 2', 3', and 4') is middle surface of thick shell element (where, z=0).

#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----|
#-  2. GEOMETRICAL SETS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/PROP/TYPE22/1/2
TSH_COMP example
#   Isolid    Ismstr                         Icstr     Inpts      Iint                            dn
        14         0                           010       333         0                             0
#                q_a                 q_b
                   0                   0
#                 Vx                  Vy                  Vz   skew_ID     Iorth      Ipos
                   1                  -1                   1         0         0         0
#             Ashear
                   0
#              PHI_I               T_I/T                  ZI     MAT_I                             
                  45                 0.3                   0         1 
                  90                 0.4                   0         2 				  
                 -45                 0.3                   0         1                              
#             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

I_solid - Solid formulation
- I_solid =14 formulation (H8 element) must use constant stress formulation (I_cstr > 0), which refers to local isoparametric orthogonolized system r-s-t. Definition of the system is described in the comments of /PROP/TYPE6 (SOL_ORTH).
- When using I_solid=15 with pentahedron elements, /PENTA6 elements are recommended but degenerated /BRICK elements can also be used.
Number of layers.
- For I_solid = 14 formulation (HA8 element), number of layers ( < 9 ) is defined as:
  - If I_cstr = 001, the number of layers in t direction is equal to k value from I_npts field.
  - If I_cstr = 010, the number of layers in s direction is equal to j value from I_npts field (I_cstr = 010; I_npts = 282; for a number of 8 layers in s direction).
  - If I_cstr = 100, the number of layers in r direction is equal to i value from I_npts field.
- For I_solid formulation (HA8 element) when the number of layers > 9 is defined as:
  - Use I_int for I_solid formulation (HA8 element) when the number of layers > 9.
    In this case, the thickness direction integration points defined by I_npts should be zero.
    
    Example, I_cstr = 010; I_npts = 202; I_int = 100 for a number of 100 layers in "s" direction.
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.

I_smstr - Small strain formulation flag.

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	Element Formulation	Strain Formulation	Negative Volume Handling Method
/BRICK	`I`_solid =14, 15	Full geometric nonlinearities `I`_smstr = 2, 4	Switch to small strain using element shape from cycle before negative volume.

Othotropy in local coordinate system.
- The thick shell orthotropy is planar and the third orthotropy direction is coincident with the normal to the shell plane.
- Global vector $V$ or skew_ID is used to define the othotropy direction. The global vector $V$ or the $X$ -axis of specified skew (in this case, global vector is ignored) is projected to the mean plane of soslid element.
- For I_solid=14, the mean plane of the element depends on I_cstr.
  - r-s for I_cstr=001
  - r-t for I_cstr=010
  - s-t for I_cstr=001.
- $ϕ$ is the angle (in degrees) between the first direction of orthotropy and projection of reference vector on the shell mean plane for layer i.
Material used for layer
- Material law type used in Mat_ID_i can be different for each layer.
- For I_solid= 15, the material law number defined in /PART will be used to compute the contact interface stiffness and the hourglass stresses.