Abstract

Spatial joints with end-plate connections show significant spatial coupling effects under spatial loading. Mechanical behaviour and failure modes of these spatial joints differ from those of planar joints. This study involved experiments and finite element analyses with respect to planar joints with end-plate connections under static load. The numerical results agreed well with the experimental data, and this verified the adequacy of the finite element analyses. Then, finite element models of the spatial interior joint, exterior joint, and corner joint were established to analyse the difference between the mechanical behaviour of spatial joints and planar joints. The component method was used to analyse components contributing to the initial stiffness of spatial joints. An initial rotation stiffness calculation model of spatial joints was proposed based on the deformation of joints. The findings indicated that results of the calculation models were in good agreement with those of the finite element analyses, and this proved that the calculation model proposed in this study could act as a reference method.

1. Introduction

Joints are an important part of an entire structure. Some previous studies [1, 2] examining beam-column connections focused on in-plane behaviour. However, the actual loading state of the joint should be spatial. The behaviour of the frame structure joints depends on the combined effect of the major-axis and minor-axis connections. Only planar connections were considered in actual structural design, and this does not reflect the actual mechanical behaviour of joints. Thus, it is necessary to separately consider coupling effects under spatial loading.

The moment-rotation behaviour of the end-plate joint with an extended stiffener was investigated by Tartaglia et al. [3, 4], and also, the design method was developed by Francavilla et al. [5]. Katalin and Miklós [6] investigated a modified component method that could be applied to the analysis of spatial joints. Cabrero and Bayo [7, 8] investigated the connection behaviour of a new spatial joint subjected to proportional spatial loading. Simoes da Silva [9] presented a conceptual design model of spatial joints. This conceptual design model was a combination of the calculation models of spatial joints based on the component method principle. Loureiro et al. [10] investigated spatial joints with end-connections on the major axis and minor axis. The results indicated that minor-axis connections offered a higher degree of improvement to the rotation stiffness of the joints when compared with that of major-axis connections. Gil et al. [11] used an experimental study and finite element analyses to investigate the interaction between major-axis and minor-axis connections of spatial joints. Zhou [12] conducted a set of experiments with respect to spatial composite joints.

Researchers also developed simplified models and calculation methods of initial rotation stiffness of planar joints, such as modifying the component method [13], application of new joints [14], and new loading condition [15]. The component method in the Eurocode was the most widely used analytical method. However, a few studies investigated the initial rotation stiffness of spatial joints. Calculation models for the initial rotation stiffness of steel joints loaded along the major axis were provided by EC3 [16], which were only applied to planar joints. The interaction stiffness effect between the major axis and minor axis was not investigated before. And, this effect could not be neglected in mechanical behaviour research of spatial joints.

This study investigated the rotation stiffness of the spatial joints considering the interaction stiffness effect between the major axis and minor axis. The accurate finite element modelling method was proposed by comparing the resultant curves with static experiments of planar joints. Finite element models of spatial joints were then constructed based on the verified modelling method. Differences in the initial rotation stiffness between spatial joints with end-plate connections and planar joints were evaluated. The theoretical rotational stiffness calculation method of the spatial joint was developed, and the result proved to be accurate.

2. Experimental Tests on the Joint Stiffness

Two kinds of specimens were designed. These joints were designed according to the code for design of steel structure [17], and the dimension was referenced from a practical engineer of an automobile repair workshop. According to the code, the column web and flange were designed to satisfy the shear deformation of the panel zone, and for the connection along the minor axis, the thickness of the column was increased to enhance the mechanical performance of the panel zone. The first specimen involved a major-axis interior joint with an extended end-plate connection that was numbered as MZ-1. The other specimen involved a minor-axis interior joint with flush end-plate connections that was numbered as MZ-2. Figure 1 shows the details of the panel zone.

The dimensions of the triangular stiffener on the right beam of MZ-1 corresponded to 75 mm × 75 mm × 10 mm (Table 1). Both columns and beams were composed of hot-rolled H-section steel with a strength grade of Q235B. The test mainly includes planar joints under symmetric loading and the material quality testing of the components. The properties of components are listed in Table 2. The high-strength bolt is of 10.9 grade, and the material property of the bolt was referenced from Guo [18] and listed in Table 3. Different connection types on the left and right of the joints were labelled as the left connection and right connection, respectively.

