Maria Bayard Dühring - Solid Mechanics
Maria Bayard Dühring - Solid Mechanics
Maria Bayard Dühring - Solid Mechanics
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
20 Chapter 3 Topology optimization applied to time-harmonic propagating waves<br />
1. Choose target frequencies ωi and initialize the design variable vector ξ, the<br />
counter iter = 0, the tolerance tol=0.01 and change=1.<br />
2. for iter = 1:1000<br />
3. Calculate filtered densities ˜ ξ by the chain rule.<br />
4. Calculate the material values and solve the FE problem based on ˜ ξ.<br />
5. Compute the objective function Ψ and the sensitivities based on ˜ ξ.<br />
6. Calculate filtered sensitivities.<br />
7. Update the design variables ξ new by the MMA-algorithm.<br />
8. Compute change=max(||ξ new -ξ||).<br />
9. If mod(iter,50) == 0 | change < tol, then update target frequencies, the material<br />
damping or the filter as appropriate and set change=0.5.<br />
10. Break if change < tol.<br />
11. end<br />
12. Postprocess and display the results.<br />
The optimized designs are in general sensitive to the choice of driving frequency<br />
or frequency interval as well as other factors such as the volume fraction β, the initial<br />
guess and the filter size. It is therefore recommended to choose the driving frequency<br />
or the frequency interval carefully and to perform more optimizations with several<br />
parameter combinations to obtain the best possible design.