Buckling of slender string in cylindrical tube under axial load: Experiments and theoretical analysis

Buckling of drill and tubing strings in drill or casing holes will affect the life of the string and cause difficulties in drilling and in oil production. A number of theoretical models have been developed to study the buckling load and the post-buckling configuration of the strings under some idealized conditions. However, verification of these theories in practice or in laboratory has not been extensively reported.In this paper work on laboratory buckling tests of strings is presented and the results are compared with theoretical formulas. The results indicate that the process of string buckling can be divided into two stages, namely plane buckling and spatial buckling. In the spatial buckling the string-buckling configuration can roughly be represented by a helix. The relation between the buckling force and the helical pitch, obtained from the laboratory tests, is similar to that predicted by helical buckling theories. The friction between the string and the wall of the confining tube causes a nonuniform buckling shape and prohibits further buckling at some critical loads. The buckling theories do not cover these effects and may thus induce significant errors in predicting the string-post-buckling configuration. An approximate formula for estimating the friction force is proposed, which provides more accurate results. Further studies of the friction effect are suggested.     

液压作用下套管柱临界屈曲载荷计算的能量法

推导了内外液压作用下套管柱微弯时的能量平衡公式,并计算了其屈曲临界载荷计算公式。得到的结果低于目前已发表的结果。内外液压对套管稳定性的影响相当于在套管底端作用一个附加的轴向力和一个沿轴线均布的线载荷。对含有内压的两端封堵的细长管柱,内压对临界屈曲载荷没有影响。结果可供油井设计及作业参考。     

The computer simulation of drill column dragging in inclined bore-holes with geometrical imperfections

The problem about identification of resistance forces acting on a drill column moving in an inclined bore-hole is stated. It is supposed that the well trajectories can have geometrical imperfections in the shape of cylindrical spiral or plane cosinusoidal curves. The system of ordinary differential equations is derived on the basis of the theory of curvilinear flexible elastic rods. It permits one to describe static effects of the drill column bending accompanying the processes of its raising, lowering and rotating inside the bore-hole. Through the use of this system the direct and inverse problems of the drill column deforming are formulated for calculation of internal and external resistance forces acting on the drill column tube. The methods for numerical solution of the constructed equations are elaborated. With their use the phenomena of the drill columns motion and their frictional seizure inside the bore-holes are simulated for different geometrical imperfections and relations between the velocities and directions of their rotation and axial motion.     

Twin-characteristic-parameter solution of axisymmetric dynamic plastic buckling for cylindrical shells under axial compression waves

In the present investigation on the dynamic plastic buckling of cylindrical shells under axial compression waves, the critical axial stress and the exponential parameter of inertia terms in stability equations are treated as a couple of characteristic parameters. The criterion of transformation and conservation of energy in the process of buckling initiation is used to derive the supplementary restraint equation of buckling deformation at the fronts of axial elastic and plastic compression waves. The supplementary restraint equation, stability equations, boundary conditions and continuity conditions constitute the necessary and sufficient conditions of determining the two characteristic parameters. Two characteristic equations are derived for the two characteristic parameters. The critical axial stress or the critical buckling time, the exponential parameter of inertia terms and the initial modes of buckling deformation are calculated quantitatively from the solution of the characteristic equations.     

New method for geometric nonlinear analysis of large displacement drill strings

Introduction Drillstringsina3Dlargedisplacementwellboreareusuallyregardedas3Dcurved beamswithlargedisplacementandlargerotationbutwithsmallstrain,thus,therelationship betweenstressandstrainislinear.Thedeformationofdrillstringscanbecomposedoftwo parts:oneiscalledtheinitialdeformationcausedwhentheaxisofthedrillstringisdeformed fromitsstraightconfigurationtothepositioncoincidentwiththewellaxis,involvinglarge rotationsandwiththeimplicationoffiniterotationalproblem.Theotheristhedeformation fromthew…     

Dynamic buckling of cylindrical shells subject to an axial impact in a symplectic system

