Skip to main content

MINI REVIEW article

Front. Built Environ., 14 November 2022
Sec. Earthquake Engineering
This article is part of the Research Topic Rising Stars in Built Environment View all 10 articles

Phenomenological rate-independent uniaxial hysteretic models: A mini-review

Raffaele Capuano
Raffaele Capuano1*Aguinaldo FraddosioAguinaldo Fraddosio2Mario Daniele PiccioniMario Daniele Piccioni2
  • 1Department of Structures for Engineering and Architecture, University of Naples Federico II, Napoli, Italy
  • 2Department of Civil Engineering Sciences and Architecture, Polytechnic University of Bari, Bari, Italy

A great variety of phenomenological models has been proposed over the years to model rate-independent hysteretic forces in structural mechanics. The classification of such models is usually based on the type of equation that needs to be solved to evaluate the output variable. In particular, we distinguish among algebraic, transcendental, differential and integral models. For algebraic (transcendental) models, an algebraic (a transcendental) equation needs to be solved to compute the output variable; conversely, differential equations are employed for differential models, whereas equations expressed in integral form characterize integral models. This paper provides a mini-review of the most adopted phenomenological rate-independent uniaxial hysteretic models. Such models are selected in order to provide a complete overview of the four types of previously mentioned models, currently available in the literature. In particular, we illustrate the fundamental characteristics of each model and discuss their peculiarities in terms of 1) number of adopted parameters and variables, 2) physical interpretation of parameters and related calibration procedures, 3) type of hysteresis loop shapes that can be simulated.

1 Introduction

The mathematical modeling of hysteresis phenomena is a challenging problem addressed during the years by different authors by means of different models. The main aim of the various authors was to obtain accurate and computationally efficient models based on a small number of parameters (Vaiana et al., 2019b). The model proposed in the literature can be classified into four main categories based on the type of equation that need to be solved to compute the output variable (Vaiana et al., 2018; Vaiana et al., 2021b). In particular, we have: 1) algebraic, 2) transcendental, 3) differential, 4) integral models.

2 Algebraic models

2.1 Ramberg-Osgood model

The Ramberg-Osgood model is an algebraic model formulated by Ramberg and Osgood (1943). In this model, generally used for simulating the behavior of steel, the generalized displacement u is assumed as the output variable and is expressed as the sum of two terms:

u=ue+up,(1)

an elastic and a plastic generalized displacement indicated, respectively, as ue and up. In particular, the generalized displacement is evaluated in a closed-form as:

u=fE+KfEn,(2)

where E is the elastic modulus, K and n are constants that depend on the considered material. The first term on the right-hand side of Eq. 2 is equal to the elastic part of the generalized displacement, whereas the second term accounts for the plastic part. The parameters K and n describe the hardening behavior of the material. Introducing the generalized yield force of the material fy, and defining a new parameter α as:

α=KfEn1,(3)

it is possible to rewrite the general expression as follows:

u=fE+αfEffyn1.(4)

Since Eq. 4 defines the skeleton curve, it is necessary to adopt the Masing rule (Masing and Mauksch, 1926) in order to simulate an entire hysteresis loop.

2.2 Giuffrè-Menegotto-Pinto model

As reported in Carreño et al. (2020), the model was first developed by Giuffrè (1970), and was based on the nonlinear stress-strain relation proposed by Goldberg and Richard (1963) and incorporates the effect of plastic deformations on the Bauschinger effect observed in steel tested experimentally. The formulation was further improved by Menegotto and Pinto (1973) and subsequently used by multiple authors due to its simplicity and accuracy in predicting the response of reinforcing steel. Filippou et al. (1983) later incorporated the effect of isotropic hardening into the constitutive law. The formulation proposed by Filippou is often preferred in the modeling of reinforced concrete systems given the consistent response and limited failures in convergence observed in the nonlinear analysis of large structural models. The constitutive relation of the model is:

f=Eh+EEh1+uuyr1ru,(5)

