Seismic Response Analysis of Prestressed Concrete Rocking Frame

: In order to investigate the seismic performance of prestressed concrete rocking frame (PCRF), a theoretical model based on rigid body is established for a one-story single-span PCRF. The PCRF studied in this paper has the connecting interfaces set at the column feet and at the inner faces of the beam–column joints, allowing the columns to be uplifted with the accompanying separation of the beam–column interface and rotation of the beam and column around the interface. The tendons are arranged along the centerline of the beam and columns. The connections between the beam and columns and the anchoring of columns are accomplished by prestressing the tendons. The theoretical model consists of a rigid beam, rigid columns and elastic tendons. The governing motion equation of the PCRF is derived based on the model and a numerical solution of the equation of motion is obtained. The energy dissipation of the PCRF is analyzed and the calculation method for the coefﬁcient of restitution is derived. Time history analysis and parameter analysis of seismic response of the PCRF are conducted and the results show that the PCRF has promising seismic behavior.


Introduction
The rocking structure concept was first proposed by Housner, based on the surveys of earthquake damage [1,2].With continuous research of resilient structures, the rocking structure has received increasing attention from researchers [3][4][5][6][7].A rocking structure generally has minor damage after an earthquake and its main components can even be undamaged [8].Frames are a very common form in concrete structures.The energy dissipation of a traditional frame during an earthquake mainly depends on its structural ductility.Such an earthquake damaged frame usually requires a long repair time with great cost.To deal with this problem, a type of rocking frame is presented based on the conceptual fusion of the rocking structure and frame.
Makris et al. [9] analyzed the rocking response and stability of the rocking frame.It was found that the stability of the rocking frame is strengthened with an increase in cap beam mass.The seismic response of a plane rocking frame with geometric asymmetry was studied by Dimitrakopoulos et al. [10].Although the rocking mechanism of a symmetrical rocking frame is very different from that of an asymmetric one, the influence of structural asymmetry on the stability of the rocking frame is minute.Vassiliou et al. [11] suggested a new finite element modeling method for a deformable rocking frame based on zero length fiber cross-section elements and Hilber-Hughes-Taylor energy consumption algorithm.
Prestressed tendons are often arranged in the column or beam of a rocking structure in actual projects to enhance the self-centering ability of the structure, to increase the structural stability and to reduce the residual deformation.Based on this measure, different forms of prestressed concrete rocking frame (PCRF) were presented by researchers.Priestley et al. [12] proposed a concrete rocking frame derived from a precast unbonded prestressed concrete frame.Beams can be connected to columns through tendons, and the disengagement of the beam from the interface for the beam-column joint can be allowed.A cyclic loading test of beam-column joints for precast PCRF was carried out by Cheok et al. [13] and the test results demonstrated that the joints had great ductility.Priestley et al. [14] completed an experimental study on the seismic behavior of joints of an unbonded PCRF.The failure of joints did not occur in the experiment and the damage of tested joints was minor, which can be quickly repaired after an earthquake.Nonlinear static pushover analysis and dynamic time history analysis of an unbonded post-tensioned PCRF were conducted by El Sheikh et al. [15].The fiber model and spring model used for such structural analysis were also presented.The analysis results showed that the strength, ductility and self-centering ability of an unbonded post-tensioned PCRF can meet the requirements for resisting rare earthquake.
The rocking response and stability of a prestressed rocking frame was studied theoretically by Makris et al. [16].It was found that the effects of prestress on rocking columns with different sizes are basically the same as that of a single independent prestressed rigid rocking column.The seismic response of a prestressed rocking frame with buckling restrained braces which has flag-shaped hysteretic behavior was analyzed by Giouvanidis et al. [17].They revealed that prestress is not always beneficial to the seismic performance of the prestressed rocking frame, but this also depends on the size of the frame columns.Lu et al. [18] proposed a kind of controlled rocking concrete frame.A reversed cyclic loading test and shaking table test were accomplished and the corresponding finite element model was built in the study.The research results proved that the frame has excellent seismic performance.A self-centering prestressed concrete frame with web friction devices was devised by Guo et al. [19,20].The result of a low-cycle reversed loading test implied that the selfcentering prestressed concrete frame has prominent seismic performance and remarkable self-centering ability.
Existing studies show that the PCRF has, in general, superior seismic performance and stability.Nevertheless, the theoretical analysis methods of PCRFs in previous studies are rare and lack depth.Hence, the theoretical model based on rigid body of a one-story singlespan PCRF is established in this paper.The model consists of a rigid beam, rigid columns and elastic prestressing tendons.The tendons are arranged along the centerline of the beam and columns.The connections between the beam and columns and the anchoring of columns are accomplished by prestressing the tendons.The governing equation of motion of the PCRF is derived based on the model and a numerical solution of the equation of motion is obtained.The energy dissipation of the PCRF is analyzed and the calculation method for the coefficient of restitution is proposed.The time history and parameter analysis of the seismic response of the PCRF are presented.

