Skip to main content

ORIGINAL RESEARCH article

Front. Built Environ., 09 November 2021
Sec. Earthquake Engineering

A Simple Critical Response Evaluation Method for Base-Isolation Building-Connection Hybrid System Under Double Impulse as Representative of Near-Fault Ground Motion

  • Department of Architecture and Architectural Engineering, Graduate School of Engineering, Kyoto University, Kyoto, Japan

We revisit a unique building system including a base-isolation, building-connection hybrid control system. The base-isolation system withstands pulse-type earthquake ground motions effectively and the building-connection system resists long-duration earthquake ground motions efficiently. A simple smart critical response evaluation method without nonlinear time-history response analysis is proposed for this hybrid building system under near-fault ground motions. An analytical expression of the maximum elastic-plastic deformation of a damped bilinear hysteretic single-degree-of-freedom (SDOF) model under critical double impulse as a representative of pulse-type ground motions derived in our previous paper plays an important role in the development of the simple critical response evaluation method. A two-step transformation procedure into an SDOF model is proposed. The first step is the transformation of the main base-isolated building into an SDOF system and the second step is the reduction of the connecting dampers supported on a sub building to a damper with a sophisticated compensation factor on an assumed rigid wall. The evaluation of damping coefficients with the consideration of yielding of the base-isolation story is a key step in this paper. Different from the previous work, the equivalent damping coefficient is derived depending on the response range before and after yielding of the base-isolation story. This treatment enhances the accuracy of the proposed method. The accuracy and reliability of the proposed response evaluation method is demonstrated by the time-history response analysis of the multi-degree-of-freedom (MDOF) model.

1 Introduction

The resilience of building structures during and after natural disasters is a central theme of late, and many investigations have been accumulated (Bruneau and Reinhorn, 2006; Takewaki et al., 2012). Bruneau and Reinhorn (2006) introduced four factors (redundancy, robustness, rapidity, and resourcefulness) to characterize structure resilience. In earthquake-prone countries and regions, it is essential to design building structures so as to resist severe earthquakes without major damage that obstructs their continuing use (Amadio et al., 2003; Takewaki et al., 2012; Takewaki, 2013; Takewaki et al., 2013). Since intrinsically uncertain characteristics of earthquake ground motions are inevitable, the reliable prediction of forthcoming events in regards to time, space, and character is extremely difficult (Takewaki, 2013; Takewaki et al., 2011; Takewaki et al., 2012; Takewaki et al., 2013). Although newer buildings systems such as passive control systems and base-isolation systems are anticipated as effective strategies for guaranteeing the structural safety of building structures, the structural properties of constituent members and elements of such innovative systems are not certain (Ben-Haim, 2006; Fujita et al., 2021). In Japan, the consideration of the variabilities of structural properties of isolators and dampers is mandatory in the design of passively controlled buildings and base-isolated buildings because the uncertainty degree (± 20–30%) of isolators and dampers is usually larger than those of building frame members (beams and columns) (Takewaki et al., 2012; Fujita et al., 2021). The concept of robustness and redundancy plays a central role in the resilient seismic-resistant design of such building structures (Ben-Haim, 2006; Takewaki et al., 2012; Fujita et al., 2021). In the field of resilience, the theme of “Community Resilience” is receiving much interest recently because of the complexity of the concept of resilience and the difficulty in its realization in a real community consisting of many built environments (Mieler et al., 2015; Masoomi and van de Lindt, 2019; You et al., 2021).

Since the 1980s, base-isolated buildings have been developed rapidly. They are effective for pulse-type ground motions with predominant periods shorter than a few seconds or random earthquake ground motions (Hall et al., 1995; Jangid, 1995; Kelly, 1999; Naeim and Kelly, 1999; Jangid and Kelly, 2001; Hino et al., 2008; Takewaki and Fujita, 2009). However, their earthquake resilience is not completely guaranteed for long-period ground motions with a predominant period of 5–8 s (Kamae et al., 2004; Ariga et al., 2006). Although it is believed in general that passive energy dissipating systems (especially connected building systems) are effective for long-duration and long-period ground motions due to allowing sufficient time for energy dissipation (Luco and Barros, 1998; Basili and Angelis, 2007; Takewaki, 2007; Patel and Jangid, 2011; Takewaki et al., 2011; Takewaki et al., 2012; Richardson et al., 2013; Fukumoto and Takewaki, 2015; Kawai et al., 2021), they are not necessarily resilient against impulsive pulse-type ground motions.

In this paper, a new kind of hybrid or dual passive control building systems is treated in which a base-isolated building model is supported by a sub building (e.g., car parking tower) through oil dampers (Murase et al., 2013; Kasagi et al., 2016; Fukumoto and Takewaki 2017). Some researchers demonstrated that this new passive control system is effectively robust for two counterpart-type earthquake ground motions, i.e., pulse-type ground motions and long-duration, long-period ground motions (Murase et al., 2013; Kasagi et al., 2016; Fukumoto and Takewaki 2017). It has also been demonstrated through the input energy analysis in the frequency domain that, although the connecting oil dampers do not withstand pulse-type ground motions effectively due to the lack of time for energy dissipation, those dampers work smartly for long-duration, long-period ground motions (Taniguchi et al., 2016; Tamura et al., 2017). However, these analyses are limited to linear models. Although Hayashi (2018) and Hayashi et al. (2018) developed a simple method using a single-degree-of-freedom (SDOF) model for simulating the earthquake response of these buildings including the base-isolation, building connection hybrid damper system of nonlinear properties, the accuracy for near-fault ground motions is relatively low because of the non-robust evaluation of damping properties (Hayashi, 2018), or the method was applied only to the response under long-duration ground motion (Hayashi et al., 2018).

