#### 1. Introduction

Hence, this paper aims to explore the thermoelastic behavior and in-plane instability of fixed arches under a linear temperature gradient field and a uniformly distributed radial load by using the method of effective centroid and principle of virtual work. The analytical solutions of the in-plane displacements, internal forces, and in-plane instability load of fixed arches are obtained while the influences of the gradient temperature differential, the included angle and the slenderness ratio of arches on the in-plane stability analysis of arches are considered. All analytical solutions are further verified by FE results.

#### 2. Model and Assumptions

The model of the fixed circular arch is shown in Figures 1(a) and 1(b). An axis system is applied to the circular arch. The origin of the axis system is the center of cross section in the arch crown, the axis is the vertical axis of cross section, and is tangent with the arch axis which is coincident with the arch axis, and pointing to the right side of arch, the axis is the horizontal axis of cross section of the arch, where R, S, L, f, ϕ, and 2Φ are the radius, length, span, rise, angular coordinates, and including angle of the circular arch, respectively. b and h are the width and height of the cross section, respectively. In addition, Tu and Tb are the cross-sectional temperatures at the top and bottom fiber, respectively; q is the uniformly distributed radial load.

The thermoelastic and in-plane instability analyses of this paper use the following assumptions:(1)There is no correlation between temperature T and thermal expansion α.(2)The deformation state and the temperature state are considered to be time-independent.(3)The deformations of the arch satisfy the theory of Bernoulli–Euler beam.(4)It is assumed that the section size of the arches is much smaller than the length and radius of the arches.(5)The expansion of the cross section is assumed to be small and can be ignored in the analysis.(6)The temperature of environment is assumed to be 20°C, and the temperature gradient is linearly distributed along the spindle axis of the cross section but evenly distributed along the spindle axis and the geometric centroid axis . Therefore, the temperature of arbitrary point of the cross section is not the function of its coordinates and but is the function of its coordinate . The temperature of arbitrary point of the cross section can be given by

#### 3. Thermoelastic Analysis

##### 3.1. Effective Centroid

For a steel arch, the relation between Young’s modulus E and temperature T is defined by [21]where E20 and are Young’s modulus of steel at the temperature 20°C and the reduction coefficient, respectively, and can be expressed as

The temperature distribution of the cross section is plotted in Figure 2(a), and typical changes in Young’s modulus of steel vs. gradient temperature differential are plotted in Figure 2(b). Figure 2(a) shows that the temperature is distributed in a linear gradient across the cross section of the arch. Figure 2(b) shows that Young’s modulus of low-carbon steel decreases obviously as the gradient temperature differential increases.

The rectangular solid section shown in Figure 3 is used by Pi and Bradford [19] to derive the effective centroid. According to theory of Pi and Bradford [19], it is not precise to conduct thermoelastic analysis of the arch if the axial system composed of the geometric spindles and of the cross section and the geometric centroid axis is used. Therefore, it is necessary to define an effective centroid axis os to accuracy and simplify the analysis. The position of the geometric centroid related to the effective axes oxy was given by [19]where is the effective axial stiffness which can be expressed as

The position of the effective centroid o can be expressed as

In addition, the axial deformation caused by thermal expansion at the effective centroid can be expressed aswhere α and To are the coefficient of linear expansion and the temperature of the effective centroid, respectively, and To can be expressed aswhere Tc is the temperature of the geometric centroid.

The difference value yc between geometric centroids and the effective centroids for arches under linear temperature gradient field with the rectangular solid section are shown in Figure 4(a), where the temperatures at the top fiber Tu = 20°C. Figure 4(a) shows that the linear gradient temperature field has a significant influence on the effective centroid of the rectangular section, and yc increases as the gradient temperature differential increases.

The differences between the effective centroid temperature To and the average temperature Tc for arches with a rectangular solid section under linear temperature gradient field are shown in Figure 4(b), where the temperatures at the top fiber Tu = 20°C. Figure 4(b) shows that the linear temperature gradient field has great influence on Tc − To of the rectangular solid section. The value of Tc − To increases as the gradient temperature differential increases.

##### 3.2. Differential Equations of Equilibrium