Rigid Body Model of the PCRF
The PCRF has the connecting interfaces set at the column feet and at the vertical inner faces of the beam-column joints, allowing the columns to be uplifted with the accompanying separation of the beam-column interface and rotation of the beam and column around the interface.The theoretical model based on rigid body of a one-story single-span PCRF is shown in Figure 1, which consists of two columns, a beam and prestressing tendons.Relative slips between the rigid members on each interface are not allowed.The tendons are arranged along the centerline of the beam and columns.The connections between the beam and columns and the anchoring of columns are accomplished by prestressing the tendons.The column has a mass m c , height 2h and width 2b.Its semi-diagonal and slenderness ratios are R c =   The motion law of the PCRF is the same as that for the planar four-bar linkage mechanism shown in Figures 2-4 when the PCRF is rocking.The rigid body model of the PCRF is a single degree of freedom system and can be described by a generalized coordinate.The rocking mechanism and planar four-bar linkage mechanism of the PCRF during counterclockwise rotation are shown in Figure 2. The dotted lines in Figures 1 and 2 are used to illustrate the four-bar mechanism presented in Figures 3 and 4. Hence, the planar fourbar linkage model can be used to obtain the relationship between the angles of the beam and columns.That is, the motion of the columns and beam when the PCRF is rocking can be represented by the motion of the bars in the four-bar linkage model.The governing equation of motion of the PCRF can be derived by the Lagrangian equation method.The angles of rotation of bars AB, BC and CD are denoted as 1 q , 2 q and 3 q , respectively.
Figure 4 illustrates the planar four-bar model for the case 1 0 q < and the case 1 0 q > .
The angle 2 q is positive and 3 q is negative in the case 1 0 q < , while the opposite is positive for the case 1 0 q > .The motion law of the PCRF is the same as that for the planar four-bar linkage mechanism shown in Figures 2-4 when the PCRF is rocking.The rigid body model of the PCRF is a single degree of freedom system and can be described by a generalized coordinate.The rocking mechanism and planar four-bar linkage mechanism of the PCRF during counterclockwise rotation are shown in Figure 2. The dotted lines in Figures 1 and 2 are used to illustrate the four-bar mechanism presented in Figures 3 and 4. Hence, the planar four-bar linkage model can be used to obtain the relationship between the angles of the beam and columns.That is, the motion of the columns and beam when the PCRF is rocking can be represented by the motion of the bars in the four-bar linkage model.The governing equation of motion of the PCRF can be derived by the Lagrangian equation method.The angles of rotation of bars AB, BC and CD are denoted as θ 1 , θ 2 and θ 3 , respectively.Figure 4 illustrates the planar four-bar model for the case θ 1 < 0 and the case θ 1 > 0. The angle θ 2 is positive and θ 3 is negative in the case θ 1 < 0, while the opposite is positive for the case θ 1 > 0.  The motion law of the PCRF is the same as that for the planar four-bar linkage mechanism shown in Figures 2-4 when the PCRF is rocking.The rigid body model of the PCRF is a single degree of freedom system and can be described by a generalized coordinate.The rocking mechanism and planar four-bar linkage mechanism of the PCRF during counterclockwise rotation are shown in Figure 2. The dotted lines in Figures 1 and 2 are used to illustrate the four-bar mechanism presented in Figures 3 and 4. Hence, the planar fourbar linkage model can be used to obtain the relationship between the angles of the beam and columns.That is, the motion of the columns and beam when the PCRF is rocking can be represented by the motion of the bars in the four-bar linkage model.The governing equation of motion of the PCRF can be derived by the Lagrangian equation method.The angles of rotation of bars AB, BC and CD are denoted as 1 q , 2 q and 3 q , respectively.
Figure 4 illustrates the planar four-bar model for the case 1 0 q < and the case 1 0 q > .
The angle 2 q is positive and 3 q is negative in the case 1 0 q < , while the opposite is positive for the case 1 0 q > .

Kinematics of the PCRF
The geometrical relationship of the planar four-bar linkage mechanism for 1 0 q < is shown in Figure 3. Firstly, 1 q is assumed as a known angle, and i φ i = are calculation auxiliary angles.From the geometry relationships, we have

