Skip to main content

ORIGINAL RESEARCH article

Front. Built Environ., 25 August 2021
Sec. Earthquake Engineering
This article is part of the Research Topic Energy-Based Seismic Engineering View all 8 articles

Explicit Overturning Limit of Rigid Block Using Triple and Pseudo-Triple Impulses Under Critical Near-Fault Ground Motions

Sae HommaSae Homma1Kunihiko NabeshimaKunihiko Nabeshima2Izuru Takewaki
Izuru Takewaki1*
  • 1Department of Architecture and Architectural Engineering, Graduate School of Engineering, Kyoto University, Kyoto, Japan
  • 2Graduate School of Advanced Science and Engineering, Hiroshima University, Hiroshima, Japan

An explicit limit for the overturning of a rigid block is derived on the input level of the triple impulse and the pseudo-triple impulse as a modified version of the triple impulse that are a substitute of a near-fault forward-directivity ground motion. The overturning behavior of the rigid block is described by applying the conservation law of angular momentum and the conservation law of mechanical energy (kinetic and potential). The initial velocity of rotation after the first impulse and the change of rotational velocity after the impact on the floor due to the movement of the rotational center are determined by using the conservation law of angular momentum. The maximum angle of rotation after the first impulse is obtained by the conservation law of mechanical energy. The change of rotational velocity after the second impulse is also characterized by the conservation law of angular momentum. The maximum angle of rotation of the rigid block after the second impulse, which is mandatory for the computation of the overturning limit, is also derived by the conservation law of mechanical energy. This allows us to prevent from computing complex non-linear time-history responses. The critical timing of the second impulse (also the third impulse timing to the second impulse) is featured by the time of impact after the first impulse. As in the case of the double impulse, the action of the second impulse just after the impact is employed as the critical timing. It is induced from the explicit expression of the critical velocity amplitude limit of the pseudo-triple impulse that its limit is slightly larger than the limit for the double impulse. Finally, it is found that, when the third impulse in the triple impulse is taken into account, the limit input velocity for the overturning of the rigid block becomes larger than that for the pseudo-triple impulse. This is because the third impulse is thought to prevent the overturning of the rigid block by giving an impact toward the inverse direction of the vibration of the rigid block.

Introduction

The rocking behavior of rigid blocks plays an important role in the evaluation of earthquake response of slender structures, monuments and furniture. Since many base-isolated high-rise buildings are being constructed in Japan and most isolators do not have tensile resistance, the seismic risk analysis of such slender structures is extremely focused recently.

The seismic rocking response of rigid blocks has been investigated for long time. The work was initiated by Milne (1885) and the detailed fundamental formulation was made by Housner (1963). Subsequently, Yim et al. (1980) did extensive investigation for versatile recorded ground motions following Housner (1963) and Ishiyama (1982) classified the seismic rocking responses into various types of nonlinear overturning response. Many investigations were conducted continuously so far (Spanos and Koh 1984, Hogan 1989, Shenton III and Jones 1991, Pompei et al., 1998, Zhang and Makris 2001, Prieto et al., 2004, Yilmaz et al., 2009, ElGawady et al., 2010, DeJong 2012, Dimitrakopoulos and DeJong 2012a; Dimitrakopoulos and DeJong 2012b). Rather recently, DeJong (2012) and Dimitrakopoulos and DeJong (2012a), Dimitrakopoulos and DeJong (2012b) investigated the rocking response of a rigid block in a more detailed manner. Casapulla et al. (2010) and Casapulla and Maione (2016) considered the multiple sequence of impulses for the rocking response of a rigid block and found the resonant response. They investigated the responses under two-type multiple impulses, i.e., with gradually increasing intervals for resonance and with equal intervals. Furthermore, Makris and Kampas (2016) derived an important fact related to the scale effect of blocks on the overturning limit level of sinusoidal inputs and earthquake ground motions.

It is well known that slender structures are affected strongly by impulsive pulse-type loading. After Northridge earthquake in 1994, Hyogoken-Nanbu earthquake in 1995 and Chi-Chi earthquake in 1999, many earthquake structural engineers and designers focused on the structural design under impulsive pulse-type earthquake ground motions (Hall et al., 1995; Sasani and Bertero 2000; Mavroeidis et al., 2004; Makris and Black 2004; Kalkan and Kunnath 2006; Xu et al., 2007; Rupakhety and Sigbjörnsson 2011; Champion and Liel 2012; Gerami and Sivandi-Pour 2014; Vafaei and Eskandari 2015; Khaloo et al., 2015). The principal part of such ground motions in the form of pulse-type waves represents the essential features of near-fault ground motions and such principal part was found to be classified into two types, i.e., the fling-step (fault-parallel) and forward-directivity (fault-normal) inputs (Mavroeidis and Papageorgiou 2003, Bray and Rodriguez-Marek 2004, Kalkan and Kunnath 2006, Mukhopadhyay and Gupta 2013a; Mukhopadhyay and Gupta 2013b; Zhai et al., 2013, Hayden et al., 2014, Yang and Zhou 2014, Kohrangi et al., 2018, Khansefid and Bakhshi 2019; Khansefid 2020). Afterward, it was made clear that the principal parts of the fling-step and forward-directivity inputs can be represented by some wavelets or a series of harmonic waves. Mavroeidis and Papageorgiou (2003) and some others pointed out the remarkable characteristics of this class of ground motions and provided some simple models. For example, the fling-step and forward-directivity inputs were characterized by two or three wavelets. It was reported that these pulse-type ground motions affect greatly the long-period building structures and large inelastic deformations could be concentrated to lower stories.

