SlideShare a Scribd company logo
1
Integration of reduced 3D models in
vibration design processes.
Examples from various industries.
Etienne Balmes
SDTools
Arts et Métiers ParisTech
NAFEMS Simulation des systèmes
3 Juin 2015
• FEM simulations
• System models (model reduction, state-space, active control, SHM)
• Experimental modal analysis
• Test/analysis correlation, model updating
Activities
2
CAD/Meshing
FEM
Simulation
Testing
CATIA, Workbench, …
NASTRAN, ABAQUS, ANSYS,...
Adams, Simulink,...
LMS TestLab, ME-Scope, …
Simulation
Validation
SDT : MATLAB based toolbox
Commercial since 1995 > 700 licenses sold
Pantograph/catenary Modal test correlation Track dynamics
Outline
• Systems, models, dynamic models
• Tools for model reduction
– Variable separation
– Parametric models
– Domain decomposition
• Conclusion
3
A system = I/O representation
Prototype Virtual prototype
 All physics (no risk on validity)  limited physics (unknown & long CPU)
 in operation response  design loads
 limited test inputs  user chosen loads
 measurements only  all states known
 few designs  multiple (but 1 hour, 1 night,
several days, … thresholds)
 Cost : build and operate  Cost : setup, manipulate
In Out
Environment/Design point
System
Meta/reduced models
5
Full numerical
model
expensive
Meta-model
acceptable cost
Learning
points
Responses
Computation
points
LearningX LearningY
X
Estimations
Yˆ
Validity ?
• Regular relation
• Band-limited
• Spatial position
of inputs
• …
Predictive monitoring
of fuel circuit
Ph.D. of B. Lamoureux
~500 parameters
~100 indicators
~20 Inputs
Data from in
operation
measurements
System models of structural dynamics
Simple linear time invariant system
Extensions
• Coupling (structure, fluid,
control, multi-body, …)
• Optimization, variability,
damping, non linearity, …
When
Where
Sensors
Large/complex FEM
Historical keywords :
Modal analysis
Superelements, CMS, …
Ingredient 1 : variable separation
• General transient but
– limited bandwidth
– time invariant system
• Modal Analysis
response well approximated
in spatial sub-space
𝑞(𝑥, 𝑡) = 𝜙𝑗 𝑥 𝑎(𝑡)
𝑁𝑀
𝑗=1
• Space shapes =modes
• Time shapes =
generalized coordinates
7
SVD on the time response
• coincides with modes if
isolated resonances
• similar info for NL systems
8
Space / Time decomposition
Squeal limit cycle
PhD Vermot (Bosch)
NL system with impacts, PhD Thénint (EDF)
Data/model reduction
9
• SVD = data reduction through variable separation
– Extension to higher dimension variable separation see Chinesta (afternoon)
• Ritz analysis : build reduced dynamic models
– Reduced model = differential/analytic equation for qR(𝑡)/qR(𝑠)
– States qR allow restitution
– Assumptions on loading : band limited 𝑢 𝑡
restricted loads in space 𝑏𝑖 𝑥
F x,t = bi x u t
NA
i=1
– Learning = full FEM static & modes (McNeal, Craig-Bampton, …)
{q}N=
qR
Nx NR
T
𝑀𝑠2
+ 𝐶𝑠 + 𝐾 𝑞(𝑠) = 𝑏 𝑢(𝑠)
𝑇 𝑇 𝑍(𝑠) 𝑇 𝑁𝑅×𝑁𝑅 𝑞 𝑅(𝑠) = 𝑇 𝑇 𝑏 𝑢(𝑠)
Validity of reduced system models
Test & FEM system models assume
• Input restrictions
– Frequency band (modes)
– Localization (residual terms)
• System
– Time invariant
– Linear
Implemented in all major FEM & Modal Testing software
10
In Out
Environment/Design point
System
qR
Nx NR
T
System=IO relation
System=modal series
Challenge :
account for environment/design change
Sample design changes
• Material changes (visco damping)
• Junctions (contact)
• Component/system
Mesh/geometry
11
12
Ingredient 2 : parametric matrices
•Viscoelastic damping
𝐾𝑣 = 𝐾 𝐸(𝜔, 𝑇)
•Rotation induced stiffening
𝐾 𝐺 = 𝐾 Ω
•Contact stiffness evolution with
operating pressure
𝐾 𝑁 = 𝐾 p(x, 𝐹𝐺𝑙𝑜𝑏𝑎𝑙)
Reduction basis T can be fixed
for range of parameters
Speedup : 10-1e5
13
• Multi-model
• Other + residue iteration
• Example : strong coupling
With heavy fluids : modes of structure & fluid give
poor coupled prediction
Bases for parametric studies
Example water filled tank
With residualWithout residual
[T(p1) T(p2) … ]
Orthogonalization
[T]
[Tk] Rd
k=K-1 R(q(Tk))
Orthog [Tk Rd
k]
1th vertical mode: Main frame and
bow moving in phase
2
Co-simulation
SDToos/OSCAR
MSC/Motion
(VSD 2014)
Ingredient 3 : domain decomposition
• 1D  models coupled by few in/out :
hydraulic circuits, shaft torsion
• 3D FEM : classical uses
– Component Mode Synthesis/ Craig-Bampton
– Multibody with flexible superelements
• For each component base assumptions
remain
– LTI, few band-limited I/O
Two challenges
• Performance problems for large interfaces
• Component/system relation
14
Concept
Requirements &
architecture
Component design
System
operation
AVL Hydsim
Basic component coupling
Start : disjoint component models
Coupling relation between disjoint states
• Continuity 𝑞𝐼1 − 𝑞𝐼2 = 0
• Energy
+
16
Coupling + reduction
Classical CMS
• Reduced independently
• All interface motion (or interface modes)
• Assembly by continuity
Difficulties
• Mesh incompatibility
• Large interfaces
• Strong coupling (reduction requires knowledge of coupling)
Physical interface coupling
• Assembly by computation of interface energy
(example Arlequin)
Difficulties
• Use better bases than independent reduction
17
Squeal example : trace of system modes
CMS with trace of system modes
• No reduction of DOFs internal to contact area
• Reduction : trace of full brake modes on
reduced area & dependent DOFs (no need for
static response at interface)
Reduced model with exact system modes
Very sparse matrix for faster for time
integration
Component mode tuning method
• Reduced model is sparse
• Component mode amplitudes are DOFs
• Reduced model has exact nominal modes
(interest 1980 : large linear solution, 2015 : enhanced
coupling)
• Change component mode frequency  change the diagonal
terms of Kel
Disc
OuterPad
Inner Pad
Anchor
Caliper
Piston
Knuckle
Hub
wj
21
[M] [Kel] [KintS] [KintU]
CMT & design studies
• One reduced model /
multiple designs
Examples
• impact of modulus change
• damping real system or component mode
19
Component redesign
Sensitivity
energy analysis
Nom
.
+10
%
+20
%
-
20%
20
Conclusion
Reduced 3D models combine
• Variable separation
• Solve using generalized/reduced DOFs
–u(t) just assumed band-limited
–Restitution is possible
• Parametric matrices
• Domain decomposition
– Craig-Bampton is very costly
– Generalized coordinates can make sense
Challenges
• Engineering time to manage experiments
• Control data volume (>1e3 of NL runs)
• Control accuracy : develop software / train
engineers
In
Out
Environment
Design point
System
qR
Nx NR
T
www.sdtools.com/publications
Ritz

