Made available in DSpace on 2014-12-17T14:58:15Z (GMT). No. of bitstreams: 1 MayaraSCFL_DISSERT.pdf: 5165682 bytes, checksum: c5b249c3b897f27db4e517452be9b9ce (MD5) Previous issue date: 2012-04-30 / Coordena??o de Aperfei?oamento de Pessoal de N?vel Superior / Composites based on PEEK + PTFE + CARBON FIBER + Graphite (G_CFRP) has increased application in the top industries, as Aerospace, Aeronautical, Petroleum, Biomedical, Mechanical and Electronics Engineering challenges. A commercially available G_CFRP was warmed up to three different levels of thermal energy to identify the main damage mechanisms and some evidences for their intrinsic transitions. An experimental test rig for systematize a heat flux was developed in this dissertation, based on the Joule Effect. It was built using an isothermal container, an internal heat source and a real-time measurement system for test a sample by time. A standard conical-cylindrical tip was inserted into a soldering iron, commercially available and identified by three different levels of nominal electrical power, 40W (manufacturer A), 40W (manufacturer B), 100W and 150W, selected after screening tests: these power levels for the heat source, after one hour of heating and one hour of cooling in situ, carried out three different zones of degradation in the composite surface. The bench was instrumented with twelve thermocouples, a wattmeter and a video camera. The twelve specimens tested suffered different degradation mechanisms, analyzed by DSC (Differential Scanning Calorimetry) and TG (Thermogravimetry) techniques, Scanning Electron Microscopy (SEM) and Energy-Dispersive X-Rays (EDX) Analysis. Before and after each testing, it was measured the hardness of the sample by HRM (Hardness Rockwell M). Excellent correlations (R2=1) were obtained in the plots of the evaporated area after one hour of heating and one hour of cooling in situ versus (1) the respective power of heat source and (2) the central temperature of the sample. However, as resulting of the differential degradation of G_CFRP and their anisotropy, confirmed by their variable thermal properties, viscoelastic and plastic properties, there were both linear and non-linear behaviour between the temperature field and Rockwell M hardness measured in the radial and circumferential directions of the samples. Some morphological features of the damaged zones are presented and discussed, as, for example, the crazing and skeletonization mechanism of G_CFRP / Comp?sitos baseados em matrizes polim?ricas de PEEK e PTFE, refor?adas com fibra de carbono e grafite (G_CFRP) apresentam crescente aplica??o e desafios ? Engenharia nas ind?strias Aeroespacial, Aeron?utica, de Petr?leo, Biom?dica, Mec?nica e Eletr?nica. Um comp?sito G_CFRP foi aquecido em tr?s n?veis de energia t?rmica para identificar os principais mecanismos de dano e algumas evid?ncias em suas transi??es de mecanismos. Uma bancada experimental foi desenvolvida para sistematizar o fluxo t?rmico com base no Efeito Joule. Foi constru?da usando-se um recipiente isot?rmico, uma fonte quente interna e um sistema de medidas em tempo real para ensaiar um corpo-de-prova (CP) de cada vez. Uma ponta c?nica-cil?ndrica foi inserida em um ferro de soldar, comercialmente dispon?vel e identificado por tr?s diferentes n?veis de pot?ncia el?trica, 40W (fabricante A), 40W (fabricante B), 100W e 150W, selecionados ap?s ensaios piloto: estes n?veis de pot?ncia para a fonte quente, ap?s uma hora de aquecimento e uma hora de resfriamento in situ, promoveu tr?s zonas diferentes de degrada??o na superf?cie do comp?sito. A bancada foi instrumentada com doze termopares, um watt?metro e uma c?mera de v?deo. Os doze C.P. ensaiados apresentaram diferentes mecanismos de degrada??o, analisados pelas t?cnicas de Calorimetria Diferencial Explorat?ria (DSC) e Termogravimetria (TG), e pelas an?lises de Microscopia Eletr?nica de Varredura (MEV) e Energia Dispersiva de Raios-X (EDS). Antes e ap?s cada ensaio, foram feitos ensaios de dureza Rockwell M (HRM). Excelentes correla??es (R2=1) foram obtidas nas curvas da ?rea evaporada ap?s uma hora de aquecimento e uma hora de resfriamento in situ versus (1) a respectiva pot?ncia da fonte quente e (2) a temperatura central do C.P. entretanto, como resultado da degrada??o diferencial do G_CFRP e da sua anisotropia, confirmadas por suas propriedades t?rmicas vari?veis, propriedades viscoel?sticas e viscopl?sticas, houve comportamentos linear e n?o-linear entre o campo de temperatura e a HRM medidos nas dire??es radial e circunferencial dos C.P. Algumas peculiaridades morfol?gicas das zonas de dano s?o apresentadas e discutidas, como, por exemplo, os mecanismos de dano por crazing e esqueletiza??o do G_CFRP