Kinematics of the PCRF
The geometrical relationship of the planar four-bar linkage mechanism for 1 0 q < is shown in Figure 3. Firstly, 1 q is assumed as a known angle, and i φ i = are calculation auxiliary angles.From the geometry relationships, we have The equations below can be obtained from sine and cosine law: The expressions for 2 q and 3 q can be derived from Equations (1)-( 5).The expressions are written as

Kinematics of the PCRF
The geometrical relationship of the planar four-bar linkage mechanism for θ 1 < 0 is shown in Figure 3. Firstly, θ 1 is assumed as a known angle, and ϕ i (i = 1, 2, 3) are calculation auxiliary angles.From the geometry relationships, we have The equations below can be obtained from sine and cosine law: The expressions for θ 2 and θ 3 can be derived from Equations (1)- (5).The expressions are written as Appl.Sci.2021, 11, 585 5 of 25 Substituting Equation (2) into Equations ( 6) and (7) yields Differentiating Equations ( 8) and ( 9) with respect to time t, the corresponding angular velocities can be obtained and expressed for simplicity as . . .
θ 1 is the angular velocity of bar AB and the upper dot denotes differentiation with respect to time t.The symbol ∂ θ 1 denotes a partial derivative with respect to θ 1 .Similarly, the second partial derivatives of θ 2 and θ 3 with respect to θ 1 can be written as If point A is the coordinate origin, the positions after the deformation of the left column mass center G AB , beam mass center G BC and right column mass center G CD are

Rocking Motion Equation
The governing equation of motion of the PCRF can be derived from the Lagrangian equation: where T and V are the kinetic energy and potential energy of the PCRF.If the mass of the tendons is neglected, the potential energy of the PCRF can be calculated by where V F is the gravitational potential energy of frame members, and V T is the elastic potential energy due to the elongation of tendons.If AD is regarded as the zero potential energy surface, then V F is expressed as During rocking motion, the elastic potential energy V T can be derived as The additional elongation e cL of the tendon at the base of the left column can be obtained from the cosine law and the geometry relationship shown in Figure 5.We have  The kinetic energy of the PCRF can be calculated by Equation (24): where bx v is the velocity along the x-axis of the beam mass center.by v is the velocity along the y-axis of the beam mass center.Substituting Equation (25) into Equation (24) yields the kinetic energy of the PCRF: During an admissible rotation 1 δθ , Equation ( 27) can be derived by the principle of virtual work The virtual work caused by the external field forces is where ( ) Then the expression of e cL can be written as Similarly, the additional elongation e cR of the tendon at the right column base is The additional elongation e b of the beam tendon can be derived from the cosine law and the geometry relationships shown in Figure 2: The expression of e b can be rewritten as For simplicity, the two column tendons have the same parameters.Assuming P c0 is the initial pretension force of the column tendons, the initial elongation e c0 of the column tendon is where k c = EA cs /2h is the stiffness of the column tendon, E is the Young's modulus of the tendon and A cs is the cross-section area of the column tendons.Assuming P b0 is the initial pretension force of the beam tendon, the initial elongation e b0 of the beam tendon is where k b = EA bs /(2L + 2b) is stiffness of the beam tendon and A bs is the cross-section area of the beam tendon.
The kinetic energy of the PCRF can be calculated by Equation (24): where v b is the velocity of the beam mass center.
is the mass moment of inertia of the column with respect to the pivot point A or D.
is the mass moment of inertia of the beam with respect to its mass center.The expression of v b can derived from Equation (25): where v bx is the velocity along the x-axis of the beam mass center.v by is the velocity along the y-axis of the beam mass center.Substituting Equation (25) into Equation (24) yields the kinetic energy of the PCRF: During an admissible rotation δθ 1 , Equation ( 27) can be derived by the principle of virtual work The virtual work caused by the external field forces is where Substituting Equations ( 28) and (29) into Equation (27) yields the generalized force Q: The substitution of Equations ( 15), ( 16), ( 26) and (30) into Equation ( 13) results in the motion equation of the PCRF: The terms I, J, G, K and S in Equation ( 31) are nonlinear functions of the generalized coordinate θ 1 and are expressed as Equation (32): By considering the inherent symmetry of the PCRF, the motion equation of the PCRF for θ 1 > 0 can be written in the same form as Equation (31), while the expressions of I, J, G, K and S in this case are instead The simplified and conservative starting condition for the calculation program can be obtained by ignoring the influence of the prestressing force on the triggering-off condition of the rocking motion for the PCRF.With this assumption, the lower bound solution of the minimum ground acceleration .. u g,min that initiates the rocking motion of the frame can be derived by substituting the initial conditions ..
The calculation of λ g can be obtained based on the expression of G and S in Equation (32) or Equation (33):

