0
Research Papers

Design Optimization Toward Alleviating Forced Response Variation in Cyclically Periodic Structure Using Gaussian Process

[+] Author and Article Information
K. Zhou, A. Hegde, P. Cao

Department of Mechanical Engineering,
University of Connecticut,
191 Auditorium Road, Unit 3139,
Storrs, CT 06269

J. Tang

Professor
Department of Mechanical Engineering,
University of Connecticut,
191 Auditorium Road, Unit 3139,
Storrs, CT 06269
e-mail: jtang@engr.uconn.edu

1Corresponding author.

Contributed by the Technical Committee on Vibration and Sound of ASME for publication in the JOURNAL OF VIBRATION AND ACOUSTICS. Manuscript received January 20, 2016; final manuscript received October 3, 2016; published online December 7, 2016. Assoc. Editor: Mahmoud Hussein.

J. Vib. Acoust 139(1), 011017 (Dec 07, 2016) (14 pages) Paper No: VIB-16-1038; doi: 10.1115/1.4035107 History: Received January 20, 2016; Revised October 03, 2016

Cyclically periodic structures, such as bladed disk assemblies in turbomachinery, are widely used in engineering systems. It is well known that small uncertainties exist among their substructures, which in certain situations may cause drastic change in the dynamic responses, a phenomenon known as vibration localization. Previous studies have suggested that the introduction of small, prespecified design modification, i.e., intentional mistuning, may alleviate vibration localization and reduce response variation. However, there has been no systematic methodology to facilitate the optimal design of intentional mistuning. The most significant challenge is the computational cost involved. The finite-element model of a bladed disk usually requires a very large number of degrees-of-freedom (DOFs). When uncertainties occur in a cyclically periodic structure, the response may no longer be considered as simple perturbation to that of the nominal structure. In this research, a suite of interrelated algorithms is proposed to enable the efficient design optimization of cyclically periodic structures toward alleviating their forced response variation. We first integrate model order reduction with a perturbation scheme to reduce the scale of analysis of a single run. Then, as the core of the new methodology, we incorporate Gaussian process (GP) emulation to conduct the rapid sampling-based evaluation of the design objective, which is a metric of response variation under uncertainties, in the parametric space. The optimal design modification can thus be directly identified to minimize the response variation. The efficiency and effectiveness of the proposed methodology are demonstrated by systematic case studies.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.

References

