/PROP/TYPE11 (SH

Block Format Keyword

/PROP/TYPE11 - Sandwich Shell Property Set

Description

This property set is used to define the sandwich shell property set. It is possible to define sandwich composite with several layers and each lay with individual material, thickness, layer position and orthotropic direction. This property is only compatible with Material Laws 15, 25, 27, 36, 60, 72 and user laws.

Format

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
/PROP/TYPE11/prop_ID/unit_ID or /PROP/SH_SANDW/prop_ID/unit_ID
prop_title
Ishell	Ismstr	Ish3n	Idrill
hm		hf		hr		dm		dn
N	Istrain	Thick		Ashear			Ithick	Iplas
VX		VY		VZ		skew_ID	Iorth	Ipos

For each layer (integration point)

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
		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 = 24: QEPH shell formulation
Ismstr	Shell small strain formulation flag (Integer) = 0: use value in /DEF_SHELL = 1: small strain from time =0 (new formulation compatible with all other formulation flags) = 2: full geometric nonlinearities with possible small strain formulation activation in RADIOSS Engine (option /DT/SHELL/CST) = 3: old small strain formulation (only compatible with hourglass type 2) = 4: full geometric nonlinearities (in RADIOSS Engine, option /DT/SHELL/CST has no effect)
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, which based on DTK12 of BATOZ (refer to the Theory Manual)
Idrill	Drilling degree of freedom stiffness flag (Comment 23) (Integer) = 0: default set to 2 = 1: yes = 2: no
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)
dn	Shell numerical damping (Real)
N	Number of layers, with Default = 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 (Comment 21) (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 taken into account = 2: thickness is constant
Iplas	Shell plane stress plasticity flag (Integer) = 0: default set to value defined with /DEF_SHELL = 1: iterative projection with three Newton iterations = 2: radial return
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)
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 X-axis of an orthonormal co-rotational element coordinate system. = 1: the first orthotropy direction is constant with respect to a non-orthonormal system of deformed element.
Ipos	Layer positioning flag for reference vector (Comment 21) 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 for layer i (Real)
ti	Thickness of layer i (Comment 22) (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)

prop_type11_example

3 layers (N=3) with different material, different material direction (fiber direction). Reference vector is taken from x-axis of skew.

#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----|

/SKEW/FIX/1

New SKEW 1

# OX OY OZ

1.0 0 100.0

# X1 Y1 Z1

0 0 1

# X2 Y2 Z2

0 -1 0

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

#- 2. GEOMETRICAL SETS:

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

/PROP/SH_SANDW/2/2

SH_SANDW example

# Ishell Ismstr Ish3n Idril

12 0 0 0

# hm hf hr dm dn

0 0 0 .1 .1

# N Istrain Thick Ashear Ithick Iplas

3 0 1.6 0 1 1

# Vx Vy Vz Skew_ID Iorth Ipos

0 0 0 1 0 0

# Phi t Z mat_ID

45 .5 0 1

90 .6 0 2

-45 .5 0 1

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

#enddata

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

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

QEPH: formulation with hourglass physical stabilization for general use (istotropic + law25 shells only).

QBAT: modified BATOZ Q4 y24 shell with four Gauss integration points and reduced integration for in-plane shear. No hourglass control is needed for this shell.

DKT18: BATOZ DKT18 thin shell with three Hammer integration points.

2.	Small strain formulation is activated from time t =0, if Ismstr = 1 or 3. It may be used for a faster preliminary analysis, but the accuracy of results is not ensured. Any shell for which can be switched to a small strain formulation by RADIOSS Engine option /DT/SHELL/CST, except if Ismstr = 4.

3.	If Ithick or Iplas are activated, the small strain option is automatically deactivated in the corresponding type of element.

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

5.	If the small strain option is set to 1 or 3, the strains and stresses which are given in material laws are engineering strains and stresses; otherwise they are true strains and stresses.

6.	Iplas = 1 is available for Material Law 27.

7.	Flag Istrain is automatically set to 1 for Material Law 25 and Law 27.

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

9.	For hourglass type 3, hourglass maximum values may be larger, default values are 0.1 for hm and hr.

10.	The default value of dm is 5% for Law 25 and Law 27.

11.	The same default value of dm is used for QBAT shells, except for:

•	the default value of dm for QEPH is 1.5% for Material Laws 2, 19, 27, 32, 36 and 43

For further information about dm coefficient, refer to the RADIOSS Theory Manual.

12.	Shell numerical damping dn is only used for Ishell = 12 and 24:

•	for Ishell = 24, dn is used for hourglass stress calculation;

•	for QBAT, dn is used for all stress terms, except transverse shear;

•	for DKT18, dn is only used for membrane.

13.	The default value of dn is:

•	1.5% for Ishell = 24

•	0.1% for QBAT

•	0.01% for DKT18

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

15.	It is recommended to use Iplas = 1, if Ithick= 1.

prop_type9_v10

16.	Input components of the global vector are defined in Line 6 with VX, VY and VZ.

17.	Projection of vector on shell element plane becomes vector

18.	Direction 1 of local coordinate system of orthotropy is defined with vector and angle (angle in degree).

19.	Input as many formats as number of layers (one format per layer, Line 7).

20.	is the angle between direction 1 of orthotropy and projection of vector on the shell for layer i.

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

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.

22.

Each layer as well as the corresponding part must use the same material law, but may have different material properties, hence material IDs. RADIOSS checks for this condition and errors out if it is not met. Global material properties (membrane stiffness, bending stiffness, mass, and inertia) are calculated based on the material properties and layup (thicknesses, …). They are used for stability, mass and interface stiffness. A material is still required at part definition level, but is only used for pre- and post- (visualization “by material”) and its physical characteristics are ignored. The previous formulation where stiffness and mass were calculated from the material associated to the part is still used, if the version number of the input file is V13 or earlier.

23.	Idrill is available for QEPH, QBAT (Ishell = 12 and 24), and standard triangle (C0) shell elements (Ish3n = 1 and 2).

24.	Drilling DOF stiffness is recommended for implicit solutions especially for Riks method and bending dominated problems.

/PROP/TYPE11 (SH_SANDW)

/PROP/TYPE11 (SH_SANDW)

/PROP/TYPE11 (SH_SANDW)

See Also: