La conception parasismique moderne est toujours en évolution afin de trouver des solutions pour améliorer la performance des structures, pendant des évènements tels que les séismes, tout en demeurant une option économique. Les diagonales ductiles confinées (DDC) sont devenus un choix intéressant comme système de résistance aux charges sismiques, en combinant la grande rigidité latérale des cadres à treillis concentriques avec une très bonne capacité de dissipation d’énergie. Toutefois, ce système est enclin à de grandes flèches latérales permanentes, causées par une faible rigidité suivant la plastification de la diagonale, ce qui pourrait probablement être résolu par l’utilisation d’un matériau possédant une bonne capacité d’écrouissage comme matériau pour le noyau ductile de la DDC. Dans cette étude, des éprouvettes machinées à partir de plaques d’acier inoxydable austénitique 304L et d’acier au carbone 350WT ont été testées sous des chargements de traction ainsi que des chargements cycliques à amplitude de déformation constante et à amplitude de déformation variable, dans le but de caractériser et de déterminer la faisabilité pour l’introduction de l’acier inoxydable 304L comme matériau de noyau ductile pour les DDC. Les résultats des essais de traction ont démontré que l’acier 304L présente un écrouissage et une ductilité plus importante que l’acier au carbone. De plus, une relation décrivant le comportement en contrainte-déformation vraie a été calée basée sur l’analyse d’image. Les résultats des essais cycliques à amplitude constante ont démontré que l’acier inoxydable 304L présente un plus grand durcissement cyclique que l’acier au carbone avec des contraintes maximales cycliques atteignant jusqu’à près de 3 fois la valeur de la limite élastique. En contrepartie, l’acier au carbone a démontré une meilleure durée de vie en fatigue oligocyclique que l’acier inoxydable. Les résultats des essais cycliques à amplitude variable ont, quant à eux, été utilisés pour caler un modèle d’écrouissage mixte représentant chacun des matériaux afin d’effectuer une simulation numérique d’un test de qualification pour DDC en utilisant le logiciel Abaqus. / Contemporary seismic design is constantly evolving in search of innovative and economic solutions to enhance the performance of structures during earthquakes. Buckling-restrained braces (BRBs) have evolved as an effective seismic load-resisting system, combining high lateral rigidity of concentrically braced frames with high energy dissipation capacity. However, BRBs are prone to large residual drifts, resulting from a low-post yield stiffness of the brace component. This could be resolved by using a material with a high strain hardening capacity, as the core plate material of the BRB or one of the core plates in a multi-plate configuration. This research aims at characterizing the austenetic stainless steel 304L in order to examine the feasibility of introducing this material as a core plate material oin BRBs. In this study, coupons machined from stainless steel 304L and carbon steel 350WT plates were tested under monotonic tensile loading and cyclic loading at both constant and variable cyclic strain. Results from the monotonic tensile tests showed that the 304L steel exhibits higher strain hardening and higher ductility character than the carbon steel. From these tests, a relation describing the true stress-strain behavior of the 304L steel was calibrated using image analysis techniques. Results from the constant amplitude cyclic tests showed that 304L stainless steel exhibits higher cyclic hardening with maximum cyclic stress values up to about three times that of the yield stress. On the other hand, the carbon steel showed higher low-cycle fatigue life. Results from the variable amplitude cyclic tests were used to calibrate the nonlinear kinematic-isotropic hardening model for both materials for use in a numerical simulation of the novel BRB in seismic vulnerability assessment studies.