A simple and sophisticated response evaluation method using an SDOF model is proposed in this paper for the above-mentioned base-isolation, building-connection hybrid structural system under a near-fault ground motion. As in the previous method (Hayashi, 2018; Hayashi et al., 2018), a two-step transformation of the overall structural system into an SDOF model is conducted. The first step is the transformation into an SDOF system and the second step is the reduction of the connecting oil dampers supported on a sub building to upper-story concentrated oil dampers with a smart compensation factor on a rigid wall. It is shown that application of the previously derived analytical expression (Akehashi et al., 2018) of the maximum deformation to the reduced SDOF model under the critical double impulse as a representative of pulse-type near-fault ground motions enables the establishment of a simplified, but rather accurate response evaluation method.

2 Double Impulse as Representative of Near-Fault Ground Motion

2.1 Transformation of the Main Part of Near-Fault Ground Motion into Double Impulse

The concept of a pair of impulses with inverse directions to each other, called a double impulse, was introduced by Kojima and Takewaki (2015a) to represent a major part of near-fault ground motions. Since the characteristic main part of near-fault ground motions is well known to be expressed by a pulse-type one-cycle or 1.5-cycle sine wave and is influential particularly for tall and base-isolated buildings, with a rather long natural period (Sasani and Bertero, 2000; Mavroeidis and Papageorgiou, 2003; Makris and Black, 2004; Mavroeidis et al., 2004; Kalkan and Kunnath, 2006). Figure 1A shows an example of a near-fault ground motion (Rinaldi FN component during the 1994 Northridge earthquake) and its modeling into a one-cycle sine wave and a double impulse. Figure 1B illustrates the velocity and displacement properties of the one-cycle sine wave and the double impulse. The red arrows indicate the Dirac delta function δ(t). V is the given velocity (the input velocity level) and t0 is the time interval of two impulses. Consequently, the double impulse can be expressed by

u¨g(t)=Vδ(t)Vδ(tt0).(1)

FIGURE 1
www.frontiersin.org

FIGURE 1. Transformation of recorded ground motion into a one-cycle sine wave and double impulse, (A) Example of Rinaldi Station FN component during 1994 Northridge earthquake and its modeling into a one-cycle sine wave and double impulse, (B) Velocity and displacement properties of a one-cycle sine wave and double impulse.

To make the double impulse a reliable substitute for near-fault ground motions, it is important to compare it with the one-cycle sine wave. The equivalent one-cycle sine wave with the circular frequency ωp (the period is Tp) and the velocity amplitude Vp1 can be expressed by

u¨gSW(t)=0.5ωpVp1sin(ωpt)(0tTp=2t0),(2)

where Vp1=1.2222V (Kojima et al., 2018). This relation can be derived from the equivalence of the maximum Fourier amplitudes of both inputs.

3 Summary of Previous Research for Accuracy Upgrade

The base-isolation, building-connection hybrid control system as shown in Figure 2A was treated by Murase et al. (2013). They demonstrated the high performance of this control system for both impulsive and long-duration ground motions. Then, Hayashi (2018) tried to propose a simple response evaluation method using an SDOF model for impulsive and long-duration ground motions (Kojima and Takewaki, 2015b). It was shown that, while his simple model is effective for long-duration ground motions (Hayashi et al., 2018), it is not for impulsive ground motions resulting from the narrow-band estimation performance of the equivalent damping for a broader range of deformation including a plastic region. In this paper, a revised transformation method effective also for impulsive ground motions is presented. For this purpose, the method by Hayashi (2018) is explained first.

FIGURE 2
www.frontiersin.org

FIGURE 2. Base-isolation, building-connection hybrid control system, (A) Overview of system (Fukumoto and Takewaki, 2017), (B) MDOF shear building model with the base-isolation, building-connection hybrid control system, (C) Restoring force-deformation relation of the base-isolation story.

Consider a multi-degree-of-freedom (MDOF) shear building model with the base-isolation, building-connection hybrid control system as shown in Figure 2B. The common mass of the main building, the common mass of the sub building, and the mass of the base-isolation story are denoted by mU, mF and mI. The damping coefficient of each connecting viscous damper is c. The base-isolation story is assumed to consist of natural rubber isolators, steel dampers, and oil dampers. By neglecting the second-branch stiffness (stiffness of natural rubber after yielding of steel dampers), the restoring-force characteristic of the base-isolation story is assumed to obey an elastic-perfectly plastic model as shown in Figure 2C. The story shear force and the deformation of the base-isolation story are denoted by Qi and ui. The initial stiffness, the yield force, and the yield deformation of the base-isolation story are expressed by kI,fy,dyI and the damping coefficient of the base-isolation story is denoted by cI.

In the method of Hayashi (2018), an elastic-perfectly plastic SDOF model with viscous damping was introduced to make use of an effective response evaluation method by the critical (resonant) double impulse (Kojima et al., 2018).

In the single base-isolated building, the superstructure is modeled into an SDOF model and a 2DOF model is constructed with the super building model mass MU and the base-isolation story mass mI. By neglecting mI compared to MU in high-rise buildings, the base-isolated building is reduced to an SDOF model. The super building model mass MU is the sum of the super building masses. The stiffness kU and the damping coefficient cU of the super building are determined so as to attain the given fundamental natural period and the specified lowest damping ratio of the MDOF model. kI,fy,dyI, and cI defined above for the MDOF model are used in the 2DOF model (kI=fy/dyI). In the SDOF model, the mass, the initial stiffness, the yield deformation, and the damping coefficient are expressed by Me, ke, dye, cmain.

By applying the above-mentioned assumption, the mass of the SDOF model is obtained as

Me=MU.(3)

The series-spring assumption leads to the following initial stiffness of the SDOF model.

1/ke=1/kI+1/kU.(4)

As for the damping, the series-spring assumption of complex springs provides the damping coefficient cmain of the SDOF model (Hayashi et al., 2018).

1ke+iωecmain=1kI+iωecI+1kU+iωecU,(5)

where i is the imaginary unit and

ωe=ke/Me(6)

From Eqs 5, 6, cmain can be derived consequently as

cmain=(kIcU+kUcI)(kI+kU)(kIkUωe2cIcU)(cI+cU)(kI+kU)2+ωe2(cI+cU)2.(7)

