Low cycle fatigue behavior of single crystal superalloy with temperature gradient

In present study, low cycle fatigue model based on rate dependent constitutive formulation with kinematic hardening and fatigue damage rule is introduced to investigate the fatigue behavior of single crystal superalloys with temperature gradient. Low cycle fatigue tests with uniform temperature and with temperature gradient are carried out to investigate the influence of temperature gradient. The microstructure reveals that the slip deformation is a principal mechanism of low cycle fatigue with temperature gradient. Complex stress experiments of low cycle fatigue are carried out to verify the life prediction rule. Bauschinger effect, ratcheting effect and fatigue damage of single crystal superalloy at different temperatures are studied. Special attention is placed on the simulation of low cycle fatigue behavior with temperature gradient. The simulation results show that the ratcheting effect at high temperature is more remarkable than that at low temperature due to the influence of temperature gradient, which results in the difference of damage evolution between high temperature zone and low temperature zone. Comparison between experiments and simulations with temperature gradient proves that the predicted model based on the damage of the second cycle is reasonable.     

高温合金材料循环相关和时间相关热机械疲劳循环性能模拟

从平衡热力学不可逆系统出发，用非线性粘弹塑性运动强化莱模拟高温合金材料的应变控制热机械疲劳循环特性。讨论了温度变化和应变循环的相位关系，循环相关和时间相关热机械疲劳损伤机制，蠕变和疲劳间的相互作用。在建立本构关系和状态方程时，均考虑了温度变化所产生的影响。     

A model for high temperature creep of single crystal superalloys based on nonlocal damage and viscoplastic material behavior

A model for high temperature creep of single crystal superalloys is developed, which includes constitutive laws for nonlocal damage and viscoplasticity. It is based on a variational formulation, employing potentials for free energy, and dissipation originating from plasticity and damage. Evolution equations for plastic strain and damage variables are derived from the well-established minimum principle for the dissipation potential. The model is capable of describing the different stages of creep in a unified way. Plastic deformation in superalloys incorporates the evolution of dislocation densities of the different phases present. It results in a time dependence of the creep rate in primary and secondary creep. Tertiary creep is taken into account by introducing local and nonlocal damage. Herein, the nonlocal one is included in order to model strain localization as well as to remove mesh dependence of finite element calculations. Numerical results and comparisons with experimental data of the single crystal superalloy LEK94 are shown.     

Lifetime simulation of thermo-mechanically loaded components

The estimation of the lifetime of thermo-mechanically loaded components by testing is very costly and time-consuming, since the high temperature cycle time in practical application dominates the test duration. Common frequencies for TMF (thermo-mechanical fatigue) tests are at about 0.01 Hz compared to 10–100 Hz at HCF (high cycle fatigue) and about 0.1–1 Hz at isothermal LCF (low cycle fatigue) tests. Therefore, the simulation of fatigue life is an important design step in the fast moving and competitive automotive industry, where the steady rise of engine power and the demand for lightweight construction concurrent with enhanced reliability require an optimised dimensioning process. Methods and models are usually derived from results made on tests with specimens, since it is possible to systematically and exactly define loading parameters and measurement categories. After an extensive test programme (tensile tests, creep tests, low cycle fatigue tests and thermo-mechanical fatigue tests with different influences on specimens) it was possible to develop material models for the simulation of the time- and temperature dependent stress–strain hystereses and damage models for the simulation of the TMF lifetime. Based on this knowledge the whole simulation chain to determine the TMF life of a component is introduced: thermal calculation, mechanical calculation and lifetime calculation. Furthermore the transferability of specimen based simulation models to real components (an alternative test piece and a cylinder head) is investigated.     

Multiaxial creep and cyclic plasticity in nickel-base superalloy C263