where E is the elastic modulus, Eh is the post-elastic stiffness, uy is the yield generalized displacement and r is a parameter which controls the rate of variation of the tangent stiffness from E to Eh. In particular, the expression for the evaluation of the parameter r is:

r=r0a1uuya2+uuy.(6)

Since Eq. 5 defines the skeleton curve, as done for the Ramberg-Osgood model, it is necessary to adopt the Masing rule in order to simulate an entire hysteresis loop.

2.3 Algebraic model

The model proposed by Vaiana et al. (2019a) is a uniaxial phenomenological model which can simulate hysteresis loops limited by two parallel straight lines or curves by adopting a set of only five parameters easy to calibrate (Sessa et al., 2020). Furthermore, the model allows one to evaluate the output variable by means of a closed form expression having an algebraic nature.

In particular, the generalized force f, during the generic loading phase (u̇>0), can be evaluated as:

fu,uj+=c+u,uj+ifu<uj+cuuifu>uj+,(7)

whereas, during the generic unloading phase (u̇<0), it can be computed as:

fu,uj=cu,ujifu>ujcluifu<uj.(8)

In Eqs. 7 and 8, c+ and c represent, respectively, the generic loading and unloading curves whose expressions are:

c+u,uj+=β1u3+β2u5+kbu+f0+kakb1α1+uuj++2u01α1+2u01α,(9)
cu,uj=β1u3+β2u5+kbuf0+kakbα11u+uj+2u01α1+2u01α,(10)

whereas cu and cl are, respectively, the upper and lower limiting curves having expressions:

cuu=β1u3+β2u5+kbu+f0,(11)
clu=β1u3+β2u5+kbuf0.(12)

The internal variable uj+ (uj), characterizing the generic loading (unloading) phase, is given by:

uj+=1+uP+2u0{1αkakb[fPβ1uP3β2uP5kbuPf0+kakb1+2u01α1α]}11α,(13)
uj=1+uP2u0+{α1kakb[fPβ1uP3β2uP5kbuP+f0+kakb1+2u01αα1]}11α,(14)

where uP and fP are the coordinates of the initial point of c+ (c).

The internal model parameters u0 and f0 can be expressed in terms of the parameters ka, kb, and α as follows:

u0=12kakbδk1α1,(15)
f0=kakb21α1+2u01α1,(16)

where δk is the difference between the two different values assumed by the transverse tangent stiffness at uj+ (uj). However, the authors suggests, based on the results of several numerical tests (Vaiana et al., 2020; Vaiana et al., 2021a; Pellecchia et al., 2022), to set δk = 10–20.

3 Transcendental models

3.1 Kikuchi-Aiken model

The transcendental model proposed by Kikuchi and Aiken (1997) was formulated with the purpose of accurately predicting the response of elastomeric bearings (Losanno et al., 2019; Losanno et al., 2022a; Losanno et al., 2022b; Losanno et al., 2022c).

The generalized force f is expressed as the sum of two terms:

fu=f1u+f2u,(17)

where

f1u=121fufmfmu+sgnu|u|n,(18)

and

f2u=fu12ea1+uum+b1+uumec1+uumifu̇>0fu12ea1uum+b1uumec1uumifu̇<0.(19)

In particular, fu is the generalized force at zero generalized displacements, fm is the generalized peak force on the skeleton curve, um is the generalized peak displacement on the skeleton curve and the parameter n specifies the stiffening phenomenon. The model parameters a, b, and c are obtained by imposing the following expressions:

1e2aa=1πheqfm2fu,(20)
b=c2πheqfmfu2+2ae2a1,(21)

where heq is the equivalent viscous damping ratio, evaluated from an empirical formula as a function of the generalized force, which is determined from the results of tests performed on individual bearings. Both Eqs. 20 and 21 are derived assuming that the analytical hysteresis loop area matches the experimental one.

Eq. 19 are derived for application to steady-state hysteresis behaviors for elastomeric bearings. Masing rule is applied to fully define the generalized force under a randomly-varying displacement, and Eq. 19 are replaced by:

