/PROP/TYPE22 (TSH

Block Format Keyword

/PROP/TYPE22 - Composite Thick Shell Property Set

Description

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

Format

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
/PROP/TYPE22/prop_ID/unit_ID or /PROP/TSH_COMP/prop_ID/unit_ID
prop_title
Isolid	Ismstr			Icstr	Inpts	Iint		dn
qa		qb
VX		VY		VZ		skew_ID	Iorth	Ipos
Ashear

For each layer (integration point):

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
i		ti/t		Zi		mat_IDi

Last card:

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)

Field	Contents	SI Unit Example
prop_ID	Property identifier (Integer, maximum 10 digits)
unit_ID	Optional unit identifier (Integer, maximum 10 digits)
prop_title	Property title (Character, maximum 100 characters)
Isolid	Solid elements formulation flag (Integer) = 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.
Ismstr	Small strain formulation flag (Integer) = 0: default, set to value defined in /DEF_SOLID = 1: small strain from time =0 = 2: full geometric nonlinearities with possible small strain formulation in RADIOSS Engine (/DT/BRICK/CST) = 3: simplified small strain formulation from time =0 (non-objective formulation) = 4: full geometric nonlinearities (/DT/BRICK/CST has no effect)
Icstr	Constant stress formulation flag (Isolid = 14 only) (Integer) = 001: reduced stress integration in t direction = 010: reduced stress integration in s direction = 100: reduced stress integration in r direction
Inpts	Number of integration points (Comment 2) (Integer) = j: 1 < j < 200 for Isolid =15 = ijk: 2 < i,j,k < 9 for Isolid =14 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
Iint	Number of layers when 9 < number of layers < 200 (Isolid = 14 only) (Comment 3) (Integer)
dn	Numerical damping for stabilization (Isolid = 15 only) Default = 0.1 (Real)
qa	Quadratic bulk viscosity Default = 1.10 (Real) Default = 0.0 for /MAT/LAW70
qb	Linear bulk viscosity Default = 0.05 (Real) Default = 0.0 for /MAT/LAW70
Ashear	Shear factor Default = 1.0 (Real)
VX	X component for reference vector Default = 1.0 (Real)
VY	Y component for reference vector Default = 0.0 (Real)
VZ	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 (VX, Vy, and VZ will be ignored). (Integer)
Iorth	Orthotropic system formulation flag for reference vector Default = 0 (Integer) = 0: 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.
Ipos	Layer positioning flag for reference vector Default = 0 (Integer) = 0: 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.
	Minimum time step Default = 106 (Real)
1	Angle for layer i (Real)
ti/t	Relative thickness of layer i. where, ti is the thickness of i_th layer, and t is the total thickness (Real)
Zi	Z position (normalized by the thickness) of layer i (-0.5 < Zi < 0.5) Default = 0.0 (Real)
mat_IDi	Material identifier for layer i (Integer)

prop_type22_example

3 layers (Inpts=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

# 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

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

#enddata

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

1.	Isolid =14 formulation (H8 element) must use constant stress formulation (Icstr > 0), which refers to local isoparametric orthogonolized system r-s-t. Definition of the system is described in the comments of /PROP/TYPE6.

2.	For Isolid = 14 formulation (HA8 element), number of layers ( < 9 ) is defined as follows:

If Icstr = 001, the number of layers in t direction is equal to k value from Inpts field

If Icstr = 010, the number of layers in s direction is equal to j value from Inpts field
(Icstr = 010; Inpts = 282; for a number of 8 layers in s direction)

If Icstr = 100, the number of layers in r direction is equal to i value from Inpts field

3.	Use Iint for Isolid formulation (HA8 element) when the number of layers > 9. In this case, the thickness direction integration points defined by Inpts should be zero. (example, Icstr = 010; Inpts = 202; Iint = 100 for a number of 100 layers in s direction).

4.	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 reference orthotropy direction. The global vector or the X-axis of specified skew (in this case global vector is ignored) is projected to the mean plane of solid element. For Isolid=14, the mean plane of the element depends on flag Icstr: r-s for Icstr=001, r-t for Icstr=010 and s-t for Icstr=001.

6.	i is the angle (in degrees) between the first direction of orthotropy this projection of reference vector on the shell mean plane (r and t) for layer i.

7.	Material law type can be different for each layer.

8.	The material law number defined in part definition will be used to compute the contact interface stiffness and the hourglass stresses (Isolid = 15 only).

9.	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

Isolid = 14, 15

Full geometric nonlinearities

Ismsr = 2, 4

Switch to small strain using element shape from cycle before negative volume.

/PROP/TYPE22 (TSH_COMP)

/PROP/TYPE22 (TSH_COMP)

/PROP/TYPE22 (TSH_COMP)