Using a General Beam Section to Define the Section Behavior

A general beam section:

  • is used to define beam section properties that are computed once and held constant for the entire analysis;

  • can be used to define linear or nonlinear section behavior;

  • for linear section behavior can be associated only with linear material behavior

  • enables the use of meshed cross-sections (Meshed Beam Cross-Sections); and

  • enables the use of tapered cross-sections (Abaqus/Standard only).

This page discusses:

Linear Section Behavior

Linear section response is calculated as follows. At each point in the cross-section the axial stress, σ, and the shear stress, τ, are given by

σ=E(θ¯,fβ)(ε-εth) and τ=G(θ¯,fβ)γ,

where

E(θ¯,fβ)

is Young's modulus (which may depend on the temperature, θ¯, and field variables, fβ, at the beam axis);

G(θ¯,fβ)

is the shear modulus (which may also depend on the temperature and field variables at the beam axis);

ε

is the axial strain;

γ

is the shear caused by twist; and

εth

is the thermal expansion strain.

The thermal expansion strain is given by

εth=α(θ¯,fβ)(θ-θ0)-α(θ¯I,fβI)(θI-θ0),

where

α(θ¯,fβ)

is the thermal expansion coefficient,

θ

is the current temperature at a point in the beam section,

fβ

are field variables,

θ0

is the reference temperature for α,

θI

is the initial temperature at this point (see Defining Initial Temperatures), and

fβI

are the initial values of the field variables at this point (see Defining Initial Values of Predefined Field Variables).

If the thermal expansion coefficient is temperature or field-variable dependent, it is evaluated at the temperature and field variables at the beam axis. Therefore, since we assume that θ varies linearly over the section, εth also varies linearly over the section.

The temperature is defined from the temperature of the beam axis and the gradients of temperature with respect to the local x1- and x2-axes:

θ=θ¯+θx1x1+θx2x2.

The axial force, N; bending moments, M1 and M2 about the 1 and 2 beam section local axes; torque, T; and bimoment, W, are defined in terms of the axial stress σ and the shear stress τ (see Beam element formulation). These terms are

N=E(A(εc-εcth)+Γ0χ),M1=E(I11(κ1-εthx2)-I12(κ2+εthx1)),M2=E(-I12(κ1-εthx2)+I22(κ2+εthx1)),T=GJϕ+GIpwp,W=E(Γ0(εc-εcth)+ΓWχ),

where

A

is the area of the section,

I11

is the moment of inertia for bending about the 1-axis of the section,

I12

is the moment of inertia for cross-bending,

I22

is the moment of inertia for bending about the 2-axis of the section,

J

is the torsional constant,

Γ0

is the sectorial moment of the section,

ΓW

is the warping constant of the section,

εc

is the axial strain measured at the centroid of the section,

εcth

is the thermal axial strain,

κ1

is the curvature change about the first beam section local axis,

κ2

is the curvature change about the second beam section local axis,

ϕ

is the twist,

χ

is the bicurvature defining the axial strain in the section due to the twist of the beam, and

wp=wf-w

is the difference between the unconstrained warping amplitude, wf, and the actual warping amplitude, w.

Γ0, ΓW, χ, and wp are nonzero only for open-section beam elements.

Defining Linear Section Behavior for Library Cross-Sections or Linear Generalized Cross-Sections

Linear beam section response is defined geometrically by A, I 11 , I 12 , I 22 , J, and (if necessary) Γ 0 and Γ W .

You can input these geometric quantities directly or specify a standard library section and Abaqus calculates these quantities. In either case define the orientation of the beam section (see Beam Element Cross-Section Orientation).

You can specify Young's modulus, the shear modulus, and the coefficient of thermal expansion as functions of temperature; and associate the section properties with a region of your model. If the thermal expansion coefficient is temperature dependent, the reference temperature for thermal expansion must also be defined as described later in this section.

Alternatively, you can associate a material definition (Material Data Definition) with the section definition. Abaqus determines the equivalent section properties. You must associate the section behavior with a region of your model.

Specifying the Geometric Quantities Directly

You can define “generalized” linear section behavior by specifying A, I11, I12, I22, J, and—if necessary—Γ0 and ΓW directly. In this case you can specify the location of the centroid, thus allowing the bending axis of the beam to be offset from the line of its nodes. In addition, you can specify the location of the shear center.

