/PROP/TYPE16 (SH

Block Format Keyword

/PROP/TYPE16 - Anisotropic Layered Shell Property Set

Description

This property set is used to define the anisotropic layered shell property set. This property is currently only compatible with Elastic Anisotropic Fabric (/MAT/LAW58 - FABR_A) and only one layer is allowed. For compatibility with shell elements, only standard and Batoz fully integrated element are available:

➢	Q4: original 4 nodes RADIOSS shell with hourglass perturbation stabilization.

➢	QBAT: modified BATOZ Q4γ24 shell with 4 Gauss integration points.

Format

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
/PROP/TYPE16/prop_ID/unit_ID or /PROP/SH_FABR/prop_ID/unit_ID
prop_title
Ishell	Ismstr	Ish3n
hm		hf		hr		dm
N	Istrain	Thick		Ashear			Ithick
VX		VY		VZ		skew_ID	Ipos
				ti		Zi		mat_IDi

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)
Ishell	Shell element formulation flag (Integer) = 0: use value in /DEF_SHELL = 1: Q4, visco-elastic hourglass modes orthogonal to deformation and rigid modes (Belytschko) = 2: Q4, visco-elastic hourglass without orthogonality (Hallquist) = 3: Q4, elasto-plastic hourglass with orthogonality = 4: Q4 with improved type 1 formulation (orthogonalization for warped elements) = 12: QBAT shell formulation
Ismstr	Shell small strain formulation flag (Integer) = 0: use value in /DEF_SHELL (only larger strain allowed) = 2: full geometric nonlinearities with possible small strain formulation activation in RADIOSS Engine (option /DT/SHELL/CST) = 4: full geometric nonlinearities (/DT/SHELL/CST) =11: total small strain from time = 0
Ish3n	3 node shell element formulation flag (Integer) = 0: use value in /DEF_SHELL = 1: standard triangle (C0) = 2: standard triangle (C0) with modification for large rotation = 30: DKT18 = 31: DKT_S3
hm	Shell membrane hourglass coefficient Default = 0.01 (Real)
hf	Shell out-of-plane hourglass Default = 0.01 (Real)
hr	Shell rotation hourglass coefficient Default = 0.01 (Real)
dm	Shell membrane damping (Real)
N	Number of layers, with (Comment 12) Default set to 1 (Integer)
Istrain	Compute strains for post-processing flag (Integer) = 0: default set to value defined with /DEF_SHELL = 1: yes = 2: no
Thick	Shell thickness (Real)
Ashear	Shear factor Default is Reissner value: 5/6 (Real)
Ithick	Shell resultant stresses calculation flag (Integer) = 0: default set to value defined with /DEF_SHELL = 1: thickness change is accounted for = 2: thickness is constant (Comment 6)
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 for reference vector If the local skew is defined, its X-axis replaces the global vector V. VX, VY, and VZ coordinates are ignored. Default = 0 (Integer)
Ipos	Layer positioning flag for reference vector (Comment 9) Default = 0 (Integer) = 0: Layer positions Zi are automatically calculated with regard to layer thicknesses. = 1: All layer positions Zi must be user-defined.
	Angle of 1st local axis for layer i (Real)
	Angle between first and second axis Default = 90.0 (Real)
ti	Thickness ti of layer i (Real)
Zi	Z position of layer i (Zi defines the position of the middle of the layer) Default = 0.0 (Real)
mat_IDi	Material identifier for layer i (Integer)

1.	Total small strain formulation is activated from time t= 0, if Ismstr = 11; this formulation has been developed to improve the robustness for Airbag models.

2.	Shell membrane damping dm is used for Material LAW58:

dm = 25% (default), if Ishell=12 (BATOZ) shell is used

otherwise, dm = 0

3.	The hourglass formulation is visco-elastic for Q4 shells.

4.	hm, hf, and hr are used only for Q4 shells. They must have a value between 0 and 0.05.

5.	For hourglass type 3 (Ishell= 3), hourglass maximum values may be larger, default values are 0.1 for hm and hr.

6.	Input components of global vector used to define direction 1 of local coordinate system. Alternatively, it may be defined by a local skew system.

prop_type9_v10

7.	is the angle between local X-axis and projection of vector on the shell for layer i.

8.	is the angle between local axes X and Y (directions of anisotropy) for layer i.

9.	Ipos = 0: layer positions are calculated automatically. If Thick is not equal to the sum of layer thicknesses, a warning message is displayed, and individual ply thickness will be adjusted such that the sum of the ply thicknesses will equal the total thickness which specified on Thick.

Ipos_0

Ipos = 1: all layer positions in the element thickness are user-defined. Multiple layers are allowed to have the same space position. In this case Thick is not checked, as it does not need to be equal to the sum of layer thickness.

For more details, refer to Layer thickness and layer position calculation in FAQs.

10.	The material law number given in the element input section will be used to define the mass and the sound speed of the composite, as well as the interface stiffness.

11.	In case of reference metrics, the orientation for directions of anisotropy must be defined with the reference geometry, not the initial one.

12.	Each i layer has only 1 integration point through its thickness.

/PROP/TYPE16 (SH_FABR)

/PROP/TYPE16 (SH_FABR)

/PROP/TYPE16 (SH_FABR)

Description