ProductsAbaqus/StandardAbaqus/ExplicitAbaqus/CAE Beam section behaviorDefining a beam section's response to stretching, bending, shear, and torsion of the beam's axis requires a suitable definition of the axial force, N; bending moments, ${M}_{11}$ and ${M}_{22}$; and torque, T, as functions of the axial strain, $\epsilon $; curvature changes, ${\kappa}_{11}$ and ${\kappa}_{22}$; and twist, $\varphi $. Here the subscripts 1 and 2 refer to local, orthogonal axes in the beam section. If opensection beam types are used, the section behavior must also define the warping bimoment, W, and the generalized strain measures include the warping amplitude, w, and the bicurvature of the beam, $\chi $, which is the gradient of the warping amplitude with respect to position along the beam: $\chi =dw/dS$. The type of section definition you choose determines whether the beam section properties are recomputed during the progression of the analysis or established in the preprocessor for the duration of the analysis. If a general beam section definition is used (see Using a general beam section to define the section behavior), the crosssection properties are computed once, during preprocessing. Alternatively, a beam section definition that is integrated during the analysis can be used (see Using a beam section integrated during the analysis to define the section behavior), in which case Abaqus will use numerical integration of the stress over the crosssection to define the beam's response as the analysis proceeds. Since planar beams deform only in the X–Y plane, only N and ${M}_{11}$, and $\epsilon $ and ${\kappa}_{11}$, contribute to the response in these elements: ${\kappa}_{22}$, $\varphi $, and w are assumed to be zero. In Abaqus bending moments in beam sections are always measured about the centroid of the beam section, while torque is measured with respect to the shear center. The beam axis (defined as the line joining the nodes that define the beam element) need not pass through the centroid of the beam section. The degrees of freedom of the beam element are at the origin of the local $\left(1,2\right)$ coordinate system defined in the crosssection of the beam; that is, the line of the element connecting the element's nodes passes through the origin of the crosssection's local coordinate system. Determining whether to use a beam section integrated during the analysis or a general beam sectionWhen a beam section integrated during the analysis is used (see Using a beam section integrated during the analysis to define the section behavior), Abaqus integrates numerically over the section as the beam deforms, evaluating the material behavior independently at each point on the section. This type of beam section should be used when the section nonlinearity is caused only by nonlinear material response. When a general beam section is used (see Using a general beam section to define the section behavior), Abaqus precomputes the beam crosssection quantities and performs all section computations during the analysis in terms of the precomputed values. This method combines the functions of beam section and material descriptions (a material definition is not needed). The precomputed section properties may be specified in a variety of ways, including quite complex geometries defined with a twodimensional finite element mesh (see Meshed beam crosssections). A general beam section should be used when the beam section response is linear or when it is nonlinear and the nonlinearity arises from more than just material nonlinearity, such as in cases when section collapse occurs. Input File Usage Use the following option to define a beam section integrated during the analysis: BEAM SECTION Use the following option to define a general beam section: BEAM GENERAL SECTION Abaqus/CAE Usage To define a beam section integrated during the analysis: Property module: Create Section: select Beam as the section Category and Beam as the section Type: Section integration: During analysis To define a general beam section: Property module: Create Section: select Beam as the section Category and Beam as the section Type: Section integration: Before analysis Geometric section quantitiesThe section quantities described below are needed to define the behavior of a general beam section. Moments of inertiaThe moments of inertia with respect to the centroid are defined as $$\begin{array}{cc}\hfill {I}_{11}& ={\int}_{A}{\left({x}_{2}{x}_{2}^{c}\right)}^{2}dA,\hfill \\ \hfill {I}_{22}& ={\int}_{A}{\left({x}_{1}{x}_{1}^{c}\right)}^{2}dA,\text{and}\hfill \\ \hfill {I}_{12}& ={\int}_{A}\left({x}_{1}{x}_{1}^{c}\right)\left({x}_{2}{x}_{2}^{c}\right)dA,\hfill \end{array}$$
where (${x}_{1},{x}_{2}$) is the position of the point in the local $\left(1,2\right)$ beam section axis system and (${x}_{1}^{c},{x}_{2}^{c}$) is the position of the centroid of the crosssectional area. Bending stiffness and rotary inertia contributions for a meshed section profile (see Meshed crosssections) are calculated using the twodimensional crosssection model. The following integrated properties are defined for the entire crosssection model meshed with warping elements: $$\begin{array}{cc}\hfill {\left(EI\right)}_{11}& ={\int}_{A}E{\left({x}_{2}{x}_{2}^{c}\right)}^{2}dA,\hfill \\ \hfill {\left(EI\right)}_{22}& ={\int}_{A}E{\left({x}_{1}{x}_{1}^{c}\right)}^{2}dA,\hfill \\ \hfill {\left(EI\right)}_{12}& ={\int}_{A}E\left({x}_{1}{x}_{1}^{c}\right)\left({x}_{2}{x}_{2}^{c}\right)dA,\hfill \\ \hfill {\left(\rho I\right)}_{11}& ={\int}_{A}\rho {\left({x}_{2}{x}_{2}^{m}\right)}^{2}dA,\hfill \\ \hfill {\left(\rho I\right)}_{22}& ={\int}_{A}\rho {\left({x}_{1}{x}_{1}^{m}\right)}^{2}dA,\text{and}\hfill \\ \hfill {\left(\rho I\right)}_{12}& ={\int}_{A}\rho \left({x}_{1}{x}_{1}^{m}\right)\left({x}_{2}{x}_{2}^{m}\right)dA,\hfill \end{array}$$
where (${x}_{1}^{m},{x}_{2}^{m}$) is the center of mass of the cross section. Torsional constantThe torsional constant, J, depends on the shape of the cross section. The torsional constant of a circular section is the polar moment of inertia, $J={I}_{11}+{I}_{22}$. The torsional constant for the rectangular and trapezoidal library sections is calculated numerically by Abaqus using the Prandtl stress function approach. A local finite element model of the crosssection is created internally for this purpose. The number of integration points selected for the crosssection determines the accuracy of this finite element model. For increased accuracy specify a higherorder rule by selecting nondefault integration. The above rule is also applied to both the thinwalled box section and the arbitrary section to increase the accuracy of the model. If the thickness for each segment making up the section varies significantly, more integration points for the box section or smaller segments for the arbitrary section should be specified in the area where the thickness varies. The torsional stiffness for a meshed section is calculated over the twodimensional region meshed with warping elements. The accuracy of the integration depends on the number of elements in the model: $$\left(GJ\right)={\int}_{A}{G}_{\alpha \beta}\left({\psi}_{,\alpha}+{\u03f5}_{\gamma}^{\alpha}\left({x}^{\gamma}{x}_{s}^{\gamma}\right)\right)\left({\psi}_{,\beta}+{\u03f5}_{\delta}^{\beta}\left({x}^{\delta}{x}_{s}^{\delta}\right)\right)dA,$$
where ${\psi}_{,\alpha}$ denotes the derivative of the warping function with respect to the crosssection (1, 2) axis and ${x}_{\alpha}^{s}$ is the position of the shear center of the crosssectional area. All indices take values 1, 2. For more details, see Meshed beam crosssections. For closed thinwalled sections the torsional constant is calculated from $$J=\frac{4{A}_{c}^{2}}{\oint 1/tds},$$
where t is the thickness of the section, ${A}_{c}$ is the area enclosed by the median line of the section, and s is the length of the median line, measured along the circumference of the section in a counterclockwise direction. For open, builtup, thinwalled sections, $$J=\int \frac{1}{3}{t}^{3}ds.$$
Abaqus will check if a builtup section is closed or not and will use the appropriate torsional constant. Sectorial moment and warping constantFor open, thinwalled sections the sectorial moment is defined as $${\mathrm{\Gamma}}_{0}={\int}_{s}{S}_{w}tds$$
and the warping constant is defined as $${\mathrm{\Gamma}}_{W}={\int}_{s}{S}_{w}^{2}tds,$$
where ${S}_{w}$ is the sectorial area at a point in the section with the shear center as its pole. Rotary inertia for Timoshenko beamsIn general, the rotary inertia associated with torsional modes is different from that of flexural modes. For unsymmetric crosssections the rotary inertia is different in each direction of bending. For crosssections where the beam node is not located at the center of mass, coupling exists between the translational and rotational degrees of freedom. By default, the exact (anisotropic and coupled) rotary inertia is used for Timoshenko beams. In Abaqus/Standard the anisotropic rotary inertia introduces unsymmetric terms in the Jacobian operator during geometrically nonlinear, transient, directintegration dynamic simulations. If the rotary inertia effects are significant in the geometrically nonlinear dynamic response and the exact rotary inertia is used, the unsymmetric solver should be used for better convergence. Optionally, an approximate isotropic and uncoupled rotary inertia can be selected. In Abaqus/Standard this means that the rotary inertia associated with the torsional mode only is used for all rotational degrees of freedom; potentially destabilizing rotary inertia effects in impact problems due to the anisotropy or displacementrotation coupling will not be introduced. In Abaqus/Explicit this means a scaled flexural inertia with a scaling factor chosen to maximize the stable element time increment is used for all rotational degrees of freedom; i.e., the stable time increment will not be determined by the flexural response of the beam. In some slender beam analyses an isotropic approximation to the rotary inertia may be accurate enough. If beam elements are used to model platetype structures in Abaqus/Explicit (i.e., if the moment of inertia about one section axis of the beam is more than a thousand times greater than the moment of inertia about the other axis), the exact rotary inertia formulation may lead to a sharp cutback in the stable time increment. In this case it is recommended that you either use the isotropic approximation or alternatively consider modeling the structure with shell elements, which might be better suited to this type of analysis. For a definition of rotary inertia for the beam's crosssection, see Mass and inertia for Timoshenko beams. Input File Usage Use the following option to specify isotropic rotary inertia for a beam section integrated during the analysis: BEAM SECTION, ROTARY INERTIA=ISOTROPIC Use the following option to specify isotropic rotary inertia for a general beam section: BEAM GENERAL SECTION, ROTARY INERTIA=ISOTROPIC Abaqus/CAE Usage Isotropic rotary inertia for beam sections is not supported in Abaqus/CAE. The default exact rotary inertia is always used. Adding inertia to the beam section behavior for Timoshenko beamsAdditional mass and rotary inertia properties for Timoshenko beams (including PIPE elements) can be defined. This added inertia defined within the crosssection per unit length along the beam contributes to the inertia response of the beam without contributing to the structural stiffness. Additional beam inertia cannot be defined for a section if isotropic rotary inertia is used. To specify additional beam inertia, you define the mass (per unit length) with the mass center positioned at point $\left({x}_{1},{x}_{2}\right)$ in the local (1, 2) beam crosssection axis system. To include rotary inertia (per unit length), you can also define the angle $\alpha $ (in degrees) within the crosssection local (1, 2) system that positions the first axis of the rotary inertia coordinate system $\left(X,Y\right)$ relative to the local 1direction in the beam crosssection axis system. See Figure 1 for an illustration. Figure 1. Beam element with added inertia.
The rotary inertia components relative to the rotary inertia coordinate system $\left(X,Y\right)$ are defined as $$\begin{array}{cc}\hfill {I}_{XX}^{\rho}& ={\int}_{A}\rho {Y}^{2}dA,\hfill \\ \hfill {I}_{YY}^{\rho}& ={\int}_{A}\rho {X}^{2}dA,\text{and}\hfill \\ \hfill {I}_{XY}^{\rho}& ={\int}_{A}\rho XYdA,\hfill \end{array}$$
where A is the area, $\rho $ is the mass density, and X and Y are the local rotary inertia system coordinates measured from $\left({x}_{1},{x}_{2}\right)$, the center of the added mass contribution. As many point masses and rotary inertia contributions as are needed to define the added inertia can be specified. Mass proportional damping associated with the added inertia can be specified by assigning a value to the mass proportional Rayleigh damping coefficient, ${\alpha}_{R}$, or the composite damping coefficient, ${\xi}_{\alpha}$. Abaqus will use the mass weighted average (based on the material damping and the added inertia damping) for the element mass proportional damping. Input File Usage Use the following option in conjunction with the beam section definition to specify additional inertia properties: BEAM ADDED INERTIA, ALPHA=${\alpha}_{R}$, COMPOSITE=${\xi}_{\alpha}$ mass per unit length, ${x}_{1}$, ${x}_{2}$, $\alpha $, ${I}_{11}$, ${I}_{22}$, ${I}_{12}$ Abaqus/CAE Usage Additional inertia properties are not supported in Abaqus/CAE. Additional inertia due to immersion in fluidWhen a beam is fully or partially submerged, the effect of the surrounding fluid can be modeled as an additional distributed inertia on the beam (see Loading due to an incident dilatational wave field). By default, the beam is assumed to be fully submerged. Alternatively, you can specify that the added inertia per unit length should be reduced by a factor of onehalf to model a partially submerged beam. You specify the fluid mass density, ${\rho}_{f}$ (per unit volume); beam local x and y coordinates of the wetted crosssection centroid; wetted section effective radius, r; and empirical drag or flow coefficients, ${C}_{A}$ and ${C}_{AE}$. The inertia added per unit length to a fully immersed beam crosssection is given by $$\pi {r}^{2}{\rho}_{f}{C}_{A}.$$
Because the beam crosssection origin may not be coincident with the centroid of the wetted crosssection, the additional fluid inertia may include rotary effects. Nonzero values for the x and yoffsets of the wetted crosssection centroid will produce rotationdisplacement coupling in the inertia formulation. The default model for the added inertia derives from inviscid flow around a cylindrical crosssection (${C}_{A}=1.0$); you can specify a coefficient, ${C}_{A}$, that models flow around a different crosssection geometry. An immersed beam also experiences an additional added mass effect at its free ends. If a beam element's end node is not attached to any other element and additional fluid inertia is defined for this element, an additional mass may be added in the form: $$\frac{8}{3}{r}^{3}{\rho}_{f}{C}_{AE}.$$
For ${C}_{AE}=1.0$ this added mass corresponds to that of a hemispherical cap; the default value is ${C}_{AE}=0.0$. The coefficient ${C}_{AE}$ can be changed to model other geometries. If the beam is partially submerged, the end inertia is automatically reduced by onehalf. However, the added mass at the free ends is always isotropic: axial and transverse motions experience the same additional inertia. The virtual mass added to a submerged or partially submerged beam is not included in the total mass, center of mass, moments, or products of inertia reported in the data (.dat) file. Input File Usage Use the following option in conjunction with the beam section definition to define a fully immersed beam: BEAM FLUID INERTIA, FULL ${\rho}_{f}$, x, y, r, ${C}_{A}$, ${C}_{AE}$ Use the following option in conjunction with the beam section definition to define a partially immersed beam: BEAM FLUID INERTIA, HALF ${\rho}_{f}$, x, y, r, ${C}_{A}$, ${C}_{AE}$ Abaqus/CAE Usage To define a fully immersed beam: Property module: beam section editor: Fluid Inertia: toggle on Specify fluid inertia effects: Fully submerged To define a partially immersed beam: Property module: beam section editor: Fluid Inertia: toggle on Specify fluid inertia effects: Half submerged References