Specifying a Standard Library Section and Allowing Abaqus to Calculate the Geometric Quantities

You can select one of the standard library sections (see Beam Cross-Section Library) and specify the geometric input data needed to define the shape of the cross-section. Abaqus then calculates the geometric quantities needed to define the section behavior automatically. In addition, you can specify an offset for the section origin.

Specifying the Linear Section Response with a Material Definition

The material definition (Material Data Definition) can contain isotropic linear elastic behavior (Linear Elastic Behavior) and isotropic thermal expansion behavior (Thermal Expansion). If both the isotropic linear elastic material behavior and the isotropic thermal expansion behavior are temperature or field variable dependent, the values of the independent variables (temperature or field variables) that you specify must be the same for both the linear elastic moduli and thermal expansion coefficient. You can specify damping behavior (Material Damping). In Abaqus/Explicit you must define the density (Density) of the material. In an Abaqus/Standard analysis the density is needed only when the mass of the beam elements is required. Any nonlinear material properties (such as plastic behavior) are ignored.

Defining Linear Section Behavior for Meshed Cross-Sections

Linear beam section response for a meshed section profile is obtained by numerical integration from the two-dimensional model. The numerical integration is performed once, determining the beam stiffness and inertia quantities, as well as the coordinates of the centroid and shear center, for the duration of the analysis. These beam section properties are calculated during the beam section generation and are written to the text file jobname.bsp. This text file can be included in the beam model. See Meshed Beam Cross-Sections for a detailed description of the properties defining the linear beam section response for a meshed section, as well as for how a typical meshed section is analyzed.

Defining Linear Section Behavior for Tapered Cross-Sections in Abaqus/Standard

In Abaqus/Standard you can define Timoshenko beams with linearly tapered cross-sections. General beam sections with linear response and standard library sections are supported, with the exception of arbitrary sections. The section parameters are defined at the two end nodes of each beam element. The effective beam area and moment of inertia for bending about the 1- and 2-axis of the section used in the calculation of the beam stiffness matrix, section forces, and stresses are

Aeff=AI+AIAJ+AJ3,I11eff=I11I+(I11I)3I11J4+I11II11J+I11I(I11J)34+I11J5,I22eff=I22I+(I22I)3I22J4+I22II22J+I22I(I22J)34+I22J5,

where the superscripts I and J refer to the two end nodes of the beam. The remaining effective geometric quantities are calculated as the average between the values at the two end nodes. This approximation suffices for mild tapering along each element, but it can lead to large errors if the tapering is not gradual. Abaqus/Standard issues a warning message during input file preprocessing if the area or inertia ratio is larger than 2.0 and an error message if the ratio is larger than 10.0.

The effective area and inertia are not used in the computation of the mass matrix. Instead, terms on the diagonal quadrants use the properties from the respective nodes, while off-diagonal quadrants use averaged quantities. For example, the axial inertia a linear element would have the diagonal term coming from node I of ρAI/3, while node J contributes with ρAJ/3 and the two off-diagonal contributions equal ρ(AI+AJ)/12. Mild tapering is assumed in this formulation, since the total mass of the element totals ρ(AI+AJ)/2.

Nonlinear Section Behavior

Typically nonlinear section behavior is used to include the experimentally measured nonlinear response of a beam-like component whose section distorts in its plane. When the section behaves according to beam theory (that is, the section does not distort in its plane) but the material has nonlinear response, it is usually better to use a beam section integrated during the analysis to define the section geometrically (see Using a Beam Section Integrated during the Analysis to Define the Section Behavior), in association with a material definition.

Nonlinear section behavior can also be used to model beam section collapse in an approximate sense: Nonlinear dynamic analysis of a structure with local inelastic collapse illustrates this for the case of a pipe section that may suffer inelastic collapse due to the application of a large bending moment. In following this approach you should recognize that such unstable section collapse, like any unstable behavior, typically involves localization of the deformation: results, therefore, are strongly mesh sensitive.

Calculation of Nonlinear Section Response

Nonlinear section response is assumed to be defined by

N=N(εc-εcth,θ¯,fβ),M1=M1(κ1,θ¯,fβ),M2=M2(κ2,θ¯,fβ),T=T(ϕ,θ¯,fβ),

