不同服役条件下FGH96合金的蠕变特征

图1 蠕变试样及其解剖分析图

Fig.1 Creep sample and analysis method (RD—rolling direction, ND—normal direction, TD—transverse direction)

表1 FGH96合金蠕变实验方案

Table 1 The schedule of creep test

T / °C	σ / MPa	ε_c
704	690	3% and fracture
	720
	750
	780
	810
650	690	1% and fracture
670
720
750

Note:T—temperature, σ—applied stress, ε_c—residual strain

新窗口打开| 下载CSV

采用DM8000光学显微镜(OM)对蠕变前的FGH96合金的显微组织进行观察,采用S-4800 SEM对蠕变前合金的γ'相组织进行表征,采用Quanta 650-Nordlys Nano EBSD对合金中的晶界特征、局域取向差分布等信息进行分析,采用Tecnai G² F20 TEM对蠕变后FGH96合金内部的位错形态进行研究。

2 实验结果

2.1 显微组织

FGH96合金的金相组织和γ'相组织如图2所示。由图可知,经过固溶热处理的FGH96合金中的一次γ'相全部回溶,平均晶粒尺寸约为20 μm。合金内部均匀分布着较大尺寸的二次γ'相和细小的三次γ'相。采用物理化学相分析^[25]的方法分析了γ'相的尺寸分布,如图2c所示。本实验中FGH96合金的γ'相尺寸呈现单峰分布状态,其平均尺寸约为93.5 nm。

图2

图2 1160℃固溶处理的FGH96合金试样的显微组织及γ'相尺寸分布

Fig.2 OM (a) and SEM (b) images of the microstructure of FGH96 superalloy with 1160^oC solution, and the size distribution of γ' precipitates (c)

2.2 蠕变实验结果

不同条件下FGH96合金蠕变过程的应力-应变曲线如图3所示。一般来说,金属材料的蠕变过程可以分为3个阶段,即蠕变起始阶段、稳态蠕变阶段以及快速蠕变阶段。对于粉末高温合金,蠕变起始阶段所经历的时间较短,几乎无法在蠕变应力-应变曲线中体现出来；随着蠕变变形的进行,FGH96合金快速进入稳态蠕变阶段,在该阶段,合金的蠕变速率保持稳定；当蠕变变形达到一定程度后,蠕变速率开始加速,此时蠕变进入快速蠕变阶段。

图3

图3 FGH96合金蠕变曲线

Fig.3 The creep curves of FGH96 superalloy tested at 704^oC, 690-810 MPa (a) and 690 MPa, 650-750^oC (b)

图3a为704℃时不同加载应力条件下FGH96合金的蠕变曲线。可知,在蠕变温度为704℃时,随着加载应力的提高,合金达到1.0%残余应变的蠕变寿命逐渐降低。当加载应力为690 MPa时,达到1.0%残余应变所需要的蠕变时间为158.5 h,而当加载应力提高至810 MPa时,达到1.0%残余应变的蠕变寿命缩减到34.8 h。

图3b是加载应力为690 MPa时,不同温度条件下FGH96合金的蠕变曲线。可知,随着蠕变温度的提高,合金的蠕变寿命显著降低。当蠕变温度为650℃时,蠕变时间达到1400 h以上时,蠕变残余应变仅为0.04%,说明在650℃、690 MPa条件下,FGH96合金几乎不发生蠕变变形；当蠕变温度提高至750℃,达到1.0%的蠕变时间仅仅8.17 h。

综上分析,当涡轮盘的服役温度在700℃时,FGH96合金的服役应力应不超过780 MPa,否则将会引起严重的蠕变变形,从而影响涡轮盘的正常服役；当涡轮盘的服役应力在690 MPa时,应避免其服役温度超过720℃,否则将会出现涡轮盘蠕变失效风险。

根据蠕变应力-应变曲线,计算得到了不同测试条件下FGH96合金的稳态蠕变速率( $\dot{γ}$ ),如图4和表2所示。如图4a所示,当蠕变温度为704℃时,FGH96合金的 $\dot{γ}$ 随着加载应力的增大而增大,当加载应力由690 MPa提高至810 MPa时, $\dot{γ}$ 由3.76 × 10^-9 s^-1提高至4.98 × 10^-8,且加载应力每增加30 MPa,蠕变速率增长1倍。图4b是蠕变加载应力为690 MPa时,FGH96合金的 $\dot{γ}$ 随蠕变温度的变化曲线。由图可知,随着蠕变温度由650℃提高至750℃,合金的 $\dot{γ}$ 由4.19 × 10^-11 s^-1提高至1.24 × 10^-7 s^-1。蠕变温度每提高30℃, $\dot{γ}$ 提高一个数量级。