Figure 2 presents the test setup and details of the material test. A reaction frame and a distributive beam were used to fix the top of the column, and a hydraulic jack was used to support the other end. It should be noted that it was necessary for the point of the resultant force of the hydraulic jack and the loading beam to coincide with the centroid of the column. Synchronized loads were applied to both ends of the beams. An oil distributor system was adopted to guarantee the synchronization of loading. Additionally, brackets were installed close to both ends of the beams to prevent lateral buckling.

The entire loading process involved unidirectional symmetric loading, and the specific loading process is as follows: (1) The column was first subjected to a prescribed axial force (axial compressive ratio = 0.3) corresponding to 0.3fyAs, where fy and As denote the design strength of steel and the cross-section area of the column, respectively. (2) The ends of the beams were then subjected to unidirectional symmetric forces. Force control was employed prior to the yielding of the specimens. Following this, displacement control was applied at the ends of beams until there was local buckling at the upper flange of the beam, and the specimens could not continue to bear the load.

Mechanical parameter results are provided in Table 4. The moment-rotation response for all specimens is illustrated in Figure 3. Figure 3(a) shows that large differences existed between the left connection and right connection. These differences were attributed to the stiffener in the end-plate of the right connection. The stiffener enhanced both bearing capacity and stiffness of the joint and guaranteed the rotation capacity of the joint. The test data of the minor-axis end-plate connection are shown in Figure 3(b). This study only presented the test results of the left connection because specimen MZ-2 was bilaterally symmetrical and under symmetric loads. Although the ends of the specimen were asymmetric due to the nonsynchronization of loading positions and loading processes, they were still controllable.

Wang and Wang [19] reported that the displacement of the beam end has a nonlinear relationship with the rotation angle. The displacement of the beam end in the elastic stage consists of the joint rotation angle and deflection of the beam itself. In the plastic stage, the effect of beam plastic deformation should be considered. Therefore, the moment-rotation relationship of the connection was calculated using the following equation:where s is the vertical displacement of the beam end, θ is the rotation angle of the beam to the column connection, θ′ is the plastic rotation angle of the beam, M is the moment of the beam end, θL is the elastic deformation at the beam end in the elastic stage, θ′L is the plastic deformation at the beam end in the plastic stage, and EI is the bending rigidity of the beam.

3. Finite Element Model

A FE modelling software, ANSYS [20], was used to develop three-dimensional FE models of planar joints. The SOLID45 element was used to model the beams, columns, end-plates and bolts of the spatial joints. The element sizes of the bolts and the joint panel zone corresponded to 8 mm, and the element sizes of other parts were in the range of 20–30 mm. A surface-to-surface contact formulation was defined between the different parts of steel and between the bolts and the steel. The target faces and contact faces were modelled by the TARGE170 element and the CONTA174 element, respectively. The friction coefficient between the contact faces was 0.45.

The grade of the steel used was Q235B in which fy = 235 MPa, according to the code for design of steel structure [17]. Tables 2 and 3 provide a summary of the material properties of the H-section steel and bolts. Considering that the high-strength bolt is less ductile as compared with the H-section steel and the yield stress of the high-strength bolt is much higher, two constitutive models were provided to describe the H-section steel and the bolt, respectively. A trilinear isotropic hardening model was used as the constitutive model of the H-section steel. Grade 10.9 (M20) high-strength bolts with a yield strength of 940 MPa were used. A bilinear kinematic hardening model was used as the constitutive model of the bolt, and the value of the modulus of the hardening stage was 10% of the elasticity modulus as shown in Figure 4. For the H-section steel, fy = 260 Mpa, fu = 436 Mpa, εy = 0.001327, and εu = 0.010306. For the bolt, fy = 940 Mpa, fu = 1040 Mpa, εy = 0.004563, and εu = 0.009417.

Figure 5 shows the FE models and mesh division. The boundary conditions of the spatial joints in the FE model included fixing the bottom of the column and constraining the horizontal movement of the top of the column. An axial force with a value that was 0.3 times the yielding bearing capacity of the column cross-section was applied on the top of the column. The loading mode corresponded to symmetric loading at the ends of beams.

3.1. Comparison of the Test Results with Finite Element Analysis

Figure 6 shows the comparisons of the moment-rotation curves of specimen MZ-1 and specimen MZ-2. The left connection curves were compared to illustrate the effectiveness of the FE simulation. It can be observed that there was good agreement between the test results and the FE simulation in the elastic stage. Additionally, in the stage after yielding, there was some disparity due to the simplification of the constitutive behaviour of the materials and test errors. As shown in Figure 7, the final deformations of the test specimens and the numerical simulation were in good agreement. According to the comparison of curves and deformations, this finite element model had been proven to be sufficiently accurate for simulating these planar joints (Figure 8).

