Failure mechanisms and reinforcement of pipe-framed solar greenhouses under snow loads
Failure mechanisms
The load-displacement curves of the structure under snow loads are shown in Fig. 7. For each example in the figure, the X-coordinate denotes the vertical displacement of the ridge node, and the Y-coordinate denotes the snow load. The maximum load on the load-displacement curve is the ultimate bearing capacity of the structure. Under the same load conditions, the ultimate bearing capacity of the greenhouse with flat elliptical hollow section is the highest, followed by the hat-shaped section and the rectangular hollow section is the smallest. The main reason for this phenomenon is that the moment of inertia of the flat elliptical section along the strong axis is greater than that of the hat-shaped section and the rectangular hollow section (as listed in Table 2). For the same reason, the ultimate bearing capacity of the greenhouse with hat-shaped section is close to that of the greenhouse with hollow rectangular section. A similar phenomenon was also founded in the auther’s previous study41. Under uniform snow load, the ultimate bearing capacities of the greenhouse skeleton with flat elliptical hollow section, hat-shaped cross-section and rectangular hollow section are 546.6 N m− 2, 491.8 N m− 2 and 485.3 N m− 2, respectively. Under non-uniform snow load, the ultimate bearing capacities of the greenhouse skeleton with flat elliptical hollow section, hat-shaped cross-section and rectangular hollow section are 178.2 N m− 2, 170.0 N m− 2 and 168.6 N m− 2, respectively. From the above results it can be concluded that when the amount of steel is the same, it is recommended to give priority to the flat elliptical hollow section. A similar conclusion was also reached by Liu et al.42, who compared the safety performances of greenhouse skeletons with circular hollow section, rectangular hollow section, hat-shaped section and flat elliptical hollow section. In addition, comparing Fig. 7a and b, we can find that the ultimate bearing capacities of the greenhouse skeleton with flat elliptical hollow section, hat-shaped cross-section and rectangular hollow section under uniform snow load are about 3.1, 2.9 and 2.9 times that under non-uniform snow load. These results mean that the greenhouse is sensitive to non-uniform snow loads. Therefore, more attention should be paid to the safety of greenhouses subjected to non-uniform snow loads.

Load-displacement curves under: (a) uniform and (b) non-uniform snow loads.
To explore the failure process of the greenhouse under snow loads, the greenhouse skeleton with flat elliptical hollow section is selected as an example. Figure 8 shows the stress development process of the greenhouse skeleton under non-uniform snow load. Based on Fig. 7b, the failure process of greenhouses under snow loads can be classified by three stages: elastic stage (OA), elastic-plastic stage (AB) and failure stage (BC). When the snow load is less than 140.9 N m− 2 (before point A), the greenhouse skeleton is in the elastic stage. The displacement increases linearly with the increase of snow load. When the snow load increases to 140.9 N m− 2 (Fig. 8a), the greenhouse skeleton begins to enter the elastic-plastic stage. The section flanges at the corner of the south roof, at the mid-span, at the junction of the north roof and the column, and at the bottom of the column begin to yield. As the snow load continues to increase, the yield zone expands from the flange to the web. When the snow load increases to 178.2 N m− 2 (Fig. 8b), the greenhouse skeleton reaches the ultimate bearing capacity. The junction between the north roof and the column, and the bottom of the column enter the full-section yield state. Once the greenhouse skeleton reaches its ultimate bearing capacity, the displacement will continue to increase even if the snow load is reduced. When the snow load is reduced to 169.3 N m− 2 (Fig. 8c), the ridge node displacement reaches 170.7 mm. The section at the corner of the south roof, at the mid-span, at the junction of the north roof and the column, and at the bottom of the column all enter the full-section yield state, leading to the greenhouse skeleton failure.

Stress development process of the greenhouse skeleton with flat elliptical hollow section under non-uniform snow load: (a) snow load s0 = 140.9 N m− 2 (point A in Fig. 7b); (b) snow load s0 = 178.2 N m− 2 (point B in Fig. 7b); (c) snow load s0 = 169.3 N m− 2 (point C in Fig. 7b).
To further explore the failure mechanism of the greenhouse under non-uniform snow loads, Fig. 9 presents the axial force and bending moment distributions of the greenhouse at its peak load. From Fig. 9a, we can see that the axial force values are all negative, indicating that the greenhouse subjected to compression force. The maximum axial force occurs at the south roof end. Compared with axial force, the bending moment distribution is relatively complex. From Fig. 9b, we can see that the bending moment values at the corner of the south roof, the mid-span, the junction of the north roof and the column, and the bottom of the column are relatively large. These locations are consistent with those in Fig. 8c where full-section yield occurs. Comparing Figs. 8c and 9b, we can conclude that the main reason of the greenhouse failures may be attributed to excessive bending moment caused by snow loads.

