On the Use of Cloud Analysis for Structural Glass Members under Seismic Events

: Current standards for seismic-resistant buildings provide recommendations for various structural systems, but no speciﬁc provisions are given for structural glass. As such, the seismic design of joints and members could result in improper sizing and non-efﬁcient solutions, or even non-efﬁcient calculation procedures. An open issue is represented by the lack of reliable and general-ized performance limit indicators (or “engineering demand parameters”, EDPs) for glass structures, which represent the basic input for seismic analyses or q -factor estimates. In this paper, special care is given to the q -factor assessment for glass frames under in-plane seismic loads. Major advantage is taken from efﬁcient ﬁnite element (FE) numerical simulations to support the local/global analysis of mechanical behaviors. From extensive non-linear dynamic parametric calculations, numerical outcomes are discussed based on three different approaches that are deeply consolidated for ordinary structural systems. Among others, the cloud analysis is characterized by high computational efﬁciency, but requires the deﬁnition of speciﬁc EDPs, as well as the choice of reliable input seismic signals. In this regard, a comparative parametric study is carried out with the support of the incremental dynamic analysis (IDA) approach for the herein called “dynamic” (M1) and “mixed” (M2) procedures, towards the linear regression of cloud analysis data (M3). Potential and limits of selected calculation methods are hence discussed, with a focus on sample size, computational cost, estimated mechanical phenomena, and predicted q -factor estimates for a case study glass frame.


Introduction
The large use of glass structures in civil engineering applications represents a challenging issue for designers. In addition to intrinsic mechanical features of the involved load-bearing materials [1,2], careful consideration should be paid in earthquake-prone regions to satisfy rigid resistance and displacement demands. This is the case of primary, stand-alone glass structures, but also secondary glass systems belonging to different primary buildings and constructional assemblies [3][4][5][6][7][8].
According to various literature studies, the seismic capacity of glass structures can benefit from innovative tools and special fasteners [9][10][11]. At the component level, refined calculation approaches and investigations of literature have been dedicated to both the pre-and post-cracked analysis of laminated glass (LG) elements [12][13][14], including considerations of their residual strength [15]. In any case, glass structures are still a rather new domain for several professional designers, and certainly require dedicated methods of analysis [16]. Among others, an open issue is represented by the seismic design of glass structures. Most of the available technical documents do not provide specific recommendations for glass [17,18], but suggest the use of "reliable calculation methods" to verify the seismic resistance/displacement capacity of glass components and restraints. Such a the seismic resistance/displacement capacity of glass components and restraints. Such a technical difficulty is further enforced by the need for a realistic calibration of the expected q-factor [19].
The main goal of present study, in this regard, is to assess the sensitivity of q-factor for glass structures based on simplified or more advanced calculation approaches. As shown in Section 2, consolidated strategies are common for conventional constructional typologies/materials. Moreover, the q-factor itself (with q > 1) is known to represent the intrinsic dissipation capacity of the structure/material to verify. At the same time, established performance indicators (or "engineering demand parameters", EDPs) in support of seismic analysis and design are available in literature for structural members and systems made of steel, reinforced concrete, timber, or masonry, while such a calibration is missing for glass. To summarize the present discussion, the numerical analysis is focused on a case study glass frame that was earlier investigated in [19]. Differing from [19], however, the attention was given to the seismic performance and capacity of the full-size frame, rather than its key base connections only. To this aim, an original finite element (FE) numerical model was developed and optimized to support the local/global analysis of the frame as a whole. Extended sets of non-linear dynamic analyses were in fact carried out for the frame under in-plane seismic lateral loads. In doing so, three selected methods of analysis that are deeply consolidated for ordinary constructions (M1 to M3 in Section 3) were adapted to the examined structural glass frame and assessed for the q-factor prediction. Basic comparative calculations were first carried out with the support of the incremental dynamic analysis (IDA) approach for the herein called "dynamic" (M1) and "mixed" (M2) procedures. The cloud analysis procedure (M3), as shown, is characterized by high efficiency compared to M1 and M2 methods, but requires the calibration of specific EDPs for glass, as well as an accurate selection of input signals for the structural system to verify. FE comparative results are thus discussed in  showing the potential and limits of selected M1 to M3 calculation methods, in support of a realistic and computationally efficient estimation of seismic behavioral trends for glass structures, thus resulting in their optimized structural design.

State-of-Art and Literature Review on q-factor Methods
Following EC8 [17], the seismic design of buildings is today conducted by using the so-called force-based design (FBD) method. The design base shear is conventionally obtained as the ratio between the elastic base shear and the q-factor of the structure to verify ( Figure 1). The q-factor introduction, as such, simplifies its complex energy dissipation capacity (by means of plastic deformations) to a linear elastic model. Due to its strategic role, the q-factor definition is thus a topic which has been deeply discussed in the seismic engineering field as a primary focus of several studies since the 1950s.  A first simple formulation was proposed for the q-factor in the 1980s [20]. Further, it was first recognized by the modern design strategy that structures able to resist severe earthquakes are expected to experience permanent damage. As a matter of fact, design Sustainability 2021, 13, 9291 3 of 25 seismic actions are scaled by taking advantage of an intrinsic plastic capacity that is correlated to irreversible deformations.
Typical q-factor values by EC8 are known to span in the range of 1.5 (inverted pendulum systems), 2.0 (torsionally flexible systems), and 3.0 (frame systems), or even higher. Due to this, these constant EC8 values should be generally treated as an upper bound, and thus moving the decision on the acceptable level of damage becomes a designer responsibility. This decision is directly affected by structural features, details, and material properties. In such a general discussion, the choice of standards to adopt constant q-factor values looks very conservative. Several literature studies proved that the dissipative capacity of a structure is generally greater than the recommended limit values. For example, as concerns steel moment resisting frames (MRFs), the EC8 prescribes different q-factors for medium (DCM, local plastic deformations) or high (DCH, global plastic deformations) ductility classes. Macedo et al. [21] evaluated the consequences of adopting the EC8 recommended q-factor and presented a more rational selection methodology based on the specific structure and the site seismic hazard. Costanzo et al. [22] discussed existing design provision for both MRFs and chevron concentrically braced frames (CCBFs), giving evidence of a large lateral overstrength due to the codified design requirements. Also for reinforced concrete frames, studies by Kappos [23], Borzi and Elnashai [24], and Chryssanthopoulos et al. [25] assessed the reliability of q-factor values by EC8, and emphasized their high conservativity. Although the cited results from [23][24][25] looked conflicting, the joint EC8 conservatism was jointly justified with either structural overstrength or ductility supply, or both the aspects. In this context, it is thus recognized that the primary goal of standardizing committees is to simplify, on the safe side, the computational burden for designers. Such a strategy makes it possible to avoid performing complex non-linear analyses, but at the same time can severely limit the actual structural plastic capacity of the examined building systems.
For glass structures, to date, legislative and research efforts have not provided a general recommendation about realistic q-factor values that could be adopted in design. Furthermore, it is already required to satisfy global and local verifications for resistance and displacement capacities in seismic conditions [3]. As such, the typical effect often takes the form of fully elastic design (q = 1). The present study aimed to investigate further the expected structural behavior trends of seismically loaded glass members, based on the observations of a case study frame. Calibrated parameters are presented to possibly support the adaptation of consolidated general procedures to glass structures. The potentials/issues of available methodologies are assessed towards the q-factor calculation for similar structural typologies.

