Au cours des 50 dernières années, l’étude du phénomène de liquéfaction de sols sableux saturés ont fait l'objet de nombreuses recherches en laboratoire. La plupart de ces recherches antérieures se sont concentrées sur la liquéfaction de sables propres, en supposant que la présence de fines limite le développement de la surpression interstitielle et donc le risque de liquéfaction. Cependant, les sables sableux naturels se trouvent généralement dans la nature sous la forme d'un mélange de sable et de fines et, l’influence de ces fines sur le risque de liquéfaction de ce type de matériau n'est toujours pas claire. En effet, on trouve des résultats contradictoires, dans la littérature, sur l'effet des particules fines sur le phénomène de liquéfaction des sables. Dans ce contexte, l’objectif principal de ce travail expérimental est de clarifier et de quantifier l’influence des particules fines présentes au sein d’une matrice sableuse sur l’initiation et le développement du phénomène de liquéfaction. On s’est intéressé, en particulier, à la nature des particules fines (plastique/non-plastique), à leur proportion au sein de la matrice, à l’état de compacité de la matrice et à l’état de contrainte initial (état de consolidation). A cet effet, à l'aide d'un nouveau dispositif triaxial servo-hydraulique, une série d'essais monotones et cycliques ont été effectués afin de clarifier cet effet. Les résultats expérimentaux montrent que l'augmentation de la proportion de fines non plastiques augmente la résistance à la liquéfaction du mélange sous chargements monotone et cyclique. Par contre, cette tendance est inversée pour les mélanges qui contiennent des fines plastiques. Enfin, une comparaison a été établie entre le comportement de ces sols sous cisaillement monotone et cyclique non drainé en termes d'évaluation du déclenchement du phénomène d'instabilité et d'autres paramètres mécaniques / Since the last 50 years, the study of the phenomenon of liquefaction of saturated sandy soils has been a topic of extensive laboratory research. Most of the earlier research has focused on the liquefaction of clean sands assuming that the presence of fines resists the development of pore water pressure as well as the risk of liquefaction. However, natural sand is found in nature under the form of a mixture of sand and fines and, the influence of these fines on the liquefaction risk of this type of material is still unclear. In fact, we could find contradictory results in the literature review of the effect of fine particles on the sand liquefaction phenomenon. In this context, the main objective of this study is to clarify and quantify the influence of fine particles (plastic and non-plastic) present in a sandy matrix on the initiation and development of the liquefaction phenomenon .We’ve focused, in particular, on the nature of fine particles (plastic/non-plastic), their proportion in the matrix, the matrix compactness condition and initial stress state (state consolidation). For this purpose, with the aid of a new triaxial servo-hydraulic device, a series of monotonic and cyclic tests were done in order to clarify this effect. Experimental results show that the increase in non-plastic fines increases the resistance to liquefaction of the mixture under both monotonic and cyclic loading. However, this trend is reversed for the mixtures containing plastic fines. Finally a comparison has been established between the behavior of these soils under undrained monotonic and cyclic shearing in terms of evaluation of the initiation of instability phenomenon and other mechanical parameters