Recently Nabeshima et al. (2016) introduced a double impulse for the rocking vibration of a rigid block in order to substitute the response of the rigid block to near-field ground motions. The double impulse was proposed by Kojima and Takewaki (2015a), Kojima and Takewaki (2016), Taniguchi et al. (2016) to express the peculiar characteristics of near-fault ground motions. It was clarified that the double impulse enables the derivation of the closed-form maximum deformation of elastic-plastic SDOF systems. Nabeshima et al. (2016) derived a closed-form solution of the maximum critical response and a closed form expression of the limit input velocity level for the overturning. They also made clear that the critical timing of the second impulse is just after the impact on the ground. Taniguchi et al. (2017) extended the work of Nabeshima et al. (2016) dealing with the critical case (critical timing of the second impulse which gives the maximum response of a rigid block) to the non-critical case which is important in view of actual situations.

In the present paper, the triple impulse as a substitute of the forward-directivity input of near-fault ground motions and its modified version called the pseudo-triple impulse are treated and an explicit expression is obtained for the critical velocity amplitude limit of the triple impulse and the pseudo-triple impulse for overturning of a rigid block. The formulation for the pseudo-triple impulse is presented first and it is extended to the triple impulse.

Triple Impulse and Corresponding 1.5-Cycle Sinusoidal Wave

In this paper, it is aimed at modeling a principal part of a near-fault ground motion into a 1.5-cycle sinusoidal wave and then simplify such 1.5-cycle sinusoidal wave into a triple impulse following the references (Kojima and Takewaki 2015b) as shown in Figure 1. This is because the triple impulse has a simple characteristic and a straightforward derivation of its maximum response can be expected even for non-linear elastic responses based on an energy approach to free vibrations.

FIGURE 1
www.frontiersin.org

FIGURE 1. Modeling of principal part of near-fault ground motion into 1.5-cycle sinusoidal wave and modeling of such 1.5-cycle sinusoidal wave into triple impulse.

Following the reference (Kojima and Takewaki 2015b), consider a ground acceleration u¨g(t) as a triple impulse, as shown in Figure 1A, expressed by

u¨g(t)=0.5Vδ(t)Vδ(tt0)+0.5Vδ(t2t0)(1)

where 0.5V is the given initial velocity and –V, 0.5V are the second and third impulse velocity amplitude. t0 is the time interval between two consecutive impulses. The time derivative is denoted by an over-dot. The comparison with the corresponding 1.5-cycle sinusoidal wave expressed by Eq. 2 is plotted in Figure 1.

