/TH/SHEL

Block Format Keyword

/TH/SHEL - Time History Output for 4-Node Shell Elements

Description

Describes the time history output for 4-node shell elements. Output parameters are components of stress and strain tensor expressed in different coordinate systems, plastic strains, strain rate, internal energy, element thickness, element deletion flag, etc.

Format

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
/TH/SHEL/thgroup_ID
thgroup_name
var_ID1	var_ID2	var_ID3	var_ID4	var_ID5	var_ID6	var_ID7	var_ID8	var_ID9	var_ID10
elem_ID	skew_ID	elem_name

Field	Contents
thgroup_ID	TH group identifier (Integer, maximum 10 digits)
thgroup_name	TH group name (Character, maximum 100 characters)
var_ID1, ...n	Variables saved for TH (see table below) (Character, maximum 8 characters)
elem_ID	Element identifier (Integer)
skew_ID	Skew coordinate system identifier for stress output. Default element coordinate system (Integer)
elem_name	Name of the element to appear in the time history (Integer, maximum 80 characters)

Keyword

Object saved

Variables

SHEL

4-node shells

F1, F2, F12, Q1, Q2, M1, M2, M12

IEM, IEB, OFF, THIC, EMIN, EMAX, EPSD, E1, E2, E12, SH1, SH2, K1, K2, K12, USRi (i=1,60), MF1, MF2, MF12, MQ1, MQ2

USRII_JJ (II=1,60; JJ=1,99), USRII_JKK (II=1,60; J=1,4; KK=1,99)

SX_JJ, SY_JJ, SXY_JJ, SYZ_JJ, SZX_JJ (JJ=1,99)

Keyword

Variable Group

Saved TH Variables

SHEL

DEF

STRESS

STRAIN

PLAS

FAILURE

WPLAY01_10

WPLAY11_20

WPLAY21_30

WPLAY31_40

WPLAY41_50

WPLAY51_60

WPLAY61_70

WPLAY71_80

WPLAY81_90

WPLAY91_99

F1, F2, F12, M1, M2, M12, IEM, IEB, OFF, EMIN, EMAX

F1, F2, F12, Q1, Q2, M1, M2, M12

E1, E2, E12, SH1, SH2, K1, K2, K12

EMIN, EMAX

NFAIL, PFAIL, FAIL_D1, FAIL_D2, FAIL_EN

WPLAY01, …, WPLAY10

WPLAY11, …, WPLAY20

WPLAY21, …, WPLAY30

WPLAY31, …, WPLAY40

WPLAY41, …, WPLAY50

WPLAY51, …, WPLAY60

WPLAY61, …, WPLAY70

WPLAY71, …, WPLAY80

WPLAY81, …, WPLAY90

WPLAY91, …, WPLAY99