where (    ) means a functional dependence on the conjugate variables: N=N(ε), M1=M1(κ1), etc. For example, N(εc-εcth,θ¯,fβ) means that N is a function of:εc-εcth; θ¯, the temperature of the beam axis; and of fβ, any predefined field variables at the beam axis. When the section behavior is defined in this way, only the temperature and field variables of the beam axis are used: any temperature or field-variable gradients given across the beam section are ignored.

These nonlinear responses may be purely elastic (that is, fully reversible—the loading and unloading responses are the same, even though the behavior is nonlinear) or may be elastic-plastic and, therefore, irreversible.

The assumption that these nonlinear responses are uncoupled is restrictive; in general, there is some interaction between these four behaviors, and the responses are coupled. You must determine if this approximation is reasonable for a particular case. The approach works well if the response is dominated by one behavior, such as bending about one axis. However, it may introduce additional errors if the response involves combined loadings.

Defining Nonlinear Section Behavior

You can define “generalized” nonlinear section behavior by specifying the area, A; moments of inertia, I11 for bending about the 1-axis of the section, I22 for bending about the 2-axis of the section, and I12 for cross-bending; and torsional constant, J. These values are used only to calculate the transverse shear stiffness; and, if needed, A is used to compute the mass density of the element. In addition, you can define the orientation and the axial, bending, and torsional behavior of the beam section (N, M1, M2, T), as well as the thermal expansion coefficient. If the thermal expansion coefficient is temperature dependent, the reference temperature for thermal expansion must also be defined as described below.

Nonlinear generalized beam section behavior cannot be used with beam elements with warping degrees of freedom.

The axial, bending, and torsional behavior of the beam section and the thermal expansion coefficient are defined by tables. See Material Data Definition for a detailed discussion of the tabular input conventions. In particular, you must ensure that the range of values given for the variables is sufficient for the application since Abaqus assumes a constant value of the dependent variable outside this range.

Defining Linear Response for N, M1, M2, and T

If the particular behavior is linear, N, M1, M2, and T should be specified as functions of the temperature and predefined field variables, if appropriate.

As an example of axial behavior, if

N=(AE)(εc-εcth),

where (AE) is constant for a given temperature, the value of (AE) is entered. (AE) can still be varied as a function of temperature and field variables.

Defining Nonlinear Elastic Response for N, M1, M2, and T

If the particular behavior is nonlinear but elastic, the data should be given from the most negative value of the kinematic variable to the most positive value, always giving a point at the origin. See Figure 1 for an example.

Example of elastic nonlinear beam section behavior definition.

Defining Elastic-Plastic Response for N, M1, M2, and T

By default, elastic-plastic response is assumed for N, M1, M2, and T.

The inelastic model is based on assuming linear elasticity and isotropic hardening (or softening) plasticity. The data in this case must begin with the ( 0 , 0 ) point and proceed to give positive values of the kinematic variable at increasing positive values of the conjugate force or moment. Strain softening is allowed. The elastic modulus is defined by the slope of the initial line segment, so that straining beyond the point that terminates that initial line segment is partially inelastic. If strain reversal occurs in that part of the response, it is elastic initially. See Figure 2 for an example.

Example of inelastic nonlinear beam section behavior definition.

Defining the Reference Temperature for Thermal Expansion

The thermal expansion coefficient may be temperature dependent. In this case the reference temperature for thermal expansion, θ0, must be defined.

Defining the Initial Section Forces and Moments

You can define initial stresses (see Defining Initial Stresses) for general beam sections that are applied as initial section forces and moments. Initial conditions can be specified only for the axial force, the bending moments, and the twisting moment. Initial conditions cannot be prescribed for the transverse shear forces.

Defining a Change in Cross-Sectional Area due to Straining

In the shear flexible elements Abaqus provides for a possible uniform cross-sectional area change by allowing you to specify an effective Poisson's ratio for the section. This effect is considered only in geometrically nonlinear analysis (see Defining an Analysis) and is provided to model the reduction or increase in the cross-sectional area for a beam subjected to large axial stretch.

The value of the effective Poisson's ratio must be between −1.0 and 0.5. By default, this effective Poisson's ratio for the section is set to 0.0 so that this effect is ignored. Setting the effective Poisson's ratio to 0.5 implies that the overall response of the section is incompressible. This behavior is appropriate if the beam is made of rubber or if it is made of a typical metal whose overall response at large deformation is essentially incompressible (because it is dominated by plasticity). Values between 0.0 and 0.5 mean that the cross-sectional area changes proportionally between no change and incompressibility, respectively. A negative value of the effective Poisson's ratio results in an increase in the cross-sectional area in response to tensile axial strains.

