matlab钢材切割,一种基于MATLAB的钢材裂纹扩展速率试验数据处理方法

第 33 卷第 6 期 武 汉 科 技 大 学 学 报 Vol. 33 , No. 6 2010 年 12 月 Journal of Wuhan University of Science and Technology Dec. 2010 　　收稿日期:20100401 　　作者简介:刘建亮(1984) ,男 ,第二炮兵工程学院硕士生. Email :chenfengsuibi @sohu. com 　　通讯作者:岳应娟(19722) ,女 ,第二炮兵工程学院副教授. E2mail :yingjuanyue @163. com 一种基于 MA TLAB 的钢材裂纹扩展速率试验数据处理方法 刘建亮 ,岳应娟 (第二炮兵工程学院 ,陕西 西安 ,710025) 摘要 :提出一种基于 MATLAB 软件设计的钢材裂纹扩展速率试验数据处理程序 ,并利用该程序对18MnMoNb 钢和 16MnR 钢裂纹扩展速率试验数据进行处理。在分析程序运行结果的基础上 ,评价七点递增法、五点递增法、三点递增法针对不同试验数据处理精度的优劣性 ,并分析 Paris 公式中常数 C、m 的相关性 , 得出了 C、m 的相关性经验公式。 关键词 :MATLAB ;裂纹扩展速率 ;试验数据处理程序 ;递增法 中图分类号 : TG405 　　文献标志码 :A 　　文章编号 :167423644(2010)0605902051 　问题的提出 疲劳裂纹的扩展大致可以分为 3 个阶段[1] , 如图 1 所示。在第一阶段内 ,裂纹前沿应力场强度因子幅ΔK 值较低 ,裂纹扩展速率 da/ dN ( a 为裂纹长度; N 为载荷循环次数) 也较低 ,当ΔK 低于疲劳裂纹扩展的临界应力场强度范围因子幅 , 即应力门槛值ΔKth时 ,da/ dN = 0 ,裂纹基本上不扩展。第二阶段是疲劳裂纹扩展的稳定阶段 ,也是决定疲劳寿命的主要阶段。一般认为 ,在双对数坐标系中 ,该阶段 da/ dN 与ΔK 的关系是一条直线。用于描述这一阶段的 da/ dN 的关系式有数十个 ,其中用得最广泛而且最简单的是 Paris 提出的表达式[1] : da dN = C(ΔK) m (1) 　　式(1) 中 , C、m 是与试验条件如环境、温度、频率、应力比等有关的材料常数 ,可由标准的疲劳小试样获得 ,在概率断裂力学的计算中 , C、m常 图 1 　 ln( da/ dN)与 lnΔK的关系曲线 Fig. 1 Relationship between ln( da/ dN) and lnΔK 视为随机变量 ,不少文献认为 lnC 与 m 具有相关 性[223] 。 在双对数坐标系中 ,式(1)的形式为 ln da dN = lnC + mlnΔK (2) 　　若ΔK 继续增大 ,则进入疲劳裂纹扩展第三阶段 ,这时最大循环应力 Kmax已接近材料的平面应变断裂韧度 K1C ,随着裂纹的扩展 ,ΔK 亦迅速增大并导致材料失稳断裂。与裂纹扩展第一阶段一样 ,第三阶段仅占疲劳寿命极少部分。 由此可见 ,在材料的疲劳分析中 ,第二阶段的分析是关键。有多种方法测定材料的 da/ dN ,具体方法可以参照文献[4 ]。该阶段疲劳分析的基本思路 :按照标准设计裂纹扩展速率试验 ,测得( ai , N i ) ,求得 ( da/ dN ) i 和 ΔKi , 最后对 ( da/ dN) i 和ΔKi 做幂乘回归 ,即可求得材料常数 C 和 m 。 要获得可靠的材料常数 C 和 m ,首先要在试验中提高测量精度 ,其次就是要有一