This paper discusses the dynamic pre-buckling of finite cylindrical shells in the propagation and reflection of axial stress waves. By introducing the Hamiltonian system into dynamic buckling of structures, the problem can be described mathematically in a symplectic space. The solutions of Hamiltonian dual equations shown in canonical variables are obtained. The problem is reduced to the determination of eigenvalues and eigensolutions, with the former indicating critical buckling loads and the latter buckling modes. Numerical example presented shows phenomena of axisymmetric and non-axisymmetric dynamic buckling subject to impacts of axial load.     

黏弹性层合圆柱壳的临界轴压分析

基于经典屈曲理论,研究了轴向受压黏弹性复合材料层合圆柱壳的临界屈曲载荷. 利用Boltzmann线性积分型本构关系描述铺设单层的各向异性黏弹性行为. 结合解析与数值 方法,由Donnell型屈曲控制方程以及边界条件的Laplace变换导出相空间的特征方程,根 据Laplace逆变换的极值定理,获得层合圆柱壳的瞬时弹性临界载荷与持久临界载荷. 针对 多组铺设方式,通过数值算例重点分析了临界载荷随铺设角的变化特征,两种临界载荷的峰 值点差异程度与铺设方式、几何参数以及材料类型的关系,得到了一些对黏弹性层合圆柱壳 的优化设计有参考价值的结论.     

SIMULATION OF TRANSVERSE VIBRATION OF SLENDER ROD STRING IN THE VERTICAL CYLINDER UNDER BUCKLING DISPLACEMENT EXCITATION1)

对于铅直圆筒内受交变拉压轴向载荷作用的细长杆柱,当杆柱底端所受到的轴向压力大于杆柱屈曲的临界载荷时,细长杆柱在圆筒内将产生螺旋屈曲,杆柱的屈曲变形将激励杆柱在圆筒内产生横向振动。以细长杆在圆筒内的瞬时屈曲构型作为杆柱横向振动的位移激励,建立了屈曲位移激励下的细长杆在圆筒内横向振动与杆管碰撞规律的仿真模型。采用有限差分法对井深进行离散,Newmark 法对时间进行离散,以恢复系数法建立了细长杆和圆筒的碰撞条件,对细长杆在圆筒内的横向振动数学模型进行了数值仿真。仿真结果表明,细长杆和圆筒内壁的碰撞现象主要发生在细长杆底端受压发生屈曲后,且几乎沿整个圆筒都有发生,从圆筒顶部至底部的碰撞力峰值逐渐增大;而在杆柱底端受拉时碰撞现象很少,碰撞力也较小。     

Spatially complex localisation in twisted elastic rods constrained to a cylinder

We consider the problem of a long thin weightless rod constrained to lie on a cylinder while being held by end tension and twisting moment. Applications of this problem are found, for instance, in the buckling of drill strings inside a cylindrical hole. In a previous paper the general geometrically exact formulation was given and the case of a rod of isotropic cross-section analysed in detail. It was shown that in that case the static equilibrium equations are completely integrable and can be reduced to those of a one-degree-of-freedom oscillator whose non-trivial fixed points correspond to helical solutions of the rod. A critical load was found where the rod coils up into a helix.Here the anisotropic case is studied. It is shown that the equations are no longer integrable and give rise to spatial chaos with infinitely many multi-loop localised solutions. Helices become slightly modulated. We study the bifurcations of the simplest single-loop solution and a representative multi-loop as the aspect ratio of the rod's cross-section is varied. It is shown how the anisotropy unfolds the `coiling bifurcation'. The resulting post-buckling behaviour is of the softening–hardening–softening type typically seen in the cellular buckling of long structures, and can be interpreted in terms of a so-called Maxwell effective failure load.     

Effects of local thickness defects on the buckling of micro-beam