Calculation of Collision Energy Dissipation
It is assumed that energy dissipation occurs only during the collisions when the PCRF is rocking, that is, there is no energy consumption in the process of rotation.Considering that the duration of impact is extremely short, it can be assumed that the PCRF is always in a horizontal state during impacts.Collisions occur at the column feet and the inner faces at the beam-column joint in the case of a one-story single-span PCRF.Regardless of the plastic deformation of the tendons, the elongation of the tendons at the moment of collision is the initial elongation.The forces in the tendons are constant during the collision, which will not affect the momentum moments of the PCRF before and after the collision.The following assumptions are introduced to simplify the analysis and calculation of the collision energy dissipation of the PCRF.It is assumed that the collision forces on each collision surface concentrate on one point which is the postimpact rotation point of the identical rigid body.Relative slips between rigid members on each collision surface are not allowed.Non-collision forces such as gravity and the inertial force of members can be ignored compared to the collision force.
The case θ 1 < 0 is taken as an example for analysis.The collision analysis model of the PCRF is shown in Figure 6 and the frame has returned to a zero-rotation state at this time.It is assumed that the angular velocity where F ci is the impact force and ∆t i is the duration of impact.The ratio of angular velocity after and before the collision is defined as the coefficient of restitution r = . • Linear momentum conservation along the x-axis for the whole frame • Linear momentum conservation along the y-axis for the whole frame ( ) ( ) • Conservation of moment of momentum about point A′ for the whole frame ( )

LI I mh mb Lbm I m h h h h mbL I mh mb Lbm I m h h h h mbL
• Conservation of moment of momentum about point B′ for the left column • Conservation of moment of momentum about point C′ for the right column where is the moment of inertia of the column with respect to its mass center.The calculated expression for the coefficient of restitution r is shown as where • Linear momentum conservation along the x-axis for the whole frame • Linear momentum conservation along the y-axis for the whole frame • Conservation of moment of momentum about point A for the whole frame • Conservation of moment of momentum about point B for the left column • Conservation of moment of momentum about point C for the right column where is the moment of inertia of the column with respect to its mass center.The calculated expression for the coefficient of restitution r is shown as where The dimensionless parameters in Equation ( 43) are the mass ratio of beam to column κ = m b /m c , span-height ratio λ = L/h and height coefficient λ h = h b /h for the beam and λ H = h /h for the column.

