书
书
书
第
53
卷
第
4
期
2018
年
8
月
西
南
交
通
大
学
学
报
JOURNAL OF SOUTHWEST JIAOTONG UNIVERSITY
Vol. 53 No. 4
Aug. 2018
收稿日期
:
20170114
基金项目
:国家自然科学基金资助项目(
51478071
)
作者简介
:渠昱(
1975
—),男,博士研究生,研究方向为桥梁及隧道工程,桥梁钢结构,电话:
15215047021
,
Email
:
5225158@ qq. com
通信作者
:曾勇(
1980
—),男,副教授,研究方向为桥梁结构理论与养护策略研究,
Email
:
zycquc@ 126. com
引文格式
:渠昱,顾安邦,曾勇,等
.
桥梁结构钢裂纹塑性区的研究及应用[
J
]
.
西南交通大学学报,
2018
,
53
(
4
):
720726.
QU Yu
,
GU Anbang
,
ZENG Yong
,
et al. Study on the crack plastic zone of bridge structure steel and its application
[
J
]
. Journal of
Southwest Jiaotong University
,
2018
,
53
(
4
):
720726.
文章编号
:
02582724
(
2018
)
04072007 DOI
:
10. 3969 / j. issn. 02582724. 2018. 04. 008
桥梁结构钢裂纹塑性区的研究及应用
渠
昱
1
,
2
,
顾安邦
1
,
2
,
曾
勇
1
,
2
,
杜柏松
1
,
2
(
1.
重庆交通大学山区桥梁与隧道工程国家重点实验室,重庆
400074
;
2.
重庆交通大学山区桥梁结构与材料
教育部工程研究中心,重庆
400074
)
摘
要
:为了研究裂纹塑性对裂纹扩展的影响,利用工程简化算法、应力函数法、扩展有限元法对桥梁钢裂尖塑
性区的尺寸和形状分别进行了计算;由于平面应力和平面应变情况下尾迹场循环塑性的特性不同,利用不连续
扩展有限元对两种情况下尾迹场的循环塑性和塑性累积进行了模拟分析,探讨了裂尖塑性区、循环塑性区的形
成和尾迹场产生压应力的机理
.
研究结果表明:裂尖塑性区尺寸与应力水平(名义应力与屈服极限的比值)的平
方成正比,当应力水平大于
0. 4
时,裂尖塑性区尺寸需要考虑应力水平的影响;裂尖塑性区的形状以蝶形向前伸
展,使裂纹尾迹场免受裂尖高应力场的拉伸作用,有利于裂纹闭合;裂尖塑性区存在材料的逆向流动,在循环塑
性区裂纹表面的塑性累积产生压应力效应有利于裂纹提前闭合,这种塑性诱发的裂纹提前闭合对研究变幅加
载、过载引起的裂纹扩展滞后有重要意义
.
关键词
:裂尖塑性区;循环塑性区;塑性累积;裂纹闭合;裂纹扩展滞后;不连续伽辽金扩展有限元法
中图分类号
:
O342
文献标志码
:
A
Study On the Crack Plastic Zone of
Bridge Structure Steel and Its Application
QU Yu
1
,
2
,
GU Anbang
1
,
2
,
ZENG Yong
1
,
2
,
DU Baisong
1
,
2
(
1. State Key Laboratory Breeding Base of Mountain Bridge and Tunnel Engineering
,
Chongqing Jiaotong University
,
Chongqing 400074
,
China
;
2. Mountain Bridge and Materials Engineering Research Center of Ministry of Education
,
Chongqing Jiaotong University
,
Chongqing 400074
,
China
)
Abstract
:
In order to study the influence of crack plasticity on crack growth
,
the shape and size of
crack tip plastic zone of bridge steel were calculated by using engineering simplification algorithm
,
the
stress function method and extended finite element method respectively. As a result of the cyclic
plasticity of wake field under plane stress and plane strain being different
,
cyclic plasticity and plastic
accumulation of crack tip and wakes were simulated and analyzed by using the discontinuous extended
finite element method. The formation of crack tip plastic zone
,
cyclic plastic zone and the mechanism
of compressive stress of wake field were studied. The results showed that the appearances and sizes of
plastic zone at the crack tip were proportional to the square of the stress level
(
the ratio of nominal
stress to the yield limit
)
. When the stress level was greater than 0. 4
,
the size of the crack tip plastic
zone had to be taken into account the effect of the stress level. The butterflyshaped crack tip plastic
zone symmetrically extended forward
,
shielding the tensile force of the crack tip high stress field on the
wake field
,
which was beneficial to crack closure. There exists reverse flow of material in plastic zone
西
南
交
通
大
学
学
报
第
53
卷
略高阶项的影响
.
由于裂尖塑性区和构件剩余部分
尺寸相比较小,状态仍然以弹性为主,线弹性分析
仍适用
.
所以,如何改进线弹性方法估计塑性区是
具有实际意义
.
图
1
平面应力
I
型裂纹裂尖塑性区
Fig. 1 Plane stress
crack tip plastic zone for mode I
1. 2
应力函数计算裂尖塑性区
Westergaard
应力函数
Z
(
z
)提供了一个严格的
从线弹性应力场计算裂尖塑性区方法
[
13
]
,可以避
免简化算法截断引起的误差
.
利用应力函数
Z
(
z
)
得到的应力场满足所有边界条件,但是弹塑性边界
并不靠近裂尖,这可以将坐标原点移动到裂尖,再
次利用平面
Ⅰ
型裂纹加载的
Irwin
解
[
12
]
.
如果
(
x
,
y
)和(
r
,
θ
)分别是直角坐标系和极坐标系,原
点设在裂尖上,
z = rexp θ = x + iy
是复函数,
i
2
=
- 1
,
Irwin
解为
σ
xx
σ
yy
τ
{ }
xy =
Re Z
(
z
)
- yIm Z′
(
z
)
Re Z
(
z
)
+ yIm Z′
(
z
)
- yRe Z′
(
z
{
)
}
,
(
6
)
Z
(
z
)
=
(
a + rcos θ + irsin θ
)
σ
n
(
a + rcos θ + irsin θ
)
2
- a
槡
2
,
Z′
(
z
)
=
- a
2
σ
n
[(
a + rcos θ + irsin θ
)
2
- a
2
]
3 / 2
.
(
7
)
首先利用
Z
(
z
)和
Z′
(
z
)函数得到应力,然后
利用等效
Mises
应力等于屈服强度
S
y
,就可以得到
平面应力情况下的塑性区方程,如式(
8
)
.
{[
Re Z -yIm Z′ -σ
n
]
2
+
[
Re Z +yIm Z′
]
2
-
[
Re Z - yIm Z′- σ
n
][
Re Z + yIm Z′
]
+
3
(
- yRe Z′
)
2
}
1 / 2
- S
y
= 0.
(
8
)
为了考虑
SIF
的影响,对于尖裂纹,椭圆裂纹
镜像半轴取裂纹开度
C
TOD
的一半,即
b = C
TOD
/ 2 =
2K
2I
/
(
πS
y
E
),计算的平面应力裂尖塑性区结果见
图
2.
尽管应力函数表达的完备应力场对平面裂纹
的线弹性解是正确的,但是,塑性区内的应力分布
与线弹性场的奇异性不匹配
.
塑性区内的能量不平
衡会改变塑性区周围的应力分布,因此塑性区的大
小和形状需要进行修正
.
图
2
由应力函数计算的
Ⅰ
型
裂纹裂尖塑性区结果
Fig. 2 Estimated results
from stress function for the mode Ⅰ
1. 3
满足平衡条件要求的塑性区修正
应力修正就是强制平面裂纹满足平衡条件
.
塑
性区内利用
Mises
应力对屈服产生的应力截断进
行平衡校正是合理的解决方案,因为它能将所有应
力都与
Mises
应力相关
.
本文只讨论(
σ
n
/ S
y
= 0. 8
)
Mises
应力对塑性区的修正
[
14
]
,按式(
9
)计算,结果
见图
3.
p
z eq M
=
∫
pz M
0
σ
M
(
r
,
θ
)
dr / S
y
,
(
9
)
式中:
p
z eq M
为
Mises
应力平衡修正后的裂尖塑性
区;
p
z M
为利用
Westergaard
函数计算的裂尖塑性
区;
σ
M
(
r
,
θ
)为修正点的
Mises
应力
.
图
3 Ⅰ
型裂纹裂尖塑性区平衡
校正前后的结果比较
Fig. 3 Comparison of crack tip plastic zone before
and after equilibrium correction for mode Ⅰ crack
2
2
7
第
4
期
渠
昱,等:桥梁结构钢裂纹塑性区的研究及应用
图
3
中:
A
表示
K
I
+ T
的裂尖塑性区;
B
表示
Mises
应力修正的
K
I
+ T
裂尖塑性区;
C
表示应力
函数
Z
计算的裂尖塑性区;
D
表示
Mises
应力修正
的应力函数
Z
计算的裂尖塑性区
.
从图
3
可以看出:高应力水平对塑性区的影响
需要修正;低应力水平对塑性区的影响很小,不需
要进行塑性区修正
.
2
数值方法
为了与理论值进行比较,利用扩展有限元
(
XFEM
)对桥梁钢裂尖塑性区进行数值计算
.
2. 1
利用
XFEM
计算
Q345qD
裂尖塑性区
XFEM
的基本想法是将不连续位移场分解为
两部分:连续部分是标准有限元插值;不连续部分
是根据局部单位分解定理
[
15
]
增加的部分,并将这
些附加信息嵌入有限元插值
.
这样的矩阵仍然是稀
疏的,仍可以用现有的有限元分析框架处理不连续
问题
.
裂纹采用两个水平集函数描述
[
15
]
:描述裂纹
路径的水平集函数和裂尖位置的水平集函数
.
附加
函数:完全被裂纹分割的单元节点附加采用有向距
离函数;裂尖单元节点附加取分支函数
[
15
]
.
为了消
除在附加和非附加单元之间存在混合(过渡)单元
产生多余项问题,利用内罚函数方法使附加和不附
加覆盖之间强制连续,称为不连续伽辽金扩展有限
元
[
16
]
,比标准的扩展有限元更精确,收敛速度得到
了优化
.
利用
Abaqus
对
Q345qD
桥梁结构钢板进行计
算,钢板 尺 寸 为 宽
1200 mm ×
长
1 600 mm ×
厚
4 mm
,中心裂纹长度
10 mm
,划分成
329 873
个
C3D8R
单元,
147 070
个节点,裂纹处网格加密尺
寸为
0. 5 mm × 0. 5 mm
,裂纹垂直方向单向加载,
应力比
σ
n
/ S
y
= 0. 8
,
0. 6
,
0. 4.
根据应力场得到的
裂尖塑性区结果见图
4.
图
4
不同应力比情况下中心裂纹塑性区
Fig. 4 Central crack plastic zone under different stress ratios
2. 2
计算结果比较
根据以上计算的结果,取裂纹长度
a = 10 mm
,
裂尖塑性区半径比较见表
1.
表
1
各种计算方法的塑性半径比较
Tab. 1 Comparison of plastic radii
of calculation methods mm
计算方法
修正前后
σn/ Sy
0. 8 0. 6 0. 4
KI+ T
前
3. 20 1. 80 0. 80
后
6. 40 3. 60 1. 60
Westergoard
函数
前
后
4. 03 2. 27 1. 01
8. 06 4. 54 2. 02
XFEM
—
4. 28 1. 72 0. 70
由表
1
可知:在低应力比情况下,几种算法基
本接近,在高应力比情况下,
Westergoard
函数与
XFEM
法计算结果很接近,说明数值方法是可行的
.
从图
1 ~ 4
可知,裂尖塑性区对称向前伸展
.
由
于塑性区不能传递拉应力,裂尖塑性区的两个“翅
膀”遮蔽了裂尖高应力对裂纹尾迹场的拉伸作用,
这有利于裂纹闭合
.
3
裂尖塑性区和尾迹场塑性区
裂纹塑性区分为裂尖塑性区(单调塑性区和
逆向塑性区)和尾迹场塑性区(遗留在尾迹场的单
调塑性变形区和循环塑性区)
[
17
]
两部分,见图
5.
图
5
过载情况下裂纹扩展时的塑性区
Fig. 5 Plastic zone induced by overload
3
2
7
西
南
交
通
大
学
学
报
第
53
卷
图
5
中:
β
为最大剪应力平面方位角,数学尖
裂纹为
70. 5°
,光滑试件为
45°
,钝裂纹为
45° ~
70. 5°
[
1718
]
;
d
为裂纹一个循环的扩展长度
.
循环载荷作用下,裂尖出现塑性变形,虽然这
是小范围屈服,但表明当应力从最小应力
σ
min
上升
到最大应力
σ
max
的过程中塑性区就产生了
.
当加载
到最大应力
σ
max
时,塑性区的尺寸与应力强度因子
K
max
和
S
y
比值的平方成正比
.
塑性区在加载方向产
生了塑性伸长,相比原来的尺寸有所增大,于是在
卸载过程中,这个塑性区受到压缩荷载作用,出现
了反向塑性
.
反向塑性的出现,要求在相反方向的
局部应力增量达到两倍
S
y
的量级,这意味着反向
塑性区的尺寸应该约为加载过程中所产生的塑性
区的
1 / 4.
所以反向塑性区比最大应力
σ
max
所产生
的塑性区要小很多
.
在反向塑性区之外,相对大得
多的塑性区变形仅发生在
σ
max
时,这称为单调塑
性,而正是这个单调塑性变形引起了材料在加载方
向上的永久性伸长
.
3. 1
长裂纹门槛值的裂尖塑性区半径
疲劳裂纹萌生和裂纹扩展是循环滑移的结果,
循环滑移意味着循环塑性变形(见图
5
)
.
在周期加
载下不可逆的位错滑移产生永久滑移带,表面晶粒
挤入或挤出,以最佳取向滑移
.
疲劳在低于屈服应
力的作用下发生,塑性变形限于材料中少量的晶
粒
.
随着滑移带增多和应力的增加,在裂尖形成塑
性区
I
型裂纹裂尖塑性区半径
r
pl
和逆向塑性区半
径
r
rpl
分别为
[
19
]
r
pl
= K
2
I
/
(
2pS
2
y
),
r
rpl
= K
2
I
/
(
8pS
2
y
)
.
(
10
)
根据
Hobbacher
[
20
]
建议,
C
和
CM
结构钢
Paris
裂纹扩展区的下确界(
SIF
门槛值)
ΔK
th
与塑性区
尺寸的关系为
ΔK
th
=
6. 0
,
R = 0
,
6. 0 - 4. 6r
,
R > 0
{
,
(
11
)
式中:
R
为加载应力比,
R = 0. 2.
根据刘艳萍等
[
21
]
的建议,对
Q370qE
钢,
ΔK
th
= 5. 556
(
1 - 0. 825R
)
1. 147
.
(
12
)
对桥梁钢
Q345qD
,
ΔK
th
= 4. 52 MPa
·
m
0. 5
,计
算得到
r
pl
= 27. 29 μm
,
r
rpl
= 6. 82 μm.
El Haddad
等
[
22
]
建议使用
ΔK
th
和疲劳极限
σ
f
表达虚拟裂纹长度
a
,即
a = ΔK
th
/
(
πσ
f
Y
),
(
13
)
式中:
Y
为几何形状系数,取决于裂纹形态,对于穿
透裂纹(长度
2a
)的无限平板,
Y = 1. 0
,其它裂纹类
型
Y
值可以从应力强度因子求解
.
正交异性钢桥面板开孔处的常幅疲劳极限为
70 MPa
[
23
]
,根据式(
13
)可以计算初始裂纹长度
a
0
= 1. 262 mm
,
a
0
/ r
pl
= 46 > 25
,符合小范围屈服假定
.
3. 2
循环塑性和裂纹提前闭合机理
图
6
是利用
DGXFEM
计算平面应力和平面
应变条件下,循环塑性过程中塑性尾迹场中的材料
以不同方式向裂尖转移的情况
.
图中:
①
平面应
力:颈缩实现材料向裂尖转移;
②
平面应变:塑性
区单元剪切、变形、旋转,体积增大,实现材料向裂
尖转移;
③
单元剪切、变形伸长、旋转锲入模型
.
在
每一个载荷循环中都会发生裂纹尖端塑性,所以裂
纹扩展时要穿过以前荷载循环所产生的塑性区,因
此,塑性变形就被残留在裂纹尾迹场
.
塑性尾迹场
大部分材料在加载过程中经受了单调塑性变形
(见图
5
),而沿着裂纹面的一个很小的边缘则经受
了循环塑性变形
.
塑性尾迹场中的材料在载荷方向
上已发生塑性伸长(见图
6
),在邻近裂尖塑性区的
裂纹面发生收缩产生压应力
.
在卸载过程中,试样
完全卸载之前裂纹却能够闭合
.
这就是塑性诱发的
裂纹提前闭合
[
17
]
.
图
6
塑性尾迹场中材料向裂尖转移
Fig. 6 Transfer of material in wake field to crack tip
平面应力情况下,试件表现为侧向颈缩实现材
料转移,平面应变则通过单元剪切滑移、变形伸长
和旋转锲入的方式实现向裂尖的转移
.
当裂纹通过
塑性区的时候,遗留或填充的材料在裂纹两侧表面
上塑性累积,伸长的塑性累积产生压缩效应的尾迹
场,使裂尖提前闭合
.
3. 3
塑性区对过载循环的影响
在变幅疲劳问题中,一个过载能够使裂纹尖端
张开并且变钝,这将降低张开应力
S
op
,并且助长裂
纹开始向过载形成新的裂尖塑性区扩展,而且由于
在反向塑性区内存在的残余拉伸应力随后可能发
生某种程度的裂纹扩展加速,残余拉伸应力促使裂
纹尖端张开
.
在一个小的裂纹长度增量范围内观察
4
2
7
第
4
期
渠
昱,等:桥梁结构钢裂纹塑性区的研究及应用
到了这种紧随过载后的裂纹扩展加速现象
.
由于受
到因过载在单调塑性区内引起的残余压缩应力的
作用,裂纹进入过载塑性区后继续扩展将遇到越来
越严重的裂纹闭合,于是裂纹扩展速率下降
.
这种裂纹塑性诱发的裂纹闭合对过载情况下
的裂纹扩展具有明显的迟滞效应,是由于在先前产
生的裂纹尖端塑性区中,已发生塑性伸长的材料残
留在裂纹尾迹上引起的
.
高峰载荷引起了较大的塑
性区,当裂纹尖端进入这些塑性区时就出现更多的
裂纹闭合,这已被试验所证实
.
在裂纹进一步扩展
后,将在裂纹尾迹上留下更多的塑性变形,这使得
S
op
升高,从而降低有效应力强度因子幅
ΔK
eff
,导致
了裂纹扩展迟滞
.
在通过了裂纹扩展速率的最小值
并且相当程度的后续扩展以后,扩展速率恢复到正
常的常幅加载水平
.
在
S
op
降低到
S
op
= S
min
后,由于
在后续的常幅循环中不再发生裂纹闭合,于是裂纹
扩展迟滞结束
.
4
结
论
(
1
)裂尖塑性区尺寸与应力水平的平方成正
比,扩展有限元法与
Westergoard
函数法的计算结
果基本一致
.
(
2
)正交异性钢桥面板开孔处,材料
Q345qD
取
R = 0. 2
,常幅疲劳极限
70 MPa
,
ΔK
th
= 4. 52 MPa
·
m
0. 5
,
计算得到
r
pl
= 27. 29 μm
,
r
rpl
= 6. 82 μm
,
a
0
=
1. 262 mm.
(
3
)对称向前伸展的裂尖塑性区形成的遮蔽
效应,使裂纹尾迹场避免了裂尖应力场的拉伸作
用,有利于裂纹闭合
.
(
4
)
DGXFEM
对裂纹表面循环塑性变形进行
模拟,表明平面应力情况下,面外以颈缩的方式实
现材料往裂尖的转移;平面应变情况下,材料以剪
切变形、旋转伸长和锲入实现材料往裂尖转移
.
当
裂纹通过塑性区的时候,残留
/
填充的材料在裂纹
表面产生塑性累积,循环塑性区内的单元伸长产生
压缩效应使裂纹提前闭合
.
(
5
)塑性区对过载循环情况下裂纹扩展滞后
有重要影响
.
致谢:山区桥梁与隧道工程国家重点实验室培育基
地开放基金资助
.
参考文献
:
[
1
]
李传习,李游,陈卓异,等
.
钢箱梁横隔板疲劳开裂原
因及 补 强 细 节 研 究 [
J
]
.
中 国 公 路 学 报,
2017
,
30
(
3
):
121131.
LI Chuanxi
,
LI You
,
CHEN Zhuoyi
,
et al. Fatigue
cracking reason and detail dimension of reinforcement
about transverse diaphragm of steel box bridge
[
J
]
.
China Journal of Highway and Transport. 2017
,
30
(
3
):
121131.
[
2
]
张清华,卜一之,李乔
.
正交异性钢桥面板疲劳问题
的研究进展[
J
]
.
中国公路学报,
2017
,
30
(
3
):
14
30.
ZHANG Qinghua
,
BU Yizhi
,
LI Qiao. Review on
fatigue problems of orthotropic steel bridge deck
[
J
]
.
China Journal of Highway and Transport
,
2017
,
30
(
3
):
1430.
[
3
]
CARROLL J D
,
ABUZAID W Z
,
LAMBROS J
,
et al.
On the interactions between strain accumulation
,
microstructure
,
and fatigue crack behavior
[
J
]
.
International Journal of Fracture
,
2013
,
180
(
2
):
223
241.
[
4
]
MIAO Jiashi
,
POLLOCK T M
,
WAYNE JONES J.
Microstructural extremes and the transition from fatigue
crack initiation to small crack growth in a polycrystalline
nickelbase superalloy
[
J
]
. Acta Materialia
,
2012
,
60
(
6 / 7
):
28402854.
[
5
]
GROSS D W
,
NYGREN K
,
PATAKY G J
,
et al. The
evolved microstructure ahead of an arrested fatigue crack
in Haynes 230
[
J
]
. Acta Materialia
,
2013
,
61
(
15
):
57685778.
[
6
]
PILCHAK A L. Fatigue crack growth rates in alpha
titanium
:
faceted vs. striation growth
[
J
]
. Scripta
Materialia
,
2013
,
68
(
5
):
277280.
[
7
]
HERBIG M
,
KINGM A
,
REISCHIG P
,
et al. 3D
growth of a short fatigue crack within a polycrystalline
microstructure studied using combined diffraction and
phasecontrast Xray tomography
[
J
]
. Acta Materialia
,
2011
,
59
(
2
):
590601.
[
8
]
WILLIAMS J J
,
YAZZIE K E
,
PADILLA E
,
et al.
Understanding fatigue crack growth in aluminum alloys
by in situ Xray synchrotron to mography
[
J
]
.
International Journal of Fatigue
,
2013
,
57
(
12
):
7985.
[
9
]
陈景杰,黄一,李玉刚
.
压应力对裂尖塑性区影响的
量化研究[
J
]
.
华中科技大学学报:自然科学版,
2016
(
2
):
102106.
CHEN Jingjie
,
HUANG Yi
,
LI Yugang. Quantitative
study on the influence of compressive stress on crack tip
plastic zone
[
J
]
. Journal of Huazhong University of
Science and Technology
:
Nature Science Edition
,
2016
(
2
):
102106.
5
2
7
西
南
交
通
大
学
学
报
第
53
卷
[
10
]
赵均海,魏雪英,马淑芬
.
混凝土结构
I
型裂纹裂尖
塑性区研究[
J
]
.
工程力学,
2006
,
23
(
9
):
141145.
ZHAO Junhai
,
WEI Xueying
,
MA Shufen. Crack tip
plastic zone of type I cracks for concrete structure
[
J
]
.
Engineering Machanics
,
2006
,
23
(
9
):
141145.
[
11
]
VASCO OLMO J M
,
JAMESM M N
,
CHRISTOPHER
C J
,
et al. Assessment of crack tip plastic zone size
and shape and its influence on crack tip shielding
[
J
]
.
Fatigue & Fracture of Engineering Materials &
Structures
,
2016
,
39
(
8
):
969981.
[
12
]
IRWIN G R. Analysis of stresses and strains near the
end of a crack traversing a plate
[
J
]
. Journal of
Applied Mechanics
,
1957
,
24
:
361364.
[
13
]
SANFORD R J. Principles of fracture mechanics
[
M
]
.
New Jersey
:
Pearson Education
,
2002
:
2944.
[
14
]
SOUSA R A
,
CASTRO J T P
,
LOPES A A O
,
et al
,
On improved crack tip plastic zone estimates based on
Tstress and on complete stress fields
[
J
]
. Fatigue &
Fracture of Engineering Materials & Structures
,
2012
,
36
(
1
):
2538.
[
15
]
BELYTSCHKO T
,
GRACIE R
,
VENTURA G. A
review of extended generalized finite element methods
for material modeling
[
EB / OL
]
.
[
20170512
]
.
http
:
/ / iopscience. iop. org / article / 10. 1088 / 0965
0393 / 17 / 4 / 043001 / meta.
[
16
]
GRACIE R
,
WANG Hongwu
,
BELYTSCHKO T.
Blending in the extended finite element method by
discontinuous Galerkin and assumed strain
methods
[
J
]
. International Journal for Numerical
Methods in Engineering
,
2008
,
74
(
11
):
16451669.
[
17
]
PIPPAN R
,
HOHENWARTER A. Fatigue crack
closure
:
a review of the physical phenomena
[
J
]
.
Fatigue & Fracture of Engineering Materials &
Structures
,
2017
,
40
(
4
):
471495.
[
18
]
XIANG Yibing
,
LU Zizi
,
LIU Yongming. Crack
growthbased fatigue life prediction using an equivalent
initial flaw model. part Ⅰ
:
uniaxial loading
[
J
]
.
International Journal of Fracture
,
2010
,
32
(
2
):
341
349.
[
19
]
IACOVIELLO F
,
COCCO V D
,
CAVALLINI M.
Fatigue crack propagation and overload damaging
micromechanisms in a ferriticpearlitic ductile cast
iron
[
J
]
. Fatigue & Fracture of Engineering Materials
& Structures
,
2016
,
39
(
8
):
9991011.
[
20
]
HOBBACHER A. XIII2151r107 / XV1254r107
Recommendation for fatigue design of welded joints and
components
[
S
]
. Paris
:
International Institute of
Welding
,
2007.
[
21
]
刘艳萍,陈传尧,李建兵,等
. 14MnNbq
焊接桥梁钢
的疲劳裂纹扩展行为研究[
J
]
.
工程力学,
2008
,
25
(
4
):
209213.
LIU Yanping
,
CHEN Chuanyao
,
LI Jianbing
,
et al.
Fracture crack growth behavior for the weld heat
affected zone of 14MnNbq bridge steel
[
J
]
.
Engineering Mechanics
,
2008
,
25
(
4
):
209213.
[
22
]
El HADDAD M H
,
TOPPER T H
,
SMITH K N.
Prediction of non propagating cracks
[
J
]
. Engineering
Fracture Mechanics
,
1979
,
11
(
3
):
573584.
[
23
]
中交公路规划设计院有限公司
. JTG D64
—
2015
公
路钢结构桥梁设计规范[
S
]
.
北京:人民交通出版
社股份有限公司,
2015.
(编辑:郭菊彬)
6
2
7