q-factor and Selected Calculation Methods
Different approaches can be used from literature to analytically or numerically predict the q-factor of a given structural system [26]. As far as the computational effort and accuracy of a method increase, and the reference EDPs are well defined, moreover, the q-factor estimation is progressively more robust and reliable. Figure 2 shows a typical push-over (PO) analysis result for a reinforced concrete building, in which the EDPs are qualitatively pointed out, depending on various performance levels and limit states. As usual, the most common EDPs are represented: • for structural components, by inter-story drift ratios (IDR), with inelastic component deformations and associated forces; • for non-structural components (and contents), by inter-story drift ratios (IDR) or peak floor accelerations (PFA), see [27][28][29].
materials and systems, but these parameters cannot be directly transferred to glass structures. It follows that secondary glass members that take place in a primary building must necessarily accommodate the seismic performance and capacity of the building itself (and thus satisfy the corresponding EDPs). For primary/stand-alone glass structures, otherwise, no recommended parameters are available, and thus the present study tries to provide some research developments in this direction.

Figure 2.
Expected building response and damage under seismic events. In evidence, the reference limit states and EDPs for design.
The above issue arises for novel structural systems and/or innovative materials (glass included), for which no or indicators are provided by design standards for earthquake resistant buildings. Relevant examples of literature can be found in [30][31][32][33][34][35].
From a practical point of view, the flowchart in Figure 3 can be adapted to general constructions/materials, once standardized procedures are established and a primary calculation method is chosen. In case of structural glass (as well as other innovative solutions), the critical step takes place in #2, as a direct/major effect of the analysis method choice and its input basic assumptions (first of all, the reference EDPs).  It is worth noting that recommended EDPs are available for traditional constructional materials and systems, but these parameters cannot be directly transferred to glass structures.
It follows that secondary glass members that take place in a primary building must necessarily accommodate the seismic performance and capacity of the building itself (and thus satisfy the corresponding EDPs). For primary/stand-alone glass structures, otherwise, no recommended parameters are available, and thus the present study tries to provide some research developments in this direction.
The above issue arises for novel structural systems and/or innovative materials (glass included), for which no or indicators are provided by design standards for earthquake resistant buildings. Relevant examples of literature can be found in [30][31][32][33][34][35].
From a practical point of view, the flowchart in Figure 3 can be adapted to general constructions/materials, once standardized procedures are established and a primary calculation method is chosen. In case of structural glass (as well as other innovative solutions), the critical step takes place in #2, as a direct/major effect of the analysis method choice and its input basic assumptions (first of all, the reference EDPs). It is worth noting that recommended EDPs are available for traditional constructional materials and systems, but these parameters cannot be directly transferred to glass structures.
It follows that secondary glass members that take place in a primary building must necessarily accommodate the seismic performance and capacity of the building itself (and thus satisfy the corresponding EDPs). For primary/stand-alone glass structures, otherwise, no recommended parameters are available, and thus the present study tries to provide some research developments in this direction.

Figure 2.
Expected building response and damage under seismic events. In evidence, the reference limit states and EDPs for design.
The above issue arises for novel structural systems and/or innovative materials (glass included), for which no or indicators are provided by design standards for earthquake resistant buildings. Relevant examples of literature can be found in [30][31][32][33][34][35].
From a practical point of view, the flowchart in Figure 3 can be adapted to general constructions/materials, once standardized procedures are established and a primary calculation method is chosen. In case of structural glass (as well as other innovative solutions), the critical step takes place in #2, as a direct/major effect of the analysis method choice and its input basic assumptions (first of all, the reference EDPs).   In the present paper, such an issue is further discussed with a focus on the structural glass frame described in Figure 4 and [19]. Three calculation methods (M1, M2, and M3 from Figure 3) are compared in terms of predicted q-factor, computational efficiency, accuracy, and sufficiency of results. In doing so, special care is taken for the detection of reliable EDPs that could be used for design, especially with regard to the key configurations of yielding and collapse. The so-called M1, M2, and M3 methods herein explored find inspiration from literature, but in the current study are specifically adapted to glass frames. Examples for traditional structures can be found in [36], as regards the M1 (Section 3.1) and M2 (Section 3.2) methods whilst, for the M3 one (Section 3.3), the procedure in use for the construction of fragility curves is adapted to glass. As such, the M3 q-factor is derived from linear regression on a cloud of points that is obtained from non-linear time-history analyses. In the present paper, such an issue is further discussed with a focus on the structural glass frame described in Figure 4 and [19]. Three calculation methods (M1, M2, and M3 from Figure 3) are compared in terms of predicted q-factor, computational efficiency, accuracy, and sufficiency of results. In doing so, special care is taken for the detection of reliable EDPs that could be used for design, especially with regard to the key configurations of yielding and collapse. The so-called M1, M2, and M3 methods herein explored find inspiration from literature, but in the current study are specifically adapted to glass frames. Examples for traditional structures can be found in [36], as regards the M1 (Section 3.1) and M2 (Section 3.2) methods whilst, for the M3 one (Section 3.3), the procedure in use for the construction of fragility curves is adapted to glass. As such, the M3 q-factor is derived from linear regression on a cloud of points that is obtained from non-linear time-history analyses.

