Les fonctions D-finies (ou holonomes) à une variable sont les solutions d'équations différentielles linéaires à coefficients polynomiaux. En calcul formel, il s'est avéré fructueux depuis une vingtaine d'années d'en développer un traitement algorithmique unifié. Cette thèse s'inscrit dans cette optique, et s'intéresse à l'évaluation numérique des fonctions D-finies ainsi qu'à quelques problèmes apparentés. Elle explore trois grandes directions. La première concerne la majoration des coefficients des développements en série de fonctions D-finies. On aboutit à un algorithme de calcul automatique de majorants accompagné d'un résultat de finesse des bornes obtenues. Une seconde direction est la mise en pratique de l'algorithme " bit burst " de Chudnovsky et Chudnovsky pour le prolongement analytique numérique à précision arbitraire des fonctions D-finies. Son implémentation est l'occasion de diverses améliorations techniques. Ici comme pour le calcul de bornes, on s'attache par ailleurs à couvrir le cas des points singuliers réguliers des équations différentielles. Enfin, la dernière partie de la thèse développe une méthode pour calculer une approximation polynomiale de degré imposé d'une fonction D-finie sur un intervalle, via l'étude des développements en série de Tchebycheff de ces fonctions. Toutes les questions sont abordées avec un triple objectif de rigueur (résultats numériques garantis), de généralité (traiter toute la classe des fonctions D-finies) et d'efficacité. Pratiquement tous les algorithmes étudiés s'accompagnent d'implémentations disponibles publiquement.

4 chapitres : 1. Quelques notions sur les modélisations de la houle. 2. Le problème élastique. 3. Stabilité d'un massif de sol cohérent. - 4. Stabilité d'un massif de sol granulaire.

Les éléments de coques multicouches ont montré une grande efficacité dans l'analyse non linéaire des structures. A l'origine, ces éléments avaient été développés pour des lois de comportement spécifiques au béton et à l'acier. Leur généralisation à l'analyse des structures dont le comportement est régi par des lois de comportement à variables internes (plasticité, élasto-visco-plasticité, endommagement) constitue le propos de ce travail de thèse. Issus de l'approche semi-globale, les éléments de coques multicouches permettent de rendre compte des phénomènes non-linéaires tant matériels que géométriques pour l'analyse des structures constituées de coques, sous divers chargements, avec un temps de calcul raisonnable. Dans une première partie, nous développons ces éléments dans le cadre de l'hypothèse des petites perturbations. Une extension est ensuite présentée en intégrant les effets non-linéaires géométriques. Une description lagrangienne actualisée avec un traitement semi-tangentiel des paramètres de rotation est utilisée à cet effet. Dans une seconde partie, nous présentons dans un cadre général, l'aspect théorique ainsi que le traitement numérique des lois de comportement à variables internes. A partir de là, un ensemble de lois de comportement à variables internes permettant la modélisation des comportements de l'acier et du béton sous divers chargements, est présenté. Le but étant d'utiliser ces modèles au sein des éléments finis de coques multicouches, pour le calcul de structures coques en béton armé soumises à des chargements statiques et dynamiques. Une attention particulière est portée sur la condition de contraintes planes d'une part dans le schéma d'intégration des contraintes, et d'autre part dans le calcul du tenseur de comportement tangent. Enfin, les développements effectués sont implantés dans le code de calcul par éléments finis CESAR-LCPC. Au travers d'exemples d'applications numériques, nous précisons le domaine d'application ainsi que les limites des outils proposés.

