Residual Capacity of Reinforced Concrete Frame With and Without Steel Plate Shear Wall Under Sequence Earthquakes By Incremental Dynamic Analysis

The residual capacity of a damaged structure after the main earthquake is equal to the smallest spectral acceleration of the first mode, which causes local or general failure during the aftershock. In this research, the effect of seismic sequence on the residual capacity of reinforced concrete frame systems and a relatively new system of reinforced concrete frames eqquipped with steel plate shear wall has been investigated. Based on this, four system of 4, 8, 12 and 24 stories, which represent short, intermediate, tall, are modeled in finite element software and subject to three sets of of single and real seismic sequence, taking into account the damage, The effects of mainshock earthquakes have been analyzed under aftershock earthquakes nonlinear increment dynamic analysis. The analysis showed that in the real seismic sequence, the residual capacity of reinforced concrete frame with steel plate shear wall in short and intermediate structures on average 3 times and in tall structures up to 8 times compared to the residual capacity of reinforced concrete frame without steel plate shear wall. Also, in the real seismic sequence, the residual capacity of the structure decreased with increasing the height of short to intermediate structures and intermediate to tall structures, so that this capacity reduction decreased by an average of 65% in reinforced concrete frame with steel plate shear wall and 75% in reinforced concrete frame without steel plate shear wall .