The problem in deriving cmain is that the effect of yielding in the base-isolation story is not reflected. This issue will be overcome in this paper. On the other hand, the reduction of connecting dampers requires a sophisticated procedure which considers 1) the effect of rigid modeling of the sub building for simple treatment, 2) the effect of the location of connecting dampers, 3) the effect of height of the sub building. Hayashi (2018) introduced a coefficient βd on the connecting damping coefficient c for compensating these effects in the model where the connecting dampers are allocated in the top nc stories and the main building is connected to the sub building at the top nc stories (this model is called the RMDOF model). Consequently, the total connecting damping coefficient is determined by

C=βdncc.(8)

The compensation coefficient βd is determined so that the fundamental damping ratios of the MDOF model and the RMDOF model by complex eigenvalue analysis coincide. In the analysis, the equivalent stiffness kIeq of the base-isolation story is to be evaluated iteratively for the convergent deformation uI of the base-isolation story.

4 A Simple Critical Response Evaluation Method for Base-Isolated Building Under Double Impulse as Representative of Near-Fault Ground Motion

As explained before, it was shown that, while Hayashi’s simple model (Hayashi, 2018) is effective for long-duration ground motions (Hayashi et al., 2018), it is not for impulsive ground motions. In this paper, a revised transformation method effective for even impulsive ground motions is presented.

4.1 Plane Frame Model of Base-Isolated High-Rise Building

Consider a 40-story base-isolated frame building model as shown in Figure 3. Unlike the base-isolation system used in Section 3, lead rubber isolators and oil dampers are employed as the constituent members of the base-isolated story. The total depth of these isolators is 200 [mm]. This base-isolated frame building model is reduced to an SDOF mass-spring model with a dashpot via a 2DOF model as illustrated in Figure 3. For the simplicity of presentation of the proposed theory, the super building frame has a constant floor mass mU and the free wall (sub building) has a constant floor mass mF. Let mI denote the mass of the base-isolation story and let c denote the constant damping coefficient of connecting oil dampers placed at several stories. The base-isolation story is assumed to have a bilinear hysteretic restoring-force characteristic as usually used for lead rubber isolators. kI1,fy,dyI indicate the initial stiffness in the elastic range, the yield force, and the yield deformation of the base-isolation story, respectively (kI1=fy/dyI). The first and second stiffnesses of each lead rubber isolator are K1=20.01 [kN/mm] and K2=1.54 [kN/mm] and the force axis percept is Qd= 360.6 [kN]. These properties are multiplied by 36 (number of lead rubber bearings) to obtain the restoring-force characteristics of the base-isolation story.

FIGURE 3
www.frontiersin.org

FIGURE 3. Reduction of base-isolated tall plane frame into SDOF model.

Tup and Tsub denote the fundamental natural period of the main frame building for the fixed base-isolation story and the fundamental natural period of the sub building. Furthermore, cI indicates the total damping coefficient of oil dampers in the base-isolation story. The total damping coefficient cI of oil dampers is to be specified so as to attain the damping ratio 0.15 for the equivalent stiffness kIeq of the base-isolation story at the deformation uI=0.4(m). The structural damping ratio of the super-structure (instantaneous stiffness-proportional damping) is set to 0.03. The strength of concrete is Fc60 (60 [N/mm2]) and the Young’s modulus is 33,500 [N/mm2]. The floor mass is 286.2 × 103 [kg] and the mass of the base-isolation story is 858.6 × 103 [kg]. The parameters of the base-isolated building are shown in Table 1.

TABLE 1
www.frontiersin.org

TABLE 1. Member parameters of plane frame.

4.2 Transformation of Base-Isolated Building Frame into SDOF Model

Consider the transformation of a base-isolated building frame into a 2DOF model.

Since the fundamental natural period of the super building with a fixed base is 3.0 (s), the stiffness of the super building can be obtained as

kU=MU(2π/3.0)2.(9)

The initial stiffness kI1 and the second stiffness kI2 of the base-isolation story are determined by considering the hysteretic characteristics and the number of isolators and dampers (Figure 4A). Let dyI and fy denote the yield deformation and the yield force of the base-isolation story. The damping coefficient of the super building as an SDOF model is determined for the fundamental damping ratio = 0.03 by

cU=2×0.03×MUkU.(10)

FIGURE 4
www.frontiersin.org

FIGURE 4. Transformation from 2DOF model into SDOF model, (A) Force-deformation relation of a base-isolation story in 2DOF model and that in SDOF model, (B) Evaluation of deformation of the base-isolation story in SDOF model.

The damping coefficient of the base-isolation story is determined by considering the number of oil dampers.

The 2DOF model set above is modeled into an SDOF model as shown in Figure 3. Since mI is rather small compared to MU in high-rise buildings, mI is neglected here. Then

Me=MU.(11)

The stiffness and damping coefficient of the SDOF model are determined differently in the elastic region and the plastic region of the base-isolation story. While Hayashi (2018) determined the equivalent damping coefficient only for the initial stiffness, that is determined differently in the elastic region and the plastic region of the base-isolation story.

In the region before the yielding of the base-isolation story, the equivalent stiffness ke1 and the equivalent damping coefficient ce1 of the SDOF model can be evaluated by the following series modeling of complex springs.

1ke1+iωe1ce1=1kI1+iω1cI+1kU+iω1cU,(12)

where ω1 is the elastic fundamental natural circular frequency of the 2DOF model and the natural circular frequency of the SDOF model is obtained from ωe1=ke1/Me. From Eq. 12, ke1 and ce1 are obtained. While Hayashi (2018) obtained the elastic stiffness and the damping coefficient from Eqs 4, 5, the proposed method provides those in a unified manner from Eq. 12.

After the yielding of the steel damper in the base-isolation story, the stiffness reduces to kI2 (see Figure 4A) and the equivalent stiffness in the base-isolation story changes gradually in time. By responding to this phenomenon, assume the deformation uIn=1 of the base-isolation story and the equivalent stiffness of the base-isolation story by

uIn=1=dyI,(13)
kIeqn=1=fy/uIn=1,(14)

where n indicates the iteration cycle number. Based on these values, the equivalent stiffness and the equivalent damping coefficient of the SDOF model are obtained from

