1 Introduction to LS-OPT – 14.05.08
DYNAmore GmbH Industriestraße 2 70565 Stuttgart http://www.dynamore.de
Heiner Müllerschön
hm@dynamore.de
Introduction to LS-OPT Heiner Mllerschn hm@dynamore.de DYNA more - - PowerPoint PPT Presentation
Introduction to LS-OPT Heiner Mllerschn hm@dynamore.de DYNA more GmbH Industriestrae 2 70565 Stuttgart http://www.dynamore.de Introduction to LS-OPT 14.05.08 1 Overview n Introduction/Features n Methodologies Optimization n
1 Introduction to LS-OPT – 14.05.08
DYNAmore GmbH Industriestraße 2 70565 Stuttgart http://www.dynamore.de
Heiner Müllerschön
hm@dynamore.de
2 Introduction to LS-OPT – 14.05.08
3 Introduction to LS-OPT – 14.05.08
stand alone optimization software
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.3 / Outlook
4 Introduction to LS-OPT – 14.05.08
variables (e.g. sheet thicknesses)
robustness of the optimum
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.3 / Outlook
5 Introduction to LS-OPT – 14.05.08
(abnormal termination)
as a “Process Manager”
HyperMorph, DEP-Morpher, SFE-Concept
any Pre-Processor
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.3 / Outlook
6 Introduction to LS-OPT – 14.05.08
(*PARAMETER_)
progress
LS-DYNA response types
(node and part selection)
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.3 / Outlook
7 Introduction to LS-OPT – 14.05.08
2 1
= P p i i i i
calibrate materials/systems
identification applications
and simulation curve to compare
individual parameters in parameter identification
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.3 / Outlook
8 Introduction to LS-OPT – 14.05.08
Objective D e s i g n V a r i a b l e 1 D e s i g n V a r i a b l e 2
Subregion (Range) Starting (base) design Response surface Response values Experimental Design points
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.3 / Outlook
9 Introduction to LS-OPT – 14.05.08
Optimization of sub-problem (response surface) using LFOPC algorithm Optimum (predicted by response surface) Optimum (computed by simulation using design variables) Starting value on response surface
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.3 / Outlook
10 Introduction to LS-OPT – 14.05.08
Design Variable 1 Design Variable 2
Region of Interest Design Space
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.3 / Outlook
11 Introduction to LS-OPT – 14.05.08
2 2 4 2 1 2 1 2 1 2 1 1 2
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.3 / Outlook
movie
12 Introduction to LS-OPT – 14.05.08
T
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.3 / Outlook
13 Introduction to LS-OPT – 14.05.08
linear polynomial neural network quadratic polynomial
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
14 Introduction to LS-OPT – 14.05.08
2 2 4 2 1 2 1 2 1 2 1 1 2
analytical function (green) global neural net approx. with 20 points (red) simulation points
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
movie
15 Introduction to LS-OPT – 14.05.08
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
16 Introduction to LS-OPT – 14.05.08
results of optimization or stochastic analysis
(property of DYNAmore)
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
feasible
response surface
17 Introduction to LS-OPT – 14.05.08
Gradient Based Methods Random Search Genetic Algorithms RSM / SRSM
§accuracy of
solution
§number of solver
calls
§very robust, can not
diverge
§easy to apply §good for problems
with many local minimas
§very effective,
particularly SRSM
§trade-off studies
§filter out noise,
smoothing of results
§can diverge §can stuck in local
minimas
§ step-size
dilemma for numerical gradients
§bad convergence,
not effective
§Chooses best
may not be representative of a good (robust) design
§many solver calls,
solver runs
§Chooses best
may not be representative of a good (robust) design
§approximation
error
§verification run
might be infeasible
§number of
variables control minimum number of required runs
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
18 Introduction to LS-OPT – 14.05.08
due to Variation of Input (Parameter)
respect to Responses
parameter parameter
response response
Feasible range
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
19 Introduction to LS-OPT – 14.05.08
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
20 Introduction to LS-OPT – 14.05.08
dimension of the problem (number of variables) and the type
contributions clearly
Multi Meta-Model exploration with D-SPEX
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
21 Introduction to LS-OPT – 14.05.08
due to Variation of Input (Parameter)
respect to Responses
parameter parameter
response response
Feasible range
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
22 Introduction to LS-OPT – 14.05.08
important variables input
input
Correlation Matrix ANOVA
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
23 Introduction to LS-OPT – 14.05.08
due to Variation of Input (Parameter)
respect to Responses
parameter parameter
response response
Feasible range
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
24 Introduction to LS-OPT – 14.05.08
1 2 3 4 5 6 7 8 1 1,05 1,1 1,15 1,2 1,25 1,3 1,35
Verteilung max RWFOCE
Probability of 8,4% for violating the FORCE-bound
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
25 Introduction to LS-OPT – 14.05.08
due to Variation of Input (Parameter)
respect to Responses
parameter parameter
response response
Feasible range
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
26 Introduction to LS-OPT – 14.05.08
High Variance
Buckling mode B
Buckling mode A
Courtesy DaimlerChrysler
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
27 Introduction to LS-OPT – 14.05.08
4 cpus, distributed memory
belted / not belted and “Hybrid III 5th Female“ / „Hybrid III 50th Male“ possible
Dummy 56 km/h – belted 40 km/h – not belted Hybrid III 5th Female H305a(ktiv) H305p(assiv) Hybrid III 50th Male H350a(ktiv) H350p(assiv)
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
28 Introduction to LS-OPT – 14.05.08
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
Adaptive Airbag Deployment (6 Variables)
H305a
5%-dummy, belted
H305p
5%-dummy, not belted
H350a
50%-dummy, belted
H350p
50%-dummy, not belted
Area Venthole1 FAB_VENT FAB_VENT FAB_VENT FAB_VENT
Lower – Upper B.
Area Venthole2 SBA_VENT SBA_VENT SBA_VENT SBA_VENT
Lower – Upper B.
Trigger Time FAB_ADT1_05a FAB_ADT1_05p FAB_ADT1_50a FAB_ADT1_50p
Lower – Upper B.
Venthole 1 FAB_VENT Venthole 2 SBA_VENT Trigger Time FAB_ADT1_05a FAB_ADT1_50a FAB_ADT1_05p FAB_ADT1_50p total Venthole Area time [t]
29 Introduction to LS-OPT – 14.05.08
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
Adaptive Steering Column (5 Variables)
H305a
5%-dummy, belted
H305p
5%-dummy, not belted
H350a
50%-dummy, belted
H350p
50%-dummy, not belted
Force Level StCo LKS_SKAL LKS_SKAL LKS_SKAL LKS_SKAL
Lower – Upper Bound
Trigger Time LKS_DOWN05a LKS_DOWN50a LKS_DOWN05p LKS_DOWN50p
Lower – Upper Bound
Force Level StCo LKS_SKAL time [t] Force Steer. Colu.
30 Introduction to LS-OPT – 14.05.08
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
Adaptive Steering Column (5 Variables)
H305a
5%-dummy, belted
H305p
5%-dummy, not belted
H350a
50%-dummy, belted
H350p
50%-dummy, not belted
Force Level StCo LKS_SKAL LKS_SKAL LKS_SKAL LKS_SKAL
Lower – Upper Bound
Trigger Time LKS_DOWN05a LKS_DOWN50a LKS_DOWN05p LKS_DOWN50p
Lower – Upper Bound
Force Level StCo LKS_SKAL time [t] Force Steer. Colu.
31 Introduction to LS-OPT – 14.05.08
–> min BrustA3ms-05a –> min BrustA3ms-50a –> min BrustA3ms-05p –> min BrustA3ms-50p
–> HIC15-05a –> HIC15-50a –> HIC15-05p –> HIC15-50p
–> FemurLi-05a –> FemurLi-50a –> FemurLi-05p –> FemurLi-50p
–> BrustSx-05a –> BrustSx-50a –> BrustSx-05p –> BrustSx-50p
–> BrustA3ms-05a –> BrustA3ms-50a –> BrustA3ms-05p –> BrustA3ms-50p
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
32 Introduction to LS-OPT – 14.05.08
H305a.pc H305p.pc H350a.pc H350p.pc
Input Files
PAM-CRASH LSF EVALUATOR
HIC15-05a HIC15-05p HIC15-50a HIC15-50p ThoraxSx-05a ThoraxSx-05p ThoraxSx-50a ThoraxSx-50p Femur-05a Femur-05p Femur-50a Femur-50p Thoraxa3ms-05a Thoraxa3ms-05p Thoraxa3ms-50a Thoraxa3ms-50p GUR_ENDE05a FABADT1_05a LKS_DOWN05a FABADT1_05p LKS_DOWN05p FABADT1_50p LKS_DOWN50p GUR_ENDE50a FABADT1_50a LKS_DOWN50a GUR_FOR1 Variables local remote remote LKS_SKAL FAB_VENT SBA_VENT
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
33 Introduction to LS-OPT – 14.05.08
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
Airbag and Seatbelt
(active=6, passive=3)
Surface Methodology (SRSM) using linear polynomial approximations
8 iterations - total runs: 276
Design Space
Optimum Start 2 3 global region of interest chosen
34 Introduction to LS-OPT – 14.05.08
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
a result which meets all requirements is gained in 8 iterations, each with 34 shots 5% aktive 50% passive 5% passive 50% aktive
History of Thorax Acceleration
35 Introduction to LS-OPT – 14.05.08
2 1 3 1
= =
j
i n j j j i α
t N m m
m
β − =
1
; N=1, 2, 3, 4, 5, 6 ;
kl t ijkl ij
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
36 Introduction to LS-OPT – 14.05.08
quasi-static curve – > Ogden fit Strain rate A –> fit of Prony terms 2 1
P p
=
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
37 Introduction to LS-OPT – 14.05.08
start design – no beads
displacement force
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
38 Introduction to LS-OPT – 14.05.08
to consider uncertainties
mean = 2.5mm; σ = 0.05mm
mean = 2.4mm; σ = 0.05mm
mean = 0.8mm; σ = 0.05mm
mean = 1.0mm; σ = 0.05mm
mean = 1.0 ; σ = 0.05
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
39 Introduction to LS-OPT – 14.05.08
internal energy sheet thickness T_1139
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
40 Introduction to LS-OPT – 14.05.08
Probability of 8,4% for violating the RWFORC-bound
constraint-bound
0,2 0,4 0,6 0,8 1 1 ,2 1 ,4 1 ,6 1 20 30 40 50 60 70
Max Min Mean
and mean history values
history values
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
41 Introduction to LS-OPT – 14.05.08
Dashboard young_alu x_transl z_transl Airbag Mass Flow scal_massflow Slip Ring Friction sfric1 Slip Ring Friction sfric2 Steering Wheel rot_stwh Pre-Tensioner preten Force Limit Retractor forcelimit Sled Acceleration scalaccel
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
42 Introduction to LS-OPT – 14.05.08
x
scalaccel sfric1 sfric2 preten forcelimit rot_stwh transl_x transl_z scalmassflow young_alu all variables residuals Total Design Variable 2,5% 25,0% 25,0% 4,4% 5,6% 4,8% 50,0% 50,0% 5,0% 5,0% S tandard Deviation
max_bf_pelvis 2,3% 1,8% 3,7% 1,1% 0,6% 0,0% 4,5% 1,6% 2,2% 0,5% 7,2% 6,0% 9,4% Standard Deviation Contribution HIC 36 3,1% 1,3% 0,5% 0,0% 1,3% 0,5% 0,1% 1,2% 1,8% 0,3% 4,3% 4,7% 6,4% max_chest_intru 1,5% 0,6% 0,6% 0,5% 0,4% 0,1% 0,1% 1,0% 1,8% 0,3% 2,8% 1,9% 3,4% max_b_f_shoulder 0,1% 4,1% 0,1% 0,0% 4,4% 0,1% 0,7% 0,3% 0,6% 0,0% 6,1% 1,8% 6,3% max_chest 1,9% 0,7% 0,1% 0,3% 1,4% 0,1% 0,5% 0,2% 0,6% 0,1% 2,6% 3,5% 4,3% max_pelvis 2,9% 0,7% 0,1% 0,2% 0,2% 0,1% 0,8% 0,9% 0,9% 0,1% 3,4% 2,3% 4,1%
Contribution of
results
2 . Determ
σ
Dσ
Dσ
Varσ
Varσ
Rσ
2 Residual
σ
Rσ
2 Total
σ
Tσ
Tσ
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.1 / Outlook
43 Introduction to LS-OPT – 14.05.08
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.2 / Outlook
44 Introduction to LS-OPT – 14.05.08
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.3 / Outlook
45 Introduction to LS-OPT – 14.05.08
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.3 / Outlook
46 Introduction to LS-OPT – 14.05.08
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.3 / Outlook
47 Introduction to LS-OPT – 14.05.08
§ Introduction/Features § Methods – Optimization § Methods - Robustness § Examples - Optimization § Examples - Robustness § Version 3.3 / Outlook
48 Introduction to LS-OPT – 14.05.08