En hydrologie, la simulation de la transformation de la pluie en débit dans les rivières constitue un axe de recherche dynamique. À la mise en oeuvre des nouveaux modèles, il faut ajouter les tentatives d'améliorer ceux existant grâce à la possibilité qu'offrent des nouveaux outils d'acquisition de données et à la puissance de calcul des ordinateurs toujours croissante. La puissance de calcul des nouveaux ordinateurs rend utilisable des algorithmes autrefois difficiles à mettre en oeuvre comme les réseaux de neurones (RN). Les réseaux de neurones ont connu un essor dans la modélisation hydrologique dans les années 1990 où ils ont été essentiellement utilisés dans la mise en oeuvre de modèles pluie-débit. Dans cette thèse l'on cherche à remplacer le module BV3C (bilan vertical 3 couches) du modèle distribué HYDROTEL par un ensemble de réseaux de neurones. BV3C divise le sol en trois couches où il simule les teneurs en eau et les débits sortant de chacune des couches. Cette démarche a pour but d'explorer l'opportunité de remplacer des modules de modèles complexes par des réseaux de neurones qui, une fois optimisés, constituent des outils de calculs très simples, rapides et transportables sur des supports informatiques simples. Le défi d'une telle démarche est de trouver une base de données représentative susceptible d'être utilisée par le module substitué. Dans le cas présent, des données provenant de zones hydro-climatologiques différentes ont été utilisées. Ces données sont utilisées comme entrées du module original extrait de l'ensemble du modèle HYDROTEL. Les résultats de simulation sont classés avant d'être utilisés en partie pour l'optimisation et le test des réseaux de neurones. Les réseaux mis en oeuvre sont testés sur une autre partie des données et dans un cadre opérationnel où les réseaux de neurones sont réintégrés dans le modèle. Les résultats des différents tests montrent tout d'abord que la substitution donne des résultats satisfaisants sur l'ensemble des données qui n'ont servi ni à l'optimisation, ni aux tests des réseaux de neurones. En plus, on enregistre un léger gain de temps. Les résultats sur les teneurs en eau sont nettement meilleurs. Cela s'explique par le fait que celles-ci connaissent de faibles variations dans le temps. Les variations plus importantes des débits des différentes couches rendent plus difficile leur modélisation mais les résultats obtenus rendent la substitution envisageable aussi bien dans le présent cas que dans des modules plus complexes.

Une configuration alternative de la chaîne de modélisation hydroclimatique est proposée contournant la rareté des observations météorologiques nécessaires au post-traitement statistique des simulations climatiques ainsi qu’au calage du modèle hydrologique. Cette configuration propose de calibrer les fonctions de transfert nécessaires au post-traitement statistique conjointement aux paramètres du modèle hydrologique et repose sur l’hypothèse que les variations de débits observées à l’exutoire d’un bassin versant constituent un proxy valable des forçages climatiques correspondant. Cette configuration repose également sur la notion de « fonction-objectif asynchrone », c’est-à-dire un critère d’erreur négligeant volontairement la corrélation entre deux variables. Parce qu’elles sont largement dépendantes des observations disponibles, les configurations conventionnelles de la chaîne de modélisation hydroclimatique sont généralement limitées au déploiement des champs précipitations et températures. Elles sont conséquemment limitées à l’usage de formulations empiriques pour évaluer les processus hydrologiques à l’échelle du bassin versant. Par la mise en place de la configuration alternative, les champs humidité de l’air, rayonnement solaire et vitesse du vent sont intégrés à une chaîne complète de modélisation hydroclimatique, permettant la construction de projections hydrologiques à partir de la formulation physique d’évapotranspiration de référence de Penman-Montheith. Les projections ainsi produites sont comparées à celles issues d’une configuration conventionnelle employant des réanalyses comme substituts aux observations. Les deux configurations sont mises en place sur le bassin versant de la rivière Du Loup, un bassin forestier de taille intermédiaire localisé dans la Vallée du Saint-Laurent, à partir des simulations issues de l’ensemble climatique CRCM5-LE et du modèle hydrologique distribué à base physique WaSiM-ETH. L’analyse des projections hydrologiques produites démontrent la capacité de la configuration alternative à produire une réponse hydrologique cohérente, généralement plus robuste que celle produite par une configuration conventionnelle. Cette dernière s’avère sensible aux biais inscrits à même les réanalyses, affectant le post-traitement des séries simulées et l’identification des paramètres calibrés. Sur la période historique, la formulation de Penman-Montheith offre une représentation plus cohérente des fluctuations d’évapotranspiration comparativement à une formulation empirique basée exclusivement sur la température de l’air. La projection de l’évapotranspiration, des teneurs en eau dans le sol ainsi que des étiages estivaux se sont avérée sensibles aux choix de la formulation d’évapotranspiration. La configuration alternative proposée répond d’abord à une incohérence structurelle inscrite dans les configurations conventionnelles de la chaîne de modélisation hydroclimatique. Elle permet également le déploiement d’une modélisation hydroclimatique à base physique cohérente pour des régions où les observations météorologiques sont rares tel que les régions montagneuses ou septentrionales, ou bien en voie de développement.

