Design Optimization of Concrete Aqueduct Structure considering Temperature Effects

An aqueduct is a water conveyance structure that enables channel flow across canals, valleys, depressions, roads, and other structures. *e optimal structural selection of the aqueduct is particularly important to ensure engineering quality and optimize project investment. To optimize the design of an aqueduct structure, this study established a mathematical model based on the three-dimensional finite element method that considers the temperature field and structural stress field coupling among its design parameters. *e model was used to optimize and design the main wall thickness and tie spacings of the aqueduct structure. *e Caohe aqueduct was considered as an example for the proposed design. *e influences of temperature-induced stress on the reinforced concrete structure of the aqueduct in winter and summer were investigated based on the actual engineering conditions of the structure, and the corresponding structural optimization was obtained. *e results showed that the optimized aqueduct can offset temperature and structural stresses, thus reducing the amount of material required. *e maximum generated stress was also lower than that of the original design. Furthermore, this study is expected to provide guidance for similar engineering construction projects.


Introduction
An aqueduct is a water-conveying building that is significantly influenced by water. In particular, the expansion of water during freezing in winter will affect the aqueduct, and floods during summer will increase the water level of the aqueduct. Temperature stresses also influence aqueducts. us, an aqueduct structure that uses less material and resists temperature stress is of great significance to water transfer projects.
ere are three main types of concrete structures: plain concrete, which has no steel reinforcements, a high compressive strength, and a low tensile strength; reinforced concrete, which can withstand significant tensile stresses due to steel reinforcements; and prestressed concrete, where the concrete is squeezed due to high-strength prestressed steel reinforcement, enabling the resulting precompressive stresses to counteract external loads and increasing the crack resistance of structural members. is study focuses on prestressed concrete aqueducts, which are capable of withstanding considerable tensile and compressive stresses with reduced displacements, thereby preventing cracking [1].
Many reported studies have investigated the optimization of the aqueduct structure. Based on the application of the genetic algorithm to the pi (π)-type arch aqueduct structure optimization design, Xia et al. [2] verified the feasibility and reliability of the genetic algorithm. Zheng et al. [3] proposed an improved genetic algorithm and applied it to optimize the design of aqueduct structures, revealing improved performance and optimization results. Based on the improved ant colony algorithm, Zhou et al. [4] developed a mixed discrete variable program for optimizing aqueduct structural design, and they demonstrated that this algorithm is reliable and has a strong global convergence ability. Liu [5] attained an economical and reasonable structural form for an aqueduct and its dimensions by optimizing its longitudinal slope, span cross-sectional form, and structural design. Furthermore, this scheme was applied to the Eastern Water Supply Works of Shenzhen city, and the economic benefits were validated. Existing reports in the literature do not consider the role of temperature stress during the optimization of aqueduct structures. In particular, for aqueducts spanning long distances and subjected to varying climate conditions, optimizing the structure while accounting for temperature stresses can help reduce the amount of concrete required and the associated costs. In this study, the subproblem method was used for the first optimization calculation of the aqueduct, and the scanning method was used for the second optimization analysis. e results will provide a valuable reference for similar projects [3]. Currently, breakthroughs in the finite element method (FEM) have enabled its widespread application in various fields, including nonlinear analyses of embedded viscoelastic micro composite beams with initial geometrical imperfections [6]. e present study uses the finite element method (FEM) to simulate an aqueduct under a variety of load conditions. Simulations were then performed to optimize the proposed aqueduct design for summer and winter conditions. is study makes a significant contribution to the literature because the optimized structure considering temperature enabled reducing the amount of material required for aqueduct construction by 5-6%.

