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

On the stress analysis of cracked bodies by means of finite elements

The work described in this thesis deals with the development and application of a finite element program for the analysis of several cracked structures. In order to simplify the organisation of the material presented herein, the thesis has been subdivided into two Sections : In the first Section the development of a finite element program for the analysis of two-dimensional problems of plane stress or plane strain is described. The element used in this program is the six-mode isoparametric triangular element which permits the accurate modelling of curved boundary surfaces. Various cases of material aniftropy are included in the derivation of the element stiffness properties. A digital computer program is described and examples of its application are presented. In the second Section, on fracture problems, several cracked configurations are analysed by embedding into the finite element mesh a sub-region, containing the singularities and over which an analytic solution is used. The modifications necessary to augment a standard finite element program, such as that developed in Section I, are discussed and complete programs for each cracked configuration are presented. Several examples are included to demonstrate the accuracy and flexibility of the technique.

An analysis of the effect of process parameters on the formability of sheet metal

This investigation is in two parts, theory and experimental verification. (1) Theoretical Study In this study it is, for obvious reasons, necessary to analyse the concept of formability first. For the purpose of the present investigation it is sufficient to define the four aspects of formability as follows: (a) the formability of the material at a critical section, (b) the formability of the material in general, (c) process efficiency, (d) proportional increase in surface area. A method of quantitative assessment is proposed for each of the four aspects of formability. The theoretical study also includes the distinction between coaxial and non-coaxial strains which occur, respectively, in axisymmetrical and unsymmetrical forming processes and the inadequacy of the circular grid system for the assessment of formability is explained in the light of this distinction. (2) Experimental Study As one of the bases of the experimental work, the determination of the end point of a forming process, which sets the limit to the formability of the work material, is discussed. The effects of three process parameters on draw-in are shown graphically. Then the delay of fracture in sheet metal forming resulting from draw-in is analysed in kinematical terms, namely, through the radial displacements, the radial and the circumferential strains, and the projected thickness of the workpiece. Through the equilibrium equation of the membrane stresses, the effect on the shape of the unsupported region of the workpiece, and hence the position of the critical section is explained. Then, the effect of draw-in on the four aspects of formability is discussed throughout this investigation. The triangular coordinate system is used to present and analyse the triaxial strains involved. This coordinate system has the advantage of showing all the three principal strains in a material simultaneously, as well as representing clearly the many types of strains involved in sheet metal work.

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.

Ductile fracture mechanisms in a structural steel

Microstructural fracture processes in a BS4360 Grade 50D structural steel with lower sulphur content were studied in smooth tensile specimen tests and Charpy-size bend bar tests. Based on the experimental analysis, an experimental void growth relation with the plastic strain and stress triaxiality and multiplying factor on void growth were determined. Experimental results show that the void growth relation can be reasonably used to estimate the constraint in the specimens containing the notch or crack, also they can be used to evaluate the variations of the stress triaxiality in front of the notch and crack tip under general yielding condition. Side-grooves obviously increase the constraint of the CVN specimens. Strain hardening leads to increasing the stress triaxiality, and decelerating the net void growth. This is especially true for the values of stress triaxiality more than about one. Additionally, the effect of the stress triaxiality on the critical void growth corresponding to the onset of ductile tearing was preliminarily investigated. In this work, a large number of smaller specimens were tested to investigate the ductile-brittle transition behaviour of the structural steel. A void growth rate explanation was suggested for evaluating the temperature transition behaviour. The elastic-plastic fracture tough-ness values based on small specimen tests, such as pre-cracked side-grooved bending specimen and short bar tensile specimen, may give large overestimates of the plane strain fracture toughness.

Shear induced toughening in bonded joints:experiments and analysis

Bonded joint specimens were fabricated from composite adherends and either an epoxy or a urethane adhesive. In mixed-mode fracture experiments, the epoxy bonded specimens generally failed by subinterfacial fracture in the composite, while specimens bonded with urethane failed very close to the adhesive/substrate interface. For the epoxy bonded specimens, fracture toughness did not change significantly with mode-mix, but for urethane bonded joints, fracture toughness increased with increasing shear load. Finite element analysis, which modeled specimens bonded with the two adhesives, showed similar trends. The different toughening behaviors for the two bonded joints can be attributed to dissipation of energy through inelastic deformation, which was insignificant in the epoxy-bonded joints but substantial when the urethane was used as the bonding agent.