Time History Analysis of Seismic Response of the PCRF
The seismic response analysis of the PCRF is based on a basic model.The sections of the beam and columns are square.The frame is in the zero-rotation state without initial angular velocity and angular acceleration.The parameters of the basic model are shown in Table 1.The term ρ c is the density of concrete, and the fracture elongation of the tendons is taken as 1%.The set of earthquake motion records recommended by the ATC-63 project [21] is selected as the source of earthquake records.The information of the earthquakes is shown in Table 2. Records ER1 to ER11 are far-field earthquakes, records ER12 to ER18 are pulse near-field earthquakes and the rest are no-pulse near-field earthquakes.The values of the peak ground acceleration (PGA) and peak ground velocity (PGV) of each earthquake record are also shown in Table 2.According to the calculation results, the time t RI of the rocking motion triggering-off, the time t EP of peak ground acceleration (PGA), the maximum rotation θ 1,max and its occurrence time t RM , as well as the maximum angular velocity .θ 1,max and its occurrence time t VM , can be obtained as shown in Table 3, where the maximum values of rotation and angular velocity are given in absolute values without the consideration of direction.The result of ER15 is illustrated here to explain the characteristics of the seismic response of the PCRF.The time history curve of ER15 is shown in Figure 7 and the calculation results are presented in Figure 8.As shown in Figure 8a,b, t EP does not exactly coincide with t RM and t VM .After the appearance of PGA and local maximum values of acceleration, the amplitudes of responses decay rapidly.At a time around 10.5 s, the frame motion approached the high-frequency vibration stage, in which the vibration amplitude is small and the peak of angular velocity is basically stable.It can be observed from Figure 8c,d that the phase orbit finally formed a spindle-shaped limit cycle which was centered on the origin in the phase diagram.The energy dissipation in the collisions of the PCRF is continuously supplied by the earthquake.Although the frame vibration in the high-frequency vibration stage is similar to steady-state vibration, it is named pseudo-steady-state vibration, considering that its external excitation is a random earthquake.It can be seen that the influence of external excitation on the frame motion is weakened significantly during this pseudo-steady-state vibration.In order to investigate the relationship between the key time points RI  In order to investigate the relationship between the key time points t RI , t EP , t RM and t VM of the basic model under different earthquakes, the t RI of each response are all set as 1 s.The corresponding unified key time points can be obtained.The unified key time points of far-field and near-field earthquake records are shown in Figures 9 and 10, respectively.
Figure 9 illustrates that, under the far-field earthquakes, t RM is very close to t VM and they generally occur shortly after t RI or near t EP .In most of the cases that t RM and t VM appeared shortly after t RI , the frame vibration has entered the pseudo-steady-state vibration stage or a stage with very small amplitudes of responses before t EP .The amplification effect of PGA on frame response is no longer significant in this situation.When the time interval between t RI and t EP is small, the peak values of frame response appear around t EP , or the t RM , and t VM are in close proximity to t RI and t EP , respectively.Figure 10 suggests that the difference between t RM and t VM is very small when the earthquakes are near-field earthquakes and they appear around t RI or t EP simultaneously.It can be found that the relationships between t RI , t EP , t RM and t VM are strongly dependent on the specific earthquake.The variety of earthquakes should be considered in the analysis of the PCRF.
The maximum elongations of prestressed tendons can be calculated based on the maximum rotation of the PCRF, as shown in Figure 11.It can be observed from Figure 11 that the maximum elongations of prestressed tendons are less than 0.5%, which have a large surplus from the 1% value of the fracture elongation of prestressed tendons.That is to say, the reserve of deformation capacity is sufficient and the seismic performance of the PCRF is prominent.It is worth emphasizing that the member connections of the PCRF proposed in this study depend mainly on the prestressed tendons and prestressing forces.Therefore, the structure failure of the PCRF will be led by the fracture of the prestressed tendons.Special attention should be given to the calculation of the deformation and elongation of prestressed tendons in analysis and design.amplification effect of PGA on frame response is no longer significant in this situa When the time interval between t and t is small, the peak values of frame resp Time interval (s) large surplus from the 1% value of the fracture elongation of prestressed tendons.T to say, the reserve of deformation capacity is sufficient and the seismic performance PCRF is prominent.It is worth emphasizing that the member connections of the proposed in this study depend mainly on the prestressed tendons and prestressing Therefore, the structure failure of the PCRF will be led by the fracture of the prest tendons.Special attention should be given to the calculation of the deformation and gation of prestressed tendons in analysis and design.

Parameter Analysis of PCRF Seismic Response
Eight independent dimensionless parameters can be proposed naturally duri derivation of the motion equation of the PCRF, which are the mass ratio of beam to co κ, span-height ratio λ , column aspect ratio tanα , height coefficient h λ , initia stressing force of beam tendon .quake record ER15 is used as a test record and a parameter analysis of the PCRF s response is conducted.
A frame with a larger κ has a smaller response in the early stage of vibrati shown in Figure 12.However, the responses of the frame are changed dramaticall the emergence of a large acceleration pulse at 2~3 s.In the meantime, the frame w larger κ has a significantly larger angle and angular velocity compared wi smaller ones.A number of acceleration pulses appeared around 4 s, which increa responses of the frame with

