There are two types of electric crossing : i ) subterranean / submarine line ii ) overhead-line crossing. We always consider the last one as a more economic option. The inconvenience of an overhead-line crossing would be the environmental constraints among which the existing obstacles, the clearance for the navigation and the aesthetics demanded by the public. The overhead-line crossings usually have conductors of long ranges which are outside of the field of application for the current transmission line codes. These are limited to reaches of a length included between 200 m and 800 m, as well as a height of support lower than 60 m. However, for reaches over 800 m and over a height over 60 m, the criteria of conception in the transmission line codes for the calculation of wind loads are not applicable. In this study we concentrate on loads on the supports owed to the limit wind applied to bare conductors and insulators chains The objective of the present study is to examine the effect of the temporal and spatial correlation of the wind load along the conductors on a finite element model. A special attention was brought to the evaluation of the importance of the dynamic load transmitted on by the conductors and the insulators chains for the case of a turbulent wind load. The numerical study on finite element model for the example of a overhead-line crossing was done with the software ADINA. The wind load for the finite element model for the example of a overhead-line crossing was generated by the software WindGen which uses the method of Simiu-Scanlan and the method of spectral representation developed by Shinozuka-Deodatis. Wind loads generated where integrated into the finite element model ADINA for a dynamic analysis of the overhead-line crossing. For the first part, the current methods are used to calculate the efforts in supports due to the wind loads with an engineering approach and a comparaison approach. The current methods are then compared with the efforts obtained from an advanced method, transient dynamic and spectral stochastic, and specifically for the case of a simple overhead-line and an overhead-line crossings. For the second part, the effect of the longitudinal correlation of the wind load on two parallel conductors was examined. Finally, dynamic experiments on an insulators chain were made to determine the variation of the damping and the rigidity of the system for different type of insulators, different speed of application of the load and the inclination of the insulator.

Dans ce travail de thèse, une stratégie de calcul multi-échelle en dynamique transitoire basée sur la méthode LATIN est proposée. Ce travail fait suite, entre autre, à la thèse de H. Lemoussu qui a appliqué la méthode LATIN dans sa version mono-échelle au cas de la dynamique, ainsi qu'aux avancées plus récentes concernant l'introduction d'une vision à deux échelles au sein de la méthode LATIN pour des calculs statiques et quasi-statiques. Notre but a été d'étendre cette vision à deux échelles au cas de la dynamique transitoire. Une écriture de la stratégie de résolution multi-échelle est proposée pour le cas de la dynamique. Ce travail a permis de mettre en évidence plusieurs particularités de la méthode au sujet des conditions d'admissibilité des quantités macroscopiques en dynamique. L'introduction de l'approche multi-échelle en dynamique a permis d'obtenir l'extensibilité numérique de la méthode de décomposition de domaine pour le cas des interfaces. L'autre volet de cette thèse concerne la prise en compte de variabilité paramétrique en dynamique transitoire. Le but était ici de s'appuyer sur le savoir faire du LMT-Cachan, en matière de technique de calcul multi-résolution pour des problèmes avec contact pour développer une stratégie de calcul multi-résolution multi-échelle en dynamique, apte à prendre en compte les incertitudes, tout en diminuant de façon drastique le coût de calcul par rapport aux approches conventionnelles.

La stratégie de calcul proposée est une approche fréquentielle dans laquelle la partie basses fréquences est traitée grâce aux éléments finis classiques, alors que les moyennes fréquences sont calculées grâce à la Théorie Variationnelle des Rayons Complexes (TVRC) initialement introduite pour des problèmes de vibration forcée.

La simulation dynamique des structures composites jusqu'à rupture pose des difficultés, notamment numériques. Dans ce cadre, les mésomodèles d'endommagement ont été développés au LMT-Cachan, tant en statique qu'en dynamique, par l'introduction de lois à effet retard pour la description objective de la rupture.Ce travail discute l'identification de ces paramètres de retard, soulevant ainsi le problème de la prise en compte de conditions expérimentales fortement incertaines.Cette étude propose une méthode d'identification inverse à partir d'essais en dynamique transitoire, basée sur le principe d'erreur en relation de comportement. Elle a été développée, dans un premier temps, en élasticité. Cela amène à résoudre des problèmes dynamiques non standards demandant des méthodes de résolution adaptées. L'étude a permis de mettre en avant les points clés de la robustesse de la méthode face aux incertitudes de mesure.

