A New Elastoplastic Time-History Analysis Method for Frame Structures

+is study aimed to analyze the formation and application of the time-domain elastoplastic response spectrum. +e elastoplastic response spectrum in the time domain was computed according to the trilinear force-restoring model. +e time-domain elastoplastic response spectrum corresponded to a specific yield strength coefficient, fracture stiffness, and yield stiffness. However, the force-restoring models corresponding to different structural systems and the states of the structural systems at different moments were not the same. +erefore, the dynamic characteristics of a particular periodic point corresponding to a particular structure weremeaningful for the elastoplastic response spectrum. In addition, the curve in the time-domain dimension along the periodic point truly reflected the real-time response of the structure when the structure encountered a seismic load.


Introduction
e elastic response spectrum provides the maximum structural response to an earthquake in the elastic stage. However, it cannot be directly used for structural seismic internal stress analysis. erefore, understanding the timedomain elastoplastic response spectrum of the structure is necessary to understand the nonlinear performance and development history of the structure [1]. Elastoplastic response analysis of structures is generally based on displacement. At present, the commonly used elastoplastic seismic response calculation methods are the dynamic elastoplastic time-history analysis method, static elastoplastic analysis method, and simplified elastoplastic analysis method [2][3][4][5][6]. e abovementioned methods are widely used as traditional methods to analyze elastoplastic deformation of structures, but the simplified elastoplastic analysis method is an empirical amplification method. Because of the complex and changeable forms of structures and the various influencing factors of ground motion, it is difficult to consider the influence of different factors by only one amplification coefficient, and the accuracy of the analysis is poor. e current spectrum curve cannot well reflect the actual cyclic stress process of the structure. erefore, this paper presents a method to solve the development process of elastoplastic seismic force of structures by means of mode decomposition method using time-domain elastoplastic response spectrum.
e force-restoring model of the structure needs to be determined first to solve the timedomain elastoplastic response spectrum of the structure. en, the elastoplastic dynamic differential equation of the structure can be derived. Finally, the structural response spectrum can be analyzed by the numerical integration method [7,8].

Elastoplastic Force-Restoring Model
e force-restoring model is a function between internal forces and deformations obtained by mathematical induction based on numerous experiments [9,10]. e force-restoring curve is the embodiment of the structure or component strength, stiffness, and comprehensive ductility. It also reflects the energy absorption ability of the structure or component. us, it is the basis for studying the elastoplastic properties of the structure [11][12][13][14][15]. For the reinforced concrete structure, the force-restoring performance can be understood from materials and components.
e force-restoring model of the material reflects the constitutive relationship between steel and concrete.
e force-restoring model of the components reflects the relationship between internal force and deformation hysteresis. Two rules are applied to select the force-restoring model of reinforced concrete structures: to ensure the accuracy of the model, which can reflect the structural stress state, and to express in a mathematical form that is able to analyze the structure using nonlinear mathematics. Recently, the most common force-restoring models for reinforced concrete structural components include the bilinear model ( Figure 1) and the degenerate trilinear model (Figure 2). e degenerate trilinear model is determined using parameters such as fracture load P c , yield load P y , elastic stiffness, fracture stiffness, and postyield stiffness. e stiffness of the model during unloading remains unchanged, which is always the tangential stiffness between the yield point and the displacement origin [16]. e degenerate trilinear model is more reasonable and accurate in simulating the variation in internal forces and displacements in reinforced concrete structural components in an elastoplastic state during large earthquakes.

