轴向载荷下周期梁超材料带隙可调工程：基于Bloch-Floquet理论的钻柱横向波传播与振动控制框架

《Next Materials》：Tunable band-gap engineering in axially loaded periodic beam metamaterials: A drill string lateral dynamics case study

　　

轴向张力（axial tension）控制横向波速、波长和频率。周期性工具接头接触（periodic tool-joint contacts）产生阻带（stop-bands），使横向扰动局域化。轴向压缩（axial compression）使零频阻带（zero-frequency stop band）变窄，指示屈曲起始（buckling onset）。Bloch-Floquet理论提供带隙调制（band-gap modulation）的闭合形式框架。该框架为自动化钻井系统（automated drilling systems）中的自适应振动控制（adaptive vibration control）提供依据。

论文解读：

一、研究背景与问题。钻柱（drill string）是地表设备与钻头之间的主要机械连接，传递旋转动力、钻井液和轴向钻压（WOB, weight-on-bit），由钻杆、重钻杆（HWDP）和井底钻具组合（BHA, bottom hole assembly）组成。钻柱动力学包括轴向、扭转和横向三种耦合模态，其中横向动力学（lateral dynamics）最复杂，涉及垂直于中性轴的四个自由度，并对轴向载荷状态敏感。WOB增大降低下部BHA净轴向张力，降低弯曲刚度并放大弯曲位移。钻杆因离散工具接头与井壁接触形成周期性，改变波传播并产生类似光子晶体的频带隙。已有连续梁模型未严格处理操作轴向力与工具接头接触空间周期性耦合。因此，研究人员采用将轴向载荷Euler–Bernoulli梁与周期接触效应集成的分析框架，以刻画钻柱横向动力相互作用和力学行为。研究目标与新颖性在于提出将Euler–Bernoulli梁与Bloch–Floquet波理论结合的闭合形式框架，证明轴向张力是控制振动阻带位置和宽度的连续调谐参数，并为自动钻井系统频率靶向减振提供物理基础。

二、关键技术方法（不超过250字）。研究人员采用解析与数值结合方法。将钻柱简化为周期铰支梁，工具接头周期接触井壁。依次建立零轴向载荷静力模型、含惯性Euler–Bernoulli动态模型、含轴向载荷与周期接触全动态模型。基于Bloch–Floquet周期条件和传递矩阵推导闭合色散关系，定义传播波数??_p与衰逝波数??_e。数值仿真在MATLAB中完成，参数依据API 5DP规范：基线跨度?=24 m，??=12 m，总长??_tot=240 m，??=200 GPa，????=1.19×10⁶ N m²，钢密度7850 kg/m³，泥浆密度1000–1800 kg/m³，单位长度质量29.02 kg/m。轴向载荷用浮力修正公式计算。以独立两节点Euler–Bernoulli梁有限元法（FEM, finite element method）验证。无样本队列，数据来自API规范和模型参数。

三、研究结果。

3.1 模型参数：基于API 5DP Range 3双接头钻杆，采用基线参数。结论是为单跨和全钻柱仿真提供代表值。

3.2 轴向载荷剖面：解析计算浮力修正轴向张力，并在离散深度点数值验证。结论是张力从表面最大208.8 kN线性降至中和点2617.6 m为零，底部为?30.5 kN压缩；深度变化轴向载荷直接改变每跨有效弯曲刚度，形成深度依赖色散关系。

3.3 结果：单跨分析。3.3.1 零频阻带（无轴向载荷）：基线仿真确认cos(θ)=?2，|cos(θ)|>1，处于阻带；横向位移逐跨指数衰减，无波传播。结论是静态横向扰动局域于起源跨。3.3.2 频率相关通/阻带结构：非零频Euler–Bernoulli动态模型显示随ω增加从阻带进入通带；四个基本横向响应形状为两个衰逝和两个传播。结论是四模态结构源于四阶控制方程。3.3.3 全动态模型：含轴向载荷与周期接触。色散结果显示轴向载荷根本移动cos(θ)。拉伸时低频cos(θ)钳制在?2，形成深阻带，约ω=4.27 rad/s穿过通带；压缩时相演化更平滑，ω=3.71 rad/s位于通带内；无载时ω=1.42 rad/s穿过cos(θ)=0。结论是轴向拉伸扩大第一阻带并移向高频，增加高张力段衰减；压缩逐步关闭阻带，促进波传播，增加振动敏感性；周期接触修改每跨有效传递矩阵，放大传播特性。3.3.4 静态屈曲阈值：闭合阈值??_cr=?????_p,cr²=???²??/(4??²)与经典Euler屈曲载荷??_cr=??²????/?²一致。对基线跨|??_cr|≈20.4 kN；轴向载荷剖面中于约2873 m超过阈值，位于最后约127 m。结论是近钻头段单跨超过静态屈曲阈值，与BHA稳定器作用一致。

