Application of the finite element method to problems in fracture mechanics

Numerical techniques have been finding increasing use in all aspects of fracture mechanics, and often provide the only means for analyzing fracture problems. The work presented here, is concerned with the application of the finite element method to cracked structures. The present work was directed towards the establishment of a comprehensive two-dimensional finite element, linear elastic, fracture analysis package. Significant progress has been made to this end, and features which can now be studied include multi-crack tip mixed-mode problems, involving partial crack closure. The crack tip core element was refined and special local crack tip elements were employed to reduce the element density in the neighbourhood of the core region. The work builds upon experience gained by previous research workers and, as part of the general development, the program was modified to incorporate the eight-node isoparametric quadrilateral element. Also. a more flexible solving routine was developed, and provided a very compact method of solving large sets of simultaneous equations, stored in a segmented form. To complement the finite element analysis programs, an automatic mesh generation program has been developed, which enables complex problems. involving fine element detail, to be investigated with a minimum of input data. The scheme has proven to be versati Ie and reasonably easy to implement. Numerous examples are given to demonstrate the accuracy and flexibility of the finite element technique.

The application of finite element analysis to fracture mechanics.

DUE TO COPYRIGHT RESTRICTIONS ONLY AVAILABLE FOR CONSULTATION AT ASTON UNIVERSITY LIBRARY AND INFORMATION SERVICES WITH PRIOR ARRANGEMENT

The relationship between fracture toughness and microstructure in titanium alloys

Linear Elastic Fracture Mechanics has been used to study the microstructural factors controlling the strength and toughness of two alpha-beta, titanium alloys. Fracture toughness was found to be independent of orientation for alloy Ti/6A1/4-V, but orientation dependent for IMI 700, bend and tension specimens giving similar toughness values. Increasing the solution temperature led to the usual inverse relationship between strength and toughness, with toughness becoming a minimum as the beta transus was approached. The production of a double heat treated microstructure led to a 100% increase in toughness in the high strength alloy and a 20% increase in alloy Ti/6A1/4V, with little decrease in strength. The double heat treated microstruoture was produced by cooling from the beta field into the alpha beta field, followed. by conventional solution treatment and ageing. Forging above the beta transus led to an increase in toughness over alpha beta forging in the high strength alloy, but had little effect on the toughness of Ti/6A1/4V. Light and electron microscopy showed that the increased toughness resulted from the alpha phase being changed from mainly continuous to a discontinuous platelet form in a transformed beta matrix. Void formation occurred at the alpha-beta interface and crack propagation was via the interface or across the platelet depending on which process required the least energy. Varying the solution treatment temperature produced a varying interplatelet spacing and platelet thickness. The finest interplatelet spacing was associated with the highest toughness, since a higher applied stress was required to give the necessary stress concentration to initiate void formation. The thickest alpha platelet size gave the highest toughness which could be interpreted in terms of Krafftt's "process zone size" and the critical crack tip displacement criterion by Hahn and Rosenfield from an analysis by Goodier and Field.

The fracture toughness of a range of cast steels

A range of plain carbon, carbon-manganese and low alloy cast steels were tested in order to determine their various fracture toughness values under elastic and elastic-plastic conditions. The main fracture toughness parameters which are considered are (1) Linear Elastic Fracture Mechanics (LEFM), (2) the J-Contour Integral, and (3) Crack Opening Displacement (COD). Results are obtained from fracture toughness specimens of various dimensions and the relevance of the validity criteria to cast steels is considered in some detail. In addition, the effect of casting position on specimen toughness values was noted. Valid KIC results according to LEFM, were obtained for three of the eight cast steels tested. Although KIC values from LEFM were not obtained from the remaining five steels, critical COD and J-integral values were determined. It is postulated that these values and particularly the critical J values can be used, with confidence for material selection or in defect tolerance calculations using these steels. Toughness values were found to vary with casting position in several of the steels tested and the possible reasons for such variations are discussed in the Thesis.

Fracture toughness of a series of high chromium abrasion resistant alloys

The fracture properties of a series of alloys containing 15% chromium and 0.8 to 3.4% carbon are investigated using strain fracture toughness testing techniques. The object of the work is to apply a quantitative method of measuring toughness to abrasion resistant materials, which have previously been assessed on an empirical basis; and to examine the relationship between microstructure and K10 in an attempt to improve the toughness of inherently brittle materials. A review of the relevant literature includes discussion of the background to the alloy series under investigation, a survey of the development of fracture mechanics and the emergence of K10 as a toughness parameter. Metallurgical variables such as composition, heat treatment, grain size, and hot working are ???? to relate microstructure to toughness, and fractographic evidence is used to substantiate the findings. The results are applied to a model correlating ductile fracture with plastic strain instability, and the nucleation of voids. Strain induced martensite formation in austenitic structures is analysed in terms of the plastic energy dissipation mechanisms operating at the crack tip. Emphasis is placed on the lower carbon alloys in the series, and a composition put forward to optimise wear resistance and toughness. The properties of established competitive materials are compared to the proposed alloy on a toughness and cost basis.

