Restoring Force Model of an Energy-Dissipation Joint in Hybrid Frames: Simplified Skeleton Curve and Hysteretic Rules

In this paper, a restoring force model, composed of a trilinear skeleton curve and hysteretic rules, is proposed based on nine pseudostatic tests of the energy-dissipation joint under horizontal low cyclic loading. *e critical points of the simplified skeleton curve are obtained via theoretical derivation and FE simulation. *e hysteretic rules for the joints are simplified as a concave hexagon, where the parameters of the critical points are optimized by the genetic algorithm (GA). Using the established trilinear skeleton curve, three different working stages, i.e., elastic, hardening, and softening, were divided by the critical points and the moment stiffness of three stages can be calculated. *e proposed hysteretic rules of each stage can reveal and explain the “pinching” in the cyclic loading, which make it easier to understand the mechanism of the energy-dissipation joint. *e comparison between the restoring force model and the tests shows that the simplified skeleton curves, the established hysteretic rules, and the ductility and the damping ratio are consistent with the experimental results. Finally, the effectiveness of the established restoring force model is verified.


Introduction
In the past few decades, engineered bamboo products (EBPs) have emerged as alternatives to traditional building materials because of their advantages, such as energy saving and environmental friendliness [1,2]. As the manufacturing technology is getting increasingly mature, the promotion of application of wood and bamboo materials in mid-and high-rise buildings has caught great attention of researchers worldwide [3][4][5][6][7][8].
However, compared with concrete and steel, there are some inherent limitations of EBPs when used as perpendicular bearing components, i.e., the accumulated creep and requirement for a larger size. To overcome the issues mentioned above, a hybrid frame made up of steel columns and engineered bamboo beams is proposed as an effective solution [9,10]. Obviously, the connection between the steel column and engineered bamboo beam becomes an essential aspect of steel-engineered bamboo frame.
Studies on variable parameters have been conducted to explore the appropriate design methods of the joint as a reliable connection. e dowel and bolt are the most commonly used joints in modern wood or EBP buildings, whose ductility is achieved by means of the plastic deformation in metallic connectors [11][12][13][14][15]. Conversely, in heavy portal frames, the energy-dissipation capacity of the connections turns out to be relatively low because of the limited number of joints. To improve the hysteretic performance of timber or EBP connections, some new joints or reinforcement methods were proposed. An innovative posttensioned beam-column timber joint, which turned out to be of good structural performances, has been developed and studied by many researchers [16,17]. Other reinforcement methods such as gluing wood-based panels on each side of the shear planes, by orienting unidirectional fibers perpendicular to the grain and inserting a box-type steel bracket between timber beam and column [10,[18][19][20][21][22][23][24][25][26][27], were tested, and satisfactory results were obtained.
Among these studies, Huang et al. [28] proposed a semirigid joint made up of a steel hinge and a pair of top and seat brackets to carry the shear force and moment from the end of the beam, as shown in Figure 1. In this setup, the joint can serve as a frame connection and an energy-dissipation device. Low cyclic loading experiments for nine joints were conducted in Nanjing Forest University, and the results revealed that the joint provided sufficient strength and stiffness to satisfy the requirements of serviceability, demonstrating a satisfactory energy-dissipation capacity. With a proper design, the ductility ratio and damping ratio of the connections can reach more than 3.0 and 30%, respectively. Huang et al. [28] also suggested a formula to predict the load-carrying capacity of the proposed semirigid joint together with a design method to fully display its load-carrying capacity and ductility. e research results and theoretical derivation verified the pronounced advantages of the joint in connecting the steel-engineered bamboo hybrid frame, which requires further study to investigate its restoring characteristics.
e cyclic restoring force model is an useful tool for performing seismic analysis of energy-dissipation joints in mid-and high-rise buildings [29,30]. Based on the hysteresis loops and skeleton curves of energy-dissipation joints, a simplified trilinear skeleton curve is constructed, where the characteristic values can be obtained via the load-carrying capacity formula and theoretic derivation. e hysteretic rules of a single loop have a reversed S-shape, which can be modeled by a concave hexagon. In this paper, the characteristic values of the concave hexagon are optimized by the genetic algorithm.

Brief Introduction of the Test
e schematic picture of the energy-dissipation joint is illustrated in Figure 1. e base panel and the stub were connected through a 20 mm dimeter bolt. e PSB panels were mounted on the two sides of steel stub of the connection, respectively, through 12 bolts of 22 mm in dimeter. e middle parts of the two energy-dissipation plates (EDPs in the following) were welded on the stub symmetrically, while one end was weld on the base panel.
e details of the test connections (i.e., the diameter and spacing of the bolts, the sizes of the base panel, and stub) were designed with a suitable over strength. In this way, the bolt connection is stronger than the EDPs and the tension and compression forces lead to the yielding or buckling of the EDPs that dissipated energy when exerted earthquake or reversed loading. eoretically, the dimension of the brackets is one of the major factors impacting on the mechanical behavior of the connection.
According to the purpose of the energy-dissipation joint design [28], the width-to-thickness ratio of the energydissipation plate should be greater than 1/5 to ensure that there is only minimal shear under seismic loading. Presimulations in ABAQUS with variable widths, thicknesses, and cross sections were carried out to study the effect of the buckling stress of the plate, which showed that the length-tothickness ratio was a major and reliable factor to explore the hysteretic behaviors of energy-dissipation plate (see Section 4.1.2). erefore, a group of tests were performed to explore the effect of the length-to-thickness ratio on the momentcarrying capacity of the energy-dissipation joint by changing thicknesses with identical cross section and length. e parameters of the tests are presented in Table 1, where l is the length of EDP, b is the width of EDP, and t is the thickness of EDP.
All specimens were loaded under a lateral cyclic force without axial pressure. For the convenience of loading, the specimens were rotated 90 degrees compared to their actual service position in the steel-engineered bamboo frame [28]. Using the setup presented in Figure 2(b), the engineered bamboo beam could be idealized as a cantilever beam that is fixed at the steel column and free at the other end. e lengths of the beam and column were 2 m and 1.2 m to avoid significant influence on the connection [31]. Hence, a similar bending moment to that in an actual steel-engineered bamboo hybrid frame was exerted on the specimens. e advantages of the setup also include the convenience for calculation and installment.
Preloading was conducted to check if all the bolt-joints were properly mounted and worked well and then unloading the actuator to zero to reset the acquisition system. Load was controlled by the movement of actuator at the speed of 4 mm/min before the EDP yielding and of 0.2 Δ mm/min after the first yielding of the EDP, where Δ represents the top displacement of first yielding in the EDP. ree cycles were carried out for each loading grade after the first yielding occurred [32,33]. e loading-unloading regime is illustrated in Figure 2(c).

Analysis of the Experimental Results
According to the nine experimental results, the failure modes can be categorized into two different types, as presented in Table 2.
Since the joint is designed to dissipate earthquake energy by the EDPs, the fatigue strength of the weld in failure mode II dominates the load-carrying capacity of the joint and the EDPs are too thick to fully display their energy-dissipation ability. According to design suggestions in [28], the lengthto-thickness ratio should be limited in the range of 9 to 16. Only the first 4 tests in Table 1 meet the requirement. Besides that, there are only two working stages in the skeleton curves of the last 5 tests, due to the experimental setups' limits. Since the analytical model is intended to reveal the whole working stages, only the results of the first four groups in Table 1 were used for evaluating these coefficients.
As summarized from the hysteresis loops and skeleton curves in [28], a typical skeleton curve can be divided into three stages, as presented in Figure 3. e characteristics of the three working stages are included in Table 3.
e initial stiffness and moment-carrying capacity of the joint were also presented by Huang et al. [28] so that the stiffness in the elastic stage (I) and the maximum moment of the joint can be calculated. However, when structures undergo strong lateral forces (i.e., earthquakes and wind), the energy is dissipated by joints when the elastic stage is 2

Advances in Civil Engineering
Energy-dissipation plate (EDP)      exceeded. erefore, the restoring force model including the nonlinear behavior of the joint should be studied. e restoring force characteristics can be represented by a simplified skeleton curve and the hysteretic rules. e method to construct the restoring force model based on the experimental results (i.e., the skeleton curve and hysteresis loops) is explained as following: (1) by simplifying the nonlinear segments (hardening and softening segments) into straight lines with constant slope, the experimental skeleton curve can be transferred as a trilinear curve with six critical points; (2) the hysteretic rules can be simplified as a hexagon concluded from every single loop, where the maximum points can be obtained from the simplified skeleton curve and the other critical points need to be determined with GA.
Using the established trilinear skeleton curve, three different working stages, namely, elastic, hardening, and softening, were identified and defined by the critical points and the moment stiffness relationship of the three stages. e proposed hysteretic rules of each stages reveal and explain the "pinching" in the cyclic loading, which make it easier to understand the mechanism of the energy-dissipation joint.

Restoring Force Model of the Energy-Dissipation Joint
Based on the analysis of the nine pseudostatic tests of the energy-dissipation joint and the extended FE simulation, the critical points of the energy-dissipation joint with different sizes can be calculated. e restoring force model of this type joint is constructed, where the model parameters obtained can work for other connections of the same type.

Simplified Trilinear Skeleton Curve.
A typical skeleton curve of the energy-dissipation joint in the steel-engineered bamboo hybrid frame is presented in Figure 3, which has three segments: a straight line, a hardening segment with a continuously changing stiffness, and a softening stage with a slight decline after the peak load. ere are two critical points after the yielding, i.e., the peak point and the ultimate point. e comparison of the negative and positive values (absolute moments and rotations) of the two critical points is presented in Table 4, which shows that the errors of the two absolute moments and rotations are less than 5%, within the acceptable range. But the size relationship of the two values cannot be determined in designing since the EDP yields in an uncertain side first. To solve the problem and ensure the design redundancy, the minimum moments and rotations of the two absolute values were selected to represent the peak/ultimate moment and rotation in positive and negative directions.
In this paper, a symmetric trilinear simplified skeleton curve is selected to relatively match with the experimental skeleton curve. In Figure 4, points Y, P, and U refer to the yield point, peak point, and ultimate point of the joint, respectively. erefore, six key values of the simplified trilinear skeleton curve should be determined based on the experimental skeleton curves. M y and θ y represent the moment and rotation of Point Y, respectively; M p and θ p are the moment and rotation of Point P, respectively; M u and θ u are the moment and rotation of Point U, respectively. K 1 , K 2 , and K 3 are the moment stiffness values of each segment belonging to the simplified skeleton curve of the energydissipation joint.

Rotations of the Skeleton Curve.
e yield point, which is also called the soft point, refers to the critical state where the load and displacement curves at the end of the beam begin to deviate from the linear variation [34]. Combined with the test phenomenon, the yield point should be considered as the mutation point where the slope of the skeleton curve begins to change and the point where the EDP begins to buckle. e initial yielding is the point when the stress of EDP reaches the yield strength of steel, f y . us, the tensile or compressive force in EDP must be F EDP � f y bt when initial yielding takes place; from the deformation of energy-dissipation joint in the elastic stage presented in Figure 5, the initial yield moment can be calculated by (1) e initial moment stiffness of the connection can be calculated by K E � f y bt(h + t)/θ E , where θ E represents the rotation of connection, which can be calculated by equation (2) as follows: where E represents Young's module of the steel; b and t are the width and the thickness of EDP, respectively; l represents the length of EDP; and h stands for the clear distance between the two EDPs. Advances in Civil Engineering e comparison of the experimental and the theoretical elastic stiffness is presented in Table 5, which shows that the errors between them are all within 5%. erefore, K E can be considered as the stiffness of stage I.
With all the rotations obtained from the tests, linear regression algorithm in ORIGIN was selected to describe the relation between each rotation, i.e., yield rotation θ y , peak rotation θ P , and ultimate rotation θ u with the initial rotation θ E , respectively. e result shows that the relationship between them are all constants. us, three coefficients (i.e., α 1 , α 2 , and α 3 ) are introduced to obtain them starting from the initial rotation θ E as follows: Table 6 shows the coefficients of rotation in the simplified trilinear skeleton curve.

Load-Carrying Capacity and Ultimate States of the Joint.
e load-carrying capacity is studied considering the following hypotheses: (1) the shear of the EDPs and the moment of the hinge can be omitted to ensure that the steel hinge and pair of EDPs can completely carry the shear force and moment, respectively; (2) all of the connections between the engineered bamboo beam and the stub, steel column, and base panel are considered as rigid connections; and (3) the friction of the hinge is ignored.
Since the buckling of the EDPs dominates the momentcarrying capacity of the joint, it is necessary to explore the mechanical behavior of EDPs. After getting the maximum reaction force of buckling, the peak moment of the energydissipation joint can be obtained by F EDP (h + t). Finite element method (FEM) in ABAQUS was utilized to obtain numerical simulation of the energy-dissipation plate under axial loading. e plate was divided in two segments, namely, (1) the real energy-dissipation segment (red part) and (2) the weld segment (blue part) in Figure 6, to better understand the mechanism. e boundary conditions were set the same as for the energy-dissipation joint: only the outof-plane ends movement was constrained. e connection between the two segments in ABAQUS was set as "tie" (a kind of constrains), which indicated the rigid connection of them. e displacement/force was exerted on the reference point, which was coupled with the top surface of the EDP.
e established 3D FE model is shown in Figure 6. e element type of C3D8R (8-node linear brick element, reduced integration with hourglass control) in ABAQUS was chosen to simulate the large deformation of the actual dissipation segment for its better precision and computing efficiency, while the other components type C3D20R (20-node linear brick element, reduced integration    Advances in Civil Engineering with hourglass control) were adopted since there is no need for large strain and complex damage evolution. In the tests, the behavior of EDPs turns out to be in the elastoplastic range, so the whole life stage, i.e., elastic stage, hardening stage, and softening stage, should be included in the strainstress relation. Figure 7 shows the stress and strain relation of steel. Formulas (4) to (8) show the transformation from nominal to real strain-stress relationship. Besides that, nonnegligible damage of the EDPs was observed in the process of cyclic loading, so the ductile damage criteria was chosen to simulate the evolution of the material damage, as presented in Figure 8. e transformation of the real and nominal strain-stress relationship of the steel can be explained as follows: considering the incompressibility of the plastic deformation, the volume of the material does not change after large deformations, and the following equations can be obtained: where l and l 0 and A and A 0 refer to the initial and final (i.e., after large deformations) lengths and cross sections of the test material. e nominal stress σ nom can be calculated as the force F divided by the initial cross section. Substituting A in the computation of strain, the real stress σ rea can be presented as

Advances in Civil Engineering
where l/l 0 is equal to 1 + ε nom ; in this way, the relationship of nominal and real strain and stress can be described as where ε nom and ε rea represent nominal and real strains, respectively. e transformation of the detailed nominal and real strain and stress of the steel is presented in Table 7.
e stress and deformation of the EDP at the maximum displacement are shown in Figures 9(a) and 9(b). Von Mises stress nephogram is applied to show the stress distribution. e maximum stress or strain occurs mainly in the actual dissipation segment, while the other parts remain in the elastic range. It corresponds well with the design aims and failure mode of the tests. Figure 9(c) shows the displacement-force relationship of the EDP, whereas the loading capacity of the energydissipation joint can be obtained from equation M E � F EDP (h + t) [28]. e result indicates that the relative errors between the tests and simulations of the four energy-dissipation joint are all less than 15%, within the acceptable range.
Considering the fact that the cross section of EDP cannot fully yield when connection failed, a coefficient η � M p / [f u bt(h + t)] is used to quantify the partial yielding effect, where f u is the ultimate strength of the steel. e size of the energy-dissipation plate has the prevalent influence on the moment-carrying capacity and ductility of the energy-dissipation joint. FE simulations were conducted by varying only width (from 50 mm to 100 mm), thickness (from 6 mm to 12 mm), cross section (from 500 mm 2 to 1000 mm 2 ), and length-to-thickness ratio (from 6.6 to 16). e correlations between loading capacity and these parameters were analyzed in SPSS, and the results are presented in Table 8. e results revealed that the Pearson correlations of the four elements are all greater than 0.6, indicating the pronounced relation of these elements with the loading capacity of energy-dissipation plate. Table 8 also shows the reliability of these parameters which equals to 1− ξ, and ξ refers to the maximum variation range when the other parameters changes with one parameter determined. e conclusion that length-to-thickness ratio can serve as a stable and reliable parameter compared with thickness and width can be drawn. In this paper, the length-to-thickness ratio was chosen as an independent variable and η(M p / [f u bt(h + t)]) as the dependent variable to explore the loading capacity of the energy-dissipation joint.
Other 88 FE simulations by varying length-tothickness ratio of the EDP (λ � l/t) were carried out based on the same method above. It turned out that the coefficients η and λ tend to have a negative exponential curve. An empirical relationship between η and λ can be obtained by data fitting which is expressed as follows [28]: e comparison of the peak moments computed by equation (10) with those obtained by testing for groups J-6 and J-8 shows good agreements, with all the relative errors less than 15% [28]. By combining the above hypotheses with finite element simulation performed in ABAQUS, the loadcarrying capacity of the energy-dissipation joint can be estimated by equation (10). In accordance with the Chinese standard GB50017-2017 [32], 0.85M p is taken as the ultimate moment in the simplified skeleton curve, shown as follows:

Method to Obtain the Simplified Skeleton
Curve. e calculation flow chart for constructing the simplified trilinear skeleton curve is presented in Figure 10. e elastic stiffness and rotation can be calculated at the first step by assigning the sizes of the joint and its mechanical properties related to Young's modulus and the yield strength. e yield, peak, and ultimate rotations of the joint can be obtained by three coefficients multiplying the elastic rotation, respectively. Hence, the coordinates of the yield point can be gotten. e other two moments of the critical points of the simplified skeleton curve can be calculated from equations (10) and (11). With all the values of the critical points, the stiffness of the three working stages can be calculated as follows:

Simplified Hysteretic Rules
To examine the theoretic hysteretic rules of the energydissipation joint, a single loop within the whole hysteresis loops of J-8-1 is highlighted, as shown in Figure 11. e shape and feature of the experimental hysteresis loops were analyzed, and some conclusions were reached: (1) the reverse S-shape or Z-shape of the hysteresis loops indicates their pronounced pinching characteristics, (2) the relation of M − θ in the loading process is a straight line, and (3) the reloading and unloading relations of M − θ are approximately straight lines, the slope of which is smaller than that of the loading process. erefore, there is a strength degradation in the reloading and unloading processes. In Figure 12, a concave hexagon made up of six lines is proposed to simulate the reloading, pinching, and unloading processes of the hysteresis loop, where different slopes 8 Advances in Civil Engineering    Advances in Civil Engineering between loading and unloading can simulate its stiffness degradation. e advantages of this type of simplified hysteretic rules are presented in the following: (1) the four movable points in the hexagon indicate that the hexagon can be changed into different shapes; i.e., there is flexibility to simulate hysteresis loops of various shapes; (2) with the critical points properly determined, the hysteretic characteristics of the loop, such as the stiffness of reloading and unloading, and the dissipated energy can be simultaneously simulated without losing precision.
In the concave hexagon, the maximum and minimum points, got from the skeleton curve, refer to the maximum/ minimum reaction forces when the corresponding displacements are exerted. e other four points are related to the pinching and reloading stiffness deterioration of the energy-dissipation joint. Point (x 1 , y 1 ) and Point (x 4 , y 4 ) serve as "breaking points" that indicate the extend of pinching. e combination of Point (x 1 , y 1 ) and Point (x 3 , y 3 ) or Point (x 2 , y 2 ) and Point (x 4 , y 4 ) can be used to calculate the reloading stiffness, which is lower than K 1 because of the deterioration of pinching.

Parametric Optimization.
e genetic algorithm (GA) is a metaheuristic inspired by the natural selection process, which belongs to the larger class of evolutionary algorithms (EA). In recent studies, GAs have proven their pronounced advantages [35,36], i.e., the simplicity, robustness, and insensitivity of initial values for global optimization problems of some large and complex nonlinear systems. ey are commonly used to generate high-quality solutions for optimization and search problems by relying on bio-inspired operators, such as mutations, crossover, and selection. e mechanism of GA, presented in Figures 13 and 14, can be explained as follows: (1) Initial value: random solution values of parameters within specified boundary are chosen in the first step to simulate the first generation of the population. lution vectors are all checked if they meet the tolerance of the fitness function, which is also the optimization objective of the algorithm. e solution vectors with less errors will become the dominant individuals, which will be passed to the next generation. In this way, the solution vector will be optimized to a fitter one. e process of selection will be repeated generation by generation until the best solution vector occurs.
In Matlab, the GA toolbox can be used to achieve the operations above. e construction of fitness function and boundary determination of the parameters are the two aspects of upmost importance.
In this section, the GA is used to identify and optimize the critical parameters of the hysteretic curve, and the fitness function of the genetic algorithm is shown in the following equation: where S and s are the areas of the test and the simplified hysteresis loop, respectively. e range of eight parameters is set as follows: x l < x < x u , where x l and x u represent the upper and lower bounds of the variable. eoretically, the critical values of the concave hexagon are in the range of the maximum and minimum, but some restrictions must be added to ensure the best fitness for the single hysteresis loop. e method to determine the lower  and upper bounds is illustrated in Figure 15 with a single hysteresis loop, where LB represents the lower bound and UB represents the upper bound of the parameter. Two main steps are included in the method: the rough values of each critical points are first found and then expanded to a certain range. Point 1 and Point 3 are the turning points of the hysteresis loop, whose rough values can be gotten directly from the picture. However, the other two points (Point 2 and Point 4) can be calculated as the cross points of the two tangent lines of the reloading and unloading curves. With the rough values of the four critical points determined, the ranges of each parameters can be set as [0.90, 1.10] times the rough values.
Other parameters set in the GA tool are presented as follows: the population size is equal to 1000, the number of generations is equal to 50, the ratio of crossover is set as 0.8, and the function tolerance of function J is set as 10 − 6 . Parameters of a single loop optimized by the GA are shown in Table 9.  Figure 16(b) shows the relative area error between the test and the simplified model with the increase in the evolution algebra of the populations. From Figure 16(a), the overall trend between the test result and the simplified model indicates their agreement. e latter figure shows that the area error between the test and the simplified model is less than 10 − 6 , so the energy dissipation of the hexagon is identical to that of the actual test when the simplified hysteretic rule is used.

Parameters Optimized by the GA.
e same method can be used to optimize the parameters of other hysteresis loops, and the data collection processes are identical to that of Section 5.1. e result reveals that the values (displacement and force) of the critical points increase with the increase in maximum displacement or force and the relationships between the two factors tend to be straight lines or nonlinear curves. To make the relationships more explicit, the ratio of θ to θ p is selected as the independent variable instead of θ.
As mentioned in Section 3, the maximum and minimum moments and rotations at each cycle can be calculated from the skeleton curve. us, eight coefficients (i.e., ξ 1 to ξ 8 ) are used to determine the ratio between x 1 , x 2 , x 3 , x 4 , y 1 , y 2 , y 3 , and y 4 to the maximum and minimum values x max , x min , y max , and y min of each hysteresis loop. Detailed definition of ξ is presented in Table 10, and μ is the ratio of θ to θ p . Linear algorithm in ORIGIN was selected to describe the relationships of μ and ξ, as shown in Figure 17.
In the analysis of the parameters optimized by the GA, the critical points of the hysteretic rules are related to the ratios of measured rotation θ to the peak rotation θ p . In the regression analysis, the four critical points considering the effect of μ (0.5 < μ < 1.5) are presented in the following equations:

Comparison of the Tests and Simplified Model
6.1. Comparison of the Skeleton Curves. Based on the flow chart shown in Figure 10, the simplified skeleton curves of the first four groups in Table 1 are obtained, and Figure 18 presents the comparison of the skeleton curve from the test and simplified skeleton curves. e results show that the two skeleton curves are notably consistent with each other and that the simplified skeleton curve can basically simulate the relation of M and θ.    Advances in Civil Engineering e comparison between the ductility calculated from the experimental and the predicated skeleton curve is presented in Table 11, which shows that the errors are less than 10%.

Comparison of Hysteresis Loops.
e hysteretic rules for the energy-dissipation joint were simplified as a concave hexagon, the critical points of which were optimized by the GA as previously mentioned. e hysteretic curves obtained from the simplified model and the tests of the first four group in Table 1 are compared in Figure 19.
e comparison between the damping ratio obtained from the experimental and the simplified model is presented in Table 12. e errors are less than 10%, within acceptable range.

Simplified Restoring Force Model.
Combining the simplified skeleton curve and hysteretic rules of the energydissipation joints, a simplified restoring force model of the energy-dissipation joint is constructed, as presented in Figure 20.
In Figure  detailed loading and unloading path to establish the restoring force model consists of the following steps: (I) When the maximum load is less than yield moment M y , the energy dissipated by the joint is negligible. Both the loading stiffness and unloading stiffness are equal to the elastic moment stiffness K E or preyield moment stiffness K 1 , and the hysteretic rule is a straight line. (II) When the maximum load is between the yield moment and peak moment, the pinching in the hysteresis loop of the energy-dissipation joints becomes prominent. e shape becomes a reversed  S-shape, and the hysteretic rule in this segment can be simplified as a concave hexagon. e critical point values of this hexagon are presented in equations (11) to (20). e unloading and reloading path is 7 ⟶ 13 ⟶ 14 ⟶ 8 ⟶ 11 ⟶ 12 ⟶ 7. (III) When the joint reaches the softening stage, the energy dissipated by the joint increases with the displacement due to the lateral loading and the hysteresis loop becomes slightly chubby so that its pinching effect is weaker than that in the stage II. e shape of the hysteresis loop becomes a hexagon. e parameters of the hysteretic rule can be obtained from equations (11) to (20). e unloading and reloading path is 9 ⟶ 17 ⟶ 18 ⟶ 10 ⟶ 15 ⟶ 16 ⟶ 9.

Conclusion
Based on the experimental skeleton curves of the energydissipation joint in steel-engineered bamboo hybrid frames, a simplified symmetric trilinear skeleton curve and a simplified concave hexagon were established to relatively match with the experimental skeleton curve and hysteretic rules. e errors of ductility and damping ratio between the tests and calculated hysteretic model are all within 10%, and the skeleton curve and hysteresis loops of tests and proposed hysteretic model are in good agreement.

Data Availability
e XLSX data used to support the findings of this study may be accessed by emailing to the corresponding author. e available time will be six months after the publication since there are still some researches and articles for submission based on this manuscript and the original data.