In this thesis, a miniaturized energy harvester system is developed. The energy harvester converts fluid kinetic energy into electrical energy without using any rotating components. The working principle of the energy harvester is based on the so called vortex-induced vibration. Such systems have the potential to provide energy for wireless sensor networks in the field of inline measurements for gas, oil or water transportation systems. The theoretical background of the vortex-induced vibration (VIV) is studied. Based on the studies, a fluid-structure interaction simulation is carried out to optimize the structure of the energy harvester. As result, the conversion efficiency is significantly improved, which is experimentally confirmed. A series of demonstrators are manufactured according to the simulation and optimization results. It is tested on a self-constructed test bench. To further improve the performance, an electromagnetic generator is proposed, and therefore, a multimethod demonstrator realized. The demonstrators are working in air flow already at a velocity of 2 m/s, and reach the maximum efficiency at 3.6 m/s. This performance ranks among the best published results and is discussed in detail. / In der vorliegenden Arbeit wird ein miniaturisiertes Energiegewinnungssystem entwickelt, das unter Verzicht auf rotierende Komponenten kinetische Strömungsenergie in elektrische Energie umwandelt. Die Funktion dieses Wandlers basiert auf der sogenannten wirbelinduzierten Vibration. Derartige Systeme besitzen unter anderem das Potenzial, drahtlose Sensornetzwerke zur Erfassung von Messdaten in Gas-, Öl- oder Wassertransportsystemen mit Energie zu versorgen zu können. In der Arbeit wird der theoretische Hintergrund der wirbelinduzierten Vibration untersucht und darauf basierend werden Fluid-Struktur-Wechselwirkungssimulationen zur Strukturoptimierung durchgeführt in deren Ergebnis eine theoretische Verbesserung der Effizienz des Wandlers um ein Mehrfaches erreicht wird, die auch praktisch bestätigt wird. Unter Berücksichtigung der Simulations- und Optimierungsergebnisse wurden eine Reihe von Demonstratoren gefertigt, die auf einem selbst konstruierten Prüfstand getestet wurden. Zur weiteren Erhöhung der Leistungsfähigkeit des Wandlers wird ein zusätzlicher elektromagnetischer Generator vorgeschlagen und damit ein Multi-Methoden-Demonstrator technisch realisiert. Die Demonstratoren arbeiten in strömender Luft bereits bei Geschwindigkeiten von 2 m/s und erreichen bei 3,6 m/s ihre maximale Effizienz. Die erreichten Ergebnisse ordnen sich im Vergleich mit denen aus entsprechenden Publikationen vorn ein und werden ausführlich diskutiert.