Icing has since long been identified as a serious issue in the aeronautical world. Ice build up, due to the presence of supercooled water droplets in clouds, leads to degradation of aerodynamic and/or air intake performances, among other undesirable consequences. Hence aircraft manufacturers must comply with certifications and regulations regarding flight safety in icing conditions. In order to do so, ice protection systems are used. Due to the multi-physical context within which these systems operate, numerical simulation can be a valuable asset. The present work deals with the numerical modelling of electro-thermal ice protection systems. It is built around the development of three modules. Two of them are dedicated to modelling heat transfer in the system and in the ice block. The other one models the mechanical behaviour of ice and fracture. Hence, the mechanical properties of atmospheric ice are reviewed in order to identify some mechanical parameters relevant to the fracture model. The fracture mechanics numerical method is used to investigate possible ice shedding mechanisms, that is to say the mechanisms leading to the detachement of ice, which are not yet well understood. The final goal of this work is to propose a completely coupled 2nd generation simulation methodology for electro-thermal ice protection systems. Hence the feasibility of a coupled thermal computation with ice shedding prediction based on the developed modules is shown. / Le givrage a depuis longtemps été identifié comme une problématique sérieuse dans le monde aéronautique.L’accrétion de givre, due à la présence de gouttelettes d’eau surfondue dans les nuages, dégrade les performances aérodynamiques et le rendement des entrées d’air parmi d’autres conséquences néfastes. Ainsi, les avionneurs sont sujets à des règles de certifications concernant la capacité à voler en conditions givrantes. Pour se faire, des systèmes de protection contre le givre sont utilisés. En raison de la complexité des phénomènes physiques mis en jeux, la simulation numérique constitue un atout lors de la phase de conception. Ce travail de thèse porte sur la modélisation et la simulation numérique des systèmes électrothermiques de protection contre le givre. Il s’articule autour de trois approches de modélisation, qui ont donné lieu au développement de trois modules. Deux d’entre eux sont dédiés à la simulation du transfert de chaleur dans le système et dans la glace (changement de phase). Le troisième est lié à la modélisation du comportement mécanique du givre atmosphérique avec fissuration. Ainsi, les propriétés mécaniques du givre atmosphérique sont revues de façon à pouvoir identifier les paramètres intervenant dans le modèle de fissuration. Ce modèle est ensuite utilisé pour étudier les mécanismes possibles de détachement du givre, qui ne sont à l’heure actuelle pas encore bien compris. Le but final de ce travail est de proposer une méthodologie de simulation couplée pour les systèmes électrothermiques de protection contre le givre. Ainsi, la faisabilité d’un calcul couplé thermique-fissuration avec prédiction de détachement de givre est présentée.

