Abstract
An integral finite element model of steel circular hollow section and welded hollow spherical joints is established to consider the restraint of welded hollow spherical joints on steel circular hollow section in this paper. The buckling modes are used as the initial geometrical imperfection forms of the steel circular hollow section. Based on the statistics method, a method for calculating the buckling capacity design value of the steel circular hollow section with uniform reliability considering the restraint of welded hollow spherical joints is established, and the formula for calculating error of this value at certain confidence is deduced. The relationships between the buckling capacity design value of steel circular hollow section and the structural parameters of the member and joint are determined quantitatively by linear regression method. A practical formula for the buckling capacity design value of steel circular hollow section considering the restraint of welded hollow spherical joints is derived. Compared to the results calculated by the design code indicates that if the difference of restraint effect which is caused by different welded hollow spherical joint is ignored, the buckling capacity value of steel circular hollow section calculated based on the unified effective length coefficient may be bigger or smaller than that of the actual case. Therefore, refined consideration should be placed on the different restraint of welded hollow spherical joints of different sizes on steel circular hollow section, and the restraint of the end joints should be accurately considered when determining the buckling capacity of steel circular hollow sections.
Introduction
Grid space structure is widely used in the roof structure of large public buildings such as airport terminal, gymnasium and auditorium due to its large span, light weight, high strength and beautiful shape. The bearing capacity of grid space structure is closely related to the mechanical properties of structural members and joints. As a large number of members in grid space structures bear axial compression loads, the member buckling capacity is always the key point in the structural design. The buckling capacity of slender members has been studied for many years,1,2 and it is also clearly specified in the grid space grid structure design codes. The stability checking calculation of compressive members is based on the stability coefficient. According to the current structure design codes, such as Chinese design code for steel structures, 3 the buckling capacity of the member bearing compression load is calculated based on the stability coefficient and effective length factor. The stability coefficient for certain section is usually determined by the column curves, and the column curves are calculated based on the free members at both ends with the consideration of initial imperfections. The effective length factor of the compression member is used to consider the influence of restraints caused by the joints. By the code method, the effective length coefficient of slender member is determined by the type of the joint. For example, the Chinese Technical Specification for Grid Space Structures 4 stipulates that for the structure using welded hollow spherical joints, the effective length coefficients are 0.9 for the chord members, 0.9 for web members near the supports and 0.8 for other web members; for the structure using bolted spherical joints, the effective length coefficient is 1.0 for each member. The effective length coefficient for compression members is related to the restraint effect caused by the joints.
In fact, even the same type of joint is used for the grid structure, the restraint effect caused by the joint to the member buckling capacity varies with the joint dimension. The restraint effect on the member is related to the type and dimension of the joint as well as the dimension of the member itself. Therefore, the effective length coefficient of the structural member should differ with the dimensions and type of the joint. If the actual effective length coefficient of members in the grid space structure is larger than that calculated by the code, the buckling capacity of the structural member will be overestimated in current design. Therefore, refined method by which the restraint of joints on the buckling capacity of the structural member can be actually considered needs to be established.
At present, scholars have carried out extensive research on the calculation method of bearing capacity for slender member. Galvão et al. 5 studied slender frames with semi-rigid connections and the results highlight the importance the stiffness of the semi-rigid connections on the buckling of these structures. But the initial geometrical imperfections of the members were not considered in the research process. There is also a calculation method for the buckling capacity of a single steel hollow section, which does not consider the joint restraint in the calculation process. Gardner and Nethercot6,7 carried out multiple sets of steel square, rectangular and circular hollow sections tests, which forms the basis for an explicit relationship between cross-sectional slenderness and cross-sectional deformation capacity The buckling capacities of steel square and rectangular hollow section were calculated using the numerical simulation method considering the initial imperfections. Characterisation of local plate imperfection amplitudes is described whereby a model originally devised for hot-rolled carbon steel cross-sections was re-calibrated and applied to stainless steel cross-sections. Of the features investigated, strain hardening characteristics of the material were identified as being primarily responsible for the differences in structural behaviour between hot-rolled and cold-formed sections.
Wang et al. 8 studied the influence of initial imperfections on the ultimate bearing capacity of axially compressed steel hollow sections, and gave a local stability calculation method for axial compression members. The test results show that the design of Q690 steel circular hollow section in the existing design standards is mostly conservative. McCann et al. 9 examined the response of EHS stub columns with parameters variation of imperfection amplitude and imperfection shape. Imperfection sensitivity was found to change with slenderness and aspect ratio. Hassan et al. 10 considered three methods of geometrical imperfections for square and rectangular structural steel hollow sections, base on which a new procedure of determining geometrical imperfection was proposed. Shen and Ahmer 11 formed a nonlinear variational approach to model local-global mode interaction. Qi et al. 12 established the member calculating model considering the buckling effect, and calculated the ultimate bearing capacity of the premier geometric faulty structures based on the probability and statistics method. It is shown that the nonlinear design bearing capacity of the lattice shell at a uniform probability level cannot be acquired by the methods suggested by the structure design code. However, the studies above did not examine the buckling capacity of the member in relation to the initial imperfection as well as the joint restraint. Therefore, refined consideration should be carried out.
In this paper, a method for calculating the buckling capacity design value of the steel circular hollow section with uniform reliability considering the restraint of welded hollow spherical joints is established, and the error formula of this design value at a certain confidence level is also derived. This method is verified by comparisons with experimental results recorded in published reference. The relationships between the buckling capacity design value of steel circular hollow section and the geometrical and material parameters of the member and joint are determined quantitatively by linear regression method. A practical formula for the buckling capacity design value of the steel circular hollow section considering the restraint of welded hollow spherical joints is derived.
Simulation of initial geometrical imperfections
Initial geometrical imperfection of the grid space structure members inevitably appears in the process of manufacture, transportation and assembly.13,14 Initial geometrical imperfections will reduce the buckling capacity of members. Therefore, the effect of initial geometrical imperfections must be taken into account when studying the buckling capacity of steel circular hollow sections with the consideration of joint restrains. The initial geometrical imperfections of the structure member must satisfy the displacement compatibility condition. In this paper, linear buckling mode is taken as the initial geometrical imperfection form of structural members, and l/300 (l is the member length) is taken as the maximum initial geometrical deviation value according to the Technical Specification for Grid Space Structures. 4 Linear buckling mode is a possible deformation form when a member buckles, and it satisfies displacement continuity conditions.
In reference, 8 a series of buckling capacity tests were carried out on steel circular hollow sections. Both the two ends of the steel circular hollow section specimen were welded on the steel plates of the upper and lower pedestals in the testing machine, which is shown in Figure 1. The numerical model of the specimen was established by the general finite element software ABAQUS. The initial geometrical imperfections of the specimen were simulated by buckling modes, and the buckling capacity of the specimen was calculated by arc length method with the consideration of geometrical nonlinearity and material nonlinearity. The finite element model is shown in Figure 2. In the finite element model, the mesh size is 1/15 of the longest edge of the steel circular hollow section, and the meshes are refined partially. The mesh element type is linear, finite film strain, complete integral quadrilateral shell element S4 in ABAQUS/Standard. The material model is nonlinear material model, and the yield strength is 690 MPa, which is consistent with the reference. 8

