Finite Element Analysis Method for Material Nonlinearity of Steel Tube Concrete Structures

Literature Overview

The paper by Wu Qingxiong, Chen Baochun, and Wei Jiangang (2008), published in Engineering Mechanics (Volume 25, Issue 6, pages 68-74), presents a novel finite element analysis method for capturing the material nonlinearity of steel tube concrete (SRC) structures. Funded by the Fujian Province Youth Science and Technology Talent Innovation Fund (2005J010), this research from Fuzhou University's College of Civil Engineering introduces a separation technique for the elastic-plastic stiffness matrix that enables efficient computation of material nonlinearity in beam elements.

Core Technical Content

The method applies the concept of generalized forces and generalized strains to beam element analysis, separating the elastic-plastic stiffness matrix into an elastic stiffness matrix and a plastic stiffness matrix. This separation allows the deformation of a beam element to be simply superimposed from elastic and plastic components. The internal forces are computed through the product of the elastic strain energy slope (elastic stiffness matrix) and displacement, enabling accurate and rapid calculation of structural unbalanced forces during incremental-iterative computation.

Methodology Framework

Key Technical Features

1. Deformation superposition: Elastic and plastic deformations are computed independently and superimposed, simplifying the nonlinear solution process.
2. Internal force computation: Internal forces are derived from the elastic stiffness matrix multiplied by displacement, rather than from the total stiffness matrix, which improves convergence.
3. Unbalanced force calculation: The method enables accurate and rapid computation of unbalanced forces after structural deformation, which is critical for incremental-iterative nonlinear analysis.
4. Double nonlinearity: The program NL_Beam3D can simultaneously consider material nonlinearity and geometric nonlinearity (P-Δ effects).

Technical Interpretation and Engineering Practice Relevance

Fiber Model and Material Behavior

The fiber model approach discretizes the cross-section into multiple material fibers, each of which follows its own constitutive law. This is particularly important for SRC structures where the steel tube and concrete core exhibit fundamentally different stress-strain behaviors:

• Steel tube: Bilinear or multilinear elastic-plastic behavior with strain hardening
• Concrete core: Nonlinear compression with confinement effect, potentially tensile cracking
• Interface: Slip and bond behavior between steel and concrete

The accuracy of the fiber model depends on the appropriate constitutive laws for each material. For the steel tube, the material model should account for the cyclic stress-strain behavior if seismic analysis is performed, including the Bauschinger effect and kinematic hardening.

Steel Pipe Manufacturing Relevance

The finite element model's accuracy is directly dependent on the material properties input, which in turn depend on the quality of the steel pipe material:

1. Yield strength variability: The actual yield strength of the steel tube may differ from the nominal grade (e.g., Q345) due to manufacturing variations. The fiber model should use the actual measured material properties, not just the nominal grade designation.
2. Strain hardening ratio: The ratio of ultimate tensile strength to yield strength (typically 1.15-1.30 for structural steels) affects the post-yield behavior and must be accurately captured in the material model.
3. Weld effects: The weld zone in a welded steel tube may have different mechanical properties from the base metal, which should be considered in the fiber model for critical applications.

Comparison with Conventional Methods

Study Insights and Reflections

This paper presents an elegant computational technique that addresses a fundamental challenge in nonlinear structural analysis: efficiently and accurately capturing material nonlinearity in beam elements. The separation of the elastic and plastic stiffness matrices is conceptually straightforward but computationally powerful, as it allows for better convergence behavior in incremental-iterative schemes. For engineers working on SRC structures, this method provides a practical tool for analyzing complex loading scenarios — including seismic loading, fire loading, and impact loading — where both material and geometric nonlinearities are significant.

The integration of the fiber model with the proposed stiffness separation technique is particularly valuable because it allows for a detailed representation of the SRC cross-section without the computational cost of a full 3D shell or solid element model. This makes it feasible to perform nonlinear analysis of large-scale SRC buildings within reasonable computational time. The program NL_Beam3D, which incorporates both material and geometric nonlinearities, represents a significant advancement over conventional linear or geometrically nonlinear-only analysis tools.

From a broader perspective, this research demonstrates the importance of computational methods in modern structural engineering. The ability to accurately predict the nonlinear behavior of SRC structures enables more efficient and safer design, reducing the need for conservative empirical design rules that may lead to over-engineering. For steel pipe manufacturers, this underscores the importance of providing detailed material property data — not just yield strength and ultimate strength, but also the full stress-strain curve, including the strain hardening region — to enable accurate finite element analysis of SRC structures. The methodology presented here is a valuable addition to the structural engineer's toolkit and should be considered for applications where the accuracy of nonlinear analysis is critical, such as the design of long-span SRC bridges, offshore platforms, and seismic isolation systems.