Fatigue in concrete is a complex phenomenon involving formation of microcracks, their coalescence into major crack and simultaneous formation of the fracture process zone ahead of the crack tip. Complex phenomena are best dealt through an energy approach and hence it is reasonable to use the theory of thermodynamics. Fracture mechanics and damage mechanics are two theories that are based on physically sound principles and are used to describe failure processes in materials. The former deals with the study of macroscopic cracks, whereas the latter defines the state of microcracking. In this study, the concepts from these theories are utilized to improve our understanding and modeling of fatigue process in concrete. In this thesis, a closed form expression for the thermodynamic function entropy is proposed and examined for its size independency and its use as a material property to characterize failure of concrete under fatigue. In the thermodynamic formalism, dissipative phenomena are described by a dissipation potential or its dual, from which evolution laws for internal variables could be defined. In this work, closed form expressions for dual of dissipation potential are derived using concepts of dimensional analysis and self-similarity within the framework of fracture mechanics and damage mechanics. Consequently, a fatigue crack propagation law and a fatigue damage evolution law are proposed respectively. A method is proposed in this study to correlate fracture mechanics and damage mechanics theories by equating the potentials obtained in each theory. Through this equivalence, a crack could be transformed into an equivalent damage zone and vice versa. Also, damage state corresponding to a given crack in a member can be quantified in terms of a damage index. An analytical way of computing size independent S-N curves is proposed, using a nonlocal damage theory by including aggregate size and specimen size in the formulation. It is realized from this study that fracture mechanics and damage mechanics theories should be used in a unified manner in order to accurately model the process of fatigue in concrete. Furthermore, based on the models developed in this study, several damage indicators for fatigue of concrete are proposed. The advantages and limitations of each of these indices are presented such that, the relevant damage index could be used, based on available parameters. Additionally, deterministic sensitivity studies are carried out to determine the most important parameters influencing fatigue life of a concrete member.

Problems in the control and identification of structural dynamic systems can lead to multimodal optimization problems, which are difficult to solve using classical gradient based methods. In this work, optimization problems pertaining to the vibration control of smart structures and the exploration of isospectral systems are addressed. Isospectral vibrating systems have identical natural frequencies, and existence of the isospectral systems proves non-uniqueness in system identification. For the smart structure problem, the optimal location(s) of collocated actuator(s)/sensor(s) and the optimal feedback gain matrix are obtained by maximizing the energy dissipated by the feedback control system. For the isospectral system problem, both discrete and continuous systems are considered. An error function is designed to calculate the error between the spectra of two distinct structural dynamic systems. For the discrete system, the Jacobi matrix, derived from the given system, is modified and the problem is posed as an optimization problem where the objective is to minimize the non-negative error function. Isospectral spring-mass systems are obtained. For the continuous system, finite element modeling is used and an error function is designed to calculate the error between the spectra of the uniform beam and the non-uniform beam. Non-uniform cantilever beams which are isospectral to a given uniform cantilever beam are obtained by minimizing the non-negative error function. Numerical studies reveal several isospectral systems, and optimal gain matrices and sensor/actuator locations for the smart structure. New evolutionary algorithms, which do not need genetic operators such as crossover and mutation, are used for the optimization. These algorithms are: Artificial bee colony (ABC) algorithm, Glowworm swarm optimization (GSO) algorithm, Firefly algorithm (FA) and Electromagnetism inspired optimization (EIO) algorithm.

Práce se zabývá tvárným lomem, který je výsledkem víceosého kvazi-statického monotónního namáhání doprovázeného rozsáhlými plastickými deformacemi, přičemž pro degradaci materiálu je uvažován lokální přístup. Ve výpočtech o rozvoji poškození rozhodují použité mezní podmínky tvárného lomu. Tyto byly teoreticky studovány v úvodu práce a po výběru vhodné mezní podmínky byl stanoven postup kalibrace. Dále byl rozpracován plán měření a realizovány zkoušky při pokojové teplotě na slitině hliníku 2024-T351, zahrnující tah, krut a tlak, pro studium rozvoje poškození a věrohodnou kalibraci vybraného fenomenologického modelu tvárného porušování, vyjádřeného pomocí lomového přetvoření a závislého na hydrostatickém tlaku a deviátoru tenzoru napětí. Mezní podmínka tvárného lomu byla posléze svázána s podmínkou plasticity. Plasticita byla pro zkoumaný materiál uvažována ve tvaru zohledňujícím i stav třetího invariantu deviátoru tenzoru napětí. Celý navržený přístup, plně aplikovatelný na víceosé úlohy, byl implementován pomocí uživatelské rutiny do komerčního programu založeného na explicitní variantě metody konečných prvků. V závěru práce je předložena aplikace navrženého přístupu k modelování tvárného porušování v podobě verifikace na vybraných zkušebních testech, z níž plynou závěry a doporučení pro další práci.