3.2. FE Models of Spatial Joints

The joints can be classified as three types, namely, middle column joints, edge column joints, and corner column joints, as shown in Figure 9. Tables 5 and 6 list the specific sizes and labels of the FE models of planar joints and spatial joints, respectively.

3.3. Influence of Different Loading Modes

Spatial loading represents the loading mode that involved concentrated loads at the ends of the beams connected with the joint. Planar loading represents the loading mode that only included the action of unidirectional loading. In order to improve computation efficiency, the FE models of the middle column joints and the edge column joints could be reduced by half based on the structural forms and the force conditions. The FE models and the mesh division of the middle column joints, edge column joints, and corner column joints are shown in Figure 9.

Figure 10 shows comparisons of the moment-rotation curves of the major-axis connections of the spatial joints under spatial loading and planar loading, and Table 7 lists the mechanical properties of the joints. The elastic ultimate moment could be defined by the yielding of components such as a beam flange, beam web, column flange, column web, end-plate, and high-strength bolt. And, the initial stiffness could be calculated by the ratio between moment and rotation at the elastic ultimate stage.

As shown in Figure 10 and Table 7, in the elastic stage, the initial rotation stiffness and elastic ultimate moments of the spatial joints under planar loading were almost identical to those under spatial loading. The maximum error of the initial rotation stiffness of the spatial corner column joints under different loading modes corresponded to 1.8%. However, the middle column joints and the edge column joints exhibited a greater difference in their bearing capacities under different loading modes. The joints under spatial loading yielded first, and then the plastic development was more rapid than that of the planar joints and the plastic region exceeded more than that of the planar joints. The bearing capacity of the joint under spatial loading was smaller than that under planar loading, and the maximum decrease of the corner column joints corresponded to 13%. The degradation of the bearing capacity of planar joints was slower than that of spatial joints. With respect to the middle column joints, the difference in the behaviour of the spatial load and that of the planar load was small.

Furthermore, as shown in Table 7, the initial stiffnesses of the middle column joints and edge column joints are larger than that of the corner column joints. With respect to the major-axis connection, the middle column joints and the edge column joints were under symmetric loading while the corner column joints were under unilateral loading. Hence, it was necessary to consider the effect of the shear deformation of the column web in the corner column joints which increased the relative rotation of the connection and reduced the rotation stiffness accordingly. The rotation stiffnesses of the middle column joints and the edge column joints were larger because it was not necessary to consider the effect of the shear deformation of the column web.

4. Initial Stiffness of Spatial Joints with End-Plate Connections

4.1. Flexural Stiffness of the Column Flange and End-Plate

An equivalent T-stub was considered to calculate the flexural stiffness of the column flange and the end-plate [21]. A very low level of prying force is produced by bolts in the elastic stage. Prying force usually occurs in the plastic stage, and it is a significant characteristic of the ultimate state of the joint [22]. Hence, it could be considered that the tightening effort of the bolt is continuous in the elastic stage, and based on this assumption, the bolted connection was assumed as a fixed connection and the prying action could be ignored. Figure 11 shows the mechanical model of the simplified equivalent T-stub.

The flexural stiffness of the column flange is calculated based on the theory of structural mechanics and material mechanics:where .

The flexural stiffness of the end-plate is expressed as follows:where leff,cf and leff,ep denote the effective lengths of the equivalent T-stubs of the column flange and the end-plate, respectively, which are calculated by a method proposed by Faella et al. [23]. The determined value is shown in Figure 12. Table 8 presents the values of leff,cf and leff,ep. Additionally, tcf and tep represent the thickness of the column flange and the end-plate, respectively.

Furthermore, dm in Figure 12 denotes the effective diameter of the bolt head, and it is calculated as dm = 1.5db where db denotes the nominal diameter of the bolt. Moreover, e denotes the distance between the centre of the bolt and the outer edge of the end-plate or the column flange, m denotes the distance between the bolt axis and the weld foot of the beam flange or the minor-axis connection plate, he denotes the effective height of the welded seam, and rc denotes the radius of the web-to-flange fillet of the column.

4.2. Tension and Compression Stiffness of the Column Web

