Abstract
For calculating the thin-walled closed curved box girder caused by the temperature gradient of the internal force and displacement, based on the fundamental differential equation of the curve beam and the principle of minimum energy, set a reverse statically indeterminate simply supported curve beam as the basic structure, consider the warping effect of the closed curve box girder, and put forward a kind of plane curve beam temperature deformation simple analytical calculation method. Compared with the finite element calculation results, the relative error of the analytical calculation results is less than 5%. It is concluded that the analytical method has sufficient accuracy in calculating the out-of-plane deformation of the thin-walled closed curved box girder under the temperature gradient.
1. Introduction
In recent years, with the construction of urban overpasses, curved beams are increasingly used in modern engineering structures. Scholars at home and abroad have studied the coupling characteristics of bending and torsion of curved beams [1–6] and the factors such as warping, distortion, and shear lag of thin-walled closed curved box girders [7–10]. Temperature effect is one of the main causes of bridge engineering problems. The research on the temperature effect of bridges mainly focuses on straight bridges [11–13], and the research on the temperature effect of curved beams is less. Literature studies [14–18] used the finite element method to analyze the stress and deformation of the curved beam under temperature load, and literature studies [19–22] used the analytical method to analyze the temperature effect of the curved beam but do not consider the warping effect of the curved beam under temperature in the out-of-plane direction.
Based on the basic differential equation of the curved beam and the principle of minimum potential energy, this paper presents a primary torsional statically indeterminate simply supported curved beam as the basic structure. A simple analytical method for calculating the vertical temperature gradient deformation of plane curved beams is presented.
2. Deformation of the Simply Supported Statically Indeterminate Curved Beam under the Vertical Temperature Gradient
2.1. Basic Assumptions
To analyze the temperature effect of the thin-walled closed curved box girder, the following basic assumptions are adopted:(1)Warping deformation is only considered comprehensively with vertical bending deformation(2)The strain and strain curvature of the curved box girder are considered as curved beams with small curvature, and the deformation is in a small deformation range(3)Distortion of the cross section is not considered(4)The temperature does not change along the axial direction of the beam, but only exists on the cross section (5)The horizontal displacement of the structure caused by torsion is not considered(6)The warping function sum and are used for theoretical analysis, and the assumption of plane section in pure torsion theoretical analysis is not used
The coordinate system of the curved beam is shown in Figure 1.
2.2. The Solution of Vertical Deflection Torsion Angle and Warping Function
According to [4, 7], the geometric equations of vertical bending and torsion of curved beams are as follows:where is the vertical bending curvature; is the converted twist rate, which is the torsional angle per unit length; is the closed section warping function corresponding to ; is the vertical deflection; is the section torsion angle; is the closed section warping function corresponding to ; is the center angle coordinate of the curved beam; and is the corresponding radius of the curved beam.
The physical equation iswhere is the bending moment in the axial direction; is the pure torsional torque; is the warping torsional torque; is the warping double moment; is the sectional bending moment of inertia in the axial direction; is the section pure torsion constant; is the polar moment of inertia of the section; and is the moment of inertia of section warping. Total torque is .
Considering a curved equal-section beam with center angle and taking the midpoint of the circular arc as the coordinate zero point, determined by formulas (1)–(7), under the temperature effect, the potential energy of the structure iswhere is the coefficient of thermal expansion. The first four terms of the potential energy expression are strain energy, and the fifth term is temperature load potential energy. According to the principle of minimum potential energy, the possible deformation is the deformation that makes the potential energy reach the stationary value. For statically indeterminate simply supported curved beams, the function curves , , and (or , , and ) are needed to make the strain energy get the minimum value. Let so that formulas (9)–(16) are obtained, where (9)–(11) are governing equations and (12)–(16) are boundary conditions. Get the governing equation:
Boundary conditions are
At the bearing, which is , because the bearing limits the torsion and vertical displacement, , , and . Boundary conditional expressions (13), (14), and (16) are naturally satisfied. However, simply supported bearings cannot limit the displacement and , so and are not 0 in , so in ,
First, it is solved by governing equations (9)–(11) and obtained by (10):
Substitute equation (19) into equation (9) to get
From (11),
Substitute formula (21) into formula (20) to get
Solving it,and among them,
Because the symmetric structure has symmetric stress and deformation under symmetric temperature load, is an even function. can be known from symmetry.
Substitute (24) into (21) to get
Solving it,
Similarly, is an even function, and are known by symmetry. Substitute equations (24) and (27) into equation (19), sort out, and solve them.
In the same way, is an even function, and are known by symmetry, obtained by sorting out (24), (26), (27), and (28).
Considering boundary conditions (17) and (18), substituting equations (18), (30), and (31) into equation (17), we get
For most sections, there are
Therefore, substitute equation (29) into equation (18), and get .
By the boundary conditions of the structure,
From equations (30), (34), and (35), we get
So,
It is obtained from formulas (30), (34), (35), and (36) that
Therefore, the solution of vertical deflection , section torsion angle , and warping function is
If the warping effect in torsion is ignored, the potential energy of the structure is
It can be obtained by a similar calculation process as above.
Because the coordinates on the symmetric box section are left and right antisymmetric about the section, for the temperature distribution of the left and right symmetric section, such as the temperature gradient along the beam height specified in General Code for Design of Highway Bridges and Culverts (JTG D60-2015)[23], there are.
At this time, equations (41) and (42) are completely consistent with equations (45) and (46). Therefore, under the assumptions mentioned above, the deformation calculation results of the thin-walled closed curve box girder with and without considering the warping effect are the same for the temperature distribution of the left and right sections. On the contrary, the deformation calculation results of the thin-walled closed curve box girder with and without the warping effect are different for the temperature distribution of the left-right asymmetric section. It can be seen that, for the thin-walled closed curve box girder, due to the existence of the warping effect, the transverse temperature gradient distribution in the plane will affect the vertical deflection, torsion angle, and warping out of plane. This property is different from that of the straight beam, which shows that the transverse temperature gradient is also of great significance when considering the temperature field distribution and temperature effect of the curved beam.
3. Finite Element Verification
The thin-walled closed curved box girder is taken as an example, and the sectional shape is shown in Figure 2. The temperature load here is the temperature gradient load specified in 《Code for Design of Concrete Structures of Railway Bridges and Culverts》(TB 10092-2017) [24], and , whose unit is °C. Section characteristics and data related to calculation are shown in Table 1.
3.1. Single-Span Statically Indeterminate Simply Supported Curved Beam
Solid45 solid element in finite element analysis software Ansys is used to establish a single-span curved beam model, and the boundary conditions are set with reference to the bearing conditions of the structure, as shown in Figure 3. The deflection calculation result diagram is shown in Figure 4, and the calculation results of the finite element and theoretically calculated deflection at some coordinate angles are shown in Table 2 and Figure 5. The relative difference in Table 2 is the comparison result between the theoretical value and finite element calculation result (average value).
It can be seen from Table 2 that when calculating the vertical deflection of a single-span statically indeterminate simply supported curved beam, the calculated results with and without warpage are smaller than those calculated by the finite element method, the difference is about 5%, and the results with warpage are slightly more accurate than those without warpage.
3.2. Multispan Continuous Curved Beam
Taking the deformation equation of a single-span statically indeterminate simply supported curved beam under temperature effect as the basic structure, the multispan continuous curved beam is solved according to the force method equation. In this structure, the generalized forces are the bending moment at the intersection point of each span and the warping double moments. Generalized displacement is the relative rotation angle and warpage at the intersection point, and the magnitudes are, respectively,
Therefore, under the action of temperature load,
If warping is not considered, the generalized displacement has only one formula, which is .
And the force method equation of the n + 1 span continuous curved beam iswhere , , , and are n × n order stiffness matrices, and each item in the matrix is given by [4].
For three-span continuous curved beams, the corresponding parameters are the same as those in Figure 2 and Table 1. Under the temperature load specified in Figure 2, the calculation results of the bearing torque reaction of the statically indeterminate structure are shown in Table 3, and the negative sign indicating direction is omitted from the calculation results of torque in the table.
A three-span continuous beam model is established with Solid45 solid element in finite element analysis software Ansys, and the boundary conditions are arranged with reference to the bearing conditions of the structure, as shown in Figure 6. Under the temperature gradient load specified in Figure 2, the calculation results of finite element deflection and bearing torque reaction are shown in Figure 7, and the comparison between theoretical calculation results and finite element calculation results of each bearing torque reaction is shown in Table 3.
Comparing the calculation results of the reaction force of each support of the statically indeterminate structure, it is found that the relative difference between the calculation results considering warpage and the finite element calculation results is about 5%, and the calculation results have sufficient accuracy, while the calculation results ignoring warpage have larger errors, and some items exceed 10%. There is a big difference between the calculation results considering warpage and ignoring warpage, and the difference between each redundant force and bearing torque reaction force exceeds 5%.
4. Parameter Analysis
The radius , cross-section bottom length , height , and flange length (for example, the cross section shown in Figure 2 is = 4 m, = 2 m, and = 1 m) are taken as parameters for parameter analysis. The theoretical calculation results of the three-span continuous curved beam with different parameters are compared with the finite element calculation results. For the convenience of analysis, the bearing torque reaction near the middle of the span is taken as the comparison item.
4.1. Radius r
On the premise of keeping the cross-section shape, temperature gradient distribution, and other parameters unchanged, three results and relative differences of the bearing torque reaction near the midspan with radius r of 60 m, 80 m, 100 m, 20 m, and 140 m are obtained, as shown in Figure 8. It can be seen from Figure 8 that, with the decrease of radius, the bearing reaction force increases. However, when the radius is different, the relative difference between the calculated results of the torque reaction of the three bearings is basically unchanged. The relative difference between the calculation results considering warpage and ignoring warpage is −5.5%, and that between the calculation results considering warpage and the finite element calculation results is 5.8%. It can be seen that the radius has little influence on the relative difference among the three calculation results. For other sections with different sizes and shapes, the same rules are also observed. Therefore, after parameter analysis, the radius r = 100 m is set.
4.2. Cross-Section Bottom Edge Length l1 and Height l2
On the premise of keeping the temperature gradient distribution and radius unchanged, three results and relative differences of the bearing torque reaction near the midspan are obtained when the length l1 of the cross-section base is 3 m, 4 m, 5 m, and 6 m and the height l2 is 1 m, 1.5 m, 2 m, 2.5 m, and 3 m, as shown in Table 4. See Figure 9 for the relationship between the relative difference between considering warpage, ignoring warpage, and finite element results and the warping torsion attenuation coefficient (), which is an important index to measure the warping torsion geometric characteristics of thin-walled curved beams, so this coefficient is used for analysis. It can be seen from Table 4 and Figure 9 that when the length l1 and height l2 of the bottom edge of the cross section change, the relative difference between the calculated result considering the warping effect and the finite element calculation result is about 5%, which has good accuracy. Under the same parameters, the relative difference between the calculated result without considering warping and the finite element calculation result is higher than the former, and the relative difference of multiple calculation results is over 10% and sometimes reaches 22%. However, there is no obvious positive or negative correlation between the two relative differences and the warping torsional attenuation coefficient.
4.3. Flange Length l3
On the premise of keeping temperature gradient distribution, radius, and other parameters unchanged, taking the three-span curved beam structure with cross sections of l1 = 4 m, l2 = 2 m, and flange length l3 of 0.5, 0.75, 1, 1.25, and 1.5 m as an example, three results and relative differences of the bearing torque reaction near the midspan are calculated, which are considered warping, not considered warping, and finite element calculation, as shown in Figure 10. It can be seen from Figure 10 that when the flange length l3 is changed, the relative difference between the calculation result considering the warping effect and the finite element calculation result is between 2% and 10%, which is larger than the previous analysis. When the flange is long, the theoretical solution proposed in this paper is not accurate enough. The reason may be that the shear lag effect is not considered, and the deformation assumption in this paper is inconsistent with the real deformation distribution of the flange. It is necessary to apply a more precise deformation assumption model to consider the calculation of the temperature effect of the curved beam structure in this case.
5. Conclusion
According to the basic differential equation of the curved beam and the principle of minimum energy, this study takes the primary torsional statically indeterminate simply supported curved beam as the basic structure and, according to whether the warping effect of the thin-walled closed curved box girder is considered, puts forward a simple analytical calculation method for the temperature deformation of the plane curved beam. Comparing the analytical calculation results with the finite element calculation results, the following conclusions are obtained:(1)When calculating the vertical deflection of a single-span statically indeterminate simply supported curved beam, the difference between the calculation results provided by this study with and without warping is small, both of which are about 5%, and the calculation results with warping are slightly more accurate than those without warping. For multispan continuous curved beams, the relative difference between analytical calculation results and finite element calculation results considering warping is about 5% in the analysis, and the calculation results have sufficient accuracy. There is a big error in the calculation results when the warpage is ignored, and the relative difference of multiple calculation results is over 10%. Therefore, considering the warping effect is helpful to calculate the temperature effect of the multispan continuous curved beam more accurately.(2)By analyzing the parameters of curved beams with different cross sections and different radii, it is found that the relative difference between the theoretical calculation results and the finite element results considering the warping effect is greatly influenced by the cross-section shape, but less influenced by the radius size. However, there is no obvious positive or negative correlation between the relative difference and the warping torsional attenuation coefficient of the thin-walled box section.(3)When the flange length is taken as the parameter, the relative difference between the calculation result considering the warping effect and the finite element calculation result reaches 10% when the flange is long. Therefore, more detailed deformation assumptions, such as shear lag deformation assumptions, are needed to more accurately consider the temperature effect calculation of the thin-walled curved box girder when the flange is long.(4)Through theoretical calculation, it is found that, for the thin-walled closed curved box girder, the transverse temperature gradient distribution in the plane will affect the vertical deflection, torsion angle, and warpage out of plane. This property is different from that of the straight beam.
Data Availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the Shandong Transportation Science and Technology Program (2013A08-03 and 2021B85).