Physically-based constitutive equations for uniaxial creep deformation in nickel alloy C263 [Acta Mater. 50 (2002) 2917] have been generalised for multiaxial stress states using conventional von Mises type assumptions. A range of biaxial creep tests have been carried out on nickel alloy C263 in order to investigate the stress state sensitivity of creep damage evolution. The sensitivity has been quantified in C263 and embodied within the creep constitutive equations for this material. The equations have been implemented into finite element code. The resulting computed creep behaviour for a range of stress state compares well with experimental results. Creep tests have been carried out on double notched bar specimens over a range of nominal stress. The effect of the notches is to introduce multiaxial stress states local to the notches which influences creep damage evolution. Finite element models of the double notch bar specimens have been developed and used to test the ability of the model to predict correctly, or otherwise, the creep rupture lifetimes of components in which multiaxial stress states exist. Reasonable comparisons with experimental results are achieved. The γ^′ solvus temperature of C263 is about 925 °C, so that thermo-mechanical fatigue (TMF) loading in which the temperature exceeds the solvus leads to the dissolution of the γ^′ precipitate, and a resulting solution treated material. The cyclic plasticity and creep behaviour of the solution treated material is quite different to that of the material with standard heat treatment. A time-independent cyclic plasticity model with kinematic and isotropic hardening has been developed for solution treated and standard heat treated nickel-base superalloy C263. It has been combined with the physically-based creep model to provide constitutive equations for TMF in C263 over the temperature range 20–950 °C, capable of predicting deformation and life in creep cavitation-dominated TMF failure.     

镍基单晶高温合金定向粗化行为及高温蠕变力学性能研究进展

镍基单晶高温合金是一种广泛应用于航空发动机和工业燃气轮机的两相叶片材料,由软的$\gamma $ 基体相和均匀镶嵌在其中的立方状 $\gamma'$ 沉淀强化相组成.它有个显著的特征,即在高温施加应力条件下, $\gamma '$沉淀相会发生定向粗化, 形成筏状.这种筏化行为直接影响了合金的蠕变疲劳寿命,是镍基单晶高温合金强化机制研究的重点. 此外,镍基单晶高温合金无晶界, 不存在高温晶界弱化、纵向晶界裂纹等问题.因此, $\gamma$/$\gamma'$相界面的位错运动、微观结构以及在载荷和温度作用下的演化决定了其蠕变力学性能.本文从镍基单晶高温合金的微观强化机制出发对定向粗化行为及蠕变力学性能进行了综述.重点介绍了定向粗化行为发生的微观机理、驱动力、影响因素和蠕变过程中界面微结构演化、蠕变力学模型以及定向粗化对高温蠕变力学性能的影响,指出了高温蠕变力学性能研究的发展方向和仍待解决的问题.      

金属材料多轴非比例低周疲劳寿命预测概述

工程中的大多数构件承受着比例或非比例多轴疲劳荷载作用,而非比例强化效应会大大影响其多轴疲劳寿命。精确预测金属材料在多轴非比例荷载下的低周疲劳寿命需要同时考虑等向强化、随动强化及非比例强化效应下的材料本构关系,并在临界面上计算出相应应力应变值,根据不同疲劳失效形式选取不同类型的失效模型来确定疲劳寿命．本文针对这一过程中重要知识点进行阐述,并介绍了相关模型与理论．     

On the evolution of lattice deformation in austenitic stainless steels—The role of work hardening at finite strains

In this work, a three dimensional crystal plasticity-based finite element model is presented to examine the micromechanical behaviour of austenitic stainless steels. The model accounts for realistic polycrystal micromorphology, the kinematics of crystallographic slip, lattice rotation, slip interaction (latent hardening) and geometric distortion at finite deformation. We utilise the model to predict the microscopic lattice strain evolution of austenitic stainless steels during uniaxial tension at ambient temperature with validation through in situ neutron diffraction measurements. Overall, the predicted lattice strains are in very good agreement with those measured in both longitudinal and transverse directions (parallel and perpendicular to the tensile loading axis, respectively). The information provided by the model suggests that the observed nonlinear response in the transverse {200} grain family is associated with a competitive bimodal evolution of strain during inelastic deformation. The results associated with latent hardening effects at the microscale also indicate that in situ neutron diffraction measurements in conjunction with macroscopic uniaxial tensile data may be used to calibrate crystal plasticity models for the prediction of the inelastic material deformation response.     