Le but de cette étude est d'étudier le renforcement des sols par inclusions rigides sous chargement complexe et cyclique. L'effet de certains paramètres liés à la définition d’un chargement complexe et cyclique sur le comportement du système est mis en évidence.Du point de vue des chargements statiques complexes, des semelles de fondation posées sur un sol compressible renforcé par inclusions rigides sans matelas soumis à des charges centrées, excentrées verticales et horizontales et à quelques cycles de charge ont été étudiées. Des approches numériques et expérimentales sont présentées. Les résultats des mesures expérimentales et numériques permettent de mettre en évidence le comportement de ces systèmes en termes de contrainte sur la tête d'inclusion et sur le sol compressible, de déplacements verticaux et latéraux de la semelle et du déplacement latéral de l'inclusion. L'efficacité de la semelle renforcée est comparée à celle d’une semelle non renforcée.Une modélisation 3D de solutions de fondations pour les éoliennes est étudiée. La combinaison d’un chargement vertical et de différents moments appliqués à la fondation est prise en compte. Le sol compressible renforcé par inclusions rigides est considéré comme une option qui est comparée à d’autres solutions plus classiques (fondation superficielle et radier sur pieux). Les résultats obtenus permettent de présenter l’impact sur le tassement du sol, la rotation de la fondation, les efforts axiaux et les moments fléchissants dans les inclusions rigides. Les résultats numériques indiquent enfin que la technique d'amélioration du sol par inclusions rigides peut être une solution appropriée pour les fondations d'éoliennes.En ce qui concerne les aspects cycliques, trois points principaux sont abordés. Dans un premier temps, la modélisation numérique d’essais en laboratoire d’un renforcement de sol par inclusions rigides soumis à des chargements monotones et des cycles limités de chargement mise en œuvre. Le modèle hypoplastique (HYP) est utilisé pour modéliser le comportement de la plate-forme de transfert de charge. Les résultats numériques sont validés à la fois par rapport aux données expérimentales et numériques de Houda (2016). L'influence des conditions aux limites et de l'état du sol compressible est mise en évidence. Les résultats numériques indiquent qu'il est possible de considérer le comportement cyclique du sol renforcé par inclusions rigides en utilisant le modèle HYP.Dans un second temps, un remblai renforcé par des inclusions rigides sous un nombre élevé de chargement cyclique est étudié. Deux niveaux de complexité différents pour le modèle constitutif (HYP et le modèle élastique linéaire parfaitement plastique avec un critère de rupture de type Mohr-Coulomb) ont été pris en compte pour étudier le comportement de la LTP et analyser le comportement cyclique du système. Le modèle HYP est proposé pour la suite des études car il permet de bien capturer la décroissance et l’accumulation des tassements avec le nombre de cycles de charge. L'effet des paramètres qui sont le nombre de cycles de charge, l'amplitude et la fréquence (induite par la vitesse du trafic) et la hauteur du remblai est également présentée.Finalement, une étude sur la réponse cyclique d'un remblai de GRPS est menée. En comparant le remblai renforcé par des géosynthétiques (GRPS) avec le remblai renforcé par inclusions (PE), le rôle du géosynthétique est mis en évidence sous des chargements statiques et cycliques. L'influence du nombre de cycles de chargement et du nombre de géosynthétiques sur l'effet de voute et les tassements cumulés est également discutée. / The aim of the study is to investigate the soil improvement by rigid inclusions under complex and cyclic loadings, and to highlight the effect of some parameters related to complex and cyclic loading on the system behavior.Concerning the static complex loading, footings over rigid inclusion-reinforced soil without mattress subjected to centered, eccentrically vertical and horizontal loads, and load cycles are first studied. Numerical and experimental approaches are presented. Monitored and numerical results permit to show the behavior of these reinforced systems in terms of stresses on the inclusion head and soft soil, vertical and lateral displacements of the footings and lateral displacement of the inclusions. The efficiency of the reinforced footing is also presented and compared to the unreinforced one.A 3D modeling of the foundation solutions for wind turbines is presented. The combination of vertical loading and different moments applied to the foundation is taken into account. The inclusion-improved soft soil under footing is considered as a foundation option and, compared to classical ones (shallow foundation and piled raft). The obtained results are illustrated concerning the ground surface settlements, the foundation rotations, the axial forces and bending moments of the reinforcements. The numerical results indicate that the soil improvement technique by rigid inclusions can be an appropriate solution for the wind turbine foundations.With regard to the cyclic aspects, three main concerns are studied. Firstly, the numerical modeling of laboratory tests on a soil improvement by rigid inclusions subjected to monotonic loading and a limited load cycles is carried out, in which the hypoplasticity (HYP) model is used to model the load transfer platform (LTP). The numerical results are validated against both the experimental data and numerical ones of Houda (2016). The influence of the boundary condition and soft soil state are figured out. The numerical results indicate that it is possible to address the cyclic behavior of the rigid inclusion-reinforced soil by using the HYP model.Secondly, a piled embankment under a high number of cyclic loadings is studied. Two different levels of complexity for the constitutive models are used (HYP and a simpler one the linear elastic perfectly plastic constitutive model with a shear criteria of Mohr-Coulomb). These models were considered to model the behavior of the LTP and analyze the cyclic behavior of the system. The HYP model is then suggested for the following studies since it can capture well the arching decrease and the cumulated settlements under the load cycles number. The effect of the parameters that are load cycles number, amplitude and frequency (induced by traffic speed), and embankment height is illustrated as well.Finally, a study on the cyclic response of a GRPS embankment is conducted. By comparing the geosynthetic-reinforced pile-supported (GRPS) embankment with the piled embankment (PE), the role of the geosynthetic is verified under static and cyclic loading aspects. The influence of the load cycles number and the geosynthetic layers number on the arching effect and cumulative settlements is shown as well.