Gesintertes Silber (SAG) stellt eines der vielversprechendsten Materialien für Hochtemperaturanwendungen in der Leistungselektronik dar. Im Vergleich zu konventionellen Loten sind die mechanischen und thermischen Vorteile enorm, allerdings hochgradig prozessabhängig. Zusammen mit den relativ zeitintensiven Ermüdungstestmethoden ist das die Ursache, dass es aktuell nur wenige Lebensdauermodelle dazu gibt. In dieser Arbeit wird am Beispiel solcher SAG-Proben ein mechanisch beschleunigter, isothermer Biegeversuch vorgestellt, welcher das Potenzial hat, die zeitkritischen Temperatur- oder Lastwechselversuche zu ersetzen. Zum Vergleich wurde ein Temperaturwechseltest als Referenzversuch durchgeführt. Hierzu wird zunächst der Stand der Technik des Silber-Sinterns aufgezeigt, wobei der Schwerpunkt auf der mechanischen Materialcharakterisierung liegt.Wo das elastische Verhalten als näherungsweise allein porositätsabhängig gelten kann, ist die inelastische Dehnung noch unzureichend untersucht. Besonders die zeitabhängige inelastische Dehnung (Kriechen) zeigt noch kein vollständig konsistentes Bild, wodurch auch die Fehlermechanismen und deren Gewichtung noch nicht grundsätzlich als geklärt gelten können. Die gängigsten Belastungstests, welche in der Literatur zu finden sind, haben schwerwiegende Nachteile. Der hohe Zeitbedarf, die teils schwer quantifizierbaren Fehlerparameter und die fehlende Einstellmöglichkeit des Verhältnisses Kriechdehnung zu plastischer Dehnung sind hier im Besonderen zu nennen. Der rein dehnungsgesteuerte Biegeversuch hat diese Nachteile nicht. Über die Biegegeschwindigkeit ließe sich der Kriechanteil nahezu beliebig erhöhen (ggf. unter Nutzung von Haltezeiten). Die Biegeversuche wurden isotherm bei fünf Temperaturen von 22◦C bis 125◦C mit je drei Amplituden und drei Biegegeschwindigkeiten durchgeführt. Schlecht gesinterte Proben machten sich reproduzierbar als Frühausfall bemerkbar, so dass sich die Methode bereits gleich zu Beginn als hervorragender Qualitätstest bewährte. In puncto Ermüdung konnte ein stabiles und reproduzierbares Ausfallverhalten in Form von vergleichbaren Weibull-Formfaktoren und Ausfallbildern in den metallografischen Schliffen gefunden werden. Mit den Daten der Biegeversuche wurde ein fehlerphysikalisches Lebensdauermodell (Coffin-Manson) aufgestellt, welches erfolgreich den Ausfall des Temperaturwechseltests als Referenzversuch vorhersagen konnte. / Sintered silver (SAG) as die attach material is one of the promising solutions to exploit the advantages of high-gap semiconductors in power electronics. The mechanical and thermal properties are far superior to solders, but severely process-dependent. Combined with the time requirements of the state of the art (SoA) fatigue test methods this is most likely the reason for the lack of profound reliability studies yet. This thesis presents an isothermal bending test, which has the ability to replace the time-consuming thermal shock test as primary fatigue experiment for physics of failure based (PoF) lifetime models. A benchmark against a conventional thermal cycling test was done. The state of the art of the silver-sintering technique will be given with focus on the mechanical material characterization. While the elastic properties are mostly porosity-dependent, the inelastic properties are insufficiently examined yet. Especially the creep does not show a consistent image, what leads to many questions regarding the failure mechanism. The most common fatigue tests in the literature do have serious disadvantages. The time-consumption is high, the failure parameter can hardly be quantified and the ratio of plasticity to creep cannot be adjusted easily. The pure mechanical bending test does not have those disadvantages. By changing the bending-speed and the addition of holding times, the creep can be adjusted almost at will. The bending-experiments were conducted at five different temperatures between 22°C and 125°C and with three bending amplitudes as well as three speeds. Insufficiently sintered samples could be identified very early. This already proofed the value of the test as a quality test. Furthermore, a stable and repeatable fatigue behaviour could be observed, what was given by stable Weibull-exponents and repeatable cross sections. A lifetime-model was established by usage of the bending-test-data, what eventually predicted successfully the lifetime of a thermal cycling reference test.