Constitutive and damage model for 63Sn37Pb solder under uniaxial and torsional cyclic loading

A constitutive model with Ohno–Wang kinematic hardening rule is developed and employed to simulate the isothermal cyclic behavior of Sn–Pb solder under uniaxial and torsional loading. An implicit constitutive integration scheme is presented for inelastic flow of solder. Then a modified low cycle fatigue life prediction model is put forward in which the sum of maximum shear strain range and normal strain range based on the critical plane concept is adopted to replace the uniaxial strain range used by Stolkarts et al. [Stolkarts, V., Keer, L.M., Fine, M.E., 1999. Damage evolution governed by microcrack nucleation with application to the fatigue of 63Sn–37Pb solder. J. Mech. Phys. Solids 47, 2451–2468]. Comparison of the experimental results and simulation verifies that the stress strain hysteresis loops and peak stress decline curve of solder can be reasonably modeled over a wide range of loading conditions with implement of damage coupled constitutive model, and the lifetime estimations of 63Sn37Pb solder based on the assumption of microcrack nucleation governed damage is effective to provide a conservative prediction.     

Strain gradient crystal plasticity analysis of a single crystal containing a cylindrical void

The effects of void size and hardening in a hexagonal close-packed single crystal containing a cylindrical void loaded by a far-field equibiaxial tensile stress under plane strain conditions are studied. The crystal has three in-plane slip systems oriented at the angle 60° with respect to one another. Finite element simulations are performed using a strain gradient crystal plasticity formulation with an intrinsic length scale parameter in a non-local strain gradient constitutive framework. For a vanishing length scale parameter the non-local formulation reduces to a local crystal plasticity formulation. The stress and deformation fields obtained with a local non-hardening constitutive formulation are compared to those obtained from a local hardening formulation and to those from a non-local formulation. Compared to the case of the non-hardening local constitutive formulation, it is shown that a local theory with hardening has only minor effects on the deformation field around the void, whereas a significant difference is obtained with the non-local constitutive relation. Finally, it is shown that the applied stress state required to activate plastic deformation at the void is up to three times higher for smaller void sizes than for larger void sizes in the non-local material.     

Dwell fatigue crack nucleation model based on crystal plasticity finite element simulations of polycrystalline titanium alloys

In this paper a crystal plasticity-based crack nucleation model is developed for polycrystalline microstructures undergoing cyclic dwell loading. The fatigue crack nucleation model is developed for dual-phase titanium alloys admitting room temperature creep phenomenon. It is a non-local model that accounts for the cumulative effect of slip on multiple slip systems, and involves evolving mixed-mode stresses in the grain along with dislocation pileups in contiguous grains. Rate dependent, highly anisotropic behavior causes significant localized stress concentration that increases with loading cycles. The crystal plasticity finite element (CPFE) model uses rate and size-dependent anisotropic elasto-crystal plasticity constitutive model to account for these effects. Stress rise in the hard grain is a consequence of time-dependent load shedding in adjacent soft grains, and is the main cause of crack nucleation in the polycrystalline titanium microstructure. CPFE simulation results are post-processed to provide inputs to the crack nucleation model. The nucleation model is calibrated and satisfactorily validated using data available from acoustic microscopy experiments for monitoring crack evolution in dwell fatigue experiments.     

Strain gradient plasticity by internal-variable approach with normality structure