Parameter Analysis of PCRF Seismic Response
Eight independent dimensionless parameters can be proposed naturally during the derivation of the motion equation of the PCRF, which are the mass ratio of beam to column κ, span-height ratio λ, column aspect ratio tan α, height coefficient λ h , initial prestressing force of beam tendon Π σb = σ pb0 A bs /m c g, initial prestressing force of column tendon Π σc = σ pc0 A cs /m c g, linear stiffness of beam tendon Π Eb = E p A bs /m c g and linear stiffness of column tendon Π Ec = E p A cs /m c g.These are kept because any reduction of independent dimensionless parameters will inevitably influence the precision of the analysis in consideration of the nonlinearity of the rocking motion for the PCRF.The values of the dimensionless parameters for the basic frame model are λ = 2, tan α = 0.2, λ h = 0.2, κ = 1.8, Π σb = 4.13, Π σc = 6.61,Π Eb = 2684.1 and Π Ec = 4294.6.Earthquake record ER15 is used as a test record and a parameter analysis of the PCRF seismic response is conducted.
A frame with a larger κ has a smaller response in the early stage of vibration, as shown in Figure 12.However, the responses of the frame are changed dramatically with the emergence of a large acceleration pulse at 2~3 s.In the meantime, the frame with a larger κ has a significantly larger angle and angular velocity compared with the smaller ones.A number of acceleration pulses appeared around 4 s, which increase the responses of the frame with κ = 1.8 further.The other frames have a smaller rotation before encountering multiple acceleration pulses and the amplification effect of multiple acceleration pulses on their response is limited.Different frames with different κ got into the pseudo-steady-state vibration eventually, while the frames with smaller κ entered earlier.
The excessive increase in the response for the frame with κ = 1.8 caused by the first large acceleration pulse is counterintuitive.The mechanism of this phenomenon can be explained by the frame with larger κ having greater inertia force and energy after the acceleration pulse, causing a sudden increase in its response.The subsequent multiple acceleration pulses made the system response increase further.This phenomenon can be considered to suggest that the acceleration pulse has amplitude sensitivity to the amplification effect of the response of the PCRF.That is, when the vibration amplitude of the frame is large enough or greater than a threshold, the acceleration pulse has an amplification effect on the system response similar to resonance.The magnification of the system response increases with the enlargement of the system amplitude.
For further study of the influence of κ and the influence of earthquake characteristics, the far-field earthquake record of the Northridge earthquake in 1994, named ER2, is selected for analysis.The time history curve of ER2 and similar calculation results are illustrated in Figure 13.It can be seen that the frame with larger κ had a smaller vibration response in the whole vibration process and entered into pseudo-steady-state vibration earlier.Based on the above analysis, cases of different ground motion types should be fully considered in the analysis of PCRFs and time history analysis of various earthquake records is necessary.
It can be found in Figure 14 that the responses of frames with different λ are substantially close at the initial stage of vibration and the differences between responses are not changed significantly with the increase in vibration duration.When multiple consecutive acceleration pulses occurred, the response of the frame with smaller λ is obviously reduced, but the response of the frame with larger λ is increased significantly.Frames with different λ have different amplitudes when multiple continuous acceleration pulses are input.The accumulation of response amplification caused by acceleration pulses results in a large difference in subsequent calculation results.Parameter κ is a constant when λ is variational so the inertia forces of beam and columns are consistent in the analysis for λ.The variation range of λ is not very large in the research.It is for the above reasons that the influence of λ on the dynamic response of the PCRF is not similar to κ.The excessive increase in the response for the frame with 1.8   caused by th large acceleration pulse is counterintuitive.The mechanism of this phenomenon c explained by the frame with larger  having greater inertia force and energy aft acceleration pulse, causing a sudden increase in its response.The subsequent mu acceleration pulses made the system response increase further.This phenomenon c considered to suggest that the acceleration pulse has amplitude sensitivity t amplification effect of the response of the PCRF.That is, when the vibration amplitu the frame is large enough or greater than a threshold, the acceleration pulse h should be fully considered in the analysis of PCRFs and time history analysis of variou earthquake records is necessary.As shown in Figure 15, the frames with different tan have different valu motion triggering-off time, which signifies that the frame with a larger column aspect is more likely to be excited.A frame with a smaller column aspect ratio tends to sto faster speed.After the occurrence of multiple continuous acceleration pulses, the resp of the frame with tan 0.175   has a very obvious amplification effect.The respon As shown in Figure 15, the frames with different tan α have different values of motion triggering-off time, which signifies that the frame with a larger column aspect ratio is more likely to be excited.A frame with a smaller column aspect ratio tends to stop at a faster speed.After the occurrence of multiple continuous acceleration pulses, the response of the frame with tan α = 0.175 has a very obvious amplification effect.The response of the frame with tan α = 0.225, which entered the state of pseudo-steady-state vibration, hardly fluctuates.Figure 16 illustrates that the frame with a larger h  tends to become motio faster than the smaller h  ones.When the first acceleration pulse is imported, it c clearly seen that the response amplification of the frame with small amplitu apparently weaker than that of the other frames.The acceleration pulse basically h response amplification effect on the frame which is in the state of pseudo-steady vibration before the appearance of the acceleration pulse.The response of the frame  Figure 16 illustrates that the frame with a larger λ h tends to become motionless faster than the smaller λ h ones.When the first acceleration pulse is imported, it can be clearly seen that the response amplification of the frame with small amplitude is apparently weaker than that of the other frames.The acceleration pulse basically has no response amplification effect on the frame which is in the state of pseudo-steady-state vibration before the appearance of the acceleration pulse.The response of the frame with λ h = 0.2 has a very obvious amplification phenomenon.It can be seen from the time history curves in Figures 17-20 that the influ dimensionless parameters related to prestressed tendons and initial prestressing fo the frame responses are relatively consistent.The increase in dimensionless para can reduce the historical maximum responses of the frame markedly, which also the earlier appearance of pseudo-steady-state vibration.The difference is t influences on the vibration response in the initial stage are miscellaneous.A frame large initial prestressing force has a relatively large response at the initial s vibration, while the cases of the linear stiffness of prestressed tendons are the opp It can be seen from the time history curves in Figures 17-20 that the influence of dimensionless parameters related to prestressed tendons and initial prestressing forces on the frame responses are relatively consistent.The increase in dimensionless parameters can reduce the historical maximum responses of the frame markedly, which also leads to the earlier appearance of pseudo-steady-state vibration.The difference is that the influences on the vibration response in the initial stage are miscellaneous.A frame with a large initial prestressing force has a relatively large response at the initial stage of vibration, while the cases of the linear stiffness of prestressed tendons are the opposite.

