Le nouvel élément de poutre est une évolution d'un élément de Timoshenko poutre avec un nœud supplémentaire situé à mi-longueur. Ce nœud supplémentaire permet l'introduction de trois composantes supplémentaires de contrainte afin que la loi constitutionnelle 3D complète puisse être utilisée directement. L'élément proposé a été introduit dans un code d'éléments finis dans Matlab et une série d'exemples de linéaires/petites contraintes ont été réalisées et les résultats sont systématiquement comparés avec les valeurs correspondantes des simulations ABAQUS/Standard 3D. Ensuite, la deuxième étape consiste à introduire le comportement orthotrope et à effectuer la validation de déplacements larges / petites contraintes basés sur la formulation Lagrangienne mise à jour. Une série d'analyses numériques est réalisée qui montre que l'élément 3D amélioré fournit une excellente performance numérique. En effet, l'objectif final est d'utiliser les nouveaux éléments de poutre 3D pour modéliser des fils dans une préforme composite textile. A cet effet, la troisième étape consiste à introduire un comportement de contact et à effectuer la validation pour un nouveau contact entre 3D poutres à section rectangulaire. La formulation de contact est dérivée sur la base de formulation de pénalité et de formulation Lagrangian mise à jour utilisant des fonctions de forme physique avec l'effet de cisaillement inclus. Un algorithme de recherche de contact efficace, qui est nécessaire pour déterminer un ensemble actif pour le traitement de contribution de contact, est élaboré. Et une linéarisation constante de la contribution de contact est dérivée et exprimée sous forme de matrice appropriée, qui est facile à utiliser dans l'approximation FEM. Enfin, on présente quelques exemples numériques qui ne sont que des analyses qualitatives du contact et de la vérification de l'exactitude et de l'efficacité de l'élément de 3D poutre proposé. / The new beam element is an evolution of a two nodes Timoshenko beam element with an extra node located at mid-length. That extra node allows the introduction of three extra strain components so that full 3D stress/strain constitutive relations can be used directly. The second step is to introduce the orthotropic behavior and carry out validation for large displacements/small strains based on Updated Lagrangian Formulation. A series of numerical analyses are carried out which shows that the enhanced 3D element provides an excellent numerical performance. Indeed, the final goal is to use the new 3D beam elements to model yarns in a textile composite preform. For this purpose, the third step is introducing contact behavior and carrying out validation for new 3D beam to beam contact with rectangular cross section. The contact formulation is derived on the basis of Penalty Formulation and Updated Lagrangian formulation using physical shape functions with shear effect included. An effective contact search algorithm is elaborated. And a consistent linearization of contact contribution is derived and expressed in suitable matrix form, which is easy to use in FEM approximation. Finally, some numerical examples are presented which are only qualitative analysis of contact and checking the correctness and the effectiveness of the proposed 3D beam element.

