A Seismic Risk Assessment of Concrete-Filled Double-Skin Steel Tube (CFDST) Frames with a Beam-Only-Connection for Reinforced Concrete Shear Walls (BRWs)

: The beam-only connected reinforced concrete shear wall (BRW) is used as a reinforcing component to enhance the seismic performance of concrete-ﬁlled, double-skin steel tube (CFDST) frames. The effects of the BRW on seismic risks of CFDST frames are investigated. Three performance levels of limit states are deﬁned and described according to the typical failure of test specimens and the existing deﬁnition of current guidance. A simpliﬁed numerical model is calibrated for CFDST frame-BRW structures, and test results validate it. Nonlinear dynamic analyses on a nine-story CFDST-BRW building are performed to investigate the effects of BRW on reducing the seismic risk of CFDST buildings. The results show that the BRW reduces the probability of collapse of the CFDST frame to 2.76% in 50 years, which can effectively reduce the probability of different degrees of damage in the service cycle of the structure. The results provide information for risk-informed decision-making on the design of CFDST frame-BRW structures.


Introduction
The traditional concrete-filled steel tubular (CFST) frame structures rely on the connection of CFST frame joints (finger beam-column joint) to achieve lateral resistance.Once the CFST structure joints are damaged, they will collapse.Therefore, the frame-shear wall structure system and seismic design method have gradually attracted attention.Many researchers have conducted experimental, theoretical and numerical studies on the design methods of different structures.
So far, many experimental and theoretical studies have been carried out on the performance enhancement of CFST columns and frames.Wang [1] designed thin steel circular CFSTs.The specimens were tested, and their seismic performances were evaluated.Skalomenos [2] proposed using the HFD method to guide the seismic assessment of CFST frames.Gan [3] studied the mechanical properties of shear connectors in CFST frame structures.Lai [4] and Wang [5] used high-strength materials to improve the seismic capacities of CFST frames.Bai [6] established a model of a CFST frame and verified the correctness of the connection parts of the numerical model through experiments.Ren [7] and Jia [8] studied the typical working relationships between energy dissipation devices, buckling restrained braces, and CFST frames.
However, because the frame structure realizes the lateral resistance through the connections of beam and column members, only one seismic defense is included, and the scope of the application is limited.People gradually advocate using shear walls to improve the performance of frame structures.
For this reason, many shear wall types, frame connection methods, and seismic theories have been proposed.Nie [9] suggested using a composite shear wall to strengthen a CFST frame.The cooperative working relationship between the composite shear wall and the CFST frame was deduced.Wang [10] conducted seismic performance tests on concretefilled steel tubular infilled walls.Guo [11] proposed a connection method between the CFST frames and partially infilled walls and conducted experimental research on their seismic performance.Hu [12] studied the cyclic performance of a CFDST frame strengthened by a BSW.Bian [13] proposed an assembled structure composed of a new frame-shear wall structure.Hu [14] proposed a fully precast CFSDT frame, which is only connected to the BSW system.An SPSW with a beam-only connection for this system was presented.Xu [15] proposed a precast CFST frame support system by using reinforced concrete beams and a shaking table test.Skalomenos [16] studied different modeling methods concerning the seismic vulnerability of CFST integral structures with or without rigid PZs.The influence of varying modeling methods on the structural vulnerability was quantified.Leon [17] proposed the concept of resilience of CFST composite structures and suggested that the resilience design of systems under different levels of earthquakes can be realized by determining the performance level indicators.Kamaris [18] provided a simplified seismic design method for CFST structures.Ahmadi [19] established a vulnerability model of CFST structures and calculated the failure probability of the structure under earthquake actions, which was used as the risk data of structural performance-based design.Hu [20] assessed the seismic risk of the CFDST frames.Six typical extreme performance levels of the CFDST structure were determined.The models for the specimens were established.The seismic responses of these models to seismic excitation were studied by nonlinear time history analysis, and the limit capacity were determined by incremental dynamic analysis.
Although many different types of shear walls have been proposed, most of these shear walls are connected to the frames with complete filling.The additional stresses by the shear walls will directly act on the column, which are unfavorable to the seismic resistance of the column.In addition, although there is already research foundations for CFST frame-shear walls, there are almost no research reports on composite CFDST frame-shear walls.Based on the above existing research, we propose a CFDST frame-BRW structure.
This study aims to define and describe three types of performance limit states applicable to this type of structure through the typical failure characteristics of frame-BRW structures.Based on the test-validating method, the numerical analysis models of specific composite frame-BRW structures are established.The regional representative seismic waves are selected, and the nonlinear dynamic analysis of the nine-story structure is carried out.Then, based on the theory of seismic risk assessment, a structural probability calculation model is constructed to evaluate the seismic risk of the structure quantitatively.The results will provide ideas and risk data for the performance-based seismic design of assembled CFDST frame-BRW structures.