图4

图4 FGH96合金稳态蠕变速率随温度和应力的变化关系

Fig.4 Relationship between steady creep strain rate and loading stress (a) and temperature (b)

表2 不同条件下FGH96合金稳态蠕变速率

Table 2 Steady creep strain rates ( $\dot{γ}$ ) at different creep conditions

T / ^oC	$σ$ / MPa	$\dot{γ}$ / s^-1	T / ^oC	$σ$ / MPa	$\dot{γ}$ / s^-1
704	690	3.76 × 10^-9	650	690	4.19 × 10^-11
	720	7.82 × 10^-9	670		6.72 × 10^-10
	750	1.76 × 10^-8	704		3.76 × 10^-9
	780	2.74 × 10^-8	720		3.01 × 10^-8
	810	4.98 × 10^-8	750		1.24 × 10^-7

新窗口打开| 下载CSV

对于FGH96合金涡轮盘,其轮缘服役温度一般处于650~750℃之间,在此温度区间内,FGH96合金的稳态蠕变速率对服役温度更为敏感。服役温度的提高将会导致FGH96合金的蠕变速率呈指数增长。尤其当服役温度超过720℃后,应严格控制合金的应力水平,否则可能造成整个发动机因涡轮盘蠕变而失效。同时,当涡轮盘轮缘服役温度达到700℃后,其服役应力应控制在780 MPa以下。

3 分析讨论

3.1 蠕变速率分析

Peng等^[20]基于Gibbs位错运动理论和晶体塑性理论,建立了粉末高温合金的蠕变速率物理模型,该模型可以准确计算粉末高温合金的在稳态蠕变阶段的蠕变速率,如式(1)所示：

\dot{γ} = ρ_{g} b^{2} ν e x p (- \frac{Δ F}{k_{b} T}) s i n h (\frac{(σ c o s θ c o s φ - τ_{c}) Δ V}{k_{b} T})

(1)

式中,ρ_g是滑移位错密度,b是Burgers矢量模,ν是位错震动频率,k_b是Boltzmann常数,T是温度,ΔF是Helmholtz自由能,ΔV是热激活体积,σ是外加载荷,θ是加载方向与滑移方向之间的夹角,φ是加载方向与滑移面法向之间的夹角,τ_c是位错运动的临界剪切应力。对于FGH96合金,当合金的显微组织确定后, 式(1)中仅有σ、ΔF和T是变量,其余的参数均是定值,ΔF则会随着加载应力和温度的变化而变化。因此,根据式(1)可以判断,随着T和σ的增大,蠕变速率均会增大。

位错运动的本质是热激活作用,位错的热激活作用主要取决于激活距离d和逃逸速率Γ (s^-1),而逃逸速率会随着温度的提高而增加。因此,当温度提高时,位错运动的逃逸速率增大,其克服析出相和溶质原子运动的速率也就相应增大,这就导致了合金中位错滑移速率的增大,最终导致了蠕变速率的增大。

3.2 应力对蠕变变形机理的影响

局域取向差(local misorientation)是指同一晶粒内某一数据点与相邻数据点之间的取向偏差。局域取向差一部分是由EBSD的角分辨率偏差导致的系统误差,另一部分则是由样品本身局部取向变化导致的,这种取向变化常常与塑性变形过程中形成的亚结构有关。Kamaya^[26~29]采用EBSD技术获取了变形金属内部的局域取向差值,发现局域取向差的大小与材料内部的局部变形程度、塑性残余变形以及累计疲劳损伤等具有一定的对应关系。

图5是实验温度为704℃时,不同加载应力条件下FGH96合金经过3.0%蠕变残余变形后局域取向差KAM (kernel average misorientation)图和反极图(IPF)分析结果。可知,不同的加载应力条件下FGH96合金的局域取向差并未出现明显差异,整个样品上的局域取向差分布较为均匀。同时,随着加载应力的增大,经过蠕变变形后合金内部的反极值未出现明显差异,且反极值均较小,说明在690~810 MPa范围内加载,并不会造成合金蠕变变形机制的改变。

