# Simultaneous Material and Topology Optimization Based on Topological Derivatives

Jannis Greifenstein, Michael Stingl

► To cite this version:

Jannis Greifenstein, Michael Stingl. Simultaneous Material and Topology Optimization Based on Topological Derivatives. 26th Conference on System Modeling and Optimization (CSMO), Sep 2013, Klagenfurt, Austria. pp.118-127, 10.1007/978-3-662-45504-3_11. hal-01286404

HAL Id: hal-01286404

https://inria.hal.science/hal-01286404

Submitted on 10 Mar 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Simultaneous material and topology
optimization based on topological derivatives*

Jannis Greifenstein and Michael Stingl

Chair of Applied Mathematics 2,
Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU),
Nägelsbachstr. 49b, 91052 Erlangen, Germany.

{greifenstein,stingl}@math.fau.de

http://www.math.fau.de

**Abstract.** We use an asymptotic expansion of the compliance cost func-
tional in linear elasticity to find the optimal material inside elliptic in-
clusions. We extend the proposed method to material optimization on
the whole domain and compare the global quality of the solutions for
different inclusion sizes. Specifically, we use an adjusted free material
optimization problem, that can be solved globally, as a global lower ma-
terial optimization bound. Finally, the asymptotic expansion is used as
a topological derivative in a simultaneous material and topology opti-
mization problem.

**Keywords:** material optimization, topology optimization, material ori-
entation, asymptotic expansions, discrete material optimization

# 1 Introduction

We investigate material and topology optimization of compliance problems in two dimensions. To this end, we first present an asymptotic formula<fnref n="1" /> of the compliance functional for the insertion of a number of ellipsoidal bodies in an elastic domain, originally derived in [1]. Later, we study numerically the feasibility of replacing all material using the same asymptotic expansion as for the ellipses and finally make use of the formula as topological derivative.

The problem described above is by far not new. There are many publications
dealing with similar types of problems. For the rotational optimization consid-
ered later, in [2] an analytical formula for the strain energy is derived, to directly
compute the optimal material orientation. In [3], this has been embedded into
a structural optimization algorithm for compliance minimization. A similar ap-
proach is discussed in [4] for a plate model. The method proposed in this article,
however, can be used for a broader spectrum in material optimization as well,

<footnote>* The authors want to thank the German Research Foundation (DFG) for funding this research work within Collaborative Research Centre 814, subproject C2.</footnote>

<footnote n="1">¹ We note that the asymptotic formulae are also available for the three-dimensional case.</footnote>
such as discrete material optimization. The algorithm for simultaneous material and topology optimization presented at the end of this article is very similar to topology gradient methods, see e.g. [5]. For more references, see [6].

## 2 Material and topology optimization

We consider a domain $\Omega \subset \mathbb{R}^2$ of isotropic elastic material. The domain is subject to exterior traction and other boundary conditions (e.g. homogeneous Dirichlet conditions). The objective is to find the optimal material $C^0$ in a set of admissible materials $C$ to insert into an inclusion, for which the compliance as defined in (1) is minimized.

The elastic body is modeled by the equation of linear elasticity

$$ \int_{\Omega} C_{ijkl}(x) \varepsilon_{ij}(u) \varepsilon_{kl}(v) \, dx = \int_{\Gamma} f u \, ds, $$

with the displacement field $u$, the linearized strains $\varepsilon_{ij}(u) = (\frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i})$ and a traction $f$. We use Voigt notation to denote the fourth-order stiffness tensor $C_{ijkl}$ by a symmetric 3 × 3-matrix and the strains and stresses by a vector

$$ C = \begin{pmatrix} C_{1111} & C_{1122} & \sqrt{2}C_{1112} \\ C_{2222} & \sqrt{2}C_{2212} & 2C_{1212} \\ \text{sym.} & 2C_{1212} & 2C_{2122} \end{pmatrix}, \quad \varepsilon = \begin{pmatrix} \varepsilon_{11} \\ \varepsilon_{22} \\ \sqrt{2}\varepsilon_{12} \end{pmatrix}, \quad \sigma = \begin{pmatrix} \sigma_{11} \\ \sigma_{22} \\ \sqrt{2}\sigma_{12} \end{pmatrix}. $$

The stresses are given by Hooke's law

$$ \sigma = C\varepsilon. $$

### 2.1 Optimal material in elliptic inclusions