Internal force distributions of the greenhouse skeleton with flat elliptical hollow section at the peak non-uniform snow load: (a) axial force distribution and (b) bending moment distribution.
Reinforcement methods
According to the results in Sect. “Failure mechanisms”, under non-uniform snow load, the ultimate bearing capacity of the greenhouse skeleton with flat elliptical hollow section is 178.2 N m− 2, which is less than the basic snow pressure of 300 N m− 2 once in 10 years in Shenyang city. Therefore, in order to ensure the safety performance of this type of greenhouse under snow loads in Shenyang city, it is necessary to strengthen its structure. In this section, taking greenhouse skeletons with flat elliptical hollow section as examples, three reinforcement methods are proposed and their effects on the structure safety are discussed.
Reinforced by one Brace
One of the most common reinforcement methods is to add one brace between the south roof and the north roof or between the south roof and the column, as shown in Fig. 10. However, the bracing arrangement is mainly based on the builder’s experience and lacks the support of design theory. To investigate the effects of bracing position on the structure safety, the bracing arrangement method is shown in Fig. 10. The brace is 48 × 2.0 mm circular hollow section (CHS) steel tube. One end of the brace is pin connected to the south roof (S1, S2, S3, S4, S5), and the other end is pin connected to the north roof (N1, N2, N3) or the column (N4, N5). And there are 25 bracing arrangements in total. For the sake of expression, the greenhouse with bracing is named after the numbers at both ends of the brace. For example, the greenhouse model in Fig. 10 is named S4N3. The model without bracing is named original model.

Bracing arrangement method of the greenhouse.
Figure 11 shows the effects of bracing position on the ultimate bearing capacity of greenhouses under snow loads. The dashed line represents the ultimate bearing capacity of the original greenhouse without bracing. It can be seen from Fig. 11 that the ultimate bearing capacity of greenhouses cannot always be improved with bracing. On the contrary, when the brace is not set appropriately, it may reduce the ability of the structure to resist snow loads. Only when one end of the brace is connected to the column (S1N4, S1N5, S2N4, S2N5, S3N4, S3N5, S4N4, S4N5, S5N4, and S5N5) can the ultimate bearing capacity of the structure be increased, regardless of uniform and non-uniform snow loads. Under uniform snow load (Fig. 11a), the optimal model is S1N5, with an ultimate bearing capacity of 670.5 N m− 2, which is 22.7% higher than the original model. Under non-uniform snow load (Fig. 11b), the optimal model is also S1N5, with an ultimate bearing capacity of 217.6 N m− 2, which is 22.1% higher than the original model.

Effect of bracing on the ultimate bearing capacity of greenhouses under (a) uniform and (b) non-uniform snow loads.
To explore the influence mechanism of the bracing on the ultimate bearing capacity of greenhouses, taking the S1N1, S1N2, S1N3, S1N4 and S1N5 models as examples, Fig. 12 compares the load-displacement curves of reinforcement models and the original model. As illustrated in Fig. 12, when the snow load is small, the displacement increases linearly as the load increases, that is, the structure is in the elastic stage. In the elastic stage, the slope of load-displacement curves of models S1N4 and S1N5 is larger than that of the original model, while the slope of models S1N1, S1N2, S1N3 is smaller than that of the original model. These phenomena indicate that the bracings of models S1N4 and S1N5 increase the initial structural stiffness of the original model, while the bracings of models of S1N1, S1N2, S1N3 reduce the initial structural stiffness of the original model. In addition, the initial structural stiffness of model S1N5 is the largest, while the initial structural stiffness of models S1N1, S1N2, S1N3 are not much different. These phenomena are consistent with the results in Fig. 11.

Effect of the bracing on the load-displacement curves of greenhouses under (a) uniform and (b) non-uniform snow loads.
To further explore the influence mechanism of the bracing on the ultimate bearing capacity of greenhouses under non-uniform snow loads, taking the S1N5 model as examples, Fig. 13 presents the internal force distributions of the greenhouse at the peak load. Comparing Figs. 9 and 13, it can be found that the existence of the brace changes the internal force distribution of the greenhouse. In terms of axial force distribution, as shown in Fig. 13a, the greenhouse frame above the brace is subjected to axial tension, while the greenhouse frame below the brace is subjected to axial compression. While for the original model, the whole greenhouse frame is subjected to axial compression, as shown in Fig. 9a. In terms of bending moment distribution, as shown in Fig. 13b, the existence of the brace reduces the bending moment value of the north roof skeleton, and the maximum bending moment appears at the joint of the column and the brace.

