/PROP/TYPE10 (SH_COMP)

Block Format Keyword This property set is used to define the composite shell property set. It is possible to define composite with several layers and each layer with individual orthotropic direction.

Format

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/PROP/TYPE10/prop_ID/unit_ID or /PROP/SH_COMP/prop_ID/unit_ID
prop_title
Ishell Ismstr Ish3n Idril     P_thickfail    
hm hf hr dm dn
N   Thick Ashear   Ithick Iplas  
VX VY VZ skew_ID     IP
ϕ1 ϕ2 ϕ3 ϕ4 ϕ5
... ϕN

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 Shell element formulation flag. 1
= 0
Use value in /DEF_SHELL.
= 1 Default, if /DEF_SHELL is not defined
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.

(Integer)

 
Ismstr Shell small strain formulation flag. 2
= -1
Automatically set the best value according element type and material law.
= 0
Use value in /DEF_SHELL.
= 1
Small strain from time = 0 (formulation compatible with all other formulation flags).
= 2 Default, if /DEF_SHELL is not defined
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).

(Integer)

 
Ish3n 3 node shell element formulation flag.
= 0
Use value in /DEF_SHELL.
= 1
Standard triangle (C0).
= 2 Default, if /DEF_SHELL is not defined
Standard triangle (C0) with modification for large rotation.
= 30
DKT18.
= 31
DKT_S3, which based on DTK12 of BATOZ (refer to Element Library in the Theory Manual).

(Integer)

 
Idril Drilling degree of freedom stiffness flag. 7
= 0
Use value in /DEF_SHELL.
= 1
Yes.
2 Default, if /DEF_SHELL is not defined.
No.

(Integer)

 
P_thickfail Percentage of layers that must fail before the element is deleted. 9 10

0.0 P _ t h i c k f a i l 1.0 (Real)

 
hm Shell membrane hourglass coefficient.

Default = 0.01 (Real)

Default = 0.1 for hourglass type 3 (Ishell =3)

 
hf Shell out-of-plane hourglass.

Default = 0.01 (Real)

 
hr Shell rotation hourglass coefficient.

Default = 0.01 (Real)

Default = 0.1 for hourglass type 3 (Ishell =3)

 
dm Shell Membrane Damping.

Only active for Material Laws 25, 32 and 36.

Default =1.5% for Ishell =24 (QEPH)+LAW25, 32 and 36

Default =5.0% for Ishell =1,2,3,4,12+LAW25

Default =0.0% for Ishell =1,2,3,4,12+LAW32 and 36

(Real)

 
dn Shell numerical damping. 4

Only used for Ishell =12 and 24.

Default =1.5% for Ishell =24 (QEPH)

Default =0.1% for Ishell =12 (QBAT)

Default =0.01% for Ishell =30 (DKT18)

(Real)

 
N Number of layers.

Layer thickness is equal to Thick/N with 0 ≤ N ≤ 100.

Default =1 (Integer)

 
Thick Shell thickness.

(Real)

[ m ]
Ashear Shear factor.

Default is Reissner value: 5/6 (Real)

 
Ithick Shell resultant stresses calculation flag. 5
= -1
Automatically set the best value according element type and material law.
= 0
Use value in /DEF_SHELL.
= 1
Thickness change is taken into account.
= 2 Default, if /DEF_SHELL is not defined
Thickness is constant.

(Integer)

 
Iplas Shell plane stress plasticity flag. 6
It is available for Material Laws 2, 22, 32, 36 and 43.
= -1
Automatically set the best value according element type and material law.
= 0
Use value in /DEF_SHELL.
= 1
Iterative projection with three Newton iterations.
= 2 Default, if /DEF_SHELL is not defined
Radial return.

(Integer)

 
VX X component for reference vector. 8

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

Default = 0 (Integer)

 
IP Reference direction in shell plane. 8
= 0 (Default)
Use 1st direction of skew_ID or vector V (if skew_ID is not defined) projected on the shell element.
= 20
Defined from element connectivity (N1,N2) of the shell element.
= 22
Defined from 1st direction of skew_ID projected on the shell element and angle phi. (Vector V is ignored).
= 23

Defined from reference vector V projected on the shell element and angle phi (skew_ID is ignored).

(Integer)

 
ϕ 1 Angle for layer 1.

(Real)

[ deg ]
ϕ 2 Angle for layer 2.

(Real)

[ deg ]
ϕ 3 Angle for layer 3.

(Real)

[ deg ]
ϕ 4 Angle for layer 4.

(Real)

[ deg ]
ϕ 5 Angle for layer 5.

(Real)