Dynamic (or PGA) Method (M1)
The dynamic (or PGA) method (M1) conventionally defines the q-factor as the ratio between PGAu and PGAy, that is the peak ground acceleration values corresponding to "collapse" or "first yielding" respectively: In accordance with Equation (1), IDA were thus carried out in this paper. Based on Figure 4a and a set of input accelerograms, sequential non-linear time-history numerical analyses were performed to estimate the PGAu and PGAy values of interest. A minimum set of 7 input signals was taken into account [17].

Mixed Method (M2)
The mixed (M2) method examined in this paper still takes advantage from efficient FE simulations. The q-factor estimation was based in this case on two different contributions, that is: In Equation (2), PGAu and PGAy values agree with the definition in Section 3.1, and can be derived from the non-linear IDA for the structural system object of analysis.  [19]) for the q-factor estimation, based on (a) non-linear incremental dynamic analyses (IDA) or (b) push-over (PO) numerical procedures.

Dynamic (or PGA) Method (M1)
The dynamic (or PGA) method (M1) conventionally defines the q-factor as the ratio between PGA u and PGA y , that is the peak ground acceleration values corresponding to "collapse" or "first yielding" respectively: In accordance with Equation (1), IDA were thus carried out in this paper. Based on Figure 4a and a set of input accelerograms, sequential non-linear time-history numerical analyses were performed to estimate the PGA u and PGA y values of interest. A minimum set of 7 input signals was taken into account [17].

Mixed Method (M2)
The mixed (M2) method examined in this paper still takes advantage from efficient FE simulations. The q-factor estimation was based in this case on two different contributions, that is: In Equation (2), PGA u and PGA y values agree with the definition in Section 3.1, and can be derived from the non-linear IDA for the structural system object of analysis. At the same time, V y and V d in Equation (2) denote the base shear load corresponding to "first significant yield strength" and "allowable design strength". These base shear values were conventionally derived from non-linear PO simulations according to Figure 4b.

Cloud Analysis (M3) with Linear Regression
The q-factor of the examined frame was finally estimated in this paper by using the inelastic response spectrum, with the support of the so-called cloud analysis and the spectral acceleration (S a ) definitions [36]. Successful cloud analysis applications can be found in [37][38][39] for various structural typologies and materials.
Differing from Sections 3.1 and 3.2 (IDA procedure), the cloud analysis is carried out with the support of a set of unscaled accelerograms. The set of input signals (60 in the present study, Table A1 [40]) must be established to ensure an appropriate distribution of cloud data [41,42]. It is in fact known that the major issue of IDA method may consist of a significant computational cost, and most often a marked scaling of original records to various intensity levels, before the desired EDPs could be achieved. This effort is not required in cloud analysis. Moreover, the use of unscaled natural accelerograms in cloud analysis allows to keep all the information related to the event, also known as "record-to-record variability". From unscaled signals and non-linear dynamic analyses, the correlation is established between selected EDPs and some intensity measure (IM) values of the imposed signals by taking advantage of linear regression [37]. As in case of IDA, however, the unscaled signals require an accurate definition of reference EDPs and are also expected to cover a useful range of values for identifying the required limit states. These signals are thus sensitive to the fundamental vibration period T 1 (to predict) and the characteristics of the structure to verify (material properties, damage mechanisms, etc.). Moreover, the input signals should be selected to be representative of the seismic hazard of the site under investigation. When appropriate signals are not available, site-specific ground motion modeling techniques can also be used [43][44][45]. Based on EC8, the q-factor can be finally calculated as: where S a (T 1 ) is the spectral ordinate corresponding to the characteristic period of the design spectrum. The subscripts "u" and "y" in Equation (3) refer to the "collapse" and "first yielding" configurations.

Geometrical and Mechanical Properties
The current research study follows and extends the analytical and numerical investigation reported in [19]. As such, some geometrical and mechanical features are summarized herein for the system in Figure 5a. Each glass frames followed the layout in Figure 5b, with H = 6 m and L = 8 m. Both the beam and column sections were composed of heatstrengthened (HS) LG members, with uniform size (h = 600 mm high × t tot = 66 mm thick) given by 5 × t g = 12 mm glass layers and t int = 1.52 mm thick ionoplast foils. The mechanical connection at each beam-column interception took the form of an ideal pin, see Figure 5b. Possible out-of-plane deformations of the frame were restrained, and the related mechanisms (including lateral-torsional buckling for beams [46,47], or coupled bending-compressive buckling for columns [48]) can be preliminarily disregarded. For the base restraints of columns, stainless steel pins pass through two holes in the glass (ϕ g = 32 mm in diameter, with ϕ = 24 mm the nominal diameter of bolts and D = 500 mm their distance).

Preliminary Elastic Seismic Design of the Frame
For calculation purposes, the glazed assembly of Figure 5 is located in a high seismicity region of Italy. Based on [17,49,50], its strength and stiffness should be verified to resist the most unfavorable expected seismic combination of actions Ed, that is: where Rd is the structural capacity. A simple design of glass members should properly verify that LG columns and beams are not subjected-due to the imposed in-plane seismic loads-to relevant stress peaks and premature fracture. As for regular structures in plan and elevation, the input seismic force Fi the frame should resist is given by:

Preliminary Elastic Seismic Design of the Frame
For calculation purposes, the glazed assembly of Figure 5 is located in a high seismicity region of Italy. Based on [17,49,50], its strength and stiffness should be verified to resist the most unfavorable expected seismic combination of actions E d , that is: where R d is the structural capacity.
A simple design of glass members should properly verify that LG columns and beams are not subjected-due to the imposed in-plane seismic loads-to relevant stress peaks and Sustainability 2021, 13, 9291 8 of 25 premature fracture. As for regular structures in plan and elevation, the input seismic force F i the frame should resist is given by: with i = 1, . . . n, m i , m j the story masses and z i , z j their height from the foundation, while the base shear F b is: and: • W the aboveground total mass of the building object of analysis, • S d (T 1 ) the design acceleration from the reference spectrum, as a function of the vibration period T 1 , with S d = 0.35 g the peak ground acceleration (high seismic region of Italy), • λ = 1 a correction factor for one-story buildings with T 1 > 2T C (otherwise 0.85), and • q ≥ 1 the behavior factor of the system.
Given the lack of more appropriate recommendations, the conventional design of the case study columns suggests the assumption that q = 1. Disregarding the vertical loads that the LG members must sustain (as a part of the framed system of Figure 5), the in-plane lateral force affects the region of glass holes at the base connections, that is: with: the magnification factor for stresses [51][52][53], while the tensile stress σ t in glass is given by: with: The resistance verification of the LG columns in seismic conditions requires that: with: the design resistance [49]. Among the coefficients in Equation (12), the short-term duration of seismic events (conventionally set in 30 s [3]) suggests k mod = 0.78. Given that the columns are composed of HS glass, Equation (12) results in f g;d ≈ 75 MPa, that is ≈ 1/3rd the maximum stress from Equation (7), due to the seismic shear from Equation (6).
To avoid the improper sizing of load-bearing glass members, the design (with given input parameters) would require the exploitation of a minimum q min ≈ 3. In other words, the ratio of Equation (7) to Equation (12) and combination with Equation (6) can be used for simple analytical estimates of (minimum) required plastic capacities of the frame, towards the seismic demand. In this regard, it is also worth noting that the analytical model developed in [19] for the ductility estimation of base angle brackets (and properly combined with the stress analysis in Equations (7)-(12)), would result in q = 4.58 (collapse governed by stress peaks in the region of glass holes). However, such an analytical prediction is not able to account for complex mechanisms in the frame as a whole, under dynamic seismic accelerations.

