Les propriétés mécaniques et fonctionnelles des matériaux métalliques sont conditionnées par leurs microstructures, qui sont elles-mêmes héritées des traitements thermomécaniques subis. Etre capable de prévoir et simuler la microstructure et ses hétérogénéités lors des procédés de mise en forme complexes est récemment devenu crucial dans l'industrie métallurgique. C'est également un véritable challenge d'un point de vue numérique qui met en évidence l'importance des matériaux numériques dans les nouvelles méthodes de modélisation. Dans ce travail, nous nous intéressons à un modèle en champ complet récent basé sur la méthode level set (LS) dans un cadre éléments finis (EF) pour la modélisation des mécanismes de recristallisation.Les points forts de cette approche par rapport à l'état de l'art ont motivé le développement d'un logiciel appelé DIGIMU® par la société TRANSVALOR avec le soutien de grandes entreprises industrielles. Toutefois, le principal inconvénient de cette approche, commun aux autres méthodes en champ complet utilisant des maillages EF non structurés, reste son coût numérique important.Le principal objectif de ce travail a donc été d'améliorer considérablement le coût numérique de la formulation LS utilisée dans le contexte de maillages EF non structurés. De nouveaux développements génériques ont été réalisés pour améliorer l'efficacité globale du modèle. La formulation 2D LS existante, déjà utilisée pour modéliser la croissance de grains, la recristallisation statique et l'effet d'ancrage de Smith-Zener, a été étendue et améliorée afin de modéliser ces mécanismes en 3D pour des polycristaux à grand nombre de grains en des temps de calcul raisonnables. / Mechanical and functional properties of metallic materials are strongly related to their microstructures, which are themselves inherited from thermal and mechanical processing. Being able to accurately predict and simulate the microstructure and its heterogeneities after complex forming paths recently became crucial for the metallurgy industry. This is also a real challenge from a numerical point of view which highlights the importance of digital materials in new modeling techniques. In this work, we focus on a recent front-capturing full field model based on the level set (LS) method within a finite element (FE) framework to model recrystallization mechanisms.The strengths of this approach comparatively to the state of the art have motivated the development of a software package called DIGIMU® by the company TRANSVALOR with the support of major industrial companies. However, the main drawback of this approach, common with other front-capturing full field approaches working on unstructured FE meshes, is its important computational cost, especially in 3D.Main purpose of this work was finally to drastically improve the numerical cost of the considered LS-FE formulation in context of unstructured FE meshes. New generic numerical developments have been proposed to improve the global efficiency of the model. The existing 2D LS formulation, already used to model grain growth, static recrystallization and the Smith-Zener pinning effect, has been extended and improved in order to model these mechanisms in 3D for large-scale polycrystals with reasonable computational costs.

La croissance dendritique est calculée en utilisant un modèle champ de phase avec adaptation automatique anisotrope et non structurées d’un maillage éléments finis. Les inconnues sont la fonction champ de phase, une température adimensionnelle et une composition adimensionnelle, tel que proposé par [KAR1998] et [RAM2004]. Une interpolation linéaire d’éléments finis est utilisée pour les trois variables, après des techniques de stabilisation de discrétisation qui assurent la convergence vers une solution correcte non-oscillante. Afin d'effectuer des calculs quantitatifs de la croissance dendritique sur un grand domaine, deux ingrédients numériques supplémentaires sont nécessaires: un maillage adaptatif anisotrope et non structuré [COU2011], [COU2014] et un calcul parallèle [DIG2001], mis à disposition de la plateforme numérique utilisée (CimLib) basée sur des développements C++. L'adaptation du maillage se trouve à réduire considérablement le nombre de degrés de liberté. Les résultats des simulations en champ de phase pour les dendrites pour une solidification d'un matériau pur et d’un alliage binaire en deux et trois dimensions sont présentés et comparés à des travaux de référence. Une discussion sur les détails de l'algorithme et le temps CPU sont présentés et une comparaison avec un modèle macroscopique sont faite. / Dendritic growth is computed using a phase-field model with automatic adaptation of an anisotropic and unstructured finite element mesh. Unknowns are the phase-field function, a dimensionless temperature and a dimensionless composition, as proposed by [KAR1998] and [RAM2004]. Linear finite element interpolation is used for all variables, after discretization stabilization techniques that ensure convergence towards a correct non-oscillating solution. In order to perform quantitative computations of dendritic growth on a large domain, two additional numerical ingredients are necessary: automatic anisotropic unstructured adaptive meshing [COU2011], [COU2014] and parallel implementations [DIG2001], both made available with the numerical platform used (CimLib) based on C++ developments. Mesh adaptation is found to greatly reduce the number of degrees of freedom. Results of phase-field simulations for dendritic solidification of a pure material and a binary alloy in two and three dimensions are shown and compared with reference work. Discussion on algorithm details and the CPU time are outlined and a comparison with a macroscopic model are made.

