28.08.2013 Views

Maria Bayard Dühring - Solid Mechanics

Maria Bayard Dühring - Solid Mechanics

Maria Bayard Dühring - Solid Mechanics

SHOW MORE
SHOW LESS

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.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!