Issue 
Int. J. Simul. Multisci. Des. Optim.
Volume 5, 2014



Article Number  A09  
Number of page(s)  9  
DOI  https://doi.org/10.1051/smdo/2013013  
Published online  06 February 2014 
Article
Structural optimization in ESAC: annals 2011
Engineering Simulation and Aerospace Computing, Northwestern Polytechnical University, 710072, Xi’an, Shaanxi, P.R. China
^{*} email: jh.zhu@nwpu.edu.cn
Received:
2
June
2012
Accepted:
23
October
2013
The purpose of this paper is to give an overall introduction of the structural optimization research works in ESAC group in 2011. Four main topics are involved, i.e. 1) topology optimization with multiphase materials, 2) integrated layout and topology optimization, 3) prediction of effective material properties and 4) composite design. More detailed techniques and some numerical results are also presented and discussed here.
Key words: ESAC annals / Topology optimization / Shape optimization / Composite
© J.H. Zhu, W.H. Zhang, Published by EDP Sciences, 2014
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The Laboratory of Engineering Simulation and Aerospace Computing (ESAC) in Northwestern Polytechnical University is an active international research team in the field of structural optimization design. Since 1990s, we have been dedicated to the theoretical study of shape and topology optimization and their applications in aeronautics and aerospace industries.
In 2011, the research works and the published papers are mainly concentrated on the research of the following four topics, i.e., (1) topology optimization with multiphase materials, (2) integrated layout and topology optimization, (3) prediction of effective material properties and (4) composite design.
2 Topology optimization with multiphase materials
This work is focused on the topology optimization of lightweight structures consisting of multiphase materials [1, 2]. Instead of adopting the common idea of using volume constraint, a new problem formulation with mass constraint is firstly proposed, which is a significant breakthrough of the topology optimization with multiphase materials. For a discretized structure with n designable finite elements, the mass constraint is stated as(1)where ρ_{i} denotes the material density related to element i, which has to be interpolated from the set of available material densities. M is the upper bound controlling the structural mass.
Meanwhile, to have a deep understanding of the mass constraint, recursive multiphase materials interpolation (RMMI) and uniform multiphase materials interpolation (UMMI) schemes are discussed and compared based on numerical tests and theoretical analysis.
RMMI interpolation scheme can be treated as an extension of the traditional SIMP/RAMP interpolation model reserved for the case of only one solid material phase and the void to the case of multiphase materials. This scheme is constructed through a recursive process. For a structure with m solid materials and n designable finite elements, the RMMI scheme interpolates Young’s modulus of element i in the following manner:(2)where(3)where E_{(j)} denotes Young’s modulus of material phase j. x _{i} denotes the set of design variables associated with element i. Clearly, each element has m design variables that have different significations. v_{ij}(x_{ij}) is the weighting function related to element i.
It is indicated that the nonlinearity of the mass constraint introduced by RMMI brings numerical difficulties to attain the global optimum of the optimization problem. On the contrary, the UMMI2 scheme makes it possible to formulate the mass constraint in a linear form with separable design variables.
UMMI originated from the socalled DMO. It means that the continuous laminate design of fiber orientations was treated as a DMO problem. In this paper, the concept of DMO is generalized for structural topology optimization with multiphase materials. As the weighting function is uniform for each material phase, the interpolation scheme is, in fact, a uniform multiphase materials interpolation (UMMI). For a problem of m solid materials and n designable elements, each element will have m design variables. The UMMI interpolation of Young’s modulus can be expressed as the weighted sum of all candidate material phases: (4) w_{ij} is the weighting function associated with the jth material phase which can be defined in UMMI2 as(5)
As shown in Figure 1, one such formulation favors very much the problem resolution by means of mathematical programming approaches, especially the convex programming methods. Moreover, numerical analysis indicates that fully uniform initial weighting is beneficial to seek the global optimum when UMMI2 scheme is used. Besides, the relationship between the volume constraint and mass constraint is theoretically revealed. The filtering technique is adapted to avoid the checkerboard pattern related to the problem with multiphase materials.
Figure 1. UMMI2 interpolation model with two solid materials (m = 2, p = 3). 
Numerical examples show that the UMMI2 scheme with fully uniform initial weighting is reliable and efficient to deal with the structural topology optimization with multiphase materials and mass constraint. Meanwhile, the mass constraint formulation is evidently more significant than the volume constraint formulation.
As shown in Figures 2a and 2b, with the RMMI scheme and uniform x_{ij} or w_{ij}, the optimal structures have configurations consisting of both VM1 and VM2. Comparatively, the result shown in Figure 2c with only one single solid material (VM2) is found to be stiffer. This implies that the RMMI scheme cannot yield the global optimum at least even uniform initial weighting are used.
Figure 2. Optimization results of a cantilever beam. 
Using the UMMI2 scheme with fully uniform initial variables, the optimal structure only consists of VM2, as shown in Figure 2d, even VM1 andVM2 are available. Although configurations given in Figures 2c and 2d are very similar, the structural compliances using UMMI2 are smaller. Besides, partial uniform initial variables of Case B are adopted and the results are shown in Figures 2e and 2f. Obviously, different initial variables lead to totally different configurations. Note that some elements at the rightbottom corner in Figure 2f consist of “mixed” materials.
Another interesting thing is to compare the optimization results obtained with mass and volume constraints. To do this, suppose the structure mass is less than 25.92 × 10^{3} kg. The upper bound of the volume constraint can be correspondingly determined for VM3 if the upper bound of the volume constraint of VM1 is fixed a priori. For different volume fractions of VM1, optimization results are plotted in Figure 3. The optimal solution with the mass constraint corresponds to the volume fractions of VM1 and VM3 being 0.463 and 0.068, respectively. It can be seen that this is the best solution among all possible ones using volume constraint in the sense of compliance minimization. In other words, the mass constraint always leads to a better solution than the volume constraint with the precondition of the same structure mass. Meanwhile, the percentage of each solid material phase will be automatically identified. In fact, the mass constraint has a more significant and important sense than the volume constraint in most engineering problems.
Figure 3. Comparisons of the optimization results under the mass and volume constraints with the equal structure mass. 
From the development of this work, it is seen that the combination of UMMI2 scheme together with fully uniform initial variables and the convex programming methods provides an efficient approach for the structural layout design of multiphase materials in static load case when the mass is controlled. Meanwhile, due to its generality, this approach may be extended to the structural design and microstructure design of multiphase materials in generalized load cases.
3 Integrated layout and topology optimization
In many engineering products, typical structural systems consists of components and structures with integrated interconnections inside a limited packing space. Global system performances significantly depend upon the layout of components and the topology of supporting structures. The design of multicomponent structural systems is a relatively new area and is still not extensively researched. Since most real life engineering design problems involve multiple components or structures, it is a subject of great relevance
Recently, further researches are carried out in our group with significant improvements and systematical methodologies [3, 4]. Packing optimization is performed by means of the FCM (Finite Circle Method), and the supporting structure is designed by means of topology optimization techniques. Because movable components inevitably change the FE mesh and locally corrupt the updating of the material distribution, the traditional topology optimization method using a fixed FE mesh is no longer applicable. To circumvent this problem, a density point technique was proposed to associate the density variables with the centroids of a fixed background mesh. As shown in Figure 4, density points are a set of predefined points in the design domain to control the material distribution. No matter how the mesh changes, each element receives the pseudodensity value from its nearest density point.
Figure 4. Element material properties received from the density points. 
In order to reduce the computing cost and simplify sensitivity analysis with respect to location design variables of each component, the superelement technique is adopted to model the movable components.
Similar to classic topology optimization, the sensitivity of C with respect to pseudodensity design variables can be analytically obtained by(6)
As to the sensitivity with respect to location design variables of components, a semianalytical scheme combined with superelement formulation is applied in the following way. The global stiffness matrix of the system is partitioned into three parts: supporting structure, mesh transition region, and components. The derivative of the system compliance with respect to s_{ε} representing one of the three location parameters x_{ε}, y_{ε}, and θ_{ε} of component ε gives rise to(7)
A design example of simple structure system with two rectangular components is considered here. The global strain energy is the objective function to be minimized and the volume fraction is set to be 50%. Assume that all components are made of steel, and the supporting structure is made of aluminum. As shown in Figure 5, the components and the supporting structure are optimized simultaneously. Finally, the components are located in the main loading path of the structure due to their high stiffness.
Figure 5. Integrated layout design with components of stiffer materials: (a) iteration 5, (b) iteration 20, and (c) final design; final strain energy C50.074 J. 
Design of multicomponent systems can also be extended to support layout of bridgelike problems. Fixations are, here, defined as movable solid support components, which are partially embedded into the design domain. Positions of these support components are optimized together with the structural topology to obtain the maximum stiffness of the global system.
A typical example is illustrated in Figure 6. Materials used for the structure and support components are aluminum and steel, respectively. The road on the bridge is a nondesignable area. Four identical support components are located symmetrically along the bottom boundary of the design domain. A uniform surface pressure is applied on the top of the bridge. The volume fraction for the design domain is set to be 40%. By optimization, it is found that a slight change of the support locations will lead to a significant variation of the structural topology.
Figure 6. Integrated support and layout design of a bridgelike structure: (a) iteration 5, (b) iteration 10, and (c) final design, C50.1298 J. 
It is possible to extend the methodology to the layout design of 3D multicomponent systems. However, the embedded meshing technique for 3D problems is much more complicated than in 2D situation. Meanwhile, the semianalytical sensitivity analysis scheme for location design variables of components is timeconsuming because each 3D component has now six location design variables. Recently, a Heaviside functionbased approach is investigated for sensitivity analysis with respect to location design variables. The material discontinuity across the component boundary is approximately smoothed by means of the Heaviside function. When a location design variable is perturbed, sensitive elements attached to the component boundary are assumed to undertake a material property change, while the mesh remains unchanged. In this way, sensitivity analysis with respect to location design variables is carried out as easily as for pseudodensity variables. Figure 7 shows the optimization result of a 3D structure, in which tetrahedral elements are used and the involved components have a greater Young’s modulus than that of the supporting structure.
Figure 7. Integrated layout design of a 3D structure with three components: (a) iteration 2, (b) iteration 5, (c) iteration 10, (d) iteration 18, and (e) iteration history. 
4 Prediction of effective material properties
This work presents a microstructural modeling method considering the oxidation effect. Further studies on the prediction of material properties were carried for the multiphase and multilayer braided composite [5]. Regarding the multiscale characteristics of the composite, the microstructure modeling is carried out sequentially from fiber to tow scale. The geometrical configuration of the microstructure is first analyzed, and mathematical relations among different geometrical parameters are derived on each scale. Second, effective elastic properties are obtained based on the sequential homogenization from fiber to tow scale. A strain energy based method is proposed to evaluate effective elastic properties with specific boundary conditions imposed on the microstructure. Numerical results obtained by the proposed method and the microstructure model show a good agreement with the results measured experimentally.
The present study aimed at the microstructure modeling and numerical prediction of elastic properties of 3D fourstep MPML braided composites. First, the microstructure modeling is carried out sequentially from fiber to tow scale. Geometrical descriptions and finite element modeling of the microstructure are presented in detail.
The RVC (Representative Volume Cell) model of the 3D fourstep MPML braided composite involves two scales: fiber and tow scales. The first scale concerns the modeling of RVC for MPML tows, and the second scale concerns the modeling of RVC for 3D fourstep MPML braided composites as shown in Figure 8. Considering the multiscale characteristics, the multiscale analysis methods are needed for the analysis of MPML plain weave composites. The multiscale analysis procedure is as follows. First, on the fiber scale, finite element model is built to obtain the effective moduli of the fiber scale RVC, and the results are used for the tows that were treated as a homogeneous transversely isotropic material. Second, based on the obtained elastic properties of tows, finite element model of the tow scale RVC is created to evaluate effective elastic properties of the composite. Here, RVC models of both scales are established using the ANSYS finite element software.
Figure 8. Images (SEM) of 3D braided MPML C/CSiC. 
On the tow scale, the finite element model of RVC is depicted in Figure 9. Here, a 10node tetrahedral solid element SOLID187 is used for the discretization, and the finite element models of both tows and matrices are shown respectively.
Figure 9. Finite element model of RVC on tow scale. (a) Tow, (b) First layer of matrix, (c) Second layer of matrix, and (d) Last layer of matrix. 
Second, a strain energy based method is applied to evaluate the effective elastic properties. Numerical results are then compared with the experimental results to validate the modeling.
In the elastic regime, the macroscopic behaviors of the RVC can be characterized by the effective stress tensor σ and strain tensor ε over the homogeneous equivalent model. They are interrelated by the effective, also termed homogenized, stiffness matrix C ^{ H } (8)
The matrix coefficient C ^{ H } _{1111} can be derived(9)
In the same way, demonstrations can be made for other coefficients. The effective stiffness matrix can be written as(10)
The elastic constants can be derived by inversing the above matrix. In practice, the considered RVC will be discretized as a finite element model on which the initial strain will be imposed to evaluate the strain energy. The comparison with the experimental results shows the validity and the rationality of the present method. Future studies will be focused upon the numerical prediction of thermal properties (the thermal expansion coefficient and thermal conductivity), progressive failure analysis and design optimization of the microstructure for 3D MPML braided composite.
Similar model incorporates the modeling of oxidized microstructure and computing of degraded elastic moduli can also be extended to simulate the oxidation behaviors of 2D C/SiC composites exposed to air oxidizing environments below 800 °C [6]. Regarding the multiscale characteristics of 2D C/SiC composite, the microstructure modeling is carried out on microscopic and macroscopic scale, respectively to compute the degraded elastic properties in terms of time duration, temperature and pressure, whose influences upon the oxidation microstructure morphology and degraded properties of 2D C/SiC composites are also investigated. The predicted microstructure morphology after oxidation is found in good agreement with experiment data as shown in Figure 10.
Figure 10. Morphologies of microscopic RVC at different oxidation time. (a) 10 h, (b) 20 h, (c) 30 h, and (d) experimental observation. 
It is found that the computed inplane extensional moduli after oxidation are in good agreement with experimental data over the oxidation time. In addition, numerical results show that the temperature and pressure have a profound influence upon the oxidation results. The ablation of carbon phases is aggravated with the increase of temperature and pressure and therefore the elastic modulus of composite decreases gradually. An immediate perspective is to address other composite architectures such as the 3D braided C/SiC composite and to extend the developed oxidation model to other oxidizing environments, like water or gaswater environments. Besides, it is worth to note that the oxidation simulation of 2D C/SiC composite above 800 °C is not taken into account in the current work because the oxidation scheme is nonuniform and much complicated.
5 Composite design
Composite design with fiber orientation optimization plays an important role in aircraft structure design. Earlier orientation optimization methods face the main problem of huge number of design variables. Recently, a patch concept was proposed to reduce the number of design variables as shown in Figure 11. The elements in one patch share identical definition of fiber orientation. However, the traditional stressbased method cannot deal with patch orientation optimization of composite structures. In this paper, we propose an extended stressbased method to deal with such problems. The considered problems are to minimize the mean compliance under multiple load cases or to maximize the eigenvalues of a composite structure. Four numerical examples are solved in order to demonstrate the efficiency of the new method. It is shown that the new method has the ability to deal with constraints on orientation angle, such as symmetric, antisymmetric and discrete orientation angle constraints. The iteration is less timeconsuming because no sensitivity analysis is needed and a quick convergence rate can be achieved.
Figure 11. Definitions of layers, patches, and finite elements. 
In this work [7], the stressbased method is extended to deal with optimal patch orientation problems for composite structures under single and multiple load conditions. Some special constraints related to patch orientation angles are applied for the symmetric and antisymmetric laminates in which patch orientation angles have to take the same and the opposite values, respectively. Besides, according to manufacturing considerations, investigations are made to the cases where the orientation angles are restricted to discrete values of 0°, 45°, −45°, and 90°. Several examples are finally given to show the validity of this method.
For thin laminated composite structures, the planestress condition may be applied due to the fact that outofplane stresses are often small compared with inplanestresses. Suppose each layer of the laminate is made of orthotropic materials. Figure 12 shows a patch over a layer where the stressstrain relationship of a point is(11)where ε = {ε_{1}, ε_{2}, γ_{12}}^{ T } and σ = {σ_{1}, σ_{2}, τ_{12}}^{T} represent the planar strains and stresses in the primed coordinate system (O123), respectively; S _{0} is the unrotated orthotropic compliance matrix.
Figure 12. Orthotropic material. 
Unlike the traditional stressbased method related to each element, the Extended stressbased method (ESBM) expands the integration over the patch domain so that(12)where θ_{p} is the orientation variable related to patch p; Ω_{p} represents the spatial domain of patch p.
Let σ_{x}, σ_{y}, and τ_{xy} denote components of the stress tensor in the unprimed coordinate system (Oxyz), equation (12) can explicitly be developed as a function of the patch orientation variable θ_{p} (13)where coefficients A, B, D and E are computed as(14)
Note that the numerical Gaussian integration can be applied for computing coefficients in equation (13). In practice, if the finite element mesh is dense enough, the stress in each element can be assumed to be a constant approximated by the value at the centroid of the element.
When θ_{p} ≠ π/2, roots of equation (13) can equivalently be obtained by solving the following fourthorder polynomial equation in terms of tan θ_{p} within the interval (−π/2, π/2).(15)
By means of equation (13), the following evaluation function can be defined to determine the optimal value of θ_{p} selected from the real roots of equation (15) and θ_{p} = π/2.
Note that if each patch has just one finite element, the ESBM is identical to the traditional stressbased method.
One example as shown in Figure 13 is to show the validity of the ESBM in dealing with constraints of patch orientation variables. The number of patches is 64 and each patch consists of 8 × 8 elements. The left side of the cylinder (Z = 0) is clamped and a bending moment Mx is applied on the right side (Z = 36) of the structure. Assume the laminated composite with four layers is used and they are numbered sequentially from the inside to the outside of the cylinder. Likewise, the objective function is to minimize the structural mean compliance.
Figure 13. Optimal patch orientation angles with antisymmetric laminate. 
Acknowledgments
This work is supported by National Natural Science Foundation for Distinguished Young Scholar (10925212), National Natural Science Foundation of China (90916027, 11002113, 51275424, 11172236), 973 Program (2011CB610304), the 111 Project (B07050), NPU Foundation for Fundamental Research (JC20120229).
References
 Gao T, Zhang WH. 2011. A mass constraint formulation for structural topology optimization with multiphase materials. International Journal for Numerical Methods in Engineering, 88(8), 774–796. [CrossRef] [Google Scholar]
 Gao T, Zhang WH, Duysinx P. 2011. Topology optimization of structures designed with multiphase materials: volume constraint or mass constraint? Chinese Journal of Theoretical and Applied Mechanics, 43(2), 296–305. [Google Scholar]
 Zhang WH, Xia L, Zhu JH, Zhang Q. 2011. Some recent advances in the integrated layout design of multicomponent systems. Journal of Mechanical Design, 133(10), 104503115. [Google Scholar]
 Xia L, Zhu JH, Zhang WH. 2011. New strategies for the efficient integrated layout design of multicomponent system. Journal of Mechanical Engineering, 46(10), 148–155. [Google Scholar]
 Xu YJ, Zhang WH, Domaszewski M. 2011. Microstructure modelling and prediction of effective elastic properties of 3D multiphase and multilayer braided composite. Material Science and TechnologyLondon, 27(7), 1213–1221. [CrossRef] [Google Scholar]
 Xu YJ, Zhang WH. 2011. Numerical modelling of oxidized microstructure and degraded properties of 2D C/SiC composites in air oxidizing environments below 800 degrees C. Material Science and Engineering A, 528(27), 7974–7982. [CrossRef] [Google Scholar]
 Zhang J, Zhang WH, Zhu JH. 2011. An extended stressbased method for orientation angle optimization of laminated composite structures. Acta Mechanica Sinica, 27(6), 977–985. [CrossRef] [MathSciNet] [Google Scholar]
