/PROP/TYPE1 (SHELL)
Block Format Keyword Describes the shell property set which is used for 3-node or 4-node shell element. Belytschko, QBAT or QEPH shell formulation are available.
Format
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
/PROP/TYPE1/prop_ID/unit_ID or /PROP/SHELL/prop_ID/unit_ID | |||||||||
prop_title | |||||||||
Ishell | Ismstr | Ish3n | Idril | Ipinch | P_thickfail | ||||
hm | hf | hr | dm | dn | |||||
N | Thick | Ashear | Ithick | Iplas |
Definition
Field | Contents | SI Unit Example |
---|---|---|
prop_ID | Property
identifier. (Integer, maximum 10 digits) |
|
unit_ID | Unit Identifier. (Integer, maximum 10 digits) |
|
prop_title | Property
title. (Character, maximum 100 characters) |
|
Ishell | 4 node shell element
formulation flag.
(Integer) |
|
Ismstr | Shell small strain
formulation flag. 4
(Integer) |
|
Ish3n | 3 node shell element
formulation flag.
(Integer) |
|
Idril | Drilling degree of freedom
stiffness flag. 3
(Integer) |
|
Ipinch | Pinching degree of freedom
flag. 12
(Integer) |
|
P_thickfail | Ratio of through thickness
integration points that must fail before the element is deleted.
11
(Real) |
|
hm | Shell membrane hourglass
coefficient. 5 Only used when Ishell=1, 2, 3, 4. Default = 0.01 for Ishell =1, 2, 4. Default = 0.1 for Ishell = 3. (Real) |
|
hf | Shell out-of-plane
hourglass. 5 Only used when Ishell=1, 2, 3, 4. Default = 0.01 (Real) |
|
hr | Shell rotation hourglass
coefficient. 5 Only used when Ishell=1, 2, 3, 4. Default = 0.01 for Ishell =1, 2, 4. Default = 0.1 for Ishell = 3. (Real) |
|
dm | Shell Membrane Damping For default values, refer to the table in Comment 7. (Real) |
|
dn | Shell numerical damping.
6 It is only used for Ishell =12, 24 and Ish3n = 30 Default =0.015 for Ishell =24 (QEPH) Default =0.001 for Ishell =12 (QBAT) Default =0.0001 for Ish3n =30 (DKT18) (Real) |
|
N | Number of integration
points through the thickness with 0 < N <
10. 10
(Integer) |
|
Thick | Shell
thickness. (Real) |
|
Ashear | Shear factor. Default is Reissner value: 5/6 (Real) |
|
Ithick | Shell resultant stresses
calculation flag. 8
(Integer) |
|
Iplas | Shell plane stress
plasticity flag. 9 It is available for material laws 2, 22, 27 and 36.
(Integer) |
Comments
- Ishell - Shell small strain
formulation flag
- Ishell=1,2,3,4 (Q4): original 4 node Radioss shell with hourglass perturbation stabilization
- Ishell=24 (QEPH): formulation with hourglass physical stabilization for general use
- Ishell=12 (QBAT): modified BATOZ Q4γ24 shell with four Gauss integration points and reduced integration for in-plane shear. No hourglass control is needed for this shell
- Ishell=2 is incompatible with one integration point (through thickness) for shell element
- Ish3n - 3 node shell element
formulation flag
- Ish3n=30 (DKT18): BATOZ DKT18 thin shell with three Hammer integration points
- Idril - Drilling degree of
freedom stiffness flag
- Drilling DOF stiffness is recommended for implicit solutions especially for Riks method and bending dominated problems
- Idril is available for QEPH, QBAT (Ishell =12, 24), and standard triangle (C0) shell elements (Ish3n = 1, 2)
- Ismstr - Small strain formulation
- 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 the 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.
- 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.
- Ismstr =10 is only compatible with these material laws 42, 62, 69, and 82 that use total strain formulation. Generally, the Left Cauchy-Green strain is used. For User laws, the deformation gradient tensor and the right stretch tensor may be used.
- hm, hf, and hr - Hourglass coefficients
- Are only used for Q4 shells (Ishell=1,2,3,4). They must have a value between 0 and 0.05.
- For Ishell=3, the default values of hm and hr are 0.1 with larger values possible.
- dn - Shell
numerical damping coefficient
- dn is only used for Ishell =12, 24 and Ish3n = 30
- Ishell =24 (QEPH) dn is used for hourglass stress calculation
- Ishell =12 (QBAT) dn is used for all stress terms, except transvers shear
- Ish3n =30 (DKT18) dn is only used for membrane
- dn is only used for Ishell =12, 24 and Ish3n = 30
- Default value for
numerical damping coefficient
dm.
dm Ishell= 1, 2, 3, 4 (Q4) 0.05 for LAW27 0.0 for LAW32, 36
Ishell=12 (QBAT) 0.0 for LAW2, 27, 32, 36 Ishell=24 (QEPH) 0.015 for all other material Ish3n= 1, 2, 31 0.0 Ish3n=30 (DKT18) 0.0 Any Ishell, Ish3n 0.05 for LAW65 - Ithick - Shell resultant stresses
calculation flag
- Ithick is automatically set to 1 for material LAW32
- Ithick =1, the small strain option is automatically deactivated in the corresponding type of element
- Iplas - Shell plane stress
plasticity flag
- It is recommended to use Iplas =1, if Ithick =1
- If Iplas =1, the small strain option is automatically deactivated in the corresponding type of element
- Global integration,
N=0
- Global integration is only compatible with material LAW1.
- Failure models are not available with global integration for shells.
- State files (/STATE/DT which creates /INISHE and /INISH3N) are not compatible with global integration.
- Automatically set to N=3 for the other material laws.
- Element deletion rules
used with /FAIL models.
- Under-integrated elements (Belytschko, QEPH, DKT18):
- If only one /FAIL model with its own P_thickfail (/FAIL/BIQUAD, /FAIL/TAB1) is applied to the shell material, and P_thickfail on the /FAIL card is set > 0, the value of P_thickfail is used regardless of P_thickfail set on the /PROP/SHELL card. If the value of P_thickfail on the /FAIL card is set at 0, the P_thickfail from the property card is used (no default value, so should be set to > 0 to avoid rupture after single integration point failure if that is not desired).
- If only one /FAIL model without its own P_thickfail (/FAIL/JOHNSON) is applied to the shell material, and the P_thickfail value from the /PROP/SHELL property is set to > 0, the property P_thickfail is taken by default and overwrites Ifail_sh settings from the failure model. If P_thickfail on the property is set = 0 (or blank), then Ifail_sh from the failure model is used (for legacy compatibility).
- If more than one /FAIL failure model is
applied to the shell material, the value of
P_thickfail
is calculated individually from each of the failure model
P_thickfail
or Ifail_sh
settings and
P_thickfail
on the property card is ignored.For Example, per failure model:
- If Ifail_sh = 1 (one integration point failure sufficient to delete element) = P_thickfail= 1.0 e-6
- If Ifail_sh = 2 (all integration points failure is necessary to delete the element) = P_thickfail = 1.0
- If
P_thickfail
(/FAIL) = 0, =
P_thickfail
= 0.0
In this case on a /FAIL card with P_thickfail, failure will occur when a single integration point fails, due to this /FAIL (set to > 0 if this is not desired).
- Fully-integrated shells (Batoz, DKT_S3)
The rule described for under-integrated shells applies to each Gauss point separately. P_thickfail ratio is checked for all integration points in thickness for each in-plane Gauss point. The element is deleted only when all Gauss points reach P_thickfail ratio criteria.
- Under-integrated elements (Belytschko, QEPH, DKT18):
- Ipinch is used to enable or disable set pinching degree of freedom. Currently, only supported for QBAT shells (Ishell=12).