Le flottement des tabliers de ponts flexibles constitue l'une des instabilités aérodynamiques les plus sévères pouvant survenir lors de tempêtes de vent. Le flottement est un phénomène dit auto-excité et se caractérise par de grandes amplitudes de déplacements dynamiques qui croissent avec la vitesse de vent et qui peuvent mener à la ruine de l'ouvrage. Avec l'évolution des matériaux, des techniques de constructions et des connaissances dans le domaine du génie civil, l'augmentation de la longueur des portées des ponts a pour effet de rendre ces structures plus flexibles. Historiquement, la participation du déplacement horizontal au flottement a été considérée comme négligeable. Toutefois, l'augmentation de la flexibilité latérale du tablier jumelée au couplage des modes dû aux caractéristiques des sections de tabliers et à l'interaction avec les câbles soulève quelques doutes concernant cette hypothèse. Le modèle mathématique actuellement utilisé pour représenter les forces auto-excitées de flottement a été développé par Scanlan et Tomko [1971]. Ce modèle définit des matrices de rigidité aérodynamique et d'amortissement aérodynamique qui varient en fonction de la vitesse réduite de l'écoulement. Ces deux matrices s'ajoutent aux matrices structurales et constituent la contribution aéroélastique de l'écoulement. Quelques recherches ont mis en évidence l'aspect non linéaire du flottement en démontrant la dépendance des coefficients instationnaires à la fréquence et à l'amplitude auxquelles ils ont été extraits. Les objectifs principaux de ce projet de recherche sont d'étudier l'effet du déplacement horizontal sur le flottement, d'étudier la dépendance des coefficients instationnaires à la fréquence et à l'amplitude d'extraction et d'élaborer une méthode expérimentale et un modèle mathématique permettant l'intégration temporelle non linéaire des forces auto-excitées. Des essais en régime libre sur 6 modèles sectionnels réguliers ont d'abord été réalisés pour évaluer l'effet du déplacement horizontal sur le flottement. Par la suite, une étude paramétrique de l'effet de la fréquence et de l'amplitude auxquelles les coefficients intationnaires sont extraits a été réalisée en régime forcé harmonique. Cette étude ayant mis en évidence la dépendance de certains coefficients à ces paramètres d'extraction, une dernière étude expérimentale a été réalisée. Cette dernière a permis l'extraction de coefficients de force dynamiques à partir de mouvements cycliques à vitesses ou accélérations absolues constantes. L'intégration de ces coefficients temporels, par opposition à fréquentiels, dans un outil de calcul dynamique transitoire non linéaire a permis la mise en évidence du flottement de sections de pont atteignant des amplitudes de vibration variant en fonction de la vitesse de l'écoulement. La vitesse de flottement obtenue concorde bien avec la moyenne des prédictions calculées à l'aide de la méthode classique. Un tel modèle permettrait à moyen terme de prendre en compte l'aspect non linéaire du flottement dans un modèle dynamique transitoire complet d'un pont, y incluant également les non-linéarités de nature structurale.

Pour les groupes turbo-alternateurs du parc nucléaire français, la situation accidentelle de référence lors du dimensionnement est définie par le départ d'une ailette terminale des corps basse pression. Après avoir détecté la perte de l'ailette, la turbine est déconnectée du réseau électrique, et la turbomachine ralentit progressivement sous l'effet de frottements d'origine fluide. Le dimensionnement consiste alors à vérifier que les structures palières sont capables de supporter l'effort résultant du balourd généré par cette perte d'ailette pendant la phase de ralentissement, notamment lors du passage des vitesses pour lequel des contacts rotor-stator peuvent se produire. A cet effet, un modèle éléments finis de turbine a été développé pour simuler des transitoires de vitesse d'un rotor en considérant la vitesse de rotation de l'arbre comme une inconnue, permettant ainsi d'estimer la décélération angulaire due au contact frottant. Trois modèles de stator ont été proposés appréhender le comportement dynamique du carter, et le couplage des structures est réalisé au moyen des forces de contact calculées par la méthode des multiplicateurs de Lagrange. La résolution des équations non-linéaires du mouvement a nécessité le développement d'un algorithme de contact, adapté à la géométrie de chacun des carters. Des simulations réalisées sur un modèle simple de turbine montrent que la vitesse de rotation et les chargements aux paliers sont bien estimés lors de l'interaction rotor-stator. Ainsi l'outil numérique permet d'étudier la réponse couplée rotor-stator et de prédire l'influence de paramètres tels que le coefficient de friction et les propriétés matérielles du stator.