More Related Content

PDF
Nafems15 systeme
PPTX
Optimization of Electrical Machines in the Cloud with SyMSpace by LCM
PPTX
11 generalsisation of fluent
PDF
NONLINEAR MODEL PREDICTIVE CONTROL FOR OPERATION OF A POST COMBUSTION ABSORPT...
PPTX
Transformation of Random variables & noise concepts
PPTX
Multi phase mixture media
PPTX
Block Diagram For Control Systems.
PPTX
Optimal real-time landing using DNN
Nafems15 systeme
Optimization of Electrical Machines in the Cloud with SyMSpace by LCM
11 generalsisation of fluent
NONLINEAR MODEL PREDICTIVE CONTROL FOR OPERATION OF A POST COMBUSTION ABSORPT...
Transformation of Random variables & noise concepts
Multi phase mixture media
Block Diagram For Control Systems.
Optimal real-time landing using DNN

What's hot (20)

PPTX
Welch Verolog 2013
PDF
Parallel Algorithms
PDF
Lecture 6 ME 176 2 Time Response
PDF
E0812730
PDF
Real Time most famous algorithms
PDF
Keep Calm and React with Foresight: Strategies for Low-Latency and Energy-Eff...
PDF
Fpga implementation of optimal step size nlms algorithm and its performance a...
PDF
Start MPC
DOCX
Matrix Multiplication Report
PDF
Fpga implementation of optimal step size nlms algorithm and its performance a...
PDF
[Capella Day 2019] Model execution and system simulation in Capella
PDF
Power and Clock Gating Modelling in Coarse Grained Reconfigurable Systems
PPT
Chap2 slides
PDF
論文紹介 Probabilistic sfa for behavior analysis
PPTX
Test different neural networks models for forecasting of wind,solar and energ...
PPT
Chap8 slides
PDF
Ch.01 inroduction
PPT
Concepts of predictive control
PPT
UIC Panella Thesis
PPT
Lecture24 clockpower routing
Welch Verolog 2013
Parallel Algorithms
Lecture 6 ME 176 2 Time Response
E0812730
Real Time most famous algorithms
Keep Calm and React with Foresight: Strategies for Low-Latency and Energy-Eff...
Fpga implementation of optimal step size nlms algorithm and its performance a...
Start MPC
Matrix Multiplication Report
Fpga implementation of optimal step size nlms algorithm and its performance a...
[Capella Day 2019] Model execution and system simulation in Capella
Power and Clock Gating Modelling in Coarse Grained Reconfigurable Systems
Chap2 slides
論文紹介 Probabilistic sfa for behavior analysis
Test different neural networks models for forecasting of wind,solar and energ...
Chap8 slides
Ch.01 inroduction
Concepts of predictive control
UIC Panella Thesis
Lecture24 clockpower routing
Ad

