CN111723536B - 一种瓦斯抽采钻机***的多体动力学分析方法 - Google Patents
一种瓦斯抽采钻机***的多体动力学分析方法 Download PDFInfo
- Publication number
- CN111723536B CN111723536B CN202010550055.5A CN202010550055A CN111723536B CN 111723536 B CN111723536 B CN 111723536B CN 202010550055 A CN202010550055 A CN 202010550055A CN 111723536 B CN111723536 B CN 111723536B
- Authority
- CN
- China
- Prior art keywords
- generalized
- equation
- independent
- speed
- coordinates
- 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
- 238000005553 drilling Methods 0.000 title claims abstract description 52
- 238000004458 analytical method Methods 0.000 title claims abstract description 50
- 238000000605 extraction Methods 0.000 title claims abstract description 39
- 230000003068 static effect Effects 0.000 claims abstract description 30
- 230000001133 acceleration Effects 0.000 claims abstract description 26
- 238000000034 method Methods 0.000 claims abstract description 25
- 230000008569 process Effects 0.000 claims abstract description 7
- 238000006467 substitution reaction Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 37
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000002940 Newton-Raphson method Methods 0.000 claims description 6
- 238000012933 kinetic analysis Methods 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 3
- 230000010354 integration Effects 0.000 abstract description 2
- 238000006073 displacement reaction Methods 0.000 description 5
- 239000003245 coal Substances 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- 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/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B44/00—Automatic control systems specially adapted for drilling operations, i.e. self-operating systems which function to carry out or modify a drilling operation without intervention of a human operator, e.g. computer-controlled drilling systems; Systems specially adapted for monitoring a plurality of drilling variables or conditions
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- General Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种瓦斯抽采钻机***的多体动力学分析方法,包括:利用牛顿‑拉夫逊Newton‑Raphson和差分法求出静力学方程的近似解;利用代入法求解运动学方程;利用直接积分和构造状态方程法求解动力学方程,并通过动力学分析得到每一时刻的广义坐标、广义速度、广义加速度,从而得到钻机钻进过程中的一系列动态响应。与现有技术相比,根据本发明的分析方法可以预先调整钻进参数来降低该时刻钻孔失效的风险,从而保证钻进过程中钻孔的稳定性。
Description
技术领域
本发明涉及刻钻孔失效的风险预防技术领域,尤其涉及一种瓦斯抽采钻机***的多体动力学分析方法。
背景技术
煤炭资源是世界主要能源之一。近年来,数次煤炭开采工程所引发的瓦斯***灾难引发了广泛关注。有效的抽采煤层瓦斯可以将瓦斯含量和压力降低到安全水平,从根本上消除瓦斯灾害。
近年来,随着煤层勘探开发规模的扩大,瓦斯抽采的成孔率问题受到了国内外学者的持续关注。但是,由于钻孔失效是由多个因素相互影响、相互促进的复杂的***性问题,涉及到摩擦学、振动学、运动学、材料科学、化学、传热学等一系列学科,存在许多不确定性和随机性,因此,关于瓦斯抽采钻进方面的研究工作还需要进一步完善。
发明内容
针对上述缺陷或不足,本发明的目的在于提供一种瓦斯抽采钻机***的多体动力学分析方法,使瓦斯抽采技术的优化或控制更容易实现,为瓦斯抽采向智能化发展提供数据支持。
为达到以上目的,本发明的技术方案为:
一种瓦斯抽采钻机***的多体动力学分析方法,包括以下步骤:
S1、静力学分析:构建瓦斯抽采钻机***的静力学方程,然后利用牛顿-拉夫逊Newton-Raphson法和差分法求出静力学方程的近似解,并通过静力学分析得到广义坐标;
S2、运动学分析:在静力学分析的前提下,构建瓦斯抽采钻机***的运动学方程,然后利用代入法求解运动学方程,并通过运动学分析得到广义速度;
S3、动力学分析:构建瓦斯抽采钻机***的动力学方程,利用直接积分和构造状态方程法求解动力学方程,并通过动力学分析得到每一时刻的广义坐标、广义速度、广义加速度,从而得到钻机钻进过程中的一系列动态响应。
具体地,所述步骤S1具体包括:
当***静止时,构建的瓦斯抽采钻机***的静力学方程:
设
Re=Qe+Fmz-Kq (2)
将广义坐标写成分块形式:
qi、qd分别表示广义独立坐标和相关坐标,求约束方程对广义坐标的偏导,并写成紧凑形式为:
将其写为独立坐标、相关坐标分块形式为:
Bdi为n×(n-nc)阶矩阵,将式(6)代入式(5),两边乘以虚位移得:
求解式(7)即得到qi;
设
则有:
具体地,所述步骤S2具体包括:
将广义速度分块为:
求约束方程对时间的偏导,并写成分块形式为:
参考静力学分析的计算方法,参考(1)和(3),由式(12)推导出相关速度与独立速度的关系为:
独立速度由静力学分析得出的独立坐标直接求导而得,将其代入式(6)求得相关速度,按照式(11)的形式整合独立速度和相关速度,即得到整个***的初始广义速度。
具体地,所述步骤S3具体包括:
构建瓦斯抽采钻机***的运动学方程:
式中,M和K分别为***的质量矩阵和刚度矩阵;和λ分别为约束的雅可比矩阵和拉格朗日乘子;Qe和Qv分别表示***的广义外力矢量和速度的二次矢量;为相互独立的约束方程组;为与广义坐标q和时间t相关的函数,nc表示约束方程的个数;
将式(13)对时间进行求导,即对约束方程组求二阶导数得:
式(15)整理为:
设:
则式(16)写为:
将式(14)中的动力学方程和式(18)整合成矩阵形式为:
式中:Ctt表示约束方程对时间的二次导数,Cq表示约束方程对广义坐标的一次导数;
设:
式(19)写成紧凑形式为:
求解式(21)得:
将式(22)写为分块形式为:
求解式(23)得加速度矢量和拉格朗日乘子为:
由式(24)知广义加速度与广义坐标、广义速度和时间相关,因此将式(24)写为更加简洁的形式为:
将广义加速度写成相应的分块形式为:
构造出状态方程:
对式(29)直接积分求得下一时刻的独立坐标和独立速度,其中初始时刻的广义坐标、速度和加速度分别由静力学分析、运动学分析和动力学分析得到;
随后,将所求得的独立坐标和独立速度分别代入式(6)和式(13),得该时刻的相关坐标和相关速度;
整合所求得的独立坐标和相关坐标即得到该时刻的整体广义坐标,整合所求得的独立速度和相关速度即得到该时刻的整体广义速度;此时,若运算时间大于所需要的终止时间,则可整合输出全部时刻的整体广义坐标、速度和加速度;若运算时间小于所需要的终止时间,则将该时刻的整体广义坐标代入步骤S2中的运动学分析中,再次对***进行运动学和动力学分析,直至运算时间大于所需要的终止时间。
与现有技术比较,本发明提供了一种瓦斯抽采钻机***的多体动力学分析方法,通过对瓦斯抽采钻机***静力学、运动学和动力学分析,提高瓦斯抽采钻机***的动力学模型求解效率,精确求解结果,使瓦斯抽采技术的优化或控制更容易实现,为瓦斯抽采向智能化发展提供数据支持;另外,通过对瓦斯抽采钻机***刚柔耦合多体动力学模型进行求解,求解结果体现了瓦斯抽采工况的复杂性,为瓦斯抽采工艺的优化提供了方法支持,进而根据本发明的分析方法可以预先调整钻进参数来降低该时刻钻孔失效的风险,从而保证钻进过程中钻孔的稳定性。
附图说明
图1是本发明一种实施例的方法流程图。
图2是本发明一种实施例的初始位置的离散图。
图3是本发明一种实施例的第1根钻杆的Z向位移的结果图。
图4是本发明一种实施例的第1根钻杆的Z向速度的结果图。
图5是本发明一种实施例的第1根钻杆的Z向加速度的结果图。
具体实施方式
下面将结合附图对本发明做详细描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明的保护范围。
为了更加确切验证本发明,本实施例以图2所示的一种瓦斯抽采钻机***为例,柔性体以第1根钻杆为例,刚性体以钻头为例,其余物体的预估值同理可得,钻杆和钻头的初始位置图如2所示,其中,将每个离散柔性物体均用4个节点来标识,将每个刚性物体力学性质均简化到质心上,以第1根钻杆为起始物体,从节点1开始标识,节点1为钻机和第1根钻杆的连接点,节点2、3为第一根钻杆离散后的标识点,节点4为第1根钻杆和第2根钻杆的连接点,以此类推,钻头利用节点0表示。为第i根钻杆以及钻头的物体坐标系,也可以用一组表示位置和方向的数组代替,其中,表示第i根钻杆和钻头的物体坐标系相对于整体坐标系的位置坐标,θi表示第i根钻杆和钻头在整体坐标系中的方位。为第i根钻杆第j个单元的单元坐标系,同样可以用一个数组(ql,qm,qn)来代替,例如(q4,q5,q6)表示节点2的三个弹性位移方向。
在上述如图2所示的瓦斯抽采钻机***的基础上,参照图1所示,本实施例的一种瓦斯抽采钻机***的多体动力学分析方法,其具体步骤如下:
1、预估每个物体初始位置的广义坐标和广义速度。以图2中的第1根钻杆为例进行数值计算,其初始广义坐标和广义速度的预估值可为表示:
***的静力学模型为:
设
Re=Qe+Fmz-Kq (4)
则式(3)的静力学方程可改写为:
广义坐标的分块形式可写为:
式中,qi和qd分别表示独立坐标和相关坐标。
将式(5)的两边分别乘以虚位移δq可得:
式(3)中约束方程对广义坐标求导的紧凑形式为:
将其按独立坐标和相关坐标进行分块可得:
由于式(9)中约的束方程均具有独立性,因此,其雅可比矩阵均为满秩矩阵,即存在逆矩阵,故式(9)可改写为:
令
则式(10)可写为:
考虑到式(12),总体虚位移可表示为:
其中,
式中,Bdi为n×(n-nc)阶矩阵。
将式(13)代入式(7)可得:
考虑到qi的独立性,式(15)可改写为:
因此,将式(17)代入式(16)可得
-ReBdi=0 (18)
将式(4)和式(14)代入式(18),可得
考虑到式(6)将广义坐标写成分块形式,并代入式(12)消去方程组中的相关坐标项,即可确定初始广义独立坐标qi,从而通过式(12)求出初始广义相关坐标qd,整合独立坐标和相关坐标得出初始广义坐标,记为q(0)。
因此,设
则有
由差分法可知,当选用hx作为等步距时,式(21)中的系数矩阵为:
将式(22)的系数矩阵代入式(21)后,即可通过求解式(21)得到初始差值Δq(0)。
分别设εe和εs为方程和解的误差范围,当且|q(i)-q(i-1)|<εs时,迭代终止。当而时,则需要重新确定***的初值,在修改好初值后,重新进行迭代,直至求得满足条件的的解。若该算法是收敛点,则该迭代的收敛阶数为2,即
|q(i+1)-q*|<β|q(i)-q*|2 (24)
式中,β为与i无关的常数;i为迭代的阶数。
式中,Ak为n×n阶的非奇异矩阵。
考虑到n维平行弦法,若选取全部Ak≡A,则式(27)可改写为:
式中,Δq(k)表示第k次迭代的矫正量,通常被称为Newton差。
将已求得的初值q(0)和Δq(0)代入式(31),通过k=0,1,2,…次迭代求解式(3.31)中的近似值q(k+1),即可确定静态***的整体广义坐标q。期间,与单变量***相似,若两次相邻迭代解的误差满足且 j=1,2,…,n,则迭代终止。否则,则需要在是奇异矩阵且的前提下重新估算初值并进行迭代。
2、在静力学分析的前提下,通过对瓦斯抽采钻机***进行运动学分析,来求解***运动的广义速度。此时约束方程组可写为:
将其对时间进行求导可得:
考虑到式(6),将广义速度写成相应的分块形式为:
将式(34)代入式(33)可得:
参考静力学分析的计算方法,按式(9)到(12)的步骤,可由式(35)推导出相关速度与独立速度的关系为:
独立速度可由静力学分析得出的独立坐标直接求导而得,将其代入式(12)可求得相关速度,按照式(34)的形式整合独立速度和相关速度,即可得到整个***的初始广义速度。
3、在静力学分析和运动学分析的前提下,通过对瓦斯抽采钻机***进行动力学分析,来求解***在运动时每一时刻的广义坐标、广义速度和广义加速度。***的动力学模型为:
式中,M和K分别为***的质量矩阵和刚度矩阵;和λ分别为约束的雅可比矩阵和拉格朗日乘子;Qe和Qv分别表示***的广义外力矢量和速度的二次矢量;为相互独立的约束方程组;为与广义坐标q和时间t相关的函数,nc表示约束方程的个数。
将式(36)对时间进行求导,即对约束方程组求二阶导数可得:
式(38)可整理为:
设
则式(39)可写为:
将式(37)中的动力学方程和式(41)整合成矩阵形式为:
式中:Ctt表示约束方程对时间的二次导数,Cq表示约束方程对广义坐标的一次导数。
设:
式(42)可写成紧凑形式为:
求解式(44)可得:
将其写为分块形式为:
求解式(46)可得加速度矢量和拉格朗日乘子为:
由式(47)可知广义加速度与广义坐标、广义速度和时间相关,因此可将式(47)写为更加简洁的形式为:
考虑到式(6)和式(34),将广义加速度写成相应的分块形式为:
因此,可将状态方程构造为:
对式(52)直接积分求得下一时刻的独立坐标和独立速度,其中初始时刻的广义坐标、速度和加速度分别由静力学分析、运动学分析和动力学分析得到。随后,将所求得的独立坐标和独立速度分别代入式(12)和式(36),可得该时刻的相关坐标和相关速度。整合所求得的独立坐标和相关坐标即可得到该时刻的整体广义坐标,整合所求得的独立速度和相关速度即可得到该时刻的整体广义速度。此时,若运算时间大于所需要的终止时间,则可整合输出全部时刻的整体广义坐标、速度和加速度。若运算时间小于所需要的终止时间,则将该时刻的整体广义坐标代入运动学分析中,再次对***进行运动学和动力学分析,直至运算时间大于所需要的终止时间。以第1根钻杆为例,利用MATLAB软件计算了其在下钻过程中的Z向位移、速度和加速度分别如图3、4和5所示。
由图可知,刚柔耦合模型中钻杆的Z向位移受钻杆重力的影响较为严重,而Z向速度和加速度受钻杆重力的影响相对较小。此外,三条曲线对于振动和摩擦的耦合影响均较为敏感,曲线的波动无明显规律,幅值的突变现象时有发生。由牛顿第二定律可知,力和加速度成正比,因此,由加速度曲线可以更直观的看出力的突变发生在哪一时刻,从而通过预先调整钻进参数来降低该时刻钻孔失效的风险,从而保证钻进过程中钻孔的稳定性。
对于本领域技术人员而言,显然能了解到上述具体事实例只是本发明的优选方案,因此本领域的技术人员对本发明中的某些部分所可能作出的改进、变动,体现的仍是本发明的原理,实现的仍是本发明的目的,均属于本发明所保护的范围。
Claims (3)
1.一种瓦斯抽采钻机***的多体动力学分析方法,其特征在于,包括以下步骤:
S1、静力学分析:构建瓦斯抽采钻机***的静力学方程,然后利用牛顿-拉夫逊Newton-Raphson法和差分法求出静力学方程的近似解,并通过静力学分析得到广义坐标;S1具体包括:
当***静止时,构建的瓦斯抽采钻机***的静力学方程:
设
Re=Qe+Fmz-Kq (2)
将广义坐标写成分块形式:
Bdi为n×(n-nc)阶矩阵,将式(6)代入式(5),两边乘以虚位移得:
求解式(7)即得到qi;
设
则有:
S2、运动学分析:在静力学分析的前提下,构建瓦斯抽采钻机***的运动学方程,然后利用代入法求解运动学方程,并通过运动学分析得到广义速度;
S3、动力学分析:构建瓦斯抽采钻机***的动力学方程,利用直接积分和构造状态方程法求解动力学方程,并通过动力学分析得到每一时刻的广义坐标、广义速度、广义加速度,从而得到瓦斯抽采钻机钻进过程中的一系列动态响应。
3.根据权利要求2所述的瓦斯抽采钻机***的多体动力学分析方法,其特征在于,所述步骤S3具体包括:
构建瓦斯抽采钻机***的运动学方程:
式中,M和K分别为***的质量矩阵和刚度矩阵;和λ分别为约束的雅可比矩阵和拉格朗日乘子;Qe和Qv分别表示***的广义外力矢量和速度的二次矢量;为相互独立的约束方程组;为与广义坐标q和时间t相关的函数,nc表示相互独立的约束个数;
将式(13)对时间进行求导,即对约束方程组求二阶导数得:
式(15)整理为:
设:
则式(16)写为:
将式(14)中的动力学方程和式(18)整合成矩阵形式为:
式中:Ctt表示约束方程对时间的二次导数,Cq表示约束方程对广义坐标的一次导数;
设:
式(19)写成紧凑形式为:
求解式(21)得:
将式(22)写为分块形式为:
求解式(23)得加速度矢量和拉格朗日乘子为:
由式(24)知广义加速度与广义坐标、广义速度和时间相关,因此将式(24)写为更加简洁的形式为:
将广义加速度写成相应的分块形式为:
构造出状态方程:
对式(29)直接积分求得下一时刻的广义独立坐标和广义独立速度,其中初始时刻的广义坐标、广义速度和广义加速度分别由静力学分析、运动学分析和动力学分析得到;
随后,将所求得的独立坐标和独立速度分别代入式(6)和式(13),得该时刻的相关坐标和相关速度;
整合所求得的独立坐标和相关坐标即得到该时刻的广义坐标,整合所求得的独立速度和相关速度即得到该时刻广义速度;此时,若运算时间大于所需要的终止时间,则可整合输出全部时刻的广义坐标、速度和加速度;若运算时间小于所需要的终止时间,则将该时刻的广义坐标代入步骤S2中的运动学分析中,再次对***进行运动学和动力学分析,直至运算时间大于所需要的终止时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010550055.5A CN111723536B (zh) | 2020-06-16 | 2020-06-16 | 一种瓦斯抽采钻机***的多体动力学分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010550055.5A CN111723536B (zh) | 2020-06-16 | 2020-06-16 | 一种瓦斯抽采钻机***的多体动力学分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111723536A CN111723536A (zh) | 2020-09-29 |
CN111723536B true CN111723536B (zh) | 2022-10-04 |
Family
ID=72567071
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010550055.5A Active CN111723536B (zh) | 2020-06-16 | 2020-06-16 | 一种瓦斯抽采钻机***的多体动力学分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111723536B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101639681A (zh) * | 2008-07-29 | 2010-02-03 | 深圳市大族激光科技股份有限公司 | 一种电子装备运动机构性能参数优化方法 |
NL2007656C2 (en) * | 2011-10-25 | 2013-05-01 | Cofely Experts B V | A method of and a device and an electronic controller for mitigating stick-slip oscillations in borehole equipment. |
CN104239599A (zh) * | 2014-07-07 | 2014-12-24 | 西安工业大学 | 一种基于多点定位柔性工装***的动力学仿真分析方法 |
CN106695793A (zh) * | 2017-01-18 | 2017-05-24 | 宁波韦尔德斯凯勒智能科技有限公司 | 一种xyzr四轴钻孔机器人主动柔顺控制装置及方法 |
CN107220421A (zh) * | 2017-05-18 | 2017-09-29 | 北京理工大学 | 一种空间复杂柔性结构多体***动力学建模与计算方法 |
CN109543264A (zh) * | 2018-11-12 | 2019-03-29 | 天津理工大学 | 一种基于多维度重构校正的柔性多体机器人建模与求解方法 |
CN110259433A (zh) * | 2019-06-28 | 2019-09-20 | 宝鸡石油机械有限责任公司 | 一种实体钻机数字化监测方法 |
CN111291499A (zh) * | 2020-03-04 | 2020-06-16 | 岭南师范学院 | 一种基于多体动力学的瓦斯抽采钻机建模方法 |
-
2020
- 2020-06-16 CN CN202010550055.5A patent/CN111723536B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101639681A (zh) * | 2008-07-29 | 2010-02-03 | 深圳市大族激光科技股份有限公司 | 一种电子装备运动机构性能参数优化方法 |
NL2007656C2 (en) * | 2011-10-25 | 2013-05-01 | Cofely Experts B V | A method of and a device and an electronic controller for mitigating stick-slip oscillations in borehole equipment. |
CN104239599A (zh) * | 2014-07-07 | 2014-12-24 | 西安工业大学 | 一种基于多点定位柔性工装***的动力学仿真分析方法 |
CN106695793A (zh) * | 2017-01-18 | 2017-05-24 | 宁波韦尔德斯凯勒智能科技有限公司 | 一种xyzr四轴钻孔机器人主动柔顺控制装置及方法 |
CN107220421A (zh) * | 2017-05-18 | 2017-09-29 | 北京理工大学 | 一种空间复杂柔性结构多体***动力学建模与计算方法 |
CN109543264A (zh) * | 2018-11-12 | 2019-03-29 | 天津理工大学 | 一种基于多维度重构校正的柔性多体机器人建模与求解方法 |
CN110259433A (zh) * | 2019-06-28 | 2019-09-20 | 宝鸡石油机械有限责任公司 | 一种实体钻机数字化监测方法 |
CN111291499A (zh) * | 2020-03-04 | 2020-06-16 | 岭南师范学院 | 一种基于多体动力学的瓦斯抽采钻机建模方法 |
Non-Patent Citations (6)
Title |
---|
Dynamics of flexible multibody systems with hybrid uncertain parameters;Wang Zhe et al.;《Mechanism and Machine Theory》;20180330;第121卷;第128-147页 * |
Flexible rod model for the rotation of a drill string in an arbitrary borehole;Alexander K. Belyaev et al.;《Acta Mechanica》;20171117;第229卷;第841-848页 * |
基于广义坐标形式牛顿-欧拉方法的空间并联机构动力学正问题分析;陈根良等;《机械工程学报》;20090715(第07期);第41-48页 * |
旋挖钻机力学特性分析与疲劳寿命研究;李康健;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20160415;第C038-485页 * |
螺旋钻机变幅时机液耦合动力学的键合图建模;胡均平等;《中南大学学报(自然科学版)》;20160226;第47卷(第2期);第495-502页 * |
钻杆自动传送***结构设计与仿真分析;朱吉良;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅰ辑》;20120915;第B019-17页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111723536A (zh) | 2020-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111136633B (zh) | 针对时变时延下柔性主-从机器人***的全状态控制方法 | |
CN106483964B (zh) | 一种基于接触力观测器的机器人柔顺控制方法 | |
CN106767780B (zh) | 基于Chebyshev多项式插值逼近的扩展椭球集员滤波方法 | |
CN112904728A (zh) | 一种基于改进型趋近律的机械臂滑模控制轨迹跟踪方法 | |
CN108628172A (zh) | 一种基于扩张状态观测器的机械臂高精度运动控制方法 | |
Piperno et al. | Design of efficient partitioned procedures for the transient solution of aeroelastic problems | |
Huang et al. | Generalization of Newton-Euler formulation of dynamic equations to nonrigid manipulators | |
CN111159636A (zh) | 基于绝对节点坐标描述的柔性多体***动力学半解析灵敏度分析方法 | |
CN103345155B (zh) | 微陀螺仪的自适应反演控制***及方法 | |
CN111291499B (zh) | 一种基于多体动力学的瓦斯抽采钻机建模方法 | |
CN110069854A (zh) | 多重分辨tweno格式对可压流场问题的模拟方法 | |
Ding et al. | Second order adjoint sensitivity analysis of multibody systems described by differential–algebraic equations | |
CN114707385A (zh) | 基于Krylov子空间的深部地层导热系数三维预测方法及装置 | |
CN111723536B (zh) | 一种瓦斯抽采钻机***的多体动力学分析方法 | |
Farhat et al. | Transient aeroelastic computations using multiple moving frames of reference | |
CN113276114B (zh) | 一种基于终端任务指派的可重构机械臂协同力/运动控制***与方法 | |
Awrejcewicz et al. | Nonlinear coupled problems in dynamics of shells | |
CN111241728B (zh) | 一种欧拉方程的间断伽辽金有限元数值求解方法 | |
CN111581789B (zh) | 一种基于matlab的平流层飞艇升空多物理场耦合的解耦方法 | |
CN112084592A (zh) | 折叠式桁架动力学分析***、方法、装置和存储介质 | |
Zhang et al. | Geometrically exact shell with drilling rotations formulated on the special Euclidean group SE (3) | |
Zhe et al. | Adaptive complementary filtering algorithm for imu based on mems | |
Gupta et al. | Dynamic programming approach to load estimation using optimal sensor placement and model reduction | |
CN114147713B (zh) | 基于自适应神经网络高阶动态滑模的轨迹跟踪控制方法 | |
Vidoni et al. | Comparison of model order reduction techniques for flexible multibody dynamics using an equivalent rigid-link system 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 |