图5

图5 704℃、不同加载应力条件下经过1.0%蠕变变形的FGH96合金局域取向差(KAM)和反极图(IPF)分析

Fig.5 Kernel average misorientation (KAM) maps and inverse pole figures (IPFs) of FGH96 alloy with 1.0% crept strain under various loading stresses at 704^oC

(a) 690 MPa (b) 720 MPa (c) 750 MPa (d) 780 MPa (e) 810 MPa

图6是实验温度为704℃时,不同加载应力条件下FGH96合金经过1.0%蠕变残余变形后的TEM分析结果。由图可知,在不同应力条件下经过蠕变变形的FGH96合金内部均出现了大量的层错,对层错进行HRTEM观察发现,这些层错均是在( $11 \bar{1}$ )面上沿着[ $2 \bar{1} 1$ ]方向形成的。研究^[21]表明,对于FGH96合金,当基体位错在滑移过程中遇到γ'相后,会分解成2个不全位错a/6[112]和a/3[ $11 \bar{2}$ ],其中a/6[112]不全位错在连续的(111)原子面上滑移就会产生层错,当层错的原子层数超过3层时,即形成了伪微孪晶,层错界面再经过原子重排,即形成了微孪晶。因此,FGH96合金在704℃、加载应力为690~810 MPa范围内,FGH96合金的蠕变变形机制仍然以位错滑移为主。

图6

图6 704℃、不同加载应力条件下经过1.0%蠕变变形的FGH96合金TEM像

Fig.6 TEM images of FGH96 alloy with 1.0% crept strain under various loading stresses at 704^oC

(a) 690 MPa (b) 720 MPa (c) 750 MPa (d) 780 MPa (e) 810 MPa

(f) HRTEM and selected area electron diffraction (SAED) pattern (inset) of the area in Fig.6e (Arrow [ $2 \bar{1} 1$ ] shows the dislocation slip direction of stacking fault, [ $\bar{1} 1 \bar{1}$ ] and [ $11 \bar{1}$ ] represent the normal directions of ( $\bar{1} 1 \bar{1}$ ) and ( $11 \bar{1}$ ) faces)

根据Schmid原理,位错在滑移方向的剪切力τ取决于加载应力σ和Schmid因子。当加载应力较小时,合金内部Schmid因子较大的晶粒内部会发生位错的滑移；随着加载应力的增大,尽管某些晶粒的Schmid因子较小,但受到加载应力增大的影响,该晶粒内部在滑移方向上的剪切力仍然达到了位错滑移所需的应力,其内部位错也会在滑移面上滑移。因此,随着加载应力的增大,合金内部会有更多的晶粒满足位错滑移的要求从而导致塑性变形的加速,也就是蠕变速率的增大。

当加载应力为690 MPa (图6a)和720 MPa (图6b)时,合金内部的微孪晶集中分布在少数晶粒内部,与该晶粒相邻的晶粒或孪晶粒内部未出现微孪晶,说明此时的应力水平较低,仅有少数晶粒内部的滑移系满足Schmid原理,发生了位错滑移；当加载应力提高至750 MPa (图6c)和780 MPa (图6d)时,在孪晶界两侧的晶粒内部均出现了微孪晶,说明随着应力水平的提高,相邻晶粒内部的滑移系均满足了Schmid原理；当加载应力进一步提高至810 MPa (图6e)时,图中显示的3个晶粒内部均形成了大量的微孪晶,且微孪晶密度非常高,说明此时合金中大量晶粒均满足了Schmid原理,位错滑移现象更为显著,导致了蠕变速率的增大。

3.3 温度对蠕变变形机理的影响

