Publications

2002

• Caractérisation Biomécanique et Modélisation 3D par Imagerie X et IRM haute résolution de l'os spongieux humain : Evaluation du risque fracturaire
  
  • Follet Hélène
  , 2002. Dans le cadre de la prévention des pathologies osseuses, l'étude présentée a pour but de tester une nouvelle méthodologie d'évaluation du risque fracturaire de l'os spongieux de calcanéum humain. En s'appuyant i) sur des données cliniques (de densité minérale osseuse, de microstructure, DXA, Scanner X, Histomorphométrie) et ii) sur des techniques d'imagerie haute définition (IRM à 78 µm et tomographie à 10 µm), l'objectif est d'estimer les propriétés mécaniques de l'os spongieux (Module d'Young et contrainte maximale de compression). Deux types d'essais mécaniques ont été mis en œuvre : essais de compression sur échantillons cubiques d'os spongieux, essai de micro-flexion sur trabécules osseuses. Différents modèles par éléments finis (MEF) de ces essais ont été construits et permettent de déterminer les propriétés mécaniques du tissu trabéculaire. Le risque fracturaire pourrait être évalué par quantification de la charge de ruine et du degré de déformation de ce tissu. Les résultats de cette méthodologie sont alors confrontés aux méthodes classiques d'évaluation clinique du risque fracturaire.
• Développement et mise en oeuvre de méthodes optiques pour la mesure de relief et de champ de déformations en vue de la modélisation d'organes biologiques
  
  • Brèque Cyril
  , 2002. Cette étude porte sur l'application des techniques de mesure optiques pour la caractérisation mécanique des organes biologiques. Deux voies ont été explorées : la réalisation de maillages en éléments finis et la mesure de déformations planes. Les maillages en éléments finis peuvent être réalisés à partir du relief des organes. Nous avons alors développé une technique de projection de réseaux de traits. Afin d'améliorer l'exactitude, l'analyse des réseaux s'effectue par des techniques de décalages de phases. Une procédure de calibrage du montage a été spécialement développée. La mesure des déformations sur de la peau a été effectuée par une extension de la technique de suivi de marqueurs. Afin de voir l'ensemble des phénomènes mécaniques qui régissent ce tissu, nous avons réalisé des essais de traction simple, de relaxation, de traction à accélération constante et de traction biaxiale. Le caractère viscoplastique de la peau et l'orthotropie de sa structure ont été ainsi mis en évidence.
• Tomographie aux rayons X de l'écoulement des matériaux granulaires ensilés
  
  • Luong Minh-Phong
  , 2002. Cette communication propose l’utilisation de la tomographie aux rayons X pour mettre en évidence l'intérêt d’un indicateur de l'effet de voûte, appelé indice de coulabilité des matériaux granulaires, facilement obtenu à partir des essais triaxiaux conventionnels bien connus en mécanique des sols. En extension, cette analyse permet d’interpréter correctement le bouchonnage des pieux tubulaires à fond ouvert observé en géotechnique pétrolière.
• Sur l'identification des paramètres des lois de comportement des roches argileuses.
  
  • Lecampion Brice
  , 2002. Ce travail a pour but de mettre au point des méthodes d'identification des paramètres constitutifs des roches argileuses. Dans le cadre du projet de laboratoire de recherche souterrain de l'ANDRA, il est nécessaire de développer ce type de méthodes en vue de l'interprétation des nombreuses mesures qui seront pratiquées sur site. On s'est focalisé sur deux aspects rhéologiques majeurs du comportement de ce type de roche: le comportement poroélastique d'une part et le comportement élasto-viscoplastique d'autre part. La première partie se focalise sur l'identification des paramètresporoélastiques. Le chapitre 2 rappelle le problème direct et discute un certain nombre de points importants concernant le problème inverse d'identification. Le troisième chapitre est dédié à la formulation des techniques de calcul de gradient pour le cas poroélastique linéaire. La résolution numérique par éléments finis est discutée. Les méthodes de Différentiation directe et de l'état adjoint sont validées sur un exemple numérique bidimensionnel à l'aide du code de calcul par éléments finis Cast3M. Uidentificationdes coefficients poroélastiques des roches argileuses de Meuse Haute-Marne à partir d'essais de laboratoire est discutée en détail dans le chapitre 4. L'utilisation de solutions semi-explicites approchées des problèmes directs permet d'obtenir une méthode d'identification rapide. La deuxième partie du mémoire porte sur l'identification des paramètres élasto-viscoplastiques. Les traits ducomportementviscoplastiques des roches argileuses de Meuse Haute-Marne sont discutés dans le chapitre 5 sur la base des résultats expérimentaux. La modélisation de ce comportement est envisagée. On propose notamment un modèle viscoplastique avececrouissage isotrope non linéaire permettant de reproduire les essais. Les paramètres de cette loi de comportement sont identifiés sur un essai de fluage unîdîmensionnelen conditions drainées. Les déformations différées proviennent alors du comportement poroélastique et viscoplastique de la roche. On montre qu'il est possible de séparer ces deux phénomènes. Tous les paramètres poro-elastoviscoplastiques sont identifiés,une solution semi-explicite de l'essai de fluage est utilisée. Le chapitre 6 présente une méthode d'identification des paramètres élasto-viscoplastiques dans le cas général. Comme déjà présenté au cours du présent chapitre, l'identification se ramène à la minimisation d'une fonctionnelle coût. Le gradient de la fonctionnelle est calculé par Différentiation Directe. La méthode de Différentiation directe est développée dans Cast3M pour la loi élasto-visco-plastique de Norton-Hoff. L'identification des paramètres de cette loi sur des configurations unidimensionelleset bidimensionnelles d'ouvrages souterrains revêtus ou non est ensuite présentée. Les mesures sont simulées par un premier calcul direct puis éventuellement perturbées. L'unicité et la stabilité de ce problèmed'identification sont discutées. La méthode est ensuite appliquée dans le chapitre 7 à un cas réel : le tunnel de la Chamoise. Des modélisations unidimensionnelles et bidimensionnelles en déformation plane sont utilisées pour résoudre le problème direct. Les hypothèses de modélisation sont justifiées. Les paramètres obtenus permettent de reproduire convenablement les données. L'étude de ce cas réel permet également d'émettre un certain nombre de recommandations sur les dispositifs de mesure in situ en vue de l'identification des paramètres des lois de comportement.
• Dynamique optimale de systèmes articulés à cinématique fermée. Application à la synthèse d'allures de marche optimales
  
  • Chessé Stéphane
  , 2002. Ce mémoire est centré sur l'optimisation dynamique de mouvements de systèmes articulés plans travaillant en boucle fermée. L'application visée est l'optimisation de pas de marche de robots bipèdes. La première partie est consacrée au développement de modèles dynamiques de systèmes bouclés adaptés au problème de l'optimisation lorsqu'on fait appel au principe du maximum de Pontryagin. Le problème d'optimisation est traité dans la deuxième partie. La démarche utilisée consiste à ouvrir la chaîne au niveau d'une liaison bien choisie et à traiter les efforts de liaison comme des efforts de commande complémentaires. La condition de fermeture de la chaîne est alors considérée comme une contrainte du problème d'optimisation et est traitée par une méthode de pénalité. Le troisième chapitre contient des simulations numériques relatives à des mécanismes plans suractionnés de type robots parallèles. Les principaux résultats concernent l'optimisation de pas de marche de bipèdes plans à pattes bisegmentaires et trisegmentaires.
• Identification of poroelastic constants of deep argillaceous rocks. II: Inverse analysis
  
  • Lecampion Brice
  • Constantinescu Andrei
  , 2002 (1), pp.245-250. This paper discusses the identification of poroelastic constants from drained isotropic confinement and pulse test. Closed form solutions of the direct problem are obtained for a very thin specimen. The validity of these solutions for a realistic aspect ratio of the core is assessed by a comparison with a finite element model. The identification problem is solved by minimizing a least square functional using an explicit gradient and a Levenberg-Marquardt algorithm. Uniqueness of this inverse problem as well as the effect of noise on input data are fully discussed. The identification procedure is then applied to measurements performed on a deep argillaceous rock (argillite of Meuse Haute-Marne). (10.1201/9781003078807-37)
  DOI : 10.1201/9781003078807-37
• Analyse des ondes dans les barres élastiques et viscoélastiques par la méthode du maximum de vraisemblance
  
  • Othman Ramzi
  • Blanc Robert
  • Bussac Marie-Noëlle
  • Collet Pierre
  • Gary Gérard
  , 2002.
• Application of the likelihood method to the analysis of waves in elastic and viscoelastic rods
  
  • Othman Ramzi
  • Gary Gérard
  • Bussac Marie-Noëlle
  • Collet Pierre
  , 2002. In this paper, we are interested in separating waves in elastic and viscoelastic rods propagating in opposite direction. N strain and P velocity measurements are taken into account. This application of the likelihood method gives a solution in the frequency domain. Using the inverse Fourier transform, one can recover the strain, stress, displacement and velocity at any section of the rod. In experimental conditions, the results are stable against noise when N+P>2 and NP ≠ 0. (10.1142/9789812704405_0035)
  DOI : 10.1142/9789812704405_0035
• Dispersion identification using the Fourier analysis of resonances in elastic and viscoelastic rods
  
  • Othman Ramzi
  • Gary Gérard
  • Blanc Robert H
  • Bussac Marie-Noëlle
  • Collet Pierre
  , 2003, pp.229-235. A new method for identifying the dispersion relation in elastic and viscoelastic rods is presented. It relies on the multiple resonance of a strain measurement. The wave velocity is related to the resonance position whereas the damping is related to the resonance bandwidth. Applied to an aluminum bar, the method provides wave dispersion for frequencies up to 60 kHz. (10.1142/9789812704405_0034)
  DOI : 10.1142/9789812704405_0034
• Détermination de la dispersion à partir de l'analyse de Fourier des résonances dans les barres élastiques et viscoélastiques
  
  • Ohtman Ramzi
  • Blanc Robert
  • Bussac Marie-Noëlle
  • Gary Gérard
  • Collet Pierre
  , 2002.
• Determination of the electronic density in a medium by an inverse method based on double-Compton scattering in transmission imaging
  
  • Truong Tuong Trong
  • Nguyen Maï K.
  • Bui Huy Duong
  • Daveau Christian
  , 2002, pp.405-412. We propose a new method for reconstructing defects in bulk materials, via the recovery of the electronic volume density ne in the material, using doubly Compton-scattered gamma photons in transmission imaging modality. We establish an integral relationship between the photon flux density after a double scattering and the product of electronic densities at two different sites. This photon flux density at different scattering angles can be measured by an energy-position detector in a fixed configuration. Thus the reconstruction of ne can be formulated as an inverse problem of the integral relation, which is bilinear in the electronic densities. This theoretical result opens the way for a new imaging principle, which exploits scattered radiation rather than discarding it as in most existing imaging procedures. In this new procedure, the motion of the detector is no longer necessary as it is the case of conventional tomography.
• Mode II fatigue crack growth in a maraging steel. Frequency and environment effects
  
  • Doquet Véronique
  Revue de Métallurgie, EDP Sciences, 2002, 99 (5), pp.441 - 448. (10.1051/metal:2002163)
  DOI : 10.1051/metal:2002163
• LOSAC: a first salt cavern abandonment software
  
  • Brouard Benoit
  • Berest Pierre
  • Durup Gérard
  , 2002, pp.29-48. Over recent years, considerable research effort, especially by the SMRI, has been invested on the subject of the long-term abandonment of salt caverns. At this point in time, the industry has still not clearly identified a universally accepted concept or approach to this major problem. One seemingly promising trend already identified is to consider the evolution of a saturated brine filled salt cavern with respect to the system’s thermal (brine expansion), geomechanical (cavern closure by creep) and hydraulic (brine micro-permeation in the salt massif) characteristics. Gaz de France has applied such a concept in the development of a computer programme named LOSAC intended to evaluate the long-term evolution of its salt caverns in different abandonment scenarios. LOSAC is financed by Gaz de France and developed with Brouard Consulting and LMS/Ecole polytechnique (France). This presentation proposes a short demonstration of the software with the aim of illustrating a practical application of this specific cavern abandonment concept in a very user-friendly computer program.
• Tightness tests in salt-cavern wells
  
  • Bérest Pierre
  • Brouard Benoit
  • Durup Gérard
  , 2002, 56 (5), pp.1-27. Thousands of caverns have been leached out from salt formations. They are used for storing a variety of fluid products ranging from compressed air and hydrogen to LPG, natural gas and crude oil, which requires that the caverns be tight. The main factors in the onset of well leakage and its prevention are discussed: fluid pressure distribution, geological environment, cementing workmanship and well architecture. The Mont Belvieu accident is described to illustrate the importance of periodic cavern testing. Test methods are discussed; apparent, corrected and actual leaks are distinguished. Factors contributing to apparent leaks are described, as are two actual in situ tests that use fuel oil and nitrogen as test fluids. It is proven that a thorough test analysis allows good estimations of actual leaks. (10.2516/ogst:2001037)
  DOI : 10.2516/ogst:2001037
• Structures sous chargements mobiles : aspects mécaniques et fatigue
  
  • Maitournam Habibou
  , 2002. Si la compréhension et la prévision théorique des diverses dégradations affectant la tenue des structures mécaniques en contact constituent des préoccupations majeures de l'industrie et de la recherche, il n'en reste pas moins que, malgré l'importance des travaux effectués, des verrous demeurent sur la voie qui mène à ces objectifs. Qu'il s'agisse de la fissuration sous ses divers aspects ou de l'usure sous toutes ses formes, douces ou sévères, contre lesquels il est important de se prémunir, un des enjeux est toujours l'établissement de modèles fiables de prédiction. Or dans la plupart des cas, cette fiabilité ne pourrait être obtenue que par l'utilisation d'une démarche globale prenant en compte la genèse depièces et leur comportement en service. Et c'est là que résident les obstacles ; il s'agit en somme de retrouver toute l'histoire thermomécanique des pièces par la modélisation de leur procédé de fabrication (générateurs des contraintes résiduelles), le suivi de l'évolution de leur état thermomécanique sous les sollicitations variables de service, et d'appliquer les critères d'apparition et de cumul des dommages. Ainsi, lorsqu'on considère les problèmes des endommagements induits par les contacts entre solides (écaillage, fretting, fatigue,etc.) l'absence d'une démarche globale prédictive tient essentiellement à deux raisons. La première est le manque de méthodes de calculsrapides et fiables permettant d'évaluer les grandeurs thermomécaniques pertinentes responsables de l'usure ou des autres dégradations au voisinage des contacts mobiles et répétés. En effet, les méthodes de calcul classiques incrémentales y sont mal adaptées à cause du grand nombre de cycles du chargement. La seconde provient des faiblesses des critères destinés à évaluer certains phénomènes d'endommagement sous sollicitations multiaxiales, comme c'est le cas des états de contraintes induits par les contacts entre solides. Devant la profusion des propositions de critères, un premier tri s'avère déjà nécessaire pour éliminer ceux établis dans un contexte particulier et souvent faisant intervenir des paramètres dont les définitions ne sont souvent pas claires, encore moins intrinsèques. Il faudrait aussi éliminer ceux, qui, de formulations complexes, ne sont pas adaptés au calcul de structure essentiellement à cause de coût "computationnel". La démarche que nous développons vise précisément à surmonter ces deux difficultés dans un certain nombre de cas spécifiques. Les axes principaux de recherche qui, naturellement s'en dégagent, sont : (i) la formulation et la mise en oeuvre numérique d'algorithmes de calcul spécifiques pour traiter des problèmes importants (faisant essentiellement intervenir des sollicitations mobiles : rail/roue, dudgeonnage, roulements à billes, fretting, etc.) difficilement abordables par les méthodes classiques ; (ii) la formulation de lois de comportement de matériaux en relationavec les structures et les processus à étudier : divers degrés de finesse ou de sophistication sont à considérés, seuls certains phénomènes sont à prendre en compte ; selon qu'il s'agit de modéliser le soudage (acier), le rochet plastique dans le rail, le modèle de comportement serait bien sûr différent ; (ii) la recherche de critères d'apparition de défauts ; (iii) la modélisation de procédés tels que ceux d'assemblage (soudage, dudgeonnage) ou ceux destinés à introduire des contraintes résiduelles (galetage, grenaillage, trempe laser, etc.) ; elle s'avère indispensable pour une prédiction quantitative du comportement à la fatigue des structures concernées. Plus précisément, le mémoire est organisé comme suit. Dans le chapitre 1, nous abordons les méthodes de calcul des structures sous chargements cycliques. Il s'agit des cas de contacts répétés et mobiles sur de grandes distances (exemple du contact rail/roue) et de ceux de contacts intermittents ou à faibles débattements (exemple du contact essieu-axe/roue) . Nous rappelons les travaux que nous avons effectués sur la mise au point de méthodes permettant d'accéder directement et manière fiable à l'état asymptotique de structures thermoélastoplastiques sous chargement cyclique. Ainsi, dans le cas de chargements mobiles et répétitifs, nous avons proposé l'algorithme de thermoélastoplasticité stationnaire. Cette procédure numérique permet de déterminer les états thermomécaniques générés par des sollicitations mobiles et répétées. Cette évaluation peut être effectuée cycle par cycle (Méthode Stationnaire Passage par Passage) ou directement à l'éventuel état stabilisé atteint généralement au bout d'un grand nombre de cycles (Méthode Stationnaire Directe). Actuellement, nos travaux sur les algorithmes de calcul des structures thermoélastoplastiques en contacts intermittents ou à faibles débattements avec comme application immédiate le fretting-fatigue. La méthode cyclique directe} que nous développons, qui s'apparente à une procédure de recherche de cycle limite à grand incrément de temps, s'appuie d'une part sur la méthode LATIN de Ladevèze et d'autre part sur lapropriété de périodicité de l'état limite. Le deuxième chapitre est consacré à l'évaluation de l'endommagement induit par les contacts mobiles. Dans le cas de la fatigue oligocyclique nous avons utilisé des critères disponibles dans la littérature (essentiellement la loi de Manson-Coffin). En ce qui concerne la fatigue polycyclique, c'est l'approche multi-échelle initiée par Dang Van qui est développée et appliquée. Les critères qui en sont issus être identifiés par des expériences simples à réaliser, d'une autre nature que des essais de roulement. Nous présentons une extension de l'approche macro-méso au cas de la fatigue polycyclique anisotherme et montrons qu'elle permet d'obtenir des formulations anisothermes cohérentes des critéres de fatigue polycyclique. Le troisième chapitre porte sur l'estimation de la dégradation des rails des chemins de fin. Ceux-ci constituent une source intéressante de problèmes des endommagements de contact : de la fatigue polycyclique s'initiant en profondeur (tache ovale) ou en surface (squat) en passantpar la fatigue oligocyclique (shelling) jusqu'à l'usure douce ou sévère, régulière ou ondulatoire, une large palette d'avaries est rencontrée. C'est donc naturellement que ce domaine illustre bien l'intérêt des outils de calcul et des modèles de fatigue que nous avions mis au point : l'utilisation de la démarche couplant calculstationnaire et critères de fatigue permet une étude systématique de tels problèmes en apportant des réponses quantitatives. Dans le dernier chapitre, on montre que le problème de fretting (endommagements induits par les contacts à faibles débattements peut être résolu par l'utilisation d'une démarche classique découplée consistant en un calcul de l'état stabilisé mécanique (par la méthode cyclique directe) et une prévision numérique d'endommagement qui est liée directement à la nature de l'état stabilisé. Ainsi la bonne concordance obtenue entre les simulations numériques et les essais dans le cas de la fissuration du fretting indique que celle-ci relève d'un problème de fatigue polycyclique classique.