A buckling model of Timoshenko micro-beam with local thickness defects is established based on a modified gradient elasticity. By introducing the local thickness defects function of the micro-beam, the variable coe-cient differential equations of the buckling problem are obtained with the variational principle. Combining the eigensolution series of the complete micro-beam with the Galerkin method, we obtain the critical load and buckling modes of the micro-beam with defects. The results show that the depth and location of the defect are the main factors affecting the critical load, and the combined effect of boundary conditions and defects can significantly change the buckling mode of the micro-beam. The effect of defect location on buckling is related to the axial gradient of the rotation angle, and defects should be avoided at the maximum axial gradient of the rotation angle. The model and method are also applicable to the static deformation and vibration of the micro-beam.     

黏弹性板的大挠度蠕变屈曲

研究了具有初始小挠度受轴向压载黏弹性板的蠕变屈曲问题,在建立控制方程时,利用了von Karman非线性应变-位移关系,并考虑了初始挠度,用标准线性固体模型描述材料的黏弹性特性,在求解非线性积分方程时,利用梯形公式计算记忆积分式,将非线性积分方程化为非线性代数方程进行数值求解,得到了结构的蠕变变形过程,又将问题退化到小挠度情况进行研究,得到了挠度随时间扩展的解析解,分析了瞬时失稳临界载荷、持久临界载荷的物理意义,讨论了考虑几何非线性对黏弹性板蠕变屈曲的影响。     

弹性圆柱壳在轴向冲击载荷和温度耦合作用下的屈曲

通过引入哈密顿体系,将临界载荷和临界温度及它们所对应的屈曲模态归结为辛体系下的广义本征值和本征解问题。根据辛本征解的正交性和完备性,给出了全部的且独立存在的屈曲模态。数值结果表明,在轴向冲击载荷和温度耦合作用下,弹性圆柱壳的屈曲呈现出复杂的模式,温度直接影响冲击临界载荷的大小。随着温度的增加,冲击临界荷载降低,最后,文中给出各种条件下的屈曲模态。     

Frequency analysis of nonlinear free vibration of a cross string

This work concerns nonlinear free vibration of a cross string under large amplitude. The equations governing the nonlinear vibration of the cross string are derived at first from the Hamilton principle, and they take the form of Duffing equation. Then the perturbation method is used to solve the nonlinear coupled natural frequency of the cross string. The nonlinear natural frequency not only has the characteristic of nonlinearity, but also reflects the coupled characteristic, i.e., the natural frequency of the cross string varying with that of its constituent strings. The results show that the overall effect on the cross string is somehow averaged due to the nonlinearity of each constituent string, i.e., the natural frequencies of the cross string contain both the linear natural frequencies of the constituent strings and the nonlinear parts that depend upon the vibration amplitude, the diameter of one constituent string, the length ratio of the two strings, etc., but the contribution of each constituent string to the natural frequency is in different proportions.     

Buckling of long thick-walled circular cylindrical shells of isotropic incompressible hyperelastic materials under uniform external pressure

For isotropic incompressible hyperelastic materials, the problem of determining the critical external pressure at which a long thick-walled circular cylindrical shell will buckle involves solving a fourth-order system of highly non-homogeneous, ordinary differential equations. Closed-form solutions of this system are derived here for plane-strain conditions and for the particular case of the Varga material. These solutions are used to derive the buckling criterion and numerical values are obtained for the resulting critical pressures. They are found to be in good agreement with existing theoretical and experimental results for the neo-Hookean material.     

轴向均布载荷下压杆稳定问题的DQ解

叙述了微分求积法(differential quadrature method)的一般方法，研究用微分求积法求解在均布轴向载荷下细长杆的稳定问题．通过Newton-Raphson法求解非线性方程组，以及对问题进行线性假设后求解广义特征值方程，得到了精度很高的后屈曲挠度数值和临界载荷数值．与解析解和其他近似解相比，微分求积法具有较高的精度和简便性．     

Controlling torsional vibrations of drill strings via decomposition of traveling waves