根据蠕变实验结果,FGH96合金的稳态蠕变速率对温度较为敏感,温度每提高30℃,合金的稳态蠕变速率提高一个数量级。图7是实验应力为690 MPa时,不同温度条件下FGH96合金经过1.0%蠕变残余变形后局域取向差的KAM和IPF图。可知,当蠕变实验温度在704℃以下时(图7a~c),FGH96合金的局域取向差并未出现明显差异,整个样品上的局域取向差分布较为均匀,说明此时合金的蠕变变形较为均匀；当实验温度提高至720℃时(图7d),FGH96合金的少数晶界附近的局域取向差较大,说明经过蠕变变形后,合金晶界位置出现了轻微的应变集中；当温度提高至750℃后(图7e),经过变形后的合金内部的局域取向差明显增大,说明此时合金内部具有较高的应变水平,且应变集中区沿晶界分布。对于金属材料,当温度超过等强温度后,合金内部的晶界逐渐变为弱化区,从而导致应变在晶界处集中。根据图7中KAM随蠕变温度的演化结果,可以在一定程度上说明FGH96合金的等强温度应该处于704~720℃范围内,因此当蠕变温度提高至720℃以上后,应变逐渐在晶界处累积。

图7

图7 690 MPa应力、不同温度条件下经过1.0%蠕变变形的FGH96合金KAM和IPF分析

Fig.7 KAM and IPF analyses of FGH96 alloy with 1.0% crept strain under the stress of 690 MPa and various temperatures

(a) 650℃ (b) 670℃ (c) 704℃ (d) 720℃ (e) 750℃

同时,随着实验温度的提高,经过蠕变变形后的合金内部的反极值未出现明显差异,且反极值均较小,说明当实验温度在750℃以下时,FGH96合金在蠕变过程中晶粒取向无明显变化,合金的蠕变变形机制仍然以位错滑移为主。

采用TEM对不同温度下蠕变后的FGH96合金中的位错进行分析,如图8所示。与不同应力状态下的蠕变类似,FGH96合金中形成了大量的微孪晶,这些微孪晶依然是由( $11 \bar{1}$ )面上的位错沿着[ $2 \bar{1} 1$ ]方向运动形成。与应力增大后的规律相同,随着温度的提高,经过蠕变变形的合金内部的微孪晶数量逐渐增多,且当温度提高至750℃后,同一晶粒内部不同滑移系同时开启滑移。FGH96合金在650~750℃条件下进行蠕变变形时,位错滑移并形成微孪晶是其主要变形机制,位错运动是一种热激活过程,随着温度的提高,位错运动的能力加强,合金中不同方向的滑移系同时开启滑移,导致出现微孪晶的交叉滑移现象,因此蠕变速率随着温度的提高而明显增大。

图8

图8 FGH96合金在690 MPa、不同温度条件下蠕变后形成的微孪晶

Fig.8 Microtwins in FGH96 superalloy after creep at 690 MPa and various temperatures

(a) 650^oC (b) 670^oC (c) 704^oC (d) 720^oC (e) 750^oC

(f) HRTEM and selected area electron diffraction (SAED) pattern of the area in Fig.8d (Arrow [ $2 \bar{1} 1$ ] shows the dislocation slip direction of stacking fault, [ $11 \bar{1}$ ] and [ $\bar{1} 1 \bar{1}$ ] represent the normal directions of ( $11 \bar{1}$ ) and ( $\bar{1} 1 \bar{1}$ ) faces)

3.4 温度和应力对蠕变断裂机理的影响

Tian等^[4]研究了FGH100L合金的蠕变断裂机理,发现该合金在不同温度和应力条件下的蠕变断裂特征均为穿沿晶混合特征。在低温低应力状态下,合金的断裂机理以沿晶断裂为主,穿晶断裂为辅；随着蠕变温度和应力水平的提高,由于较高的应变速率,导致合金的断裂机理以穿晶断裂为主,沿晶断裂为辅。

Peng等^[22]对不同γ'相状态的FGH96合金的蠕变断裂机理进行了研究,发现随着γ'相尺寸的细化,FGH96合金的蠕变断裂机理由穿晶-沿晶断裂混合特征转变为沿晶断裂特征。图9是不同温度和应力条件下,FGH96合金蠕变断口源区特征图。由图可知,当温度处于650~750℃范围、应力处于690~810 MPa范围内时,FGH96合金的蠕变断裂均呈现典型的沿晶断裂特征。FGH96合金在蠕变过程中,合金表面被氧化后,O原子会沿着晶界向内扩散,使晶界逐渐氧化,导致晶界弱化,最终变现为沿晶断裂特征。

图9

图9 不同温度和应力条件下FGH96合金的蠕变断裂源区形貌

Fig.9 Fractographs of primary stages of FGH96 superalloy tested at different temperatures and stresses