Output for Shell 4-Nodes	Available for Material Law
OFF: element flag for deactivation (see Comment 7)	1, 2, 3, 4, 5, 6, 10, 11, 15, 16, 18, 19, 20, 21, 22, 23, 24, 25, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 69, 72, 76, 78, 79, 80, 82, 87, 92
F1: stress in direction 1	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 52, 66, 72, 76, 78, 80, 82, 87, 92
F2: stress in direction 2	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 52, 66, 72, 76, 78, 80, 82, 87, 92
F12: stress in direction 12	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 52, 66, 72, 76, 78, 80, 82, 87, 92
Q1: mean stress in direction 13	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 52, 66, 72, 76, 78, 80, 82, 87, 92
Q2: mean stress in direction 23	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 52, 66, 72, 76, 78, 80, 82, 87, 92
M1: moment per unit length per unit thickness square in direction 1	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 72, 76, 78, 80, 82, 87, 92
M2: moment per unit length per unit thickness square in direction 2	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 72, 76, 78, 80, 82, 87, 92
M12: moment per unit length per unit thickness square in direction 12	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 72, 76, 78, 80, 82, 87, 92
IEM: internal membrane energy per volume unit	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 72, 76, 78, 80, 82, 87, 92
IEB: internal bending energy per volume unit	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 72, 76, 78, 80, 82, 87, 92
THIC: thickness	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 72, 76, 78, 80, 82, 87, 92
EMIN: minimum equivalent plastic strain over integration point	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 48, 66, 72, 76, 78, 80, 82, 87, 92
EMAX: maximum equivalent plastic strain over integration point	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 72, 76, 78, 80, 87
EPSD: equivalent strain rate	2, 15, 25, 27, 36, 44, 48, 76
E1: membrane strain in direction 1	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 48, 66, 72, 76, 78, 80, 82, 87, 92
E2: membrane strain in direction 2	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 72, 76, 78, 80, 82, 87, 92
E12: membrane strain in direction 12	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 72, 76, 78, 80, 82, 87, 92
SH1: shear strain in direction 1	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 72, 76, 78, 80, 82, 87, 92
SH2: shear strain in direction 2	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 72, 76, 78, 80, 82, 87, 92
K1: curvature in direction 1	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 72, 76, 78, 80, 82, 87, 92
K2: curvature in direction 2	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 72, 76, 78, 80, 82, 87, 92
K12: curvature in direction 12	0, 1, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 72, 76, 78, 80, 82, 87, 92
USR1 ..., USR5	29, 30, 31, 52
USR1 ..., USR5, ..., USR60	29, 30, 31
USRII_JJ, USII_JKK	29, 30, 31, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 66, 72, 76
SX_JJ, SY_JJ, SXY_JJ, SYZ_JJ, SZX_JJ (JJ=1,99)	0, 2, 3, 19, 22, 25, 27, 28, 29, 32, 36, 37, 38, 40, 42, 43, 44, 46, 48, 52, 66, 72, 76, 78, 80, 82, 87, 92
NFAIL	25
PFAIL	25
FAIL_D1	25
FAIL_D2	25
FAIL_EN	25
WPLAYJJ (JJ=1,99)	25

1.	Available names are listed in the two tables above.

In the first table, TH-variables are listed. If a TH-variable is requested, this variable is written.

In the second table, other variable groups are listed. If a variable group is requested, all the associated TH-variables are written at once.

2.	The value for EPSD (equivalent strain rate) is only computed in case of the strain rate filtering is asked for the material law. It is available for Laws 2, 15, 25, 27, 36, 44 and 48.

3.	Shear strains are .

4.	For Material Law 25, the output EMIN and EMAX corresponds to minimum plastic work and respectively maximum plastic work.

5.	The value of EPSD will be computed only in case of strain rate filtering is asked for in the material law.

6.	For Material Law 25:

NFAIL is the total number of failed layers

PFAIL is the percentage of failed layers

FAIL_D1 is the number of layers which reached the failure level in direction 1

FAIL_D2 is the number of layers which reached the failure level in direction 2

FAIL_EN is the number of layers which reached the failure level in plastic work

WPLAYJJ is the plastic work for layer JJ

OFF:

0.0 deleted element
1.0 active element
2.0 (shells / solids): active element using small strain
Negative value -1/-2 when the element is on stand-by due to rigid body

Some intermediate values can be found during few cycles before element deletion (progressive reduction of off during a few cycles applies as a multiplicative factor to the stress, so that the stress is progressively vanished and element is deleted after ~ 10 cycles). This process applies to some element types, but not all.

Skew system is used to determine coordinate system to output generalized (averaged through the thickness) stress tensor components and moments F1, F2, F12, Q1, Q2, M1, M2, and M12. By default the skew system is coincident with elemental coordinate system. When the skew system is specified, then X- axis of the system is projected to the element plane to obtain direction 1 of the local coordinate system. Direction 2 is obtained by rotation of the first axis by 90 degrees. The angle between direction 1 and X-direction of the elemental coordinate system is fixed and does not change during simulation. Components F1, F2, F12, Q1, Q2, M1, M2, and M12 are projected to this local coordinate system.

9.	SX_JJ, SY_JJ, SXY_JJ, SYZ_JJ, SZX_JJ (JJ=1,99) are the stress components of the stress matrix expressed in the local skew per layer through element thickness.

The local skew is elemental coordinate system for any isotropic material law as well as material laws 19, 27, and 32. Local skew is the orthotropic material coordinate system in case of orthotropic material laws 15, 25, and all orthotropic material laws higher than 28. LAW1 has only global integration or membrane formulations (with NP=1 in /PROP/SHELL). Therefore stress tensor is calculated and only global force output (F1, F2, F12, M1, M2 and M12) is available.