f2u=f2i+fu22eauuium+buuiumecuuiumifu̇>0f2ifu22eauuiumbuuiumecuuiumifu̇<0,(22)

where

f2i=fif1,(23)

and (ui, fi) is the most recent point of load reversal.

3.2 Vaiana-Rosati model

The model proposed by Vaiana and Rosati (2023) is a new uniaxial phenomenological model which can simulate a great number of rate-independent hysteretic responses characterized by symmetric, asymmetric, pinched, S-shaped, flag-shaped hysteresis loops or by a combination of them. The model is based on two sets of eight parameters controlling the loading and unloading phase in a separate way. Furthermore, the model allows one to evaluate the output variable by means of a closed form expression having an exponential nature.

In particular, the generalized force f, during the generic loading phase (u̇>0), can be evaluated as:

fu,uj+=c+u,uj+ifu<uj+cuuifu>uj+,(24)

whereas, during the generic unloading phase (u̇<0), it can be computed as:

fu,uj=cu,ujifu>ujcluifu<uj.(25)

In Eqs. 24 and 25, c+ and c represent, respectively, the generic loading and unloading curves whose expressions are:

c+u,uj+=β1+eβ2+uβ1++4γ1+1+eγ2+uγ3+2γ1++kb+u+f0+1α+eα++uuj++ū+eα+ū+,(26)
cu,uj=β1eβ2uβ1+4γ11+eγ2uγ32γ1+kbuf0+1αeαu+uj+ūeαū,(27)

whereas cu and cl are, respectively, the upper and lower limiting curves having expressions:

cuu=β1+eβ2+uβ1++4γ1+1+eγ2+uγ3+2γ1++kb+u+f0+,(28)
clu=β1eβ2uβ1+4γ11+eγ2uγ32γ1+kbuf0.(29)

The internal variable uj+ (uj), characterizing the generic loading (unloading) phase, is given by:

uj+=uP+ū++1α+ln{+α+[β1+eβ2+uPβ1++4γ1+1+eγ2+uPγ3+2γ1++kb+uP+f0++1α+eα+ū+fP]},(30)
uj=uPū1αln{α[β1eβ2uPβ1+4γ11+eγ2uPγ32γ1+kbuPf01αeαūfP]},(31)

where uP and fP are the coordinates of the initial point of c+ (c).

4 Differential models

4.1 Bouc-Wen model

The Bouc-Wen model is a uniaxial rate-independent hysteretic model, having a differential nature, that has been originally formulated by Bouc (1967) and subsequently extended by Wen (1976). The model is capable of simulating symmetric generalized force-displacement hysteresis loops (Visintin, 2013; Carboni et al., 2018). In the formulation of the model proposed by Wen (1980), the differential equation that defines the model and that needs to be solved is:

ḟ=βu̇fnγu̇fn1f+Au̇,(32)

where β, γ, A and n influence the generalized force-displacement hysteresis loop shape. In particular, parameters γ and β control the shape of the hysteresis loop, A the restoring force amplitude, and n the rate of variation between the initial and the asymptotic tangent stiffness, a large value of n corresponds to an almost bilinear hysteresis loop.

4.2 Ozdemir model

The model was proposed by Ozdemir (1976) and illustrated in more detail in Graesser and Cozzarelli (1991). The model proposed by Ozdemir is a modified form of the Bouc-Wen model. This model is general enough to include a variety of different types of hysteretic behavior. In addition, the material model is linked to experimental studies of the properties of a nickel-titanium shape-memory alloy (SMA). The equations that describe the model are:

ḟ=Eu̇u̇fβYnβ̇=αEu̇fβYn,(33)

where f is the generalized force, u is the generalized displacement, β is the generalized back-force, E is the elastic modulus, Y is the generalized force correspondent to the point obtained as the intersection of the upper limiting straight line and the tangent line to the origin, n is a constant parameter controlling the rate of variation of the tangent stiffness from E to the asymptotic value Ey, and α is a constant controlling the tangent stiffness of the hysteresis loop and is given by:

α=EyEEy.(34)

5 Integral models

5.1 Preisach model

The Preisach model, first derived in Preisach (1935) and than described in Kuczmann (2010) and Mayergoyz (1986), is an integral model given by:

ft=αβμα,βγα,β,utdαdβ.(35)

The model is obtained considering an infinite set of simplest hysteresis operators γ. These operators can be represented by rectangular loops on the input-output plane. Numbers α and β correspond to up and down switching values of input. The outputs of these operators may assume only two values, + 1 and −1.

Along with the set of operators γ, the model considers an arbitrary weight function μ(α, β). To determine μ(α, β), the set of first order reversal curves should be experimentally found. This can be done by bringing first the input to such a value that outputs of all operators γ are equal to −1. Now if we gradually increase the input value, then we will follow along a limiting ascending branch. The notation fα will be used for the output value on this branch corresponding to the input u = α. The notation fαβ will be used for the output value on the transition curve attached to the limiting ascending branch at the point fα. This output value corresponds to the input u = β. At this point we can define the function:

Fα,β=fαfαβ2.(36)

Finally it can be proved that:

μα,β=2Fα,βαβ.(37)

6 Discussion

In this review, we made an overview of some celebrated and recently formulated models to illustrate their fundamental characteristics with particular attention to the number of parameters that need to be set (Table 1) and, if present, their physical interpretation. We also show the different kinds of hysteresis loop shapes, shown in Figure 1, that can be simulated by using the different models.

TABLE 1
www.frontiersin.org

TABLE 1. Uniaxial phenomenological models.

FIGURE 1
www.frontiersin.org

FIGURE 1. Hysteresis loops shapes.

In particular, from Table 1 it can be noted that:

• differential models require a small number of parameters but can simulate only simple hysteresis loops shape;

• models such as the algebraic Ramberg-Osgood, the Giuffrè-Menegotto-Pinto, and the transcendental Kikuchi-Aiken, define only the skeleton curve of a loop. Therefore, it is necessary to use the Masing rule to simulate an entire hysteresis loop, which involves the addition of a series of internal variables;

• the Vaiana-Rosati model, despite the number of parameters, allows for the simulation of the vast majority of hysteresis loop shapes through closed-form expressions;

• the hysteresis loop shape in Figure 1I, typical of the magnetic hysteresis, can be simulated only by the Preisach model, which presents significant difficulties in the evaluation of the output variable.

Author contributions

RC: Writing—original draft, Writing—review and editing. AF: Supervision. MP: Supervision.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer GM declared a past co-authorship with the authors AF, MP to the handling editor.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Bouc, R. (1967). “Forced vibration of mechanical systems with hysteresis,” in Proceedings of the 4th Conference on Nonlinear Oscillation, Prague, Czechoslovakia, September 5–9, 1967.

Google Scholar

Carboni, B., Lacarbonara, W., Brewick, P. T., and Masri, S. F. (2018). Dynamical response identification of a class of nonlinear hysteretic systems. J. Intelligent Material Syst. Struct. 29, 2795–2810. doi:10.1177/1045389x18778792

CrossRef Full Text | Google Scholar

Carreño, R., Lotfizadeh, K., Conte, J., and Restrepo, J. (2020). Material model parameters for the giuffrè-menegotto-pinto uniaxial steel stress-strain model. J. Struct. Eng. (N. Y. N. Y). 146, 04019205. doi:10.1061/(asce)st.1943-541x.0002505

CrossRef Full Text | Google Scholar

Filippou, F. C., Popov, E. P., and Bertero, V. V. (1983). Effects of bond deterioration on hysteretic behavior of reinforced concrete joints. EERC Report 83/19, Berkeley: Earthquake Engineering Research Center, University of California.

Google Scholar

Giuffrè, A. (1970). Il comportamento del cemento armato per sollecitazioni cicliche di forte intensità. G. del Genio Civ. 28, 1–20.

Google Scholar

