会员注册 | 登录 | 微信快捷登录 支付宝快捷登录 QQ登录 微博登录 | 帮助中心 人人文库renrendoc.com美如初恋!
站内搜索 百度文库

热门搜索: 直缝焊接机 矿井提升机 循环球式转向器图纸 机器人手爪发展史 管道机器人dwg 动平衡试验台设计

水冷盘式制动器热疲劳失效有限元分析.pdf水冷盘式制动器热疲劳失效有限元分析.pdf -- 1 元

宽屏显示 收藏 分享

资源预览需要最新版本的Flash Player支持。
您尚未安装或版本过低,建议您

水冷盘式制动器热疲劳失效有限元分析蔡运迪唐文献黄秋芸王玲玲江苏科技大学,镇江,212003摘要针对海洋钻井绞车在连续下钻过程中的制动工况,建立了水冷盘式制动器三维热机耦合分析模型,运用大型有限元分析软件ABAQUS数值模拟了制动器的制动过程,获得了制动盘表面及内部温度场与应力场的分布特征,并以此为基础分析了制动盘热疲劳失效的机理。研究结果表明连续制动工况下,制动盘的温度场与应力场相互耦合,两者具有相似的变化规律周向热应力是形成制动盘表面初始裂纹的主要应力分量,在热应力反复作用下,该初始裂纹发展为粗大的裂纹,最终导致制动盘的断裂。分析结果与实际情况吻合较好,从而证明了该分析方法的正确性和可行性。关键词水冷盘式制动器连续制动温度场应力场热疲劳裂纹中图分类号TE951TK121DOI10.3969/j.issn.1004-132X.2012.22.016FiniteElementAnalysis(FEM)onThermalFatigueFailureofWater-coolingDiscBrakeCaiYundiTangWenxianHuangQiuyunWangLinglingJiangsuUniversityofScienceandTechnology,Zhenjiang,Jiangsu,212003AbstractAimingattherepeatedbrakingconditionsofoffshoredrilingdrawworksduringtheroundtrip,athree-dimensionalFEMofthermo-mechanicalcouplinganalysisofthewater-coolingbrakediscwasestablished.ByusingtheFEAsoftwareABAQUS,therepeatedbrakingprocesswassimulated.Distributioncharacteristicsofthesurfaceandinnertemperaturefieldandstressfieldwereobtained,basedonwhich,thethermalfatiguefailuremechanismofbrakediscwasanalyzed.There-sultsshowthat,intherepeatedbraking,thetemperatureandstressofbrakediscarecoupled,andtheyhavethesimilarvariationdisciplinarian.Thermalstressinthecirculardirectionisthemainstressforcausingtheinitialcracksonthesurfaceofbrakedisc.Thesecracksaregettinglargerundertheeffectofrepeatedthermalstress,eventualy,inducingthefractureofbrakedisc.Theanalysisresultsareingoodagreementwiththepracticalsituation,whichprovethevalidityandfeasibilityoftheanalysismethodusedherein.Keywordswater-coolingdiscbrakerepeatedbrakingtemperaturefieldstressfieldthermalfa-tiguecracks收稿日期20110801基金项目江苏省普通高校研究生科研创新计划资助项目0引言盘式制动器是海洋钻井绞车的关键装置,其工作可靠性直接关系到钻井作业的安全问题。在制动系统承担的诸多作业中,下放钻柱作业是载荷最大的一种[1]。在下放钻柱过程中,制动器需在相当长的时间内周期性地连续制动,每次制动过程中,制动器都要吸收由钻柱和游动系统动能转化而来的大量热能。实践表明,连续、重载的周期性制动所形成的热疲劳裂纹是海洋钻井绞车盘式制动器的主要失效形式,因此,有必要对海洋钻井绞车盘式制动器进行热-机耦合分析,研究制动盘的热疲劳裂纹的形成机理。关于盘式制动器的热-机耦合分析已有大量的研究[2-4]。Aviles等[5]针对制动盘建立了三维热-机耦合的有限元模型,并与台架试验作了比较,其研究表明,制动盘翘曲是产生热点和导致热弹性不稳定的主要原因。Hwang等[6]通过建立部分三维热-机耦合模型,研究了通风式制动盘的温度场和热变形情况。吕振华等[7]建立了蹄-鼓式制动器热弹性耦合分析模型,研究了摩擦副间接触力分布、制动鼓瞬态温度场、应力场和变形场。Gao等[8]通过建立实心制动盘三维热-机耦合模型,分析了制动盘瞬态温度场的分布特征和热疲劳裂纹形成的原因。杨智勇等[9]建立了制动盘、闸片和盘毂之间的三体接触弹塑性热-机耦合模型,分析了制动盘的温度场和应力应变场的分布情况。虽然国内外学者在盘式制动器的热-机耦合分析方面做了许多工作,但对于海洋钻井绞车用水冷盘式制动器热-机耦合的研究尚处于空白。本文以某海洋钻井绞车水冷盘式制动器为研究对象,建立制动器的三维热-机耦合有限元分析模型,运用有限元软件ABAQUS分析连续起下钻6272中国机械工程第23卷第22期2012年11月下半月工况下水冷盘式制动器热-机耦合特性,并以此为基础对制动盘热疲劳裂纹的形成机理进行分析。1热-机耦合分析模型1.1制动器有限元模型海洋钻井绞车盘式制动器结构见图1,它主要由制动盘、钳架、制动钳、液压缸等组成[10]。制动钳分为工作钳和安全钳,在正常工作情况下,使用工作钳实现送钻、起下钻等工况下的制动,而安全钳不参与制动。制动钳为一钳式杠杆机构,液压油使液压缸活塞杆伸出,推动杠杆使摩擦片压紧在与绞车滚筒连为一体的制动盘上,从而实现绞车的制动。活塞杆收回时,摩擦片离开制动盘,制动力矩消失。图1盘式制动器结构示意图图2a为水冷盘式制动器的简化模型,仅保留了参加制动的工作钳摩擦片和制动盘。制动盘内部设有水冷通道,起下钻作业时,给制动盘通冷却循环水,以平衡摩擦副产生的热量。由于制动器结构及载荷均具有对称性,在数值计算中,为节省计算资源和时间,取其一半建立有限元计算模型,如图2b所示。(a)三维物理模型(b)有限元模型图2水冷盘式制动器三维有限元模型1.2热传导数学模型图3所示为水冷盘式制动器在柱坐标系(r,θ,z)下的三维计算模型,图中S1~S5分别代表制动盘的摩擦表面、外圆侧面、内圆侧面、水冷通道内侧面和盘的中心平面,P1~P3分别代表摩擦片的侧面、接触摩擦表面和背面,rd、Rd分别代表制动盘的内外圆半径,rp、Rp分别代表摩擦片的内外圆半径。图3水冷盘式制动器三维计算模型假设摩擦副为各向同性材料,根据热传导理论,摩擦副在柱坐标系下的热传导微分方程为2Tr2+1rTr+1r22Tθ2+2Tz2=ρcλTt(1)式中,λ为材料的导热系数ρ为材料的密度c为材料的质量热容T为摩擦界面对应点的瞬态温度t为制动时间r为盘、片有效作用半径。下面给出三维温度场计算的边界条件。在摩擦界面(S1和P2)上λ(Trnr+1rTθnθ+Tznz)=q(2)在所有换热界面(S2~S5、P1和P3)上λ(Trnr+1rTθnθ+Tznz)=h(T0-T)(3)三维温度场计算的初始条件为T(r,θ,z,t)=T0t=0(4)式中,T0为环境温度q为边界的热流密度h为边界上的对流换热系数nr、nθ、nz为边界外法线的方向余弦。在制动过程中,不考虑材料的磨损,认为绞车的起升系统动能全部转化为热能,则摩擦副间产生的摩擦热流为q(r,t)=μp(r,t)v(t)=μp(r,t)ω(t)r(5)式中,μ为摩擦副间的摩擦因数p为摩擦副间的接触压力v为盘、片相对移动速度ω为制动盘的角速度。摩擦热流在摩擦副间按照热传导和材料物理特性自然分配,摩擦副间的热流分配系数为[11]κ=qdqp=(ρdcdλdρpcpλp)1/2(6)式中,下标d、p分别代表制动盘和摩擦片。1.3热-机耦合计算模型盘式制动器在制动过程中处在温度场和力场的共同作用下,这是一个非常复杂的多物理场非线性耦合作用的问题,应同时考虑温度和应力的相互影响。制动器在制动时产生的摩擦热流q(r,θ,t)使摩擦副间产生了不均匀的温度场T(t),在温度场7272水冷盘式制动器热疲劳失效有限元分析蔡运迪唐文献黄秋芸等T(t)的作用下摩擦副发生机械应变εe和热应变εT,从而改变了摩擦副间的接触关系和接触压力p(r,θ,t),接触关系和接触压力的改变反过来影响到摩擦热流q(r,θ,t)的大小,由此可见,制动过程是一个复杂的热-机耦合的过程。图4为制动器的热-机耦合关系图。图4制动器热-机耦合关系其中,盘、片的总应力为[12]σ=|E|(εe+εT)(7)式中,E为材料的弹性矩阵。单元应力σ与节点上的接触压力p(r,θ,t)之间的平衡关系可表示为∫V|B|σdV=p(r,θ,t)(8)式中,B为节点位移和单元应变之间的转换矩阵。本文采用更新Lagrange方法求解此热-机耦合问题,热-机耦合有限元方程如下KuMTMuC[]uu(t)T(t熿燀燄燅)=F(t)KTT(t)+Q+[]H(9)式中,Ku为力学刚度矩阵MT为热学刚度矩阵F为载荷矢量Cu为热容矩阵u为节点位移矢量T为节点温度矢量KT为热传导矩阵Mu为热力耦合矩阵Q为热载荷矢量H为耗散矢量。在每个增量步开始时将几何形状更新,在新的Lagrange坐标下分析温度场方程,采用非线性方程迭代解法求解热传导方程的等效温度场递推关系式。收敛后,在同一增量步中更新温度值,评价材料的力学性质和热应变,迭代求解力平衡方程,收敛后进行下一增量步的分析。2热边界条件制动器的散热部件主要是制动盘,因此主要分析制动盘热边界条件,而对于固定不动的摩擦片,在其非摩擦表面施加自然对流换热边界条件。制动盘外表面换热由对流换热和辐射换热两部分组成。对流换热分为两个阶段,即在制动过程中的强制对流换热和制动结束后的自然对流换热,因此,制动盘外表面综合换热系数为h=hd+hr制动过程中h′d+hr{制动结束后(10)式中,hd为制动过程中制动盘外表面对流换热系数hr为制动盘外表面辐射换热系数h′d为制动结束后制动盘外表面对流换热系数。对于制动盘内部水冷通道表面,则是通过冷却循环水将热量带出制动盘,为强制对流换热。2.1对流换热2.1.1强制对流换热(1)制动盘外表面。制动盘外表面的对流换热系数可以采用有限空间自然对流换热系数的经验公式[13],即hd=0.70Re0.55ka/DRe≤2.41050.04Re0.8ka/DRe>2.410{5(11)Re=ωρaD2/(2ua)式中,ka为空气的导热系数D为制动盘的外直径Re为空气的雷诺数ρa为空气的密度ua为空气的动力黏度。(2)制动盘内部水冷通道表面。对于制动盘内部水冷通道中冷却循环水的强制对流换热,使用最广的关联式为迪图斯-贝尔特(Dittus-Boelter)公式[13],即Nuf=0.023Re0.8fPrnf=hfdeλf(12)式中,Nuf为水的努塞尔数Ref为水的雷诺数Prf为水的普朗特数hf为冷却循环水与制动盘内部水冷通道间的对流换热系数de为水冷通道的当量直径,对于截面为矩形的水冷通道,de=2ab/(a+b),其中a、b分别为水冷通道截面矩形的长和宽λf为水的导热系数。加热流体时n=0.4,式(12)中采用流体平均温度为定性温度,由此可得hf=0.023λfdeRe0.8fPr0.4f(13)2.1.2自然对流换热由于制动器被制动机构外的安全外壳罩住,因此,制动结束后,制动盘外表面对流换热按有限空间的自然对流换热计算,可知其对流换热系数为[13]h′d=NuakaLc(14)Nua=0.197(GrδPra)1/4(H/δ)-1/9Grδ=gαVΔTδ3/ν2式中,Nua为空气的努塞尔数Lc为特征长度,即制动盘的直径Grδ为格拉晓夫数g为重力加速度αV为体胀系数,αV≈1/TΔT为表面温度与环境温度的差ν为空气的运动黏度Pra为空气的普朗特数H、δ分别为有限空间的高度和宽度。2.2辐射换热热辐射能的计算是一个高度非线性的问题,处理辐射换热问题的思路是将辐射换热边界条件转化为物体表面的对流换热边界。热辐射过程中,辐射热流密度为qr=ψεd(T4d-T40)=hr(Td-T0)(15)式中,ψ为斯蒂芬-波尔兹曼(Stefan-Boltzmann)常数8272中国机械工程第23卷第22期2012年11月下半月εd为黑体辐射系数Td、T0分别为制动盘瞬时温度和周围环境温度。则辐射换热系数折算成的对流换热系数为hr=ψεdT4d-T40Td-T0=ψεd(T3d+T30)(Td-T0)(16)3材料特性参数制动盘的内外半径分别为rd=580mm和Rd=990mm,厚度为90mm,摩擦片的内外半径分别为rp=745mm和Rp=975mm。制动盘的材料为35CrMo中碳调质钢,摩擦片为无石棉有机摩擦材料,摩擦副材料热物性参数如表1所示。制动盘材料35CrMo的热物性随温度变化,其参数如表2所示,温度点之间的数值可通过线性插值得到。表1摩擦副材料热物性参数制动盘摩擦片密度ρ(kg/m3)78662206弹性模量E(GPa)206.05.4导热系数λ(W/(mK))1.7质量热容c(J/(kgK))1204热膨胀系数α(10-6K-1)30.0泊松比0.280.30表2制动盘材料35CrMo热物性参数温度(℃)导热系数λ(W/(mK))质量热容c(J/(kgK))热膨胀系数α(10-6K-1)2048.052012.010047.756112.520047.759913.130044.061113.640041.065714.050038.171614.34仿真计算与结果分析下钻作业的制动周期一般为90~150s,本文取制动周期为123s,其中制动时间3s,松闸停留时间120s。本文将模拟绞车在6次连续下钻工作中的制动过程,取钻柱下放速度为2m/s,并将该制动过程近似为匀减速。随着下放钻柱数量的增大,制动载荷也增大,为了模拟这一工况,6个周期内作用于摩擦片背部的制动压力依次在0.9~1.2MPa范围内等间距递增取值,并认为在各自制动周期内制动压力值恒定不变。摩擦副初始温度、冷却循环水温度和环境温度均为20℃。4.1制动盘温度场分布特征分析图5为第一个制动周期内制动盘温度场分布云图。制动结束时,制动盘温度场分布很不均匀,有较浅的热影响层,在摩擦区域形成很大的温度梯度。高温区主要分布在表面摩擦环附近,呈现出带状,最高温度出现在制动盘表面的摩擦区域外沿。当第一个制动周期结束时,制动盘在热传导、热对流及冷却循环水的共同作用下,温度场变得更加均匀,整个制动盘的温度比较低。(a)制动结束时(t=3s)(b)制动周期结束时(t=123s)图5第一个制动周期内制动盘温度场分布在摩擦半径R=900mm、θ=45°处的表面节点7406为6次循环制动周期内温度值最大节点,该节点温度随时间变化的曲线如图6所示。从图6可以看出,制动开始时,节点温度急剧上升制动结束后,节点温度快速下降,在后期下降趋势变得平缓。随着制动次数和制动载荷的增大,节点温度的最大值也由第一个制动周期的186.5℃逐渐增大到第六个制动周期的295.0℃在每个制动周期结束的时刻,节点温度为最低值,均维持在23.6~24.8℃范围之内,结合实际工程经验可知,这主要是由于水冷通道内冷却循环水对制动盘强制冷却。图6节点7406在6个制动周期内温度变化曲线图7所示为前三个制动周期内,半径R=900mm、θ=132°处节点在制动盘厚度方向上的温度变化曲线,可以看出,制动盘表面节点的温度波动较大,也是最先达到最大值的。随着厚度的增9272水冷盘式制动器热疲劳失效有限元分析蔡运迪唐文献黄秋芸等
编号:201311161814335016    大小:589.89KB    格式:PDF    上传时间:2013-11-16
  【编辑】
1
关 键 词:
行业资料 机械制造 精品文档 水冷盘式
温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 人人文库网仅提供交流平台,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
  人人文库网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

当前资源信息

4.0
 
(2人评价)
浏览:11次
zhike上传于2013-11-16

官方联系方式

客服手机:13961746681   
2:不支持迅雷下载,请使用浏览器下载   
3:不支持QQ浏览器下载,请用其他浏览器   
4:下载后的文档和图纸-无水印   
5:文档经过压缩,下载后原文更清晰   

相关资源

相关资源

相关搜索

行业资料   机械制造   精品文档   水冷盘式  
关于我们 - 网站声明 - 网站地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5