Because aftershocks usually occur shortly after the main earthquake, and some of these aftershocks cause damage to the structure with a destructive force equal to that of the main earthquake, strengthening the damaged structure of the main earthquake in this short time interval is not possible and in order to reduce the level of hazards, it is necessary to evaluate the performance of the damaged structure in general under the main earthquake and aftershocks [1].
In 2011, about 100 aftershocks of magnitude six or more occurred four days after the Tohoku earthquake in Japan, causing structural damage and infrastructure damage [2].
In recent years, due to the increasing development of design methods and the tendency of codes to design based on performance, seismic assessment and performance of structures during earthquakes has become one of the most important issues in earthquake engineering. In the performance-based seismic design method, the structure is designed for different levels of expected performance, related to different levels of earthquake risk. An important step in performance-based design is to evaluate the nonlinear seismic response of structures. Accordingly, it is necessary to study the performance of structures during an earthquake and determine the capacity required to resisting possible aftershocks in response to the structure and the possibility of collapse of the building.
The amount of damage to the structure is directly related to the demand for ductility [3]. Therefore, evaluating the relationship between these two parameters is very important. Current codes only evaluate the ductility demands of a structure in a main earthquake. Research in recent years has shown that the occurrence of aftershocks can significantly increase the demand for ductility of structures and lead to an increase in earthquake risk [4].
The main study of sequence earthquakes since the end of the 19th century has been done by Omori [5]. Omori concluded that the amount of deterioration during an aftershock was inversely related to the deterioration during the main earthquake. One of the first studies by Mahin was conducted in 1980 [6]. Mahin realized that severe aftershocks could double the ductility and displacement demands of many elastoplastic one-degree-of-freedom systems. In addition to Mahin work, there are several other studies to investigate the seismic performance of one-degree-of-free inelastic systems under a variety of ground motions in an earthquake sequence or various ground motions. These include the research of Elnath et al. In 1998 [7], Sanasaka and Kirmidjian in 2002 [8], Amadio et al. In 2003 [9], Das et al. In 2007 [10] and Iancovivi in 2007 [11].
Nicholas and Cornell in 2004 [12] conducted studies to calculate the residual capacity of structures damaged by the main earthquake using nonlinear static and dynamic analyzes of time history. The results of this study showed that in the elastic the capacity calculated by two methods is the same but in the inelastic range the residual capacity curve (for aftershock tolerance) calculated by the static method compared to the dynamic method has a small estimate of the residual capacity. They had a static method in correcting the answer; The capacity calculated dynamically was the main earthquake using the same aftershocks as the main earthquake 1 . On the other hand, since selecting the main earthquake as the aftershock is the worst possible scenario, modifying the capacity curve of the static method is not very correct.
Hotzigeorgio in 2009 [13] proposed nonlinear displacement values suitable for Sequence earthquakes, in addition to Hotzigeorgio in 2010 [14] Studied on the demand for ductility and coefficient of behavior of the nonlinear free degrees system under sequence earthquakes near and far field. Hotzigeorgio in 2010 [15] examined the response of one degree of freedom structures under sequence earthquakes near and far field using artificial earthquakes and concluded that the use of the main earthquake in evaluating the inelastic response of structures contrary to certainty and designers must pay attention to the effects of sequence earthquakes.
Mustafa and Takewaki in 2011 [16] examined simple random models and showed the repetition of sequence earthquakes. In addition, they proposed the effects of sequence earthquakes on freedom multi degree systems. Hotzigeorgio and Liolios in 2010 [17] investigated the effect of seismic sequence on the ductility demand of reinforced concrete frames. According to this research, seismic sequence has an important effect on the response of the structure and ductility demands in seismic sequence can be calculated with high accuracy using a combination of single earthquake demands such as seismic sequence. Garcia et al. in 2011 [18] assessed the demand for drift in steel frames under sequence seismic records in the near and near field. The main purpose of this study was to investigate the relationship between the frequency content of the main seismic and aftershock, the effects of aftershocks on the drift and residual drift of freedom multi degree systems and to evaluate the difference in response between real and artificial sequence earthquakes.
Research by Efraimiadou et al. in 1997 [19], Di Sarno in 2013 [20], and Abdelnaby and Elnashai in 2014 [21] focused on multistory reinforced concrete frames. Almost all of this research has been done on the response of two-dimensional free degree systems.
Zhai et al. in 2012 [3] studied the response of inelastic single degree of freedom systems under a variety of ground motions. They compared the different responses of structural demand parameters in aftershocks, such as maximum acceleration, maximum velocity, maximum displacement and residual displacement with their corresponding values in the main earthquake, and concluded that in potential aftershocks there is an increase in the response of the structure and aftershocks can not always increase both components of maximum acceleration and maximum velocity in the structure in the same way. In other words, for a fixed structure, aftershocks may increase the maximum acceleration, while there is no maximum velocity. Aftershock may also increase or decrease the residual displacement of structural due to the effects of the main earthquake. Therefore, aftershock can change the accumulation of damage to structures by affecting the displacement and sometimes reverse the results. Shin et al in 2014 [22], by analyzing the fragility of reinforced concrete frames, they concluded the change structural performance under sequence earthquake. Hatzivassiliou et al. in 2015 [23] investigated the effects of real seismic sequence with vertical earthquake component on regular and irregular threedimensional concrete buildings. The most important result of this research is that the maximum story displacement of the structure under seismic sequence is more than single earthquakes and the displacement and residual drift increase due to seismic sequence compared to single earthquake.
Raghunandan et al. In 2017 [24] investigated the change in the lateral strength of the structure in helping its residual capacity and according to his research based on incremental dynamic analysis of seven structural models with different lateral design resistances, under a set of seismic sequences Structures designed for maximum lateral forces were more resistant to failure than structures designed for lower seismic levels. Increasing the lateral strength at the design (e.g. using a lower strength reduction coefficient in the US design) can help ensure adequate residual strength in the structure after an earthquake exposure. The results of this study also helped to evaluate how the residual capacity of reinforced concrete frames with lateral strength less than the required regulations for different levels of earthquake damage.
Salim Bahrami and Gholhaki in 2018 [25] conducted an analytical study to evaluate the effects of higher modes of reinforced concrete frames with thin steel plate shear walls under simple pulse. According to their research, the demands of maximum displacement and force for different structures with the ratio of pulse period to the first mode period of vibration is equal to one and with the increase of stories, the effects of higher modes decrease and the most important effects of higher modes are roof shear force and base shear. In 2019, they also investigated the effect of higher modes and degrees of freedom on the coefficient of reduction of resistancein reinforced concrete frames with steel plate shear walls and the coefficient of reduction of earthquake resistance of far and near field, which can be different, for a multi degree of freedom have corrected.
Alberto CUEVAS in 2018 [26] in his P.h d dissertation entitled Seismic residual capacity of reinforced concrete frames, while pointing to the lack of studies on assessing the residual capacity of damaged buildings to withstand future aftershocks, comprehensive information Evidence-based reinforcement to assess the residual damage capacity of reinforced concrete buildings and provide a feasibility study. Analytical relationships of the remaining crack width are also prepared and presented in this research. Kusunoki et al. in 2018 [27] conducted research on a new method for assessing the residual seismic capacity of existing structures with a few degrees of freedom using an accelerometer.  [28] investigated post-earthquake damage and the remaining seismic capacity of structures. In his research, the remaining seismic 2 Back to Back capacity of the structure was evaluated. Damage was assessed based on the level of damage, taking into account the remaining cracks, crushing of concrete.
Recently, the residual seismic capacity has been evaluated by these researchers based on the response spectrum method. Steel plate shear wall is a type of system suitable for earthquake and wind lateral loads consists of a series of separate panels, each panel is enclosed inside two beams and columns and a steel plate is connected to the surrounding elements [29]. Use of steel plate shear wall due to stiffness, strength and significant energy dissipation of seismic, as an effective lateral system in seismic rehabilitation to increase the lateral strength and stiffness of buildings against earthquake (technical rehabilitation strategies) in steel structures is considered. Recently, according to the analysis and design of steel plate shear walls, which requires of environmental elements (beams and columns) with high rigidity, to increase the lateral strength and stiffness of concrete buildings with a moment resisting frame that inherently has such elements is used. Therefore, recently, the reinforced concrete frame with steel plate shear wall has been proposed as a new system against lateral loads [30].
Due to the newness of this system, this system needs more numerical and laboratory studies. This research is also part of the researches to better understand this system.
Numerous studies have been done on a type of structural systems, including resisting moment frames, braced systems, and concrete shear walls under sequence earthquakes. However, the reinforced concrete frames with steel plate shear wall has been studied from various aspects recently, and therefore the effect of sequence earthquakes (near and far field) on this system has not been studied.
In this research, the effect of seismic sequence on the residual capacity of reinforced concrete frame systems and a relatively new system of reinforced concrete frames eqquipped with steel plate shear wall has been investigated and compared. For this purpose, a number of sequence real and artificial records of far and near field, including the main earthquake and aftershock records have been used. The effect of main shocks and aftershocks of real earthquakes on the residual capacity of structural models has been investigated and the results of residual capacity in the two systems have been compared. .

