盲孔法中释放系数的数值计算方法

JournalofMechanicalStrength

2008,30(1):033-036

盲孔法中释放系数的数值计算方法

*

NUMERICALMETHODFORCALCULATIONOFRELEASE

COEFFICIENTSINBLINDHOLEMETHOD

**刘一华

贝晖 詹春晓 李 昊 贺斌

(合肥工业大学土木建筑工程学院,合肥230009)LIUYiHua HeYunHui ZHANChunXiao LIHao

(SchoolofCivilEngineering,HefeiUniversityofTechnology,Hefei230009,China)

摘要 盲孔法测量残余应力时,钻孔会在盲孔附近产生加工硬化。由于硬度与弹性模量的平方有关,在进行释放系数A、B的数值计算时,通过增大盲孔周围材料的弹性模量来反映该加工硬化。应用有限元软件ANSYS,对304不锈钢的释放系数A、B进行三维有限元数值计算,得到释放系数随加工硬化程度增大而增大的变化规律,同时给出考虑加工硬化后释放系数随盲孔深度的变化规律。通过与实验测定结果的比较,发现在对304不锈钢进行释放系数的有限元数值计算时,当取盲孔附近的加工硬化层为0.2mm厚、弹性模量增大37%时,有限元的数值计算值与实验测定值最接近,释放系数A、B的数值计算误差分别为-1.8%和2.4%。

关键词 盲孔法 释放系数 有限元法 数值计算 加工硬化 弹性模量中图分类号 O348.2

Abstract Work-hardeningneartheblindholeiscausedbyholedrillingwhenresidualstressesaremeasuredbytheblindholemethod.Ashardeningisrelatedtothesquareofelasticmodulus,thework-hardeningisconsideredbyincreasingtheelasticmodulusofthematerialneartheholewhenreleasecoefficientsAandBarenumericallyanalyzed.ThereleasecoefficientsAandBofstainlesssteel304arenumericallycalculatedbythethreedimensionalfiniteelementmethodprogramANSYS.Therelationshipthatthereleasecoeffi-cientsAandBareincreasingwiththeincreasingofdegreeofwork-hardeningisobtained.Meanwhile,paringthecalculatingresultstotheexperimentalones,itcanbefoundthatthecalculatingvaluesisproximaltotheexperimentaloneswhilethethicknessofwork-hardeninglayeris0.2mmandtheelasticmodulusis37%higherthantheoriginalone.Onthiscondition,thecalculatingerrorsofthereleasecoefficientsAandBareonly-1.8%and2.4%,respectively.

Keywords Blindholemethod;Releasecoefficient;Finiteelementmethod;Numericalcalculation;Work-hardening;Elasticmodulus

Correspondingauthor:LIUYiHua,E-mail:liuyihua@,Tel:+86-551-2902547;Fax:+86-551-2902066TheprojectsupportedbytheNaturalScienceKeyFoundationofEducationDepartmentofAnhuiProvince(No.2003kj045zd),andtheResearchProgramforYoungTeachersofHigherEducationofAnhuiProvince(No.2006jgl005)andtheResearchFoundationofHefeiUniversityofTechnology(No.060803F).

Manuscriptreceived20060406,inrevisedform20060717.

1 引言

盲孔法是由德国学者Mathar

[2]

[1]

中得到广范应用。

于1934年提出的,

在用盲孔法测量残余应力时,确定释放系数A、

B致关重要,它直接影响残余应力的测量精度。释放系数A、B通常采用实验测定

[7-8]

[5-6]

后经Soete和Vancombrugge等学者发展起来的一种测量残余应力的半无损方法,它具有对工件损伤小、不影响工件的性能、简单易行、测量精度高等特点。为此,

1981年美国材料试验协会将其制定为ASTM标准E837)81,1992年我国船舶工业总公司也将其制定为我国船舶行业标准CB3395)92,使盲孔法在工程

[4]

[3]

,因实验测定颇费人

力物力,有人提出基于通孔理论的近似计算公式;为提高确定释放系数的精度和效率,近年有