Numerical Model
The reference numerical analysis was carried out in ABAQUS/Standard [54] on a full three-dimensional model representative of the glass frame object of analysis, inclusive of LG members, and reproducing the geometrical details for base connections ( Figure 6). For symmetry, 1/4th the geometry was taken into account.

Numerical Model
The reference numerical analysis was carried out in ABAQUS/Standard [54] on a full three-dimensional model representative of the glass frame object of analysis, inclusive of LG members, and reproducing the geometrical details for base connections ( Figure 6). For symmetry, 1/4th the geometry was taken into account.
Differing from [19], the seismic response of the frame as a whole was explored for the purpose of this study. To this aim, a novel optimized FE model was developed to maximize its computational efficiency.
Solid brick elements (C3D8R type from ABAQUS library) were used for rigid base support, angular members, frictionless foils, and bolts. For the glass plate, a mix of brick solid elements and shell elements was used to preserve the accuracy of stress distributions in the regions of holes. Finally, the column in elevation (and top beam) were described as pinned rigid beams. The overall symmetry assumption resulted in 18,000 solid/shell elements and 75,000 DOFs for the frame model in Figure 6. After preliminary validation, such a solution was used to replace the FE assembly from [19], in which a total of 45,000 solid elements and 170,000 DOFs were used for half geometry of the base connection only.  Differing from [19], the seismic response of the frame as a whole was explored for the purpose of this study. To this aim, a novel optimized FE model was developed to maximize its computational efficiency.
Solid brick elements (C3D8R type from ABAQUS library) were used for rigid base support, angular members, frictionless foils, and bolts. For the glass plate, a mix of brick solid elements and shell elements was used to preserve the accuracy of stress distributions in the regions of holes. Finally, the column in elevation (and top beam) were described as pinned rigid beams. The overall symmetry assumption resulted in 18,000 solid/shell elements and 75,000 DOFs for the frame model in Figure 6. After preliminary validation, such a solution was used to replace the FE assembly from [19], in which a total of 45,000 solid elements and 170,000 DOFs were used for half geometry of the base connection only.

Materials and Contact Interactions
Key mechanical assumptions for glass and steel members were derived from [19]. An elastic-perfectly plastic law was used for mild steel, with E s = 210 GPa the modulus of elasticity, ν s = 0.3 the Poisson' ratio, and σ s,y = σ s,u = 235 MPa the yielding/failure strength, with corresponding strain values equal to ε y = 0.112% and ε u = 25%. The ductile damage material option was also accounted for in FE analyses to detect the possible initiation of ductile failure mechanism in angle brackets.
An elasto-plastic law was also used for steel bolts, with σ b,y = σ b,u = 1000 MPa the yielding/ultimate resistance (8.8 resistance class). Finally, the rubber layers were described in the form of an equivalent elastic-perfectly plastic material, with E r = 30 GPa and ν r = 0.3. The yielding/ultimate stress was conventionally set at 2.4 MPa [19].
The tensile brittleness of glass was included with the concrete damaged plasticity (CDP) model from the ABAQUS library [19]. While the CDP model was primarily developed for concrete, literature studies show that the same model can be efficiently used for structural glass members (under specific loading/boundary conditions), as in the present application. From the post-processing of FE results, fracture initiation in glass was in fact assumed as a reference for "failure". This means that the overall post-cracked stage was disregarded but the simulation was prevented from additional uncertainties that are typical of the post-cracked response of glass under cyclic loads. In doing so, the nominal mechanical properties for HS were taken into account (E g = 70 GPa, ν g = 0.23 and σ tk =70 MPa). Further, the characteristic compressive strength was set to σ ck = 300 MPa (350-500 MPa [19] the reference strength).
A set of surface-to-surface contacts at the interface of adjacent FE components allowed to reproduce the in-plane lateral response of the frame under seismic loads (with "penalty" tangential characteristic (friction µ = 0.3) and "hard" normal features). "Tie" mechanical constraints were also used to rigidly connect some FE components (i.e., the head of each bolt and the corresponding angle bracket, or the frictionless layer and the adjacent angle bracket).

Loading Strategy
The frame was investigated by taking into account the presence of in-plane seismic loads and dead loads due to constructional members, plus a vertical accidental load Q k = 3 kN/m 2 (with i = 1 m). The FE assembly of Figure 6 was used for both the required non-linear static PO and time-history dynamic analyses. As such, two different solving procedures were taken into account, based on two separate steps representative of: • S1 = an initial stage for introduction of dead and accidental loads (5 s), followed by • S2 = seismic analysis of the pre-loaded glass frame (60 s).
In case of IDA, the main seismic input consisted in the selected accelerogram in Figure 7 (acceleration-time history at the base of the frame for M1 and M2 methods). The used earthquake records were derived from [55], that is considering a PGA of 0.35 g, with type A soil (rock soil), topographic category T1, and a reference nominal life of 50 years. A maximum lower and upper tolerance of 10% was also considered.
For the PO analyses (M2 approach), otherwise, the FE system of Figure 6 was subjected to a linear increasing in-plane shear force in accordance to Figure 7 (base connection rigidly fixed to ground). Finally, in case of cloud analysis (M3), a total of 60 unscaled signals (Table A1 [40]) was taken into account to replace the 7 input signals of M1 and M2 procedures. Sustainability 2021, 13, x FOR PEER REVIEW 11 of 25 Figure 7. Reference set of time-acceleration histories, as derived from REXEL [55] for the non-linear incremental dynamic analysis (IDA) of the case study glass frame.

