[go: up one dir, main page]

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 PDF

Info

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
Application number
CN202210635214.0A
Other languages
Chinese (zh)
Other versions
CN114936436A (en
Inventor
周已
王儒晨
蒋阳
胡家豪
雷鑫林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN202210635214.0A priority Critical patent/CN114936436B/en
Publication of CN114936436A publication Critical patent/CN114936436A/en
Application granted granted Critical
Publication of CN114936436B publication Critical patent/CN114936436B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/04Ageing analysis or optimisation against ageing
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling 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

The invention discloses a method for establishing a spiral seal abrasion model of a roller bit under high-temperature and dynamic load working conditions, which is applied to the field of abrasion models. Analyzing friction and wear properties of the friction pair material of the sealing structure at different temperatures and loads through a ball-block friction experiment, and selecting a basic model under a high-temperature dynamic load working condition; based on a basic model under a high-temperature working condition, temperature influence parameters are introduced to correct and fit, data of a ball-block friction experiment are input to obtain an optimal solution, and a spiral seal abrasion model of the roller bit under the high-temperature working condition is established; and (3) based on improvement of a basic model under a dynamic load working condition, carrying out formula deduction on a wear mechanism between ball test pieces under a dynamic load, deducing a relation between a wear volume and time according to a Hertz point contact theory, and establishing a spiral seal wear model of the roller bit under the dynamic load working condition. The built model can accurately simulate the abrasion mechanism under the working condition of high temperature and dynamic load, and provides theoretical support for the abrasion-proof improved design scheme.

Description

高温、动载工况下牙轮钻头螺旋密封磨损模型的建立方法Establishment of wear model of spiral seal of cone drill bit under high temperature and dynamic load conditions

技术领域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时刻下的磨损体积微元dV0Among 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)