Cette thèse porte sur l'étude théorique et expérimentale des structures en béton armé sous sollicitations statique puis dynamique de type séismique. Une importante recherche bibliographique a été effectuée en vue d'analyser au mieux le comportement dynamique des matériaux ainsi que celui des éléments structuraux. Cette recherche est complétée par l'étude des structures essayées sur table séismique. Les enseignements qui ont été tirés de cette étude préalable, notamment ceux concernant la modification de la rigidité des structures Béton Armé sous chargement alterné et séismique, ont permis d'établir une méthode de calcul (dite ''statique équivalente") pour l'analyse de la tenue aux séismes des structures en Béton Armé. Cette méthode tient compte, en particulier de la non-linéarité du comportement des matériaux. Les réponses au séisme qui ont été obtenues par cette méthode sur des portiques essayés sur table vibrante, concordent de façon très satisfaisante avec les valeurs d'essais.

Cette étude est consacrée à l'optimisation géométrique du contrôle actif dans les gaines de ventilation. La première partie de ce travail concerne l'étude du contrôle hybride (actif+passif) dans les gaines de ventilation et de ses avantages suivant la géométrie de la gaine et la présence ou non de revêtements absorbants. Les effets du contrôle passif du bruit (dissipatif, réflexif et diffusif) dans une gaine de ventilation améliorent l'efficacité du contrôle actif. Un contrôle hybride (actif + passif) présente donc des avantages, même aux basses fréquences. La seconde partie de ce mémoire s'intéresse au placement des microphones d'erreur et des sources secondaires. Pour une excitation harmonique, une méthode de calcul par programmation linéaire et entière détermine un nombre suffisant de microphones d'erreur ainsi qu'une position optimale des microphones d'erreur et des sources secondaires. Cette méthode de placement est appliquée pour la réduction du bruit dans les gaines de ventilation.

Ce travail présente l'application de concepts modernes du génie logiciel au développement d'un environnement de simulation modulaire de systèmes thermiques complexes. Une application stricte des idées de la Conception Orientée Objet permet la réalisation d'un environnement modulaire, souple et efficace. Une analyse hiérarchique du système à étudier nous fournit un graphe de dépendance entre les composants du système. Chaque composant est confiné dans un objet informatique et a son propre modèle de calcul. Le travail essentiel du simulateur est la résolution des conditions de connexion entre les composants. La structure hiérarchique est conservée et exploitée pour une parallélisation des calculs pendant la simulation numérique. Pour une meilleure capitalisation des modèles développés, nous nous reposons sur une « modélothèque » à partir de laquelle l'utilisateur peut récupérer des modèles pour construire son propre système de simulation. Des exemples d'application de notre environnement de simulation pour des problèmes thermiques de bâtiment sont présentés.

Le but de cette thèse est l'étude prédictive par calcul aux éléments finis des effets de la fissuration à l'échelle de structures en béton armé et précontraint soumises à des chargements statiques, cycliques et dynamiques. Il s'agit d'une part de préciser les lois de comportement des matériaux constitutifs de ces structures, et d'autre part de développer des outils numériques adaptés. Pour la modélisation du comportement du béton, un modèle du comportement macroscopique est présenté, couplant la plasticité à l'endommagement. Pour la partie plastique, il s'agit d'adapter la loi élastoplastique au comportement spécifique du béton. La (micro)fissuration du matériau est représentée à l'échelle macroscopique en termes de déformations plastiques, et l'apparition et l'évolution du dommage par l'évolution de variables plastiques. En particulier, la porosité plastique modélise la variation irréversible de l'espace poreux connecté créé par (micro)fissuration. Cette signification physique est à la base du couplage phénoménologique de la plasticité et de l'endommagement : les effets du dommage sont modélisés par une variation des caractéristiques élastiques fonction de la variable d'endommagement choisie : la porosité plastique. Cette loi de comportement est utilisée au sein d'un élément fini particulier : l'élément poutre multifibre, issu de l'extension des approches multicouches au cas tridimensionnel. L'outil numérique développé permet de prendre en compte des phénomènes non-linéaires tant matériels que géométriques dans l'analyse de structures à poutres sous des chargement divers avec un temps de calcul raisonnable. Pour l'application au cas des structures en béton précontraint, la formulation prend en compte des déplacements relatifs (glissement) entre câble précontraint et béton : la précontrainte est traitée comme un problème de condition aux limites à l'interface acier-béton. Ces développements sont mis en oeuvre dans un code de calcul par éléments finis. Au travers d'exemples d'applications numériques, le domaine d'application et les limites des outils proposés sont précisés. Ces outils se veulent une aide à la conception pour les ingénieurs, dans les études prédictives des effets du dommage à l'échelle des structures.

La méthode d'homogénéisation en calcul à la rupture permet de calculer les propriétés de résistance à l'échelle macroscopique d'un matériau hétérogène périodique (matériaux composites, sols renforcés) à partir de la solution d'un problème de calcul à la rupture défini sur la période de base. Une méthode numérique spécifique destinée à résoudre un tel problème est mise au point. Elle est fondée sur la mise en oeuvre de l'approche cinématique du calcul à la rupture et conduit : 1) à la recherche du minimum d'une fonction convexe d'un nombre fini de paramètres scalaires pour laquelle un algorithme de résolution est proposé. 2) à une évaluation par l'extérieur du convexe de résistance macroscopique qui se révèle d'autant plus précise que le nombre de paramètres de minimisation est élevé. On procède à de nombreuses applications de cette méthode à des problèmes de "contrainte plane" (plaques renforcées ou perforées chargées dans leur plan) et son extension prometteuse à des problèmes en "déformation plane" est esquissée.

A l'intérieur d'un solide, l'endommagement est provoqué par des mouvements microscopiques. Nous décidons de tenir compte de ces mouvements dans le principe des puissances virtuelles. Les vitesses microscopiques étant reliées à la vitesse d'endommagement, la puissance des efforts intérieurs que nous retenons dépend de cette vitesse et aussi de la vitesse du gradient de l'endommagement pour décrire les interactions. L'utilisation de cette approche permet de représenter le comportement non linéaire du béton, qui est en grande partie dû à l'endommagement, par des modèles adoucissants. Les applications aux calculs de structures donnent des résultats objectifs surmontant la difficulté liée à la dépendance du maillage lors de la discrétisation par éléments finis. La prise en compte du gradient de l'endommagement dans la formulation permet également de rendre compte de l'effet d'échelle de structures. Le phénomène unilatéral lié à la refermeture des microfissures, qui conduit à la restauration de la raideur lors du passage d'un chargement endommageable en traction à un chargement en compression, est pris en compte par l'utilisation de deux variables scalaires d'endommagement et par une partition du tenseur des déformations dans l'expression de l'énergie libre. Les modèles que nous présentons sont implantés dans un code de calcul par éléments finis et les résultats obtenus permettent une analyse prévisionnelle des structures.

L'étude de l'équilibre à long terme du tunnel revêtu suscite un intérêt important. D'un point de vue mécanique, il s'agit d'un problème d'interaction complexe entre deux structures de comportement très différent : le massif avec la cavité d'une part, le soutènement d'autre part. Les effets différés qui apparaissent dans un massif sont principalement liés à deux phénomènes qui peuvent être couplés : - la viscosité du squelette solide du massif ; - la diffusion du fluide interstitiel au sein du massif. On commence par présenter une étude bibliographique assez générale concernant les différentes étapes d'étude, de conception et de réalisation d'un projet de tunnel profond. Dans la seconde partie du mémoire, on étudie l'influence de l'histoire du chargement (vitesse de creusement et instant de pose du soutènement) sur l'équilibre d'un tunnel soutenu creusé dans un massif monophasique à comportement viscoplastique. Une méthode simplifiée de dimensionnement des tunnels dans de tels milieux a été élaborée. On a également abordé les effets de couplage entre comportements hydraulique et mécanique dans le cadre des tunnels profonds creusés dans un milieu poreux saturé isotherme ayant un comportement poro-élastique ou poro-élasto-plastique. Une étude analytique du tunnel 1D révèle l'influence de l'histoire du chargement mécanique et des conditions aux limites hydrauliques imposées en paroi. Des simulations numériques de creusement de tunnel en géométrie 2D axisymétrique sont menées. L'étude analytique permet, grâce à la généralisation de la notion de contrainte fictive de soutènement, de localiser le point d'équilibre du tunnel soutenu dans le plan (convergence en paroi - contrainte appliquée en paroi). Ces résultats permettent d'aborder le calcul en déformation plane des tunnels creusés en milieux biphasiques avec plus de rigueur.