• Modélisation thermomécanique des disques de frein par une approche eulérienne.
  
  • Nguyen-Tajan Thi Mac-Lan
  , 2002. L'objectif de cette thèse est de proposer une modélisation du comportement thermomécanique des disques de frein automobiles pour la prédiction de leur tenue en fatigue. L'approche adoptée associe notamment : des stratégies numériques nouvelles qui reposent sur des algorithmes eulériens adaptés aux problèmes de pièces sous chargement mobile et répété ; une modélisation pertinente du comportement du matériau constituant les disques ; des critères de tenue adaptés aux différents endommagements subis par le disque, qui prennent notamment en compte le caractère multiaxial et anisotherme des sollicitations. Dans un premier temps, est présentée une analyse des phénomènes physiques opérant dans un disque de frein en service (fortes amplitudes de température, gradients thermiques élevés, déformations plastiques, frottement disque /plaquette, etc. ) et de leurs couplages (couplages thermomécaniques, couplages contact disque/plaquettes et élévation de température, etc.). A partir de cette analyse, nous pouvons proposer une modélisation numérique du disque qui se fonde sur des hypothèses pertinentes. Cette modélisation utilise des algorithmes nouveaux qui permettent de diminuer considérablement les temps de calcul. Leur principe fondamental est de se placer dans un repère non plus lié à la structure mais au chargement et d'utiliser une approche eulérienne pour prendre en compte le mouvement de la structure dans ce repère. Ces algorithmes sont décrits dans cette thèse. Enfin, on propose une stratégie de dimensionnement qui repose sur une approche dite globale. Cette dernière consiste à proposer une modélisation des sollicitations, une modélisation du comportement du matériau, une stratégie de calcul de structure performant et des critères de tenue adaptés aux phénomènes d'endommagement observés. L'objectif est de donner un même « degré de précision » dans chacun des modèles pour assurer au global une filière robuste qui répond au besoin de l'ingénieur en charge de la conception des disques de frein.
