Scholarly article on topic 'On cold dwell facet fatigue in aero-engine components'

On cold dwell facet fatigue in aero-engine components Academic research paper on "Materials engineering"

Share paper
Academic journal
International Journal of Fatigue
OECD Field of science
{"Cold dwell fatigue" / "Crystal plasticity" / "Titanium alloys" / "Aero-engine discs"}

Abstract of research paper on Materials engineering, author of scientific article — M.A. Cuddihy, A. Stapleton, S. Williams, F.P.E. Dunne

Abstract This paper investigates the mechanisms of facet nucleation through combining aero-engine manufacturer disc component test data with microstructure-sensitive crystal plasticity finite element (CPFE) models. Full-scale component testing has been carried out in a manner representative of in-service conditions. Elastic FE analyses of discs under these conditions and fully accounting for thermal and residual processing strains have also been carried out. Disc facet nucleation sites have been identified and the local stress states evaluated in order to establish crystal plasticity oligocrystal sub-models. The oligocrystal RVE models provide knowledge of hard-soft grain stresses under dwell loading, and the consequent load shedding in order to provide stresses required for the facet nucleation. The disc component facet observations together with the crystal plasticity sub-model oligocrystal approach provide persuasive evidence that a hard-soft grain combination is required for facet formation, that the remote stress state influences the resolved shear stress on the soft grain initiating slip (with tensile uniaxial stress state more damaging than a tension-tension biaxial stress state), and that the load shedding which results is essential in pushing up the hard-grain basal stress to nucleate facets.

Academic research paper on topic "On cold dwell facet fatigue in aero-engine components"

Accepted Manuscript

On cold dwell facet fatigue in aero-engine components

M.A. Cuddihy, A. Stapleton, S. Williams, F.P.E. Dunne



S0142-1123(16)30405-4 JIJF 4149

To appear in:

International Journal of Fatigue

International Journal of


Received Date: Revised Date: Accepted Date:

19 August 2016 25 November 2016 28 November 2016

Please cite this article as: Cuddihy, M.A., Stapleton, A., Williams, S., Dunne, F.P.E., On cold dwell facet fatigue in aero-engine components, International Journal of Fatigue (2016), doi: 2016.11.034

This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

On cold dwell facet fatigue in aero-engine components

M.A. Cuddihy1*, A. Stapleton2, S. Williams2, F.P.E. Dunne1

1. Department of Materials, Imperial College London, United Kingdom

2. Rolls-Royce plc, Derby, United Kingdom


This paper investigates the mechanisms of facet nucleation through combining aero-engine manufacturer disc component test data with microstructure-sensitive crystal plasticity finite element (CPFE) models. Full-scale component testing has been carried out in a manner representative of inservice conditions. Elastic FE analyses of discs under these conditions and fully accounting for thermal and residual processing strains have also been carried out. Disc facet nucleation sites have been identified and the local stress states evaluated in order to establish crystal plasticity oligocrystal sub-models. The oligocrystal RVE models provide knowledge of hard-soft grain stresses under dwell loading, and the consequent load shedding in order to provide stresses required for the facet nucleation.

th the cry

The disc component facet observations together with the crystal plasticity sub-model oligocrystal approach provide persuasive evidence that a hard-soft grain combination is required for facet formation, that the remote stress state influences the resolved shear stress on the soft grain initiating slip (with tensile uniaxial stress state more damaging than a tension-tension biaxial stress state), and that the load shedding which results is essential in pushing up the hard-grain basal stress to nucleate facets.

esults is <

Keywords: cold dwell fatigue; crystal plasticity; titanium alloys; aero-engine discs

Email address:

1 Introduction

In-service conditions of an operating gas turbine engine represent some of the most complex loading regimes of any thermo-mechanical system. Titanium alloys are widely used in gas turbine engine often under extremes of loading, in components such as discs and blades. They are ut extensively in the aerospace industry due to their low density, excellent corrosion resistance and high fatigue strength [1]. However, it has been well recognised over the past 40 years that titanium alloys suffer from a significant failure mode known as cold dwell fatigue, often characterised by facet crack nucleation and subsequent growth. The lifetime reduction resulting from the inclusion of load holds at maximum stress, and at ambient temperature (<200°C), during each loading cycle is termed the dwell fatigue life debit. It continues to be a significant industrial concern because its non-destructive evaluation management is hugely costly, and it remains a safety-critical issue. The compressor discs used in the intermediate pressure section of a jet engine are among the components deemed to be 'life-critical', meaning that rupture or failure is not an acceptable outcome due to risk to the airframe integrity; the required reliability for these components is 1 event in 108 flights [2].

Considered fundamental to cold dwell fatigue is the formation of facets which are always associated with the failure initiation site. The commonly utilised model for facet nucleation is an adaption of the classical Stroh model for crack initiation due to dislocation pile-up [3], introduced by Evans and Bache [4]. An important subsequent modification to this model included the incorporation of the dwell component, by Hasija and colleagues [5]. The time-dependent stress redistribution, or load shedding, is now believed to be key to the formation of facets and this is recognised in this work. The facet nucleation process occurs on basal planes and is argued to preferentially develop on a grain badly orientated for slip (hard grain), with c-axis near or parallel to the local maximum principal stress direction, adjacent to a soft grain well orientated for slip, be that occurring on prismatic or basal systems. The stress acting normal to the basal plane (and, as some authors argue, the shear stress) during the dwell period, leading to load shedding, is thought to be necessary (but not sufficient) for facet nucleation [6], [7].

Evans and Bache [4] were the first to explicitly attempt an investigation of a multiaxial dwell loading regime. They performed tension, torsion, and combined tension-torsion experiments and were able to make a comparison of effective stress at failure and lives-to-failure (Nf). For equal values of Von Mises effective stress, it was observed that tension-torsion tests (i.e. the inclusion of an in-phase shear loading component) often had a higher lives-to-failure value than simple uniaxial tension. However, as Doquet and De Greef point out [8], shear loading resulting from torsion is not particularly representative of aero engine in-service conditions so these findings are not directly comparable with

the dwell regime. A study carried out by Doquet and De Greef [8] consisted of a series of tests comparing uniaxial and tensile equibiaxial loadings in low cycle fatigue (LCF) and dwell fatigue. The results obtained, while providing valuable insight, were subject to quite considerable scatter, but the indication from the biaxial loading testing in low cycle fatigue was that tensile equibiaxial loading conditions have an ameliorating effect, resulting in extended fatigue life. In many of the dwell fatigue tests performed, the specimens ruptured prematurely, precluding conclusive interpretation. Doquet and De Greef attempted to rationalise their findings in light of the classic multiaxial theories (such as Fatemi-Socie, Sines and Crossland etc.) but found that none of these was able to accurately capture the behaviour observed.