This effective Poisson's ratio is not available for use with Euler-Bernoulli beam elements.

Defining Damping

When the beam section and material behavior are defined by a general beam section, you can include mass and viscous stiffness proportional damping in the dynamic response (calculated in Abaqus/Standard with the direct time integration procedure, Implicit Dynamic Analysis Using Direct Integration).

See Material Damping for more information about the material damping types available in Abaqus.

Specifying Temperature and Field Variables

Define temperatures and field variables by giving the values at the origin of the cross-section as either predefined fields or initial conditions (see Predefined Fields or Initial Conditions). Temperature gradients can be specified in the local 1- and 2-directions; other field-variable gradients defined through the cross-section are ignored in the response of beam elements that use a general beam section definition.

Output

Only the section forces, moments, and transverse shear forces and section strains, curvatures, and transverse shear strains can be output (see Element Output and Writing Element Output to the Output Database).

You can output stress and strain at particular points in the section. For linear section behavior defined using a standard library section or a generalized section, only axial stress and axial strain values are available. For linear section behavior defined using a meshed section, axial and shear stress and strain are available. For nonlinear generalized section behavior, axial strain output only is provided.

Specifying the Output Section Points for Standard Library Sections and Generalized Sections

To locate points in the section at which output of axial strain (and, for linear section behavior, axial stress) is required, specify the local (x1,x2) coordinates of the point in the cross-section: Abaqus numbers the points 1, 2, … in the order that they are given.

The variation of ε over the section is given by

ε=εc+κ1(x2-x2c)-κ2(x1-x1c),

where (x1c,x2c) are the local coordinates of the centroid of the beam section and κ1 and κ2 are the changes of curvature for the section.

For open-section beam element types, the variation of ε over the section has an additional term of the form ψ(x1,x2)χ, where ψ(x1,x2) is the warping function. The warping function itself is undefined in the general beam section definition. Therefore, Abaqus will not take into account the axial strain due to warping when calculating section points output. Axial strains due to warping are included in the stress/strain output if a beam section integrated during the analysis is used.

Abaqus uses St. Venant torsion theory for noncircular solid sections. The torsion function and its derivatives are necessary to calculate shear stresses in the plane of the cross-section. The function and its derivatives are not stored for a general beam section. Therefore, you can request output of axial components of stress/strain only. A beam section integrated during the analysis must be used to obtain output of shear stresses.

Requesting Output of Maximum Axial Stress/Strain in Abaqus/Standard

If you specify the output section points to obtain the maximum axial stress/strain (MAXSS) for a linear generalized section, the output value is the maximum of the values at the user-specified section points. You must select enough section points to ensure that this is the true maximum. MAXSS output is not available for nonlinear generalized sections or for an Abaqus/Explicit analysis.

Specifying the Output Section Points for Meshed Cross-Sections

For meshed cross-sections you can indicate in the two-dimensional cross-section analysis the elements and integration points where the stress and strain will be calculated during the subsequent beam analysis. Abaqus will then add the section points specification to the resulting jobname.bsp text file. This text file is then included as the data for the general beam section definition in the subsequent beam analysis. See Meshed Beam Cross-Sections for details.

The variation of the axial strain ε over the meshed section is given by

ε=εc+κ1(x2-x2c)-κ2(x1-x1c),

where (x1c,x2c) are the local coordinates of the centroid of the beam section and κ1 and κ2 are the changes of curvature for the section.

The variations of shear components γ1 and γ2 over the meshed section are given by

γ1=γs1+ϕ(Ψx1-(x2-x2s)),
γ2=γs2+ϕ(Ψx2+(x1-x1s)),

where (x1s,x2s) are the local coordinates of the shear center of the beam section, ϕ is the twist of the beam axis, Ψ(x1,x2) is the warping function, and γs1 and γs2 are shear strains due to the transverse shear forces.

For the case of an orthotropic composite beam material, the axial stress σ and the two shear components τ1 and τ2 are calculated in the beam section (1, 2) axis as follows:

{στ1τ2}=[E00G1(cosα)2+G2(sinα)2(G1-G2)cosαsinαsymG1(sinα)2+G2(cosα)2]{εγ1γ2},

where α determines the material orientation.