[论文解读] Computational mechanics of soft filaments

Mattia Gazzola, Levi H. Dudte+2 more|arXiv (Cornell University)|Jul 1, 2016

Advanced Materials and Mechanics被引用 5

一句话总结

本文提出了一套鲁棒且完全符合物理规律的计算框架，基于SE(3)上的完整Cosserat杆理论，用于模拟软丝状物在六自由度变形（弯曲、扭转、拉伸和剪切）下的动力学行为。该方法能够精确模拟与环境之间的复杂相互作用，包括摩擦力、流体动力学效应以及主动肌肉力，并在正向问题与反向问题中均表现出色，例如无肢运动与扭结结构的形成。

ABSTRACT

Soft slender structures are ubiquitous in natural and artificial systems and can be observed at scales that range from the nanometric to the kilometric, from polymers to space tethers. We present a practical numerical approach to simulate the dynamics of filaments that, at every cross-section, can undergo all six possible modes of deformation, allowing the filament to bend, twist, stretch and shear, while interacting with complex environments via muscular activity, surface contact, friction and hydrodynamics. We examine the accuracy of our method by means of several benchmark problems with known analytic solutions. We then demonstrate the capabilities and robustness of our approach to solve forward problems in physics and mechanics related to solenoid and plectoneme formation in twisted, stretched filaments, and inverse problems related to active biophysics of limbless locomotion on solid surfaces and in bulk liquids. All together, our approach provides a robust computational framework to characterize the mechanical response and design of soft active slender structures.

研究动机与目标

• 开发一种数值鲁棒且高效的计算方法，用于模拟可发生全部六种形变模式（弯曲、扭转、拉伸与剪切）的软丝状物，突破传统不可伸长杆模型的局限。
• 在单一统一框架中整合真实物理相互作用，包括各向异性摩擦、流体阻力、表面接触以及主动肌肉力。
• 通过解析基准测试验证该方法，并展示其在软活性力学中复杂正向与反向问题求解中的能力。
• 通过模拟由非对称摩擦与肌肉激活模式驱动的陆地与水下运动，实现对无肢生物运动策略的反向设计。

提出的方法

• 该方法基于完整欧氏群SE(3)的离散Cosserat杆公式，允许每个截面的中线位置与局部取向独立演化。
• 采用哈密顿结构与快速离散算子，以保持能量与动量守恒，确保数值稳定与精确。
• 模型采用广义本构律，分别使用弯曲、扭转、拉伸与剪切的刚度矩阵，材料参数包括弹性模量与剪切模量。
• 通过基于惩罚的方法实现接触力的稳健处理，包含各向同性与各向异性摩擦模型，涵盖静摩擦与动摩擦，并设置与速度相关的临界阈值。
• 肌肉激活通过B样条参数化扭矩分布，实现空间与时间上变化的内力矩与内力建模。
• 流体阻力与地面接触力通过黏性阻尼与刚度项建模，地面摩擦力遵循依赖于速度与方向的分段线性力律。

实验结果

研究问题

• RQ1 全柔顺杆模型在模拟具有复杂形变模式的扭转、拉伸丝状物时，其数值结果对已知解析解的逼近精度如何？
• RQ2 在受控加载条件下，该模型能否稳健模拟扭转、拉伸丝状物中扭结结构与螺旋结构的形成？
• RQ3 当与对称肌肉激活模式结合时，各向异性摩擦如何实现无肢丝状物的净运动？
• RQ4 在物理上真实的仿真环境中，实现陆地与水下无肢运动所需的最小机械与控制参数是什么？
• RQ5 该框架能否用于求解反向问题，例如识别在不同环境中最大化前进速度的最优肌肉激活模式？

主要发现

• 该数值方法能精确复现静态与动态杆构型的解析解，包括大变形与复杂边界条件的情形。
• 该模型成功捕捉了扭转、拉伸丝状物中扭结结构与螺旋结构的形成，其形变模式与理论预期一致。
• 各向异性摩擦使得在对称肌肉激活下仍能实现净前向运动，证明通过摩擦实现对称性破缺可有效驱动运动。
• 即使肌肉激活中无行波分量，系统仍因摩擦各向异性而实现净位移，验证了模型的物理解释力。
• 在对称情况下，仿真未出现虚假运动或数值噪声，证实该方案在各向同性与各向异性摩擦条件下均具有鲁棒性。
• 该方法在小时间步长下实现稳定且精确的长时间积分，始终维持能量守恒与物理保真度，适用于多样的机械与环境条件。