Partial areas between the minor-axis bolts are subject to minor tension and compression forces and involve low elastic strain in the compressive and tensile regions. These components were ignored when compared with the other components that are in high stress. Additionally, the decrease in the calculation height of the column web led to the increase in the initial stiffness of these joints. With respect to the effect of the thickness of the column web, the spread angle of the stress that is transferred from the bolt-hole edge to the centre of the column web was assumed as θ, and its values are presented in equation (4). The areas of the ignored parts of the compressive and tensile regions were then calculated, and the corresponding illustration of the calculation is shown in Figure 13. Table 9 presents the calculated results where db in the table denotes the nominal diameter of the bolt, and the definitions of other symbols can be found in Figure 13.where trep is the thickness of the minor-axis end-plate and tcw is the thickness of the web.

Where c is the transverse bolt pitch of the up flange along the minor-axis, d is the longitudinal bolt pitch of the up flange along the minor-axis, s is the distance from the bolt-hole center of the bottom flange to the bottom flange center.

In the table, if the column is made of hot-rolled H steel, then beff,c = tbf + 2he,ep + 2tep + 2(tcf + rc). However, if the column is made of welded H steel, beff,c = tbf + 2he,ep + 2tep + 2(tcf + he,c), where he,ep and he,c denote the effective heights of welded seams of the end-plate and column, respectively. Additionally, rc denotes the radius of the web-to-flange fillet of the column. The definitions of he,ep, he,c, tcf, rc, and tbf are shown in Figure 14.

The principle of equivalent volume states that the ignored volume of the column web equals the decreased volume of the tensile and compressive stiffness components. According to this principle, the equivalent calculation height of the column web in tension or compression is calculated as follows:where hcw denotes the calculation height of the column web that is given as follows:

The tensile and compressive stiffness of the column web can be calculated by simplifying the column web as axial tension and axial compression plates. Furthermore, it is necessary to consider the effect of different loading modes on the stiffness of the joint.

The compressive stiffness of the column web is calculated as follows:where hcw,c denotes the equivalent calculation height of the column web in compression; tcw denotes the thickness of the column web; and beff,c denotes the effective height of the column web in compression. According to Eurocode 3 [16], φ denotes the influence coefficient of load and corresponds to 0.5 for symmetric loads, 1 for antisymmetry loads, and 1 for other loads.

The compressive stiffness of the column web is expressed as follows:where hcw,t denotes the equivalent calculation height of the column web in tension. This is equal to λ(p/)3 based on the principle that states that the distribution of the tension of the bolts depends on the stiffness of the column web and the minor-axis connection plate in conjunction with the principle that states that the column web and the minor-axis connection plate as shown in the present study have the same equivalent tensile deformation. Additionally, λ denotes the influence coefficient of bolt spacing. Each parameter is illustrated in Figure 12. The effective tensile width is shown in Figure 15, and its values are listed in Table 10.

4.3. Tensile Stiffness of Bolts

The study ignored the effect of the prying force of the bolts in the rotation stiffness calculation. Additionally, it was necessary to modify the suggested formula of the Eurocode to calculate the tensile stiffness of the bolts. Furthermore, the coefficient γ was introduced into the calculation model for high-strength bolts based on the effect of prestress. The value of γ is based on a previous study [24] and generally corresponded to 10. The tensile stiffness of the bolts can then be derived as follows:where As denotes the effective area of the bolt rod and is typically (in engineering) equal to 0.8 times the nominal area of the rod. Lb denotes the calculation length of the bolt. The meaning of thickness variation t can be found in the appendix.

4.4. Modification of the Rotation Stiffness Calculation Model of the Spatial Joints

The deformation of the upper connection includes the deformation of the bolts, end-plates, the bending deformation of column flanges, and tensile deformation of the column web. The deformation of the lower flanges could be reflected by the compressive deformation of the column web. The connection was in accordance with the plane section assumption of the steel beam and the assumption that the rotational centre of the joint was located in the center of the lower flanges along the thickness direction. The initial rotation stiffness calculation model of the spatial middle column joint is illustrated in Figure 16.

The total deformation of the spatial middle column joints can be divided into three parts, namely, Δt that corresponds to the tensile deformation of the column web and bolts, Δb that corresponds to the bending deformation of the end-plate and column flange, and Δc that corresponds to the compressive deformation of the column web.