(a) 650^oC, 690 MPa (b) 750^oC, 690 MPa (c) 704^oC, 690 MPa (d) 704^oC, 810 MPa

4 结论

(1) 当蠕变温度为704℃时,FGH96合金的蠕变性能随着应力水平的提高而降低；当加载应力为690 MPa时,FGH96合金的蠕变性能随着温度提高而显著降低,且FGH96合金的稳态蠕变速率对服役温度更为敏感,服役温度每提高30℃,将会导致蠕变速率提高一个数量级。

(2) 当FGH96合金涡轮盘在700℃以下服役时,最大应力不超过780 MPa即可保证涡轮盘安全性；当FGH96合金涡轮盘的服役应力不超过690 MPa时,其服役温度可以提高至720℃,同时应尽量避免FGH96合金在720℃以上长时使用。

(3) 在涡轮盘主要的服役条件范围内,即温度处于650~750℃范围,应力处于690~810 MPa范围时,FGH96合金的蠕变变形均以位错滑移为主,且位错在滑移过程中,会在( $11 \bar{1}$ )原子面上形成大量的微孪晶,并导致蠕变变形的累积。位错运动一种热激活运动,随着温度和应力的提高,合金中更多的滑移系开动,从而导致蠕变速率的增大。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

Zheng

X L

. Mechanical Properties of Materials [M]. 2nd Ed., Xi'an: Northwestern Polytechnical University Press, 1999: 114

郑修麟

. 材料的力学性能 [M]. 第2版. 西安: 西北工业大学出版社, 1999: 114

[2]

J S

, Che

X L

, Song

Effect of temperature and stress on creep behavior of welding joint in stainless steel

[J]. Weld. Joining, 2003, (1): 17

屈金山, 车小莉, 宋兵.

温度和应力对不锈钢焊接接头蠕变行为的影响

[J]. 焊接, 2003, (1): 17

[3]

Azadi

, Azadi

Evaluation of high-temperature creep behavior in Inconel-713C nickel-based superalloy considering effects of stress levels

[J]. Mater. Sci. Eng., 2017, A689: 298

[4]

Tian

, Hao

Z B

, Ge

C C

, et al.

Effects of stress and temperature on creep behavior of a new third-generation powder metallurgy superalloy FGH100L

[J]. Mater. Sci. Eng., 2020, A776: 139007

[本文引用: 2]

[5]

Hou

Z P

, Zhang

, et al.

High temperature creep damage behavior of a novel Co-Cr-Mo-Ni alloy

[J]. J. Iron Steel Res., 2019, 31: 683

侯智鹏, 张姝, 张鹏等.

新型Cr-Co-Mo-Ni合金的高温蠕变损伤

[J]. 钢铁研究学报, 2019, 31: 683

[6]

Condat

, Decamps

Shearing of γ' precipitates by single a/2<110> matrix dislocations in a γ/γ' Ni-based superalloy

[J]. Scr. Metall., 1987, 21: 607

[7]

Décamps

, Morton

A J

, Condat

On the mechanism of shear of γ' precipitates by single (a/2) <110> dissociated matrix dislocations in Ni-based superalloys

[J]. Philos. Mag., 1991, 64A: 641

[8]

Décamps

, Raujol

, Coujou

, et al.

On the shearing mechanism of γ' precipitates by a single (a/6)<112> Shockley partial in Ni-based superalloys

[J]. Philos. Mag., 2004, 1A: 91

[9]

Knowles

D M

, Chen

Q Z

Superlattice stacking fault formation and twinning during creep in γ/γ' single crystal superalloy CMSX-4

[J]. Mater. Sci. Eng., 2003, A340: 88

[10]

Chen

Q Z

, Knowles

D M

Mechanism of <112> /3 slip initiation and anisotropy of γ' phase in CMSX-4 during creep at 750^oC and 750 MPa

[J]. Mater. Sci. Eng., 2003, A356: 352

[11]

Kolbe

The high temperature decrease of the critical resolved shear stress in nickel-base superalloys

[J]. Mater. Sci. Eng., 2001, A319-321: 383

[12]

Viswanathan

G B

, Sarosi

P M

, Whitis

D H

, et al.

Deformation mechanisms at intermediate creep temperatures in the Ni-base superalloy René 88 DT

[J]. Mater. Sci. Eng., 2005, A400-401: 489