We insert a finite number of inclusions $\omega_i$, $i = 1, \dots, n_{\text{ell}}$ with $n_{\text{ell}} > 1$ and centers $z_1, \dots, z_{n_{\text{ell}}}$ into the domain $\Omega$ and search for the optimal material to be used in the inclusions. For an exemplary setup of boundary conditions and loads, a sketch of a possible problem specification is shown in Fig. 1. We place the elliptic inclusions on a regular grid within the domain $\Omega$, so that the center points of the inclusions are distributed equidistantly. In the sketch in Fig. 1, we have $n_{\text{ell}} = 100$ disjoint elliptic inclusions. In the remaining domain

$$ \Omega^1 := \Omega \setminus \bigcup_{i=1}^{n_{\text{ell}}} \omega_i, $$

we insert an isotropic matrix material $C^1$, and, in each inclusion, a material $C_i$, $i = 1, \dots, n_{\text{ell}}$.
<figure><img src="image_2.png" /><figcaption>Fig. 1. Example for placing elliptic inclusions $\omega_i, i = 1, \dots, n_{\text{ell}}$, in the domain $\Omega$.</figcaption></figure>

Now, the compliance function of the elastic body $\Omega$ (without inclusions) with respect to a set of given load functions $f_k \in [L^2(\Gamma)]^2$, $k \in \mathcal{K} = \{1, 2, \dots, n_{\text{loads}}\}$ applied to a part $\Gamma$ of $\partial\Omega$ is defined as

$$ \mathcal{J}_c(f) = \sum_{k \in \mathcal{K}} \int_{\Gamma} f_k(x) u_k(x) \, ds. \qquad (1) $$

By virtue of the asymptotic expansion (cf. [1]) for the two-dimensional case, the compliance function $\mathcal{J}_{c(C^0)}(f)$ of a body with a single inclusion of material $C^0$ can be approximated as

$$ \left| \mathcal{J}_{c(C^0)}(f) - \left( \mathcal{J}_c(f) - \frac{h^2}{2} \varepsilon(u; z)^\top P_{(C^0)} \varepsilon(u; z) \right) \right| \le c h^{5/2} (1 + |\ln h|) \|f\|^2, \quad (2) $$

where $h$ is a dimensionless scaling parameter for the elliptic inclusion $\omega$, $\varepsilon(u; z)$ is the strain corresponding to the displacement of the domain $\Omega$ without inclusion evaluated at the center $z$ of the ellipse, and $P_{(C^0)}$ is the so called polarization matrix given by

$$ P_{(C^0)} = -|\omega| \left( C^1 - C^0 + (C^1 - C^0) (\mathbb{I}_3 - \Psi(C^1 - C^0))^{-1} \Psi(C^1 - C^0) \right). \quad (3) $$

The matrix $\Psi$ in formula (3) depends solely on the isotropic matrix material and is, for the stretched coordinates $\omega = \{\xi : \xi_1^2/a^2 + \xi_2^2/b^2 < 1\}$, $\xi = x/h$, computed by

$$ \Psi = \int_{\partial \omega} D(\nu(\xi)) (D(\nabla_{\xi})\Phi(\xi))^{\top} \, ds_{\xi}, \quad D(x)^{\top} = \begin{pmatrix} x_1 & 0 & 2^{-1/2}x_2 \\ 0 & x_2 & 2^{-1/2}x_1 \end{pmatrix}, \quad (4) $$
where $\nu$ denotes the exterior normal of the ellipse $\omega$ and $a, b$ are the semi-major and semi-minor axes of the ellipse, respectively, cf. Fig. 1. Finally, $\Phi$ in (4) is the fundamental solution given by

$$ \Phi(x) = c_1 \left( \begin{matrix} -(\lambda + 3\mu) \ln(r) - (\lambda + \mu) \frac{x_2^2}{r^2} & (\lambda + \mu) \frac{x_1 x_2}{r^2} \\ (\lambda + \mu) \frac{x_1 x_2}{r^2} & -(\lambda + 3\mu) \ln(r) - (\lambda + \mu) \frac{x_1^2}{r^2} \end{matrix} \right), \quad (5) $$

