The modified mapping-collocation (MMC) method is applied to analyze a circumferential crack in an isotropic curved beam. Based on the method a MATLAB code was developed to obtain the stress field. Incorporating the stress correlation technique, the opening and sliding fracture mode stress intensity factors (SIF)s of the crack for the beam under pure bending moment load case are calculated. The MMC method aims to solve two-dimensional problems of linear elastic fracture mechanics (LEFM) by combining the power of analytic tools of complex analysis (Muskhelishvili formulation, conformal mapping, and extension arguments) with simplicity of applying the boundary collocation method as a numerical solution approach. Qualitatively, a good agreement between the computed stress contours and the fringe shapes obtained from the photoelastic experiment on a plexiglass specimen is observed. Quantitatively, the results are compared with that of ANSYS finite element analysis software. The effect of crack size, crack position and beam thickness variation on SIF values and mode mixity is investigated.

Probabilistic assessment and life cycle management of engineering components and systems in a nuclear power plant is intended to ensure safe and efficient operation of energy generation over its entire life. The CANDU reactor core consists of 380-480 pressure tubes, which are like miniature pressure vessels that contain natural uranium fuel. Pressure tubes operate under severe temperature and radiation conditions, which result in degradation with ageing. Presence of flaws in a pressure tube makes it vulnerable to delayed hydride cracking (DHC), which may lead to rupture or break-before-leak situation. Therefore, assessment of flaws in the pressure tubes is considered an integral part of a reactor core assessment program. The main objective of the thesis is to develop advanced probabilistic and mechanical stress field models for the assessment of flaws. The flaw assessment models used by the industries are based on deterministic upper/lower bound values for the variables and they ignore uncertainties associated with system parameters. In this thesis, explicit limit state equations are formulated and first order reliability method is employed for reliability computation, which is more efficient than simulation-based methods. A semi-probabilistic approach is adopted to develop an assessment model, which consists of a mechanics-based condition (or equation) involving partial factors that are calibrated to a specified reliability level. This approach is applied to develop models for DHC initiation and leak-before-break assessments. A novel feature of the proposed method is that it bridges the gap between a simple deterministic analysis and complex simulations, and it is amenable to practical applications. The nuclear power plant systems are not easily accessible for inspection and data collection due to exposure to high radiation. For this reason, small samples of pressure tubes are inspected at periodic intervals and small sample of data so collected are used as input to probabilistic analysis. The pressure tube flaw assessment is therefore confounded by large sampling uncertainties. Therefore, determination of adequate sample size is an important issue. In this thesis, a risk informed approach is proposed to define sample size requirement for flaw assessment. Notch-tip stress field is a key factor in any flaw assessment model. Traditionally, linear elastic fracture mechanics (LEFM) and its extension, serves the basis for determination of notch-tip stress field for elastic and elastic-perfectly-plastic material, respectively. However, the LEFM solution is based on small deformation theory and fixed crack geometry, which leads to singular stress and strain field at the crack-tip. The thesis presents new models for notch and crack induced stress fields based on the deformed geometry. In contrast with the classical solution based on small deformation theory, the proposed model uses the Cauchy's stress definition and boundary conditions which are coupled with the deformed geometry. This formulation also incorporates the rotation near the crack-tip, which leads to blunting and displacement of the crack-tip. The solution obtained based on the final deformed configuration yields a non-singular stress field at the crack-tip and a non-linear variation of stress concentration factor for both elastic and elastic-perfectly-plastic material. The proposed stress field formulation approach is applied to formulate an analytical model for estimating the threshold stress intensity factor (KIH) for DHC initiation. The analytical approach provides a relationship between KIH and temperature that is consistent with experimental results.