• Une méthode de calcul de structures soumises a des chargements mobiles
  
  • Nguyen-Tajan Mac-Lan
  • Maitournam Habibou
  • Thomas Jean Jacques
  Revue Européenne des Éléments Finis, HERMÈS / LAVOISIER, 2002, 11, pp.247-261. Le travail présenté dans cet article s'inscrit dans une étude globale de la tenue en service d'un disque de frein automobile. Celui-ci subit des chargements thermomécaniques responsables de divers types d'endommagement : faïençage, fissures radiales sur les pistes de frottement, rupture du bol. La prédiction numérique de ces avaries exige d'abord la mise au point de méthodes numériques adaptées aux problèmes de structures soumises à des chargements thermomécaniques mobiles. Celles-ci consistent en des algorithmes eulériens qui sont présentés et décrits. (10.3166/reef.11.247-261)
  DOI : 10.3166/reef.11.247-261
• Nondestructive evaluation of material instability of a deep argillaceous rock
  
  • Luong M. P.
  • Eytard J. C.
  , 2002, pp.67-76. No abstract provided
• A new viscoelastic law used for thermomechanical analysis in the automotive industry
  
  • Verger Loick
  • Constantinescu Andrei
  • Pommier Benjamin
  , 2002.
• Approche expérimentale et modélisation micromécanique du comportement des matériaux
  
  • Equipe Micromécanique Xx
  , 2002.
