裂纹尖端断裂力学参数计算
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
❖ 20世纪50年代，美国北极星导弹固体燃料发动机 发射时发生低应力脆断。
❖ 1965年，英国某大型合成塔在水压试验时断裂成 两段。
事故调查发现 →断裂起源于构件中裂纹
1、裂纹尖端断裂力学参数研究意义
❖ 传统的强度理论
缺陷：传统强度理论并没有考虑材料中是否有缺陷， 对有缺陷的材料，对其安全可靠性不能做出正确的判 断。
14 13.5
位移外推 KI
解析KI
h/b
13
h:b
0.75 1 2 4
结论
❖ 验证了1/4节点处理裂纹尖端奇异性是可以的。 ❖ 在数值法计算中，随着平板尺寸的增大，KI的值
逐渐接近于解析值。
谢谢大家！
基于ANSYS的裂纹尖端应力强度 因子KI的计算
❖ 裂纹尖端应力强度KI研究的意义 ❖ 裂纹尖端KI的计算方法 ❖ 裂纹尖端应力奇异性处理 ❖ Ansys计算过程及结果
1、裂纹尖端断裂力学参数研究意义
❖ 随着现代高强材料和大型结构的广泛应用，一些 按传统强度理论和常规方法设计、制造的产品， 发生了不少重大断裂事故。
solid95 远离裂纹采用实体 建模 solid45
4、ansys计算过程及结果
❖ 载荷及约束
4、ansys计算过程及结果
❖ 应力云图
最大应力420MPA
4、ansys计算过程及结果
不同载荷下的应力强度因子
KI 16 MPa m 14
12 10
8 6 4 2 0
10 20 30 40 50
位移外推KI 解析法KI
位移外推法
裂纹前缘 u K I 2r (1 v)
(1)
G
变换(1)得
KI
lim
r 0
G u 2(1- ) r
(2)
KI
4 3 KIB
1 - 3 KIC
2 E 4(1- 2 )
(8uB 3
uc ) L
(3)
缺点：取决于节点的位置和其位移值精度
Chen LS and Kuang JH. A modified linear extrapolation formula for determination of stress intensity
2、裂纹尖端KI的计算方法
KI的计算方法 解析法 数值法 位移外推法 J积分法
2、裂纹尖端KI的计算方法
解析法
K I a f (a, w, )
为裂纹尖端无无裂纹时的应力
f(a,w,…)为几何修正系数 缺陷：适用于几何简单的板类，杆类，梁类构 件；对于较复杂得构件，无法得到正确的解析解 。
2、裂纹尖端KI的计算方法
2、裂纹尖端KI的计算方法
J积分法
J积分定义为与路径无关的曲线积分
tx xnx xyny ty yny xynx
tx, ty 分别为X,Y轴的引力分量 n为积分路劲上的单位法向量
间接求得
KI
JE
1- 2
缺点：只能应用于穿透性裂纹，对于表面椭 圆裂纹，剪切滑移等裂纹根本无法计算。
3、裂纹尖端应力奇异性的处理
载荷：MPa
4、ansys计算过程及结果
改变平板的宽度使得 b:a=2.5，3.3，5，10
KI 16
15.5
15 14.5
14
位移外推 KI
解析解KI
b/a
13.5
13 2.5 3.3 5
b:a
10
4、ansys计算过程及结果
改变平板长度h:b=0.75，1，2，4
KI 16
15.5
15 14.5
K I a
ij K I K I 越大， ij就越大
K反映了裂纹尖端应力场的强弱程度
1、裂纹尖端断裂力学参数研究意义
❖ K断裂准则
K KIC
Fra Baidu bibliotek
K IC 为材料的断裂韧性
(1)确定含裂纹构件的临界载荷。G，a，KIC → Fc
(2) 确定裂纹的极限尺寸。G，F，KIC → a
(3) 确定带裂纹构件的安全性。
3、裂纹尖端应力奇异性的处理
❖ 经奇异化处理的裂纹尖端单元划分
4、ansys计算过程及结果
❖ 中心穿透裂纹
几何参数：a=25.2mm；b=126mm； t=6.3mm，2h=2b
材料参数：E 206 Gpa ; 0.3
载荷： 0 50MPa
4、ansys计算过程及结果
❖ 模型建立
(1) 逐节点直接建模方法 (2)实体建模方法 本文：裂纹尖端采用逐节点
1、裂纹尖端断裂力学参数研究意义
❖ 工程中常见的几种裂纹
ij K I
1、裂纹尖端断裂力学参数研究意义
y
y
x
K I cos 1 sin sin 3
2r 2
2 2
2a
xy
dy
x
r dx
x
y
K I cos 1 sin sin 3
2r 2
2 2
xy
K I cos sin cos 3 2r 2 2 2
3、裂纹尖端应力奇异性的处理
L
H/4 B
A L/4
当 xB L / 2
H
C x 0.5(1)L u (1)[(1)uB 0.5uC ] (8)
xx
u x
u
x
1 L
[-4u
B
(2
1)u C ]
(9)
由(9)式知 应力，应变不具有奇异性 因此，裂纹尖端采用1/4节点等参二次单元是合理的。
u (1)[(1)uB 0.5uC ]
(5)
3、裂纹尖端应力奇异性的处理
位移u对坐标x进行微分：
xx
u x
u
x
1 x
1 L
[2uB
(
1 2
)uC
]
(6)
在极坐标中 x r, 0 故(6)变为
xx
u x
u
x
1 r
1 L
[2uB
(
1 2
)uC
]
(7)
根据材料的本构关系，应力与应变成正比，故应力也与 1/ r 项成比例
❖ 由上知，裂纹尖端各点应力表达式中均含有1/ r 项 r→0，尖端应力∞，裂纹尖端具有奇异性。
对这种奇异性问题，用常规的单元，势必得不到 精确的解；一般采用1/4节点的二次等参单元。
L
1 2 H/4 L/4
H 3
3、裂纹尖端应力奇异性的处理
L
A B H/4 L/4
H C
XA,XB,XC为A，B，C节点的坐标 ， uA,uB,uC分别为三节点的水平位移。
裂纹线上任意一点的坐标x和 位移u都可以用形函数插值为：
x 0.5(1 )xA (1 )(1 )xB 0.5(1 )xC
u 0.5(1 )uA (1 )(1 )uB 0.5(1 )uC
(4)
1 1
其中 xA 0 xB L / 4 xC L uA 0 代入(4)式化简
x 0.25(1)2 L
❖ 1965年，英国某大型合成塔在水压试验时断裂成 两段。
事故调查发现 →断裂起源于构件中裂纹
1、裂纹尖端断裂力学参数研究意义
❖ 传统的强度理论
缺陷：传统强度理论并没有考虑材料中是否有缺陷， 对有缺陷的材料，对其安全可靠性不能做出正确的判 断。
14 13.5
位移外推 KI
解析KI
h/b
13
h:b
0.75 1 2 4
结论
❖ 验证了1/4节点处理裂纹尖端奇异性是可以的。 ❖ 在数值法计算中，随着平板尺寸的增大，KI的值
逐渐接近于解析值。
谢谢大家！
基于ANSYS的裂纹尖端应力强度 因子KI的计算
❖ 裂纹尖端应力强度KI研究的意义 ❖ 裂纹尖端KI的计算方法 ❖ 裂纹尖端应力奇异性处理 ❖ Ansys计算过程及结果
1、裂纹尖端断裂力学参数研究意义
❖ 随着现代高强材料和大型结构的广泛应用，一些 按传统强度理论和常规方法设计、制造的产品， 发生了不少重大断裂事故。
solid95 远离裂纹采用实体 建模 solid45
4、ansys计算过程及结果
❖ 载荷及约束
4、ansys计算过程及结果
❖ 应力云图
最大应力420MPA
4、ansys计算过程及结果
不同载荷下的应力强度因子
KI 16 MPa m 14
12 10
8 6 4 2 0
10 20 30 40 50
位移外推KI 解析法KI
位移外推法
裂纹前缘 u K I 2r (1 v)
(1)
G
变换(1)得
KI
lim
r 0
G u 2(1- ) r
(2)
KI
4 3 KIB
1 - 3 KIC
2 E 4(1- 2 )
(8uB 3
uc ) L
(3)
缺点：取决于节点的位置和其位移值精度
Chen LS and Kuang JH. A modified linear extrapolation formula for determination of stress intensity
2、裂纹尖端KI的计算方法
KI的计算方法 解析法 数值法 位移外推法 J积分法
2、裂纹尖端KI的计算方法
解析法
K I a f (a, w, )
为裂纹尖端无无裂纹时的应力
f(a,w,…)为几何修正系数 缺陷：适用于几何简单的板类，杆类，梁类构 件；对于较复杂得构件，无法得到正确的解析解 。
2、裂纹尖端KI的计算方法
2、裂纹尖端KI的计算方法
J积分法
J积分定义为与路径无关的曲线积分
tx xnx xyny ty yny xynx
tx, ty 分别为X,Y轴的引力分量 n为积分路劲上的单位法向量
间接求得
KI
JE
1- 2
缺点：只能应用于穿透性裂纹，对于表面椭 圆裂纹，剪切滑移等裂纹根本无法计算。
3、裂纹尖端应力奇异性的处理
载荷：MPa
4、ansys计算过程及结果
改变平板的宽度使得 b:a=2.5，3.3，5，10
KI 16
15.5
15 14.5
14
位移外推 KI
解析解KI
b/a
13.5
13 2.5 3.3 5
b:a
10
4、ansys计算过程及结果
改变平板长度h:b=0.75，1，2，4
KI 16
15.5
15 14.5
K I a
ij K I K I 越大， ij就越大
K反映了裂纹尖端应力场的强弱程度
1、裂纹尖端断裂力学参数研究意义
❖ K断裂准则
K KIC
Fra Baidu bibliotek
K IC 为材料的断裂韧性
(1)确定含裂纹构件的临界载荷。G，a，KIC → Fc
(2) 确定裂纹的极限尺寸。G，F，KIC → a
(3) 确定带裂纹构件的安全性。
3、裂纹尖端应力奇异性的处理
❖ 经奇异化处理的裂纹尖端单元划分
4、ansys计算过程及结果
❖ 中心穿透裂纹
几何参数：a=25.2mm；b=126mm； t=6.3mm，2h=2b
材料参数：E 206 Gpa ; 0.3
载荷： 0 50MPa
4、ansys计算过程及结果
❖ 模型建立
(1) 逐节点直接建模方法 (2)实体建模方法 本文：裂纹尖端采用逐节点
1、裂纹尖端断裂力学参数研究意义
❖ 工程中常见的几种裂纹
ij K I
1、裂纹尖端断裂力学参数研究意义
y
y
x
K I cos 1 sin sin 3
2r 2
2 2
2a
xy
dy
x
r dx
x
y
K I cos 1 sin sin 3
2r 2
2 2
xy
K I cos sin cos 3 2r 2 2 2
3、裂纹尖端应力奇异性的处理
L
H/4 B
A L/4
当 xB L / 2
H
C x 0.5(1)L u (1)[(1)uB 0.5uC ] (8)
xx
u x
u
x
1 L
[-4u
B
(2
1)u C ]
(9)
由(9)式知 应力，应变不具有奇异性 因此，裂纹尖端采用1/4节点等参二次单元是合理的。
u (1)[(1)uB 0.5uC ]
(5)
3、裂纹尖端应力奇异性的处理
位移u对坐标x进行微分：
xx
u x
u
x
1 x
1 L
[2uB
(
1 2
)uC
]
(6)
在极坐标中 x r, 0 故(6)变为
xx
u x
u
x
1 r
1 L
[2uB
(
1 2
)uC
]
(7)
根据材料的本构关系，应力与应变成正比，故应力也与 1/ r 项成比例
❖ 由上知，裂纹尖端各点应力表达式中均含有1/ r 项 r→0，尖端应力∞，裂纹尖端具有奇异性。
对这种奇异性问题，用常规的单元，势必得不到 精确的解；一般采用1/4节点的二次等参单元。
L
1 2 H/4 L/4
H 3
3、裂纹尖端应力奇异性的处理
L
A B H/4 L/4
H C
XA,XB,XC为A，B，C节点的坐标 ， uA,uB,uC分别为三节点的水平位移。
裂纹线上任意一点的坐标x和 位移u都可以用形函数插值为：
x 0.5(1 )xA (1 )(1 )xB 0.5(1 )xC
u 0.5(1 )uA (1 )(1 )uB 0.5(1 )uC
(4)
1 1
其中 xA 0 xB L / 4 xC L uA 0 代入(4)式化简
x 0.25(1)2 L