This paper presents a constitutive formulation for materials with strain gradient effects by internal-variable approach with normality structure. Specific micro-structural rearrangements are assumed to account for the inelasticity deformations for this class of materials, and enter the constitutive formulations in form of internal variables. It is further assumed that the kinetic evolution of any specific micro-structural rearrangement may be fully determined by the thermodynamic forces associated with that micro-structural rearrangement, by normality relations via a flow potential. Macroscopic gradient-enhanced inelastic behaviours may then be predicted in terms of the microscopic internal variables and their conjugate forces, and thus a micro–macro bridging formulation is available for strain-gradient-characterised materials. The obtained formulations are first applied to crystallographic materials, and a crystal gradient plasticity model is developed to account for the influence of microscopic slip rearrangements on the macroscopic gradient-dependent mechanical behaviour for this class of materials. Micro-cracked geomaterials are also treated with these formulations and a gradient-enhanced damage constitutive model is developed to address the impacts of the evolutions of micro-cracks on the macroscopic inelastic deformations with strain gradient effects for these materials. The available formulations are further compared with other thermodynamic approaches of constitutive developing.     

考虑塑性应变记忆恢复的循环棘轮本构模型研究

在统一粘塑性循环本构理论框架下,以Ohno-Abdel-Karim非线性随动硬化模型为基础,建立了一个循环本构模型。模型通过引入塑性应变幅值记忆效应,并在塑性应变记忆项中加入恢复系数,提高了对循环硬化材料单轴棘轮行为的预言能力。将模型应用于316L不锈钢单轴棘轮行为的描述中,模拟不同平均应力、应力幅值下的棘轮应变,均与实验数据吻合较好,证明本文改进的本构模型能合理地描述循环硬化材料的单轴棘轮行为。     

Constitutive modelling of anisotropic creep deformation in single crystal blade alloys SRR99 and CMSX-4

A damage mechanics based model has been developed to model stress rupture and creep behaviour of the first and second generation single crystal superalloys SRR99 and CMSX-4. In this article the creep behaviour of CMSX-4 in several different orientations at 950°C is simulated using finite elements, these simulations are compared with the results of creep tests. In order that the effects of rotation and specimen bending can be accounted for in the analysis the entire creep specimen is modelled. The FE program ABAQUS has been used and the slip system model is written using a User MATerial subroutine (UMAT). EBSD (electron back scattered diffraction) measurements of the lattice rotations occurring during creep indicate that the active slip systems at 950°C are the <101>{111} and <112>{111} systems, our simulations show that the creep results can be explained by activating these two families of slip system. There is strong microstructural evidence that the significant components of the hardening matrix should be those causing self and latent hardening of the <101>{111} systems and latent hardening by the <101>{111} systems on the <112>{111} systems.     

An elasto-plastic damage constitutive theory and its prediction of evolution of subsequent yield surfaces and elastic constants

Based on pair functional potentials, Cauchy-Born rule and slip mechanism, a material model assembling with spring-bundle components, a cubage component and slip components is established to describe the elasto-plastic damage constitutive relation under finite deformation. The expansion/shrink, translation and distortion of yield surfaces can be calculated based on the hardening rule and Bauschinger effect defined on the slip component level. Both kinematic and isotropic hardening are included. Numerical simulations and predictions under tension, torsion, and combined tension-torsion proportional/non-proportional loading are performed to obtain the evolution of subsequent yield surfaces and elastic constants and compare with two sets of experimental data in literature, one for a very low work hardening aluminum alloy Al 6061-T6511, and another for a very high work hardening aluminum alloy annealed 1100 Al. The feature of the yield surface in shape change, which presents a sharp front accompanied by a blunt rear under proportional loading, is described by the latent hardening and Bauschinger effect of slip components. Further, the evolution law of subsequent yield surfaces under different proportional loading paths is investigated in terms of their equivalence. The numerical simulations under non-proportional loading conditions for annealed 1100 Al are performed, and the subsequent yield surfaces exhibit mixed cross effect because the kinematic hardening and isotropic hardening follow different evolution tendency when loading path changes. The results of non-proportional loading demonstrate that the present model has the ability to address the issue of complex loading due to the introduction of state variables on slip components. Moreover, as an elasto-plastic damage constitutive model, the present model can also reflect the variation of elastic constants through damage defined on the spring-bundle components.     