1keeq+iωe2ce2=1kIeq+iω2cI+1kU+iω2cU,(15)

where keeq denotes the equivalent stiffness of the SDOF model and is defined by

keeqn=1=fy/uen=1,(16)

ue is an assumed maximum deformation of the SDOF model and is determined by uI. In this computation, the lowest mode of the 2DOF model is used so that the equivalent stiffness of the base-isolation story in the 2DOF model attains kIeq and defines the deformation of the super building as ue for the deformation uI of the base-isolation story.

In Eq. 15, ω2 is the fundamental natural circular frequency of the 2DOF model with the equivalent stiffness kIeq of the base-isolation story. In addition, ωe2=keeq/Me. From Eq. 15, keeq and ce2 are consequently obtained.

Once ue and keeq are obtained for the SDOF model, the second stiffness ke2 of the SDOF model is determined by using the force-deformation relation of the 2DOF model. Since the equivalent damping coefficient of the SDOF model changes depending on the yielding of the base isolation story, ce2 for the post-yield case is used as the equivalent damping coefficient cmain in place of ce1 for the elastic case.

Through the above-mentioned procedure, the parameters Me, ke1, ke2, cmain for the SDOF model are obtained.

The maximum deformation of the bilinear hysteretic SDOF model to the critical double impulse was derived by Akehashi et al. (2018). In this paper, this closed-form expression is used. The relation between umax and uiso is derived as follows.

uiso={fmax/kI1(fmaxbI)/kI2(umax<dye)(umaxdye),(17)
fmax={ke1umaxke2umax+be(umax<dye)(umaxdye),(18)

where dye is the yield displacement of the SDOF model and bI, be are the force-axis intercept in the 2DOF model and the SDOF model. Because the SDOF model is modeled from the 2DOF model, the maximum deformation uiso of the base-isolation story has to be transformed from the maximum deformation umax of the SDOF model. In this process, the maximum story shear force of the SDOF model is defined by fmax , and the deformation uiso of the base-isolation story of the 2DOF model is computed for the same maximum story shear force fmax (see Figure 4B).

In Eq. 14, the parameters of the SDOF model are determined for an assumed equivalent stiffness kIeq of the base-isolation story of the 2DOF model. The assumed maximum deformation uI of the base-isolation story is determined iteratively so that the difference from the value uiso obtained from Eqs 17, 18 is minimized. The flowchart is shown in Figure 5.

FIGURE 5
www.frontiersin.org

FIGURE 5. Flow for determining equivalent stiffness.

4.3 Accuracy Check of Proposed Simple Response Evaluation Method Using SDOF Model Through Comparison With Time-History Response Analysis Result for Plane Frame Model

The parameters of the 2DOF model and the SDOF model obtained by the above-mentioned procedure are shown in Table 2.

TABLE 2
www.frontiersin.org

TABLE 2. Parameters of 2DOF model and SDOF model of base-isolated building.

Since the second stiffness ke2 and the damping coefficient cmain of the SDOF model after yielding in the base-isolation story are computed iteratively, those values change depending on the input velocity level of the critical double impulse. The relations of ke2, cmain with the input velocity V are shown in Figures 6A,B. It can be observed that, as V increases and the plastic deformation becomes larger, the second stiffness ke2 and the damping coefficient cmain of the SDOF model become larger and those values converge to constant values finally. It should be remarked that, while the damping coefficient cmain is constant at the lower level around 0.7 × 107 (N/(m/s)) in the previous method (Hayashi, 2018; Hayashi et al., 2018), the proposed method provides a larger value depending on the input velocity level. Figure 6C indicates the maximum deformation comparison of the base-isolation story among the proposed one, the 2DOF model under the one-cycle sine wave, the plane frame under the one-cycle sine wave, and the previous one by Hayashi (2018). It can be seen that, while the previous method (Hayashi, 2018) exhibits an inaccurate result, especially for a larger input level due to the factors explained in Figure 6B (the estimation of a smaller damping coefficient for the SDOF model), the proposed method for the SDOF model can evaluate the maximum deformation of the base-isolation story of the plane frame accurately. In addition, it is noted that the present model employs a bilinear hysteretic restoring-force model (second stiffness ratio = 0.077) for the base-isolation story while an elastic-perfectly plastic model was used in the previous method mentioned above. This may slightly affect the difference in Figure 6C.

FIGURE 6
www.frontiersin.org

FIGURE 6. Parameters of reduced SDOF model with respect to input velocity level and accuracy check of proposed simple response evaluation method for SDOF model with respect to input velocity level, (A) Second-branch stiffness of SDOF model, (B) Damping coefficient, (C) Maximum deformation comparison of base-isolation story among proposed one, 2DOF model under one-cycle sine wave, plane frame under one-cycle sine wave and previous work by Hayashi (2018).

5 A Simple Critical Response Evaluation Method for Base-Isolation Building-Connection Hybrid System Under Double Impulse as Representative of Near-Fault Ground Motion

Consider a base-isolation, building-connection hybrid control system as shown in Figure 7. By extending the previous method by Hayashi (2018) to a more sophisticated method including an accurate damping evaluation, a more reliable response evaluation method for a simple SDOF model is presented. This method enables the evaluation of the maximum deformation of the base-isolation story under a double impulse as a representative of near-fault ground motions without laborious time-history response analysis.

FIGURE 7
www.frontiersin.org

FIGURE 7. Plane frame including base-isolation building-connection hybrid control system and its reduction to SDOF model supported on rigid wall by compensated oil dampers.

5.1 Base-Isolation Building-Connection Hybrid System

The member parameters of this system are the same as those in Table 1. The connection dampers are placed on stories 4, 8, 12, 16, 18, 20, 22, 24, and 26. While Hayashi (2018) treated a shear building model, the present paper deals with a reinforced concrete plane frame. As in Section 4.1, the structural damping ratio of the super-structure (instantaneous stiffness-proportional damping) is set to 0.03. The strength of concrete is Fc60 (60 [N/mm2]) and the Young’s modulus is 33,500 [N/mm2]. In addition, the restoring-force characteristic of the base-isolation story is extended from the elastic-perfectly plastic model to the bilinear hysteretic model for the purpose of practicality.

5.2 Transformation of Base-Isolation Building-Connection Hybrid System into SDOF Model

The simplified SDOF model of a base-isolated building is further reduced to another SDOF model by connecting that SDOF model to a rigid wall as a representative of a stiff sub-building as shown in Figure 7. In this simplification procedure, a compensation coefficient βd for connecting dampers is introduced as in the method of Hayashi (2018) to include 1) the effect of rigid modeling of the sub building for simple treatment, 2) the effect of the location of connecting dampers, and 3) the effect of height of the sub building. The novel point is in the damping evaluation in the simplification of the base-isolated building explained in Section 4.