L’objectif de la thèse est de développer un modèle de prédiction de l’uni des chaussées flexibles. La réalisation de cet objectif a impliqué quatre principales étapes. Premièrement, une revue de connaissance relative à la thématique a été faite. Il ressort de cette étude que l’Indice de Rugosité Internationale (IRI) est l’indice qui représente le mieux le phénomène d’uni. Cependant, la prédiction de celui-ci est encore essentiellement faite par des modèles empiriques. Deuxièmement, les deux principales dégradations des chaussées qui peuvent être les sources de problème d’uni pendant leur durée de vie utile sont notamment la déformation permanente différentielle (ornière) et le soulèvement différentiel. Ils ont été analysés et la nécessité de leur prédiction a été établie. Les modèles de prédiction recherchés doivent intégrer dans leur structure mathématique les paramètres géotechniques de sol qui sont les sources de la variabilité le long d’un tracé. Pour le soulèvement au gel, le modèle de Konrad (1981 et 2005) a été jugé approprié pour être intégré dans le modèle de prédiction de l’uni. Cependant, pour la déformation permanente, un nouveau modèle a été développé pour les besoins du projet. Troisièmement, un modèle mathématique de prédiction de l’uni des chaussées flexibles en termes d’IRI a été dérivé. Il a été bâti sur les hypothèses rigoureuses et comporte les paramètres complexes suivants : la déformation permanente ( ε p ), le soulèvement au gel (S), la fissuration transversale ( Ft ) et le coefficient représentant l’effet d’amplification ou d’atténuation des longueurs d’onde (K). Quatrièmement, le modèle dérivé a été étalonné et validé avec de données obtenues sur des chaussées réelles. Pour l’étalonnage, les données des bases de données du LTTP et du MTQ, d’une part, et, d’autre part, les données de la cueillette des paramètres géotechniques issus des essais de caractérisations des échantillons prélevés à chaque 5 m le long de 5 différents sites au Québec, ont été utilisées. Les coefficients d’étalonnage requis ont été déterminés avec succès. Pour la validation du modèle, par manque de paramètres géotechniques à chaque 5 m le long des sections identifiées pour réaliser cet objectif, le principe de niveau d’utilisation du modèle d’uni a été développé. Trois niveaux d’utilisation du modèle ont été développés et le modèle a été validé au troisième niveau d’utilisation avec un niveau de succès acceptable. / The objective of the thesis is to develop a model of prediction of the flexible pavements roughness. This goal was achieved through a rigorous process that involved four steps. First, a review of knowledge was made. It appears from this study that the International Roughness Index (IRI) is the index that better represents the roughness phenomenon, but the prediction is still empirical. Second, the two main pavement degradations that could cause roughness problems during the useful life are the permanent deformation and frost heave. After analysis of the existing models the need of the prediction model was established. The desired prediction models should have in their mathematical structure the geotechnical parameters of soil, which are the origin of variability along a road. For frost heave, the Konrad model (1981 and 2005) was found to be appropriate. However, for permanent deformation there was no adequate model. Thus, a permanent deformation model was developed. Third, a mathematical model for the prediction of pavements roughness in terms of IRI was derived. It has been built on rigorous assumptions and includes the following parameters: permanent deformation ( ε p ), frost heave (S), transverse crack ( Ft ) and coefficient representing the effect of amplification or attenuation of wavelengths (K). Fourth, the derived model was calibrated and validated with data obtained on actual pavements. For calibration, data from the LTTP and MTQ databases and the geotechnical parameters resulting from the characterization tests of the samples, taken at every 5 m along 5 different sites in Quebec Province were used. The required calibration coefficients have been determined successfully. For the validation of the model, due to lack of geotechnical parameters at every 5 m along the identified sections to achieve this objective, the principle of the model use level has been developed. Three levels of the model use were developed and the model was validated at the third level with an acceptable level of success.

