30.07.2013 Views

WAVES AND VIBRATIONS IN INHOMOGENEOUS STRUCTURES ...

WAVES AND VIBRATIONS IN INHOMOGENEOUS STRUCTURES ...

WAVES AND VIBRATIONS IN INHOMOGENEOUS STRUCTURES ...

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.

<strong>IN</strong>TERNATIONAL JOURNAL FOR NUMERICAL METHODS <strong>IN</strong> ENG<strong>IN</strong>EER<strong>IN</strong>G<br />

Int. J. Numer. Meth. Engng 2007; 72:1605–1630<br />

Published online 10 April 2007 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/nme.2065<br />

Topology optimization of dynamics problems<br />

with Padé approximants<br />

Jakob S. Jensen ∗,†<br />

Department of Mechanical Engineering, Solid Mechanics, Technical University of Denmark,<br />

Nils Koppels Allé, Building 404, Denmark<br />

SUMMARY<br />

An efficient procedure for topology optimization of dynamics problems is proposed. The method is based<br />

on frequency responses represented by Padé approximants and analytical sensitivity analysis derived using<br />

the adjoint method. This gives an accurate approximation of the frequency response over wide frequency<br />

ranges and a formulation that allows for design sensitivities to be computed at low computational cost<br />

also for a large number of design variables. Two examples that deal with optimization of forced vibrations<br />

are included. Copyright q 2007 John Wiley & Sons, Ltd.<br />

Received 30 October 2006; Revised 1 February 2007; Accepted 2 March 2007<br />

KEY WORDS: Padé approximants; topology optimization; forced vibrations<br />

1. <strong>IN</strong>TRODUCTION<br />

The solution to a discretized linear vibration problem with harmonic excitation can be expressed<br />

in the exact form<br />

u =<br />

2Nd−2 i=0<br />

ãi( − 0) i<br />

2Nd i=0 ˜bi( − 0) i<br />

in which Nd is the number of degrees of freedom in the model, ãi and ˜bi are expansion coefficients,<br />

is the excitation frequency, and 0 is an arbitrary expansion frequency. This representation of<br />

the solution is rarely used except for systems with only few degrees of freedom for which the<br />

corresponding small number of expansion coefficients can be found analytically. Solutions for larger<br />

problems are usually obtained numerically with a factorization method or an iterative method.<br />

∗ Correspondence to: Jakob S. Jensen, Department of Mechanical Engineering, Solid Mechanics, Technical University<br />

of Denmark, Nils Koppels Allé, Building 404, Denmark.<br />

† E-mail: jsj@mek.dtu.dk<br />

Copyright q 2007 John Wiley & Sons, Ltd.<br />

(1)

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

Saved successfully!

Ooh no, something went wrong!