u¨g(t)={0.5Apsin(ωpt)(0t<0.5Tp,Tpt1.5Tp)Apsin(ωpt)(0.5Tpt<Tp)(2)

In Eq. 2, Ap,Tp,ωp=2π/Tp denote the amplitude of the 1.5-cycle sinusoidal acceleration wave, its period and its circular frequency. In this case, the relation Tp=2t0 holds. Due to Takewaki and Kojima (2021), the relations Vp=Ap/ωp and Vp/V=0.62235722... hold.

The velocity and displacement of the triple impulse and the corresponding modified sinusoidal wave can also be found in the reference (Kojima and Takewaki 2015b). It has been confirmed that the triple impulse is a good approximation of the corresponding 1.5-cycle sinusoidal wave even in the form of velocity and displacement.

The Fourier transform of the acceleration u¨g(t) of the triple impulse can be expressed as

U¨g(ω)={0.5Vδ(t)Vδ(tt0)+0.5Vδ(t2t0)}eiωtdt=V(0.5eiωt0+0.5eiω2t0)(3)

While most of the previous methods utilize the equivalent linearization of the oscillating model for the unchanged input (see Figure 2A including an equivalent linear stiffness) or introduce detailed nonlinear structural models for which time-consuming response analysis is necessary, the method proposed in the works (Kojima and Takewaki 2015a; Kojima and Takewaki 2015b; Kojima and Takewaki 2016) and in this paper transforms the input into the triple or pseudo-triple impulse for the unchanged oscillating model (see Figure 2B). It should be pointed out that the rocking response of a rigid block includes a complicated phenomenon, e.g., strong nonlinearity, sudden energy loss at impact, the equivalent linearization may be difficult to apply.

FIGURE 2
www.frontiersin.org

FIGURE 2. Comparison of new method of input transformation with conventional method of oscillating model transformation. (A) Conventional method: equivalent linearization of oscillating model for unchanged input, (B) New method: transformation of input into triple impulse for unchanged oscillating model.

Maximum Rotation Angle of Rigid Block Subjected to Critical Pseudo-Triple Impulse

Consider the rocking vibration of a rigid block (mass: m, width: 2b, height: 2h) subjected to a base acceleration u¨g(t) as shown in Figure 3 (input of u¨g(t) does not change its essence). Slipping of the block is assumed to be ignored. The geometrical properties are denoted by R=b2+h2 and α as shown in Figure 3. I(=(4/3)mR2) and g denote the mass moment of inertia around the edge of bottom right (also bottom left) and the acceleration of gravity, respectively. It was pointed out that this rocking vibration can be represented by the mechanical behavior of a rigid bar, as shown in Figure 3, which has the same mass moment of inertia and is supported by a non-linear elastic rotational spring with rigid initial stiffness and negative second slope. The moment-rotation relation of the non-linear elastic rotational spring is shown in Figure 4.

FIGURE 3
www.frontiersin.org

FIGURE 3. Modeling of rocking rigid block by rigid bar supported by non-linear elastic rotational spring with rigid initial stiffness and negative second slope (Nabeshima et al., 2016).

FIGURE 4
www.frontiersin.org

FIGURE 4. Moment-rotation relation for rocking response of rigid block and timing of impulses (red solid circle is the first impulse and red solid triangle is the second impulse) (Nabeshima et al., 2016).

In this paper, a pseudo-triple impulse as shown in Figure 5A is introduced where the third impulse in the triple impulse as shown in Figure 5B is neglected. This new model is treated because the third impulse may prevent the overturning of the rigid block by giving an impact toward the inverse direction of the vibration of the rigid block and the derivation of the explicit overturning limit for the triple impulse seems to be complicated and difficult. A scenario that the overturning occurs after the second impulse in the pseudo-triple impulse is employed here. This scenario seems valid because the input limit on the overturning corresponding to this scenario provides a lower limit in general. As for more detailed scenarios, see Ishiyama (1982) and Dimitrakopoulos and DeJong (2012b).

FIGURE 5
www.frontiersin.org

FIGURE 5. Input motions in numerical examples, (A) Pseudo-triple impulse, consisting of the first two impulses in triple impulse, and the corresponding equivalent modified 1.0-cycle sinusoidal wave, (B) Triple impulse and the corresponding equivalent 1.5-cycle sinusoidal wave.

Figure 6 shows the overview of the rocking vibration under the pseudo-triple impulse, consisting of the first two impulses in triple impulse. The critical acting timing of the second impulse is at the instant of impact where the velocity of rotation attains the maximum (see Nabeshima et al., 2016). Furthermore, it can be shown that critical timing is just after the impact because the velocity of rotation is reduced greatly at the impact. More detailed verification of the critical timing of the second impulse is shown in Nabeshima et al. (2016).

FIGURE 6
www.frontiersin.org

FIGURE 6. Overview of rocking vibration under pseudo-triple impulse, consisting of the first two impulses in triple impulse, (A) Overview of rocking vibration, (B) Mechanical energy conservation and angular momentum conservation.

When the angle of rotation of the rigid block is denoted by θ(t), the equation of motion can be expressed by

Iθ¨(t)+mgRsin{αθ(t)}=mu¨g(t)Rcos{αθ(t)}θ(t)<0(4a)
Iθ¨(t)+mgRsin{αθ(t)}=mu¨g(t)Rcos{αθ(t)}θ(t)>0(4b)

The block rotates around an edge after the uplift from the ground. During the rotation, the position (potential) energy increases and the block resists the input acceleration.

Let θ1max denote the maximum angle of rotation after the first impulse and θ2maxdenote that after the second impulse. Since the formulation is similar to Nabeshima et al. (2016) except the amplitude (0.5V) of the first impulse, the formulation is explained simply. The flowchart presented in Nabeshima et al. (2016) may be useful.

The first conservation law of angular momentum just after the first impulse and the first conservation law of mechanical energy after the first impulse can provide the initial rotational velocity θ˙1 after the first impulse as θ˙1=3Vcosα/8R and 0.5Iθ˙12=A1 where A1 is the area of the trapezoid-like zone corresponding to θ1max shown in Figure 4. The expression of A1 in terms of θ1max, and use of the linear approximation sin(αθ)(αθ),sinαα lead to θ1maxαα2{3V2cos2α/(16gR)}.

The second conservation law of angular momentum at the impact (θ˙1: rotational velocity just before the impact) provides Iθ˙12mRbsinαθ˙1=Iθ˙2. According to Housner (1963), the rotational velocity θ˙2 just after the impact to the floor may be expressed by θ˙2=rθ˙1 where r={(1+3cos2α)/4}2. This parameter r was defined by Housner (1963). When the surface boundary condition between the block and the base is provided, e.g. the surface material properties (ElGawady et al., 2010) or the rocking of tall buildings (limited contact area), another coefficient could be added on or incorporated into the parameter r.

The conservation law of angular momentum just after the second impulse provides the rotational velocity change Δθ˙ by the second impulse as IΔθ˙=mVRcosα. This leads to Δθ˙=2θ˙1 and θ˙2=θ˙2+Δθ˙=(2+r)θ˙1. The rotation angle θ˙2 is obtained following the Housner’s formulation (conservation law of angular momentum at the impact) and Δθ˙ is derived from the transformation of the horizontal impulse into rotation. The second conservation law of mechanical energy provides 0.5Iθ˙22=A2 where A2 is the area of the trapezoid-like zone corresponding to θ2max shown in Figure 4.

Limit Input Level of Critical Pseudo-Triple Impulse Characterizing Overturning of Rigid Block

The overturning of the rigid block can be characterized by the attainment of θ2max to α as shown in Figure 7. Application of this condition to 0.5Iθ˙22=A2 provides the following limit input level for the pseudo-triple impulse.

Vc1over=2R(1+0.5r)h2(Rh)3g(5)

FIGURE 7
www.frontiersin.org

FIGURE 7. Moment-rotation relation for rocking response of rigid block and limit of overturning (Nabeshima et al., 2016).

The subscript “c1” indicates “critical” for the pseudo-triple impulse corresponding to the “1.0”-cycle sine wave. This limit input velocity is slightly larger than the limit for the double impulse derived in Nabeshima et al. (2016) (the coefficient 0.5 is a new addition). The critical timing Δtcr can be obtained approximately (linear approximation: see the reference (Housner 1963)) by solving the equation of free-rocking motion from the first impulse to the second impulse.

Δtcr2pcosh1(11θ1maxα)(6)

where p=mgR/I=3g/(4R). While the overturning limit input velocity expressed above was derived by the energy balance approach without linearizing the equation of motion, the critical timing expressed above was instead obtained approximately by solving the linearized equation of free rocking motion as in Housner (1963), Nabeshima et al. (2016) and Casapulla (2015).

The limit critical timing Δtc1over characterizing the just overturning of the block after the second impulse under the pseudo-triple impulse can be derived by substituting V=Vc1over into θ1maxαα2{3V2cos2α/(16gR)} and using Eqs. 5, 6 and an approximation cosα=1(1/2)α2.

Δtc1over(2/p)cosh1{(2+r)/(2+r)21}(7)

The above formulation was made for the case where the third impulse was disregarded. Figure 9 shows the overview of the rocking vibration under the triple impulse. The overturning limit for the complete triple impulse will be derived in the later section.

Numerical Examples and Discussion

Pseudo-triple Impulse and 1.0-Cycle Sinusoidal Wave

Numerical examples are introduced here to demonstrate the accuracy and reliability of the proposed method. Three numerical examples of rectangular columns are considered for the models with width (2b) = 1, 2, 4 m corresponding to the references (Makris and Kampas 2016; Nabeshima et al., 2016). The column height is changed parametrically.

Figure 8 shows the comparison of time-history responses θ(t)/α at the overturning limit between the pseudo-triple impulse, consisting of the first two impulses in the triple impulse, and the corresponding equivalent modified 1.0-cycle sinusoidal wave magnified by a coefficient about 1.31–1.32. The “equivalent” means that the maximum Fourier amplitudes of the triple impulse (not the pseudo-triple) and the 1.5-cycle sinusoidal wave have the same value. Vc indicates the value of Eq. 20 and Asin presents the acceleration amplitude of the modified 1.0-cycle sinusoidal wave which gives the overturning limit. Ap presents the acceleration amplitude of the equivalent 1.0-cycle sinusoidal wave transformed using the relation Vp=Ap/ωp=ApTp/(2π). It can be understood that the pseudo-triple impulse, consisting of the first two impulses in triple impulse, attains the overturning limit accurately and the corresponding equivalent modified 1.0-cycle sinusoidal wave magnified by a coefficient about 1.31–1.32 indicates a rather accurate overturning limit.

FIGURE 8
www.frontiersin.org

FIGURE 8. Comparison of time-history responses θ(t)/α at overturning limit between the pseudo-triple impulse, consisting of the first two impulses in triple impulse, and the corresponding equivalent modified 1.0-cycle sinusoidal wave magnified by a coefficient about 1.31–1.32.

Triple Impulse and 1.5-Cycle Sinusoidal Wave

It may be useful to investigate the response of a rigid block under the triple impulse to study the role of the third impulse which was disregarded in the pseudo-triple impulse. Figure 9 shows the overview of rocking vibration under the triple impulse from the viewpoint of mechanical energy conservation and angular momentum conservation.

FIGURE 9
www.frontiersin.org

FIGURE 9. Overview of rocking vibration under triple impulse from the viewpoint of mechanical energy conservation and angular momentum conservation.

Figure 10 shows the comparison of time-history responses θ(t)/α at the overturning limit between the triple impulse and the corresponding equivalent 1.5-cycle sinusoidal wave magnified by a coefficient about 1.38–1.39. Vc indicates the limit velocity amplitude of the triple impulse and Asin presents the acceleration amplitude of the 1.5-cycle sinusoidal wave which gives the overturning limit. It can be understood that the triple impulse attains the overturning limit accurately and the corresponding equivalent 1.5-cycle sinusoidal wave magnified by a coefficient about 1.38–1.39 indicates a rather accurate overturning limit. It can also be found that the third impulse in the triple impulse plays a role of preventing the overturning of the rigid block and larger amplification coefficients are required to overturn the rigid block. This phenomenon will be investigated further in the following part.

FIGURE 10
www.frontiersin.org

FIGURE 10. Comparison of time-history responses θ(t)/α at overturning limit between the triple impulse and the corresponding equivalent 1.5-cycle sinusoidal wave magnified by a coefficient about 1.38–1.39.

Comparison of Overturning Limit Acceleration Amplitude and Overturning Limit Velocity Amplitude

Figure 11 shows the comparison of the overturning limit acceleration amplitude AP among the closed-form limit and the limit simulated by the time-history response analysis (THRA) using the pseudo-triple impulse (consisting of the first two impulses in the triple impulse), the limit simulated by THRA for the triple impulse, the limits simulated by THRA for the equivalent modified 1.0-cycle sinusoidal wave and for the 1.5-cycle sinusoidal wave for 2b = 1, 2, 4 (m). AP for the pseudo-triple impulse and the triple impulse was transformed using the relation Vp=Ap/ωp=ApTp/(2π). It can be found that the simulation by THRA for the pseudo-triple impulse demonstrates the accuracy of the closed-form overturning limit. Furthermore, the limit AP for the equivalent modified 1.0-cycle sinusoidal wave corresponds to the level of the closed-form limit multiplied by about 1.32 and the limit AP for the 1.5-cycle sinusoidal wave is slightly larger than the above-mentioned limit (modified 1.0-cycle sine). It seems that the additional 0.5-cycle sine wave of the 1.5-cycle sine wave prevents the overturning a little bit. In addition, the limit for the triple impulse is slightly smaller than the closed-form limit. This may result from the criticality of the triple impulse (input timing) related to the relation Vp=Ap/ωp=ApTp/(2π). This will be clarified in the following investigation on the velocity limit.

FIGURE 11
www.frontiersin.org

FIGURE 11. Comparison of overturning limit acceleration amplitude AP among closed-form limit using pseudo-triple impulse (consisting of the first two impulses in triple impulse), triple impulse limit, equivalent modified 1.0-cycle sinusoidal wave, 1.5-cycle sinusoidal wave and Rinaldi Station FN.

In Figure 11, the plot of the overturning limit acceleration of Rinaldi Station FN shown in Figure 1 is also presented. Only the amplitude modulation was conducted. Since the resonant phenomenon can be found only around 2[m]<R<3[m] for 2b = 1 (m), other plots are far from “1.5 sine”. It can be observed that the present formulation for the critical resonant case can capture approximately the essential feature of near-field recorded ground motions in the critical case.

Figure 12 presents the comparison of the overturning limit velocity amplitude VP among the closed-form limit and the limit simulated by THRA using the pseudo-triple impulse, the limit simulated by THRA for the triple impulse, the limits simulated by THRA for the equivalent modified 1.0-cycle sinusoidal wave and for the 1.5-cycle sinusoidal wave for 2b = 1, 2, 4 (m). VP for the equivalent modified 1.0-cycle sinusoidal wave and the 1.5-cycle sinusoidal wave was transformed using the relation Vp=Ap/ωp=ApTp/(2π). It can be found that the simulation by THRA for the pseudo-triple impulse demonstrates the accuracy of the closed-form overturning limit as in Figure 11. Furthermore, the limit VP for the equivalent modified 1.0-cycle sinusoidal wave corresponds to the level of the closed-form limit multiplied by about 1.32 and the limit VP for the 1.5-cycle sinusoidal wave is larger than the above-mentioned limit (modified 1.0-cycle sine). It seems that the additional 0.5-cycle sine wave of the 1.5-cycle sine wave prevents the overturning as in Figure 11. In addition, the limit for the triple impulse is larger than the closed-form limit by approximately 15%. This means that the third impulse in the triple impulse plays a role to prevent the overturning compared to the pseudo-triple impulse. The comparisons in terms of AP and VP are useful in understanding the features of the overturning limits by various inputs.

FIGURE 12
www.frontiersin.org

FIGURE 12. Comparison of overturning limit velocity amplitude VP among closed-form limit using pseudo-triple impulse (consisting of the first two impulses in triple impulse), triple impulse limit, equivalent modified 1.0-cycle sinusoidal wave limit and 1.5-cycle sinusoidal wave limit.

It should be remarked that, since the input periods of the 1.0-cycle and 1.5-cycle sinusoidal waves used in Figures 11, 12 are based on Eq. 23 for the pseudo-triple impulse, their criticality is not necessarily clear. If true critical periods of the 1.0-cycle and 1.5-cycle sinusoidal waves are sought and used, the corresponding limit input levels may become slightly smaller ones.

Additional Theoretical Investigation for Triple Impulse

In this section, a closed-form overturning velocity limit for the triple impulse is derived. It is assumed here again that the second impulse acts just after the impact of the rigid block on the floor. In this case, the timing of the third impulse is prescribed.

The equation of free-vibration motion of the rigid block after the second impulse can be described by

Iθ¨+mgRsin(αθ(t))=0    (θ(t)<0)(8)

Using the approximation sin(αθ(t))αθ(t) as in the previous sections, Eq. 8 leads to

Iθ¨+mgR(αθ(t))=0    (θ(t)<0)(9)

The general solution of Eq. 9 can be expressed by

θ(t)=α+C1ep(tΔtcr)+C2ep(tΔtcr)(10)

Using the rotation condition θ(Δtcr)=0,θ˙(Δtcr)=(r+2)θ˙1 at t=Δtcr (the impact of the rigid block on the floor), the constants C1,C2 can be determined as

C1=12((r+2)θ˙1pα)(11)
C2=12((r+2)θ˙1pα)(12)

The rotation during the free vibration after the second impulse can then be expressed by

θ(t)=α12((r+2)θ˙1pα)ep(tΔtcr)+12((r+2)θ˙1pα)ep(tΔtcr)     (Δtcrt<2Δtcr)(13)

Here, θ˙1 and the critical timing Δtcr were obtained in the previous section.

From these conditions, the rotation θ3=θ(2Δtcr) and the rotation velocity θ˙3=θ˙(2Δtcr) just before the third impulse can be expressed as

θ3=θ(2Δtcr)=α12((r+2)θ˙1pα)epΔtcr+12((r+2)θ˙1pα)epΔtcr(14)
θ˙3=θ˙(2Δtcr)=12p((r+2)θ˙1pα)epΔtcr12p((r+2)θ˙1pα)epΔtcr(15)

Let Δθ˙3 denote the change of the rotation velocity at the impact of the third impulse. Then the conservation law of the angular momentum leads to

IΔθ˙3=m(0.5V)Rcos(αθ3)(16)

Therefore, Δθ˙3 can be obtained as

Δθ˙3=3Vcos(αθ3)8R=cos(αθ3)cosαθ˙1(17)

The rotation θ˙3 just after the action of the third impulse can be expressed by

θ˙3=θ˙3+Δθ˙3=12p((r+2)θ˙1pα)epΔtcr12p((r+2)θ˙1pα)epΔtcr+{cos(αθ3)cosα}θ˙1(18)

Let E1 and E2 denote the energy (kinetic plus position energy) just after the action of the third impulse and the position energy at the maximum rotation θ3max after the third impulse. These quantities can be obtained as

E1=12I(θ˙3+Δθ˙3)2+12mgRα2{1(1+θ3α)2}(19)
E2=12mgRα2{1(1+θ3maxα)2}(20)

From the energy conservation law, the following relation holds.

E1=E2(21)

θ3max can be derived from Eqs 1921.

θ3maxα=1+11p2(θ˙3α+Δθ˙3α)2{1(1+θ3α)2}(22)

Substitution of the overturning condition θ3max=α into Eq. 22 provides

1p2(θ˙3α+Δθ˙3α)2=(1+θ3α)2(23)

Substitution of all the derived quantities into Eq. 23 and solution of the resulting equation for the velocity lead to the final form of the overturning limit velocity Vc1_TIover.

Figure13 shows the overturning limit velocity Vc1_TIover using the relation Vp/V=0.62235722....... among the closed-form limit using the pseudo-triple impulse (consisting of the first two impulses in the triple impulse), the triple impulse limit (simulation), the equivalent modified 1.0-cycle sinusoidal wave limit (simulation), the 1.5-cycle sinusoidal wave limit (simulation) and the closed-form limit obtained in this section using the triple impulse. It can be observed that the derived closed-form expression corresponds well to the value obtained through the time-history response analysis simulation conducted in the previous section.

FIGURE 13
www.frontiersin.org

FIGURE 13. Comparison of overturning limit velocity amplitude VP among closed-form limit using pseudo-triple impulse (consisting of the first two impulses in triple impulse), triple impulse limit, equivalent modified 1.0-cycle sinusoidal wave limit, 1.5-cycle sinusoidal wave limit and closed-form limit using triple impulse.

Conclusion

An approximate closed-form limit has been derived on the input level of the pseudo-triple impulse as a representative of the principal part of a near-fault ground motion for the overturning of a rigid block. The conclusions may be summarized as follows.

1) As in the double impulse case, the rocking vibration of a rigid block can be formulated by incorporating the conservation law of angular momentum and the conservation law of mechanical energy. The conservation law of angular momentum enables the determination of the initial rotational velocity just after the first impulse and the rotational velocity changes at the impact of the rigid block on the floor and the impact due to the second impulse. On the other hand, the conservation law of energy enables the determination of the maximum rotational angle after the first impulse and that after the second impulse. These maximum rotational angles are required for the computation of the overturning limit. These two conservation laws are quite effective and enables us to avoid the complicated non-linear time-history response analysis.