Detection of First Yielding, Collapse, and Allowable Strength Parameters
Given the general definitions of calculation approaches in Sections 4.1 to 4.3 and the structural system object of study, special care must be taken for the definition of the key design configurations. Major design challenges derive from the lack of explicit recommendations for glass structures in current design standards for seismic resistant buildings. On the other hand, the in-plane seismic response of the frame is strongly affected by the intrinsic mechanical properties of its basic components, namely: • brittle elastic glass panels (with holes), and • flexible angle brackets at the base of the frame.
In other words, the "first yielding" condition of the system of Figure 6 was defined in this project as the first plastic deformation of angle brackets in tension (with σs,y = 235 MPa the reference strength and δy = 0.677 mm the corresponding vertical deformation [19]). Regarding the "collapse" damage state for the frame, maximum drift amplitudes (or column rotations, or even vertical deformations of the steel angle brackets) should be checked. Globally, for the parametric investigation herein summarized, the control of local and global critical conditions for the frame was primarily based on local stress and displacement controls, namely representing a potential: In addition, for comparative purposes, conventional deformation limits available in design standards were also taken into account. For the case study frame, as far as the glass holes are properly protected from potential local damage, the seismic analysis could take advantage of the intrinsic flexibility and dissipative capacity of angle brackets. In this Figure 7. Reference set of time-acceleration histories, as derived from REXEL [55] for the non-linear incremental dynamic analysis (IDA) of the case study glass frame.

Detection of First Yielding, Collapse, and Allowable Strength Parameters
Given the general definitions of calculation approaches in Sections 4.1 and 4.2 and the structural system object of study, special care must be taken for the definition of the key design configurations. Major design challenges derive from the lack of explicit recommendations for glass structures in current design standards for seismic resistant buildings. On the other hand, the in-plane seismic response of the frame is strongly affected by the intrinsic mechanical properties of its basic components, namely: • brittle elastic glass panels (with holes), and • flexible angle brackets at the base of the frame.
In other words, the "first yielding" condition of the system of Figure 6 was defined in this project as the first plastic deformation of angle brackets in tension (with σ s,y = 235 MPa the reference strength and δ y = 0.677 mm the corresponding vertical deformation [19]). Regarding the "collapse" damage state for the frame, maximum drift amplitudes (or column rotations, or even vertical deformations of the steel angle brackets) should be checked. Globally, for the parametric investigation herein summarized, the control of local and global critical conditions for the frame was primarily based on local stress and displacement controls, namely representing a potential: In addition, for comparative purposes, conventional deformation limits available in design standards were also taken into account. For the case study frame, as far as the glass holes are properly protected from potential local damage, the seismic analysis could take advantage of the intrinsic flexibility and dissipative capacity of angle brackets. In this sense, the reliability of limit values of Table 1 [58], EC8 [17], and NTC2018 [50] documents and their applicability to the examined frame were taken into account in this study. Finally, the "allowable strength" condition required by the M2 approach should also be defined. As far as the brackets are assumed responsible of the overall in-plane seismic performance of the frame, the yielding stress of steel suggests that: with γ M the partial safety factor, thus a minimum: to account in Equation (2). Figure 8 shows the base shear-lateral deformation of the frame from PO analysis. It is clear that the above assumption can strongly penalize the frame response, and the critical glass members prove to offer a safety factor in the order of ≈2.1 against potential tensile cracks. Moreover, it is possible to see that the compressive limit in glass holes is not achieved, neither under large displacements. This results from the high deformation capacity of the frame, thanks to detailing of base connections explored in [19].  [17], and NTC2018 [50] documents and their applicability to the examined frame were taken into account in this study. Finally, the "allowable strength" condition required by the M2 approach should also be defined. As far as the brackets are assumed responsible of the overall in-plane seismic performance of the frame, the yielding stress of steel suggests that: with γM the partial safety factor, thus a minimum: to account in Equation (2). Figure 8 shows the base shear-lateral deformation of the frame from PO analysis. It is clear that the above assumption can strongly penalize the frame response, and the critical glass members prove to offer a safety factor in the order of ≈2.1 against potential tensile cracks. Moreover, it is possible to see that the compressive limit in glass holes is not achieved, neither under large displacements. This results from the high deformation capacity of the frame, thanks to detailing of base connections explored in [19].

Seismic Performance Assessment
The seismic response of the frame was found to agree with Figure 9, where the typical IDA deformed shape (detail) is proposed for brackets under large in-plane lateral displacements.