Measurement of Strain Localization Resulting from Monotonic and Cyclic Loading at 650 <Superscript>∘</Superscript>C in Nickel Base Superalloys

A methodology is presented for the use of the oxide scale that develops in polycrystalline Ni-base superalloys at service temperature, as a speckle pattern for μm-scale resolution strain measurements. Quantitative assessment of the heterogeneous strain field at the grain scale is performed by high-resolution SEM digital image correlation under monotonic and cyclic loading in polycrystalline Ni-base superalloys up to 650 ^°C. In the René 88DT superalloy, strain localization is observed near twin boundaries during low cycle fatigue (LCF) at intermediate temperatures, correlating with activation of {111} 〈110〉 and {111} 〈112〉 slip systems. A strong correlation between the microstructural configuration that promotes strain localization during monotonic loading and crack initiation at 650 ^°C in low cycle fatigue was observed.     

粉末冶金涡轮盘简单模拟件低周疲劳寿命研究

采用光滑圆棒试样和带孔平板试样,对不同温度下的镍基粉末高温合金(FGH95)的低周疲劳(LCF)寿命进行了试验研究和有限元分析。在详细分析试验和有限元计算结果的基础上,提出了复杂应力状态下的低周疲劳寿命模型。模型寿命表达为真实应力幅的函数,模型参数由不同应力水平加载作用下的光滑圆棒试样试验结果给定,进一步采用涡轮盘简单模拟件即带孔平板试样对比验证LCF寿命模型的有效性。有限元计算结果显示,理论预测寿命与试验结果能很好地吻合。     

沉淀相的尺寸与形态对镍基合金蠕变行为的影响

镍基合金具有优良的高温力学性能,广泛应用于涡轮叶片等热端部件。沉淀相的尺寸和形态是影响镍基合金力学性能的重要因素。本文在考虑应变梯度的镍基合金晶体塑性本构模型的基础上,引入了各向异性损伤张量,研究了包含两种不同尺寸和三种不同长细比的沉淀相形态的镍基合金蠕变行为。结果表明,该模型能够很好地反映沉淀相的尺寸对镍基合金蠕变行为的影响,与实验结果符合较好。同时,沉淀相的形态也对镍基合金的力学性能产生重要影响,随着沉淀相长细比的增加,镍基合金的蠕变寿命延长,这体现了粗化和形态对镍基合金蠕变行为影响的一种竞争的机制。     

基于微结构动态演化机制的单晶镍基高温合金晶体塑性本构及其有限元模拟

因其优异的高温力学性能,镍基单晶高温合金在航空航天和能源等领域得到了广泛的应用.镍基单晶高温合金优异的高温性能来源于其特有的两相微结构.基于代表体胞模型及分块均匀化方法,以位错密度为主要内变量,发展了一个包含两相微结构和位错演化信息的单晶镍基高温合金塑性行为的本构模型.该本构模型充分考虑了镍基单晶合金中位错在基体相和沉淀增强相中的多种演化机制,例如,基体位错八面体滑移、立方滑移、位错攀移、交滑移、位错弓出、位错切过沉淀增强相以及位错Kear-Wilsdolf(K-W)锁形成与解锁等.在商用有限元软件ABAQUS的框架下,编制了UMAT用户材料子程序.利用该用户子程序,对单晶和多晶镍基高温合金在不同温度、不同加载方向下的单调塑性、循环塑性、蠕变等典型行为进行了计算模拟.结果表明:该晶体塑性本构模型能"统一地"刻画镍基高温合金在不同温度、不同方向下的多种变形行为,并与实验结果具有良好的一致性.