刘倩等：基于Lemaitre模型的辊式冲裁工艺断面质量分析 ·1201· (a损伤参数 )损伤参数 (c)损伤参数 图1损伤分布与断面质量.(a)p=12.16MPa:(b)p=12.32MPa:(c)断后形貌 Fig.1 Damage distribution and sheared edge quality:(a)p=12.16 MPa:(b)p=12.32 MPa;(c)cross-section shape (a损伤参数 b)损伤参数 c)损伤参数 t6.256x10- +4.260x10- +1.007×10 +5.370x10-2 +3.652×10 +8.621×10- +4.475×10-2 +3.043×10-1 +7.189×10- +350x10-- +2.435×10-1 +5.746×10-1 +2.685x10-2 +1.826x10 +4.309x10- +1790x10- 41217×10 +2.873×10- 8951×10-2 0只 43610-1 9.098×10-24 -9.098×10-4 .538×10-34 图2文献中的损伤分布与断面质量o],(a)p=12.16MPa;(b)p=12.32MPa;(c)断后形貌 Fig.2 Damage distribution and sheared edge quality in literature[o]:(a)p=12.16 MPa;(b)p=12.32 MPa;(c)cross-section shape a)◆载荷 (2)根据弹性变形阶段内应力与应变的关系，计 算屈服前初始杨氏模量E. (3)根据卸载时应力应变的变化，得到各个残余 应变对应的杨氏模量E,计算出每次卸载时损伤变量 D: D=1-E/E. (31) (4)假定损伤变量与真应变间为线性关系，对损 截荷 伤变量D与真应变ε进行线性拟合，如图4所示： 应变 D=bo+b18. (32) 图3弹性模量的测量方法.(a)拉伸试验：(b)应力-应变 式中：b。与b,为拟合系数 曲线 (5)当D=0时，计算损伤发生时的对数应变值 Fig.3 Measure of Young's modulus:(a)tensile experiment;(b) E。,即拟合直线与应变横坐标轴的交点： stress-strain curve Ep=-00/61. (33) 载拉伸试验，直至试样发生断裂，记录应力-应变关系 (6)对拉伸试验实体建模，进行数值模拟直至达 曲线. 到断裂时最大伸长量，见图5，由此获得变形区内最大刘 倩等: 基于 Lemaitre 模型的辊式冲裁工艺断面质量分析 图 1 损伤分布与断面质量. (a) p = 12郾 16 MPa; (b) p = 12郾 32 MPa; (c) 断后形貌 Fig. 1 Damage distribution and sheared edge quality: (a) p = 12郾 16 MPa; (b) p = 12郾 32 MPa; (c) cross鄄section shape 图 2 文献中的损伤分布与断面质量[10] . (a) p = 12郾 16 MPa; (b) p = 12郾 32 MPa; (c) 断后形貌 Fig. 2 Damage distribution and sheared edge quality in literature [10] : (a) p = 12郾 16 MPa; (b) p = 12郾 32 MPa; (c) cross鄄section shape 图 3 弹性模量的测量方法. ( a) 拉伸试验; ( b) 应力鄄鄄 应变 曲线 Fig. 3 Measure of Young爷s modulus: (a) tensile experiment; ( b) stress鄄鄄strain curve 载拉伸试验,直至试样发生断裂,记录应力鄄鄄应变关系 曲线. (2)根据弹性变形阶段内应力与应变的关系,计 算屈服前初始杨氏模量 E. (3)根据卸载时应力应变的变化,得到各个残余 应变对应的杨氏模量 E寛,计算出每次卸载时损伤变量 D: D = 1 - E寛/ E. (31) (4)假定损伤变量与真应变间为线性关系,对损 伤变量 D 与真应变 着 进行线性拟合,如图 4 所示: D = b0 + b1着. (32) 式中:b0 与 b1 为拟合系数. (5)当 D = 0 时,计算损伤发生时的对数应变值 着D ,即拟合直线与应变横坐标轴的交点: 着D = - b0 / b1 . (33) (6)对拉伸试验实体建模,进行数值模拟直至达 到断裂时最大伸长量,见图 5,由此获得变形区内最大 ·1201·