Cite this article as: Zhu JH & Zhang WH: Structural optimization in ESAC: annals 2011. Int. J. Simul. Multisci. Des. Optim., 2014, 5, A09.
All Figures
Figure 1. UMMI2 interpolation model with two solid materials (m = 2, p = 3). 

In the text 
Figure 2. Optimization results of a cantilever beam. 

In the text 
Figure 3. Comparisons of the optimization results under the mass and volume constraints with the equal structure mass. 

In the text 
Figure 4. Element material properties received from the density points. 

In the text 
Figure 5. Integrated layout design with components of stiffer materials: (a) iteration 5, (b) iteration 20, and (c) final design; final strain energy C50.074 J. 

In the text 
Figure 6. Integrated support and layout design of a bridgelike structure: (a) iteration 5, (b) iteration 10, and (c) final design, C50.1298 J. 

In the text 
Figure 7. Integrated layout design of a 3D structure with three components: (a) iteration 2, (b) iteration 5, (c) iteration 10, (d) iteration 18, and (e) iteration history. 

In the text 
Figure 8. Images (SEM) of 3D braided MPML C/CSiC. 

In the text 
Figure 9. Finite element model of RVC on tow scale. (a) Tow, (b) First layer of matrix, (c) Second layer of matrix, and (d) Last layer of matrix. 

In the text 
Figure 10. Morphologies of microscopic RVC at different oxidation time. (a) 10 h, (b) 20 h, (c) 30 h, and (d) experimental observation. 

In the text 
Figure 11. Definitions of layers, patches, and finite elements. 

In the text 
Figure 12. Orthotropic material. 

In the text 
Figure 13. Optimal patch orientation angles with antisymmetric laminate. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.