Conclusions
The theoretical model based on rigid body of a one-story single-span P established.The corresponding governing motion equation of the PCRF is derive on the model and the numerical solution of the motion equation is obtained.The dissipation of the PCRF is analyzed and the calculation method for the coeffi restitution is proposed.Time history analyses of the seismic response of the P carried out with different types of earthquake records.The analysis results show maximum rotation of the frame is small and the capacity of structural deform sufficient, which means the PCRF has promising seismic performance.The connections of the PCRF proposed in this study depend on the prestressed tend prestressing forces.The structure failure of the PCRF will be led by the fra prestressed tendons so the calculation of the deformation and elongation of pre tendons should be cautious in analysis and design.
A pseudo-steady-state vibration of the PCRF may appear after some dur vibration under earthquake excitation, in which case a steady spindle-shaped lim is finally formed by phase orbit in the phase diagram.The influence of external ex on the vibration of the PCRF is weakened significantly in the state of pseudo-stea vibration.The parameter analysis of the PCRF is conducted based on eight dimen parameters.The results reveal the characteristic influences of the dimen

Conclusions
The theoretical model based on rigid body of a one-story single-span PCRF is established.The corresponding governing motion equation of the PCRF is derived based on the model and the numerical solution of the motion equation is obtained.The energy dissipation of the PCRF is analyzed and the calculation method for the coefficient of restitution is proposed.Time history analyses of the seismic response of the PCRF are carried out with different types of earthquake records.The analysis results show that the maximum rotation of the frame is small and the capacity of structural deformation is sufficient, which means the PCRF has promising seismic performance.The member connections of the PCRF proposed in this study depend on the prestressed tendons and prestressing forces.The structure failure of the PCRF will be led by the fracture of prestressed tendons so the calculation of the deformation and elongation of prestressed tendons should be cautious in analysis and design.
A pseudo-steady-state vibration of the PCRF may appear after some duration of vibration under earthquake excitation, in which case a steady spindle-shaped limit cycle is finally formed by phase orbit in the phase diagram.The influence of external excitation on the vibration of the PCRF is weakened significantly in the state of pseudo-steady-state vibration.The parameter analysis of the PCRF is conducted based on eight dimensionless parameters.The results reveal the characteristic influences of the dimensionless parameters on the seismic response of the PCRF.The parameters related to prestressed tendons and initial prestressing forces have a consistent influence on the seismic response of the PCRF, and maximum values of seismic response decrease with the increase in the respective parameters.The influences of dimensionless parameters of mass and geometry on the seismic response of the PCRF are complicated, in which case the type of earthquake also needs to be considered.It is worth noting that when the rotational earthquake component and the translational earthquake component are considered simultaneously in the analysis, the motion of the rocking frame will be a three-dimensional motion.The influence of the rotational earthquake component on the seismic response of the PCRF needs further research.

√ h 2 +
b 2 and b/h = tan α, respectively.The beam has a mass m b , height 2h b and length 2L b .Its half-diagonal length is R b = h 2 b + L 2 b and its height span ratio is h b /L b = tan γ.The inclined angle between line AB and AD is β.The coordinate system and other geometric parameters are shown in Figure 1 and the clockwise rotation of the beam and column are positive.The mass centers of the columns and beam are G AB , G CD and G BC .
line AB and AD is β .The coordinate system and other geometric parameters are shown in Figure1and the clockwise rotation of the beam and column are positive.The mass centers of the columns and beam are AB G , C D G and BC G .