Les bases de fonctions sont des outils indispensables de la géométrie numérique puisqu'ils permettent de représenter des fonctions comme des vecteurs, c'est à dire d'appliquer les outils de l'algèbre linéaire à l'analyse fonctionnelle. Dans cette thèse, nous présentons plusieurs constructions de bases de fonctions sur des surfaces pour la géométrie numérique. Nous commençons par présenter les bases de fonctions usuelles des éléments finis et du calcul extérieur discret, leur théorie et leurs limites. Nous étudions ensuite le Laplacien et sa discrétisation, ce qui nous permettra de construire une base de fonctions particulière~: les fonctions propres de l'opérateur de Laplace-Beltrami, ou harmoniques variétés. Celles-ci permettent de généraliser la transformée de Fourier et le filtrage spectral aux fonctions définies sur des surfaces. Nous présentons ensuite des applications de cette base de fonction à la géométrie numérique. En particulier, nous montrons qu'une fois calculée, cette base de fonction permet de filtrer la géométrie en temps interactif. Pour pouvoir définir des bases de fonctions de façon plus indépendante du maillage de la surface, nous nous intéressons ensuite aux paramétrisations globales, et en particulier aux champs de directions à symétries qui permettent de les définir. Ainsi, dans la dernière partie, nous étudions ces champs de directions à symétries, et en particulier leur géométrie et leur topologie. Nous donnons alors des outils pour les construire, les manipuler et les visualiser / Function bases are fundamental objects in geometry processing as they allow to represent functions as vectors, that is to apply tools from linear algebra to functional analysis. In this thesis, we present various constructions of useful functions bases for geometry processing. We start by presenting usual function bases, their theory and limits. We then study the Laplacian operator and its discretization, and use it to define a particular function basis: Laplacian eigenfunctions or Manifold harmonics. The Manifold Hamonics form a function basis that allows to generalize the Fourier transform and spectral filtering on a surface. We present some applications and extensions of this basis for geometry processing. To define function bases in a mesh-independant manner, we need to build a global parameterization, and especially the direction fields required to define them. Thus, in the last part of this thesis we study N-symmetry direction fields on surfaces, and in particular their geometry and topology. We then give tools to build, edit, control and visualize them