research method
The method used in this research is presented with a flowchart of Figure 1.

Numerical modeling validation 2.1. Numerical modeling validation of reinforced concrete frame equipped with steel plate shear wall
Validation of analytical models is one of the steps of a research. In this research, in order to validate the model, the laboratory study of Choi and Park in 2011 [31] according to Figure 2a has been used. He conducted a laboratory study to investigate the cyclic behavior of walls consisting of boundary elements of reinforced concrete frames and thin steel sheets. In order to ensure the accuracy of the modeling, the numerical model of the laboratory sample was modeled and analyzed in OpenSees finite element software [32]. For modeling, the NonlinearBeamColumn element has been used for the beam and column elements with deformation control, which has the ability to take into account the P-Δ effect and large deformations. The strip method has been used to model the steel sheet [33]. In this method, truss element is used to model the tensile strips. To model the wide plasticity of the elements in the program, the cross sections of the beam and column elements are divided into a number of fibers. Concrete01 and Steel02 materials have been used for modeling concrete and steel materials of reinforcements, respectively. To model the actual behavior of the strips that should not react when pressed, Hysteretic materials are used, which with the three-line behavior in tension and pressure gives the strips the property that does not show resistance when under pressure and allows that the diagonal tensile field of a steel shear wall is well modeled. Also, the discussion of concrete confinement of columns is seen in the model. Numerical results from cyclic loading are compared with laboratory results (Figure 2b). The values of load bearing capacity, initial stiffness and energy absorption determined from the experiment and the corresponding simulated model are presented in Table 1. The comparison between the two diagrams in Figure 1b shows the acceptable accuracy in the modeling phase of this research.

Table 1
Comparison of Finite Element Analysis Results [27] and Choi and Park Model Test [31].
(a) (b) Fig. 2. (a) Geometric and reinforcing details of the concrete frame equipped with thin steel plate infill wall [31] (b) Hysteresis curve of numerical model [32]and Choi and Park laboratory sample [31].

Numerical modeling validation of reinforced concrete frame
Validation of analytical models is one of the steps of a research. In this research, in order to validate the model, the laboratory study of Choi and Park in 2011 [31] according to Figure 3a has been used [31]. He conducted a laboratory study to investigate the cyclic behavior of reinforced concrete frames. In order to ensure the accuracy of the modeling, the numerical model of the laboratory sample was modeled and analyzed in OpenSees finite element software [32]. For modeling, the NonlinearBeamColumn element has been used for the beam and column elements with deformation control, which has the ability to take into account the P-Δ effect and large deformations. To model the wide plasticity of the elements in the program, the cross sections of the beam and column elements are divided into a number of fibers. Concrete01 and Steel02 materials have been used for modeling concrete and steel materials of reinforcements, respectively. Also, the discussion of concrete confinement of columns is seen in the model. Numerical results from cyclic loading are compared with laboratory results (Figure 3b). The values of load bearing capacity, initial stiffness and energy absorption determined from the experiment and the corresponding simulated model are presented in Table 2. The comparison between the two diagrams in Figure 3b shows the acceptable accuracy in the modeling phase of this research.