L’hydrolyse de la fraction organique des matières résiduelles municipales (FOMRM) est une approche peu étudiée comme alternative au tri mécanique des traitements mécano-biologiques (TMB) appliqués aux matières résiduelles municipales (MRM) issues de la collecte à deux voies. L’objectif est d’étudier ce tri biologique pour extraire la FOMRM en un lixiviat riche en acides organiques et valorisable par biométhanisation. Les MRM sont représentées par une matrice artificielle composée de résidus de fruits et légumes et d’inertes. L’impact de plusieurs traitements (immersion, recirculation du lixiviat, pré-compostage, prétraitement thermique), seuls ou combinés, sur les performances de tri et le potentiel de biométhanisation (PBM) est étudié. Comparativement à une dégradation passive, l’immersion et le pré-compostage augmentent respectivement : 1) le tri de 1 à 5%, séparant 65% de la FOMRM en près d’une semaine (2.5 fois plus rapidement), et 2) le PBM de 10 à 65%, pouvant générer de 10 à 33 Nm³ CH₄/tonne STV.

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2010-2011 / La sédimentation des particules contenues dans les eaux de ruissellement est reconnue comme étant le phénomène le plus important pour l'amélioration de la qualité de l'eau en bassin de rétention avant son rejet au milieu récepteur. Ce projet de maîtrise vise à apporter des éléments clés dans sa compréhension. L'objectif est de mettre en place un protocole de mesure de vitesses de chute adapté aux besoins du projet. La combinaison de ces vitesses à l'évaluation de la distribution de la taille des particules fournira des informations sur l'influence des temps de rétention en bassin autant du point de vue qualité de l'eau que sur le gain en termes de décantation physique des particules. L'objectif est donc d'obtenir des résultats décrivant les vitesses de chute des particules contenues dans les eaux de ruissellement et les diamètres de ces particules. Un protocole devra aussi être élaboré afin de mesurer cette dernière caractéristique. D'un point de vue gestion du bassin, il y a aussi un objectif visant à mettre en lumière l'existence ou non de relations entre le type d'événement pluvieux et les vitesses de chute enregistrées. À l'été 2008 et 2009, des campagnes de mesures ont été menées afin de récolter des données sur la qualité des eaux de ruissellement dans une zone résidentielle typique de Québec et des eaux à l'intérieur du bassin. L'évaluation des différents protocoles existants permettant de mesurer la vitesse de sédimentation a mené à l'emploi du protocole ViCAs élaboré à Lyon en France. Les résultats obtenus à partir du ViCAs ont permis de fractionner sous différentes classes de vitesse de chute les échantillons composites récupérés à l'entrée et à la sortie du bassin de rétention étudié. Les proportions de chaque classe varient grandement d'un événement à l'autre. Leur mise en relation avec les caractéristiques des événements pluvieux ayant généré ces eaux de ruissellement a permis de conclure que les vitesses de chute sont proportionnelles à l'intensité de pluie maximale, à la hauteur de précipitation et au volume ruisselé. La durée de l'événement a aussi un effet notable, mais non linéaire. Les essais sur le terrain combinés à ceux du laboratoire permettent de croire qu'une période de décantation de deux jours et demi entraîne l'enlèvement maximal atteignable des matières en suspension. Par contre, la comparaison des mesures de distribution de taille de particules, effectuées sur les échantillons composites des eaux de ruissellement, avec les distributions de vitesses de chute n'a pu mener à des conclusions. Il ne semble pas possible de faire un lien direct entre ces deux variables.