3.4 结果：全钻柱分析。3.4.1 波传播的空间变化：顶部最大张力时波传播快、??_p小、波长长、局部阻带最宽且移向高频；向BHA张力降低，相速度下降、??_p增大、波长缩短、阻带变窄，接近通带边界。结论是上部高张力段如高通滤波器，只允许足够高频横向扰动，排除低频内容。3.4.2 非零频响应（无轴向载荷）：展示无轴向载荷、零频下周期铰支钻柱横向位移模态。3.4.3 频率相关通/阻带结构：展示由有效质量导致的有限频率效应。3.4.4 轴向载荷下横向位移：仅分析轴向载荷作用的横向位移模态。3.4.5 轴向载荷与周期接触下横向位移：对240 m钻柱在拉伸、压缩、无载下计算归一化横向位移模态，工具接头处??=0。拉伸使振荡周期延长、接触间幅值受抑，与阻带扩大一致；压缩使振荡周期收缩、幅值升高，与阻带变窄和波传播增强一致；无载振荡最显著。色散结果中，拉伸ω?3.6 rad/s时cosθ=?2，约4.27 rad/s进入通带；压缩通带起始移向低频，ω=3.71 rad/s时cosθ=0.317；无载ω=1.42 rad/s时cosθ=0。结论是轴向载荷和周期边界共同控制空间波长和接触间振幅。

3.5 独立有限元验证：用两节点Euler–Bernoulli梁单元建模单元胞，含一致刚度、质量和几何刚度，施加Bloch周期性，缩聚动态刚度矩阵，由cos(θ)=?(S₁₁+S₂₂)/(2S₁₂)求色散。结论是FEM随网格加密单调收敛到闭合色散关系，n=16时误差低于10^?2%，在ω∈[0.05,10] rad/s范围吻合，确认闭合色散关系为独立数值模型的精确连续极限。

四、讨论部分总结。4.1 关键结果的物理解释：无轴向载荷时零频宽阻带cosθ=?2，横向波在工具接头间衰减、扰动局域；轴向拉伸使阻带加宽，顶部波传播快、近中和点慢；轴向压缩使阻带变窄或消失，ω=0处非零??_p指示屈曲失稳起始；高频时轴向载荷影响可忽略；周期接触产生通/阻带结构，可由轴向载荷调谐，与BHA设计相关。4.2 分岔机制：轴向载荷是带隙结构的连续分岔参数；从高张力宽阻带经逐步变窄到完全波传播，遵循色散关系平滑曲线，原则上可由地面可测大钩载荷数据预测，但需扩展理想假设。4.3 周期接触：周期井壁接触打破铰支梁平移对称性，产生轴向载荷可调通带结构，并显著放大传播模态分量，与BHA密集稳定器段横向振动增强数值观察一致。4.4 对振动控制的意义：阻带频率范围可作为天然振动滤波器；若主导横向激励频率落入阻带，横向波不能传播，能量局域。调整WOB或RPM移动激励频率相对带隙可借助此过滤机制；模型提供操作参数与带隙边界关系，可作为控制器设计输入之一，但需解决理想化问题。4.5 有限元验证：闭合色散关系与独立FEM比较，随网格从n=1到n=16单调收敛，n=16误差低于10^?2%，确认精确连续极限；拉伸n=4轻微非单调偏差归因于工具接头附近曲率梯度粗网格，而非模型差异。4.6 局限：假设理想无摩擦铰支，排除流固耦合，限于垂直井轴向载荷剖面和Euler–Bernoulli理论；这些理想化会定量移动带隙边界，但保持定性分岔行为。

五、研究结论翻译。本研究对可变轴向载荷下周期铰支钻柱横向动力学进行严格解析和数值研究。主要结论：轴向张力控制横向波传播速度：顶部张力较高产生更长波长、更快传播和更宽零频阻带，近中和点张力较低产生较慢传播和较窄阻带。在零轴向载荷、零频极限，模型一致确认完全阻带，无实波数。轴向压缩使零频阻带变窄或消除，提供即将屈曲失稳的直接波传播特征。阻带关闭的闭合压缩阈值与等效铰支-铰支柱的经典Euler屈曲载荷完全一致，严格连接Bloch-Floquet框架与经典屈曲理论。四个基本横向响应模态共存：两个衰逝和两个传播。周期工具接头-井壁接触显著放大传播分量，是连续梁模型缺失的关键效应。闭合色散关系由周期单元胞有限元梁模型独立验证，两种方法在所有轴向载荷情况和所考虑全频范围内接近收敛，确认色散关系是独立数值模型的精确连续极限，而非传递矩阵推导的假象。Bloch-Floquet色散关系提供闭合形式、计算高效框架，预测带隙结构随轴向载荷、频率、弯曲刚度和跨长的函数。轴向载荷调制的带隙构成分岔机制，其边界原则上可与地面可测钻井参数关联；在采用的理想无摩擦无阻尼公式内，这为频率靶向振动缓解策略建立物理基础，而非已证明的实时控制能力。未来扩展包括用真实摩擦接触力学替代理想周期铰支，纳入环形钻井液附加质量和黏性阻尼，考虑钻柱旋转和间歇工具接头接触，扩展至斜井和水平井几何及非线性轴向载荷剖面，采用Timoshenko梁公式用于高频区，并用井下MWD振动数据验证。这些扩展是将带隙调制框架从解析和FEM验证基础转化为实时振动控制策略所必需。