The in-plane thermoelastic analysis of fixed arches under a linear temperature gradient field and a uniformly distributed radial load is presented in this section. It is known that the axial compressive force and the bending moment of arches under a linear temperature gradient field and a uniformly distributed radial load are variable. In addition, before carrying out the in-plane instability of an arch, it is essential to accurately determine the axial compressive force of the arch. Therefore, the preinstability analysis of the arch is needed. Based on the energy method, the energy differential equation for thermoelastic analysis of the arch under a linear temperature gradient field and a uniformly distributed radial load can be expressed aswhere εss0 is the linear preinstability longitudinal normal strain, which can be expressed as

σss0 is the preinstability longitudinal normal stress, which can be expressed aswhere and are defined as the dimensionless axial displacement and radial displacement, respectively.

Substituting equations (10) and (11) into equation (9) can be obtained as

The initial expression of axial compressive force N and bending moment M are given bywhere is the effective stiffness, which can be expressed as

Integrating by parts of equation (12), the differential equations can be obtained as

##### 3.3. Displacements

The static boundary conditions for fixed archesare required. The dimensionless displacement and for arches under linear temperature gradient field and a uniformly distributed radial load can be solved by substituting equations (13), (14), (18), and (19) into equations (16) and (17):where .

The distributions of the dimensionless displacement and along the arch axis given by equations (20) and (21) are plotted in Figures 5(a) and 5(b) for arches with bottom fiber temperature of cross section Tb = 50°C, 100°C, and 200°C, respectively, where the top fiber temperature of cross section Tu = 20°C, the included angle 2Φ = 60°, and with .

Figures 5(a) and 5(b) show that the dimensionless radial displacement of arches under a linear temperature gradient field and a uniformly distributed radial load are symmetrically distributed along the arch axis, while the dimensionless axial displacement of arches are antisymmetrically distributed along the arch axis. In addition, Figure 5(a) shows that the dimensionless radial displacement increases as the bottom fiber temperature of cross section Tb increases. Figure 5(b) shows that the axial displacements of the left half of the arch increases as the bottom fiber temperature of cross section Tb increases, while the dimensionless axial displacements of the right half of the arch decreases as the bottom fiber temperature of cross section Tb increases. Figure 5(b) also shows that the axial displacements at the crown and two ends are equal to zero.

The distributions of the dimensionless displacement and along the arch axis given by equations (20) and (21) are plotted in Figures 6(a) and 6(b) for arches with the included angle 2Φ = 30°, 60°, and 90°, respectively, where the top fiber temperature of cross section Tu = 20°C, the bottom fiber temperature of cross section Tb = 100°C, and .

Similarly, Figures 6(a) and 6(b) show that the dimensionless radial displacement of arches under a linear temperature gradient field and a uniformly distributed radial load are symmetrically distributed along the arch axis, while the dimensionless axial displacement of arches are antisymmetrically distributed along the arch axis. In addition, Figure 6(a) shows that the radial displacements increase as the included angle 2Φ increases. Figure 6(b) shows that the axial displacements of the left half of the arch increase as the included angle 2Φ increases, while the axial displacements of the right half of the arch decrease as the included angle 2Φ increases.

##### 3.4. Internal Forces

Substituting equations (24)–(26) into equations (14) and (15) leads to the solution of axial compressive force N and bending moment M for arches under a linear temperature gradient field and a uniformly distributed radial load, which can be expressed as

To demonstrate the influences of the linear temperature gradient field on N and M given by equations (23) and (24). The variations of Nc with the included angle 2Φ for arches with various bottom fiber temperatures of cross section Tb (Tb = 20°C, 200°C, and 400°C) are plotted in Figure 7(a), where Nc is the central axial compressive force, and the top fiber temperature of cross section Tu = 20°C, q = 2000 kN/m. In addition, the variations of Mc with the included angle 2Φ for arches having different bottom fiber temperatures of cross section Tb (20°C, 200°C, and 400°C) are plotted in Figure 7(b), where Mc is the central bending moment and the top fiber temperature of cross section Tu = 20°C, q = 2000 kN/m.

Figure 7(a) shows that Nc increases as the bottom fiber temperature of cross section Tb increases. Figure 7(a) also shows that Nc increases as the included angle 2Φ increases initially and, after that, it decreases as the included angle 2Φ increases in case the included angle of arches reaches to a certain value. Figure 7(b) shows that Mc decreases as the bottom fiber temperature of cross section Tb increases. Figure 7(b) also shows that Mc decreases as the included angle 2Φ increases initially and, after that, it increases as the included angle 2Φ increases once included angle of arches attach a certain value.

