弹塑性材料中裂纹的仿真研究.pdfVIP

第l9卷 5期 应用基础与工程科学学报 V0l_19．No．5 2011年10月 JOURNALOFBASICSCIENCEANDENGINEERING October2011 文章编号：1005-0930(2011)05-0810-007 中图分类号 ：034 文献标识码：A doi：10．3969／j．issn．1005-0930．2011．05．014 弹塑性材料中裂纹的仿真研究 李 成，铁 瑛， 郑艳萍 (郑州大学机械工程学院，河南 郑州 450001) 摘要：针对含裂纹的弹塑性材料结构，根据弹塑性断裂理论分析了裂纹增长的条 件．建立了基于裂纹尖端存在塑性区的计算模型，得到了描述含裂纹弹塑性材料 在外荷载作用下的塑性区尺寸及塑性区前端裂纹张开位移的解析解．按照所建 立的计算模型对不同的裂纹长度、不同的外荷载对塑性区尺寸及张开位移的影 响进行探讨．并给出了材料在外荷载作用下破坏的临界应力强度因子曲线． 关键词：弹塑性材料；裂纹；仿真；塑性区尺寸；张开位移；临界应力强度因子 线弹性断裂力学是把材料当作完全弹性体，运用弹性力学的理论方法研究裂纹的开 裂和扩展规律，从而提出裂纹临界扩展的判据．不仅如此，由于裂纹尖端的高度应力集中， 在裂纹尖端附近存在着塑性区，但如果这个塑性区的尺寸和裂纹尺寸相比小得多，这个微 小的塑性区被大范围的弹性区所包围，这时仍可以用线弹性断裂力学来研究裂纹的扩展 情况．但对于一些所谓的韧性材料，在裂纹的尖端往往会形成明显的塑性区，随着外力的 增大，裂纹尖端不断钝化，出现裂纹逐渐张开的现象，同时塑性区也在逐渐增大，对于这种 情况，线弹性断裂理论将不再适用．用应力强度因子来描述裂纹的扩展也将失去意义．这 种塑性变形占较大比重的断裂问题就要用弹塑性断裂理论来解决．目前应用较为广泛的 是J积分理论、COD原理…等方法．理论上，结构的裂纹问题的求解、分析，数学上要归结 为边值问题的建立和求解 J．近几十年来，弹塑性裂纹问题得到了一定程度的研究，一些 研究结果已经发表 引．但弹塑性断裂理论远不如线弹性断裂理论成熟，还在不断的发展 中．已有的文献主要在其它条件已知的情况下对塑性区的尺寸和裂纹张开位移进行解析 计算，而对裂纹的初始长度和外荷载对塑性区的影响，以及外荷载对临界应力强度因子的 影响的解析研究较少．关于这方面的研究主要是用实验方法和数值方法来进行 引．但在 采用有限元数值方法对裂纹进行分析时，由于裂纹尖端存在奇异单元，所得到的结果往往 误差较大．采用弹粘塑性力学模型，文献[9．10]分别对稳恒扩展裂纹尖端的弹粘塑性场和 动态裂纹尖端弹粘塑性场进行了研究，获得了平面应变动态扩展裂纹尖端的指数奇异性 场和对数奇异性场，取得了较好的分析结果． 本文根据弹塑性断裂理论建立计算模型．采用相应的解析计算方法对不同的初始裂 纹长度、不同的外荷载对塑性区尺寸及张开位移的影响进行全面的仿真研究．得到了塑性 收稿 日期：2009-03-09；修订 日期：2011-05—17 基金项 目：国家 自然科学基金资助项 目；教育部留学回国人员科研启动基金资助项 目(208082) 作者简介：李 成 (1962一)，男，博士，教授．TelE—mail：chengli@ZZU．edu．cn 李 成等：弹塑性材料中裂纹的仿真研究 811 区尺寸随外荷载变化的情况．给出了临界荷载和裂纹初始长度的关系以及材料在外荷载 作用下破坏的临界应力强度因子曲线． 1 力学模型与基本方程 在对带有裂纹的这类韧性材料试验时发现，裂纹尖端的塑性区呈扁平带状 [1]，如图1 所示，这就是所谓的D-M模型．其中，d为裂纹塑性区长度、0为裂纹初始长度，单位 mm． or，为塑性区上所作用的分布应力．基于这个模型可以较好地处理具有穿透裂纹的板的弹 塑性问题．所谓张开位移，是指裂纹体受载后，在裂纹尖端沿垂直于裂纹的方向所产生的 位移(COD)，如上图中的 ．其临界值为 是材料的断裂韧性表征参数，为一常数，由试 验得到．在弹塑性断裂理论中，将张开位移6作为一个参量，而将临界张开位