Cyclic Tests and Performance Limits
A total of two CFDST-BRW specimens were designed and tested.The BRW is the reinforced concrete shear wall.There were two test specimens investigated in this test.
These specimens were scaled down to half and were fabricated using the prefabrication method proposed in Ref. [20].Details of these specimens are shown in Figure 1.
In the elastic stage, the CFDST frame and the BRW did not show obvious deformation.The state changed when the story drift ratio (SDR) was close to 1.0%, and the experimental curves showed an inflection point.Then, test specimens entered the elastic-plastic state.The cracks occurred on the BRW at this stage, resulting in the nonlinear behaviors of the test specimens.The specimens were significantly deformed as the SDR reached 3.0% and entered the failure stage.The loading capacity of all CFDST frame-BRW specimens began to degrade significantly, and the BRW gradually failed.The weld at the CFDST joint was torn, the obvious buckling deformation occurred on the stiffener of the CFDST joint and the plate at the beam flange, and the plate at the bottom of the column sustained local buckling.When the SDR reached 6.0%, the ultimate loading capacity of all specimens was obviously degraded, and all of the BRW failed.The CFTSD joints and column bottoms in the CFDST frame form obvious plastic hinges, and the structure was about to lose its bending resistance and cause overall deformation.The envelop curves of the specimens are shown in Figure 2. In the elastic stage, the CFDST frame and the BRW did not show obvious d mation.The state changed when the story drift ratio (SDR) was close to 1.0%, and experimental curves showed an inflection point.Then, test specimens entered the ela plastic state.The cracks occurred on the BRW at this stage, resulting in the nonlinea haviors of the test specimens.The specimens were significantly deformed as the reached 3.0% and entered the failure stage.The loading capacity of all CFDST frame-B specimens began to degrade significantly, and the BRW gradually failed.The weld a CFDST joint was torn, the obvious buckling deformation occurred on the stiffener o CFDST joint and the plate at the beam flange, and the plate at the bottom of the col sustained local buckling.When the SDR reached 6.0%, the ultimate loading capacity specimens was obviously degraded, and all of the BRW failed.The CFTSD joints and umn bottoms in the CFDST frame form obvious plastic hinges, and the structure about to lose its bending resistance and cause overall deformation.The envelop curv the specimens are shown in Figure 2.There are three typical mechanical phases.
(1) Stage I: from 0 to 1% drift.The CFDST frame-BRW structure can be approxim considered to be in the elastic stage.Although the BRW has small cracks, the d In the elastic stage, the CFDST frame and the BRW did not show obvious deformation.The state changed when the story drift ratio (SDR) was close to 1.0%, and the experimental curves showed an inflection point.Then, test specimens entered the elasticplastic state.The cracks occurred on the BRW at this stage, resulting in the nonlinear behaviors of the test specimens.The specimens were significantly deformed as the SDR reached 3.0% and entered the failure stage.The loading capacity of all CFDST frame-BRW specimens began to degrade significantly, and the BRW gradually failed.The weld at the CFDST joint was torn, the obvious buckling deformation occurred on the stiffener of the CFDST joint and the plate at the beam flange, and the plate at the bottom of the column sustained local buckling.When the SDR reached 6.0%, the ultimate loading capacity of all specimens was obviously degraded, and all of the BRW failed.The CFTSD joints and column bottoms in the CFDST frame form obvious plastic hinges, and the structure was about to lose its bending resistance and cause overall deformation.The envelop curves of the specimens are shown in Figure 2.There are three typical mechanical phases.
(1) Stage I: from 0 to 1% drift.The CFDST frame-BRW structure can be approximately considered to be in the elastic stage.Although the BRW has small cracks, the deformation has little effect on their structural performance.There are three typical mechanical phases.
(1) Stage I: from 0 to 1% drift.The CFDST frame-BRW structure can be approximately considered to be in the elastic stage.Although the BRW has small cracks, the deformation has little effect on their structural performance.(2) Stage II: from 1 to 3% drift.Many obvious cracks appeared on the BRW, and the bolt connection between the BRW and CFDST frame was damaged to a certain extent.
Local buckling of column and beam flanges of the CFDST frames was obvious, and the weld connections and flange plates near the high-strength bolts were torn.(3) Stage III: from 3 to 6% drift.The BRW was obviously deformed, and it had basically separated from the CFDST frame and lost its loading capacity.The steel beam plate and high-strength bolts near the CFDST column were torn, and the CFDST frame joints and the bottom of the column formed a plastic hinge, which has basically lost its loading capacity.
The American design codes FEMA 356 [21], HAZUS [22], and SEAOC [23] have not defined the values of the limit state of the CFDST frame-BRW structure, and other design codes have not provided relevant reference information.We used the maximum inter-story drift ratio (ISDAmax) as the performance index to describe the limit state of the building structure.Combined with the definition of the CFDST frame defined by Ref. [20], and the meaning of the performance limits of the CFDST structure in FEMA 356 performance limit of the CFDST frame-BRW structure is defined as follows.
(1) Immediate Occupancy (IO): the structure needs to be used immediately after the earthquake.This stage is described for the CFDST frame-BRW structure, i.e., the BRW has small cracks, the steel beam flange or web exhibits a slight yield or buckling phenomenon, and no obvious fracture phenomenon is observed.The performance limit is defined as 1.0%.(2) Structural Damage (SD): the structural or nonstructural components may be seriously damaged or partially collapsed, but the building still has a certain ability to resist collapse, which can ensure that the structure will not collapse.For the CFDST frame-BRW structure, this stage is described as follows.The BRW has many cracks, and the concrete at the bolt connection between the BRW and CFDST frame exhibits spalling phenomena.There are many steel plates yielding, buckling, or fracturing on the steel beam near the column, but they do not completely fail, and some members have obvious local buckling or plastic hinges.The performance limit is defined as 3.0%.(3) Structural Damage (SD): the structure collapses due to its inability to withstand its own gravity or insufficient lateral stiffness.For the CFDST frame-BRW structure, this stage is described as follows.The BRW has basically separated from the CFDST frame and has lost its loading capacity.The steel beam plate in the CFDST frame is torn, the joint concrete is crushed, and the outer steel tube is severely buckled.The performance limit is defined as 6.0%.

CFDST Framework Numerical Model and Verification
Hu established the test-validated numerical model of the CFDST framework [20].Based on this model, a nine-story CFDST frame-BRW was constructed, and the influence of a BRW on the seismic risk of a nine-story CFDST frame was studied.A seismic risk assessment was carried out.To ensure the integrity of the paper, the modeling and verification of the CFDST frame structure are still explained.

