/PROP/TYPE17 (STACK)

Block Format Keyword

/PROP/TYPE17 - Stacking information for ply-based Sandwich Shell Property Set

Description

This property set is used to define the sandwich shell property set using the stack and ply approach.

Two possible input:

•	All elements may have different number of plies by associating a group of elements to the ply in a way similar to OptiStruct.

•	The layout is described by a sequence of substacks (“SUB”), which are comprised of several plies. The interplay between substacks is defined with the keyword INT.

The input can be either “substack” or by ply, but not both.

Format

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

By ply

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

Or by substack

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
Sub	Nsub	Sub-plyn
Substack Name
Pply_IDi			Zi
Minterply
Sub	Nsub	Sub-plyn
Substack Name
Pply_IDj			Zj
Minterply
INT	Pply_IDt	Pply_IDb

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 (Comment 4) (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 (Integer) = 0: no = 1: yes
Plyxfem	Activates plyxfem formulation (beta version) (Integer) = 0: no = 2: yes
Z0	Distance from the shell element reference plane to the bottom surface of the shell (Real or blank). It should be considered, if Ipos =2 Default = 0.0 (Real)
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)
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 4) (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 34) Default = 0 (Integer) = 0: layer positions are automatically calculated with regard to layer thicknesses. 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. Global thickness is not checked in this case, since it needs to not be equal to sum of layer thicknesses. =2: the shell element reference plane is at Z0 from the bottom surface of the shell. =3: the top surface of the shell is considered as element reference plane. =4: the bottom surface of the shell is considered as element reference plane.
Pply_IDi Pply_IDj	Ply property identifier for layer i (j) (Integer) i = 1, 2, 3 …N j = 1, 2, 3 …N
Minterply	Interply material identifier (needed only if plyxfem formulation is activated). This material is ignored, if the plyxfem formulation is not activated. (Integer)
	Angle for layer i (j) (Comment 23) (Real)
Zi Zj	Z position of layer i (j) (Zi or Zj defines the position of the middle of the layer) Default = 0.0 (Real)
Sub	=SUB: Indicates the beginning of a substack definition (must be left justified) (Integer)
Nsub	Substack identification number (Integer)
Sub-plyn	Number of plies in substack (Integer)
Substack Name	Name of the substack (Maximum 100 characters)
Int	=INT: Indicates the connection between 2 substacks (Optional, if used must be left justified) (Integer)
Pply_IDt	(Optional) Ply identification number of the top ply of the bottom substack (Integer)
Pply_IDb	(Optional) Ply identification number of the bottom ply of the top substack (Integer)

There are three plies in this example. Use /DRAPE to change the ply orientation so that each ply has a different material direction (fiber direction).And different material (which is defined in each ply with /PROP/TYPE19) is also applied.

prop_type17_ply_example

#RADIOSS STARTER

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

/PROP/TYPE17/2/2

composite combine by ply

# Ishell Ismstr Ish3n Idrill Plyxfem Z0

12 0 0 1 0 0

# hm hf hr dm dn

0 0 0 .1 .1

# Istrain Thick Ashear Ithick Iplas

0 0 0 1 1

# VX VY VZ skew_ID Iorth Ipos

1 0 1 0 0 0

# Pply_IDi PHIi Zi

11 0 0

#Minterply

# Pply_IDi PHIi Zi

12 90 0

#Minterply

# Pply_IDi PHIi Zi

13 0 0

#Minterply

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

/PROP/TYPE19/11/2

PROP number 11

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 0 0 0 0

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

/PROP/TYPE19/12/2

PROP number 12

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

2 .6 0 0 0 0 0

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

/PROP/TYPE19/13/2

PROP number 12

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 -45 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----|

Substack definition for the ply layout.

prop_type17_example

#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/TYPE17/2/2

composite combine by substack

# Ishell Ismstr Ish3n Idrill Plyxfem Z0

12 0 0 1 0 0

# hm hf hr dm dn

0 0 0 .1 .1

# Istrain Thick Ashear Ithick Iplas

0 0 0 1 1

# VX VY VZ skew_ID Iorth Ipos

0 1 0 0 0 0

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

# Sub Nsub Sub-plyn

SUB 1 4

Top substack

# Pply_IDi PHIi Zi

11 0 0

#Minterply

12 90 0

#Minterply

13 0 0

#Minterply

14 90 0

#Minterply

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

# Sub Nsub Sub-plyn

SUB 2 4

Left substack

# Pply_IDi PHIi Zi

21 90 0

#Minterply

22 0 0

#Minterply

