CN103698188B - 慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法 - Google Patents
慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法 Download PDFInfo
- Publication number
- CN103698188B CN103698188B CN201410012971.8A CN201410012971A CN103698188B CN 103698188 B CN103698188 B CN 103698188B CN 201410012971 A CN201410012971 A CN 201410012971A CN 103698188 B CN103698188 B CN 103698188B
- Authority
- CN
- China
- Prior art keywords
- load
- ing
- point displacement
- loading
- crack
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法,包括步骤:一、紧凑拉伸试样预处理;二、采用慢拉升试验机进行慢应变速率拉伸试验,并绘制出加载力P与加载点位移ΔL之间的实验数据关系曲线图;三、对步骤二进行的慢应变速率拉伸试验进行模拟,并绘制出加载力P与加载点位移ΔL之间的理论分析数据关系曲线图;四、获取裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,所对应的试验时间t0、t1、t2、…、tn;五、获得裂纹扩展速率。本发明能将裂纹扩展造成的加载力P和加载点位移ΔL的变化分离出来,能够实时预测Δai的变化,获得裂纹的扩展速率,且能有效减少环境影响带来的测试数据失真,提高数据精度。
Description
技术领域
本发明属于断裂力学技术领域,尤其是涉及一种慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法。
背景技术
裂纹萌生、扩展和断裂是影响重要机械结构寿命和安全服役的重要因素之一,涉及到航空、航天、交通运输、化工、机械、材料、能源等多个工程领域。实验测试结构材料裂纹扩展特性是本项研究的重要手段之一;
为了提高核电设备的耐腐蚀性,轻水核反应堆的结构材料大量选用奥氏体不锈钢和镍基合金。实践和研究表明,奥氏体不锈钢和镍基合金材料在高温高压水环境中会产生以应力腐蚀开裂(SCC)为代表的环境致裂(EAC)。为了保证核电设备的安全运行,实验测定核电结构材料在高温水环境中环境致裂裂纹扩展速率是该领域的一项重要工作;
紧凑拉伸(CT)试样是断裂力学的标准试样之一,也是高温高压水环境环境中测定核电结构材料环境致裂裂纹扩展速率的最常用试样,该实验是将一个CT试样放置在一个具有模拟核压力容器高温高压水环境的高压釜中,并施加一定的载荷,实时动态观测该材料环境致裂裂纹扩展速率,整个实验过程可持续数个月或数年之久,耗时、耗资均较大;
为了全面评价材料的性质与行为,除须知道裂纹的起始长度外,精确测定在缓慢或稳定增长的动过程中裂纹长度是十分必要的。通常,在实验室中用显微镜跟踪或照相的光学方法确定裂纹增长,有它的局限性,这就是记录到的裂纹很可能不代表试件内部的实际情形,因为它是在自由表面上观察到的。另一方面,在高温高压水环境中核电结构材料环境致裂实验中,试样置于高压釜中,并不外露,用光学仪器做观察则是不可能的,而且,材料的环境致裂裂纹扩展速率试验要持续较长的时间,直接观察与记录也是不现实的;
电位差法,包括交流电位差法(ACPD)和直流电位差法(DCPD),是目前实时动态观测核电结构材料环境致裂裂纹扩展速率的最主要方法。其基本原理是一定值的电流通过含裂纹试样时,裂纹长度不同会造成电位分布的不同,通过微伏计的测量与换算,可以得到裂纹的扩展状况。但由于该方法是利用裂纹周围电位的微小变化来测定裂纹的扩展情况,影响因素较多,任何干扰都可能造成监测结果的误差;
慢应变速率拉伸(SSRT)为评价材料在特定介质中的应力腐蚀敏感性而设计的拉伸速率小于某一临界值的拉伸试验方法。它是以一个恒定不变的或相当缓慢的应变速率对置于腐蚀环境的试样施加应力,目的是为了得到试样的在特定环境下的应力腐蚀敏感性,它可作为高温高压水环境中测定核电结构材料环境致裂裂纹扩展速率的加速试验,把自然环境条件下长时间的裂纹扩展过程在较短时间内的实验环境条件下呈现出来。
综合分析显示,高温水环境下核电结构材料环境致裂的慢应变速率拉伸试验中实时动态观测裂纹扩展速率是一项及其重要而艰难的工作,为了有效地提高观测精度,必须对现有的测试方法进行改进。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足,提供一种慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法,其方法简单,实现方便,能够获得裂纹的扩展速率dΔai/dti,且能够有效地减少环境影响因素带来的测试数据失真,提高获取数据的精度。
为解决上述技术问题,本发明采用的技术方案是:一种慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法,其特征在于该方法包括以下步骤:
步骤一、紧凑拉伸试样预处理:根据ASME标准选取一个用于慢应变速率拉伸试验的标准紧凑拉伸试样,并采用线切割和预制疲劳裂纹的方法在标准紧凑拉伸试样上设置用于初期引导裂纹扩展的预裂裂纹,所述预裂裂纹的长度为a;
步骤二、采用慢拉伸试验机进行慢应变速率拉伸试验,并绘制出加载力P与加载点位移ΔL之间的实验数据关系曲线图,其具体过程为:
步骤201、在慢拉伸试验机上设定慢应变速率拉伸试验的试验参数,包括加载点位移的变化速度S和试验总时间T;
步骤202、慢拉伸试验机开始对标准紧凑拉伸试样进行加载,且在加载过程中,慢拉伸试验机上的力传感器探头对加载力P进行实时检测并将所检测到的信号实时传输给与所述慢拉伸试验机连接的计算机,所述计算机采用Origin软件,以加载点位移ΔL为横坐标,加载力P为纵坐标,绘制出加载力P与加载点位移ΔL之间的实验数据关系曲线图并显示;同时,测量记录裂纹实时的扩展长度A;其中,ΔL=St且ΔL≦L,t为试验时间且t≦T,L为在试验总时间T内加载点位移总变化长度;
步骤三、对步骤二中进行的慢应变速率拉伸试验进行模拟,并绘制出加载力P与加载点位移ΔL之间的理论分析数据关系曲线图:在计算机上采用有限元分析软件对步骤二中进行的慢应变速率拉伸试验进行模拟,并分析慢应变速率拉伸试验中裂纹扩展长度递增变化的过程中,加载力P与加载点位移ΔL的关系,且在裂纹扩展长度分别为Δa0、Δa1、Δa2、…、Δan时,所述计算机采用Origin软件,以加载点位移ΔL为横坐标,加载力P为纵坐标,绘制出加载力P与加载点位移ΔL之间的理论分析数据关系曲线图并显示在同一幅图中;其中,Δan不大于步骤202中测量记录的裂纹实时的扩展长度A的最大值;
步骤四、获取裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,所对应的试验时间t0、t1、t2、…、tn:在计算机上采用Origin软件将步骤202中绘制的加载力P与加载点位移ΔL之间的实验数据关系曲线图和步骤三中绘制的在裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,加载力P与加载点位移ΔL之间的理论分析数据关系曲线图整合到同一幅图中,得到整合后的加载力P与加载点位移ΔL之间的关系曲线图,且在整合后的加载力P与加载点位移ΔL之间的关系曲线图中,提取出裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,加载力P与加载点位移ΔL之间的实验数据关系曲线与n条加载力P与加载点位移ΔL之间的理论分析数据关系曲线的交点Qi的坐标数据,得到了裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,所对应的加载点位移ΔL0、ΔL1、ΔL2、…、ΔLn,由进而计算得到了裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,所对应的试验时间t0、t1、t2、…、tn;其中,i的取值为0、1、2、…、n;
步骤五、获得裂纹扩展速率:在计算机上采用Origin软件,以试验时间ti为横坐标,裂纹扩展长度Δai为纵坐标,绘制出步骤四中得到的实验时间t0、t1、t2、…、tn和裂纹扩展长度Δa0、Δa1、Δa2、…、Δan对应的各个点,并拟合得到裂纹扩展长度Δai与试验时间ti之间的关系曲线,裂纹扩展长度Δai与试验时间ti之间的关系曲线上任意一点的切线斜率k=dΔai/dti即为当时的裂纹扩展速率,其中,i的取值为0、1、2、…、n。
上述的慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法,其特征在于:步骤三中所述有限元分析软件为ABAQUS软件。
上述的慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法,其特征在于:步骤三中在计算机上采用有限元分析软件对步骤二中进行的慢应变速率拉伸试验进行模拟,并分析慢应变速率拉伸试验中裂纹扩展长度递增变化的过程中,加载力P与加载点位移ΔL的关系的具体过程为:
步骤301、在ABAQUS软件中,建立带有长度为a的预裂裂纹的标准紧凑拉伸试样的三维模型;
步骤302、在步骤301中建立的标准紧凑拉伸试样的三维模型中,添加长度为Δa1的裂纹;
步骤303、在ABAQUS软件中,建立慢拉伸试验机对标准紧凑拉伸试样进行加载时所用的两个加载销的刚体模型;
步骤304、赋予标准紧凑拉伸试样的三维模型材料参数,所述材料参数包括材料的类型、材料的杨氏模量、材料的泊松比、材料的屈服强度、材料的硬化指数和材料的偏移指数;
步骤305、将标准紧凑拉伸试样和两个加载销生成装配体,并移动两个加载销的位置,使得两个加载销完全对称地装配在所述标准紧凑拉伸试样上,且使得两个加载销伸出到所述标准紧凑拉伸试样左右两侧的长度相等;
步骤306、创建第一分析步和第二分析步两个分析步,与ABAQUS软件自带的初始分析步构成三个分析步;其中,初始分析步用于建立两个加载销和标准紧凑拉伸试样之间的初始接触关系;第一分析步用于施加一个载荷使初始接触关系平稳建立,便于计算的收敛;第二分析步用于加载;
步骤307、在初始分析步中,设置两个加载销和标准紧凑拉伸试样之间的初始接触关系为硬接触关系,主接触表面选择为销的外表面,从接触表面选择为孔的内表面;并设置两个加载销除拉伸方向外的其余两个直线自由度和三个旋转自由度均为0,对两个加载销除拉伸方向的直线自由度外所有的自由度进行约束;
步骤308、在第一分析步中,对两个加载销沿拉伸方向分别施加一个大小为0.02mm~0.08mm的位移载荷;在第二分析步中,对两个加载销沿拉伸方向分别施加线性变化的大小为2mm~8mm位移载荷;
步骤309、对由标准紧凑拉伸试样和两个加载销生成的装配体进行网格划分,在裂纹尖端划分密致网格,其余位置划分稀疏网格,并将网格的单元类型选择为三维应力六面体一次单元;
步骤310、提交计算,ABAQUS软件按照步骤301~309设置的参数进行分析计算,模拟慢应变速率拉伸试验中的裂纹扩展过程;
步骤311、裂纹扩展过程中,进行数据提取,先提取加载销的加载点位移ΔL随时间变化的数据,再提取加载力P随时间变化的数据,得到了加载力P与加载点位移ΔL的关系。
上述的慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法,其特征在于:步骤三中裂纹扩展长度Δai=iΔx,其中,Δx为裂纹单次扩展长度,i的取值为0、1、2、…、n,n为整数。
上述的慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法,其特征在于:步骤五中在计算机上采用Origin软件,获得裂纹扩展长度Δai与试验时间ti之间的关系曲线上任意一点的切线斜率k=dΔai/dti的具体过程为:首先,打开分析analysis菜单,选用微积分calculus,并选其中的微分Derivative;然后,在微分Derivative菜单下方的文件列表中,找到Derivative文件,打开后查看微分结果,即得到了裂纹扩展长度Δai与试验时间ti之间的关系曲线上任意一点的切线斜率k=dΔai/dti。
本发明与现有技术相比具有以下优点:
1、本发明方法简单,实现方便。
2、本发明针对单一采用电位差法,包括交流电位差法(ACPD)和直流电位差法(DCPD),实时动态观测核电结构材料环境致裂裂纹扩展速率可能会受到的干扰和精度不足等问题,利用弹塑性条件下,标准断裂力学试样加载点位移、加载力以及裂纹长度可以通过理论分析或有限元分析获得的特点,建立起的一种理论分析、数值分析与试验观察相融合的高温水环境下核电结构材料环境致裂的慢应变速率拉伸试验裂纹扩展速率中实时动态观测和分析方法,能够有效地减少环境影响因素带来的测试数据失真,提高获取数据的精度。
3、在慢应变速率试验中,加载点位移ΔL是一个基本的监测数据,如果能够获得裂纹增长长度Δai与加载点位移ΔL的关系,就能够知道在测试过程裂纹增长的情况;而在实际过程中两个原因会促使加载点位移增加,一是由试样弹塑性变形引起的加载点位移ΔL1,二是由裂纹扩展Δai引起的加载点位移ΔL2,即ΔL=ΔL1+ΔL2,本发明有效地区分出了试样弹塑性变形和裂纹增长长度Δai对加载点位移ΔL的影响,获得了加载点位移ΔL与裂纹扩展长度Δai之间的关系,最终获得了裂纹的扩展速率dΔai/dti,是应力腐蚀开裂慢应变速率拉伸试验物理实测方式的一致有效补充。
综上所述,本发明方法简单,实现方便,能够获得裂纹的扩展速率dΔai/dti,且能够有效地减少环境影响因素带来的测试数据失真,提高获取数据的精度。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明的方法流程图。
图2为加载力P与加载点位移ΔL之间的实验数据关系曲线图。
图3为在裂纹扩展长度为Δai条件下的加载力P与加载点位移ΔL之间的理论分析数据关系曲线图。
图4为整合后的加载力P与加载点位移ΔL之间的关系曲线图。
图5为裂纹扩展长度Δai与试验时间ti之间的关系曲线图。
图6为裂纹扩展速率dΔai/dti与试验时间ti之间的关系曲线图。
具体实施方式
如图1所示,本发明所述的慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法,包括以下步骤:
步骤一、紧凑拉伸试样预处理:根据ASME标准选取一个用于慢应变速率拉伸试验的标准紧凑拉伸试样,并采用线切割和预制疲劳裂纹的方法在标准紧凑拉伸试样上设置用于初期引导裂纹扩展的预裂裂纹,所述预裂裂纹的长度为a;
本实施例中,所述a的取值为10mm。
步骤二、采用慢拉伸试验机进行慢应变速率拉伸试验,并绘制出加载力P与加载点位移ΔL之间的实验数据关系曲线图,其具体过程为:
步骤201、在慢拉伸试验机上设定慢应变速率拉伸试验的试验参数,包括加载点位移的变化速度S和试验总时间T;
本实施例中,S的取值为0.005mm/h,T的取值为900h;
步骤202、慢拉伸试验机开始对标准紧凑拉伸试样进行加载,且在加载过程中,慢拉伸试验机上的力传感器探头对加载力P进行实时检测并将所检测到的信号实时传输给与所述慢拉伸试验机连接的计算机,所述计算机采用Origin软件,以加载点位移ΔL为横坐标,加载力P为纵坐标,绘制出加载力P与加载点位移ΔL之间的实验数据关系曲线图并显示;同时,测量记录裂纹实时的扩展长度A;其中,ΔL=St且ΔL≦L,t为试验时间且t≦T,L为在试验总时间T内加载点位移总变化长度;
本实施例中得到的加载力P与加载点位移ΔL之间的实验数据关系曲线图如图2所示。由于以上获得加载点位移ΔL和加载力P数据的过程中,裂纹扩展长度Δai是不断变化的,因此得到的是在裂纹动态扩展条件下的加载力P与加载点位移ΔL之间的关系,还不能够区分得到加载点位移ΔL和影响因素加载力P、裂纹扩展长度Δai之间唯一确定的关系。
实质上,加载点位移ΔL是由两部分组成的,一是由试样弹塑性变形引起的加载点位移ΔL1,二是由裂纹扩展长度Δai时引起的加载点位移ΔL2即ΔL=ΔL1+ΔL2;利用弹塑性有限元分析方法,将裂纹扩展长度Δai作为限定的定值条件,分析计算在此条件限制下加载力P与加载点位移ΔL之间的关系,此时得到的关系是在裂纹扩展长度为Δai条件下的加载力P与加载点位移分量ΔL2的理论分析数据关系;
步骤三、对步骤二中进行的慢应变速率拉伸试验进行模拟,并绘制出加载力P与加载点位移ΔL之间的理论分析数据关系曲线图:在计算机上采用有限元分析软件对步骤二中进行的慢应变速率拉伸试验进行模拟,并分析慢应变速率拉伸试验中裂纹扩展长度递增变化的过程中,加载力P与加载点位移ΔL的关系,且在裂纹扩展长度分别为Δa0、Δa1、Δa2、…、Δan时,所述计算机采用Origin软件,以加载点位移ΔL为横坐标,加载力P为纵坐标,绘制出加载力P与加载点位移ΔL之间的理论分析数据关系曲线图并显示在同一幅图中;其中,Δan不大于步骤202中测量记录的裂纹实时的扩展长度A的最大值;
本实施例中,步骤三中裂纹扩展长度Δai=iΔx,i的取值为0、1、2、…、n,n为整数。其中,Δx为裂纹单次扩展长度,Δx的取值为0.375mm,n的取值为8,即在裂纹扩展长度分别为0mm、0.375mm、0.75mm、1.125mm、1.5mm、1.875mm、2.25mm、2.625mm、3mm时,裂纹长度分别为10mm、10.375mm、10.75mm、11.125mm、11.5mm、11.875mm、12.25mm、12.625mm、13mm时,所述计算机采用Origin软件,以加载点位移ΔL为横坐标,加载力P为纵坐标,绘制出加载力P与加载点位移ΔL之间的理论分析数据关系曲线图并显示在同一幅图中;
本实施例中,步骤三中所述有限元分析软件为ABAQUS软件。步骤三中在计算机上采用有限元分析软件对步骤二中进行的慢应变速率拉伸试验进行模拟,并分析慢应变速率拉伸试验中裂纹扩展长度递增变化的过程中,加载力P与加载点位移ΔL的关系的具体过程为:
步骤301、在ABAQUS软件中,建立带有长度为a的预裂裂纹的标准紧凑拉伸试样的三维模型;
步骤302、在步骤301中建立的标准紧凑拉伸试样的三维模型中,添加长度为Δa1的裂纹;即添加长度为0.375mm的裂纹,使得裂纹总长度为10.375mm;
步骤303、在ABAQUS软件中,建立慢拉伸试验机对标准紧凑拉伸试样进行加载时所用的两个加载销的刚体模型;
步骤304、赋予标准紧凑拉伸试样的三维模型材料参数,所述材料参数包括材料的类型、材料的杨氏模量、材料的泊松比、材料的屈服强度、材料的硬化指数和材料的偏移指数;本实施例中,材料的类型为镍基合金、材料的杨氏模量为210GPa、材料的泊松比为0.3、材料的屈服强度为436MPa、材料的硬化指数为6.495,材料的偏移指数为1;
步骤305、将标准紧凑拉伸试样和两个加载销生成装配体,并移动两个加载销的位置,使得两个加载销完全对称地装配在所述标准紧凑拉伸试样上,且使得两个加载销伸出到所述标准紧凑拉伸试样左右两侧的长度相等;
步骤306、创建第一分析步和第二分析步两个分析步,与ABAQUS软件自带的初始分析步构成三个分析步;其中,初始分析步用于建立两个加载销和标准紧凑拉伸试样之间的初始接触关系;第一分析步用于施加一个载荷使初始接触关系平稳建立,便于计算的收敛;第二分析步用于加载;
步骤307、在初始分析步中,设置两个加载销和标准紧凑拉伸试样之间的初始接触关系为硬接触关系,主接触表面选择为销的外表面,从接触表面选择为孔的内表面;并设置两个加载销除拉伸方向外的其余两个直线自由度和三个旋转自由度均为0,对两个加载销除拉伸方向的直线自由度外所有的自由度进行约束;其中,硬接触关系的特点是:接触面能传递的接触压力大小不受限制,当接触压力变为零或负值时,两个面分离,不存在接触约束;
步骤308、在第一分析步中,对两个加载销沿拉伸方向分别施加一个大小为0.02mm~0.08mm的位移载荷;在第二分析步中,对两个加载销沿拉伸方向分别施加线性变化的大小为2mm~8mm位移载荷;优选地,在第一分析步中,对两个加载销沿拉伸方向分别施加一个大小为0.05mm的位移;在第二分析步中,对两个加载销沿拉伸方向分别施加线性变化的大小为5mm位移载荷;其中,施加的线性变化的位移载荷的大小需要保证慢应变速率拉伸试验能够进行至少900h。
步骤309、对由标准紧凑拉伸试样和两个加载销生成的装配体进行网格划分,在裂纹尖端划分密致网格,其余位置划分稀疏网格,并将网格的单元类型选择为三维应力六面体一次单元;本实施例中,在裂纹尖端划分的密致网格,1立方米有106个网格;其余位置划分的稀疏网格,1立方米有103个网格;
步骤310、提交计算,ABAQUS软件按照步骤301~309设置的参数进行分析计算,模拟慢应变速率拉伸试验中的裂纹扩展过程;
步骤311、裂纹扩展过程中,进行数据提取,先提取加载销的加载点位移ΔL随时间变化的数据,再提取加载力P随时间变化的数据,得到了加载力P与加载点位移ΔL的关系。
本实施例中得到的加载力P与加载点位移ΔL之间的理论分析数据关系曲线图如图3所示,图3中,曲线L0、L1、L2、L3、L4、L5、L6、L7和L8分别对应裂纹扩展长度分别为0mm、0.375mm、0.75mm、1.125mm、1.5mm、1.875mm、2.25mm、2.625mm、3mm时,即裂纹长度分别为10mm、10.375mm、10.75mm、11.125mm、11.5mm、11.875mm、12.25mm、12.625mm、13mm时,加载力P与加载点位移ΔL之间的理论分析数据关系曲线。这里模拟得到的是静态裂纹条件的情况下,即限制了裂纹扩展过程中裂纹扩展长度的这一影响因素的情况下,得到的单一影响因素加载力P与加载点位移ΔL的影响关系,即确定出了在某一固定裂纹扩展长度条件下的加载力P和加载点位移ΔL的关系,这样可以消除裂纹扩展长度对于加载点位移ΔL的影响。
步骤四、获取裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,所对应的试验时间t0、t1、t2、…、tn:在计算机上采用Origin软件将步骤202中绘制的加载力P与加载点位移ΔL之间的实验数据关系曲线图和步骤三中绘制的在裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,加载力P与加载点位移ΔL之间的理论分析数据关系曲线图整合到同一幅图中,得到整合后的加载力P与加载点位移ΔL之间的关系曲线图,且在整合后的加载力P与加载点位移ΔL之间的关系曲线图中,提取出裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,加载力P与加载点位移ΔL之间的实验数据关系曲线与n条加载力P与加载点位移ΔL之间的理论分析数据关系曲线的交点Qi的坐标数据,得到了裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,所对应的加载点位移ΔL0、ΔL1、ΔL2、…、ΔLn,由进而计算得到了裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,所对应的试验时间t0、t1、t2、…、tn;其中,i的取值为0、1、2、…、n;
本实施例中,为了合理区分数据类型,将步骤202中绘制的加载力P与加载点位移ΔL之间的实验数据关系曲线设置为实线并加粗,将步骤三中绘制的加载力P与加载点位移ΔL之间的理论分析数据关系曲线设置为虚线,得到的整合后的加载力P与加载点位移ΔL之间的关系曲线图如图4所示,图4中,加载力P与加载点位移ΔL之间的实验数据关系曲线与裂纹扩展长度为0mm时加载力P与加载点位移ΔL之间的理论分析数据关系曲线L0没有交点,即交点Q0不存在;加载力P与加载点位移ΔL之间的实验数据关系曲线与裂纹扩展长度Δa1为0.375mm时加载力P与加载点位移ΔL之间的理论分析数据关系曲线L1的交点为Q1,从Q1的坐标可以得到,裂纹扩展长度Δa1为0.375mm时,所对应的加载点位移ΔL为0.826mm;加载力P与加载点位移ΔL之间的实验数据关系曲线与裂纹扩展长度Δa2为0.75mm时加载力P与加载点位移ΔL之间的理论分析数据关系曲线L2的交点为Q2,从Q2的坐标可以得到裂纹扩展长度Δa2为0.75mm时所对应的加载点位移;加载力P与加载点位移ΔL之间的实验数据关系曲线与裂纹扩展长度Δa3为11.125mm时加载力P与加载点位移ΔL之间的理论分析数据关系曲线L3的交点为Q3,从Q3的坐标可以得到裂纹扩展长度Δa3为1.125mm时所对应的加载点位移;加载力P与加载点位移ΔL之间的实验数据关系曲线与裂纹扩展长度Δa4为1.5mm时加载力P与加载点位移ΔL之间的理论分析数据关系曲线L4的交点为Q4,从Q4的坐标可以得到裂纹扩展长度Δa4为1.5mm时所对应的加载点位移;加载力P与加载点位移ΔL之间的实验数据关系曲线与裂纹扩展长度Δa5为1.875mm时加载力P与加载点位移ΔL之间的理论分析数据关系曲线L5的交点为Q5,从Q5的坐标可以得到裂纹扩展长度Δa5为1.875mm时所对应的加载点位移;加载力P与加载点位移ΔL之间的实验数据关系曲线与裂纹扩展长度Δa6为2.25mm时加载力P与加载点位移ΔL之间的理论分析数据关系曲线L6的交点为Q6,从Q6的坐标可以得到裂纹扩展长度Δa6为2.25mm时所对应的加载点位移;加载力P与加载点位移ΔL之间的实验数据关系曲线与裂纹扩展长度Δa7为2.625mm时加载力P与加载点位移ΔL之间的理论分析数据关系曲线L7的交点为Q7,从Q7的坐标可以得到裂纹扩展长度Δa7为2.625mm时所对应的加载点位移;加载力P与加载点位移ΔL之间的实验数据关系曲线与裂纹扩展长度Δa8为3mm时加载力P与加载点位移ΔL之间的理论分析数据关系曲线L8的交点为Q8,从Q8的坐标可以得到裂纹扩展长度Δa8为3mm时所对应的加载点位移;由于已知S=0.005mm/h,因此由公式能够得到裂纹扩展长度Δai分别为0mm、0.375mm、0.75mm、1.125mm、1.5mm、1.875mm、2.25mm、2.625mm和3mm时,所对应的试验时间为t0、t1、t2、t3、t4、t5、t6、t7、t8。
具体实施时,为了能够得到更准确的裂纹扩展长度Δai和试验时间ti之间的关系,可以在步骤三中减小Δx的取值,即增大n。
步骤五、获得裂纹扩展速率:在计算机上采用Origin软件,以试验时间ti为横坐标,裂纹扩展长度Δai为纵坐标,绘制出步骤四中得到的实验时间t0、t1、t2、…、tn和裂纹扩展长度Δa0、Δa1、Δa2、…、Δan对应的各个点,并拟合得到裂纹扩展长度Δai与试验时间ti之间的关系曲线,裂纹扩展长度Δai与试验时间ti之间的关系曲线上任意一点的切线斜率k=dΔai/dti即为当时的裂纹扩展速率,其中,i的取值为0、1、2、…、n。本实施例中,裂纹扩展长度Δai与试验时间ti之间的关系曲线图如图5所示。
本实施例中,步骤五中在计算机上采用Origin软件,获得裂纹扩展长度Δai与试验时间ti之间的关系曲线上任意一点的切线斜率k=dΔai/dti的具体过程为:首先,打开分析analysis菜单,选用微积分calculus,并选其中的微分Derivative;然后,在微分Derivative菜单下方的文件列表中,找到Derivative文件,打开后查看微分结果,即得到了裂纹扩展长度Δai与试验时间ti之间的关系曲线上任意一点的切线斜率k=dΔai/dti。本实施例中,得到的裂纹扩展速率dΔai/dti与试验时间ti之间的关系曲线图如图6所示。
以上所述,仅是本发明的较佳实施例,并非对本发明作任何限制,凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化,均仍属于本发明技术方案的保护范围内。
Claims (5)
1.一种慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法,其特征在于该方法包括以下步骤:
步骤一、紧凑拉伸试样预处理:根据ASME标准选取一个用于慢应变速率拉伸试验的标准紧凑拉伸试样,并采用线切割和预制疲劳裂纹的方法在标准紧凑拉伸试样上设置用于初期引导裂纹扩展的预裂裂纹,所述预裂裂纹的长度为a;
步骤二、采用慢拉伸试验机进行慢应变速率拉伸试验,并绘制出加载力P与加载点位移ΔL之间的实验数据关系曲线图,其具体过程为:
步骤201、在慢拉伸试验机上设定慢应变速率拉伸试验的试验参数,包括加载点位移的变化速度S和试验总时间T;
步骤202、慢拉伸试验机开始对标准紧凑拉伸试样进行加载,且在加载过程中,慢拉伸试验机上的力传感器探头对加载力P进行实时检测并将所检测到的信号实时传输给与所述慢拉伸试验机连接的计算机,所述计算机采用Origin软件,以加载点位移ΔL为横坐标,加载力P为纵坐标,绘制出加载力P与加载点位移ΔL之间的实验数据关系曲线图并显示;同时,测量记录裂纹实时的扩展长度A;其中,ΔL=St且ΔL≦L,t为试验时间且t≦T,L为在试验总时间T内加载点位移总变化长度;
步骤三、对步骤二中进行的慢应变速率拉伸试验进行模拟,并绘制出加载力P与加载点位移ΔL之间的理论分析数据关系曲线图:在计算机上采用有限元分析软件对步骤二中进行的慢应变速率拉伸试验进行模拟,并分析慢应变速率拉伸试验中裂纹扩展长度递增变化的过程中,加载力P与加载点位移ΔL的关系,且在裂纹扩展长度分别为Δa0、Δa1、Δa2、…、Δan时,所述计算机采用Origin软件,以加载点位移ΔL为横坐标,加载力P为纵坐标,绘制出加载力P与加载点位移ΔL之间的理论分析数据关系曲线图并显示在同一幅图中;其中,Δan不大于步骤202中测量记录的裂纹实时的扩展长度A的最大值;
步骤四、获取裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,所对应的试验时间t0、t1、t2、…、tn:在计算机上采用Origin软件将步骤202中绘制的加载力P与加载点位移ΔL之间的实验数据关系曲线图和步骤三中绘制的在裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,加载力P与加载点位移ΔL之间的理论分析数据关系曲线图整合到同一幅图中,得到整合后的加载力P与加载点位移ΔL之间的关系曲线图,且在整合后的加载力P与加载点位移ΔL之间的关系曲线图中,提取出裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,加载力P与加载点位移ΔL之间的实验数据关系曲线与n条加载力P与加载点位移ΔL之间的理论分析数据关系曲线的交点Qi的坐标数据,得到了裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,所对应的加载点位移ΔL0、ΔL1、ΔL2、…、ΔLn,由进而计算得到了裂纹扩展长度为Δa0、Δa1、Δa2、…、Δan时,所对应的试验时间t0、t1、t2、…、tn;其中,i的取值为0、1、2、…、n;
步骤五、获得裂纹扩展速率:在计算机上采用Origin软件,以试验时间ti为横坐标,裂纹扩展长度Δai为纵坐标,绘制出步骤四中得到的实验时间t0、t1、t2、…、tn和裂纹扩展长度Δa0、Δa1、Δa2、…、Δan对应的各个点,并拟合得到裂纹扩展长度Δai与试验时间ti之间的关系曲线,裂纹扩展长度Δai与试验时间ti之间的关系曲线上任意一点的切线斜率k=dΔai/dti即为当时的裂纹扩展速率,其中,i的取值为0、1、2、…、n。
2.按照权利要求1所述的慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法,其特征在于:步骤三中所述有限元分析软件为ABAQUS软件。
3.按照权利要求2所述的慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法,其特征在于:步骤三中在计算机上采用有限元分析软件对步骤二中进行的慢应变速率拉伸试验进行模拟,并分析慢应变速率拉伸试验中裂纹扩展长度递增变化的过程中,加载力P与加载点位移ΔL的关系的具体过程为:
步骤301、在ABAQUS软件中,建立带有长度为a的预裂裂纹的标准紧凑拉伸试样的三维模型;
步骤302、在步骤301中建立的标准紧凑拉伸试样的三维模型中,添加长度为Δa1的裂纹;
步骤303、在ABAQUS软件中,建立慢拉伸试验机对标准紧凑拉伸试样进行加载时所用的两个加载销的刚体模型;
步骤304、赋予标准紧凑拉伸试样的三维模型材料参数,所述材料参数包括材料的类型、材料的杨氏模量、材料的泊松比、材料的屈服强度、材料的硬化指数和材料的偏移指数;
步骤305、将标准紧凑拉伸试样和两个加载销生成装配体,并移动两个加载销的位置,使得两个加载销完全对称地装配在所述标准紧凑拉伸试样上,且使得两个加载销伸出到所述标准紧凑拉伸试样左右两侧的长度相等;
步骤306、创建第一分析步和第二分析步两个分析步,与ABAQUS软件自带的初始分析步构成三个分析步;其中,初始分析步用于建立两个加载销和标准紧凑拉伸试样之间的初始接触关系;第一分析步用于施加一个载荷使初始接触关系平稳建立,便于计算的收敛;第二分析步用于加载;
步骤307、在初始分析步中,设置两个加载销和标准紧凑拉伸试样之间的初始接触关系为硬接触关系,主接触表面选择为销的外表面,从接触表面选择为孔的内表面;并设置两个加载销除拉伸方向外的其余两个直线自由度和三个旋转自由度均为0,对两个加载销除拉伸方向的直线自由度外所有的自由度进行约束;
步骤308、在第一分析步中,对两个加载销沿拉伸方向分别施加一个大小为0.02mm~0.08mm的位移载荷;在第二分析步中,对两个加载销沿拉伸方向分别施加线性变化的大小为2mm~8mm位移载荷;
步骤309、对由标准紧凑拉伸试样和两个加载销生成的装配体进行网格划分,在裂纹尖端划分密致网格,其余位置划分稀疏网格,并将网格的单元类型选择为三维应力六面体一次单元;
步骤310、提交计算,ABAQUS软件按照步骤301~309设置的参数进行分析计算,模拟慢应变速率拉伸试验中的裂纹扩展过程;
步骤311、裂纹扩展过程中,进行数据提取,先提取加载销的加载点位移ΔL随时间变化的数据,再提取加载力P随时间变化的数据,得到了加载力P与加载点位移ΔL的关系。
4.按照权利要求1所述的慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法,其特征在于:步骤三中裂纹扩展长度Δai=iΔx,其中,Δx为裂纹单次扩展长度,i的取值为0、1、2、…、n,n为正整数。
5.按照权利要求1所述的慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法,其特征在于:步骤五中在计算机上采用Origin软件,获得裂纹扩展长度Δai与试验时间ti之间的关系曲线上任意一点的切线斜率k=dΔai/dti的具体过程为:首先,打开分析analysis菜单,选用微积分calculus,并选其中的微分Derivative;然后,在微分Derivative菜单下方的文件列表中,找到Derivative文件,打开后查看微分结果,即得到了裂纹扩展长度Δai与试验时间ti之间的关系曲线上任意一点的切线斜率k=dΔai/dti。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410012971.8A CN103698188B (zh) | 2014-01-11 | 2014-01-11 | 慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410012971.8A CN103698188B (zh) | 2014-01-11 | 2014-01-11 | 慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103698188A CN103698188A (zh) | 2014-04-02 |
CN103698188B true CN103698188B (zh) | 2016-04-13 |
Family
ID=50359792
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410012971.8A Expired - Fee Related CN103698188B (zh) | 2014-01-11 | 2014-01-11 | 慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103698188B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105092457B (zh) * | 2015-07-09 | 2018-03-13 | 中国石油天然气股份有限公司 | 交变载荷和腐蚀介质共同作用下的注采管柱腐蚀评价方法 |
CN105571945B (zh) * | 2015-12-18 | 2018-03-27 | 中国科学院地质与地球物理研究所 | 岩石原位微拉伸试样及测试方法 |
US10775286B2 (en) * | 2016-01-28 | 2020-09-15 | The Boeing Company | Method for test data reduction of composite intralaminar failure mode |
CN109724872A (zh) * | 2018-12-11 | 2019-05-07 | 江苏亚星锚链股份有限公司 | 一种r6级系泊链模拟工况下抗环境脆化性能评价方法 |
CN110411819B (zh) * | 2019-07-12 | 2022-04-08 | 中国船舶重工集团公司第七二五研究所 | 一种用光滑圆棒拉伸试样测试应力腐蚀门槛值kiscc的方法 |
CN110455627A (zh) * | 2019-08-26 | 2019-11-15 | 中国特种设备检测研究院 | 基于恒位移加载的材料与高压氢气相容性评价方法及系统 |
CN110763575B (zh) * | 2019-11-07 | 2022-01-25 | 建龙北满特殊钢有限责任公司 | 一种快速检测残余延伸强度的方法 |
CN111044367B (zh) * | 2019-12-31 | 2020-12-18 | 成都理工大学 | 一种基于三轴应力-应变曲线的岩石裂缝亚临界扩展速率实验测试方法 |
CN112557229B (zh) * | 2020-12-08 | 2024-03-19 | 中国石油化工股份有限公司 | 一种评价金属材料慢拉伸应力腐蚀敏感性的方法 |
CN112801981B (zh) * | 2021-01-28 | 2022-06-17 | 中国科学院武汉岩土力学研究所 | 确定混合压剪裂纹尖端断裂过程区传播速度的方法及设备 |
CN113283137B (zh) * | 2021-05-24 | 2023-03-03 | 中国兵器工业第五九研究所 | 一种有机防护涂层体系下结构件腐蚀效应控制仿真方法 |
CN114414364A (zh) * | 2021-11-30 | 2022-04-29 | 温州电力建设有限公司 | 基于2d激光的导地线单丝力学强度测试系统及方法 |
CN116399694A (zh) * | 2023-03-10 | 2023-07-07 | 北京建筑大学 | 一种直接拉伸脆性固体材料长期寿命评价方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0545269A (ja) * | 1991-08-12 | 1993-02-23 | Japan Steel Works Ltd:The | 引張試験における降伏点の検出方法および引張試験装置 |
CN100510696C (zh) * | 2006-12-01 | 2009-07-08 | 中国科学院海洋研究所 | 动载条件下研究大气中应力腐蚀敏感性的装置及方法 |
CN202757839U (zh) * | 2012-08-16 | 2013-02-27 | 中国石油化工集团公司 | 一种慢应变速率拉伸设备 |
CN102928625B (zh) * | 2012-10-22 | 2015-07-15 | 山东大学 | 用于板状载荷试样腐蚀研究的扫描电化学显微镜电解池及其应用 |
CN203216827U (zh) * | 2013-05-15 | 2013-09-25 | 国家电网公司 | 一种应力腐蚀慢应变速率拉伸实验用极化实验装置 |
-
2014
- 2014-01-11 CN CN201410012971.8A patent/CN103698188B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN103698188A (zh) | 2014-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103698188B (zh) | 慢应变速率拉伸应力腐蚀裂纹扩展速率测量方法 | |
Li et al. | A statistical damage constitutive model for softening behavior of rocks | |
CN103344705B (zh) | 一种应用声发射能量值测定岩石脆性指数的方法 | |
Wang et al. | Numerical investigation on the creep crack-tip constraint induced by loading configuration of specimens | |
Yang et al. | A damage identification approach for plate structures based on frequency measurements | |
CN102865952B (zh) | 一种混凝土工作应力的无损测试方法 | |
CN104990793B (zh) | 火灾后木结构建筑的承载能力评估方法 | |
Li et al. | A 3D grain-based creep model (3D-GBCM) for simulating long-term mechanical characteristic of rock salt | |
He et al. | Characterization of 3-D creep constraint and creep crack growth rate in test specimens in ASTM-E1457 standard | |
Mu et al. | Three‐dimensional analyses of in‐plane and out‐of‐plane crack‐tip constraint characterization for fracture specimens | |
CN103471940A (zh) | 面内及面外统一拘束与材料延性断裂韧性关系的确定方法 | |
CN106596255A (zh) | 一种基于光纤光谱图像分析的孔边裂纹诊断方法 | |
CN105136596B (zh) | 具有恒位移载荷的裂尖应力腐蚀开裂状态测试系统及方法 | |
Van Paepegem | Fatigue damage modelling of composite materials with the phenomenological residual stiffness approach | |
Wu et al. | Crack diagnosis method for a cantilevered beam structure based on modal parameters | |
CN105352433A (zh) | 一种船体典型焊接结构表面裂纹深度及形状测量装置及其方法 | |
Lai et al. | A new method to measure crack extension in nuclear graphite based on digital image correlation | |
Zhang et al. | Dynamic shear fracture behaviour of granite under axial static pre-force by 3D high-speed digital image correlation | |
He et al. | Effects of side‐groove depth on creep crack‐tip constraint and creep crack growth rate in C (T) specimens | |
Lindström et al. | Non-linear fracture mechanics in LS-DYNA and LS-PrePost | |
Nor et al. | An overview of fracture mechanics with ANSYS | |
Yang et al. | Computation and analysis of high rocky slope safety in a water conservancy project | |
Ayatollahi et al. | Computation of V-notch shape factors in four-point bend specimen for fracture tests on brittle materials | |
Garcia et al. | Back-face strain compliance relation for SEN (B) specimens for wide range in crack lengths | |
Maier et al. | Inverse analysis problems in structural engineering of concrete dams |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160413 Termination date: 20170111 |
|
CF01 | Termination of patent right due to non-payment of annual fee |