CFDT Column, Steel Beam, and CFDST Joints
The CFDST frame structure model is established by using the OpenSees.Nonlinear elements are used to simulate the performances of CFDST columns and beams, and the CFDST column filled with concrete uses the Concrete 02 material, while the steel tube and beam use the Steel 02 material.The Kent-Scott-Park damage model is used for concrete materials to simulate the hysteretic characteristics and material degradation properties accurately.Three nonlinear spring elements, i.e., a horizontal K H, a vertical K V , and a rotating K θ , are used to simulate the CFDST joint model in Figure 3, and the mechanical behaviors of the nonlinear springs used to simulate the CFDST joints and nonlinear springs are simulated by the Zerolength element in OpenSees.In the numerical model, the model of the horizontal spring corresponds to two nodes at the connection position of the beam and column, and the modeling is achieved by maintaining the same horizontal displacement of the two nodes.
According to the failure modes observed during the test, the local buckling and fracture behavior of the steel beams occurred, resulting in strength and stiffness degradations of CFDST joints.Therefore, a simple nonlinear model cannot accurately simulate the mechanical behaviors of steel beams.An energy-based improved IK steel damage model developed by Lignos [24] was selected to simulate the intricate mechanical behaviors of steel beams accurately.The model is based on the Ibarra-Krawinkler (IK) [24] and more than 300 wide-flange steel beam tests.
Three nonlinear spring elements, i.e., a horizontal KH, a vertical KV, and a rotating Kθ, are used to simulate the CFDST joint model in Figure 3, and the mechanical behaviors of the nonlinear springs used to simulate the CFDST joints and nonlinear springs are simulated by the Zerolength element in OpenSees.In the numerical model, the model of the horizontal spring corresponds to two nodes at the connection position of the beam and column, and the modeling is achieved by maintaining the same horizontal displacement of the two nodes.According to the failure modes observed during the test, the local buckling and fracture behavior of the steel beams occurred, resulting in strength and stiffness degradations of CFDST joints.Therefore, a simple nonlinear model cannot accurately simulate the mechanical behaviors of steel beams.An energy-based improved IK steel damage model developed by Lignos [24] was selected to simulate the intricate mechanical behaviors of steel beams accurately.The model is based on the Ibarra-Krawinkler (IK) [24] and more than 300 wide-flange steel beam tests.
The monotonic curve model of the improved IK model was adopted, as shown in Figure 4.The curve contains the yield moment (My), capping moment (Mc), residual moment (Mr), yield rotation (θy), pre-capping plastic rotation (θp), post-capping plastic rotation (θpc) and ultimate rotation capacity (θu).The strength and nonlinear deformation values should be appropriately reduced since the specimen is subjected to a reciprocating load.The values that need to be modified are shown as follows: (1) Mc is defined as 0.9 times the initial value but not less than My; (2) θp is determined to be 0.7 times the initial value; (3) θpc is determined to be 0.5 times the initial value.The monotonic curve model of the improved IK model was adopted, as shown in Figure 4.The curve contains the yield moment (M y ), capping moment (M c ), residual moment (M r ), yield rotation (θ y ), pre-capping plastic rotation (θ p ), post-capping plastic rotation (θ pc ) and ultimate rotation capacity (θ u ).The strength and nonlinear deformation values should be appropriately reduced since the specimen is subjected to a reciprocating load.The values that need to be modified are shown as follows: (1) M c is defined as 0.9 times the initial value but not less than M y ; (2) θ p is determined to be 0.7 times the initial value; (3) θ pc is determined to be 0.5 times the initial value.
column, and the modeling is achieved by maintaining the same horizontal displacement of the two nodes.

Two nodes at this position
Three springs in the connection K H , K V and K θ According to the failure modes observed during the test, the local buckling and fracture behavior of the steel beams occurred, resulting in strength and stiffness degradations of CFDST joints.Therefore, a simple nonlinear model cannot accurately simulate the mechanical behaviors of steel beams.An energy-based improved IK steel damage model developed by Lignos [24] was selected to simulate the intricate mechanical behaviors of steel beams accurately.The model is based on the Ibarra-Krawinkler (IK) [24] and more than 300 wide-flange steel beam tests.
The monotonic curve model of the improved IK model was adopted, as shown in Figure 4.The curve contains the yield moment (My), capping moment (Mc), residual moment (Mr), yield rotation (θy), pre-capping plastic rotation (θp), post-capping plastic rotation (θpc) and ultimate rotation capacity (θu).The strength and nonlinear deformation values should be appropriately reduced since the specimen is subjected to a reciprocating load.The values that need to be modified are shown as follows: (1) Mc is defined as 0.9 times the initial value but not less than My; (2) θp is determined to be 0.7 times the initial value; (3) θpc is determined to be 0.5 times the initial value.The three-linear shear deformation model established by Skalomenos [16] modeled the PZ in the CFT column connection.In this idealized model, the total shear strength is the sum of the strengths of the steel tube and the concrete in the panel under the same deformation.To be sure of the shear yield strength and post-yield stiffness of the ideal model, this study considers the performances of the joint zone of the CFT column-wide flange beam bolted bending moment connection and the inner end plate and the CFT joint of the through-beam [25].The PZ model of CFDST column joints is established, as shown in Figure 3.A rotating spring between the beam-column represents the relative rotation.The rotating spring is modeled by the Zerolength element.Two stiffness bar elements are used to simulate the relative rotation between the two motions of beams and columns and are modeled by stiffness-link elements.The rotational spring represents the mechanical behaviors of the joint, including rotational stiffness and shear capacity.Therefore, it is essential to model the stiffness and shear capacity of the spring and then simulate the nonlinear behaviors of the PZ.By determining the stiffness and shear loading capacity of the spring, the complex force relationship of the joint can be simplified into a functional relationship of multiple line segments.Then, it can be simplified into a relatively simple mathematical calculation problem.
The horizontal shear capacity of the PZ includes the shear contributions provided by steel tube and concrete.The shear capacity Vu is the sum of the shear capacities of steel tube and concrete: where A sv and A cv are the effective shear extents of the steel tube and concrete, respectively, and f y and f c are the steel tube yield strength and the core concrete compressive strength, respectively.The core concrete is in a three-dimensional compressive stress state under the action of gravity load.Figure 5 is the mechanical model of core concrete without considering the internal and external steel pipes under gravity load.