[ deg ]
ϕ N Angle for layer N with 0 ≤ N ≤ 100 (5 angles per Line).

(Real)

[ deg ]

Example


prop_type10_example
Figure 1. 6 Layers with Different Material Direction (Fiber Direction)
#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/SH_COMP/2/2
SH_COMP example
#   Ishell    Ismstr     Ish3n    Idrill                             Pthick_fail
        12         0         0         0                                       0
#                 hm                  hf                  hr                  dm                  dn
                   0                   0                   0                  .1                  .1
#        N                         Thick              Ashear              Ithick     Iplas
         6                           1.8                   0                   1         1
#                 Vx                  Vy                  Vz   skew_ID                            Ip
                   1                   0                   1         0                             0
#              Phi_1               Phi_2               Phi_3               Phi_4               Phi_5
                 -60                 -30                   0                  30                  60
                  90
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

Comments

  1. Ishell, Ish3n – 4-node and 3-node shell 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.
    • Ish3n =30 (DKT18): BATOZ DKT18 thin shell with three Hammer integration points.
    • Flag Ishell =2 is incompatible with one integration point for shell element.
  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 Δ t< Δ t min can be switched to a small strain formulation by Radioss Engine option /DT/SHELL/CST, except if Ismstr = 4.
    • If Ismstr = 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.
  3. hm, hf, and hr - Hourglass coefficients
    • hm, hf, and hr are only used for Q4 shells (Ishell =1,2,3,4). They must have a value between 0 and 0.05.
    • For Ishell =3, default values of hm and hr are 0.1 with larger values possible.
  4. dn - Shell numerical damping coefficient
    • dn is only used for Ishell =12, 24 and Ish3n = 30:
      • for Ishell =24 (QEPH) dn is used for hourglass stress calculation.
      • for Ishell =12 (QBAT) dn is used for all stress terms, except transvers shear
      • for Ish3n =30 (DKT18) dn is only used for membrane
  5. Ithick - Shell resultant stresses calculation flag
    • Flag Ithick is automatically set to 1 for /MAT/LAW32 (HILL).
    • If Ithick =1, the small strain option is automatically deactivated in the corresponding type of element.
  6. 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.
  7. Idril
    • Drilling DOF stiffness is recommended for implicit solutions especially for Riks method and bending dominated problems.
    • Idril is available for QEPH, QBAT (Ishell =12 and 24), and standard triangle (C0) shell elements (Ish3n = 1 and 2).
  8. Anisotropy direction definition.
    The reference vector V MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOvaaaa@36B5@ is defined as following according flag IP:
    • If IP=0 and skew_ID = 0, the reference vector V MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOvaaaa@36B5@ is defined with VX, VY and VZ.
    • If IP=0 and skew_ID ≠ 0, the reference vector V MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOvaaaa@36B5@ is the first direction (local X) of the local coordinate system skew_ID.
    • If IP = 20, the reference vector V MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOvaaaa@36B5@ is defined with the node N1 and N2 of the shell elements.
    • If IP = 22, the reference vector V MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOvaaaa@36B5@ is the first direction (local X) of the local coordinate system skew_ID. Vector components VX, VY and VZ are ignored.
    • If IP = 23, the reference vector V MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOvaaaa@36B5@ is defined with VX, VY and VZ. Local coordinate system skew_ID is ignored.
    The reference vector V MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOvaaaa@36B5@ is projected on the shell element plane and becomes the vector V . Then for each layer, the 1st material direction (m1) is vector V turned ϕ i degrees (turns positive direction around shell normal n ).

    prop_type9_v10
    Figure 2.
    The hierarchy order to define the reference vector V is:
    • initial state card (/INISHE/ORTHO)
    • shell property

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

    The 2nd material direction m2 is derived from direction m1 rotated 90 degrees (orthotropic).

    Layers are overlaying one by one from bottom to top.

    Ipos_0
    Figure 3.
  9. P_thickfail parameter is not compatible with failure defined within the material law itself, such as plastic failure strain in LAW36.
  10. Element deletion rules used with /FAIL models:
    • Under-integrated elements (Belytchko, QEPH, DKT18):
      • f more than one failure model is applied to the shell material, or P_thickfail=0 (blank), the value of P_thickfail is calculated individually from each failure model settings.
        Example:
        Ifail_sh = 1
        One integration point failure sufficient to delete element.
        P_thickfail = 1.0 e-6
        Ifail_sh = 2
        All integration points failure is necessary to delete the element.
        P_thickfail = 1.0
      • If only one failure model is applied to the material, the P_thickfail value from the property is taken by default and overwrites local Ifail_sh settings from the failure model.
    • 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.