A partir dune approche du type modélisation physique, on sintéresse à la quantification des interactions, appelées couramment « effets de groupe », qui peuvent se développer entre les micropieux rapprochés sous sollicitations monotones et cycliques verticales.<br />Dans le cas monotone, leffet de groupe est étudié à partir du comportement dune inclusion générique au sein dun groupe en interaction avec des inclusions adjacentes. Les résultats montrent un effet de groupe positif sur le frottement latéral et négatif sur la résistance en pointe. Dans le cas des chargements cycliques, deux cas sont considérés à savoir les chargements à force contrôlée destinés à établir des diagrammes de stabilité cyclique et les chargements à déplacement contrôlé pour les mécanismes de dégradation. Pour les premiers, laugmentation du nombre dinclusions montre un effet positif et pour les seconds, les résultats montrent un effet de groupe positif sur le frottement latéral et négatif sur la pointe.

Les travaux de recherche menés dans le cadre de cette thèse concernent l'étude du comportement des groupes et des réseaux de micropieux. L'interaction entre les micropieux d'un groupe, ainsi que l'effet d'inclinaison des micropieux au sein d'un réseau sont traités à travers les trois parties de cette thèse.<br />La première partie est consacrée à la définition, la classification et les domaines d'application des micropieux. Les différentes méthodes usuelles de calcul des micropieux (isolés, en groupes ou en réseaux) sont également présentées.<br />La deuxième partie présente la synthèse de divers essais expérimentaux réalisés, dans le cadre du projet national FOREVER, sur des groupes et/ou des réseaux de micropieux. Ces essais sont réalisés en vraie grandeur ou en modèle réduit (centrifugeuse, chambre d'étalonnage ou cuve expérimentale).<br />La troisième partie traite de deux approches numériques différentes pour l'étude du comportement des micropieux. La méthode des fonctions de transfert de charge est utilisée à travers le programme de calcul des groupes de pieux (GOUPEG) qui est développé pour tenir compte de l'effet de l'inclinaison des micropieux. Le cas d'un réseau élémentaire de 4 micropieux et d'un groupe équivalent est étudié. Une étude paramétrique sur l'effet de l'inclinaison des micropieux au sein d'un réseau élémentaire est présentée. Enfin, la méthode des éléments finis est utilisée (le logiciel CESAR-LCPC) pour l'étude du comportement d'un réseau élémentaire de 4 micropieux (en double chevalet). Cette dernière approche fournit des résultats intéressants sur le comportement du massif du sol sous les micropieux inclinés.

Les nids d'abeille d'aluminium combinent légèreté et grande capacité d'absorption d'énergie. Ils sont alors de plus en plus utilisés dans les secteurs du transport (automobile, aéronautique ...) pour contribuer conjointement à l'allègement structural et à la sécurité. Dans cette thèse, le comportement à l'écrasement des nids d'abeille est étudié en tenant compte de l'effet combiné de l'angle d'orientation dans le plan des cellules, de l'angle de chargement et de la vitesse de sollicitation, que la littérature ne relate pas. Un dispositif de chargement mixte compression/cisaillement est conçu pour mener l'étude expérimentale. L'analyse des résultats porte sur le pic initial d'effort, le plateau d'effort, ainsi que sur les modes de déformation. Les résultats montrent une augmentation de la résistance sous sollicitation dynamique dépendante de l'angle de chargement Ψ. Elle devient moins significative quand l'angle de chargement augmente jusqu'à atteindre un angle critique. Pour Ψ > Ψcritique, les réponses quasi-statiques sont même plus élevées que les réponses dynamiques. Une étude numérique est alors entreprise. Elle permet de comprendre ce phénomène qui est imputé aux mécanismes de déformation locaux des cellules. Les résultats numériques montrent également que l'effet de l'angle d'orientation □ dans le plan est plus prononcé sur la force tangentielle que sur la force normale, que cela influence également les modes d'effondrement et donc la réponse mécanique. Ces simulations numériques, couplées aux résultats expérimentaux, permettent alors de dissocier les composantes normale et tangentielle de la réponse des nids d'abeille et d'identifier les paramètres d'un critère macroscopique de résistance exprimé en fonction de la vitesse d'impact, de l'angle de chargement et de l'angle d'orientation dans le plan. Finalement, dans le but de réduire le coût des simulations numériques, un modèle élément fini (EF) réduit basé sur un critère de périodicité tenant compte de l'angle d'orientation dans le plan est proposé et son domaine de validité est évalué.