• Crack shape sensitivity by the adjoint variable method using a boundary-only formula
  
  • Rus Carlborg G.
  • Bonnet Marc
  • Gallego Roberto
  , 2002. No abstract provided
• A micromechanical model for the TiAl lamellar grain
  
  • Gélébart Lionel
  • Bornert Michel
  • Crépin Jérôme
  • Bretheau Thierry
  , 2002.
• Computational basis for elastodynamic identification in a semi-infinite solid
  
  • Nintcheu Fata S.
  • Guzina B. B.
  • Bonnet Marc
  , 2002. No abstract provided
• Sensitivity analysis for shape perturbation of cavity or internal crack using BIE and adjoint variable approach
  
  • Bonnet Marc
  • Burczynski T.
  • Nowakowski M.
  International Journal of Solids and Structures, Elsevier, 2002, 39, pp.2365-2385. This paper deals with the application of the adjoint variable approach to sensitivity analysis of objective functions used for defect detection from knowledge of supplementary boundary data, in connection with the use of BIE/BEM formulations for the relevant forward problem. The main objective is to establish expressions for crack shape sensitivity, based on the adjoint variable approach, that are suitable for BEM implementation. In order to do so, it is useful to consider first the case of a cavity defect, for which such boundary-only sensitivity expressions are obtained for general initial geometry and shape perturbations. The analysis made in the cavity defect case is then seen to break down in the limiting case of a crack. However, a closer analysis reveals that sensitivity formulas suitable for BEM implementation can still be established. First, particular sensitivity formulas are obtained for special shape transformations (translation, rotation or expansion of the crack) for either two- or three-dimensional geometries which, except for the case of crack expansion together with dynamical governing equations, are made only of surface integrals (three-dimensional geometries) or line integrals (two-dimensional geometries). Next, arbitrary shape transformations are accommodated by using an additive decomposition of the transformation velocity over a tubular neighbourhood of the crack front, which leads to sensitivity formulas. This leads to sensitivity formulas involving integrals on the crack, the tubular neighbourhood and its boundary. Finally, the limiting case of the latter results when the tubular neighbourhood shrinks around the crack front is shown to yield a sensitivity formula involving the stress intensity factors of both the forward and the adjoint solutions. Classical path-independent integrals are recovered as special cases. The main exposition is done in connection with the scalar transient wave equation. The results are then extended to the linear time-domain elastodynamics framework. Linear static governing equations are contained as obvious special cases. Numerical results for crack shape sensitivity computation are presented for two-dimensional time-domain elastodynamics. (10.1016/S0020-7683(02)00131-2)
  DOI : 10.1016/S0020-7683(02)00131-2
• Parameter Identification for lined tunnels in a viscoplastic medium
  
  • Lecampion Brice
  • Constantinescu Andrei
  • Nguyen Minh Duc
  International Journal for Numerical and Analytical Methods in Geomechanics, Wiley, 2002, 26 (12), pp.1191-1211. This paper is dedicated to the identification of constitutive parameters of elasto‐viscoplastic constitutive law from measurements performed on deep underground cavities (typically tunnels). This inverse problem is solved by the minimization of a cost functional of least‐squares type. The exact gradient is computed by the direct differentiation method and the descent is done using the Levenberg–Marquardt algorithm. The method is presented for lined or unlined structures and is applied for an elastoviscoplastic constitutive law of the Perzyna class. Several identification problems are presented in one and two dimensions for different tunnel geometries. The used measurements have been obtained by a preliminary numerical simulation and perturbed with a white noise. The identified responses match the measurements. We also discuss the usage of the sensitivity analysis of the system, provided by the direct differentiation method, for the optimization of in situ monitoring. The sensitivity distribution in space and time assess the location of the measurements points as well as the time of observation needed for reliable identification. (10.1002/nag.241)
  DOI : 10.1002/nag.241