Analysis of fracture behaviour of the cement mantle of reconstructed acetabulum

In this study, the finite element method was used to analyse the crack behaviour in the cement of a reconstructed acetabulum by computing the stress intensity factors at the crack tip. Three loading cases were examined (Fig. 6). These cases present the different human body postures. Both positions and orientations of crack effect on the SIF variation were analysed. When valuating the crack position effect, one notices no risk of crack propagation under the load type 1; however, under the load type 2 and the load type 3 this risk is more important. Load type 3 is the most dangerous loading condition. When computing crack orientation, one noted that the risk of crack propagation is higher when the crack inclination is 20° and 100°.     

Numerical analysis of the notch effect and the behaviour of notch crack in adhesively bonded composite laminates

In this study the finite element method is used to analyse the notch effect and the behaviour of notch cracks in adhesively composite laminate under tension by computing respectively the stress concentration factor at the notch tip which characterize the notch strength and the stress intensity factor at the crack tip which characterize the resistance to the crack propagation. The effects of the adhesive properties and fiber orientation on the variation of both stress concentration and stress intensity factors are highlighted. The obtained results show that the notch strength is reduced in the layer of the laminate of which the fiber orientation is in the applied load direction; the resistance to the crack propagation is also reduced in this type of layer. The stress intensity factor at the tip of notch crack exhibits an asymptotic behaviour as the crack length increases.     

Internal–external circumferential crack behaviour in the cement layer of total hip replacement

This study aimed to investigate crack behaviour at the internal and external surfaces of the cement layer in total hip replacement. A three‐dimensional model of the femur with the cemented prosthesis was developed and analysed. Cracks were placed on the internal, external and both internal and external surfaces of the cement layer. Stress intensity factors were measured during gait. Results revealed that the stress intensity factors modes I and III were the most dominant in the crack propagation in the cement layer. The domain of mode I was the medial and lateral sides of the cement layer. Meanwhile, the domain of mode III was the anterior and posterior sides of the cement layer. The stress intensity factor and distance from the distal end indicated an inverse relationship. The internal and external cracks had no significant interaction. Moreover, stress intensity factors at the external surface of the cement layer were higher than those on the internal surface.     

Numerical analysis of the crack growth path in the cement mantle of the reconstructed acetabulum

In this study, we use the finite element method to analyze the propagation's path of the crack in the orthopedic cement of the total hip replacement. In fact, a small python statement was incorporated with the Abaqus software to do in loop the following operations: extracting the crack propagation direction from the previous study using the maximal circumferential stresses criterion, drawing the new path, meshing and calculating again (stresses and fracture parameters). The loop is broken when the user's desired crack length is reached (number of propagations) or the value of the mode I stress intensity factor is negative. Results show that the crack propagation's path can be influenced by human body posture. The existing of a cavity in the vicinity of the crack can change its propagation path or can absolutely attract it enough to meet it. Crack can propagate in the outward direction (toward the acetabulum bone) and cannot propagate in the opposite direction, the mode I stress intensity factor increases with the crack length and that of mode II vanishes.     

FE analysis of the behaviour of octagonal bonded composite repair in aircraft structures

Bonded composite repair has been recognized as an efficient and economical method to extend the fatigue life of cracked aluminium components. In this work, the finite element method is applied to analyze the central crack’s behaviour repaired by a boron/epoxy composite patch. The knowledge of the stress distribution in the neighbourhood of cracks has an importance for the analysis of their repair according to the patch geometry. The effects of mechanical and geometrical properties of the patch on the variation of the stress intensity factor at the crack tip were highlighted. The obtained results show that the stress intensity factor at the repaired crack with composite patch of height 2c/3 is reduced about 5% compared to cracks repaired with octagonal patch of size c. For patch height of c/3 the reduction is about 7%. The adhesive properties must be optimised in order to increase the repair performances and to avoid the adhesive failure.     

Calculation of stress intensity factors by the force method

The force method is a simple and accurate technique for calculating stress intensity factors (SIFs) from finite element (FE) models, but it has been scarcely used. This paper shows three important advantages of the force method, which make it particularly attractive for designers and researchers. First, it can be employed without special singular quadratic finite elements at the crack tip. Actually, linear reduced integration elements may be used. Second, the force method can be applied to highly anisotropic materials without requiring knowledge of complicated elasticity relations for the stress field around the crack tip. Third, it can handle mixed-mode fracture problems.     

On fracture analysis using an element overlay technique

In this paper, an element overlay technique (s-FEM [Comput. Struct. 43 (1992) 539]) is applied to various two dimensional linear fracture problems. When s-FEM is adopted, local finite element model concerning cracks can be built independently from the global finite element mesh for modeling overall structure. The local model is superposed on the global one. Therefore, it is tractable to introduce cracks in an existing finite element model. The accuracy of s-FEM is critically examined and it is found that the size of local mesh region needs to be larger than or roughly equal to that of an element in the global mesh.     