[13]

Viswanathan

G B

, Karthikeyan

, Sarosi

P M

, et al.

Microtwinning during intermediate temperature creep of polycrystalline Ni-based superalloys: Mechanisms and modelling

[J]. Philos. Mag., 2006, 86A: 4823

[14]

Sarosi

P M

, Viswanathan

G B

, Mills

M J

Direct observation of an extended complex stacking fault in the γ' phase of a Ni-base superalloy

[J]. Scr. Mater., 2006, 55: 727

DOI URL

[15]

Unocic

R R

, Viswanathan

G B

, Sarosi

P M

, et al.

Mechanisms of creep deformation in polycrystalline Ni-base disk superalloys

[J]. Mater. Sci. Eng., 2008, A483-484: 25

[16]

Kovarik

, Unocic

R R

, Li

, et al.

Microtwinning and other shearing mechanisms at intermediate temperatures in Ni-based superalloys

[J]. Prog. Mater. Sci., 2009, 54: 839

[17]

Zou

J W

, Wang

W X

Development and application of P/M superalloy

[J]. J. Aeronaut. Mater., 2006, 26: 244

邹金文, 汪武祥.

粉末高温合金研究进展与应用

[J]. 航空材料学报, 2006, 26: 244

[18]

Zhang

G Q

, Zhang

Y W

, Zheng

, et al.

Research progress in powder metallurgy superalloys and manufacturing technologies for aero-engine application

[J]. Acta Metall. Sin., 2019, 55: 1133

张国庆, 张义文, 郑亮等.

航空发动机用粉末高温合金及制备技术研究进展

[J]. 金属学报, 2019, 55: 1133

[19]

Wang

X Q

, Luo

X J

, Zou

J W

Effects of HIP temperature on microstructure of FGH96 superalloy

[J]. J. Aeronaut. Mater., 2006, 26: 293

王旭青, 罗学军, 邹金文.

热等静压温度对FGH96粉末高温合金显微组织的影响

[J]. 航空材料学报, 2006, 26: 293

[20]

Peng

Z C

, Tian

G F

, Jiang

, et al.

Mechanistic behaviour and modelling of creep in powder metallurgy FGH96 nickel superalloy

[J]. Mater. Sci. Eng., 2016, A676: 441

[本文引用: 2]

[21]

Peng

Z C

, Zou

J W

, Wang

X Q

Microstructural characterization of dislocation movement during creep in powder metallurgy FGH96 superalloy

[J]. Mater. Today Commun., 2020, 25: 101361

[22]

Peng

Z C

, Zou

J W

, Yang

, et al.

Influence of γ' precipitate on deformation and fracture during creep in PM nickel-based superalloy

[J]. Prog. Nat. Sci., 2021, 31: 303

DOI URL [本文引用: 2]

[23]

M Z

, Pham

M S

, Peng

Z C

, et al.

Creep deformation mechanisms and CPFE modelling of a nickel-base superalloy

[J]. Mater. Sci. Eng., 2018, A718: 147

[24]

Feng

Y F

, Zhou

X M

, Zou

J W

, et al.

Effect of cooling rate during quenching on the microstructure and creep property of nickel-based superalloy FGH96

[J]. Int. J. Miner. Metall. Mater, 2019, 26: 493

[25]

Liu

Q B

, Lu

C F

, Yan

Physicochemical phase analysis of a directionally solidified Ni-based cast superalloy

[J]. Metall. Anal., 2006, 26(2): 9

刘庆斌, 卢翠芬, 燕平.

一种定向凝固镍基铸造高温合金的物理化学相分析

[J]. 冶金分析, 2006, 26(2): 9

[26]

Kamaya

Measurement of local plastic strain distribution of stainless steel by electron backscatter diffraction

[J]. Mater. Charact., 2009, 60: 125

[27]

Kamaya

A smoothing filter for misorientation mapping obtained by EBSD

[J]. Mater. Trans., 2010, 51: 1516

DOI URL

[28]

Kamaya

Assessment of local deformation using EBSD: Quantification of accuracy of measurement and definition of local gradient

[J]. Ultramicroscopy, 2011, 111: 1189

DOI URL

[29]

Kamaya

Characterization of microstructural damage due to low-cycle fatigue by EBSD observation

[J]. Mater. Charact., 2009, 60: 1454