Die vorliegende Arbeit umfasst die Entwicklung und Implementierung eines nichtlokalen duktilen Schädigungsmodells, die Anwendung von Verfahren zur Bestimmung von Materialparametern aus Versuchsdaten sowie die Berechnung von bruchmechanischen Kennwerten im spröd-duktilen Übergangsbereich durch die numerische Simulation von Bruchmechanikversuchen unter Verwendung der ermittelten Parameter. Das entwickelte nichtlokale duktile Schädigungsmodell basiert auf dem Modell nach Gurson-Tvergaard-Needleman (GTN). Durch die Einführung eines zusätzlichen Längenparameters mittels einer impliziten Gradientenformulierung wird die Netzabhängigkeit des GTN-Modells eliminiert. Zur Lösung des nunmehr gekoppelten Feldproblems wird das nichtlokale Schädigungsmodell in Form eines benutzerdefinierten Elements in ein Finite-Elemente-Programm implementiert. Zur Parameterbestimmung werden Kraft-Verschiebungs-Kurven des Small-Punch-Tests (SPT), einem Kleinstprobenversuch, ausgewertet. Aufgrund des sehr geringen Materialbedarfs für die Herstellung der benötigten SPT-Proben werden Rest- bzw. Bruchstücke von Charpy-Tests weiterverwendet. Es werden zwei ferritische Reaktorstähle im bestrahlten und unbestrahlten Zustand untersucht. Die Versuchsreihen decken den vollständigen Bereich der Zähigkeit beider Stähle ab. Die Bestimmung von Bruchzähigkeiten erfolgt nach dem Konzept des Local Approach allein durch numerische Berechnung des Spannungs- und Verformungszustandes in Bruchmechanikproben. Hierbei werden die zuvor aus dem SPT bestimmten Fließkurven und Schädigungsparameter verwendet. Die berechneten Bruchzähigkeiten werden mit experimentellen Ergebnissen verglichen. / This work comprises the development and implementation of a non-local ductile damage model, the application of methods for the identification of material parameters from experimental data as well as the calculation of fracture mechanics parameters in the brittle-ductile transition zone through numerical simulations of fracture mechanical tests using the identified parameters. The developed non-local ductile damage model is based on the Gurson-Tvergaard-Needleman model (GTN). The pathological mesh sensitivity of the GTN model is eliminated by introducing an additional length parameter by means of an implicit gradient formulation. To solve the coupled field problem, the non-local damage model is implemented in a finite element program in the form of a userdefined element. Force-displacement-curves of the small punch test (SPT), a miniaturised test, are evaluated for the determination of material parameters. Given the modest material requirements for the preparation of the required samples remnants of Charpy-specimens are reused. Two ferritic reactor steels, both irradiated and unirradiated, are examined. The experiments cover the full brittle-ductile transition region of the steels. Following the concept of the Local Approach, fracture toughness values are determined by numerical calculation of the stress and deformation state in fracture mechanics specimen only. Here, the yield curves and damage parameters previously determined from the SPT are used. The calculated fracture toughness values are compared with experimental results.