Multiple penny-shaped cracks interaction in a finite body and their effect on stress intensity factor

In this paper we investigate the stress intensity factors (SIFs) of multiple penny-shaped cracks in an elastic solid cylinder under mode I (axial tension) loading. The cracks are located symmetrically and in parallel to one another in the isotropic cylinder. The fractal-like finite element method (FFEM) is employed to study the interaction of multiple cracks and to demonstrate the efficiency of the FFEM for multiple crack problems. The results show that the SIF values of the inner cracks, which are denoted as crack number 1,2,3,…,(n+1)/2 of a stack of n parallel cracks, are lower than the SIF values of a single crack by between 16% and 48%. Also, the outermost crack, that is the crack closest to the boundaries of a multiple cracked body, has the highest SIF values and is, therefore, likely to fail first.     

Numerical investigation of through crack behaviour under welding residual stresses

Numerous engineering structures operate under the presence of residual stresses resulting from welding or other manufacturing processes. In the present work, the effect of typical residual stress fields on stress intensity factors and crack propagation angle of cracks developing into the residual stress field under mixed mode loading conditions is studied. For the calculations a numerical methodology based on linear elastic finite element analysis is used. The presented results provide a useful tool for an efficient assessment of the influence of residual stress field on the crack evolution behaviour.     

Stress intensity factor error index for finite element analysis with singular elements

An error index for the stress intensity factor (SIF) obtained from the finite element analysis results using singular elements is proposed. The index was developed by considering the facts that the analytical function shape of the crack tip displacement is known and that the SIF can be evaluated from the displacements only. The advantage of the error index is that it has the dimension of the SIF and converges to zero when the actual error of the SIF by displacement correlation technique converges to zero. Numerical examples for some typical crack problems, including a mixed mode crack, whose analytical solutions are known, indicated the validity of the index. The degree of actual SIF error seems to be approximated by the value of the proposed index.     

Numerical methods for determination of stress intensity factors of singular stress field

The asymptotic solution of the singular stress field near a singular point is generally comprised of one or more singular terms in the form of Kr^λ-1f_ij(θ). Based on the asymptotic solution of the singular stress field and the common numerical solution (stresses or displacements) obtained by an ordinary tool such as the finite element method or boundary element method, a simple and effective numerical method is developed to calculate stress intensity factors for one and two singularities. Three examples show that the stress intensity factors evaluated using the method proposed in this paper are very accurate.     

Influence of disbond on notch crack behaviour in single bonded lap joints

The behaviour of the adhesive bonded joints due to the imposed eccentric loading generates a very complex distribution of the stress in the structure. Good adhesion between substrate and adhesive ensures a successful and lasting assembly. In this study the finite element method is used to analyze the behaviour of a bonded lap joint of dissimilar materials. The effects of the mechanical properties of the joints on the shear stress variation with and without presence of a circular notch are investigated. The results show that the maximum shear stresses are located at a distance of about 18% that of the lap length whatever the type of material used. In addition, the stress intensity factor is amplified by the presence of the negative effect of disband whose increase is linearly proportional to the square of the stress intensity factor. It reached its maximum value for a crack length equal to two-fifths of the notch radius.     

Stress intensity factors in spot welds

This paper looks at stress intensity factors of cracks in resistance spot welded joints. Stress intensity factors have been used in the past to predict fatigue crack propagation life of resistance spot welds. However, the stress intensity factors from all previous work was based on assumed initial notch cracks at the nugget, parallel to the sheets. Physical evidence shows, however, that fatigue cracks in spot welds propagate through the thickness of the sheets rather than through the nugget. In this work, stress intensity factors of assumed notch cracks and through thickness cracks in tensile shear (TS) and modified coach peel (MCP) specimens were determined by the finite element method. The finite element results from the assumed notch cracks were compared with the results in the literature and were found to be in agreement with the results from Zhang’s equations [Int. J. Fract. 88 (1997) 167]. The stress intensity factors of assumed notch cracks were found to be different from those of through thickness cracks. To date, no analytic equations for stress intensity factors of through thickness cracks in spot welds have been published. In the current work, simple equations are proposed to estimate the K_I and K_II values of through thickness cracks in TS and MCP specimens.     

Computation of the stress intensity factors for repaired cracks with bonded composite patch in mode I and mixed mode

In this study, the finite element method is used to analyse the behaviour of repaired cracks with bonded composite patches in mode I and mixed mode by computing the stress intensity factors at the crack tip. The effects of the patch size and the adhesive properties on the stress intensity factors variation were highlighted. The plot of the stress intensity factors according to the crack length in mode I, shows that the stress intensity factor exhibits an asymptotic behaviour as the crack length increases. In mixed mode, the obtained results show that the Mode I stress intensity factor is more affected by the presence of the patch than that of mode II.     