Numerical Expression of Degenerate Trilinear Restoration Force Model
Not only the mechanical properties of structural components but also the deformation properties are included in the force-restoring curve model [17]. is study used a degenerate trilinear force-restoring model to describe the structural characteristics such as stiffness degradation and ductility. e force-restoring equation is assumed to be where kis the slope of the force-restoring equation and brepresents different possible stages of the restoring force and displacement of the particle during earthquakes. Different stages represented the initial stiffness as k e , the fracture stiffness as kc, and the postyield stiffness as k y . e ratio of the fracture stiffness kc to the initial stiffness ke was denoted as α, and the ratio of the postyield stiffness k y to the fracture stiffness k c was denoted as β. e load when the structural component was fractured was denoted as P c , and the load at the time of yield was denoted as P y . For a structural component whose initial stiffness k e , fracture stiffness k c , postyield stiffness k y , fracture load P c , and the yielding load P y were all known, its degeneration trilinear force-restoring model was then determined. As shown in Figure 2, the strain or displacement corresponding to the fracture of the component was assumed as u1, and the strain or displacement when the component yielded was u 2 . e numerical expressions of different stages were given as follows: (vii) Stages 9-10: f(x) � (αk e P y /(P y − (1 − α)P c )) · x + (P y − αk e P y u 9 /(P y − (1 − α)P c )) e reactions of the structural components in different stages can be simplified using the mathematical expressions of the aforementioned stages, and the stress state during the whole process can be expressed. In this study, the maximum seismic response of the single-degree-of-freedom system with different periods was calculated based on the aforementioned force-restoring model when calculating the elastoplastic response spectrum. In the computing  procedure, the maximum elastic response of the structure was first calculated. en, the force-restoring curve was determined using the yield strength coefficient of the structure and setting different fracture strength ratios and yield strength ratios. us, the corresponding reaction analysis was performed according to the recordings of different earthquake ground dynamics.

Elastoplastic Dynamic Differential Equation
e elastoplastic dynamic differential equation for a singledegree-of-freedom system during horizontal earthquakes is as follows: where m is the mass; c is the damping; x is the particle displacement; t is the particle displacement; u g is the ground motion;f(x) is the elastoplastic restoring force; dx/dtis the relative velocity of the system;d 2 x/dt 2 is the relative acceleration of the system; and d 2 u g (t)/dt 2 is the ground movement acceleration.
During the deformation of single-degree-of-freedom systems, the variation in restoring force is nonlinear [18].
is study used the trilinear hysteresis curve for the forcerestoring model and assumed that the damping remained constant during the deformation of the system. After substituting the restoring force into the single-degree-offreedom elastoplastic dynamic differential equation of the system, the following expression was obtained: Dividing both sides of the equation by m yielded the following: erefore, the dynamic differential equation of the single-degree-of-freedom system was established based on the degenerate trilinear model: where ξ is the damping ratio; ω is the natural frequency of the single-degree-of-freedom system; and r is 1 in stages 0-1, αin stages 1-2, αβ in stages 2-9, αP y /(P y − (1 − α)P c ) in stages 3-4, and stages 9-10, and