[9]632-636

人提出数值计算方法。由于在用实验方法测定释放系数时需要钻孔,钻孔后会在孔边附近产生加工硬化

[10]139

,从而降低应变的释放。为了

真实模拟释

*20060406收到初稿,20060717收到修改稿。安徽省教育厅自然科学重点科研计划项目(2003kj045zd)、安徽省高校青年教师科研资助计划项目

(2006jgl005)和合肥工业大学科研发展基金(060803F)资助。

**刘一华,男,1959年8月生,湖南桃源人,汉族。合肥工业大学土木建筑工程学院教授,博士生导师,博士。研究方向为残余应力和断裂力学

034机 械 强 度2008年

放系数的测定实验,本文提出一种考虑钻孔后在孔边附近产生的加工硬化、应用有限元法确定释放系数的数值计算方法。

系数A和B,为监控实验加载,在试样的两侧粘贴2个单向应变计(应变计3、4)。残余应变计的敏感栅尺寸为1mm@1mm。用钻孔前由2个残余应变计测得的应变值可得到材料的弹性模量E和泊松比M。钻孔时先用Á1的钻头快速钻孔,再用Á1.5的钻头慢速扩孔,钻孔深度h=2mm。弹性常数和释放系数的实验测定结果如表1所示。

表1 弹性常数E、M和释放系数A、B的实验测定值Tab.1 ExperimentalvaluesofelasticconstantsE、M

andreleasecoefficientsA、B

测定点Measuredpoint

1

EPGPa210213212

M0.2750.2800.278

A@106PMPa-1-0.3649-0.3937-0.3793

B@106PMPa-1-0.8353-0.8257-0.8305

2 盲孔法的基本原理

若构件中一点O的主残余应力分别为R1和R2,在O点粘贴一残余应变计,如图1所示,当在O点钻一直径为d的盲孔后,由敏感栅R1、R2和R3所测得的释放应变分别为E1、E2和E3,则O点的主残余应力及其方向分别为

R1

E1+E3=4AºR2

1

4B

(E1-E3)+(E1+E3-2E2)

(1)(2)

2平均值Averagevalue

E1+E3-2E2

tan(2H)=1-E3

式中,A和B称为释放系数,由实验测定。

释放系数A、B通常在均匀单向拉伸应力场R1=R、R2=0中测定,在此情况下释放系数A、B分别为E1+E3E1-E

3

B=(3)2R2R

式中,E1和E3分别为沿R1

方向和垂直于R1方向的线

A=

应变。

4 释放系数的数值计算

4.1 计算模型及其网格划分

为了模拟上述释放系数的测定实验,取有限元分

析模型如图3所示,在模型的中间有一盲孔,孔的直径d=1.5m,孔深h=2mm,应用ANSYS有限元软件进行三维数值模拟。由于对称性,取1P4进行计算,有限元网格划分如图4所示。为了模拟真实情况,在粘贴应变计的位置处,按照敏感栅的形状和大小划分网格(图4b中的A、B矩形区域),矩形区域尺寸为1mm@1mm,其靠近孔口一侧到孔中心的距离为1.6mm。由于钻孔后在孔边附近会产生加工硬化,在钢材上钻孔及扩孔时的平均加工硬化深度为0.18mm~0.2mm,为此,在盲孔边缘取0.2mm厚作为加工硬化层,划分环形单元(图4b中的C区域)。

[10]139

图1 残余应变计的布置

Fig.1 Arrangementofresidualstraingauge

3 释放系数的实验测定结果

用304不锈钢板加工成图2所示的试样,在试样的中间粘贴2个残余应变计(应变计1、2),用于测定释放

图3 有限元分析模型Fig.3 Finiteelementanalysismodel

4.2 释放系数随硬化程度的变化情况

[11]

研究表明,硬度H与弹性模量E有如下关系

H=KE

2

(4)

式中,K为常数。

根据式(4),加工硬化程度N(即加工后的表面显

图2 测定释放系数的试样

Fig. formeasuringrelease微硬度与原始显微硬度之比@100%)

