Carbon Fiber-Reinforced Polymer Composites Integrated Beam–Column Joints with Improved Strength Performance against Seismic Events: Numerical Model Simulation

: Strength enhancement of non-seismic concrete beam–column joints (NSCBCJs) via carbon fiber-reinforced polymer composites (CFRPCs) integration has become a viable strategy. However, the implementation of these NSCBCJs without transverse reinforcement shows poor performance during earthquakes in seismic locations. Thus, strengthening the anti-seismic performance of NSCBCJs to meet the acceptance criteria of ACI 374.1-05 is fundamentally significant. Yet, in addition to limited experimental results, only a few numerical studies based on the finite element model have been performed to determine the anti-seismic behavior of NSCBCJs. Consequently, the stress contribution of CFRPCs to NSCBCJs is not clearly understood. Therefore, we used a finite element model to examine the strength contribution of CFRPCs to NSCBCJs. The performance of the proposed finite element model was validated using the experimental results, demonstrating a good agreement between them. It was shown that the strength of NSCBCJs was improved due to CFRPC incorporation, thereby achieving compliance with the seismic requirements of ACI 374.1-05. In addition, CFRPCs presence could enhance the confinement, reduce the deformation of the NSCBCJs and, thus, decrease their stiffness and strength degradation, while simultaneously improving the energy dissipation.


Introduction
In the context of construction technology, the term "strengthening" refers to the retrofit, rehabilitation, and repair of building materials.In retrofitting, the structural components are installed with the goal of improving their overall performance in terms of shear, flexural, and compressive strength [1][2][3].This improvement in the structural integrity becomes more significant than the initial performance.Rehabilitation is the process of reparation and restoration of the performance and strength lost in the structures under varied situations [4][5][6][7].In the repairing process, the structural members' quality are improved to meet the desired requirements.Over the decades, numerous experimental and analytical studies have been conducted to gain an in-depth understanding regarding the behavior of beam-to-column joints (BCJs).In the reinforced concrete structures, the critical zones are identified at the BCJ wherein the vertical and lateral loads meet and are transferred to the footing [8][9][10].Consequently, the risk of failure at the joints is much higher compared to other structural components.During an earthquake, the plastic hinge mechanism (PHM) can potentially occur near the joints.Essentially, BCJs are a segment of columns shared by the beams at their intersections in reinforced concrete buildings [11][12][13].Needless to say that with the increase in the frequency of earthquakes worldwide, the sustainability of building structures in the civil engineering sectors has become a major concern unless steps are taken to inhibit this damage.It has been realized that the implementation of stable and durable BCJs with high performance is essential in the construction process [14].To determine the legitimacy of such joints, in-depth experimental, analytical, and numerical modelling and simulation studies are vital.
In recent years, many details of the BCJs have been examined to protect the BCJs from failure during seismic vibration.In this regard, impacts of several critical factors, such as the degree of steel reinforcement in various types of concretes [15], changes in the longitudinal reinforcement [16] and shear reinforcement [17][18][19][20], and the use of rebar at different diameters [21] on the BCJs were evaluated.It is established that diverse strengthening methods, like the incorporation of different carbon fiber sheets (CFRP and glass fiber-reinforced polymer (GFRP)) [22][23][24], steel plate jacketing [25], and external post-tension rods [26], can effectively be used to improve the anti-seismic behavior of BCJs.The results of several experimental and theoretical studies on the BCJs concerning influencing factors and strengthening methods have significantly contributed towards their anti-seismic performance.
The concept of using fiber-reinforced polymer (FRP) confinement involves bonding the FRP to the concrete surface with a bonding agent, such as epoxy resin.This bond allows the FRP and concrete to function as a composite structure, enhancing the BCJs' seismic performance when properly designed and constructed [27,28].FRP is an anisotropic material with superior engineering properties along the direction of the reinforcing fiber, but it is weak in the transverse direction.Therefore, it is highly recommended to confine the structure in the principal fiber direction of the FRP [29][30][31].The application process of FRP begins with cleaning the concrete surface, applying a layer of primer and epoxy, and then laying up the FRP sheets with epoxy until the required number of layers is achieved.A variety of FRP materials have been used to retrofit BCJs, including CFRP [32][33][34], GFRP [35,36], aramid fiber-reinforced polymer [37], basalt fiber-reinforced polymer (BFRP) [38], FRP hybrid composites [39], quasi-isotropic laminates [40], and sprayed FRP [41].Among these composite materials, CFRP is the most common due to its superior mechanical properties [42].
The application of CFRP has become a recognized method for strengthening or rehabilitating existing reinforced concrete elements (RCEs).Many studies have reported the excellent efficacy of CFRP in enhancing the flexural performance of RC beams under monotonic loading [29,43,44].The main elements of the CFRP-concrete bonded system used for external reinforcement of RC structures include CFRP dry fiber sheets, adhesive, and anchors when necessary [45].It is generally assumed that complete composite action occurs between the concrete elements, like beams, and the CFRP materials.However, the bonding quality relies heavily on the shear stiffness and strength of the adhesive, which must be sufficient to transfer shear forces between the CFRP and the concrete substrate [46][47][48].Furthermore, the bond quality is notably affected by the properties of the concrete and the preparation of the surface.
In study conducted by Laseima et al. [49], the authors have been examined the seismic behavior of exterior RC-BCJs strengthened with CFRP.From the obtained results and compared to control specimens, the authors reported that the strengthened system improved by 62%, 209%, and 61.9% in terms of strength, energy dissipation, and ductility, respectively.In [50], the authors examined the impact of using CFRP with different anchors on the seismic performance of RC-BCJs assemblies.The results showed that the ultimate load-carrying capacity increased by between 19 and 43.6%, depending on the anchoring configurations, compared to the control specimens.Additionally, all strengthened specimens demonstrated greater total energy dissipation and damping ratios than the control specimens.
In a recent study conducted by Karayannis and Golias [51], the authors investigated full-scale RC-BCJs that were externally reinforced with CFRP ropes.The joints were strengthened using X-shaped ropes, either singly or doubly, on each side of the joint, while single or double straight ropes were applied to each side of the beam.The strengthened specimens demonstrated an improved hysteretic response compared to the unstrengthened ones.Enhancements were observed in maximum loads per loading step, stiffness, and energy dissipation.The maximum strength of the reinforced samples increased by 20% to 57%, depending on the number of CFRP ropes used.
In experiments conducted by Lim et al. [52], four RC-BCJs were retrofitted utilizing two different techniques involving CFRP grids covered with engineered cementitious composite materials and high-strength mortar.The dimensions of the retrofitted specimens were preserved by removing the concrete cover of the joint and applying a cementitious matrix.The results showed that the failure mode of the specimens could be modified by redistributing the load concentration from the joint to the beam.Additionally, the overall performance of the specimens in terms of strength improved, exhibiting greater ductility and, therefore, delaying failure.
For instance, Atari et al. [53] strengthened ten BCJs through different types of fibers and evaluated their influence on the ductility, durability, strength performance, and energy dissipation capacity against earthquakes.Golias et al. [54] examined six external joints subjected to seismic loads, where three of them were strengthened with carbon fiber ropes of different shapes.In addition, Mady et al. [55] designed five T-shaped joints to withstand seismic loads wherein the effects of some vital factors, like the details of the reinforcement with carbon fiber GFRP in strengthening the performance were determined.Bsisu and Hiari et al. [56] theoretically studied the behavior and strength of CFRPGreinforced joints to examine the feasibility of improving their performance against seismic loading.Ercan et al. [57] used four different strengthening patterns for the joints to examine their anti-seismic response.These joints were strengthened from the inside by changing the reinforcement details and from the outside with carbon fiber CFRPCs.
In early structural designs, the significance of BCJs is frequently overlooked, making the retrofitting of these joints crucial for structural stability.Over time, external confinement has proven to be one of the most effective retrofit methods.Nevertheless, its practical application is hindered by several factors, such as the absence of comprehensive guidelines for analysis, design, and confinement schemes.Few studies have examined the seismic performance of reinforced concrete structures enhanced with carbon fiber-reinforced polymers.
Considering the basic significance of and applied interest in high-performance CBCJs, we performed a numerical study to determine the possibility of strengthening NSCBCJs by integrating CFRPCs.The results obtained from the numerical model simulation were compared with the state-of-the-art experimental findings available in the literature.In addition, more details about the selected specimens and evaluation methods are presented.The results revealed NSCBCJ strength enhancement because of CFRPC integration, fulfilling the compliance with the seismic requirements standard ACI 374.1-05.

