/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}$		`Vdef_min`		`Vdef_max`		`APS_max`		`COL_min`

Definition

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)
$ϕ {}_{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)
$Δ t_{\min}$	Minimum time step. Default = 10⁶ (Real)	$[s]$
`Vdef_min`	Minimum volume ratio (V/Vo) to delete solid element. Default = 0
`Vdef_max`	Maximum volume ratio (V/Vo) to delete solid element. Used only if different from 0. Default = 0
`ASP_max`	Maximum aspect ratio to delete solid element. Used only if different from 0. Default = 0
`COL_min`	Minimum collapse value to delete solid element. Default = 0

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                              
#         deltaT_min            vdef_min            vdef_max             ASP_max             COL_min
                   0                   0                   0                   0                   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 (refer to /DEF_SOLID) 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.