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

水冷盘式制动器热疲劳失效有限元分析蔡运迪唐文献黄秋芸王玲玲江苏科技大学，镇江，２１２００３摘要针对海洋钻井绞车在连续下钻过程中的制动工况，建立了水冷盘式制动器三维热机耦合分析模型，运用大型有限元分析软件ＡＢＡＱＵＳ数值模拟了制动器的制动过程，获得了制动盘表面及内部温度场与应力场的分布特征，并以此为基础分析了制动盘热疲劳失效的机理。研究结果表明连续制动工况下，制动盘的温度场与应力场相互耦合，两者具有相似的变化规律；周向热应力是形成制动盘表面初始裂纹的主要应力分量，在热应力反复作用下，该初始裂纹发展为粗大的裂纹，最终导致制动盘的断裂。分析结果与实际情况吻合较好，从而证明了该分析方法的正确性和可行性。关键词水冷盘式制动器；连续制动；温度场；应力场；热疲劳裂纹中图分类号ＴＥ９５１；ＴＫ１２１ＤＯＩ１０．３９６９／Ｊ．ＩＳＳＮ．１００４－１３２Ｘ．２０１２．２２．０１６ＦＩＮＩＴＥＥＬＥＭＥＮＴＡＮＡＬＹＳＩＳ（ＦＥＭ）ＯＮＴＨＥＲＭＡＬＦＡＴＩＧＵＥＦＡＩＬＵＲＥＯＦＷＡＴＥＲ－ＣＯＯＬＩＮＧＤＩＳＣＢＲＡＫＥＣＡＩＹＵＮＤＩＴＡＮＧＷＥＮＸＩＡＮＨＵＡＮＧＱＩＵＹＵＮＷＡＮＧＬＩＮＧＬＩＮＧＪＩＡＮＧＳＵＵＮＩＶＥＲＳＩＴＹＯＦＳＣＩＥＮＣＥＡＮＤＴＥＣＨＮＯＬＯＧＹ，ＺＨＥＮＪＩＡＮＧ，ＪＩＡＮＧＳＵ，２１２００３ＡＢＳＴＲＡＣＴＡＩＭＩＮＧＡＴＴＨＥＲＥＰＥＡＴＥＤＢＲＡＫＩＮＧＣＯＮＤＩＴＩＯＮＳＯＦＯＦＦＳＨＯＲＥＤＲＩＬＩＮＧＤＲＡＷＷＯＲＫＳＤＵＲＩＮＧＴＨＥＲＯＵＮＤＴＲＩＰ，ＡＴＨＲＥＥ－ＤＩＭＥＮＳＩＯＮＡＬＦＥＭＯＦＴＨＥＲＭＯ－ＭＥＣＨＡＮＩＣＡＬＣＯＵＰＬＩＮＧＡＮＡＬＹＳＩＳＯＦＴＨＥＷＡＴＥＲ－ＣＯＯＬＩＮＧＢＲＡＫＥＤＩＳＣＷＡＳＥＳＴＡＢＬＩＳＨＥＤ．ＢＹＵＳＩＮＧＴＨＥＦＥＡＳＯＦＴＷＡＲＥＡＢＡＱＵＳ，ＴＨＥＲＥＰＥＡＴＥＤＢＲＡＫＩＮＧＰＲＯＣＥＳＳＷＡＳＳＩＭＵＬＡＴＥＤ．ＤＩＳＴＲＩＢＵＴＩＯＮＣＨＡＲＡＣＴＥＲＩＳＴＩＣＳＯＦＴＨＥＳＵＲＦＡＣＥＡＮＤＩＮＮＥＲＴＥＭＰＥＲＡＴＵＲＥＦＩＥＬＤＡＮＤＳＴＲＥＳＳＦＩＥＬＤＷＥＲＥＯＢＴＡＩＮＥＤ，ＢＡＳＥＤＯＮＷＨＩＣＨ，ＴＨＥＴＨＥＲＭＡＬＦＡＴＩＧＵＥＦＡＩＬＵＲＥＭＥＣＨＡＮＩＳＭＯＦＢＲＡＫＥＤＩＳＣＷＡＳＡＮＡＬＹＺＥＤ．ＴＨＥＲＥ－ＳＵＬＴＳＳＨＯＷＴＨＡＴ，ＩＮＴＨＥＲＥＰＥＡＴＥＤＢＲＡＫＩＮＧ，ＴＨＥＴＥＭＰＥＲＡＴＵＲＥＡＮＤＳＴＲＥＳＳＯＦＢＲＡＫＥＤＩＳＣＡＲＥＣＯＵＰＬＥＤ，ＡＮＤＴＨＥＹＨＡＶＥＴＨＥＳＩＭＩＬＡＲＶＡＲＩＡＴＩＯＮＤＩＳＣＩＰＬＩＮＡＲＩＡＮ．ＴＨＥＲＭＡＬＳＴＲＥＳＳＩＮＴＨＥＣＩＲＣＵＬＡＲＤＩＲＥＣＴＩＯＮＩＳＴＨＥＭＡＩＮＳＴＲＥＳＳＦＯＲＣＡＵＳＩＮＧＴＨＥＩＮＩＴＩＡＬＣＲＡＣＫＳＯＮＴＨＥＳＵＲＦＡＣＥＯＦＢＲＡＫＥＤＩＳＣ．ＴＨＥＳＥＣＲＡＣＫＳＡＲＥＧＥＴＴＩＮＧＬＡＲＧＥＲＵＮＤＥＲＴＨＥＥＦＦＥＣＴＯＦＲＥＰＥＡＴＥＤＴＨＥＲＭＡＬＳＴＲＥＳＳ，ＥＶＥＮＴＵＡＬＹ，ＩＮＤＵＣＩＮＧＴＨＥＦＲＡＣＴＵＲＥＯＦＢＲＡＫＥＤＩＳＣ．ＴＨＥＡＮＡＬＹＳＩＳＲＥＳＵＬＴＳＡＲＥＩＮＧＯＯＤＡＧＲＥＥＭＥＮＴＷＩＴＨＴＨＥＰＲＡＣＴＩＣＡＬＳＩＴＵＡＴＩＯＮ，ＷＨＩＣＨＰＲＯＶＥＴＨＥＶＡＬＩＤＩＴＹＡＮＤＦＥＡＳＩＢＩＬＩＴＹＯＦＴＨＥＡＮＡＬＹＳＩＳＭＥＴＨＯＤＵＳＥＤＨＥＲＥＩＮ．ＫＥＹＷＯＲＤＳＷＡＴＥＲ－ＣＯＯＬＩＮＧＤＩＳＣＢＲＡＫＥ；ＲＥＰＥＡＴＥＤＢＲＡＫＩＮＧ；ＴＥＭＰＥＲＡＴＵＲＥＦＩＥＬＤ；ＳＴＲＥＳＳＦＩＥＬＤ；ＴＨＥＲＭＡＬＦＡ－ＴＩＧＵＥＣＲＡＣＫＳ收稿日期２０１１０８０１基金项目江苏省普通高校研究生科研创新计划资助项目０引言盘式制动器是海洋钻井绞车的关键装置，其工作可靠性直接关系到钻井作业的安全问题。在制动系统承担的诸多作业中，下放钻柱作业是载荷最大的一种［１］。在下放钻柱过程中，制动器需在相当长的时间内周期性地连续制动，每次制动过程中，制动器都要吸收由钻柱和游动系统动能转化而来的大量热能。实践表明，连续、重载的周期性制动所形成的热疲劳裂纹是海洋钻井绞车盘式制动器的主要失效形式，因此，有必要对海洋钻井绞车盘式制动器进行热－机耦合分析，研究制动盘的热疲劳裂纹的形成机理。关于盘式制动器的热－机耦合分析已有大量的研究［２－４］。ＡＶＩＬＥＳ等［５］针对制动盘建立了三维热－机耦合的有限元模型，并与台架试验作了比较，其研究表明，制动盘翘曲是产生热点和导致热弹性不稳定的主要原因。ＨＷＡＮＧ等［６］通过建立部分三维热－机耦合模型，研究了通风式制动盘的温度场和热变形情况。吕振华等［７］建立了蹄－鼓式制动器热弹性耦合分析模型，研究了摩擦副间接触力分布、制动鼓瞬态温度场、应力场和变形场。ＧＡＯ等［８］通过建立实心制动盘三维热－机耦合模型，分析了制动盘瞬态温度场的分布特征和热疲劳裂纹形成的原因。杨智勇等［９］建立了制动盘、闸片和盘毂之间的三体接触弹塑性热－机耦合模型，分析了制动盘的温度场和应力应变场的分布情况。虽然国内外学者在盘式制动器的热－机耦合分析方面做了许多工作，但对于海洋钻井绞车用水冷盘式制动器热－机耦合的研究尚处于空白。本文以某海洋钻井绞车水冷盘式制动器为研究对象，建立制动器的三维热－机耦合有限元分析模型，运用有限元软件ＡＢＡＱＵＳ分析连续起下钻６２７２中国机械工程第２３卷第２２期２０１２年１１月下半月工况下水冷盘式制动器热－机耦合特性，并以此为基础对制动盘热疲劳裂纹的形成机理进行分析。１热－机耦合分析模型１．１制动器有限元模型海洋钻井绞车盘式制动器结构见图１，它主要由制动盘、钳架、制动钳、液压缸等组成［１０］。制动钳分为工作钳和安全钳，在正常工作情况下，使用工作钳实现送钻、起下钻等工况下的制动，而安全钳不参与制动。制动钳为一钳式杠杆机构，液压油使液压缸活塞杆伸出，推动杠杆使摩擦片压紧在与绞车滚筒连为一体的制动盘上，从而实现绞车的制动。活塞杆收回时，摩擦片离开制动盘，制动力矩消失。图１盘式制动器结构示意图图２Ａ为水冷盘式制动器的简化模型，仅保留了参加制动的工作钳摩擦片和制动盘。制动盘内部设有水冷通道，起下钻作业时，给制动盘通冷却循环水，以平衡摩擦副产生的热量。由于制动器结构及载荷均具有对称性，在数值计算中，为节省计算资源和时间，取其一半建立有限元计算模型，如图２Ｂ所示。（Ａ）三维物理模型（Ｂ）有限元模型图２水冷盘式制动器三维有限元模型１．２热传导数学模型图３所示为水冷盘式制动器在柱坐标系（Ｒ，Θ，Ｚ）下的三维计算模型，图中Ｓ１～Ｓ５分别代表制动盘的摩擦表面、外圆侧面、内圆侧面、水冷通道内侧面和盘的中心平面，Ｐ１～Ｐ３分别代表摩擦片的侧面、接触摩擦表面和背面，ＲＤ、ＲＤ分别代表制动盘的内外圆半径，ＲＰ、ＲＰ分别代表摩擦片的内外圆半径。图３水冷盘式制动器三维计算模型假设摩擦副为各向同性材料，根据热传导理论，摩擦副在柱坐标系下的热传导微分方程为２ＴＲ２＋１ＲＴＲ＋１Ｒ２２ＴΘ２＋２ＴＺ２＝ΡＣΛＴＴ（１）式中，Λ为材料的导热系数；Ρ为材料的密度；Ｃ为材料的质量热容；Ｔ为摩擦界面对应点的瞬态温度；Ｔ为制动时间；Ｒ为盘、片有效作用半径。下面给出三维温度场计算的边界条件。在摩擦界面（Ｓ１和Ｐ２）上Λ（ＴＲＮＲ＋１ＲＴΘＮΘ＋ＴＺＮＺ）＝Ｑ（２）在所有换热界面（Ｓ２～Ｓ５、Ｐ１和Ｐ３）上Λ（ＴＲＮＲ＋１ＲＴΘＮΘ＋ＴＺＮＺ）＝Ｈ（Ｔ０－Ｔ）（３）三维温度场计算的初始条件为Ｔ（Ｒ，Θ，Ｚ，Ｔ）＝Ｔ０Ｔ＝０（４）式中，Ｔ０为环境温度；Ｑ为边界的热流密度；Ｈ为边界上的对流换热系数；ＮＲ、ＮΘ、ＮＺ为边界外法线的方向余弦。在制动过程中，不考虑材料的磨损，认为绞车的起升系统动能全部转化为热能，则摩擦副间产生的摩擦热流为Ｑ（Ｒ，Ｔ）＝ΜＰ（Ｒ，Ｔ）Ｖ（Ｔ）＝ΜＰ（Ｒ，Ｔ）Ω（Ｔ）Ｒ（５）式中，Μ为摩擦副间的摩擦因数；Ｐ为摩擦副间的接触压力；Ｖ为盘、片相对移动速度；Ω为制动盘的角速度。摩擦热流在摩擦副间按照热传导和材料物理特性自然分配，摩擦副间的热流分配系数为［１１］Κ＝ＱＤＱＰ＝（ΡＤＣＤΛＤΡＰＣＰΛＰ）１／２（６）式中，下标Ｄ、Ｐ分别代表制动盘和摩擦片。１．３热－机耦合计算模型盘式制动器在制动过程中处在温度场和力场的共同作用下，这是一个非常复杂的多物理场非线性耦合作用的问题，应同时考虑温度和应力的相互影响。制动器在制动时产生的摩擦热流Ｑ（Ｒ，Θ，Ｔ）使摩擦副间产生了不均匀的温度场Ｔ（Ｔ），在温度场７２７２水冷盘式制动器热疲劳失效有限元分析蔡运迪唐文献黄秋芸等Ｔ（Ｔ）的作用下摩擦副发生机械应变ΕＥ和热应变ΕＴ，从而改变了摩擦副间的接触关系和接触压力Ｐ（Ｒ，Θ，Ｔ），接触关系和接触压力的改变反过来影响到摩擦热流Ｑ（Ｒ，Θ，Ｔ）的大小，由此可见，制动过程是一个复杂的热－机耦合的过程。图４为制动器的热－机耦合关系图。图４制动器热－机耦合关系其中，盘、片的总应力为［１２］Σ＝｜Ｅ｜（ΕＥ＋ΕＴ）（７）式中，Ｅ为材料的弹性矩阵。单元应力Σ与节点上的接触压力Ｐ（Ｒ，Θ，Ｔ）之间的平衡关系可表示为∫Ｖ｜Ｂ｜ΣＤＶ＝Ｐ（Ｒ，Θ，Ｔ）（８）式中，Ｂ为节点位移和单元应变之间的转换矩阵。本文采用更新ＬＡＧＲＡＮＧＥ方法求解此热－机耦合问题，热－机耦合有限元方程如下ＫＵＭＴＭＵＣ［］ＵＵ（Ｔ）Ｔ（Ｔ熿燀燄燅）＝Ｆ（Ｔ）ＫＴＴ（Ｔ）＋Ｑ＋［］Ｈ（９）式中，ＫＵ为力学刚度矩阵；ＭＴ为热学刚度矩阵；Ｆ为载荷矢量；ＣＵ为热容矩阵；Ｕ为节点位移矢量；Ｔ为节点温度矢量；ＫＴ为热传导矩阵；ＭＵ为热力耦合矩阵；Ｑ为热载荷矢量；Ｈ为耗散矢量。在每个增量步开始时将几何形状更新，在新的ＬＡＧＲＡＮＧＥ坐标下分析温度场方程，采用非线性方程迭代解法求解热传导方程的等效温度场递推关系式。收敛后，在同一增量步中更新温度值，评价材料的力学性质和热应变，迭代求解力平衡方程，收敛后进行下一增量步的分析。２热边界条件制动器的散热部件主要是制动盘，因此主要分析制动盘热边界条件，而对于固定不动的摩擦片，在其非摩擦表面施加自然对流换热边界条件。制动盘外表面换热由对流换热和辐射换热两部分组成。对流换热分为两个阶段，即在制动过程中的强制对流换热和制动结束后的自然对流换热，因此，制动盘外表面综合换热系数为Ｈ＝ＨＤ＋ＨＲ制动过程中Ｈ′Ｄ＋ＨＲ｛制动结束后（１０）式中，ＨＤ为制动过程中制动盘外表面对流换热系数；ＨＲ为制动盘外表面辐射换热系数；Ｈ′Ｄ为制动结束后制动盘外表面对流换热系数。对于制动盘内部水冷通道表面，则是通过冷却循环水将热量带出制动盘，为强制对流换热。２．１对流换热２．１．１强制对流换热（１）制动盘外表面。制动盘外表面的对流换热系数可以采用有限空间自然对流换热系数的经验公式［１３］，即ＨＤ＝０．７０ＲＥ０．５５ＫＡ／ＤＲＥ≤２．４１０５０．０４ＲＥ０．８ＫＡ／ＤＲＥ＞２．４１０｛５（１１）ＲＥ＝ΩΡＡＤ２／（２ＵＡ）式中，ＫＡ为空气的导热系数；Ｄ为制动盘的外直径；ＲＥ为空气的雷诺数；ΡＡ为空气的密度；ＵＡ为空气的动力黏度。（２）制动盘内部水冷通道表面。对于制动盘内部水冷通道中冷却循环水的强制对流换热，使用最广的关联式为迪图斯－贝尔特（ＤＩＴＴＵＳ－ＢＯＥＬＴＥＲ）公式［１３］，即ＮＵＦ＝０．０２３ＲＥ０．８ＦＰＲＮＦ＝ＨＦＤＥΛＦ（１２）式中，ＮＵＦ为水的努塞尔数；ＲＥＦ为水的雷诺数；ＰＲＦ为水的普朗特数；ＨＦ为冷却循环水与制动盘内部水冷通道间的对流换热系数；ＤＥ为水冷通道的当量直径，对于截面为矩形的水冷通道，ＤＥ＝２ＡＢ／（Ａ＋Ｂ），其中Ａ、Ｂ分别为水冷通道截面矩形的长和宽；ΛＦ为水的导热系数。加热流体时Ｎ＝０．４，式（１２）中采用流体平均温度为定性温度，由此可得ＨＦ＝０．０２３ΛＦＤＥＲＥ０．８ＦＰＲ０．４Ｆ（１３）２．１．２自然对流换热由于制动器被制动机构外的安全外壳罩住，因此，制动结束后，制动盘外表面对流换热按有限空间的自然对流换热计算，可知其对流换热系数为［１３］Ｈ′Ｄ＝ＮＵＡＫＡＬＣ（１４）ＮＵＡ＝０．１９７（ＧＲΔＰＲＡ）１／４（Ｈ／Δ）－１／９ＧＲΔ＝ＧΑＶΔＴΔ３／Ν２式中，ＮＵＡ为空气的努塞尔数；ＬＣ为特征长度，即制动盘的直径；ＧＲΔ为格拉晓夫数；Ｇ为重力加速度；ΑＶ为体胀系数，ΑＶ≈１／Ｔ；ΔＴ为表面温度与环境温度的差；Ν为空气的运动黏度；ＰＲＡ为空气的普朗特数；Ｈ、Δ分别为有限空间的高度和宽度。２．２辐射换热热辐射能的计算是一个高度非线性的问题，处理辐射换热问题的思路是将辐射换热边界条件转化为物体表面的对流换热边界。热辐射过程中，辐射热流密度为ＱＲ＝ΨΕＤ（Ｔ４Ｄ－Ｔ４０）＝ＨＲ（ＴＤ－Ｔ０）（１５）式中，Ψ为斯蒂芬－波尔兹曼（ＳＴＥＦＡＮ－ＢＯＬＴＺＭＡＮＮ）常数；８２７２中国机械工程第２３卷第２２期２０１２年１１月下半月ΕＤ为黑体辐射系数；ＴＤ、Ｔ０分别为制动盘瞬时温度和周围环境温度。则辐射换热系数折算成的对流换热系数为ＨＲ＝ΨΕＤＴ４Ｄ－Ｔ４０ＴＤ－Ｔ０＝ΨΕＤ（Ｔ３Ｄ＋Ｔ３０）（ＴＤ－Ｔ０）（１６）３材料特性参数制动盘的内外半径分别为ＲＤ＝５８０ＭＭ和ＲＤ＝９９０ＭＭ，厚度为９０ＭＭ，摩擦片的内外半径分别为ＲＰ＝７４５ＭＭ和ＲＰ＝９７５ＭＭ。制动盘的材料为３５ＣＲＭＯ中碳调质钢，摩擦片为无石棉有机摩擦材料，摩擦副材料热物性参数如表１所示。制动盘材料３５ＣＲＭＯ的热物性随温度变化，其参数如表２所示，温度点之间的数值可通过线性插值得到。表１摩擦副材料热物性参数制动盘摩擦片密度Ρ（ＫＧ／Ｍ３）７８６６２２０６弹性模量Ｅ（ＧＰＡ）２０６．０５．４导热系数Λ（Ｗ／（ＭＫ））１．７质量热容Ｃ（Ｊ／（ＫＧＫ））１２０４热膨胀系数Α（１０－６Ｋ－１）３０．０泊松比０．２８０．３０表２制动盘材料３５ＣＲＭＯ热物性参数温度（℃）导热系数Λ（Ｗ／（ＭＫ））质量热容Ｃ（Ｊ／（ＫＧＫ））热膨胀系数Α（１０－６Ｋ－１）２０４８．０５２０１２．０１００４７．７５６１１２．５２００４７．７５９９１３．１３００４４．０６１１１３．６４００４１．０６５７１４．０５００３８．１７１６１４．３４仿真计算与结果分析下钻作业的制动周期一般为９０～１５０Ｓ，本文取制动周期为１２３Ｓ，其中制动时间３Ｓ，松闸停留时间１２０Ｓ。本文将模拟绞车在６次连续下钻工作中的制动过程，取钻柱下放速度为２Ｍ／Ｓ，并将该制动过程近似为匀减速。随着下放钻柱数量的增大，制动载荷也增大，为了模拟这一工况，６个周期内作用于摩擦片背部的制动压力依次在０．９～１．２ＭＰＡ范围内等间距递增取值，并认为在各自制动周期内制动压力值恒定不变。摩擦副初始温度、冷却循环水温度和环境温度均为２０℃。４．１制动盘温度场分布特征分析图５为第一个制动周期内制动盘温度场分布云图。制动结束时，制动盘温度场分布很不均匀，有较浅的热影响层，在摩擦区域形成很大的温度梯度。高温区主要分布在表面摩擦环附近，呈现出带状，最高温度出现在制动盘表面的摩擦区域外沿。当第一个制动周期结束时，制动盘在热传导、热对流及冷却循环水的共同作用下，温度场变得更加均匀，整个制动盘的温度比较低。（Ａ）制动结束时（Ｔ＝３Ｓ）（Ｂ）制动周期结束时（Ｔ＝１２３Ｓ）图５第一个制动周期内制动盘温度场分布在摩擦半径Ｒ＝９００ＭＭ、Θ＝４５处的表面节点７４０６为６次循环制动周期内温度值最大节点，该节点温度随时间变化的曲线如图６所示。从图６可以看出，制动开始时，节点温度急剧上升；制动结束后，节点温度快速下降，在后期下降趋势变得平缓。随着制动次数和制动载荷的增大，节点温度的最大值也由第一个制动周期的１８６．５℃逐渐增大到第六个制动周期的２９５．０℃；在每个制动周期结束的时刻，节点温度为最低值，均维持在２３．６～２４．８℃范围之内，结合实际工程经验可知，这主要是由于水冷通道内冷却循环水对制动盘强制冷却。图６节点７４０６在６个制动周期内温度变化曲线图７所示为前三个制动周期内，半径Ｒ＝９００ＭＭ、Θ＝１３２处节点在制动盘厚度方向上的温度变化曲线，可以看出，制动盘表面节点的温度波动较大，也是最先达到最大值的。随着厚度的增９２７２水冷盘式制动器热疲劳失效有限元分析蔡运迪唐文献黄秋芸等