Open Access
ARTICLE
Xudong Jiang^{1,*}, Jiaqi Ma^{1}, Xiaoyan Teng^{2}
1 School of Mechanical and Power Engineering, Harbin University of Science and Technology, Harbin, 150080, China
2 College of Mechanical and Electrical Engineering, Harbin Engineering University, Harbin, 150001, China
* Corresponding Author: Xudong Jiang. Email:
(This article belongs to this Special Issue: Computer Modeling in Ocean Engineering Structure and Mechanical Equipment)
Computer Modeling in Engineering & Sciences 2023, 135(1), 511-526. https://doi.org/10.32604/cmes.2022.022785
Received 26 March 2022; Accepted 24 May 2022; Issue published 29 September 2022
Topology optimization aims to produce design solutions of high performance by finding the optimal structural layout in structural design. It has the computational ability to solve great problems and produce reliable solutions to various engineering problems, therefore creating an immense opening for research in applied mechanics. During the last three decades, many gradient- or heuristic based optimization methods have been extensively investigated, according to [1], including density approach, topological derivatives, level set approach, phase field approach and evolutionary approaches. Among them, the convergent and mesh-independent bi-directional evolutionary structural optimization (BESO) developed by Huang et al. [2], iteratively removing inefficient material in addition to add material to the most demanding places, has become a widely adopted methodology for both academic research and engineering application [3,4].
In the classic BESO procedure, the number of elements to be rejected at the current iteration is prescribed by its ratio to the total number of elements at the previous iteration. This ratio is called as Evolutionary Volume Ratio (ER). The rejected elements are removed due to representing unnecessary regions so that their stiffness matrices are deleted. Based on the philosophy, the conventional BESO procedure can produce the following numerical issue due to improper ER. The smaller the value of the evolutionary volume ratio used, the more accurate is the final design, at the expense of larger computation time. The use of larger evolutionary volume ratio will reduce the number of elements of the resulting design more rapidly, whereas it may result in removing efficient elements incorrectly and fail to evolve to optimal solution. To address this issue, SESO (Smoothing Evolutionary Structural Optimization) procedure presents an organization of elements where a defined p% of rejected elements at each iteration are removed and (1 − p%) of them are returned to the structure [5,6]. For an unnecessary element for the structure, its corresponding structural stiffness gradually reduces till it no more contributes to the structure. This procedure can provide a typical characteristic of the continuous optimization. However, the key parameter of p% is selected according to a priori knowledge of optional solution without deterministic formula [7–9]. In this context, this paper applies an alternative topological optimization technique that provides a variable evolutionary volume ratio along the iterative process inside an extended fixed domain for a structure. It is suggested that one could use a high evolutionary volume ratio to sharply decrease the number of elements at early iterations while a low evolutionary volume ratio to obtain the sufficiently precious optimum at final iterations. This technique is an extension of the Morphing ESO methodology by Luo et al. [10]. The present variant of BESO presents advantages over the classical BESO method and these are demonstrated in this work. The numerical examples demonstrate that the improvements included in the formulation provide a compromise between sufficient accuracy and significant time saving.
Attenuation of unwanted vibrations is important in engineering structures as they could have detrimental effects on structural performances. Precious topology optimization focuses mainly on maximizing single dynamic performance like natural frequency [11,12], modal damping [13], or frequency response [14]. Topology optimization minimizing dynamic responses in time domain is another case of more difficulties. The equivalent static load (ESL) method is an increasingly popular approach to solve dynamic response structural optimization problems. By ESL method the dynamic topology optimization of a structure can be transformed into a static one under multiple loading cases [15,16]. Sun et al. [17] performed topology optimization of a three-dimensional flexible multi-body system via equivalent static load in the moving morphable components (MMC) based frame. Xu et al. [18] proposed a method for the concurrent topology optimization of macro-and micro-structural material distribution under dynamic loading based on ESL in the time domain. Up to the present, dynamic stiffness and frequency behavior of structures as two important factors in real-design problems have been deeply studied as separate objective functions. Therefore, achieving a trade-off between these two performances can have vital importance during structural design and analysis related to shipbuilding, automotive and aviation industries.
Over the last decades, a considerable effort has been devoted to single-objective optimization problems. But there may be several objective functions, usually conflicting in many real cases. As such, it is more appropriate to generate a set of optimal solutions which constitute the so-called Pareto set. For example, Simonetti et al. [19] explored the application of an evolutionary optimization technique for multi-objective optimization problems using the stress and strain energy criteria. Xu et al. [20] performed a mixed integer linear programming for multi-objective optimization of tensegrity structures using the ground structure method. Sleesongsom et al. [21] presented a multi-objective reliability-based topology optimization considering uncertain structural parameters using a fuzzy set model. Recently, Teimouri et al. [22] and Zhu et al. [23] have performed multi-objective optimization of continuum structures considering static stiffness and natural frequency. Crescenti et al. [24] generated minimal Pareto sets in multi-objective topology optimization of the wing box structural layout using smart normal constraint method. Lim et al. [25] obtained the optimal topology of a periodic spaceframe structure for maximal effective flexural and torsional stiffnesses along with minimal mass by Genetic Algorithm. Simonetti et al. [26] explored the application of the SESO technique to implement a parallel optimization to minimize the Von Mises stress and the internal strain energy growth.
In the conventional optimal design procedure, the natural frequency and the dynamic stiffness are exclusively considered as two independent factors for practical design problems. Consequently, it is of significance to establish a balance between these two indexes at the conceptual design stage. In this article, a multi-objective optimization scheme is implemented to obtain optimal topologies of a structure considering the natural frequency and the dynamic stiffness. The modified BESO with variable EVR is also introduced to improve the efficiency and stability during optimization.
The remainder of the article is organized as follows. Section 2 introduces the variable evolutionary volume ratio representing the dynamic removal behavior during optimization. Section 3 formulates the mono-objective optimization problem in the natural frequency and the dynamic stiffness. A weighted sum technique according to the sensitivity information obtained from mono-objective optimization is employed to carry out the multi-objective optimization. Section 4 presents several numerical examples to verify the proposed algorithm. Concluding remarks are made in Section 5.
2 Modified BESO with Variable Evolutionary Volume Ratio
Huang et al. [2] treated material removal and addition with a fixed evolutionary volume ratio to impose modifications on the topology of a structure, using the optimality criterion for the gradual removal of the finite elements in the mesh, which do not effectively contribute to a better performance of the structure. According to [2], the target volume for the next iteration (
where
The evolutionary volume ratio in the BESO method plays the similar role as the move limit in mathematical programming and the step size in optimality criteria methods. Therefore, it is an important parameter impacting on evolution process and optimal solution. It is expected that a high evolutionary volume ratio is used to significantly decrease the number of inefficient elements at early iterations while a low evolutionary volume ratio to obtain the sufficiently accurate optimum at final iterations. Consequently, a reasonable evolutionary strategy should gradually reduce the evolutionary volume ratio as the iteration proceeds. The evolutionary procedure proposed can be performed by
where
These functions are continuous and differentiated with an image varying from 0 to 1, as shown in Fig. 1. Nonetheless, here the sinusoidal function is used to regulate the evolution process. It is noteworthy that for larger value of
3 Multi-Objective Design Optimization
Topological design considering multiple dynamic characteristics is of great importance for a real engineering structure. A multi-objective topology optimization technique is implemented to obtain optimal topology for Maximizing dynamic stiffness and natural frequency. In this case, weighting factors are imposed on multi-objective sensitivity numbers to reflect the importance of both dynamic stiffness and natural frequency. As a result, using several combinations of weighting factors, the resulting different topology dependent on the level of importance is obtained representing a Pareto-optimal solution.
3.1 Sensitivity Number for Dynamic Stiffness
To attenuate the unwanted vibration, the averaging summation of the dynamic strain energy during the time domain is defined as the objective function for dynamic stiffness in this research. The general formulation can be stated as follows:
In the above formulation,
Using the vibration theory with the finite element method, the dynamic behavior of a continuum structure is expressed by the following differential equations:
where
Based on equivalent static loads method (ESLM) proposed by Jang et al. [15] for dynamic response topology optimization and verified by Stolpe [16], an ESL set,
where
Then the present topology optimization problem can be transformed into that for multiple load cases, which can be rewritten as
According to [2,4], the sensitivity number for static stiffness with single load case is generally defined by the following equation:
where
With the static multiple load case in Eq. (9) equivalent to the original one in Eq. (6), the sensitivity of the dynamic strain energy with the design variable can be written as
3.2 Sensitivity Number for Natural Frequency
For a solid-void design, the topological optimization problem of maximizing the l-th natural frequency
Using finite element analysis, the dynamic behavior of a structure may be expressed by the following eigenvalue equation:
where
The sensitivity number of the j-th element for natural frequency can be formulated as [8]
where
3.3 Multi-Objective Sensitivity Number
Since the obtained multi-objective sensitivity number for dynamic stiffness should be combined with the sensitivity number for natural frequency, it is normalized as follows:
where
In the same way, the sensitivity numbers for natural frequency can be normalized by the following equation:
where
The multi-objective sensitivity number for both dynamic stiffness and natural frequency is defined by the following Eq. (17), where the weighting factors of both dynamic stiffness and natural frequency are imposed on sensitivity numbers for them to reflect the importance of dynamic stiffness and natural frequency, respectively.
where
To circumvent the checkerboard pattern, a smoothing filter scheme [4] is implemented to blur the element sensitivities using a low-pass filter of radius
where
Then the above nodal sensitivity numbers will be converted to smoothed elemental sensitivity numbers by summing up weighted
where
where
Therefore, in terms of the above discussion, the flowchart of multi-objective topology optimization problem is shown in Fig. 2. To verify the modified BESO method, the multi-objective optimization problem is solved using the classic BESO method compared with the resulting optimal designs generated by the present methodology.
The design domain and loading condition of a long slender beam is illustrated in Fig. 3. The beam is 140 mm long and 20 mm high. All the degrees of freedom at the left and right sides are fixed. A dynamic force with amplitude 2.0 kN and duration 0.2 s, is applied at the center of the bottom edge. The structure is discretized into 2800 four nods plane stress elements. The Young’s modulus, the Poisson’s ratio, and the density are ρ = 7860 Kg/m3, E = 724 GPa, μ = 0.3, respectively. The final volume is restricted as 50% of the initial design.
According to the comprehensive sensitivity number in Eq. (17), the five weighted factors
Figs. 4–8 illustrate the evolution histories of dynamic compliance, natural frequency as well as the volume fraction for various weighted factors using the modified BESO method. In terms of Figs. 4–8, large weighted factor
Figs.10–14 show the evolution histories of dynamic compliance, natural frequency as well as the volume fraction for various weighted factors using the classic BESO method with constant EVR. Fig. 15 depicts the resulting optimal topology corresponding to various weighted factors. Table 1 compares the results obtained by the classic BESO with those obtained by the modified BESO. It is evident that the modified BESO method can reproduce the optimal solutions obtained from the classic BESO. It is verified that the dynamic removal strategy with variable EVR substantially saves the computational time when compared with that with constant EVR. It attributes to the fact that high EVR at early iterations is favorable to significantly remove inefficient elements while low EVR at final iterations to obtain the sufficiently accurate optimum.
The dynamic stiffness and natural frequency are parallelly optimized as a two-objective function in the structural topology by modified BESO method. A weighted sum method is introduced to establish a balance between these two objectives. The numerical results reveal that the optimal topology from the multi-objective optimization problem is provided with a large dynamic stiffness compared with that exclusively from maximization of natural frequency, and a high natural frequency compared with that exclusively from minimization of dynamic stiffness. Any improvement in one objective performance requires a certain amount of compensation with the other objective performance. This type of topology optimization scheme is suitable for structural design in vehicle and aerospace industries where dynamic stiffness and natural frequency are equally important.
The present examples demonstrate that the modified BESO method with variable EVR can produce similar Pareto-optimum sets to those generated by classic BESO method with constant EVR. However, the dynamic removal strategy with the variable EVR presents a low computational cost since it consumes a small number of iterations to capture the accurate optimal topology using the constant EVR. Although only two objective functions are considered in this study, the scheme can be extended to other multiple objectives such as stress, displacement and frequency or even to thermal, fluidic and acoustic applications. These will be investigated and reported in the near future.
Funding Statement: This study was funded by the National Natural Science Foundation of China (Grant No. 51505096), and the Natural Science Foundation of Heilongjiang Province (Grant No. LH2020E064).
Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.
References
View
Download
Like
Optimal Mass Distribution Prediction for Human Proximal Femur with Bi-modulus Property
Jiao Shi, Kun Cai, Qing H. QinA Multi-Criteria Topology Optimization for Systematic Design of Compliant Mechanisms
Zhen Luo, Nong Zhang,Interval-Based Uncertain Multi-Objective Optimization Design of Vehicle Crashworthiness
F.Y.Li, G.Y.LiMulti-Disciplinary Optimization for Multi-Objective Uncertainty Design of Thin Walled Beams
Fangyi Li, Guangyao Li, Guangyong...