A general methodology for calculating mixed mode stress intensity factors and fracture toughness of solder joints with interfacial cracks

Solder joints in electronic packages undergo thermo-mechanical cycling, resulting in nucleation of micro-cracks, especially at the solder/bond-pad interface, which may lead to fracture of the joints. The fracture toughness of a solder joint depends on material properties, process conditions and service history, as well as strain rate and mode-mixity. This paper reports on a methodology for determining the mixed-mode fracture toughness of solder joints with an interfacial starter-crack, using a modified compact mixed mode (CMM) specimen containing an adhesive joint. Expressions for stress intensity factor (K) and strain energy release rate (G) are developed, using a combination of experiments and finite element (FE) analysis. In this methodology, crack length dependent geometry factors to convert for the modified CMM sample are first obtained via the crack-tip opening displacement (CTOD)-based linear extrapolation method to calculate the under far-field mode I and II conditions (f(1a) and f(2a)), (ii) generation of a master-plot to determine a(c), and (iii) computation of K and G to analyze the fracture behavior of joints. The developed methodology was verified using J-integral calculations, and was also used to calculate experimental fracture toughness values of a few lead-free solder-Cu joints. (C) 2014 Elsevier Ltd. All rights reserved.

A meshless adaptive multiscale method for fracture

The paper presents a multiscale method for crack propagation. The coarse region is modelled by the differential reproducing kernel particle method. Fracture in the coarse scale region is modelled with the Phantom node method. A molecular statics approach is employed in the fine scale where crack propagation is modelled naturally by breaking of bonds. The triangular lattice corresponds to the lattice structure of the (111) plane of an FCC crystal in the fine scale region. The Lennard-Jones potential is used to model the atom-atom interactions. The coupling between the coarse scale and fine scale is realized through ghost atoms. The ghost atom positions are interpolated from the coarse scale solution and enforced as boundary conditions on the fine scale. The fine scale region is adaptively refined and coarsened as the crack propagates. The centro symmetry parameter is used to detect the crack tip location. The method is implemented in two dimensions. The results are compared to pure atomistic simulations and show excellent agreement. (C) 2014 Elsevier B. V. All rights reserved.

Fracture in metallic glasses: mechanics and mechanisms

Significant progress in understanding the mechanical behavior of metallic glasses (MGs) was made over the past decade, particularly on mechanisms of plastic deformation. However, recent research thrust has been on exploring the mechanics and physics of fracture. MGs can be very brittle with K-Ic values similar to silicate glasses and ceramics or very tough with K-Ic akin to high toughness crystalline metals. Even the tough MGs can become brittle with structural relaxation following annealing at temperatures close to glass transition temperature (T-g). Detailed experimental studies coupled with complementary numerical simulations of the recent past have provided insights on the micromechanisms of failure as well as nature of crack tip fields, and established the governing fracture criteria for ductile and brittle glasses. In this paper, the above advances are reviewed and outstanding issues in the context of fracture of amorphous alloys that need to be resolved are identified.