L’allègement des structures imposé par les réductions de coût se traduit par des structures de plus en plus souples qui les rendent de plus en plus sensibles aux vibrations. Le contrôle des vibrations devient donc un enjeu majeur dans de nombreuses applications industrielles et les limites des matériaux imposent maintenant un recours au contrôle actif de plus en plus fréquent. L’évolution des structures au cours du temps (viellisement, conditions aux limites, architecture, …) pose le problème de la robustesse du contrôle. Par ailleurs, l’actionnement de plus en plus présent dans le domaine mécanique constitue à la fois une source supplémentaire de vibrations, mais aussi de contrôle et d’évolution d’architecture des structures. La thèse s’intéresse au contrôle actif autoadaptatif des vibrations permettant de maintenir automatiquement la performance et la stabilité des structures évolutives. Il s’agit donc de s’affranchir de la connaissance des causes et des informations sur les évolutions. La méthode proposée s’appuie sur un développement modal permettant de limiter le nombre de composants de contrôle et de cibler les modes à contrôler en limitant l’énergie de contrôle. Ainsi, il est nécessaire de reconstruire les caractéristiques du modèle modal indispensables pour réactualiser le contrôle en figeant seulement une structure de modèle. S’affranchissant à la fois des causes d’évolution de la structure et utilisant seulement une structure de modèle, la méthode est généralisable à toute application en mécanique des structures. La méthode proposée, basée sur l’utilisation d’un identificateur exploitant à la fois excitation et réponse de la structure, prend en compte les limites imposées par le contrôleur. Le modèle constitue le lien qui doit être établi entre identificateur et contrôle pour permettre la réactualisation. Par ailleurs, un compromis entre l’objectif d’atténuation des vibrations et les performances de l’identification est alors nécessaire du fait du couplage identification/contrôle apparaissant dans la boucle fermée. Ce compromis est également conditionné par le matériel utilisé. La méthode proposée est exploitée sur une structure discrète mettant en évidence une inversion de formes modales au cours de son évolution qui déstabilise un contrôle figé. Le choix opéré pour répondre aux différents compromis cités ci dessus a conduit à l’utilisation d’un contrôleur classique (LQG) et un identificateur basé sur la méthode des sous-espaces (N4SID). Cette application sur une structure simple a permis de caractériser un certain nombre de limites physiques : la bande passante, densité modale, vitesse d’évolution, Le contrôle modal autoadaptatif proposé s’avère robuste en performance et efficace lorsque la réactualisation est systématique. Une variante conditionnelle, toujours basée sur l’analyse de la réponse de la structure, est enfin proposée pour optimiser le processus de réactualisation afin de suivre plus efficacement les évolutions. / The lightness of structure due to the reduction of cost results in some structures which are more and more flexible. This flexibility makes these structures more sensitive to vibrations. The vibration control becomes an important issue in lots of industrial applications, and now the limitation of materials imposes a requirement of active control more and more frequently.The change of time-varying structure(ageing effect, boundaries conditions, architecture of structure etc)brings the robust problem of control.Further more,the action of device which emerges more and more frequently in mechanical fields introduces not only an additional cause of vibrations,but also a source of control and a source for changing the architecture of structures.The thesis focuses on self-adaptive active control of vibration which permits to keep up automatically the performance and stability of the time-varying structures.So it needs to overcome the knowing about cause and information on the changes.The proposed method relies on a development of modal technology which permits to limit the amount of component in control system and to target on the modes which need to be controlled.So the energy of control is limited. Further more,it needs to reconstruct the characteristics of modal model which are indispensable for updating the control.In this case, only the structure of model is fixed.Overcoming the knowing about cause of change in the structure and using only the structure of model, this method can be generalized for all applications in mechanical structures.The proposed method is based on the utilization of an identifier which uses both the excitation and response of the structure.And this method considers the limitations induced by the controller.The model forms le link which should be established between the identifier and the controller for allowing the updating. Further more, a compromise between the objective of reducing vibrations and the performance of identification is necessary due to the coupling effect of identification/control which appears in the closed-loop. This compromise is also conditioned by the used equipments.The proposed method is carried out on a discrete time-varying structure for showing an inversion of mode shape during its change. This inversion of mode shape destabilises a fixed control system. The operated choices for responding the different previous quoted compromise lead to a classic controller (LQG) and an identifier based on the subspace method (N4SID).This application on a simple structure permitted to characterise some physical limitation: the bandwidth, the modal density and the velocity of change…The proposed self-adaptive modal control is proved to be robust in terms of performance and be efficient when the updating is systematical. Always based on the analysis of the response of the structure, a conditional variant is finally proposed for optimizing the process of updating in order to follow the change more efficiently.