where $c_1 = \frac{1}{4\pi} \frac{1}{\mu(\lambda+2\mu)}$, $r = \sqrt{x_1^2 + x_2^2}$ and $\lambda$ and $\mu$ are the Lamé parameters corresponding to the isotropic material $C^1$, see e.g. [7, Chap. 3]. A more detailed explanation may be found in [1, Sec. 3.1].

In order to minimize the compliance for a single inclusion with center $z_i$ we can now use the asymptotic expansion (2) and The inserted material will be chosen out of a set of admissible materials $C$

$$ \begin{aligned} \min_{C^0 \in C} \mathcal{J}_{c(C^0)}(f) &\approx \min_{C^0 \in C} \mathcal{J}_c(f) - \frac{h^2}{2} \varepsilon(u; z_i)^\top P_{(C^0)} \varepsilon(u; z_i) \\ &= \mathcal{J}_c(f) + \frac{h^2}{2} \min_{C^0 \in C} (-\varepsilon(u; z_i)^\top P_{(C^0)} \varepsilon(u; z_i)). \end{aligned} $$

Thus, taking into account that $J_c(f)$ is independent of the rotation angle and noting that $h^2$ is a constant scaling parameter, the functional

$$ \mathcal{D}_c(C^0, z_i) := -\varepsilon(u; z_i)^\top P_{(C^0)} \varepsilon(u; z_i) \quad (6) $$

can be used as an approximate model to find the optimal rotation angle resulting in the optimization problem

$$ \min_{C^0 \in C} \mathcal{D}_c(C^0, z_i). \quad (7) $$

We note that the functional $\mathcal{D}_c$ depends on the inserted material $C^0$ only via the polarization matrix $P_{(C^0)}$ and on the displacement field $u(x)$ only locally through the evaluation of the strain at center $z_i$ of the ellipse. Furthermore, only a single evaluation of the state problem for the unperturbed domain is required.

Due to the local character of (6) the optimal orientation of $n_{ell} < \infty$ non-intersecting inclusions at once can be approximated by the solution of $n_{ell}$ optimization problems of type (7) or equivalently by the solution of the problem

$$ \min_{(C_0^0, \dots, C_{n_{ell}}^0) \in C^{n_{ell}}} \mathcal{J}_{c(C_0^0, \dots, C_{n_{ell}}^0)}^{\text{asymp}} := \sum_{i=1}^{n_{ell}} \min_{C_i^0} D_c(C_i^0, z_i), \quad (8) $$

which is separable in terms of the different materials $C_0^0, \dots, C_{n_{ell}}^0$ used in the different inclusions. The latter is certainly only an approximation of the original simultaneous material optimization problem

$$ \min_{(C_0^0, \dots, C_{n_{ell}}^0) \in C^{n_{ell}}} \mathcal{J}_c(C_0^0, \dots, C_{n_{ell}}^0), \quad (9) $$
however it is shown in [6] for the optimal rotation of an orthotropic material, that the results for this approximation are close to the solution of the original problem. Moreover, we will later on compare the quality of the approximated solution to rigorous lower bounds. We want to stress that still only a single evaluation of the state problem for the unperturbed domain is required, which allows for a highly efficient numerical solution.

The optimization procedure for solving (8) is performed by the following algorithm:

(S1) choose matrix material $C^1$ and admissible materials for inclusions $C$;

(S2) compute $\Psi$ from (4);

(S3) define loads and boundary conditions;

(S4) solve state problem without inclusions for isotropic material;

foreach inclusion $\omega_i$, $i = 1, \dots, n_{ell}$ do

(S5) solve $\min_{C^0 \in C} \mathcal{D}_c(C^0, z_i)$ to global optimality;

end

**Algorithm 1:** Basic algorithm for minimization of the compliance based on the asymptotic model.

## 2.2 Admissible material choices

In order to avoid local minima, the set of admissible materials $C$ should allow for a global solution of the local material optimization problem (7) in a reasonable time. Interestingly, the easiest choice here would be a set of discrete materials. In the following, however, we concentrate on parametric material formulations. Using the properties of the asymptotic expansion, we can give a wider class of parametrizations that lead to globally solvable problems. According to [6] (Theorem 2.8, p. 14), the polarization matrix (3) is positive definite if $(C^0)^{-1} - (C^1)^{-1}$ is negative definite. If the isotropic material $C^1$ is then chosen s.t. every $C^0 \in C$ is strictly stiffer, it can then be shown that the functional (6) is convex for any linear material parametrization. An example of this would be the so-called free material optimization (FMO, see e.g. [8, 9]):

