#### Abstract

Topology optimization has developed rapidly in the past three decades; as a creative and efficient optimization technique, it has been applied in engineering fields of aerospace and mechanical. However, there are a few attempts in bridge form design. In this paper, the parametric level set method is utilized to solve the form finding of arch bridges. The optimization model for minimizing the structural compliance under the volume constraint is built. Three numerical examples of form finding of arch bridges are studied. Results show that the optimal structures which have well-distributed stress and smooth force transmission are almost identical with the actual forms of arch bridges. The optimal forms can be treated as alternatives in the preliminary design stage, and topology optimization has a bright prospect in form finding of arch bridges.

#### 1. Introduction

In conventional process of engineering design, there is blindness that the initial form of structure is achieved according to the knowledge and experience of engineers. Topology optimization refers to achieving optimal distribution of a material in a given design domain and also obtaining the optimal mechanical performance under specified constraints simultaneously. It can provide a novel and reasonable design in the preliminary design stage. Based on these advantages, topology optimization is widely applied in different fields of engineering, for instance, aerospace, mechanical, electromagnetic, and electrochemical [1–4].

Topology optimization has been developed rapidly in the past three decades. Numerous topology optimization methods including homogenization method (HM) [5, 6], density method (DM) [7–9], evolutionary structural optimization (ESO) [10, 11], and level set method (LSM) [12, 13] have been proposed. Compared with the previous three methods, LSM has a smooth boundary which contains a lot of geometry information, and it is easier to convert the optimal configuration into the CAD model for manufacturing. In addition, LSM handles complex boundary change easily. LSM refers to using of the zero isosurfaces of the level set functions to represent the structural boundaries. In conventional LSM, the evolution of structural boundaries is tracked by solving a Hamilton–Jacobi partial differential equation (H-J PDE). An upwind scheme was employed by Mei et al. [14] to solve the H-J PDE using the discrete solution. In this method, the accuracy and rate of convergence are affected as they are often too steep or too flat near the boundaries during the evolution of structural boundaries. Periodically, reinitializing the signed distance function is needed to alleviate numerical instability. And the optimal result generally depends on the initial topology because the new holes cannot be created in the design domain in the most conventional LSM.

Parametric level set method (PLSM) solves these problems effectively without directly solving the H-J PDE. In this method, the implicit level set function is interpolated by radial basis functions (RBFs), the H-J PDE is converted into a series of ordinary differential equations which are easier to handle, and the propagation of boundaries is driven by updating interpolated coefficients using mathematical programming methods. The PLSM not only inherits the favorable features of the conventional LSM but also avoids the difficulties in solving the HJ-PDE. Wang et al. [15] firstly introduced globally supported radial basis functions (GSRBFs) into topology optimization of minimizing the structural compliance under the volume constraints. Luo et al. [16] applied the compactly supported radial basis functions (CSRBFs) [17, 18] into topology optimization of compliant mechanisms and established the optimization model with a nonconvex objective function and specified constraints. The evolution of the boundary was the process of updating the interpolated coefficients using the method of moving asymptotes (MMA) [19]. Luo et al. [20, 21] also proposed an effective method for structural stiffness design which combines the CSRBFs with well-established optimization algorithms, such as MMA and OC.

Up to date, topology optimization has been widely used in fields of aerospace, mechanical, etc. However, it has seen little use in form finding of bridge. Chen et al. [22] obtained the optimal form of arch bridge using BESO and discussed the optimization results with different rise-to-span ratios. Felicetti et al. [23] applied BESO to bridge design in consideration of different requirements, for instance, support types, constraint types, and selections of span. However, this method has the numerical problem of checkerboard. Furthermore, the process of converting the optimized result to the CAD model needs a lot of postprocessing efforts.

In this paper, PLSM which has a lot of geometry information is applied into form finding of arch bridges. The optimization model of minimizing the structural compliance under the volume constraint is proposed. And three examples of arch bridges are studied. The paper is composed as follows. In Section 2, level set representation and parameterization using CSRBF are presented. After that, the topology optimization model is formulated. The sensitivity analysis and the flowchart of optimization are given in Section 3. Section 4 discusses the results of three examples in form finding of arch bridges. Finally, some conclusions are presented in Section 5.

#### 2. Parametric Level Set Method

##### 2.1. Level Set Representation

In the level set method, the structural boundaries are represented by the zero-level set of a higher dimensional level set function, which is usually defined as follows:where is the level set function, is the design domain, is the positive value of the level set function, and denotes the zero-level set.

A pseudotime is introduced to manage the evolution of structural boundaries; differentiating both the sides of with respect to , we can obtain thatwhere is the normal velocity of the level set function.

##### 2.2. Parameterization Using CSRBF

