Probabilistic Investigation of Sensitivities of Advanced Test- Analysis Model Correlation Methods

Similar documents
Multivariate problems and matrix algebra

4.5 JACOBI ITERATION FOR FINDING EIGENVALUES OF A REAL SYMMETRIC MATRIX. be a real symmetric matrix. ; (where we choose θ π for.

1 Online Learning and Regret Minimization

Thermal Diffusivity. Paul Hughes. Department of Physics and Astronomy The University of Manchester Manchester M13 9PL. Second Year Laboratory Report

Review of Calculus, cont d

Driving Cycle Construction of City Road for Hybrid Bus Based on Markov Process Deng Pan1, a, Fengchun Sun1,b*, Hongwen He1, c, Jiankun Peng1, d

13: Diffusion in 2 Energy Groups

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

8 Laplace s Method and Local Limit Theorems

System Identification with Noisy Data

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Purpose of the experiment

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

Chapter 9: Inferences based on Two samples: Confidence intervals and tests of hypotheses

u( t) + K 2 ( ) = 1 t > 0 Analyzing Damped Oscillations Problem (Meador, example 2-18, pp 44-48): Determine the equation of the following graph.

Student Activity 3: Single Factor ANOVA

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

Lecture 3. In this lecture, we will discuss algorithms for solving systems of linear equations.

Factors affecting the phonation threshold pressure and frequency

The Properties of Stars

Geometric Sequences. Geometric Sequence a sequence whose consecutive terms have a common ratio.

Integral equations, eigenvalue, function interpolation

Department of Electrical and Computer Engineering, Cornell University. ECE 4070: Physics of Semiconductors and Nanostructures.

A Signal-Level Fusion Model for Image-Based Change Detection in DARPA's Dynamic Database System

4.1. Probability Density Functions

Matrix Eigenvalues and Eigenvectors September 13, 2017

1 The fundamental theorems of calculus.

Tests for the Ratio of Two Poisson Rates

CHAPTER 4a. ROOTS OF EQUATIONS

Information processing via physical soft body. Supplementary Information

INTRODUCTION. The three general approaches to the solution of kinetics problems are:

Tutorial 4. b a. h(f) = a b a ln 1. b a dx = ln(b a) nats = log(b a) bits. = ln λ + 1 nats. = log e λ bits. = ln 1 2 ln λ + 1. nats. = ln 2e. bits.

Numerical Analysis: Trapezoidal and Simpson s Rule

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

Acceptance Sampling by Attributes

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

New data structures to reduce data size and search time

For the percentage of full time students at RCC the symbols would be:

Math 61CM - Solutions to homework 9

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS

Designing Information Devices and Systems I Spring 2018 Homework 7

99/105 Comparison of OrcaFlex with standard theoretical results

Chapter 4: Techniques of Circuit Analysis. Chapter 4: Techniques of Circuit Analysis

CS 188: Artificial Intelligence Spring 2007

Module 6: LINEAR TRANSFORMATIONS

SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD. Mehmet Pakdemirli and Gözde Sarı

Bridging the gap: GCSE AS Level

Reinforcement learning II

Simulation of Eclipsing Binary Star Systems. Abstract

Lecture 3 Gaussian Probability Distribution

Scientific notation is a way of expressing really big numbers or really small numbers.


Non-Linear & Logistic Regression

Industrial Electrical Engineering and Automation

Problem Set 3 Solutions

LECTURE NOTE #12 PROF. ALAN YUILLE

FEM ANALYSIS OF ROGOWSKI COILS COUPLED WITH BAR CONDUCTORS

Multi-Armed Bandits: Non-adaptive and Adaptive Sampling

A027 Uncertainties in Local Anisotropy Estimation from Multi-offset VSP Data

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3

Terminal Velocity and Raindrop Growth

Ordinary differential equations

CS 188 Introduction to Artificial Intelligence Fall 2018 Note 7

ELE B7 Power System Engineering. Unbalanced Fault Analysis

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

Emission of K -, L - and M - Auger Electrons from Cu Atoms. Abstract

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

12. DYNAMIC ANALYSIS. Force Equilibrium is Fundamental in the Dynamic Analysis of Structures 12.1 INTRODUCTION

Monte Carlo method in solving numerical integration and differential equation

Genetic Programming. Outline. Evolutionary Strategies. Evolutionary strategies Genetic programming Summary

Chapter 3 MATRIX. In this chapter: 3.1 MATRIX NOTATION AND TERMINOLOGY

DETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING MOMENT INTERACTION AT MICROSCALE

On the Uncertainty of Sensors Based on Magnetic Effects. E. Hristoforou, E. Kayafas, A. Ktena, DM Kepaptsoglou

Direct Design of Orthogonal Filter Banks and Wavelets by Sequential Convex Quadratic Programming

CALCULATED POWDER X-RAY DIFFRACTION LINE PROFILES VIA ABSORPTION

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

Introduction to Determinants. Remarks. Remarks. The determinant applies in the case of square matrices

LECTURE 14. Dr. Teresa D. Golden University of North Texas Department of Chemistry

Using QM for Windows. Using QM for Windows. Using QM for Windows LEARNING OBJECTIVES. Solving Flair Furniture s LP Problem

p-adic Egyptian Fractions

Lecture 14: Quadrature

Lecture 1. Functional series. Pointwise and uniform convergence.

Chapter 3 Polynomials

Quantum Physics II (8.05) Fall 2013 Assignment 2

Data Assimilation. Alan O Neill Data Assimilation Research Centre University of Reading

Multiscale Fourier Descriptor for Shape Classification

Numerical integration

Chapter 1: Fundamentals

We partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.

SOLVING SYSTEMS OF EQUATIONS, ITERATIVE METHODS

Derivations for maximum likelihood estimation of particle size distribution using in situ video imaging

New Expansion and Infinite Series

Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 25, March 29, 2006

The Islamic University of Gaza Faculty of Engineering Civil Engineering Department. Numerical Analysis ECIV Chapter 11

Math 426: Probability Final Exam Practice

Golden Section Search Method - Theory

3.1 Review of Sine, Cosine and Tangent for Right Angles

Sufficient condition on noise correlations for scalable quantum computing

INVESTIGATION OF MATHEMATICAL MODEL OF COMMUNICATION NETWORK WITH UNSTEADY FLOW OF REQUESTS

Transcription:

Probbilistic Investigtion of Sensitivities of Advnced Test- Anlysis Model Correltion Methods Liz Bergmn, Mtthew S. Allen, nd Dniel C. Kmmer Dept. of Engineering Physics University of Wisconsin-Mdison Rndll L. Myes Sndi Ntionl Lbortories

Test Anlysis Correltion 7,146 DOF FEM Updte FEM fils TAM/Test Correltion psses FEM vlidted Select sensor loctions nd reduce FEM 108 DOF TAM 108 Sensor Test fils TAM/FEM Correltion psses FEM mss mtrix must be reduced to test degrees of freedom (TAM) in order to compute modl orthogonlity. Orbitl Sciences 2

The Controversy Current FEM reduction lgorithms Sttic TAM: fils for hevy, soft structures. My be difficult to chieve good TAM/FEM correltion Fundmentl FEM propellnt mode (left) nd fundmentl FEM propellnt mode predicted by Sttic TAM (right) Improved Reduced Sttic (IRS) TAM: ill-conditioned under certin circumstnces Modl TAM: Trivil to chieve perfect TAM/FEM correltion, however it hs reputtion of being highly sensitive to experimentl or modl-mismtch errors 3

Purpose of Reserch Study the sensitivity of vrious TAMs to gin insight into fctors tht strongly ffect sensitivity A probbilistic nlysis will be used to chrcterize the effect of mesurement errors on TAM sensitivity 4

Relevnt Literture Guyn Freed, AM nd Flnign, CC (1990): Modl TAM most sensitive, sensors plced using modl kinetic energy Avitbile, P, Pechinsky, F, nd O Cllhn, J (1992): Sensor plcement is vitl to TAM performnce, SEREP nd Hybrid perform better thn Sttic TAM for smll sensor sets Chung, YT (1998): Sensor plcement ws not discussed nd no significnt difference could be seen between the TAMs Modl Lrge Reduced Cross Orthogonlity of Test Tower (Chung 1998) Avitbile, P, Pechinsky, F, nd O Cllhn, J (1992) 5

Relevnt Literture Gordis, JH (1992), Blelloch, P nd Vold, H (2005) : Notes ill-conditioning in dynmic reduction eqution: Proposes tht IRS TAM will be ill conditioned if the nturl frequencies of the structure with the o-set DOF pinned re similr to the frequencies of the structure of interest. Recently, this theory seems to hve been pplied to other TAM techniques such s the Modl TAM. 6

Model Generic Stellite 7,146 DOF Trget modes: first 18 consecutive flexible modes (0.3-11.8 Hz) 108 sensors Trget Mode 5 Trget Mode 6 Trget Mode 7 Trget Mode 8 2.7 Hz 2.8 Hz 3.5 Hz 3.7 Hz 7

Test Anlysis Models Sttic TAM = sensor loction o = omitted DOF Eigenvlue problem 2 M ωi M o M M o oo φ i + φio K K o K K o oo φ i = φio 0 2 2 Lower prtition eqution [ Ko ωi M o ]{ φi} + [ Koo ωi M oo ]{ φio} = 0 Neglect the mss of the o-set DOF 0 0 { } 2 1 φ = K ω M K ω 2 M { φ } io oo i oo o i o i Sttic Trnsformtion Mtrix (ech column represents constrint mode) I K K [ ] TS = 1 oo o 8

Test Anlysis Models IRS TAM { } 2 1 φ = K ω M K ω 2 M { φ } io oo i oo o i o i Ill-conditioned when ω i2 is ner ny of the eigenvlues of the Koo, Moo system Approximte the frequency terms Clculte the IRS trnsformtion mtrix ω 2 ~ ~ i i S S 1 { φ } = M K { φ } [ T ] = [ T ] + [ T ] IRS S i i M 0 0 ~ 1 1 1 0 K oo M o M oo Koo K o o [ Ti ] = M S K S M I ~ 9

Test Anlysis Models IRS TAM O-set system Mode 1 16.8 Hz FEM Trget Mode 18 11.8 Hz 10

Test Anlysis Models Sttic nd IRS TAM Mss weighted effective independence did not select the lumped msses (the lumped msses were essentil to TAM- FEM correltion) Modl kinetic energy pplied to ll 18 trget modes ws not sufficient A significnt mount of hnd selection nd engineering judgment ws used (modified modl kinetic energy method) 5 core lumped msses 11

Test Anlysis Models Modl TAM Physicl coordintes in terms of modl coordintes Prtitioned Equtions Solve for modl coordintes in terms of the sensor DOF Modl trnsformtion mtrix x x o φ = φo {} q { x } = [ φ ]{ q} { x } = [ φ ]{ q} o {} [ ] [ T 1 T q = φ φ φ ]{ x } [ T ] M o I T φo ( φ φ ) = 1 T φ 12

Test Anlysis Models Modl TAM Sensor plcement chieved with Effective Independence Mximize the determinnt of the Fisher informtion mtrix mx Q = mx φ T φ Effective Independence E Di 1 = φ Q φ 0.0 E 1. 0 i it Di 13

Test Anlysis Models Modl TAM Modl TAM o-set frequencies re similr to the FEM frequencies, so the theory of Gordis suggests tht this TAM will be sensitive. O-set system Mode 2 FEM Trget Mode 4 O-set system Mode 5 FEM Trget Mode 7 1.2 Hz 1.8 Hz 3.2 Hz 3.5 Hz 14

Test Anlysis Models Modl using Condition Number Sensor Plcment Modl coordintes in terms of the sensor DOF {} [ ] [ T T q = φ φ φ ]{ x } 1 Solution is more sensitive if the condition number of, is lrge. φ Begin with visuliztion set, nd dd sensors tht minimize the condition number of φ 15

Test Anlysis Correltion 7,146 DOF FEM Updte FEM fils TAM/Test Correltion psses FEM vlidted Select sensor loctions nd reduce FEM 108 DOF TAM 108 Sensor Test fils TAM/FEM Correltion psses Orbitl Sciences 16

Orthogonlity Criteri: O = φ Cross Orthogonlity Criteri: CO = φfem Frequency Comprison Criteri: Correltion Metrics 0 off digonl term 0.1 [ ] [ T M ][ φ ] FEM ~ TAM 0 off digonl term 0.1 0.95 digonl term 1.0 f error ~ FEM [ ] [ T M ][ φ ] = f FEM f FEM f TAM TAM TAM * 100 3% 17

TAM-FEM Correltion Mx off digonl term: 0.05 Mx off digonl term: 6e-4 *Modl TAM lwys produces perfect orthogonlity for TAM-FEM correltion 18

Test Anlysis Correltion 7,146 DOF FEM Updte FEM fils TAM/Test Correltion psses FEM vlidted Select sensor loctions nd reduce FEM 108 DOF TAM 108 Sensor Test fils TAM/FEM Correltion psses Orbitl Sciences 19

{} i φ FEM { U}= Noise Model nd Simulted Test Mode Shpes i φ1 i φ2 = i φ n i φn 1 Mx vlue *2%* {} i φ {} { } i Test = φ noise + φ FEM { U} = { φ} noise column vector of uniformly distributed rndom numbers between -1 nd 1 FEM trget mode 15 Test trget mode 15 20

Noise Model FEM ssumed to be perfect Noise contminted mode shpe Noise vector models the net effect of ll errors tht cuse the FEM mode shpes to disgree with the test mode shpes. Noise Distribution: Uniform no ssumption is mde bout the distribution of noise Noise Amplitude: Sensors with the smllest motion hve the lrgest noise to signl rtio Noise is smll on verge: ± 2% t sensor loctions with the lrgest motion. 21

TAM-Test Test Correltion Results (1 cse of Rndom Noise) Mx off digonl term: 0.27 Mx off digonl term: 0.79 22

TAM-Test Test Correltion Results (1 cse of Rndom Noise) Mx off digonl term: 0.05 Mx off digonl term: 0.04 23

Monte Crlo Simultion Thus fr, TAM-Test correltion hs been studied using only one noise profile Rndom noise dded in 10,000 itertions Orthogonlity computed for ech itertion Mximum off-digonl term of orthogonlity ws stored 24

TAM-Test Test Correltion Results Pss Fil Despite its low o-set frequencies, Modl TAM does not show high sensitivity! 25

TAM-Test Test Correltion Results If Orthogonlity > 0.1 one might Refine FEM before exiting test Repet test nd/or look for errors Updte the FEM In this cse, the FEM ws perfect (errors in test modes were purely rndom) Note: The specific rnking of different TAM methods my depend on: The structure of interest The chrcteristics of the noise Systemtic errors between the test nd FEM 26

TAM-Test Test Correltion Results Pss Fil Sensor selection is criticl to the performnce of ech TAM Most previous studies used the sme sensor set, usully optimized for the Sttic TAM 27

Predicting Stndrd Devition Recently, we hve developed formuls to nlyticlly predict sensitivity of TAM bsed on simple metrics For exmple, for the noise model used in this study: T Oij = [ φi + ni ] M TAM [ φi + ni ] ni = noise ( 2 T ) ( 2 ) ( ) 2 2 2 σ ( O ) = M φ σ + φ M σ + M σσ ij TAM j i i TAM j TAM mn i j m m m m m n Sttic Modl EFI Modl C#TAM Modl with Sttic DOF Mximum Orthogonlity Off-Digonl Predicted STD 0.03 0.009 0.006 0.02 Actul STD 0.03 0.01 0.006 0.02 28

Conclusions nd Future Work Conclusions IRS TAM ws ill-conditioned, s predicted by Gordis Modl TAM did not show high sensitivity even though its o-set frequencies were ner those of the trget modes Probbilistic nlysis more fully explins TAM sensitivity One cn even predict the sensitivity of the TAMs nlyticlly given the TAM Mss mtrix, mode shpes nd noise model. Future Work Develop more ccurte noise models Study the effect of systemtic mismtch between FEM nd test due to modeling errors. My need the Hybrid TAM in these cses Apply these methods to other physicl systems, nlyticlly nd experimentlly. Investigte systems with non-consecutive trget modes 29