Numerical Integration Method of Elastoplastic Response Spectrum
e earthquake acceleration and structure response of the elastoplastic system are both nonlinear [19,20], which cannot be described by an analytical solution solved directly from the dynamic differential equation, such as an elastic system. erefore, it can only be computed by the step-by-step integration method. is study used the Newmark-β method to solve the increments based on the quasistatic equation.
e Newmark-β method is used to calculate the structural dynamic response according to the assumed acceleration variation within the time increment. It is a stepwise integration method. Its advantage is that it avoids any superposition application and can be well adapted to nonlinear analysis. Owing to various ways of acceleration variation in the time increment, the basic assumptions of the Newmarkβ method are as follows: is study used the average acceleration method (i.e., β is 0.5 and c is 0.25).
e acceleration was assumed as the average value during time interval t i -t i +1 : us, the speed and displacement were as follows: Because of a sudden change in stiffness in the force-restoring model, a cumulative error occurs during the reciprocal integral calculation if the breakpoint is not accurately processed, reducing the reliability of the computing result. e breakpoint often appears within a certain integration step. A basic assumption of the Newmark-β method is that the stiffness of the structure in each calculation section is linear. erefore, extra efforts are required to reach the stiffness breakpoint. is study used the precise breakpoint processing method [21]. e polygonal line in the degenerate trilinear model in Figure 2 indicates two types of breakpoints in the model: the loading point ( Figure 3), in which the symbol of speed does not change on both sides of the breakpoint, and the unloading point (Figure 4), in which the sign of speed changes on different sides of the breakpoint. e aforementioned characteristics were used to determine the loading point based on the increase in absolute displacement and determine the unloading point using the sign change of speed. e elastoplastic dynamic differential equation in increments was expressed as follows: Assuming that the breakpoint appeared at a certain time t + Δt 0 in the (t, t + Δt) time period, the displacement at the loading breakpoint was expressed as follows: Advances in Civil Engineering which was converted into the following incremental form: Equation (13) could be used to obtain the following expression: e incremental form of the elastoplastic dynamic differential equation (11) was used to obtain the following: By substituting equation (16) into equation (13), the displacement increment in t to t + △t 0 period was obtained: By substituting equation (15) into equation (13), the following expression was obtained: Substituting equation (18) into equation (17) yielded the following expression: Equation (19) was organized into the standard form as follows: e coefficients of A 1 , B 1 , C 1 , and D 1 in the aforementioned equation were calculated as follows: For the unloading breakpoint, the velocity at the breakpoint was zero: which was converted into the following incremental form: Equation (23) yielded the following:  Advances in Civil Engineering (du(t + Δt 0 ))/dt was expressed using the incremental Newmark-β method: Equation (25) yielded the following: By substituting equations (13), (24), and (26) into equation (11), the following expression was obtained: e equation (27) is organized into the standard form as follows: e coefficients of A 2 , B 2 , C 2 , and D 2 in the aforementioned equation were calculated as follows: e coefficients of the loading and unloading breakpoints were related to displacement, velocity, acceleration, and structural stiffness at time t, as shown by equations (20) and (28). erefore, four coefficient values were determined based on the state of the structure at time t and the breakpoint type. en, the time of breakpoint was accurately solved using equations (20) and (28).
For any seismic wave, the elastic absolute acceleration response spectra were first analyzed to get the maximum absolute acceleration s a of the elastic system. en, the elastic seismic response analysis of the structure and actual reinforcement structure of each component was performed by calculating the elastic seismic response of the structure. After that, the yield strength coefficient ζ of each floor was obtained by comparing the floor shear-bearing capacity calculated using real component reinforcement and standard material strength with floor elastic seismic shear stress calculated using rare earthquake interaction (i.e., the ratio of floor elastic seismic shear stress calculated using yield shearbearing capacity of different floors and rare earthquake seismic interaction). e corresponding yield load was obtained as follows: Using the initial stiffness of the structure and setting the ratio of the postfracture stiffness to the elastic stiffness of the structure as β and the ratio of stiffness after yielding to fracture stiffness as p, the postfracture stiffness k c and stiffness after yielding k y were obtained as follows: (31) Using the aforementioned known structural floor yield load, elastic stiffness, and postfracture stiffness and by setting the ratio of the fracture displacement and yield displacement of the structure as α, the calculation formula of the floor displacement when the structure was fractured was obtained as follows: e corresponding yield displacement calculation formula is as shown in equation (33): e corresponding trilinear force-restoring model was determined by substituting equation (30) into equation (33). en, the displacement, velocity, and acceleration of the system at any time under seismic wave excitation were determined by numerical integration of the elastoplastic dynamic differential equation using the Newmark-β method [22] according to the displacement u 0 , speed du 0 /dt, and acceleration d 2 u 0 /dt 2 of the single-degree-of-freedom system for any seismic wave recording. e seismic wave is divided as Figure 5 to the same time interval △T, and then the full period seismic response spectrum is calculated for the time period after each segmentation. e maximum value of seismic influence coefficient of each period point in t m− 1 -t m time slot is a m1 a m2 . . . a mn . When the natural vibration period of the single-degree-of-freedom system is T n , the maximum value of seismic influence coefficient in each time slot (0-t 1 , t 1 -t 2 , . . ., t m− 1 -t m ) is a 1n a 2n . . . a mn . e maximum value of seismic influence coefficient corresponding to natural vibration period T n in t m− 1 -t m time slot is a mn , thus the time domain response spectrum matrix of each period point in each time period is obtained. Hence, the seismic force of the time domain and the elastoplastic displacement of the layer of weak floors of the structural system were analyzed using the aforementioned data.
Advances in Civil Engineering e elastoplastic acceleration response spectra were computed using the lab-made computational programs based on the earthquake recordings of the level 9.0 EW direction earthquake in East Japan AOM021 on March 11, 2011 (300-s seismic recording length; recording seismic peak acceleration amplitude was 220 cm/s 2 during 7 degrees 0.1 g earthquake; damping ζ was 0.05). e calculation example used in this paper is a frame with two spans and three floors with both directions of 6 m, and floor height is 3.2 m. e frame column section is 400 mm × 400 mm, the frame beam section is 200 mm × 500 mm, the floor dead load is 2.0 kN/ m 2 , the live load is 3.5 kN/m 2 , the roof dead load is 2.0 kN/ m 2 , the live load is 0.5 kN/m 2 , and the line load of side beam is 6.0 kN/m). e elastic analysis in this study showed that the yield strength of the elastic factor was set as 0.8 and the ratio of fracture stiffness to initial elastic stiffness of the hysteretic curve of concrete structures was set as 0.5. e ratio of postyielding stiffness to fracture stiffness was 0.1, and the ratio of fracture displacement to yield displacement was 0.5. Figure 6 is a time-domain elastoplastic response spectrum curve generated by calculating the time-domain seismic response spectrum with a selected time interval of 10 s and a period interval of 0.02 s. e corresponding matrix comprised 600 rows and 30 columns. Figure 7 is a timedomain elastoplastic seismic impact coefficient surface with a time interval of 0.01 s and a period interval of 0.02 s. e corresponding seismic impact coefficient matrix comprised 600 rows and 30,000 columns. Figures 6 and 7 show the time-domain absolute acceleration seismic influence coefficient surface generated using the same seismic acceleration record but different time intervals. A comparison of the two figures showed that the period and time of the maximum peak of seismic influence coefficient were basically the same, and the calculated peaks are all 0.38. However, more peaks appeared on the coefficient surface in the time domain due to the very short time interval selected in Figure 7. Horizontal seismic force analysis in the time domain using this response spectrum surface was accurate enough to capture the change in 0.01 s time interval, but the corresponding computation cost was quite large. For example, if the response spectrum curve generated using this seismic record was used for calculating mode decomposition calculation, the number of time-domain seismic forces under each mode reached 30,000; the total number of timedomain seismic forces reached 30,000 n (n is the mode number). Moreover, the vibration modes of 30,000 n seismic forces were recombined and calculated. Such large-scale computations were obviously quite inefficient; even using computers to perform the computation took a long time. erefore, the time-domain response spectrum matrix generated at intervals of 1, 5, and 10 s was recommended in this study to be used to analyze the time-domain seismic force according to the stage of the project and required calculation accuracy, so as to meet the design requirements of the corresponding project stage. For example, as shown in Figure 6, the time interval is 10 s and the corresponding elastoplastic response spectrum matrix has only 30 seismic forces in each mode. Although it could not reflect the trend of seismic force and deformation in the 10 s time period, the peak value of seismic force in the 10 s time period was not lost, roughly reflecting the change in seismic force in the time domain of 300 s seismic record. e premise of studying the elastoplastic response spectrum is that the load corresponds to the yield strength of the structure system. When the load response of the particle exceeds its cracking load for the first time, the stiffness of the system will be reduced, and the natural vibration period will be correspondingly extended. erefore, the peak acceleration corresponding to each period of the elastoplastic response spectrum curve is smaller than that of the elastic response. As the seismic vibration intensity increased, the stiffness of the system underwent a second reduction when the load response of the particle exceeded the set yield load f y . en, the stiffness of the system was reduced for the second time, and the corresponding natural vibration period was extended for the second time, which again reduced the acceleration response of the system subjected to a seismic dynamic load. Using EPDA (Elastoplastic Dynamic Analysis) software developed by the Chinese Academy of Architectural Sciences and using the same seismic records and structures as those of Figures 6 and 7, the peak value of elastoplastic dynamic time-history analysis is 0.39, which shows that the time-domain elastoplastic seismic response surface calculated in this paper is effective and practical.
e aforementioned analysis showed that the elastoplastic response spectrum curve corresponded to a specific yield strength coefficient, fracture stiffness, and yield stiffness reduction factor. erefore, the corresponding periodic point of the specific structure was meaningful for the elastoplastic response spectrum curve. In addition, the curve corresponding to the time-domain dimension along the periodic point truly reflected the time-dependent reaction when the specific structure encountered the seismic dynamic load. With response spectrum curves corresponding to time-domain dimensions, the development history of structural seismic force can be obtained by mode decomposition response spectrum method, and then the time-domain values of interfloor shear force and displacement can be obtained. Because the maximum elastoplastic seismic influence coefficients of interfloor shear force and displacement calculated by ordinary mode decomposition response spectrum method are adopted under each mode of vibration, while the elasticolastic seismic shear force and displacement mining in time-domain are adopted. e maximum elasticolastic seismic impact coefficient in each period is used, and the corresponding calculation accuracy is greatly improved. e maximum elastoplastic interstory displacement angle of the frame structure is 1/78. e maximum interstory displacement angle of the structure is 1/76 obtained by elastoplastic time history analysis under the same seismic wave. erefore, the calculation results of the two methods are basically consistent. e calculation theory in this paper meets the needs of engineering practice.