L’effet Biefeld-Brown, du nom de ses découvreurs dans les années 1920, désigne la force électrohydrodynamique (EHD) s’appliquant sur deux électrodes sous haute tension dans l’air. Si l’origine de cette force a pu faire l’objet de certaines spéculations, il est aujourd'hui admis qu’elle repose sur l’accélération par un fort champ électrique d’un volume d’air partiellement ionisé. Cet effet aussi appelé vent ionique intéresse diverses applications : contrôle actif d’écoulement, augmentation du transfert de chaleur par convection forcée, séchage de denrées alimentaires ou encore la propulsion. Cette thèse, présente une étude expérimentale, théorique et numérique du vent ionique dans une configuration modèle à deux électrodes parallèles. Le faible rendement du vent ionique l’a écarté des applications à la propulsion mais des expériences récentes menées en 2013 montrent qu’il permet d’atteindre un rapport poussée/puissance étonnement élevé. Nous montrons dans une première partie, à partir de mesures et de considérations aérodynamique générales que la poussée générée pourrait suffire à contrebalancer la force de traînée pour certains aéronefs ultra-légers. Ces mesures ont permis de quantifier la force EHD et sa dépendance avec la géométrie des électrodes. En outre, la meilleure configuration à deux collecteurs peut produire une poussée presque deux fois plus importante qu’une configuration avec un seul collecteur, à tension fixée. Ces premiers résultats ont été affinés dans un second temps par les mesures PIV qui ont permis la reconstruction de l’écoulement et du champ de force entre les électrodes. Les vitesses mesurées dépassent rarement 3 m/s, et la force volumique est de l’ordre de 10 N/m 3. L’origine physique de la configuration optimale à deux collecteurs a été éclaircie par la mise en évidence des structures de sillages et de leurs effets instationnaires. Par ailleurs, une analyse théorique générale de la force propulsive nous a permis de confirmer sa dépendance explicite avec le rapport courant sur mobilité ionique. Le courant étant directement lié à la physique de la décharge couronne, la seconde partie de la thèse s’est concentrée sur son analyse théorique et numérique. Une analyse asymptotique a ainsi permis de trouver une expression analytique du champ électrique critique et de la caractéristique courant-tension permettant de connaître l’influence de la densité du gaz et de sa composition sur le courant produit dans des électrodes concentriques. Cette approche asymptotique a été associée à une formulation de décomposition de domaine dans le cadre d’une discrétisation par éléments finis pour analyser des configurations plus générales. Une résolution itérative du système d’équations stationnaires non-linéaire couplées par méthode de Newton est proposée, testée et validée. Cette méthode peut être étendue à des géométries plus complexes, permettant ainsi d’obtenir une condition d’injection des charges prenant en compte la physique complexe de la décharge. / The Biefeld-Brown effect, named after its discoverers in 1920s, stands for the electrohydrodynamic (EHD) force applied on two high voltage electrodes in air. The origin of this force has been subject to controverse, but it is establised that it relies on the acceleration of ionized air by a strong electric field. Numerous applications are associated with ionic wind : active flow control, heat transfer enhancement, food drying or even propulsion. At the first glance, the low efficiency is unattractive for propulsion. However recent experiments highlighted a surprisingly high thrust/power ratio. This PhD research aimed to better understand the phenomena, through experiments and theoretical or numerical analysis. First, an experimental study was carried out to quantify the EHD force and its variations with the geometry of the electrodes. For instance, the best position of two collecting electrodes could produce nearly twice more thrust than the one collector configuration, for a given applied voltage. Considering the mass and the aerodynamic of some already existing very light aircrafts, it is shown that the produced thrust could at most balance the aerodynamic drag. This first results were enhanced by PIV measurements, which gave deeper insight into the flow and the force field between the electrodes. The air speed recorded was at most around 3 m/s, while the volumetric force of the order of 10 N/m3. The physical explaination of the optimal two collectors configuration relied partially on the wake flow structures and their unsteady effects. Then, a theoretical analysis of the propulsive force confirmed its explicit dependence on the current to mobility ratio. The discharge current being determined by the corona discharge physics, the second part of this work focuses on its theoretical and numerical analysis. An asymptotic approach of the corona discharge for concentric cylindrical electrodes led to an explicit expression of both the onset surface electric field and the current-voltage law as functions of the gas density, the effective ionization coefficient and the electrodes size. This asymptotic approach was reformulated in the frame of a domain decomposition method, implemented numerically with a finite elements discretization, in order to generalize the asymptotic approach. The iterative algorythm for the steady non-linear coupled system of equations is based on Newton method. This method provides a physically relevant boundary condition for the charge injection and can be applied to more complex geometries.

