Toward Predictive Digital Twins
via component-based reduced-order models and interpretable machine learning
Michael Kapteyn*, Dr. David Knezevic, Prof. Karen Willcox
AIAA SciTech | Paper AIAA-2020-0418 | January 6, 2020
Toward Predictive Digital Twins via component-based reduced-order - - PowerPoint PPT Presentation
Toward Predictive Digital Twins via component-based reduced-order models and interpretable machine learning Michael Kapteyn*, Dr. David Knezevic, Prof. Karen Willcox AIAA SciTech | Paper AIAA-2020-0418 | January 6, 2020 Outline 1 Motivation
AIAA SciTech | Paper AIAA-2020-0418 | January 6, 2020
sense structural data estimate structural state predict flight capabilities dynamically replan mission
1
sense structural data estimate structural state predict flight capabilities dynamically replan mission
2
Customized 12ft Telemaster aircraft:
Complex structure with multiple materials Custom wing sets: pristine & damaged Custom sensor suite
3
3 axis accelerometer 3 axis gyro Dual high-frequency dynamic strain and vibration sensors Temperature, pressure and humidity sensors
*One of the authors has a family member who is co-founder of Divinio. Purchase of the sensors for use in the research was reviewed and approved in compliance with all applicable MIT policies and procedures.
4
Use model library to train a classifier that predicts asset state based on sensor data Construct library of reduced-order models representing different asset states
sensor data
Analysis, Prediction, Optimization
updated digital twin current digital twin 5
component interior component port 𝑑 governing PDE (in our case linear elasticity) computational mesh damage parameters
reduced stiffness material loss crack length delamination size … Young’s modulus Poisson’s ratio number of plies ply angles …
geometric parameters non-geometric parameters
6
system parameters 𝜈 = [𝜈%, 𝜈', 𝜈(] component parameters 𝜈% Instantiate and Assemble Apply Loads + assembly parameters 𝜈' + load parameters 𝜈( =
7
Start with the usual finite element problem statement: Find 𝑣, ∈ 𝑊
, such that 𝑏 𝑣,, 𝑤 ; 𝜈 = 𝑔 𝑤; 𝜈 ∀ 𝑤 ∈ 𝑊 ,
𝐵5,5 𝐵5,67 𝐵5,68 𝐵5,67
9
𝐵67,67 𝐵5,68
9
𝐵68,68 𝕍 𝑣67 𝑣68 = 𝑔
5
𝑔
67
𝑔
68
Express interior DOFs in terms of port DOFs 𝐵6<,6<𝑣67 = 𝑔
67 − 𝐵5,67 9
𝕍 Substitute to get a system involving only port DOFs: 𝕋 𝜈 𝕍 𝜈 = 𝔾(𝜈) Issue: Schur complement 𝕋(𝜈) is large (𝐍×𝐍), and expensive to compute
M port DOFs N interior DOFs ΩG ΩH 𝑄
Solve on each component independently
port DOFs Interior DOFs
8
[Huynh 2013]
i. Port Reduction: Retain only the first 𝑛 dominant modes at component ports ▸ Reduces the size of 𝕋: M × M 𝑛 × 𝑛 ii. Component Interior Reduction: Replace the finite element space inside each component with a reduced basis (RB) space of dimension 𝑜 ▸ Reduces the size of matrices required to compute entries of 𝕋: N × N 𝑜 × 𝑜
M port DOFs N interior DOFs ΩG ΩH 𝑄
9
▸ Never have to solve full-system FE model
▸ System may have many spatially distributed parameters
▸ Allows for expressive adaptation: changes to topology, meshes etc.
10
root top skin top skin bottom skin spar caps shear web ribs flaps aileron linkages circular rods
11
increasing effective damage (reduction in stiffness)
0% 80% 20% 40% 60%
damage region
12
asset state noisy sensor data Forward (predictive) model Inverse (reactive) model
13
From component-based model to digital twin: Interpretable machine learning
0% 80% 20% 40% 60%
60 70 80 90 100 80 90 100 110 120
Sensor 16 Sensor 24
80% 60% 40% 20% 0%
Component 1 Component 2
300 350 400 450 500 550 600 650 80 90 100 110
Sensor 22 Sensor 8
80% 60% 40% 20% 0%
sensor 22 < 429? 80% 𝑜 𝑧 20% 0% 𝑜 𝑧 sensor 22 < 495? sensor 22 < 383? sensor 22 < 351? 𝑜 𝑧 60% 40% 𝑜 𝑧 sensor 24 < 103? 80% 𝑜 𝑧 20% 0% 𝑜 𝑧 sensor 16 + 1.6365*(sensor 24) < 237? sensor 16 +0.3675*(sensor 24) < 112? sensor 16 < 73? 𝑜 𝑧 60% 40% 𝑜 𝑧
14
sensor selection
standard neural networks
15
From component-based model to digital twin: Interpretable machine learning
Use model library to train a classifier that predicts asset state based on sensor data Construct library of reduced-order models representing different asset states
sensor data
Analysis, Prediction, Optimization
updated digital twin current digital twin 16
Test with experimental data Incorporate multimodal observations Flight demonstration
18
For a project overview, slides, and the full paper, visit https://kiwi.oden.utexas.edu/research/digital-twin
Funding acknowledgements:
19
[Kordonowy 2011] [Singh 2017] [Vallaghé 2015] [Huynh 2013] [Bertsimas 2019] Kordonowy, D., and Toupet, O., "Composite airframe condition-aware maneuverability and survivability for unmanned aerial vehicles." Infotech@ Aerospace 2011, 2011-1496. Singh, V., and Willcox, K.E., "Methodology for Path Planning with Dynamic Data-Driven Flight Capability Estimation." AIAA Journal (2017): 2727-2738. Vallaghé, S., et al. "Component-based reduced basis for parametrized symmetric eigenproblems."Advanced Modeling and Simulation in Engineering Sciences 2.1 (2015): 7. Huynh, D.B.P., D.J. Knezevic, and A.T. Patera. "A static condensation reduced basis element method: approximation and a posteriori error estimation." ESAIM: Mathematical Modelling and Numerical Analysis 47.1 (2013): 213-251. Bertsimas, D., and Dunn, J., "Machine Learning under a Modern Optimization Lens." Dynamic Ideas (2018).