Cette thèse est consacrée à la modélisation de la ventilation mécanique chez l'humain et à l'analyse numérique des systèmes en découlant. Des simulations directes d'écoulement d'air dans l'ensemble des voies aériennes étant impossibles (maillages indisponibles et géométrie trop complexe), il est nécessaire de considérer un domaine d'intérêt réduit, qui implique de travailler dans une géométrie tronquée, comportant des frontières artificielles ou encore de considérer des modèles réduits simples mais représentatifs. Si on cherche à effectuer des simulations numériques 3D où l'écoulement du fluide est décrit par les équations de Navier-Stokes, différentes problématiques sont soulevées :- Si on considère que la ventilation est la conséquence de différences de pression, les conditions aux limites associées sont des conditions de type Neumann. Cela aboutit à des questions théoriques en terme d'existence et d'unicité de solution et à des questions numériques en terme de choix de schémas et de méthodes adaptées.- Lorsque l'on travaille dans un domaine tronqué, il peut être nécessaire de prendre en compte les phénomènes non décrits grâce à des modèles réduits appropriés. Ici nous considérons des modèles 0D. Ces couplages 3D/0D sont à l'origine d'instabilités numériques qu'on étudie mathématiquement et numériquement dans ce manuscrit. Par ailleurs, lorsqu'on s'intéresse à des régimes de respiration forcée, les modèles usuels linéaires sont invalidés par les expériences. Afin d'observer les différences entre les résultats expérimentaux et numériques, il est nécessaire de prendre en compte plusieurs types de non linéarités, comme la déformation du domaine ou les phénomènes de type Bernoulli. Une approche par modèles réduits est adoptée dans ce travail.Pour finir, on a cherché à valider les modèles obtenus en comparant des résultats numériques et des résultats expérimentaux dans le cadre d'un travail interdisciplinaire.Parvenir à modéliser et simuler ces écoulements permet de mieux comprendre les phénomènes et paramètres qui entrent en jeu lors de pathologies (asthme, emphysème...). Un des objectifs à moyen terme est d'étudier l'influence du mélange hélium-oxygène sur le dépôt d'aérosol, toujours dans le cadre du travail interdisciplinaire. A plus long terme, l'application de ces modèles à des situations pathologiques pourrait permettre de construire des outils d'aide à la décision dans le domaine médical (compréhension de la pathologie, optimisation de thérapie...). / In this thesis, we study the modelling of the human mecanical ventilation and the numerical analysis of linked systems. Direct simulations of air flow in the whole airways are impossible (complex geometry, unavailable meshes). Then a reduced area of interest can be considered, working with reduced geometries and artificial boundaries. One can also use reduced models, simple but realistic. If one try to make 3D numerical simulations where the fluid flow is described by the Navier-Stokes equations, various issues are raised:- If we consider that ventilation is the result of pressure drops, the associated boundary conditions are Neumann conditions. It leads to theoretical questions in terms of existence and uniqueness of solution and numerical issues in terms of scheme choice and appropriate numerical methods.- When working in a truncated domain, it may be necessary to take into account non-described phenomena with appropriate models. Here we consider 0D models. These 3D/0D couplings imply numerical instabilities that we mathematically and numerically study in this thesis.Furthermore, when we focus on forced breathing, linear usual models are invalidated by experiments. In order to observe the differences between the experimental and numerical results, it is necessary to take into account several types of non-linearities, such as deformation of the domain or the Bernoulli phenomenon. A reduced model approach is adopted in this work. Finally, we sought to validate the obtained models by comparing numerical and experimental results in the context of interdisciplinary work.Achieving model and simulate these flows allow to better understand phenomena and parameters that come into play in diseases (asthma, emphysema ...). A medium-term objective is to study the influence of helium-oxygen mixture in the aerosol deposition. In the longer term, the application of these models to pathological situations could afford to build decision support tools in the medical field (understanding of pathology, therapy optimization ...).

O presente trabalho consiste no uso do Método dos Elementos Finitos (MEF) para a análise de interação fluido-estruturas de barras com foco em problemas transientes envolvendo explosões ou outras ações com propagação de ondas de choque. Para isso é necessário o estudo de três diferentes aspectos: a dinâmica das estruturas computacional, a dinâmica dos fluidos computacional e o problema do acoplamento. No caso da dinâmica das estruturas computacional deve-se identificar em função da cinemática de deformações, quais são os requisitos para que um elemento seja adequado para analisar tais problemas, tendo em vista que a formulação deve admitir grandes deslocamentos. Para evitar problemas relacionados com aproximações de rotações finitas, opta-se por empregar uma formulação descrita em termos de posições e que leva em consideração os efeitos de empenamento da seção transversal. No caso da dinâmica dos fluidos computacional, busca-se uma formulação para escoamentos compressíveis que seja estável e ao mesmo tempo sensível ao movimento da estrutura, sendo empregado um algoritmo de integração temporal explícito baseado em características com as equações governantes descritas na forma Lagrangeana-Euleriana Arbitrária (ALE). No que se refere ao acoplamento, busca-se modularidade e versatilidade, empregando-se um modelo particionado fraco (explícito) de acoplamento e técnicas de transferência das condições de contorno (Dirichlet-Neummann), sendo estudados os efeitos de utilizar transferência bidirecional ou unidirecional dessas condições de contorno. / This work consists in the use of the Finite Element Method (FEM) for numerical analysis of fluid-bar structures, focusing on transient problems involving explosions or other actions with shock waves propagation. For this purpose, one needs to study three different aspects: the computational structural dynamics, the computational fluid dynamics and the coupling problem. Regarding computational structural dynamics, one need firstly to identify the requirements for an element to be adequate to analyze such problems, taking into account the fact that such element should admit large displacements. In order to avoid problems related to finite rotation approximations and to give a realist representation of a 3D bar structure, we chose a formulation defined in terms of positions and that considers the cross-section warping effects. Regarding computational fluid dynamics, we seek for a stable formulation for compressible flows, and at same time, sensitive to the movement of the structure, leading to an explicit time integration algorithm based on characteristics with governing equations described in the Arbitrary Lagrangian-Eulerian (ALE) form. Regarding to coupling, we chose to use a weak (explicit) partitioning coupling model in order to ensure modularity and versatility. The developed coupling scheme is bases on boundary conditions transfer techniques (Dirichlet-Neummann), and we study the effects of using bidirectional or unidirectional boundary conditions transfers.