Conclusions
is study analyzed the formation and application of timedomain elastoplastic seismic response spectrum in depth. First, the elastoplastic force-restoring model, the numerical expression of the degenerate trilinear force-restoring model, and the treatment of the inflection point for reinforced concrete structures were analyzed. en, the elastoplastic response spectrum curve of any seismic wave for a certain yield strength coefficient was solved using the Newlab-β numerical integration method using Matlab according to the basic principles of the elastoplastic dynamic differential equation. e practicality and effectiveness of time-domain elastoplastic seismic response spectrum were illustrated using computation examples. e main conclusions of this study were as follows: (1) For the time-domain elastoplastic seismic response spectrum surface generated by different time intervals using the same seismic acceleration record, the period and time of occurrence of the maximum peak of the response spectrum were basically the same. For a shorter time interval, the time-domain range displayed a more common trend, and more peaks appeared in the surface figure. However, the minimum time interval could not be less than the recording step size of the seismic acceleration. (2) e elastoplastic response spectrum curve corresponded to the specific yield strength coefficient, fracture stiffness, and yield stiffness reduction factor. erefore, the corresponding periodic point of the specific structure was meaningful for the elastoplastic response spectrum curve. e curve corresponding to the time-domain dimension along the periodic point truly reflected the time-dependent reaction when the specific structure encountered the seismic dynamic load.
(3) e time-domain values of shear force and displacement between floors were analyzed using the time-domain elastoplastic response spectrum matrix. e maximum elastoplastic seismic influence coefficient of each mode was used for calculating the interlayer shear force and displacement by the ordinary mode decomposition reaction spectrum method. Also, the maximum elastoplastic seismic influence coefficient in each period was used to calculate the time-domain elastoplastic seismic shear force and displacement. Hence, the corresponding computation accuracy was greatly improved. e accuracy of the new method meets the needs of engineering. e calculation efficiency is obviously higher than that of the elastoplastic dynamic timehistory analysis method.

Data Availability
All the data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare no conflicts of interest.