Internal force distributions of the model S1N5 at the peak non-uniform snow load: (a) axial force distribution and (b) bending moment distribution.
According to the above results, we can conclude that: (1) Setting up one brace does not always improve the structure’s ability to resist snow loads. Improper bracing arrangement will reduce the ultimate bearing capacity of the greenhouse, thereby accelerate failure progress of the structure; (2) When arranging one brace, it is recommended to connect one end of the brace to the south roof and the other end to the column, and avoid connecting the brace to the north roof.
Reinforced by one lattice column
Another reinforcement method is to replace the single tube column with a lattice column on the basis of the original model. The schematic of the lattice column model of the greenhouse is shown in Fig. 14. The lattice column is CHS steel tube with the size of 48 × 2.0 mm. To facilitate the expression, the greenhouse reinforced by one lattice columns is referred as the lattice column model.

The schematic of the lattice column model of the greenhouse.
Figure 15 compares the load displacement curves of the lattice column model and the original model under uniform and non-uniform snow loads. The ultimate bearing capacity of the lattice column model is obviously greater than that of the original model under both uniform and non-uniform snow loads. Under uniform snow load, the ultimate bearing capacity of the lattice column model is 984.7 N m− 2, representing an increase of 80.2% compared to the original model. Under non-uniform snow load, the ultimate bearing capacity of the lattice column model is 326.9 N m− 2, representing an increase of 83.4% compared to the original model.

Comparison of load-displacement curves between the lattice column model and the original model under (a) uniform and (b) non-uniform snow loads.
To explore the influence of the lattice column on the stiffness of the greenhouse, Fig. 16 shows the deformed shapes of the original model and the lattice column model under non-uniform snow load. As shown in Fig. 16, under the same load, the deformation of the lattice column model is significantly smaller than that of the original model. This phenomenon is consistent with the results in Fig. 15, in which the slope of the load displacement curve of the lattice column model is clearly greater than that of the original model. In particular, there is almost no deformation in the north roof and the column of the lattice column model. This is due to the fact that the stiffness of the greenhouse structure significantly improves when the lattice column replaces the single tube column. As a result, smaller structural displacements and improved ultimate bearing capacity are achieved.

Deformed shapes of the original model and the lattice column model under non-uniform snow loads (The snow load is 150 N m− 2).
To further explore the failure mechanism of the lattice column model under non-uniform snow load, Fig. 17 presents the von Mises stress distribution of the greenhouse at the peak load. As depicted in Fig. 17, under the peak load, the section flanges at the corner of the south roof, at the mid-span, and at the bottom of the column entered to yield. While fully sectional yielding occurred at the junction of the column and the north roof. Therefore, more attention should be paid for junction of the column and the north roof.

The von Mises stress distribution of the lattice column model under the peak non-uniform snow load.
Reinforced by one temporary column under the South roof
Before extreme snow loads come, setting up temporary columns is also a reinforcement method to improve the ultimate bearing capacity of the greenhouse. Figure 18 shows a schematic diagram of a temporary column arrangement. In Fig. 18, the horizontal distance from the temporary column to the south roof end is defined as HD. To explore the effect of the temporary column position on the ultimate bearing capacity of the greenhouse, the temporary column moved from the south roof end to the ridge in steps of 0.5 m. The temporary column is CHS steel tube with the size of 48 × 2.0 mm. Considering the detachability of the temporary column, the connection between the upper end of the column and the greenhouse skeleton was simplified to a hinge. And the column was assumed to be hinge-connected with the ground.

Diagram of the temporary column arranged on the greenhouse.
Figure 19 shows the effects of the temporary column position on the ultimate bearing capacity of the greenhouse under snow loads. Under uniform snow load (Fig. 19a), as HD increase, the ultimate bearing capacity of the greenhouse increases at first and then decreases. When HD = 4.5 m, the ultimate bearing capacity reaches the maximum value of 2660.1 N m− 2, which is about 4.9 times that of the original model. Under non-uniform snow load (Fig. 19b), when HD is less than 4.5 m, the ultimate bearing capacity increases with the increase of HD. When HD = 4.5 m, the ultimate bearing capacity reaches the maximum value of 995.0 N m− 2. There is almost no change in the ultimate bearing capacity when HD increases from 4.5 m to 5.0 m. When HD is greater than 5.0 m, the ultimate bearing capacity decreases with the increase of HD. Based on the above results, we can conclude that for a greenhouse with a span of 10 m, the optimal position of the temporary column is 4.5 m away from the south roof end.