Steel circular hollow section specimen. 8

Finite element model.
The test result and numerical result of the buckling capacity value of the specimen are listed in Table 1.
Buckling capacity value of steel circular hollow section specimen.
From Table 1, it can be seen that the buckling capacity value of numerical result is very close to that of the test result. The maximum relative error is 5.86%, and the minimum relative error is only −0.81%. The initial geometrical imperfections of the steel circular hollow section specimen could not be accurately known, so the trial calculation needs to be carried out using multi-order buckling modes as the geometrical imperfection forms of the numerical model. For one steel circular hollow section there are different numerical buckling capacity values corresponding to different buckling modes. Generally there would be one value that is most close to that of the test result, and the buckling mode corresponding to this value best fits the actual geometrical imperfection of the very specimen. This indicates that it is feasible to use buckling mode as the initial geometrical imperfection distribution form of steel circular hollow section.
The formation of initial geometrical imperfections of space structural members is a stochastic process, so which order of buckling mode can be used to simulate the initial geometrical imperfections of the members is also a stochastic process.15–17 In order to accurately determine the buckling capacity value of steel circular hollow section in space structures considering the joint restraint, it is necessary to simulate the initial geometrical imperfections by buckling modes of sufficient orders. Therefore, parametric analysis needs to be carried out to derive a practical formula for the buckling capacity value of steel circular hollow section.
The parametric analysis is carried out by finite element program ABAQUS. The finite element numerical model is illustrated in Figure 3. As a commonly used type of joint for grid space structure, the welded hollow spherical joint is contained in the numerical model to consider its restraint on steel circular hollow section. The two hemispheres of welded hollow spherical joint are welded to each other, so an equivalent model can be established with the steel circular hollow section and two half joints which are rigidly fastened to rigid plates at both ends. In the finite element model, shell element model is used for the steel circular hollow section and the joint, and both the geometrical nonlinearity and material nonlinearity are taken into account.18,19 The complete axial load-displacement curve of steel circular hollow section member with geometrical imperfection considering restraint of welded hollow spherical joints is calculated.

