Fundamentally, design is the practice of defining a distribution of some resource(s) in a space to accomplish an objective. Many approaches have been developed to help enhance and streamline structural design over countless decades, but none do so with the directness of topology optimization (TO).
This approach to structural design was pioneered in the 1980’s by Bendsoe and Kikuchi [1] and considers the design problem in terms of:
This set of definitions is paired with an optimization algorithm to iteratively optimize the material distribution within the design domain. The resulting optimization problem can be expressed as
|
minχ f(ρ(χ))s.t.g(ρ(χ))≤0h(ρ(χ))=0 |
( 1 ) |
A common variation of this problem to minimize compliance by finding an optimal distribution of finite element densities, ρe, with a volume fraction constraint
|
minρe f(ρe)=F⊤d(ρe)s.t.vmax−∑eρeve≤0F=K(ρe)d0≤ρe≤1∀e |
( 2 ) |
In this problem formulation, an element stiffness of ρe=1 indicates that the element has the same stiffness as the material that it represents, while an element with ρe=0 has no stiffness and represent a region that is void. However, this problem formulation does not preclude the formation of regions with intermediate densities of 0<ρe<1, which calls for a definition of a function for the constitutive properties as a function of ρe. For example, Figure 1 shows a topology optimization result for a case where
|
E(ρe)=ρeE0+Emin |
For compliance minimization, the sensitivity of the objective with respect to element density is
|
∂f∂ρe=−d⊤e(∂∂ρeKe)de |
The element stiffness matrix, Ke, can be expressed as
|
Ke=(ρeE0+Emin)K0e |
where K0e is the stiffness matrix for a totally solid element. Thus, the objective sensitivity becomes
|
∂f∂ρe=−E0d⊤eK0ede |
( 6 ) |
The sensitivity of the constraint function is
|
∂g∂ρe=−ve |
( 7 ) |
where ve is finite element volume.
If intermediate density is allowed within the solution, then it is also necessary to define a physical interpretation of intermediate material density. One such interpretation is that intermediate density regions are regions populated by architected materials with volume fractions equal to ρe, and this interpretation is discussed in section 1.1. Alternatively, intermediate density can be assumed to be undesirable, and measures can be added to the problem formulation to avoid it in topology optimization solutions, as discussed in section 1.2.
Figure 2 – Boundary conditions for the MBB beam problem with the design domain marked in gray |
An Architected Material Interpretation of Topology Optimization
Architected Materials Background
An “architected material” is a material that consists of small cellular “microstructures”. Due to the geometrical properties of these microstructures, such architected materials exhibit properties that are different from the properties of the material(s) from which they are made. Two examples of naturally occurring examples of architected materials that are bones and wood (Figure 3) [2]. Research into manmade architected microstructure materials began to gain traction in the latter half of the 20th century [3]. These manmade architected materials can broadly be separated into two categories: random and structured (Figure 4), the latter of which can often be described in terms of “unit cells”. The impetus for exploring such materials is that they provide design engineers with access to material property combinations that are not attainable with naturally existing or conventional materials (Figure 5).
References
[1] Bendsøe, M. P., and Kikuchi, N., 1988, “Generating Optimal Topologies in Structural Design Using a Homogenization Method,” Comput. Methods Appl. Mech. Eng., 71(2), pp. 197–224.
[2] Schaedler, T. A., and Carter, W. B., 2016, “Architected Cellular Materials,” Annu. Rev. Mater. Res., 46(1), pp. 187–210.
[3] Fleck, N. A., Deshpande, V. S., and Ashby, M. F., 2010, “Micro-Architectured Materials: Past, Present and Future,” Proc. R. Soc. Math. Phys. Eng. Sci., 466(2121), pp. 2495–2516.