Le Service de l’environnement de la Ville de Québec a opté pour la protection du milieu par la mise en place de projets pilotes évaluant des pratiques de gestion optimale (PGO). Dans le cas actuel, on s’interroge sur le rendement de la noue engazonnée sur le territoire de la Ville de Québec. L’étude de cas de la noue est déterminante pour l’utilisation future de cette phytotechnologie. Le suivi a mis en évidence plusieurs défauts de construction, mais surtout des problèmes de croissance des végétaux dans le chenal principal. Conséquemment après deux ans d’implantation, l’ouvrage n’était toujours pas arrivé à maturité. La cueillette de données n’a pas permise d’estimer l’efficacité de l’ouvrage en 2010 et 2011. Plusieurs questions demeuraient sans réponse. La modélisation de la noue a permis de conclure que la conception, bien qu’imparfaite, n’était pas la cause des problèmes rencontrés puisque le projet est adapté à son milieu. / Quebec City’s « Service de l’environnement » has invested in watershed protection by doing several pilot studies to evaluate promising best management practices (BMP). This study focuses on the efficiency of a grass swale located on its main drinking water intake watershed. In addition to several construction shortcomings, vegetative growth problems in the main channel were observed: the grass swale did not reach maturity after two summer seasons (2010 and 2011). A hydraulic model of the grass swale was developed to better understand the hydraulic behaviour of the swale and to evaluate the impact of various changes in design. Results show that water velocities under 0.05 m/s were optimal for treatment and that problems met were not directly attributed to the design but mostly to construction issues. Finally, analysis of field data and modelling results led to several recommendations to improve both the design and implementation of this phytotechnology.

Tableau d'honneur de la Faculté des études supérieures et postdorales, 2015-2016 / Les structures composites bois-béton sont largement utilisées aujourd’hui grâce à l’optimisation des performances et du coût de la structure qu’elles permettent. Le comportement de ces structures dépend fortement de la connexion entre les deux matériaux. La conception de ce type d’ouvrage peut rapidement devenir complexe, car il existe une multitude de connecteurs qui ont tout un chacun des comportements qui peuvent être très différents des uns des autres. Dans ce mémoire de maîtrise, un modèle numérique qui permet de calculer avec précision les fréquences naturelles d’une poutre composite est présenté et implémenté. De plus, un modèle qui permet de prédire le comportement statique d’un connecteur ponctuel et d’une poutre composite en tenant compte du comportement non linéaire de chacun des matériaux est présenté et implémenté. Ces modèles reproduisent très bien les essais expérimentaux. Avec les simulations effectuées, il a pu être conclu qu’un connecteur dimensionnable serait l’idéal pour achever une performance ductile d’une poutre composite bois-béton. L’idée de ce type de connecteur est donc explorée et testée en laboratoire. Garantir une performance ductile rend la structure plus sécuritaire et permet ainsi d’augmenter l’utilisation du bois dans le secteur non résidentiel de la construction. / The timber-concrete composite structures are widely used today because they allow the performance optimization and the structure cost. The behaviour of timber concrete composite (TCC) structures depends heavily on the connection between the two materials. The design of this type of structure can quickly become complex because there are a host of connectors whose behaviour may be very different from each other. In this master thesis, a numerical model that can accurately calculate the natural frequencies of a composite beam is introduced and implemented. In addition, a model that predicts the static behaviour of a discrete connector and of a composite beam by taking into account the non linear behaviour of each material is presented and implemented. These models reproduce very well the experimental tests.With simulations, it has been found that an engineered connector would be ideal to provide for a ductile performance of a TCC beam. The idea of this type of connector is explored and tested in a laboratory. Ensure ductile performance makes the structure more secure and thus can increase the use of wood in non-residential construction sector.

Les bâtiments actuels faits de pièces en béton préfabriqué utilisent généralement des assemblages qui n’offrent que peu ou pas de continuité à la structure. Il en résulte une sous utilisation importante des matériaux par rapport à une même structure se comportant de façon monolithique. Plusieurs études portant sur des assemblages rigides, soudés ou assemblés par post-tension, ont déjà été réalisées. Bien que structuralement viable, aucun de ces modèles ne s'est avéré capable de répondre pleinement aux besoins de l'industrie. Les tolérances de construction trop serrées et la complexité de réalisation ont souvent posé problème. Les travaux présentés ici portent donc sur le développement d'un assemblage rigide boulonné, utilisable à grande échelle, destiné à l'industrie du bâtiment multiétagé en béton préfabriqué. L'étude entreprise comporte deux phases, soit la conception du système structural et la détermination de la rigidité des assemblages. Les résultats obtenus sont utilisés afin de comparer le nouveau concept aux charpentes actuelles. La démarche d’analyse s’appuie sur le logiciel VisualDesign pour élaborer la structure du bâtiment, et sur le logiciel d’analyse par éléments finis ANSYS pour la caractérisation du comportement des assemblages. Les efforts de recherche ont été orientés de manière à établir la faisabilité du modèle, afin de permettre le développement futur de projets de recherche en partenariat avec l’industrie. / To this day, the connections of precast concrete buildings offer little or no continuity to the structure. Consequently, when compared to monolithic structures, the material is clearly not used to its full capacity. Studies on rigid connections, welded or assembled by post tension, have been carried out by a few investigators. Although structurally valid, none of the proposed connections proved able to meet the needs of the construction industry, due to stringent erection tolerances and complexity of assembly. The work presented here concerns the development of a bolted rigid connection system likely to be used on a large scale by the precast concrete construction industry. The study is separated in two parts. The first is concerned with the development and tuning of the new structural system while the second one deals with the determination of the rigidity of the proposed connection. The connection behaviour is then used to compare the amount of material used in the new system and the actual precast system. The software Visual Design was used to analyse the building structure, and the finite element software ANSYS to determine the characteristic behaviour of the new connection. The intent of the research is to establish the feasibility of the proposed model, so that future research and development efforts can be carried out in partnership with the precast industry.