Finite element model.
Considering the restraint of welded hollow spherical joints at the ends, the factors that may affect the buckling capacity of steel circular hollow section in the space structure are:
Parameters of the finite element models.
Based on the previous analysis, buckling modes of sufficient orders should be used as the initial geometrical imperfection forms of steel circular hollow section. For each of the 54 numerical models listed in Table 2, the 1st to 100th buckling modes are all used to create imperfect steel circular hollow sections. With the consideration of the joint restraint, the buckling capacity value of the steel circular hollow section with geometrical imperfection of certain buckling mode can be obtained as the peak value of axial load-displacement curve which is calculated by the arc length method. For example, the load-displacement curve of the model no. 4 with geometrical imperfection of its first buckling mode is shown in Figure 4. The peak value of this curve can be used as the buckling capacity value of this imperfect structural member.

Axial load-displacement curve.
Buckling capacity of steel circular hollow section
For each of the 54 models listed in Table 2, the 1st to 100th buckling modes are all used to create imperfect steel circular hollow sections. There would be one buckling capacity value corresponding to one imperfect model. So for one steel circular hollow section of the given size, it has 100 possible buckling capacity values. Theoretically there are countless imperfect forms for one steel circular hollow section, so these 100 values need to be analysed statistically to establish a reasonable method for the buckling capacity value of steel circular hollow sections with initial imperfection and joints restraint.
The initial geometrical imperfection of steel circular hollow section is unpredictable, so the buckling capacity value is a stochastic variable essentially. The 100 buckling capacity values calculated with the 100 buckling modes are just samples of this stochastic variable. The population of this stochastic variable can be analysed based on the sample using statistical method. Then the method for calculating the buckling capacity design value of the steel circular hollow section with uniform reliability considering the restraint of welded hollow spherical joints can be established accordingly.
Basic assumptions
Multiple independent factors cause initial geometrical flaw of steel circular hollow section in space structures. As each one of the factors plays a less prominent part, initial geometrical imperfection of steel circular hollow section is normally distributed as the probability theory describes.
Initial geometrical imperfection of steel circular hollow section can be represented by deviation vectors of finite element joint coordinate. According to Design Specification for Grid Space Structures, 4 some assumptions are made: The mean value of the finite element joint coordinate deviation of the member is 0; the tolerance of the joint coordinate deviation is l/300 (where l is the length of the steel circular hollow section), and the standard deviation of the joint coordinate deviation is a half of its tolerance. Accordingly the element joint coordinate deviation of steel circular hollow section distributes in normal fashion of N(0, (l/600)2), and range from −l/300 to l/300. The geometrical initial imperfections of the steel circular hollow section are independent multiple random variables, because deviation of the element joint coordinate is thought to be independent.
Distribution of buckling capacities of steel circular hollow sections
Take the fourth model of steel circular hollow section in Table 2 as an example. One hundred buckling modes are used to create initial geometrical imperfect steel circular hollow section. A 100 sized sample of the buckling capacity of steel circular hollow section can be obtained by FEM analysis. Based on this sample, the population distribution of the buckling capacity of this steel circular hollow section considering the joint restraint can be estimated.
Based on the analysis in Section ‘Basic assumptions’, it can be assumed that the bearing capacity values of imperfect steel circular hollow section considering the joint restraint are normally distributed, and a further nonparametric K-S (Kolmogorov-Smirnov) test is carried out for the verification.
Let the cumulative frequency distribution function of the sample be
Where
The inspection level of 0.05 is used to test the population distribution. Tn = 0.050, p = 0.675 and p > 0.050, it indicates that the buckling capacity values of the imperfect steel circular hollow section considering the joint restraint are normally distributed.
The mean µ and variance
Where
Create a likelihood function as following:
Constitute equations set as follows:
Use
The estimators of the mean µ and standard variance σ of the population of the example imperfect steel circular hollow section considering the joint restraint are 533 and 36 kN respectively.
Buckling capacity design value with uniform reliability
Use
Where
The buckling capacity design value of steel circular hollow section is defined by:
The buckling capacity design value
Three statistical samples of buckling capacity values of the No. 4 steel circular hollow section are calculated using 1st~100th, 101th ~200th and 201th ~300th buckling modes as the initial geometrical imperfection forms. The minimum values of the three samples are listed in Table 3 respectively, and the corresponding reliabilities calculated based on equation (6) are also listed in the fourth column.
Minimum value and its reliability of statistical samples of buckling capacity values.
It can be seen from Table 3 that different initial geometrical imperfections of steel circular hollow section will lead to different buckling capacity samples. If the minimum value of the sample is used as the buckling capacity design value, it depends on the selected initial geometrical imperfection of the steel circular hollow section, and a stable buckling capacity design value of steel circular hollow section with uniform reliability cannot be obtained. It is also not economical to just take the minimum value of the sample as the buckling capacity design value without consideration of reliability.
Error of the buckling capacity design value
Theoretically, if buckling capacity values of infinite imperfect steel circular hollow sections are calculated, the accurate value of the bearing capacity design value
The buckling capacity values of steel circular hollow sections distribute in the normal fashion of
Where
According to equation (8), the error of the buckling capacity design value is relative to the reciprocal of the arithmetic square root of the sample size.
According to equations (7) and (8), buckling capacity design values and relative errors of 54 steel circular hollow sections considering welded hollow spherical joint restraints listed in Table 2 can be calculated. The results are shown in Table 4.
Buckling capacity design values and relative errors of steel circular hollow sections.
Buckling capacity of steel circular hollow section considering the restraint of welded hollow spherical joint
It can be seen from the analysis in Section ‘Buckling capacity of steel circular hollow section’ that it is necessary to calculate the sample of plenty of data by statistical method when determining the buckling capacity design value of steel circular hollow sections considering the joints restraint. Although this method is fully justified by the stochastic characteristic, it is not practical due to the massive calculation work. A formula for the buckling capacity design value of steel circular hollow sections considering the joints restraint needs to be derived. Therefore, the influence of the steel circular hollow section length
The buckling capacity design value