Seismic Performance Assessment
The seismic response of the frame was found to agree with Figure 9, where the typical IDA deformed shape (detail) is proposed for brackets under large in-plane lateral displacements.
A total of 140 non-linear analyses was carried out with the imposed scaled accelerograms from Figure 7 (with an average of ≈20 differently scaled simulations for each accelerogram). The IDA results still confirmed the close correlation with PO results in Figure 8, with a qualitative agreement of damage phenomena and maximum effects due to the imposed design loads. Figure 9b A more detailed analysis of IDA results can be found in Figures 10 and 11, in terms of relevant EDPs, as a function of the imposed PGA for each one of the input scaled signals. It is worth noting that the parametric analysis was carried out in the ideal PGA range of 0-50 m/s 2 to address the performance of structural components. In this regard, typical PGA values can be seen as associated to limited stress levels in the structure, as is expected due to the limited structural mass and high flexibility of the system. Key benefits derive also from gaps in the region of glass holes to prevent premature stress peaks at the edges. Such an approach is also in line with other studies on the seismic performance of structural systems with flexible joints (see for example [59]). A total of 140 non-linear analyses was carried out with the imposed scaled accelerograms from Figure 7 (with an average of ≈20 differently scaled simulations for each accelerogram). The IDA results still confirmed the close correlation with PO results in Figure 8, with a qualitative agreement of damage phenomena and maximum effects due to the imposed design loads. Figure 9b, in this regard, presents the evolution of maximum stress peaks in the region of holes, while Figure 9c,d focus on the bracket and frame responses, respectively.
A more detailed analysis of IDA results can be found in Figures 10 and 11, in terms of relevant EDPs, as a function of the imposed PGA for each one of the input scaled signals. It is worth noting that the parametric analysis was carried out in the ideal PGA range of 0-50 m/s 2 to address the performance of structural components. In this regard, typical PGA values can be seen as associated to limited stress levels in the structure, as is expected due to the limited structural mass and high flexibility of the system. Key benefits derive also from gaps in the region of glass holes to prevent premature stress peaks at the edges. Such an approach is also in line with other studies on the seismic performance of structural systems with flexible joints (see for example [59]). and vertical loads. Otherwise, it is also interesting to notice that the LG members can sustain relatively strong earthquake motions, before glass could fracture. A relevant aspect is hence represented, in both figures, by the non-linear evolution of stress peaks with the imposed PGA. A first linear trend of the charts can be observed for PGA up to ≈6 m/s 2 , and such a slope change coincides with first yielding (and progressive plastic deformation) of angle brackets. This limit condition was generally achieved for PGA in the order of ≈4 m/s 2 (Figure 10c,d). Compared to the stress evolution in the holes region, similar trends can also be observed for the deformations of the frame in Figure 11.
The vertical displacement δ, in-plane lateral drift u / H and base rotation θ are proposed, as obtained from IDA and maximum envelopes of selected EDPs. Under the input assumptions of this study, the collapse condition is never achieved on the side of angle brackets (Figure 11a). The limit drift of 2% or 5% is exceeded for PGA in the order of ≈10 m/s 2 and 20 m/s 2 (average value), see Figure 11b. The 2% drift, finally, is mostly in line with the 0.03 rad rotation of the frame (Figure 11c). As shown in Figure 10a,b, the compressive stress peaks were mostly observed to double the corresponding tensile peaks in glass, due to a combination of in-plane lateral and vertical loads. Otherwise, it is also interesting to notice that the LG members can sustain relatively strong earthquake motions, before glass could fracture. A relevant aspect is hence represented, in both figures, by the non-linear evolution of stress peaks with the imposed PGA. A first linear trend of the charts can be observed for PGA up to ≈6 m/s 2 , and such a slope change coincides with first yielding (and progressive plastic deformation) of angle brackets. This limit condition was generally achieved for PGA in the order of ≈4 m/s 2 (Figure 10c,d).
Compared to the stress evolution in the holes region, similar trends can also be observed for the deformations of the frame in Figure 11.
The vertical displacement δ, in-plane lateral drift u / H and base rotation θ are proposed, as obtained from IDA and maximum envelopes of selected EDPs. Under the input assumptions of this study, the collapse condition is never achieved on the side of angle brackets (Figure 11a). The limit drift of 2% or 5% is exceeded for PGA in the order of ≈10 m/s 2 and 20 m/s 2 (average value), see Figure 11b. The 2% drift, finally, is mostly in line with the 0.03 rad rotation of the frame (Figure 11c).

q-factor Estimates
The analysis was first focused on the so-called M1 method. For the M2 case, the M1 value was adapted with the magnification factor in Equation (14). The so-calculated IDA results are proposed in Figure 12.
Note that the attention was focused on the most unfavorable collapse mechanism for the frame as a whole. This was generally observed to coincide with tensile glass cracking ("glass"), while in two cases, only the deformability of the base joint allowed to reach a lateral drift of 2% ("Drift 2%"). The corresponding q-factor values are presented for the M1 method (Equation (1)), in the range from 1 (a#1) to 5 (a#6). The average value of q = 2.59 is compared with the corresponding M2 estimate from Equation (2), q = 2.72. The preliminary analytical requirement (Equations (7)-(12) combined with Equation (6)) is also highlighted (qmin = 3), while the analytical value based on local analysis (q = 4.58 from [19]) gives evidence of intrinsic limits due to simple predictions carried out for the base connection only (with collapse governed by stress peaks in the region of glass holes).

q-factor Estimates
The analysis was first focused on the so-called M1 method. For the M2 case, the M1 value was adapted with the magnification factor in Equation (14). The so-calculated IDA results are proposed in Figure 12.
Note that the attention was focused on the most unfavorable collapse mechanism for the frame as a whole. This was generally observed to coincide with tensile glass cracking ("glass"), while in two cases, only the deformability of the base joint allowed to reach a lateral drift of 2% ("Drift 2%"). The corresponding q-factor values are presented for the M1 method (Equation (1)), in the range from 1 (a#1) to 5 (a#6). The average value of q = 2.59 is compared with the corresponding M2 estimate from Equation (2), q = 2.72. The preliminary analytical requirement (Equations (7)-(12) combined with Equation (6)) is also highlighted (q min = 3), while the analytical value based on local analysis (q = 4.58 from [19]) gives evidence of intrinsic limits due to simple predictions carried out for the base connection only (with collapse governed by stress peaks in the region of glass holes). Sustainability 2021, 13, x FOR PEER REVIEW 16 of 25

Input Records and EDPs
The maximum inter-story drift (IDR) was chosen as reference EDP for the frame (T1 = 0.3 s), with PGA and pseudo-spectral acceleration Sa(T1) being selected as IM parameters. While the PGA value is only related to the seismic ground motion, the Sa value depends on the dynamic behavior of the structure; therefore, the performance to different IMs depends on the type of structure and the governing failure mechanism.
The preliminary analysis was focused on the distribution of tensile stress peaks in a glass column (with nominal height H) deprived of the angle brackets (rigid base connection). Given the characteristic tensile resistance of HS glass (σtk = 70 MPa), an inter-story displacement uu= 0.048 m was calculated as in Figure 13. The reference EDP corresponds to contour plots in Figure 13a,b), while Figure 13c shows the data trend obtained at the column base and in terms of maximum envelope.
It is worthy of interest that the so-calculated value corresponds to u/H ≈ 0.007 and is in close correlation with consolidated limit values for constructional materials characterized by typical brittle behavior in tension, such as, for example, masonry [50]. At the same time, the calculated value significantly minimizes the expected seismic capacity of the frame, thus confirming the key role of its base connections.
The total set of 60 unscaled ground motion records in Table A1 were chosen from [40], depending on the possible collapse mechanism of the frame. According to procedures for general buildings, special attention was paid to cover a wide range of spectral accelerations, but also to respect the consistency between the characteristics of selected records and the supposed classification for the site of interest. As a result, the selected accelerograms were characterized by a moment magnitude (Mw) between 5.6 and 7.6, an epicentral distances (R) ranging between 3.5 km and 62.9 km, and a soil class type A or B (EC8 classification).

Input Records and EDPs
The maximum inter-story drift (IDR) was chosen as reference EDP for the frame (T 1 = 0.3 s), with PGA and pseudo-spectral acceleration S a (T 1 ) being selected as IM parameters. While the PGA value is only related to the seismic ground motion, the S a value depends on the dynamic behavior of the structure; therefore, the performance to different IMs depends on the type of structure and the governing failure mechanism.
The preliminary analysis was focused on the distribution of tensile stress peaks in a glass column (with nominal height H) deprived of the angle brackets (rigid base connection). Given the characteristic tensile resistance of HS glass (σ tk = 70 MPa), an inter-story displacement u u = 0.048 m was calculated as in Figure 13. The reference EDP corresponds to contour plots in Figure 13a,b), while Figure 13c shows the data trend obtained at the column base and in terms of maximum envelope.
It is worthy of interest that the so-calculated value corresponds to u/H ≈ 0.007 and is in close correlation with consolidated limit values for constructional materials characterized by typical brittle behavior in tension, such as, for example, masonry [50]. At the same time, the calculated value significantly minimizes the expected seismic capacity of the frame, thus confirming the key role of its base connections.
The total set of 60 unscaled ground motion records in Table A1 were chosen from [40], depending on the possible collapse mechanism of the frame. According to procedures for general buildings, special attention was paid to cover a wide range of spectral accelerations, but also to respect the consistency between the characteristics of selected records and the supposed classification for the site of interest. As a result, the selected accelerograms were characterized by a moment magnitude (M w ) between 5.6 and 7.6, an epicentral distances (R) ranging between 3.5 km and 62.9 km, and a soil class type A or B (EC8 classification).

