/PROP/TYPE18 (INT_BEAM)

Block Format Keyword Describes the integrated beam property set. This beam model is based on Timoshenko theory and takes into account transverse shear strain without warping in torsion.

It can be used for deep beam cases (short beams). Beam section and position of integration points can be either used as predefined or prescribed directly.

Format

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/PROP/TYPE18/prop_ID/unit_ID or /PROP/INT_BEAM/prop_ID/unit_ID
prop_title
Isect Ismstr
dm df
NIP Iref Y0 Z0
If NIP > 0, add NIP cards defining the subsection parameters (each integration point per line)
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Yi Zi Area
If Isect > 0, add following 2 lines
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
NITR L1 L2 L3 L4
L5 L6
Add flag for rotational DOF for the beam nodes.
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
ω DOF

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)

Isect Section type. 5
= 0 (Default)
Integrated beam.
= 1
Predefined rectangular section.
=2
Predefined circular section.
= 3
Predefined rectangular section with Gauss-Lobatto quadrature.
= 4
Predefined circular section with Gauss-Lobatto quadrature.
= 5
Predefined circular section.
= 10
Predefined I-shape section.
= 11
Predefined channel section.
= 12
Predefined L-shape section.
= 13
Predefined T-shape section.
= 14
Predefined tubular box section.
= 15
Predefined Z-shape section.
= 16
Predefined trapezoidal section.
= 17
Predefined circular section.
= 18
Predefined tubular section.
= 19
Predefined I-shape type 2 section.
= 20
Predefined solid box section.
= 21
Predefined cross-shape section.
= 22
Predefined H-shape section.
= 23
Predefined T-shape type 2 section.
= 24
Predefined I-shape type 3 section.
= 25
Predefined channel type 2 section.
= 26
Predefined channel type 3 section.
= 27
Predefined T-shape type 3 section.
= 28
Predefined box shape section.
= 29
Predefined hexagonal section.
= 30
Predefined hat shape section.
= 31
Predefined closed hat shape section.

(Integer)

Ismstr Small strain option flag.
= 0 (Default)
Set to 4.
= 1
Small strain formulation from t = 0.
= 4
Full geometric nonlinearities.

(Integer)

dm Beam membrane damping.

Default = 0.00 (Real)

df Beam flexural damping.

Default = 0.01 (Real)

NIP Number of integration points (subsections).

Only for Isect =0; otherwise, NIP=0.

(Integer)

Iref Subsection center reference flag.

Only for Isect =0.

= 0 (Default)
Subsection center is calculated as a barycenter of the integration points.
= 1
Subsection center is defined by using local coordinates (Y0 and Z0).

(Integer)

Y0 Local Y coordinate of the section center.

Only for Isect =0.

(Real)

[ m ]
Z0 Local Z coordinate of the section center.

Only for Isect =0.

(Real)

[ m ]
Yi Local Y coordinate of the integration point.

(Real)

[ m ]
Zi Local Z coordinate of the integration point.

Only for Isect =0.

(Real)

[ m ]
Area Area of the subsection.

Only for Isect =0.

(Integer)

[ m 2 ]
NITR Option for integration points in predefined section for Isect > 0. 5.

Default = 2 if (1Isect5)

Default = 0 if (Isect10).

(Integer)

L1 First size of the predefined section for Isect > 0. 5

(Real)

[ m ]
L2 Second size of the predefined section for Isect > 0. 5

Default = L1 for Isect = 1 or 3 (Real)

[ m ]
L3 Third size of the predefined section for Isect > 0. [ m ]
L4 Fourth size of the predefined section for Isect > 0. [ m ]
L5 Fifth size of the predefined section for Isect > 0. [ m ]
L6 Sixth size of the predefined section for Isect > 0. [ m ]
ω D O F Rotation DOF code of nodes 1 and 2 (see detail input below).

(6 Booleans)

Detail of Rotation DOF Input Fields for Nodes 1 and 2

(1)-1 (1)-2 (1)-3 (1)-4 (1)-5 (1)-6 (1)-7 (1)-8 (1)-9 (1)-10
ω X 1 ω Y 1 ω Z 1 ω X 2 ω Y 2 ω Z 2

Definition

Field Contents SI Unit Example
ω X 1 = 1 Rotation DOF about X at node 1 is released.

(Boolean)

ω Y 1 = 1 Rotation DOF about Y at node 1 is released.

(Boolean)

ω Z 1 = 1 Rotation DOF about Z at node 1 is released.

(Boolean)

ω X 2 = 1 Rotation DOF about X at node 2 is released.

(Boolean)

ω Y 2 = 1 Rotation DOF about Y at node 2 is released.

(Boolean)

ω Z 2 = 1 Rotation DOF about Z at node 2 is released.

(Boolean)