In further LCF (as opposed to dwell) work, Bonnand et al. [9] tested a range of Ti-6Al-4V samples in combined tension-torsion and equibiaxial tension under LCF and found that few of the theories predict, with confidence, observed behaviour. They developed an ad-hoc theory which represents a compromise between the two major categories of multiaxial fatigue criteria; namely critical plane theory and effective stress amplitude. A review of the myriad multiaxial stress theories has been presented by Kallmeyer et al. [10] who conclude that the Findley parameter and the Manson-McKnight model present the most convincing fit in a survey of a large dataset of experimental LCF findings. These data are, however, for low cycle fatigue only and for non- hexagonal close packed (hcp) metals, so not directly applicable in the context of cold dwell fatigue.

In passing, we note that multiaxial stress states are often associated with ductile fracture. Rice and Tracey [11] considered spherical voids and developed a fundamental description of ductile fracture based on void growth and coalescence. The Rice and Tracey model was employed by Helbert et al. [12] who undertook a detailed investigation of the effects of stress triaxiality on the damage mechanisms of an equiaxed a/p Ti-6Al-4V. They found that in the low triaxiality range (values less than 0.66) no significant growth of voids is observed, even up to high plastic strains (sp < 0 . 3 ) . This point has been demonstrated again in more recent work on notched samples of IMI834 by Kumar et al. [13].

The Stroh model is not the only proposed explanation for dwell fatigue in the literature. Lefranc et al. [14] investigated dwell fatigue in Ti-6242 with a bimodal a/p microstructure and argued that crack nucleation was caused by the coalescence of cavities or voids, induced by shear stresses. These voids nucleated at the a/p interfaces and were observed to be greater in size and number under dwell fatigue loading, as opposed to static creep or normal fatigue conditions. Shear induced porosities causing micro cracks under dwell loading were also observed by Gerland et al. [15]. These authors assert that it is the porosities developed under creep conditions of the load hold that are responsible for dwell

fatigue, which is in keeping with classic ductile fracture theories [11]. The loading applied to generate this failure mechanism was close to or exceeding macroscopic yield (00.2%) and the failure may therefore not be mechanistically the same as that occurring at lower in-service stresses in aero-engine discs.

It is known from industrial experience that under in-service conditions the stresses carried by the titanium alloy components (compressor discs) are far lower than yield, at approximately 50-70% of the yield stress. Small-scale laboratory dwell testing is typically found to require much higher stresses in order to induce dwell fatigue failure. Importantly, it is noted that in the majority of cases of failure resulting from dwell fatigue observed in industrial testing, the "worst case" grain orientations are always present (see later). Laboratory samples may consist of relatively small numbers of grains across the gauge section and may not capture the statistics of full-scale components. In such samples, the worst case (rogue) combination of grain crystallographic orientations [16] may simply not be present; a compressor disc contains a volume of material several orders of magnitude larger than that for a simple tensile specimen. The high stresses required to generate facets in the absence of a rogue pair may also be the cause of other undesirable failure states. Bache et al. demonstrate that in the absence of the worst case grain orientation, crack nucleation may be forced onto other combinations, again at high applied stresses [17]. Pilchak et al. [18] found that facets may occur on basal planes, mis-orientated with the macroscopically applied load by approximately 43°. However, this test was also carried out at the relatively high stress of 0.95 cy.

This paper presents an investigation on facet nucleation in a-titanium alloy discs from a comprehensive industry spin test programme in such a way to unambiguously avoid the small volume, high stress problem associated with laboratory dwell fatigue testing. We begin by considering the stress states in a rotating disc, establishing the basis for a further in depth study of multiaxial stress loading. This is followed by a detailed examination of stress states and facet nucleation results from disc spin tests carried out by an aero-engine manufacturer. This continuum level study indicates the necessity of a detailed crystal plasticity sub-modelling approach which is used to illustrate the effects of multiaxial stress state on the rogue grain combination and facet nucleation. The nature of the hard grain crystallographic orientation with respect to the principal stress is briefly assessed in the context of experimental observations reported in the literature. The macro-level stress state analysis of the disc is then coupled with the crystal plasticity analyses which allow us to investigate more fully the connection between load shedding and facet nucleation in discs under dwell fatigue. This facilitates prediction of the likely sites of facet nucleation and the stress levels and states necessary to cause them.

2 Disc spin test studies

2.1 Introduction to spin tests

In the absence of predictive capabilities due to the complexities of structure, texture and properties, aero-engine manufactures have had to rely on an expensive (> £150,000 each) set of time-consuming (5-8 weeks per test) spin tests to establish a reliable applied stress-to-life (S-N) curve [2]. Spin tests are a form of quality control which involve taking full sample discs manufactured as standard product hardware with fully representative geometric features, which are cycled in a manner similar to inservice use, often until failure.

Rolls-Royce have undertaken a comprehensive spin test programme for a variety of in-service components and one particular set of results is examined here. This set is defined as a series of spin

Figure 1 Compressor disc geometry from spin tests, (a) cutaway of a 3D rendering of a disc, and (b) axisymmetric slice of

disc highlighting the key bore region with facet locations

The illustration in Figure 1(a) shows a 3D rendering of the disc using the geometry supplied by the manufacturer (with a 1/4 cutaway for clarity). The inset figures illustrate an axisymmetric slice of the disc and show the facet locations. These facet locations, identified by the black spots in the bore

region, have been measured and quantified as the failure origins for the various dwell cycle tests. These sites are collectively similar, such that they represent the qualities associated with classical facet behaviour seen in dwell fatigue, i.e. basal plane orientation approximately perpendicular to the principal stress direction (hoop), and quasi-cleavage fracture planes.

2.2 Stress distributions in compressor disc spin tests

An axisymmetric finite element (FE) model with isotropic, elastic material properties representing the compressor disc material (an a-Ti variant, Ti-834) is utilised to determine stress states for disc spin tests. The model replicates the spin test temperature (<200oC, known to be the range conducive for dwell) and assumes an isothermal state. Interestingly, a temperature of 120oC is argued to be the worst case for dwell fatigue [19]. The model is loaded to replicate the typical cycling of the spin test, i.e., rotational loading at minimum RPM and spinning up to a peak rotational speed which is held for 2

minutes (dwell condition), after which the speed is reduced back down to the idle condition.

