Bayesian operational modal analysis: Difference between revisions

Revision as of 23:26, 5 August 2020

Latest revision as of 21:18, 28 January 2023

(9 intermediate revisions by 7 users not shown)

Line 1:

{{refimprovemore citations needed|date=December 2013}}

'''Bayesian operational modal analysis (BAYOMA)''' adopts a [[Bayesian inference|Bayesian]] [[system identification]] approach for [[Operationaloperational Modalmodal Analysisanalysis]] (OMA). [[Operational Modalmodal Analysis]] (OMA)analysis aims at identifying the modal properties ([[natural frequency|natural frequencies]], [[damping ratio]]s, [[mode shape]]s, etc.) of a constructed structure using only its (output) vibration response (e.g., velocity, acceleration) measured under operating conditions. The (input) excitations to the structure are not measured but are assumed to be '[[Ambient vibrations|ambient]]' ('broadband random'). In a Bayesian context, the set of modal parameters are viewed as uncertain parameters or random variables whose probability distribution is updated from the prior distribution (before data) to the posterior distribution (after data). The peak(s) of the posterior distribution represents the most probable value(s) ('''MPV''') suggested by the data, while the spread of the distribution around the MPV reflects the remaining uncertainty of the parameters.

==Pros and cons==

Line 27:

|author2=Katafygiotis, L.S.

|title=Bayesian spectral density approach for modal updating using ambient data

|journal=Earthquake Engineering and& Structural Dynamics

|year=2001

|volume=30

|issue=8

|pages=1103–1123

|doi=10.1002/eqe.53|s2cid=110355068
}}</ref> and [[fast Fourier transform]] (FFT)<ref>

{{cite journal

|last=Yuen

Line 43 ⟶ 44:

|issue=2

|pages=81–95

|doi=10.1260/136943303769013183|s2cid=62564168
}}</ref> of ambient vibration data. Based on the formulation for FFT data, fast algorithms have been developed for computing the posterior statistics of modal parameters .<ref name="bayomabook" />. Recent developments based on [[EM algorithm]] <ref>

{{cite journal

|last=Li

Line 51 ⟶ 53:

|journal=Mechanical Systems and Signal Processing

|year=2019

|volume=132

|pages=490–511

|doi=10.1016/j.ymssp.2019.06.036

|bibcode=2019MSSP..132..490L

|url=}}</ref> show promise for simpler algorithms and reduced coding effort. The fundamental precision limit of OMA has been investigated and presented as a set of '''uncertainty laws''' which can be used for planning ambient vibration tests <ref name=ulaw2018>

|hdl=10356/149983

|s2cid=199124928

|hdl-access=free

|url=}}</ref> show promise for simpler algorithms and reduced coding effort. The fundamental precision limit of OMA has been investigated and presented as a set of '''uncertainty laws''' which can be used for planning ambient vibration tests .<ref name=ulaw2018>

{{cite journal

|last=Au

Line 61 ⟶ 69:

|journal=Mechanical Systems and Signal Processing

|year=2018

|volume=102

|pages=139–157

|doi=10.1016/j.ymssp.2017.09.017

|bibcode=2018MSSP..102..139A

|url=}}</ref>.

|hdl=10871/30384

|hdl-access=free

|url=}}</ref>.

==Connection with [[maximum likelihood method]]==

Bayesian method and [[maximum likelihood method]] (non-Bayesian) are based on different philosophical perspectives but they are mathematically connected; see, e.g., <ref>

{{cite journal

|last=Au

Line 73 ⟶ 86:

|journal = Mechanical Systems and Signal Processing

|year=2017

|volume=25

|issue=3

|pages=e2113

|doi=10.1002/stc.2113

|s2cid=55868193

|doi-access=free

}}

</ref> and Section 9.6 of .<ref name="bayomabook" />. For example,

*Assuming a uniform prior, the most probable value (MPV) of parameters in a Bayesian method is equal to the ___location where the likelihood function is maximized, which is the estimate in Maximum Likelihood Method

*Under a Gaussian approximation of the posterior distribution of parameters, their covariance matrix is equal to the inverse of Hessian of the negative log of likelihood function at the MPV. Generally, this covariance depends on data. However, if one assumes (hypothetically; non-Bayesian) that the data is indeed distributed as the likelihood function, then for large data size it can be shown that the covariance matrix is asymptotically equal to the inverse of the [[Fisher information]] matrix (FIM) of parameters (which has a non-Bayesian origin). This coincides with the [[Cramer–Rao bound]] in classical statistics, which gives the lower bound (in the sense of matrix inequality) of the ensemble variance of any unbiased estimator. Such lower bound can be reached by maximum-likelihood estimator for large data size.

*In the above context, for large data size the asymptotic covariance matrix of modal parameters depends on the 'true' parameter values (a non-Bayesian concept), often in an implicit manner. It turns out that by applying further assumptions such as small damping and high signal-to-noise ratio, the covariance matrix has mathematically manageable asymptotic form, which provides insights on the achievable precision limit of OMA and can be used to guide ambient vibration test planning. This is collectively referred as 'uncertainty law'.<ref name="ulaw2018"/>.

==See also==

*[[Operational Modalmodal Analysisanalysis]]

*[[Bayesian inference]]

*[[Ambient Vibrationsvibrations]]

*[[Microtremor]]

*[[Modal analysis]]

*[[Modal testing]]

==Notes==

Line 89 ⟶ 115:

|publisher=Kluwer Academic Publisher

|___location=Boston

|url=}}</ref> <ref>

|isbn=

|url=}}</ref> <ref>

{{cite book

|first=M. |last=Schipfors

Line 97 ⟶ 122:

|year= 2014

|publisher=Springer

|isbn=

|url=https://www.springer.com/gp/book/9781493907663}}

</ref> <ref>

{{cite book

|first=R. |last=Brincker

Line 106 ⟶ 130:

|year= 2015

|publisher=John Wiley & Sons

|doi=10.1002/9781118535141

|isbn=9781118535141

|url=https://onlinelibrary.wiley.com/doi/book/10.1002/9781118535141}}

</ref> and Bayesian OMA <ref name=bayomabook>

Line 127 ⟶ 153:

|publisher=Cambridge University Press

|___location=United Kingdom

}}

|isbn=

|url=}}

</ref> and Cox<ref name=cox>

{{cite book

Line 136 ⟶ 161:

|publisher=Johns Hopkins University Press

|___location=Baltimore

}}

|isbn=

|url=}}

</ref> for Bayesian inference in general.

*See Beck<ref>

Line 150 ⟶ 174:

|pages=825–847

|doi=10.1002/stc.424

|s2cid=122257401

|url=}}</ref> for Bayesian inference in structural dynamics (relevant for OMA)

|doi-access=free

|url=}}</ref> for Bayesian inference in structural dynamics (relevant for OMA)

*The uncertainty of the modal parameters in OMA can also be quantified and calculated in a non-Bayesian manner. See Pintelon et al.<ref>

Line 164 ⟶ 190:

|pages=2359–2373

|doi=10.1016/j.ymssp.2006.11.007

|url=|bibcode=2007MSSP...21.2359P

}}</ref>

==See also==

*[[Operational Modal Analysis]]

*[[Bayesian inference]]

*[[Ambient Vibrations]]

*[[Microtremor]]

*[[Modal analysis]]

*[[Modal testing]]

==References==