3D mechanic compliance minimization

The simulations in 3d look much better, but still need some hours to calculate. I hope to speed this up by a factor of 5-10, I will need it for real piezo optimization.

This is a bridge, fixed at the four lower corner points and three load points in the lower plane. I allowed half of the total volume. It is quite symmetric. The optimization is done with the optimiality criteria and the iterations for all this examples are between 15 and 40.


Here the bridge has only a single lower load point.


Here the single load point example but a volume fraction of 10%.


Here Tom (thanks again :)) suggested to optimize the eiffel tower. The top center force in positive z-direction, the lower four corners are fixed.


This is the most challenging example. For this bridge I apply a downwards force on every node on the lower plane. Currently it is still running, so this is a snapshot of iteration 15. The volume fraction is 20% and the penalization is 5 instead of 3. Still I don't think that there will ever be a nice black and white result. The green is for a pseudo density of 0.5. I'm curious if you have any idea to solve this optimization problem. I thought about adaptive penalization (increasingly higher depending on the greyness) - does anyone have experience in this field?


Finally the aqueduct which looks much better in 3d :) Here I use a penalization of 4. With the standard 3 I get nice big holes but I cannot get rid of the greyness. At leas here higher penalization works. The load applied on the upper long central line (a "pipe"). I fix the upper four corners and the front end back point of the feet. Note that there is a bead!