Energy-based damage and fracture framework for viscoelastic asphalt concrete

A framework based on the continuum damage mechanics and thermodynamics of irreversible processes using internal state variables is used to characterize the distributed damage in viscoelastic asphalt materials in the form of micro-crack initiation and accumulation. At low temperatures and high deformation rates, micro-cracking is considered as the source of nonlinearity and thus the cause of deviation from linear viscoelastic response. Using a non-associated damage evolution law, the proposed model shows the ability to describe the temperature-dependent processes of micro-crack initiation, evolution and macro-crack formation with good comparison to the material response in the Superpave indirect tensile (IDT) strength test.

Effects of truck traffic on top-down fatigue cracking performance of flexible pavements using a new mechanics-based analysis framework

The mechanics-based analysis framework predicts top-down fatigue cracking initiation time in asphalt concrete pavements by utilising fracture mechanics and mixture morphology-based property. To reduce the level of complexity involved, traffic data were characterised and incorporated into the framework using the equivalent single axle load (ESAL) approach. There is a concern that this kind of simplistic traffic characterisation might result in erroneous performance predictions and pavement structural designs. This paper integrates axle load spectra and other traffic characterisation parameters into the mechanics-based analysis framework and studies the impact these traffic characterisation parameters have on predicted fatigue cracking performance. The traffic characterisation inputs studied are traffic growth rate, axle load spectra, lateral wheel wander and volume adjustment factors. For this purpose, a traffic integration approach which incorporates Monte Carlo simulation and representative traffic characterisation inputs was developed. The significance of these traffic characterisation parameters was established by evaluating a number of field pavement sections. It is evident from the results that all the traffic characterisation parameters except truck wheel wander have been observed to have significant influence on predicted top-down fatigue cracking performance.

Creep and fatigue of a single crystal nickel base superalloy

The introduction of single crystal casting techniques has led to the development of existing nickel-base superalloys to produce materials with optimum mechanical properties in the single crystal condition. As single crystals are known to be anisotropic, a study is needed to determine the general mechanical properties of these materials, and determine the effects of crystal orientation upon them. A study has been carried out to identify the effect of orientation and temperature on the creep and fatigue properties of a development single crystal superalloy, SRR 99. Creep testing and crystal rotation experiments have been made on SRR 99 and an earlier development alloy, SRR 9. Fatigue experiments at elevated temperatures have been carried out on both notched and un-notched specimens of alloy SRR 99. To aid in this analysis, several analytical techniques have been employed including Laue x-ray orientation analysis, measurement of strain by photographic methods and microstructural examination. Crystal rotation experiments have indicated that shear of 1 precipitates by lbrace111rbrace< 112> slip systems is operative during primary creep deformation at temperatures of 750oC and 850oC. The effect of orientation variation obtained by standard casting practices was not found to be significant. Creep rupture was found to be associated with multiple crack initiation from micropores. Fatigue crack initiation in un-notched specimens was found to be related to microporosity and microstructural defects. Failure was predominantly by crystallographic crack growth on lbrace111rbrace planes. The use of linear elastic fracture mechanics to describe fatigue crack propagation in alloy SRR 99 was found to be acceptable at temperatures up to 850oC. Variation of temperature, frequency and crystal orientation was found to have only moderate effect upon crack propagation rates.

Micromechanics of agglomerate damage processes

This thesis reports a detailed investigation of the micromechanics of agglomerate behaviour under free-fall impact, double (punch) impact and diametrical compression tests using the simulation software TRUBAL. The software is based on the discrete element method (DEM) which incorporates the Newtonian equations of motion and contact mechanics theory to model the interparticle interactions. Four agglomerates have been used: three dense (differing in interface energy and contact density) and one loose. Although the simulated agglomerates are relatively coarse-grained, the results obtained are in good agreement with laboratory test results reported in the literature. The computer simulation results show that, in all three types of test, the loose agglomerate cannot fracture as it is unable to store sufficient elastic energy. Instead, it becomes flattened for low loading-rates and shattered or crushed at higher loading-rates. In impact tests, the dense agglomerates experience only local damage at low impact velocities. Semi-brittle fracture and fragmentation are produced over a range of higher impact velocities and at very high impact velocities shattering occurs. The dense agglomerates fracture in two or three large fragments in the diametrical compression tests. Local damage at the agglomerate-platen interface always occurs prior to fracture and consists of local bond breakage (microcrack formation) and local dislocations (compaction). The fracture process is dynamic and much more complex than that suggested by continuum fracture mechanics theory. Cracks are always initiated from the contact zones and propagate towards the agglomerate centre. Fracture occurs a short time after the start of unloading when a fracture crack "selection" process takes place. The detailed investigation of the agglomerate damage processes includes an examination of the evolution of the fracture surface. Detailed comparisons of the behaviour of the same agglomerate in all three types of test are presented. The particle size distribution curves of the debris are also examined, for both free-fall and double impact tests.

