二维裂纹在线弹性及弹塑性应力场下J积分计算.pdfVIP

二二维裂纹在线弹性及弹塑性应力场下的J积分计算 李强 (中困原iF能科学甜f究院，102413) N^ 摘要：本文按照断裂力学中埘J积分的定义．为数值”算韭{j全面的理睦推导．编写程辟．与^DI 程厅接口’实现’}t面问题中、线弹{叫发弹塑忡J租分的王出潋值”锋，以裂纹拉仲武村为算例对 程片进行校核，用本文的程序研究小同的私!分路释列J扒分汁算结米的影响，对小H的汁算成力强 度压子的方值逃{J比较． 关键词：弹塑性，J积分，数值计算，AD川A 0引言 工程结构的破坏、故障和事故，有很大一部分是由于结构部件的断裂引起的。据 我【目劳动部统计，我国在80年代发牛的锅炉和』土力容器的爆炸事故约五千起人员累 计伤亡近力人，居国内劳动安全事故的第二位“。断裂力学的研究可对防范上述行为产 生萤大作用a对于由韧P芏和塑性较好的材料制造的各种设备，基于线弹性断裂力学的 应，勺强度因子理论的应用受到限制，『『iiJ积分理论日益被广泛应用口l。 1 J积分的定义 不，：。般性，如图l i所示，设 有一平面穿透裂纹体，裂纹表面无力 ／，一 作用，外加载荷使此裂纹周围产生一 ■≮nY 维应力(应变)场。若围绕裂纹尖端 任意取凹路I，始j：裂纹下表面，终 于裂纹上表向，按逆时针转功。于是 J积分的定义如下： ＼ ／j’ (Ⅱ，p=Y．z) 。，=I(¨”}一门。盯。俨?√卢J)(如， 其中w是应变能密度， cr,,／，dE“㈦胪Y，， (I2) w=¨(r。，)=r -维裂纹在线弹性及弹塑性应力场p的J积分刊算 !o 2 J积分的数值计算 2．1 J积分路径的选取 本文在数值计算th取积分路径沿单元高斯积分点，即选定的遗书高斯积分^ 连线构成了积分路径I，这样呵利崩高斯积分公式刈式(【l1求拟分 2．2高斯积分点上u。值的求解 2 2 1坐标变换借助于形函数‘jI来确定星本单7i内仟何一电(：，1 1 R0i体坐标值： = A f!Ia ，， 声一 咖 ，， ， 1 = ! Hy一=∑= Ⅳ f( f， 呜 2 b) 2 2 ， 2位移根据单元节点位移，fⅢ约J：杉幽数来确；之基本单7c内任何幢 的f々移： ”(锄)=∑．A’胪，咖i， ㈦训 _ ”’(锄)=∑．Ⅳ∥，，小tj 猢 2 2 3 LI¨的计算从式(二2)得- 尘：；，0i。 一Ⅻ 4 i：I毋’‘ 窘=蓦》 ∞㈨ 和j