Les nids d'abeille d’aluminium combinent légèreté et grande capacité d’absorption d'énergie. Ils sont alors de plus en plus utilisés dans les secteurs du transport (automobile, aéronautique …) pour contribuer conjointement à l’allègement structural et à la sécurité. Dans cette thèse, le comportement à l’écrasement des nids d'abeille est étudié en tenant compte de l'effet combiné de l'angle d'orientation dans le plan des cellules, de l’angle de chargement et de la vitesse de sollicitation, que la littérature ne relate pas. Un dispositif de chargement mixte compression/cisaillement est conçu pour mener l’étude expérimentale. L’analyse des résultats porte sur le pic initial d’effort, le plateau d’effort, ainsi que sur les modes de déformation. Les résultats montrent une augmentation de la résistance sous sollicitation dynamique dépendante de l’angle de chargement Ψ. Elle devient moins significative quand l’angle de chargement augmente jusqu’à atteindre un angle critique. Pour Ψ > Ψcritique, les réponses quasi-statiques sont même plus élevées que les réponses dynamiques. Une étude numérique est alors entreprise. Elle permet de comprendre ce phénomène qui est imputé aux mécanismes de déformation locaux des cellules. Les résultats numériques montrent également que l’effet de l’angle d’orientation □ dans le plan est plus prononcé sur la force tangentielle que sur la force normale, que cela influence également les modes d’effondrement et donc la réponse mécanique. Ces simulations numériques, couplées aux résultats expérimentaux, permettent alors de dissocier les composantes normale et tangentielle de la réponse des nids d’abeille et d’identifier les paramètres d'un critère macroscopique de résistance exprimé en fonction de la vitesse d'impact, de l'angle de chargement et de l'angle d'orientation dans le plan. Finalement, dans le but de réduire le coût des simulations numériques, un modèle élément fini (EF) réduit basé sur un critère de périodicité tenant compte de l'angle d'orientation dans le plan est proposé et son domaine de validité est évalué. / Aluminium honeycombs combine lightweight with an efficient energy absorption capability (specific energy). They are widely used as crash energy absorbing and protective structures in a wide range of transport applications (automotive, aircraft …) to reduce energy consumption and greenhouse gas emission. Reducing vehicle mass has however to be done while at least maintaining the same safety levels. In this thesis, the honeycomb behaviour is investigated under mixed shear-compression loadings taking into account the combined effect of the in-plane orientation angle and the impact velocity, which has not been deeply investigated in the literature. Experimental study based on an improvement of a mixed shear-compression loading device is realised. Experimental analysis focuses on the initial peak and average crushing forces as well as the deforming pattern modes. Comparing quasi-static and dynamic results, a dynamic enhancement depending of the loading angle Ψ is observed under mixed shear-compression loading until a critical loading angle (Ψcritical). Beyond, a negative enhancement is observed. Thus, a numerical study is carried out. The negative enhancement phenomenon is attributed to the collapse mechanisms which are affected by the loading angle Ψ. Numerical results also highlight that the in-plane orientation angle has an effect on the collapse mechanisms and consequently on the mechanical response. This effect is more pronounced on the tangential force than the normal force. The combined effect of the in-plane orientation angle and the loading angle is analysed on the three identified deforming pattern modes. Combining numerical and experimental results, the average crushing normal and shear forces are dissociated. Therefore, the parameters of a macroscopic yield criterion for the mixed shear-compression honeycomb behaviour depending of the impact velocity, the loading angle and the in-plane orientation angle are identified. Finally, in order to optimise the cost in CPU-time of the numerical simulation, a reduced FE model based on the periodicity procedure taking into account the in-plane orientation angle is proposed and its validity range is evaluated.

