一种振荡器的受力与运动分析方法技术

本发明专利技术公开了一种振荡器的受力与运动分析方法，首先分析转子的重力做功与动能、振荡器定子的重力功与动能、振荡器的重力元功与动能、套管内泥浆中振荡器上力系(包括泥浆阻尼力)功，最后得到振荡器的运动微分方程及运动方程。本发明专利技术的有益效果是能够准确定义振荡器的受力与主要运动。

【技术实现步骤摘要】
一种振荡器的受力与运动分析方法
本专利技术属于钻井
，涉及一种振荡器的受力与运动分析方法。
技术介绍
在钻井完井过程中，灌注水泥固实套管的质量将直接影响着完井后套管的工作质量。若固井水泥分布不均匀不密实，就会造成套管受力不均匀而易于发生破损。因此，灌注水泥如何达到均匀密实的效果，一直是完井工作者关心的问题。为此，人们设计了不同的振荡器下入套管中，以图通过振荡器的摆动冲击套管激起套管振动进而促使套管外水泥波动，达到固井水泥分布均匀密实效果。但对此缺乏理论分析验证支撑，资料少见。本文试图通过对振荡器在直井套管中的运动和对套管的作用主要行为，进行模拟分析，探索分析振荡器对套管水泥固实效果的评估方案。振荡器在套管内的泥浆中在偏心转子的周期惯性力激励下运动，为了简化问题而反映其主要运动影响，视振荡器的运动为绕悬挂点的定点运动，但其运动受到套管内壁的限制，并且何时何处受到运动约束，是未知的，构成边界非线性问题。振荡器除了受重力作用外，还受到泥浆运动阻尼力作用。
技术实现思路
本专利技术的目的在于提供一种振荡器的受力与运动分析方法，本专利技术的有益效果是能够准确定义振荡器的受力与主要运动。本专利技术所采用的技术方案是首先分析转子的重力做功与动能、振荡器定子的重力功与动能、振荡器的重力元功与动能、套管内泥浆中振荡器上力系(包括泥浆阻尼力)功，最后得到振荡器的运动微分方程及运动方程。进一步，转子的重力做功与动能如下，设振荡器转子角速度为ω，将所有转子组成部分模拟为扇形截面偏心块，任一瞬时其对称面与x1的夹角为ωt，转子上一点距z1轴为r，与x1的极角为ψ，则该点相对运动方程：x1＝rcosψ,y1＝rsinψ,z1＝z1(1)经坐标变换可得转子上任一点绝对运动方程：考虑直井中，振荡器较长外加更长的悬挂绳索，与套管之间的间隙相对小得多，即则s≈1，并略去高阶小量，转子上任一点绝对坐标、速度分别为设转子密度为ρ，内外半径分别为r1与r2，长为lz，上端距悬挂点为lz1，两侧面所夹圆心角为α，任一瞬时其对称面与x1的夹角为ωt时，两侧面与x1的夹角分别为为ωt-0.5α与ωt+0.5α。围绕转子任一点取微块质量为dm＝ρrdψddz1dr，则由扇形形心位置公式由式(3)知转子的重心坐标振荡器的转子的重力元功与动能分别为转子由电机轴、电机转子、偏心块、偏心块轴等四部分组成，设振荡器顶端距O点为lzo，钢的密度ρ＝7.85*103kg/m3，偏心块α＝0.5π，其余部分α＝2π，四部分其它参数分别为lz＝0.495m，lz1＝lzo+0.787，r1＝0，r2＝0.0085m；lz＝0.35m，lz1＝lzo+0.846m，r1＝0.0085m，r2＝0.0225m；lz＝0.47m，lz1＝lzo+1.455m，r1＝0.0135m，r2＝0.0358m；lz＝0.725m，lz1＝lzo+1.29m，r1＝0，r2＝0.0135m，由式(4)可确定各部分重力元功动能，并分别叠加获得整个转子部分重力元功与动能分别为进一步，振荡器定子的重力功与动能，由于定子是由圆筒或圆柱组成，将所有定子组成部分模拟为圆筒，任一瞬时其上一点距z1为r，与x1的极角为ψ，则该点在动系Ox1y1z1下：x1＝rcosψ,y1＝rsinψ,z1＝z1在静系Oxyz下其坐标与速度分别为其重心在z1轴上，即rc＝0,z1c＝lz1+0.5lz，由此知振荡器定子的重力元功与动能由式(7)计算各部分重力元功及动能，并分别叠加获得整个定子部分重力元功与动能分别为进一步，振荡器的重力元功与动能：略去相对小量或高阶小量，有进一步，套管内泥浆中振荡器上力系功：考虑泥浆对振荡器的影响，振荡器在泥浆中运动，将受到泥浆的阻尼作用，将水泥浆阻尼对振荡器横向运动所起的作用视为粘性阻尼，设单位表面积泥浆粘性阻尼系数为c，取振荡器外直径为D，长为l，顶端距悬挂点为l1，表面某点速度为v，围绕该点微小表面积dS上受到的泥浆阻尼力dFf＝-cvdS泥浆中振荡器的泥浆运动阻尼力功考虑泥浆对振荡器的影响主要反映在振荡器在泥浆中运动时迎泥浆半圆柱表面上的阻尼力，由式(6)和dS＝0.5Ddψdz1，并考虑则s≈1，且略去高阶小量，其元功为由式(11)分别计算其阻尼力元功并叠加它们获得振荡器泥浆阻尼力元功套管内泥浆中振荡器的力系功由式(9)、(12)，泥浆中振荡器系统的各力总元功进一步，振荡器的运动微分方程及运动方程；将(10)、(13)代入拉氏方程对于广义坐标θ分别有这是一个二自由度非线性微分方程组，考虑研究直井中振荡器对套管的作用在各过井轴的铅锤面内的几率是相同的，即主要考察振荡器对套管某铅锤面内的作用带来的影响即可，故将参量θ固定，设为θ＝0，即只考察θ＝0铅锤平面内的作用，则上述方程第一式有引入振荡器系统的阻尼比、固有频率、荷载系数方程(14)改写为方程(15)为是一个常系数二阶线性非齐次微分方程，表示振荡器主要为简谐激扰下的强迫振动，对应的齐次微分方程，设并考虑一般泥浆阻尼较小，即ξ<1，对应齐次微分方程特征根为共轭复根，即通解为复指数函数，根据复指数函数与三角函数的关系，通解解可写为其中，Ф1与δ为由初始条件确定的积分常数。为振荡器的衰减固有频率，方程(15)的特解，取为代入方程(15)，解得待定常数Ф与α于是方程(15)通解为式中，第一项初始条件引起的自由振动为瞬态振动，第二项振荡器偏心转子引起的强迫振动部分，为等幅的简谐振动即稳态振动。为了研究的方便，引入力幅静变形、频率比和振幅放大系数分别为则系统的响应相位差、振幅放大系数和振幅为Φ＝βΦ0,初始条件取为可解得Ф1与δ，从而方程(15)的解为由于振荡器开动前，即从而方程(15)的解为附图说明图1是振荡器转子模拟为扇形截面坐标示意图；图2是振荡器定子模拟为圆筒截面坐标示意图。具体实施方式下面结合具体实施方式对本专利技术进行详细说明。1、转子的重力做功与动能设振荡器转子角速度为ω，将所有转子组成部分模拟为图1所示扇形截面偏心块。任一瞬时其对称面与x1的夹角为ωt，转子上一点距z1为r，与x1的极角为ψ，则该点相对运动方程即动系Ox1y1z1下x1＝rcosψ,y1＝rsinψ,z1＝z1(1)经坐标变换可得转子上任一点绝对运动方程即任一瞬时静系Oxyz下坐标位置：考虑直井中，振荡器较长外加更长的悬挂绳索，与套管之间的间隙相对小得多，即则s≈1，并略去高阶小量，转子上任一点绝对坐标、速度分别为设转子密度为ρ，内外半径分别为r1与r2，长为lz，上端距悬挂点为lz1，两侧面所夹圆心角为α，任一瞬时其对称面与x1的夹角为ωt时，两侧面与x1的夹角分别为为ωt-0.5α与ωt+0.5α。围绕转子任一点取微块质量为dm＝ρrdψddz1dr，则由扇形形心位置公式由式(3)知转子的重心坐标振荡器的重力元功与动能分别为转子由电机轴、电机转子、偏心块、偏心块轴等四部分组成，设振荡器顶端距O点为lzo，钢的密度ρ＝7.85*103kg/m3，偏心块α＝0.5π，其余部分α＝2π，四部分其它参数分别为lz＝0.495m，lz1＝lzo+0.787，r1＝0，r2＝0.0085m；lz＝0.35m，lz1＝lzo+0.846m，r1＝0.00本文档来自技高网...