Les stratifiés composites sont de plus en plus utilisés dans les pièces de structures des aéronefs ce qui fait émerger de nouvelles problématiques comme celle des Impacts de Faible Energie (IFE). En effet, bien qu’ils possèdent des propriétés rapportées à leur masse très intéressantes ces matériaux peuvent être vulnérables aux petits chocs. Or, compte tenu des nombreux paramètres influents lors d’un tel impact (énergie, vitesse, stratification...), les essais actuellement majoritairement privilégiés à l’échelle industrielle sont long et coûteux. Ainsi, l’apport de la simulation numérique pourrait être d’une grande aide pour les constructeurs. La pratique du « virtual testing », en particulier, permettrait d’aller dans cette direction ce qui aurait pour effet de rationaliser les campagnes d’essais et les coûts financier qui en découlent. Cependant, elle peine à être mise en place ici car le temps CPU nécessaire pour la simulation fine des ndommagements induits par les IFE est trop important avec les méthodes actuelles. Partant de ce constat, ce travail a consisté à tirer avantageusement partie de la localisation spatiale et temporelle des délaminages, fissurations matricielles et ruptures de fibres qui peuvent apparaître pendant l’impact pour diminuer le coût de calcul. Ainsi une méthode multiéchelle en espace et en temps a été mise en place. Elle consiste à découper la structure impactée en deux zones. L’une est située autour du point d’impact, elle contient l’ensemble des non-régularités du problème (contact, loi adoucissante, modèle de zone cohésive). Elle est traitée avec le code de dynamique explicite Europlexus. L’autre correspond à la partie complémentaire. Le problème mécanique y est beaucoup plus régulier et il est traité avec le code de dynamique implicite Zset/Zébulon. Un couplage peu intrusif basé sur la méthode GC est donc réalisé entre ces deux codes. Il permet d’utiliser une modélisation adaptée dans chacune des deux régions ce qui permet en particulier d’utiliser des pas de temps différents. Un rapport supérieur à 1000 peut ainsi être obtenu entre celui du code explicite fixé par la condition de stabilité et celui utilisé dans la partie complémentaire. Un gain de temps CPU significatif confirmé par la simulation d’un impact réalisé sur un panneau composite raidi est ainsi obtenu. Il est également montré que la répartition implicite/explicite peut évoluer au cours du calcul. Pour cela un mécanisme de bascule a été mis en place. Il permet ainsi de faire transiter la résolution d’une partie de la structure initialement traitée dans le code Zebulon dans Europlexus. Un gain de temps supplémentaire est alors obtenu grâce à cette méthode sur le même cas d’application. / The composite laminates are increasingly used in aircraft structural parts which lead to new issues such as the Low Energy Impacts (LEI). Indeed, although they have well mechanical properties relative to their mass, small shocks may be very harmfull for laminates. Controlling such situations is essential for manufacturers that why lot of testing campaigns are currently performed. Yet, they are time consuming and expensive considering the many influential parameters (energy, speed, layup...). Numerical simulations of this phenomenon by practicing the so called “virtual testing” process could be really helpfull to rationalize testing campaigns in order to save money. Yet, this practice remain currently hard to do at the industrial scale due to the excessive CPU time required for fine simulation of damages induced by the LEI. Based on this observation, this work has consisted in taking advantage of the spatial and temporal location of delamination, matrix cracking and fiber breakage that can occur during impact in order to reduce the computational cost. Thus, a space and time multiscale method has been put in place. The impacted structure is split into two areas. One is located around the impacted point, it contains all the non-regularities of the problem (contact, softening law, cohesive zone model). This domain is treated with the explicit dynamics code Europlexus. The other one corresponds to the complementary part. The mechanical problem is much more regular and it is treated with the implicit dynamics code Zset / Zebulon. A low intrusive coupling based on the GC method is carried out between these two codes. It allows to use an adapted model in both regions different time step are in particular used. A time step ratio upper to 1000 can be reach between the one of the explicit code set by the stability condition and the one used in the complementary part. As a results, significant CPU time is saved. This is confirmed by the simulation of a stiffened composite panel impacted. It is also shown that the implicit / explicit allocation can change over the calculation. To do that, a switch mechanism has been established. It thus makes it possible to transit the resolution of a portion of the structure initially solved in the code Zebulon to Europlexus. As a results, further gain is obtained.

