Exact analytical solutions for some popular benchmark problems in topology optimization III: L-shaped domains revisited

The paper provides the corrections and extensions of the paper of the same title, published in vol. 35, 2008, pp. 165–174 of the current journal. The analytical solutions of Michell trusses constructed within the L-shaped domains are confirmed by the ground structure method of the truss optimization, thus giving the upper bounds of the exact solutions.

domains of polygons, allowing also for non-convex shapes. The current version of the program developed makes use of the adaptive ground structure approach with the selective subsets of active bars (Sokół 2011b). This approach is similar in spirit to that proposed by Gilbert and Tyas (2003) but applies a different strategy of adding-removing of the active bars in the fully connected ground structure. The program is capable of solving large-scale problems with the number of potential members exceeding one billion (10 9 ) which provides highly accurate results both for volume and layout of fully stressed truss of minimal weight, designed within a given feasible domain, supported along given parts of the boundary or at fixed supports and subjected to a loading applied at fixed points. In the present paper we shall confine attention to the case of the absolute value of the lower and upper limits for the stress being equal, denoted by σ 0 . Equipped with this new software one can solve numerically a broad class of the optimum design problems of Michell type for the L-shaped domains thus checking the analytical solutions predicted by Lewiński and Rozvany (2008). The noted discrepancies between the numerical predictions and the analytical results mentioned were the impulse for the present authors to look at this class of problems more closely. The result of the current verification is the confirmation of all the layouts predicted in Lewiński and Rozvany (2008) in the case of the point load being directed vertically (or parallel to the vertical edges of the L-shaped domain). However, the formulae for the optimal weight should be corrected. The source of the mistake was neglecting the rigid rotation ω at the vicinity of the re-entrant corner. The present paper makes up for this shortcoming, by deriving the new and now fully correct formulae for the virtual (adjoint) displacement fields within the feasible L-shaped domains, leading to analytical formulae for the optimum weight. The new formulae provide the results which compare favorably with the numerical predictions involving up to 10 9 truss members.
The cases of the force being non-vertical need a special care. The numerical predictions reveal new, till now never published layouts. Some of them lie within the scope of the formulae derived for the case of the vertical forces, while other are completely different, hence their analytical discussion is planned to be the subject of a future research.
For brevity, Eq.n, Fig.m and Sec.p of Lewiński and Rozvany (2007) (or Lewiński and Rozvany 2008) will be referred to as Eq.II.n, Fig.II.m and Sec.II.p (or Eq.III.n, Fig.III.m and Sec.III.p). The notation used is adopted after the papers referred. The present paper can be viewed as the addendum and extension of the publication by Lewiński and Rozvany (2008).

Optimal topology for an inclined line support and point load
We revisit the Section III.2. Note that Eq.II.78 refers to the circular domain in Fig where x(α, β), y(α, β) are given by Eqs.III.3. Thus the position of point P with coordinates (α p , β p ) is determined by with ω 0 = −ω(0, 0) or ω 0 = −d 1 and since d 1 = −1 − π , see the result below Eqs III.6, we have while u 0 , v 0 are defined by Eq.III.11. The results III.12, III.13 should be corrected to the form where k = π . The formula III.15 is replaced by where k = π and the formula III.18 should be corrected accordingly. The rest of Sec.III.2 is correct.

L-shaped domain with a horizontal support
The Sec.III.3 is correct, yet the Table III.1, comprising the results of Sec.III.4 should assume the new form (Table 1). The volume of an optimal truss is equal to V = V V 0 , The

Extended L-shaped domains with an inclined lower boundary
The formula III.25 holds, where and u P is the virtual displacement along the force P , as in Fig. 3. The displacement u P is computed by with Fig. 1 Case of θ = 6π/10, V = 19.1096, and α P = π/6. The analytical layout Substitution of (9) into (8) gives and u, v are given by (5), where one should put k = 2μ. These equations for u and v can be derived from the integral formulae (3), where one should put ω 0 = 1 + 2μ, cf. Eqs.III.24 and II.78.