Engineering Profiles.
In the simulation, the main riverbed slot side uses a 30-meter-long three-way prestressed reinforced concrete structure with multiple sidewalls (Figure 1). e single slot section is 6.0 × 5.4 m 2 . e thickness of the sidewall is 0.6 m, and the pavement plate at the top of the structure is 2 m wide. e thickness of the middle wall is 0.7 m, and the pavement plate at the top of the structure is 2.7 m wide. Tie bars are set between the two adjacent walls, and each bar is 2.5 m long, 0.3 m wide, and 0.4 m high. e bars are spaced 2.5 m apart. e two sidewalls of the water flume part of the aqueduct (without the support section) are equipped with side ribs, and each rib is 0.5 m wide and 0.7 m high. e bottom rib is placed at the bottom of the aqueduct. e ribs near the piers are each 0.7 m wide and 2.0 m high, and the other ribs are each 0.5 m wide and 0.9 m high. Each cross end between two slots has a reserved 0.55-meter-long postcast strip to allow the workspace for prestressed tensioning. e piers of the structure are of the hollow gravity type, and the foundation is an end-bearing pile; the platform is 26.8 × 8.5 × 2.0 m 3 (length × width × height). Twelve root end-bearing piles are arranged in a double row under the platform.
e prestressed water flume of the aqueduct (without the support section) mainly adopts straight and curved structures; its diameter is 15.2 mm. ere are 980 rectilinear steel stranded wires that are of 10 different kinds with 150 channels and a total length of 22,771 cm ( Figure 2). Furthermore, there are 328 curved steel stranded wires (Figure 3), which are of nine different kinds of components with 26 channels and a total length of 25,417.1 cm (Figure 4).

Finite Element Models.
Element selection for the water flume of the aqueduct (without the support section) is the main analytical object of the aqueduct structure. e second object of importance is the steel stranded wire. In this study, the water flume of the aqueduct structure (without support) was regarded as a hydraulic concrete structure for calculations; the Caohe aqueduct is a kind of large-span, highcarrying-capacity aqueduct, in which the simulation of the prestressed steel strands must be considered. is also involves considering the influence of the steel elements in the concrete structure. Two types of elements in the Ansys simulation were introduced as follows [7]. e structural analysis uses Link8 (three-dimensional pole element) elements to simulate the steel strands and Solid45 (three-dimensional element) elements to simulate the concrete components in the finite element software Ansys [8]; this analysis is converted to a thermal analysis after applying the gravity load, water load, wind load, and other loads. e steel strands are changed to the Link33 (three-dimensional pole element) elements for this simulation, and the concrete is changed into the Solid70 element (three-dimensional body element). e temperature is applied in a steady-state thermal convection mode using Ansys to apply the convection temperature from outside the solid. In order to obtain the internal temperature of the solid, a thermal analysis of the aqueduct is performed using the temperature field results, and the node temperature of the thermal analysis is obtained as a body load; when the Ansys thermal analysis is converted into a structural analysis, the body load is read as a load. Finally, the temperature field and stress field of the structure are simultaneously simulated to investigate the stress and strain of the entire structure.

Establishment of the Aqueduct Model.
e structure of the aqueduct is regular and symmetrical, and the design requirement is significant, so it is appropriate to use FEM meshes. e FEM model is shown in Figures 5-8 below. Figure 6 is an FEM model of the solid aqueduct structure in Figure 5. After applying the mesh, the model contains 50,240 elements and 62,273 nodes. Figure 7 is the FEM model of the steel strand Cache aqueduct.
e steel strand simulation theory shows that curved steel parts in the Caohe aqueduct have complicated conditions. An appropriate constraint equation is used to model the real physical reinforcement. us, concrete elements and reinforced node groups were used, and the results are shown in Figure 8. Figure 8 illustrates that the constraint equation is applied around the steel strands, and the method does not change the node locations; this calculation method is much simpler than the node-coupled method, and the calculation accuracy is higher. erefore, this method is better for solving the structural force line problem of a complex shape. e modeling problem of the steel-reinforced concrete structure for complex calculation was thus solved. e water flume part of the stress calculation for the aqueduct structure (without support) considered the following eight combinations. According to the design data, the main aqueduct loads are weight, water, wind, static ice pressure, prestressed force, and temperature action. e combinations of these conditions are summarized in Table 1.