Similar to Nafems15 Technical meeting on system modeling (20)

PDF
Areva10 Technical Day
PDF
Imac10 Component
PDF
ECL15_Damping_SDTools
PPTX
SDT experimental modal analysis capabilities
PDF
Efficiency enhancement techniques in finite element analysis navigating compl...
PPT
Introduction to Modeling and Simulations.ppt
PPTX
Ansys Stimulation Study
PDF
Flexible Multibody Dynamics
PDF
Finite Element Analysis Made Easy Lr
PPTX
Fem utkarsh
PDF
LANLSS__LRS_JMT_AJW_JDB
PDF
Towards Analysis Led Design
PPTX
cae ppt about its features and the applications
PPTX
Flow and heat transfer analysis through a brake disc: a CFD approach
PDF
Introduction LS-DYNA Session delivered by Mr.Suman M.L.J
PDF
Samtech Group
PDF
Samtech Group
PDF
CM2_Mode.pdf
PDF
Taguchi and the tablecloth trick | Insight, issue 5
PDF
MEMS Extraction & Verification
Areva10 Technical Day
Imac10 Component
ECL15_Damping_SDTools
SDT experimental modal analysis capabilities
Efficiency enhancement techniques in finite element analysis navigating compl...
Introduction to Modeling and Simulations.ppt
Ansys Stimulation Study
Flexible Multibody Dynamics
Finite Element Analysis Made Easy Lr
Fem utkarsh
LANLSS__LRS_JMT_AJW_JDB
Towards Analysis Led Design
cae ppt about its features and the applications
Flow and heat transfer analysis through a brake disc: a CFD approach
Introduction LS-DYNA Session delivered by Mr.Suman M.L.J
Samtech Group
Samtech Group
CM2_Mode.pdf
Taguchi and the tablecloth trick | Insight, issue 5
MEMS Extraction & Verification
Ad

Recently uploaded (20)