Modeling of models 3.1. Modeling models in Etabs software and the basics of designing steel plate shear walls
Regarding the classification of structural systems, some have considered the ratio of the height to the smallest horizontal dimension of the structure as a criterion for classifying buildings, and height to the smallest dimension ratios greater than 1.5 π, between π and 1.5 π, between π and 0.5 π and less than 0.5 π are known as high tall, tall, intermediate and short buildings, respectively [34]. Accordingly, in this research, four models of 4, 8, 12 and 24 story with height to the smallest dimension ratios of 0.54, 1.09, 1.63 and 3.26 in the classification of short, short, intermediate and tall structures with a rectangular plan according to Fig. 4 a, is selected with a reinforced concrete frame with steel plate shear wall and high ductility. The height of the stories of the models is 3.4 meters and the roof is considered as a block joist. The construction site of the structures is considered to be a high earthquake risk and soil type III. Concrete used in C22 grade concrete structures has a characteristic strength of 220 kg/cm2 and rebars of type A3 with a yield stress of 4000 kg/cm2. The steel used for equivalent braces is construction soft steel with a yield stress of 2400 kg/cm2. In the analysis and design of the studied structures, the sixth [35] and ninth [36] national building regulations and the Iranian 2800 earthquake standard, fourth edition [37] have been used, according to sixth of the National Building Regulations, the dead weight load of stories and roof is 640 kg/m2, live load of stories and roof is 200 kg/m2 and the load of walls of stories is 600 kg/m. For design of structures, according to the results of laboratory work performed on a reinforced concrete frame with steel plate shear wall and the criteria set by the American loading code ASCE07-2010 [38] for steel plate shear wall, the behavior factor 7 for the reinforced concrete equipped with steel plate shear wall lateral system and , the behavior factor 7.5 for the reinforced concrete frame system [37] was used. In addition, the criteria for designing dual systems under sixth of the National Regulations have been controlled. In order to design thin steel plate shear walls, according to Canadian and American steel regulations, an equivalent brace is considered instead of each steel plate, and after calculating the cross-sectional area of each brace, the thickness of the steel plate is calculated from equation (1).
Where θ is the angle between the brace and the column, L the span of the frame, Ab is the equivalent cross sectional area of the brace, α is the angle of formation of the diagonal tensile field in the steel plate After determining the thickness, each plate is converted into a number of diagonal strips, the cross section of each strip is obtained from equation 2:   Where Mfpb is the plastic anchor of the cross-section of the beam and σty is the final stress of the diagonal tensile field, which for thin plates is equal to their yield stress. Due to the small difference in the intensity of the diagonal tensile field between two adjacent stories, control of this relationship is required only for the end beam, but if the difference between the diagonal tensile fields between two adjacent stories is large, the relationship should be controlled for the middle beams.
In order to ensure that the perimeter columns can withstand the stresses due to the environmental loads along with the stresses due to the tensile field effect, it is necessary to check equation 5 for the plastic anchor of the columns: fpc M Equivalent tensile bracing was used to model of the steel plate shear wall in the design stages due to the impossibility of modeling the steel plate in the model geometry and its analysis by ETABS software. In the last two stories of the 24 story model, due to the negative shear in the braces, the use of a steel plate shear wall in these stories was avoided in the structure and only the special reinforced concrete frame system was used to deal with the lateral force.