【技术保护点】
1.一种振荡器的受力与运动分析方法，其特征在于：首先分析转子的重力做功与动能、振荡器定子的重力功与动能、振荡器的重力元功与动能、套管内泥浆中振荡器上力系功，最后得到振荡器的运动微分方程及运动方程。

【技术特征摘要】
1.一种振荡器的受力与运动分析方法，其特征在于：首先分析转子的重力做功与动能、振荡器定子的重力功与动能、振荡器的重力元功与动能、套管内泥浆中振荡器上力系功，最后得到振荡器的运动微分方程及运动方程。2.按照权利要求1所述一种振荡器的受力与运动分析方法，其特征在于：所述转子的重力做功与动能如下，设振荡器转子角速度为ω，将所有转子组成部分模拟为扇形截面偏心块，任一瞬时其对称面与x1的夹角为ωt，转子上一点距z1为r，与x1的极角为ψ，则该点在动坐标系Ox1y1z1下的相对运动方程：x1＝rcosψ,y1＝rsinψ,z1＝z1(1)经坐标变换可得转子上任一点在静坐标系Oxyz下的绝对运动方程式中，[A]为静坐标系Oxyz和动坐标系Ox1y1z1间的坐标变换矩阵，s为y轴单位正向矢量和z1轴单位正向矢量矢乘的模，为振荡器轴线与z轴的夹角，θ为z1轴在xy面内投影与x轴夹角。考虑直井中，振荡器较长外加更长的悬挂绳索，与套管之间的间隙相对小得多，即则并略去高阶小量，转子上任一点绝对坐标、速度分别为设转子密度为ρ，内外半径分别为r1与r2，长为lz，上端距悬挂点为lz1，两侧面所夹圆心角为α，任一瞬时其对称面与x1的夹角为ωt时，两侧面与x1的夹角分别为为ωt-0.5α与ωt+0.5α。围绕转子任一点取微块质量为dm＝ρrdψddz1dr，则由扇形形心位置公式由式(3)知转子的重心坐标振荡器的转子的重力元功与动能分别为转子由电机轴、电机转子、偏心块、偏心块轴等四部分组成，设振荡器顶端距O点为lzo，钢的密度ρ＝7.85*103kg/m3，偏心块α＝0.5π，其余部分α＝2π，四部分其它参数分别为lz＝0.495m，lz1＝lzo+0.787，r1＝0，r2＝0.0085m；lz＝0.35m，lz1＝lzo+0.846m，r1＝0.0085m，r2＝0.0225m；lz＝0.47m，lz1＝lzo+1.455m，r1＝0.0135m，r2＝0.0358m；lz＝0.725m，lz1＝lzo+1.29m，r1＝0，r2＝0.0135m，由式(4)可确定各部分重力元功动能，并分别叠加获得整个转子部分重力元功与动能分别为。3.按照权利要求1所述一种振荡器的受力与运动分析方法，其特征在于：所述振荡器定子的重力功与动能，由于定子是由圆筒或圆柱组成，将所有定子组成部分模拟为圆筒，任一瞬时其上一点距z1为r，与x1的极角为ψ，则该点在动系Ox1y1z1下：x1＝rcosψ,y1＝rsinψ,z1＝z1在静系Oxyz下其坐标与速度分别为其重心在z1轴上，即rc＝0,z1c＝lz1+...

【专利技术属性】
技术研发人员：于桂杰，刘延强，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：山东,37