As in the method by Hayashi (2018), the RMDOF model with the compensation factor βd for the connecting dampers is considered as shown in Figure 8A to reflect the above three factors (rigid modeling of sub building, location of connecting dampers, height of the sub building).

FIGURE 8
www.frontiersin.org

FIGURE 8. Introduction of RMDOF model and procedure for evaluating compensation factor βd, (A) MDOF plane frame model with hybrid control system and RMDOF plane frame model, (B) Procedure for evaluating compensation factor βd without repetition.

The compensation factor βd for the connecting dampers is determined so that the lowest damping ratio hR (by complex eigenvalue analysis) of the RMDOF model coincides with that hM of the hybrid-controlled plane building frame.

hM=hR(19)

Equation 19 requires repetition for obtaining βd. To avoid this repetition, it is assumed that hR is linear to the total damper quantity. Let hR denote the lowest damping ratio of the RMDOF model with the original total damper quantity c×nc as shown in Figure 8B.

In the complex eigenvalue analysis, it is necessary to determine the equivalent stiffness of the base-isolation story. In this paper, the equivalent stiffness is to be obtained from the analysis for the SDOF model with βd=1 (initial assumption) where the equivalent stiffness of the base-isolation story is evaluated using the 2DOF model. While the Hayashi’s model and method (Hayashi, 2018) requires repetition for evaluating βd due to the inaccuracy in the evaluation of the equivalent damping coefficient, the proposed method highly enhances the reliability and accuracy because of the upgrade of the accuracy in the response evaluation in terms of the SDOF model.

Figure 9A shows the relation of the lowest-mode damping ratio of the MDOF model and the RMDOF model with the equivalent stiffness of the base-isolation story where five damping coefficients of the connecting damper are taken. The following observation can be drawn from Figure 9A

*The lowest-mode damping ratio of the MDOF model exhibits different values depending on the equivalent stiffness. The variation is large in the range of smaller equivalent stiffness.

*The lowest-mode damping ratio of the RMDOF model indicates larger values than the MDOF model. The difference is larger in the range of larger equivalent stiffness.

FIGURE 9
www.frontiersin.org

FIGURE 9. Evaluation of compensation factor, (A) Relation of lowest-mode damping ratio of MDOF model and RMDOF model with equivalent stiffness of base-isolation story where five damping coefficients of connecting damper, (B) Relation of equivalent stiffness of base-isolation story with input velocity level, (C) Relation of compensation factor with input velocity level.

Once such a relationship is obtained, the compensation factor βd can be determined by using the equivalent stiffness corresponding to the input velocity level. Figures 9B,C present the relation of the equivalent stiffness of the base-isolation story with the input velocity level and the relation of the compensation factor with the input velocity level. It can be observed from Figures 9B,C that, as the damping coefficients of connecting dampers becomes large, the compensation factor becomes smaller for a smaller input velocity level. It is noted that only eigenvalue analysis is needed for MDOF and RMDOF models and the response evaluation is conducted only for the SDOF model by using the closed-form expression (Akehashi et al., 2018).

Figure 10 presents the comparison of the maximum deformation of the base-isolation story by the proposed method for the SDOF model (with and without compensation factor) with that by time-history response analysis for the plane frame model for various connection damper levels. It can be observed that, while the proposed method exhibits a good correspondence to the plane frame model regardless of the compensation factor βd in the small connection damper level, the compensation factor plays an important role in the larger connection damper level. It seems that this results from the fact that, when the connection damper level becomes larger, the effect of the assumption of the rigid sub-building arises strongly. This fact can also be understood from the result that, as the connection damper level becomes larger, βd becomes smaller (see Figure 9C). This indicates the compensation effect clearly.

FIGURE 10
www.frontiersin.org

FIGURE 10. Comparison of the maximum deformation of the base-isolation story by the proposed method for the SDOF model (with and without compensation factor) with that by time-history response analysis for the plane frame model for various connection damper levels, (A) c = 0.208 [kN/(mm/s)], (B) c = 0.417 [kN/(mm/s)], (C) c = 0.521 [kN/(mm/s)], (D) c = 0.625 [kN/(mm/s)], (E) c = 0.833 [kN/(mm/s)].

6 Verification of Validity of Using One-Cycle Sine Wave Equivalent to Critical Double Impulse for SDOF Model as Critical Input for Plane Frame

In Sections 4, 5, the one-cycle sine wave equivalent to the critical double impulse for the SDOF model was used to verify the accuracy of the proposed response evaluation method. However, it seems necessary to check whether the critical double impulse for the SDOF model is actually critical for the plane frame model. In this section, this issue will be investigated.

6.1 Critical Double Impulse for Plane Frame

Recently, many multi-purpose structural analysis programs have become available. In this paper, a general-purpose structural analysis program SNAP (Kozo System Co., 2019) is used. In this program, the accuracy check is included automatically, e.g., the processing of unbalanced forces.

The double impulse u¨gDI is to be expressed by a one-point (triangular) acceleration input u¨gTRI where the triangle area of this input corresponds to the impulse VI,N. This relation can be expressed by

AI,NΔt=VI,N.(20)

Kawai et al. (2021) showed that Δt=0.01 [s] is sufficient enough for the accuracy. It is noted that, since a rapid change of input acceleration may cause the occurrence of unbalanced forces, division of the time interval into short durations is often conducted.

