طراحی چیدمان جامع از سیستم های چند جزئی با استفاده از XFEM و تجزیه و تحلیل حساسیت تحلیلی

This study presents the integrated layout optimization of multi-component systems using a fixed mesh. The optimization formulation is established under the framework of the extended finite element method (XFEM). The level set method is used to represent components and is combined with the XFEM to describe material discontinuities across elements. Sensitivity analysis is proposed with respect to geometric variables of components and pseudo-densities of the basic structure. An analytical shape sensitivity analysis method with respect to positions and shapes of components is developed. Both solid and void components are considered to show the efficiency and accuracy of the proposed shape sensitivity analysis method. Furthermore, a revised finite circle method that adapts shape changes of elliptical components is proposed for the definition of non-overlapping constraints. Finally, numerical examples of maximizing the structural stiffness are tested to demonstrate the proposed method.

Structural optimization, including topology, shape and size optimization, has been widely used in industrial and civil engineering to improve the performance of mechanical designs. As illustrated in Fig. 1, the integrated layout optimization studied here aims to find a proper position and shape design of components and the topology of the base structure in a limited base domain simultaneously. Embedded components may be considered as solids or substructures to afford loads or as voids to reduce the overall structure weight. Three kinds of interfaces thus exist between the base structure and the components, i.e., material–material (m–m) interface, material–void (m–v) interface and void–void (v–v) interface. Full-size image (20 K) Fig. 1. Illustration of the integrated layout optimization. Figure options In this paper, as the void part of the base structure is considered to be a special domain filled with materials of very weak stiffness as in traditional topology optimization [1], attention is paid only to the first kind of interface (m–m) that can be regarded as a generalization of the rest two ones. Historically, components were designed as a packing optimization problem in the previous work [2], [3], [4], [5] and [6] where only their positions were changeable without shape variation. The current work is an extended approach that combines the packing optimization and shape optimization together for more design flexibilities of integrated layout optimization. Clearly, there exist two main concerns in the integrated layout design of multi-component systems as discussed by Zhang et al. [6]. First, the avoidances of overlap between components and overlap between the component and the domain boundary. Second, shape sensitivity analysis with respect to geometric design variables related to positions and shapes of components. The overlapping avoidance is a very important issue for optimization of multi-component systems. In the earlier work of Qian and Ananthasuresh [2], each component was approximated by one circle so that the overlapping judgment between components turned out to determine the distance between two circles. However, when the component has a complex shape, especially concave, this approximation is very rough. To improve the approximation accuracy, Zhang and Zhang [7] proposed the so-called finite circle method (FCM) by which each component is approximated with a set of circles and non-overlapping constraints between components are thus converted into simple constraints between two sets of circles. Another underlying difficulty concerns shape sensitivity analysis related to the shape representation of the component. Explicit parametric representations, e.g., Bézier or NURBS [8], implicit representations, e.g., the signed distance function [9] or R-function [10] and combined representations [11] can be used. To the authors’ knowledge, there exist mainly three sensitivity analysis methods related to different shape representations and different mesh discretizations. Fig. 2 illustrates the case of a movable void component inside a meshed plate. Full-size image (58 K) Fig. 2. Three methods for shape sensitivity analysis of a void component. (a) Grid perturbation method with Lagrangian mesh. (b) Material projection method with Eulerian mesh. (c) XFEM based method with Eulerian mesh. Figure options (a) Grid perturbation method using Lagrangian mesh. As shown in Fig. 2a, the FE mesh will move in accordance with the variation of the component boundary that is often parameterized to define geometric design variables. This implies that the model is based on a Lagrangian mesh that varies when geometric design variables are perturbed or modified in shape sensitivity analysis and optimization process, respectively. As a result, the computational domain has to be remeshed for any change of geometric design variables [12]. This kind of method always works in conjunction with a mesh morphing technique [13] to relocate the nodes in sensitivity analysis while the mesh topology is kept unchanged. (b) Material projection method using Eulerian mesh. In order to avoid the sophisticated remeshing procedure, a fixed Eulerian mesh is used for the base domain independently of the shape variation of the component boundary. In this case, an implicit shape representation approach, e.g., level set method [9] is often used for the geometric description of the component boundary. To characterize the effect of boundary variation, the fixed mesh will, however, undergo a change of attributed material properties depending upon the location of each element. As illustrated in Fig. 2b, material attributes stay unchanged for inner elements that do not share the component boundary, while intermediate values of material attributes have to be assigned to boundary elements. A simple way is to take the area fraction of the solid part in the boundary elements for the correction of material properties, as discussed in the work of Allaire et al. [14], Kim and Chang [15] and Luo et al. [16] and [17]. However, the stress sensitivity is still questionable. Another way to make one such correction is the so-called filtering technique as used in the work of Norato et al. [18]. In the work of Belytschko et al. [19] and Wang et al. [20], an approximate Heaviside function was used for the correction. However, this kind of method might cause low accuracy for analysis results near the component boundary because the analysis model is not the same as the exact boundary shape representation [21]. Meanwhile, a rather fine mesh is sometimes needed. (c) The XFEM based method using Eulerian mesh. This method is based on a fixed mesh that works with the level set method firstly proposed by Osher and Sethian [9] to represent moving interfaces and has been extensively applied to structural optimization [14], [20], [16], [17] and [22]. The level set-based structural optimization method is essentially a type of shape optimization method but able to achieve topological changes during shape variations, which can flexibly describe complex geometries via topological shape evolvements by remaining a concise and smooth boundary. Until now, several level set methods have been proposed for the layout description of multi-phase materials. In the standard level set method [9], components of the same material are represented using one level set function. In the work of Wang and Wang [23], a color level set method which represents n different material phases uses log2nlog2n level set functions by allowing overlap of the level set functions. Luo et al. [24] proposed a piecewise level set method using one level set function of constant integer value for each material phase. The representation of multi-phase materials is fairly easy by making use of discontinuities of the piecewise level set function. Notice that the piecewise level set function is non-differentiable at material interfaces. As the overlap of components is avoided by FCM and the differentiability of the level set function is needed in sensitivity analysis, the standard level set method is employed in this work and combined with the XFEM [25], [26] and [27] for computing structural responses and sensitivities. As shown in Fig. 2c, whether an embedded void or solid component is concerned, the whole base domain is completely meshed. Material discontinuities within each element are handled by the XFEM based on the partition of unity method (PUM) [28]. As a result, the shape representation of the component boundary and the analysis model are unified to increase the accuracy of analysis results. As the finite element approximation is locally enriched by the XFEM, special treatments in shape sensitivity analysis should be considered correspondingly. Miegroet and Duysinx [29] adopted the semi-analytical method (SAM). The dimension change of the stiffness matrix of an enriched element is ignored when the design boundary is very close to a node. To some extent, this treatment might degrade the accuracy of sensitivity results. In the work of Wei et al. [21] and Shu et al. [30], shape sensitivity is calculated by using the shape derivative method. When the element stiffness matrix is calculated by the XFEM, a simplified strategy is used by removing the integration points on the void part of an enrichment element without additional degrees of freedom. As indicated by Daux et al. [31] and Miegroet and Duysinx [29], the inconvenience lies in that this simplification is only suitable to material–void interface problems but unsuitable to the considered material–material interface problems studied here. In this paper, the integrated layout optimization of multi-component systems is addressed within the XFEM framework. An analytical sensitivity analysis approach oriented to material–material interface problems is developed as a general scheme. It can also be used to material–void interface problems by setting the void phase to be a relatively weak material. Besides, a revised version of the FCM is proposed to follow the shape changes of elliptical components adaptively and avoid the overlap. This paper is organized as follows: Section 2 describes the optimization problem formulation in conjunction with some special treatments of the XFEM. Section 3 gives the discrete form of the problem as well as the XFEM computing procedure. The procedure of sensitivity analysis is presented in Section 4. Numerical optimization results are given in Section 5 to demonstrate the proposed method. Finally, a conclusion is drawn out in Section 6.

In this paper, the integrated layout optimization of multi-component system is investigated. An analytical shape sensitivity analysis approach is developed within the XFEM framework. Compared with the existing work, the main difference is that both positions and shapes of components are optimized. Although the sensitivity analysis approach is derived from the material–material interface model, it can also be used in material–void interface problems by setting a relatively weak material to the void without loss of the accuracy. Unlike the existing semi-analytical method, our approach makes it possible to improve the computing efficiency without remeshing and to eliminate the confusion of dimension change of the element stiffness matrix. Numerical tests including solid and void components are studied. It is shown that the revised FCM is able to avoid the overlap and to follow adaptively the shape changes of elliptical components. Stiffness layout designs are achieved with a stable convergence for a variety of geometries of components and base structures. Further investigations can be made to 3D cases, thermo-elastic and dynamic optimization problems of multi-component systems.