$$ C_{\text{FMO}} := \left\{ \begin{pmatrix} e_1 & e_2 & e_3 \\ e_2 & e_4 & e_5 \\ e_3 & e_5 & e_6 \end{pmatrix} \succeq \bar{\tau} \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix}, e_1 + e_4 + e_6 = \bar{\rho} \right\}, \quad (10) $$

where $\bar{\tau}$ is a fixed lower eigenvalue bound and $\bar{\rho}$ bounds the total stiffness of the material tensor. The positivity constraint is still linear as a semidefinite program, s.t. the resulting problem stays convex.

Furthermore, we will now consider the usually difficult problem of rotational optimization that is also vastly simplified by the breakdown to local minimization problems. An example of this would be the variation of angle and stiffness

$$ C^0(\theta, s) := \Theta(\theta)^{\top} C(s) \Theta(\theta), \quad C(s) := \begin{pmatrix} 1.0 + s & 0.5 & 0 \\ 0.5 & 15 - s & 0 \\ 0 & 0 & 1.0 \end{pmatrix} \quad (11) $$
where we rotate an orthotropic material by an angle $\theta$ using the orthogonal rotation matrix $\Theta$ defined as

$$ \Theta(\theta) = \begin{pmatrix} \cos(\theta)^2 & \sin(\theta)^2 & -\frac{\sqrt{2}}{2}\sin(2\theta) \\ \sin(\theta)^2 & \cos(\theta)^2 & \frac{\sqrt{2}}{2}\sin(2\theta) \\ \frac{\sqrt{2}}{2}\sin(2\theta) & -\frac{\sqrt{2}}{2}\sin(2\theta) & \cos(2\theta) \end{pmatrix}. \quad (12) $$

In the results section, we will consider two different admissible material sets based on this parametrization, namely

$$ C_{\theta,s} := \{C^0(\theta, s) : \theta \in [0, \pi], s \in [0, 14]\} \quad (13) $$

and a pure rotational optimization

$$ C_{\theta,s=0} := \{C^0(\theta, s) : \theta \in [0, 2\pi], s = 0\}. \quad (14) $$

Now this parametrization fails the proposed linearity. However, for any fixed rotation angle $\theta$, the local minimization problems (8) are still strictly convex. It follows, that the problem can be solved as a bilevel problem

$$ \min_{\theta_i, s_i} \mathcal{J}_c(C_i^0(\theta_i, s_i)) = \min_{\theta_i} G(\theta_i) $$

$$ G(\theta_i) := \min_{s_i} \mathcal{J}_c(C^0(\theta_i, s_i)), $$

which still allows for a global solution with moderate cost as there is only a single primary variable left. Note that this could be done similarly for more complicated linear parametrizations with rotation or in 3D with 2-3 rotation angles.

Lastly, we consider an orthotropic material parametrization using engineering constants with fixed Poisson's ratio $\nu_{xy} = \nu_{yx} = 0.3$:

$$ C(\theta, E_x, E_y, G_{xy}) = \Theta(\theta)^{\top} \begin{pmatrix} \frac{E_x}{1-0.09} & \frac{\sqrt{0.09 E_x E_y}}{1-0.09} & 0 \\ \frac{\sqrt{0.09 E_x E_y}}{1-0.09} & \frac{E_y}{1-0.09} & 0 \\ 0 & 0 & 2G_{xy} \end{pmatrix} \Theta(\theta), \quad (15) $$

with elasticity moduli $E_x, E_y$, shear modulus $G_{xy}$ and $\Theta$ as in (12). We define the admissible material set corresponding to (13) as

$$ C_{\text{Eng}} := \{C(\theta, E_x, E_y, G_{xy}) : \theta \in [0, \pi], E_x, E_y \in [1, 15], G_{xy} \in [0.5, 7.5]\}. $$

## 2.3 From elliptic inclusions to material optimization