Figure 2 .
Figure 2. Rocking mechanism and planar four-bar linkage mechanism of PCRF during counterclockwise rotation.
Appl.Sci.2020, 10, x FOR PEER REVIEW 3 of 26 line AB and AD is β .The coordinate system and other geometric parameters are shown in Figure 1 and the clockwise rotation of the beam and column are positive.The mass centers of the columns and beam are AB G , C D G and BC G .

Figure 2 .
Figure 2. Rocking mechanism and planar four-bar linkage mechanism of PCRF during counterclockwise rotation.Figure 2. Rocking mechanism and planar four-bar linkage mechanism of PCRF during counterclockwise rotation.

Figure 2 .
Figure 2. Rocking mechanism and planar four-bar linkage mechanism of PCRF during counterclockwise rotation.Figure 2. Rocking mechanism and planar four-bar linkage mechanism of PCRF during counterclockwise rotation.

)
Appl.Sci.2020, 10, x FOR PEER REVIEW 7 the beam tendon and bs A is the cross-section area of the beam tendon.

Figure 5 .
Figure 5.The additional elongation cL e of the tendon at the base of the left column.

v
is the velocity of the beam mass center.mass moment of inertia of the column with respect to the pivot point A or D. moment of inertia of the beam with respect to its mass center.The expression of b v can derived from Equation (25):

Figure 5 .
Figure 5.The additional elongation e cL of the tendon at the base of the left column.

θ − 1 1
just before the collision is known of the frame column represented by bar AB.There are five unknowns: impulse I Ax and I Ay at collision point A', I Dx and I Dy at collision point D' and the postimpact angular velocity .θ + Appl.Sci.2021, 11, 585 9 of 25 of the frame column represented by bar C'D'.The impulse I i at an arbitrary collision point i is defined as

Figure 8 .
Figure 8. Calculation results of basic frame model under ER15: (a) time history of rocking rotation; (b) time history of angular velocity; (c) phase orbit; (d) limit cycle.

15 Figure 9 .
Figure 9. Unified key time points in the case of far-field earthquakes.

Figure 10 .
Figure 10.Unified key time points in the case of near-field earthquakes.

Figure 9
Figure 9 illustrates that, under the far-field earthquakes, RM

Figure 9 . 15 Figure 9 .
Figure 9. Unified key time points in the case of far-field earthquakes.

Figure 10 .
Figure 10.Unified key time points in the case of near-field earthquakes.

Figure 9
Figure 9 illustrates that, under the far-field earthquakes, RM t is very close to

Figure 10 .
Figure 10.Unified key time points in the case of near-field earthquakes.

. 2 λ
These are kept because any red of independent dimensionless parameters will inevitably influence the precision analysis in consideration of the nonlinearity of the rocking motion for the PCRF.Th ues of the dimensionless parameters for the basic frame model are

Figure 12 .
Figure 12.Effect of variation in  : (a) time history of rocking rotation; (b) time history of angular velocity.

Figure 12 .
Figure 12.Effect of variation in κ: (a) time history of rocking rotation; (b) time history of angular velocity.

Figure 13 .
Figure 13.Effect of variation in  for the case of ER2: (a) time history of earthquake record ER2; (b) time history of rocking rotation; (c) time history of angular velocity.Figure 13.Effect of variation in κ for the case of ER2: (a) time history of earthquake record ER2; (b) time history of rocking rotation; (c) time history of angular velocity.

Figure 13 .Figure 14 .
Figure 14.Effect of variation in λ: (a) time history of rocking rotation; (b) time history of angular velocity.

Figure 14 .
Figure 14.Effect of variation in λ: (a) time history of rocking rotation; (b) time history of angular velocity.

Figure 15 .
Figure 15.Effect of variation in tanα: (a) time history of rocking rotation; (b) time history of angular velocity.

Figure 15 .
Figure 15.Effect of variation in tanα: (a) time history of rocking rotation; (b) time history of angular velocity.

Figure 16 .
Figure 16.Effect of variation in λh: (a) time history of rocking rotation; (b) time history of angular velocity.

Figure 16 .
Figure 16.Effect of variation in λ h : (a) time history of rocking rotation; (b) time history of angular velocity.

Figure 17 .
Figure 17.Effect of variation in σb Π : (a) time history of rocking rotation; (b) time history of angular velocity.

Figure 19 .
Effect of variation in Π Eb : (a) time history of rocking rotation; (b) time history of angular velocity.

Figure 20 .
Figure 20.Effect of variation in

Figure 20 .
Figure 20.Effect of variation in Π Ec : (a) time history of rocking rotation; (b) time history of angular velocity.

Table 1 .
Parameters of the basic model.

Table 3 .
Time history analysis results of seismic response.