PDF
TFEC-4-2020-Design-Guide-for-Timber-Roof-Trusses.pdf
PPTX
M Tech Sem 1 Civil Engineering Environmental Sciences.pptx
PPTX
Engineering Ethics, Safety and Environment [Autosaved] (1).pptx
PPTX
IOT PPTs Week 10 Lecture Material.pptx of NPTEL Smart Cities contd
PPTX
CARTOGRAPHY AND GEOINFORMATION VISUALIZATION chapter1 NPTE (2).pptx
PPTX
CH1 Production IntroductoryConcepts.pptx
PDF
R24 SURVEYING LAB MANUAL for civil enggi
PDF
Operating System & Kernel Study Guide-1 - converted.pdf
PPTX
MCN 401 KTU-2019-PPE KITS-MODULE 2.pptx
PPTX
CYBER-CRIMES AND SECURITY A guide to understanding
PPTX
Construction Project Organization Group 2.pptx
PDF
Mitigating Risks through Effective Management for Enhancing Organizational Pe...
PDF
BMEC211 - INTRODUCTION TO MECHATRONICS-1.pdf
PPTX
Welding lecture in detail for understanding
DOCX
573137875-Attendance-Management-System-original
PDF
composite construction of structures.pdf
PDF
PRIZ Academy - 9 Windows Thinking Where to Invest Today to Win Tomorrow.pdf
PDF
Evaluating the Democratization of the Turkish Armed Forces from a Normative P...
PDF
SM_6th-Sem__Cse_Internet-of-Things.pdf IOT
PPTX
bas. eng. economics group 4 presentation 1.pptx
TFEC-4-2020-Design-Guide-for-Timber-Roof-Trusses.pdf
M Tech Sem 1 Civil Engineering Environmental Sciences.pptx
Engineering Ethics, Safety and Environment [Autosaved] (1).pptx
IOT PPTs Week 10 Lecture Material.pptx of NPTEL Smart Cities contd
CARTOGRAPHY AND GEOINFORMATION VISUALIZATION chapter1 NPTE (2).pptx
CH1 Production IntroductoryConcepts.pptx
R24 SURVEYING LAB MANUAL for civil enggi
Operating System & Kernel Study Guide-1 - converted.pdf
MCN 401 KTU-2019-PPE KITS-MODULE 2.pptx
CYBER-CRIMES AND SECURITY A guide to understanding
Construction Project Organization Group 2.pptx
Mitigating Risks through Effective Management for Enhancing Organizational Pe...
BMEC211 - INTRODUCTION TO MECHATRONICS-1.pdf
Welding lecture in detail for understanding
573137875-Attendance-Management-System-original
composite construction of structures.pdf
PRIZ Academy - 9 Windows Thinking Where to Invest Today to Win Tomorrow.pdf
Evaluating the Democratization of the Turkish Armed Forces from a Normative P...
SM_6th-Sem__Cse_Internet-of-Things.pdf IOT
bas. eng. economics group 4 presentation 1.pptx

