engineering optimization concepts and applications fred van keulen matthijs langelaar cla h21.1...

49
Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 [email protected]

Post on 20-Dec-2015

247 views

Category:

Documents


5 download

TRANSCRIPT

Page 1: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Engineering Optimization

Concepts and Applications

Fred van Keulen

Matthijs Langelaar

CLA H21.1

[email protected]

Page 2: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Delft in The Netherlands

Delft

Page 3: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Background

Page 4: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Overview of research projects• Optimization with Uncertainties

• Approximate optimization

• Topology Optimization

• Multilevel optimization

• Fast reanalysis

• Buckling of submarine

• Impregnation

• Shoulder endoprosthesis

• SMA actuators

• Microactuation for Butterfly

• Microactuators (them./electr)

• MEMS packaging

• MEMS surface effects

• MEMS measurement structures

• Electronic interface modeling

• Modeling of MEMS

• MEMS optimization

Page 5: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

5

Man-made Insect

Page 6: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Topology

Optimization

Page 7: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Submarines

Page 8: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Micro actuator● 13 μm, ie 2.5% longitudinal strain

● at 2 V, 27 mW, Tmax = 200C

60 um

530 um

Page 9: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Who are you?

Page 10: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Course Objectives

● Understanding of principles and possibilities of

optimization

● Knowledge of optimization algorithms, ability to choose

proper algorithm for given problem

● Practical experience with optimization algorithms

● Practical experience in application of optimization to

design problems

Page 11: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Course overview

● General introduction, problem formulation, design space / optimization terminology

● Modeling, model simplification

● Optimization of unconstrained / constrained problems

● Single-variable, zeroth-order and gradient-based optimization algorithms

● Design sensitivity analysis (FEM)

● Topology optimization

Page 12: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Course material

● Main text: “Principles of Optimal Design

– Modeling and Computation”, P.Y.

Papalambros & D.J. Wilde, Cambridge

University Press

● Selected topics: “Elements of Structural

Optimization”, R.T. Haftka & Z. Gurdal, Kluwer

Academic Publishers

● Exercises and references

Page 13: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Examination

a) Report on practical exercises using Matlab and Optimization Toolbox (individual or in groups of 2 students)

b) Report on optimization project (individual or in groups of 2 students):

Definition of problem, approach (ca. 1 page A4, Deadline March 28, via email)

Final report

c) Oral exam (individual)

Page 14: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Course Schedule

● No lectures on: 19-2, 11-3 and 1-4

● How to find alternative time slots?

● Training lectures?

Page 15: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

What is optimization?

● “Making things better”

● “Generating more profit”

● “Determining the best”

● “Do more with less”

● Papalambros: “The determination of values for design

variables which minimize (maximize) the objective,

while satisfying all constraints”

Page 16: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Historical perspective

● Ancient Greek philosophers: geometrical optimization

problems

Zenodorus, 200 B.C.:

“A sphere encloses the greatest

volume for a given surface area”

? g

● Newton, Leibniz, Bernoulli, De l’Hospital (1697):

“Brachistochrone Problem”:

Page 17: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Historical perspective (cont.)

● Lagrange (1750): constrained minimization

● Cauchy (1847): steepest descent

● Dantzig (1947): Simplex method (LP)

● Kuhn, Tucker (1951): optimality conditions

● Karmakar (1984): interior point method (LP)

● Bendsoe, Kikuchi (1988): topology optimization

                         

Page 18: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

What can be achieved?

● Optimization techniques can be used for:

– Getting a design/system to work

– Reaching the optimal performance

– Making a design/system reliable and robust

● Also provide insight in

– Design problem

– Underlying physics

– Model weaknesses

Page 19: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Optimization problem

● Design variables: variables with which the design

problem is parameterized:

● Objective: quantity that is to be minimized (maximized)

Usually denoted by:

( “cost function”)

● Constraint: condition that has to be satisfied

– Inequality constraint:

– Equality constraint:

( ) 0g x

( ) 0h x

( )f x

1 2, , , nx x xx

Page 20: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Optimization problem (cont.)

● General form of optimization problem:

xxx

x

xh

xg

xx

nX

f