L’apparition des fissures macroscopiques peut manifestement se déclencher à plusieurs étapes du processus de la vie d’une cathode dans une cuve d’électrolyse. Diverses causes peuvent engendrer un tel surgissement. Au préchauffage électrique, par exemple, la répartition spatiale irrégulière des gradients thermiques peut conduire à la naissance d’une ou plusieurs fissures. L’objectif principal de ce travail se résume à la modélisation et la simulation de la propagation de fissures dans un contexte multiphysique, thermo-électro-mécanique. L’approche se base sur une philosophie de développement numérique monolithique, du couplage fortement non linéaire avec la prise en compte des échanges thermoélectriques à travers les lèvres de fissures, par la méthode des éléments finis étendus (XFEM). Ce projet a été dressé sous forme de trois sous-projets. En premier lieu, un nouveau modèle de la propagation de fissures en contexte thermomécanique dans un domaine fissuré a été développé. Une stratégie d’enrichissement géométrique par XFEM a été élaborée avec succès pour modéliser la propagation de fissures ainsi que pour atteindre le taux optimal de convergence espéré sur les réponses physiques et sur le calcul des facteurs d’intensité de contraintes. Trois benchmarks ont été examinés et validés à partir des résultats existants dans la littérature. Ensuite, un nouvel exemple de propagation de fissures et défauts multiples a été conçu. La comparaison des résultats obtenus montre la robustesse ainsi que la précision du développement numérique. En second lieu, une approche originale de la modélisation du couplage fort thermoélectrique, en présence d’une discontinuité matérielle forte statique et en tenant compte des échanges thermique et électrique à travers ses lèvres, a été développée. Tout d’abord, un modèle a été élaboré et validé dans le cas d’une fissure capacitive, où les propriétés du diélectrique, entre les lèvres de la fissure, jouent un rôle déterminant pour la quantification de son intensité d’échange. De même, un modèle a été confirmé pour le choix d’un milieu gazeux remplissant la région entre les deux lèvres, assurant l’aspect résistif d’échange thermique. Le modèle thermoélectrique a été mis en œuvre par XFEM, en implémentant la version améliorée, et en tenant compte du traitement des parasites numériques dus aux éléments de transition. Trois cas d’études ont été proposés, le premier a été appliqué pour valider la mise en œuvre numérique du modèle proposé via un patch test développé. Une bonne convergence et précision des résultats numériques ont été observées. Le second a porté sur le cas d’une fissure courbée, moins fréquente dans la littérature, qui peut être rencontré dans de nombreuses applications, et représente un défi numérique notable. Le troisième est une étude élaborée sur la cathode. L’analyse de l’impact de la présence des fissures sur l’efficacité énergétique a été soulevée aussi pour le cas de fissures multiples. En troisième lieu, une originalité numérique de la propagation de fissures en couplage thermo-électro-mécaniques (TEM) avec les échanges thermique et électrique entre les lèvres de la fissure a été exposée. La fissure n’est pas statique, cette fois-ci. L’aire générée par le déplacement des lèvres est prise en considération dans les expressions physiques des coefficients d’échanges thermique et électrique. Cela est transcrit par la quantification du saut de déplacement à travers les lèvres. Cet aspect rend le problème mécanique non linéaire, aussi, par le biais des échanges, et par conséquent le système global TEM est fortement non linéaire. Une approche originale du point de vue de la technicité d’intégration a été implémentée. Elle est fondée sur une technique de sous-triangulation barycentrique par une ’toile d’araignée’ pour les éléments de surface coupés par la fissure. Une autre technique basée sur l’intégration par ’bras’ de fissure réservée pour les éléments d’échanges thermique et électrique. Deux cas d’études essentiels ont été envisagés. Le premier est un exemple de propagation en mode mixte d’une fissure inclinée en TEM avec les échanges thermique et électrique en fonction du saut de déplacement. Le deuxième a été réservé pour un cas d’étude industrielle d’une cuve en opération, en contexte TEM. L’impact de la présence de la fissure sur les différentes réponses physiques est analysé. Comme constatation, en expansion thermoélectrique du bloc cathodique, la fissure a plus tendance à rejoindre la surface supérieure pour mettre, probablement, fin à la vie de la cathode. / The outbreak of visible cracks can be triggered at several stages in the life of a cathode block in an electrolysis cell. Various matters can prompt such an upsurge. Under electrical pre-heating, for example, the random spatial distribution of thermal gradients can lead to the rise of one or more cracks. The main objective of this work is to model and simulate the crack propagation in a multi-physical, thermo-electro-mechanical, context. The approach is based on a monolithic numerical development philosophy of a strongly nonlinear coupling, with the consideration of thermoelectric exchanges through the crack lips using the extended finite element method (XFEM). This project was sketched essentially on three sub-projects. Firstly, a new model of crack propagation in a thermomechanical fashion over a cracked domain has been developed. A geometrical enrichment strategy by XFEM has been successfully performed to model crack growth as well as to achieve the expected optimal convergence rate in physical responses and the computation of stress intensity factors. Three benchmarks were examined and validated from existing results in the literature. Then, anew example of the propagation of multiple cracks and multiple defects was designed. The comparison of the results obtained shows a good agreement with the reported works as well as the robustness and the accuracy of the numerical development. Secondly, an original approach to the modeling of full thermoelectric coupling, in the presence of a strong static material discontinuity and taking into account the thermal and electrical exchanges through its lips, has been elaborated. First, a model has been established and validated in the case of a capacitive crack, where the properties of the dielectric, between the inter-crack-lips, play a determining role for the ascertainment of its exchange intensity. Similarly, a model has been settled for the choice of a gaseous medium filling the gap between the two rims, ensuring the resistive heat exchange aspect. The thermoelectric model has been implemented via XFEM by performing the enhanced version and taking into account the treatment of numerical noise due to the blending elements. Three case studies were intended, the first one was employed to validate the numerical implementation of the stated model via a developed patch test. Good convergence and accuracy of numerical outcomes have been noted. The second focused on the case of a curved crack, less considered in the literature, which can be encountered in many applications, and represents a significant numerical challenge. The third is an elaborate study on the cathode. The analysis of the impact of cracks on energy efficiency was also raised for the case of multiple cracks. Thirdly, numerical originality of crack propagation in the context of the full thermo-electro-mechanical (TEM) coupling combined with thermal and electrical exchanges between the crack lips has been displayed. The area formed by the displacement of the two crack bank sis carried out in the physical expressions of the heat and electrical coefficients. The quantification of this gap is rendered by the resulting displacement jump across the lips. This aspect makes the mechanical problem nonlinear as well through exchanges, and therefore the global TEM system is strongly nonlinear. An original integration strategy, from a technical point of view, has been realized. It is based on a technique of barycentric sub-triangulation through a ’spider-web’ for the surface elements cut by the crack. Another procedure based on integration by crack ’arms’ reserved for heat and electrical exchanges elements. Two case studies have been examined. The first one is an example of the mixed-mode growth of a sloped crack in TEM context with thermal and electrical exchanges as a function of the displacement jump. The second was held for an industrial case of an electrolysis cell in operation and TEM circumstances. The impact of the presence of crack on several physical responses is analyzed. As an outcome, due to the thermoelectric expansion of the cathode block, the crack is expected to join the upper surface and lead to the end of the cathode’s life.