Nafems15 Technical meeting on system modeling

  • 1. 1 Integration of reduced 3D models in vibration design processes. Examples from various industries. Etienne Balmes SDTools Arts et Métiers ParisTech NAFEMS Simulation des systèmes 3 Juin 2015
  • 2. • FEM simulations • System models (model reduction, state-space, active control, SHM) • Experimental modal analysis • Test/analysis correlation, model updating Activities 2 CAD/Meshing FEM Simulation Testing CATIA, Workbench, … NASTRAN, ABAQUS, ANSYS,... Adams, Simulink,... LMS TestLab, ME-Scope, … Simulation Validation SDT : MATLAB based toolbox Commercial since 1995 > 700 licenses sold Pantograph/catenary Modal test correlation Track dynamics
  • 3. Outline • Systems, models, dynamic models • Tools for model reduction – Variable separation – Parametric models – Domain decomposition • Conclusion 3
  • 4. A system = I/O representation Prototype Virtual prototype  All physics (no risk on validity)  limited physics (unknown & long CPU)  in operation response  design loads  limited test inputs  user chosen loads  measurements only  all states known  few designs  multiple (but 1 hour, 1 night, several days, … thresholds)  Cost : build and operate  Cost : setup, manipulate In Out Environment/Design point System
  • 5. Meta/reduced models 5 Full numerical model expensive Meta-model acceptable cost Learning points Responses Computation points LearningX LearningY X Estimations Yˆ Validity ? • Regular relation • Band-limited • Spatial position of inputs • … Predictive monitoring of fuel circuit Ph.D. of B. Lamoureux ~500 parameters ~100 indicators ~20 Inputs Data from in operation measurements
  • 6. System models of structural dynamics Simple linear time invariant system Extensions • Coupling (structure, fluid, control, multi-body, …) • Optimization, variability, damping, non linearity, … When Where Sensors Large/complex FEM Historical keywords : Modal analysis Superelements, CMS, …
  • 7. Ingredient 1 : variable separation • General transient but – limited bandwidth – time invariant system • Modal Analysis response well approximated in spatial sub-space 𝑞(𝑥, 𝑡) = 𝜙𝑗 𝑥 𝑎(𝑡) 𝑁𝑀 𝑗=1 • Space shapes =modes • Time shapes = generalized coordinates 7
  • 8. SVD on the time response • coincides with modes if isolated resonances • similar info for NL systems 8 Space / Time decomposition Squeal limit cycle PhD Vermot (Bosch) NL system with impacts, PhD Thénint (EDF)
  • 9. Data/model reduction 9 • SVD = data reduction through variable separation – Extension to higher dimension variable separation see Chinesta (afternoon) • Ritz analysis : build reduced dynamic models – Reduced model = differential/analytic equation for qR(𝑡)/qR(𝑠) – States qR allow restitution – Assumptions on loading : band limited 𝑢 𝑡 restricted loads in space 𝑏𝑖 𝑥 F x,t = bi x u t NA i=1 – Learning = full FEM static & modes (McNeal, Craig-Bampton, …) {q}N= qR Nx NR T 𝑀𝑠2 + 𝐶𝑠 + 𝐾 𝑞(𝑠) = 𝑏 𝑢(𝑠) 𝑇 𝑇 𝑍(𝑠) 𝑇 𝑁𝑅×𝑁𝑅 𝑞 𝑅(𝑠) = 𝑇 𝑇 𝑏 𝑢(𝑠)
  • 10. Validity of reduced system models Test & FEM system models assume • Input restrictions – Frequency band (modes) – Localization (residual terms) • System – Time invariant – Linear Implemented in all major FEM & Modal Testing software 10 In Out Environment/Design point System qR Nx NR T System=IO relation System=modal series Challenge : account for environment/design change
  • 11. Sample design changes • Material changes (visco damping) • Junctions (contact) • Component/system Mesh/geometry 11
  • 12. 12 Ingredient 2 : parametric matrices •Viscoelastic damping 𝐾𝑣 = 𝐾 𝐸(𝜔, 𝑇) •Rotation induced stiffening 𝐾 𝐺 = 𝐾 Ω •Contact stiffness evolution with operating pressure 𝐾 𝑁 = 𝐾 p(x, 𝐹𝐺𝑙𝑜𝑏𝑎𝑙) Reduction basis T can be fixed for range of parameters Speedup : 10-1e5
  • 13. 13 • Multi-model • Other + residue iteration • Example : strong coupling With heavy fluids : modes of structure & fluid give poor coupled prediction Bases for parametric studies Example water filled tank With residualWithout residual [T(p1) T(p2) … ] Orthogonalization [T] [Tk] Rd k=K-1 R(q(Tk)) Orthog [Tk Rd k]
  • 14. 1th vertical mode: Main frame and bow moving in phase 2 Co-simulation SDToos/OSCAR MSC/Motion (VSD 2014) Ingredient 3 : domain decomposition • 1D  models coupled by few in/out : hydraulic circuits, shaft torsion • 3D FEM : classical uses – Component Mode Synthesis/ Craig-Bampton – Multibody with flexible superelements • For each component base assumptions remain – LTI, few band-limited I/O Two challenges • Performance problems for large interfaces • Component/system relation 14 Concept Requirements & architecture Component design System operation AVL Hydsim
  • 15. Basic component coupling Start : disjoint component models Coupling relation between disjoint states • Continuity 𝑞𝐼1 − 𝑞𝐼2 = 0 • Energy +
  • 16. 16 Coupling + reduction Classical CMS • Reduced independently • All interface motion (or interface modes) • Assembly by continuity Difficulties • Mesh incompatibility • Large interfaces • Strong coupling (reduction requires knowledge of coupling) Physical interface coupling • Assembly by computation of interface energy (example Arlequin) Difficulties • Use better bases than independent reduction
  • 17. 17 Squeal example : trace of system modes CMS with trace of system modes • No reduction of DOFs internal to contact area • Reduction : trace of full brake modes on reduced area & dependent DOFs (no need for static response at interface) Reduced model with exact system modes Very sparse matrix for faster for time integration
  • 18. Component mode tuning method • Reduced model is sparse • Component mode amplitudes are DOFs • Reduced model has exact nominal modes (interest 1980 : large linear solution, 2015 : enhanced coupling) • Change component mode frequency  change the diagonal terms of Kel Disc OuterPad Inner Pad Anchor Caliper Piston Knuckle Hub wj 21 [M] [Kel] [KintS] [KintU]
  • 19. CMT & design studies • One reduced model / multiple designs Examples • impact of modulus change • damping real system or component mode 19 Component redesign Sensitivity energy analysis Nom . +10 % +20 % - 20%
  • 20. 20 Conclusion Reduced 3D models combine • Variable separation • Solve using generalized/reduced DOFs –u(t) just assumed band-limited –Restitution is possible • Parametric matrices • Domain decomposition – Craig-Bampton is very costly – Generalized coordinates can make sense Challenges • Engineering time to manage experiments • Control data volume (>1e3 of NL runs) • Control accuracy : develop software / train engineers In Out Environment Design point System qR Nx NR T www.sdtools.com/publications Ritz