Ce travail tente de développer des outils numériques efficaces pour une approche plus rationnelle et moins empirique du dimensionnement à la ruine des ouvrages de génie civil. Contrairement aux approches traditionnelles reposant sur une combinaison de calculs élastiques, l'adoption de coefficients de sécurité et une vérification locale des sections critiques, la théorie du calcul à la rupture nous semble être un outil prometteur pour une évaluation plus rigoureuse de la sécurité des ouvrages. Dans cette thèse, nous proposons de mettre en œuvre numériquement les approches statique par l'intérieur et cinématique par l'extérieur du calcul à la rupture à l'aide d'éléments finis dédiés pour des structures de plaque en flexion et de coque en interaction membrane-flexion. Le problème d'optimisation correspondant est ensuite résolu à l'aide du développement, relativement récents, de solveurs de programmation conique particulièrement efficaces. Les outils développés sont également étendus au contexte de l'homogénéisation périodique en calcul à la rupture, qui constitue un moyen performant de traiter le cas des structures présentant une forte hétérogénéité de matériaux. Des procédures numériques sont spécifiquement développées afin de déterminer puis d'utiliser dans un calcul de structure des critères de résistance homogènes équivalents. Enfin, les potentialités de l'approche par le calcul à la rupture sont illustrées sur deux exemples complexes d'ingénierie : l'étude de la stabilité au feu de panneaux en béton armé de grande hauteur ainsi que le calcul de la marquise de la gare d'Austerlitz / This work aims at developping efficient numerical tools for a more rational and less empirical assessment of civil engineering structures yield design. As opposed to traditionnal methodologies relying on combinations of elastic computations, safety coefficients and local checking of critical members, the yield design theory seems to be a very promising tool for a more rigourous evaluation of structural safety. Lower bound static and upper bound kinematic approaches of the yield design theory are performed numerically using dedicated finite elements for plates in bending and shells in membrane-bending interaction. Corresponding optimization problems are then solved using very efficient conic programming solvers. The proposed tools are also extended to the framework of periodic homogenization in yield design, which enables to tackle the case of strong material heterogeneities. Numerical procedures are specifically tailored to compute equivalent homogeneous strength criteria and to use them, in a second step, in a computation at the structural level. Finally, the potentialities of the yield design approach are illustrated on two complex engineering problems : the stability assessment of high-rise reinforced concrete panels in fire conditions and the computation of the Paris-Austerlitz railway station canopy

Depuis plusieurs dizaines d’années, les cristaux photoniques et phononiques font l’objet d’études poussées notamment enoptique, électromagnétisme et en acoustique. Ces Métamatériaux, constitués de diffuseurs périodiques, ont despropriétés impossibles à observer pour des matériaux usuels et peuvent par exemple, courber les rayons où interdire latransmission des ondes sonores sur certaines gammes de fréquences (bandes interdites). En agissant sur lescaractéristiques géométriques du cristal il est possible de combiner les pertes en transmission liées à la période, à deseffets de résonances plus basses fréquences liés aux diffuseurs (rigide, résonnant...) et obtenir des coefficients de transmission quasi nuls, où d’absorption quasi totale sur de larges bandes de fréquences. Deux métamatériaux sont étudiés, visant à trouver des solutions alternatives à des problématiques rencontrées en acoustique et utilisant un réseau périodique d’inclusions résonantes. Le premier est un cristal sonique utilisé comme barrière acoustique et créé à l’aide de cannes de bambou percées comportant des pertes en transmission basses fréquences. Le second est un panneau de matériau poreux enfermant des inclusions résonantes et offrant une absorption acoustique quasi totale pour des longueurs d’ondes jusqu’à 10 fois supérieures à l’épaisseur du matériau. Les comportements de ces deux dispositifs ont été étudiés théoriquement, expérimentalement et numériquement via plusieurs méthodes qui ont permis de mettre en évidence leurs excellentes performances pour des applications acoustiques dans l’audible. / Since several decades, photonic and phononic crystals are the center of numerous studies and in particular in the optics,electromagnetism and acoustics fields. These metamaterials, created by a periodic array of inclusions, have propertiesimpossible to obtain with usual materials. They can, for example, bend the waves or stop the waves for some frequencyranges (band gap). By changing the characteristic of the unit cell, it is possible to combine transmission losses linked to theperiodicity, with low frequency resonances linked to the type of scatterer (rigid, resonator...) and obtain very low transmissioncoefficient or very high absorption coefficient on very large frequency ranges depending on the device. Two metamaterialsdevices are studied to find alternative solutions, for acoustics problems, by using periodic array of scatterers. The first deviceis a sonic crystal used has an noise barrier and built with drilled bamboo rods, that have low frequency transmission losses(around 300 Hz and around 2000 Hz). The second device is a periodic array of resonant inclusions embedded in a porousplate that can absorb almost all the waves for a wide frequency range that correspond to wavelength up to 10 times bigger than de thickness of the plate. The behavior of this two devices are studied theoretically, experimentally and numerically by using several methods (Plane Waves Expansion, Multiple Scattering Theory for the first device and finite element method for the second). All this methods allow to bring out the very good performances of this metamaterials devices in audiblefrequency range.