Torsional vibrations in drill strings, especially stick-slip vibrations, are detrimental to the drilling process as they slow down the rate of penetration and may lead to failure of the drilling equipment. We present a method for controlling these vibrations by exactly decomposing the drill string dynamics into two traveling waves traveling in the direction of the top drive and in the direction of the drill bit. The decomposition is derived from the wave equation governing the string vibrations and is achieved with only two sensors that can be placed directly at the top drive and at a short distance below the top drive (e.g., 5 m). Therefore, downhole measurements along the string and at the bit are not necessary, which is a major advantage compared to other control concepts for drill string dynamics. The velocity of the top drive is then controlled in order to absorb the wave traveling in the direction of the top drive, thus achieving a reflection coefficient of zero for the frequency range of the undesired torsional vibrations. The proposed algorithm is implemented for both a numerical example and an experimental setup; results show that the control concept works very effectively.     

Post-buckling analysis of slender elastic vertical rods subjected to terminal forces and self-weight

This article presents the behavior of slender elastic rods subjected to axial terminal forces and self-weight. The mathematical formulation is presented, a solution is sought for a double-hinged boundary condition and the analysis is carried out for different values of non-dimensional weight. The formulation derives from geometrical compatibility, equilibrium of forces and moments and constitutive relations yielding a set of six first order non-linear ordinary differential equations with boundary conditions specified at both ends, which characterizes a complex two-point boundary value problem. Furthermore, a perturbation method is used to find the critical buckling loads and initial post-buckling solutions. A numerical integration scheme based on a three parameter shooting method is employed in the post-buckling solutions.     

直杆塑性动力屈曲的双特征参数解

在关于轴向弹塑性压缩波下直杆塑性动力屈曲的研究中,将临界应力和屈曲惯性项指数参数作为一对特征参数求解．由动力屈曲发生瞬间的能量转换和守恒准则,导出波阵面上的屈曲变形补充约束条件．失稳控制方程、边界条件、塑性波阵面上的连续条件和补充约束条件构成定量求解两个特征参数和动力屈曲模态的完备条件。     

Method of equilibrium differential equation for analysis of strength of large deflection drill string

ＩｎｔｒｏｄｕｃｔｉｏｎＤｒｉｌｌｓｔｒｉｎｇｉｓｍｏｒｅｅａｓｉｌｙｄａｍａｇｅｄｓｏｔｈａｔａｃｃｉｄｅｎｔｉｎｐｒｏｄｕｃｔｉｏｎｏｃｃｕｒｓｍｏｒｅｏｆｔｅｎｉｎｌａｔｅｒａｌｈｏｒｉｚｏｎｔａｌｗｅｌｌｔｈａｎｉｎｇｅｎｅｒａｌｗｅｌｌｂｅｃａｕｓｅｔｈｅｄｉａｍｅｔｅｒｏｆｄｒｉｌｌｓｔｒｉｎｇｉｓｓｍａｌｌａｎｄｔｈｅｌａｒｇｅｃｏｍｐｌｅｍｅｎｔａｒｙｄｅｆｌｅｃｔｉｏｎｌｏａｄｓｏｎｄｒｉｌｌｓｔｒｉｎｇａｒｉｓｅｓｂｙｔｈｅｗｅｌｌｂｏｒｅｏｆｖｅｒｙｌａｒｇｅｃｕｒｖａｔｕ…     

The effect of variations in the creep exponent on the buckling of circular cylindrical shells

The senior author solved the problem of axially symmetrical creep buckling of thin circular cylindrical shells subjected to uniform axial compression. In that analysis the constitutive equation was a power law, and the exponent was taken to be equal to three. The purpose of this work was to extend the solution to a range of values of the creep exponent, n. To cope with the increasing algebraic complexity, a digital computer was employed in two ways: to generate the set of equations symbolically, and then to solve these equations. The computer programs were used to generate numerical solutions for the cases in which n was equal to 3, 5, 7 and 9. Two simple extrapolation techniques were then employed to obtain approximate solutions to the critical time problem for values of n up to 29.