[10]139

和弹性模量

增大程度M(即加工后的表面弹性模量与原始弹性模

第30卷第1期刘一华等:盲孔法中释放系数的数值计算方法035

图4 有限元网格划分Fig.4 FEMmesh

图5 不同弹性模量增大程度下释放系数的计算误差Fig.5 Calculatingerrorsofreleasecoefficientsindifferent

increasingdegreeofelasticmodulus

量之比@100%)之间有下列关系

N=M

2

(5)

为此,在取加工硬化层厚度为0.2mm的情况下,通过增大硬化层的弹性模量,计算出不同弹性模量增大程度M(即不同加工硬化程度N)下释放系数A、B的数值计算值,并与表1所示的实验测定值进行比较,如表2所示,图5给出不同弹性模量增大程度M下释放系数的数值计算值与实验测定值的误差关系。由图5可见,数值计算误差与M成线性关系,用最小二乘法进行线性拟合,可以得到释放系数A、B的数值计算值的计算误差DDB随M的线性变化关系为A、

DA=0.5205-0.3934MDB=0.3941-0.2704M

表2 不同加工硬化程度时释放系数的计算值及其误差Tab.2 Calculatingvaluesofreleasecoefficientsandtheir

errorsindifferentdegreeofwork-hardening

A@106PMPa-1

M(%)

N(%)

计算值

误差

B@106PMPa-1计算值

误差

(6)

图6 释放系数随孔深的变化趋势

Fig.6 Trendtoreleasecoefficientswithdepthofhole

CalculatingvalueError(%)CalculatingvalueError(%)

100120130140150160

100144169196225256

-0.4300-0.3967-0.3813-0.3668-0.3529-0.3397

13.44.60.5-3.3-7.0-10.4

-0.9370-0.8859-0.8629-0.8414-0.8211-0.8021

12.86.73.91.3-1.1-3.4

图7 表面自熔焊试件

Fig.7 Sampleforsurfacesel-ffusiblewelding

由式(6)可得,当M=137%时,释放系数的数值计算值A=-0.3710@10MPa、B=-0.8469@10MPa2.4%。

4.3 释放系数随孔深的变化情况

除了孔口表面的加工硬化对释放系数会产生影响外,孔的深度对释放系数的影响也很明显

[9]634

-6

-1

-6

-1

的误差最小,分别为DDA=-1.8%、B=

图6中还给出未考虑加工硬化时的释放系数A、B的数值计算值。由图6可见,当考虑钻孔产生的加工硬化后,释放系数随孔深的变化趋势与未考虑加工硬化时

的变化趋势基本相同,但释放系数的数值增大(绝对值减小)。

。为此,

取加工硬化层厚度为0.2mm、弹性模量增大程度M=

137%(即硬化程度N=188%)时,对不同孔深情况下

6,

5 实验检验

A

036机 械 强 度2008年

有效性,用盲孔法对图7所示的304不锈钢表面自熔焊试件焊缝表面附近沿焊缝方向的焊接残余应力Rx的分布情况进行测量,焊缝宽8mm,深2mm~3mm,采用两根试件,共布置7个测点。分别用实验测定的释放系数和上述数值计算的释放系数计算出的焊接残余应力如表3所示,由表3可见,用数值计算的释放系数计算出的残余应力与用实验测定的释放系数计算出的残余应力很接近,最大误差0.72%。

表3 残余应力Rx随坐标y的变化情况Tab.3 VaryingofresidualstressRxwithcoordinatey

RxPMPa

xPmm

yPmm

误差

实验测定A、BExperimentalA,B

2222-44-220-44-22

-11-8-5-3016

152382412479442434319

数值计算A、BNumericalA,B

153383413480443435321

Error(%)0.580.330.360.160.200.130.72

参考文献(References)

[1] [2] [3]

MatharJ.Determinationofinitialstressesbymeasuringthedeformationsarounddrilledholes[J].Trans.ASME,1934,4:249-254.