Prior to spin loading, an initial step is carried out to include the residual stresses developed by the post-forging heat treatment, aging, and quenching operations. The thermal strains generated by the quenching process are calculated in a separate model of the bulk forged element, and the resulting strains are mapped onto the rig test model. The calculation of the thermal strains induced by the quenching process is verified against a measurement of the strains by an incremental hole-drilling test1. This involves bonding standard copper strain gauges to the specimen, then drilling a hole through the strain gauge and the material. As the material relaxes around the hole, the displacements measured by the strain gauge can be used to back calculate the surface residual strains [20]. This technique is only suitable for near-surface measurements of strain. Alternative methods of measuring the process residual strains at depth have been developed using synchrotron x-ray diffraction; this technique has been used to measure residual strains in nickel superalloy turbine discs and may probe up to ~50mm beneath the surface [21].

The contour plots in Figure 2 show the principal normalised stress distributions, corresponding to the hoop, radial and axial stress directions respectively (there are minimal amounts of shear observed in these models) at full spin loading, and including the initial residual stresses. These operating stresses have been normalised by a representative tensile strength. Figure 2 also shows three representative facet locations.

1 The residual strain calculation was carried out at Rolls-Royce plc

Figure 2 Distributions of normalised hoop (1), radial (2) and axial (3) principal stresses; values indicated on the isolines (except for G3). Example facet nucleation sites shown by blue circles

Evidently, the hoop stresses dominate over the other two directions, and the peak concentrations are localised in the lower left region of the disc bore. However, relatively large radial stresses are carried by the diaphragm section of the disc where the radial stress is almost 40% of the hoop stress for that region.

Figure 3 shows a close up and more detail of the hoop stresses with respect to the 22 experimentally observed facet nucleation sites. An immediate observation is the strong dependence on the hoop stress i.e. no faceting is observed to occur when the normalised hoop stress is less than 0.58 with the majority of facets occurring at hoop stress values greater than 0.66.

Figure 3 Detailed assessment of facet location sites wit]

ormalised hoop stress distribution

A simple design paradigm may be to assume that peak hoop stress is the dominant driver of facet nucleation, but the facet outliers (e.g. a13, a1) cannot be explained in this way. In addition, why is there just the one observed facet location on the 0.58 stress contour? For completeness, this is examined further at the continuum level by determining the stress triaxiality ratio for the bore region, which is shown in Figure 4. The stress triaxiality x, is the ratio of the hydrostatic mean stress and the

von Mises effective stress given by

Figure 4 shows that the facet nucleation locations fall within a quite tight banding of triaxiality ratio, with average triaxiality of ~0.34, which represents an approximately uniaxial tensile stress state. The facet point a13, previously considered an outlier with respect to hoop stress, is now observed to have the same stress triaxiality as the majority of the other facet locations.

= V3 Ol +

h) +{a2~ X = Он/°е

г) +{a3~ai)

Figure 4 Bore section contour plot showi

triaxiality and facet nucleation locations

However, other bore locations with equivalent triaxiality values on the right-hand side of the bore section (indicated in Figure 4) do not show any incidences of facet nucleation such that a triaxiality of X = 0. 3 4 is clearly not a sufficient condition for faceting. Finally, in contrast with the earlier discussion on stress triaxiality related to laboratory specimens, high (tensile) triaxiality is definitively not associated with dwell facets. The highest value of X, 0.55, is far removed even from the closest facet location, a13. It is important to note that the operating stresses in the bore region are much lower than the ca2% yield stress, such that classical multiaxial continuum stress theory or ductile fracture models are simply not relevant. Hence, the next section seeks to provide mechanistic understanding through the use of crystal plasticity sub-models which replicate the location-specific crystallography, stress states and load shedding resulting from the spin test loading.

Crystal plasticity modelling of biaxial loading in dwell fatigue

The previous section demonstrated that multiaxial stress states develop within the disc, and that while the hoop stress due to rotational speed is large, it is not the sole component of stress. Doquet and DeGreef s work [8] indicates that equibiaxial tensile loading may have a positive effect on fatigue life, as opposed to uniaxial tension. Secondly, the torsion testing work by Evans and Bache [4] demonstrates that shear loading affects fatigue life, however it is not a directly representative loading for disc components, so it will not be considered further. Finally, the examination of disc spin data above suggests a link to the multiaxial nature of the stress state. Hence, an approach is presented based on crystal plasticity sub-modelling which provides the link between facet nucleation at the microstructural level with the continuum disc stress states developed under in-service loading. With knowledge of the facet nucleation sites within the spin tested discs, the macroscale stress states determined above at these locations are imposed at the sub-model crystal plasticity level, as indicated schematically in Figure 5.

Figure 5 Schematic of sub-modelling procedure to calculate facet stress state

For each facet nucleation site, a multigrain representative volume element (RVE) is developed, each of which contains the classical hard-soft grain orientation combination conducive to dwell load shedding. The macroscale stress state at each facet location is extracted from the disc analysis above

and imposed locally at the sub-model crystal plasticity level. However, before presenting disc facet nucleation results, we briefly discuss the crystal plasticity framework utilised in the analyses, and carry out systematic studies of the role of biaxial stress states in load shedding, and the effect of the orientation of the hard grain.

3.1 Crystal plasticity framework

The rate-dependent crystal plasticity framework utilised in this study is based on that originally developed by Dunne et al. [22]. We begin by stating that the multiplicative kinematic decomposition of the deformation gradient (F), into elastic (Fe) and plastic (Fp) tensors is given by Lee [23] such that:

given b

F = — = FeFp . ^ 4

Assuming crystallographic slip accounts for plastic deformation, and for now considering just single slip,

Fp = / + ^P = / + y( s®n ) 5

- = / + r(»®n)

normal res

in which s and n are slip direction and plane normal respectively, and y is the magnitude of the slip. Since the material properties at a point are time-dependent, it is convenient to write the plastic deformation gradient in rate form as

Fp = Y ( s ® n ) 6

For a spatially varying velocity field, the velocity gradient is defined and decomposed into the symmetric and anti-symmetric parts as

L = F (F ) " 1 = sym (L) + asym (L) = D + W 7

where D and W are the deformation rate and spin tensors, respectively. Both D and W may be further decomposed into elastic and plastic parts as:

D = D e + D p 8

W = We + Wp 9

The rate of crystal plastic deformation is

D p = sym (L p) 10

and the plastic spin

Wp = asym (Lp) 11

The plastic velocity gradient is that associated with plastic flow through a fixed lattice, including contributions from all slip systems, and is given by

LP = ls l0n l 12