Modeling of models in OpenSees
For modeling of beam and column elments (NonlinearBeamColumn )the ability to take into account the effect P-Δ and large deformation, the beams and columns have been used. For modeling of steel plate shear wall strip method is used. In this method for the of the tensile strip of the truss element is used. For the concrete and steel reinforcement materials respectively of materialsConcrete01 and Steel02 has been used. For the actual behavior of strips is used the materials Hysteretic which allows the use of diagonal tension steel plate shear. For example, the bracing of the four stories in table 3 is presented. In the eight-story model, the thickness of the plate is calculated to be 2 mm for stories one to six, 1.6 mm for the seventh story and 1.1 mm for the eighth story. In the twelve-story model, the thickness of the plate is 1.1 mm for the first to fourth stories, 1.8 mm for the fifth to eighth stories, 1.6 mm for the ninth and tenth stories, and 1.4 mm for the eleventh story. And for the twelfth story is 1.1 mm. In the twenty-four story model, the thickness of the sheet in the first story is 4.4 mm, in the second to eighth stories, 3.2 mm, in the ninth to fifteenth stories, 2.9 mm, in the sixteenth story, 2.6 mm, in the seventeenth to eighteenth stories, 2.1 mm, in the nineteenth to twentieth stories, 1.8 mm, in the twenty-first and twenty-second stories, 1.6 mm, in the simplest and third stories, 4.4. 1 mm and on the twenty-fourth story, 1.1 mm was calculated and used in the analysis. Also concrete confinement of columns is cosidered in the model. OpenSees uses the distributed plasticity by the fiber element.
Regarding geometric nonlinearity, it should be said that the effects of geometric nonlinearity are defined by the transfer matrices that are a feature of OpenSees. In the mentioned program, after defining the geometry of the model, the gravitational analyzes are gravitationally analyzed (non-linear static) and by setting the time in the amplitude of the problem to zero before performing the nonlinear dynamic analysis, the gravitational load values remain constant in subsequent dynamic analyzes. The effect of P-Delta is considered in the analyzes. Due to the fact that in the case of peripheral frames with considerable stiffness and random rotation relative to the middle frames, for the analysis of the structure designed in Etabs, peripheral frames and their specifications have been used for two-dimensional nonlinear analysis. In order to take into account the effects of P-Delta, the Leaning Column method has been used. In fact, the concept of a Leaning Column is used to consider the effect of P-Delta due to gravitational loads. In this method, an additional bay is defined. Columns between which there are rotating springs with low stiffness are connected to the main frame with high-rigidity beams. The purpose of defining rotary springs with negativity is negligible so that the leaning columns do not absorb momentum. Gravity loads are then applied to these columns. In the nonlinear dynamic analysis, the selection of a time interval of 30 seconds between the main earthquakes and aftershocks with zero acceleration amplitude is considered in order to stabilize the frame under the effect of the main earthquake excitation. 5% damping was applied for all models.
In order to evaluate the seismic demands of reinforced concrete frames equipped with steel plate shear walls, the parameters of residual capacity of models have been selected to assess the seismic bahavior of models.

Selected seismic records in nonlinear dynamic analysis
Sequence seismic records have been applied to the models in two ways, including real seismic records and artificial seismic records. In the real seismic records method, using the sequence seismic records recorded in the recorded stations and in the artificial method due to the unavailability or insufficient number of sequence real seismic records an artificial aftershock have been applied to the models. Due to the recording of a small number of sequence seismic records in different seismic regions of the world, to investigate the effects of this event on the performance of structures with different construction conditions and also to more realistically investigate the effect of sequence seismic on the performance of structures demand to use artificial records is more felt [34]. According to studies, the most important approaches to the production of artificial seismic records using the real seismic record are the back to back (repetitive) approach and the randomized approach [39][40].
In this study, in the method of real seismic records, three groups of seismic records including records of real main earthquakes with effective peak acceleration and their aftershock records, other real main seismic record with nearly effective peak acceleration and their aftershock records, records of a number of real main earthquakes used in the research of researchers and their aftershock record have been used and the method of artificial seismic records are some of the important seismic records in the far and near field is studied in the form of seismic scenarios.