Dans ce travail, nous nous sommes intéressés à la mise en place de schémas semi-implicites pour l’amélioration des simulations numériques du déplacement d’un globule rouge dans le sang. Nous considérons la méthode levelset, où l’interface fluide-structure est représentée comme la ligne de niveau 0 d’une fonction auxiliaire et le couplage est effectué en ajoutant un terme source dans l’équation fluide.Le principe de ces schémas semi-implicites est de prédire la position et la forme de la structure par une équation de la chaleur et d’utiliser cette prédiction pour obtenir un terme de force dans l’équation fluide plus précis. Ce type de schémas semi-implicites a d’abord été mis en place dans le cadre d’un système diphasique ou d’une membrane élastique immergée afin d’utiliser un plus grand pas de temps que pour un couplage explicite. Cela a permis d’améliorer les conditions sur le pas de temps et ainsi augmenter l’efficacité globale de l’algorithme complet par rapport à un schéma explicite classique.Pour étendre ce raisonnement au cas d’un globule rouge, nous proposons un algorithme pour simuler le flot de Willmore en dimension 2 et 3. Notre méthode s’inspire des méthodes de mouvements d’interface générés par diffusion et nous arrivons à obtenir un flot non linéaire d’ordre 4 uniquement avec des résolutions d’équations de la chaleur. Pour assurer la conservation du volume et de l’aire d’un globule rouge, nous proposons ensuite une méthode de correction qui déplace légèrement l’interface afin de recoller aux contraintes.La combinaison des deux étapes précédentes décrit le comportement d’un globule rouge laissé au repos. Nous validons cette méthode en obtenant une formed’équilibre d’un globule rouge. Nous proposons enfin un schéma semi-implicite dans le cas d’un globule rouge qui ouvre la voie vers l’utilisation de cette méthode comme prédicteur de l’algorithme de couplage complet. / In this work, we propose new semi-implicit schemes to improve the numerical simulations of the motion of an immersed red blood cell. We consider the levelset method where the interface is described as the 0 isoline of an auxiliary function and the fluid-structure coupling is done by adding a source term in the fluid equation.The idea of these semi-implicit scheme is to predict the position and the shape of the structure through a heat equation and to use this prediction to improve the accuracy of the source term in the fluid equation. This type of semi-implicit scheme has firstly been implemented in the case of a multiphase flow and a immersed elastic membrane and has shown better temporal stability than an explicit scheme, resulting in an improved global efficiency.In order extend this method to the case of a red blood cell, we propose an algorithm to compute the Willmore flow in dimenson 2 and 3. In the spirit of the diffusion generated motion methods, our method simulate a non linear four order flow by only solving heat equations. To ensure the conservation of the volume and area of the the vesicle, we add to the method a correction step that slightly moves the interface so that we recover the constraints.Combnation of these two steps allows to compute the behavior of a red blood cell left at rest. We validate this method by obtaining the convergence to an equilibrium shape in both 2D and 3D. Finaly we introduce a semi-implicit scheme in the case of a red blood cell that shows how we can use this method as a prediction in the complete coupling model.