Le but de cette thèse pluridisciplinaire est d’étudier le problème de l’interaction fluide-structure à partir du point de vue mathématique et physique. Des problèmes d’interaction d’un fluide visqueux avec une structure élastique décrivent, par exemple, des interactions entre le manteau terrestre et de la croûte terrestre, le sang et la paroi vasculaire dans un vaisseau sanguin, etc. En génie l’interaction fluide visqueux-structure apparaît lors de la formation de solution colloïdale quand un laser passe à travers le fluide influençant le substrat (ablation laser dans un liquide). Fusion sélective au laser (FSL) est utilisée pour étudier le comportement des contraintes résiduelles en dépendance des propriétés thermoélastiques et mécaniques du matériau et des formes variées des cordons rechargés. A partir du point de vue mathématique le système couplé “flux fluide visqueux – plaque mince élastique” en 3D lorsque l’épaisseur de la plaque, E, tend vers zéro, tandis que la densité et le module de Young du matériau élastique sont d’ordre 1 et E-3, respectivement, est considéré. Le solide est couché par le fluide qui occupe un domaine épais. La modélisation multi-échelle est effectuée pour la partie élastique. Le développement asymptotique complet est construit lorsque E tend vers zéro. L’existence, la régularité et l’unicité de la solution pour le problème initial sont étudiées au moyen de techniques variationnelles. La méthode de décomposition asymptotique partielle du domaine est appliquée pour le système couplé. L’erreur de la méthode est évaluée / The goal of this multi-disciplinary thesis is to study the fluid-structure interaction problem from mathematical and physical viewpoints. Viscous fluid-structure interaction problems describe, for example, interactions between the Earth mantle and the Earth crust, the blood and the vascular wall in a blood vessels, etc. In engineering viscous fluid-structure interaction appears during colloidal solution formation when a laser pierce through the fluid influencing the substrate (laser ablation in a liquid). Selective laser melting (SLM) is used to study the behavior of residual stresses depending on the thermoelastic and mechanical properties of the material and on various forms of reloaded beads. From mathematical point of view the coupled system “viscous fluid flow-thin elastic plate” in 3D when the thickness of the plate, E, tends to zero, while the density and the Young’s modulus of the plate material are of order 1 and E-3, respectively, is considered. The plate lies on the fluid which occupies a thick domain. The multi-scale modeling is performed for the elastic part. The complete asymptotic expansion is constructed when E tends to zero. The existence, the regularity and the uniqueness of the solution for the original problem are studied by means of variational techniques. The method of asymptotic partial domain decomposition is applied for the coupled system. The error of the method is evaluated