On the variability in fracture toughness of `ductile' bulk metallic glasses

The mode I fracture toughness, K-Ic, of ductile bulk metallic glasses (BMGs) exhibits a high degree of specimen-to-specimen variability. By conducting fracture experiments in modes I and II, we demonstrate that the observed high variability in mode I, vis-a-vis mode II, is a result of highly variable propensity for the conversion of shear bands into cracks in mode I whereas in mode II, crack growth direction is fixed. Thus, the measured variability in K-Ic is intrinsic to the nature of BMGs. (C) 2015 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Optimization of clamped beam geometry for fracture toughness testing of micron-scale samples

Fracture toughness measurements at the small scale have gained prominence over the years due to the continuing miniaturization of structural systems. Measurements carried out on bulk materials cannot be extrapolated to smaller length scales either due to the complexity of the microstructure or due to the size and geometric effect. Many new geometries have been proposed for fracture property measurements at small-length scales depending on the material behaviour and the type of device used in service. In situ testing provides the necessary environment to observe fracture at these length scales so as to determine the actual failure mechanism in these systems. In this paper, several improvements are incorporated to a previously proposed geometry of bending a doubly clamped beam for fracture toughness measurements. Both monotonic and cyclic loading conditions have been imposed on the beam to study R-curve and fatigue effects. In addition to the advantages that in situ SEM-based testing offers in such tests, FEM has been used as a simulation tool to replace cumbersome and expensive experiments to optimize the geometry. A description of all the improvements made to this specific geometry of clamped beam bending to make a variety of fracture property measurements is given in this paper.

Fracture behavior of magnesium alloys - Role of tensile twinning

In this work, Mode-I fracture experiments are conducted using notched compact tension specimens machined from a rolled AZ31 Mg alloy plate having near-basal texture with load applied along rolling direction (RD) and transverse direction (TD). Moderately high notched fracture toughness of J(C) similar to 46 N/mm is obtained in both RD and TD specimens. Fracture surface shows crack tunneling at specimen mid-thickness and extensive shear lips near the free surface. Dimples are observed from SEM fractographs suggesting ductile fracture. EBSD analysis shows profuse tensile twinning in the ligament ahead of the notch. It is shown that tensile twinning plays a dual role in enhancing the toughness in the notched fracture specimens with reduced triaxiality. It provides significant dissipation in the background plastic zone and imparts hardening to the material surrounding the fracture process zone via operation of several mechanisms which retards micro-void growth and coalescence. (C) 2015 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

A thermodynamic correlation between damage and fracture as applied to concrete fatigue

Fatigue damage in concrete is characterized by the simultaneous presence of micro and macrocracks. The theory of fracture mechanics conveniently handles the propagation of macrocracks, whereas damage mechanics precisely describes the state of microcracking. This paper provides a platform to correlate fracture mechanics and damage mechanics theories through an energy equivalence within a thermodynamic framework by equating the energy dissipated according to each theory. Through this correlation, damage corresponding to a given crack length could be obtained, and alternatively a discrete crack could be transformed into an equivalent damage zone. The results are validated using available experimental data on concrete fatigue including stiffness degradation and acoustic emission. (C) 2015 Elsevier Ltd. All rights reserved.

A thermodynamic correlation between damage and fracture as applied to concrete fatigue

Fatigue damage in concrete is characterized by the simultaneous presence of micro and macrocracics. The theory of fracture mechanics conveniently handles the propagation of macrocracks, whereas damage mechanics precisely describes the state of microcracking. This paper provides a platform to correlate fracture mechanics and damage mechanics theories through an energy equivalence within a thermodynamic framework by equating the energy dissipated according to each theory. Through this correlation, damage corresponding to a given crack length could be obtained, and alternatively a discrete crack could be transformed into an equivalent damage zone. The results are validated using available experimental data on concrete fatigue including stiffness degradation and acoustic emission. (C) 2015 Elsevier Ltd. All rights reserved.

Twinning induced enhancement of fracture toughness in ultrafine grained Hydroxyapatite-Calcium Titanate composites

Despite being highly bioactive and biocompatible, the limitations of monolithic hydroxyapatite (HA) include extremely low fracture toughness, poor electrical conductivity. While addressing these issues, the present study demonstrates how CaTiO3 (CT) addition to HA can be utilized to obtain a combination of long crack fracture toughness (1.7 MPa m(1/2) SEVNB technique) and flexural strength of 98-155 MPa (3-point bending) and a moderate tensile strength (diametral compression) of 17-36 MPa. The enhancement in fracture resistance in spark plasma sintered HA-CT composites has been explained in reference to the observed twin morphology. TEM reveals the presence of twins in CT grains due to 1800 rotation about 101]. The measured properties along with our earlier reports on biocompatibility and electrical properties make HA-CT suitable for bone tissue engineering applications. When compared with other competing HA-based biocomposites, HA-CT composites are found to have a better combination of properties useful for medium load bearing implant applications. (C) 2015 Elsevier Ltd. All rights reserved.

Fracture Testing at Small-Length Scales: From Plasticity in Si to Brittleness in Pt

The field of micro-/nano-mechanics of materials has been driven, on the one hand by the development of ever smaller structures in devices, and, on the other, by the need to map property variations in large systems that are microstructurally graded. Observations of `smaller is stronger' have also brought in questions of accompanying fracture property changes in the materials. In the wake of scattered articles on micro-scale fracture testing of various material classes, this review attempts to provide a holistic picture of the current state of the art. In the process, various reliable micro-scale geometries are shown, challenges with respect to instrumentation to probe ever smaller length scales are discussed and examples from recent literature are put together to exhibit the expanse of unusual fracture response of materials, from ductility in Si to brittleness in Pt. Outstanding issues related to fracture mechanics of small structures are critically examined for plausible solutions.

In-situ study of microscale fracture of diffusion aluminide bond coats: Effect of platinum