Pour être commercialisés et installés à bord des navires, les panneaux structuraux tels que les cloisons et les pontsdoivent passer avec succès un essai normalisé de résistance au feu de type ISO 834. De tels essais sont longs et coûteux,les constructeurs de panneaux souhaitent par conséquent maximiser les chances de succès lors du passage de leur produitau test en particulier pour les designs alternatifs que constituent les sandwichs composite.La simulation par éléments finis est un outil pour modéliser le comportement thermomécanique de telles structures.Les codes industriels comme SAMCEF qui a été utilisé dans ce travail sont capables de réaliser des analyses thermiquesavec dégradation et des analyses mécaniques mettant en jeu la gestion du contact, la dégradation des propriétésmécaniques et la perte des structures par ruine.L’enjeu de cette étude est d’abord identifier les spécificités de telles structures particulières par leur taille, leurconception, de déterminer quelles sont les données thermiques et mécaniques nécessaires pour alimenter le modèlenumérique et le cas échéant de les construire à partir des résultats de la réaction au feu des matériaux. Des modèlesnumériques que l’on souhaite robustes et utilisables dans un contexte industriel, sont développés pour déterminer lecomportement thermomécanique de tels panneaux. Ils prennent en compte l’évolution des propriétés thermiques etmécaniques des matériaux en cours de dégradation. Ces modèles doivent permettre par la suite l’estimation desperformances de nouveaux designs lors d’un essai de certification ISO 834. / Structural panels used in naval industry such as bulkheads and decks must succeed in standard certification testssuch as ISO 834 to be commercialized and settled on board. As these tests are long and expensive, panel manufacturerswish to maximize chances of success for their panels when submitted to certification tests especially when it comes toalternative designs such as composite sandwich panels.Finite elements analyses are used to model thermo-mechanical behavior. Industrial software such as SAMCEF,which was used to conduct this work, are able to solve thermal analyses with degradation and mechanical analyses involvingcontact conditions, degradation of mechanical properties and loss of structures due to failure.The objective in this study is to first identify characteristics of these structures. They are special due to their size andtheir manufacturing. This study aims also to determine thermal and mechanical data required for numerical modeling.When necessary some data can be computed from results coming from the results of the materials’ reaction to fire.Numerical models are developed to determine thermo-mechanical behavior and are designed to be robust and used inindustrial context. They include the evolution of thermal and mechanical properties during the degradation process. Thesemodels must enable to estimate the performances of innovative designs during an ISO 834 certification test.

Le système ostéo-articulaire a un rôle de soutien essentiel pour l’ensemble du corps humain. Il supporte les forces de gravité et les efforts produits par les activités quotidiennes. Ainsi l’os optimise et adapte sa masse et sa géométrie à travers le processus de remodelage osseux. Cette régulation spatio-temporelle peut subir des déséquilibres métaboliques comme l’ostéoporose qui conduisent à la survenue de fracture de l’extrémité supérieure du fémur lors d’une chute latérale ou de traumatismes divers. En effet l’ostéoporose se traduit par une diminution au niveau de la densité osseuse et des détériorations au niveau de la microarchitecture du tissu osseux augmentant ainsi le risque de fractures.Dans le cadre de cette thèse une modélisation mécanique du comportement de l’os ainsi qu’une modélisation biologique des activités cellulaires ont été proposées dans un premier temps. Dans un second temps, la pharmacocinétique de quatre médicaments a été modélisé soient l’Alendronate, le Denosumab,le Romosozumab et l’Odanacatib. Finalement, le couplage de ces modèles a permis d’avoir un modèle mécano-biologique couplé aux effets des médicaments contre l’ostéoporose qui permet de prédire l’évolution de la densité minérale osseuse et celle de l’endommagement par fatigue permettant ainsi d’analyser l’évolution de la qualité osseuse. Ce modèle a été implémenté au code de calcul par éléments finis ABAQUS/standard à travers sa routine utilisateur UMAT. Le modèle a été appliqué pour simuler différents scénarii de remodelage sur des fémurs humains (2D et 3D). Différents facteurs ont été analysés tels que l’amplitude des activités physiques, les doses de médicaments injectées, la durée des traitements, etc. les résultats obtenus sont cohérents (qualitativement) avec les études cliniques existantes. En conclusion, le modèle mécano-biologique couplé aux effets des médicaments proposés contribue à l’analyse fine du comportement de l’os et l’application des algorithmes a permis d’effectuer des essais virtuels permettant d’analyser les effets combinés de nombreux facteurs pluridisciplinaires caractérisant la qualité osseuse. / The osteo-articular system plays the role of crucial support for the whole human body. It supports the gravity forces as well as the efforts generated by daily activities. Thus the bone optimizes and adapts its mass and its geometry through the process of bone remodeling. This spatio temporal regulation can undergo metabolic imbalances such as osteoporosis which lead to the occurrence of the upper end of the femur fracture during a fall side or various traumas. In fact, osteoporosis is reflected into a decrease in bone density and level of damage at the level of the microarchitecture of bone tissue increasing therefore the risk of fractures. In the framework of this thesis, a mechanical modelling of the behaviour of the os as well as a biological modelling of cellular activities were firstly proposed. Secondly, the pharmacokinetics of four drugs were modeled which are Alendronate, Denosumab, the Romosozumaband the Odanacatib. Finally, the coupling of these models allowed us to have a mechanic-biological model coupled with the effects of drugs against osteoporosis that can predict the evolution of bone mineral density and of damage by fatigue allowing to analyze the evolution of the bone quality. This model has been implemented by finite elements ABAQUS/standard through its user routine UMAT.The model has been applied to simulate different scenarios of remodeling on human femurs (2D and3D). Different factors were analysed such as the range of physical activities, the doses of injected drugs,the duration of treatment, etc. Obtained results are consistent (qualitatively) with existing clinical studies.As a conclusion, the mechanic-biological model coupled with the effects of the proposed drugs contributes to the fine analysis of the bone behaviour and the application of algorithms allowed to conduct virtual tests in order to analyze the combined effects of many multidisciplinary factors characterizing the bone quality.