A prevalent morphology in the microscopic world of artificial microswimmers, bacteria and viruses is that of a helix. The intriguingly different physics at play at the small scale level make it necessary for bacteria to employ swimming strategies different from our everyday experience, such as the rotation of a helical filament. Bio-inspired microswimmers that mimic bacterial locomotion achieve propulsion at the microscale level using magnetically actuated, rotating helical filaments. A promising application of these artificial microswimmers is in non-invasive medicine, for drug delivery to tumours or microsurgery. Two crucial features need to be addressed in the design of microswimmers. First, the ability to selectively control large ensembles and second, the adaptivity to move through complex conduit geometries, such as the constrictions and curves of the tortuous tumour microvasculature. In this dissertation, a mechanics-based selective control mechanism for magnetic microswimmers is proposed, and a model and simulation of an elastic helix passing through a constricted microchannel are developed. Thereafter, a theoretical framework is developed for the propulsion by stiff elastic filaments in viscous fluids. In order to address this fluid-structure problem, a pertubative, asymptotic, elastohydrodynamic approach is used to characterise the deformation that arises from and in turn affects the motion. This framework is applied to the helical filaments of bacteria and magnetically actuated microswimmers. The dissertation then turns to the sub-bacterial scale of bacteriophage viruses, 'phages' for short, that infect bacteria by ejecting their genetic material and replicating inside their host. The valuable insight that phages can offer in our fight against pathogenic bacteria and the possibility of phage therapy as an alternative to antibiotics, are of paramount importance to tackle antibiotics resistance. In contrast to typical phages, flagellotropic phages first attach to bacterial flagella, and have the striking ability to reach the cell body for infection, despite their lack of independent motion. The last part of the dissertation develops the first theoretical model for the nut-and-bolt mechanism (proposed by Berg and Anderson in 1973). A nut being rotated will move along a bolt. Similarly, a phage wraps itself around a flagellum possessing helical grooves, and exploits the rotation of the flagellum in order to passively travel along and towards the cell body, according to this mechanism. The predictions from the model agree with experimental observations with respect to directionality, speed and the requirements for succesful translocation.

L’objectif de ce travail porte sur le développement d’un modèle d’interaction fluide-structure adapté à la dynamique des éoliennes de grandes tailles avec des pales flexibles qui se déforment de manière significative sous l’effet de la pression exercée par le vent. Le modèle développé est basé sur une approche efficace d’IFS partitionnée pour un fluide incompressible et non visqueux en interaction avec une structure flexible soumise a des grandes transformations. Il permet de fournir une meilleure estimation de la charge aérodynamique et de la réponse dynamique associée du système (pales, mat, attachements, câbles) avec un temps de calcul raisonnable et pour des simulations sur des longues périodes. Pour la modélisation structurale, un élément fini de type solide 3D est développé pour l’étude dynamique des pales d’éolienne soumises à des grands déplacements et des grandes rotations. Une amélioration du comportement en flexion est proposée par l’introduction des degrés de liberté en rotation et l’enrichissement du champ de déplacements afin de décrire plus précisément la flexibilité des pales. Cet élément solide est apte de capter des modes de hautes fréquences qui peuvent s’avérer néfastes pour la stabilité du calcul. Deux techniques sont donc proposées pour les contrôler : la régularisation de la matrice masse et le développement des schémas d’intégration robustes de conservation et de dissipation d’énergie. Les chargements aérodynamiques sont modélisés en utilisant la Panel Method. Il s’agit d’une méthode aux frontières, relativement rapide par rapport à la CFD mais suffisamment précise pour calculer la distribution de la pression exercée sur la pale. Les modèles fluide et structure interagissent via un algorithme de couplage partitionné itératif dans lequel des considérations particulières sont prises en compte dans le contexte des grandes transformations. Dans un effort visant à instaurer un indicateur de fatigue dans la méthodologie proposée, des câbles précontraints sont introduits reliant le mat de l’éolienne au support. Une nouvelle formulation complémentaire en termes de contraintes est ainsi développée pour l’analyse dynamique des câbles 3D en comportement élasto-visco-plastique. Chaque méthode proposée a été d’abord validée sur des cas tests pertinents. Par la suite, des simulations numériques d’éoliennes avec des pales flexibles sont effectuées en vue d’affiner la compréhension de leur comportement dynamique et l’intérêt que la flexibilité des pales peut apporter à leur fonctionnement. / In this work, a numerical model of fluid-structure interaction is developed for dynamic analysis of giant wind turbines with flexible blades that can deflect significantly under wind loading. The model is based on an efficient partitioned FSI approach for incompressible and inviscid flow interacting with a flexible structure undergoing large transformations. It seeks to provide the best estimate of true design aerodynamic load and the associated dynamic response of such system (blades, tower, attachments, cables). To model the structure, we developed a 3D solid element to analyze geometrically nonlinear statics and dynamics of wind turbine blades undergoing large displacements and rotations. The 3D solid bending behavior is improved by introducing rotational degrees of freedom and enriching the approximation of displacement field in order to describe the flexibility of the blades more accurately. This solid iscapable of representing high frequencies modes which should be taken under control. Thus, we proposed a regularized form of the mass matrix and robust time-stepping schemes based on energy conservation and dissipation. Aerodynamic loads are modeled by using the 3D Vortex Panel Method. Such boundary method is relatively fast to calculate pressure distribution compared to CFD and provides enough precision. The aerodynamic and structural parts interact with each other via a partitioned coupling scheme with iterative procedure where special considerations are taken into account for large overall motion. In an effort to introduce a fatigue indicator within the proposed framework, pre-stressed cables are added to the wind turbine, connecting the tower to the support and providing more stability. Therefore, a novel complementary force-based finite element formulation is constructed for dynamic analysis of elasto-viscoplastic cables. Each of theproposed methods is first validated with differents estexamples.Then,several numerical simulations of full-scale wind turbines are performed in order to better understand its dynamic behavior and to eventually optimize its operation.