23 90 0

#Minterply

24 0 0

#Minterply

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

# Sub Nsub Sub-plyn

SUB 3 4

Right substack

# Pply_IDi PHIi Zi

31 90 0

#Minterply

32 0 0

#Minterply

33 90 0

#Minterply

34 0 0

#Minterply

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

# Sub Nsub Sub-plyn

SUB 4 3

Middle substack

# Pply_IDi PHIi Zi

41 90 0

#Minterply

42 0 0

#Minterply

43 90 0

#Minterply

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

# Int Pply_IDt Pply_IDb

INT 14 21

INT 14 31

INT 43 21

INT 34 41

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

/PROP/TYPE19/11/2

PROP number 11

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 45 0 0 0

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

/PROP/TYPE19/12/2

PROP number 12

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 45 0 0 0

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

/PROP/TYPE19/13/2

PROP number 13

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 45 0 0 0

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

/PROP/TYPE19/14/2

PROP number 14

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 45 0 0 0

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

/PROP/TYPE19/21/2

PROP number 21

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 46 0 0 0

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

/PROP/TYPE19/22/2

PROP number 22

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 46 0 0 0

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

/PROP/TYPE19/23/2

PROP number 23

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 46 0 0 0

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

/PROP/TYPE19/24/2

PROP number 24

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 46 0 0 0

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

/PROP/TYPE19/31/2

PROP number 31

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 47 0 0 0

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

/PROP/TYPE19/32/2

PROP number 32

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 47 0 0 0

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

/PROP/TYPE19/33/2

PROP number 33

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 47 0 0 0

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

/PROP/TYPE19/34/2

PROP number 34

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 47 0 0 0

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

/PROP/TYPE19/41/2

PROP number 41

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 44 0 0 0

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

/PROP/TYPE19/42/2

PROP number 42

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 44 0 0 0

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

/PROP/TYPE19/43/2

PROP number 43

# mat_ID_i t delta_phi grsh4n_ID grsh3n_ID Npt_ply A_i

1 .5 45 44 0 0 0

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

#enddata

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

1.	The stack property is used in combination with /PROP/PLY (/PROP/TYPE19) to create composites properties through the ply-based definition.

2.	Only compatible with Material Laws 25, 27, 36, 60 and user laws.

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

QEPH: formulation with hourglass physical stabilization for general use.

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.

The input thickness may be equal to zero as the true shell thickness is always calculated from the plies thicknesses and positions. If the input thickness is different from 0.0, a warning will be displayed if the calculated thickness from the plies is different from the input thickness, but the element thickness will be calculated from the plies thicknesses and positions as above.

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

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

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

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

9.	Flag Iplas is available for Material Law 27.

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

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

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

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

14.	The same default value of dm is used in 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.

15.	Shell numerical damping dn is only used for Ishell = 12, 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.

16.	The default value of dn is:

•	1.5% for Ishell = 24

•	0.1% for QBAT

•	0.01% for DKT18

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

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

prop_type9_v10

19.	Input components of global vector are defined in Line 6.

20.	Projection of vector on shell element plane becomes vector

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

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

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

24.	Layer positions must be defined if the Ipos flag is active. The Zi values are real layer positions in the local Z axis (negative and positive value are allowed and Zi defines the position of the middle of the layer).

25.	Material law type must be the same for each layer.

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

27.	The material law type input in the element definition must be identical to the material law type used in Line 7.

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

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

30.	If the substack definition is used, each substack definition begins with the keyword “SUB” and the substack identification number.

All the plies until the next keyword “SUB” or the keyword “INT” is encountered belong to the same substack. The substack definition can only be used if N (the input number of layer) is equal to 0.

31.	All substacks must be defined before defining the INT connections.

32.	If an element contains plies belonging to several susbtacks, the connection between the substacks is defined with an INT connection: the ply PPLY_IDt of the first substack is connected to the ply PPLY_IDb of the second substack.

33.	The old input is still supported, but considered as obsolete.

34.	Ipos = 1 is not recommended unless the number of ply is constant for all elements.

The mid-surface of the shell is the surface defined by the nodes connected to the element. For an exact definition of this surface, refer to the Theory Manual.

35.

Ipos = 2:

prop_type17_ipos2

36.

Ipos = 3:

prop_type17_ipos3

37.

Ipos = 4:

prop_type17_ipos4

38.	The plyxfem formulation is activated using the Plyxfem flag. An interply material must be defined between each ply. The material must use LAW1 and be compatible with Ladeveze failure criteria for delamination. This option is in the beta stage.