To demonstrate the influences of the slenderness ratio on N and M given by equations (23) and (24), the variations of Nc with 2Φ for arches with different slenderness ratio () are plotted in Figure 8(a), where the top fiber temperature of cross section Tu = 20°C and the bottom fiber temperature of cross section Tb = 100°C, q = 2000 kN/m. In addition, the variations of Mc with 2Φ for arches with different slenderness ratio () are plotted in Figure 8(b), where the top fiber temperature of cross section Tu = 20°C and the bottom fiber temperature of cross section Tb = 100°C, q = 2000 kN/m.

Figure 8 shows that both the central axial force Nc and the central bending moment Mc increase as the slenderness ratio increases.

#### 4. Instability Analysis

According to the principle of virtual work, using the condition of equation (25) and substituting the nonlinear strain given by equation (26) and the bending moment M given by equation (14) lead to the virtual work statement for instability equilibrium as

Integrating equation (27) into parts results in the differential equation for the dimensionless radial instability displacement aswhere

The general solution of equation (28) is obtained as

Substituting the boundary conditions of fixed arches at into equation (30) leads to a group of four homogeneous equations for C1, C2, C3, and C4 aswhere the coefficient matrix is given by

The existence of nontrivial solutions for C1–4 requires the determinant of the matrix of equation (32) to vanish, resulting in

The antisymmetric solutions obtained from the second factor of equation (33) satisfy the axial inextensibility condition for bifurcation instability of an arch given by equation (25), which can represent bifurcation instability modes of the arch as pointed out in [15]. The lowest antisymmetric solutions of fixed one can be obtained from the second factor of equation (33) aswhere the coefficient η is related to the included angle 2Φ as shown in Figure 9.

In addition, substituting the solution (34) into equation (29) leads to lowest antisymmetric bifurcation instability load for fixed arches as

Substituting the axial force N obtained from equation (23) into equation (35) leads to an equation as

The lower bound and upper bound of the antisymmetric bifurcation instability loads can be given by equation (36) by setting and . In addition, the average antisymmetric bifurcation instability loads can be used and given by averaging equation (23) over the arch length as

Typical changes of the lower and upper bounds of the dimensionless in-plane antisymmetric bifurcation instability loads and the average dimensionless in-plane antisymmetric bifurcation instability loads qcrR/NE2 vs. the included angle 2Φ are plotted in Figure 10. Figure 10 shows that the lower and upper bounds of the critical instability loads are very close to the critical average instability load given by equation (37) for fixed arches.

Typical changes of the dimensionless in-plane antisymmetric bifurcation instability load qcrR/NE2 vs. 2Φ obtained from equation (37) for fixed arches having different bottom fiber temperatures Tb are plotted in Figure 11, where

In Figure 11, arches with rectangular solid section have its dimensions as follows: the width b = 2000 mm and the height h = 500 mm. Besides, Poisson’s ratio ν and Young’s modulus E at T = 20°C of the material are assumed as ν = 0.32 and E20 = 200 GPa, respectively.

Figure 11 shows that the gradient temperature differential influences the in-plane antisymmetric bifurcation instability loads of fixed arches significantly. As expected, the dimensionless instability load qcrR/NE2 decreases as the gradient temperature differential increases, and it is also shown from Figure 11 that as the included angle 2Φ increases, the critical instability load decreases rapidly for fixed shallow arches with 2Φ < 40° and the critical instability load becomes nearly a constant when fixed nonshallow arches 2Φ > 40°.

Typical changes of the dimensionless in-plane antisymmetric bifurcation instability load qcrR/NE2 vs. 2Φ obtained from equation (37) for fixed arches with different slenderness ratios are plotted in Figure 12, where the dimensions of the rectangular solid section and the properties of the materials are the same as those in Figure 11.

Figure 12 shows that the slenderness ratio influences the in-plane antisymmetric bifurcation instability loads of fixed arches significantly. As expected, the dimensionless instability load qcrR/NE2 decreases obviously as the slenderness ratio increases for shallow arches about 2Φ < 60° and decreases slightly as the slenderness ratio increases for nonshallow arches about 2Φ > 60°.

#### 5. Comparisons with FE Results