The cases of the non-vertical point load
The optimal trusses equilibrating given forces and reactions are composed of two families of bars, one of which in tension and the second in compression. All the nodes, including the supports should be equilibrated. Thus the layout in figure III.12 is valid only for the case of the vertical force-parallel to the bars AG 2 and A 1 B. The layouts in Fig. 2 The numerically determined layout for problem of Fig. 1; the volume of the structure equals V num = 19.12494V 0 , cf. Table 1. The result was obtained using the ground structure with 150 mln potential members of The analytical layouts for these two cases (Figs. 4 and 5) will not be constructed, as they lie outside the scope of the layouts discussed here. Instead, let us consider the problem of Fig. 6.
In this problem the reactions H A = P 3 and P 1 are positive. Thus the bars which meet at the node A will have different signs. The analytical layout shown in Fig. 6 will be now explained and parameterized, see Fig. 7. Our main task is to match two virtual displacement fields within the domains: ABG 2 and BG 2 D along the radius BG 2 . Both The conditions: assure that the formulae (11, 12) describe a continuous displacement field. Note that because the sign conventions for the rigid rotation are different for both the parameterizations. The choice d 1 = −1 − 2μ will assure that v A = 0, where v A is denoted in Fig. 7, while u A does not vanish. This choice does not satisfy the kinematic conditions of the problem, but the correct formulae for the optimal weight will be achieved by taking into account the virtual work of the reactions.
Our aim is now to find both virtual displacements (u, v) at arbitrary point (ξ * , η * ) of the domain G 0 AG 2 . Let us compute the free terms in (3): and now compute the displacements: u(ξ * , η * ), v(ξ * , η * ) by using (3), where we use (15) and put ω * 0 = −1. Rather lengthy derivation, with using the integration formulae for Lommel functions lead to and v ξ * , η * = av 1 ξ * , η * , Let us note that where α 0 = ∠G 0 BG 2 = 50.309 • . Note that α 0 is the root of the equation Note, also that the displacements (16), (18) can be equivalently found from (5) by specifying k = −2, because this choice corresponds to a correct value of ω 0 for the parameterization in the domain G 1 BPE 2 G 1 . The weight of the optimal structure is given by V = V 1 + V , where σ 0 V 1 is equal to the virtual work of the forces P 1 , P 3 , P (see Fig. 6) on the virtual displacement field and V is the weight of the bar BA 1 : V = P 2 P V 0 . Thus where u A , v A are given by (20) and where u is defined by (6) for k = 2μ or Let us note that v A = 0. Indeed, by (19) but G 0 = F 0 + F 2 , G 1 = F 1 + F 3 (cf. (7) in Lewiński et al. 1994) hence Recalling the identity (cf. (8) in Lewiński et al. 1994) and taking into account (21) we conclude that v A = 0. The equilibrium equations for the structure of Fig. 6 read P 1 = P [L/a + cos μ], P 2 = P [L/a + sin μ + cos μ], and L = |G 1 P| or, see Eq.III.26 Let us compute u A /a by (23), (17), (21): We compute V by (22), (28), (29), taking into account that v A = 0: Let us define the function We rearrange (31) to the form The above equation is applicable if the net of bars is composed of one family of bars in tension and the second family of bars in compression, mutually orthogonal. Consider now the problem of Fig. 6 for the data: where l is the horizontal coordinate of point P, see Fig. 6. From the equation: we find L = 3.91995a, and by solving (29) we find θ = 113.4456 • . We remember that α 0 = 50.309 • . The volume of the optimal truss is computed by (33), hence V /V 0 = 42.85857. The analytical layout is shown in Fig. 6, while the numerical layout is given in Fig. 8. The numerically

Final remarks
The present paper provides new analytical and numerical solutions to the Michell problem, of transmitting a given point load to given supporting nodes, corresponding to the feasible domain being L-shaped, thus being concave, with one reentrant corner. The solutions displayed are characterized by the following features: -the optimal structures do not make use of the whole feasible domain; some empty domains, around the boundary as well as within the feasible domain appear, see Fig. 5; -the solutions turn out to be highly sensitive to the load direction, see Fig. 1; the domain between the vertical bars will be (partly) filled up if the load deviates slightly from the vertical direction, see Fig. 6.
All analytical solutions presented in this paper have been confirmed with four-digits precision by the numerical solutions corresponding to highly dense ground structures.