0)(

0)(

)(

:to subject

min

Page 21: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Solving optimization problems

● Optimization problems are typically solved using an

iterative algorithm:

Model

Optimizer

Designvariables

Constants Responses

Derivatives ofresponses(design sensi-tivities)

hgf ,,

iii x

h

x

g

x

f

,,

x

Page 22: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Curse of dimensionality

Looks complicated … why not just sample the design

space, and take the best one?

● Consider problem with n design variables

● Sample each variable with m samples

● Number of computations required: mn

                                                                    

Take 1 s per computation,

10 variables, 10 samples:

total time 317 years!

Page 23: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Parallel computing

● Still, for large problems,

optimization requires lots

of computing power

● Parallel computing

Page 24: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Optimization in the design process

Conventional design process:

Collect data to describe the system

Estimate initial design

Analyze the system

Check performance criteria

Is design satisfactory?

Change design based on experience /

heuristics / wild guesses

Done

Optimization-based design process:

Collect data to describe the system

Estimate initial design

Analyze the system

Check the constraints

Does the design satisfy convergence criteria?

Change the design using an optimization

method

Done

Identify:1. Design variables2. Objective function3. Constraints

Page 25: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Optimization popularity

● Increasing availability of numerical modeling techniques

● Increasing availability of cheap computer power

● Increased competition, global markets

● Better and more powerful optimization techniques

● Increasingly expensive production processes (trial-and-error approach too expensive)

● More engineers having optimization knowledge

Increasingly popular:

Page 26: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Optimization pitfalls!

● Proper problem formulation critical!

● Choosing the right algorithm

for a given problem

● Many algorithms contain lots

of control parameters

● Optimization tends to exploit

weaknesses in models

● Optimization can result in very sensitive designs

● Some problems are simply too hard / large / expensive

Page 27: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Structural optimization

● Structural optimization = optimization techniques

applied to structures

● Different categories:

– Sizing optimization

– Material optimization

– Shape optimization

– Topology optimization

t

E, R

r

L

h

Page 28: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Shape optimization

Yamaha R1

Page 29: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Topology optimization examples

Page 30: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Classification● Problems:

– Constrained vs. unconstrained

– Single level vs. multilevel

– Single objective vs. multi-objective

– Deterministic vs. stochastic

● Responses:

– Linear vs. nonlinear

– Convex vs. nonconvex (later!)

– Smooth vs. nonsmooth

● Variables:

– Continuous vs. discrete (integer)

Page 31: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Practical example: Airbus A380

● Wing stiffening ribs

of Airbus A380:

● Objective: reduce weight

● Constraints: stress, bucklingLeading edge ribs

Page 32: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Airbus A380 example (cont.)

● Topology and shape optimization

Page 33: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Airbus A380 example (cont.)

● Topology optimization:

● Sizing / shape

optimization:

Page 34: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Airbus A380 example (cont.)

● Result: 500 kg weight savings!

Page 35: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Other examples

● Jaguar F1 FRC front wing:reduce weight

constraints on

max. displacements

5% weight saved

Page 36: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Other examples (cont.)

● Design optimization of packaging products

(Van Dijk & Van Keulen):

● Objective: minimize

material used

● Constraints:

stress, buckling

● Result: 20% saved

Page 37: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

SMA active catheter optimization

Page 38: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

But also …● Optimization is also applied in:

– Protein folding

– System identification

– Financial market forecasting (options pricing)

– Logistics (traveling salesman problem),route planning, operations research

– Controller design

– Spacecraft trajectory planning

● This course: focus on (structural) design optimization

Page 39: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

What makes a design

optimization problem interesting?

● Good design optimization problems often show a

conflict of interest / contradicting requirements:

– Aircraft wing: stiffness vs. weight

– F1 car: idem

– Oil bottle: stiffness / buckling load vs. material usage

● Otherwise the problem could be trivial!

Page 40: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

The optimization model

Model

Optimizer

Designvariables

Constants Responses

Derivatives ofresponses(design sensi-tivities)

hgf ,,

iii x

h

x

g

x

f

,,

x

Page 41: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Systems approach

● Systematic way of thinking:

– What is input / output?

– What belongs to system / environment?

– What level of detail?

– Distinguish sub-systems, hierarchies

System functionInput Output

Environment

Page 42: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Example: cantilever beamh E,

F, U

U(t) F(t)

E, , h, L i

F(t) U(t)Etc.

Page 43: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Model example

Mathematical model:

123

3 3

33

bhE

FL

EI

FLU

Finite element model:FhbLEKU 1),,,(

F, Uh, b

h

b

E, L

Steel

U(x), M(x), V(x)

Page 44: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Model example (2)

● System (state) variables: U(x), M(x), V(x)

● System parameters: h, b, L

● System constants: E,

F, Uh, b

h

b

E, L

Steel

U(x), M(x), V(x)

Page 45: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Features of computer models

● Finite accuracy due to:

– Discretization in time and space

– Finite number of iterations

(eigenvalues, nonlinear models)

– Numerical round-off errors, ill-conditioning

● Responses can be “noisy”:

– Due to different discretization in space and/or time

(e.g. remeshing)

                              

                   

Page 46: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Noisy response

● Example: effect of remeshing

Normalized stress constraint

Hole radius

Page 47: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Features of computer models (cont.)

● Computational models are (very) time consuming

● Often design sensitivities can be calculated

– Cost of design sensitivity analysis?

– Accuracy / consistency of sensitivities

Response

Design variable

ExactNumericalmodel

Page 48: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Finite difference sensitivities

● Straightforward way to compute sensitivities:

finite differences

● More later!

( ) ( )df f x x f x

dx x

Small!

f

x

( )

( )

f x x

f x

x

Page 49: Engineering Optimization Concepts and Applications Fred van Keulen Matthijs Langelaar CLA H21.1 A.vanKeulen@tudelft.nl

Einstein’s advice

“Everything should be made as simple as

possible, but not simpler”

● Model simplification important for optimization!

More in next lectures.