随钻扩眼工具的动力学仿真研究.pdfVIP

马清明等：随钻扩眼工具的动力学仿真研究 5 为对应于非完整约束的拉式乘子列阵；咖：为完整约 (图2中实线，约为30kN)，直至第2个刀翼开始与 束矩阵；彰为非完整约束矩阵。 井壁接触，第1个刀翼的受力逐渐减小直至等于零， 与井壁脱离，工具开始绕第2个刀翼与井壁接触形 采用拉格朗日(Lagmnge)方程建立系统动力学 方程的优点是程式化程度较高，适用于通过计算机 成的新的瞬时转动中心转动，第2个刀翼开始与井 高效求解。6一? 壁接触后，其受力开始逐渐增大(图2中虚点线)， 1．2．2动力学方程的求解对式(1)～(3)进行改达到最大值(大约50kN)后迅速稳定(约为30 写，首先将A和肚合成同一列阵，仍用A表示。将 kN)，上述过程不断重复，即形成扩大的井眼。结合 实际分析扩眼刀翼的动态响应，不难发现扩眼刀翼 非完整约束方程p(g，i，￡)=o对时间求导，并令u =西(即M为广义速度矩阵，相应地五为广义加速度 上的切削元件和钻柱在这种周期性冲击载荷的作用 矩阵)，然后代人式(1)，并写成新的函数F。构造下，寿命和可靠性将被逐渐削弱。 函数G(M，i)=配一百并令其等于o，即得动力学求解 方程如下 F(q，M，五，A，￡)=0 (4) G(M，香)=M一香=o (5) 中(q，￡)=0 (6) 式中，F为系统动力学微分方程；西为描述约束的代 数方程列阵。 采用Gear预估一校正算法求解式(4)～(6)所图2扩眼刀翼受力的变化规律(1个周期) 示的微分方程㈠引。 1．2．3边界条件基于法向接触约束的惩罚准则， 计算扩眼刀翼与岩石之问的接触力一一1。 0 g=0 F，=培 + 一 2蔓“ 堑出 0gd。。 无凶次半径 培 + 嵯亟出卫“建如 g≥d。、 (7) 式中，七是刚度系数，为罚参数；g代表2个几何体之 间的穿透值；dg／df代表在接触点的穿透率；e表示 力的指数，为正实数；d。。、代表当阻尼系数为最大值 c。Ⅲ、时，边界的最大穿透值，为正实数；c；。。。为最大阻 尼系数。 2仿真结果及分析 图4扩眼刀翼最外缘的运动轨迹包络线 2．1 扩眼工具在井眼内的运动 图3是单个刀翼最外缘的单圈运动轨迹，综合 设定扩眼工具有3个刀翼，扩眼刀翼的半径为 3个刀翼最外缘的连续运动轨迹，即得到运动轨迹 118．5 mm，工具的偏心率为5％，钻具的转速为100包络线如图4所示，因为实际刀翼外缘是安装有切 r／min，仿真结果见图2至图4(图中无因次半径为削元件的，因此图4也是井眼的形状。图4显示的 井眼半径与扩眼刀翼半径118．5mm的比值)。分 井眼内部均匀分布有多条“棱柱”，结合上述扩眼工 析图2和图3可发现，由于工具偏离井眼轴心，只有 具在井眼内的运动过程分析可知，“棱柱”的形成是 1个刀翼与井壁接触，在上部钻柱扭矩的作用下，导 由于扩眼工具偏心，扩眼刀翼间断性地切削井壁造 致扩眼工具以该刀翼与井壁的接触点为瞬时转动中 成的，在文中条件下，所形成井眼外轮廓的半径为 124．425 心转动，在转动的过程中，该刀翼的受力保持不变 mm(图4中实线)，图4中虚线为包络线内 马清明等：随钻扩眼工具的动力学仿真研究 7 3 结论