Load Distribution on Item.
e listed loads were all standard values according to the design specification for hydraulic concrete structures [9]; when calculating the bearing capacity, these values should be multiplied by the corresponding subcoefficients of the load components. e self-weight of the structure as a permanent load corresponds to a subcoefficient of 1.05. e subcoefficients for variable loads include a subcoefficients of 1.1 for water pressure, 1.1 for ice pressure, 1.3 for wind pressure, and 1.1  When water flows through the aqueduct in summer (according to the working conditions of uniform temperature warming), the temperature calculation was based on the outside surface temperature of the aqueduct (41°C), the external nighttime temperature was 35°C, and the tank water temperature was 28°C. When the water flows through the aqueduct in winter (according to the working conditions of uniform temperature cooling), the temperature calculations were based on an aqueduct external surface temperature of −10°C and a tank water temperature of 4°C. e ice thickness was 0.4 m; this impacts the unit length of the building for the static ice pressure standard value of

Structure Assumption.
In the x-direction, the aqueduct was considered a simply supported beam in the flow direction, and the cross direction is positive. In the y-direction, the aqueduct was considered a cantilever beam in the vertical flow direction, and the downward direction is positive. In the z-direction, the aqueduct was considered a three-span continuous beam in the transverse view, and the left side of the front view is considered positive.

Temperature Calculation
Assumption. e temperature action is steady-state thermal convection modeled using  is kind of calculation method is simple, and the accuracy meets the previously reported requirements [3].   (3) Reinforced Concrete. A comprehensive elastic modulus was adopted: E c (concrete) + k * E c (elasticity modulus of rebar multiplied by its volume ratio in concrete) � 3.71 × 104 N/mm 2 .

Boundary Conditions.
According to the conditions of the support constraints (horizontal arrows represent the bearing movement in the vertical flow direction, and vertical arrows represent the bearing movement in the direction of flow), the boundary conditions of the water flume of the aqueduct (without the support) are three constraints in the x-, y-, and z-directions in the end-circle bearing of the water transfer part of the aqueduct (without the support part), and the remaining three bearings were fixed in the vertical ground and flow directions. e corresponding circle bearing at the other end was fixed in the vertical ground and vertical flow directions, and the remaining three bearings were fixed in the vertical direction ( Figure 9).

Applied Loads
(1) Self-weight for C50 concrete density: c � 25.0 kN/m 3 as the applied inertia load. (5) e temperature was chosen from among conduction, convection, or radiation heat transfer. is can be used to transform analysis element types into the structural analysis and to apply the winter and summer temperature actions.
e above loads were considered with subcoefficients.
(6) Prestress is applied by the initial strain method.

Calculation Analysis.
By using the Ansys software, boundary constraints and loads for the working condition were applied to calculate the stress and strain. e component stress and maximum deformation values of the control section were found on the stress cloud map. e calculation results are shown in Table 2.

Stress Analysis of Key
Parts. e calculations show that most parts of the Caohe aqueduct structure were in a state of compression under various conditions and had a partial section in a smaller tensile state. e stress ranges for all parts are listed in Table 3.
Stress Results Summary. Stress in the water flume of the aqueduct (without support) was analyzed under eight working conditions; working conditions 1, 3, and 5 are of the same regularity, and working conditions 2, 4, and 6 are similar as well. e stresses in the water flume of the aqueduct (without support) are satisfactory for all parts; the tensile stress is less than that specified in the design specification for hydraulic concrete structures (DL/T5057-2009) for the allowed range; the compressive stress is not strong, which means that the structural design is reasonable [11]. e control condition is a short-term load condition. With a full tank of water pressure and ice load in the tie bar, the vertical stress in the direction of flow is 1.26 MPa [12].

