Next Article in Journal
Wind Speed Measurement via Visual Recognition of Wind-Induced Waving Light Stick Target
Next Article in Special Issue
A Fracture Model for Dynamic Sliding Safety Evaluation of a Concrete Dam Subjected to Seismic Excitation
Previous Article in Journal
Pareto-Optimized AVQI Assessment of Dysphonia: A Clinical Trial Using Various Smartphones
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Seismic Performance of Frame Structure with Hysteretic Intermediate Discontinuity

by
Angelo Di Egidio
1,
Stefano Pagliaro
1 and
Alessandro Contento
2,*
1
Department of Civil, Construction-Architectural and Environmental Engineering, University of L’Aquila, 67100 L’Aquila, Italy
2
College of Engineering, Fuzhou University, Fuzhou 350108, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(9), 5373; https://doi.org/10.3390/app13095373
Submission received: 20 March 2023 / Revised: 18 April 2023 / Accepted: 24 April 2023 / Published: 25 April 2023
(This article belongs to the Special Issue Structural Seismic Design and Evaluation)

Abstract

:
The introduction of an intermediate discontinuity in frame structures is commonly named inter-storey isolation. Inter-storey isolation is an effective technique for the seismic protection of new or existing frame structures. The devices that are used to perform the discontinuity, mainly of the structural stiffness, are placed at a higher storey level of the structure and not at the base level. In the latter years, this technique has gained increasing interest because, especially for existing buildings, it is cheaper and technically easier to implement than base isolation. In this paper, the attention is focused on the effects on a frame structure of an intermediate elasto-plastic discontinuity that can be described by a Bouc-Wen model. The frame structure with the intermediate discontinuity is modelled with a 3-DOF reduced model. Its dynamical behaviour is investigated by considering both harmonic and seismic external excitation. The results are summarized in gain maps aimed at finding the parameters that optimize the seismic behaviour of the structure.

1. Introduction

The seismic protection of structures is one of the biggest challenges in civil engineering. In recent years, the building technical codes increased the requirements in terms of structural capacity that a structure needs to fulfil, both for new buildings and the seismic retrofit of existing buildings. New buildings can meet the needed requirements by means of a proper design. Instead, improving the seismic response of existing buildings may require invasive interventions that may modify the serviceability of such buildings. There is a wide literature available about the seismic retrofit of structures. Base isolation (BI) and Tuned Mass Damper (TMD) are two widely known techniques. The first studies on base isolation and tuned mass damper systems were presented in [1,2], respectively. Tsai, in [3], and Taniguchi et al., in [4], showed that a TMD placed on the floor immediately above the isolation layer in a base isolated structure can reduce the base displacement. A more complex version of TMD, named Tuned Mass Damper Inerter (TMDI), was presented in [5] and was obtained by coupling a TMD with an inerter. In [5], such a device was intended to reduce the displacement demand of base-isolated structures, but in [6], instead, the same device was used to reduce the vibrations in tall buildings.
Protection systems based on the TMD have been used with limited success in the seismic protection of structures because of the difficulty in tuning the TMD to multiple frequencies, as it is needed for seismic excitations. De Angelis et al. [7] showed that it is possible to decrease a TMD tendency to detuning by designing a TMD with a high ratio between the mass of the TMD and the mass of the structure to be protected. To achieve such a high mass ratio, some authors purposed the combined use of TMD and BI. This combination is commonly named inter-storey isolation: in inter-storey isolation, base-isolated and tuned-mass-damper behaviours coexist.
The inter-storey isolation has become widely used for the seismic protection of structures also because it is usually characterized by a low impact on the aesthetic of the building. Even though there are many works about the use of inter-storey isolation, most of them are related to case studies, and only a few of them are devoted to clarify the general aspects of its mechanical behaviour. For example, Chey et al. added new storeys on the rooftop of an existing 12-storey building [8] and connected the two parts of the structure with an isolation layer, thus making the upper part of the structure work as a TMD to reduce the vibration of the lower part of the structure. Wang et al. [9,10,11] presented analytical and experimental analyses of the response of building structures with inter-storey isolation. In [12,13,14], the roof of the building was used as a vibrating damper to protect the structure, whereas in [15,16], entire sections of the studied structures were converted into vibrating masses by means of a uniform disconnection along the height of the building. An analogous technique for the reduction of the seismic displacements is the partial mass isolation (PMI) technique presented in [17], where the performance of the PMI is compared to those of conventional TMD and BI.
In order to highlight the general mechanisms behind inter-storey isolation, a suitable option is to analyse its effects on an archetype model of a structure. In the literature, archetype models are often used to highlight the main characteristics of multi-Degree-Of-Freedom (DOF) models. For instance, in [1,2,18,19], the studied systems are represented by archetype 2-DOF models. A similar approach, aimed at studying the general behaviour of the inter-storey isolation with a 2-DOF archetype model of the structure, has been presented in [20,21,22]. Specifically, in [20,21] and in this paper, the inter-storey isolation is named intermediate discontinuity, since it represents a sudden change in the stiffness of the structure. A 2-DOF model is suitable to describe systems where the isolation is placed at higher or lower levels, but cannot properly describe systems where the isolation is placed at middle-floors. In [23], a third degree of freedom has been added to the model to improve the describing capacity of the inter-storey isolation. In that paper, the disconnection devices are modelled by linear links with a non-classical Rayleigh damping model. However, the choice of a linear model to describe the behaviour of the disconnection could prevent the correct evaluation of the effectiveness of the introduction of the intermediate discontinuity. Hence, a non-linear modelling of the discontinuity, with constitutive laws closer to the real behaviour of the disconnection devices, could improve the ability of the model to capture the behaviour of real structures.
In this paper, a system with intermediate discontinuity is studied by means of a 3-DOF model, which represents a medium-rise frame structure (i.e., an archetype model), and a Bouc–Wen hysteresis model representing the disconnection devices, which are hypothetically assumed to be lead-rubber-bearing isolators [24,25]. The main purpose of the paper is to define the influence of the hysteretic characteristics of the discontinuity devices on the seismic response of the structure. Since in the paper the numerical simulations refer only to a six-storey frame structure, the results cannot be generalised, but are intended to provide preliminary information about the hysteretic discontinuity. Although the results are not generalizable, the method used to investigate the behaviour of the structure with intermediate discontinuity provides information for most of the medium-rise frame structures with hysteretic inter-storey isolation. Additionally, the paper does not provide technical specifications about the practical realization of an intermediate discontinuity, but it has to be framed as an analytical/numerical study, which precedes any practical implementation.
The first part of the paper is focused on the description of the mechanical behaviour of the system. In particular, Section 2 presents the mechanical model and the equations of motion. The second part is focused on the dynamical response. Specifically, Section 3 studies the response to harmonic excitation, whereas Section 4 presents the response to seismic excitation. The results are summarized in gain maps used to find the optimal parameters for the inter-storey isolation of the archetype model.

2. The Mechanical Model

2.1. Archetype Model