SoeteW,VancombruggeR.Anindustrialmethodforthedeterminationofresidualstresses[J].Proc.SESA,1950,8(1):17-26.

ASTMCommitteeE28onMechanicalTesting.ASTME837)81.Stan-dardtestmethodfordeterminingresidualstressesbythehole-drillingstrain-gagemethod.Philadelphia:AmericanSocietyforTestingMaterials(ASTM),1981:1-10.

[4]

中国船舶总公司.CB3395)92残余应力测试方法钻孔应变释放法[S].北京:中国船舶工业总公司,1992.1-10.

ChinashipIndustryCorporation.CB3395)92.Methodformeasuringresidualstressesbythehole-drillingstrainreleasemethod[S].Beijing:GeneralShipCompanyofChina,1992:1-10(InChinese).

[5]

袁发荣,伍尚礼.残余应力测试与计算[M].长沙:湖南大学出版社,1987:55-64.

YUANFaRong,WUShangLi.Themeasurementandcalculationofresid-ualstress[M].Changsha:HunanUniversityPress,1987:55-46(InChinese).

[6]

裴 怡,包亚峰,唐慕尧.盲孔法测定时计算公式中A、B值的研究[J].机械强度,1997,19(1):18-21.

PEIYi,BAOYaFeng,TANGMuYao.StudyontheA、Bvaluesinca-lculatingformulaofblind-holetechnique[J].JournalofMechanicalStrength,1997,19(1):18-21(InChinese).

[7]

KabiriM.Towardmoreaccurateresidua-lstressmeasurementbythehole-drillingmethod:analysisofrelieved-straincoefficients[J].ExperimentalMechanics,1986,25(1):14-21.[8]

陈惠南.盲孔法测量残余应力A、B系数计算公式讨论[J].机械强度,1989,11(2):31-36.

CHENHuiNan.AdiscussionofcalculationofcoefficientsAandBbyblind-hole-drillingmethod[J].JournalofMechanicalStrength,1989,11(2):31-36(InChinese).

[9]

侯海量,朱 锡,刘润泉.盲孔法测量焊接残余应力应变释放系数的有限元分析[J].机械强度,2003,25(6):632-636.

HOUHaiLiang,ZHUXi,LIURunQuan.FEManalysisonthestrainre-leasefactorsformeasuringweldingresidualstressesusingblindholemethod[J].JournalofMechanicalStrength,2003,25(6):632-636(InChinese).

[10] 上海市金属切削技术协会.金属切削手册[M].3版.上海:上

海科学技术出版社,2000:170-174.

ShanghaiSocietyofMetalCuttingTechnology.Handbookofmetalcutting[M].3thed.Shanghai:ShanghaiScienceandTechnologyPress,2000.170-174(InChinese).

[11] BaoYW,WangW,ZhouYC.Investigationoftherelationshipbetween

elasticmodulusandhardnessbasedondepth-sensingindentationmea-surements[J].ActaMaterialia,2004,52:5397-5404.

6 结论

本文通过增大盲孔周围区域材料的弹性模量E来考虑因钻孔引起的加工硬化,应用有限元软件

ANSYS对304不锈钢的释放系数A、B的测定实验进行三维有限元数值模拟,并与实验测定结果进行比较,得出以下结论:

1)在释放系数的有限元数值计算中,必须考虑钻孔引起的加工硬化。对304不锈钢进行释放系数的有限元数值计算时,与实验测定值相比,当不考虑盲孔周围的加工硬化时,释放系数A、B的有限元计算值的误差分别高达13.4%和12.8%;当取盲孔附近的加工硬化层厚为0.2mm、弹性模量增大程度为M=137%时,释放系数A、B的有限元计算值的误差仅分别为-1.8%和2.4%。

2)释放系数随加工硬化程度的增大而增大。

3)考虑加工硬化后释放系数随孔深的变化规律与未考虑加工硬化时释放系数随孔深的变化规律相同。

本文来源：https://www.bwwdw.com/article/kco1.html