基于轴向力的钻柱瞬态动力学建模.docxVIP

基于轴向力的钻柱瞬态动力学建模 单元动力学模型 在钻井过程中，钻头（尤其是牙轮钻头）的牙齿交替地接触不均匀的井底，导致钻头柱垂直滑动和垂直交变强度。这种交变应力的存在使钻柱的纵向振动、横向振动、扭转振动之间相互耦合,形成一种非线性的、复杂的且不断变化的振动。基于钻柱的这种运动形式,钻柱的破坏主要是在动载荷作用下发生的。关于钻柱瞬态动力学的研究目前已有大量的报道,但大部分处于理论阶段,计算工作量大,而且针对钻柱的动态失效分析还比较少。为此,笔者利用大型有限元软件的非线性瞬态动力学分析功能,分析钻柱承受任意随时间变化的钻压的动力响应,提取钻柱上的危险点,确定钻柱的失效位置,这对钻柱的优化设计具有一定的意义。 根据钻柱的实际工作情况,建立钻柱的动力学模型,如图1所示。由于在ANSYS有限元软件中,过于细长杆的非线性问题分析容易造成求解过程不收敛,得不到计算结果,况且,钻柱上部主要承受拉应力的作用,出现异常失效的情况很少,因此,笔者只研究下部钻柱组合的动力学特性,把握住钻柱的主要运动特征,通过对钻柱边界条件的合理设置,使其能真实反应钻柱受力振动的情况。由于钻柱类似于一根细长杆,结构简单,在这里采用直接建模法,即通过节点直接生成有限单元模型。现做如下简化和假设: (1)钻柱为均质弹性直杆; (2)考虑下部钻柱在钻压作用下发生屈曲而与井壁的接触; (3)考虑钻柱旋转产生的向心力对变形的影响; (4)没有考虑内外钻井液流动过程中对钻柱动力学产生的扰动影响。 考虑钻柱与井壁的随机接触碰撞,在钻柱每个单元节点上设置一个间隙元,对所有间隙元进行坐标转化和一系列对号入座拼装过程,形成总体刚度矩阵,得到钻柱在总体坐标系下的动力学平衡方程为 [M]{u¨}+[C]{u˙}+([K]+[KG]){u}={F(t)}(1)[Μ]{u¨}+[C]{u˙}+([Κ]+[ΚG]){u}={F(t)}(1) 式中 [M]、 [C]——整体钻柱的质量矩阵和阻尼矩阵; [K]、[KG]——整体钻柱的刚度矩阵和所有动力间隙元刚度矩阵; {u¨u¨}、{u˙u˙}、{u}、{F(t)}——钻柱节点的加速度向量、速度向量、位移向量和等效节点力向量。 根据钻柱的动力学特性,选用Newmark方法进行求解,该方法中近似假设 {u˙t+Δt}={u˙t}+[(1?β){u¨t}+β{u¨t+Δt})]Δt(2){u˙t+Δt}={u˙t}+[(1-β){u¨t}+β{u¨t+Δt})]Δt(2) 经过一系列推导,由t时刻的已知参数求t+Δt的计算公式为 (1αΔt2[M]+βαΔt[C]+[K]+[KG]){ut+Δt}={Ft+Δt}+[M][1αΔt2{ut}+1αΔt{u˙t}+(12α?1){u¨t}]+[C][βαΔt{ut}+(βα?1){u˙t}+(β2α?1)Δt{u¨t}](3){u¨t+Δt}=1αΔt2({ut+Δt}?{ut})?(1αΔt2[Μ]+βαΔt[C]+[Κ]+[ΚG]){ut+Δt}={Ft+Δt}+[Μ][1αΔt2{ut}+1αΔt{u˙t}+(12α-1){u¨t}]+[C][βαΔt{ut}+(βα-1){u˙t}+(β2α-1)Δt{u¨t}](3){u¨t+Δt}=1αΔt2({ut+Δt}-{ut})- 1αΔt{u˙t}?(12α?1){u¨t}(4){u˙t+Δt}=βαΔt({ut+Δt}?{ut})+1αΔt{u˙t}-(12α-1){u¨t}(4){u˙t+Δt}=βαΔt({ut+Δt}-{ut})+ (1?βα){u˙t}+(1?β2α)Δt{u¨t}(5){ut+Δt}={ut}+{u˙t}Δt+12×(1-βα){u˙t}+(1-β2α)Δt{u¨t}(5){ut+Δt}={ut}+{u˙t}Δt+12× (1?2α){u¨t}+2α(u¨t+Δt)(1-2α){u¨t}+2α(u¨t+Δt)Δt2(6) 式中α、β——Newmark积分参数(0≤2α≤1,0≤β≤1)。 1. 钻柱的高应力交变分析 根据上述动力学模型和理论方法,对新疆焉耆盆地山前构造推覆体城二井所采用的某钻柱组合进行非线性瞬态动力学分析,该钻柱的钻具结构为: ?311.2 mm(121/4英寸)钻头+转换接头+?203.2 mm(8英寸)减振器1根+?203.2 mm(8英寸)钻铤1根+扶正器+?203.2 mm(8英寸)钻铤1根+扶正器+转换接头+?203.2 mm(8英寸)钻铤3根+转换接头+?177.8 mm(7英寸)钻铤4根+转换接头+?158.8 mm(61/4英寸)钻铤6根+转换接头+?127 mm(5英寸)加重钻杆18根+?127 mm(5英寸)钻杆至井口。 相关参数为:井斜角1°;钻压137.2 kN;转速65 r/min;