Probabilistic assessment and life cycle management of engineering components and systems in a nuclear power plant is intended to ensure safe and efficient operation of energy generation over its entire life. The CANDU reactor core consists of 380-480 pressure tubes, which are like miniature pressure vessels that contain natural uranium fuel. Pressure tubes operate under severe temperature and radiation conditions, which result in degradation with ageing. Presence of flaws in a pressure tube makes it vulnerable to delayed hydride cracking (DHC), which may lead to rupture or break-before-leak situation. Therefore, assessment of flaws in the pressure tubes is considered an integral part of a reactor core assessment program. The main objective of the thesis is to develop advanced probabilistic and mechanical stress field models for the assessment of flaws. The flaw assessment models used by the industries are based on deterministic upper/lower bound values for the variables and they ignore uncertainties associated with system parameters. In this thesis, explicit limit state equations are formulated and first order reliability method is employed for reliability computation, which is more efficient than simulation-based methods. A semi-probabilistic approach is adopted to develop an assessment model, which consists of a mechanics-based condition (or equation) involving partial factors that are calibrated to a specified reliability level. This approach is applied to develop models for DHC initiation and leak-before-break assessments. A novel feature of the proposed method is that it bridges the gap between a simple deterministic analysis and complex simulations, and it is amenable to practical applications. The nuclear power plant systems are not easily accessible for inspection and data collection due to exposure to high radiation. For this reason, small samples of pressure tubes are inspected at periodic intervals and small sample of data so collected are used as input to probabilistic analysis. The pressure tube flaw assessment is therefore confounded by large sampling uncertainties. Therefore, determination of adequate sample size is an important issue. In this thesis, a risk informed approach is proposed to define sample size requirement for flaw assessment. Notch-tip stress field is a key factor in any flaw assessment model. Traditionally, linear elastic fracture mechanics (LEFM) and its extension, serves the basis for determination of notch-tip stress field for elastic and elastic-perfectly-plastic material, respectively. However, the LEFM solution is based on small deformation theory and fixed crack geometry, which leads to singular stress and strain field at the crack-tip. The thesis presents new models for notch and crack induced stress fields based on the deformed geometry. In contrast with the classical solution based on small deformation theory, the proposed model uses the Cauchy's stress definition and boundary conditions which are coupled with the deformed geometry. This formulation also incorporates the rotation near the crack-tip, which leads to blunting and displacement of the crack-tip. The solution obtained based on the final deformed configuration yields a non-singular stress field at the crack-tip and a non-linear variation of stress concentration factor for both elastic and elastic-perfectly-plastic material. The proposed stress field formulation approach is applied to formulate an analytical model for estimating the threshold stress intensity factor (KIH) for DHC initiation. The analytical approach provides a relationship between KIH and temperature that is consistent with experimental results.