In the following, two archetype systems are used to compare results with and without inter-storey isolation. The Main System archetype (MS) represents the structure with the introduction of an intermediate discontinuity (Figure 1a), whereas the Auxiliary System archetype (AS) represents the same structure without the intermediate discontinuity (Figure 1b). In the case of the MS, the part of the structure below the discontinuity is called substructure, whereas the part above the discontinuity is called superstructure.
The MS and AS are 3-DOF and 2-DOF equivalent models, respectively (Figure 2). In such models, m 1 and k 1 represent the mass and the stiffness of the substructure, whereas k 2 and m 2 refer to the superstructure. In the MS, the constitutive behaviour of the disconnection devices is modelled by means of the Bouc–Wen model. There are several variations of such a model, for instance, Song and Der Kiureghian [26] proposed a generalized Bouc–Wen model to describe asymmetric hysteresis, whereas Aloisio [27] extended the six-parameters formulation to an eight-parameter formulation to include the cyclical degradation of the material properties. In this paper, the classical Bouc–Wen formulation [28] is deemed fit to accurately describe the constitutive behaviour of the disconnection devices.
It is worth observing that no specific reference is made to any new or existing structures. Moreover, it is assumed that the intermediate discontinuity keeps the structure to be protected in the elastic field. To evaluate the effectiveness of the intermediate discontinuity, the (linear elastic) displacements provided by the MS are compared with those of the AS.
For the MS, the equations of motion, obtained by a direct approach, are written as
c 11 u ˙ 1 ( t ) + c 12 u d ˙ F h ( t ) + k 1 u 1 ( t ) + m 1 u ¨ 1 ( t ) = m 1 a g ( t ) ,
c 21 u 1 ˙ + c 22 u ˙ d ( t ) + c 23 u ˙ 2 ( t ) + F h ( t ) + k 2 u d ( t ) u 2 ( t ) + m d u ¨ d ( t ) = m d a g ( t ) ,
c 32 u d ˙ ( t ) + c 33 u 2 ˙ ( t ) + k 2 u 2 ( t ) k 2 u d ( t ) + m 2 u ¨ 2 ( t ) = m 2 a g ( t ) ,
z ˙ ( t ) = u ˙ d ( t ) u ˙ 1 ( t ) { A ˜ z ( t ) n { β 1 + γ 1 sgn { z ( t ) u ˙ d ( t ) u 1 ˙ ( t ) } } } u y 1 ,
where the restoring force F h ( t ) is defined as
F h ( t ) = Ψ k d ( u d ( t ) u 1 ( t ) ) + k d u y ( 1 Ψ ) z ( t ) ,
z ( t ) is a hysteretic variable that describes the post-yielding behaviour, whereas A ˜ , β 1 , n, and γ 1 are the Bouc–Wen parameters defined in [28]. The Bouc–Wen parameters are redundant; therefore, both Charalampakis [29] and Ma et al. [30] suggested assigning specific values to some of them. In details, they suggested using A ˜ = 1 , so that k d = F y / u y acquires the meaning of pre-yielding (elastic) stiffness. Additionally, some researchers [29,31] suggested using β 1 + γ 1 = 1 , where β 1 and γ 1 are parameters that control the shape and dimension of the hysteresis cycle. The choice to adopt these two constraints, A ˜ = 1 and β 1 + γ 1 = 1 , provides a physical meaning to F y , u y , and Ψ . Specifically, F y is the yielding force, u y is the yielding displacement, and Ψ is the ratio between post-yielding and pre-yielding (elastic) stiffness. Moreover, the values of the hysteretic dimensionless parameter z ( t ) remain inside the range [ 1 , 1 ] . Therefore, in this study, the values adopted for the parameters are A ˜ = 1 , β 1 = γ 1 = 1 / 2 , and n = 2 .
The damping coefficients c i j in Equations (1a)–(1c) can be grouped in matrix form as
C = c 11 c 12 0 c 21 c 22 c 23 0 c 32 c 33 ,
and Matrix C is derived by means of the classical Rayleigh formulation where C = α R M + β R K . In this equation, α R and β R are the Rayleigh coefficients, and it is assumed that the modal damping of the first two modes is ξ = 0.05 . The elastic stiffness matrix K and the mass matrix M are obtained for Ψ = 1 .
The equations of motion of the AS are described those of a shear-type 2-DOF model. The mechanical behaviour of both MS and AS are considered elastic and linear because the aim of the proposed protection strategy is to maintain the linear elastic behaviour of the structure to be protected. The linear elastic behaviour of the AS is only considered as a reference for the comparison because the AS would conceivably undergo plastic deformations under severe earthquakes. The proposed analytical model does not account for the P-Delta effects that may arise in the columns of the substructure due to the lateral displacement of the superstructure because such displacement is considered negligible, as in [9,10,11].

2.2. Dynamic Equivalence between Multiple-DOF Systems and Archetype Models