Cette thèse de doctorat porte sur le calcul par la méthode des éléments finis du comportement dynamique de réservoirs élastiques précontraints contenant un liquide interne à surface libre. Nous considérons que la pression hydrostatique exercée par le fluide interne incompressible sur les parois flexibles du réservoir est à l’origine de grands déplacements, conduisant ainsi à un état d’équilibre non-linéaire géométrique. Le changement de raideur lié à cet état précontraint induit un décalage des fréquences de résonances du problème de vibrations linéaires couplées.L’objectif principal du travail est donc d’estimer, par des approches numériques précises et efficaces, l’influence des non-linéarités géométriques sur le comportement hydroélastique du système réservoir/liquide interne autour de différentes configurations d’équilibre. La méthodologie développée s’effectue en deux étapes. La première consiste à calculer l’état statique non-linéaire par une approche éléments finis lagrangienne totale. L’action du fluide sur la structure est ici modélisée par des forces suiveuses hydrostatiques. La deuxième étape porte sur le calcul des vibrations couplées linéarisées. Un modèle d’ordre réduit original est notamment proposé pour limiter les coûts de calcul associés à l’estimation de l’effet de masse ajoutée. Enfin, divers exemples sont proposés et comparés à des résultats de la littérature (issus de simulations numériques ou d’essais expérimentaux) pour montrer l’efficacité et la validité des différentes approches numériques développées dans ce travail. / This doctoral thesis focuses on the calculation by the finite element method of the dynamic behavior of prestressed elastic tanks containing an internal liquid with a free surface. We consider that the hydrostatic pressure exerted by the incompressible internal fluid on the flexible walls of the tank causes large displacements, thus leading to a geometric non-linear equilibrium state. The change of stiffness related to this prestressed state induces a shift in the resonance frequencies of the coupled linear vibration problem. The main objective of the work is therefore to estimate, through precise and efficient numerical approaches, the influence of geometric nonlinearities on the hydroelastic behavior of the reservoir/internal liquid system around different equilibrium configurations. The methodology developed is carried out in two stages. The first one consists in calculating the non-linear static state by a total Lagrangian finite element approach.The action of the fluid on the structure is modelled here by hydrostatic following forces. The second step is the calculation of linearized coupled vibrations. In particular, an original reduced order model is proposed to limit the calculation costs associated with the estimation of the added mass effect. Finally, various examples are proposed and compared with results from the literature (from numerical simulations or experimental tests) to show the effectiveness and validity of the different numerical approaches developed in this work.