Experimental Work
Following the recommendation of Trung et al. [58], this work selected four concrete specimens to examine the anti-seismic behavior of the proposed NSCBCJ-integrated CBCJs.The first specimen was designed using non-seismic details.The second specimen was made based on the first specimen details of CFRPC to strengthen the joint, beam area, and volume adjacent to the joint.Those specimens had identical dimensions and degrees of reinforcement. Figure 1 illustrates the reinforcement details of specimens for the nonseismic (NS) and seismic design (SD).The designed specimens (NS and SD) had identical dimensions and the same number of longitudinal reinforcements.The NS specimen differed from the SD specimen in that it lacked transverse reinforcements in the beam-column joint region, had relatively large stirrup spacing in the beam, and featured down-bending anchorage of beam bars away from the joints.In contrast, the SD specimen was designed to meet the requirements for an intermediate moment frame as specified by the ACI 318-02 code [59], which included confinement reinforcement in the beam-column joint region and relatively small stirrup spacing in the beam.Comparing the behaviors of these two specimens under cyclic loading can reveal the impact of adequate transverse reinforcement on the shear capacity and ductility of the joints.Figure 2a shows the experimental design for testing the specimens.The joint supports were fixed to the columns and beams at both ends to indicate the inflection points in the element's central span.Next, the joint support was attached under the column to the strong floor.The roller support was connected at the end of the beam to enable rotation and longitudinal movements only.In contrast, the column base was connected to a hinge support.The load on the specimens (at the end of the column) was applied using a 500 kN actuator.The cyclic load test was conducted under displacement control, as described by Trung [58] in the figure below (Figure 2b).A displacement-controlled quasi-static analysis was performed in ABAQUS using the dynamic implicit method, which is optimal for simulating the nonlinear behavior of beamcolumn joints.The default full Newton solving algorithm was utilized.The mesh was refined to increase node points, ensuring a reasonable mesh density in the joint areas with an acceptable aspect ratio of the elements.The two materials shared the same nodes because the steel reinforcing truss element was connected between the nodes of each concrete element to create an embedded bond.The embedded bond assumption used in the structural modeling did not result in significant errors in the predicted load and displacement response.In the CFRP specimen, full contact was simulated at the interface between the CFRP shell elements and the solid concrete elements.
Eng 2024, 5, FOR PEER REVIEW 4 these two specimens under cyclic loading can reveal the impact of adequate transverse reinforcement on the shear capacity and ductility of the joints.Figure 2a shows the experimental design for testing the specimens.The joint supports were fixed to the columns and beams at both ends to indicate the inflection points in the element's central span.Next, the joint support was attached under the column to the strong floor.The roller support was connected at the end of the beam to enable rotation and longitudinal movements only.
In contrast, the column base was connected to a hinge support.The load on the specimens (at the end of the column) was applied using a 500 kN actuator.The cyclic load test was conducted under displacement control, as described by Trung [58] in the figure below (Figure 2b).A displacement-controlled quasi-static analysis was performed in ABAQUS using the dynamic implicit method, which is optimal for simulating the nonlinear behavior of beam-column joints.The default full Newton solving algorithm was utilized.The mesh was refined to increase node points, ensuring a reasonable mesh density in the joint areas with an acceptable aspect ratio of the elements.The two materials shared the same nodes because the steel reinforcing truss element was connected between the nodes of each concrete element to create an embedded bond.The embedded bond assumption used in the structural modeling did not result in significant errors in the predicted load and displacement response.In the CFRP specimen, full contact was simulated at the interface between the CFRP shell elements and the solid concrete elements.Eng 2024, 5, FOR PEER REVIEW 4 these two specimens under cyclic loading can reveal the impact of adequate transverse reinforcement on the shear capacity and ductility of the joints.Figure 2a shows the experimental design for testing the specimens.The joint supports were fixed to the columns and beams at both ends to indicate the inflection points in the element's central span.Next, the joint support was attached under the column to the strong floor.The roller support was connected at the end of the beam to enable rotation and longitudinal movements only.
In contrast, the column base was connected to a hinge support.The load on the specimens (at the end of the column) was applied using a 500 kN actuator.The cyclic load test was conducted under displacement control, as described by Trung [58] in the figure below (Figure 2b).A displacement-controlled quasi-static analysis was performed in ABAQUS using the dynamic implicit method, which is optimal for simulating the nonlinear behavior of beam-column joints.The default full Newton solving algorithm was utilized.The mesh was refined to increase node points, ensuring a reasonable mesh density in the joint areas with an acceptable aspect ratio of the elements.The two materials shared the same nodes because the steel reinforcing truss element was connected between the nodes of each concrete element to create an embedded bond.The embedded bond assumption used in the structural modeling did not result in significant errors in the predicted load and displacement response.In the CFRP specimen, full contact was simulated at the interface between the CFRP shell elements and the solid concrete elements.