Ce mémoire présente les résultats de calculs analytiques de prédiction du comportement de dispositifs de retenue de niveau PL-2 armé d'armature en polymères renforcés de fibres de verre (PRFV) soumis aux essais de chargement quasi statique horizontal. L'un des modèles analytiques à l'étude est le calcul de la charge ultime que peut reprendre un dispositif de retenue en fonction des dimensions de la fissure au niveau de la connexion. Le deuxième modèle analytique à l'étude est le calcul de prédiction de la relation charge-déflexion en tête de muret. La collaboration entre le Ministère des Transports du Québec (MTQ), le Ministère des Transports de l'Alberta et l'Université de Sherbrooke a permis d'effectuer des essais quasi statiques sur plusieurs configurations durant les dernières années. Les configurations de dispositifs provenant du MTQ sont des barrières de type 210 et 311 dont les murets et les dalles sont entièrement armés de barres en PRFV. Les configurations des dispositifs de l'Alberta sont quant à eux des dispositifs de type S-1642 et S-1650 dont les murets sont armés de barres en PRFV et les dalles sont armées de barres d'acier enduit d'époxy. Grâce à cette banque de données, il est possible d'effectuer l'étude des modèles sur différentes configurations et ainsi déterminer certaines de leurs limites. De plus, une modélisation par éléments finis de la configuration de dispositif de retenue de type 311 (PL-2) a été effectuée. Le dispositif a été modélisé en trois dimensions avec le logiciel ADINA. De ce modèle, une comparaison avec les données expérimentales a été effectuée afin de le valider. De plus, une étude paramétrique a été effectuée sur l'influence des résistances du béton de la dalle et du muret et du module d'élasticité des barres en polymères renforcés de fibres (PRF). En effet, des bétons de 35 MPa et 50 MPa ont été utilisés afin de représenter trois cas de combinaison de résistance en compression des bétons. Finalement, trois grades de modules élastiques ont été utilisés sur le dispositif ayant un béton de 35 MPa au niveau de la dalle et un béton de 50 MPa au niveau du muret (combinaison la plus courante) afin d'étudier l'effet du module élastique des barres.

Depuis le milieu des années 90, avec la revalorisation des déplacements urbains par la marche et le vélo, les passerelles piétonnières sont fortement demandées dans nos villes. Étant également très convoitées par les architectes et les artistes pour exprimer leurs idées, ces structures deviennent de plus en plus élancées, légères, flexibles et dotées de fréquences propres qui peuvent facilement être mises en résonance par la marche normale des piétons. Les premiers ouvrages ayant manifesté ces problèmes de façon remarquable sont la passerelle du Millénium de Londres et la passerelle de Solferino de Paris qui ont dû être fermées suite à leur inauguration en raison de vibrations latérales excessives. Ces structures ambitieuses créent donc de nouveaux défis qui ne sont pas adéquatement considérés dans plusieurs codes de pratique, et ce, en raison d'un mauvais contrôle des vibrations induites par des charges de groupes et de foules de piétons. De plus, les piétons ont tendance à se synchroniser au mouvement latéral de la structure, ce qui résulte en des mouvements forts qui compromettent grandement le confort des usagers. L'objet de ce projet de recherche vise à mieux définir et caractériser le chargement des groupes et des foules de piétons sur les passerelles piétonnières afin de permettre une meilleure compréhension et un meilleur contrôle des vibrations de ces structures dès leur phase de conception. Pour ce faire, les vibrations induites par des chargements de groupes et de foules de piétons ont été mesurées sur trois structures existantes soit deux travées indépendantes de la Passerelle de Sherbrooke et le pont du Cosmos Montréal. Les critères de confort et les niveaux de vibration prédits selon différentes normes ont aussi été calculés pour ces trois structures. Également, un programme a été créé afin de simuler des groupes ou des foules de piétons de taille variable circulant sur une passerelle piétonnière et de calculer la réponse de la structure selon différents paramètres. Des simulations numériques de chargement de 1 piéton, de groupes et de foules de piétons ont été effectuées l'aide de ce programme sur les trois structures afin de simuler les vibrations induites par les piétons. Différents modèles de chargements piétonniers ont été utilisés. Finalement, une comparaison des différents niveaux de vibration mesurés, calculés et simulés permet une caractérisation du chargement piétonnier et de l'influence du nombre de piétons dans le chargement sur la réponse de la structure, et ce, tant pour les vibrations verticales que latérales de la structure. Ce projet aura pour résultat de fournir une meilleure définition du chargement piétonnier et d'assurer ainsi un meilleur contrôle des vibrations des passerelles piétonnières dès les phases de conception de la structure. Ainsi, ce projet favorisera le développement des passerelles piétonnières et permettra aux architectes et artistes une plus grande liberté pour l'édification de ces sculptures urbaines sans compromettre le confort des usagers, tout en limitant les coûts par un meilleur contrôle des vibrations.