La thèse présente la modélisation, l'étude expérimentale et le contrôle actif des instabilités aéroacoustiques rencontrées en cavité sous écoulement turbulent à faible nombre de Mach. On propose une formulation problème de stabilité de l'interface fluide séparant deux écoulements uniformes de vitesse différente en integrant les effets acoustiques. Les modes d'instabilité de l'interface sont alors étudiés en fonction du nombre de Mach et de la configuration géométrique. Une maquette comportant une cavité de dimensions réglables est ensuite étudiée en soufflerie à l'aide de mesures de pression. Ces données valident en partie l'approche analytique adoptée. On conçoit alors un dispositif de contrôle des modes d'instabilité, appliqué en particulier dans le cas de leur couplage avec l'acoustique de la veine d'essais. Enfin, un système de vélocimétrie par images de particules synchronisé sur les modes d'oscillation permet de valider l'étude théorique et la stratégie de contrôle.

The need to protect marine structures from the high-intensity impulsive loads created by underwater explosions has stimulated renewed interest in the mechanical response of sandwich structures. The objective of this combined numerical and experimental study is to analyze the dynamic response of composite sandwich structures and develop material-structure-property relations and design criteria for improving the blast-resistance of marine structures. Configurations analyzed include polymer foam core structures with planar geometries. A novel experimental facility to generate high-intensity underwater impulsive loads and carry out in-situ measurements of dynamic deformations in marine structures is developed. Experiments are supported by fully dynamic finite-element simulations which account for the effects of fluid-structure interaction, and the constitutive and damage response of E-glass/polyester composites and PVC foams. Results indicate that the core-density has a significant influence on dynamic deformations and failure modes. Polymeric foams experience considerable rate-effects and exhibit extensive shear cracking and collapse under high-magnitude multi-axial underwater impulsive loads. In structures with identical masses, low-density foam cores consistently outperform high-density foam cores, undergoing lesser deflections and transmitting smaller impulses. Calculations reveal a significant difference between the response of air-backed and water-backed structures. Water-backed structures undergo much greater damage and consequently need to absorb a much larger amount of energy than air-backed structures. The impulses transmitted through water-backed structures have significant implications for structural design. The thickness of the facesheets is varied under the conditions of constant material properties and core dimensions. The results reveal an optimal thickness of the facesheets which maximizes energy absorption in the core and minimizes the overall deflection of the structure.