Effect of temporary column on the ultimate bearing capacity of greenhouses under (a) uniform and (b) non-uniform snow loads.
Reinforcement efficiency analysis
As discussed in Sect. “Reinforcement methods”, although the ultimate bearing capacity of the greenhouse structure may be improved by setting bracings, lattice columns and temporary columns, the increase in the steel amount cannot be ignored. To identify the most suitable reinforcement method, the concept of reinforcement efficiency P is introduced here:
$$P=\Delta F_{cr}-\Delta S$$
(3)
where \(\Delta F_{cr}\) is the ultimate bearing capacity improvement factor, \(\Delta S\) is the steel amount increase factor for the greenhouse. \(\Delta F_{cr}\) is defined as:
$$\Delta F_{cr}=(F_{cr}^r-F_{cr}^0)/F_{cr}^o$$
(4)
where \(F_{cr}^r\) is the ultimate bearing capacity of the reinforcement model. \(F_{cr}^o\) is the ultimate bearing capacity of the original model.
\(\Delta S\) is defined as:
$$\Delta S=(S^r-S^0)/S^o$$
(5)
where \(S^r\) is the steel amount of the reinforcement model. \(S^o\) is the steel amount of the original model. In this study, a positive P value signifies that the reinforcement method is economical and effective. Conversely, a negative P value indicates that the reinforcement method is uneconomical. The situation, particularly in the case of a negative \(\Delta F_{cr}\) and positive \(\Delta S\), is most uneconomical, as an increase in the steel amount results in a reduction in ultimate bearing capacity.
Figure 20 shows the effect of the bracing position on reinforcement efficiency. It can be clearly seen from Fig. 20 that most reinforcement efficiencies are negative, regardless of uniform or non-uniform snow loads. Under uniform snow load (Fig. 20a), only models S1N5, S2N5 and S3N5 have positive reinforcement efficiencies. The maximum reinforcement efficiency value appears in model S1N5, which is 7.0%. The worst model is S5N3 with a reinforcement efficiency value of -23.1%. Similar with the results under uniform snow load, the maximum and minimum reinforcement efficiencies under non-uniform snow load (Fig. 20b) also appear in model S1N5 and model S3N5, which are 6.4% and − 22.5%, respectively. Based on the above results, we can conclude that S1N5 is the optimal bracing reinforcement model.

Effect of bracing on the reinforcement efficiency of greenhouses under (a) uniform and (b) non-uniform snow loads.
For the lattice column model, the ultimate bearing capacity improvement factors under uniform and non-uniform snow loads are 80.2% and 83.4%, respectively. The steel amount increase factor is 20.3%. After calculation, the reinforcement efficiency of the lattice column model under uniform snow and non-uniform snow load is 59.9% and 63.1% respectively. Compared with the S1N5 model, the reinforcement efficiency of the lattice column model under uniform and non-uniform snow loads is about 8.6 and 9.9 times that of the S1N5 model. Therefore, under the condition that crop cultivation operation is not significantly affected, the lattice column model is preferred to improve the snow resistance ability of greenhouse.
Figure 21 shows the effect of the temporary column position on the reinforcement efficiency. Comparing Figs. 19 and 21, it can be found that the change trend of reinforcement efficiency with the temporary column position is similar to that of the ultimate bearing capacity with the temporary column position. Under uniform snow load (Fig. 21a), the reinforcement efficiency P increases firstly and then decreases at the range form HD = 0.5 m to HD = 8.0 m. As for HD = 4.5 m, the reinforcement efficiency P reaches the maximum value of 365.3%. Under non-uniform snow load (Fig. 21b), the reinforcement efficiency P also reaches its maximum value when HD = 4.5 m. The maximum reinforcement efficiency is 437.1%. Based on the above results, we can conclude that for a greenhouse with a span of 10 m, the most cost-effective position of the temporary column is 4.5 m away from the south roof end.

Effect of temporary column position on the reinforcement efficiency of greenhouses under (a) uniform and (b) non-uniform snow loads.
link