1. The method for establishing the spiral seal abrasion model of the roller bit under the high-temperature and dynamic load working condition is characterized by comprising the following steps of:
Step (1): obtaining friction and wear properties of a friction pair material of a sealing structure under different temperatures and different loads through a ball-block friction experiment, and selecting a basic model for spiral seal wear of the roller bit under a high-temperature working condition and a basic model for spiral seal wear of the roller bit under a dynamic load working condition;
step (2) based on the basic model of the spiral seal abrasion of the roller bit under the high-temperature working condition, introducing temperature influence parameters to correct, and obtaining a preliminary model of the spiral seal abrasion of the roller bit under the high-temperature working condition; fitting the initial model of the spiral seal abrasion of the roller bit under the high-temperature working condition, inputting data of the ball-block friction experiment to obtain an optimal solution, and establishing the model of the spiral seal abrasion of the roller bit under the high-temperature working condition;
step (3): based on the basic model of spiral seal abrasion of the roller bit under the dynamic load working condition, carrying out formula deduction on the abrasion mechanism between ball test pieces under the dynamic load, deducing the relation between abrasion volume and time according to the Hertz point contact theory, and establishing a spiral seal abrasion model of the roller bit under the dynamic load working condition;
In the step (1), a basic model of spiral seal abrasion of the roller bit under the high-temperature working condition and a basic model of spiral seal abrasion of the roller bit under the dynamic loading working condition are Archard abrasion models;
In step (2), the temperature influencing parameters include: temperature coefficient, pressure index and velocity index;
Based on an Archard abrasion model, introducing the temperature coefficient, and correcting the pressure index and the speed index to obtain a preliminary model of spiral seal abrasion of the roller bit under a high-temperature working condition;
The initial model of the spiral seal abrasion of the roller bit under the high-temperature working condition is as follows:
Wherein k s is the temperature coefficient; t is the temperature; t 0 is room temperature, and 25 ℃; m is the pressure index; n is a velocity index; v is the sliding speed in mm/s; is the depth wear rate; p is the pressure of the contact wear area in MPa; k' is the wear coefficient;
In the step (3), based on the Archard abrasion model, the abrasion mechanism between the ball test pieces under dynamic load is deduced by a formula, and the abrasion mechanism is obtained:
Wherein dh is the wearing depth infinitesimal in mm; k 1 is the improvement in the wear coefficient in the equation as a dimensionless constant; v is the relative sliding speed of the ball and block test pieces, and the unit is m/s; σ 0 (t) is a function of the change in contact stress at the contact point over time, and the expression is rewritten to σ 0(t)=σs S (t) according to the contact principle, where the S (t) contact area function is related to the change in dynamic load amplitude; h m is the hardness parameter of the block test piece material, and is the attribute of the material;
integrating the time t of the improved Archard abrasion model, wherein k 1, v and Hm are constants, and obtaining:
In the step (3), after integrating the improved expression based on the Archard abrasion model to obtain the time t, according to the hertz point contact theory, the relation between the abrasion volume and the time is deduced, and the abrasion volume V 0 is written as the relation of the abrasion volume infinitesimal dV 0 at the time t, as follows:
Wherein K is the introduced wear coefficient; 1 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 test piece; theta is a cutting-in half angle, and the value range of theta is between 0 and pi/8;
The normal load applied to the ball test piece is a sinusoidal dynamic load, and the function of the dynamic load changing along with time is as follows: f (t) =x 0 +10sin ωt, integrating the relation of the wear volume infinitesimal dV 0 when the worn volume V 0 is written as t to obtain a wear volume improvement formula when the wear duration is t 0 under dynamic load, and building a spiral seal wear model of the roller bit under dynamic load working condition as follows:
2. The method for establishing the spiral seal abrasion model of the roller bit under the high-temperature dynamic load working condition according to claim 1, wherein in the step (1), a control variable method is adopted in the ball-block friction experiment, the change condition of friction coefficients between friction pairs under different temperatures and different loads, the abrasion loss and the surface morphology of a worn block test piece are measured, and the friction and abrasion performance of the friction pair material with the seal structure under different temperatures and different loads is obtained.
3. The method for building a spiral seal abrasion model of a roller bit under a high-temperature and dynamic load working condition according to claim 1, wherein in the step (2), a preliminary model of the spiral seal abrasion of the roller bit under the high-temperature working condition is fitted, data of the ball-block friction experiment is input to obtain an optimal solution, and the spiral seal abrasion model of the roller bit under the high-temperature working condition is built specifically as follows: inputting a preliminary model of spiral seal abrasion of the roller bit under the high-temperature working condition into 1st Opt software, defining solving parameters and a value range, wherein the solving parameters comprise the temperature coefficient k s, the abrasion coefficient k', the pressure index m and the speed index n to obtain a fitting code, inputting data of the ball-block friction experiment into the fitting code to obtain an optimal solution, and establishing a spiral seal abrasion model of the roller bit under the high-temperature working condition;
the spiral seal abrasion model of the roller bit under the high-temperature working condition is as follows:
wherein, 25 ℃ is taken.
4. The method for building a spiral seal wear model of a roller bit under high-temperature and dynamic load conditions according to claim 1, further comprising: after the spiral sealing abrasion model of the roller bit under the high-temperature working condition and the spiral sealing abrasion model of the roller bit under the dynamic loading working condition are established, simulation verification is carried out on the spiral sealing abrasion model of the roller bit under the high-temperature working condition and the spiral sealing abrasion model of the roller bit under the dynamic loading working condition;
simulation verification of the spiral seal abrasion model of the roller bit under the high-temperature working condition is specifically as follows:
static contact simulation is carried out by adopting finite element software, so that the contact stress distribution condition of the friction pair at the moment of loading process is obtained;
by adopting dynamic friction heat generation simulation, the friction heat generation quantity between the ball blocks is not obvious under the test working condition, and the influence of the friction heat generation on the temperature is ignored;
adopting dynamic friction and wear simulation, using Abaqus finite element simulation software, importing a secondary development user subroutine Umeshmotion, calling a self-adaptive grid division technology to realize the process of loss of materials of the model due to friction, obtaining the change condition of the wear volume under different temperatures and loads, and comparing with test results to obtain more accurate wear volume;
calculating the change amount of abrasion to the groove depth of the spiral ring according to the geometric relationship through the abrasion volume, and calculating an optimal groove depth value;
Simulation verification of the spiral seal abrasion model of the roller bit under the dynamic load working condition is specifically as follows:
And simulating ball-block abrasion under dynamic load and static load through finite element simulation, and comparing a simulation result with a test result.
CN202210635214.0A 2022-06-07 2022-06-07 Method for establishing spiral seal abrasion model of roller bit under high-temperature and dynamic load working conditions Active CN114936436B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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