The PLSM with compactly supported radial basis function (CSRBF) [20, 21] which has sufficiently smooth and accurate to approximate the level set function is applied into topology optimization of arch bridge. The level set function can be expressed aswhere is the total number of computational points, denotes the CSRBF at the th point, and is the interpolated coefficient at the th point.

In this work, CSRBF with C2 continuity [17, 18] which has strictly positive definiteness and the sparseness of collection matrices is employed to approximate the level set function, and it can be expressed aswhere is the radius of support, and it is usually defined in a 2D Euclidean space aswhere is the coordinates of the th computational point and is the radius of support domain of the computational point , which is defined as .

By substituting Equation (3) into Equation (2), the HJ-PDE in Equation (2) can be rewritten as follows:

The normal velocity of the level set function can be calculated bywhere .

Equation (7) shows that is obtained by the variation of interpolated coefficients, and the shape change problem is transformed to a parameterization problem.

#### 3. Topology Optimization for Form Finding of Arch Bridges

##### 3.1. Optimization Model Using PLSM

In this optimization model, the interpolated coefficients of CSRBFs are adopted as design variables, optimization model of minimizing the structural compliance under the volume constraint is proposed, and it can be mathematically expressed as follows:where denotes the objective function, is the displacement field, and is the virtual temperature field. is the strain, is the volume of design domain, is the volume of structure, denotes the volume constraint, and are the lower and the upper bounds of the design variables, and is the elastic stiffness which can be calculated as follows:where is the material elastic stiffness and is the Heaviside function. In order to avoid regenerating the element mesh, we employ a smoothed approximation of , which is defined aswhere and are small positive numbers to avoid singularity of stiffness matrix and is the width of the numerical approximation.

is the weak form of the state equation, and the energy bilinear form and the load linear form can be expressed as follows:where is the body force and is the boundary traction.

##### 3.2. Sensitivity Analysis

The optimization problem in Equation (8) can be solved by sensitivity-based optimization algorithms. Thus, the derivatives of the objective function and volume constraints with respect to the design variables are required. The Lagrangian method is employed to transform the constrained optimization problem into an unconstrained problem, which can be expressed aswhere is the Lagrangian multiplier.

The material derivative is used to obtain the shape derivative of . For more detailed theoretical proofs, the readers can refer to [20, 21]. The shape derivative of can be obtained byandwhere is the mean curvature.

The derivative of Lagrangian with respect to can also be obtained by the chain rule as follows:

By comparing the corresponding terms in Equations (14) and (16), the derivatives of the objective function and the volume constraint with respect to the interpolated coefficients can be obtained by

When Equations (17) and (18) are obtained, the optimization problem in Equation (8) can be solved by sensitivity-based optimization algorithms.

##### 3.3. Optimization Procedure

In optimization problem, with the sensitivities of objective function and volume constraint with respect to interpolated coefficients obtained, the sensitivity-based optimization algorithms such as optimality criteria and mathematical programming can be applied straightforwardly. Then, the shape and topology optimization problem converts into the parameter optimization problem. In this work, the method of moving asymptotes (MMA) which is a widely used mathematical programming algorithm is employed to solve the optimization problem. The process of optimization is shown in Figure 1, and the corresponding introduction is presented as follows: Step 1: initialize the design domain of arch bridge. Define the load, boundary conditions, and set the volume constraint. Step 2: FE analysis and calculate the real and virtual displacement fields. Step 3: calculate the derivatives of objective function and the volume constraint with respect to the interpolated coefficients in Equations (17) and (18). Step 4: calculate the values of objective function and the volume constraint. Update the interpolated coefficients with MMA optimization algorithm, and then, the update level set function. Step 5: check whether the convergence condition is shown in Equation (19). If the terminal condition is met, stop, and the optimal structure is found. Otherwise, go back to Step 2. The terminal condition is defined aswhere denotes the value of the objective function at th iteration, and is a small positive value.

#### 4. Numerical Examples

In this section, form finding of three kinds of arch bridges is presented to illustrate the effectiveness of the proposed method. The conventional level set method (CLSM), BESO, and SIMP are applied for comparison. In these examples, the computational points in PLSM are corresponding to the finite element nodes and four-node plane stress elements are used for finite element analysis (FEA). The parameters in Equation (8) are and , where denotes the initial values of interpolation coefficients, and the support radius for computational points is set to be . , and are the parameters in Equations (10) and (11). The parameters in BESO are , and , while in SIMP, . The Lenovo personal computer equipped with 1 Intel^{@} Core^{TM} i5-3230M processor with a clock speed of 2.6 GHz and 4.00 GB RAM is employed to calculate central processing unit (CPU) times.

*Example 1. *In this deck arch bridge design example, the admissible design domain, boundary conditions, and load conditions are shown in Figure 2. The design domain is modeled into elements, a uniform load is applied in the top of the design domain, and the prescribed volume fraction is .

