[论文解读] Low-Mach-number--slenderness limit for elastic Cosserat rods

Franziska Baus, Axel Klar+2 more|arXiv (Cornell University)|Jul 13, 2015

Micro and Nano Robotics参考文献 5被引用 4

一句话总结

本文建立了动态弹性Cosserat杆模型在低马赫数-细长比极限下的理论，表明当细长比和马赫数同时趋于零且其比值保持恒定（即 Ma/ε = μ = 常数）时，系统渐近收敛至无转动惯量的Kirchhoff梁模型。分析揭示了小参数至四阶的统一数学结构，从而支持一种渐近保持的数值格式，该格式在过渡区域中保持能量守恒与数值鲁棒性。

ABSTRACT

This paper deals with the relation of the dynamic elastic Cosserat rod model and the Kirchhoff beam equations. We show that the Kirchhoff beam without angular inertia is the asymptotic limit of the Cosserat rod, as the slenderness parameter (ratio between rod diameter and length) and the Mach number (ratio between rod velocity and typical speed of sound) approach zero, i.e. low-Mach-number--slenderness limit. The asymptotic framework is exact up to fourth order in the small parameter and reveals a mathematical structure that allows a uniform handling of the transition regime between the models. To investigate this regime numerically, we apply a scheme that is based on a spatial Gauss-Legendre collocation and an $α$-method in time.

研究动机与目标

• 严格建立在低马赫数与细长比极限下，动态弹性Cosserat杆模型与Kirchhoff梁模型之间的渐近关系。
• 识别极限系统的数学结构及其一阶校正项，其有效性可达 O(ε⁴) 阶，从而实现在过渡区域中的一致性数值处理。
• 基于Gauss-Legendre配点法与时间方向的α-方法，开发一种渐近保持的数值格式，确保在小 ε 与 Ma 条件下具有鲁棒性与能量守恒性。
• 证明当Euler-Bernoulli关系控制接触力偶且转动惯量项消失时，极限模型即为Kirchhoff-Love方程。
• 通过数值验证，展示该格式在半离散化产生的DAE索引（0 至 3）范围内均具有收敛性与稳定性。

提出的方法

• 在小参数 ε（细长比）下进行渐近展开，其中 Ma/ε = μ = 常数，以推导极限系统及其一阶校正项。
• 采用幂级数展开至 O(ε⁴) 阶，以确保在此阶次内渐近框架的精确性。
• 应用空间Gauss-Legendre配点法对空间导数进行离散化，形成微分代数方程组（DAE）。
• 在时间方向采用α-方法，通过可调节参数（λ = 0.5 时保证能量守恒，λ = 1 时引入耗散），以处理极限区域中的刚性问题。
• 对Cosserat杆（ε > 0）、Kirchhoff梁（ε = 0）及校正项（O(ε²)）系统统一处理DAE，其索引范围为 0 至 3。
• 通过二维Euler-Bernoulli悬臂梁的数值收敛性研究验证该格式，分析空间与时间方向的收敛阶数。

实验结果

研究问题

• RQ1 当细长比 ε 与马赫数 Ma 同时趋于零且 Ma/ε = μ = 常数时，动态弹性Cosserat杆模型的渐近极限为何？
• RQ2 极限系统及其一阶校正项的数学结构如何支持从Cosserat杆到Kirchhoff梁过渡区域中的一致性数值处理？
• RQ3 所提出的渐近保持格式在半离散化产生的索引为 0 至 3 的DAE系统中，其收敛性与稳定性特性如何？
• RQ4 在超弹性本构关系与适当边界条件下，该数值格式如何保持能量守恒？
• RQ5 该数值格式在 ε 极小时（系统趋于刚性）是否仍能保持精度与鲁棒性？

主要发现

• 无转动惯量的Kirchhoff梁模型是动态Cosserat杆模型在低马赫数-细长比区域下的精确渐近极限，其有效性可达 O(ε⁴) 阶。
• 极限系统及其一阶校正项构成一个在数学上统一且精确至 ε 四阶的框架，支持一致的数值处理。
• 所提出的渐近保持格式在 Δt = Δs → 0 时达到 p = 1 的收敛阶，当 λ = 0.5 时在 ε 相关情形下达到 p = 2，与理论预期一致。
• 数值收敛性研究确认空间离散化收敛阶为 p_s = 2，时间离散化收敛阶为 p_t ≈ 1（或更低），即使在 DAE 索引为 3 的情况下也成立，尽管缺乏理论收敛性保证。
• 在使用超弹性本构关系、外力势能及适当边界条件时，该格式可保持能量守恒，且可通过调节 λ 参数切换至耗散模式。
• 极限系统的雅可比结构在不同模型变体（M）、（T）、（S）中均表现出鲁棒性，所有情况下均观察到一致的数值行为。