2) The critical timing of the second impulse after the action of the first impulse has been characterized by the time of impact of the rigid block on the floor. The action of the second impulse just after the impact corresponds to the critical timing as in the case of the double impulse. Although, in most cases under the triple impulse, the third impulse acts before the overturning occurs, the neglect of the third impulse leads to a closed-form expression of the overturning limit for the pseudo-triple impulse where the third impulse in the triple impulse is disregarded.

3) By introducing the condition that the attainment of the maximum rotational angle θ2max after the second impulse at the limit value α of rotation characterizes the overturning of the rigid block under the pseudo-triple impulse, the critical velocity amplitude of the pseudo-triple impulse just inducing the overturning of the rigid block was derived. The closed-form expression of the area A2 in the restoring-force characteristic in the negative side enabled the derivation of the closed-form velocity amplitude limit of the pseudo-triple impulse for overturning with the aid of the conservation law of energy.

4) The limit input velocity for the pseudo-triple impulse is slightly larger than the limit for the double impulse derived in Nabeshima et al. (2016). This finding was enabled through the closed-form expression of the limit input velocity.

5) Numerical examples including the comparison with the numerical simulation results by the Runge-Kutta method demonstrated the accuracy and reliability of the proposed method. However, as for the comparison of the response to the pseudo-triple impulse with that to the equivalent modified 1.0-cycle sinusoidal wave, a magnification coefficient (about 1.31–1.32 in this case) should be introduced for guaranteeing the correspondence of the responses to the pseudo-triple impulse and to the equivalent modified 1.0-cycle sinusoidal wave.

