content variations on cusum tests for flutter monitoring
play

Content Variations on CUSUM tests for flutter monitoring - PowerPoint PPT Presentation

Content Variations on CUSUM tests for flutter monitoring Introduction Mich` ele Basseville, Laurent Mevel, Rafik Zouari Subspace-based residual for modal monitoring IRISA (CNRS & INRIA), Rennes, France CUSUM test for monitoring a scalar


  1. Content Variations on CUSUM tests for flutter monitoring Introduction Mich` ele Basseville, Laurent Mevel, Rafik Zouari Subspace-based residual for modal monitoring IRISA (CNRS & INRIA), Rennes, France CUSUM test for monitoring a scalar instability index Variations on the CUSUM test Eurˆ eka project no 3341 FliTE2 Experimental results michele.basseville@irisa.fr -- http://www.irisa.fr/sisthem/ Conclusion 1 2 Introduction - (1) Introduction - (2) • Flutter: critical aircraft instability phenomenon unfavorable interaction of aerodynamic, elastic and inertial • Aim of in-flight online flutter monitoring: forces; may cause major failures Early detection of a deviation in the aircraft modal parameters • Flight flutter testing, very expensive and time consuming : before it develops into flutter. Design the flutter free flight envelope • Flutter clearance techniques: • Change-point detection: natural approach In-flight identification: output-only, or using input excitations Data processing: time-frequency, wavelet, envelope function • For a scalar instability criterion ψ and a critical value ψ c , Flutter prediction based on model-based approaches: online hypotheses testing: flutterometer ( µ -robustness), physical model updating H 0 : ψ > ψ c and H 1 : ψ ≤ ψ c • Some challenges: Real time on-board monitoring, • CUSUM test as an approximation to the optimal test Handling transients between steady flight test points • Our approach: • Variations on the CUSUM test Statistical detection for monitoring instability indicators 3 4

  2.   Λ modes θ ∆   Subspace-based residual for modal monitoring Canonical parameter : =     vec Φ mode shapes   Φ   X k +1 = F X k + V k F φ λ = λ φ λ       Φ∆       O p +1 ( θ ) =   Observability in modal basis :  .   .   .        ∆  Φ∆ p   Y k = H X k ϕ λ = H φ λ     Given: R 0 R 1 R 2 . . .       � � R 1 R 2 R 3 . . .   H ⋆ ∆ ∆ • a reference parameter θ ⋆ , by SVD of � Y k Y T p +1 ,q (reference data)   H R i = E , =   k − i   R 2 R 3 R 4 . . .       . . ... .  . . .  . . . U ( θ ⋆ ) T H ⋆ � p +1 ,q = 0 parameter estimating function R i = H F i G = ⇒ H = O C U ( θ ⋆ ) T O p +1 ( θ ⋆ ) = 0 , U ( θ ⋆ ) T U ( θ ⋆ ) = I   � G H � C ∆ F 2 G   = F G . . .   HF   O ∆   = ,   • a n -size sample of new data; H p +1 ,q � HF 2   � �   G ∆ X k Y T  .  = E . . k For testing θ = θ ⋆ , statistics (residual) : Output-only covariance-driven subspace identification = √ n vec � � ζ n ( θ ⋆ ) ∆ U ( θ ⋆ ) T H p +1 ,q � SVD of H − → O − → ( H, F ) − → ( λ, ϕ λ ) 5 6 Local approach to testing H 1 : θ = θ ⋆ + Υ / √ n Data-driven computation for online detection � � H 0 : θ = θ ⋆ and n − p k = q Z k ( θ ⋆ ) / √ n Mean sensitivity and covariance matrices: � ζ n ( θ ⋆ ) = = 1 / √ n ∂/∂ ˜ � � � J n ( θ ⋆ , θ ) ∆ � θ = θ ⋆ , Σ n ( θ ⋆ , θ ) ∆ ζ n ( θ ⋆ ) ζ n ( θ ⋆ ) T θ E θ ζ n (˜ θ ) � = E θ � � ˜ � � T Y + Z k ( θ ⋆ ) ∆ k,p +1 Y − T = K k ( θ ⋆ , θ ) vec U ( θ ⋆ ) k,q If Σ n ( θ ⋆ , θ ) is positive definite, and for all Υ , under both hypoth: Σ n ( θ ⋆ , θ ) − 1 / 2 ( ζ n ( θ ⋆ ) − J n ( θ ⋆ , θ ) Υ) Another approximation → N (0 , I ) n → ∞ For n large enough, and k = 1 , . . . , n , Normalized residual: Z k ( θ ⋆ ) ≈ Gaussian i.i.d., mean 0 before change and � = 0 after. ζ n ( θ ⋆ ) ∆ = K n ( θ ⋆ , θ ) ζ n ( θ ⋆ ) K n ( θ ⋆ , θ ) ∆ , Σ n ( θ ⋆ , θ ) ∆ = Σ − 1 / 2 J T n Σ − 1 = J T n Σ − 1 n J n Monitoring any function ψ ( θ ) n n Replace J n ( θ ⋆ , θ ) with J n ( θ ⋆ , θ ) J ⋆ θψ , where J ⋆ θψ = ∂θ/∂ψ | θ = θ ⋆ . � � ζ n ( θ ⋆ ) − Σ n ( θ ⋆ , θ ) 1 / 2 Υ → N (0 , I ) n → ∞ 7 8

  3. CUSUM test for monitoring a scalar index The crossing of a critical ψ c by ψ is reflected into a change with the same sign in the mean ν of the i.i.d. Gaussian Z k ( θ ⋆ ) . Variations on the CUSUM test - (1) The CUSUM test may be used for testing between: H 1 : ν ≤ 0 H 0 : ν > 0 and For detecting aircraft instability precursors, select: Procedure for unknown sign and magnitude of change in ψ : An instability criterion ψ and a critical value ψ c ; a) i) Set a min. change magnitude ν m > 0 , and test between: H 0 : ν > ν m / 2 and H 1 : ν ≤ − ν m / 2 b) A left kernel matrix U ( . ) ; n − p S n ( θ ⋆ ) ∆ k = q ( Z k ( θ ⋆ ) + ν m ) , T n ( θ ⋆ ) ∆ � = = k = q,...,n − p S k ( θ ⋆ ) max A reference θ ⋆ for estimating J n ( θ ⋆ ) and Σ − 1 c) n ( θ ⋆ ) ; H 1 g n ( θ ⋆ ) ∆ > = T n ( θ ⋆ ) − S n ( θ ⋆ ) < ̺ threshold d) A min. change magnitude ν m and a threshold ̺ . H 0 ii) Run 2 tests in parallel, for decreasing and increasing ψ ; iii) Make a decision from the first test which fires; iv) Reset all sums and extrema to 0, switch to the other test. 9 10 Solution 2. - Details Variations on the CUSUM test - (2) i) Compute the critical eigenvalues λ c at flight point t using identified modal signatures ( θ 1 , . . . , θ t ) Three solutions for b)-c): and extrapolation of the characteristic polynomial associated with the quasi-steady aeroelastic model 1. θ ⋆ ∆ q + ( K + V 2 C ) q = 0 ; = θ 0 identified on reference data for the stable system; M ¨ q + ( D + V B ) ˙ U ( θ ⋆ ) computed, J n ( θ 0 ) , Σ − 1 ii) Build the critical modal signature θ c from λ c and n ( θ 0 ) estimated recursively with the test data. the mode-shapes ϕ λ identified at flight point t ; 2. θ ⋆ ∆ = θ c , critical parameter closer to instability, computed iii) Compute the Z k ( θ c ) ’s and S n ( θ c ) ; at each flight point using θ 0 and an aeroelastic model; Σ − 1 J n ( θ c ) and ˆ ˆ Compute n ( θ c ) with the test data; U ( θ ⋆ ) computed, J n ( θ c ) , Σ − 1 n ( θ c ) estimated recursively with the test data. iv) Run the CUSUM test between flight points t and t + 1 ; v) Repeat these steps for flight point t + 1 : 3. U ( . ) ∆ � = U n estimated on test data, modal identification of θ t +1 to update the prediction of θ c , J n ( θ 0 ) , Σ − 1 n ( θ 0 ) estimated recursively with the test data. J n ( θ c ) and ˆ ˆ Σ − 1 computation of n ( θ c ) , CUSUM test between t + 1 and t + 2 . 11 12

  4. Solution 3. - Details Solution 3. - Details (Contd.) i) Initialization For an initial airspeed: ˆ J n ; Estimate a reference θ 0 and compute the constant terms in Select data sample size L , lag τ , block size K , ν m , and ̺ ; Σ − 1 Compute ˆ ˆ J L + τ with the first L + τ samples; L + τ and Compute ˆ U L + τ with ( Y 1 , ..., Y L ) ; Compute the Z k ’s and S L + τ . ii) Recursive loop Running the CUSUM test: for each n ≥ L + τ : � Compute recursively U n with ( Y n − τ − L +1 , ..., Y n − τ ) ; Σ − 1 � � J n to compute S n and g n until g n ≥ ̺ . � Use U n with and n Σ − 1 � J n every K samples. � Update recursively and n 13 14 Example - Aeroelastic Hancock wing model Example - Numerical results Rigid wing with constant chord; 2 d.o.f. in bending and torsion. CUSUM test run with ν m = 0 . 1 , ̺ = 100 , and the damping as ψ . J n , ˆ ˆ Solution 1. with θ ⋆ = θ 0 at V = 20 m/s , online recursive Σ n . J n , ˆ ˆ Matrix F , and eigenvalues λ : functions of airspeed V . Solution 2. with θ ⋆ = θ c at V = 85 m/s , online recursive Σ n . Solution 3. with online recursive ˆ U n , ˆ J n , ˆ Flutter airspeed: V f = 88 . 5 m/s . Σ n . Stability indicator ψ : Damping coefficient Alarm onset times depend on threshold; ˆ V f is more important. Modal frequencies variation with airspeed Modal damping coefficients variation with airspeed 0.2 13 Solution 1. θ ⋆ far from instability, alarm at V =67 m/s , ˆ V f =65 m/s . Bending mode Bending mode Torsional mode Torsional mode 12 0.15 The test detects that torsional damping decreases under ψ c . Damping coefficient 11 0.1 F(Hz) 10 0.05 Solution 2. θ ⋆ close to instability, alarm at V =88 m/s , ˆ 9 V f =85 m/s . 0 8 The test detects that flutter is happening between two steady 7 −0.05 0 10 20 30 40 50 60 70 80 90 0 10 20 30 40 50 60 70 80 88.5 V(m/s) V(m/s) points, and confirms the flutter prediction. Frequencies Damping coefficients Solution 3. Alarm at V =88 m/s , ˆ V f =78 m/s much closer to flutter. Bending & torsion modes Good behavior for light damping decrease before alarm. 20700 -size 2D-samples simulated (300 for each V =20:1:88 m/s ). Detection (before flutter) of torsional damping drop. 15 16

Download Presentation
Download Policy: The content available on the website is offered to you 'AS IS' for your personal information and use only. It cannot be commercialized, licensed, or distributed on other websites without prior consent from the author. To download a presentation, simply click this link. If you encounter any difficulties during the download process, it's possible that the publisher has removed the file from their server.

Recommend


More recommend