Les dépôts effectués sur substrats de silicium sont très courants notamment dans le domaine de la microélectronique. Les propriétés physiques recherchées pour ce type de structures dépendent fortement de celles de la couche. Il apparaît donc essentiel de connaître les paramètres élastiques ainsi que l’épaisseur des films considérés. De plus, la détection de certains défauts concernant la couche est souvent recherchée. L’objectif de ce travail a été de contribuer à la caractérisation de structures du type couche sur substrat. Pour cela, les ultrasons-lasers présentent de nombreux avantages puisqu’ils autorisent entre autres leur contrôle non destructif sans contact. Les ondes acoustiques de surface dans une gamme de fréquence s’étendant jusqu’à 45 MHz ont été utilisées. Nous avons développé différents modèles analytiques et les résultats expérimentaux ont aussi été comparés à certaines simulations par éléments finis. Plus particulièrement, nous avons montré qu’il était possible d’obtenir l’ensemble des paramètres élastiques du substrat et de la couche ainsi que l’épaisseur de cette dernière. Par ailleurs, nous nous sommes aussi intéressés à la détection de certains défauts en régime impulsionnel mais aussi quasi-monochromatique. Des résultats originaux concernant l’effet d’une absence de couche de forme déterminée sur le premier mode de Rayleigh ou bien encore de problèmes d’adhésion ont été présentés. Sur ce dernier point, une méthode innovante permettant de distinguer un fort niveau d’adhésion d’un faible a aussi été introduite. / Thin films deposited on silicon substrates are very common especially in microelectronic applications. The physical properties expected for these types of structures depend on the properties of the layer. Therefore, it is essential to know elastic parameters and thickness of the films considered. Furthermore, some layer defects detection is often required. The aim of this work was to contribute to the characterization of structures composed of a single layer deposited on a substrate. For this, laser ultrasonics offers many advantages since it is nondestructive and non-contact method. Surface acoustic waves excited in a frequency range up to 45 MHz are used. Analytical models have been developed and experimental results have also been compared with some finite element simulations. More specifically, we have shown that it was possible to obtain the thickness of the layer and all elastic parameters of both substrate and layer. In addition, we have been interested in the detection of various defects with a broadband and quasi-monochromatic excitation. Original results concerning the effect of a lack of layer with a specific geometry on the first Rayleigh mode or even adhesion problems were presented. Concerning this last point, an innovative method to tell difference between high and low adhesion was introduced.