La thèse présente une méthodologie pour construire des éléments finis de type « solide-coque » avec intégration réduite en vue des applications à la simulation de la mise en forme tel que l’emboutissage des tôles où ces éléments finis doivent présenter de bonnes aptitudes à modéliser la flexion mais également les situations de laminage de la tôle. A partir des éléments volumiques à 8 nœuds et 3 degrés de liberté par nœud (les 3 composantes du déplacement), un neuvième nœud est rajouté au centre de l’élément. Ce neuvième nœud n’est pourvu que d’un seul degré de liberté, le déplacement le long de la direction de l’épaisseur. Cette direction privilégiée a un nombre de points d’intégration supérieur ou égal à 3 mais l’intégration est réduite au centre de l’élément diminuant très sensiblement les temps CPU par rapport à une intégration complète. Un soin particulier a été pris pour contrôler tous les modes à énergie nulle dus à l’intégration réduite. Ce nœud supplémentaire permet une distribution linéaire de la déformation normale. Avec les lois de comportement complètement 3D ces nouveaux éléments solide-coque donnent des résultats similaires en flexion à ceux obtenus avec des éléments coques et état plan de contrainte. Le neuvième nœud joue le rôle d’un paramètre supplémentaire pour l’interpolation quadratique du déplacement dans la direction de l’épaisseur. Ce degré de liberté a une signification physique et un effort équivalent à une pression normale peut être prescrit. Dans les situations de pression normale et dans le cas du contact, la contrainte normale obtenue est physique ce qui n’est pas le cas de nombreux éléments solide-coque de la littérature. Le pincement ou le laminage des tôles est correctement modélisé. Pour valider ces éléments, un module d’emboutissage en U avec passage et laminage de la bande de tôle sur des rouleaux a été construit au laboratoire. La comparaison entre les efforts d’emboutissage calculés et mesurés est très bonne ainsi que la géométrie des bandes de tôle obtenue après retour élastique. / This thesis presents a methodology for developing under-integrated “solid-shell” finite elements for sheet forming simulations like deep drawing where these elements must offer a bending capability and sheet thinning conditions as well. Starting from 8-node elements endowing three degrees of freedom per node (three displacement components), a ninth node is added at the centre of the element. This extra node has just one degree of freedom: a displacement along the ‘thickness’ direction. Several integration points are distributed along this privileged direction (5 points, generally) but the in-plane reduced integration at the centre of the element decreases CPU costs compared to a full integration. A special care has been taken to control all zero-energy modes due to the reduced integration. This additional node allows a linear distribution of the normal strain. With fully-3D constitutive laws, these new solid-shell elements give similar bending results as those obtained with shell elements and a plane stress state hypothesis. This ninth node acts as an additional parameter for the quadratic interpolation of the displacement in the ‘thickness’ direction. The corresponding degree of freedom has a physical meaning and a force, equivalent to a normal pressure for instance, may be prescribed. In situations of a normal pressure and in the case of contact, the obtained normal stress is physically defined, which is not the case for many solid-shell elements found in the literature. The pinching (or the thinning) of sheets is properly modelled. To validate these elements, an apparatus for U-drawing tests with ironing or thinning on strip sheets has been built in the laboratory. The comparison between numerical and experimental punch force during sheet forming is pretty good as well as the geometry of blank after springback.