While the asymptotic model (2) rigorously holds only for elliptic inclusions of small size, in the following we will also numerically investigate the behavior when replacing the material inside squared patches of finite elements. Choosing the elements properly, the material in the whole domain can be replaced this way with the FE patches still being disjoint. Using a large number of patches, the size of the inclusion stays small compared to the domain size. Thus, we will study increasingly bigger ellipses and compare the compliance values of the different parametrizations to global lower material optimization bounds computed with an FMO solver.
**Validation methods** For the numerical evaluation, we discretize the domain $\Omega$ using rectangular finite elements. This discretization is necessary to compute the displacements used in the asymptotic expansion. The elliptic inclusions are approximated by those finite elements, for which the coordinates of their center point are contained in the inclusion $\omega_i$. In order to obtain the actual compliance value for the optimization result, the material used in those elements is then replaced by the optimal value of $C_i^0$. When replacing all material, we use equally sized squared FE patches that are uniformly distributed, disjoint and cover the whole domain.

Furthermore, we compare the compliance values to a global lower material optimization bound determined by solving a modified FMO problem. Specifically, we solve

$$ \min_{(C_0^0, \dots, C_{n_{\text{ell}}}^0) \in \mathcal{C}_{\text{FMO}}^{n_{\text{ell}}}} \mathcal{J}_c(C_0^0, \dots, C_{n_{\text{ell}}}^0) \quad (16) $$

with $C_{\text{FMO}}$ as in (10) and the bounds $\underline{\tau}$ and $\bar{\rho}$ chosen as close as possible to the ones used in the specific material parametrization, s.t. all possible tensors are a subset of $C_{\text{FMO}}$. For a more detailed description of the method, see [6]. Note that within these bounds, any physically admissible material may be used and that this problem is convex for the compliance cost functional. Thus, the problem is solved globally using the algorithm described in [9] and we obtain a global lower compliance bound.

**Numerical results** We consider the example from Fig. 1 with 10 × 10 ellipses and discretize the domain $\Omega$ using 100 × 100 finite elements. We compare the different admissible material sets as defined in Section 2.2. For $C_{\text{FMO}}$ we choose $\underline{\tau} = 1$ as lower eigenvalue bound and $\bar{\rho} = 17$ as upper trace bound both in the asymptotic material optimization as in the FMO solver. The results are shown in Table 1 and a visualization in Fig. 2. Although the error compared to the exact FMO result increases heavily with the size of the inclusions, this is largely due to the decrease of the overall compliance value. The absolute value does not increase much from the largest ellipses to the squared FE patches.

<table><caption>**Table 1.** Compliance values and FMO comparison for increasing ellipse size.</caption><thead><tr><th>a = b:</th><th>0.02</th><th>0.04</th><th>0.05</th><th colspan="2">squared patch</th></tr></thead><tbody><tr><td>C<sub>θ,s=0</sub></td><td>15.420</td><td>1.0%</td><td>11.053</td><td>6.15%</td><td>8.3173</td><td>14.2%</td><td>4.3730</td><td>42.3%</td></tr><tr><td>C<sub>θ,s</sub></td><td>15.332</td><td>0.42%</td><td>10.738</td><td>3.12%</td><td>7.7140</td><td>5.93%</td><td>3.5668</td><td>16.1%</td></tr><tr><td>C<sub>Eng</sub></td><td>15.328</td><td>0.39%</td><td>10.647</td><td>2.25%</td><td>7.6135</td><td>4.55%</td><td>3.5014</td><td>13.9%</td></tr><tr><td>C<sub>FMO</sub></td><td>15.289</td><td>0.14%</td><td>10.551</td><td>1.33%</td><td>7.4840</td><td>2.77%</td><td>3.2912</td><td>7.10%</td></tr><tr><td>FMO</td><td>15.268</td><td></td><td>10.413</td><td></td><td>7.2821</td><td></td><td>3.0730</td><td></td></tr></tbody></table>

For the squared FE patches, we furthermore study the different parametrizations separating the domain into 50 × 50 patches. The results are found in Table 2
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Optimization for $C_{\theta,s}$ with similar results. The arrows denote the principal stiffness directions and the gray value the magnitude (black: low, white: high).</figcaption></figure>

and Fig. 3(a). We can see, that for the parametrization with nonlinear subproblems $C_{Eng}$, apparently a local minimum was found, which stresses the importance of a global solution of the local material optimization problems.

<table><caption>**Table 2.** Compliance values and FMO comparison for 50x50 squared FE patches.</caption><tr><td>$C_{\theta,s=0}$</td><td>3.6850</td><td>38.89%</td></tr><tr><td>$C_{\theta,s}$</td><td>2.9963</td><td>12.93%</td></tr><tr><td>$C_{Eng}$</td><td>3.0535</td><td>15.08%</td></tr><tr><td>$C_{FMO}$</td><td>2.8964</td><td>9.16%</td></tr><tr><td>FMO</td><td>2.6533</td><td></td></tr></table>