Les structures aéronautiques peuvent être soumises à des sollicitations sévères telles que les collisions, les impacts de volatiles, etc … Sous l’effet de ces sollicitations rapides, qui du fait de leurs temps caractéristiques très courts limitent les transferts thermiques, le matériau peut dissiper l’énergie dans des zones de déformation localisée qui peuvent conduire à une ruine prématurée de la structure. Le travail de la thèse porte sur la définition d’une méthodologie expérimentale destinée à étudier les conditions de rupture de matériaux à haute résistance à vocation aéronautique consécutives à un endommagement dynamique. Ce travail comprend : •la mise au point d’essais rapides de cisaillement ; •des observations microstructurales des matériaux avant et après sollicitation ; •la simulation numérique des essais. / Aeronautical structures may be submitted to severe loading such as collisions, bird strike, etc. Under dynamic loading, involving quasi adiabatic conditions, the material may dissipate energy within zones of localised deformation wich may lead to the premature failure of the structure. The PhD work aims at defining an experimental methodology devoted to study the conditions of fracture of aeronautical, high strength materials intervening after dynamic damage. Tasks include notably: * definition of dynamic shear tests * microstructural observation of the material before and after loading * numerical simulation * Development of fracture criterion

L'objectif de cette thèse est de démontrer comment faire une analyse théorique de fatigue et de mécanique de la rupture d'une structure éolienne à l'axe horizontal. La chaîne des calculs nécessaires pour atteindre cet objectif s'avère être particulièrement longue pour deux raisons : d'abord, la vitesse du vent varie aléatoirement avec le temps ; deuxièmement, l'amplitude de vibration du mât est amplifié en raison des ses fréquences naturelles de vibration. Un chapitre entier est consacré à la modélisation de la vitesse du vent dans l'espace et dans le temps. Ce même chapitre démontre comment synthétiser un signal aléatoire à partir d'une fonction de densité spectrale de puissance (DSP). La force axiale du rotor est le chargement le plus important sur une structure éolienne à l'axe horizontal. Cette force a un rapport non linéaire avec la vitesse du vent. Cela implique la nécessité de déterminer la DSP de la force axiale à partir de son signal, en se servant d'une technique d'estimation spectrale. La méthode Thomson Multitaper s'est avéré la plus satisfaisante pour cette application. La DSP des déplacements du mât est déterminée en associant la réceptance du système structurel avec la DSP de la force qui représente tous les chargements. Un signal de contrainte peut finalement être synthétisé à partir de sa DSP. La technique de comptage de cycles de chargement connue sous le nom de rainflow est abordée et appliquée. Le fait que le signal de contraintes a une amplitude variable implique la nécessité d'employer une technique plus avancée de simulation de propagation de fissures. La technique choisie pour cette thèse est connue sous le nom de strip-yield (bande d'écoulement). / The objective of this thesis is to demonstrate how to do theoretical analyses of fatigue and fracture mechanics in a structure for horizontal axis wind turbine. The chain of calculations required to reach this objective is particularly long for two reasons : firstly, the wind speed varies randomly with time , secondly, the vibration amplitude of the mast is amplified due to its natural frequencies of vibration. A whole chapter is dedicated to modeling the wind speed in space and time. The same chapter shows how to synthesize a random signal by employing a power spectral density function (PSD). The axial force of the rotor is the most important loading on a structure for horizontal axis wind turbine. This force has a non linear relation with the wind speed. This implies the need to determine the PSD of the axial force from its signal, by employing a spectral estimation method. The Thomson Multitaper method revealed to be the most satisfactory for this application. The PSD of displacement of the mast is determined by associating the receptance of the structural system and the PSD of the force representing all loadings. Finally, a signal of stress can be synthesized from its PSD. The fatigue cycle counting method known as rainflow is discussed and employed. The fact that the signal of stress has a variable amplitude implies the need of a more sophisticated method to simulate a crack propagation. The method chosen in this thesis is called strip-yield.