Fracture and viscoelastic properties of asphalt binders during fatigue and rest periods

Viscoelastic asphalt binder plays an important role in bonding individual aggregate particles and contributes to the durability and stability of asphalt pavement. When asphalt binder is subjected to cyclic loading, deformation and fracture may develop simultaneously within it, leading to the deterioration of material properties and eventually fatigue failure. Research has found that some degree of recovery may develop if rest periods are applied after fatigue deterioration. However, it is not clear whether such recovery is caused by fracture healing, viscoelastic recovery, or both. This paper presents an analysis to differentiate the contributions of fracture healing and viscoelastic recovery to the asphalt binder during rest periods. It also evaluates the damage caused by deformation and fracture during a fatigue process. It is found that viscoelastic recovery plays an important role in the instant increase in the dynamic shear modulus at the beginning of the rest period. The effect of fracture healing on dynamic shear modulus recovery is more dominant in the long term. A healing index is developed based only on the fracture healing of asphalt binder, excluding the effect of viscoelastic recovery. It can be used to evaluate the true healing properties of different asphalt binders. Copyright © 2014 by ASTM International.

Using overlay test to evaluate fracture properties of field-aged asphalt concrete

Field material testing provides firsthand information on pavement conditions which are most helpful in evaluating performance and identifying preventive maintenance or overlay strategies. High variability of field asphalt concrete due to construction raises the demand for accuracy of the test. Accordingly, the objective of this study is to propose a reliable and repeatable methodology to evaluate the fracture properties of field-aged asphalt concrete using the overlay test (OT). The OT is selected because of its efficiency and feasibility for asphalt field cores with diverse dimensions. The fracture properties refer to the Paris’ law parameters based on the pseudo J-integral (A and n) because of the sound physical significance of the pseudo J-integral with respect to characterizing the cracking process. In order to determine A and n, a two-step OT protocol is designed to characterize the undamaged and damaged behaviors of asphalt field cores. To ensure the accuracy of determined undamaged and fracture properties, a new analysis method is then developed for data processing, which combines the finite element simulations and mechanical analysis of viscoelastic force equilibrium and evolution of pseudo displacement work in the OT specimen. Finally, theoretical equations are derived to calculate A and n directly from the OT test data. The accuracy of the determined fracture properties is verified. The proposed methodology is applied to a total of 27 asphalt field cores obtained from a field project in Texas, including the control Hot Mix Asphalt (HMA) and two types of warm mix asphalt (WMA). The results demonstrate a high linear correlation between n and −log A for all the tested field cores. Investigations of the effect of field aging on the fracture properties confirm that n is a good indicator to quantify the cracking resistance of asphalt concrete. It is also indicated that summer climatic condition clearly accelerates the rate of aging. The impact of the WMA technologies on fracture properties of asphalt concrete is visualized by comparing the n-values. It shows that the Evotherm WMA technology slightly improves the cracking resistance, while the foaming WMA technology provides the comparable fracture properties with the HMA. After 15 months aging in the field, the cracking resistance does not exhibit significant difference between HMA and WMAs, which is confirmed by the observations of field distresses.

Improved methodology to evaluate fracture properties of warm-mix asphalt using overlay test

This study developed a reliable and repeatable methodology to evaluate the fracture properties of asphalt mixtures with an overlay test (OT). In the proposed methodology, first, a two-step OT protocol was used to characterize the undamaged and damaged behaviors of asphalt mixtures. Second, a new methodology combining the mechanical analysis of viscoelastic force equilibrium in the OT specimen and finite element simulations was used to determine the undamaged properties and crack growth function of asphalt mixtures. Third, a modified Paris's law replacing the stress intensity factor by the pseudo J-integral was employed to characterize the fracture behavior of asphalt mixtures. Theoretical equations were derived to calculate the parameters A and n (defined as the fracture properties) in the modified Paris's law. The study used a detailed example to calculate A and n from the OT data. The proposed methodology was successfully applied to evaluate the impact of warm-mix asphalt (WMA) technologies on fracture properties. The results of the tested specimens showed that Evotherm WMA technology slightly improved the cracking resistance of asphalt mixtures, while foaming WMA technology provided comparable fracture properties. In addition, the study found that A decreased with the increase in n in general. A linear relationship between 2log(A) and n was established.

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.