The curve of
Each type of the parameter of model 11–model 20 has the same value except for the cross section wall thickness

The curve of
Each type of the parameter of model 21–model 30 has the same value except for the cross section outer diameter

The curve of
Each type of the parameter of model 31–model 40 has the same value except for the wall thickness

The curve of
Each type of the parameter of model 41–model 50 has the same value except for the outer diameter

The curve of
Each type of the parameter of model 51–model 60 has the same value except for the material design strength

The curve of
Construct dimensionless variables
Therefore, the expression of
By equation (10) the restraint of welded hollow spherical joint on steel circular hollow section is finely considered. It corresponds with the actual case. Based on the current design code, 3 the restraint of welded hollow spherical joint on buckling capacity of steel circular hollow section is considered by a unified effective length coefficient. The different restraint of welded hollow spherical joints of different sizes on steel circular hollow section is not finely taken into account. The buckling capacity values of model 2, model 13, model 26, model 36, model 43 and model 53 in Table 2 are respectively calculated based on the code and equation (10). The calculation results are shown in Table 5.
Buckling capacity values of steel circular hollow section based on the code[3] and Eq. (10).

Relationship curve of
It can be seen from Table 5 that the buckling capacity values of model 2, model 36, model 43 and model 53 calculated by the code are smaller than those calculated by equation (10). The buckling capacity values of model 43 calculated by the code is 29.35% smaller than that calculated by equation (10). The buckling capacity values of model 13 and model 26 calculated by the code are bigger than those calculated by equation (10). The buckling capacity values of model 13 calculated by the code is 18.93% bigger than that calculated by equation (10).
Therefore, when the restraint of welded hollow spherical joints is finely considered, the buckling capacity value of steel circular hollow section could be quite different from that calculated based on the current code. If the difference of restraint effect which is caused by different welded hollow spherical joint is ignored, the buckling capacity value of steel circular hollow section calculated based on the unified effective length coefficient may be bigger or smaller than that of the actual case. If the buckling capacity value calculated by the code is smaller than that calculated by equation (10), steel circular hollow section of bigger size would be used in the design, which makes the design uneconomical. If the buckling capacity value calculated by the code is bigger than that calculated by equation (10), steel circular hollow section of smaller size would be used in the design, which makes the design unsafe. Therefore, refined consideration should be placed on the different restraint of welded hollow spherical joints of different sizes on steel circular hollow section in the design, and the restraint of the end joints should be accurately considered when determining the buckling capacity of steel circular hollow sections.
Conclusions
In this paper, buckling mode is used as the initial geometrical imperfection form of steel circular hollow section. There are different numerical buckling capacity values corresponding to different buckling modes for one steel circular hollow section. Generally, there would be one value that is most close to that of the test result. This indicates that the buckling mode corresponding to this value best fits the actual geometrical imperfection, and it is feasible to use buckling mode as the initial geometrical imperfection distribution form of steel circular hollow section.
An integral finite element model of steel circular hollow section and welded hollow spherical joints is established to consider the restraint of welded hollow spherical joints on steel circular hollow section. The buckling modes are used as the initial geometrical imperfection forms of the steel circular hollow section when parametric analysis of the integral finite element model is carried out.
By nonparametric test, the buckling capacity value of geometrical imperfect steel circular hollow section considering the restraint of welded hollow spherical joint is verified to distribute normally. A new way of calculating the buckling capacity design value of steel circular hollow section at a uniform probability with the consideration of the restraint of welded hollow spherical joint is inferred, and its error expression at certain confidence level is also proposed.
The relationships between the buckling capacity design value of steel circular hollow section considering the joint restraint and the geometrical and material parameters of the member and joint are determined quantitatively by linear regression method. A practical formula for the buckling capacity design value of the steel circular hollow section considering the restraint of welded hollow spherical joints is derived.
When the restraint of welded hollow spherical joints is finely considered, the buckling capacity value of steel circular hollow section could be quite different from that calculated based on the current code. If the difference of restraint effect which is caused by different welded hollow spherical joint is ignored, the buckling capacity value of steel circular hollow section calculated based on the unified effective length coefficient may be bigger or smaller than that of the actual case. According to the calculation example of this paper, the relative difference of the stability bearing capacity of steel circular hollow section calculated by these two methods range between −29.35% and 18.93%. Refined consideration should be placed on the different restraint of welded hollow spherical joints of different sizes on steel circular hollow section, and the restraint of the end joints should be accurately considered when determining the buckling capacity of steel circular hollow sections.
One concern about the findings is that there are many types of joints applied in spatial latticed structures, and they offer different restraint. This paper focuses on welded hollow spherical joints, and future research should be undertaken to study the design stability bearing capacity of steel circular hollow section considering the restraint effect of other kinds of joints, such as bolt-ball joint, intersecting joint and cast steel joint.
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The authors are grateful for the financial support from National Key Research and Development Plan of China (2021YFB2600500).
