Automatic Differentiation for a Structural Optimization Solver.pdf

Automatic Differentiation for a Structural Optimization Solver Presented to EUCCO 2004 - European Conference on Computational Optimization, Minisymposia on Automatic Differentiation and Large-Scale Optimization, March 29-31 2004, Technical University of Dresden, Germany. Shaun Forth Mohamed Tadjouddine Andy Keane Applied Maths. OR Group Computational Engineering Design Centre Engineering Systems Department School of Engineering Sciences Cranfield University (Shrivenham Campus) University of Southampton, Highfield Swindon SN6 8LA, UK Southampton SO17 1BJ, UK Work funded by EPSRC under Grant GR/R85358/01 AD2CompEng - Automatic Differentiation and Adjoints Applied to Computational Engineering Structural Design – p.1/19 Introduction Lightweight cantilever structure for suspending scientific instruments away from satellite. Wish to minimise trans- mission of vibration through structure from satellite to instrument NASA Photo ID: STS61B-120-052 Structural Design – p.2/19 Model Using Receptance Theory [SK95] Structure of long, thin beams executing transverse purely longitudinal Euler-Bernoulli modes. Beams coupled at their ends or ends are hinged/clamped. Formulate beam end displacements in terms of applied forces via Green functions. Leads to large linear system whose solution gives forces at end of each beam. Vibrations- Consider forcing x1,1 = ?x1,1e iωt at first end of beam 1 Solve complex linear system for amplitude and phase of structure response xn,j at end j of beam n. From force fn,j at end j of node n calculate power P == 1 2 Re(fn,j · xn,j) Structural Design – p.3/19 Optimisation by Genetic Algorithms 6 ? 1) Periodic forcing at this end 2) Measure transmitted power at this end        * 3) Use GA to minimise power Structural Design – p.4/19 Reduction in Transmitted Power (2D) Structural Design – p.5/19 Optimisation by Genetic Algorithm Not a huge problem n = 453 Use of GA with population of 100 gives slow convergence Run times of 83 CPU days for 10 ge