Panel Zone
This paper establishes the panel zone (PZ) model to simulate the nonlinear behaviors of the CFDST joint panel.The PZ model should consider the influence of multiple factors, such as concrete cracking, concrete crushing, and steel yield, and the model should be based on the existing members of OpenSees.
The three-linear shear deformation model established by Skalomenos [16] modeled the PZ in the CFT column connection.In this idealized model, the total shear strength is the sum of the strengths of the steel tube and the concrete in the panel under the same deformation.To be sure of the shear yield strength and post-yield stiffness of the ideal model, this study considers the performances of the joint zone of the CFT column-wide flange beam bolted bending moment connection and the inner end plate and the CFT joint of the through-beam [25].The PZ model of CFDST column joints is established, as shown in Figure 3.A rotating spring between the beam-column represents the relative rotation.The rotating spring is modeled by the Zerolength element.Two stiffness bar elements are used to simulate the relative rotation between the two motions of beams and columns and are modeled by stiffness-link elements.The rotational spring represents the mechanical behaviors of the joint, including rotational stiffness and shear capacity.Therefore, it is essential to model the stiffness and shear capacity of the spring and then simulate the nonlinear behaviors of the PZ.By determining the stiffness and shear loading capacity of the spring, the complex force relationship of the joint can be simplified into a functional relationship of multiple line segments.Then, it can be simplified into a relatively simple mathematical calculation problem.
The horizontal shear capacity of the PZ includes the shear contributions provided by steel tube and concrete.The shear capacity Vu is the sum of the shear capacities of steel tube and concrete: ' ( ) 1.99 3 where Asv and Acv are the effective shear extents of the steel tube and concrete, respectively, and fy and f ' c are the steel tube yield strength and the core concrete compressive strength, respectively.
The core concrete is in a three-dimensional compressive stress state under the action of gravity load.Figure 5 is the mechanical model of core concrete without considering the internal and external steel pipes under gravity load.The functional relationship between the compressive strength of core concrete and the standard strength of concrete is as follows: where f c is the standard compressive strength of concrete, k c is the strength improvement coefficient of concrete under lateral restraint, and ζ is the constraint reduction coefficient, which can be calculated by Equation (3).P and P 0 are the pressures caused by the outer and inner steel tubes, respectively.B and t are the width and thickness of the outer steel tube, respectively, and D is the diameter of the inner steel tube.Figure 6 shows the idealized trilinear model of the PZ of CFDST frame joints.The yield point of the PZ is the yield point of the steel tube.
where fc is the standard compressive strength of concrete, kc is the strength improvement coefficient of concrete under lateral restraint, and ζ is the constraint reduction coefficient, which can be calculated by Equation (3).P and P0 are the pressures caused by the outer and inner steel tubes, respectively.B and t are the width and thickness of the outer steel tube, respectively, and D is the diameter of the inner steel tube.
Figure 6 shows the idealized trilinear model of the PZ of CFDST frame joints.The yield point of the PZ is the yield point of the steel tube.The calculation of the shear deformation is shown as follows: where κs1 and κs2 are the shear coefficients of the square and circular steel tubes, which are equal to 1.2 and 10/9, respectively.Asv1 and Asv2 are the shear extents of the outer and inner tubes, respectively.Gs is the shear modulus of the steel tube.When the shear force of the PZ reaches 60% of the ultimate loading capacity Vu, the PZ yields.The inelastic stiffness of the second stage K2 is 20% of the initial stiffness of the first stage K1.The bending moment-rotation (M-θ) relationship of the rotating spring can be converted from the sheardeformation relationship mentioned above.
The bending moment M and the rotation angle θ can be calculated by the following formulas as:

CFDST Frame-BRW Structure Numerical Model and Verification
In the CFDST frame-BRW structure, the BRW and CFDST frame are connected by high-strength bolts.Therefore, many nonlinear behaviors are involved in seismic analysis.The calculation of the shear deformation is shown as follows: where κ s1 and κ s2 are the shear coefficients of the square and circular steel tubes, which are equal to 1.2 and 10/9, respectively.A sv1 and A sv2 are the shear extents of the outer and inner tubes, respectively.G s is the shear modulus of the steel tube.When the shear force of the PZ reaches 60% of the ultimate loading capacity V u , the PZ yields.The inelastic stiffness of the second stage K 2 is 20% of the initial stiffness of the first stage K 1 .The bending moment-rotation (M-θ) relationship of the rotating spring can be converted from the shear-deformation relationship mentioned above.
The bending moment M and the rotation angle θ can be calculated by the following formulas as:

CFDST Frame-BRW Structure Numerical Model and Verification
In the CFDST frame-BRW structure, the BRW and CFDST frame are connected by high-strength bolts.Therefore, many nonlinear behaviors are involved in seismic analysis.Considering the use of classical finite element software ANSYS or ABAQUS to finely model the CFDST frame-BRW structure, it will not only cause the numerical model to be complicated and the calculation cost to be too high, but also may cause the numerical model to be difficult to converge due to various nonlinearities.The simplified analysis method is used to establish the numerical analysis model.The simplified method utilized in Ref. [20] is used to model the CFDST framework.We recommend using two cross-arranged nonlinear springs to consider the nonlinear behaviors of the BRW.The seismic performance simulation analysis of the CFDST frame-BRW structure is still based on the OpenSees.
The two-node link element is used to simulate the nonlinear springs, and the hysteresis parameters of the BRW are input through the Pinching04 material model.The length of the two-node link element can be zero or nonzero.The element can simulate the deformation of the structure with 1 to 6 degrees of freedom.The two-node link element used in this paper only considers the axial deformation, that is, the contribution of the BRW to the lateral stiffness of the CFDST frame-BRW structure is simulated by the axial deformation.Pinching04 is a uniaxial material model defined by 39 parameters in Figure 7.
resis parameters of the BRW are input through the Pinching04 material model.The length of the two-node link element can be zero or nonzero.The element can simulate the deformation of the structure with 1 to 6 degrees of freedom.The two-node link element used in this paper only considers the axial deformation, that is, the contribution of the BRW to the lateral stiffness of the CFDST frame-BRW structure is simulated by the axial deformation.Pinching04 is a uniaxial material model defined by 39 parameters in Figure 7.The element can characterize the deformation response of the structural member under uniaxial tension-compression loads.They can also represent the loading capacity and stiffness degradation behavior of the member under cyclic loading.The Pinching04 material model uses 16 parameters to define the skeleton curve of the members (ePf1, ePd1, ePf2, ePd2, ePf3, ePd3, ePf4, ePd4, eNf1, eNd1, eNf2, eNd2, eNf3, eNd3, eNf4 and eNd4) and uses 6 parameters to define the hysteretic characteristics of the members under a reciprocating load (uForceP, uForceN, rDispP, rDispN, rForceP and rForceN).The loading capacity degradation behaviors of the member under cyclic loading is defined by five parameters (gF1, gF2, gF3, gF4, and gF5).The unloading stiffness degradation of the member is controlled by five parameters (gK1, gK2, gK3, gK4, and gKLim).The loading stiffness degradation of the member is controlled by five parameters (gD1, gD2, gD3, gD4, and gDLim).The maximum degradation degree of each cycle is limited by two parameters (gE and dmgType).In addition, the Pinching04 material model can consider positive and negative loadings, so the Bauschinger effect of the model can be considered.For each loading direction, the skeleton curve can be defined by four parameters.Similarly, each direction can also define the hysteretic parameters, loading capacity, and stiffness degradation parameters.Because the Pinching04 material model contains many parameters, it can simulate the hysteretic behavior of complex structures under an earthquake.The element can characterize the deformation response of the structural member under uniaxial tension-compression loads.They can also represent the loading capacity and stiffness degradation behavior of the member under cyclic loading.The Pinching04 material model uses 16 parameters to define the skeleton curve of the members (ePf1, ePd1, ePf2, ePd2, ePf3, ePd3, ePf4, ePd4, eNf1, eNd1, eNf2, eNd2, eNf3, eNd3, eNf4 and eNd4) and uses 6 parameters to define the hysteretic characteristics of the members under a reciprocating load (uForceP, uForceN, rDispP, rDispN, rForceP and rForceN).The loading capacity degradation behaviors of the member under cyclic loading is defined by five parameters (gF1, gF2, gF3, gF4, and gF5).The unloading stiffness degradation of the member is controlled by five parameters (gK1, gK2, gK3, gK4, and gKLim).The loading stiffness degradation of the member is controlled by five parameters (gD1, gD2, gD3, gD4, and gDLim).The maximum degradation degree of each cycle is limited by two parameters (gE and dmgType).In addition, the Pinching04 material model can consider positive and negative loadings, so the Bauschinger effect of the model can be considered.For each loading direction, the skeleton curve can be defined by four parameters.Similarly, each direction can also define the hysteretic parameters, loading capacity, and stiffness degradation parameters.Because the Pinching04 material model contains many parameters, it can simulate the hysteretic behavior of complex structures under an earthquake.
Figure 8 is a geometric transformation relationship between the load-displacement curve of the element and the load-displacement curve of the BRW, including the force and displacement borne by the element: where F and F are the horizontal forces borne by the wall and the axial force borne by a two-node link element, respectively; ∆ w and ∆ w are the horizontal displacements of the BRW and the axial deformation of a two-node link element, respectively; θ is the angle between the two-node link element and the steel beam and represents the angle between the force and deformation direction of the two-node link and the horizontal direction.
where F and F′ are the horizontal forces borne by the wall and the axial force bor two-node link element, respectively; Δw and Δw′ are the horizontal displacements BRW and the axial deformation of a two-node link element, respectively; θ is th between the two-node link element and the steel beam and represents the angle b the force and deformation direction of the two-node link and the horizontal direct After selecting the element and material model of the simulated BRW, it is ne to ensure the parameters in the material model, in which the skeleton curve is rel the loading capacity and deformation capacity of the BRW.To match the Pinching terial model, this paper recommends a fivefold line model for the skeleton curve BRW connected to the beam end.There are four key points and five stages in Fi The stage is known as the elastic stage, in which the BRW concrete does not crack, high-strength bolted connection between the BRW and the CFDST frame is not tor The loading capacity Py and the displacement Dy at point A are defined as the l capacity corresponding to the cracking of concrete in the BRW, and the displace the ratio of the loading capacity to the elastic lateral stiffness of the BRW.Because After selecting the element and material model of the simulated BRW, it is necessary to ensure the parameters in the material model, in which the skeleton curve is related to the loading capacity and deformation capacity of the BRW.To match the Pinching04 material model, this paper recommends a fivefold line model for the skeleton curve of the BRW connected to the beam end.There are four key points and five stages in Figure 9.The stage is known as the elastic stage, in which the BRW concrete does not crack, and the high-strength bolted connection between the BRW and the CFDST frame is not torn.
two-node link element, respectively; Δw and Δw′ are the horizontal displacements o BRW and the axial deformation of a two-node link element, respectively; θ is the a between the two-node link element and the steel beam and represents the angle betw the force and deformation direction of the two-node link and the horizontal direction After selecting the element and material model of the simulated BRW, it is neces to ensure the parameters in the material model, in which the skeleton curve is relate the loading capacity and deformation capacity of the BRW.To match the Pinching04 terial model, this paper recommends a fivefold line model for the skeleton curve o BRW connected to the beam end.There are four key points and five stages in Figu The stage is known as the elastic stage, in which the BRW concrete does not crack, an high-strength bolted connection between the BRW and the CFDST frame is not torn.The loading capacity Py and the displacement Dy at point A are defined as the loa capacity corresponding to the cracking of concrete in the BRW, and the displaceme the ratio of the loading capacity to the elastic lateral stiffness of the BRW.Because the The loading capacity P y and the displacement D y at point A are defined as the loading corresponding to the cracking of concrete in the BRW, and the displacement is the ratio of the loading capacity to the elastic lateral stiffness of the BRW.Because the test model involves the initial defects of the BRW, the uncertainty of the material properties, and the bolt connection between the BRW and the CFDST frame, the ideal state theory cannot sufficiently reflect the stress state of the CFDST frame-BRW structure.To this end, this paper draws on the simplified processing method adopted by the energy-based improved IK damage model proposed by Lignos [26].It uses the test data to determine the skeleton curve's loading capacity and deformation parameters, represented by a simple constant coefficient.According to the analysis of the results of the two test models, the following key points are defined.
(1) Point A is defined as follows, i.e., the loading capacity P y of the BRW is 0.59 times the peak loading capacity P u1 of the BRW.The SDR is 0.6%, that is, the displacement is 9.0 mm.(2) Point B is defined as follows, i.e., the peak loading capacity P u1 of the BRW is calculated according to the simplified calculation method.The SDR of the BRW is 1.3%, and the displacement is 19.5 mm.(3) Point C is defined as follows, i.e., the loading capacity P u2 of the BRW is 0.75 times (BF-BRW-1) and 0.85 times (BF-BRW-2) of the peak loading capacity P u1 of the BRW, and the SDRs are 2.6% (BF-BRW-1) and 3.3% (BF-BRW-2), that is, the displacements are 39.0 mm (BF-BRW-1) and 49.5 mm (BF-BRW-2), respectively.(4) Point D is defined as follows, i.e., the loading capacity P m of the BRW is 0.27 times that of the peak loading capacity P u1 of the BRW.The SDR is 5.7%, that is, the displacement is 85.5 mm.
When the skeleton curve parameters of the BRW are determined, the skeleton curve parameters of the BRW are converted into the skeleton curve parameters of the two-node link element according to Equations ( 8) and ( 9) and are input into the simplified model of the structure to realize the simulation analysis of the test model.Among them, ePf1, ePf2, ePf3, and ePf4 correspond to the forces at points A, B, C, and D after the geometric relationship is transformed; ePd1, ePd2, ePd3, and ePd4 correspond to the displacements at points A, B, C and D after geometric transformation, respectively.
where P i max represents the peak load of the specimen under the ith reciprocating load; P 1 max denotes the peak load of the specimen in the 1 st loading stage; K i is the secant stiffness of the test specimen in the ith loading stage; P i j indicates the peak load capacity of the test specimen in the jth direction (forward loading and reverse loading) when the ith level is loaded; and ∆ i j indicates the displacement corresponding to the peak loading capacity of the specimen in the j direction when the ith stage is loaded.
In addition to the forces and displacements at points A, B, C, and D, the hysteretic parameters, loading capacity, and stiffness degradation parameters of the two-node link element need to be determined.Analysis of the hysteretic curves measured in the test, the hysteretic parameters, loading capacity, and stiffness degradation parameters of the two-node link element can be obtained as follows.
The above parameters are input into the numerical analysis model, and the CFDST frame-BRW structure is numerically simulated and analyzed.The two CFDST frame-BRW structure tests completed are used to verify the numerical model as shown in Figure 10.
Since the CFDST frame connection and the BRW are represented by multiline segment functions, the predicted hysteresis curves show the prominent line segment shape characteristics.The simplified numerical model not only simulates the ultimate loading capacity of the test model but also has high for the degradation behaviors of the stiffness and loading capacity of the test model under cyclic loading.This shows that the model can accurately predict the seismic performances of the CFDST frame-BRW structure.In addition, compared with traditional finite element analysis software, such as ANSYS and ABAQUS, the calculation time of a single model using the simplified analysis method is less than 5 min, the calculation efficiency is high, and the convergence is good.Therefore, this method can be considered as an effective numerical method for seismic analysis and dynamic calculation of CFDST frame-BRW structures.
The above parameters are input into the numerical analysis model, and th frame-BRW structure is numerically simulated and analyzed.The two CFDST fra structure tests completed are used to verify the numerical model as shown in Fig Since the CFDST frame connection and the BRW are represented by multi ment functions, the predicted hysteresis curves show the prominent line segme characteristics.The simplified numerical model not only simulates the ultimate capacity of the test model but also has high accuracy for the degradation behavio stiffness and loading capacity of the test model under cyclic loading.This shows model can accurately predict the seismic performances of the CFDST frame-BR ture.In addition, compared with traditional finite element analysis software, suc SYS and ABAQUS, the calculation time of a single model using the simplified method is less than 5 min, the calculation efficiency is high, and the convergence Therefore, this method can be considered as an effective numerical method fo analysis and dynamic calculation of CFDST frame-BRW structures.