Le cadre général de ce travail concerne la prédiction des ouvrages en béton armé soumis à des charges dynamiques sévères de type impact. Le dimensionnement de ces structures sous de telles conditions doit être capable de modéliser de manière fiable et efficace le comportement discontinu local ainsi que la réponse globale de l'ouvrage. Localement, la méthode des éléments discrets est utilisée pour représenter correctement les phénomènes discontinus. Cette méthode a été validée sur des essais quasi-statiques et dynamiques réalisés sur des structures de béton et de béton armé. <br /><br />Ce travail s'attache en particulier à développer le processus d'identification des paramètres locaux du modèle discret afin de le rendre plus prédictif. L'application de la méthode des éléments discrets sur des structures de grandes dimensions est limitée par le nombre d'éléments discrets nécessaire qui les rendent très peu efficaces. Pour garder la qualité de prédiction local et rendre cette méthode adaptée à l'échelle d'une structure, elle est couplée à la méthode des éléments finis dans les parties qui ne sont pas soumises à d'importantes discontinuités. Cette approche couplée permet des gains de temps importants, tant pour la modélisation que la pour la simulation. <br /><br />La méthode multi-domaine couplée ED/EF proposée utilise une zone de recouvrement sur laquelle les relations de couplage cinématique sur les déplacements et les rotations sont assurées par multiplicateurs de Lagrange. L'énergie totale du modèle couplé est partagée par une combinaison linéaire des énergies de la partie discrète et de la partie continue. <br /><br />Les réflexions d'ondes parasites liées à la différence de taille de discrétisation entre les deux approches sont atténuées par une méthode de relaxation. La validation de la méthode couplée porte sur des simulations d'impact sur des dalles en béton et en béton armé.

L'objectif de la thèse est de proposer une stratégie robuste permettant l'identification en dynamique transitoire de loi de comportement gouvernant la réponse du matériau jusqu'à la rupture. Il s'agit d'être à même d'identifier les paramètres matériaux supposés gouverner la rupture dans un contexte où les données expérimentales sont fortement incertaines. Faisant suite à un premier travail où la stratégie avait été élaborée dans un cadre élastique, le travail s'est concentré sur l'extension de la méthode d'identification pour les cas non linéaires tout d'abord la viscoplasticité puis les modéles d'endommagement à taux limités.<br />Les difficultés rencontrées dans ces cas résident dans la non-linéarité et le caractère instable du problème de minimisation sous contraintes non linéaires auquel la formulation nous amène. Une extension de la méthode LATIN aux problèmes mal posés a été proposée et développée afin de permettre la résolution itérative de ce type de problèmes d'optimisation. La résolution de ces derniers fait appel à une méthode de traitement robuste issue du contrôle optimal et basée sur l'équation de Riccati.<br />Une fois ces difficultés résolues et dans les cas simples unidimensionnels traités pour le moment, la stratégie d'identification proposée s'avère très robuste face aux perturbations des mesures même dans le cas très sévère de la localisation et de la rupture.