Element Type
The steel reinforcing bar and concrete are the two main components of any reinforced CBCJs.Moreover, the key components of the retrofitted CBCJs are the concrete, steel reinforcement bars, and external CFRP sheet.To introduce a realistic model under cyclic force, it is necessary to accurately replicate the real material characteristics of every component.The ABAQUS 19 library (computer software, user's manual, Providence, RI, Rising Sun Mills) offers effective material models for three-dimensional (3D) continuum elements that can simulate the actual behavior of each component with acceptable accuracy.Using three different element types, including solid, shell, and link, various geometries were modeled.We used a 3D eight-node linear brick element with reduced integration and hourglass control (C3D8R) for the solid section, like concrete (Figure 3a).The linear brick components were chosen because they can be utilized with contact, unlike the quadratic brick elements that require more time to compute the constant nodal loads over the surface.The truss elements T3D2 (Figure 3b) used for the steel bars were embedded in the concrete.The discrete reinforcement bars in the CBCJ are modeled as T3D2 elements.Four-node, doubly curved, linear, three-dimensional general-purpose shell elements with decreased integration and hourglass control (S4R5) were used in the CFRP plate (Figure 3c).Since the nodes of a conventional shell element are situated in a well-defined planar dimension, the thickness of this shell element can be determined by the section properties.The conventional shell elements are regarded to be more accurate in contact modelling than continuum shell elements because they can measure strain or slip without affecting the thickness of the CFRP composite plates.Furthermore, the conventional shell elements have exceptional computational efficiency.

Material Properties
In this study, the material properties used to prepare the specimens, such as concrete compressive strength and splitting tensile strength, longitudinal reinforcement and stirrups' yield strength, as well as the tensile strength and the elastic modulus of CRFP and epoxy, are summarized in Table 1.Concrete is known to exhibit a high compressive strength (CS) and low tensile strength (TS), and it is prone to fracture when subjected to tensile loads.This is mainly due to the susceptibility of aggregates in concrete that can endure compressive stress and, thus, allow it to carry compressive loads.However, cracks that separate the cement particles in the aggregates are generated due to tensile stresses.The separation of cement particles causes the entire structure to fail as cracks start propagating.The issue with concrete can be addressed via the reinforcing of elements, such as metallic bars or fibers.These components serve as the framework of the entire construction and can support the aggregates when subjected to tensile forces.This process is referred to as concrete reinforcing.The term "brittle material" might apply to concrete.The difference in behavior between concrete and ductile materials, like steel, is due to their distinct load responses.Brittle materials tend to develop tensile fracture perpendicular to the direction of the largest tensile strain.Thus, the cracks tend to develop parallel to the maximum compressive stress when a concrete specimen is subjected to a uniaxial compressive load.In compression, the stress-strain curve of concrete being elastic can gradually reach up to the maximum CS followed by a decrease into a softening region and eventually crushing failure at an ultimate strain.Tao et al. [60] and Chen et al. [61] modeled the normal concrete stress-strain relationship (Figure 4).The stress-strain characteristics of concrete can be expressed as follows: where f c , ε c , f p , ε p , and E 0 are the concrete's corresponding stress (in MPa), strain, maximum experimental stress (in MPa), strain at maximum experimental stress (taken as 0.002), and elastic modulus (in MPa), E 0 = 4700 f c , respectively.
difference in behavior between concrete and ductile materials, like steel, is due to their distinct load responses.Brittle materials tend to develop tensile fracture perpendicular to the direction of the largest tensile strain.Thus, the cracks tend to develop parallel to the maximum compressive stress when a concrete specimen is subjected to a uniaxial com pressive load.In compression, the stress-strain curve of concrete being elastic can gradu ally reach up to the maximum CS followed by a decrease into a softening region and even tually crushing failure at an ultimate strain.Tao et al. [60] and Chen et al. [61] modeled the normal concrete stress-strain relationship (Figure 4).The stress-strain characteristics of concrete can be expressed as follows: where f c , ε c , f p , ε p , and  0 are the concrete's corresponding stress (in MPa), strain maximum experimental stress (in MPa), strain at maximum experimental stress (taken as 0.002), and elastic modulus (in MPa),  0 = 4700 √ f c , respectively.The stress-strain curve of concrete under tension is approximately linearly elastic up to the maximum TS.After this point, the cracks and strength of concrete decrease gradu ally to zero [26,27].The uniaxial TS of concrete can be modeled via the following equation (2 The damaged plasticity model in ABAQUS is characterized by the degradation of the compression and tension.When plasticizers are introduced, the elastic stiffness of the ma terial decreases due to its deterioration, making it unable to restore the original elastic stiffness.This is crucial for cyclic loading, as the two damage parameters are considered The stress-strain curve of concrete under tension is approximately linearly elastic up to the maximum TS.After this point, the cracks and strength of concrete decrease gradually to zero [26,27].The uniaxial TS of concrete can be modeled via the following equation: The damaged plasticity model in ABAQUS is characterized by the degradation of the compression and tension.When plasticizers are introduced, the elastic stiffness of the material decreases due to its deterioration, making it unable to restore the original elastic stiffness.This is crucial for cyclic loading, as the two damage parameters are considered to be dependent on the plastic strains, indicating a reduction in the elastic stiffness.Figure 5 displays the stress-strain curve for tension and compression modeled using the damage plasticity model.to be dependent on the plastic strains, indicating a reduction in the elastic stiffness.Figure 5 displays the stress-strain curve for tension and compression modeled using the damage plasticity model.The element develops cracks and eventually damages the material when it experiences a tension that surpasses its TS.This behavior can be expressed by a variable called the damage parameter (dt).Furthermore, the behavior of the element under CS is represented by the variable damage parameters (dc).Generally, the developed cracks do not affect the stiffness in compression unless the degradation and compression stiffness reach the same level as the stiffness in tension, in which the dc becomes zero.Following the experimental results of Trung [58], this study used (in ABAQUS for the concrete material) a CS value of 33.8 MPa, a Young's modulus of 27,324 MPa, and a Poisson's ratio of 0.2.Five parameters were considered to describe the Drucker-Prager flow potential yield surface that was first developed by Lubliner et al. [62] and then upgraded by Lee et al. [63].To accurately determine the values of the selected parameters, several tests were conducted in the model wherein the default parameters of ABAQUS were used.These parameters included (1) the dilation angle equal to 31 that denoted the volume change to shear strain ratio, (2) the eccentricity equal to 0 when the flow potential approaches a straight line or The element develops cracks and eventually damages the material when it experiences a tension that surpasses its TS.This behavior can be expressed by a variable called the damage parameter (dt).Furthermore, the behavior of the element under CS is represented by the variable damage parameters (dc).Generally, the developed cracks do not affect the stiffness in compression unless the degradation and compression stiffness reach the same level as the stiffness in tension, in which the dc becomes zero.Following the experimental results of Trung [58], this study used (in ABAQUS for the concrete material) a CS value of 33.8 MPa, a Young's modulus of 27,324 MPa, and a Poisson's ratio of 0.2.Five parameters were considered to describe the Drucker-Prager flow potential yield surface that was first developed by Lubliner et al. [62] and then upgraded by Lee et al. [63].To accurately determine the values of the selected parameters, several tests were conducted in the model wherein the default parameters of ABAQUS were used.These parameters included (1) the dilation angle equal to 31 that denoted the volume change to shear strain ratio, (2) the eccentricity equal to 0 when the flow potential approaches a straight line or 0.1, (3) the initial uniaxial to biaxial CS ratio of 1.16 (Figure 6), ( 4) the viscosity that represented the relaxation time of the visco-plastic system that was assumed to be zero because the BCJ model caused a severe convergence difficulty, enabling us to improve the convergence of the model in the softening region, and ( 5) the yield surface on the deviatoric plane defined by the ratio of the second stress invariant on the tensile to the compressive meridian was equal to 2/3.
Eng 2024, 5, FOR PEER REVIEW 8 0.1, (3) the initial uniaxial to biaxial CS ratio of 1.16 (Figure 6), ( 4) the viscosity that represented the relaxation time of the visco-plastic system that was assumed to be zero because the BCJ model caused a severe convergence difficulty, enabling us to improve the convergence of the model in the softening region, and ( 5) the yield surface on the deviatoric plane defined by the ratio of the second stress invariant on the tensile to the compressive meridian was equal to 2/3.The linear isotropic and bilinear kinematic model was applied to simulate the elastic and inelastic behavior of the steel components in the reinforcement bar.The material model in ABAQUS used the Young's modulus value of 200,000 MPa and Poisson's ratio of 0.3 to accurately represent the behavior of steel.Following the tensile test conducted by Trung [58], the steel material parameters utilized in the finite element models were derived.These properties were then incorporated into the model (Table 2 and Figure 7).An orthotropic elastic material was used to represent CFRPCs.Following the report of Narayanan [65], the material properties were determined.The required material properties for the 3D analysis of CFRP are displayed in Table 1.The model assumed full bonding between CFRP and concrete.The linear isotropic and bilinear kinematic model was applied to simulate the elastic and inelastic behavior of the steel components in the reinforcement bar.The material model in ABAQUS used the Young's modulus value of 200,000 MPa and Poisson's ratio of 0.3 to accurately represent the behavior of steel.Following the tensile test conducted by Trung [58], the steel material parameters utilized in the finite element models were derived.These properties were then incorporated into the model (Table 2 and Figure 7).An orthotropic elastic material was used to represent CFRPCs.Following the report of Narayanan [65], the material properties were determined.The required material properties for the 3D analysis of CFRP are displayed in Table 1.The model assumed full bonding between CFRP and concrete.

Geometry
As mentioned previously, to accurately represent the geometrical configuration and dimensions of the test specimens', 3D numerical models were developed.The mesh for the beam and columns was chosen so that the solid elements' node points could line up

Geometry
As mentioned previously, to accurately represent the geometrical configuration and dimensions of the test specimens', 3D numerical models were developed.The mesh for the beam and columns was chosen so that the solid elements' node points could line up with the real positions of the reinforcement.The mesh was subdivided to obtain more node points, resulting in a reasonable mesh density in the joint areas with the acceptable aspect ratio of the elements.The two materials shared the same nodes because the steel reinforcing truss element was attached between the nodes of each concrete element to create a perfect bond.It is worth noting that the perfect bond assumption used in the structural model did not cause a significant error in the predicted load-displacement response.In the CFRPC, the full contact was simulated at the interface between the CFRPC shell and the concrete elements.Values for the cross-sectional area and initial strain of reinforcement as well as the CFRP were entered.A value of zero was entered for the initial strain due to the absence of initial stress on the reinforcement.Figure 8 illustrates the finite element mesh and geometry of the joints.Four joints were modeled including non-seismic detail, seismic detail, CFRPC to strengthen the NSCBCJs, and two layers of CFRPC specimens for strengthening the NSCBCJs.

Geometry
As mentioned previously, to accurately represent the geometrical configuration and dimensions of the test specimens', 3D numerical models were developed.The mesh for the beam and columns was chosen so that the solid elements' node points could line up with the real positions of the reinforcement.The mesh was subdivided to obtain more node points, resulting in a reasonable mesh density in the joint areas with the acceptable aspect ratio of the elements.The two materials shared the same nodes because the steel reinforcing truss element was attached between the nodes of each concrete element to create a perfect bond.It is worth noting that the perfect bond assumption used in the structural model did not cause a significant error in the predicted load-displacement response.In the CFRPC, the full contact was simulated at the interface between the CFRPC shell and the concrete elements.Values for the cross-sectional area and initial strain of reinforcement as well as the CFRP were entered.A value of zero was entered for the initial strain due to the absence of initial stress on the reinforcement.Figure 8 illustrates the finite element mesh and geometry of the joints.Four joints were modeled including non-seismic detail, seismic detail, CFRPC to strengthen the NSCBCJs, and two layers of CFRPC specimens for strengthening the NSCBCJs.

Boundary Condition
The loads and limitations applied in the finite element models were identical to the conducted test.In addition, the lateral cyclic loading was carried out at each drift ratio in the finite element models.Herein, U X , U Y , and U Z are the displacement in the direction of x (beam axis), y (column axis), and z (displacement perpendicular to the x and y axis) axis, respectively.The z-direction displacements were completely restrained at the top of the column (Uz = 0) and the lateral displacements were applied in the x-axis direction.To simulate the hinge support, the displacements in the column bottom at each direction were restrained (Ux = Uy = Uz = 0).Additionally, the displacements in the y and z directions were restricted (Uy = Uz = 0) at the end of the beam to indicate roller support.Figure 9 shows the boundary condition used in finite element analysis.

Boundary Condition
The loads and limitations applied in the finite element models were identical to the conducted test.In addition, the lateral cyclic loading was carried out at each drift ratio in the finite element models.Herein, UX, UY, and UZ are the displacement in the direction of x (beam axis), y (column axis), and z (displacement perpendicular to the x and y axis) axis, respectively.The z-direction displacements were completely restrained at the top of the column (Uz = 0) and the lateral displacements were applied in the x-axis direction.To simulate the hinge support, the displacements in the column bottom at each direction were restrained (Ux = Uy = Uz = 0).Additionally, the displacements in the y and z directions were restricted (Uy = Uz = 0) at the end of the beam to indicate roller support.Figure 9 shows the boundary condition used in finite element analysis.

Force versus Displacement Curves
To confirm the viability of the simulated finite element models, the force-displacement curves of four specimens obtained in the current study were compared with the experimental findings of Trung [58] as shown in Figure 10.The finite element results were observed to be consistent with the experimental ones (Table 3).It was asserted that the simulated finite element model is able to predict the inelastic response of the test specimens.The observed slight difference between the finite element model simulations and experimental results may be due to (1) unexpected occurrences in the experimental work, (2) variables, like the concrete's tensile and compressive properties selection, (3) the refinement of meshes, (4) the use of idealized boundary conditions in the finite element models, and (5) nonlinearity in the materials' behavior (materials' properties modeled in ABAQUS).
Instead of comparing the results of every load cycle, the backbone curve should be developed from the experimental and numerical data for each direction of loading with unique behavior.The curve should be plotted in a single quadrant (positive force versus positive displacement or negative force versus negative displacement; consequently, the backbone curve created by connecting the first cycle peak points with the related displacement is shown in Figure 11.The results indicated that the backbone curve (envelope curve) obtained from the finite element analysis was in good agreement with the experimental results, indicating a considerable correlation at all stages of the lateral cycling loading.

. Force versus Displacement Curves
To confirm the viability of the simulated finite element models, the force-displacement curves of four specimens obtained in the current study were compared with the experimental findings of Trung [58] as shown in Figure 10.The finite element results were observed to be consistent with the experimental ones (Table 3).It was asserted that the simulated finite element model is able to predict the inelastic response of the test specimens.The observed slight difference between the finite element model simulations and experimental results may be due to (1) unexpected occurrences in the experimental work, (2) variables, like the concrete's tensile and compressive properties selection, (3) the refinement of meshes, (4) the use of idealized boundary conditions in the finite element models, and (5) nonlinearity in the materials' behavior (materials' properties modeled in ABAQUS).
Instead of comparing the results of every load cycle, the backbone curve should be developed from the experimental and numerical data for each direction of loading with unique behavior.The curve should be plotted in a single quadrant (positive force versus positive displacement or negative force versus negative displacement; consequently, the backbone curve created by connecting the first cycle peak points with the related displacement is shown in Figure 11.The results indicated that the backbone curve (envelope curve) obtained from the finite element analysis was in good agreement with the experimental results, indicating a considerable correlation at all stages of the lateral cycling loading.
The results of the finite element studies showed that a total of 30 cycles at 12 different drift levels varied from 0 to 10% applied to the specimens.The compression axial load on the column was not considered based on the experimental work.The total deformations that occurred near the top end of the columns were measured with a Linear Variable Differential Transformer as mentioned by Trung [58].The hysteretic curves were initially linear and elastic.With the increase in displacement, the specimens were converted into an inelastic response due to the occurrence of pinch effects.The experimental and finite element results proved that the NSCBJ specimens were strengthened due to the integration of CFRPCs, displaying a good seismic behavior with endurance strength at the test end.In short, the specimens showed maximum lateral strength and no crushing of the concrete occurred.In addition, the observed difference in the main reinforcement bar that caused the force-displacement curves was different in the positive and negative direction of loading therefore it differed in the lateral resistance of the peaks.
The ideal bilinear drawing demonstrated by Paulay and Priestley and described in Section 3.4 was used to obtain the first yield and ultimate strengths, respectively, from the experimental ( Q ye ) and finite element models ( Q yFEM ).The ultimate strengths from the finite element and experimental model are Q ue and Q uFEM , respectively.The experimental and numerical results showed that the capacity of the seismic design specimen is improved by 32% compared to the non-seismic specimen.Also, the strengthened specimen with one layer of CFRPC showed a maximum ultimate strength of specimen greater than the non-seismic one because of the rise in the strength capacity of the BCJ by 28%.Finally, the specimen strengthened with two layers of CFRPC revealed better performance when the ultimate load reached to 11.2 kN, with an increase in ultimate load capacity by 41% compared to the NS specimen.
Eng 2024, 5, FOR PEER REVIEW 11 The results of the finite element studies showed that a total of 30 cycles at 12 different drift levels varied from 0 to 10% applied to the specimens.The compression axial load on the column was not considered based on the experimental work.The total deformations that occurred near the top end of the columns were measured with a Linear Variable Differential Transformer as mentioned by Trung [58].The hysteretic curves were initially linear and elastic.With the increase in displacement, the specimens were converted into an inelastic response due to the occurrence of pinch effects.The experimental and finite element results proved that the NSCBJ specimens were strengthened due to the integration of CFRPCs, displaying a good seismic behavior with endurance strength at the test end.In short, the specimens showed maximum lateral strength and no crushing of the concrete occurred.In addition, the observed difference in the main reinforcement bar that caused the force-displacement curves was different in the positive and negative direction of loading therefore it differed in the lateral resistance of the peaks.
The ideal bilinear drawing demonstrated by Paulay and Priestley and described in Section 3.4 was used to obtain the first yield and ultimate strengths, respectively, from the experimental (  ) and finite element models (  ).The ultimate strengths from the finite element and experimental model are   and   , respectively.The experimental and numerical results showed that the capacity of the seismic design specimen is improved by 32% compared to the non-seismic specimen.Also, the strengthened specimen with one layer of CFRPC showed a maximum ultimate strength of specimen greater than the non-seismic one because of the rise in the strength capacity of the BCJ by 28%.Finally, the specimen strengthened with two layers of CFRPC revealed better performance when the ultimate load reached to 11.2 kN, with an increase in ultimate load capacity by 41% compared to the NS specimen.

Failure Modes
To further verify the finite element models, the failure modes of four joint specimens simulated by ABAQUS were validated by the experimental work (Figures 12-15).The model simulation results of the failure modes were in good agreement with the experimental findings, indicating the precision of the established model calculation.According to the experimental results of Trung [58], the failure mode of the NS specimen was the joint shear failure with the beam and column bending cracks.Shear failure of the joint happened because of the transverse reinforcement deficiency in the joint.The specimen designed according to the seismic details showed the beam flexural cracks near the joint.Compared to the non-seismic specimen, the seismic specimen developed more surface cracks in the beam generated from the bending.The shear cracks at the joint for the specimen made with the seismic detail were fewer, due to sufficient joint transverse reinforcement.Thus, it can be concluded that the transverse reinforcement restricted the joint shear failure in the beam.Rao [66] acknowledged that the failure mode in the specimen strengthened by CFRPC is mainly due to the following factors: (1) concrete crushing in compression before the yield of the longitudinal steel bar, (2) the steel yield in tension with the debonding of the CFRPC sheet, (3) the steel yield in tension with the concrete crushing, and (4) the spalling of the concrete cover.
The single-layer CFRPC-integrated specimen showed the improved strength of the beam, column, and joint.In addition, experimental results of the CFRPC-strengthened specimen exhibited fewer cracks because the restraint provided by the two strips [Ushaped] reduced the T-shaped sheets' debonding in the column.Finally, the presence of CFRPC could debond at a large lateral displacement, classifying the failure mode as steel yield (Figure 14c) in tension with debonding of the CFRPC sheet.The behavior of the

Failure Modes
To further verify the finite element models, the failure modes of four joint specimens simulated by ABAQUS were validated by the experimental work (Figures 12-15).The model simulation results of the failure modes were in good agreement with the experimental findings, indicating the precision of the established model calculation.According to the experimental results of Trung [58], the failure mode of the NS specimen was the joint shear failure with the beam and column bending cracks.Shear failure of the joint happened because of the transverse reinforcement deficiency in the joint.The specimen designed according to the seismic details showed the beam flexural cracks near the joint.Compared to the non-seismic specimen, the seismic specimen developed more surface cracks in the beam generated from the bending.The shear cracks at the joint for the specimen made with the seismic detail were fewer, due to sufficient joint transverse reinforcement.Thus, it can be concluded that the transverse reinforcement restricted the joint shear failure in the beam.Rao [66] acknowledged that the failure mode in the specimen strengthened by CFRPC is mainly due to the following factors: (1) concrete crushing in compression before the yield of the longitudinal steel bar, (2) the steel yield in tension with the debonding of the CFRPC sheet, (3) the steel yield in tension with the concrete crushing, and (4) the spalling of the concrete cover.
The single-layer CFRPC-integrated specimen showed the improved strength of the beam, column, and joint.In addition, experimental results of the CFRPC-strengthened specimen exhibited fewer cracks because the restraint provided by the two strips [Ushaped] reduced the T-shaped sheets' debonding in the column.Finally, the presence of CFRPC could debond at a large lateral displacement, classifying the failure mode as steel yield (Figure 14c) in tension with debonding of the CFRPC sheet.The behavior of the specimen made with two layers of CFRPC integration was same as the single-layer specimen, except that there was an improvement in the performance.This was mainly due to the increased confinement of the joint and the reduced shear cracks of the joint, meaning that the specimen attained larger lateral strength compared to the others.To understand the joint performance more precisely, the failure mechanisms of the specimens were assessed while taking the parameters, like concrete tensile and compressive damage, Von-Mises's stress of the reinforcement, and maximum principal plastic strain.In brief, the finite element models simulated the joint nonlinear behavior appropriately.Figure 12a clearly indicated that the NS beam-column developed joint, column, and beam tensile damage.From the joint to the beam, the concrete exhibited tensile damage.The longitudinal bars in the beam were found to be in tension, while yielding of the transverse joint bars occurred.Concrete damage developed in the joint.The NS specimen failed due to extensive cracking of the concrete at the joint as its ultimate strength was attained.Figure 12b shows the concrete damage under compression.The concrete damage that originated from the repeated loading was cumulative, leading to the strength and stiffness degradation.The principal plastic strain value served as an indicator for the occurrence of cracks in the concrete damage model.The cracks formed for the positive value of the principal plastic strain.Figure 12d shows that the joint failed in the concrete diagonal crack region, which refers to the same crack pattern obtained by the finite element model (Figure 12e).
Figure 13a revealed that the SD beam-column suffered from beam tensile damage wherein the concrete's tensile damage mostly spread in the beam.The yielding of the longitudinal bars occurred without yielding of the transverse joint bars (Figure 13c).In addition, the concrete's damage developed in the joint.The failure of this specimen occurred due to cracking of the concrete beam.The principal plastic strain showed the formation of cracks in the beam which is very similar to the crack pattern obtained by the finite element model (Figure 14d).
Figures 14 and 15 elucidate the finite element output and comparison between the failure mode at ultimate strength for the specimens strengthened by one and two layers of CFRPC.The stress contribution of the CFRPC sheet to the concrete was observed, damaging the concrete under tension.The compression value of the CFRPC specimen was lower compared to the NS specimen.This finding indicated that the cracks in the CFRPCintegrated specimen were lower than the NS specimen, as mentioned above.In addition, the maximum plastic strain distribution for the specimen was mainly concentrated at the beam end and joint, indicating the appearance of same crack pattern as observed in the experiment.
Eng 2024, 5, FOR PEER REVIEW 14 specimen made with two layers of CFRPC integration was same as the single-layer specimen, except that there was an improvement in the performance.This was mainly due to the increased confinement of the joint and the reduced shear cracks of the joint, meaning that the specimen attained larger lateral strength compared to the others.To understand the joint performance more precisely, the failure mechanisms of the specimens were assessed while taking the parameters, like concrete tensile and compressive damage, Von-Mises's stress of the reinforcement, and maximum principal plastic strain.In brief, the finite element models simulated the joint nonlinear behavior appropriately.Figure 12a clearly indicated that the NS beam-column developed joint, column, and beam tensile damage.From the joint to the beam, the concrete exhibited tensile damage.The longitudinal bars in the beam were found to be in tension, while yielding of the transverse joint bars occurred.Concrete damage developed in the joint.The NS specimen failed due to extensive cracking of the concrete at the joint as its ultimate strength was attained.Figure 12b shows the concrete damage under compression.The concrete damage that originated from the repeated loading was cumulative, leading to the strength and stiffness degradation.The principal plastic strain value served as an indicator for the occurrence of cracks in the concrete damage model.The cracks formed for the positive value of the principal plastic strain.Figure 12d shows that the joint failed in the concrete diagonal crack region, which refers to the same crack pattern obtained by the finite element model (Figure 12e).Figure 13a revealed that the SD beam-column suffered from beam tensile damage wherein the concrete's tensile damage mostly spread in the beam.The yielding of the longitudinal bars occurred without yielding of the transverse joint bars (Figure 13c).In addition, the concrete's damage developed in the joint.The failure of this specimen occurred due to cracking of the concrete beam.The principal plastic strain showed the formation of cracks in the beam which is very similar to the crack pattern obtained by the finite element model (Figure 14d).
Figures 14 and 15 elucidate the finite element output and comparison between the failure mode at ultimate strength for the specimens strengthened by one and two layers of CFRPC.The stress contribution of the CFRPC sheet to the concrete was observed, damaging the concrete under tension.The compression value of the CFRPC specimen was lower compared to the NS specimen.This finding indicated that the cracks in the CFRPCintegrated specimen were lower than the NS specimen, as mentioned above.In addition, the maximum plastic strain distribution for the specimen was mainly concentrated at the beam end and joint, indicating the appearance of same crack pattern as observed in the experiment.

Requirement of ACI 374.1-05
To achieve a joint that can be used in the high seismic risk area, the test results obtained by Trung were checked with the ACI 374.1-05.According to this standard, three criteria including energy dissipation, strength, and stiffness must be checked.

Relative Energy Dissipation
To ensure that the joints have a good damping after seismic action, relative energy dissipation at 3.5% of the drift ratio must exceed 0.125.If the beam-column joint fails to meet this condition, the structure will experience excessive oscillation after a seismic event.The dissipation energy in the cycle represents the elastic-plastic behavior.ACI 374.1-05 stipulates that the energy is calculated from the hysteric curve area at drift ratio of 3.5% (area parallelogram).The parallelogram was created by two horizontal lines and four parallel lines, with first positive stiffness on the upper parallel lines and initial negative stiffness on the bottom parallel lines.Figure 16 shows the calculated relative energy dissipation obtained from the definition of the parallelograms.Table 4 displays the specimens' relative energy dissipation ratios.

Requirement of ACI 374.1-05
To achieve a joint that can be used in the high seismic risk area, the test results obtained by Trung were checked with the ACI 374.1-05.According to this standard, three criteria including energy dissipation, strength, and stiffness must be checked.

Relative Energy Dissipation
To ensure that the joints have a good damping after seismic action, relative energy dissipation at 3.5% of the drift ratio must exceed 0.125.If the beam-column joint fails to meet this condition, the structure will experience excessive oscillation after a seismic event.The dissipation energy in the cycle represents the elastic-plastic behavior.ACI 374.1-05 stipulates that the energy is calculated from the hysteric curve area at drift ratio of 3.5% (area parallelogram).The parallelogram was created by two horizontal lines and four parallel lines, with first positive stiffness on the upper parallel lines and initial negative stiffness on the bottom parallel lines.Figure 16 shows the calculated relative energy dissipation obtained from the definition of the parallelograms.Table 4 displays the specimens' relative energy dissipation ratios.The obtained relative energy ratio values in the table are above 0.125, which clearly shows that the specimens could meet the requirement of the standard.This can be assigned to the design of the specimens that follow the criteria of a weak beam/strong column.Consequently, plastic hinges can occur at the end of beam near joint, making the energy dissipation in plastic hinge region.In addition, transverse joint reinforcement was observed to reduce the joint shear distortion and joint cracks, resulting in higher performance of the specimens made using seismic details compared to the non-seismic specimens.The better energy dissipation obtained in the specimen strengthened by two layers of CFRPC can be ascribed to the increase in the confinement of the joint by double layers of CFRPC integration.

Strength Criteria
In the standard mandates, the test specimen model must adhere to the principle of a strong column and weak beam.The maximum lateral load of the model specimen (E max ) must be lower than λE n (where E n represents the nominal strength of the specimen calculated using a strength reduction factor of 1 and where λ is the percentage ratio of column strength to beam strength).Trung's test result revealed that the lateral load varies in the positive and negative loading directions, thus necessitating separate checks for E max in each direction.Table 5 shows the evaluation of the strength criteria.The specimen meets this requirement since the value of E max was below λE n .The failure mode of NS is the flexural crack of the beam followed by shear failure in the joint.For the CFTPCstrengthened specimen, the crack was smaller in the retrofitted area.Thus, the obtained values of λ indicated the occurrence of beam flexural failure in the specimen.The joints must meet the requirement of the strengths degradation as specified by the standard.The strength degradation of the joints is assessed at a drift level of 3.5% (Table 6) and should exceed 0.75E max for the same loading direction.The results in Table 6 showed that the NSCBCJs failed to satisfy the strength criteria.Conversely, in the specimens strengthened by CFRPC and seismic details, the beam-column joint met the strength criteria.The deficiency in transverse reinforcement in the joint led to a high loss of lateral resistance capacity at a later stage.In the NS specimen, the value strength dropped dramatically by about 58% at a drift ratio of 3.5% due to the cracking of the concrete joint.The reason for an abrupt reduction in the strength value might be due to the bond failure of the longitudinal reinforcement.On the contrary, the CFRPC-strengthened specimens achieved excellent strength values due to the increased confinement provided by the CRFPC to the joint, thereby leading to a reduction in the joint cracks and bond failure.

Stiffness Criteria
After the seismic action, damage occurred in the frame structure during a large lateral displacement under a small lateral load.To overcome this effect, the standard requires that the stiffness of the joints be high (0.05) compared to the initial stiffness.The standard requires that the criterion be checked at a drift range of −0.35% to 0.35% using the secant stiffness, and it must meet the prescribed limitation.The initial positive and negative stiffnesses were calculated by dividing the load per associated displacement.Table 7 shows the results for the stiffness criterion validation of the joints, wherein the specimens were observed to satisfy this criterion.Comparing the seismic and strengthened specimens against the specimen made with non-seismic details against the requirement of ACI 374.1-05 for stiffness, it was found that the non-seismic joint had lower initial and secant stiffness due to the lack of transverse reinforcement in the joint.The specimen with seismic details demonstrated a good performance, since this specimen was designed by providing the reinforcement according to the requirement of the ACI 318-19, wherein the transverse reinforcement reduces the cracks in the joints area leading to an increase in the stiffness.This improvement was attributed to the presence of CFRPC, which enabled us to increase the strength of the beam, column, and joint.Consequently, the stiffness was increased.In short, the CFRMC-strengthened specimens demonstrated improvements in their initial and secant stiffness values.The main remarks obtained from the data analysis can be described as follows: (1) the CFRPC-strengthened specimens showed an increase in lateral resisting force without any spalling of concrete; (2) compared to the non-seismic specimens, both CFRPC-strengthened specimens displayed a good amount of energy dissipation, which was increased from 40 to 57%; (3) at the end of loading, the secant stiffness of the CFRPC-strengthened specimens was higher than that of the non-seismic specimen, indicating an increase in their initial stiffness; (4) a comparative evaluation of the acceptance criteria of ACI 374-05 showed that the specimens with seismic details achieved superior seismic performances.

Effective Elastic Stiffness and Displacement Capacity
The backbone curve of all specimens was developed depending on the data obtained from the experiential study by Trung.The results were verified using the finite element model calculation for each loading direction, revealing a unique behavior.The curve was plotted in a single quadrant (positive force with associated positive displacement, and negative force with associated negative displacement) and the backbone curve was created by connecting the first cycle peak points from each drift level (Figure 17).The backbone curve provided the information related to the peak load (P y ) in each loading direction, effective elastic stiffness (K e ), yield displacement (∆ y ), ultimate displacement (∆ u ), and displacement ductility factor (µ).The determination of the yield displacement was challenging due to the absence of a clear yield displacement in the load-displacement curve of a joint.To begin with, the nonlinear capacity curve was simplified by representing it as a bilinear curve that included a specified yield displacement.The values of ∆ y for various specimens were calculated following the recommendation of Paulay and Priestley [67]. Figure 18 illustrates an idealized bilinear curve for the joints, which is then utilized to derive the ductility parameters.The values of K e were determined using the tangent line of the secant slope from the original point and the point when the line had 0.75 of the yield load, where ∆ yi is the displacement corresponding to 0.75 P y in the ascending part of the tangent line.The following relation must satisfy the intersection of the horizontal by trial and error of the line drawn in the failure part to the line drawn from the original by slope provided the ∆ y .
A1 + A2 = A3 + A4 Eng 2024, 5, FOR PEER REVIEW 21 negative force with associated negative displacement) and the backbone curve was created by connecting the first cycle peak points from each drift level (Figure 17).The backbone curve provided the information related to the peak load (P y ) in each loading direction, effective elastic stiffness (Ke), yield displacement (Δ y ), ultimate displacement (Δ u ), and displacement ductility factor (µ ).The determination of the yield displacement was challenging due to the absence of a clear yield displacement in the load-displacement curve of a joint.To begin with, the nonlinear capacity curve was simplified by representing it as a bilinear curve that included a specified yield displacement.The values of (Δ y ) for various specimens were calculated following the recommendation of Paulay and Priestley [67]. Figure 18 illustrates an idealized bilinear curve for the joints, which is then utilized to derive the ductility parameters.The values of Ke were determined using the tangent line of the secant slope from the original point and the point when the line had 0.75 of the yield load, where Δ yi is the displacement corresponding to 0.75 P y in the ascending part of the tangent line.The following relation must satisfy the intersection of the horizontal by trial and error of the line drawn in the failure part to the line drawn from the original by slope provided the Δ y .
The specimen's displacement ductility factor (µ ) was calculated by dividing its ultimate displacement by the yield displacement via the following equation:  ture, leading to an enhancement in the strength performance of the specimen.The specimen containing two layers of CFRPC showed higher effective stiffness than the specimen designed with one layer of CFRPC.The effective stiffness of the specimen designed with seismic details was higher than that of the non-seismic specimen.This strength enhancement can be ascribed to the seismic reinforcement details of the specimen.The specimen's displacement ductility factor (µ) was calculated by dividing its ultimate displacement by the yield displacement via the following equation: The bilinear diagram was used to represent the backbone curve.The effective stiffness ( K e ) and the deformation capacity (∆ y , ∆ u , and µ) in the negative and positive loading directions of the specimens are enlisted in Table 7.The results clearly indicated that the ductility factor ranged from 3.10 to 6.16, which was larger than the suggested ductility for reinforced concrete structures (µ ≥ 2) as obtained by Peng, Zhong, et al. [68].It can be asserted that the ductility of the specimens was higher than the suggested ductility for the reinforced concrete type of structures.Also, the CFRPC-strengthened specimens and the specimen with seismic details had higher ductility factors compared to the non-seismic specimen.Table 8 shows the values for the effective stiffness (K e ) and the displacement capacity (∆ y , ∆ u ) in both the negative and positive loading directions for the specimen.The effective stiffness of the CFRPC-reinforced specimens was observed to be higher than that of the non-seismic specimen, which was mainly due to the presence of CFRPC in the structure, leading to an enhancement in the strength performance of the specimen.The specimen containing two layers of CFRPC showed higher effective stiffness than the specimen designed with one layer of CFRPC.The effective stiffness of the specimen designed with seismic details was higher than that of the non-seismic specimen.This strength enhancement can be ascribed to the seismic reinforcement details of the specimen.Figure 19 shows the drift ratio-dependent variation in the cumulative dissipated energy of various specimens.The energy dissipation capacity was determined using the enclosed loop area of each specimen in the hysteresis curve.The non-seismic specimen experienced pinching behavior during testing.This observation can be mainly attributed to the occurrence of cracks in the joint, leading to a decrease in the energy dissipation in the specimen.The presence of CFRPC significantly increased the energy dissipation.At the final loading, the energy dissipation capacity of the two-layer CFRPC-integrated specimen was approximately 43% higher than the one obtained for the non-seismic specimen.In addition, the one-layer CFRPC-strengthened specimen showed an increase in energy dissipation of 25% compared to the non-seismic specimen.Generally, the specimens reinforced with CFRPC exhibited a greater ability to dissipate the energy compared to the non-seismic specimen.In brief, the specimens reinforced with CFRPC displayed higher flexural strength due to the presence of CFRPC material at the beam-column junction, producing an enhanced energy dissipation.It is worth noting that the energy dissipation was directly impacted by the presence of transverse reinforcement in the joint.During the final loading stage, the seismic specimen exhibited an energy dissipation capacity approximately 20% more than that of the non-seismic specimen.Essentially, the confinement effect on the core concrete was enhanced by the transverse reinforcement, leading to a reduction in joint cracking.In contrast, the non-seismic specimen showed a smaller energy dissipation compared to the others until the end of the test.Figure 19 shows the drift ratio-dependent variation in the cumulative dissipated energy of various specimens.The energy dissipation capacity was determined using the enclosed loop area of each specimen in the hysteresis curve.The non-seismic specimen experienced pinching behavior during testing.This observation can be mainly attributed to the occurrence of cracks in the joint, leading to a decrease in the energy dissipation in the specimen.The presence of CFRPC significantly increased the energy dissipation.At the final loading, the energy dissipation capacity of the two-layer CFRPC-integrated specimen was approximately 43% higher than the one obtained for the non-seismic specimen.In addition, the one-layer CFRPC-strengthened specimen showed an increase in energy dissipation of 25% compared to the non-seismic specimen.Generally, the specimens reinforced with CFRPC exhibited a greater ability to dissipate the energy compared to the non-seismic specimen.In brief, the specimens reinforced with CFRPC displayed higher flexural strength due to the presence of CFRPC material at the beam-column junction, producing an enhanced energy dissipation.It is worth noting that the energy dissipation was directly impacted by the presence of transverse reinforcement in the joint.During the final loading stage, the seismic specimen exhibited an energy dissipation capacity approximately 20% more than that of the non-seismic specimen.Essentially, the confinement effect on the core concrete was enhanced by the transverse reinforcement, leading to a reduction in joint cracking.In contrast, the non-seismic specimen showed a smaller energy dissipation compared to the others until the end of the test.Figure 20 shows the simplified hysteresis loop that consists of the lateral load versus the displacement.The equivalent damping coefficient (ℎ  ) of the studied specimens was calculated to evaluate the actual energy dissipation ability via the following equation: S ABC +S CDA S OBF +S ODE (5) where   is the area of the curve ABC (Figure 20), wherein similar definitions were used for   ,   , and   .Figure 20 shows the simplified hysteresis loop that consists of the lateral load versus the displacement.The equivalent damping coefficient (h e ) of the studied specimens was calculated to evaluate the actual energy dissipation ability via the following equation: where S ABC is the area of the curve ABC (Figure 20), wherein similar definitions were used for S CDA , S OBF , and S ODE .
Figure 21 shows the drift ratio-dependent variation in the equivalent damping coefficient.The value of (h e ) is about 0.1 for the normal reinforcement concrete joints and 0.3 for the joints with steel beam reinforcement concrete columns joints [68].The obtained results indicated that the joints' capability for energy dissipation was within the range of reinforcement concrete joints and that it was lower than that of steel beam reinforcement concrete column joints.In the present case, the non-seismic specimen showed lower values of (h e ) compared to other specimens, while the two-layer CFRPC-strengthened sample showed the best performance.Figure 21 shows the drift ratio-dependent variation in the equivalent damping coefficient.The value of (ℎ  ) is about 0.1 for the normal reinforcement concrete joints and 0.3 for the joints with steel beam reinforcement concrete columns joints [68].The obtained results indicated that the joints' capability for energy dissipation was within the range of reinforcement concrete joints and that it was lower than that of steel beam reinforcement concrete column joints.In the present case, the non-seismic specimen showed lower values of (ℎ  ) compared to other specimens, while the two-layer CFRPC-strengthened sample showed the best performance.Figure 21 shows the drift ratio-dependent variation in the equivalent damping coefficient.The value of (ℎ  ) is about 0.1 for the normal reinforcement concrete joints and 0.3 for the joints with steel beam reinforcement concrete columns joints [68].The obtained results indicated that the joints' capability for energy dissipation was within the range of reinforcement concrete joints and that it was lower than that of steel beam reinforcement concrete column joints.In the present case, the non-seismic specimen showed lower values of (ℎ  ) compared to other specimens, while the two-layer CFRPC-strengthened sample showed the best performance.

Stiffness Degradation
Figure 22 displays the drift ratio-dependent variation in the secant stiffnesses for all the studied specimens.The degradation of the secant stiffness for the beam-column joint in each loading cycle was computed using the following equation:

Stiffness Degradation
Figure 22 displays the drift ratio-dependent variation in the secant stiffnesses for all the studied specimens.The degradation of the secant stiffness for the beam-column joint in each loading cycle was computed using the following equation: where +F ij , −F ij is the positive (negative) peak load and +∆ ij , −∆ ij is the displacement associated with the positive (or negative) ultimate lateral strength [34].
specimens.Specimen strengthened by two layers of CFRPC showed higher stiffness, since all the specimens had the same value in terms of the beam reinforcement ratio.In essence, the confinement given by the presence of CFRPC to the beam-column joint could mitigate the distortion of the joint panel, hence enhancing its strength and reducing crack formation in the surrounding area.In short, CFRPC-integrated NSCBCJs showed improved performance against seismic events.

Conclusions
Based on the obtained results ,the following conclusions can be drawn: i.The strengthening of NSCBCJs using CFRPC reinforcement enabled us to reduce crack formation in the plastic hinge location, thus leading to a reduced risk of building collapse during seismic events.The plastic hinges for the frame of the reinforced concrete building were placed near the column and beam joint.ii.
The strengthening of the NSCBCJs was appropriate when assessing the behavior of non-seismic details due to the overall lateral displacement increase that produces an increase in ductility demand.iii.
The NSCBCJs being non-ductile meant that they could display large lateral displacement which can be attributed to the significant cracks and deformation of the joint area.Consequently, the observed large lateral displacement can cause severe failure in the frame building during strong seismic vibration.iv.
The proposed NSCBCJs showed significant joint shear failure together with the beam and column flexural failure.However, the seismic joint failure was beam flexural, while CFRPC-integrated joints showed beam flexural failure with the debonding of the CFRPC sheet.Since the studied joints had beam flexural strength, the lateral load capacity difference can be mainly ascribed to the transverse reinforcement and strengthening of the joint.v.
Both seismic and CFRPC-strengthened joints revealed good strength, stiffness, ductility, and energy-dissipation capabilities.The degradation of the secant stiffness values for all specimens increased with the increase in drift ratio (lateral displacement) at various load steps, significantly impacting the specimens' stiffness degradation.The stiffness of the specimens was reduced with the increase in the load step.The non-seismic specimen showed lower stiffness than other specimens.Specimen strengthened by two layers of CFRPC showed higher stiffness, since all the specimens had the same value in terms of the beam reinforcement ratio.In essence, the confinement given by the presence of CFRPC to the beam-column joint could mitigate the distortion of the joint panel, hence enhancing its strength and reducing crack formation in the surrounding area.In short, CFRPC-integrated NSCBCJs showed improved performance against seismic events.

Conclusions
Based on the obtained results, the following conclusions can be drawn: i The strengthening of NSCBCJs using CFRPC reinforcement enabled us to reduce crack formation in the plastic hinge location, thus leading to a reduced risk of building collapse during seismic events.The plastic hinges for the frame of the reinforced concrete building were placed near the column and beam joint.ii The strengthening of the NSCBCJs was appropriate when assessing the behavior of non-seismic details due to the overall lateral displacement increase that produces an increase in ductility demand.iii The NSCBCJs being non-ductile meant that they could display large lateral displacement which can be attributed to the significant cracks and deformation of the joint area.Consequently, the observed large lateral displacement can cause severe failure in the frame building during strong seismic vibration.iv The proposed NSCBCJs showed significant joint shear failure together with the beam and column flexural failure.However, the seismic joint failure was beam flexural, while CFRPC-integrated joints showed beam flexural failure with the debonding of the CFRPC sheet.Since the studied joints had beam flexural strength, the lateral load capacity difference can be mainly ascribed to the transverse reinforcement and strengthening of the joint.v Both seismic and CFRPC-strengthened joints revealed good strength, stiffness, ductility, and energy-dissipation capabilities.vi Both seismic and CFRPC-strengthened joints achieved the required seismic code provisions of ACI 374.1-05.Conversely, the NSCBCJs failed to satisfy the desired requirement due to a failure in achieving the strength criteria.vii The finite element model simulation results for the stress and strain response demonstrated the considerable impact of CFRPC integration into NSCBCJs, leading to a reasonable strength capacity enhancement in the joints, thereby reducing the joint cracks.
viii It is established that the present comprehensive numerical modeling and simulation study may contribute to the development of CFRPC-reinforced high-performance anti-seismic CBCJs, thus contributing towards the sustainable construction sector.

Figure 3 .
Figure 3. Elements used in ABAQUS to model: (a) C3D8R for concrete, (b) T3D2 for longitudinal and transverse reinforcement, and (c) the S4R5 four-node shell for FRP.

Figure 5 .
Figure 5. Stiffness recovery curve of concrete under tension and compression.

Figure 5 .
Figure 5. Stiffness recovery curve of concrete under tension and compression.

Figure 8 .
Figure 8. 3D finite element model for the analysis of NSCBCJ-and CFRPC-integrated NSCBCJs.(a) NSCBCJ specimen, (b) seismic specimen, (c) one and two layers of the CFRPC specimen, (d) longitudinal and transverse reinforcement of the NSCBCJ specimen, (e) longitudinal and transverse reinforcement of the seismic specimen, and (f) mesh of all models.

Figure 8 .
Figure 8. 3D finite element model for the analysis of NSCBCJ-and CFRPC-integrated NSCBCJs.(a) NSCBCJ specimen, (b) seismic specimen, (c) one and two layers of the CFRPC specimen, (d) longitudinal and transverse reinforcement of the NSCBCJ specimen, (e) longitudinal and transverse reinforcement of the seismic specimen, and (f) mesh of all models.

Figure 9 .
Figure 9.The boundary conditions used in the model simulation.

Figure 9 .
Figure 9.The boundary conditions used in the model simulation.

4 .
Results and Discussion 4.1.Verification of the Finite Element Models with Experiments 4.1.1

Figure 10 .
Figure 10.Comparison force-displacement curves for the non-seismic specimen (NS), seismic specimen (SD), strengthened one-layer CFRPC specimen, and strengthened two-layer CFRPC specimen.(a) Experimental results of Trung [58] on NSCBCJ; (b) comparison of the proposed model results with experiments with NSCBCJ; (c) experimental results of Trung [58] with SD; (d) comparison of the proposed model results with experiments with SD; (e) experimental results of Trung [58] with one-layer CFRPC; (f) comparison of the proposed model results with experiments with one-layer CFRPC; (g) experimental results of Trung [58] on two-layer CFRPC; (h) comparison of the proposed model results with experiments of two-layer CFRPC.

Figure 18 .
Figure 18.Ideal bilinear load-displacement behavior of the specimen.

Figure 18 .
Figure 18.Ideal bilinear load-displacement behavior of the specimen.

Figure 19 .
Figure 19.Drift ratio-dependent variation in the cumulative dissipated energy.

Figure 19 .
Figure 19.Drift ratio-dependent variation in the cumulative dissipated energy.

Figure 21 .
Figure 21.Drift ratio-dependent variation in the equivalent damping coefficient.

4. 3 .
Figure22displays the drift ratio-dependent variation in the secant stiffnesses for all the studied specimens.The degradation of the secant stiffness for the beam-column joint

Figure 21 .
Figure 21.Drift ratio-dependent variation in the equivalent damping coefficient.

Figure 21 .
Figure 21.Drift ratio-dependent variation in the equivalent damping coefficient.

Table 2 .
Model properties of reinforcement material.

Table 2 .
Model properties of reinforcement material.

Table 3 .
Comparison between experimental and finite element results.

Table 3 .
Comparison between experimental and finite element results.

Table 4 .
Validation of relative energy dissipation ratios of various specimens.

Table 5 .
Validation of maximum lateral load resistance for various specimens.

Table 6 .
Assessment of joint strength degradation at a drift level of 3.5%.
* Less than the ratio of the strength.

Table 7 .
Validation of the secant stiffness of various specimens.

Table 8 .
Effective stiffness and ductility of various specimens.

Table 8 .
Effective stiffness and ductility of various specimens.