6) It was demonstrated that, when the third impulse is taken into account, the limit input velocity for the overturning of the rigid block becomes larger. This is because the third impulse may prevent the overturning of the rigid block by giving an impact toward the inverse direction of the vibration of the rigid block. In this case, a magnification coefficient (about 1.38–1.39 in the present case) should be introduced for guaranteeing the correspondence of the responses to the triple impulse and to the equivalent 1.5-cycle sinusoidal wave. The closed-form overturning limit velocity for the triple impulse was also derived by using the conservation law of angular momentum and energy.

7) Since the overturning limit is sensitive to the input frequency (period) of the related sinusoidal waves, the comparisons of the overturning limits in terms of AP and VP are useful in understanding the features of the overturning limits by various inputs.

The realization of the critical case (combination of the critical impulse timing and the critical velocity amplitude) is very rare and it may be difficult to simulate this critical case using actual ground motions without amplitude modulation. The principal purpose of this paper is to derive the explicit overturning limit even for a simple input. The detailed application of the present result to actual ground motions will be investigated in the future.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author Contributions

SH formulated the problem, conducted the computation, and wrote the paper. KN conducted the computation and discussed the results. IT supervised the research and wrote the paper.

Funding

Part of the present work is supported by the Grant-in-Aid for Scientific Research (KAKENHI) of Japan Society for the Promotion of Science (No.18H01584). This support is greatly appreciated.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Bray, J. D., and Rodriguez-Marek, A. (2004). Characterization of Forward-Directivity Ground Motions in the Near-Fault Region. Soil Dyn. Earthquake Eng. 24 (11), 815–828. doi:10.1016/j.soildyn.2004.05.001