Based on the above assumptions, the rotation of the joint in the elastic range is given by the following expression:where heq denotes the equivalent calculation height with a value that is equal to the distance between the centre of the upper flange of the beam and the rotational centre, and it can be expressed as hep = hb − tbf where hb denotes the height of the beam.

The calculation of the initial rotation stiffness of the joint is as follows:

4.5. Verification of the Rotation Stiffness Calculation Model of the Spatial Joints

Several FE models were built to verify the accuracy of equation (11). The specific forms and sizes of joints are specified in Table 11. For spatial joints of C series, the form of the major-axis connection remains unchanged; the variable parameters along the minor axis include the thickness of the end-plate, the bolt pitch, and the bolt diameter. The minor axis connections of D series are in different forms. The end-plate was extended at both ends, while the variable parameters along the minor axis include end-plate thickness, bolt-hole diameter, and bolt pitch. Comparisons between results by the theoretical analysis and results by FEM analysis are listed in Table 12.

As shown in Table 12, the results between theoretical calculation and FEM analysis results are in good accordance for most of the spatial joints, with errors less than ±5%. While for few of the joints, the relatively larger errors occurred due to the poor design of bolt spacing, such as for joint D04, there is a large bolt spacing difference between the vertical and horizontal axis. So generally, the theoretical equation proposed in this article is of good applicability for most reasonably designed joints.

5. Conclusion

This study involved analyzing the interaction between the major axis and the minor axis of spatial joints and the differences in the performances of the spatial joints and planar joints. Finally, a calculation equation for the initial rotation stiffness of spatial joints was presented based on the component method. The key conclusions of the study include the following:(i)The rotation stiffness of spatial joints under spatial loading typically exceeded that under planar loading. However, the bending moments of the spatial joints under spatial loading were lower than those under symmetric planar loading.(ii)The initial rotation stiffness and bearing capacity of the spatial joints were significantly higher than those of the planar joints. Mechanical performances between the spatial joints and planar joints are different. The components that contributed to the stiffness of the spatial joints were analysed to derive the stiffness calculation formula of each component based on elastic theory. Additionally, the stiffness of the specific components of the spatial joints was calculated, and compared with the FEM results, a relatively accurate formula of stiffness was presented.(iii)The calculation model and the theoretical equation of the initial rotation stiffness of joints were rebuilt based on the sources of the deformation of the joints. Furthermore, the accuracy of the calculation formula presented in this study was verified based on calculation examples and several FE models.

Nomenclature

As:The effective area of the bolt rod
beff,c:The effective height of the column web in compression
db:The nominal diameter of the bolt
dm:The effective diameter of the bolt head
E:Elasticity modulus
fu:Tensile stress
fy:Yield stress
hb:The height of the beam
hcw:The calculation height of the column web
hcw,c:The equivalent calculation height of the column web in compression
hcw,t:The equivalent calculation height of the column web in tension
he:The effective height of the welded seam
he,c:The effective heights of welded seams of the column
he,ep:The effective heights of welded seams of the end-plate
IT:The cross-sectional moment of inertia
Ki:Initial rotation stiffness
ktn:The tensile stiffness in the height of each row of bolts
Lb:The calculation length of the bolt
leff,cf:The effective lengths of the equivalent T-stubs of the column flange
leff,ep:The effective lengths of the equivalent T-stubs of the end-plate
lT:The calculation span of T-stub
M:Elastic ultimate moments
Mu:Plastic ultimate moments
m:The distance between the bolt axis and the weld foot of the beam flange
rc:The radius of the web-to-flange fillet of the column
tcf:Thickness of the column flange
tcw:The thickness of the web
tep:The thickness of the end-plate
th:Thickness of the bolt head
tn:The thickness of the bolt nut
trep:The thickness of the minor-axis end-plate
twh:The thickness of the gasket
Δb:The bending deformation of the end-plate and column flange
Δc:The compressive deformation of the column web
Δt:The tensile deformation of the column web and bolts
εy:Yield strain
λ:The influence coefficient of bolt spacing
φ:The influence coefficient of the load.

Data Availability

All the data used to support the findings of this study have been deposited in the figshare repository (https://figshare.com/s/da80ea5e386b41b4c092). They are available from the first author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (grant nos. 51778241; 51708226; 51638009; and 51978279), by the State Key Laboratory of Subtropical Building Science of South China University of Technology (grant nos. 2018ZB35, 2017ZB28, and 2017KD22), and by the Fundamental Research Funds for the Central Universities (Grant nos. 2017BQ086 and 2019ZD47).