O trabalho trata da utilização de formulações não-convencionais de elementos finitos na obtenção de fatores de intensidade de tensão associados a múltiplas fissuras distribuídas num domínio bidimensional. A formulação do problema de múltiplas fissuras baseia-se numa abordagem de sobreposição proposta pelo Método da Partição (\"Splitting Method\"). Segundo essa abordagem a solução do problema pode ser encontrada a partir da sobreposição de três subproblemas combinados de tal forma que o fluxo de tensão resultante nas faces das fissuras seja nulo. O uso do Método dos Elementos Finitos (MEF) em sua forma convencional pode requerer um refinamento excessivo da rede nesse tipo de problema, aumentando o custo computacional da análise. Objetivando reduzir este custo, empregam-se duas formulações não-convencionais, de forma independente, num dos subproblemas, dito local: a formulação híbrido-Trefftz de tensão e o Método dos Elementos Finitos Generalizados (MEFG). Na formulação híbrido-Trefftz é adotado o recurso do enriquecimento seletivo mediante o refrno- p na aproximação dos campos de deslocamento no contorno do elemento. Já com relação ao MEFG, empregam-se funções polinomiais e a solução analítica da mecânica da fratura como funções enriquecedoras. Exemplos de simulação numérica são apresentados no sentido de comprovar que a utilização dessas formulações não-convencionais juntamente com o Método da Partição viabiliza a obtenção de resultados com boa aproximação com recurso a redes pouco refinadas, reduzindo significativamente o custo computacional de toda a análise. / This paper treats with the use of non-conventional finite element formulations to obtain the stress intensity factor of multiple cracks located in a two-dimensional domain. The formulation of the multiple cracks problem is based on an overlapping approach suggested by the Splitting Method. Accordingly, the solution of the problem can be achieved by dividing the problem in three steps, combined so that the resulting stress flux is zero on the cracks face. The use of the Finite Element Method (FEM) in its conventional formulation requires a mesh refinement in this kind of problem, then increasing the computational cost. Aiming to reduce this cost, two non-conventional formulations are used independently to solve the local problem: the Hybrid-Trefftz stress formulation and the Generalized Finite Elements Method (GFEM). The Hybrid-Trefftz formulation is applied with selective enrichment using p-refinement in the displacements field on the element boundaries. The GFEM employs polynomial functions and analytical solutions of the fracture mechanics as enrichment functions. Examples of numerical simulations are presented in order to show that non- conventional formulations and the Splitting Method can provide accurate results with coarse mesh, thus reducing the computational cost.

La compréhension du comportement de structures jusqu'à leur ruine est nécessaire pour concevoir au mieux ces structures. Selon le matériau et les sollicitations considérées, les mécanismes physiques à l'origine de la rupture changent. Nous nous intéresserons à des matériaux homogènes pour lesquels la ruine passe par le développement de fissures autour desquelles les non-linéarités de comportement n'ont pas un rôle dominant. Ces conditions sont réunies pour les matériaux fragiles pour lesquels la source principale de dissipation est la génération non réversible d'une surface libre, et pour certaines fissures de fatigue. Sur un cycle de chargement, il existe de nombreuses applications pour lesquelles les non-linéarités restent confinées. La théorie de la mécanique linéaire élastique de la rupture est alors un modèle pertinent pour approcher le comportement de la structure. Sous ces hypothèses, le front de la fissure introduit une singularité. L'étude asymptotique de cette singularité dans des situations plane et anti-plane permet de définir les séries de Williams. La singularité est alors d'ordre un demi et elle est quantifiée par les facteurs d'intensité des contraintes (FIC) pour chacun des trois modes de sollicitations. En 3D, la fissure peut avoir une géométrie complexe, et aucune expression générale de la singularité n'existe. Dans cette thèse, les séries de Williams en déplacements sont utilisées et régularisées le long du front au sens des éléments finis. À partir de cette définition 3D des séries asymptotiques en pointe de fissure, une méthode d'extraction directe des FIC (DEK-FEM) est étendue au cas 3D. Le domaine est décomposé en deux domaines, raccordés en moyenne sur l'interface. Au voisinage du front, les champs mécaniques sont approchés par une troncature des champs asymptotiques. La singularité est donc traitée avec des champs adaptés, et les degrés de liberté associés sont directement les coefficients asymptotiques. Parmi ces coefficients asymptotiques, on retrouve les FIC et les T-stresses. Pour des raisons d'efficacité numérique et pour pouvoir relier l'échelle de la structure à l'échelle de la fissure, cette méthode est intégrée dans un contexte multigrilles localisées X-FEM. Ainsi nous montrons que cette approche permet une bonne évaluation des évolutions des FIC et du T-stress. Cette méthode est développée en parallèle d'une stratégie de post-traitement expérimental (mesure de champs de déplacements par corrélation d'images) basée sur les mêmes séries asymptotiques. Les images tridimensionnels d'un essai de fatigue in situ sont obtenues par micro-tomographie à rayons X et reconstruction. La corrélation et la régularisation basées sur les séries asymptotiques permettent d'obtenir la géométrie de la fissure et les FIC pour pouvoir identifier des lois de propagation de fissures 3D en fatigue. L'efficacité de cette méthode en parallèle d'une simulation DEK-FEM est illustrée en 2D. / It is necessary to understand the behavior of structures up to their failure to enhance their design. The mechanisms and phenomena undergoing failure vary according to the considered material and boundary conditions. We consider homogeneous materials for which cracks propagate in a context where behavior nonlinearities are not dominants. These conditions are matched for brittle and quasi-brittle materials and for some fatigue cracks. For the former, the main source of dissipation is the crack propagation which can be seen as the generation of a new free-surface. For the later, there is many applications where, in one loading cycle, the nonlinearities remains confined around the crack tip. The linear elastic fracture mechanics theory is then a pertinent model to approximate the structure behavior. Under such hypotheses, a singularity appears in the crack tip vicinity. The Williams' series expansion is computed from the asymptotic study of plane and anti-plane states. The stress is singular at the crack tip and the order of this singularity is one out of two. The singularity amplitude is quantified by the stress intensity factors (SIF), one for each of the three loading modes. In 3D, the crack shape is potentially complex (front curvature and non-planar crack), and no general asymptotic series expansion exists. In this PhD thesis, the 2D Williams' series in displacements are used and regularized with a finite element evolution along the front. From this 3D definition of the asymptotic fields in the crack tip vicinity, a numerical method for direct estimation of the SIF (DEK-FEM) is extended to 3D. This method is based on domain decomposition, the two domains are bounded in a weak sense on their interface. In the crack tip vicinity, the mechanical fields are approximated by a truncation of the asymptotic series expansion. Therefore, appropriate fields are used to deal with the singularity, and the associated degrees of freedom are directly the asymptotic coefficients. Among these coefficients are the SIF and the T-stresses. To bridge the scales between the structure and the crack front singularity and to increase the numerical efficiency, this method is embedded in a localized X-FEM multigrids approach. The proposed method is shown to provide an accurate evaluation of the SIF and T-stresses evolution. This approach has been developed in combination of an experimental post-processing method (full field displacement measurement through image correlation) based on the same asymptotic series expansion. The 3D images can be obtained for in situ fatigue experiments by X-ray microtomography and reconstruction. The crack geometry and the SIF are then provided by image correlation and regularization based on Williams series expansion. These data can be used for identifying a 3D fatigue crack growth law. The efficiency of the method is illustrated in 2D.

Matériaux incontournables des matériaux fonctionnels, les alliages à mémoire de forme(AMF) peuvent présenter de très larges déformations réversibles. La Transformation de Phase Martensitique (TPM), ayant lieu lorsqu’il est soumis à une action mécanique ou thermique, lui confère des caractéristiques particulières. Le comportement thermomécanique des AMF est à présent bien maîtrisé. Cependant la connaissance de leur comportement `a la rupture reste un enjeu majeur pour leur dimensionnement dans le cadre de leur industrialisation pérenne. Ces travaux de recherche se sont attachés `a la connaissance, la description et la quantification du phénomène de localisation en pointe de fissure liée à la TPM induite sous contrainte, au travers du développement d’un modèle prédictif et de sa corrélation expérimentale par mesures de champs simultanées lors d’essais de rupture sur des éprouvettes fissurées de NiTi. Deux modèles analytiques basés sur la mécanique linéaire de la rupture, intégrant le caractère dissymétrique du comportement des AMF en traction/compression, ont été développés pour la prédiction des zones de transformation au voisinage de la pointe de fissure en tenant compte des différents modes de rupture ( élémentaires et mixtes I+II) et du rayon de courbure en pointe de fissure. Un banc de caractérisation par mesures simultanées de champs cinématiques par corrélation d’images (DIC) et thermique par thermographie infrarouge a été développé pour cartographier les champs expérimentaux d’essais de rupture en mode I sur des éprouvettes pré-fissurées. Cette bonne corrélation des modèles analytiques ouvre de nombreuses perspectives concernant l’analyse du couplage thermo mécanique associé à la TPM en pointe de fissure, l’enrichissement des modèles analytiques initiaux, et la confrontation avec les résultats expérimentaux pour des modes de rupture plus complexes (II et mixte I+II). / Major player among functional materials, Shape Memory Alloys (SMA) may undergo verylarge reversible strain. SMA exhibit a Martensitic Phase Transformation (MPT) when they aresubmitted to mechanical or thermal actions, and that gives them some specific characteristics.The thermomechanical behavior of SMA is now well controlled. However, the knowledge of theSMA fracture behavior is a major challenge for their design and sizing for their sustainableindustrialization. This research project has focused on the understanding, describing and quantifyingof the phenomenon of localization at the crack tip due to stress-induced MPT. The study includestwo main aspects: the development of an analytical model and its experimental correlation bysimultaneous field’s measurements during tests on cracked NiTi specimens. Two analytical modelsbased on the linear fracture mechanics and those introduce the asymmetrical nature of the SMAbehavior in tension/compression, were developed for the prediction of transformation zones in thevicinity of the crack tip, taking into account the fracture mode (elementary and mixed ones)and the radii of curvature of the crack tip. A testbench with the measurement of simultaneouskinematic field with Digital Image Correlation (DIC) and thermal field with infrared thermographywas designed for mapping the experimental fields during fracture tests in mode I on pre-crackedspecimen. This good correlation of analytical models opens up many perspectives on the analysisof thermomechanical coupling associated with the MPT at the crack tip, the enrichment of the initialanalytical models, and comparison with experimental results for more complex failure modes (II andmixed I+II).

The presented diploma thesis concerns linear fracture mechanics and deals with determination of the stress intensity factor of the finite crack, which is located in the vicinity of the bimaterial interface, solved by the distributed dislocation technique and theory of complex potencials. The work is possible to devide into three parts. The first part includes basic concepts of the linear fracture mechanics and is also dedicated to the mechanics of composite materials. The second part deals with the determination of the stress intensity factor from solving singular integral equation formulated by Bueckner's principle and the distributed dislocation technique. The third part includes the specific configuration of the crack with respect to the bimaterial interface and the solution, which is compared with results obtained from the FE analysis.

The objective of this diploma thesis is the stress-strain analysis of the crack terminating at the orthotropic bi-material interface suggested as the plane problem of the linear fracture mechanics. The first part is engaged in basic relations of the linear fracture mechanics. The second part is focused on the singularity exponent evaluation for the crack impinging and generally inclined with respect to the bi-material interface. It follows the determination of the generalized stress intensity factors applying the analytical-numerical approach represented by the finite element analysis. The last part of this work is focused on the testing of algorithms applied to the specific crack and bi-material interface configurations. A conclusion discusses the influence of the bi-material mechanical properties and the angel of the crack inclination to the obtained numerical results.

A new methodology of polymer pipe lifetime estimation taking into account residual stresses is described in this thesis. Engineering equations derived based on numerical simulations of a hydrostatic pressure test are proposed. Residual lifetime calculations were performed for different loading conditions using experimental data of a creep crack propagation in studied material and stress distribution in the pipe wall. The effects which significantly influence lifetime estimation were quantified with special focus on residual stresses.

The presented diploma thesis deals with a problem of a generalized stress intensity factor determination and a consecutive study of stress distribution around the bimaterial notch tip, combining analytical and numerical methods. This task is possible to sectionalize into three parts. The first part is dedicated to the fundamentals of the linear fracture mechanics and the mechanics of composite materials. The second part deals with methods of anisotropic plane elasticity solution. Pursuant to the solution the computational models in the third part are created. The first model makes for determination of a singularity exponent eigenvalue by dint of Lekhnitskii-Eshelby-Stroh formalism. The second model makes for determination of the generalized stress intensity factor using psi-integral method, which is based on the Betti reciprocal theorem. All needed calculation are performed in the software ANSYS 12, Maple 12 and Silverforst FTN95. Results will be compared with the values obtained from a direct method of the generalised stress intensity factor determination.