CrossRef Full Text | Google Scholar

Casapulla, C., Jossa, P., and Maione, A. (2010). Rocking Motion of a Masonry Rigid Block under Seismic Actions: a New Strategy Based on the Progressive Correction of the Resonance Response. Ingegneria Sismica 27 (4), 35–48.

Google Scholar

Casapulla, C., and Maione, A. (2017). Free Damped Vibrations of Rocking Rigid Blocks as Uniformly Accelerated Motions. Int. J. Str. Stab. Dyn. 17, 1750058. doi:10.1142/S0219455417500584

CrossRef Full Text | Google Scholar

Casapulla, C. (2015). On the Resonance Conditions of Rigid Rocking Blocks. Int. J. Eng. Tech. 7 (2), 760–771.

Google Scholar

Champion, C., and Liel, A. (2012). The Effect of Near-Fault Directivity on Building Seismic Collapse Risk. Earthquake Engng Struct. Dyn. 41 (10), 1391–1409. doi:10.1002/eqe.1188

CrossRef Full Text | Google Scholar

DeJong, M. J. (2012). Amplification of Rocking Due to Horizontal Ground Motion. Earthquake Spectra 28 (4), 1405–1421. doi:10.1193/1.4000085

CrossRef Full Text | Google Scholar

Dimitrakopoulos, E. G., and DeJong, M. J. (2012a). Overturning of Retrofitted Rocking Structures under Pulse-type Excitations. J. Eng. Mech. 138 (8), 963–972. doi:10.1061/(asce)em.1943-7889.0000410

CrossRef Full Text | Google Scholar

Dimitrakopoulos, E. G., and DeJong, M. J. (2012b). Revisiting the Rocking Block: Closed-form Solutions and Similarity Laws. Proc. R. Soc. A. 468, 2294–2318. doi:10.1098/rspa.2012.0026

CrossRef Full Text | Google Scholar

ElGawady, M. A., Ma, Q., Butterworth, J. W., and Ingham, J. (2011). Effects of Interface Material on the Performance of Free Rocking Blocks. Earthquake Engng. Struct. Dyn. 40, 375–392. doi:10.1002/eqe.1025

CrossRef Full Text | Google Scholar

Gerami, M., and Sivandi-Pour, A. (2014). Performance-based Seismic Rehabilitation of Existing Steel Eccentric Braced Buildings in Near Fault Ground Motions. Struct. Des. Tall Spec. Build. 23 (12), 881–896. doi:10.1002/tal.1088