Goldberg, J. E., and Richard, R. M. (1963). Analysis of nonlinear structures. J. Struct. Div. 89, 333–351. doi:10.1061/jsdeag.0000948

CrossRef Full Text | Google Scholar

Graesser, E., and Cozzarelli, F. (1991). Shape-memory alloys as new materials for aseismic isolation. J. Eng. Mech. 117, 2590–2608. doi:10.1061/(asce)0733-9399(1991)117:11(2590)

CrossRef Full Text | Google Scholar

Kikuchi, M., and Aiken, I. D. (1997). An analytical hysteresis model for elastomeric seismic isolation bearings. Earthq. Eng. Struct. Dyn. 26, 215–231. doi:10.1002/(sici)1096-9845(199702)26:2<215:aid-eqe640>3.0.co;2-9

CrossRef Full Text | Google Scholar

Kuczmann, M. (2010). Dynamic preisach hysteresis model. J. Adv. Res. Phys. 1.

Google Scholar

Losanno, D., Hadad, H., and Serino, G. (2019). Design charts for eurocode-based design of elastomeric seismic isolation systems. Soil Dyn. Earthq. Eng. 119, 488–498. doi:10.1016/j.soildyn.2017.12.017

CrossRef Full Text | Google Scholar

Losanno, D., Calabrese, A., Madera-Sierra, I., Spizzuoco, M., Marulanda, J., Thomson, P., et al. (2022a). Recycled versus natural-rubber fiber-reinforced bearings for base isolation: review of the experimental findings. J. Earthq. Eng. 26, 1921–1940. doi:10.1080/13632469.2020.1748764

CrossRef Full Text | Google Scholar

Losanno, D., De Domenico, D., and Madera-Sierra, I. (2022b). Experimental testing of full-scale fiber reinforced elastomeric isolators (freis) in unbounded configuration. Eng. Struct. 260, 114234. doi:10.1016/j.engstruct.2022.114234

CrossRef Full Text | Google Scholar

Losanno, D., Palumbo, F., Calabrese, A., Barrasso, T., and Vaiana, N. (2022c). Preliminary investigation of aging effects on recycled rubber fiber reinforced bearings (rr-frbs). J. Earthq. Eng. 26, 5407–5424. doi:10.1080/13632469.2021.1871683

CrossRef Full Text | Google Scholar

Masing, G., and Mauksch, W. (1926). “Über das verhalten von kalt gerecktem messing bei zug-und stauchbelastung,” in Wissenschaftliche Veröffentlichungen aus dem Siemens-Konzern (Springer), 142–155.

CrossRef Full Text | Google Scholar

Mayergoyz, I. (1986). Mathematical models of hysteresis. IEEE Trans. Magn. 22, 603–608. doi:10.1109/tmag.1986.1064347

PubMed Abstract | CrossRef Full Text | Google Scholar

Menegotto, M., and Pinto, P. (1973). “Method of analysis for cyclically loaded r. c. plane frames including changes in geometry and non-elastic behavior of elements under combined normal force and bending,” in Proc.of IABSE Symposium on Resistance and Ultimate Deformability of Structures Acted on by Well Defined Repeated Loads, Lisbona, 15–22.

Google Scholar

Ozdemir, H. (1976). Nonlinear transient dynamic analysis of yielding structures. Berkeley: University of California.

Google Scholar

Pellecchia, D., Lo Feudo, S., Vaiana, N., Dion, J.-L., and Rosati, L. (2022). A procedure to model and design elastomeric-based isolation systems for the seismic protection of rocking art objects. Comput. aided. Civ. Eng. 37, 1298–1315. doi:10.1111/mice.12775

CrossRef Full Text | Google Scholar

Preisach, F. (1935). Über die magnetische nachwirkung. Z. Phys. 94, 277–302. doi:10.1007/bf01349418

CrossRef Full Text | Google Scholar

Ramberg, W., and Osgood, W. R. (1943). Description of stress-strain curves by three parameters. Tech. rep.

Google Scholar