To verify the validity and accuracy of analytical solutions of axial force and bending moment and antisymmetric bifurcation instability loads for fixed arches under a linear temperature gradient field and a uniformly distributed radial load, the scheme of analytical results comparison is adopted with the FE results simulated by ANSYS. In the FE analysis, a rectangular solid section was used and its dimensions are the width b = 2000 mm and the height h = 500 mm. Besides, Poisson’s ratio ν and Young’s modulus E with T = 20°C of the steel are set as ν = 0.3 and E20 = 200 GPa, respectively. Beam 3 of ANSYS is a linear 2-node beam element in 2-D, which is used to formulate the FE model. Convergence studies were demonstrated to show that 300 beams and 3 elements can get accurate converged results. Hence, the arch is divided into 300 elements. In addition, according to equation (2), Young’s modulus of the steel decreases as temperature increases, and so the relation of Young’s modulus and temperature needs to be set in ANSYS based on equation (2).

The analytical solutions of the axial force and the bending moment obtained from equations (23) and (24) for arches under linear temperature gradient field and a uniformly distributed radial load are compared with the FE results in Figure 13(a) as the variations of Nc with the included angle 2Φ for arches with the slenderness ratio , respectively, and in Figure 13(b) as the variations of Mc with 2Φ for arches having the slenderness ratio , respectively, where Tu = 20°C and Tb = 200°C. Figure 13 shows that the analytical solutions of the axial force and the bending moment agree with the FE results very well.

The analytical solutions of the in-plane antisymmetric bifurcation instability load obtained from equation (37) for arches under a linear temperature gradient field and uniformly distributed radial load are compared with the FE results as the variations of the dimensionless instability load qcrR/NE2 with 2Φ in Figure 14(a) for arches with the slenderness ratio and in Figure 14(b) for arches with the slenderness ratio , where Tu = 20°C and Tb = 200°C.

Above all, the analytical solutions of the critical instability load agree with the FE results very well for fixed shallow arches (2Φ < 90°), and the analytical solutions of the instability load agree with the FE results not excellent for fixed nonshallow arches (2Φ > 90°) because the assumption of the axial inextensibility condition is used in instability analysis. For the instability analysis of shallow arches (2Φ < 90°), the assumption is perfectly satisfied, but for the instability analysis of nonshallow arches (2Φ > 90°), the assumption is not perfectly satisfied, and the partial extensibility of the section of nonshallow arches exists. Hence, the analytical solutions can be used to predict the in-plane antisymmetric bifurcation instability load accurately for shallow fixed arches and approximately for nonshallow fixed arches under a linear temperature gradient field and a uniformly distributed radial load.

#### 6. Conclusion

This paper presented the in-plane antisymmetric bifurcation instability of fixed circular arches under a linear temperature gradient field and a uniformly distributed radial load. To derive the critical instability load, accurate preinstability thermoelastic analyses for solutions of the displacements and the internal forces were conducted. The following was found:(1)The radial displacements of arches under a linear temperature gradient field and a uniformly distributed radial load are symmetrically distributed along the arch axis, while the axial displacements are antisymmetrically distributed along the arch axis. In addition, the radial displacement increases as the bottom fiber temperature of cross section increases. The axial displacements of the left half of the arch increases as the bottom fiber temperature of cross section increases, while the axial displacements of the right half of the arch decreases as the bottom fiber temperature of cross section increases.(2)The central axial force increases as the bottom fiber temperature of cross section Tb increases, while the central bending moment decreases as the bottom fiber temperature of cross section Tb increases.(3)The central axial force increases as the slenderness ratio increases, while the central bending moment increases as the slenderness ratio increases.

In addition, the analytical solutions of the axial force and the critical instability load agreed with the FE results very well. It was indicated that the analytical solutions for the in-plane antisymmetric bifurcation instability load of fixed arches derived in this paper can be used to predict the in-plane antisymmetric bifurcation instability load of both shallow and nonshallow arches under a linear temperature gradient field and a uniformly distributed radial load.

This paper mainly deepens the academic understanding of the in-plane thermoelastic instability of fixed arches under linear temperature gradient field and a uniformly distributed radial load. In engineering practice, arches are usually influenced by the gradient temperature field, so it is worth investigating the in-plane instability of arches considering the effect of the gradient temperature field. Of course, in addition to the linear temperature gradient field, there are also nonlinear temperature gradient fields on hollow sections such as box section and circular pipe section. The in-plane thermal instability of such arches considering the influences of shear deformation is worth to be researched in the future.

#### Data Availability

The data used to support the findings of this study are included within the article.

#### Conflicts of Interest

The authors declare there are no conflicts of interest.

#### Acknowledgments

This work was financially supported by the National Natural Science Foundation of China (no. 51578166) and Technology Planning Project of Guangzhou City (no. 201807010021).