CrossRef Full Text | Google Scholar

Hall, J. F., Heaton, T. H., Halling, M. W., and Wald, D. J. (1995). Near-source Ground Motion and its Effects on Flexible Buildings. Earthquake Spectra 11 (4), 569–605. doi:10.1193/1.1585828

CrossRef Full Text | Google Scholar

Hayden, C. P., Bray, J. D., and Abrahamson, N. A. (2014). Selection of Near-Fault Pulse Motions. J. Geotechnical Geoenvironmental Eng. ASCE 140 (7). doi:10.1061/(asce)gt.1943-5606.0001129

CrossRef Full Text | Google Scholar

Hogan, S. J. (1989). On the Dynamics of Rigid-Block Motion under Harmonic Forcing. Proc. R. Soc. Lond. A. 425, 441–476. doi:10.1098/rspa.1989.0114

CrossRef Full Text | Google Scholar

Housner, G. W. (1963). The Behavior of Inverted Pendulum Structures during Earthquakes. Bull. Seismol. Soc. Am. 53 (2), 404–417. doi:10.1785/bssa0530020403

CrossRef Full Text | Google Scholar

Ishiyama, Y. (1982). Motions of Rigid Bodies and Criteria for Overturning by Earthquake Excitations. Earthquake Engng. Struct. Dyn. 10, 635–650. doi:10.1002/eqe.4290100502

CrossRef Full Text | Google Scholar

Kalkan, E., and Kunnath, S. K. (2006). Effects of Fling Step and Forward Directivity on Seismic Response of Buildings. Earthquake Spectra 22 (2), 367–390. doi:10.1193/1.2192560

CrossRef Full Text | Google Scholar

Khaloo, A. R., Khosravi, H., and Jamnani, H. H. (2015). Nonlinear Interstory Drift Contours for Idealized Forward Directivity Pulses Using "Modified Fish-Bone" Models. Adv. Struct. Eng. 18 (5), 603–627. doi:10.1260/1369-4332.18.5.603

CrossRef Full Text | Google Scholar

Khansefid, A., Bakhshi, A., and Ansari, A. (2019). Development of Declustered Processed Earthquake Accelerogram Database for the Iranian Plateau: Including Near-Field Record Categorization. J. Seismol 23, 869–888. doi:10.1007/s10950-019-09839-w

CrossRef Full Text | Google Scholar

Khansefid, A. (2020). Pulse-like Ground Motions: Statistical Characteristics, and GMPE Development for the Iranian Plateau. Soil Dyn. Earthquake Eng. 134, 106164. doi:10.1016/j.soildyn.2020.106164

CrossRef Full Text | Google Scholar

Kohrangi, M., Vamvatsikos, D., and Bazzurro, P. (2018). Pulse-like versus Non-pulse-like Ground Motion Records: Spectral Shape Comparisons and Record Selection Strategies. Earthquake Engng Struct. Dyn. 48 (1), 46–64. doi:10.1002/eqe.3122

CrossRef Full Text | Google Scholar

Kojima, K., and Takewaki, I. (2016). Closed-form Dynamic Stability Criterion for Elastic-Plastic Structures under Near-Fault Ground Motions. Front. Built Environ. 2, 6. doi:10.3389/fbuil.2016.00006

CrossRef Full Text | Google Scholar

Kojima, K., and Takewaki, I. (2015a). Critical Earthquake Response of Elastic-Plastic Structures under Near-Fault Ground Motions (Part 1: Fling-step Input). Front. Built Environ. 1, 12. doi:10.3389/fbuil.2015.00012

CrossRef Full Text | Google Scholar

Kojima, K., and Takewaki, I. (2015b). Critical Earthquake Response of Elastic-Plastic Structures under Near-Fault Ground Motions (Part 2: Forward-Directivity Input). Front. Built Environ. 1, 13. doi:10.3389/fbuil.2015.00013

CrossRef Full Text | Google Scholar

Makris, N., and Kampas, G. (2016). Size versus Slenderness: Two Competing Parameters in the Seismic Stability of Free-Standing Rocking Columns. Bull. Seismol. Soc. Am. 106 (1), 104–122. doi:10.1785/0120150138

CrossRef Full Text | Google Scholar

Makris, N., and Black, C. J. (2004). Dimensional Analysis of Rigid-Plastic and Elastoplastic Structures under Pulse-type Excitations. J. Eng. Mech. 130, 1006–1018. doi:10.1061/(asce)0733-9399(2004)130:9(1006)

CrossRef Full Text | Google Scholar

Mavroeidis, G. P., Dong, G., and Papageorgiou, A. S. (2004). Near-fault Ground Motions, and the Response of Elastic and Inelastic single-degree-of-Freedom(SDOF) Systems. Earthquake Engng. Struct. Dyn. 33, 1023–1049. doi:10.1002/eqe.391

CrossRef Full Text | Google Scholar

Mavroeidis, G. P., and Papageorgiou, A. S. (2003). A Mathematical Representation of Near-Fault Ground Motions. Bull. Seismological Soc. America 93 (3), 1099–1131. doi:10.1785/0120020100

CrossRef Full Text | Google Scholar

Milne, J. (1885). Seismic Experiments. Trans. Seism. Soc. Jpn. 8, 1–82.

Google Scholar

Mukhopadhyay, S., and Gupta, V. K. (2013a). Directivity Pulses in Near-Fault Ground Motions-I: Identification, Extraction and Modeling. Soil Dyn. Earthquake Eng. 50, 1–15. doi:10.1016/j.soildyn.2013.02.017

CrossRef Full Text | Google Scholar

Mukhopadhyay, S., and Gupta, V. K. (2013b). Directivity Pulses in Near-Fault Ground Motions-II: Estimation of Pulse Parameters. Soil Dyn. Earthquake Eng. 50, 38–52. doi:10.1016/j.soildyn.2013.02.019

CrossRef Full Text | Google Scholar

