[论文解读] A new computational method for a model of C. elegans biomechanics: Insights into elasticity and locomotion performance

Netta Cohen, Thomas Ranner|arXiv (Cornell University)|Feb 16, 2017

Genetics, Aging, and Longevity in Model Organisms参考文献 40被引用 14

一句话总结

本文提出一种新颖的有限元方法，用于模拟秀丽隐杆线虫（C. elegans）的运动，通过将蠕虫建模为具有局部长度约束的不可伸长弹性梁，实现了在不同环境中对多种波状步态和体形的精确模拟。其主要贡献在于一种稳定且半隐式的数值格式，预测出有效运动所需的杨氏模量范围为10–1000 kPa，证实了通过步态调节实现速度控制是穿越阻力介质时的主要调控机制。

ABSTRACT

An organism's ability to move freely is a fundamental behaviour in the animal kingdom. To understand animal locomotion requires a characterisation of the material properties, as well as the biomechanics and physiology. We present a biomechanical model of C. elegans locomotion together with a novel finite element method. We formulate our model as a nonlinear initial-boundary value problem which allows the study of the dynamics of arbitrary body shapes, undulation gaits and the link between the animal's material properties and its performance across a range of environments. Our model replicates behaviours across a wide range of environments. It makes strong predictions on the viable range of the worm's Young's modulus and suggests that animals can control speed via the known mechanism of gait modulation that is observed across different media.

研究动机与目标

• 开发一种稳健的计算框架，用于模拟秀丽隐杆线虫在多种物理环境（包括牛顿流体和粘弹性凝胶）中的运动。
• 通过非线性初值边值问题建模波状运动的生物力学，以捕捉任意体形和步态动力学。
• 通过将材料特性与在阻力介质中的运动表现关联，估算蠕虫体部的有效弹性（杨氏模量）。
• 研究步态调节如何在不改变肌肉功率的前提下实现不同介质中的速度优化。
• 提供一种数值稳定且高效的模拟方法，适用于软体非分节运动器，具有神经-机械整合的应用前景。

提出的方法

• 该模型基于弹性梁力学，构建为非线性初值边值问题，通过基于质量守恒推导的局部长度约束引入不可伸长性。
• 开发了一种新型参数化有限元方法，采用中线曲线的分段线性逼近和分段常数张力，以实现局部不可伸长性。
• 采用半隐式时间推进策略，在每个时间步长生成线性系统，从而提升数值稳定性和效率。
• 引入曲率作为新的独立变量，使连续分段线性有限元得以应用，改善了收敛性。
• 通过借鉴Barrett等人（2010）的思路，引入人工切向运动，避免了网格畸变，确保在大变形下仍具鲁棒性。
• 该公式支持内压载荷的自然施加，并可适用于具有局部非保守距离的系统（如蠕动推进器）。

实验结果

研究问题

• RQ1 秀丽隐杆线虫在琼脂样介质和流体环境中实现有效前向推力的杨氏模量可行范围是什么？
• RQ2 在不改变肌肉功率的前提下，步态调节如何影响不同阻力介质中的运动速度？
• RQ3 计算模型能否在真实生物力学约束下准确模拟多种体形，包括卷曲和Ω形转弯？
• RQ4 局部不可伸长性与基于曲率的有限元方法在模拟软体波状运动器时如何提升数值稳定性？
• RQ5 在仅使用预设肌肉激活函数而无显式神经控制的情况下，该模型在多大程度上可预测运动表现？

主要发现

• 模型预测秀丽隐杆线虫的杨氏模量范围为10–1000 kPa，下限由产生推力的能力决定，上限由在阻力环境中避免失效的能力决定。
• 蠕虫可通过调节步态在不同介质中优化速度，证实步态是运动性能的主要调控机制。
• 该数值方法实现了对复杂体形（包括卷曲虫和环形蠕虫）的稳定模拟，这些形态对传统模拟器而言极具挑战。
• 模型预测结果与实验观测及先前的定性数据一致，同时对Sznitman等人（2010a）和Fang-Yen等人（2010）的弹性估算进行了优化。
• 该方法可准确模拟在牛顿流体和粘弹性环境中的运动，未来可扩展至内部粘度和游泳过渡的模拟。
• 该公式支持未来与神经肌肉控制系统集成，实现复杂导航与环境适应的闭环模拟。