/PROP/TYPE16 (SH_FABR)

Block Format Keyword This property set is used to define the anisotropic layered shell property set. This property is currently only compatible with elastic anisotropic fabric (/MAT/LAW58 (FABR_A)).

Format

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/PROP/TYPE16/prop_ID/unit_ID or /PROP/SH_FABR/prop_ID/unit_ID
prop_title
Ishell Ismstr Ish3n       P_thickfail    
hm hf hr dm    
N   Thick Ashear   Ithick    
VX VY VZ skew_ID Ipos   IP
ϕi α 1 ti Zi mat_IDi  
Each layer (integration point) per line:
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
ϕ i α 1 ti Zi mat_IDi  

Definitions

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.
= 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.
= -1
Automatically set the best value according to element type and material law.
= 4 (Default)
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

(Integer)

 
P_thickfail Element suppression criterion. It is ratio of the sum of failed integration points thickness over total shell thickness necessary to delete the element. 7

This option used for /FAIL/FABRIC.

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

(Real)

 
hm Shell membrane hourglass coefficient.

Only used when Ishell=1, 2, 3, 4.

Default = 0.01

Default = 0.1 for Ishell =3

(Real)

 
hf Shell out-of-plane hourglass.

Only used when Ishell=1, 2, 3, 4.

Default = 0.01 (Real)

 
hr Shell rotation hourglass coefficient.

Only used when Ishell=1, 2, 3, 4.

Default = 0.01

Default = 0.1 for Ishell =3

(Real)

 
dm Shell membrane damping.

It is used for material LAW58.

Default =0.075 for Ishell ≠12

Default =0.2 for Ishell =12

Default = 0.0 for all 3 node shells (Real)

 
N Number of layers, with 1 N 100 . 6

Default = 1 (Integer)

 
Thick Shell thickness.

(Real)

[ m ]
Ashear Shear factor.

Default is Reissner value: 5/6 (Real)

 
Ithick Shell resultant stresses calculation flag.
= -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)

 
VX X component for reference vector. 3

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

Default = 0 (Integer)

 
Ipos Layer positioning flag for reference vector. 4
= 0 (Default)
Layer positions Zi are automatically calculated with regard to layer thicknesses.
= 1
All layer positions Zi must be user-defined.

Default = 0 (Integer)

 
IP Reference direction in shell plane. 3
= 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)

 
ϕ i Angle ϕ i of 1st local axis for layer i. 3

(Real)

[ deg ]
α i MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySde2aaS baaSqaaiaadMgaaeqaaaaa@388F@ Angle between first and second axis. 3

Default = 90.0 (Real)

[ deg ]
ti Thickness ti of layer i.

(Real)

[ m ]
Zi Z position of layer i (Zi defines the position of the middle of the layer).

Default = 0.0 (Real)

[ m ]
mat_IDi Material identifier for layer i. Same material should be used for all layers.

(Integer)

 

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
                  Mg                  mm                   s
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  2. GEOMETRICAL SETS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/PROP/TYPE16/2/2
prop type 16
#   Ishell    Ismstr     Ish3n                                       Pthick_fail
         0         0         0                                                 0
#                 hm                  hf                  hr                  dm
                   0                   0                   0                   0
#        N                         Thick              Ashear              Ithick
         3                           1.6                   0                   0
#                 Vx                  Vy                  Vz   skew_ID      Ipos                  Ip
                   1                   0                   1         0         0                   0
#                Phi               Alpha                  Ti                  Zi   mat_IDi
                  45                  90                  .5                             1          
                  90                  90                  .6                             2          
                 -45                  90                  .5                             1          
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

Comments

  1. Ishell – 4-node shell formulation flag

    The hourglass formulation is visco-elastic for Q4 shells (Ishell=1,2,3,4).

  2. hm, hf, and hr - Hourglass coefficients
    • hm, hf, and hr are used only for Q4 shells (Ishell=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.
  3. 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 1.
    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 with the angle α i . If α i = 90 ° MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySde2aaS baaSqaaiaadMgaaeqaaOGaeyypa0JaaGyoaiaaicdacqGHWcaSaaa@3D29@ , the layer is orthotropic.

  4. Ipos – Layer position

    Ipos= 0: layer positions are calculated automatically.

    If (1)
    T h i c k i N t i
    • A warning message is displayed
    • Individual layer thickness will be adjusted to new layer thickness t i n e w with:(2)
      T h i c k = i N t i n e w
      Here "Thick" and t i are the shell thickness and layer thickness which specified in input.

      Ipos_0
      Figure 2.
    • Ipos = 1: all layer positions in the element thickness are user defined (with t i and Z i ).
      • Thick” is not checked, as it does not need to be equal to the sum of layer thickness.
      • Multiple layers are allowed to have the same space position.

    For more details, refer to Layer thickness and position calculation. in FAQs.

  5. The material law number given in /PART will be used to define the density and Young’s Modulus which then is used to calculate the time step and interface stiffness.
  6. Each i layer has only 1 integration point through its thickness.
  7. Option P_thickfail is used by /FAIL/FABRIC criterion. This is the only criterion which is compatible with /PROP/TYPE16 and /MAT/LAW58.