Nabeshima, K., Taniguchi, R., Kojima, K., and Takewaki, I. (2016). Closed-form Overturning Limit of Rigid Block under Critical Near-Fault Ground Motions. Front. Built Environ. 2, 9. doi:10.3389/fbuil.2016.00009

CrossRef Full Text | Google Scholar

Pompei, A., Scalia, A., and Sumbatyan, M. A. (1998). Dynamics of Rigid Block Due to Horizontal Ground Motion. J. Eng. Mech. 124 (7), 713–717. doi:10.1061/(asce)0733-9399(1998)124:7(713)

CrossRef Full Text | Google Scholar

Prieto, F., Lourenço, P. B., and Oliveira, C. S. (2004). Impulsive Dirac-delta Forces in the Rocking Motion. Earthquake Engng. Struct. Dyn. 33, 839–857. doi:10.1002/eqe.381

CrossRef Full Text | Google Scholar

Rupakhety, R., and Sigbjörnsson, R. (2011). Can Simple Pulses Adequately Represent Near-Fault Ground Motions? J. Earthquake Eng. 15, 1260–1272. doi:10.1080/13632469.2011.565863

CrossRef Full Text | Google Scholar

Sasani, M., and Bertero, V. V. (2000). Proc. Of the Twelfth World Conf. Auckland, New Zealand: on Earthquake Eng.Importance of Severe Pulse-type Ground Motions in Performance-Based Engineering: Historical and Critical Review

Google Scholar

Shenton III, H. W., and Jones, N. P. (1991). Base Excitation of Rigid Bodies. I: Formulation. J. Eng. Mech. 117 (10), 2286–2306. doi:10.1061/(asce)0733-9399(1991)117:10(2286)

CrossRef Full Text | Google Scholar

Spanos, P. D., and Koh, A. S. (1984). Rocking of Rigid Blocks Due to Harmonic Shaking. J. Eng. Mech. 110 (11), 1627–1642. doi:10.1061/(asce)0733-9399(1984)110:11(1627)

CrossRef Full Text | Google Scholar

Takewaki, I., and Kojima, K. (2021). An Impulse and Earthquake Energy Balance Approach in Nonlinear Structural Dynamics. Boca Raton, FL: CRC Press.

Taniguchi, R., Kojima, K., and Takewaki, I. (2016). Critical Response of 2DOF Elastic-Plastic Building Structures under Double Impulse as Substitute of Near-Fault Ground Motion. Front. Built Environ. 2, 2. doi:10.3389/fbuil.2016.00002

CrossRef Full Text | Google Scholar

Taniguchi, R., Nabeshima, K., Kojima, K., and Takewaki, I. (2017). Closed-form Rocking Vibration of Rigid Block under Critical and Non-critical Double Impulse. Int. J. Earthq. Impact Eng. 2 (1), 32–45. doi:10.1504/ijeie.2017.083708

CrossRef Full Text | Google Scholar

Vafaei, D., and Eskandari, R. (2015). Seismic Response of Mega Buckling-Restrained Braces Subjected to Fling-step and Forward-Directivity Near-Fault Ground Motions. Struct. Des. Tall Spec. Build. 24, 672–686. doi:10.1002/tal.1205

CrossRef Full Text | Google Scholar

Xu, Z., Agrawal, A. K., He, W.-L., and Tan, P. (2007). Performance of Passive Energy Dissipation Systems during Near-Field Ground Motion Type Pulses. Eng. Structures 29, 224–236. doi:10.1016/j.engstruct.2006.04.020

CrossRef Full Text | Google Scholar

Yang, D., and Zhou, J. (2014). A Stochastic Model and Synthesis for Near-Fault Impulsive Ground Motions. Earthquake Engng Struct. Dyn. 44, 243–264. doi:10.1002/eqe.2468

CrossRef Full Text | Google Scholar

Yilmaz, C., Gharib, M., and Hurmuzlu, Y. (2009). Solving Frictionless Rocking Block Problem with Multiple Impacts. Proc. R. Soc. A. 465, 3323–3339. doi:10.1098/rspa.2009.0273

CrossRef Full Text | Google Scholar

Yim, C.-S., Chopra, A. K., and Penzien, J. (1980). Rocking Response of Rigid Blocks to Earthquakes. Earthquake Engng. Struct. Dyn. 8 (6), 565–587. doi:10.1002/eqe.4290080606

CrossRef Full Text | Google Scholar

Zhai, C., Chang, Z., Li, S., Chen, Z., and Xie, L. (2013). Quantitative Identification of Near-Fault Pulse-like Ground Motions Based on Energy. Bull. Seismological Soc. America 103 (5), 2591–2603. doi:10.1785/0120120320

CrossRef Full Text | Google Scholar

Zhang, J., and Makris, N. (2001). Rocking Response of Free-Standing Blocks under Cycloidal Pulses. J. Eng. Mech. 127 (5), 473–483. doi:10.1061/(asce)0733-9399(2001)127:5(473)

CrossRef Full Text | Google Scholar

Keywords: earthquake response, near-fault ground motion, forward-directivity input, triple impulse, rigid block, rocking, critical response, overturning

Citation: Homma S, Nabeshima K and Takewaki I (2021) Explicit Overturning Limit of Rigid Block Using Triple and Pseudo-Triple Impulses Under Critical Near-Fault Ground Motions. Front. Built Environ. 7:731670. doi: 10.3389/fbuil.2021.731670

Received: 27 June 2021; Accepted: 10 August 2021;
Published: 25 August 2021.

Edited by:

Ehsan Noroozinejad Farsangi, Graduate University of Advanced Technology, Iran

Reviewed by:

Abbas Sivandi-Pour, Graduate University of Advanced Technology, Iran
Ali Khansefid, K. N. Toosi University of Technology, Iran

Copyright © 2021 Homma, Nabeshima and Takewaki. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Izuru Takewaki, dGFrZXdha2lAYXJjaGkua3lvdG8tdS5hYy5qcA==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.