Analysis of M3 Results with Linear Regression
The results of the cloud analysis method related to PGA and Sa(T1) are separately collected in Figure 14, with attention to the measured inter-story displacement.
Differing from IDA, one of the potential intrinsic limits of the M3 method can manifest in the availability of natural seismic records that possess sufficiently high accelerations to reach the desired EDPs.
In this regard, Figure 15 gives evidence of the typical observed response for the case study frame. As shown, the imposed records are able to lead the angle brackets to yielding (Figure 15a), but still relatively smooth stress peaks are achieved in glass (holes), with tensile and compressive stress peaks in Figure 15b,c. In the same way, the measured in-plane deformations of the frame are still lower than the % limit values earlier discussed for IDA.

Analysis of M3 Results with Linear Regression
The results of the cloud analysis method related to PGA and S a (T 1 ) are separately collected in Figure 14, with attention to the measured inter-story displacement.
Differing from IDA, one of the potential intrinsic limits of the M3 method can manifest in the availability of natural seismic records that possess sufficiently high accelerations to reach the desired EDPs.
In this regard, Figure 15 gives evidence of the typical observed response for the case study frame. As shown, the imposed records are able to lead the angle brackets to yielding (Figure 15a), but still relatively smooth stress peaks are achieved in glass (holes), with tensile and compressive stress peaks in Figure 15b,c. In the same way, the measured in-plane deformations of the frame are still lower than the % limit values earlier discussed for IDA.   By considering the entire set of available data from the cloud analysis of the frame, an ordinary least-square linear regression was thus performed. The analysis was carried out in the logarithmic space, given that the use of logarithm of variables improves the fit of the model by transforming the distribution of the features to a more normally shaped bell curve. In order to control the skew and counter problems in heteroskedasticity, both the dependent variable (IM) and the independent variable (EDP) were log-transformed. The final result is proposed in Figure 16, where the linear fits of cloud data are obtained from the least squares method. By considering the entire set of available data from the cloud analysis of the frame, an ordinary least-square linear regression was thus performed. The analysis was carried out in the logarithmic space, given that the use of logarithm of variables improves the fit of the model by transforming the distribution of the features to a more normally shaped bell curve. In order to control the skew and counter problems in heteroskedasticity, both the dependent variable (IM) and the independent variable (EDP) were log-transformed. The final result is proposed in Figure 16, where the linear fits of cloud data are obtained from the least squares method. The thresholds of "yielding" (EDPY,50 = 0.037 m) and "collapse" (EDPCP,50) performance levels are indicated in Figure 16 by vertical lines. Two different thresholds (noted as "A" and "B") were used to identify the collapse prevention limit and to quantify further the influence of the base steel connection in the seismic response of the frame, namely: In this regard, it should be noted that the regression line was assumed to be valid for "B" collapse value of displacement even if it is outside the available data cloud. Once the regression line is found, the IM characterizing yielding and collapse prevention were obtained in Figure 16 using the following relations: The q-factor estimation can thus be based on Figure 16, for PGA and Sa(T1), respectively. Certainly, the obtained results are affected by base steel joints and thus by EDPCP,50. As such, the above outcome should be taken into account as a general approach for basic design considerations of similar structures, given that the resistance and stiffness of joints are strictly responsible for the final ductility of the frame, and thus for the possible fracture initiation in glass. The thresholds of "yielding" (EDP Y,50 = 0.037 m) and "collapse" (EDP CP,50 ) performance levels are indicated in Figure 16 by vertical lines. Two different thresholds (noted as "A" and "B") were used to identify the collapse prevention limit and to quantify further the influence of the base steel connection in the seismic response of the frame, namely: In this regard, it should be noted that the regression line was assumed to be valid for "B" collapse value of displacement even if it is outside the available data cloud. Once the regression line is found, the IM characterizing yielding and collapse prevention were obtained in Figure 16 using the following relations: The q-factor estimation can thus be based on Figure 16, for PGA and S a (T 1 ), respectively. Certainly, the obtained results are affected by base steel joints and thus by EDP CP,50 . As such, the above outcome should be taken into account as a general approach for basic design considerations of similar structures, given that the resistance and stiffness of joints are strictly responsible for the final ductility of the frame, and thus for the possible fracture initiation in glass.

