Abstract
This article presents an analysis procedure for evaluation of load–deformation behavior of reinforced concrete shear walls with continuous or lap-spliced bar connections in plastic hinge zones under horizontal loads. For the shear walls with continuous bars, the lateral deformations caused by flexure, shear, and reinforcement slip are evaluated by considering their interaction. The flexural deformation is calculated by conventional fiber model. The shear mechanism is based on modified compression field theory with a softened smeared cracked reinforced concrete membrane element. Both the flexural and shear deformations are estimated separately in the plastic hinge and non-plastic hinge regions. In addition, an approach is proposed for analysis of plastic hinge length based on fracture energies of materials. For the shear walls with lap-spliced bars, due to its complicated behavior and mechanism, a simple way to deal with the lap splice is proposed. The equations regarding bond-slip of the lap splice with minimum spliced length are established and the stress and strain states of lap splices with different spliced lengths are analyzed on the basis of equilibrium of forces with a mean bond stress model. Finally, the validity of the proposed analysis procedure is confirmed by comparing the analytical results with previous experimental data.
Introduction
Due to the high lateral stiffness of shear walls to resist wind and earthquake action, it is effective to use shear walls for buildings as the first lateral load resisting system. There are many experimental studies performed to investigate the seismic performance and deformational characteristics of reinforced concrete (RC) shear walls (Almeida et al., 2016; Dazio et al., 2009; Kang et al., 2013; Syed, 2013; Zhang et al., 2011). The shear walls usually provide a good response under seismic action due to their dimensions and in-plane position arrangement. The axial, flexural, shear, strain penetration of the longitudinal reinforcing bars into the anchorage and sliding shear mechanisms are often used to explain the deformation behavior of shear walls with different detailing and aspect ratios under axial and horizontal loads, especially for the combined flexural and shear behavior. Shear walls with different sections, reinforcement layouts, and aspect ratios lead to different failure modes, namely, bending failure and shear failure. To analyze a flexure-dominant shear wall (slender member), many fiber beam-column elements have been developed in the past several decades. The fiber method is mainly based on the traditional Euler–Bernoulli beam theory that plane sections remain plane. However, it may not be appropriate when axial, bending, and shear actions are coupled in RC members with low or middle aspect ratios.
Many studies have attempted to cope with the axial-bending-shear coupling by adopting the Timoshenko beam theory or a series of uniaxial concrete elements in a so-called multiple vertical line element model (MVLEM; Orakcal et al., 2004). Other models have been used to determine the shear behavior of concrete members by developing rotating (or fixed) angle softened truss model (RA-STM or FA-STM; Hsu, 1988; Hsu and Zhang, 1996), distributed stress field model (Vecchio, 2000), and axial-shear-flexure interaction (ASFI) method (Mostafaei and Vecchio; Mostafaei and Kabeyasawa, 2007). The developed ASFI method was derived through both equilibrium and compatibility of flexure and shear mechanisms that employed the plane section assumption and the modified compression field theory (MCFT; Vecchio and Collins, 1986). However, this method was mainly employed in RC columns with an inflection point, and the plastic hinge length
where
where
For fixing the reinforcement conveniently, a lap-spliced joint can be used in the concrete members. Biskinis and Fardis (2010a; 2010b) developed ultimate and yield deformation models of flexure-controlled members with continuous bars and lap-spliced bars, respectively, in the plastic hinge region. The models were mainly derived with empirical expressions according to a large databank of cyclic and monotonic tests. Tastani et al. (2015) established a series of field equations governing lap-spliced behavior. The presented model is mainly used when the bars are elastic (namely,
Shear walls with continuous bars at the connection
For estimating the response of a typical RC shear wall, the total displacement
The components of the total response are represented by several springs and are illustrated in Figure 1. At the pre-yield stage, three springs in series are used to model the three components of lateral response. At the post-yield stage, the shear wall is divided into two sections that are plastic hinge region and elastic region (non-plastic hinge region). And five springs in series located at characteristic positions along the two sections are used to model the three components of lateral response. The shear and flexural deformations are calculated by the MCFT and fiber models. Their combined interaction in the analytical procedure is based on the concept of the ASFI approach (Mostafaei and Kabeyasawa, 2007). For the flexural deformation analysis, the inclined cracked concrete in the web of shear wall will change the concrete strength which is considered by employing compression softening factor based on MCFT for an element within or outside the plastic hinge region. For the shear deformation analysis, both equilibriums of axial stresses caused by flexure, shear and applied axial load, and compatibility of strains will be taken into account. Then, the flexural and shear behavior analyses can be performed separately in a simple way.

Components of the total deformation response represented by springs: (a) a test specimen in the previous experiment, (b) a flexure-type deformation of shear wall, (c) components of the total response at pre-yield stage with three springs, and (d) components of the total response at the post-yield stage with five springs.
Flexural deformation
In the fiber model analysis, flexural deformations are evaluated through integrating the curvature along the height of cantilever shear wall. Generally, the shear wall curvature distribution is idealized and illustrated in Figure 2. Before yielding of reinforcement, the curvature shows a linear and elastic curvature distribution. After yielding, the curvature shows a nonlinear curvature distribution, which assumes that the regions of inelastic curvature are concentrated at the end of shear wall with a fixed length. The yield curvature can be calculated as follows
where

Section stress and strain distributions of concrete and reinforcement and curvature distribution assumed along the shear wall axis.
After yielding, the lateral load is evaluated as
where
where

Constitutive laws of the materials: (a) confined and unconfined concrete including tensile and compressive constitutive laws and (b) reinforcement with a bilinear stress–strain relationship.
The stress–strain relationship of compressed concrete according to the pre-peak branch of a Hognestad (1951) parabola
where
The stress–strain relationship of concrete in tension is given by
where
After yielding, the flexural deformation can be calculated as follows
where H is the height of shear wall,
The plastic hinge length is considered to be constant after yielding of reinforcement. In fact, this is an approximate method, because the plastic hinge length usually increases with an increase of applied load at early stage, and then tends to constant when the reinforcement strains are large and are developed in a certain concentrated length (Dazio et al, 2009). There are several reasons to adopt the constant plastic length for this study. First, constant length plastic hinge models are currently available in commercial and research software and, thus, can be used immediately by practicing and research engineers (Berry et al., 2008). Second, according to the vertical strains along the height of the shear walls in relevant literatures (Dazio et al., 2009; Tran and Wallace, 2015), the major deformations were often lumped in a certain region at the plastic stage. Although the plastic region will spread along the height of the specimens, the value of the increasing plastic strain at the top of plastic hinge region was not large. Third, the subsequent section of this study will show that the fixed length plastic hinge model can be used to predict shear wall deformation behavior accurately for a range of demand levels, which is the primary objective. Various well-known empirical models for
In this study, it is assumed that the energy of the applied lateral load in the member, which is produced by flexural behavior when the member is in the post-yield stage, is equal to the sum of fracture energies of materials including reinforcement and concrete in the plastic hinge region. In other words, the energy from the external applied load producing flexural deformation is dissipated by the materials of plastic hinge region. This assumption is reasonable because most of the external energy is dissipated by flexural behavior of member (Eom and Park, 2013), and the flexural deformations are mainly attributed to the plastic hinge region at the post-yield stage
where

Energy dissipation in the plastic hinge region of a member: (a) a plastic hinge energy dissipation mechanism, (b) fracture energy of reinforcement, and (c) fracture energy of concrete in compression.
The integral result of equation (11) is the area of post-yield stress–strain curve multiplied by
where
The result of equation (13) is the area of the post-peak concrete stress–strain curve multiplied by
where
where
Then, an upper-bound evaluation (α = 0.5 and β = 0.1) resulted from Joint ACI-ASCE Committee 428 (1968) recommendations and a lower-bound evaluation (α = 0.25 and β = 0.044) resulted from Paulay and Priestley’s recommendations (Paulay and Priestley, 1993) are adopted in this study
Shear deformation
After selection and calculation of the constitutive models and flexural deformations, respectively, the next steps are estimation of shear deformation and implementation of a combined flexure-shear model into the proposed procedure for evaluation of the load–deformation behavior.
The shear mechanism is based on the MCFT model that utilizes softened smeared cracked RC membrane element. The shear deformation mechanism is shown in Figure 5. The shear deformations of the shear wall are concerned with average strains (axial strain
where
where
where
where
where

Shear deformation mechanism: (a) shear deformations in the plastic hinge zone and the top section of the shear wall, (b) average stresses and strains in cracked element, and (c) Mohr’s circle for average strains.
Based on the compatibility of the membrane element, the shear deformation Δ s can be determined as follows
where
or
where
Deformation due to bond-slip of the reinforcement at the fixed end
When a shear wall is subjected to lateral load, the tensile reinforcement will be pulled out from the footing due to the bond-slip of steel bar so that it will cause a fixed end rotation at the wall-footing interface associated with strain penetration (Tastani and Pantazopoulou, 2013, 2015). The fixed rotation leads to additional lateral deformation

Bond slip at the end of the shear wall: (a) strains and stresses along the bond anchorage length, (b) average bond strength, and (c) rotation caused by reinforcement pullout (
The fixed end rotation deformation caused by bond-slip is related to the bond stress along the steel bar, bond anchorage length in the footing, and the reinforcement properties. In this study, the bond stresses
where
where
Analysis procedure
The load–deformation relationship of a shear wall with continuous bars can be evaluated by the following several steps. Considering the interaction between flexure and shear behavior in a shear wall, the procedure for calculation requires iteration. The analysis procedure is based on a displacement controlled approach. In addition, the procedure will terminate when the concrete or reinforcement reaches the limit state, whichever occurs first. Thus, the major calculated steps are introduced as follows:
Define and input structural geometries, reinforcement layout, and all material properties of concrete and steel bars. Calculate the stirrup ratio per unit volume and obtain K and
Create a fiber model and differentiate the regions of confined and unconfined concrete; apply axial load and calculate the axial strain
Assume ultimate lateral deformation (or drift ratio) used to calculate
Assume plastic hinge length
Consider c,
Check the force balances of the two cross sections. If the force balances are not satisfied, then go to step 5. Compute the shear forces
Calculate
Check the convergence of
Calculate
Compute the fracture energies of concrete and reinforcement (
Impose an increment of lateral deformation (or drift ratio) which is required to do the calculations, repeat steps 5 to 9 with calculated
Shear walls with lap-spliced bars at the plastic hinge zone
The behavior and deformational mechanisms of lap splice are complicated in the plastic hinge zone and are different from that of continuous bar. In some cases, it is not allowed to adopt a lap splice in the plastic hinge zone of concrete members due to the seismic defects of a lap splice with poor detailing (Kim et al., 2006). However, lap-spliced joints are still used in RC and precast concrete (PC) shear walls owing to its good economy and convenient construction (Zhi et al., 2017). Consequently, a simple way to deal with a lap splice in the plastic hinge region is proposed in this study.
A common way to cope with a lap splice is to consider the joint as a double anchorage (see Figure 7). So the mechanism of the lap splice is affected by the bond state along the spliced length. The development capacities, including bond strength and deformation capacities, of a lap splice are related to many factors such as spliced length, transverse reinforcement, cover thickness, and concrete strength. Many test studies on the bond behavior of lap splices concerned four-point beam testing where the spliced bars were placed in the constant moment zone (Chun, 2015). Therefore, the two spliced bars are subjected to equal tensile force. Other test studies investigated the mechanics of the lap splice with short length (Cho and Pincheira, 2006; Kim et al., 2006). The full tensile strength of the spliced bars cannot be developed in a lap splice with short length and it does not meet the modern design specifications. Therefore, this study mainly discusses the lap splice with minimum spliced length that can transfer the full strength of the bars at the end of the spliced length. In other words, members will not be failing due to bond failure (such as cover splitting) of the splice in the plastic hinge zone.

Double anchorage and force equilibrium within the spliced length.
The deformation mechanisms of shear walls with continuous and lap-spliced bars are not the same though their lateral load capacities are identical according to limit state analysis of flexural sections and even their load–deformation behavior is somewhat similar. Considering a lap splice in the plastic hinge region of a shear wall, its bond states are illustrated in Figure 7. The initial two flexural cracks are usually formed at the opposite ends of the lap splice with increasing lateral loads due to abrupt changes of reinforcement stiffness. According to the spliced length, the stress and strain states of two spliced bars are classified into three types that are stress and strain distributions with
where
For the time being, the lap splice is viewed as a double anchorage, then the stress states of the two spliced bars after yielding are illustrated in Figure 7 which is similar to the stress distribution of anchor bar in footing in Figure 6. The distributions of stresses along each spliced bar are also indicated by gray (bar 2) and magenta (bar 1) regions in Figure 8, respectively, when the reinforcement has entered the yield phase. Because of the moment gradient along the wall height, the stresses of the two spliced bars are usually not the same. However, for simplification of analysis, still assume

Stress and strain states of two spliced bars: (a) normal strain and stress distributions
It can also be found that
where
where
Due to the complicated stress situation in the plastic lap-spliced region,
Verification of the proposed analytical procedure
The validity of the proposed analytical procedure for the evaluation of the deformational behavior of shear walls with continuous and lap-spliced bars is investigated in the present section. The test data are obtained from a series of experiments conducted in a previous study (Zhi et al., 2017) on cast-in-place (CIP) and precast (PC) concrete specimens with different dimensions, stirrups and reinforcement ratios, and steel bar connection ways. The information on the tested specimens is listed in Table 1. The geometry and reinforcement layout at the joint of several selected test units are shown in Figures 1 and 9. The test setup is drawn in Figure 10. The lateral loads were applied by a 1000 kN hydraulic servo control system that was mounted to a reaction wall.
List of the tested specimens.

Geometry and reinforcement layout at the joint of selected test units (all dimensions in mm).

Test setup of the tested specimens.
All the specimens were tested under low reversed cyclic loading and the proposed analysis approach is based on a monotonic loading approach. The effects of the hysteretic behavior are neglected for the analysis of the specimens. And the effects of the unavoidable natural construction joints between the precast elements are also neglected. The dimensions were 3450 mm tall, 1600 mm wide, and 200 mm thick for specimens SW1-1–6 (Group 1) and 3650 mm tall, 1700 mm wide, and 200 mm thick for SW2-1–3 (Group 2). The axial loads for specimens in Group 1 were 600 kN and for specimens in Group 2 were 750 kN, respectively. It was kept constant till specimens were failure. PC specimens SW2-2 and SW2-3 had higher stirrup ratios than CIP specimens SW2-1. Inversely, CIP specimens SW1-1 had a higher stirrup ratio than PC specimens SW1-2 to SW1-6. The main failure mechanism for these specimens was flexural failure with reinforcement fracture and concrete crushing. The proposed analysis process that has been based on displacement control is implemented for these specimens. The estimated load–deformation (load-drift ratio) relationships compared with the tested results are illustrated in Figure 11, resulting in consistent correlations. Consequently, it might be concluded that the proposed analysis process can be used for the estimation of load–deformation behavior of RC shear walls with continuous or lap-spliced bar in plastic hinge zones based on a displacement-controlled calculation method.

Experimental and analytical results of tested members SW1-1 to SW1-6 and SW2-1 to SW 2-3: (a) SW1-1, (b) SW1-2, (c) SW1-3, (d) SW1-4, (e) SW1-5, (f) SW1-6, (g) SW2-1, (h) SW2-2, (i)SW2-3.
For the analysis of specimens SW2-1–3, their footings were fixed to the strong ground with high strength fining twisted steel bars hoping to prevent the specimens from sliding horizontally. Due to no prestress applied to the footing by four high strength fining twisted steel bars with diameter of 32 mm, the footing were rotated slightly with the lateral loading like the upper specimen. Therefore, to consider the actual situation of experiment and predict the tested results accurately, the lateral deformation resulted from footing rotation and horizontal slip of SW2-1 were also considered in this analysis, which can be calculated by strain compatibility of the flexural foundation beam according to measured data of footing and section analysis. The compression region of the foundation beam can be approximately taken as a quarter of the beam length.
In the analysis process, the ultimate strains of confined concrete and reinforcement were taken as 0.008 and 0.03, respectively. The compression softening factor
Conclusion
In this study, a procedure is proposed for evaluation of the load–deformation behavior of RC shear walls with continuous or lap-spliced bars in the plastic hinge zone under horizontal loading. Determinations of three components of the total lateral deformation that are deformations due to flexure, shear, and reinforcement slip (fixed end rotation) were implemented individually by considering their interaction. The flexural behavior was analyzed by conventional fiber section analysis. The shear response was obtained from MCFT. The influence of shear on flexural behavior was achieved by incorporating a compression softening factor based on the concrete tensile strength into the flexural analysis. The shear model included the impact of the flexural cross section strains and the plastic hinge length on the shear behavior.
In addition, an approach was proposed for analysis of the plastic hinge length based on fracture energies of the materials. The energy of the lateral load applied in the member, when the member is in the post-yield stage, is assumed to be equal to the sum of the fracture energies of materials, including reinforcement and concrete in the plastic hinge region. Thus, the plastic hinge length is derived through the energy balance between the energies caused by the external lateral load and dissipated as inside fracture energies of the materials.
For the shear walls with lap-spliced bars in the plastic hinge region, a simple way to deal with the lap splice was proposed in this study due to its complicated behavior and mechanism. The stress and strain states of a lap splice with different spliced length were analyzed by force equilibrium with an average bond stress model. The established equations regarding bond-slip of the lap splice with minimum spliced length were used to calculate the total response of a lap-spliced shear wall. The validity of proposed analysis procedure was confirmed by comparing the analytical results with previously conducted experiments.
Footnotes
Acknowledgements
The authors would like to thank the anonymous reviews for their constructive comments which help the authors in depth to improve the quality of the article.
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: This work was supported by the National Key Research and Development Program of China (2016YFC 0701703) and National Natural Science Foundation of China (51678136).
