3. Rappels#
3.1. Problème d’optimisation associé à une fonctionnelle énergétique#
Cet opérateur résout le problème suivant :
Trouver le triplet de champs admissibles \((u,v,w)\text{}\) qui minimisent la fonctionnelle:
\({e}_{\omega}^{2}(u,v,w)=\frac{\gamma}{2}{(u-v)}^{T}[K](u-v)+\frac{1-\gamma }{2}{(u-w)}^{T}{\omega}^{2}[M](u-w)+\frac{1-\alpha }{\alpha}{(\mathit{Hu}-\widehat{u})}^{T}\text{}[\mathit{Gr}](\mathit{Hu}-\widehat{u})\)
sous la contrainte:
\(\left[K\right]v-{\omega}^{2}\left[M\right]\text{w=}0\)
Où :
\(K\) représente une matrice de rigidité réelle
\(M\) représente une matrice de masse
\(H\) représente une matrice d’observation
\(\mathit{Gr}\) représente une matrice symétrique définie positive servant de norme des erreurs dans l’espace d’observation
\(\omega =2\pi f\) : pulsation d’excitation
\(\widehat{u}\) observation des déplacements à la pulsation \(\omega\)
\(\gamma\) paramètre de pondération des erreurs \((u-v)\) et \((u-w)\)
\(\alpha\) paramètre de pondération de la fonctionnelle assimilable à un terme de régularisation
3.2. Équations de résolution du problème#
L’obtention du triplet \((u,v,w)\text{}\) associé au problème sous contrainte plus haut amène, à chaque fréquence considérée, à la résolution du système d’équations linéaire suivant:
\(Al=b\)
avec, pour chaque pulsation propre \({\omega}_{i}\) :
\({A}_{i}=\left(\begin{array}{cc}\gamma \left(\text{K+γ}/\left(1-\gamma \right){\omega}_{i}^{2}M\right)& -\gamma \left(K-{\omega}_{i}^{2}M\right)\\ -\gamma \left(K-{\omega}_{i}^{2}M\right)& \left(-\text{2α}/\left(1-\alpha \right)\right){H}^{T}{G}_{r}H\end{array}\right)\) ; et \({b}_{i}=\left(\begin{array}{c}{0}_{n}\\ \left(-\text{2α}/\left(1-\alpha \right)\right){H}^{T}{G}_{r}{\widehat{u}}_{i}\end{array}\right)\)
Pour davantage de détails sur la formulation, on se reportera au document [R4.10.07].