L’analyse isogéométrique pour la résolution de problèmes de la mécanique du solide suscite de vifs intérêts depuis une dizaine d’année. En effet, cette méthode de discrétisation autorise la description exacte des géométries étudiées permettant ainsi de supprimer les erreurs dues à une mauvaise description du domaine spatial étudié. Cependant elle pose un problème théorique de propagation de raffinement lors de la localisation de maillage. Des méthodes pour contourner ce problème ont été proposée dans la littérature mais complexifient grandement la mise en œuvre de cette stratégie de résolution. Cette thèse propose une stratégie de raffinement localisé adaptatif en espace pour les problèmes de statique et en espace temps pour les problèmes de dynamique transitoire dans le cadre de l’analyse isogéométrique. Pour cela une méthode de localisation pour l’analyse isogéométrique en statique basée sur une résolution multigrille est tout d’abord développée pour des problèmes en deux dimensions. Elle présente l’avantage de contourner la problématique de propagation de raffinement de maillage due à l’analyse isogéométrique tout en étant plus simple à mettre en œuvre que les méthodes déjà existantes. De plus, l’utilisation de l’analyse isogéométrique permet de simplifier les procédures de raffinement lors de l’adaptation de maillage qui peuvent être complexes lors de l’utilisationd’éléments finis classiques. Une méthode de raffinement adaptatif espace temps basée sur une résolution multigrille est ensuite développée pour des problèmes en une dimension. Une étude sur la structure des opérateurs est proposée afin de choisir un intégrateur temporel adapté. Les performances de cette stratégies sont mises en évidence, puis une modification de la méthode de résolution est proposée afin de diminuer significativement les coûts de calculs associées à cette résolution. La méthode de raffinement adaptatif espace temps est appliquée à quelques exemples académiques afin de valider son bon comportement lors de la localisation. / Isogeometric analysis applied to structural mechanics problems is a topic of intense concerns for a decade. Indeed, an exact description of geometries studied is allowed by this discretization method suppressing errors due to a bad description of the spatial domain considered. However, a theoretical problem of refinement propagation appears during mesh localization. Local refinement methods for isogeometric analysis has been developed and implied a complexification of the implementation of such a resolution strategy. This PhD thesis expose a space adaptative refinement strategy for linear elastic problems and a space-time one for transient dynamic using isogeometric analysis. For this purpose, a localization method for isogeometric analysis based on a multigrid resolution is developed for 2D linear elastic problems. This method allow to circumvent mesh refinement propagation inherent to isogeometric analysis, and is easier to implement than existing methods. Moreover, the use of isogeometric analysis simplifies refinement procedures occuring during mesh adaptation and which can be really complex using classical finite element analysis. Then, a space-time adaptative refinement based on a multigrid resolution is developed for one dimensional in space problems. A study on operators structure is exposed in order to choose a well suited time integrator. This strategy's performances are highlighted, then an evolution of this method is set up in order to lower computational costs. The space-time adaptaptive refinement is applied to some academical examples to show it good behavior during localization.

Les problèmes de machines tournantes présentant un contact rotor-stator sont généralement traités dans la littérature à partir de modèles localement 1D. Ceci conduit à un temps de simulation acceptable, mais les approximations correspondantes sont difficiles à évaluer. Or le contact est limité en espace et en temps, il peut donc en être de même pour l'utilisation d'une modélisation 3D. Une stratégie, nommée bascule, permettant d'utiliser un modèle poutre et un autre 3D (ou un modèle mixte poutre-3D), pendant deux phases différentes de la même simulation, est donc développée. Elle est mise en oeuvre de façon non-intrusive, en dynamique transitoire, pour des systèmes présentant ou non une rotation d'ensemble, avec une résolution par intégration temporelle implicite. La bascule est validée par comparaison avec une solution 3D de référence obtenue en effectuant la simulation entière sur le modèle 3D. On analyse ensuite, sur des problèmes de contact rotor-stator, les résultats des modèles 1D et 3D, dans différentes situations de contact plus ou moins sévères. Les vibrations du rotor sont dues à un balourd et sa vitesse de rotation est imposée constante. Il apparaît que les écarts entre les résultats des modèles 1D et 3D sont peu visibles sur les orbites du rotor. En revanche les limitations du modèle 1D sont mises en évidence. En effet, l'hypothèse de section rigide du rotor entraine des approximations dans la distribution spatiale et l'intensité des efforts de contact. La modélisation 3D en revanche permet de représenter certains effets locaux au voisinage de la zone de contact.