Sessa, S., Vaiana, N., Paradiso, M., and Rosati, L. (2020). An inverse identification strategy for the mechanical parameters of a phenomenological hysteretic constitutive model. Mech. Syst. Signal Process. 139, 106622. doi:10.1016/j.ymssp.2020.106622

CrossRef Full Text | Google Scholar

Vaiana, N., and Rosati, L. (2023). Classification and unified phenomenological modeling of complex uniaxial rate-independent hysteretic responses. Mech. Syst. Signal Process. 182, 109539. doi:10.1016/j.ymssp.2022.109539

CrossRef Full Text | Google Scholar

Vaiana, N., Sessa, S., Marmo, F., and Rosati, L. (2018). A class of uniaxial phenomenological models for simulating hysteretic phenomena in rate-independent mechanical systems and materials. Nonlinear Dyn. 93, 1647–1669. doi:10.1007/s11071-018-4282-2

CrossRef Full Text | Google Scholar

Vaiana, N., Sessa, S., Marmo, F., and Rosati, L. (2019a). An accurate and computationally efficient uniaxial phenomenological model for steel and fiber reinforced elastomeric bearings. Compos. Struct. 211, 196–212. doi:10.1016/j.compstruct.2018.12.017

CrossRef Full Text | Google Scholar

Vaiana, N., Sessa, S., Paradiso, M., Marmo, F., and Rosati, L. (2019b). “An efficient computational strategy for nonlinear time history analysis of seismically base-isolated structures,” in Conference of the Italian association of theoretical and applied Mechanics (Springer), 1340–1353. doi:10.1007/978-3-030-41057-5_108

CrossRef Full Text | Google Scholar

Vaiana, N., Marmo, F., Sessa, S., and Rosati, L. (2020). “Modeling of the hysteretic behavior of wire rope isolators using a novel rate-independent model,” in Nonlinear dynamics of structures, systems and devices (Springer), 309–317. doi:10.1007/978-3-030-34713-0_31

CrossRef Full Text | Google Scholar

Vaiana, N., Capuano, R., Sessa, S., Marmo, F., and Rosati, L. (2021a). Nonlinear dynamic analysis of seismically base-isolated structures by a novel opensees hysteretic material model. Appl. Sci. 11, 900. doi:10.3390/app11030900

CrossRef Full Text | Google Scholar

Vaiana, N., Sessa, S., and Rosati, L. (2021b). A generalized class of uniaxial rate-independent models for simulating asymmetric mechanical hysteresis phenomena. Mech. Syst. Signal Process. 146, 106984. doi:10.1016/j.ymssp.2020.106984

CrossRef Full Text | Google Scholar

Visintin, A. (2013). Differential models of hysteresis, 111. Berlin, Heidelberg: Springer Science & Business Media.

Google Scholar

Wen, Y.-K. (1976). Method for random vibration of hysteretic systems. J. Engrg. Mech. Div. 102, 249–263. doi:10.1061/jmcea3.0002106

CrossRef Full Text | Google Scholar

Wen, Y. (1980). Equivalent linearization for hysteretic systems under random excitation. J. Appl. Mech. 47, 150–154. doi:10.1115/1.3153594

CrossRef Full Text | Google Scholar

Keywords: phenomenological models, hysteresis, structural dynamics, structural mechanics, structural engineering

Citation: Capuano R, Fraddosio A and Piccioni MD (2022) Phenomenological rate-independent uniaxial hysteretic models: A mini-review. Front. Built Environ. 8:1048533. doi: 10.3389/fbuil.2022.1048533

Received: 19 September 2022; Accepted: 26 October 2022;
Published: 14 November 2022.

Edited by:

Yongjia Xu, Tsinghua University, China

Reviewed by:

Gabriele Milani, Politecnico di Milano, Italy
Biagio Carboni, Sapienza University of Rome, Italy
Lili Hu, Shanghai Jiao Tong University, China

Copyright © 2022 Capuano, Fraddosio and Piccioni. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Raffaele Capuano, raffaele.capuano@unina.it

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.