Zyklische Kohäsivzonenmodelle beschreiben irreversibles Separationsverhalten und Schädigungsakkumulation unter zyklischer Belastung. In der vorliegenden Arbeit wird die Formulierung zyklischer Kohäsivzonenmodelle systematisiert und ihr Potenzial zur Simulation von Ermüdungsvorgängen analysiert. Die Kohäsivspannungs-Separations-Beziehungen werden auf Basis etablierter thermodynamischer Konzepte der Schädigungsmechanik aufgestellt. Zyklische Schädigungsakkumulation wird über die Entwicklungsgleichung der Schädigungsvariablen unter Berücksichtigung einer zustandsabhängigen Dauerfestigkeit beschrieben. Das Kohäsivzonenmodell wird erfolgreich für die Simulation von Werkstoffermüdung angewandt. Numerisch mithilfe der Methode der finiten Elemente erzeugte Rissfortschrittskurven bilden das experimentell beobachtete Ermüdungsrisswachstumsverhalten in allen Bereichen ab. Über Parameterstudien wird der Einfluss der einzelnen Modellparameter ermittelt. Darüber hinaus wird die Anwendung des zyklischen Kohäsivzonenmodells auf die Simulation von Wöhler-Versuchen vorgestellt und der Probengrößeneffekt auf das Ermüdungsverhalten untersucht. Der Zusammenhang zwischen den lokalen Beanspruchungszuständen in der Kohäsivzone und dem vorhergesagten globalen Versagensverhalten wird aufgeklärt. Die gewonnenen Erkenntnisse bilden die Grundlage für ein Konzept zur Identifikation der Kohäsivparameter, das auf der Auswertung von Wöhler- und Rissfortschrittskurven beruht. / Cyclic cohesive zone models describe irreversible separation behaviour and damage accumulation under cyclic loading. In the present thesis, the formulation of cyclic cohesive zone models is systemised and their potential to simulate fatigue processes is analysed. The relation between traction and separation is described based on established thermodynamical concepts of damage mechanics. Cyclic damage accumulation is controlled by a damage evolution equation taking into account a state-dependent endurance limit. The cohesive zone model is applied successfully to the simulation of material fatigue. Fatigue crack growth rate curves, which were obtained numerically by means of the finite element method, reproduce the experimentally observed behaviour in all stages. The influences of the particular parameters of the model are determined by parametric studies. In addition, simulations of uniaxial fatigue tests using the cyclic cohesive zone model are presented. Furthermore, the size effect on the fatigue behaviour is investigated. The relation between the local states within the cohesive zone and the predicted global failure modes is explained. These findings form the foundation for a concept of parameter identification which bases on the evaluation of Wöhler-curves and fatigue crack growth rate curves.