The optimal layouts obtained by PLSM, CLSM, BESO, and SIMP are shown in Figures 3(a)–3(d), respectively, and it can be seen that the optimal layouts of different methods are almost identical and they are similar to the actual layout of deck arch bridge, which shows the effectiveness of PLSM. However, they also have some differences. The optimal layouts of PLSM and CLSM in Figures 3(a) and 3(b) have smooth boundaries, which are easier to transform into CAD model for bridge design. Figure 3(c) shows the optimal layout of BESO with the zigzag boundary. It is hard for follow-up design. Figure 3(d) shows the optimal layout of SIMP. There are a lot of gray elements, and the bars have split forks in the junction point of the bars and the bridge deck, which shows the optimal transmission path of force. In addition, the relatively smooth boundary can be produced by BESO and SIMP when the design space is modeled into a large number of finite elements. In this way, the computational efficiency is greatly reduced.

Figure 4 gives the iteration history of the objective function and several intermediate results. The value of the objective function is 25.2 after 63 iterations. Corresponding to the PLSM topological designs in Figure 4, the level set surfaces at different iterations are displayed in Figure 5. Table 1 provides the results of different methods in form finding of deck arch bridge for comparison. It can be found that the compliance of PLSM is almost the same with that of other methods. And BESO produces the lowest compliance of 24.6 as it utilizes discrete variables. The PLSM costs the shortest computation time of 55.3 s which shows the high computational efficiency.

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

* Example 2. *In this half-through arch bridge design example, the admissible design domain, boundary conditions, and load conditions are shown in Figure 6. The design domain is modeled into elements, a uniform load is applied in the height of , and the prescribed volume fraction is .

Figure 7 shows the optimal layouts obtained by different optimization methods. It can be seen that the optimal layouts of different methods are similar to the actual layout of half-through arch bridge. The optimal configuration in PLSM is almost identical with the optimal configuration in CLSM as they have the same number of bars, which shows the validity of PLSM. Both of them have smooth boundaries. The optimal configurations of BESO and SIMP have more bars than that of PLSM and CLSM. The optimal configuration obtained by BESO in Figure 7(c) has inclined supports under the deck and also has a zigzag boundary. The configuration of SIMP in Figure 7(d) has a lot of gray elements.

The iteration history of the objective function and several intermediate results are shown in Figure 8. The objective function is convergent to 19.9 after 37 iterations. For comparison, the results of different methods in form finding of deck arch bridge are presented in Table 2. It can be found that the compliance of PLSM is less than that of CLSM and SIMP, and BESO produces the lowest compliance of 18.6. The PLSM costs the shortest computation time of 100.9 s which is much less than other methods in calculation.

**(a)**

**(b)**

**(c)**

**(d)**

* Example 3. *In this through arch bridge design example, the admissible design domain, boundary conditions, and load conditions are shown in Figure 9. The design domain is modeled into elements, a uniform load is applied in the bottom of the design domain, and the prescribed volume fraction is .

The optimal results of various topology optimization methods are illustrated in Figures 10(a)–10(d), respectively. It can be seen that the optimal layouts of different methods are similar to the actual layout of through arch bridge. The optimal configurations in PLSM and CLSM have smooth boundaries, and the optimal configuration in PLSM is almost identical with the optimal configuration in CLSM. So, the effectiveness of the PLSM can be verified. The optimal configurations of BESO and SIMP have more bars than that of PLSM and CLSM. The configuration of BESO in Figure 10(c) has a zigzag boundary. The configuration of SIMP in Figure 10(c) has a lot of gray elements. And, it should be noticed that the bars have split forks in the junction point of the bars and the bridge deck.

Figure 11 gives the iteration history of the objective function and several intermediate results. The value of objective function is 24.3 after 80 iterations. Table 3 summarizes the structural compliances and computation time of different methods in form finding of through arch bridge. It can be found that the compliance of PLSM is almost the same with that of other methods, and BESO produces the lowest compliance of 24.6. The PLSM costs the shortest computation time of 186.2 s. So, the accuracy and efficiency of PLSM are proved.

**(a)**

**(b)**

**(c)**

**(d)**

#### 5. Conclusions

In this study, a form finding method for arch bridge is developed using the PLSM-based topology optimization, the solution procedure of the topology optimization method is described, and three numerical examples are employed to illustrate the feasibility and validity of PLSM in form finding of arch bridges. Numerical examples indicate that the configuration obtained by PLSM has smooth boundary. It is easier to convert the configuration to CAD model for manufacturing. In addition, the computational efficiency has been greatly improved.

This method also can be extended to form finding of arch bridges with different rise-to-span, main tower of cable stayed bridge, and main cable of suspension bridge. Due to the effectiveness in form finding of arch bridges, it is expected that this method will be applied in bridge design in the future.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This work was supported by the National Natural Science Foundation of China (11601037), Shanxi Province Natural Science Fund (2018JQ1027), and the Research Funds for the Central Universities of China (310812171002).