Alain Cimetière

Franck Delvare

Frédéric Pons

Comptes Rendus Mécanique (CRM)

/uccessivesité àmscy ofsune partie).C. R. Mecanique 333 (2005) 123–126http://france.elsevier.com/direct/CRAS2BUne méthode inverse d’ordre un pour les problèmesde complétion de donnéesAlain Cimetièrea,∗, Franck Delvareb, Frédéric Ponsca LMP, université de Poitiers, boulevard Marie et Pierre Curie, téléport 2, BP 30179, 86962 Futuroscope Chasseneuil cedex, Franceb LEES, université d’Orléans et ENSI de Bourges, 10, boulevard Lahitolle, 18020 Bourges cedex, Francec LEA, université de Poitiers et ENSMA, boulevard Marie et Pierre Curie, téléport 2, BP 30179, 86962 Futuroscope Chasseneuil cedex, FranceReçu le 7 juillet 2004 ; accepté après révision le 30 novembre 2004Disponible sur Internet le 7 janvier 2005Présenté par Huy Duong BuiRésuméUne méthode itérative est présentée pour la résolution des problèmes inverses de type Cauchy. Les itérations svérifient exactement les équations d’équilibre. Les tests numériques attestent de la précision de la méthode et de sa capacrésoudre efficacement ces problèmes lorsque le bord du domaine n’est pas régulier.Po r citer cet article : A. Cimetière et al.,C. R. Mecanique 333 (2005). 2004 Académie des sciences. Publié par Elsevier SAS. Tous droits réservés.AbstractA first order inverse method for the data completion problem.An iterative resolution method for inverse Cauchy probleis presented. The successive iterations satisfy the equilibrium equations exactly. Numerical simulations prove the accurathe method and its ability to solve Cauchy problems when thedomain boundary is not regular.To cite this article: A. Cimetièreet al., C. R. Mecanique 333 (2005). 2004 Académie des sciences. Publié par Elsevier SAS. Tous droits réservés.Mots-clés : Solides et structures ; Problème inverse; Équation de Laplace ; Complétion de donnéesKeywords: Solids and structures; Inverse problem; Laplace’s equation; Data completion1. IntroductionCertaines techniques d’identification nécessitent de connaître à la fois déplacementset forces ou températureet flux de chaleur sur toute la frontière d’un domaine, alors que ces quantités ne sont accessibles que sur* Auteur correspondant.Adresses e-mail : alain.cimetiere@univ-poitiers.fr (A. Cimetière), franck.delvare@ensi-bourges.fr (F. Delvare), pons@ensma.fr (F. Pons1631-0721/$ – see front matter 2004 Académie des sciences. Publié par Elsevier SAS. Tous droits réservés.doi:10.1016/j.crme.2004.11.009124 A. Cimetière et al. / C. R. Mecanique 333 (2005) 123–126plétional posé deur ont étée d’identi-hodesée àctementmes peut êtresente unurs sontinuités deméthodeéliorationdomainededretionclassede celle-ci. Voir par exemple [1], en ce qui concerne l’identification de fissures. Les problèmes de comde données sont des problèmes de Cauchy pour des équations aux dérivées partielles. Le caractère mces problèmes rend leur résolution très délicate [2]. Déjà, dans les années 60, de nombreux travaux leconsacrés [4–7].Une technique de complétion de données a été développée dans [8] et mise en oeuvre dans un problèmfication de fissures [9]. L’opérateur modèle choisi était le laplacien. Cette technique, appelée dans la suite métd’ordre zéro, s’apparente à une méthode de Tikhonov itérative [3], où la connaissance a priori est réactualichaque itération. La particularité de la technique est qu’à chaque itération la solution déterminée vérifie exa« l’équilibre ». Lorsque les données sont compatibles, la suite des itérations converge vers la solution du problède Cauchy. Des niveaux de bruit importants sur les données sont très bien tolérés et le nombre de donnéetrès inférieur au nombre de quantités à déterminer [8].La reconstruction de la dérivée normale sur la frontière du domaine par la méthode d’ordre zéro préniveau de précision bien inférieur à celui obtenu pour la reconstruction de la fonction elle-même (les erreenviron dix fois plus grandes). On note également que la méthode d’ordre zéro gère assez mal les discontdérivée normale, surtout dans les coins.L’objet de cette note est de présenter une évolution de la méthode d’ordre zéro, que l’on appellerad’ordre un. Ses deux principaux apports sont une amélioration sensible de la convergence et une amtrès nette de la reconstruction des solutions, en particulier de la dérivée normale, lorsque la frontiére duprésente des points anguleux. La méthode de quasi-réversibilité [5] est inadaptée pour gérer les problèmescomplétion de données, mais on s’inspire de l’idée que les équilibres sont aussi solutions d’un problème d’orsupérieur, pour construire la méthode d’ordre un.2. La méthode d’ordre unSoitΩ un domaine borné deR2, de frontièreΓ = Γd ∪ Γi . Soientφd etψd deux fonctions données surΓd . Onse propose de résoudre le problème de Cauchy suivant :u = 0 dansΩu = φd surΓdu′ = ψd surΓd(1)où u′ désigne la dérivée normale deu. Le problème (1) admet au plus une solution. L’existence d’une solun’est assurée que si les donnéesφd etψd sont compatibles, c’est à dire sont respectivement les traces surΓd d’unefonction harmonique et de sa dérivée normale. Dans la suite, cette fonction harmonique est supposée deC2sur Ω . On introduit maintenant le problème suivant :u = 0, u1 = 0, u2 = 0 dansΩu = φd, u1 = φ1d, u2 = φ2d surΓdu′ = ψd, u′1 = ψ1d , u′2 = ψ2d surΓdu′ = u1n1 + u2n2 surΓ∫Γ (u1n1 + u2n2)ds = 0(2)où n1 et n2 sont les composantes de la normale unitaire extérieure àΩ . Les quantitésφ1d , φ2d , ψ1d , ψ2d sontévaluées à partir des donnéesφd , ψd , en écrivant queu1 et u2 sont les deux dérivées partielles deu t en utilisantl’équationu = 0 surΓd . Ainsi les valeurs deu1, u2, u′1, u′2 se trouvent données par les relations :φ1d = ψdn1 − dφddsn2, φ2d = ψdn2 + dφddsn1ψ1d = −dφ2dds , ψ2d =dφ1dds(3)A. Cimetière et al. / C. R. Mecanique 333 (2005) 123–126 125risants.l’algo-l-tion depastire deLeiretpartière nonnnées est lectionionsruit. DansLorsque les donnéesφd et ψd sont compatibles, toutes les données du problème (2) sont compatibles etu1, u2sont alors les dérivées partielles deu. Le problème (2) peut être reformulé en distinguant les relations caractél’équilibre :u = 0, u1 = 0, u2 = 0 dansΩ,∫Γ(u1n1 + u2n2)ds = 0 (4)et les relations traduisant les conditions aux limites :u = φd, u1 = φ1d, u2 = φ2d surΓdu′ = ψd, u′1 = ψ1d, u′2 = ψ2d surΓdu′ = u1n1 + u2n2 surΓ(5)La solution de ce problème estrecherchée dans l’espaceH(Ω) défini par :H(Ω) = {(v, v1, v2) ∈ [H 1(Ω)]3 vérifiant les équations d’équilibre (4)}On définit ensuiteH(Γ ), l’espace composé des traces surΓ des éléments deH(Ω) et de leurs dérivées normalePourΦd = (φd,φ1d,φ2d,ψd,ψ1d,ψ2d ), une formulation équivalente au problème (2) est :TrouverU = (u,u1, u2, u′, u′1, u′2) ∈ H(Γ ) tel que :U = Φd surΓdu′ = u1n1 + u2n2 surΓ(6)L’algorithme itératif suivant est introduit pour résoudre (2). Il est appelé méthode d’ordre un et, commerithme proposé dans [8], reprend l’esprit de la méthode itérative de Tikhonov : étant donnéc > 0 etU0 ∈ H(Γ ),trouverUk+1 ∈ H(Γ ) tel que :J kc (Uk+1)  J kc (V ) ∀V ∈ H(Γ ) avec :J kc (V ) =∥∥v′ − (v1n1 + v2n2)∥∥2−1/2,Γ + ‖V − Φd‖2Γd + c‖V − Uk‖2H(Γ )(7)Supposant que le problème (1) possède une solution de classeC2 ur Ω , la preuve de la convergence de l’agorithme est similaire à celle établie dans [8] pour l’algorithme de la méthode d’ordre zéro. La discrétisal’algorithme conduit à chaque itération à la résolution d’un système linéaire pour lequel la matrice ne dépenddek. La condition d’équilibre global, traduite par la dernière égalité de (4), joue un rôle capital. Si on la rela définition deH(Ω), la méthode d’ordre un donne alors des résultats bien moins précis.3. Tests numériquesLes simulations numériques ont pour but de comparer lméthode d’ordre un à la méthode d’ordre zéro [8].domaine considéré est donné sur la Fig. 1, la fonction à reconstruire est :u(x1, x2) = x2(x1 − 1/2)2 + x22Le supportΓd des données est composé de l’arc de cercle de rayon 1. Les donnéesφd etψd sont générées à partde la fonction solution. Les données concernantu1 et u2 sont évaluées numériquement à partir des valeurs dφdet ψd , ces valeurs pouvant être bruitées. Les dérivées tangentiell s figurant dans (3) ontété calculées en utilisanune technique de régression parabolique sur 11 points. LaFig. 2 donne la reconstruction de la dérivée normalela méthode d’ordre un. On constate que la reconstruction devient précise lorsque le domaine est à fronrégulière. On s’intéresse désormais au cas de données entachées d’un bruit de mesure. Le support des domême que précédemment. Les donnéesφd etψd sont entachées d’un bruit de 1%. La Fig. 3 donne la reconstrude la dérivée normale par la méthode d’ordre zéro et la Fig. 4 celle par la méthode d’ordre un. Ces reconstructmontrent que la méthode d’ordre un fournit des résultats beaucoup plus précis, surtout en présence de bces tests, le nombre d’inconnues est le double de celui des données.126 A. Cimetière et al. / C. R. Mecanique 333 (2005) 123–126ormaleméthodefrontièreérovailler enode à710Fig. 1. Domaine d’étude. Fig. 2. Reconstruction deu′ surΓi (ordre un).Fig. 3. Reconstruction deu′ surΓi (ordre zéro), bruit suru etu′ de 1%.Fig. 4. Reconstruction deu′ surΓi (ordre un), bruit suru etu′de 1%.4. ConclusionLa méthode d’ordre un permet d’obtenir un niveau de précision pour la reconstruction de la dérivée ncomparable à celui obtenu pour la reconstruction de la fonction elle-même par la méthode d’ordre zéro. Lad’ordre un est capable de reconstruire des fonctions et leurs dérivées normales lorsque le domaine est ànon régulière, en restituant bien les discontinuités de dérivée normale. La comparaison à la méthode d’ordre zatteste de l’efficacité et de la robustesse de la méthode d’ordre un, en particulier de sa capacité à traprésence de bruit. Par ailleurs, ce test ainsi que d’autres essais numériques ont montré la capacité de la méthdébruiter.5. RemerciementsLes auteurs expriment leur gratitude à H.D. Bui ainsi qu’aux rapporteurs pour leurs suggestions.Références[1] S. Andrieux, A. Ben Abda, H.D. Bui, Reciprocity principle and crack identification, Inverse Probl. 15 (1999) 59–65.[2] H.D. Bui, Introduction aux problèmes inverses en mécanique des matériaux, Eyrolles, 1993.[3] H.W. Engl, M. Hanke, A. Neubauer, Regularisation of Inverse Problems, Kluwer Academic, 1996.[4] M.M. Laurentiev, Some Improperly Posed Problems of Mathematical Physics, Springer-Verlag, 1967.[5] R. Lattès, J.L. Lions, Méthode de quasi-reversibilité et applications, Dunod, 1967.[6] A. Bryson, Y.C. Ho, Applied Optimal Control, Blaisdell, 1969.[7] A. Tikhonov, V. Arsénine, Méthode de résolution de problèmes mal posés, Mir, 1976.[8] A. Cimetière, F. Delvare, M. Jaoua, F. Pons, Solution of the Cauchy problem using iterated Tikhonov regularization, Inverse Probl. 1(2001) 553–570.[9] A. Cimetière, F. Delvare, M. Jaoua, M. Kallel, F. Pons, Reconstruction of cracks from incompleteboundary data, Inverse Probl. Eng.(2002) 377–392.

Une méthode inverse d'ordre un pour les problèmes de complétion de données

https://comptes-rendus.academie-sciences.fr/mecanique/item/10.1016/j.crme.2004.11.009.pdf

Une méthode inverse d'ordre un pour les problèmes de complétion de données

Abstract

Similar works

Full text

Available Versions

Comptes Rendus Mécanique (CRM)