L'objet de ce projet de recherche est de prédire la durée de vie d'éléments mécaniques qui sont soumis à des phénomènes de fatigue cyclique. L'idée est de développer un schéma numérique novateur pour prédire la rupture de structures sous de tels chargements. Le modèle est basé sur la mécanique des milieux continus qui introduit des variables internes pour décrire l'évolution de l'endommagement. Le défi repose dans le traitement des cycles de chargement pour la prédiction de la durée de vie, particulièrement pour la prédiction de la durée de vie résiduelle de structures existantes. Les approches traditionnelles de l'analyse de la fatigue sont basées sur des méthodes phénoménologiques utilisant des relations empiriques. De telles méthodes considèrent des approximations simplificatrices et sont incapables de prendre en compte aisément des géométries ou des charges complexes associées à des problèmes d'ingénierie réels. Une approche basée sur la description de l'évolution thermodynamique d'un milieu continu est donc utilisée pour modéliser le comportement en fatigue. Cela permet de considérer efficacement des problèmes d'ingénierie complexe et la détérioration des propriétés du matériau due à la fatigue peut être quantifiée à l'aide de variables internes. Cependant, cette approche peut être numériquement coûteuse et, par conséquent, des approches numériques sophistiquées doivent être utilisées.La stratégie numérique sur laquelle ce projet est basé est singulière par rapport aux schémas incrémentaux en temps usuellement utilisés pour résoudre des problèmes élasto-(visco)plastique avec endommagement dans le cadre de la mécanique des milieux continus. Cette stratégie numérique appelée méthode LATIN (Large Time Increment method) est une méthode non-incrémentale qui recherche la solution de manière itérative sur l'ensemble du domaine spacio-temporel. Une importante innovation de la méthode LATIN est d'incorporer une stratégie de réduction de modèle adaptative pour réduire de manière très importante le coût numérique. La Décomposition Propre Généralisée (PGD) est une stratégie de réduction de modèle a priori qui sépare les quantités d'intérêt spacio-temporelles en deux composantes indépendantes, l'une dépendant du temps, l'autre de l'espace, et estime itérativement les approximations de ces deux composantes. L'utilisation de l'approche LATIN-PGD a montré son efficacité depuis des années pour résoudre des problèmes élasto-(visco)plastiques. La première partie de ce projet vise à étendre cette approche aux modèles incorporant de l'endommagement.Bien que l'utilisation de la PGD réduise les coûts numériques, le gain n'est pas suffisant pour permettre de résoudre des problèmes considérant un grand nombre de cycles de chargement, le temps de calcul peut être très conséquent, rendant les simulations de problèmes de fatigue intraitables même en utilisant les techniques LATIN-PGD. Cette limite peut être dépassée en introduisant une approche multi-échelle en temps, qui prend en compte l'évolution rapide des quantités d'intérêt lors d'un cycle et leur évolution lente au cours de l'ensemble des cycles. Une description type « éléments finis » en temps est proposée, où l'ensemble du domaine temporel est discrétisé en éléments temporels, et seulement les cycles nodaux, qui forment les limites des éléments, sont calculés en utilisant la technique LATIN-PGD. Puis, des fonctions de forme classiques sont utilisées pour interpoler les quantités d'intérêt à l'intérieur des éléments temporels. Cette stratégie LATIN-PGD à deux échelles permet de réduire le coût numérique de manière significative, et peut être utilisée pour simuler l'évolution de l'endommagement dans une structure soumise à un chargement de fatigue comportant un très grand nombre de cycles. / The motivation of the research project is to predict the life time of mechanical components that are subjected to cyclic fatigue phenomena. The idea herein is to develop an innovative numerical scheme to predict failure of structures under such loading. The model is based on classical continuum damage mechanics introducing internal variables which describe the damage evolution. The challenge lies in the treatment of large number of load cycles for the life time prediction, particularly the residual life time for existing structures.Traditional approaches for fatigue analysis are based on phenomenological methods and deal with the usage of empirical relations. Such methods consider simplistic approximations and are unable to take into account complex geometries, and complicated loadings which occur in real-life engineering problems. A thermodynamically consistent continuum-based approach is therefore used for modelling the fatigue behaviour. This allows to consider complicated geometries and loads quite efficiently and the deterioration of the material properties due to fatigue can be quantified using internal variables. However, this approach can be computationally expensive and hence sophisticated numerical frameworks should be used.The numerical strategy used in this project is different when compared to regular time incremental schemes used for solving elasto-(visco)plastic-damage problems in continuum framework. This numerical strategy is called Large Time Increment (LATIN) method, which is a non-incremental method and builds the solution iteratively for the complete space-time domain. An important feature of the LATIN method is to incorporate an on-the-fly model reduction strategy to reduce drastically the numerical cost. Proper generalised decomposition (PGD), being a priori a model reduction strategy, separates the quantities of interest with respect to space and time, and computes iteratively the spatial and temporal approximations. LATIN-PGD framework has been effectively used over the years to solve elasto-(visco)plastic problems. Herein, the first effort is to solve continuum damage problems using LATIN-PGD techniques. Although, usage of PGD reduces the numerical cost, the benefit is not enough to solve problems involving large number of load cycles and computational time can be severely high, making simulations of fatigue problems infeasible. This can be overcome by using a multi-time scale approach, that takes into account the rapid evolution of the quantities of interest within a load cycle and their slow evolution along the load cycles. A finite element like description with respect to time is proposed, where the whole time domain is discretised into time elements, and only the nodal cycles, which form the boundary of the time elements, are calculated using LATIN-PGD technique. Thereby, classical shape functions are used to interpolate within the time element. This two-scale LATIN-PGD strategy enables the reduction of the computational cost remarkably, and can be used to simulate damage evolution in a structure under fatigue loading for a very large number of cycles.

No description available.