Pierre, C. , 1990, “ Weak and Strong Vibration Localization in Disordered Structures: A Statistical Investigation,” J. Sound Vib., 139(1), pp. 111–132. [CrossRef]
Slater, J. C. , Minkiewicz, G. R. , and Blair, A. J. , 1999, “ Forced Response of Bladed Disk Assemblies—A Survey,” Shock Vib. Dig., 31(1), pp. 17–24. [CrossRef]
Lim, S. H. , Pierre, C. , and Castanier, M. P. , 2006, “ Predicting Blade Stress Levels Directly From Reduced-Order Vibration Models of Mistuned Bladed Disks,” ASME J. Turbomach., 128(1), pp. 206–210. [CrossRef]
Tang, J. , and Wang, K. W. , 2003, “ Vibration Delocalization of Nearly Periodic Structures Using Coupled Piezoelectric Networks,” ASME J. Vib. Acoust., 125(1), pp. 95–108. [CrossRef]
Yu, H. , Wang, K. W. , and Zhang, J. , 2006, “ Piezoelectric Networking With Enhanced Electromechanical Coupling for Vibration Delocalization of Mistuned Periodic Structures-Theory and Experiment,” J. Sound Vib., 295(1–2), pp. 246–265. [CrossRef]
Yu, H. , and Wang, K. W. , 2007, “ Piezoelectric Networks for Vibration Suppression of Mistuned Bladed Disks,” ASME J. Vib. Acoust., 129(5), pp. 559–566. [CrossRef]
Castanier, M. P. , and Pierre, C. , 1997 , “ Consideration on the Benefits of Intentional Blade Mistuning for the Forced Response of Turbomachinery Rotors,” American Society of Mechanical Engineers, Aerospace Division (Publication) AD, Vol. 55, pp. 419–425.
Castanier, M. P. , and Pierre, C. , 2002, “ Using Intentional Mistuning in the Design of Turbomachinery Rotors,” AIAA J., 40(10), pp. 2077–2086. [CrossRef]
Choi, B. K. , Lentz, J. , Rivas-Guerra, A. J. , and Mignolet, M. P. , 2003, “ Optimization of Intentional Mistuning Patterns for the Reduction of the Forced Response Effects of Unintentional Mistuning: Formulation and Assessment,” ASME J. Eng. Gas Turbines Power, 125(1), pp. 131–140. [CrossRef]
Yu, C. B. , Wang, J. J. , and Li, Q. H. , 2011, “ Investigation of the Combined Effects of Intentional Mistuning, Damping and Coupling on the Forced Response of Bladed Disks,” J. Vib. Control, 17(8), pp. 1149–1157. [CrossRef]
Craig, R. R. , and Kurdila, A. J. , 2006, Fundamentals of Structural Dynamics, 2nd ed., Wiley, Hoboken, NJ.
Castanier, M. P. , and Pierre, C. , 2006, “ Modeling and Analysis of Mistuned Bladed Disk Vibration: Status and Emerging Directions,” J. Propul. Power, 22(2), pp. 384–395. [CrossRef]
Masson, G. , AitBrik, B. , Cogan, S. , and Bouhaddi, N. , 2006, “ Component Mode Synthesis (CMS) Based on an Enriched Ritz Approach for Efficient Structural Optimization,” J. Sound Vib., 296(4–5), pp. 845–860. [CrossRef]
Castanier, M. P. , Ottarsson, G. , and Pierre, C. , 1997, “ A Reduced Order Modeling Technique for Mistunied Bladed Disks,” ASME J. Vib. Acoust., 119(3), pp. 439–446. [CrossRef]
Rivas-Guerra, A. J. , and Mignolet, M. P. , 2003, “ Maximum Amplification of Blade Response Due to Mistuning: Localization and Mode Shape Aspects of the Worst Disks,” ASME J. Turbomach., 125(3), pp. 442–454. [CrossRef]
Marinescu, O. , Epureanu, B. I. , and Banu, M. , 2011, “ Reduced Order Models of Mistuned Cracked Bladed Disks,” ASME J. Vib. Acoust., 133(5), p. 051014. [CrossRef]
D'Souza, K. , Saito, A. , and Epureanu, B. I. , 2012, “ Reduced-Order Modeling for Nonlinear Analysis of Cracked Mistuned Multistage Bladed-Disk Systems,” AIAA J., 50(2), pp. 304–312. [CrossRef]
Williams, C. K. I. , and Barber, D. , 1998, “ Bayesian Classification With Gaussian Processes,” IEEE Trans. Pattern Anal. Mach. Intell., 20(12), pp. 1342–1351. [CrossRef]
Lebarbier, E. , 2005, “ Detecting Multiple Change-Points in the Mean of Gaussian Process by Model Selection,” Signal Process., 85(4), pp. 717–736. [CrossRef]
Rasmussen, C. E. , and Williams, C. K. I. , 2006, Gaussian Process for Machine Learning, MIT Press, Cambridge, MA.
Kennedy, M. C. , and O'Hagan, A. , 2000, “ Predicting the Output From a Complex Computer Code When Fast Approximation Are Available,” Biometrika, 87(1), pp. 1–13. [CrossRef]
O'Hagan, A. , 2006, “ Bayesian Analysis of Computer Code Outputs: A Tutorial,” Reliab. Eng. Syst. Saf., 91(10–11), pp. 1290–1300. [CrossRef]
DiazDelaO, F. A. , and Adhikari, S. , 2010, “ Structural Dynamic Analysis Using Gaussian Process Emulators,” Eng. Comput., 27(5), pp. 580–605. [CrossRef]
DiazDelaO, F. A. , and Adhikari, S. , 2011, “ Gaussian Process Emulators for the Stochastic Finite Element Method,” Int. J. Numer. Methods Eng., 87(6), pp. 521–540. [CrossRef]
Xia, Z. , and Tang, J. , 2013, “ Characterization of Dynamic Response of Structures With Uncertainty by Using Gaussian Processes,” ASME J. Vib. Acoust., 135(5), p. 051006. [CrossRef]
Zhou, K. , and Tang, J. , 2016, “ Highly Efficient Probabilistic Finite Element Model Updating Using Intelligent Inference With Incomplete Modal Information,” ASME J. Vib. Acoust., 138(5), p. 051016. [CrossRef]
Zhou, K. , Liang, G. , and Tang, J. , 2016, “ Component Mode Synthesis Order-Reduction for Dynamic Analysis of Structure Modeled With NURBS Finite Element,” ASME J. Vib. Acoust., 138(2), p. 021016. [CrossRef]
Chan, Y. J. , and Ewins, D. J. , 2010, “ Management of Variability of Vibration Response Levels in Mistuned Bladed Discs Using Robust Design Concept. Part 2: Tolerance Design,” Mech. Syst. Signal Process., 24(8), pp. 2792–2806. [CrossRef]
Happawana, G. S. , Nwokah, O. D. I. , Bajaj, A. K. , and Azene, M. , 1998, “ Free and Forced Response of Mistuned Linear Cyclic Systems: A Singular Perturbation Approach,” J. Sound Vib., 211(5), pp. 761–789. [CrossRef]
Kaminski, M. , 2002, “ Stochastic Perturbation Approach to Engineering Structure Vibrations by the Finite Difference,” J. Sound Vib., 251(4), pp. 651–670. [CrossRef]
Kaplunov, J. D. , Nolde, E. V. , and Shorr, B. F. , 2005, “ A Perturbation Approach for Evaluating Natural Frequencies of Moderately Thick Elliptic Plates,” J. Sound Vib., 281(3–5), pp. 905–919. [CrossRef]
Wei, S. T. , and Pierre, C. , 1988, “ Localization Phenomena in Mistuned Assemblies With Cyclic Symmetry Part I: Free Vibrations,” J. Vib., Acoust. Stress Reliab. Des., 110(4), pp. 429–438. [CrossRef]
Wei, S. T. , and Pierre, C. , 1988, “ Localization Phenomena in Mistuned Assemblies With Cyclic Symmetry Part II: Forced Vibrations,” J. Vib., Acoust. Stress Reliab. Des., 110(4), pp. 439–449. [CrossRef]
Zhang, Y. Q. , and Wang, W. L. , 1995, “ Eigenvector Derivatives of Generalized Non-Defective Eigenproblems With Repeated Eigenvalues,” ASME J. Eng. Gas Turbines Power, 117(1), pp. 207–212. [CrossRef]
Feiner, D. M. , and Griffin, J. H. , 2002, “ A Fundamental Model of Mistuning for a Single Family of Modes,” ASME J. Turbomach., 124(4), pp. 597–605. [CrossRef]
Feiner, D. M. , and Griffin, J. H. , 2004, “ Mistuning Identification of Bladed Disks Using a Fundamental Mistuning Model—Part I: Theory,” ASME J. Turbomach., 126(1), pp. 150–158. [CrossRef]
Feiner, D. M. , and Griffin, J. H. , 2004, “ Mistuning Identification of Bladed Disks Using a Fundamental Mistuning Model—Part II: Application,” ASME J. Turbomach., 126(1), pp. 159–165. [CrossRef]
Yoo, H. H. , Kim, J. Y. , and Inman, D. J. , 2003, “ Vibration Localization of Simplified Mistuned Cyclic Structures Undertaking External Harmonic Force,” J. Sound Vib., 261(5), pp. 859–870. [CrossRef]
Kenyon, J. A. , and Griffin, J. H. , 2003, “ Experimental Demonstration of Maximum Mistuned Bladed Disk Forced Response,” ASME J. Turbomach., 125(4), pp. 673–681. [CrossRef]
Tournour, M. A. , Atalla, N. , Chiello, O. , and Sgard, F. , 2001, “ Validation, Performance, Convergence and Application of Free Interface Component Mode Synthesis,” Comput. Struct., 79(20–21), pp. 1861–1876. [CrossRef]

