CN114936436B - Method for establishing spiral seal abrasion model of roller bit under high-temperature and dynamic load working conditions - Google Patents
Method for establishing spiral seal abrasion model of roller bit under high-temperature and dynamic load working conditions Download PDFInfo
- Publication number
- CN114936436B CN114936436B CN202210635214.0A CN202210635214A CN114936436B CN 114936436 B CN114936436 B CN 114936436B CN 202210635214 A CN202210635214 A CN 202210635214A CN 114936436 B CN114936436 B CN 114936436B
- Authority
- CN
- China
- Prior art keywords
- wear
- model
- abrasion
- working condition
- temperature
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000005299 abrasion Methods 0.000 title claims abstract 47
- 238000007789 sealing Methods 0.000 claims abstract description 33
- 238000002474 experimental method Methods 0.000 claims abstract description 26
- 239000000463 material Substances 0.000 claims abstract description 21
- 238000012360 testing method Methods 0.000 claims abstract description 21
- 230000007246 mechanism Effects 0.000 claims abstract description 14
- 230000006872 improvement Effects 0.000 claims abstract description 7
- 238000004088 simulation Methods 0.000 claims description 56
- 230000008859 change Effects 0.000 claims description 20
- 230000003068 static effect Effects 0.000 claims description 19
- 230000020169 heat generation Effects 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 8
- 238000012795 verification Methods 0.000 claims description 7
- 238000011161 development Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 4
- 238000005553 drilling Methods 0.000 description 11
- 238000006073 displacement reaction Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 238000012937 correction Methods 0.000 description 3
- 230000005489 elastic deformation Effects 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 1
- 101000827703 Homo sapiens Polyphosphoinositide phosphatase Proteins 0.000 description 1
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 1
- 102100023591 Polyphosphoinositide phosphatase Human genes 0.000 description 1
- 101100233916 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) KAR5 gene Proteins 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000007373 indentation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
Description
技术领域Technical Field
本发明涉及磨损模型领域,特别涉及一种高温、动载工况下牙轮钻头螺旋密封磨损模型的建立方法。The invention relates to the field of wear models, and in particular to a method for establishing a wear model of a spiral seal of a roller drill under high temperature and dynamic load conditions.
背景技术Background technique
2016年,周已等人提出基于“排砂”思路的螺旋组合密封结构,设计了该密封结构及密封装置,并采用仿真和实验方法研究其密封性能及磨损性能,将其应用于8 1/2SLM537G-1型号和12 1/4 SLM517GK-1型号的牙轮钻头实际钻井中,取得了较好效果。In 2016, Zhou Yi et al. proposed a spiral combined sealing structure based on the "sand removal" idea, designed the sealing structure and sealing device, and used simulation and experimental methods to study its sealing performance and wear performance. It was applied to the actual drilling of 8 1/2SLM537G-1 and 12 1/4 SLM517GK-1 roller bits, and achieved good results.
对于牙轮钻头螺旋密封,密封失效导致钻头失效,钻井成本增加,钻井效率降低,起下钻次数增加。每年因密封失效造成的钻头失效等直接经济损失不可估量,因密封失效导致的安全生产事故等间接经济损失更大。For the spiral seal of the roller bit, seal failure leads to drill bit failure, increased drilling costs, reduced drilling efficiency, and increased tripping times. The direct economic losses caused by seal failure such as drill bit failure are incalculable every year, and the indirect economic losses such as safety accidents caused by seal failure are even greater.
因此,建立牙轮钻头螺旋密封磨损模型,预测牙轮钻头螺旋密封的磨损程度,对提高钻头密封性能,提高钻头寿命,降低钻井成本,提高钻井效率具有重要意义。Therefore, establishing a wear model for the spiral seal of a roller drill bit and predicting the wear degree of the spiral seal of a roller drill bit are of great significance for improving the sealing performance of the drill bit, increasing the life of the drill bit, reducing drilling costs, and improving drilling efficiency.
而现有磨损模型中主要是针对静载、常温工况,高温、动载工况下的磨损模型的研究几乎很少。由于高温、动载作用下,材料的物理、化学及力学性能都要发生改变,其摩擦磨损性能也会受到很大的影响。因此,静载、常温工况磨损模型的研究思路并不适合高温、动载工况下磨损模型的研究,并不能满足牙轮钻头螺旋密封在日常工作中高温、动载工况下的要求。The existing wear models are mainly for static load and normal temperature conditions, and there are few studies on wear models under high temperature and dynamic load conditions. Because the physical, chemical and mechanical properties of materials will change under high temperature and dynamic load, their friction and wear properties will also be greatly affected. Therefore, the research ideas of static load and normal temperature wear models are not suitable for the research of wear models under high temperature and dynamic load conditions, and cannot meet the requirements of roller drill bit spiral seals under high temperature and dynamic load conditions in daily work.
为此,如何提供能够满足牙轮钻头螺旋密封在日常工作中高温、动载工况环境下的要求,对牙轮钻头螺旋密封的磨损程度进行合理预测的高温、动载工况下牙轮钻头螺旋密封磨损模型的建立方法是本领域技术人员亟需解决的问题。Therefore, how to provide a method for establishing a wear model of the spiral seal of a roller drill bit under high temperature and dynamic load conditions that can meet the requirements of the spiral seal of a roller drill bit in daily work under high temperature and dynamic load conditions and reasonably predict the wear degree of the spiral seal of the roller drill bit is an urgent problem that technical personnel in this field need to solve.
发明内容Summary of the invention
有鉴于此,本发明提出了高温、动载工况下牙轮钻头螺旋密封磨损模型的建立方法。首先,通过球-块摩擦实验,得出密封结构摩擦副材料在不同温度和不同载荷下的摩擦磨损性能,选择Archard磨损模型作为高温工况下牙轮钻头螺旋密封磨损的基础模型和动载工况下牙轮钻头螺旋密封磨损的基础模型;然后,通过在Archard磨损模型的基础上引入温度影响参数,并进行拟合,输入球-块摩擦实验的数据得到最优解,建立高温工况下牙轮钻头螺旋密封磨损模型;通过对Archard磨损模型进行改进,对动态载荷下的球块试件之间的磨损机理进行公式推导,再根据赫兹点接触理论,推导出磨损体积与时间的关系式,建立动载工况下牙轮钻头螺旋密封磨损模型;最后,对高温工况下牙轮钻头螺旋密封磨损模型进行仿真验证,包括:静态接触仿真,动态摩擦生热仿真,动态摩擦磨损仿真,得到较准确的磨损体积,通过磨损体积,据几何关系计算磨损对螺旋圈槽深的改变量,基于最优槽深值分析,得到当前高温工况基本不会对密封效果造成影响;对动载工况下牙轮钻头螺旋密封磨损模型仿真,采用有限元仿真模拟动载和静载下的球-块磨损,将仿真结果与试验结果对比,验证理论模型和试验的可靠性与正确性,得出了动态载荷峰值对于密封圈螺纹上的接触应力影响极大,且成线性关系;动态载荷峰值越大密封圈磨损也会越来越剧烈;频率对于密封圈表面的接触应力几乎没有影响。上述建立的高温、动载工况下牙轮钻头螺旋密封磨损模型能够准确模拟出高温、动载工况下牙轮钻头螺旋密封的磨损机理,为实际牙轮钻头螺旋密封结构的防磨改进设计方案提供理论支撑,对提高牙轮钻头螺旋密封性能,提高牙轮钻头螺旋密封寿命,降低钻井成本,提高钻井效率具有重要意义。In view of this, the present invention proposes a method for establishing a wear model of the spiral seal of a roller drill bit under high temperature and dynamic load conditions. First, through the ball-block friction experiment, the friction and wear performance of the sealing structure friction pair material under different temperatures and different loads is obtained, and the Archard wear model is selected as the basic model of the spiral seal wear of the roller drill bit under high temperature conditions and the basic model of the spiral seal wear of the roller drill bit under dynamic load conditions; then, by introducing temperature influencing parameters on the basis of the Archard wear model and fitting, the data of the ball-block friction experiment is input to obtain the optimal solution, and the wear model of the spiral seal of the roller drill bit under high temperature conditions is established; by improving the Archard wear model, the wear mechanism between the ball-block specimens under dynamic load is derived, and then according to the Hertz point contact theory, the relationship between the wear volume and time is derived, and the wear model of the spiral seal of the roller drill bit under dynamic load conditions is established; finally, the high temperature The wear model of the spiral seal of the roller drill bit under working conditions is simulated and verified, including: static contact simulation, dynamic friction heat generation simulation, dynamic friction wear simulation, to obtain a more accurate wear volume. Through the wear volume, the change in the wear on the spiral groove depth is calculated according to the geometric relationship. Based on the analysis of the optimal groove depth value, it is obtained that the current high temperature working condition will basically not affect the sealing effect; for the simulation of the wear model of the spiral seal of the roller drill bit under dynamic load conditions, finite element simulation is used to simulate the ball-block wear under dynamic and static loads, and the simulation results are compared with the experimental results to verify the reliability and correctness of the theoretical model and the experiment. It is concluded that the dynamic load peak has a great influence on the contact stress on the sealing ring thread, and the relationship is linear; the larger the dynamic load peak, the more severe the wear of the sealing ring will be; the frequency has almost no effect on the contact stress on the surface of the sealing ring. The above-mentioned wear model of the spiral seal of the roller drill bit under high temperature and dynamic load conditions can accurately simulate the wear mechanism of the spiral seal of the roller drill bit under high temperature and dynamic load conditions, and provide theoretical support for the anti-wear improvement design scheme of the spiral seal structure of the actual roller drill bit. It is of great significance to improve the performance of the spiral seal of the roller drill bit, increase the life of the spiral seal of the roller drill bit, reduce drilling costs, and improve drilling efficiency.
为了实现上述目的,本发明采用如下技术方案:In order to achieve the above object, the present invention adopts the following technical solution:
高温、动载工况下牙轮钻头螺旋密封磨损模型的建立方法,包括:The method for establishing the wear model of the spiral seal of a roller drill bit under high temperature and dynamic load conditions includes:
步骤(1):通过球-块摩擦实验得出密封结构摩擦副材料在不同温度和不同载荷下的摩擦磨损性能,选择高温工况下牙轮钻头螺旋密封磨损的基础模型和动载工况下牙轮钻头螺旋密封磨损的基础模型。Step (1): The friction and wear performance of the friction pair material of the sealing structure under different temperatures and different loads is obtained through ball-block friction experiments, and the basic model of spiral seal wear of the roller bit under high temperature conditions and the basic model of spiral seal wear of the roller bit under dynamic load conditions are selected.
步骤(2)基于高温工况下牙轮钻头螺旋密封磨损的基础模型,引入温度影响参数进行修正,得到高温工况下牙轮钻头螺旋密封磨损的初步模型;对高温工况下牙轮钻头螺旋密封磨损初步模型进行拟合,输入球-块摩擦实验的数据得到最优解,建立高温工况下牙轮钻头螺旋密封磨损模型。Step (2) is based on the basic model of spiral seal wear of the roller drill bit under high temperature conditions, and the temperature influencing parameters are introduced for correction to obtain a preliminary model of spiral seal wear of the roller drill bit under high temperature conditions; the preliminary model of spiral seal wear of the roller drill bit under high temperature conditions is fitted, and the data of the ball-block friction experiment is input to obtain the optimal solution, and the spiral seal wear model of the roller drill bit under high temperature conditions is established.
步骤(3):基于动载工况下牙轮钻头螺旋密封磨损的基础模型进行改进,对动态载荷下的球块试件之间的磨损机理进行公式推导,再根据赫兹点接触理论,推导出磨损体积与时间的关系式,建立动载工况下牙轮钻头螺旋密封磨损模型。Step (3): Based on the basic model of spiral seal wear of roller drill bits under dynamic load conditions, the wear mechanism between ball block specimens under dynamic load is improved, and then the relationship between wear volume and time is derived based on Hertz point contact theory to establish a spiral seal wear model of roller drill bits under dynamic load conditions.
可选的,步骤(1)中,在球-块摩擦实验中采用控制变量法,测得不同温度和不同载荷下摩擦副之间的摩擦系数变化情况,以及磨损后块试件的磨损量和表面形貌,得出密封结构摩擦副材料在不同温度和不同载荷下的摩擦磨损性能。Optionally, in step (1), a controlled variable method is used in the ball-block friction experiment to measure the change in the friction coefficient between the friction pairs at different temperatures and different loads, as well as the wear amount and surface morphology of the block specimen after wear, to obtain the friction and wear performance of the sealing structure friction pair material at different temperatures and different loads.
可选的,步骤(1)中,高温工况下牙轮钻头螺旋密封磨损的基础模型和动载工况下牙轮钻头螺旋密封磨损的基础模型均为Archard磨损模型。Optionally, in step (1), the basic model of spiral seal wear of the roller drill bit under high temperature conditions and the basic model of spiral seal wear of the roller drill bit under dynamic load conditions are both Archard wear models.
可选的,步骤(2)中,温度影响参数包括:温度系数,压力指数和速度指数。Optionally, in step (2), the temperature influencing parameters include: temperature coefficient, pressure index and velocity index.
基于Archard磨损模型,引入温度系数,压力指数和速度指数进行修正,得到高温工况下牙轮钻头螺旋密封磨损的初步模型。Based on the Archard wear model, the temperature coefficient, pressure index and speed index were introduced for correction, and a preliminary model of spiral seal wear of roller drill bits under high temperature conditions was obtained.
高温工况下牙轮钻头螺旋密封磨损的初步模型如下:The preliminary model of spiral seal wear of roller bit under high temperature conditions is as follows:
其中,ks为温度系数;T为温度;T0为室温,取25℃;m为压力指数;n为速度指数;v为滑动速度,单位为mm/s;为深度磨损率;P为接触磨损区域的压力,单位为MPa;k'为磨损系数。Where, k s is the temperature coefficient; T is the temperature; T 0 is the room temperature, which is 25°C; m is the pressure index; n is the velocity index; v is the sliding velocity, in mm/s; is the depth wear rate; P is the pressure in the contact wear area, in MPa; k' is the wear coefficient.
可选的,步骤(2)中,对高温工况下牙轮钻头螺旋密封磨损的初步模型进行拟合,输入球-块摩擦实验的数据得到最优解,建立高温工况下牙轮钻头螺旋密封磨损模型,具体为:将高温工况下牙轮钻头螺旋密封磨损的初步模型输入1stOpt软件中并定义求解参数和取值范围,包括温度系数ks、磨损系数k'、压力指数m和速度指数n,得到拟合代码,输入球-块摩擦实验的数据至拟合代码中,得到最优解,建立高温工况下牙轮钻头螺旋密封磨损模型。Optionally, in step (2), a preliminary model of spiral seal wear of a roller drill bit under high temperature conditions is fitted, data from a ball-block friction experiment is input to obtain an optimal solution, and a spiral seal wear model of a roller drill bit under high temperature conditions is established. Specifically, the preliminary model of spiral seal wear of a roller drill bit under high temperature conditions is input into 1stOpt software and solution parameters and value ranges are defined, including temperature coefficient k s , wear coefficient k', pressure index m and speed index n, to obtain a fitting code, data from a ball-block friction experiment is input into the fitting code, an optimal solution is obtained, and a spiral seal wear model of a roller drill bit under high temperature conditions is established.
高温工况下牙轮钻头螺旋密封磨损模型如下:The wear model of the spiral seal of the roller drill bit under high temperature conditions is as follows:
其中,为深度磨损率;T为温度;T0为室温,取25℃;P为接触磨损区域的压力,单位为MPa;v为滑动速度,单位为mm/s。in, is the depth wear rate; T is the temperature; T0 is the room temperature, which is 25°C; P is the pressure in the contact wear area, in MPa; v is the sliding velocity, in mm/s.
可选的,步骤(3)中,基于Archard磨损模型进行改进,对动态载荷下的球块试件之间的磨损机理进行公式推导,得出:Optionally, in step (3), based on the Archard wear model, the wear mechanism between the ball block specimens under dynamic load is improved and the formula is derived to obtain:
其中,dh是磨损深度微元,单位为mm;k1是改进方程中的磨损系数为无量纲常数;v是球、块试件的相对滑动速度,单位为m/s;σ0(t)是接触点处接触应力随时间的变化函数,根据接触原理可将表达式改写为σ0(t)=σsS(t),其中S(t)接触面积函数其大小和动态载荷幅值变化有关;Hm是块试件材料的硬度参数,是材料本身的属性。Among them, dh is the wear depth element, in mm; k1 is the wear coefficient in the improved equation, which is a dimensionless constant; v is the relative sliding velocity of the ball and block specimen, in m/s; σ0 (t) is the function of the change of contact stress at the contact point with time. According to the contact principle, the expression can be rewritten as σ0 (t)= σsS (t), where S(t) is the contact area function and its size is related to the change of dynamic load amplitude; Hm is the hardness parameter of the block specimen material, which is a property of the material itself.
对基于Archard磨损模型的改进式求时间t的积分,k1、v、Hm为常数,得出:The improved formula based on Archard wear model is integrated over time t, with k 1 , v, and Hm as constants, and we get:
可选的,步骤(3)中,在对基于Archard磨损模型的改进式求时间t的积分后,根据赫兹点接触理论,推导出磨损体积与时间的关系,将被磨损体积V0写为t时刻下的磨损体积微元dV0的关系式,如下:Optionally, in step (3), after integrating the improved formula based on Archard wear model with respect to time t, the relationship between the wear volume and time is derived according to Hertz point contact theory, and the wear volume V 0 is written as the relationship of the wear volume microelement dV 0 at time t, as follows:
其中,K为引入的磨损系数;l为摩擦副相对滑动距离;E*为两物体的等效弹性模量;r为球试件半径;θ为切入半角,θ的取值范围为(0-π/8)之间。Among them, K is the introduced wear coefficient; l is the relative sliding distance of the friction pair; E * is the equivalent elastic modulus of the two objects; r is the radius of the ball specimen; θ is the cutting half angle, and the value range of θ is between (0-π/8).
球试件受到的法向载荷为正弦型动态载荷,动态载荷随时间的变化函数为:F(t)=x0+10sinωt,对将被磨损体积V0写为t时刻下的磨损体积微元dV0的关系式求积分,得到动态载荷下磨损时长为t0时的磨损体积改进公式,建立动载工况下牙轮钻头螺旋密封磨损模型如下:The normal load on the ball specimen is a sinusoidal dynamic load. The function of the dynamic load changing with time is: F(t) = x 0 +10sinωt. The relationship between the wear volume V 0 and the wear volume microelement dV 0 at time t is integrated to obtain the improved formula of the wear volume when the wear time is t 0 under dynamic load. The wear model of the spiral seal of the cone drill bit under dynamic load conditions is established as follows:
其中,K为引入的磨损系数,l为摩擦副相对滑动距离,E*为两物体的等效弹性模量,r为球试件半径,θ为切入半角。Among them, K is the introduced wear coefficient, l is the relative sliding distance of the friction pair, E * is the equivalent elastic modulus of the two objects, r is the radius of the ball specimen, and θ is the cutting half angle.
可选的,还包括:在建立高温工况下牙轮钻头螺旋密封磨损模型和动载工况下牙轮钻头螺旋密封磨损模型后,对高温工况下牙轮钻头螺旋密封磨损模型和动载工况下牙轮钻头螺旋密封磨损模型进行仿真验证。Optionally, it also includes: after establishing the spiral seal wear model of the roller drill bit under high temperature conditions and the spiral seal wear model of the roller drill bit under dynamic load conditions, simulation verification is performed on the spiral seal wear model of the roller drill bit under high temperature conditions and the spiral seal wear model of the roller drill bit under dynamic load conditions.
对高温工况下牙轮钻头螺旋密封磨损模型进行仿真验证具体为:The simulation verification of the spiral seal wear model of the cone drill bit under high temperature conditions is as follows:
采用有限元软件进行静态接触仿真,得到加载过程时刻摩擦副的接触应力分布情况。Finite element software is used to perform static contact simulation to obtain the contact stress distribution of the friction pair during the loading process.
采用动态摩擦生热仿真,得出在试验工况下球块之间的摩擦生热量并不明显,忽略摩擦生热对温度的影响。By using dynamic frictional heat generation simulation, it is found that the frictional heat generation between the balls is not obvious under the test conditions, and the effect of frictional heat generation on temperature can be ignored.
采用动态摩擦磨损仿真,运用Abaqus有限元仿真软件,导入二次开发用户子程序Umeshmotion,调用自适应网格划分技术实现模型因摩擦而损失材料的过程,得到不同温度和载荷下磨损体积变化情况,并与试验结果作对比,得到较准确的磨损体积。Dynamic friction and wear simulation was adopted, Abaqus finite element simulation software was used, the secondary development user subroutine Umeshmotion was imported, and the adaptive meshing technology was called to realize the process of material loss due to friction in the model. The wear volume change under different temperatures and loads was obtained, and compared with the test results to obtain a more accurate wear volume.
通过磨损体积,据几何关系计算磨损对螺旋圈槽深的改变量,计算最优槽深值。The change in the spiral groove depth caused by wear is calculated based on the wear volume and geometric relationship, and the optimal groove depth value is calculated.
对动载工况下牙轮钻头螺旋密封磨损模型进行仿真验证具体为:The simulation verification of the wear model of the spiral seal of the cone drill bit under dynamic load conditions is as follows:
通过有限元仿真模拟动载和静载下的球-块磨损,将仿真结果与试验结果对比。The ball-block wear under dynamic and static loads was simulated by finite element simulation, and the simulation results were compared with the experimental results.
经由上述的技术方案可知,与现有技术相比,本发明提出了高温、动载工况下牙轮钻头螺旋密封磨损模型的建立方法。首先,通过球-块摩擦实验,得到密封结构摩擦副材料在不同温度和不同载荷下的摩擦磨损性能,选择Archard磨损模型作为高温工况下牙轮钻头螺旋密封磨损的基础模型和动载工况下牙轮钻头螺旋密封磨损的基础模型;然后,通过在Archard磨损模型的基础上引入温度影响参数,并进行拟合,输入球-块摩擦实验的数据得到最优解,建立高温工况下牙轮钻头螺旋密封磨损模型;通过对Archard磨损模型进行改进,对动态载荷下的球块试件之间的磨损机理进行公式推导,再根据赫兹点接触理论,推导出磨损体积与时间的关系式,建立动载工况下牙轮钻头螺旋密封磨损模型;最后,对高温工况下牙轮钻头螺旋密封磨损模型进行仿真验证,包括:静态接触仿真,动态摩擦生热仿真,动态摩擦磨损仿真,得到较准确的磨损体积,通过磨损体积,据几何关系计算磨损对螺旋圈槽深的改变量,基于最优槽深值分析,得到当前高温工况基本不会对密封效果造成影响;对动载工况下牙轮钻头螺旋密封磨损模型仿真,采用有限元仿真模拟动载和静载下的球-块磨损,将仿真结果与试验结果对比,验证理论模型和试验的可靠性与正确性,得出了动态载荷峰值对于密封圈螺纹上的接触应力影响极大,且成线性关系;动态载荷峰值越大密封圈磨损也会越来越剧烈;频率对于密封圈表面的接触应力几乎没有影响。上述建立的高温、动载工况下牙轮钻头螺旋密封磨损模型能够准确模拟出高温、动载工况下牙轮钻头螺旋密封的磨损机理,为实际牙轮钻头螺旋密封结构的防磨改进设计方案提供理论支撑,对提高牙轮钻头螺旋密封性能,提高牙轮钻头螺旋密封寿命,降低钻井成本,提高钻井效率具有重要意义。It can be seen from the above technical solutions that, compared with the prior art, the present invention proposes a method for establishing a wear model of the spiral seal of a roller drill bit under high temperature and dynamic load conditions. First, through the ball-block friction experiment, the friction and wear performance of the friction pair material of the sealing structure under different temperatures and different loads is obtained, and the Archard wear model is selected as the basic model of the spiral seal wear of the roller drill bit under high temperature conditions and the basic model of the spiral seal wear of the roller drill bit under dynamic load conditions; then, by introducing temperature influencing parameters on the basis of the Archard wear model and fitting, the data of the ball-block friction experiment is input to obtain the optimal solution, and the wear model of the spiral seal of the roller drill bit under high temperature conditions is established; by improving the Archard wear model, the wear mechanism between the ball-block specimens under dynamic load is derived, and then according to the Hertz point contact theory, the relationship between the wear volume and time is derived, and the wear model of the spiral seal of the roller drill bit under dynamic load conditions is established; finally, the high temperature The wear model of the spiral seal of the roller drill bit under working conditions is simulated and verified, including: static contact simulation, dynamic friction heat generation simulation, dynamic friction wear simulation, to obtain a more accurate wear volume. Through the wear volume, the change in the wear on the spiral groove depth is calculated according to the geometric relationship. Based on the analysis of the optimal groove depth value, it is obtained that the current high temperature working condition will basically not affect the sealing effect; for the simulation of the wear model of the spiral seal of the roller drill bit under dynamic load conditions, finite element simulation is used to simulate the ball-block wear under dynamic and static loads, and the simulation results are compared with the experimental results to verify the reliability and correctness of the theoretical model and the experiment. It is concluded that the dynamic load peak has a great influence on the contact stress on the sealing ring thread, and the relationship is linear; the larger the dynamic load peak, the more severe the wear of the sealing ring will be; the frequency has almost no effect on the contact stress on the surface of the sealing ring. The above-mentioned wear model of the spiral seal of the roller drill bit under high temperature and dynamic load conditions can accurately simulate the wear mechanism of the spiral seal of the roller drill bit under high temperature and dynamic load conditions, and provide theoretical support for the anti-wear improvement design scheme of the spiral seal structure of the actual roller drill bit. It is of great significance to improve the performance of the spiral seal of the roller drill bit, increase the life of the spiral seal of the roller drill bit, reduce drilling costs, and improve drilling efficiency.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying creative work.
图1为本发明的流程示意图。FIG. 1 is a schematic diagram of the process of the present invention.
图2为本发明的Archard磨损模型的粘着磨损过程示意图。FIG. 2 is a schematic diagram of the adhesive wear process of the Archard wear model of the present invention.
图3为本发明的基于Archard磨损模型改进的动态载荷下的球块试件的摩擦副接触区域放大示意图。FIG3 is an enlarged schematic diagram of the contact area of the friction pair of the ball block specimen under dynamic load based on the improved Archard wear model of the present invention.
图4为本发明的通过静态接触仿真得到加载完成时刻摩擦副的接触应力为10N的分布情况示意图。FIG4 is a schematic diagram of the distribution of the contact stress of the friction pair of 10N when loading is completed through static contact simulation according to the present invention.
图5为本发明的通过静态接触仿真得到加载完成时刻摩擦副的接触应力为20N的分布情况示意图。FIG5 is a schematic diagram of the distribution of the contact stress of the friction pair of 20N obtained by static contact simulation at the moment when loading is completed according to the present invention.
图6为本发明的通过静态接触仿真得到加载完成时刻摩擦副的接触应力为30N的分布情况示意图。FIG6 is a schematic diagram of the distribution of the contact stress of the friction pair of 30N obtained by static contact simulation at the moment of loading completion according to the present invention.
图7为本发明的通过静态接触仿真得到加载完成时刻摩擦副的接触应力为40N的分布情况示意图。FIG. 7 is a schematic diagram showing the distribution of the contact stress of the friction pair of 40N when loading is completed through static contact simulation according to the present invention.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.
本发明实施例公开了高温、动载工况下牙轮钻头螺旋密封磨损模型的建立方法,如图1所示,包括:The embodiment of the present invention discloses a method for establishing a wear model of a spiral seal of a roller drill bit under high temperature and dynamic load conditions, as shown in FIG1 , comprising:
步骤(1):通过在球-块摩擦实验中采用控制变量法,分别改变块试件温度和施加的载荷,测得不同温度和不同载荷下摩擦副之间的摩擦系数变化情况,以及磨损后块试件的磨损量和表面形貌,得出密封结构摩擦副材料在不同温度和不同载荷下的摩擦磨损性能,选择Archard磨损模型作为高温工况下牙轮钻头螺旋密封磨损的基础模型和动载工况下牙轮钻头螺旋密封磨损的基础模型,如图2所示。Step (1): By using the control variable method in the ball-block friction experiment, the temperature of the block specimen and the applied load are changed respectively, and the change of the friction coefficient between the friction pairs under different temperatures and different loads, as well as the wear amount and surface morphology of the block specimen after wear are measured. The friction and wear performance of the sealing structure friction pair material under different temperatures and different loads is obtained, and the Archard wear model is selected as the basic model for the wear of the spiral seal of the roller drill bit under high temperature conditions and the basic model for the wear of the spiral seal of the roller drill bit under dynamic load conditions, as shown in Figure 2.
步骤(2)基于Archard磨损模型,引入温度影响参数:温度系数,压力指数和速度指数进行修正,得到高温工况下牙轮钻头螺旋密封磨损的初步模型,具体为:Step (2) is based on the Archard wear model, and the temperature influencing parameters are introduced: temperature coefficient, pressure index and speed index for correction, so as to obtain a preliminary model of spiral seal wear of roller bit under high temperature conditions, which is specifically:
Archard模型通常在局部区域上应用,设磨损区域面积为S1,将Archard模型两边同时除以S1可得:The Archard model is usually applied to a local area. Assuming the area of the wear area is S1, dividing both sides of the Archard model by S1 yields:
其中FN为法向载荷,L为球块试验的滑动距离,W'为总磨损量,k为磨损过程中产生磨屑的概率,P为接触磨损区域的压力(MPa),H为布氏硬度值,磨损区域深度h由式(1)可得。Where F N is the normal load, L is the sliding distance of the ball block test, W' is the total wear, k is the probability of generating wear debris during the wear process, P is the pressure of the contact wear area (MPa), H is the Brinell hardness value, and the wear area depth h can be obtained by formula (1).
联立式(1)和式(2)可得:Combining equation (1) and equation (2), we can get:
式(3)对时间t是可导的,因此对时间t进行求导可得:Formula (3) is differentiable with respect to time t, so taking the derivative with respect to time t yields:
进一步整理式(4),令为深度磨损率,/>为滑动速度(mm/s),将k与H合并为新的磨损系数k',在考虑温度的影响下进行修正,得高温工况下牙轮钻头螺旋密封磨损初步模型,如式(5):Further rearrange equation (4) and let is the depth wear rate, /> is the sliding speed (mm/s), k and H are combined into a new wear coefficient k', and the wear coefficient is corrected under the influence of temperature to obtain the preliminary wear model of the spiral seal of the cone drill bit under high temperature conditions, as shown in formula (5):
其中,ks为温度系数,T0为室温,取25℃,m为压力指数,n为速度指数。Among them, ks is the temperature coefficient, T0 is the room temperature, which is 25℃, m is the pressure index, and n is the velocity index.
对高温工况下牙轮钻头螺旋密封磨损初步模型进行拟合,输入球-块摩擦实验的数据得到最优解,建立高温工况下牙轮钻头螺旋密封磨损模型,具体为:The preliminary model of spiral seal wear of cone drill bit under high temperature conditions was fitted, and the data of ball-block friction experiment was input to obtain the optimal solution. The spiral seal wear model of cone drill bit under high temperature conditions was established, which is as follows:
在1stOpt中用Function语句将所需拟合公式设定为式(5),用Parameters语句定义求解参数和取值范围,包括温度系数ks、磨损系数k'、压力指数m以及速度指数n,按照变量定义的顺序数据实现拟合,拟合代码及注释如下:In 1stOpt, the required fitting formula is set to Formula (5) using the Function statement, and the solution parameters and value ranges are defined using the Parameters statement, including the temperature coefficient k s , the wear coefficient k', the pressure index m, and the speed index n. The fitting is performed according to the sequence data defined by the variables. The fitting code and comments are as follows:
输入球-块摩擦实验的数据至拟合代码中,求解算法设置为通用全局优化算法(UGO),收敛判据为1E-10,最大迭代数为1000,并行数30,控制数50,收敛判断数15,运行后整理得到拟合结果如表1所示:The data of the ball-block friction experiment was input into the fitting code, and the solution algorithm was set to the universal global optimization algorithm (UGO), the convergence criterion was 1E-10, the maximum number of iterations was 1000, the number of parallel operations was 30, the number of controls was 50, and the number of convergence judgments was 15. After running, the fitting results were sorted out as shown in Table 1:
表1Table 1
拟合结果显示R方为0.94,说明拟合结果较好,将求解出的参数代入到式(5)中,最终得到适用于20CrNiMo-40Cr摩擦副的摩擦磨损模型如式(6):The fitting result shows that the R-square is 0.94, indicating that the fitting result is good. Substituting the solved parameters into formula (5), the friction and wear model suitable for the 20CrNiMo-40Cr friction pair is finally obtained as shown in formula (6):
其中,为深度磨损率;T为温度;T0为室温,取25℃;P为接触磨损区域的压力,单位为MPa;v为滑动速度,单位为mm/s。in, is the depth wear rate; T is the temperature; T0 is the room temperature, which is 25°C; P is the pressure in the contact wear area, in MPa; v is the sliding velocity, in mm/s.
步骤(3):基于Archard磨损模型进行改进,对动态载荷下的球块试件之间的磨损机理进行公式推导,在试验过程中球试件与块试件进行相对滑动,球试件的切入深度会不断的加深,而在经典Archard磨损模型并不能直接的将其反应出来。由于摩擦副接触形式为点接触,在进行相对滑动时虽然会在接触点产生大量的摩擦热,但块试件整体尺寸较大,摩擦热逸散极快,因此在对经典Archard磨损模型进行改进时不考虑温度变化,且将摩擦副材料的屈服极限σs和摩擦副的相对滑动速度v视为固定常量,得到以下改进公式:Step (3): Based on the improvement of Archard wear model, the wear mechanism between ball and block specimens under dynamic load is derived. During the test, the ball specimen and the block specimen slide relative to each other, and the penetration depth of the ball specimen will continue to deepen, which cannot be directly reflected in the classic Archard wear model. Since the contact form of the friction pair is point contact, although a large amount of friction heat will be generated at the contact point during relative sliding, the overall size of the block specimen is large and the friction heat dissipates very quickly. Therefore, when improving the classic Archard wear model, temperature changes are not considered, and the yield limit σs of the friction pair material and the relative sliding speed v of the friction pair are regarded as fixed constants, and the following improved formula is obtained:
其中,dh是磨损深度微元,单位为mm;k1是改进方程中的磨损系数为无量纲常数;v是球、块试件的相对滑动速度,单位为m/s;σ0(t)是接触点处接触应力随时间的变化函数,根据接触原理可将表达式改写为σ0(t)=σsS(t),其中S(t)接触面积函数其大小和动态载荷幅值变化有关;Hm是块试件材料的硬度参数,是材料本身的属性。Among them, dh is the wear depth element, in mm; k1 is the wear coefficient in the improved equation, which is a dimensionless constant; v is the relative sliding velocity of the ball and block specimen, in m/s; σ0 (t) is the function of the change of contact stress at the contact point with time. According to the contact principle, the expression can be rewritten as σ0 (t)= σsS (t), where S(t) is the contact area function and its size is related to the change of dynamic load amplitude; Hm is the hardness parameter of the block specimen material, which is a property of the material itself.
对基于Archard磨损模型的改进式求时间t的积分,k1、v、Hm为常数,得出:The improved formula based on Archard wear model is integrated over time t, with k 1 , v, and H m as constants, and we get:
将摩擦副接触区域无限放大如图3所示,将块试件表面视为x-y平面,建立O-xyz空间直角坐标系,且取任意时间的接触状态分析,在法向载荷F的作用下h无限趋近于0时,此时为极限情况下的椭圆赫兹接触。The contact area of the friction pair is infinitely enlarged as shown in Figure 3. The surface of the specimen is regarded as the x-y plane, and the O-xyz space rectangular coordinate system is established. The contact state at any time is analyzed. When h approaches 0 infinitely under the action of the normal load F, it is an elliptical Hertzian contact in the limit case.
根据赫兹点接触理论,推导出磨损体积与时间的关系式,建立动载工况下牙轮钻头螺旋密封磨损模型,具体为:According to Hertz point contact theory, the relationship between wear volume and time is derived, and the wear model of the spiral seal of the cone drill bit under dynamic load conditions is established, which is as follows:
根据赫兹点接触理论,试验中的球试件与块试件之间的接触由于试验材料的弹性形变会形成一个椭圆接触区域,切入试块的面积即为平面受球体挤压产生的弹性形变,可通过几何关系表示为:According to Hertz's point contact theory, the contact between the ball specimen and the block specimen in the test will form an elliptical contact area due to the elastic deformation of the test material. The area cut into the test block is the elastic deformation of the plane caused by the compression of the sphere, which can be expressed by the geometric relationship as follows:
S=θr2-arcosθ#(9);S = θr 2 -arcosθ# (9);
其中,θ为切入半角,a为宽度,r为球试件半径,假设摩擦副相对滑动距离为l,则被磨损体积V0为:Among them, θ is the cutting half angle, a is the width, r is the radius of the ball specimen, and assuming that the relative sliding distance of the friction pair is l, the worn volume V0 is:
V0=Sl=θr2l-arlcosθ#(10);V 0 =Sl=θr 2 l-arlcosθ#(10);
其中,θ为切入半角,a为宽度,r为球试件半径,l为摩擦副相对滑动距离,S为切入试块的面积,根据球体与平面的赫兹接触理论,块试件表面产生弹性位移的点在Z轴负方向的位移δ可表示为:Among them, θ is the cut-in half angle, a is the width, r is the radius of the ball specimen, l is the relative sliding distance of the friction pair, S is the area of the cut-in test block, and according to the Hertz contact theory of sphere and plane, the displacement δ of the point on the surface of the block specimen where elastic displacement occurs in the negative direction of the Z axis can be expressed as:
其中,h为切入试块件的深度,r为球试件半径,δ为位移,赫兹压力在接触区域的分布,以及在z轴方向上产生的垂直位移为:Where h is the depth of the specimen, r is the radius of the spherical specimen, δ is the displacement, the distribution of the Hertzian pressure in the contact area, and the vertical displacement in the z-axis direction are:
其中,σ0为接触中心的压力,r为球试件半径,a为宽度,E*为两物体的等效弹性模量,δ为位移,块试件表面在接触区域内所有发生弹性形变的点其垂直位移都相等,说明产生这样的压力分布时,压痕是由刚性圆柱压头在弹性半空间产生,此时接触区域的合力为:Among them, σ0 is the pressure at the contact center, r is the radius of the spherical specimen, a is the width, E * is the equivalent elastic modulus of the two objects, δ is the displacement, and the vertical displacement of all points on the surface of the block specimen that undergo elastic deformation in the contact area is equal, indicating that when such a pressure distribution is generated, the indentation is generated by a rigid cylindrical indenter in the elastic half space, and the resultant force in the contact area is:
其中,σ0为接触中心的压力,a为宽度,δ为位移,将式(13)代入式(10)中得:Where σ0 is the pressure at the contact center, a is the width, and δ is the displacement. Substituting equation (13) into equation (10), we get:
其中,σ0为接触中心的压力,a为宽度,E*为两物体的等效弹性模量,δ为位移,r为球试件半径,h为切入试块件的深度,在式(15)中对位移产生影响的参数为变量h和a,因此需满足以下要求:Among them, σ0 is the pressure at the contact center, a is the width, E * is the equivalent elastic modulus of the two objects, δ is the displacement, r is the radius of the ball specimen, and h is the depth of the cut into the specimen. In formula (15), the parameters that affect the displacement are variables h and a, so the following requirements must be met:
并且变量h和a还应该符合接触半径条件:And the variables h and a should also meet the contact radius conditions:
a2=rh#(18); a2 = rh#(18);
以及接触区域的最大压力条件:And the maximum pressure conditions in the contact area:
将式(18)和式(19)代入式(14)中得:Substituting equation (18) and equation (19) into equation (14), we obtain:
其中,E*为两物体的等效弹性模量,r为球试件半径,h为切入试块件的深度,根据式(20)和式(19)可推出,赫兹接触中心的压力和法向载荷及接触半径的关系:Where, E * is the equivalent elastic modulus of the two objects, r is the radius of the spherical specimen, and h is the depth of the cut into the specimen. According to equations (20) and (19), the relationship between the pressure at the Hertz contact center and the normal load and the contact radius can be deduced as follows:
其中,动态载荷F是关于时间的函数,E*为两物体的等效弹性模量,r为球试件半径,因此可将式(10)的被磨损体积V0写为t时刻下的磨损体积微元dV0:Among them, the dynamic load F is a function of time, E * is the equivalent elastic modulus of the two objects, and r is the radius of the spherical specimen. Therefore, the worn volume V 0 in equation (10) can be written as the wear volume element dV 0 at time t:
式中,K为引入的磨损系数,l为摩擦副相对滑动距离,E*为两物体的等效弹性模量,r为球试件半径,θ为切入半角,因为摩擦副在发生磨损时其接触面积上不是所有的粗糙峰和微凸起都会导致磨粒的产生;随着磨损深度的变化θ的取值是关于时间的函数,但θ的取值范围经过实际测量得出(0-π/8)之间。Where K is the introduced wear coefficient, l is the relative sliding distance of the friction pair, E * is the equivalent elastic modulus of the two objects, r is the radius of the ball specimen, and θ is the cut-in half angle, because not all rough peaks and micro-protrusions on the contact area of the friction pair will lead to the generation of abrasive particles when wear occurs; the value of θ is a function of time as the wear depth changes, but the value range of θ is obtained through actual measurement and is between (0-π/8).
由于球试件受到的法向载荷为正弦型动态载荷,动态载荷随时间的变化函数为:F(t)=x0+10sinωt,参考改进的Archard磨损模型式,对式(22)求积分,可将式子改写为本文动态载荷下磨损时长为t0时的磨损体积改进公式,建立动载工况下牙轮钻头螺旋密封磨损模型如下:Since the normal load on the ball specimen is a sinusoidal dynamic load, the function of the dynamic load changing with time is: F(t) = x 0 +10sinωt. Referring to the improved Archard wear model, the integral of equation (22) can be rewritten as the improved wear volume formula when the wear time is t 0 under dynamic load in this paper. The wear model of the spiral seal of the cone drill bit under dynamic load conditions is established as follows:
其中,K为引入的磨损系数,l为摩擦副相对滑动距离,E*为两物体的等效弹性模量,r为球试件半径,θ为切入半角,结合Archard磨损模型式即可用于ABAQUS仿真中。Among them, K is the introduced wear coefficient, l is the relative sliding distance of the friction pair, E * is the equivalent elastic modulus of the two objects, r is the radius of the ball specimen, and θ is the cutting half angle. Combined with the Archard wear model, it can be used in ABAQUS simulation.
还包括,对高温工况下牙轮钻头螺旋密封磨损模型和动载工况下牙轮钻头螺旋密封磨损模型进行仿真验证具体为:It also includes simulation verification of the wear model of the spiral seal of the cone drill bit under high temperature conditions and the wear model of the spiral seal of the cone drill bit under dynamic load conditions. Specifically:
高温工况下牙轮钻头螺旋密封磨损模型仿真具体为:The simulation of the wear model of the spiral seal of the roller drill bit under high temperature conditions is as follows:
运用Abaqus等有限元仿真软件,进行静态接触仿真,通过建模、设定相互作用、分析步及加载、网格划分及边界条件设置、计算及后处理及仿真结果及分析等步骤,得到加载完成时刻摩擦副的接触应力分布情况,如图4-图7所示。Abaqus and other finite element simulation software are used to perform static contact simulation. Through the steps of modeling, setting interaction, analysis step and loading, meshing and boundary condition setting, calculation and post-processing, simulation results and analysis, the contact stress distribution of the friction pair at the moment of loading is obtained, as shown in Figures 4 to 7.
由于温度有两个来源:一个是热流道系统按照设定温度给系统加温,此温度为已知量,在仿真中以预定义场的方式施加到球块模型上;另一个则是因摩擦产生的热量使温度上升,这部分温度需要通过仿真求解得到,因此采用动态摩擦生热仿真,摩擦生热仿真基本步骤大致与静态接触相同,最终得出在试验工况下球块之间的摩擦生热量并不明显,可以忽略摩擦生热对温度的影响。Since temperature has two sources: one is that the hot runner system heats the system according to the set temperature. This temperature is a known quantity and is applied to the ball block model in the form of a predefined field in the simulation; the other is that the heat generated by friction causes the temperature to rise. This part of the temperature needs to be solved through simulation. Therefore, dynamic friction heat generation simulation is adopted. The basic steps of friction heat generation simulation are roughly the same as those of static contact. Finally, it is concluded that the friction heat generation between the balls is not obvious under the test conditions, and the influence of friction heat generation on temperature can be ignored.
采用动态摩擦磨损仿真,运用Abaqus有限元仿真软件,导入二次开发用户子程序Umeshmotion如下:Dynamic friction and wear simulation is adopted, Abaqus finite element simulation software is used, and the secondary development user subroutine Umeshmotion is imported as follows:
调用自适应网格划分技术实现模型因摩擦而损失材料的过程,得到不同温度和载荷下磨损体积变化情况,并与试验结果作对比,得到较准确的磨损体积。Adaptive meshing technology is used to realize the process of material loss due to friction in the model, and the change of wear volume under different temperatures and loads is obtained. The wear volume is compared with the test results to obtain a more accurate wear volume.
高温工况下牙轮钻头螺旋密封磨损模型仿真研究通过摩擦仿真和磨损仿真得到的磨损体积,据几何关系计算磨损对螺旋圈槽深的改变量,基于最优槽深值分析,得到当前高温工况基本不会对密封效果造成影响。Simulation study on the wear model of spiral seal of roller drill bit under high temperature conditions. The wear volume is obtained through friction simulation and wear simulation, and the change of wear on the spiral groove depth is calculated according to the geometric relationship. Based on the analysis of the optimal groove depth value, it is found that the current high temperature condition will basically not affect the sealing effect.
动载工况下牙轮钻头螺旋密封磨损模型仿真具体为:The simulation of the wear model of the spiral seal of the cone drill bit under dynamic load conditions is as follows:
通过有限元仿真模拟动载和静载下的球-块磨损,将仿真结果与试验结果对比,得出动态载荷下的磨损体积与试验值具有很好的一致性,不同峰值载荷下和不同频率下的磨损体积变化趋势与试验所得结果相同,验证理论模型和试验的可靠性与正确性。Finite element simulation was used to simulate the ball-block wear under dynamic and static loads. The simulation results were compared with the experimental results. It was found that the wear volume under dynamic load was very consistent with the experimental value. The wear volume change trend under different peak loads and different frequencies was the same as the experimental results, verifying the reliability and correctness of the theoretical model and experiment.
动载工况下牙轮钻头螺旋密封磨损模型仿真研究得出了动态载荷峰值对于密封圈螺纹上的接触应力影响极大,且成线性关系;动态载荷峰值越大密封圈磨损也会越来越剧烈;频率对于密封圈表面的接触应力几乎没有影响。The simulation study of the wear model of the spiral seal of a roller drill under dynamic load conditions showed that the dynamic load peak has a great influence on the contact stress on the sealing ring thread, and the relationship is linear; the larger the dynamic load peak, the more severe the wear of the sealing ring will be; and the frequency has almost no effect on the contact stress on the sealing ring surface.
本发明公开了高温、动载工况下牙轮钻头螺旋密封磨损模型的建立方法。首先,通过球-块摩擦实验,得到密封结构摩擦副材料在不同温度和不同载荷下的摩擦磨损性能,选择Archard磨损模型作为高温工况下牙轮钻头螺旋密封磨损的基础模型和动载工况下牙轮钻头螺旋密封磨损的基础模型;然后,通过在Archard磨损模型的基础上引入温度影响参数,并进行拟合,输入球-块摩擦实验的数据得到最优解,建立高温工况下牙轮钻头螺旋密封磨损模型;通过对Archard磨损模型进行改进,对动态载荷下的球块试件之间的磨损机理进行公式推导,再根据赫兹点接触理论,推导出磨损体积与时间的关系式,建立动载工况下牙轮钻头螺旋密封磨损模型;最后,对高温工况下牙轮钻头螺旋密封磨损模型进行仿真验证,包括:静态接触仿真,动态摩擦生热仿真,动态摩擦磨损仿真,得到较准确的磨损体积,通过磨损体积,据几何关系计算磨损对螺旋圈槽深的改变量,基于最优槽深值分析,得到当前高温工况基本不会对密封效果造成影响;对动载工况下牙轮钻头螺旋密封磨损模型仿真,采用有限元仿真模拟动载和静载下的球-块磨损,将仿真结果与试验结果对比,验证理论模型和试验的可靠性与正确性,得出了动态载荷峰值对于密封圈螺纹上的接触应力影响极大,且成线性关系;动态载荷峰值越大密封圈磨损也会越来越剧烈;频率对于密封圈表面的接触应力几乎没有影响。上述建立的高温、动载工况下牙轮钻头螺旋密封磨损模型能够准确模拟出高温、动载工况下牙轮钻头螺旋密封的磨损机理,为实际牙轮钻头螺旋密封结构的防磨改进设计方案提供理论支撑,对提高牙轮钻头螺旋密封性能,提高牙轮钻头螺旋密封寿命,降低钻井成本,提高钻井效率具有重要意义。The present invention discloses a method for establishing a wear model of a spiral seal of a roller drill bit under high temperature and dynamic load conditions. First, through a ball-block friction experiment, the friction and wear performance of the friction pair material of the sealing structure under different temperatures and different loads is obtained, and the Archard wear model is selected as the basic model of the spiral seal wear of the roller drill bit under high temperature conditions and the basic model of the spiral seal wear of the roller drill bit under dynamic load conditions; then, by introducing temperature influencing parameters on the basis of the Archard wear model and fitting, the data of the ball-block friction experiment is input to obtain the optimal solution, and the wear model of the spiral seal of the roller drill bit under high temperature conditions is established; by improving the Archard wear model, the wear mechanism between the ball-block specimens under dynamic load is derived, and then according to the Hertz point contact theory, the relationship between the wear volume and time is derived, and the wear model of the spiral seal of the roller drill bit under dynamic load conditions is established; finally, the high temperature The wear model of the spiral seal of the roller drill bit under working conditions is simulated and verified, including: static contact simulation, dynamic friction heat generation simulation, dynamic friction wear simulation, to obtain a more accurate wear volume. Through the wear volume, the change in the wear on the spiral groove depth is calculated according to the geometric relationship. Based on the analysis of the optimal groove depth value, it is obtained that the current high temperature working condition will basically not affect the sealing effect; for the simulation of the wear model of the spiral seal of the roller drill bit under dynamic load conditions, finite element simulation is used to simulate the ball-block wear under dynamic and static loads, and the simulation results are compared with the experimental results to verify the reliability and correctness of the theoretical model and the experiment. It is concluded that the dynamic load peak has a great influence on the contact stress on the sealing ring thread, and the relationship is linear; the larger the dynamic load peak, the more severe the wear of the sealing ring will be; the frequency has almost no effect on the contact stress on the surface of the sealing ring. The above-mentioned wear model of the spiral seal of the roller drill bit under high temperature and dynamic load conditions can accurately simulate the wear mechanism of the spiral seal of the roller drill bit under high temperature and dynamic load conditions, and provide theoretical support for the anti-wear improvement design scheme of the spiral seal structure of the actual roller drill bit. It is of great significance to improve the performance of the spiral seal of the roller drill bit, increase the life of the spiral seal of the roller drill bit, reduce drilling costs, and improve drilling efficiency.
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。In this specification, each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the device disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part.
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。The above description of the disclosed embodiments enables one skilled in the art to implement or use the present invention. Various modifications to these embodiments will be apparent to one skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present invention. Therefore, the present invention will not be limited to the embodiments shown herein, but rather to the widest scope consistent with the principles and novel features disclosed herein.
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210635214.0A CN114936436B (en) | 2022-06-07 | 2022-06-07 | Method for establishing spiral seal abrasion model of roller bit under high-temperature and dynamic load working conditions |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210635214.0A CN114936436B (en) | 2022-06-07 | 2022-06-07 | Method for establishing spiral seal abrasion model of roller bit under high-temperature and dynamic load working conditions |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114936436A CN114936436A (en) | 2022-08-23 |
CN114936436B true CN114936436B (en) | 2024-06-11 |
Family
ID=82866811
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210635214.0A Active CN114936436B (en) | 2022-06-07 | 2022-06-07 | Method for establishing spiral seal abrasion model of roller bit under high-temperature and dynamic load working conditions |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114936436B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116579236B (en) * | 2023-05-06 | 2024-12-10 | 上海交通大学 | Friction and wear test method and system for typical moving pairs of aviation tail nozzles under high temperature and heavy load conditions |
CN117007429B (en) * | 2023-08-07 | 2024-03-22 | 长沙理工大学 | Test method for dynamic performance of subgrade soil considering static deflection stress level |
CN117521284B (en) * | 2023-11-07 | 2024-04-30 | 兰州理工大学 | Construction method of accelerated life prediction model for friction and wear of high-temperature floating ring seals |
CN117828946A (en) * | 2024-01-05 | 2024-04-05 | 武汉理工大学 | Bucket wear resistance evaluation method, device, equipment and storage medium |
CN118036380B (en) * | 2024-02-21 | 2024-09-27 | 江苏科技大学 | Cam profile rolling wear volume simulation method based on multi-subroutine combination |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107740679A (en) * | 2017-10-06 | 2018-02-27 | 西南石油大学 | A kind of coiled tubing drilling combination sealing drill bit |
CN109033590A (en) * | 2018-07-13 | 2018-12-18 | 北京工业大学 | A kind of process parameter optimizing method based on stamp forging mould abrasion |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110127311A1 (en) * | 2009-11-02 | 2011-06-02 | Jeremy Peterson | Out of position friction stir welding of casing and small diameter tubing or pipe |
CN105386972B (en) * | 2015-12-09 | 2017-05-17 | 合肥工业大学 | Screw vacuum pump of motive seal structure |
CN106989915B (en) * | 2017-05-31 | 2020-02-11 | 西南石油大学 | Vertical high-speed roller bit bearing spiral sealing test device |
CN107091736A (en) * | 2017-05-31 | 2017-08-25 | 西南石油大学 | A kind of method of testing of rock bit helixseal sediment outflow performance |
CN106989874B (en) * | 2017-05-31 | 2020-02-07 | 西南石油大学 | Horizontal high-speed roller bit bearing spiral seal test device |
CN107542408A (en) * | 2017-09-18 | 2018-01-05 | 杨龙飞 | Underground liquid drives cumulative drilling tool and application method |
-
2022
- 2022-06-07 CN CN202210635214.0A patent/CN114936436B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107740679A (en) * | 2017-10-06 | 2018-02-27 | 西南石油大学 | A kind of coiled tubing drilling combination sealing drill bit |
CN109033590A (en) * | 2018-07-13 | 2018-12-18 | 北京工业大学 | A kind of process parameter optimizing method based on stamp forging mould abrasion |
Non-Patent Citations (2)
Title |
---|
Sand removal mechanism of a high-speed roller bit with helical sealing;Yi Zhou 等;Section Mechanical Engineering;20190912;第9卷(第18期);1-15 * |
稠油开采中常规螺杆泵定子衬套磨损研究;韩传军 等;润滑与密封;20180515;第43卷(第05期);25-29 * |
Also Published As
Publication number | Publication date |
---|---|
CN114936436A (en) | 2022-08-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114936436B (en) | Method for establishing spiral seal abrasion model of roller bit under high-temperature and dynamic load working conditions | |
Anderson et al. | Experimental and numerical investigations of single abrasive-grain cutting | |
Nelias et al. | Elastic-plastic contact between rough surfaces: proposal for a wear or running-in model | |
Huangfu et al. | Investigation on meshing and dynamic characteristics of spur gears with tip relief under wear fault | |
Rüttimann et al. | Simulation of hexa-octahedral diamond grain cutting tests using the SPH method | |
Afrasiabi et al. | GPU-accelerated meshfree simulations for parameter identification of a friction model in metal machining | |
Agmell et al. | The link between plasticity parameters and process parameters in orthogonal cutting | |
Dehmani et al. | Numerical study of residual stress induced by multi-steps orthogonal cutting | |
Ghaedi et al. | Flexible foundation effect on seismic analysis of Roller Compacted Concrete (RCC) dams using finite element method | |
CN106844849A (en) | A kind of numerical analysis method of soil body orthotropic constitutive model | |
Rodriguez et al. | Implicit or explicit time integration schemes in the PFEM modeling of metal cutting processes | |
CN110263476A (en) | A kind of diamond abrasive grain abrasion wear prediction method based on finite element dynamic delineation emulation | |
Wen et al. | A method for calculating the real contact area of blade-disk structure combined with truncation theory between rough surfaces | |
Hu et al. | Development of BB model and investigation of P-wave propagation across jointed rock masses using CDEM | |
Bose et al. | Finite element sliding wear simulation of 2D steel-on-steel pin-on-disc tribometer | |
Karpat et al. | A novel method for calculation gear tooth stiffness for dynamic analysis of spur gears with asymmetric teeth | |
CN118070618A (en) | A constitutive model reverse identification method and system for polycrystalline diamond sintering process | |
Lin et al. | An elliptical elastic-plastic microcontact model developed for an ellipsoid in contact with a smooth rigid flat | |
Armand et al. | Numerical simulation of partial slip contact using a semi-analytical method | |
CN109359358B (en) | A side groove geometry optimization design method for compact tensile specimens with side grooves | |
CN115406776A (en) | Construction Method of Nonlinear Unified Constitutive Model for Shear Failure of Rock Structural Plane | |
Singh et al. | Analysis of wear phenomena in sliding contact surfaces | |
CN111209698B (en) | A Prediction Method of Frictional Vibration Noise Considering Interface Uncertainty and Time Variation | |
Sanjay et al. | Analysis of crack Initiation in fretting fatigue specimen | |
Ye et al. | Calculation of residual stress in ultrasonic vibration assisted grinding considering thermal-mechanical coupling: a numerical-analytical hybrid prediction approach |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |