Modeling and parameter identification of asymmetric restoring force of reinforced concrete shear walls

Abstract

Under the action of the seismic waves, reinforced concrete shear walls (RC shear walls) often exhibit hysteresis characteristics such as strength and stiffness degradation, and pinching effects due to factors such as slippage, cumulative damage of concrete, and non-uniform corrosion of steel bars or reinforcements. The hysteresis curve is asymmetric. In this paper, a nonlinear restoring force mechanical model, based on the Bouc-Wen-Baber-Noori (BWBN) model, that can describe the hysteresis characteristics of the RC shear wall is obtained by introducing the asymmetry factor. Particular attention is the solution of equation that the hysteretic energy ε at different loading stages in this model is regarded as parameters to be identified. The ordinary differential equations for parameter identification are derived into four segments. The improved particle swarm organization (IPSO) algorithm is used for parameters identification, which superiority is verified by comparison to the genetic algorithm. The accuracy of this model is verified by using RC shear wall finite element simulation data.

Keywords

RC shear wall finite element simulation BWBN model hysteresis characteristics parameter identification

Introduction

RC shear walls have obvious asymmetric hysteresis characteristics due to variable axial force loading, reinforcement, non-uniform corrosion of steel bars, and other factors.¹ Affected by factors such as the slip between steel and concrete, concrete damage, RC shear walls always exhibit obvious hysteresis characteristics such as degradation in strength and stiffness, and pinching effects under the action of seismic waves.² Therefore, the establishment of an inelastic restoring force mechanical model that can effectively consider the typical hysteretic characteristics should be the basis for the seismic dynamic response analysis of RC shear walls.

The restoring force model of RC shear wall mainly includes the ideal elastoplastic model, the bilinear model, the clough model, and the Bouc-Wen-Baber-Noori (BWBN) model.³ Among these models, the BWBN model is widely used because it belongs to the smooth curve model, which can well describe the hysteretic characteristics such as strength and stiffness degradation, and pinching effect.⁴ Bouc^5,6 proposed a smooth hysteretic model for the analysis of the inelastic dynamic response of a single degree of freedom system. Wen⁷ improved the Bouc model by equivalent linearization technology, and proposed a smooth curve model, called the Bouc-Wen model for structural random vibration analysis. Subsequently, Baber and Wen⁸ and Baber and Noori⁹ defined the stiffness and strength degradation parameters through the cumulative hysteretic energy criterion, which improved the Bouc-Wen model and applied it to the analysis of the inelastic dynamic response of a single degree of freedom system. In addition, the pinch-up effect function was further introduced, resulting in a smooth curve model that can comprehensively consider the hysteresis characteristics of strength and stiffness degradation, and pinch-up effect. In the differential relationship between relative lateral displacement and hysteresis displacement, Kwok et al.¹⁰ corrected the velocity term by introducing adjustment coefficients, and proposed an asymmetric Bouc-Wen model suitable for hysteresis damper. Similarly, Zhu and Wang¹¹ corrected the differential relationship between relative lateral displacement and hysteretic displacement by introducing correction terms related to adjustment coefficients, displacement magnitude, and speed direction, and thus proposed an asymmetric Bouc-Wen model, which is suitable for describing the hysteresis characteristics of piezoelectric ceramics. Wang and Wen¹² corrected the differential relationship between the relative lateral displacement and the hysteretic displacement by introducing correction terms related to the adjustment coefficient, the direction of the hysteretic displacement and the direction of velocity, and proposed an asymmetrical Bouc-Wen model suitable for steel structures. Song and Der Kiureghian¹³ corrected the differential relationship between relative lateral displacement and hysteretic displacement by introducing six adjustment coefficients related to the direction of displacement, hysteretic displacement and speed, and proposed an asymmetric Bouc-Wen model which could depict the hysteretic characteristic of flexible connector for substations. Pain et al.¹⁴ proposed a calibration scheme for a generalized Bouc-Wen model to capture the asymmetric behavior of soil at the high cyclic shear strain. The degradation parameters and pore–water pressure are expressed as a function of dissipated hysteretic energy. Besides, this model is also applied to structures that usually occur under significant cyclic excitation (e.g. earthquake, wind, or sea wave).^15–17 For example, Lukassen et al.¹⁸ proposed a method to calculate deformation characteristics and hysteresis effects of tensile armored steel wire in marine flexible pipes under tensile and bending loads, and studied the effects of cyclic bending loads on laying angles and friction coefficients of different tensile armored layers. The results show that the steel wire will gradually move to the compression side of the bending surface under the combined action of tension and cyclic bending. Although the above model can consider the influence of asymmetric characteristics to a certain extent, it cannot involve the influence of accumulated damage on the degradation effect. There is a need for a restoring force model of RC shear walls under the action of seismic waves that can effectively consider the strength and stiffness degradation, pinching effect, and asymmetric hysteresis characteristics.

Regarding the parameter identification of the model, because of the high nonlinearity involved in the model, parameter identification of hysteretic model represented by Bouc-Wen model is a difficult task. Up to now, the deterministic identification techniques used for nonlinear applications include least squares estimation (LSE),^19–21 extended Kalman filter (EKF),²¹ Bayes algorithm,²² etc. Compared with the above deterministic algorithm, the stochastic algorithm usually converges to the global minimum more easily. By defining the fitness function of the hysteresis model, Zhang and Xia²³ and Zhao et al.²⁴ proposed an improved particle swarm optimization algorithm for parameter identification of nonlinear dynamic hysteresis model, which improved the global search ability and avoided the precocious convergence of traditional particle swarm optimization algorithms. The Leishman-Beddoe dynamic stall model of rotor blade and the elastic displacement field model of elastomer damper are applied. Pain et al.¹⁴ expressed degradation parameters and pore water pressure as functions of dissipated energy, and used particle swarm optimization algorithm to identify parameters of the asymmetric Bouc-Wen model of sandstone. The adjusted model can successfully describe the stiffness degradation, strength degradation, excess pore water pressure, and asymmetric behavior of saturated sandstone and siltstone under high cyclic shear strain. Manikantan et al.²⁵ proposed a non-random algorithm for identifying parameters of Bouc-Wen model, namely homotopy algorithm, in order to solve the shortcomings of random methods such as simulated annealing and genetic algorithm, which have high computational cost and are not easy to converge to the global minimum. Numerical results verified the advantages of homotopy method in terms of computation amount and convergence efficiency.

In this paper, an inelastic restoring force model was established based on the BWBN model, which is effective is considering strength and stiffness degradation, pinching effects, and asymmetric hysteresis characteristics. This was done by introducing shape control parameters describing the asymmetric hysteresis properties in the first-order ordinary differential equations of hysteretic displacement and lateral displacement. The improved particle swarm optimization (IPSO) and genetic algorithm (GA) were used to identify the parameters of the model and compare the performance of the algorithms. Concurrently, the finite element simulation data of the pseudo-static loading of the RC shear wall is used to verify the accuracy of the inelastic restoring force model introduced this paper.

Asymmetric BWBN model of RC shear wall

Affected by factors such as concrete damage, cracking, and steel slippage, RC shear walls usually exhibit obvious hysteretic characteristics such as strength degradation, stiffness degradation, and pinching effects under strong earthquakes.^1,2 Based on the BWBN model,²⁶ the inelastic restoring force of the RC column can be described as:

R (x, z) = α k_{0} x + (1 - α) k_{0} z

(1)

where R(x, z) is the inelastic restoring force; k₀ is the initial linear elastic stiffness; α is the ratio of stiffness after yielding to stiffness before yielding; x is the lateral displacement, and; z is the hysteretic displacement.

Due to variable axial force loading, asymmetric cross-section or reinforcement, and uneven corrosion of steel bars, the hysteresis curve of the RC shear wall shows obvious asymmetry. The aforementioned restoring force-displacement hysteresis curve is symmetrical about the origin. Therefore, it is difficult to describe the asymmetric hysteresis characteristics of RC shear walls. In view of this, a shape control parameter, τ, describing the asymmetric characteristics is introduced into the traditional BWBN model, and the relationship of derivative of time between the relative lateral displacement, x, and the hysteretic displacement, z, is:

\overset{\cdot}{z} = h (z) \overset{\cdot}{x} \frac{1 - ν [β sign (z \overset{\cdot}{x}) + τ ξ + γ] {| z |}^{n}}{η}