Analysis of the distribution of thermal residual stresses in bonded composite repair of metallic aircraft structures

In this study, the distribution of the thermal residual stresses due to the adhesive curing in bonded composite repair is analysed using the finite element method. The computation of these stresses comprises all components of the structures: cracked plate, composite patch and adhesive layer. In addition, the influence of these residual stresses on the repair performance is highlighted by analysing their effect on the stress intensity factor at the crack tip. The obtained results show that the normal thermal stresses in the plate and the patch are important and the shear stresses are less significant. The level of the adhesive thermal stresses is relatively high. The presence of the thermal stresses increases the stress intensity factor at the crack tip, what reduce the repair performance.     

Numerical analysis of brittle materials fractured by sharp indenters

Indentation can be used to determine the fracture properties of materials. A detailed investigation is here presented on the reliability of finite element simulations of sharp indentations on cracked specimens. Elastic analyses of the stress and deformation fields arising from Vickers, Berkovich and cube-corner indenters generated the stress intensity factor along the edge of penny-shaped or elliptical cracks. Various materials with a range of properties were analysed and the results compared with published experimental data. Additional measurements from tests on soda-lime glass provided further opportunities for experimental validation of numerical predictions. The variation of the stress intensity factor indicated trends for crack growth patterns, which were consistent with experimental observations. Particular attention was given to the evaluation of the geometry-dependent parameter appearing in the relation yielding the fracture toughness of cracked indented materials. The numerical predictions of this parameter were remarkably consistent with experimental data and results from other approximate methods.     

Radial locations of strain gages for accurate measurement of mode I stress intensity factor

Strain gage methods are popular in experimental determination of stress intensity factors (SIFs). Radial location of gages with respect to the crack tip plays an important role in accuracy of strain measurements and thus accurate determination of SIFs. The present work proposes a finite element based simple, accurate and consistent method for determination of the limiting value of the radial distance (r_max) of a strain gage. This parameter is in turn useful in deciding the valid strain gage location for accurate measurement of opening mode SIF. The results obtained from the present investigation agree well with the theoretical predictions and could be used for experimental determination of SIFs for both single ended and double ended cracked specimens. The r_max values of center cracked and edge cracked plates with different crack length to width ratio are estimated. The results of the present investigation show that the relative size of the crack length and net ligament length strongly influences the r_max value and the effect of Poisson’s ratio is marginal on the r_max value.     

Numerical analysis of interaction and coalescence of numerous microcracks

The deformation and failure behaviors of brittle or quasi-brittle solids are closely related to interaction and propagation of stochastically distributed microcracks. The influence of microcrack interaction and evolution on the mechanical properties of materials presents a problem of considerable interest, which has been extensively argued but has not been resolved as yet. In the present paper, a novel numerical method is used to calculate the effective elastic moduli and the tensile strength, and to simulate the failure process of brittle specimens containing numerous microcracks. The influences of some crack distribution parameters reflecting the non-uniform spatial concentration, size and orientation distributions are examined. The effective elastic moduli and the tensile strength of brittle materials exhibit different dependences on microcrack interaction. For example, two microcrack distributions that lead to the identical effective elastic moduli may cause a pronounced difference in the tensile strengths and failure behaviors of materials. By introducing two criteria for microcrack growth and coalescence in terms of Griffith’s energy release rate, the above numerical method is extended to simulate the coalescence process of microcracks that results in a fatal crack and the final rupture of a specimen.     

FE analysis of stresses and stress intensity factors of interfacial cracks in a CTS specimen

A plane stress finite element analysis was implemented to understand the stress fields for a crack lying at an aluminium/epoxy interface of a compact tension and shear specimen. The interaction integral method was used to separate the mixed-mode stress intensity factors at the interfacial crack-tip under different loading modes, which can have important implications for characterisation of interfacial crack growth.     

Accuracy and limit of a least-squares method to calculate 3D notch SIFs

To evaluate the three-dimensional (3D) stress intensity factors (SIFs) of a sharp V-notch using the finite element result is limited in the literature. Thus, this study developed a least-squares method to solve this problem as well as study its restriction and accuracy. First, the William’s eigenfunction and complex stress function approach are deduced into a least-squares form, and then stress field from the finite element analysis is substituted into the least-squares equation to evaluate the 3D SIFs. Numerical simulations in this article show that the least-squares method can be used to calculate SIFs accurately if more than two stress terms are included. The calculated SIFs of this least-squares method are not sensitive to the maximum and minimum radiuses of the area from which data are included. The major advantage of the proposed method is that the procedure is simple and systematic, so it can be applied to any finite element code without difficulties.