## 2.4 Material and topology optimization

The term (6) in the asymptotic expansion actually corresponds to a topological derivative and can hence be used for topology optimization. We propose an algorithm, where we do not only use the information of “drilling a hole”, but use the values obtained in material optimization instead. The idea is, that those finite element patches, which provide the smallest gain, when substituting the optimal material $C^0$, are left out in the following iteration:
(S1) choose matrix material $C^1$ and admissible materials for inclusions $C$;

(S2) compute $\Psi$ from (4);

(S3) define loads and boundary conditions;

while volume > bound do

(S4) solve state problem without inclusions for current topology;
foreach inclusion $\omega_i$, $i = 1, \dots, n_{ell}$ do
| (S5) solve $\min_{C_i^0 \in C} D_c(C_i^0, z_i)$ to global optimality;
end
(S6) remove FE-patches with $\min_{C_i^0 \in C} D_c(C_i^0, z_i)$ largest;

end

**Algorithm 2:** Basic algorithm for simultaneous material and topology optimization based on the asymptotic model.

Note that, again, in each iteration the state problem only needs to be solved once.

In Fig. 3(b), the result of the algorithm for the admissible material set $C_{\theta,s}$ and $50 \times 50$ squared FE patches is shown. In this experiment, we removed in the iteration $k$ a total of $200 \times 0.83^k$ FE patches until 3 FE patches or less where removed, which lead to a final volume fraction of 0.542. The strategy for the removal of ellipses can be varied, however a decreasing volume fraction should be removed in order to obtain a smoother convergence.

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Results of the simultaneous material and topology optimization for $C_{\theta,s}$.</figcaption></figure>

Conclusion

We proposed an efficient algorithm for material optimization on multiple elliptic inclusions. The numerical evidence suggests that the accuracy of the proposed
method decreases only slightly, when replacing all material instead of just the
material in elliptic inclusions. A major advantage of this algorithm is the possi-
bility of avoiding local minima, however at the cost of only having an approx-
imate solution. The total error in the studied example for the finer resolution
was about 10%. From experience in practice, the proposed algorithm for simul-
taneous material and topology optimization seems to work well in large parts,
however oftentimes small bars are left over and if a hole happens to be drilled
in a “bad” position the algorithm struggles, as material is not reintroduced.

Nevertheless, both algorithms allow for a very efficient solution of usually quite complicated problems, such as discrete material optimization and rotational optimization. When the accuracy provided in this method does not suffice or the optimized topology appears flawed, the optimization result may still be used as a high quality initial design for other solution schemes, such as fully parametrized approaches.

References

1. Leugering, G., Nazarov, S., Schury, F., Stingl, M.: The Eshelby theorem and application to the optimization of an elastic patch. SIAM J. Appl. Math. **72**(2) (2012) 512–534

2. Pedersen, P.: On optimal orientation of orthotropic materials. Struct. Multidiscip. Opt. **1** (1989) 101–106

3. Pedersen, P.: On thickness and orientational design with orthotropic materials. Struct. Multidiscip. Opt. **3** (1991) 69–78

4. Thomsen, J.: Optimization of composite discs. Struct. Multidiscip. Opt. **3** (1991) 89–98

5. Lewiński, T., Sokołowski, J.: Energy change due to the appearance of cavities in elastic solids. Int. J. Solids Struct. **40**(7) (2003) 1765–1803

6. Schury, F., Greifenstein, J., Leugering, G., Stingl, M.: On the efficient solution of a patch problem with multiple elliptic inclusions. Optim. Eng. (2014) Accepted for publication.

7. Love, A.E.H.: A Treatise on the Mathematical Theory of Elasticity. Dover Publications, New York (1944)

8. Haslinger, J., Kočvara, M., Leugering, G., Stingl, M.: Multidisciplinary free material optimization. SIAM J. Appl. Math. **70**(7) (2010) 2709–2728

9. Stingl, M., Kočvara, M., Leugering, G.: A sequential convex semidefinite programming algorithm with an application to multiple-load free material optimization. SIAM J. Opt. **20**(1) (2009) 130–155