The influence of Pt layer thickness on the fracture behavior of PtNiAl bond coats was studied in situ using clamped micro-beam bend tests inside a scanning electron microscope (SEM). Clamped beam bending is a fairly well established micro-scale fracture test geometry that has been previously used in determination of fracture toughness of Si and PtNiAl bond coats. The increasing amount of Pt in the bond coat matrix was accompanied by several other microstructural changes such as an increase in the volume fraction of alpha-Cr precipitate particles in the coating as well as a marginal decrease in the grain size of the matrix. In addition, Pt alters the defect chemistry of the B2-NiAl structure, directly affecting its elastic properties. A strong correlation was found between the fracture toughness and the initial Pt layer thickness associated with the bond coat. As the Pt layer thickness was increased from 0 to 5 mu m, resulting in increasing Pt concentration from 0 to 14.2 at.% in the B2-NiAl matrix and changing alpha-Cr precipitate fraction, the initiation fracture toughness (K-IC) was seen to rise from 6.4 to 8.5 MPa.m(1/2). R-curve behavior was observed in these coatings, with K-IC doubling for a crack propagation length of 2.5 mu m. The reasons for the toughening are analyzed to be a combination of material's microstructure (crack kinking and bridging due to the precipitates) as well as size effects, as the crack approaches closer to the free surface in a micro-scale sample.

Fracture mechanics of piezoelectric materials

This paper presents an analysis of crack problems in homogeneous piezoelectrics or on the interfaces between two dissimilar piezoelectric materials based on the continuity of normal electric displacement and electric potential across the crack faces. The explicit analytic solutions are obtained for a single crack in an infinite piezoelectric or on the interface of piezoelectric bimaterials. For homogeneous materials it is found that the normal electric displacement D-2, induced by the crack, is constant along the crack faces which depends only on the remote applied stress fields. Within the crack slit, the perturbed electric fields induced by the crack are also constant and not affected by the applied electric displacement fields. For bimaterials, generally speaking, an interface crack exhibits oscillatory behavior and the normal electric displacement D-2 is a complex function along the crack faces. However, for bimaterials, having certain symmetry, in which an interface crack displays no oscillatory behavior, it is observed that the normal electric displacement D-2 is also constant along the crack faces and the electric field E-2 has the singularity ahead of the crack tip and has a jump across the interface. Energy release rates are established for homogeneous materials and bimaterials having certain symmetry. Both the crack front parallel to the poling axis and perpendicular to the poling axis are discussed. It is revealed that the energy release rates are always positive for stable materials and the applied electric displacements have no contribution to the energy release rates.

Dislocation theory of the fracture criterion for anisotropic solids

A dislocation theory of fracture criterion for the mixed dislocation emission and cleavage process in an anisotropic solid is developed in this paper. The complicated cases involving mixed-mode loading are considered here. The explicit formula for dislocations interaction with a semi-infinite crack is obtained. The governing equation for the critical condition of crack cleavage in an anisotropic solid after a number dislocation emissions is established. The effects of elastic anisotropy, crack geometry and load phase angle on the critical energy release rate and the total number of the emitted dislocations at the onset of cleavage are analysed in detail. The analyses revealed that the critical energy release rates can increase to one or two magnitudes larger than the surface energy because of the dislocation emission. It is also found elastic anisotropy and crystal orientation have significant effects on the critical energy release rates. The anisotropic values can be several times the isotropic value in one crack orientation. The values may be as much as 40% less than the isotropic value in another crack orientation and another anisotropy parameter. Then the theory is applied to a fee single crystal. An edge dislocation can emit from the crack tip along the most highly shear stressed slip plane. Crack cleavage can occur along the most highly stressed slip plane after a number of dislocation emissions. Calculation is carried out step by step. Each step we should judge by which slip system is the most highly shear stressed slip system and which slip system has the largest energy release rate. The calculation clearly shows that the crack orientation and the load phase angle have significant effects on the crystal brittle-ductile behaviours.

Micro-Fracture Toughness in a Nanoceramic Composite

Micro-fracture toughness in one nanoceramic composite, Al_2O_3/SiC_{nano}, is measured by Nanoindentation tests. Microscopy is adopted to probe that all nanoindentations are local at Al_2O_3, matrix component of the current composite and measure length of radial micro-cracks at corriers of nanoindentation. The measured fracture toughness at the shallow indentation is higher than those at deeper ones. Relations between micro-toughness and peak loads as well as peak depths are discussed. And an analysis to the energy consumed during the indentation and fracture toughness explained that, the lower fracture toughness measured in the deeper impression probably resulted from its fracture modes yielded by higher load and larger depth.

Numerical Analysis of Rainfall Infiltration in the Slope with a Fracture

With the finite volume method, a 2D numerical model for seepage in unsaturated soil has been established to study the rainfall infiltration in the fractured slope.The result shows that more rain may infiltrate into the slope due to existing fracture and then the pore pressure rises correspondingly. Very probably, it is one of the crucial factors accounting for slope failure. Furthermore a preliminary study has been conducted to investigate the influence of various fracture and rainfall factors such as the depth, width and location of a crack, surface condition, rainfall intensity and duration. Pore pressure and water volumetric content during the transient seepage are carefully examined to reveal the intrinsic mechanism.