Real seismic records (seismic records with effective peak acceleration)
In order to investigate the effect of seismic sequence on the response of reinforced concrete frame with steel plate shear wall, real and recorded Sequence records have been used. The mentioned records have been selected based on the Garcia in 2012 [39] and Ghodrati and Manouchehri Dana study in 2005 [41] in which the effective peak acceleration (EPA) parameter is included. The definition of the EPA parameter, which is similar to the maximum PGA ground acceleration of the acceleration type, is equal to the devided of the average spectral acceleration of an earthquake with 5% damping over a period of time of 0.1 second to 0.5 seconds by the standard magnification factor of 2.5. EPA is like the peak ground acceleration PGA of the acceleration type. According to the Ghodrati study and Manouchehri Dana [41], among the proposed parameters for selecting the critical seismic scenario such as peak ground acceleration (PGA), effective peak acceleration (EPA), peak ground velocity (PGV), peak ground displacement (PGD), duration Earthquake and Arias intensity 2 criterion, due to having the highest correlation with other parameters, has been introduced as the most appropriate criterion for selecting the most critical scenario. In such a way that the mentioned records are more critical in terms of the impact on the behavior of the structure and its response than other sequence seismic records recorded in other stations. In fact, the reason for using the records presented in this research is to avoid creating uncertainties in introducing the characteristics of the following earthquakes by using artificial methods and sometimes proposed models based on Gutenberg-Richter, Amory and Beth's laws. Thus, in all selected scenarios from the PEER 3 database, the magnitude, duration, time interval between sequence seismic records, and other seismic properties are real. The conditions for selecting sequence real earthquakes used in this paper, in addition to the fact that the recording stations of the first and second earthquakes are the same and the time interval between them is less than ten minutes, is the discussion of the effective peak acceleration of these earthquakes. So that in the sequence seismic records shown in Table 4, both the first earthquake and the second earthquake have a EPA among the other records recorded for the desired earthquake. This group of sequence earthquakes is called real seismic scenarios with EPA.

Table 4
Result for scaling the first category records of seismic scenario (with maximum EPA).

Real seismic records (seismic records with nearly maximum EPA)
Due to the limited number of selected records, a sequential seismic scenario called with nearly maximum EPA (Table 4) is considered, in which the first or second earthquake is ranked second or finally third in terms of the amount of EPA relative to the maximum.With these considerations, 11 sequence records introduced in tables 3 and 4 records have been selected for nonlinear dynamic analysis. According to the unit design for records with different PGA intensities according to the Iranian earthquake standard 2800 [37] regarding the comparison of the values of the acceleration response spectra of the earthquakes used in nonlinear dynamic analysis with the values of the standard design spectrum, each sequence to reach 0.3g were scaled. The calculated scale coefficients of earthquakes for each model are presented in the last column of Tables 4 and 5. Table 5 Result for scaling the second category of seismic scenario (with nearly maximum EPA).

Real seismic records (seismic records used by researchers)
Four records of earthquakes used by researchers [17,42] as their main earthquake and aftershock are presented in table 6.

Table 6
Records with after-earthquake acceleration to main earthquake ratio greater than one [17,42].

Effective characteristics of earthquakes in the response of structures
In nonlinear time history analysis, the behavior of the structure is observed during the earthquake time and this analysis shows the more real behavior of the structure during the earthquake compared to other analyzes. In this type of analysis, the effect of peak ground acceleration, frequency content and time of earthquake (effective duration of earthquake) is well observed.

4.4.1.Peak ground acceleration
In this study, to investigate the effect of peak ground acceleration or PGA, a relatively wide range of the ratio of peak ground acceleration aftershock to a mainshock of strong earthquakes has been considered. These ratios are: 0.04, 0.18, 0.21, 0.28, 0.32, 0.38, 0.43, 0.46, 0.74, 0.97, 1.23 for real seismic records of groups 1 and 2 and ratios of 1.06, 1.15, 1.30 and 1.89 for accelerometers of group 3.

Frequency content
peak ground motions have a greater effect on the amplitude of vibrations. If the frequency content of earthquake and vibrational frequencies of the structure are close, earthquake will cause the greatest intensification in the structure. The frequency content is reflected in the acceleration response spectrum or the Fourier spectrum of its acceleration amplitude. The predominant period of earthquakes Tg was obtained using SeismoSignal [43] software.

Duration of strong ground
Research shows that the peak ground acceleration of earthquake alone is not a good quantity to determine the earthquake destruction. Because the damage caused by an earthquake is related to the amount of energy of the earthquake movements in a building. Accordingly, the duration of earthquake based on the amount of input energy to the structure is defined as an important parameter in determining the response of the structural.
The duration of earthquake in most cases has a significant effect on the inelastic deformation of the structure and therefore in determining the amount of input energy to the structure, linear and nonlinear analysis of the structure response and also in the statistical method of determining the structure's response to earthquake is widely used. There are different methods for determining the duration of earthquake using the characteristics and effective parameters of seismic acceleration. One of these methods is the Trifunac-Brady method, which has defined the period of strong vibration as the time interval in which a significant contribution is added to the square integral of the accelerations, called the acceleration intensity. They have chosen the time interval between the 5% and 95% stocks as the period of strong. In this research, the effective duration of earthquakes was obtained by this method with the SeismoSignal software [43] and was used in nonlinear dynamic analysis of models.