Fatigue crack initiation and propagation from defects in aluminium welds

Fatigue crack initiation and propagation in aluminium butt welds has been investigated. It is shown that the initiation of cracks from both buried defects and. from the weld reinforcement may be quantified by predictive laws based on either linear elastic fracture mechanics, or on Neuber's rule of stress and strain ooncentrations. The former is preferable on the grounds of theoretical models of crack tip plasticity, although either may be used as the basis of an effeotive design criteria against crack initiation. Fatigue lives fol1owing initiation were found to follow predictions based on the integration of a Paris type power law. The effect of residual stresses from the welding operation on both initiation and propagation was accounted for by a Forman type equation. This incorporated the notional stress ratio produced by the residual stresses after various heat treatments. A fracture mechanics analysis was found to be useful in describing the fatigue behaviour of the weldments at increased temperatures up to 300°C. It is pointed out, however, that the complex interaction of residual stresses, frequency, and changes in fracture mode necessitate great caution in the application of any general design criteria against crack initiation and growth at elevated. temperatures.

Fatigue cracking from stress concentration in mild steel

Pulsating; tension fatigue tests have been carried out on edge notched specimens of a mild steel. An electrical potential drop technique was used to determine the number of cycles taken to initiate cracks and the rate at which the cracks grew across the specimen. The results could be described by the range of stress intensity factor, which for crack initiation was modified to take account of the notch root radius. Analysis of elastic stress distributions at cracks and notches and models of plasticity at crack tips are used to discuss the results. Limited evidence in the literature indicates that the fracture mechanics approach may provide a general description of crack initiation and growth in notched specimens, and a simple graphical method of calculating fatigue lives is described. The results are used to illustrate the effects of specimen size and geometry on the fatigue life of notched specimens. The relevance of the work to the assessment of the significance of defects in welds is discussed.

Implementation of pseudo J-integral based Paris’ law for fatigue cracking in asphalt mixtures and pavements

Pavement analysis and design for fatigue cracking involves a number of practical problems like material assessment/screening and performance prediction. A mechanics-aided method can answer these questions with satisfactory accuracy in a convenient way when it is appropriately implemented. This paper presents two techniques to implement the pseudo J-integral based Paris’ law to evaluate and predict fatigue cracking in asphalt mixtures and pavements. The first technique, quasi-elastic simulation, provides a rational and appropriate reference modulus for the pseudo analysis (i.e., viscoelastic to elastic conversion) by making use of the widely used material property: dynamic modulus. The physical significance of the quasi-elastic simulation is clarified. Introduction of this technique facilitates the implementation of the fracture mechanics models as well as continuum damage mechanics models to characterize fatigue cracking in asphalt pavements. The second technique about modeling fracture coefficients of the pseudo J-integral based Paris’ law simplifies the prediction of fatigue cracking without performing fatigue tests. The developed prediction models for the fracture coefficients rely on readily available mixture design properties that directly affect the fatigue performance, including the relaxation modulus, air void content, asphalt binder content, and aggregate gradation. Sufficient data are collected to develop such prediction models and the R2 values are around 0.9. The presented case studies serve as examples to illustrate how the pseudo J-integral based Paris’ law predicts fatigue resistance of asphalt mixtures and assesses fatigue performance of asphalt pavements. Future applications include the estimation of fatigue life of asphalt mixtures/pavements through a distinct criterion that defines fatigue failure by its physical significance.

Effects of crack length and shape on the fracture toughness of a high strength steel 300m

The results of fracture toughness tests on a high strength steel 300m are presented. These results show (i) that in the presence of through-thickness cracks the toughness remains constant down to (a/W)-ratios as low as 0.01 and failure loads up to 0.85σy, and (ii) that the material is more resistant to crack growth when the cracks are semi-elliptical in shape, giving a toughness value which is almost 25 per cent higher than the through-thickness one. Three independent stress analyses are used to obtain stress intensity values for the semi-elliptical cracks and additional confirmation of the increase in toughness comes from stretch zone measurements.