Example 1 (Integrated Beam)

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/2
unit for prop
                  Mg                  mm                   s
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/PROP/TYPE18/4/2
Integrated beam - bXh=10X10 with 4 integration points (subsections)
#    Isect   Ismstmr
         0         0
#                 dm                  df
                   0                   0
#      NIP      Iref                  Y0                  Z0
         4         1                   0                   0
#                  Y                   Z                Area
                 2.5                 2.5                  25
                 2.5                -2.5                  25
                -2.5                 2.5                  25
                -2.5                -2.5                  25
# OmegaDOF                                                                        
   000 000
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

Example 2 (Integrated Beam)

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/2
unit for prop
                  Mg                  mm                   s
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/PROP/TYPE18/4/2
Integrated beam  - 4 integration points in predefined section bXh=10X10
#    Isect   Ismstmr
         1         0
#                 dm                  df
                   0                   0
#      NIP      Iref                  Y0                  Z0
         0         1                   0                   0
#     NITR                            L1                  L2                  L3                  L4
         2                            10                  10                   0                   0 
#                 L5                  L6
                   0                   0
# OmegaDOF                                                                        
   000 000
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

Comments

  1. Small strain formulation is activated from time t=0, if Ismstr =1. It may be used for a faster preliminary analysis because Δ t is constant, but the accuracy of results is not ensured.
  2. If Ismstr =1, the strains and stresses which are given in material laws are engineering strains and stresses. Time history output returns true strains and stresses.
    Figure 1.

    clip0093
  3. The cross-section of the element is defined using up to 100 integration points (Figure 2). The element properties of the cross-section, that is, area moments of inertia and area, are computed by Radioss as:
    A = A i = ( d y i d z i )
    I Z = A i ( y i 2 + 1 12 d y i 2 )
    I Y = A i ( z i 2 + 1 12 d z i 2 )
  4. It can be used for deep beam cases (short beams). The use of several integration points in the section allows to get an elasto-plastic model in which von Mises criteria is written on each integration point and the section can be partially plastified contrary to the classical beam element (TYPE3). Compatible with material LAW1, LAW2, LAW36 and LAW44. However, as the element has only one integration point in its length, it is not recommended to use a single beam element per line of frame structure in order to take into account the plasticity progress in length, as well as in depth.
    Figure 2. Cross-section Definitions in the Integrated Beam

    cross-section_def
  5. Predefined cross-sections are available (circular or rectangular). Number of integration points in the section is prescribed via NITR depending on Isect and the chosen quadrature:
    • For Isect = 1 and 2. Integration points are distributed uniformly across the section according to the section type and NITR.
      Figure 3.


    • For Isect = 3. The section is rectangular, the distribution of the integration points corresponds to the Gauss-Lobatto quadrature with points on the edge. IP = NITR*NITR. The maximum NITR possible is 9 corresponding to 81 integration points.
      Figure 4.


    • For Isect = 4. The section is circular. The distribution of the integration points radially corresponds to Gauss-Lobatto quadrature with points on edge. Three options: NITR = 1, NITR = 17 and NITR = 25. Here the number of integration points is equal to NITR.
      Figure 5.


    • For Isect = 5. The section is circular. Three options: NITR = 1, NITR = 9 and NITR = 17. Here the number of integration points is equal to NITR.
      Figure 6.


    • For Isect10. The predefined section is defined with dimensions L1 to L6 and the number of integration points NITR defined between 0 to maximum value NITR_max (number of integrations points is limited to 100).
      Here is the list of required dimensions and maximum value of NITR for each predefined section:
      Table 1.
      Isect Shape (presented with NITR = 2) NITR_max Number of integration points (≤ 100)
      Isect = 10
      15 6 · N I T R + 9
      Isect = 11
      30 3 · N I T R + 9
      Isect = 12
      47 2 · N I T R + 5
      Isect = 13
      22 4 · N I T R + 9
      Isect = 14
      23 4 · N I T R + 8
      Isect = 15
      30 3 · N I T R + 9
      Isect = 16
      7 N I T R + 3 2
      Isect = 17
      2 4 · N I T R + 3 2
      Isect = 18
      2 4 · N I T R + 3 2
      Isect = 19
      15 6 · N I T R + 9
      Isect = 20
      7 N I T R + 3 2
      Isect = 21
      8 8 · N I T R + 12
      Isect = 22
      14 6 · N I T R + 13
      Isect = 23
      8 8 · N I T R + 16
      Isect = 24
      10 8 · N I T R + 14
      Isect = 25
      22 4 · N I T R + 12
      Isect = 26
      15 6 · N I T R + 17
      Isect= 27
      11 8 · N I T R + 12
      Isect= 28
      23 4 · N I T R + 8
      Isect= 29
      4 2 · N I T R + 3 2
      Isect= 30
      8 8 · N I T R + 14
      Isect= 31
      5 14 · N I T R + 22