CFDST Frame-BRW for Seismic Risk Assessment
The earthquake risk assessment theory includes the following step: (1) collec motion data; (2) dynamically analyze the target building and obtain its response; ( ure the structural performance and determine the degree of damage; (4) establi nerability curve; and (5) calculate the seismic risk and make design decisions.
According to Ref. [20], the probability of the limit state in the seismic risk ass generated the following formula:

CFDST Frame-BRW for Seismic Risk Assessment
The earthquake risk assessment theory includes the following step: (1) collect ground motion data; (2) dynamically analyze the target building and obtain its response; (3) measure the structural performance and determine the degree of damage; (4) establish a vulnerability curve; and (5) calculate the seismic risk and make design decisions.
According to Ref. [20], the probability of the limit state in the seismic risk assessment generated the following formula: where P[LS i ] is the occurring probability of each limit state in seismic evaluation; LS i is the ith performance limit state of the structure; k and k 0 are the shape parameter and proportional parameter of the seismic hazard curve, respectively; and m R β R are the median and logarithmic standard deviation of seismic demand, respectively.

Earthquake Ground Motions
The seismic load and structural resistance of existing buildings are uncertain.The uncertainty of the seismic demand is dominant in the overall response.Therefore, The parameters of the numerical model are constants such as the yield stress and modulus.As China's earthquake disaster map and regional representative seismic wave selection are still in the research and development stage, this paper takes the Los Angeles area of the US as an example and selects 40 ground motions.These ground motions are representative ground motions set by the SAC project [22] for the Los Angeles (LA) area for seismic risk analysis of steel building structures on site.These records represent two risk levels: the excess probability within 50 years is 10% (la01-20) and 2% (la21-40), each group containing US as an example and selects 40 ground motions.These ground motions are representative ground motions set by the SAC project [22] for the Los Angeles (LA) area for seismic risk analysis of steel building structures on site.These records represent two risk levels: the excess probability within 50 years is 10% (la01-20) and 2% (la21-40), each group containing 20 ground motions.The response spectrum of ground motions la01 to la40 in Figure 11 represents two seismic risk levels in the Los Angeles area.

Prototype Structures
In this paper, the nine-story CFDST frame designed in Ref. [20] is used as the base model of the CFDST frame-BRW structure, as shown in Figure 12.For the CFDST frame-BRW structure, only a BRW is filled in the middle span, and the width of CFDST-BRW model is 3 m.The specific parameters are in Table 1.

Prototype Structures
In this paper, the nine-story CFDST frame designed in Ref. [20] is used as the base model of the CFDST frame-BRW structure, as shown in Figure 12.For the CFDST frame-BRW structure, only a BRW is filled in the middle span, and the width of CFDST-BRW model is 3 m.The specific parameters are in Table 1.
Buildings 2023, 13, x FOR PEER REVIEW 12 of 18 US as an example and selects 40 ground motions.These ground motions are representative ground motions set by the SAC project [22] for the Los Angeles (LA) area for seismic risk analysis of steel building structures on site.These records represent two risk levels: the excess probability within 50 years is 10% (la01-20) and 2% (la21-40), each group containing 20 ground motions.The response spectrum of ground motions la01 to la40 in Figure 11 represents two seismic risk levels in the Los Angeles area.

Prototype Structures
In this paper, the nine-story CFDST frame designed in Ref. [20] is used as the base model of the CFDST frame-BRW structure, as shown in Figure 12.For the CFDST frame-BRW structure, only a BRW is filled in the middle span, and the width of CFDST-BRW model is 3 m.The specific parameters are in Table 1.The gravity load adopts the analysis example of Skalomenos [16], which is assumed to be G + 0.3Q = 27.5 kN/m.The yield strength of the CFDST column and the steel strength in the BRW is 275 MPa.The compressive strength of concrete and filled concrete in the BRW is 20 Mpa.

Probabilistic Seismic Demands
The nonlinear dynamic analysis is carried out to establish the numerical analysis model.The selected representative seismic waves (la01-la40) are input into the model one by one according to 10%/50 year and 2%/50 year groups, and the seismic demand results of the model can be obtained.The ISDA max is used as a parameter to evaluate the structural performance.
The functional relationship between the ISDA max of each floor and the spectral acceleration S a (T 1 ) [27] of the CFDST frame-BRW model obtained by nonlinear dynamic analysis is recorded as follows: where ε is a random variable; θ max is the ISDA max obtained by NTHA; a and b are constants; S a is the spectral acceleration corresponding to the first natural vibration period of the structure when the damping ratio is 5%.

Probabilistic Seismic Demands
The seismic vulnerability curve is the main part of structural seismic risk assessment, representing the functional relationship between the maximum dynamic response of the structure under an earthquake and the intensity of ground motion.The incremental dynamic analysis (IDA) of the CFDST frame-BRW is carried out using the selected seismic waves.Then, the curve relationship between the ground motion intensity and the dynamic response of each model is obtained.To reduce the calculation amount of dynamic analysis, each ground motion gradually increases the input ground motion intensity at a multiple of 100 gal and is gradually applied to the structure to record the relationship between the ISDA max and the input seismic intensity.The median value mC of the seismic demand of the three structural models under different performance states is obtained by counting these IDA curves.In addition, the seismic vulnerability curve also considers many uncertain factors, including the uncertainty of earthquakes, the uncertainty of structures, and the uncertainty of modeling, which correspond to the logarithmic standard deviations, including β C for the uncertainty of earthquakes, β SC for the uncertainty of structures and β u for the uncertainty of modeling; β SC can be combined with β u to form β C = β SC 2 + β u 2 .It is called the logarithmic standard deviation, corresponding to the uncertainty of seismic demand.
According to these median values and uncertainty parameters, the seismic vulnerability of each model to various performance states can be calculated by: where the demand uncertainty β D is determined by the dynamic response of a series of ground motions, which is equal to σ lnε ; m C is the median of earthquake demand obtained by IDA; and β C is the logarithmic standard deviation of the earthquake demand group.The seismic vulnerability of the CFDST frame-BRW model calculated in this paper only considers the influence of seismic uncertainty.It does not consider the impact of other uncertain factors, such as structural uncertainty and modeling uncertainty.These uncertainties need to be completed by carrying out many random numerical analyses.The authors will study the influence of the BRW on the seismic risk of CFDST frames.Therefore, only the uncertainty factors of earthquakes that have the greatest impact on them are selected.In the future, the authors will further study the seismic vulnerability and seismic risk of CFDST structures under multiple uncertainties and other damping structures [28].
The CFDST-BRW model's seismic vulnerability is calculated using Equation (12).The fragility curves corresponding to each limit state in Figure 13 and the corresponding fragility parameters are shown in Table 2.
uncertain factors, such as structural uncertainty and modeling uncertainty.These uncertainties need to be completed by carrying out many random numerical analyses.The authors will study the influence of the BRW on the seismic risk of CFDST frames.Therefore, only the uncertainty factors of earthquakes that have the greatest impact on them are selected.In the future, the authors will further study the seismic vulnerability and seismic risk of CFDST structures under multiple uncertainties and other damping structures [28].Because multiple uncertainties are essential for accurately assessing structures' seismic risk, many studies have discussed this issue [29][30][31][32][33][34][35][36][37][38].
The CFDST-BRW model's seismic vulnerability is calculated using Equation (12).The fragility curves corresponding to each limit state in Figure 13 and the corresponding fragility parameters are shown in Table 2.

Probability Estimation of Annual Performance Limits
The HAZUS software released by the USGS of the US Geological Survey was used to calculate the seismic disaster parameters of the model in the Los Angeles area of 10%/50 year and 2%/50 year earthquakes, as shown in Table 3.According to these earthquake disaster parameters, the shape parameter k and the proportional parameter k0 of CFDST frame-BRW structures in Los Angeles can be obtained.The shape parameter k of the CFDST frame-BRW structure in this area is not largely different from the shape parameter k of the CFDST frame.The shape parameter k and the proportional parameter k0 of the model are brought into the seismic risk function, and the seismic disaster curve of the model can be obtained, as shown in Figure 14.The annual average occurrence probability of the CFDST frame-BRW model can be calculated by convolution calculation of the seismic risk function and seismic vulnerability function.The results are shown in Table 4.

Probability Estimation of Annual Performance Limits
The HAZUS software released by the USGS of the US Geological Survey was used to calculate the seismic disaster parameters of the model in the Los Angeles area of 10%/50 year and 2%/50 year earthquakes, as shown in Table 3.According to these earthquake disaster parameters, the shape parameter k and the proportional parameter k 0 of CFDST frame-BRW structures in Los Angeles can be obtained.The shape parameter k of the CFDST frame-BRW structure in this area is not largely different from the shape parameter k of the CFDST frame.The shape parameter k and the proportional parameter k 0 of the model are brought into the seismic risk function, and the seismic disaster curve of the model can be obtained, as shown in Figure 14.The annual average occurrence probability of the CFDST frame-BRW model can be calculated by convolution calculation of the seismic risk function and seismic vulnerability function.The results are shown in Table 4.As shown in Figure 15, combined with Ref. [20], the BRW reduces the annual probability and collapse probability of the CFDST frame structure to achieve various performance states within the service life of 50 years.The BRW reduces the annual failure probability and 50-year collapse probability of the CFDST frame to a reasonable extent.The BRW reduces the 50-year collapse probability of the CFDST frame to 2.76%, significantly reducing the CFDST frame's collapse probability.The BRWs can effectively reduce the probability of different degrees of damage to the structure during the service life.As shown in Figure 15, combined with Ref. [20], the BRW reduces the annual probability and collapse probability of the CFDST frame structure to achieve various performance states within the service life of 50 years.The BRW reduces the annual failure probability and 50-year collapse probability of the CFDST frame to a reasonable extent.The BRW reduces the 50-year collapse probability of the CFDST frame to 2.76%, significantly reducing the CFDST frame's collapse probability.The BRWs can effectively reduce the probability of different degrees of damage to the structure during the service life.As shown in Figure 15, combined with Ref. [20], the BRW reduces the annual probability and collapse probability of the CFDST frame structure to achieve various performance states within the service life of 50 years.The BRW reduces the annual failure probability and 50-year collapse probability of the CFDST frame to a reasonable extent.The BRW reduces the 50-year collapse probability of the CFDST frame to 2.76%, significantly reducing the CFDST frame's collapse probability.The BRWs can effectively reduce the probability of different degrees of damage to the structure during the service life.

Conclusions
The seismic risk assessment of a nine-story CFDST frame-BRW structure was carried out to study the BRW on the seismic risk of the CFDST frame.The purpose was to calculate the exceedance probability of this kind of structure to reach various performance limit states during the design life period (usually 50 years).The ground motion obtained by the SAC project represented the uncertainty of seismic demand.The following conclusions can be drawn.

1.
Based on the unified theory, the model considered the mechanical properties of concrete materials under a three-dimensional stress state under compression.A simplified analysis model of a shear wall with beam-end connections of different types of materials was established.The hysteretic characteristics of a BRW under earthquake action were simulated by two nonlinear springs with Pinching04 material properties.The correctness of the adopted simplified analysis model was verified by the seismic performance test results of two CFDST frame-BRW structures.

2.
A numerical analysis model of a nine-story CFDST frame-BRW structure was established, and ground motion records were selected from the SAC project to describe the seismic risks in the Los Angeles area of the US States, which were divided into two risk levels of 2%/50 years and 10%/50 years.This model uses nonlinear time history analysis and incremental dynamic analysis on the selected records.

3.
The analysis model's fragility curves were established using three determined performance limits and structural responses.The annual probability of these three limits and the probability of collapse within 50 years were determined by convoluting the fragility with the seismic risk specified by the USGS.4.
The BRW reduced the probability of collapse of the CFDST frame to 2.76% in 50 years, which indicates the proposed BRW could effectively reduce the probability of different degrees of damage in the service cycle of the structure.

Figure 2 .
Figure 2. Envelop curves of the test specimens.

Figure 2 .
Figure 2. Envelop curves of the test specimens.

Figure 2 .
Figure 2. Envelop curves of the test specimens.

Figure 3 .
Figure 3.The PZ model of CFDST frame beam-to-column intermediate joints.

Figure 3 .
Figure 3.The PZ model of CFDST frame beam-to-column intermediate joints.

Figure 3 .
Figure 3.The PZ model of CFDST frame beam-to-column intermediate joints.

Figure 4 .
Figure 4.The monotone curve of the improved IK model.According to the CFST frame-BRW structure test, The main failure behavior of the beam is located near the column.An improved IK model simulates the degradation behavior of steel beams by using a rotating spring.The vertical spring is modeled as being fully elastic.Therefore, the steel beam is transformed into a bilinear strengthening model.3.1.2.Panel Zone This paper establishes the panel zone (PZ) model to simulate the nonlinear behaviors of the CFDST joint panel.The PZ model should consider the influence of multiple factors, such as concrete cracking, concrete crushing, and steel yield, and the model should be based on the existing members of OpenSees.

Figure 5 .
Figure 5.The stress state of concrete filled in the CFDST column.

Figure 8 .
Figure 8. Relationship between the force and deformation.

Figure 9 .
Figure 9. Definition of the skeleton curve for the BRW.

Figure 8 .
Figure 8. Relationship between the force and deformation.

Figure 8 .
Figure 8. Relationship between the force and deformation.

Figure 9 .
Figure 9. Definition of the skeleton curve for the BRW.

Figure 9 .
Figure 9. Definition of the skeleton curve for the BRW.

2 Figure 10 .
Figure 10.Comparison of FEM and test results.

Figure 10 .
Figure 10.Comparison of FEM and test results.
. The response spectrum of ground motions la01 to la40 in Figure11represents two seismic risk levels in the Los Angeles area.Buildings 2023, 13, x FOR PEER REVIEW 12 of 18

Figure 11 .
Figure 11.Response spectra in the Los Angeles region.

Figure 11 .
Figure 11.Response spectra in the Los Angeles region.

Figure 11 .
Figure 11.Response spectra in the Los Angeles region.

Figure 12 .
Figure 12.Schematic diagram of the CFDST frame-BRW structure.Figure 12. Schematic diagram of the CFDST frame-BRW structure.

Figure 12 .
Figure 12.Schematic diagram of the CFDST frame-BRW structure.Figure 12. Schematic diagram of the CFDST frame-BRW structure.
thickness of reinforced concrete is 120 mm, and the spacing between steel bars is 206.130mm.It is distributed in two layers along the thickness direction of the plate.The diameter of the steel bars is 10 mm.The thickness of the protective layer is 20 mm 1.415

Figure 13 .
Figure 13.Fragility curves of the models associated with the LA ensemble.

Figure 13 .
Figure 13.Fragility curves of the models associated with the LA ensemble.
Probability of collapse in 50 years CFDST (a) Annual probability of occurrence (b) The probability of collapse occurring in 50 years

Figure 15 .
Figure 15.Seismic risk probability of the CFDST frame shear wall model.

Figure 14 .
Figure 14.Seismic hazard curves for Los Angeles.

Figure 14 .
Figure 14.Seismic hazard curves for Los Angeles.
Annual probability of occurrence (b) The probability of collapse occurring in 50 years

Figure 15 .
Figure 15.Seismic risk probability of the CFDST frame shear wall model.Figure 15.Seismic risk probability of the CFDST frame shear wall model.

Figure 15 .
Figure 15.Seismic risk probability of the CFDST frame shear wall model.Figure 15.Seismic risk probability of the CFDST frame shear wall model.

Table 1 .
First-order natural vibration period of the BRW model.

Table 2 .
Fragility parameters for the models in the Los Angeles region.

Table 2 .
Fragility parameters for the models in the Los Angeles region.

Table 3 .
Seismic hazard parameters for the Los Angeles region.

Table 3 .
Figure 14.Seismic hazard curves for Los Angeles.Seismic hazard parameters for the Los Angeles region.

Table 4 .
Annual probability of limit state occurrence for the models in Los Angeles.

Table 4 .
Annual probability of limit state occurrence for the models in Los Angeles.

Table 3 .
Seismic hazard parameters for the Los Angeles region.

Table 4 .
Annual probability of limit state occurrence for the models in Los Angeles.