Furthermore, Akehashi and Takewaki (2019) made clear that the critical input of the second impulse can be characterized by the criterion that the critical timing of the second impulse is the time of the zero story shear in the first story. For the present model, the following relation holds.

F=FM1+FF1=0,(21)

where FM1 is the base-isolation story shear force of the main building and FF1 is the first-story shear force of the sub building. Figure 11 shows the critical timing of the second impulse and the simulation of the impulses by one-point (triangular) acceleration inputs.

FIGURE 11
www.frontiersin.org

FIGURE 11. Critical timing of second impulse and simulation of impulse by one-point acceleration input.

6.2 Comparison of the Response of Plane Frame to One-Cycle Sine Wave With Response to Impulsive Acceleration Input

Figure 12A shows the comparison of the maximum deformation of the base-isolation story of the plane frame model under the one-cycle sine wave equivalent to the critical double impulse for the SDOF model with that under the critical double impulse (a pair of inverse-direction one-point acceleration inputs) for the plane frame model. The connection damping coefficient is given by c = 0.417 [kN/(mm/s)] and the time increment is specified as Δt=0.01 [s]. It can be observed that both plots correspond well. This indicates the reliability and accuracy of using the SDOF model and the critical double impulse for the SDOF model.

FIGURE 12
www.frontiersin.org

FIGURE 12. Comparison of deformation of the base-isolation story of plane frame model under one-cycle sine wave equivalent to the critical double impulse for SDOF model with that under critical double impulse (a pair of inverse-direction one-point acceleration inputs) for plane frame model, (A) Maximum deformation of base-isolation story, (B) Time history of deformation of the base-isolation story (input velocity level V = 1.0 [m/s]).

Figure 12B illustrates the time histories of the above-mentioned two inputs. The phases of the time histories are adjusted. It can be understood that the responses under both inputs coincide well and the simple response evaluation method proposed in the paper is reliable.

7 Investigation on Performance of Proposed Evaluation Method Through Time-History Response Analysis for Recorded Ground Motion

It seems important to demonstrate the applicability of the proposed simple response evaluation method to actual recorded ground motions.

Figure 13 shows the applicability of the proposed evaluation method to a recorded ground motion through the time-history response analysis (Rinaldi Station FN comp. during1994 Northridge). Figure 13A presents the recorded ground motion acceleration and the corresponding one-cycle sinusoidal wave where the acceleration amplitude and the period of the extracted one-cycle sinusoidal wave are 7.85 [m/s] and 0.8 [s]. As for connecting dampers, the damping coefficient c = 2.5 × 106 [Ns/m] is used. The other parameters are the same as in the previous section. Figure 13B illustrates the maximum displacement of the SDOF model by the proposed method under the critical double impulse and by the time-history response analysis under the recorded ground motion. Since the amplitude and the period of the one-cycle sinusoidal wave are fixed, the model parameters Vye,ω1,dye are changed with the relation Vye=ω1dye where ω1 is the elastic natural circular frequency and dye is the elastic limit deformation of the SDOF model. This procedure may be similar to the work by Veletsos et al. (1965) for the inelastic response spectrum where the yield displacement or yield strength is changed to attain the specified ductility factor. It can be observed that the proposed method exhibits an accurate result. Figure 13C shows the time-history response for the SDOF model under the critical double impulse corresponding to the one-cycle sinusoidal wave and the recorded ground motion for the input velocity level V/Vye=4.19. It can be seen that the maximum displacement exhibits a good correspondence.

FIGURE 13
www.frontiersin.org

FIGURE 13. Comparison of proposed evaluation method with time-history response analysis for recorded ground motion, (A) Rinaldi Station FN comp. (1994 Northridge) and one-cycle sinusoidal wave, (B) Maximum displacement of SDOF model by the proposed method under critical double impulse and by time-history response analysis under recorded ground motion, (C) Time-history response analysis for SDOF model under critical double impulse and recorded ground motion for V/Vye=4.19.

8 Conclusion

An innovative building system including a base-isolation, building-connection hybrid control system proposed in the previous study has been revisited. This building system has advantageous features to withstand pulse-type earthquake ground motions through the base-isolation system and to resist long-duration earthquake ground motions through the building-connection system.

1) While the previous response evaluation method with an amplitude-narrow-band damping evaluation procedure exhibited a rather inaccurate result, a simple and more accurate response evaluation method without nonlinear time-history response analysis has been proposed for this hybrid-controlled building system under near-fault ground motions.

2) It was demonstrated that the previously derived analytical expression of the elastic-plastic deformation of a bilinear hysteretic single-degree-of-freedom (SDOF) model with lead rubber isolators and oil dampers can be used effectively under the double impulse as a representative of pulse-type ground motions.

3) A two-step transformation from this innovative hybrid structural system into an SDOF model via a 2DOF model has been proposed. The first step is the transformation of the main base-isolated building into an SDOF system via a 2DOF model and the second step is the reduction of the connecting oil dampers supported on a sub building to upper-story concentrated oil dampers with a revised compensation factor on a rigid wall. The compensation factor reflects 1) the effect of rigid modeling of the sub building, 2) the effect of the location of connecting dampers, 3) the effect of height of the sub building.

4) It was made clear that the evaluation of damping coefficients is a key step in the upgrade of accuracy. Different from the previous work (Hayashi, 2018; Hayashi et al., 2018), the equivalent damping coefficient was obtained depending on the response range before and after yielding. Effective use of the closed-form expression of the elastic-plastic deformation to the reduced SDOF model enables the development of a simple and sophisticated response evaluation method.

5) The time-history response analysis of the plane frame model under the critical double impulse and a one-cycle sine wave equivalent to the critical double impulse for an SDOF model demonstrated the accuracy of the proposed response evaluation method.

6) It was demonstrated that the proposed method exhibits a reliable performance for recorded ground motions.

The behavior of the sub-building seems to influence the total response of this base-isolation, building-connection hybrid control system because the sub-building is usually made of reinforced-concrete wall structures whose structural properties have never been made clear. The reflection of this property on the total response evaluation system will be desired for a more accurate and reliable design.