with normal vectors ri and slip direction vectors s' corresponding to the i' slip system. The slip rate, yl, is computed according to a slip rule, described next.

Dislocation pinning is taken to occur through the presence of lattice obstacles, which may include solute atoms, the sessile statistically stored dislocations (SSDs) and their associated structures, as well as geometrically necessary dislocations (GNDs), incorporated in the slip rule as an overall obstacle density. The resultant slip rule (defined explicitly in Dunne et al. [22]) determining the slip rate of a given slip system takes the form:

0, It'Ic-To

Yl = { -2 I AF\ (( I tl I -t0-xl) AVl\ . .. 13

" ! p™bl v exp (-- —j s inh i I 1 1 00- J-1 I t l j>(x + to )

in which p™ is the mobile SSD density, p0 the initial sessile dislocation density, b' the Burgers vector magnitude for slip system ', v the frequency of attempts (successful or otherwise) by dislocations to jump the energy barrier, A F the Helmholtz free energy or activation energy, k the Boltzmann constant, T the temperature in Kelvin (K), t' the resolved shear stress, xl the back-stress on slip system ', r0 the

size-independent critical resolved shear stress, y0 the shear strain that is work conjugate to the resolved

■2 1

shear stress, and AVthe activation volume, A Vl = y0 , whe re I = -¡=.

The density of sessile statistically stored dislocations is allowed to accumulate in proportion to the accumulated slip, p with hardening factor, y', chosen to ensure the experimentally observed hardening is reproduced.

Pssd = y'p

t+At t 14

Pssd = Pssd + Psso^t

This simple relation, with only one parameter, is chosen as it allows for adequate representation of experimental data. It is well understood that SSD evolution is a function of plastic slip, however, the precise nature of this relationship is less clear (and likely varies with temperature, strain rate, and material system). There remain limited quantitative experimental techniques and studies of the SSD evolution with strain, as opposed to GNDs where EBSD based methods have shown excellent results [24], [25]. We argue that incorporation of a more complicated evolution rule (such as the Kocks-

shown ant for

Mecking rule) with additional fitting factors for dislocation density without further mechanistic understanding does not add much value. GNDs are not directly considered in this work as the average grain size in this study is considered to be ~500^m likely representative of the effective structural unit size of the industrial material with respect to dwell fatigue. Previous experimental studies have show that length scale effects and strain gradient associated dislocations are likely to be significan much smaller length scales (<50^m) [26].

Each slip system becomes active when the resolved shear stress is equal to or greater than the combined effects of the back stress and the size-independent critical resolved shear stress ( | t 1 | > t0 ). The back stress on a slip system is described by the following relation:

............. . „ 15

Note that: a =


xl = aGbjptotal = aGbJpcND+PssD , and which is approximately 0.5 for titaniu

]us user-d

The crystal model is implemented within a combined Abaqus user-defined element (UEL) [27], which facilitates strain gradient calculations, and user defined material property (UMAT) subroutine.

3.2 Material calibration of slip rule

This work focuses on the cold dwell susceptible, near-alpha phase, titanium alloy Ti-6Al which is hexagonal close packed (hcp), and has 24 slip systems. Due to the level of aluminium in this alloy it is assumed that there is no twinning and that crystallographic slip is the only source of inelastic deformation [1]. The <a> and <c+a> type systems are differentiated (c/a ratio of 1.593 is used), and the slip systems for a single hcp crystal are shown in Figure 6.

Figure 6 Hexagonal close packed slip systems for a-phase titanium

Hasija et al. [5] conducted mechanical tests with samples of single-phase a-Ti-6Al, for both single crystal and polycrystalline samples. The single crystal tests were carried out at constant strain rate, for basal and prismatic slip <a>, and pyramidal slip <c+a>, and the stress strain curves generated by these tests when coupled with a CPFE single crystal model allow us to establish the base critical resolved shear stresses, see Figure 7(a). Also, as hardening is hardly observed in the tests, the hardening coefficient y' in the SSD hardening rule is set to a small non-zero value of 0.05. Equal strengths are assumed here for basal and prismatic <a> slip systems for simplicity, since the difference between them has been observed to be small in the experiments of Hasija et al. The slip strength of pyramidal <c+a> systems is often observed to be about three times that for basal or prismatic systems, recently demonstrated experimentally by Gong et al [28].

G 1600


X Pyramidal - Experimental A Prism - Experimental

-Pyramidal - Simulation

-Prism - Simulation







0.002 0.004 0.006

Plastic Strain (e^) (a)

ntal and C

0.000 0.002 0.004 0.006 0.008 Plastic Strain (ep^)

Figure 7 (a) Experimental and CPFE model curves for prismatic <a> and pyramidal <c+a> slip system response based on single crystal creep experimental data, and (b) Strain rate sensitivity calibration results, (ABQ refers to Abaqus simulations

and OSU refers to experimental data

The polycrystal experiments considered a compression test of three samples at three strain rates, 8.4 x

10-4, 1.5 x 10-4 and 1.68 x 10-5 s-1, and have been used in combination with a representative polycrystal CPFE model [19] to determine the activation energy, AF, relating to strain rate sensitivity (see Eqn. 13). The polycrystal model with a chosen activation energy of AF = 9.913x10-20 J/atom is shown alongside the three experimental strain rates indicated in Figure 7(b), providing a good representation of the polycrystal average stress-strain response. The activation volume is established to be ~18b3, which is consistent with Conrad's [29] range of values for a-Ti; ~8 to 80b3.

The small amount of strengthening observed is captured simply by means of the contribution from grains badly orientated for slip within the polycrystal, such that the hardening factor, y ' remains as would be expected from the single crystal tests. Latent hardening is considered for the polycrystal

models such that hardening on the f slip system increases the slip strength of all other systems (incorporated in the slip rule via Eqn. 15). Elastic moduli (from Hasija et al. [5]) and all the derived slip rule properties for the representative Ti-6Al alloy are given in Table 1.

Slip rule properties Value Elastic constants Value

pmo b í le 5.0 |im"2 E Exy 85 GPa

p i n itial 0.01 |m-2 E Eyz 120 GPa

p i n itial 0.0 |m-2 Vxy 0.46

V 1.0 x 1011 Hz Vyz 0.22

b <a> 3.2 x 10-4 |m Gxy 29 GPa

b <c+a> 5.1 x 10"4 |m Gyz 40 GPa

AF 9.913 x 10-20 J

k 1.381x1023 J K-1 L*1

T 293 K

AF 18 b3

v' 0.05

Table 1 Material property data

These derived slip rule constants have been utilised in a previous study, carried out by some of the authors, on the temperature sensitivity of cold dwell fatigue [19]. CPFE polycrystal models employing these slip rule parameters have been shown to accurately predict the well-known "switching off' of the dwell effect for temperatures in excess of ~200°C.

3.3 CPFE modelling of biaxial dwell stress states

We consider a representative hcp polycrystal sample with fixed in-plane ligament width of 8000^m and height of 7500^m, comprising 240 grains. This model is subject to a combination of loads used to represent the two major stress components as observed from the disc operating stress model and indicated in Figure 8(a). The lateral load, referred to as the secondary load, is representative of the radial stress and is applied in-phase with the primary vertical load - representing the hoop stress. By testing a range of loadings from uniaxial to full equibiaxial states, we can demarcate precisely the influence of the additional load component on the dwell condition.

Figure 8 CPFE polycrystal model: (a) polycrys

(b) the region surrounding the rogue grain pair

This study addresses microstructural quantities thought to be the most sensitive indicators of local behaviour. Particularly, we consider a combination of two grains buried within the polycrystal orientated such that one is favourable for slip and the other is not; this grain couple is indicated in Figure 8(b) and is where key stress magnitudes are extracted for analysis. Furthermore, grains have a directionally solidified structure and rectangular morphology, giving prismatic grains, and all grains have the same depth in the out of plane direction. Consideration is given to mesh refinement in the region local to the hard/soft pair with 768 elements per grain, whereas the mesh for the other grains is relatively coarse.

The nature of 'realistic' grain morphology has been addressed by Sauzay and co-workers [30], who have found that it contributes little to the overall mechanistic behaviour with respect to the enormous computational expense in mesh generation, thus motivating the use of simple geometry. Apart from the hard / soft grain couple with specific misorientation, all other grains in the model are assigned pseudo-random orientations, but with the imposed constraint of having c-axes inclined perpendicularly to the primary remote loading direction (within +/-10°) in order to favour easy slip2.

2 Euler angles will be in the set: (p ±,(p,0 2 ) = (p 1 = 0 ,- 1 0 < p < 1 0,- 3 60 < p2 < 3 60 )

Pseudo-random here refers to the generation of Euler angles via a random number generator. The material properties and slip rule constants are those described in §3.2.

We examine the results of a systematic study of the effects of the addition of a secondary applied stress, increasing from 0 to 700MPa; that is, from purely uniaxial to equibiaxial loading. Therefore, each analysis has a fixed primary load of 700MPa, in the yy-direction, with a differing secondary load, in the xx-direction, from 0 to 700MPa. These loads are reached after 12s and held for a dwell period of 120s. It is during the dwell that stress redistribution potentially occurs from the soft grain to the adjacent hard grain (the load shedding), and in this study, the biaxiality of the stress state and its influence on the load shedding is examined. One typical dwell cycle is chosen for analysis, largely for reasons of computational efficiency, while noting however that there is anecdotal and some physical evidence, which posits that damage may occur in the first cycle of loading, and evolve slowly in subsequent cycles [19].

We choose a point in the hard grain, just inside the boundary with the soft grain and examine the stress in the yy-direction. In this study the hard grain is orientated such that the c-axis is parallel with the primary load and the hcp unit's basal plane is therefore perpendicular to this stress. Therefore we

m-ivit* +/-» +Ulf1 n n on Q n n ' (T1

Figure 9 (a) Range of various load histories with in-phase primary (ap) and secondary (as) loads, (b) abasal for a point inside the grain boundary of the soft / hard grain across loading period (as indicated in the schematic)

0 6 12 18 24 Time {s)

The results shown in Figure 9(b) demarcate clearly the differences arising from the various load histories for the differing stress states. Firstly, we note the disparity between the two extreme cases for which it is seen that after a two minute dwell period, the difference in stress accumulated by the hard grain is almost 350MPa. Focusing on the uniaxial load, we see that at the end of the dwell period the basal stress is approximately 1150MPa, which is close to the commonly accepted value for tensile failure in representative titanium alloy with strong hcp c-axis crystallographic texture. The results compare qualitatively with the work of Doquet and DeGreef which suggests that equibiaxial loading in dwell fatigue would show longer fatigue life. This may be understood from the crystal plasticity model through the apparently far lower accumulation of residual stresses in the hard unit from cycles of equibiaxial dwell loads, as opposed to uniaxial loading. The results in Figure 9(b) are utilised further in the disc spin test discussion presented later.

3.4 Effect of hard grain orientation with respect to principal loading axis

It has been well established, through experimental observation, that facet fatigue nucleation is almost always associated with fracture on a basal plane which is normally orientated perpendicular to the applied loading. However, several studies have shown that facets can often be upwards of 10° to 15° in misalignment with the principal stress direction [31]. The figures below show a series of polycrystal analyses in which the hard grain c-axis is rotated from its hardest orientation, parallel to the applied stress, to perpendicular to the applied stress. These analyses are stress hold tests, with

ider the distribution of accumulation of plastic slip in the region

uniaxial loading only. We first cons close to the soft / hard pairing.

Rotating hard grain to soft orientation

Accumulated plastic strain

000e-02 583e-02 167ê-02 750e-02 333e-02 917e-02 500e-02 083ê-02 667e-02 2S0e-02 334e-03 168e-03 000e-06


Figure 10 Accumulated plastic strain contour plots, (a) to (d), for various rotations of the hard grain, from 0° to 90° with the

direction of loading indicated

When the hard grain is aligned parallel with the applied stress as in Figure 10(a), we see the development of strong slip bands around the hard unit and even a shielding effect is provided for by the hard grain, just below and to the right. Plastic strain development is progressive and preferential, and this is demonstrated by the intensity of the slip bands. The hard grain is rotated by 30° in (b) and by 45° in (c) and with this rotation, a diminution of the intensity of the slip bands around the grain may be seen with some slip activity occurring at the grain boundary of the now former hard unit.

When the hard grain is rotated through to 90° the former direction, it is seen to carry plastic strain. This allows for slip bands to form through the rotated hard unit and for plastic strain to be generated in the formerly shielded region.

Figure 11 Typical hexagonal close packed unit under arbitrary traction

As dwell facets are always characterised by the stark anisotropy between soft and hard grains, it is of interest to investigate the stress on the basal plane as a function of hard grain orientation, which is considered next. The stress normal to the basal plane illustrated in Figure 11 is given in terms of the traction and the normal direction by:

°b asa I = ( t -n) = ( a 7 ) .7 16

in which is the slip plane normal vector for the particular orientation of the hcp crystal. of the hard grain from the 'hard' orientation (parallel to the principal tensile orientation (perpendicular to the load) is indicated by the diagram in Figure 12

. The rota load) to th

ation e 'soft'

Figure 12(a) shows the effect of hard grain orientation with increasing biaxial stress ratio (a2 / -> 1). These stress values are for the consistent location in the hard grain and are taken at the end of a 12s load hold. It is observed that increasing the biaxiality reduces the stress not only for the hardest orientation, but also for the intermediate configurations (i.e. ) also.

Figure 12 Hard grain orientation and stress state effect for (a) biaxial loading and, (b) crystal c-axis parallel with principal

uniaxial stress loading

Figure 12(b) shows the effect on basal stress of rotating the hard grain through 90° under uniaxial loading conditions; interestingly the stress on the basal plane remains relatively flat until an orientation of about 15-20°, after which it falls off considerably. This finding may explain why facets may be found on basal planes which are not uniquely normal to the principal stress direction but often

orientated away from this axis, as significant normal basal stresses may be developed for c-axis orientations up to 20o off the principal stress direction.

Analysis of disc component results: combining CPFE modelling with spin test studies

In this section, we return to the sub-modelling approach developed above and shown in Figure 5 in order to link microstructural aspects of facet nucleation with the disc component spin test results. For context, the dependence of single crystal resolved shear stresses on biaxial stress states is considered first and is followed by an assessment of the resolved shear stresses acting on prism slip systems in a hard-soft grain combination located at each of the experimentally observed facet nucleation sites in the spin test data. The sub-modelling approach is presented next in which each spin disc facet location is represented by an oligocrystal containing the hard-soft grain combination subjected to the local disc multiaxial stress state and dwell loading. The consequences of the stress state, the dwell, and the presence of a hard-soft grain combination on stress redistribution, load shedding, and basal stresses on the hard grain are all assessed in the light of the nucleation of facets at these corresponding locations in the disc spin data.

spect to t

A single crystal is considered first, orientated with respect to the biaxial stress state shown in Figure 13. The resolved shear stress is calculated (simply by applying Schmid rule ) on all active slip systems and the highest value (normalised with respect to the threshold value for Ti-6Al) is shown against the biaxial stress ratio (i.e. the ratio of the two orthogonal stresses, c2/ci=0 corresponding to uniaxial loading and c2/ci=1 to equibiaxial loading).

Figure 13 Normalised single crystal resolved shear stress as a function of principal stress ratio

The single crystal is loaded with an increasing lateral load, and a fixed primary load. From Figure 13 it is clear that with increasing biaxiality the resolved shear stress on the soft grain prism system decreases. Naturally, the imposed biaxial stress state affects the plasticity of initiation of slip in the soft grain, which in a soft-hard grain combination is what is required to drive load shedding under dwell fatigue conditions. The absence of slip in the soft grain precludes stress redistribution to the hard grain.

Returning to the compressor disc geometry in Figure 1, the stress state at each facet nucleation site is extracted and used to determine the resolved shear stress which would develop on a prism slip system occurring at that same location, on the basis that it is known that a soft-hard grain combination (with respect to the disc hoop stress direction) exists at each facet site in the disc. For the purposes of visualisation, the same calculation is carried out everywhere in the disc section, and the corresponding resolved shear stresses developed are shown in Figure 14.

igure 14 Normalised resolved shear stresses at the experimentally observed facet sites (and calculated everywhere with respect to a soft grain prism slip system) for the disc bore section

This rather simple approach, which recognises the role of the hard-soft grain combination for facet nucleation, and takes full account of the stress state in determining the resolved shear stress on the soft grain prism slip system, provides a rational demarcation of the disc regions in which facets are seen to nucleate and those in which they are not, and the demarcation is given by a normalised resolved shear stress of trss/tcrss=1.0 where TCRss=280MPa. Compared directly with the disc bore triaxiality distribution of Figure 4, it may be seen that the region to the right hand side of the bore

section, whilst having areas of low triaxiality comparable to that at the facet sites, has very low normalised resolved shear stress. This reinforces the importance of the hard-soft grain combination and the role of slip in the soft grain in dwell facet nucleation and demonstrates why triaxiality in its own right is not a good indicator of propensity for facet nucleation. We note that all disc facet locations are located in regions where the normalised resolved shear stress is 1.0 or higher, indicating that slip on softly orientated grains is a prerequisite for facet nucleation. Disc facet sites are found to cluster in the region of high resolved shear stress. So far, the role of dwell and load shedding in relation to the observed facet sites has not been addressed, and this is the crucial aspect of this study and is addressed next.

4.1 Assessment of soft-hard grain stress state, slip, and load shedding at disc facet sites

The oligocrystal sub-models subjected to the appropriate stress state dwell loading at each of the disc facet sites are utilised in order to extract out the peak values of hard-grain basal stress at the end of the dwell period (as demonstrated and developed in Figure 9). During the dwell, load shedding potentially occurs (depending on the stress state and hence the soft grain prism system resolved shear stress). Peak hard-grain basal stresses are plotted against the associated normalised shear stress for the site-specific load history (which Figure 14 shows to be equivalent to stress state) as shown in Figure 15(a) for loading conditions in which dwell is both included and excluded.

Figure 15 Peak hard-grain stresses versus normalised re:

(b) predicted facet site

;d facet site f/

hows aga

esolved shear str : stress state, as ir

stress (a) from sub-model crystal plasticity study, and interpolated from (a)

Firstly, the graph in Figure 15 (a) shows again the role of the biaxiality stress state, or equivalently, the soft-grain prism system resolved shear stress on the resulting hard-grain basal stresses. More significantly, it makes very clear that the inclusion of the dwell leads to very strong stress redistribution and load shedding, pushing up the hard-grain basal stresses greatly. The magnitude of the load shedding, interestingly, is also seen to be dependent on the stress state (and hence equivalently, the resolved shear stress). The increase in basal stress in the hard grain, due to load shedding, is observed to be ~350MPa after dwell for a uniaxial stress state, but less than 50MPa for the equibiaxial loading case.

The results in Figure 15 (a) allow us to examine the load shedding and the hard-grain basal stresses at the disc sites at which facets are known to nucleate. With use of the map in Figure 14, the prism system resolved shear stresses at each disc facet site are extracted and combined with the results of Figure 15(a) in order to establish the peak basal stresses after load shedding at each facet site. These are shown in Figure 15(b) against the disc facet site resolved shear stress. This is reproduced in tabular format in Table 2.

Predicted Predicted

Facet Trss / T exPt a Facet Trss / Tc exPt. a

number number

(MPa) (MPa)

a1 a2 1.180 1 185 1027 1039

a3 a4 1.213 1 227 1105 1135

a5 a6 1.131 1 286 963 1287

a7 a8 1.274 1 199 1254 1076

a9 aio 1.221 1.180 1122 1039

ail ai2 1.223 1.241 1126 1168

a13 1.058 896

a14 1.231 1144

a15 1.209 1097

a16 1.272 1249

a17 1.286 1287

a18 1.174 1029

a19 1.272 1249

a20 1.286 1287

a21 1.174 1029

a22 1.091 920

a23 1.205 1088

a24 1.245 1177

a25 1.219 1118

a26 1.203 1084

Table 2 Predicted disc hard-grain basal stresses for disc facet sites from CPFE model

We note that the majority of facet sites cluster around a peak basal stress of 1100-1150MPa with some outliers on either side, both higher and lower, and are close to the measured c-axis tensile strength of titanium (~1200MPa) [32].

A clear message from these results is that the disc facet nucleation sites are in locations for which high resolved shear stresses exist on the soft grain of the hard-soft combination, for which considerable load shedding occurs, significantly pushing up the hard-grain basal stresses. In summary, all the disc facet locations generate high load shedding and high basal stresses. The primary significance of this is that load shedding is an essential factor in identifying those disc locations where faceting is observed (where also, the hard-soft grain combination also has to exist). This provides persuasive evidence that facet nucleation in discs is caused by load shedding; what is clear is that were the load shedding not to occur, the stresses at the facet sites would remain unremarkable.

The results also provide persuasive evidence that the rogue, or hard-soft grain combination, with respect to the primary (hoop) loading, is also an essential feature of facet formation in discs. Finally, it is also clear that the stress state local to the hard-soft grain combination is also a key factor in facet nucleation, since it determines the magnitude of the resolved shear stress on the soft grain prism system, thus determining the slip activation and hence the load shedding during dwell. A uniaxial stress state at the microstructure level transpires as a result to be more damaging than a biaxial (or more generally triaxial) stress state.

We have presented a multi-scale linkage of the details of microstructure using crystal plasticity (which is the scale that must be considered for cold dwell) right up to the component in-service stress scale relevant to aero-engines in operation. A methodology has been established which potentially provides a strong criterion for dwell failure. The macroscale operating stresses of the discs are far

below macroscopic yield stresses, and discs may be adequately and appropriately modelled using elastic analyses. With knowledge of the hard-soft grain combination, resolved shear stresses for any location in the disc geometry may be extracted (assuming a rogue grain combination is present) from which we can ultimately infer likely hard grain basal stresses, which provides the criterion for failure.

5 Conclusions

This paper establishes the link between disc component level testing and crystal plasticity finite element modelling, in the context of cold dwell fatigue. Furthermore this paper has presented a detailed examination of the contribution of load shedding to facet nucleation in dwell fatigue in compressor disc geometries through investigation of the local stress state at the rogue grain combination and how it relates to the macroscopically applied boundary conditions.

The hard grain orientation with respect to the stress state is a key requirement for faceting; misalignment of up to 15° with the principal tensile loading axis shows equivalent basal stresses whereas with further deviation from the principal tensile axis, the stress on the basal plane drops off significantly. This explains why most facets observed in dwell fatigue experiments have been found to be orientated to the principal loading direction within +/-15o.

The role of the biaxial remote stress state on the rogue pair has been investigated, and results have been found to be commensurate with the findings of Doquet and DeGreef [8], who hypothesise that biaxiality may have an ameliorating effect on fatigue life (with the caveat of a limited experimental dataset). The detailed micro-level study presented here provides the mechanistic basis for this phenomenon.

The disc component facet observations together with the crystal plasticity sub-model oligocrystal approach provide persuasive evidence that a hard-soft grain combination is required for facet formation, that the remote stress state influences the resolved shear stress on the soft grain initiating slip (with uniaxial stress state more damaging than biaxial), and that the load shedding which results is essential in pushing up the hard-grain basal stress to nucleate facets. The absence of any one of the rogue grain combination, the stress state, or the load shedding renders the observed disc facet locations unexplainable. However, together, they provide a persuasive mechanistic basis and quantitatively predictive methodology for facet nucleation in aero-engine discs.


The authors wish to acknowledge input and perspectives from Prof David Rugg, Dr Adrian Walker, and Dr Kate Fox at Rolls-Royce plc. We wish to thank Dr Zhen Zhang, Imperial College, for helpf discussions.


MC and FPED gratefully acknowledge financial support from Rolls-Royce plc and from the Engineering and Physical Sciences Research Council funded HexMat programme grant EP/K034332/1 ( FPED acknowledges Royal Academy of Engineering/Rolls-Royce research chair funding.


> Roya ;r, 2007.

[1] G. Lutjering and J. J. C. Williams, Titanium : Springer,

[2] D. Rugg, M. Dixon, and F. P. Dunne, 'Effective structural unit size in titanium alloys', J. Strain Anal. Eng. Des., vol. 42, no. 4, pp. 269-279, 2007.

[3] A. N. Stroh, 'The Formation of Cracks as a Result of Plastic Flow', Proc. R. Soc. Lond. Ser. Math. Phys. Sci., vol. 223, no. 1154, pp. 404-414, May 1954.

[4] W. J. Evans and M. R. Bache, 'Dwell-sensitive fatigue under biaxial loads in the near-alpha titanium alloy IMI685', Int. J. Fatigue, vol. 16, no. 7, pp. 443-452, 1994.

[5] V. Hasija, S. Ghosh, M. J. Mills, and D. S. Joseph, 'Deformation and creep modeling in polycrystalline Ti-6Al alloys', Acta Mater., vol. 51, no. 15, pp. 4533-4549, Sep. 2003.

[6] F. P. E. Dunne and D. Rugg, 'On the mechanisms of fatigue facet nucleation in titanium alloys', Fatigue Fract. Eng. Mater. Struct., vol. 31, no. 11, pp. 949-958, 2008.

M. Anahid, M. K. Samal, and S. Ghosh, 'Dwell fatigue crack nucleation model based on crystal plasticity finite element simulations of polycrystalline titanium alloys', J. Mech. Phys. Solids, vol. 59, no. 10, pp. 2157-2176, Oct. 2011.

[8] V. Doquet and V. De Greef, 'Dwell-fatigue of a titanium alloy at room temperature under uniaxial or biaxial tension', Int. J. Fatigue, vol. 38, pp. 118-129, May 2012.

[9] V. Bonnand, J. L. Chaboche, P. Gomez, P. Kanoute, and D. Pacou, 'Investigation of multiaxial fatigue in the context of turboengine disc applications', Int. J. Fatigue, vol. 33, no. 8, pp. 10061016, Aug. 2011.

[10] A. R. Kallmeyer, A. Krgo, and P. Kurath, 'Evaluation of Multiaxial Fatigue Life Prediction Methodologies for Ti-6Al-4V', J. Eng. Mater. Technol., vol. 124, no. 2, pp. 229-237, Mar. 2002.

[11] J. R. Rice and D. M. Tracey, 'On the ductile enlargement of voids in triaxial stress fields Mech. Phys. Solids, vol. 17, no. 3, pp. 201-217, Jun. 1969.

>*', J.

he dam

[12] A. L. Helbert, X. Feaugas, and M. Clavel, 'The influence of stress trlaxiality on the damage mechanisms in an equiaxeda/p Ti-6AI-4V alloy', Metall. Mater. Trans. A, vol. 27, no. 10, pp. 3043-3058, Oct. 1996.

[13] J. Kumar, B. Srivathsa, and V. Kumar, 'Stress triaxiality effect on fracture behavior of IMI-834 titanium alloy: A micromechanics approach', Mater. Des., vol. 30, no. 4, pp. 1118-1123, Apr. 2009.

[14] P. Lefranc, V. Doquet, M. Gerland, and C. Sarrazin-Baudoux, 'Nucleation of cracks from shear-induced cavities in an a/p titanium alloy in fatigue, room-temperature creep and dwell-fatigue', Acta Mater., vol. 56, no. 16, pp. 4450-4457, Sep. 2008.

[15] M. Gerland, P. Lefranc, V. Doquet, and C. Sarrazin-Baudoux, 'Deformation and damage mechanisms in an a/p 6242 Ti alloy in fatigue, dwell-fatigue and creep at room temperature. Influence of internal hydrogen', Mater. Sci. Eng. A, vol. 507, no. 1-2, pp. 132-143, May 2009.

[16] F. P. E. Dunne, A. Walker, and D. Rugg, 'A systematic study of hcp crystal orientation and morphology effects in polycrystal deformation and fatigue', Proc. R. Soc. Math. Phys. Eng. Sci., vol. 463, no. 2082, pp. 1467-1489, Jun. 2007.

[17] M. R. Bache, F. P. E. Dunne, and C. Madrigal, 'Experimental and crystal plasticity studies of deformation and crack nucleation in a titanium alloy', J. Strain Anal. Eng. Des., vol. 45, no. 5, pp. 391-399, Jul. 2010.

i] c L.

[18] A. L. Pilchak and J. C. Williams, 'Observations of Facet Formation in Near-a Titanium and iments on the Role of Hydrogen', Metall. Mater. Trans. A, vol. 42, no. 4, pp. 1000-1027, Oct. 2010.

[19] Z. Zhang, M. A. Cuddihy, and F. P. E. Dunne, 'On rate-dependent polycrystal deformation: the temperature sensitivity of cold dwell fatigue', Proc R Soc A, vol. 471, no. 2181, Sep. 2015.

[20] ASTM International, 'Standard Test Method for Determining Residual Stresses by the Hole-Drilling Strain Gage Method'. 2013.

[21] D. Dye, B. Roder, S. Tin, M. Rist, J. James, and M. Daymond, 'Modeling and measurement of residual stresses in a forged IN718 superalloy disc', in Superalloys 2004, 2004, pp. 315-322.

[22] F. P. E. Dunne, D. Rugg, and A. Walker, 'Lengthscale-dependent, elastically anisotropic, physically-based hcp crystal plasticity: Application to cold-dwell fatigue in Ti alloys', Int. J. Plast, vol. 23, no. 6, pp. 1061-1083, Jun. 2007.

[23] E. H. Lee, 'Elastic-Plastic Deformation at Finite Strains', J. Appl. Mech., vol. 36, no. 1, pp. 1-6, Mar. 1969.

[24] P. D. Littlewood and A. J. Wilkinson, 'Local deformation patterns in Ti-6Al-4V under tensile, fatigue and dwell fatigue loading', Int. J. Fatigue, vol. 43, pp. 111-119, Oct.


■ain depenc

[25] J. Jiang, T. B. Britton, and A. J. Wilkinson, 'The orientation and strain dependence of dislocation structure evolution in monotonically deformed polycrystalline copper', Int. J. Plast., vol. 69, pp. 102-117, Jun. 2015.

[26] J. Gong, T. Benjamin Britton, M. A. Cuddihy, F. P. E. Dunne, and A. J. Wilkinson, ' <a> Prismatic, <a> basal, and <c+a> slip strengths of commercially pure Zr by micro-cantilever tests', Acta Mater., vol. 96, pp. 249-257, Sep. 201

Dunne, an )f commerc '

tion in bi-

[27] H. Liang and F. P. E. Dunne, 'GND accumulation in bi-crystal deformation: Crystal plasticity analysis and comparison with experiments', Int. J. Mech. Sci., vol. 51, no. 4, pp. 326-333, Apr. 2009.

[28] J. Gong and A. J. Wilkinson, 'Anisotropy in the plastic flow properties of single-crystal a titanium determined from micro-cantilever beams', Acta Mater., vol. 57, no. 19, pp. 5693-5705, Nov. 2009.

llkinson, Ani

[29] H. Conrad, 'Thermally activated deformation of a titanium below 0.4 TM', Can. J. Phys., vol. 45, no. 2, pp. 581-590, Feb. 1967.

[30] M. Sauzay, J. Liu, F. Rachdi, L. Signor, T. Ghidossi, and P. Villechaise, 'Physically-based simulations of the cyclic behavior of FCC polycrystals', AdvMat Res, pp. 891-892, 2014.

I] V. Sir si

[31] V. Sinha, J. E. Spowart, M. J. Mills, and J. C. Williams, 'Observations on the faceted initiation site in the dwell-fatigue tested Ti-6242 alloy: Crystallographic orientation and size effects', Metall. Mater. Trans. A, vol. 37A, no. 5, pp. 1507-1518, 2006.

[32] I. P. Jones and W. B. Hutchinson, 'Stress-state dependence of slip in Titanium-6Al-4V and other H.C.P. metals', Acta Metall., vol. 29, no. 6, pp. 951-968, Jun. 1981.


Combining industrial scale testing of aero engine components with elastic finite element analysis and crystal plasticity sub models

Investigations into the multi axial stress state dependence on facet nucleation in cold dwel fatigue

Provides insight into the mechanism of facet nucleation through load shedding in tit alloys