Comparative q-factor Predictions
In conclusion, Figure 17 shows the M3 calculated q-factor and a comparison of selected methods (average). As expected, q significantly decreases as far as the M3 approach at collapse disregards the beneficial effect of brackets in the post-yielded stage (collapse "A"). At the same time, as far as the real ultimate inter-story displacement is considered for the tensile fracture of the column (collapse "B"), Figure 17 proves a stable q-factor estimation from M3 or M1-M2 methods. Equations (7)-(12) combined with Equation (6)), so as to preserve the glass columns from fracture. At the same time, it is necessary to highlight the relatively stable trend for the qfactor numerical estimates from the M1 to M3 selected approaches. Most importantly, this finding seems to confirm the potential of linear regression method based on cloud analysis, thanks to the computational efficiency of the M3 method. The M1 and M2 procedures, while limited in number of signals, are univocal in EDPs detection but could require major calculation efforts compared to M3 (140 scaled simulations, in the present study). Furthermore, the IDA calculated average q-factor can be highly sensitive to input signals (7 minimum). While the present investigation suggests a very good correlation of M1 to M3 average q-factor predictions, this could not be the case of different structural members, thus requiring even more pronounced calculation efforts from IDA (M1 or M2). Finally, compared to simple analytical estimates that are not able to account for complex mechanical phenomena of the frame as a whole (i.e., q = 4.58 from [19]), all the numerical estimates in Figure 17 are on the conservative side, thus confirming the need for refined models and non-linear dynamic procedures in support of seismic design. Figure 17. Calculated q-factor for the examined frame, based on M1 to M3 methods and EDPs. Figure 17. Calculated q-factor for the examined frame, based on M1 to M3 methods and EDPs.
The positive outcome is confirmation of intrinsic ductility and post-yielding capacity for the frame as a whole. This is in line with preliminary results from [19]. Moreover, the local analysis of angle bracket ductility and stress peak estimation in the region of glass holes was quantified in [19] up to q = 4.58 for the frame (collapse governed by tensile fracture of glass). The present study, consequently, confirms the need for full-size structural analyses for special structures and joint details.
Comparative data in Figure 17 are also a confirmation of simple analytical expectations about the minimum plastic capacity of the frame from Section 4.2 (with q min ≈ 3 from Equations (7)-(12) combined with Equation (6)), so as to preserve the glass columns from fracture. At the same time, it is necessary to highlight the relatively stable trend for the q-factor numerical estimates from the M1 to M3 selected approaches. Most importantly, this finding seems to confirm the potential of linear regression method based on cloud analysis, thanks to the computational efficiency of the M3 method. The M1 and M2 procedures, while limited in number of signals, are univocal in EDPs detection but could require major calculation efforts compared to M3 (140 scaled simulations, in the present study). Furthermore, the IDA calculated average q-factor can be highly sensitive to input signals (7 minimum). While the present investigation suggests a very good correlation of M1 to M3 average q-factor predictions, this could not be the case of different structural members, thus requiring even more pronounced calculation efforts from IDA (M1 or M2). Finally, compared to simple analytical estimates that are not able to account for complex mechanical phenomena of the frame as a whole (i.e., q = 4.58 from [19]), all the numerical estimates in Figure 17 are on the conservative side, thus confirming the need for refined models and non-linear dynamic procedures in support of seismic design.

Conclusions
Available design standards for seismic-resistant buildings provide various recommendations in support of analysis and safe design of several structures subjected to earthquakes, but no specific details are given for glass systems. Among others, major uncertainties derive from the reliable calculation of the seismic performance and dissipation capacity of glass structures, thus their q-factor.
In this paper, attention was focused on the local/global seismic analysis of a structural glass frames under in-plane lateral seismic loads. Careful consideration was paid for the development of efficient finite element (FE) numerical models in support of extended parametric non-linear dynamic analyses that could be used to adapt/assess for glass some consolidated procedures in use for structural systems composed of ordinary materials. For most traditional materials and systems, reference engineering demand parameters (EDPs) are recommended by standards or literature documents. On the other hand, reliable EDPs are still lacking for the methods' adaptation to glass structures.
Three numerical calculation methods were taken into account for q-factor estimates, based on the parametric incremental dynamic analysis (IDA; dynamic "M1" and mixed "M2" methods), and the cloud analysis based on linear regression ("M3"). Numerical calculations were also compared to simple analytical estimates.
From the FE parametric outcomes, more in detail, it was proven that the metal joints in use for structural glass applications were the major source of possible critical failure mechanisms, but also a key source of enhanced ductility performances for glass members. Such a finding was confirmed in line with ductility and flexibility capacities discussed in [19], based on local analysis of the base connection of the frame. In addition, the present study also confirmed the need of full-size FE models and non-linear dynamic procedures.
In terms of calculated q-factor values, more in detail, it was shown that: • IDA-based approaches (M1 or M2) are univocal in damage detection, thus in the corresponding estimation of reliable EDPs; • Both M1 and M2 procedures are indeed strongly expensive in computational cost. The present study, for example, was based on a minimum of 7 accelerograms and required up to 140 non-linear dynamic analyses; and • High sensitivity was observed for the predicted average q-values from M1 or M2, thus recommending a careful selection of input signals, but also the possible use of largest sets of scaled records.
At the same time, the adaptation of M3 method with linear regression to structural glass frames:

•
Confirmed the reduced computational cost of the approach, compared to M1 or M2 methods (60 unscaled signals and analyses in total for the present study, compared to 140 simulations); and • Confirmed that reliable EDPs for special structures should be properly calculated, with the support of refined numerical models or even experimental tests. Existing consolidated EDPs of literature and standards for seismic-resistant structures can hardly adapt to special glass systems and members. • However, the M3 procedure also gave evidence of some difficulties of dataset interpretation (due to limited stress/deformation levels in the load-bearing members, for some simulations). The reason was found in the set of unscaled input accelerograms that sometimes (when applied to structures characterized by limited self-weight and high flexibility as in the present study) can hardly achieve the desired EDPs at collapse; and • Furthermore, the FE parametric study proved that-once EDPs and damage mechanisms for relevant limit states are established-the M3 approach can offer rather accurate predictions for glass structures under seismic loads, and thus support as an efficient tool the estimation of q-factor for the seismic design of similar structural systems. For the present case-study frame, the calculated q-factor was in fact in line, but on the conservative side compared to simple analytical predictions from [19], based on the local analysis of bracket ductility and stress peaks in the region of glass holes. Such a finding also confirms the need for complex numerical models able to capture dynamic mechanical phenomena in similar systems, as a more detailed investigation to combine with simplified analytical procedures. Data Availability Statement: Data will be made available upon request.

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A