Aqueduct Structure Optimization Design
e optimization mainly aims to meet aqueduct safety conditions and conserve engineering materials [13]. e objective function is the minimum total volume of the solid model, and the constraints are that the optimization results are no more than the maximum stress value of the original design scheme, not exceeding the maximum displacement of the original plan [14]. e most representative aqueduct section was selected for optimization, and the subproblem method was chosen for the first optimization, with the scanning method being chosen for the second time [15,16]. e aqueduct section size parameters are shown in Figure 10. e optimization mathematical model was set up to consider the following objective function and constraint conditions: where svolume is the total volume. DMIN is the smallest value of displacement, SMAX is the maximum stress, and SMIN is the minimum stress. e variables were not optimized, which used the original design values [17]. e water load was imposed in accordance with the full tank water depth [18]. e optimized design variable values are shown in Table 4. is optimization mainly satisfies three constraints: first, the maximum absolute value of displacement does not exceed the original design value of 4.78e-4 m; second, the maximum tensile stress value does not exceed the original design value of 4.53 MPa; and third, the maximum compressive stress value does not exceed the original design value of 19.47 MPa. In case the abovementioned conditions are satisfied, a design solution that reduces the volume can be obtained. e volume of the aqueduct is calculated by stretching its cross section; therefore, in this study, the crosssectional area of the optimized aqueduct is used instead of the optimized aqueduct volume. In other words, to calculate the optimal solution, it is necessary to ensure that the maximum values of displacement and strain are less than the maximum values in the original design. Subsequently, the minimum cross-sectional area of the aqueduct needs to be determined. Reducing the cross-sectional area of the aqueduct is equivalent to reducing the volume of the aqueduct. e summer and winter project body temperature distributions are shown in Figures 11 and 12, respectively. e top part of the aqueduct temperature in contact with the water remained the same as the water temperature, and the temperature change was mainly in the sidewalls. e volume of the original design project is 73.55, and the summer optimization solution is 69.82; i.e., it reduced by 3.73 or about 5.07%. us, this design can save about 5.07% of the C50 concrete dosage. e winter optimization solution is 68.95, reduced by 4.6 or about 6.25%, so this design can save about 6.25% of the C50 concrete dosage. For different regional climate conditions, the summer project is more suitable for regions with a long summer, and the winter project is more suitable for those with a long winter [19,20].

Summary
e summer optimization solution is used for areas experiencing long periods of hot weather, whereas the winter optimization solution is used for areas experiencing long periods of cold weather. For the summer optimization solution, the left and right sides of the aqueduct are considered to have different temperatures, considering that only one side faces the sun at most times. e high specific heat capacity of water is not easily affected by temperature; hence, water temperature is mainly regarded as a constant temperature in the calculation. e results show that the temperature of the inner side of the aqueduct in contact with water is basically the same as the water temperature, which is approximately 28°C; on moving from the inside of the aqueduct to its outside, an increasing temperature gradient is observed. In the winter optimization project, the temperature of the left and right sides is considered to be −10°C because the influence of the sun is limited during winter. e water temperature inside the aqueduct is 4°C, and there is no ice. e calculations show that the point of contact between water and the aqueduct is maintained at the same temperature of 4°C, and a decreasing temperature gradient is observed from the inside of the aqueduct to its outside.
Winter optimization project used 1.18% less material than the summer project. Considering the heat bilges cold shrink phenomenon of the material, the aqueduct volume in summer is larger than that in winter, so the temperature stress is larger than that in the winter. us, the summer structure needs to bear more temperature stress, and it requires more material. is optimization project conforms to the actual engineering situation. For areas with different temperature conditions, optimization considers the    Mathematical Problems in Engineering temperature that can conserve more materials.
e optimization in the high-temperature region uses a high-temperature action. In a region of low temperature, the optimization runs under a low-temperature action. is work mainly studied the optimization of the maximum water load under different temperature conditions because the water load widely influences the structure of the maximum stress and displacement in the optimization is no more than the original design project, and the concrete dosage was optimized to be 5.07% and 6.25% less that the original material usage, respectively. Under the summer and winter temperature conditions, this study can provide theoretical references for similar projects.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.