La compréhension du comportement de structures jusqu'à leur ruine est nécessaire pour concevoir au mieux ces structures. Selon le matériau et les sollicitations considérées, les mécanismes physiques à l'origine de la rupture changent. Nous nous intéresserons à des matériaux homogènes pour lesquels la ruine passe par le développement de fissures autour desquelles les non-linéarités de comportement n'ont pas un rôle dominant. Ces conditions sont réunies pour les matériaux fragiles pour lesquels la source principale de dissipation est la génération non réversible d'une surface libre, et pour certaines fissures de fatigue. Sur un cycle de chargement, il existe de nombreuses applications pour lesquelles les non-linéarités restent confinées. La théorie de la mécanique linéaire élastique de la rupture est alors un modèle pertinent pour approcher le comportement de la structure. Sous ces hypothèses, le front de la fissure introduit une singularité. L'étude asymptotique de cette singularité dans des situations plane et anti-plane permet de définir les séries de Williams. La singularité est alors d'ordre un demi et elle est quantifiée par les facteurs d'intensité des contraintes (FIC) pour chacun des trois modes de sollicitations. En 3D, la fissure peut avoir une géométrie complexe, et aucune expression générale de la singularité n'existe. Dans cette thèse, les séries de Williams en déplacements sont utilisées et régularisées le long du front au sens des éléments finis. À partir de cette définition 3D des séries asymptotiques en pointe de fissure, une méthode d'extraction directe des FIC (DEK-FEM) est étendue au cas 3D. Le domaine est décomposé en deux domaines, raccordés en moyenne sur l'interface. Au voisinage du front, les champs mécaniques sont approchés par une troncature des champs asymptotiques. La singularité est donc traitée avec des champs adaptés, et les degrés de liberté associés sont directement les coefficients asymptotiques. Parmi ces coefficients asymptotiques, on retrouve les FIC et les T-stresses. Pour des raisons d'efficacité numérique et pour pouvoir relier l'échelle de la structure à l'échelle de la fissure, cette méthode est intégrée dans un contexte multigrilles localisées X-FEM. Ainsi nous montrons que cette approche permet une bonne évaluation des évolutions des FIC et du T-stress. Cette méthode est développée en parallèle d'une stratégie de post-traitement expérimental (mesure de champs de déplacements par corrélation d'images) basée sur les mêmes séries asymptotiques. Les images tridimensionnels d'un essai de fatigue in situ sont obtenues par micro-tomographie à rayons X et reconstruction. La corrélation et la régularisation basées sur les séries asymptotiques permettent d'obtenir la géométrie de la fissure et les FIC pour pouvoir identifier des lois de propagation de fissures 3D en fatigue. L'efficacité de cette méthode en parallèle d'une simulation DEK-FEM est illustrée en 2D. / It is necessary to understand the behavior of structures up to their failure to enhance their design. The mechanisms and phenomena undergoing failure vary according to the considered material and boundary conditions. We consider homogeneous materials for which cracks propagate in a context where behavior nonlinearities are not dominants. These conditions are matched for brittle and quasi-brittle materials and for some fatigue cracks. For the former, the main source of dissipation is the crack propagation which can be seen as the generation of a new free-surface. For the later, there is many applications where, in one loading cycle, the nonlinearities remains confined around the crack tip. The linear elastic fracture mechanics theory is then a pertinent model to approximate the structure behavior. Under such hypotheses, a singularity appears in the crack tip vicinity. The Williams' series expansion is computed from the asymptotic study of plane and anti-plane states. The stress is singular at the crack tip and the order of this singularity is one out of two. The singularity amplitude is quantified by the stress intensity factors (SIF), one for each of the three loading modes. In 3D, the crack shape is potentially complex (front curvature and non-planar crack), and no general asymptotic series expansion exists. In this PhD thesis, the 2D Williams' series in displacements are used and regularized with a finite element evolution along the front. From this 3D definition of the asymptotic fields in the crack tip vicinity, a numerical method for direct estimation of the SIF (DEK-FEM) is extended to 3D. This method is based on domain decomposition, the two domains are bounded in a weak sense on their interface. In the crack tip vicinity, the mechanical fields are approximated by a truncation of the asymptotic series expansion. Therefore, appropriate fields are used to deal with the singularity, and the associated degrees of freedom are directly the asymptotic coefficients. Among these coefficients are the SIF and the T-stresses. To bridge the scales between the structure and the crack front singularity and to increase the numerical efficiency, this method is embedded in a localized X-FEM multigrids approach. The proposed method is shown to provide an accurate evaluation of the SIF and T-stresses evolution. This approach has been developed in combination of an experimental post-processing method (full field displacement measurement through image correlation) based on the same asymptotic series expansion. The 3D images can be obtained for in situ fatigue experiments by X-ray microtomography and reconstruction. The crack geometry and the SIF are then provided by image correlation and regularization based on Williams series expansion. These data can be used for identifying a 3D fatigue crack growth law. The efficiency of the method is illustrated in 2D.

Ces travaux de thèse sont le fruit de la collaboration entre la société Hispano-Suiza et le LaMCoS de l’INSA de Lyon. Dans le cadre du développement de nouveaux systèmes de propulsion, l’implantation d’un train planétaire / épicycloïdal entre la turbine et l’hélice semble être une voie intéressante pour atteindre les performances souhaitées en terme de rendement. L’augmentation des puissances transmises et la réduction des masses embarquées dans les applications aéronautiques tendent à rendre les composants de plus en plus déformables. Lors de ces travaux de recherche, un modèle dynamique de trains planétaires a été développé, incorporant les effets des erreurs de montage, des écarts de forme et des sous-ensembles flexibles. Une approche à paramètres concentrés est utilisée, intégrant notamment des éléments spécifiques d’engrenage et des éléments d’arbre. Pour les éléments d’engrenage, le formalisme des fines tranches juxtaposées est employé pour représenter les dentures. Une raideur élémentaire et un écart normal sont attribués à chacune des tranches et sont réactualisés à chaque pas de temps en fonction de la cinématique des composants et des déviations du profil des dentures. Les déformations de la couronne sont introduites à l’aide d’un anneau discret composé de poutres droites couplé aux éléments d’engrenage. Les dentures double-hélice sont modélisées par deux éléments d’engrenage d’angles d’hélice opposés liés par une poutre de Timoshenko. Finalement, les équations du mouvement sont résolues pas à pas dans le temps par un schéma de Newmark combiné à un algorithme de contact normal, permettant de prendre en compte les pertes de contact partielles ou complètes. Dans un premier temps, un certain nombre d’éléments de validation est présenté et comparé à des résultats tirés de la littérature. Afin d’asseoir notre modélisation, l’influence des erreurs de positionnement des satellites, du décalage des hélices, des erreurs de pas et des déformations de la couronne sur les distributions de charge est abordée pour différentes configurations de train planétaire. Dans un deuxième temps, l’optimisation des corrections de denture dans le but de réduire les vibrations est investiguée. Les corrections de profil sont introduites sur les engrènements de façon (i) linéaire et symétrique en tête de dents et (ii) identique pour tous les satellites mais (iii) différente selon le flanc actif. Dans ce contexte, les corrections sont tout d’abord optimisées vis-à-vis des efforts dynamiques d’engrènement à l’aide d’un algorithme génétique. Puis, leurs performances sont analysées en fonction du couple transmis et de la vitesse de rotation. Par la suite, un critère « équivalent » est dérivé, vérifié et utilisé pour étudier l’influence du décalage des hélices et d’une correction longitudinale parabolique sur ces corrections de profil optimales. Enfin, une sous-structure du porte-couronne est introduite et son impact sur les distributions de charge est exploré. / This research work was conducted at the Contact and Structural Mechanics Laboratory (LaMCoS) of LaMCoS - INSA Lyon (UMR CNRS 5259) in partnership with Hispano-Suiza (SAFRAN group). In the context of new turbo jet engine developments, a promising technological solution consists in inserting a planetary / epicyclic gear train between the turbine and the propeller which, in theory, can improve the system performance, especially in terms of efficiency. Increasing power densities and mass reduction constraints lead to more compliant structures which need to be analyzed from a dynamic viewpoint. The present work deals therefore with the dynamic modelling of planetary / epicyclic gears and the effects of assembly errors, tooth shape deviations and deformable structural components. A lumped parameter approach has been favored which combines rigid-body gear elements, beam and lumped parameters elements. A thin-slice model has been used to simulate the time-varying elastic properties of gear teeth with an elemental stiffness and a normal deviation functions attributed to every discrete cell on the contact lines (thin slice) and updated at each time step with respect to the meshing course and the instant positions of the teeth. Ring-gear deformations are introduced via a model of elastic annulus discretized into straight beam elements and connected to the gear elements. Double-helical gears are simulated by linking two gear elements of opposite hands by Timoshenko beam elements. The possibility of helix stagger is implemented by shifting the helix positions in the base plane. Finally, the equations of motion are solved step by step in time by combining a Newmark scheme and a normal contact algorithm which makes it possible to account for partial and total instant contact losses. A number of comparisons with benchmark results from the literature are presented which prove that the proposed theoretical and numerical developments are sound and can actually be used to simulate the influence of planet position errors, helix stagger, pitch errors and ring-gear deformations. The optimization of tooth shape modifications, i.e. profile and lead modifications, re dynamic mesh forces in planetary gears is tackled. Using a genetic algorithm, optimum profile modifications are derived and compared with some analytical results in the literature. Their performance over a range of loads and speeds is assessed for helical and double helical gears with rigid and flexible ring-gears. A quasi-static “equivalent” criterion based on local transmission errors is presented and commented upon. Having proved its relevance, a number of results are derived concerning the influence of helix stagger and lead crowning superimposed on optimum profile modifications. Finally, a deformable ring-gear support is introduced using a sub-structuring technique and its contribution in terms of tooth load distribution is examined.

Avant la généralisation des toiles vers le XVI-XVIIème siècle, le bois a servi de support à d'innombrables peintures qui constituent aujourd'hui une part importante, dans les musées et les églises notamment, de notre patrimoine culturel. Après plusieurs siècles d'existence, les planches servant de support aux panneaux peints sont très souvent courbées : cette courbure est généralement imputée à la présence unilatérale de la couche picturale, imposant des échanges asymétriques d'humidité entre le bois, matériau hygroscopique, et son environnement de conservation. Par ailleurs il existe dans ces déformations une forte contribution de l'orthotropie cylindrique du bois qui est une conséquence de la croissance des arbres. Enfin, l'historique des variations hygrométriques à proximité d'un panneau peint renvoie à l'étude plus générale du comportement thermo-hygro-mécanique différé du bois. Une compréhension d'ensemble de ces phénomènes doit permettre d'orienter des décisions difficiles de conservation, notamment celles relatives au déplacement des œuvres ou à la maîtrise des environnements dans les musées. Nous proposons dans ce travail une approche générale de mécanique du matériau et des structures bois appliquée aux panneaux peints. La méthode utilisée intègre des moyens expérimentaux en laboratoire et sur des œuvres in-situ, des méthodes numériques, et l'accent est mis sur le lien fort existant entre les aspects " comportement " propres au bois (anisotropie, couplages hygromécaniques, etc.) et les aspects de structure relatifs aux panneaux peints (débit des planches, efforts extérieurs dus à leur assemblage, etc.). Un outil de calcul basé sur la modélisation mécanique des plaques orthotropes est développé. Il prend appui sur un code préexistant de transfert de masse et de chaleur décrivant les mouvements d'eau dans le support. Ce code de calcul complet est utilisé notamment sur un cas d'étude : un panneau peint de 500 ans environ, en situation d'exposition dans une église.

Cette HDR s'inscrit dans le cadre des actions de recherche pour la simulation des transitoires brutaux pour les structures et les fluides en interaction menées au Laboratoire d'Etudes de Dynamique du CEA, relatives à la définition de méthodes numériques pour la modélisation de systèmes mécaniques complexes et la résolution parallèle sur les supercalculateurs de problèmes de taille industrielle. Une particularité des approches proposées est la limitation à son minimum du nombre de paramètres non-physiques dans une simulation, pour s'accommoder des contraintes de maîtrise de la solution qu'impose le périmètre d'utilisation des concepts : sûreté nucléaire (CEA, EDF) ou aéronautique (ONERA, protection du citoyen (EC/JRC), en particulier. Ainsi, les contraintes cinématiques couplant fortement les structures entre elles (contact unilatéral par exemple) ou les fluides et les structures (avec des maillages conformes ou topologique déconnectés en fonction des situations géométriques) sont majoritairement traitées par l'intermédiaire de multiplicateurs de Lagrange, assurant la vérification exacte des équations de liaison au prix de la résolution d'un système additionnel variable dans le temps. Ce dernier aspect fait d'EPX (http://www-epx.cea.fr), le logiciel servant de réceptacle pour les méthodes, un outil à part dans la communauté du calcul en dynamique rapide. Le mémoire repose principalement sur une description des besoins en matière de modélisation pour la simulation de transitoires de référence, en particulier dans le monde du nucléaire, et des réponses apportées dans le cadre de la collaboration entre le CEA, EDF (via le LaMSID) et le LaMCoS. Sont ainsi considérés par exemple la déchirure d'un réservoir sous impact, l'accident de dimensionnement du confinement pour un réacteur de IVème génération ou la ruine d'une structure en béton armée sous impact. Sont ainsi proposés des modélisations innovantes et des algorithmes de résolution parallèles permettant de mettre en oeuvre avec efficacité les simulations correspondantes sur des calculateurs composés de noeuds multi-coeurs interconnectés sans jamais dégrader la qualité de la solution, ce qui a fait en particulier l'objet du projet ANR RePDyn (2010-2013), piloté par le CEA, à l'origine d'une collaboration étroite et en cours avec l'INRIA (Laboratoire d'Informatique de Grenoble).

Cette thèse est une contribution au développement d'un procédé tactile acoustique basé sur la perturbation des ondes de Lamb, permettant de détecter un toucher simple ou des touchers multiples sur des objets minces de forme quelconque. Après avoir étudié la propagation des modes de Lamb dans les objets minces par exemple les plaques de cuivre d'épaisseur 450 µm, nous avons observé, que l'énergie des ondes de Lamb antisymétriques était essentiellement due à la composante normale de déplacement et qu'il était aisé de perturber un mode de flexion par un toucher surfacique réalisé avec le doigt ou ponctuel avec la pointe d'un stylo, tandis que dans le cas d'un mode symétrique, l'énergie était essentiellement distribuée sur la composante de déplacement tangentielle et en conséquence plus difficilement perturbée par une interaction tactile. Nous avons en conséquence exploité le toucher sur une plaque mince perturbant le rayonnement d'une onde acoustique dans l'objet et réalisé en pratique une mesure de la perturbation d'un signal de diffraction des ondes de Lamb en un ou deux points récepteurs d'une plaque. Ainsi, une méthode de localisation basée sur la diffraction des ondes a été proposée, pour obtenir une interface tactile compacte, performante et faible coût. Cette méthode nommée méthode de la Figure de Diffraction Perturbée en Amplitude (FDPA) peut rendre une surface tactile avec une résolution spatiale de l'ordre du millimètre et un temps de réponse inférieur à 10 ms. Deux méthodes d'optimisation sont proposées pour obtenir un taux de localisation correcte proche de 100%. Le procédé tactile et les méthodes d'optimisation présentés dans cette thèse ont été testés sur plusieurs types d'objets, réalisés dans des matériaux différents et sous différentes formes géométriques. Ces réalisations ont montré des performances satisfaisantes en termes d'encombrement, de consommation et de fiabilité de localisation pour des touchers simples et multiples.

Ces travaux de thèse sont composés principalement de deux grandes parties de volumes inégaux dont voici le résumé des contenus : Partie I : Elle concerne les problèmes de stabilité et de bifurcation rencontrés dans les matériaux indépendants du temps physique (i.e., élastiques ou élasto-plastiques). Une analyse approfondie et détaillée des différentes études dans ce domaine nous permet de faire une présentation originale et synthétique de la théorie de la stabilité et de la bifurcation. Une illustration de cette théorie est ensuite donnée à travers une collaboration industrielle avec la SNECMA. On y effectue des simulations numériques du formage des aubes de réacteurs d'avions accompagnées d'analyses des états critiques de flambage. Les conclusions de ces simulations montrent clairement la nécessité de tenir compte du comportement visqueux du matériau et motivent ainsi la seconde partie de la thèse. Partie II : La deuxième partie des travaux de thèse est la plus importante quant au volume et au fond. Elle concerne les problèmes de stabilité des matériaux dépendant du temps physique (i.e., visco-élastiques ou visco-plastiques). L'objectif de cette étude est double : D'une part, une modélisation théorique du problème de stabilité des solides visqueux est proposée. Cette modélisation a pour but de poser le problème étudié comme un problème de stabilité de l'évolution quasi statique du solide visqueux considéré. En effet, l'absence d'équilibres pour de telles structures visqueuses nous amène naturellement à étudier la stabilité de leurs évolutions quasi statiques. Ces dernières dépendent bien entendu de la vitesse du chargement, ce qui constitue une source de difficultés dans une telle analyse. Cette façon de poser le problème en termes de stabilité de trajectoires est, à notre connaissance, originale et généralise la notion plus classique de la stabilité d'un équilibre. Elle consiste, en quelque sorte, à mesurer la sensibilité des mouvements aux perturbations extérieures. D'autre part, la recherche de critères assurant la stabilité de l'évolution visqueuse est menée en étroite liaison avec les résultats mathématiques relatifs aux équations différentielles non autonomes. Il est mis en évidence, notamment, que la principale difficulté dans cette analyse vient du caractère non autonome des équations différentielles régissant l'évolution quasi statique visqueuse. Les conditions suffisantes de stabilité auxquelles on aboutit sont obtenues principalement par deux méthodes. La première consiste à construire des fonctionnelles de Lyapunov adéquates à la physique du problème ou à utiliser la méthode de linéarisation autour de la réponse quasi statique. L'autre démarche est plus directe et permet d'établir le critère de stabilité de l'évolution par des considérations de continuité de la réponse par rapport aux données initiales. Cette deuxième démarche a l'avantage de s'appliquer aussi bien aux évolutions visco-plastiques qu'aux évolutions élasto-plastiques. Enfin, l'application de ces différents résultats aux modèles de tiges de Shanley et pour le cas des poutres permet d'illustrer notre méthodologie sur des exemples concrets et de comprendre la portée des critères de stabilité ainsi obtenus.