Superalloys have performed well as blade and disc materials in turbine engines due to their exceptional elevated temperature strength, high resistance to creep, oxidation and corrosion as well as good fracture toughness. This study explores the use of a relatively new technique of strain measurement, high resolution electron backscatter diffraction (HR-EBSD) to measure local deformation fields. The heart of the HR-EBSD technique lies in comparing regions in EBSD patterns from a strained region of a sample to those in a pattern from an unstrained region. This method was applied to study the elastic strain fields and geometrically necessary dislocation density (GND density) distribution near hard carbide particles in a nickel-based superalloy MAR-M-002. Significant thermal strains were initially induced by thermal treatment, which included a final cooling from the ageing temperature of 870°C. Elastic strains were consistent with a compressive radial strain and tensile hoop strain that was expected as the matrix contracts around the carbide. The mismatch in thermal expansion coefficient of the carbide particles compared to that of the matrix was sufficient to have induced localized plastic deformation in the matrix leading to a GND density of 3 x 10¹³ m^–2 in regions around the carbide. These measured elastic strain and GND densities have been used to help develop a crystal plasticity finite element model in another research group and some comparisons under thermal loading have also been examined. Three-point bending was then used to impose strain levels within the range ±12% across the height of a bend bar sample. GND measurements were then made at both carbide-containing and carbide-free regions at different heights across the bar. The average GND density increases with the magnitude of the imposed strain (both in tension and compression), and is markedly higher near the carbide particles. The higher GND densities near the carbides (order of 10¹⁴ per m²) are generated by the large strain gradients produced around the plastically rigid inclusion during monotonic mechanical deformation with some minor contribution from the pre-existing residual deformation from thermal loading. A method was developed of combining the local EBSD measurements with FE modelling to set the average residual strains within the mapped region even when a good strain free reference point was unavailable. Cyclic loading was then performed under four point loading to impose strain levels of about ±8% across the height of bend bar samples. Similar measurements as in the case of monotonic deformation were made at several interruptions to fatigue loading. Observations from the cyclic loading such as slip features, carbide cracking, GND density accumulation have been explored around carbide particles, at regions away from them and near a grain boundary.

The factors controlling the reliability of Pb-free solders when subject to thermomechanical regimes relevant to the harsh aerospace environment have been studied. Ball grid array (BGAs) typical of microelectronic devices have been manufactured in-house and subjected to isothermal ageing and thermal cycling. The BGAs comprised both Cu and Ni-Au metallizations, Pb-free Sn-Ag-Cu 400 and 600μm solder balls, FR4 and Al₂O₃ boards, and included circuits to measure resistance changes due to damage in the joints during thermal cycling. Microstructural evolution within the solders balls and complex interfacial reactions were studied in all configurations using various types of electron microscopy. The mechanical properties of the different phases formed within solder joints were studied using nanoindentation at room and elevated temperatures up to 175°C for the first time. Intermetallic compounds (IMCs) were stiff, hard and brittle with very low creep rates, while the softer primary Sn, eutectic regions and Cu metallization readily underwent creep. Two-dimensional finite element analysis (FEA) of nanoindentation was used to understand better the physical meaning of nanoindentation creep data. Reliability experiments comprised both thermal cycling and FEA of BGAs. The difference in coefficient of thermal expansion (CTE) in the BGA materials caused interfacial fatigue damage in the solder joints, which was detected primarily at the solder/metallization interface of the outermost, most strained solder joint. Accumulated creep strain per cycle at this interface was evaluated using 3D FEA of the stress-strain state of the BGA and results calibrated against experimental BGA mean lifetimes using the Coffin-Mason relationship. Nanoindentation combined with FEA has been shown to be a viable route for the rapid assessment of creep performance and lifetime in lead-free solders under aerospace thermal cycles.

This thesis has two parts. In the first part, we use the results of new fatigue experiments conducted with variable load levels as well as variable stress ratios to critically assess three (two old and one relatively new) cumulative fatigue damage models. These models are deterministic. Such models are usually tested using multiple blocks of periodic loading with differing amplitudes. However, available data pertains to zero-mean loading, and does not investigate the role of variable stress ratio (Smin/Smax). Here, we present experimental results for variable stress ratios. Two specimen geometries and two materials (Al 2014and Al 2024)are tested. Manson’s double linear damage rule (DLDR)gives the highest accuracy in predicting the experimental outcome, even in the presence of variable stress ratios, whereas predictions of the newer model (“A constructive empirical theory for metal fatigue under block cyclic loading,” Proceedings of the Royal Society A, 464 (2008), 1161-1179) are slightly inferior. The widely used Miner’s rule is least accurate in terms of prediction. The merits and drawbacks of these models, in light of the experimental results, are as follows. The DLDR, though accurate, has minor scientific inconsistencies and no clear generalization. The constructive model has possible generalizability and more appealing scientific consistency, but presently has poorer accuracy. Miner’s rule, least accurate, lies within the constructive approach for special parameter values. The DLDR can guide the new (constructive)approach through new parameter fitting criteria. In the second part of this thesis, we consider the scatter in fatigue life and use the Weibull distribution to describe ‘S-N-P’ curves. We first assume homoscedasticity (load-independent or constant variance) and present a way to draw a p-percentile line on a log-log load-life plot. Then heteroscedasticity (load-dependent variance) in fatigue life is incorporated and a simple statistical model is proposed, to obtain a straight line percentile plot at a pre-specified probability of survival ps. The proposed method is illustrated for Al 2014-T6 and Al 2024-T4 data sets (extracted manually) from MMPDS-01 (a data handbook).

The ability to predict fatigue crack growth rates is essential in safety critical systems. The discovery of fatigue crack closure in 1970 caused a flourish of research in attempts to simulate this behaviour, which crucially affects crack growth rates. Historically, crack tip plasticity models have been based on one-dimensional rays of plasticity emanating from the crack tip, either co-linear with the crack (for the case of plane stress), or at a chosen angle in the plane of analysis (for plane strain). In this thesis, one such model for plane stress, developed to predict fatigue crack closure, has been refined. It is applied to a study of the relationship between the apparent stress intensity range (easily calculated using linear elastic fracture mechanics), and the true stress intensity range, which includes the effects of plasticity induced fatigue crack closure. Results are presented for all load cases for a finite crack in an infinite plane, and a method is demonstrated which allows the calculation of the true stress intensity range for a growing crack, based only on the apparent stress intensity range for a static crack. Although the yield criterion is satisfied along the plastic ray, these one-dimensional plasticity models violate the yield criterion in the area immediately surrounding the plasticity ray. An area plasticity model is therefore required in order to model the plasticity more accurately. This thesis develops such a model by distributing dislocations over an area. Use of the model reveals that current methods for incremental plasticity algorithms using distributed dislocations produce an over-constrained system, due to misleading assumptions concerning the normality condition. A method is presented which allows the system an extra degree of freedom; this requires the introduction of a parameter, derived using the Prandtl-Reuss flow rule, which relates the magnitude of slip on complementary shear planes. The method is applied to two problems, confirming its validity.