Στην εργασία αυτή διερευνάται η δυναμική απόκριση κρυογενικών δεξαμενών υγροποιημένου φυσικού αερίου υποκείμενων σε δράση σεισμού. Επειδή οι εν λόγω δεξαμενές αποτελούν κρίσιμα στοιχεία για τη λειτουργία ενός συστήματος διανομής φυσικού αερίου και επειδή ενδεχόμενη καταστροφική αστοχία τους θα είχε δυσβάστακτες κοινωνικές επιπτώσεις, ιδιαίτερα αυστηρές απαιτήσεις ισχύουν σχετικά με τον αντισεισμικό σχεδιασμό τους. Για το σχεδιασμό έναντι σεισμού προδιαγράφονται συνήθως ένα σεισμικό γεγονός με μέση περίοδο επαναφοράς της τάξης των 500 ετών, κατά το οποίο οι δεξαμενές απαιτείται να παραμένουν πλήρως λειτουργικές, και ένα σεισμικό γεγονός με μέση περίοδο επαναφοράς της τάξης των 5000 έως 10000 ετών, κατά το οποίο απαιτείται η εξασφάλιση της ασφαλούς διακοπής της λειτουργίας των δεξαμενών, ενώ ελάχιστη μόνο βλάβη είναι αποδεκτή στα δομικά στοιχεία τους. Οι εν λόγω δεξαμενές αποτελούνται από δύο κελύφη: Ένα εξωτερικό κυλινδρικό κέλυφος από προεντεταμένο σκυρόδεμα με θολωτή στέγη από οπλισμένο σκυρόδεμα και ένα εσωτερικό κέλυφος από κρυογενικό χάλυβα, εντός του οποίου αποθηκεύεται το προϊόν. Τα δύο κελύφη εδράζονται σε κοινή πλάκα βάσης. Στην εργασία αυτή εξετάζονται δύο εναλλακτικοί τρόποι έδρασης της πλάκας βάσης: α) απλή έδραση επί του διαμορφωμένου εδάφους και β) εφαρμογή κάποιας διάταξης σεισμικής μόνωσης επί της οποίας εδράζεται η πλάκα βάσης. Ειδικότερα, εξετάζονται τα ακόλουθα συστήματα σεισμικής μόνωσης: α) μόνωση με ελαστομεταλλικά εφέδρανα με ελαστομερές υψηλής απόσβεσης, β) μόνωση με ελαστομεταλλικά εφέδρανα με πυρήνα μόλυβδου, και γ) μόνωση με συνδυασμό γραμμικών ελαστομεταλλικών εφεδράνων και μη γραμμικών αποσβεστήρων. Στις ιδιαιτερότητες του δυναμικού προβλήματος περιλαμβάνονται, εκτός από την επίδραση της μόνωσης, η δυναμική αλληλεπίδραση των δύο κελυφών, η δυναμική αλληλεπίδραση υγρού-δεξαμενής, ο κυματισμός της ελεύθερης επιφάνειας του υγρού και η δυναμική αλληλεπίδραση εδάφους-δεξαμενής. Τα φαινόμενα αυτά λαμβάνονται υπόψη στην ανάλυση. Για τη μελέτη του κυματισμού της ελεύθερης επιφάνειας χρησιμοποιείται η γραμμική θεωρία κυμάτων. Για τη μελέτη της αλληλεπίδρασης υγρού-δεξαμενής εφαρμόζεται μια προσέγγιση τύπου Lagrange με χρήση πεπερασμένων στοιχείων υγρού. Για τη μελέτη της αλληλεπίδρασης εδάφους-δεξαμενής χρησιμοποιούνται συγκεντρωμένα στοιχεία δυσπαραμορφωσιμότητας και απόσβεσης. Επίσης, συγκεντρωμένα στοιχεία δυσπαραμορφωσιμότητας και απόσβεσης χρησιμοποιούνται για την προσομοίωση των συστημάτων μόνωσης. Τα χαρακτηριστικά των παραπάνω συγκεντρωμένων στοιχείων για την προσομοίωση της αλληλεπίδρασης εδάφους-δεξαμενής και των συστημάτων μόνωσης επιλέγονται με βάση ρεαλιστικές υποθέσεις για τις δυναμικές ιδιότητες του εδάφους θεμελίωσης και τα χαρακτηριστικά των συστημάτων μόνωσης. Δύο περιπτώσεις υφιστάμενων δεξαμενών υιοθετούνται προς ανάλυση. Οι δύο δεξαμενές διαφέρουν κυρίως ως προς τη χωρητικότητα και το λόγο ύψους προς ακτίνα. Και οι δύο δεξαμενές έχουν μελετηθεί και κατασκευαστεί χωρίς σεισμική μόνωση. Ένα συνθετικό επιταχυνσιογράφημα συμβατό με το φάσμα απόκρισης σχεδιασμού του Ευροκώδικα 8 για λόγο απόσβεσης 5%, για την αντίστοιχη κατηγορία εδάφους και για μέγιστη εδαφική επιτάχυνση λίγο μεγαλύτερη από αυτή για την οποία είναι σχεδιασμένες οι δεξαμενές περιγράφει τη διέγερση της κατασκευής στην οριζόντια διεύθυνση. Για την ανάλυση αναπτύσσεται ένα συνολικό παραμετρικό προσομοίωμα πεπερασμένων στοιχείων το οποίο περιλαμβάνει τα δύο κελύφη και το αποθηκευμένο υγρό. Στο προσομοίωμα συμπεριλαμβάνονται, ανάλογα με την περίπτωση που κάθε φορά αναλύεται, τα κατάλληλα στοιχεία για την προσομοίωση της αλληλεπίδρασης εδάφους-δεξαμενής ή των διαφόρων συστημάτων μόνωσης. Το προσομοίωμα χρησιμοποιείται για τον υπολογισμό των ιδιομορφών και την πραγματοποίηση δυναμικών γραμμικών και μη γραμμικών αναλύσεων με ολοκλήρωση στο πεδίο του χρόνου. Οι μη γραμμικές αναλύσεις είναι αναγκαίες για να ληφθεί υπόψη με ρεαλιστικό τρόπο η έντονα μη γραμμική συμπεριφορά των συστημάτων μόνωσης. Στην περίπτωση αυτή η μη γραμμικότητα είναι συγκεντρωμένη στα στοιχεία που προσομοιώνουν τη μόνωση ενώ για την υπόλοιπη κατασκευή γίνεται σε κάθε περίπτωση η παραδοχή γραμμικής συμπεριφοράς. Η παραδοχή αυτή προκύπτει από το γεγονός ότι για τις δεξαμενές αυτού του τύπου η αποδεκτή βλάβη, και επομένως μη γραμμικότητα στη συμπεριφορά, των δομικών στοιχείων είναι ελάχιστη. Παρουσιάζονται επιλεγμένα αποτελέσματα των αναλύσεων που πραγματοποιήθηκαν με χρήση του προσομοιώματος πεπερασμένων στοιχείων. Τα αποτελέσματα αυτά αναφέρονται κυρίως στην τέμνουσα βάσης, τη ροπή ανατροπής, τις οριζόντιες μετακινήσεις της κατασκευής σε διάφορες στάθμες και το ύψος του κυματισμού. Έμφαση δίνεται στην ποσοτικοποίηση της επίδρασης της αλληλεπίδρασης εδάφους-δεξαμενής και της σεισμική μόνωσης στα ανωτέρω μεγέθη. Η εργασία τελειώνει με γενικά συμπεράσματα. / In the present work the dynamic response of liquefied natural gas cryogenic tanks subjected to earthquake action is investigated. Since these tanks are critical elements for the function of a natural gas distribution system and because a potential failure of them could lead to a major disaster, very severe requirements concerning the seismic design of these tanks are imposed. For the design against earthquake action two seismic events are generally considered: an event with a mean return period of 475 years during which the tanks are expected to remain fully operational and an event with a mean return period of the order of 5000 to 1000 years during which the safe shut down of the tanks is to be ensured, while minimum damage is accepted at the structural parts of the tanks. The type of tanks under consideration is constituted of two shells: an outer one made of prestressed concrete and an inner one made of cryogenic steel in which the product is stored. Both shells rest on a common base slab. In these work two alternative ways for the support of this slab are examined: a) the slab lays directly on the ground and b) the slab lays on a number of devices which provide seismic isolation. The following isolation systems are examined. a) high damping rubber bearings, b) lead core rubber bearings and c) low damping rubber bearing in conjunction with non-linear viscous dampers. The examined structural dynamic problem is characterised by certain particularities, besides the one of the implementation of an isolation system. The most important of them are sloshing of the free surface, dynamic fluid-structure interaction, dynamic soil-structure interaction and interaction between the inner and outer shell. These phenomena are taken into consideration in the preformed analyses. For the modeling of the free surface effect the realistic assumption of small wave height is made and the linear wave theory is employed. For the modeling of the fluid-structure interaction a Lagrangian approach is applied using finite elements for the modeling of the liquid and the solid domain. The soil-structure interaction is simulated by using concentrated stiffness and damping elements. Concentrated stiffness and damping elements are used for the modeling of the isolation systems as well. Two cases of existing tanks are adopted for analysis. The two tanks differ mainly in capacity and in the height to radius ratio. Both tanks are designed and constructed without seismic isolation. A global parametric finite element model is developed for the analyses. In this model the appropriate elements are incorporated for the modeling of the soil-structure interaction and the seismic isolation system, depending on the analysed case. The model is used for the calculation of eigenmodes and eigenfrequencies and for performing linear and non-linear transient analyses in time domain. Non-linear analyses are necessary in order for the highly non-linear behavior of the isolation devices to be properly simulated. In these cases the non-linearity is concentrated at the elements which simulate the seismic isolation system while the rest of the structure is considered elastic in all cases. This consideration is justified by the fact that for these tanks the acceptable damage, and consequently non-linearity in the behavior of the structural part, is minimum. Selective results of the performed analyses are presented. These results refer mainly to base shear forces, overturning moments, horizontal displacements at different levels of the tank and sloshing heights. Emphasis is on the quantification of the impact of the soil-structure interaction and the implementation of the examined seismic isolation systems.