This subsection presents the equivalence procedure used to derive the mechanical characteristics of the AS. Such a procedure was developed mainly for low- and medium-rise frame structures and aims at obtaining the masses m 1 and m 2 and stiffnesses k 1 and k 2 of the AS starting from the geometrical and mechanical characteristics of the multiple-DOF (M-DOF) frame structure without intermediate discontinuity. If the response of the M-DOF frame has a predominant vibration mode, then the AS can correctly represent the behaviour of such a M-DOF system. The equivalence procedure is described in the details in [21]. It is assumed that m 1 , m 2 , k 1 , and k 2 remain the same for the MS. The procedure is based on the following two assumptions:
  • The frequency ω m of the main mode of the M-DOF system without discontinuity and of the first mode of the AS (Figure 2b) are the same;
  • The modal displacements of the storey of the discontinuity and top storey (i.e., L-th and N-th storeys in Figure 1a of the M-DOF system without discontinuity are equal to the modal displacements u 1 and u 2 of the AS.
The masses m 1 and m 2 of the AS correspond to the total mass of the substructure and superstructure, respectively. Assuming storeys of equal mass m entails
m 1 = m L ; m 2 = m ( N L ) ,
where m is the mass of a storey, L is the level where the intermediate discontinuity is introduced, and N is the number of storeys of the frame. The masses of the MS are m ˜ 1 , m d , and m 2 . In particular, the mass of the substructure reads
m ˜ 1 = m ( L 1 ) + m 0 ; m d = m ,
where m 0 = γ m ( γ = 0.5 ) is the mass of additional stiffening elements. Such stiffening elements are those that should be added to the level just below the discontinuity to maintain the same lateral stiffness available before the introduction of the discontinuity (Figure 1a). The value γ = 0.5 is chosen because only horizontal structural elements are needed to restore the stiffness of the disconnected storey, which, consequently, has a mass smaller than that of an entire storey.
The equivalence procedure requires the identification of the main mode of the M-DOF frame structure without discontinuity. It is assumed that the main mode of the structure has the highest modal mass participation factor L . Let Φ be the vector that collects the modal displacements along a reference vertical line and along the direction to which corresponds the highest factor L . It reads
Φ = ϕ 1 , ϕ 2 , , ϕ L , , ϕ N 1 , ϕ N T
The following relationships for equivalent stiffnesses k 1 and k 2 of the AS can be obtained as in [21]:
k 1 = ω m 2 ( m 1 ϕ L + m 2 ϕ N ) ϕ L ; k 2 = ω m 2 m 2 ϕ N ϕ N ϕ L ,
where ϕ L and ϕ N are the components of the vector Φ that refer to the L and N storeys of the M-DOF system.
It is worth remarking that the stiffness values of the superstructure and substructure are assumed to remain unchanged after the discontinuity is introduced (i.e., k 1 and k 2 obtained by Equation (7) are the same for both AS and MS). The ability of both MS and AS to describe a generic 3D frame structure, without and with intermediate discontinuity, and the limitations of such archetype models were discussed in [23].

2.3. Characteristics of the Frame

The reference building is a regular six-storeys frame structure. The storey mass is m p = 301.5 kg × 10 3 , that correspond to a storey area of 300 m 2 . The viscous damping ratio of the frame structure is ξ = 5 % and the overall storey stiffness is k p = 476 , 779 KN / m . The fundamental period of such a structure is T = 0.65 s. For a six-storey building, the intermediate discontinuity can be placed at five different levels. The characteristics of the archetype models depend on the discontinuity level L and are reported in Table 1.

2.4. Gain Coefficients

The results of the MS are compared with those of the AS to verify if the introduction of the intermediate discontinuity reduces the displacements and drifts of the structure. The comparison between the displacements and drifts of the two models is made through three coefficients that are named gain coefficients and are defined as
α 1 = max [ u 1 ( t ) ] max [ u 1 A ( t ) ] ; α 2 = max [ u 2 ( t ) u d ( t ) ] max [ u 2 A ( t ) u 1 A ( t ) ] ; α 3 = max [ u d ( t ) u 1 ( t ) ] u 0 ,
where the displacements u 1 , u 2 , and u d refer to the MS (Figure 2a), u 1 A and u 2 A refer to the AS (Figure 2b), and u 0 is the threshold value of the acceptable displacements. The proposed gain coefficients represent a relative evaluation of the effectiveness of the hysteretic intermediate discontinuity. They evaluate the performances of the system with the discontinuity with respect to a reference system, represented by the corresponding linear elastic system without intermediate discontinuity.
The gain coefficient α 1 refers to the substructure and is the ratio between the displacements of the substructures of the MS and AS. Such a coefficient evaluates the performance of the superstructure as TMD for the substructure. The gain coefficient α 2 refers to the superstructure and is the ratio between the drifts of the superstructure of the MS and AS. In this case, α 2 evaluates the behaviour of the intermediate discontinuity as base isolation for the superstructure. Finally, α 3 is the ratio between the displacement of the intermediate discontinuity and the maximum threshold value u 0 = 15 cm. Such value was chosen in order to limit excessive displacements on the discontinuity layer that may jeopardize the safety of the structure. Each gain coefficient evaluates a different capability of the intermediate discontinuity. If a gain coefficient is less than unity, the introduction of the discontinuity is beneficial for a specific part of the structure. Therefore, the improvement of the dynamical behaviour of the frame means obtaining gain coefficients lower than unity (i.e., a reduction of displacements and drifts). A simultaneous control on these coefficients leads to a unitary approach to the dynamic enhancement of the structure following from the introduction of the discontinuity.
Two additional gain coefficients are introduced with the purpose to evaluate how the intermediate discontinuity affects the accelerations of the substructure and superstructure. They read
α 1 A = max [ u ¨ 1 ( t ) ] max [ u ¨ 1 A ( t ) ] ; α 2 A = max [ u ¨ 2 ( t ) u ¨ d ( t ) ] max [ u ¨ 2 A ( t ) u ¨ 1 A ( t ) ] ,
where u ¨ 1 ( t ) , u ¨ 2 ( t ) , and u ¨ d ( t ) are the accelerations of the MS. Similarly, u ¨ 2 A ( t ) and u ¨ 1 A ( t ) are the accelerations of the AS.

2.5. Variable Parameters

The parametric analysis is performed by varying the following parameters:
  • Post-yielding to pre-yielding stiffness ratio Ψ .
  • Yield displacement ratio η = u y u 0 .
  • Pre-yielding to storey stiffness ratio ρ = k d k p .
  • Ratio between m 2 (Equation (4)) and m ˜ 1 (Equation (5)) μ = m 2 m ˜ 1 .
  • Discontinuity level of the frame structure L.
The results presented in the following sections are obtained by numerically integrating the equations of motion. For this purpose, an original code is developed in the Mathematica® software. The equations of motion are integrated in an implicit form and the adopted step size is δ t = 0.01 s.

3. Harmonic Analysis

The harmonic analyses are conducted on the frame whose characteristics are reported in Table 1. Parametric analyses are carried out to identify how the parameters influence the response of the system with the discontinuity. In the analyses, the external excitation is defined as
a g ( t ) = A Sin Ω t ,
with A amplitude and Ω circular frequency of the excitation. Since the system is non-linear, A is considered as a variable parameter.

3.1. Frequency-Response Curves

A first analysis compares the frequency-response curves of the MS and AS. Such curves represent the dimensionless maximum displacement of the substructure, χ 1 = Max ( | u 1 ( t ) | / u 0 ), and the dimensionless drift of the superstructure, χ 2 = Max ( | u 2 ( t ) u d ( t ) | / u 0 ), as functions of the frequency ratio β = Ω / ω 1 , where ω 1 = k 1 / m 1 . In the analyses where the value of η is fixed, such value is assumed as η = 0.2 .
Figure 3 compares the frequency-response curves of the MS obtained for Ψ = 1.0, and the frequency-response curves of the AS. Assuming Ψ = 1.0 entails a linear behaviour of the discontinuity, the comparison shows that the discontinuity introduces a second peak in the frequency-response curves of the MS. Such a peak corresponds to a second resonant frequency, whereas the AS is characterized by a single resonant frequency in the range of the observed frequencies. The existence of a second resonant frequency is a consequence of the presence of a TMD in the system, meaning that, due to the discontinuity, the superstructure works as a TMD for the substructure. The frequency-response curve of χ 2 shows a shift of the main frequency of the MS to lower values and highlights a behaviour as a base isolated structure. Thus, the superstructure acts as a TMD for the substructure and the superstructure behaves as if it was base-isolated. In both cases, it can be asserted that the distance between the two resonant frequencies of the MS is a qualitative measure of the effectiveness of the discontinuity.
A second analysis investigates the sensitivity of the system to the values of the parameters characterizing the hysteretic behaviour of the disconnection devices. Figure 4 shows the frequency-response curves obtained for several values of the parameters characterizing the Bouc–Wen constitutive law. Specifically, Figure 4a shows the frequency-response curves of the MS for L = 3 , A = 0.6 g, ρ = 0.2 , and several values of Ψ increasing from Ψ = 0.10 to Ψ = 1.00 . The highest resonant frequency is not significantly sensitive to Ψ ; on the contrary, the smallest resonant frequency significantly depends on Ψ . The values of both χ 1 and χ 2 in correspondence of the resonant peaks strongly depend on Ψ . As Ψ decreases, the hysteresis cycles become larger, and the values of the dimensionless displacements in correspondence of the peak decrease. Moreover, the distance between the peaks increases as Ψ decreases. Figure 4b shows the frequency-response curves for L = 5 , A = 0.6 g, ρ = 0.2 , and the same values of Ψ of the case presented in Figure 4a. When the discontinuity is placed at higher levels, the values of χ 1 and χ 2 at the resonant peak that refers to the lower frequency reduce as Ψ decreases. The value of μ obtained for L = 5 is smaller than the value obtained for L = 3 , and the distance between the two resonant peaks decreases, thus reducing the frequency range where the discontinuity is effective. Figure 4c shows the frequency-response curves obtained for different values of ρ increasing from ρ = 0.02 to ρ = 0.30 , L = 3 , A = 0.6 g, and Ψ = 0.25 . As ρ increases, both the first and second resonant frequencies shift toward higher values. Moreover, χ 1 manifests higher values at the resonant peaks for smaller values of ρ . The same increase occurs for the values at the left peak of the χ 2 curves. Instead, the value of χ 2 at the higher resonant frequency decreases until it almost vanishes for decreasing values of ρ , thus confirming that, for small values of ρ , the discontinuity acts on the superstructure as a base isolation. The last sub-figure, Figure 4d, shows the frequency-response curves referring to the MS that are obtained for different values of A and assuming L = 3 , Ψ = 0.25 , and ρ = 0.2 . For increasing values of A the position of the peaks slightly shifts leftwards on the β -axis (i.e., the two resonant frequencies are only slightly sensitive to the value of A).

3.2. Gain Maps

This section analyses how Ψ affects the dynamical behaviour of the system with intermediate discontinuity. The effects of increasing values of Ψ are assessed by evaluating the gain coefficient Equations (8) and (9). The results of the parametric analysis are summarised in gain maps that represent the gain coefficients, α 1 and α 2 , plotted in the parameter plane η ρ . Each map refers to a specific value of Ψ . The maps on the left column of Figure 5 are obtained for Ψ = 0.25 , whereas the maps on the right column refer to Ψ = 0.75 . Moreover, the maps of Figure 5a show the α 1 coefficient, whereas the maps of Figure 5b show the α 2 coefficient. The maps are obtained for β = 0.75 and A = 0.1 g; they show a strong dependence on ρ and, up to a threshold value, on η . At higher values of η (where the value of η depends on the considered map), the gain coefficients become independent of η (i.e., in both the α 1 and α 2 maps, the contour levels tend to become horizontal). It is useful to remark that η directly depends on the yielding displacement; hence, high values of η refer to high values of u y . A frame with intermediate discontinuity and small values of η manifests a stronger non-linear behaviour of the disconnection devices than systems with higher η . For very high value of η , the disconnection devices tend to manifest a linear elastic behaviour since the high-yielding displacement u y cannot be reached during the motion. The impossibility to reach u y during the motion is more evident at higher values of the stiffness ratio Ψ , due to the higher post-elastic stiffness of the disconnection devices. For the chosen values of β , A, and Ψ , both α 1 and α 2 are smaller than unity in the whole parameter plane, which means that the intermediate discontinuity reduces the displacement of the substructure and the drift of the superstructure with respect to those of the original frame structure.
Under a harmonic excitation, after a transient time that depends on the damping of the structure, the system undergoes stationary dynamics. To understand the behaviour of the disconnection devices for increasing values of η , the stationary hysteretic cycles are shown for points P 1 and P 3 in Figure 5c, and for points P 2 and P 4 in Figure 5d. The hysteretic cycles show the value of the dimensionless force σ h = F h / k d u y as a function of the dimensionless displacement χ d = u d ( t ) / u y . By comparing the cycles in the first column of Figure 5c,d, which refer to Ψ = 0.25 , with the cycles in the second column, which are obtained for Ψ = 0.75 , it can be observed that a higher value of Ψ corresponds to smaller hysteresis cycles. Similarly, an increase of η also leads to a reduction of the area of the hysteretic cycles (compare the cycles of Figure 5c,d).
Figure 6 shows the gain maps obtained for β = 1 and Ψ = 0.25 . The maps in Figure 6a and Figure 6b are obtained for A = 0.1 g and A = 0.6 g, respectively, assuming that the discontinuity is located at the third level ( L = 3 ). Figure 6c, instead, shows maps obtained for A = 0.1 g and a discontinuity located at the fifth level ( L = 5 ). The comparison of Figure 6a and Figure 6b highlights the dependence of the gain maps on A. Specifically, in the gain maps of Figure 6b, the contour levels do not tend to a horizontal asymptote because the displacement of the disconnection devices exceeds the yield displacement u y due to the high value of A, thus assuring wide hysteresis cycles.
The α 1 gain maps of Figure 6a,b show that there are large regions of the parameter where the gain coefficients are higher than unity (dark grey regions). In these regions, the discontinuity is not able to reduce the displacement of the substructure. On the contrary, α 2 is smaller than unity in the whole parameter plane, thus assuring a reduction of the superstructure drift for any combination of the parameters. Nevertheless, it is always possible to choose a combination of the design parameters ρ and η inside regions where both the gain coefficients α 1 and α 2 are smaller than unity.
The comparison among the left maps of Figure 5a,b and the maps of Figure 6a (all obtained for Ψ = 0.25 and A = 0.1 g) shows that the effectiveness of the discontinuity (i.e., the smallness of the gain coefficients) strongly depends on the frequency ratio β .
Finally, Figure 6c shows the gain maps for L = 5 . In this case, the dark grey regions where α 1 is higher than unity reduce compared to the case with L = 3 , and, generally, smaller values of both α 1 and α 2 are obtained, thus showing the higher effectiveness of the intermediate discontinuity.

4. Seismic Analysis

This section shows the results of a seismic analysis. The first part of the section aims at providing an initial understanding of the seismic behaviour of the MS. For this purpose, the gain maps show the results of simulations performed using single seismic records as excitation. The two different seismic records considered in this analysis are:
(a)
El Centro, CA, Array Station 9, Imperial Valley Irrigation District, component 180;
(b)
Parkfield, CO2-065 ground motion recorded during the California earthquake 1966.
Figure 7 shows the time-histories (left graphs) and the pseudo-response spectra (right graphs) of the two seismic records. In the following, each seismic record is called with the corresponding earthquake name (i.e., the names underlined in the list above).
Since gain maps referred to a single registration are not useful for a complete assessment of the effectiveness of the intermediate discontinuity, the second part of the section shows the results of an analysis carried out considering n e = 305 seismic records. The complete list of seismic records used in the analysis can be found at http://ing.univaq.it/diegidio/Earthq/Earthquakes_Table.pdf. Additionally, in this case, the results are presented in gain maps. Such gain maps show the mean values of the gain indexes obtained with the different seismic records.

4.1. Gain Maps from a Single Seismic Record

Figure 8, Figure 9 and Figure 10 show the gain maps obtained with the seismic records of El Centro and Parkfield for L = 1 , L = 3 , and L = 5 , respectively. Specifically Figure 8a–c, Figure 9a–c and Figure 10a–c are obtained for Ψ = 0.25 , whereas Figure 8b–d, Figure 9b–d, and Figure 10b–d for Ψ = 0.75 . The regions where the gain coefficients α 1 and α 2 are smaller than unity are coloured, those where these gain coefficients are higher than unity are displayed in dark grey colour. Both the α 1 and α 2 maps show that, as Ψ increases, the gain coefficients increase, thus reducing the effectiveness of the intermediate discontinuity. The α 3 maps show the feasibility of the protection strategy. When α 3 is higher than unity (dark grey zones of the α 3 maps), the displacements of the disconnection devices are higher than the permissible threshold value. The size of the dark grey zones in the α 3 maps increases when Ψ decreases, thus even if lower values of Ψ assure a higher reduction of displacements and drifts, the displacements of the discontinuity level may not be compatible with the capability of the disconnection devices.
In Figure 8, Figure 9 and Figure 10, the thick black lines on the maps represent the admissible limit of α 3 . In each map, only in the region above such lines, the displacement of the discontinuity is below the threshold value u 0 . While the α 3 gain maps obtained from the seismic record of El Centro show only small dark grey regions, the α 3 gain maps obtained from the Parkfield seismic record show large regions where α 3 > 1 maps. The size of these regions significantly reduces as Ψ increases. Moreover, when the discontinuity level L increases, the size of the dark grey regions tends to reduce.
Figure 11 and Figure 12 show the α 1 A and α 2 A gain maps under El Centro and Parkfield seismic records, respectively. These maps evaluate the effects of the discontinuity on the accelerations of the structure. The black lines on the maps show the range highlighted in the α 3 maps of Figure 8, Figure 9 and Figure 10. Only in the regions above such lines are the displacements of the disconnection devices below the admissible threshold. In Figure 11 and Figure 12, the maps in the first two rows refer to α 1 A , whereas the maps in the third and fourth rows refer to α 2 A . The maps in the first and third rows are obtained for Ψ = 0.25 , the maps in the second and fourth rows refer to Ψ = 0.75 .
The α 1 A maps show that as Ψ increases, α 1 A decreases. Similarly, α 2 A decreases as Ψ increases. Therefore, a lower value of Ψ may increase the accelerations and, at the same time, reduce displacements and drifts. However, both the gain coefficients based on the displacements and drifts (Equation (8)) and those based on the accelerations (Equation (9)) are smaller than unity in a large area of the parameters plane. It is worth observing that α 1 A and α 2 A are lower than α 1 and α 2 , thus the use of α 1 and α 2 provides a more severe assessment of the effectiveness of the intermediate discontinuity.

4.2. Gain Maps Obtained from Multiple Seismic Records

The seismic analysis is carried out by using several seismic records obtained from [32]. The magnitude varies from 4 to 8 and the epicentral distance varies from 10 to 30 km. In the analysis, the target parameters are the mean displacements u ¯ i that are defined as the mean of the maximum displacements u i , j obtained from the integration of the equations of motion with each seismic record; they read
u ¯ i = 1 n e j = 1 n e max [ u i , j ( t ) ]
with ( i = 1 , d , 2 ) . The displacements u 1 , j , u d , j , and u 2 , j are those of the substructure, discontinuity level, and superstructure obtained from the jth seismic record. The total number of seismic records used in the simulations is n e . The number of seismic records used is sufficient to reach a coefficient of variation (CoV) of u ¯ i , CoV u ¯ i 0.2 . Therefore, the predictability of the model to identify the average displacement is higher than 80%. The CoV u ¯ i is defined as
C o V u ¯ i = σ u ¯ i μ ( u ¯ i )
where σ ( u ¯ i ) and μ ( u ¯ i ) are the standard deviation and mean value of u ¯ i , respectively. To understand the effectiveness of the disconnection, three coefficients (called mean-gain coefficients in the following) similar to those in Equation (8) are introduced:
α 1 ¯ = j = 1 n e max [ u 1 , j ( t ) ] j = 1 n e max [ u 1 , j A ( t ) ] ; α 2 ¯ = j = 1 n e max [ u 2 , j ( t ) u d , j ( t ) ] j = 1 n e max [ u 2 , j A ( t ) u 1 , j A ( t ) ] ; α 3 ¯ = j = 1 n e max [ u d , j ( t ) u 1 , j ( t ) ] n e u 0
Figure 13 shows the gain maps of the mean displacements. Such maps represent the mean-gain coefficients in the parameters plane P G A - ρ , where P G A is the peak ground acceleration and ρ is the stiffness ratio previously defined. All the gain maps are obtained for η = 0.2 . Similarly to previous figures, the regions where α ¯ i < 1 , ( i = 1 , 2 ) are named gain regions and are coloured, whereas those where α ¯ i > 1 are identified by a dark grey colour. The first, second, and third rows of the figure refer to different levels of the discontinuity, L = 1 , L = 3 , and L = 5 , respectively. All maps are obtained for Ψ = 0.25 .
The maps show wide gain regions. As can be observed, L = 1 provides the best performances in the reduction of the displacement of the substructure because the values of α ¯ 1 obtained for L = 1 are in general smaller than those obtained for other values of L. On the contrary, L = 5 provides the largest reduction of the drift of the superstructure since, in this case, α ¯ 2 is smaller than those in the other cases. The maps of α ¯ 3 show values higher than unity only in the zone with high values of P G A and low values of ρ , specifically in the region where P G A > 1.1 g and ρ < 0.1 . In this region, α ¯ 1 and α ¯ 2 have small values, but the displacement of the discontinuity level exceeds the threshold value u o (i.e., α ¯ 3 is higher than unity). Nevertheless, such a region is very small and a slight increase of ρ is sufficient to move in a region of the maps where α ¯ 3 1 and at the same α ¯ 1 1 and α ¯ 2 1 . Since all the three coefficients are less than unity the disconnection can be considered a successful protection strategy for the original frame structure. The map obtained for L = 3 shows an intermediate behaviour between those obtained for L = 1 and L = 5 , and the map of α ¯ 3 does not show any dark grey regions.
To summarize, the α ¯ 1 maps exhibit a decrease in the performance from L = 1 to L = 5 . This can be explained by assuming that the superstructure acts as a TMD for the substructure and consequently the performance decreases because the mass ratio between the superstructure and substructure decreases [7,22]. Instead, α ¯ 2 decreases from L = 1 to L = 5 . In this case, when L = 5 , the superstructure is very stiff and the superstructure acts as if it was base-isolated. On the contrary, when L = 1 , the superstructure has lower stiffness and the performance decreases since the period of the superstructure can approach that of the discontinuity layer. Similar results were already presented in [33] for tower buildings.
Figure 14 shows the mean-gain maps obtained for L = 3 . The maps in the first row refer to Ψ = 0.1 , whereas the maps in the second row are obtained for Ψ = 0.75 . A direct comparison with the maps in the second row of Figure 13 ( L = 3 ) allows evaluating the effect of Ψ on the performances of the intermediate discontinuity. The dark grey region of the α ¯ 1 map enlarges with respect to that of Figure 13 for increasing Ψ . Therefore, lower values of Ψ allow a wider range of use. Moreover, when Ψ decreases, the performance of the system increases since there is a higher reduction of the displacement of the substructure (i.e., a decrease of α 1 ). On the contrary, the coefficients α ¯ 2 and α ¯ 3 are not significantly affected by Ψ . Finally, α ¯ 3 is always smaller than unity when Ψ < 0.25 , whereas for Ψ = 0.75 there is a small region where α ¯ 3 1 .

5. Conclusions

This paper proposes a modelling approach for the study of the intermediate discontinuity. The structure with the discontinuity layer has been modelled by a 3-DOF model and the disconnection devices have been described by using the Bouc–Wen model. In order to evaluate the improvement of the dynamical behaviour obtained by the introduction of the discontinuity, two different systems have been analysed. The first system is the frame after the introduction of the discontinuity, called the main system (MS) and the second system is the reference frame, called the auxiliary system (AS).
The first part of the paper is focused on the dynamical behaviour of the MS that is analysed by using harmonic excitation. The results of the analysis have been organized in frequency-response curves and gain maps. It has been found that the value of the post-yielding to pre-yielding stiffness ratio Ψ in the Bouc–Wen model significantly affects the value of the first frequency of the MS. Moreover, by comparing the frequency-response curves of the AS and MS for different values of Ψ , it has been found that lower values of Ψ may extend the range of frequency where the intermediate discontinuity is effective. The harmonic maps highlighted that small values of Ψ may reduce the displacement of the superstructure, but also increase the displacement of the substructure.
The second part of the paper is focused on the seismic response of the models. A preliminary analysis has been carried out considering single seismic records. Five gain coefficients have been defined to evaluate the effects of the introduction of the intermediate discontinuity. Such coefficients have been defined as ratios between the displacements (or the accelerations) of the MS and those of the AS. The results of the analysis have been summarized in gain maps. Such maps show the gain coefficients in a specific parameters plane. The results also show that an increase of Ψ reduces the displacement of the substructure and increases that of the superstructure. Moreover, higher values of Ψ reduce the displacement of the discontinuity layer, thus reducing the likelihood of such displacement exceeding a threshold value. It has been found that another possibility to reduce the displacement of the discontinuity layer is to change the level of the discontinuity L. As L increases, the displacement of the discontinuity layer decreases due to the filtering effect of the substructure. Additionally, the effects of the intermediate discontinuity on the accelerations acting on the structure have been evaluated, and it has been found that when Ψ decreases, the accelerations may increase.
Finally, a seismic analysis is carried out considering several seismic records within definite ranges of magnitude and epicentral distance. The results have been organized in gain maps that show the gain coefficients in different parameters planes. Such maps may be a useful auxiliary tool in the understanding of the behaviour resulting from the introduction of an intermediate discontinuity in a frame structure.

Author Contributions

Conceptualization, A.D.E., S.P. and A.C.; Methodology, A.D.E. and A.C.; Software, S.P.; Formal analysis, S.P.; Investigation, A.D.E., S.P. and A.C.; Writing—original draft, S.P.; Writing—review & editing, A.D.E. and A.C.; Supervision, A.D.E. and A.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by CARISPAQ Foundation, year 2022.

Data Availability Statement

Data, models, and code used during the study are available from the corresponding author by request. Specifically, the following data are available: (i) data of the results, and ( i i ) the source code for the numerical integration procedure.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kelly, J. Base isolation: Linear theory and design. Earthq. Spectra 1995, 6, 223–244. [Google Scholar] [CrossRef]
  2. Den Hartog, J.P. Mechanical Vibrations, 4th ed.; McGraw-Hill: New York, NY, USA, 1956. [Google Scholar]
  3. Tsai, H.H. The effect of tuned-mass dampers on the seismic response of base-isolated structures. Int. J. Solids Struct. 1995, 32, 1195–1210. [Google Scholar] [CrossRef]
  4. Taniguchi, T.; Kiureghian, A.D.; Melkumyan, M. Effect of tuned mass damper on displacement demand of base-isolated structures. Eng. Struct. 2008, 30, 3478–3488. [Google Scholar] [CrossRef]
  5. De Domenico, D.; Ricciardi, G. An enhanced base isolation system equipped with optimal Tuned Mass Damper Inerter (TMDI). Earthq. Eng. Struct. Dyn. 2018, 47, 1169–1192. [Google Scholar] [CrossRef]
  6. Giaralis, A.; Taflanidis, A. Optimal tuned mass-damper-inerter (TMDI) design for seismically excited MDOF structures with model uncertainties based on reliability criteria. Struct. Control Health Monit. 2018, 25, e2082. [Google Scholar] [CrossRef]
  7. De Angelis, M.; Perno, S.; Reggio, A. Dynamic response and optimal design of structures with large mass ratio TMD. Earthq. Eng. Struct. Dyn. 2012, 41, 41–60. [Google Scholar] [CrossRef]
  8. Chey, M.; Chase, J.; Mander, J.B.; Carr, A. Innovative seismic retrofitting strategy of added stories isolation system. Front. Struct. Civ. Eng. 2013, 7, 13–23. [Google Scholar] [CrossRef]
  9. Wang, S.J.; Chang, K.C.; Hwang, J.S.; Lee, B.H. Simplified analysis of mid-story seismically isolated buildings. Earthq. Eng. Struct. Dyn. 2011, 40, 119–133. [Google Scholar] [CrossRef]
  10. Wang, S.J.; Chang, K.C.; Hwang, J.S.; Hsiao, J.Y.; Lee, B.H.; Hung, Y.C. Dynamic behavior of a building structure tested with base and mid-story isolation systems. Eng. Struct. 2012, 42, 420–433. [Google Scholar] [CrossRef]
  11. Wang, S.J.; Hwang, J.S.; Chang, K.C.; Lin, M.H.; Lee, B.H. Analytical and experimental studies on midstory isolated buildings with modal coupling effect. Earthq. Eng. Struct. Dyn. 2013, 42, 201–219. [Google Scholar] [CrossRef]
  12. Villaverde, R.; Mosqueda, G. Aseismic roof isolation system: Analytic and shake table studies. Earthq. Eng. Struct. Dyn. 1999, 28, 217–234. [Google Scholar] [CrossRef]
  13. Villaverde, R.; Aguirre, M.; Hamilton, C. Aseismic Roof Isolation System Built with Steel Oval Elements. Earthq. Spectra 2005, 21, 225–241. [Google Scholar] [CrossRef]
  14. Matta, E.; De Stefano, A. Robust design of mass-uncertain rolling-pendulum TMDs for the seismic protection of buildings. Mech. Syst. Signal Process. 2009, 23, 127–147. [Google Scholar] [CrossRef]
  15. Feng, M.; Mita, A. Vibration Control of Tall Buildings Using Mega SubConfiguration. J. Eng. Mech. 1995, 121, 1082–1088. [Google Scholar] [CrossRef]
  16. Sadek, F.; Mohraz, B.; Taylor, A.W.; Chung, R.M. A method of estimating the parameters of tuned mass dampers for seismic applications. Earthq. Eng. Struct. Dyn. 1997, 26, 617–635. [Google Scholar] [CrossRef]
  17. Anajafi, H.; Medina, R. Comparison of the seismic performance of a partial mass isolation technique with conventional TMD and base-isolation systems under broad-band and narrow-band excitations. Eng. Struct. 2018, 158, 110–123. [Google Scholar] [CrossRef]
  18. Rana, R.; Soong, T. Parametric analysis and simplified design of tuned mass damper. Eng. Struct. 1998, 20, 193–204. [Google Scholar] [CrossRef]
  19. Miranda, J.C. On tuned mass dampers for reducing the seismic response of structures. Earthq. Eng. Struct. Dyn. 2005, 34, 847–865. [Google Scholar] [CrossRef]
  20. Fabrizio, C.; Di Egidio, A.; de Leo, A. Top disconnection versus base disconnection in structures subjected to harmonic base excitation. Eng. Struct. 2017, 152, 660–670. [Google Scholar] [CrossRef]
  21. Fabrizio, C.; de Leo, A.; Di Egidio, A. Tuned mass damper and base isolation: A unitary approach for the seismic protection of conventional frame structures. J. Eng. Mech. 2019, 145, 04019011. [Google Scholar] [CrossRef]
  22. Reggio, A.; De Angelis, M. Optimal energy-based seismic design of non-conventional Tuned Mass Damper (TMD) implemented via inter-story isolation. Earthq. Eng. Struct. Dyn. 2015, 44, 1623–1642. [Google Scholar] [CrossRef]
  23. Pagliaro, S.; Di Egidio, A. Archetype Dynamically Equivalent 3-DOF Model to Evaluate Seismic Performances of Intermediate Discontinuity in Frame Structures. J. Eng. Mech. 2022, 148, 04022004. [Google Scholar] [CrossRef]
  24. De Domenico, D.; Ricciardi, G. Optimal design and seismic performance of tuned mass damper inerter (TMDI) for structures with nonlinear base isolation systems. Earthq. Eng. Struct. Dyn. 2018, 47, 2539–2560. [Google Scholar]
  25. Wang, S.J.; Lin, M.H.; Chiang, Y.S.; Hwang, J.S. Mechanical behavior of lead rubber bearings under and after nonproportional plane loading. Earthq. Eng. Struct. Dyn. 2019, 48, 1508–1531. [Google Scholar] [CrossRef]
  26. Song, J.; Der Kiureghian, A. Generalized Bouc-Wen model for highly asymmetric hysteresis. J. Eng. Mech. 2006, 132, 610–618. [Google Scholar] [CrossRef]
  27. Aloisio, A.; Alaggio, R.; Kohler, J.; Fragiacomo, M. Extension of Generalized Bouc-Wen Hysteresis Modeling of Wood Joints and Structural Systems. J. Eng. Mech. 2020, 146, 04020001. [Google Scholar] [CrossRef]
  28. Wen, Y.K. Method for random vibration of hysteretic systems. J. Eng. Mech. 1976, 102, 249–263. [Google Scholar] [CrossRef]
  29. Charalampakis, A. Parameters of Bouc-Wen hysteretic model revisited. In Proceedings of the 9th HSTAM International Congress on Mechanics, Limassol, Cyprus, 12–14 July 2010. [Google Scholar]
  30. Ma, F.; Zhang, H.; Bockstedte, A.; Foliente, G.; Paevere, P. Parameter Analysis of the Differential Model of Hysteresis. J. Eng. Mech. 2004, 71, 342–349. [Google Scholar] [CrossRef]
  31. Constantinou, M.; Adnane, M. Dynamics of Soil-Base-Isolated Structure Systems: Evaluation of Two Models for Yielding Systems; Report to NSAF; Department of Civil Engineering, Drexel University: Philadelphia, PA, USA, 1987. [Google Scholar]
  32. Pacific Earthquake Engineering Research Center (PEER). Strong Ground Motion Database. Available online: https://peer.berkeley.edu. (accessed on 6 November 2021).
  33. Ziyaeifar, M.; Noguchi, H. Partial Mass Isolation in Tall Buildings. Earthq. Eng. Struct. Dyn. 1998, 27, 49–65. [Google Scholar] [CrossRef]
Figure 1. Structural scheme of the (a) frame structure with intermediate discontinuity; (b) original frame structure (without intermediate discontinuity).
Figure 1. Structural scheme of the (a) frame structure with intermediate discontinuity; (b) original frame structure (without intermediate discontinuity).
Applsci 13 05373 g001
Figure 2. Archetype models: (a) MS; (b) AS (HD: Elasto-plastic behaviour of hysteretic devices).
Figure 2. Archetype models: (a) MS; (b) AS (HD: Elasto-plastic behaviour of hysteretic devices).
Applsci 13 05373 g002
Figure 3. Frequency-response curves: comparison between curves of AS and MS ( L = 3 and Ψ = 1.00 ).
Figure 3. Frequency-response curves: comparison between curves of AS and MS ( L = 3 and Ψ = 1.00 ).
Applsci 13 05373 g003
Figure 4. Frequency-response curves: (a) L = 3 , A = 0.6 g, ρ = 0.2 , and different values of the parameter ψ ; (b) L = 5 , A = 0.6 g, ρ = 0.2 , and different values of the parameter ψ ; (c) L = 3 , A = 0.6 g, Ψ = 0.25 , and different values of ρ ; and (d) L = 3 , Ψ = 0.25 , ρ = 0.2 , and different values of A.
Figure 4. Frequency-response curves: (a) L = 3 , A = 0.6 g, ρ = 0.2 , and different values of the parameter ψ ; (b) L = 5 , A = 0.6 g, ρ = 0.2 , and different values of the parameter ψ ; (c) L = 3 , A = 0.6 g, Ψ = 0.25 , and different values of ρ ; and (d) L = 3 , Ψ = 0.25 , ρ = 0.2 , and different values of A.
Applsci 13 05373 g004
Figure 5. Gain maps for fixed β = 0.75 and A = 0.1 g: (a) α 1 maps for two different values of Ψ ; (b) α 2 maps for two different values of Ψ ; (c) hysteretic cycles in P 1 ( 0.025 , 0.15 ) , Ψ = 0.25 and P 3 ( 0.025 , 0.15 ) , Ψ = 0.75 ; (d) hysteretic cycles in P 2 ( 0.15 , 0.15 ) , Ψ = 0.25 and P 4 ( 0.15 , 0.15 ) , Ψ = 0.75 .
Figure 5. Gain maps for fixed β = 0.75 and A = 0.1 g: (a) α 1 maps for two different values of Ψ ; (b) α 2 maps for two different values of Ψ ; (c) hysteretic cycles in P 1 ( 0.025 , 0.15 ) , Ψ = 0.25 and P 3 ( 0.025 , 0.15 ) , Ψ = 0.75 ; (d) hysteretic cycles in P 2 ( 0.15 , 0.15 ) , Ψ = 0.25 and P 4 ( 0.15 , 0.15 ) , Ψ = 0.75 .
Applsci 13 05373 g005
Figure 6. Gain maps for β = 1 : (a) L = 3 , A = 0.1 g, and Ψ = 0.25 ; (b) L = 3 , A = 0.6 g, and Ψ = 0.25 ; (c) L = 5 , A = 0.1 g, and Ψ = 0.25 .
Figure 6. Gain maps for β = 1 : (a) L = 3 , A = 0.1 g, and Ψ = 0.25 ; (b) L = 3 , A = 0.6 g, and Ψ = 0.25 ; (c) L = 5 , A = 0.1 g, and Ψ = 0.25 .
Applsci 13 05373 g006
Figure 7. Recorded earthquakes: (a) El Centro seismic record; and (b) Parkfield seismic record.
Figure 7. Recorded earthquakes: (a) El Centro seismic record; and (b) Parkfield seismic record.
Applsci 13 05373 g007
Figure 8. Gain maps of α 1 , α 2 , and α 3 for L = 1 : (a) El Centro seismic record and Ψ = 0.25 (b) El Centro seismic record and Ψ = 0.75 ; (c) Parkfield Seismic record and Ψ = 0.25 ; and (d) Parkfield Seismic record and Ψ = 0.75 .
Figure 8. Gain maps of α 1 , α 2 , and α 3 for L = 1 : (a) El Centro seismic record and Ψ = 0.25 (b) El Centro seismic record and Ψ = 0.75 ; (c) Parkfield Seismic record and Ψ = 0.25 ; and (d) Parkfield Seismic record and Ψ = 0.75 .
Applsci 13 05373 g008
Figure 9. Gain maps of α 1 , α 2 , and α 3 for L = 3 : (a) El Centro seismic record and Ψ = 0.25 (b) El Centro seismic record and Ψ = 0.75 ; (c) Parkfield Seismic record and Ψ = 0.25 ; and (d) Parkfield Seismic record and Ψ = 0.75 .
Figure 9. Gain maps of α 1 , α 2 , and α 3 for L = 3 : (a) El Centro seismic record and Ψ = 0.25 (b) El Centro seismic record and Ψ = 0.75 ; (c) Parkfield Seismic record and Ψ = 0.25 ; and (d) Parkfield Seismic record and Ψ = 0.75 .
Applsci 13 05373 g009
Figure 10. Gain maps of α 1 , α 2 , and α 3 for L = 5 : (a) El Centro seismic record and Ψ = 0.25 (b) El Centro seismic record and Ψ = 0.75 ; (c) Parkfield Seismic record and Ψ = 0.25 ; and (d) Parkfield Seismic record and Ψ = 0.75 .
Figure 10. Gain maps of α 1 , α 2 , and α 3 for L = 5 : (a) El Centro seismic record and Ψ = 0.25 (b) El Centro seismic record and Ψ = 0.75 ; (c) Parkfield Seismic record and Ψ = 0.25 ; and (d) Parkfield Seismic record and Ψ = 0.75 .
Applsci 13 05373 g010
Figure 11. Gain maps of α 1 A and α 2 A for L = 1 , L = 3 and L = 5 considering the El Centro seismic record as external excitation: (a) maps of α 1 A obtained for Ψ = 0.25 ; (b) maps of α 1 A obtained for Ψ = 0.75 ; (c) maps of α 2 A obtained for Ψ = 0.25 ; and (d) maps of α 2 A obtained for Ψ = 0.75 .
Figure 11. Gain maps of α 1 A and α 2 A for L = 1 , L = 3 and L = 5 considering the El Centro seismic record as external excitation: (a) maps of α 1 A obtained for Ψ = 0.25 ; (b) maps of α 1 A obtained for Ψ = 0.75 ; (c) maps of α 2 A obtained for Ψ = 0.25 ; and (d) maps of α 2 A obtained for Ψ = 0.75 .
Applsci 13 05373 g011
Figure 12. Gain maps of α 1 A and α 2 A for L = 1 , L = 3 and L = 5 considering the Parkfield seismic record as external excitation: (a) maps of α 1 A obtained for Ψ = 0.25 ; (b) maps of α 1 A obtained for Ψ = 0.75 ; (c) maps of α 2 A obtained for Ψ = 0.25 ; and (d) maps of α 2 A obtained for Ψ = 0.75 .
Figure 12. Gain maps of α 1 A and α 2 A for L = 1 , L = 3 and L = 5 considering the Parkfield seismic record as external excitation: (a) maps of α 1 A obtained for Ψ = 0.25 ; (b) maps of α 1 A obtained for Ψ = 0.75 ; (c) maps of α 2 A obtained for Ψ = 0.25 ; and (d) maps of α 2 A obtained for Ψ = 0.75 .
Applsci 13 05373 g012
Figure 13. Mean gain maps for different values of L = 1 , 3 , 5 and Ψ = 0.25 .
Figure 13. Mean gain maps for different values of L = 1 , 3 , 5 and Ψ = 0.25 .
Applsci 13 05373 g013
Figure 14. Mean gain maps for L = 3 and (a) Ψ = 0.10 ; (b) Ψ = 0.75 .
Figure 14. Mean gain maps for L = 3 and (a) Ψ = 0.10 ; (b) Ψ = 0.75 .
Applsci 13 05373 g014
Table 1. Mechanical characteristics of MS and AS archetypes.
Table 1. Mechanical characteristics of MS and AS archetypes.
L m 1 [kg × 10 3 ] m 1 ˜ [kg × 10 3 ] m 2 [kg × 10 3 ] k 1 [kN/m] k 2 [kN/m]
1301.5150.751507.50602,403182,552
2603.0425.251206.00292,175208,839
3904.5735.75904.50207,567250,374
412061055.25603.00177,681324,132
51507.51356.75301.50167,962476,779
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Di Egidio, A.; Pagliaro, S.; Contento, A. Seismic Performance of Frame Structure with Hysteretic Intermediate Discontinuity. Appl. Sci. 2023, 13, 5373. https://doi.org/10.3390/app13095373

AMA Style

Di Egidio A, Pagliaro S, Contento A. Seismic Performance of Frame Structure with Hysteretic Intermediate Discontinuity. Applied Sciences. 2023; 13(9):5373. https://doi.org/10.3390/app13095373

Chicago/Turabian Style

Di Egidio, Angelo, Stefano Pagliaro, and Alessandro Contento. 2023. "Seismic Performance of Frame Structure with Hysteretic Intermediate Discontinuity" Applied Sciences 13, no. 9: 5373. https://doi.org/10.3390/app13095373

APA Style

Di Egidio, A., Pagliaro, S., & Contento, A. (2023). Seismic Performance of Frame Structure with Hysteretic Intermediate Discontinuity. Applied Sciences, 13(9), 5373. https://doi.org/10.3390/app13095373

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop