Material Inelasticity

Distributed inelasticity elements are becoming widely employed in earthquake engineering applications, either for research or professional engineering purposes. Whilst their advantages in relation to the simpler lumped-plasticity models, together with a concise description of their historical evolution and discussion of existing limitations, can be found in e.g. Filippou and Fenves [2004] or Fragiadakis and Papadrakakis [2008], here it is simply noted that distributed inelasticity elements do not require (not necessarily straightforward) calibration of empirical response parameters against the response of an actual or ideal frame element under idealized loading conditions, as is instead needed for concentrated-plasticity phenomenological models.

In SeismoBuild, use is made of the so-called fibre approach to represent the cross-section behaviour, where each fibre is associated with a uniaxial stress-strain relationship; the sectional stress-strain state of beam-column elements is then obtained through the integration of the nonlinear uniaxial stress-strain response of the individual fibres (by default 150) in which the section has been subdivided (the discretisation of a typical reinforced concrete cross-section is depicted, as an example, in the figure below).

Such models feature additional assets, which can be summarized as: no requirement of a prior moment-curvature analysis of members; no need to introduce any element hysteretic response (as it is implicitly defined by the material constitutive models); direct modelling of axial load-bending moment interaction (both on strength and stiffness); straightforward representation of biaxial loading, and interaction between flexural strength in orthogonal directions.

Distributed inelasticity frame elements can be implemented with two different finite elements (FE) formulations: the classical displacement-based (DB) ones [e.g. Hellesland and Scordelis 1981; Mari and Scordelis 1984], and the more recent force-based (FB) formulations [e.g. Spacone et al. 1996; Neuenhofer and Filippou 1997].

In a DB approach the displacement field is imposed, whilst in a FB element equilibrium is strictly satisfied and no restraints are placed to the development of inelastic deformations throughout the member; see e.g. Alemdar and White [2005] and Freitas et al. [1999] for further discussion. In the DB case, displacement shape functions are used, corresponding for instance to a linear variation of curvature along the element.

In contrast, in a FB approach, a linear moment variation is imposed, i.e. the dual of previously referred linear variation of curvature. In the case of material inelasticity, imposing a displacement field does not enable to capture the real deformed shape since the curvature field can be, in a general case, highly nonlinear. For linear elastic material behaviour, the two approaches obviously produce the same results, provided that only nodal forces act on the element. On the contrary, in case of material inelasticity, imposing a displacement field does not enable to capture the real deformed shape since the curvature field can be, in a general case, highly nonlinear. In this situation, with a DB formulation a refined discretisation (meshing) of the structural element (typically 4-5 elements per structural member) is required for the computation of nodal forces/displacements, in order to accept the assumption of a linear curvature field inside each of the sub-domains. Still, in the latter case users are not advised to rely on the values of computed sectional curvatures and individual fibre stress-strain states.

Instead, a FB formulation is always exact, since it does not depend on the assumed sectional constitutive behaviour. In fact, it does not restrain in any way the displacement field of the element. In this sense this formulation can be regarded as always "exact". In this sense this formulation can be regarded as always "exact", the only approximation being introduced by the discrete number of the controlling sections along the element that are used for the numerical integration. A number of 4 Gauss-Lobatto integration sections are required to avoid under-integration, option which will in general simulate the spread of inelasticity in an acceptable way. Such feature enables to model each structural member with a single FE element, therefore allowing a one-to-one correspondence between structural members (beams and columns) and model elements. In other words, no meshing is theoretically required within each element, even if the cross section is not constant. This is because the force field is always exact, regardless of the level of inelasticity.

In SeismoBuild, both aforementioned DB and FB element formulations are implemented, with the latter being typically recommended, since, as mentioned above, it does not in general call for element discretisation, thus leading to considerably smaller models, with respect to when DB elements are used, and thus much faster analyses, notwithstanding the heavier element equilibrium calculations. An exception to this non-discretisation rule arises when localisation issues are expected, in which case special cautions/measures are needed, as discussed in Calabrese et al. [2010].

Finally, it is noted that, for reasons of higher accuracy, the Lobatto quadrature is used. The approximate coordinates along the element's length (measured from its baricentre) of the integration sections are for 4 integration sections: [-1  -0.447  0.447  1] x L/2.

Notes

1. Users are also invited to read the NEHRP Seismic Design Technical Brief No. 4 (i.e. Deierlein G.G., Reinhorn A.M., and Willford M.R. [2010]), in which the nonlinear modelling is well covered.
2. The Inelastic displacement-based frame element type - infrmDB is used only for members of small length, e.g. to model short columns.