(2)

ξ = sign (z) + sign (\overset{\cdot}{x})

(3)

η = 1 + δ_{η} ε; ν = 1 + δ_{ν} ε

(4)

h (z) = 1 - ζ_{1} \exp [- {(\frac{z \cdot sign (\overset{\cdot}{x})}{ζ_{2}})}^{2}]

(5)

ζ_{1} = ζ_{s} [1 - \exp (- p ε)]

(6)

ζ_{2} = (ψ + δ_{ψ} ε) (λ + ζ_{1})

(7)

z_{max} = {[ν (β + γ \pm 2 τ)]}^{- \frac{1}{n}}

(8)

where, β, γ, and n are the shape parameters of the hysteresis curve; η and ν are the stiffness and strength degradation parameters respectively; δ_η and δ_ν are the stiffness and strength degradation rate respectively; sign( $\overset{\cdot}{x}$ ) is the sign function related to speed $\overset{\cdot}{x}$ , when $\overset{\cdot}{x}$ > 0, $\overset{\cdot}{x}$ = 0, and $\overset{\cdot}{x}$ < 0, sign( $\overset{\cdot}{x}$ ) takes 1, 0, −1 respectively; h(z) is the pinching effect function, where ζ_s controls the total slip, ψ controls the pinch amplitude, δ_ψ is the pinch rate, λ_s the interaction coefficient between the pinch severity and the pinch rate, q controls the starting position of the pinch, p is the pinch slope. Furthermore, ε is the cumulative hysteretic energy consumption in a loading cycle (T):

ε = (1 - α) k_{0} \int_{0}^{T} \overset{\cdot}{x} z d t

(9)

By dividing both sides of equation (2) by $\overset{\cdot}{x}$ , the partial derivative of z to x can be obtained:

\frac{d z}{d x} = h (z) \frac{1 - ν [β sign (z \overset{\cdot}{x}) + τ ξ + γ] {| z |}^{n}}{η}

(10)

Parameter identification of the asymmetric BWBN model

Objective function

To identify the parameters of the asymmetric BWBN restoring force model, the goal is to find a set of model parameter vectors $\hat{Q}$ substituted into the asymmetric BWBN model to minimize the error between the theoretical restoring force R(t| $\hat{Q}$ ) and the restoring force R(t) of the finite element simulation. Therefore, the parameter identification of the inelastic restoring force model of the RC column can be described as a nonlinear optimization problem:

{\begin{matrix} min f (\hat{Q}) \\ s . t . {\hat{Q}}_{L} \leq \hat{Q} \leq {\hat{Q}}_{U} \end{matrix}

(11)

where, ${\hat{Q}}_{L}$ and ${\hat{Q}}_{U}$ are the upper and lower bounds of the model parameter vector $\hat{Q}$ , respectively; f( $\hat{Q}$ ) is the objective function or fitness function, which reflects the error between the measured value and the calculated value, defined as:

f (\hat{Q}) = \sum_{i = 1}^{N} [R (t_{i}) - R (t_{i} | \hat{Q})]

(12)

where, N is the number of loading steps; R(t_i) is the finite element simulation data of the RC shear wall at step t_i; R(t_i| $\hat{Q}$ ) is the restoring force calculated by substituting the model parameters $\hat{Q}$ into the BWBN model for the t_ith step. According to the difference of loading and unloading, formulas (2) through (8) can be divided into four parts:

{\begin{matrix} \begin{matrix} d z_{i} = h (z_{i}) \frac{d x_{i} - [1 + ε_{i} δ_{v}] [β sign (z_{i} {\overset{\cdot}{x}}_{i}) + τ ξ_{i} + γ] {| z_{i} |}^{n} d x_{i}}{1 + ε_{i} δ_{η}} \\ h (z_{i}) = 1 - ζ_{1} \exp [- \frac{z_{i} \cdot sign ({\overset{\cdot}{x}}_{i}) - q z_{max}}{(ψ + δ_{ψ} ε_{i}) (λ + ζ_{1})}] \\ ζ_{1} = ζ_{s} [1 - \exp (- p ε_{i})] i = 1, 2, 3, 4 \\ z_{max} = {[(1 + \sum_{i}^{4} ε_{i} δ_{v}) (β + γ)]}^{- \frac{1}{n}} \\ ξ_{i} = sign (z_{i}) + sign ({\overset{\cdot}{x}}_{i}) \\ R = α k_{0} x + (1 - α) k_{0} z \end{matrix} \end{matrix}

(13)

When $z > 0, \overset{\cdot}{u} > 0$ ,

d z_{1} = h (z_{1}) \frac{(d x_{1} - [1 + ε_{1} δ_{v}] (β + 2 τ + γ) {z_{1}}^{n} d x_{1})}{1 + ε_{1} δ_{η}}

(14)

When $z > 0, \overset{\cdot}{u} < 0$ ,

d z = h (z_{2}) \frac{(d x_{2} - [1 + ε_{2} δ_{v}] (- β + γ) {z_{2}}^{n} d x_{2})}{1 + ε_{2} δ_{η}}

(15)

When $z < 0, \overset{\cdot}{u} > 0$ ,

d z_{3} = h (z_{3}) \frac{(d x_{3} - [1 + ε_{3} δ_{v}] (- β + γ) {(- 1)}^{n} {z_{3}}^{n} d x_{3})}{1 + ε_{3} δ_{η}}

(16)

When $z < 0, \overset{\cdot}{u} < 0$ ,

d z_{4} = h (z_{4}) \frac{(d x_{4} - [1 + ε_{4} δ_{v}] (β - 2 τ + γ) {(- 1)}^{n} {z_{4}}^{n} d x_{4})}{1 + ε_{4} δ_{η}}

(17)

To facilitate the solution, the hysteretic energy consumption of each stage is set as the parameter to be identified, namely ε_i (i = 1, 2, 3, 4). Therefore, the original asymmetric BWBN model is transformed into a system of ordinary differential equations without an integral. The fitness function can be solved numerically using the fourth-order Runge-Kutta method.

The asymmetric BWBN model has 18 model parameters, including five shape parameters (k₀, α, β, γ, τ, and n), two degradation parameters (δ_η and δ_ν), six pinching effect parameters (ζ_s, p, q, ψ, δψ, and λ), and four hysteretic energy dissipation parameters, ε_i (i = 1, 2, 3, 4). Although these model parameters have clear physical meanings, they lack quantitative calculation formulas except for hysteretic energy consumption.²⁷ Therefore, it is necessary to identify the parameters of the asymmetric BWBN model based on the pseudo-static reciprocating load that comes from finite element simulation of the RC shear wall.

Improved particle swarm optimization algorithm

Particle Swarm Optimization (PSO) is a new evolutionary algorithm. The algorithm starts from a random solution and changes the historical movement state information of itself and the group through mutual cooperation and information sharing between particles, so as to achieve the purpose of gradually increasing the chance of discovery of better solutions during the solution process. PSO retains the same global search strategy and parallel processing of population as the genetic algorithm, but adopts a velocity position model that is simple to operate, which avoids crossover and mutation operations. The particles follow the optimal particles in the solution space to search the global optimal solution. This algorithm has the characteristics of easy implementation, high precision, and rapid convergence, which has been widely used in application fields such as function optimization, system identification, and fuzzy control.^23,28–31

During the flight of the particle, the velocity v_i determines the direction and position of the particle. The speed of the particle is continuously updated by the two optimal values of the local optimal solution: P_i found by the particle itself, and the global optimal solution P_g found by the entire population. The updated formula of particle velocity and position is as follows:

\begin{matrix} v_{i} (t + 1) = ω \cdot v_{i} (t) + c \cdot ran d_{1} (p_{i} (t) - x_{i} (t)) \\ + c_{2} \cdot ran d_{2} (p_{g} (t) - x_{i} (t)) \end{matrix}

(18)

x_{i} (t + 1) = x_{i} (t) + v_{i} (t + 1)

(19)

where, v_i is the velocity of the particle; ω is the inertia weighting factor; P_i is the optimal value of the individual; P_g is the optimal value of the group; x_i is the current particle position; rand is a random function, between (0, 1), and; c₁and c₂ are learning factors.

The global convergence of the algorithm is affected by the flight speed of a single particle. The function ω·v_i(t) can be understood as the previous velocity or inertia of the particle. The function c₁·rand₁(p_i(t)-x_i(t)) can be explained as the cognitive behavior of particles, representing the information sharing and mutual cooperation between particles. Throughout the solution process, the inertia weight ω, acceleration factors c₁ and c₂, and the maximum velocity v_max jointly maintain the balance of the particle’s global and local search capabilities.

Due to the large number of parameters to be identified, the solution space is complicated. In order to improve calculation efficiency, the standard particle swarm algorithm is improved by the method of linearly decreasing the inertia weight factor (IPSO). The linearly decreasing weight method can be expressed as²³:

w = w_{max} - (w_{max} - w_{min}) \frac{run}{ru n_{max}}

(20)

where, w_max is the maximum inertia weight; w_min is the minimum inertia weight; run is the current iteration number; and run_max is the total number of algorithm iterations.

A larger w has better global convergence ability, and a smaller w has stronger local convergence ability. Therefore, as the number of iterations increases, the inertia weight w should continue to decrease, so that the particle swarm algorithm has a strong global search ability in the early stage, and a strong local convergence ability in the late stage.

The steps of the algorithm are as follows:

Step 1: Initialize the population size n, the maximum number of iterations t_max, the maximum velocity v_max, the inertia weight ω, the initial particle position x, the velocity v, and other parameters. Set the initial position as the optimal position P_i experienced by each particle in the current situation, and the global optimal position of the current particle swarm as the initial P_g.

Step 2: Calculate the fitness value of each particle. Store and compare the usage value corresponding to the optimal position of each particle. Select the position of the particle with the best fitness in the entire population as the optimal position P_g of the current population.

Step 3: After the calculation of the inertia weight factor is completed in a linearly decreasing mode, the position and speed of the particles in the current population are updated according to the updated formula of the particle speed and position.

Step 4: Calculate the fitness of each particle after the position is updated and compare the fitness of the updated particle with the fitness of all the best positions that are previously updated. If the updated fitness is better, update the current particle position as the optimal position P_i of the particle.

Step 5: The fitness of each particle after the update is compared with the fitness corresponding to the best position P_g experienced by all particles. If the fitness of the current particle position is better, the value of P_g is updated.

Step 6: Judge whether the termination condition is met. If the termination condition has been met, the iterative calculation is stopped and the optimal solution is output. If the termination conditions are not met, return to Step 3 to continue the iterative calculation.

The IPSO algorithm process is shown in the Figure 1.

Figure 1.

Improved PSO algorithm.

Finite element simulation

Many scholars obtain experimental data through pseudo-static reciprocating loading tests of RC columns or RC shear walls.^32,33 However, the cost of testing is relatively high. Based on testing and numerical calculation techniques, a numerical calculation model for RC shear walls can be established,^34–36 as shown in Figure 2. The wall size is 1400 mm × 160 mm × 2800 mm, with 200 mm × 160 mm × 2800 mm post-cast concealed columns on both sides of the wall. The compressive strength of the concrete cube is 43.1 MPa, the density is 2500 kg/m³, the modulus of elasticity is 3.15 × 10⁴ MPa, and the Poisson’s ratio is 0.2. The concealed column is provided by 4 HRB335 longitudinal steel bars, with a diameter of 16 mm, and HPB300 stirrups, with a diameter of 8 mm and a spacing of 150 mm. The HPB300 steel bars were installed with a spacing of 200 mm and a diameter of 8 mm in the horizontal and vertical directions of the wall, respectively. The density of the steel bar is 7.85g/cm³, with a modulus of elasticity of 2 × 10⁵ MPa, and a Poisson’s ratio of 0.2. The plastic mechanical properties of the steel bars used in this study are presented in Table 1.

Figure 2.

Reinforced concrete shear wall.

Table 1.

The plastic mechanical properties of the steel bars used in this study.

Rebar grade	Yield stress (MPa)	Ultimate stress (MPa)
HPB300	398.14	475.27
HPB335	503.9	603

The load axial compression of concrete ratio is 0.15. Horizontal load is applied reciprocally, using force-displacement mixed-control loading. The first five levels (50, 100, 150, 200, 300 kN) are load-controlled and each level is cyclically loaded once. This is followed by displacement-controlled loading, with magnitudes of 8, 16, 24, 32, 40, 50, and 60 mm, of which each 16, 32, and 40 mm load cycle is cycled twice, for a total of 11 cycles. The force-controlled loading is equivalently converted into displacement-controlled loading, for a total of 16 loading cycles, as shown in the Figure 3.

Figure 3.

Loading method.

Due to the non-linear contact between the steel bar and the concrete, the C3D8R solid element is used for the shear wall and the T3D2 truss element is used for the steel bar. The plastic damage model is adopted to the concrete. The relationship between compression and tensile yield stress and strain is shown in Figure 4. The mechanical properties of the steel bars adopt an ideal elastoplastic model.

Figure 4.

Plastic damage model of concrete.

When the loading displacement is too large, the concrete and steel bars in the shear wall will slip. If the contact relationship is set, the calculation amount will be difficult to converge. Therefore, the ABAQUS subroutine should be used to set up the clough rebar model which depicts the hysteresis characteristic of RC walls and so on to weaken its unloading stiffness and achieve the effect of sliding equivalent. Using the static implicit algorithm, the hysteretic curve of the RC shear wall is calculated and shown in the Figure 5.

Figure 5.

Hysteresis curve of RC shear wall.

Verification and discussion

Algorithm optimization

In order to verify the superiority of the improved PSO algorithm, the genetic algorithm^37,38 was used for comparison. The data of the first and last loading cycles in the finite element simulation were selected for calculation. The population numbers of GA and IPSO algorithm were both set to 180. According to previous experience of several trial calculations, the parameter range setting is shown in Table 2.

Table 2.

Range of parameters to be identified.

Parameter	k ₀	ε ₁	ε ₂	ε ₃	ε ₄	δ _v	α	β	γ	n	δ
Range	(100, 1900)	(10³, 10⁴)	(10³, 10⁴)	(10³, 10⁴)	(10³, 10⁴)	(0, 0.1)	(0, 0.5)	(−50, 50)	(−10, 10)	(1, 5)	(−10, 10)

Parameter	ζ _s	ψ	δ _ψ	λ	q	p	δ _η
Range	(0, 1)	(1, 10)	(0, 1)	(1, 10)	(1, 10)	(1, 10)	(−0.1, 0.1)

The fitness function values after iterative convergence are presented in Table 3. It can be seen from the table that whether it is the first load cycle or the last load cycle. The fitness function value of IPSO algorithm is smaller than that of GA, meaning the IPSO algorithm obtains higher accuracy.

Table 3.

Fitness function values.

	PSO	GA
The first cycle loading	69.587	396.732
The last cycle loading	2986.7829	7890.7682

Figure 6 shows the change rule of the fitness function value and the number of iterations for the first (a) and last (b) loading cycles. It can be seen from Figure 6(a) and (b) that the fitness function value of IPSO algorithm is much larger than the fitness value of GA at a high probability at the beginning. This is because IPSO algorithm starts from a random solution, and GA requires an initial solution. As a result, when convergence is reached, the number of iterations of IPSO algorithm is greater than that of GA. However, the IPSO algorithm does not require crossover and mutation operations, so the total iteration time of IPSO algorithm (289 s) is less than that of GA (572 s).

Figure 6.

Convergence process of the fitness value of PSO and GA.

Table 4 and Figure 7 show the convergence value of each identified parameter and the iteration process under the last loading cycle with IPSO algorithm, respectively. It should be noted that k₀, ε₁, ε₂, β, γ, and δ_ψ are all at the boundary of the corresponding parameter range, indicating that there is a certain error when the parameter boundary is initially given. Therefore, the optimal parameter range can be determined through multiple trial calculations to further improve the calculation accuracy.

Table 4.

Parameter values at convergence.

Hysteretic parameters	Value
k ₀	100
ε ₁	1000
ε ₂	10,000
ε ₃	5716.520
ε ₄	9604.658
δ _v	4.4979 × 10⁻⁵
α	0.437
β	50
γ	10
n	1
τ	3.7376
ζ _s	0.0202
ψ	9.998
δ _ψ	0
λ	4.112
q	9.998
p	0
δ _η	0.04

Figure 7.

Convergence process of each parameter value of PSO.

Model accuracy verification

Based on the finite element simulation data, the accuracy of the asymmetric restoring force model in this paper was verified. As shown in Figure 8, the comparison of hysteresis loops between simulation and asymmetric restoring force model in different loading circulations. It can be seen from Figure 8(a) and (b) that the restoring force model described in this article coincides better with the finite element simulation data. In particular, in the unloading stage, the hysteresis curve based on the PSO algorithm almost coincides with the hysteresis curve of the finite element simulation, indicating that PSO algorithm can more effectively identify the parameters describing the pinch effect and curve asymmetry.

Figure 8.

Comparison of hysteresis loops between simulation and asymmetric restoring force model in different loading circulations.

Conclusions

(1) Under the action of the seismic waves, reinforced concrete shear walls (RC shear walls) often exhibit hysteresis characteristics such as strength and stiffness degradation, and pinching effects due to factors such as slippage, cumulative damage of concrete, and non-uniform corrosion of steel bars or reinforcements. Based on the Bouc-Wen-Baber-Noori (BWBN) model, by introducing the asymmetry factor, the asymmetric restoring force model is described in this article to depict the hysteresis characteristics of RC walls. In particular, we propose a different solution, that is, the hysteretic energy ε in this model is regarded as a parameter to be identified, which simplifies the solution of differential equations and solving the equation by splitting it into four sections.

(2) The method of decreasing inertia weight is used to improve the PSO algorithm (IPSO) which is used for the parameter identification of the asymmetric restoring force model. By comparing with the finite element simulation data, the model can describe the hysteresis characteristics of the RC shear wall under lateral load. At the same time, the GA is compared to verify the superiority of the IPSO algorithm in terms of calculation accuracy and speed. However, the parameters obtained by IPSO algorithm are more accurate in describing the hysteretic characteristics of RC shear walls.

Footnotes

Handling Editor: Chenhui Liang

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The authors thank for the project support of “(Sichuan university-Luzhou) Innovation Mode, Key Technologies and Application Demonstration of High-performance Intelligent Electro-hydraulic Products.”

ORCID iD

Zhou Lina

References

Foliente

. Hysteresis modeling of wood joints and structural systems. J Struct Eng 1995; 121: 1013–1022.

Ajavakom

. On system identification and response prediction of degrading structures. Struct Control Health Monit 2006; 13: 347–364.

Sengupta

. Hysteresis behavior of reinforced concrete walls. J Struct Eng 2014; 140: 1–18.

Sengupta

. Modified Bouc–Wen model for hysteresis behavior of RC beam–column joints with limited transverse reinforcement. Eng Struct 2013; 46: 392–406.

Bouc

. Forced vibration of mechanical systems with hysteresis. In: Proceedings of the fourth conference on nonlinear oscillation, Prague, Czechoslovakia, 5–9 September 1967, pp.315–321.

Bouc

. A mathematical model for hysteresis. Acta Acust United Acust 1971; 24: 16–25. (in French)

Wen

. Method for random vibration of hysteretic systems. J Eng Mech 1976; 102: 249–263.

Baber

Wen

. Random vibration of hysteretic, degrading systems. J Eng Mech Div 1981; 107: 1069–1087.

Baber

Noori

. Random vibration of degrading pinching systems. J Eng Mech 1985; 111: 1010–1026.

10.

Kwok

Nguyen

, et al. Bouc–Wen model parameter identification for a MR fluid damper using computationally efficient GA. ISA Trans 2007; 46: 167–179.

11.

Zhu

Wang

. Non-symmetrical Bouc–Wen model for piezoelectric ceramic actuators. Sens Actuators A Phys 2012; 181: 51–60.

12.

Wang

Wen

. Evaluation of pre-Northridge low-rise steel buildings. II: reliability. J Struct Eng 2000; 126: 1169–1176.

13.

Song

Der Kiureghian

. Generalized Bouc–Wen model for highly asymmetric hysteresis. J Eng Mech 2006; 132: 610–618.

14.

Pain

Nimbalkar

Hussain

. Applicability of Bouc-Wen model to capture asymmetric behavior of sand at high cyclic shear strain. Int J Geomech 2020; 20: 1–10.

15.

Gong

Bao

, et al. Asymmetric buckling of offshore pipelines under combined tension, bending and external pressure. Ships Offshore Struct 2015; 10: 162–175.

16.

Chatzi

Smyth

Masri

. Experimental application of on-line parametric identification for nonlinear hysteretic systems with model uncertainty. Struct Saf 2010; 32: 326–337.

17.

Luo

Zhou

, et al. Seismic performance of lead-filled steel tube damper: laboratory test, parameter identification and application. Eng Struct 2020; 219: 26–39.

18.

Lukassen

Krenk

Glejbøl

, et al. Non-symmetric cyclic bending of helical wires in flexible pipes. Appl Ocean Res 2019; 92: 101876.

19.

Smyth

Masri

Kosmatopoulos

, et al. Development of adaptive modeling techniques for non-linear hysteretic systems. Int J Non Linear Mech 2002; 37: 1435–1451.

20.

Wang

Ding

Ren

, et al. Structural dynamic nonlinear model and parameter identification based on the stiffness and damping marginal curves. Struct Control Health Monit 2020; 27: 56–68.

21.

Mariani

Corigliano

. Impact induced composite delamination: state and parameter identification via joint and dual extended Kalman filters. Comput Methods Appl Mech Eng 2005; 194: 5242–5272.

22.

Worden

Becker

. On the identification of hysteretic systems. Part II: Bayesian sensitivity analysis and parameter confidence. Mech Syst Signal Process 2012; 29: 213–227.

23.

Zhang

Xia

. An improved PSO algorithm for parameter identification of nonlinear dynamic hysteretic models. J Sound Vib 2017; 389: 153–167.

24.

Zhao

Zhou

. Dynamic modeling of hysteresis nonlinearity based on improved Particle Swarm Optimization algorithm. In: Proceedings of the 31st Chinese control conference, Hefei, China, 25–27 July 2012. New York: IEEE Computer Society.

25.

Manikantan

Ghosh Mondal

Suriya Prakash

, et al. Parameter identification of Bouc–Wen type hysteresis models using homotopy optimization. Mech Based Des Struct Mach 2022; 50: 26–47.

26.

Goda

Hong

Lee

. Probabilistic characteristics of seismic ductility demand of SDOF systems with Bouc-Wen hysteretic behavior. J Earthq Eng 2009; 13: 600–622.

27.

Kunnath

Mander

Fang

. Parameter identification for degrading and pinched hysteretic structural concrete systems. Eng Struct 1997; 19: 224–232.

28.

Quaranta

Marano

Greco

, et al. Parametric identification of seismic isolators using differential evolution and particle swarm optimization. Appl Soft Comput 2014; 22: 458–464.

29.

Tang

Xue

Fan

. Differential evolution strategy for structural system identification. Comput Struct 2008; 86: 2004–2012.

30.

Lei

Zhou

Tan

, et al. Optimal and secure audio watermarking scheme based on self-adaptive particle swarm optimization and quaternion wavelet transform. Signal Process 2015; 113: 80–94.

31.

Guan

Dai

Xia

. Pretension optimization and verification test of double-ring deployable cable net antenna based on improved PSO. Aerosp Sci Technol 2014; 32: 19–25.

32.

Storn

Price

. Differential evolution - a simple and efficient adaptive scheme for global optimization over continuous spaces. J Glob Optim 1997; 11: 341–359.

33.

Atalay

Penzien

. The seismic behavior of critical regions of reinforced concrete components as influenced by moment, shear and axial force. Report no. EERC 75-19, University of California, Berkeley, 1975, pp.112–113.

34.

Sengupta

. Hysteresis modeling of reinforced concrete structures: state of the art. ACI Struct J 2017; 114: 25–38.

35.

Lan

Xiao

, et al. Macro-modeling of reinforced concrete structural walls: state-of-the-art. J Earthq Eng 2017; 21: 652–678.

36.

Ouzaa

Benmansour

. Cracks in base-restrained plain and reinforced concrete walls. Turk J Eng Environ Sci 2010; 34: 215–230.

37.

Marano

Quaranta

Monti

. Modified genetic algorithm for the dynamic identification of structural systems using incomplete measurements. Comput Aided Civ Infrastruct Eng 2011; 26: 92–110.

38.

Holland

. Adaptation in natural and artificial systems: an introductory analysis with applications to biology, control, and artificial intelligence. Ann Arbor, MI: University of Michigan Press, 1975.