The results of the analysis of the main single earthquakes
After modeling, the models are subjected to earthquake of the main single earthquakes of group one seismic scenario and the results of the analysis are presented for seismic evaluation of the models. 5. 1. Assessment of residual capacity of models under real seismic sequence for group I earthquakes with effective peak acceleration.
The models under the main single-scale acceleration of this group (main earthquake Figure 7) are nonlinear dynamic analysis and their residual capacity is calculated and presented in Table 7. b a As can be seen from Table 7, the residual capacity of tall structures with reinforced concrete frame system equipped with steel shear walls is significantly higher than the residual capacity with reinforced concrete frame system without steel shear walls.

Table 7
Residual capacity of models under real seismic sequence for group I earthquakes with effective peak acceleration. In table7, Frame equal reinforced concrete frame without steel shear wall, SPSW reinforced concrete frame equipped with steel plate shear wall and Res cap, residual capacity and numbers before Freme and SPSW indicate the number of stories of model.

2. Assessment of residual capacity of models under real seismic sequence for group II earthquakes with nearly maximum EPA.
From the results of incremental dynamic analysis under the group II seismic sequence, the residual capacity of the models is calculated and presented in Table 8.

Table 8
Residual capacity of models under real seismic sequence for group II earthquakes with nearly maximum EPA. Table 8 Continous of Table 8 5. 3. Assessment of residual capacity of models under real seismic sequence for group III used by researchers.

Continous of
From the results of incremental dynamic analysis under the group III seismic sequence, the residual capacity of the models is calculated and presented in Table 9.

Table 9
Residual capacity of models under real seismic sequence for group III earthquakes Continous of Table 9 6.Comparing the residual capacity of reinforced concrete frame systems with and without steel plate shear wall.
In this part of the research, the residual capacity of the models under the earthquakes of groups I, II, and III is presented comparatively. Table  10 shows that the average residual capacity of models with reinforced concrete frame system with steel plate shear wall is much larger than the residual capacity of models with reinforced concrete frame system without steel plate shear wall. According to this table, the real seismic sequence, in short structures, the average woven capacity of models with reinforced concrete frame system with steel shear walls is 4.60 and in reinforced concrete frame system without steel plate shear walls is equal to 1.53 Which shows a 201% increase, or in other words, the capacity of the reinforced concrete frame system with steel plate shear walls is 3 times compared to the reinforced concrete frame system without steel plate shear walls.
In intermediate structures, on average, the residual capacity of models with reinforced concrete frame system with steel plate shear walls is 1.26 and in reinforced concrete frame system without steel plate shear walls is 0.43, which shows an increase of 193% or In other words, the capacity of the reinforced concrete frame system with steel plate shear walls is equal to 2.9 times compared to the reinforced concrete frame system without steel plate shear walls. In tall structures, on average, the residual capacity of models with reinforced concrete frame system with steel shear walls is 0.55 and in reinforced concrete frame system without steel shear walls is 0.07, which shows an increase of 686% or In other words, the capacity of the reinforced concrete frame system with steel plate shear walls is 7.86 times equal to that of the reinforced concrete frame system without steel plate shear walls.
In total, according to the analysis, the residual capacity of the reinforced concrete frame system with steel plate shear walls is on average 4.5 times larger than the reinforced concrete frame system without steel plate shear walls. It seems that in seismic sequence, the performance of reinforced concrete frame system with steel plate shear walls is much better than the performance of reinforced concrete frame system without steel plate shear walls. Table 10 Comparison of residual capacity of models under earthquakes of groups I, II and III 7. Strategies to deal with serious damage to structures under sequence earthquakes.
The use of steel plate shear walls in the reinforced concrete frame, due to the high residual capacity compared to the reinforced concrete frame can significantly reduce the structural damage of the reinforced concrete frame. Especially in tall structures, the use of steel plate shear walls is doubly important. Figure 1 Flowchart of methodology (a) Geometric and reinforcing details of the concrete frame equipped with thin steel plate in ll wall [26] (b) Hysteresis curve of numerical model [27]and Choi and Park laboratory sample [26]. (a) Reinforced Concrete Frame [26] (b) Hysteresis curve of numerical model [27] and Choi and Park laboratory sample [26].