Data Availability Statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Author Contributions

TN formulated the problem, conducted the computation, and wrote the paper. KF 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 JSPS KAKENHI (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

Akehashi, H., Kojima, K., and Takewaki, I. (2018). Critical Response of Single-Degree-Of-freedom Damped Bilinear Hysteretic System under Double Impulse as Substitute for Near-Fault Ground Motion. Front. Built Environ. 4, 5. doi:10.3389/fbuil.2018.00005

CrossRef Full Text | Google Scholar

Akehashi, H., and Takewaki, I. (2019). Optimal Viscous Damper Placement for Elastic-Plastic MDOF Structures under Critical Double Impulse. Front. Built Environ. 5, 20. doi:10.3389/fbuil.2019.00020

CrossRef Full Text | Google Scholar

Amadio, C., Fragiacomo, M., and Rajgelj, S. (2003). The Effects of Repeated Earthquake Ground Motions on the Non-linear Response of SDOF Systems. Earthquake Engng. Struct. Dyn. 32, 291–308. doi:10.1002/eqe.225

CrossRef Full Text | Google Scholar

Ariga, T., Kanno, Y., and Takewaki, I. (2006). Resonant Behaviour of Base-Isolated High-Rise Buildings under Long-Period Ground Motions. Struct. Des. Tall Spec. Build. 15 (3), 325–338. doi:10.1002/tal.298

CrossRef Full Text | Google Scholar

Basili, M., and Angelis, M. D. (2007). Optimal Passive Control of Adjacent Structures Interconnected with Nonlinear Hysteretic Devices. J. Sound Vib. 301 (1-2), 106–125. doi:10.1016/j.jsv.2006.09.027

CrossRef Full Text | Google Scholar

Ben-Haim, Y. (2006). Information-gap Decision Theory: Decisions under Severe Uncertainty. London: Academic Press.

Google Scholar

Bruneau, M., and Reinhorn, A. (2006). “Overview of the Resilience Concept,” in Proc. of the 8th US Nat. Conf. Earthq. Eng.

Google Scholar

Fujita, K., Wataya, R., and Takewaki, I. (2021). Robust Optimal Damper Placement of Nonlinear Oil Dampers with Uncertainty Using Critical Double Impulse. Front. Built Environ. 7, 744973. doi:10.3389/fbuil.2021.744973

CrossRef Full Text | Google Scholar

Fukumoto, Y., and Takewaki, I. (2015). Critical Demand of Earthquake Input Energy to Connected Building Structures. Earthquakes and Structures 9 (6), 1133–1152. doi:10.12989/eas.2015.9.6.1133

CrossRef Full Text | Google Scholar

Fukumoto, Y., and Takewaki, I. (2017). Dual Control High-Rise Building for Robuster Earthquake Performance, Front. Built Environ., 3: 12. doi:10.3389/fbuil.2017.00012

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, 569–605. doi:10.1193/1.1585828

CrossRef Full Text | Google Scholar

Hayashi, K., Fujita, K., Tsuji, M., and Takewaki, I. (2018). A Simple Response Evaluation Method for Base-Isolation Building-Connection Hybrid Structural System under Long-Period and Long-Duration Ground Motion. Front. Built Environ. 4, 2. doi:10.3389/fbuil.2018.00002

CrossRef Full Text | Google Scholar

Hayashi, K. (2018). Simple Response Evaluation Methods for Base-Isolation Building-Connection Hybrid Structural System under Long-Period, Long-Duration Ground Motions and Near-Fault Ground Motions. Thesis of Master of Engineering (in Japanese: Kyoto University).

Hino, J., Yoshitomi, S., Tsuji, M., and Takewaki, I. (2008). Bound of Aspect Ratio of Base-Isolated Buildings Considering Nonlinear Tensile Behavior of Rubber Bearing. Struct. Eng. Mech. 30 (3), 351–368. doi:10.12989/sem.2008.30.3.351

CrossRef Full Text | Google Scholar

Jangid, R. S., and Kelly, J. M. (2001). Base Isolation for Near-Fault Motions. Earthquake Engng. Struct. Dyn. 30, 691–707. doi:10.1002/eqe.31

CrossRef Full Text | Google Scholar

Jangid, R. S. (1995). Optimum Isolator Damping for Minimum Acceleration Response of Base-Isolated Structures. Aust. Civil Eng. Trans. 37, 325–331.

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

Kamae, K., Kawabe, H., and Irikura, K. (2004). “Strong Ground Motion Prediction for Huge Subduction Earthquakes Using a Characterized Source Model and Several Simulation Techniques,” in Proc. of the13th WCEE (Vancouver.

Google Scholar

Kasagi, M., Fujita, K., Tsuji, M., and Takewaki, I. (2016). Automatic Generation of Smart Earthquake-Resistant Building System: Hybrid System of Base-Isolation and Building-Connection. Heliyon 2, e00069. doi:10.1016/j.heliyon.2016.e00069

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawai, A., Maeda, T., and Takewaki, I. (2021). Critical Response of High-Rise Buildings with Deformation-Concentration Seismic Control System under Double and Multi Impulses Representing Pulse-type and Long-Duration Ground Motions. Front. Built Environ. 7, 649224. doi:10.3389/fbuil.2021.649224

CrossRef Full Text | Google Scholar

Kelly, J. M. (1999). The Role of Damping in Seismic Isolation. Earthquake Engng. Struct. Dyn. 28, 3–20. doi:10.1002/(sici)1096-9845(199901)28:1<3:aid-eqe801>3.0.co;2-d

CrossRef Full Text | Google Scholar

Kojima, K., Saotome, Y., and Takewaki, I. (2018). Critical Earthquake Response of a SDOF Elastic-Perfectly Plastic Model with Viscous Damping under Double Impulse as a Substitute for Near-Fault Ground Motion. Jpn. Archit Rev. 1 (2), 207–220. doi:10.1002/2475-8876.10019

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 Input and Response of Elastic-Plastic Structures under Long-Duration Earthquake Ground Motions. Front. Built Environ. 1, 15. doi:10.3389/fbuil.2015.00015

CrossRef Full Text | Google Scholar

Kozo System Co (2019). SNAP: Software for Elastic-Plastic Analysis of Three-Dimensional Frames of Arbitrary Shapes. Version 7.

Google Scholar

Luco, J. E., and De Barros, F. C. P. (1998). Optimal Damping between Two Adjacent Elastic Structures. Earthquake Engng. Struct. Dyn. 27 (7), 649–659. doi:10.1002/(sici)1096-9845(199807)27:7<649:aid-eqe748>3.0.co;2-5

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

Masoomi, H., and van de Lindt, J. W. (2019). Community-resilience-based Design of the Built Environment. Asce-asme J. Risk Uncertainty Eng. Syst. Part. A: Civ. Eng. 5 (1), 04018044. doi:10.1061/ajrua6.0000998

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

Mieler, M., Stojadinovic, B., Budnitz, R., Comerio, M., and Mahin, S. (2015). A Framework for Linking Community-Resilience Goals to Specific Performance Targets for the Built Environment. Earthquake Spectra 31 (3), 1267–1283. doi:10.1193/082213eqs237m

CrossRef Full Text | Google Scholar

Murase, M., Tsuji, M., and Takewaki, I. (2013). Smart Passive Control of Buildings with Higher Redundancy and Robustness Using Base-Isolation and Inter-connection. Earthquakes and Structures 4 (6), 649–670. doi:10.12989/eas.2013.4.6.649

CrossRef Full Text | Google Scholar

Naeim, F., and Kelly, J. M. (1999). Design of Seismic Isolated Structures. New York: Wiley.

Google Scholar

Patel, C. C., and Jangid, R. S. (2011). Dynamic Response of Adjacent Structures Connected by Friction Damper. Earthquakes and Structures 2 (2), 149–169. doi:10.12989/eas.2011.2.2.149

CrossRef Full Text | Google Scholar

Richardson, A., Walsh, K. K., and Abdullah, M. M. (2013). Closed-form Design Equations for Controlling Vibrations in Connected Structures. J. Earthquake Eng. 17 (5), 699–719. doi:10.1080/13632469.2013.771590

CrossRef Full Text | Google Scholar

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

Google Scholar

Takewaki, I., Moustafa, A., and Fujita, K. (2012). Improving the Earthquake Resilience of Buildings: The Worst Case Approach. London: Springer.

Google Scholar

Takewaki, I. (2013). Critical Excitation Methods in Earthquake Engineering. 2nd edition. Amsterdam: Elsevier Science.

Google Scholar

Takewaki, I. (2007). Earthquake Input Energy to Two Buildings Connected by Viscous Dampers. J. Struct. Eng. 133 (5), 620–628. doi:10.1061/(asce)0733-9445(2007)133:5(620)

CrossRef Full Text | Google Scholar

Takewaki, I., and Fujita, K. (2009). Earthquake Input Energy to Tall and Base-Isolated Buildings in Time and Frequency Dual Domains. Struct. Des. Tall Spec. Build. 18 (6), 589–606. doi:10.1002/tal.497

CrossRef Full Text | Google Scholar

Takewaki, I., Fujita, K., and Yoshitomi, S. (2013). Uncertainties in Long-Period Ground Motion and its Impact on Building Structural Design: Case Study of the 2011 Tohoku (Japan) Earthquake. Eng. Structures 49, 119–134. doi:10.1016/j.engstruct.2012.10.038

CrossRef Full Text | Google Scholar

Takewaki, I., Murakami, S., Fujita, K., Yoshitomi, S., and Tsuji, M. (2011). The 2011 off the Pacific Coast of Tohoku Earthquake and Response of High-Rise Buildings under Long-Period Ground Motions. Soil Dyn. Earthquake Eng. 31 (11), 1511–1528. doi:10.1016/j.soildyn.2011.06.001

CrossRef Full Text | Google Scholar

Tamura, G., Taniguchi, M., Fujita, K., Tsuji, M., and Takewaki, I. (2017). Optimal Damper Placement in Hybrid Control System of Multiple Isolation and Building Connection. Int. J. Earthq. Impact Eng. 2 (1), 67–87. doi:10.1504/ijeie.2017.083718

CrossRef Full Text | Google Scholar

Taniguchi, M., Fujita, K., Tsuji, M., and Takewaki, I. (2016). Hybrid Control System for Greater Resilience Using Multiple Isolation and Building Connection. Front. Built Environ. 2, 26. doi:10.3389/fbuil.2016.00026

CrossRef Full Text | Google Scholar

Veletsos, A. S., Newmark, N. M., and Chelapati, C. V. (1965). “Deformation Spectra for Elastic and Elasto-Plastic Systems Subjected to Ground Shock and Earthquake Motions,” in Proc. of the Third World Conference on Earthquake Engineering (New Zealand.

Google Scholar

You, T., Wang, W., and Chen, Y. (2021). A Framework to Link Community Long-Term Resilience Goals to Seismic Performance of Individual Buildings Using Network-Based Recovery Modeling Method. Soil Dyn. Earthquake Eng. 147, 106788. doi:10.1016/j.soildyn.2021.106788

CrossRef Full Text | Google Scholar

Keywords: base-isolation, building-connection, hybrid control, passive control, near-fault ground motion, double impulse, critical response

Citation: Nakamura T, Fujita K and Takewaki I (2021) A Simple Critical Response Evaluation Method for Base-Isolation Building-Connection Hybrid System Under Double Impulse as Representative of Near-Fault Ground Motion. Front. Built Environ. 7:790584. doi: 10.3389/fbuil.2021.790584

Received: 07 October 2021; Accepted: 22 October 2021;
Published: 09 November 2021.

Edited by:

Ehsan Noroozinejad Farsangi, Graduate University of Advanced Technology, Iran

Reviewed by:

Alireza Tabrizikahou, Poznań University of Technology, Poland
Vahidreza Gharehbaghi, University of Southern Queensland, Australia

Copyright © 2021 Nakamura, Fujita 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, takewaki@archi.kyoto-u.ac.jp

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.