Figures

Grahic Jump Location
Fig. 1

(a) Schematic diagram of bladed disk assembly and its substructure (N substructures in total) and (b) finite-element model of the substructure (n DOFs)

Grahic Jump Location
Fig. 2

Natural frequencies of first modal family: (a) validation of order reduction technique by comparison and (b) frequency bandwidth demonstration (zoom-in view)

Grahic Jump Location
Fig. 3

Illustration of design optimization: the response variance is reduced while the response mean does not increase (:original; – –with design modification)

Grahic Jump Location
Fig. 4

Conventional brute-force Monte Carlo-based design flow

Grahic Jump Location
Fig. 5

Complete Monte Carlo analysis result: (a) design objective (variance) distribution and (b) mean summation (constraint) distribution

Grahic Jump Location
Fig. 6

Enhanced design flow: (a) with single Gaussain process emulator and (b) with nested Gaussian process emulator

Grahic Jump Location
Fig. 7

Gaussian process emulation result (590 training data under single GP): (a) design objective variance) distribution and (b) mean summation (constraint) distribution

Grahic Jump Location
Fig. 8

Comparison of design objective evaluation versus design variable A in projection view: (a) complete Monte Carlo analysis and (b) Gaussian process emulator (590 training data under single GP)

Grahic Jump Location
Fig. 9

Comparison of design objective evaluation versus design variable τ in projection view: (a) complete Monte Carlo analysis and (b) Gaussian process emulator (590 training data under single GP)

Grahic Jump Location
Fig. 10

Error distribution of Gaussian process emulator (590 training data under single GP)

Grahic Jump Location
Fig. 11

(a) Design objective evaluation versus design variable A in projection view and (b) design objective evaluation versus design variable τ in projection view (120 training data under single GP)

Grahic Jump Location
Fig. 12

Error distribution of Gaussian process emulator (120 training data under single GP)

Grahic Jump Location
Fig. 13

Error distribution of Gaussian process emulator (51 training data under single GP)

Grahic Jump Location
Fig. 14

Prediction of response distribution under uncertainties at one specific excitation frequency (——: benchmark distribution by MC;·– –·–: inferred distribution by GP; *: training data)

Grahic Jump Location
Fig. 15

Statistical distribution of response: (a) complete Monte Carlo analysis and (b) Gaussian process emulator (67 training data points)

Grahic Jump Location
Fig. 16

(a) Comparison of prediction accuracy and (b) comparison of computational efficiency

Grahic Jump Location
Fig. 17

Scenario: A (0–4%) (a) design objective evaluation versus design variable A in projection view and (b) design objective evaluation versus design variable τ in projection view (590 training data under single GP)

Grahic Jump Location
Fig. 18

Scenario: A (0–5%) (a) design objective evaluation versus design variable A in projection view and (b) design objective evaluation versus design variable τ in projection view (590 training data under single GP)

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In