中空圆柱体降解一维建模误差量化：薄壁几何适用性准则与查找表构建新方法

《Results in Surfaces and Interfaces》：One-Dimensional Modeling of Hollow Cylinder Degradation: A Novel Approach with Quantified Error Bounds

本文针对空心圆柱体材料降解研究中一维模型适用性边界不清的问题，首次系统量化了有限尺寸圆柱体应用一维模型引入的误差。研究通过牛顿-拉夫森法求解控制立方方程，分析了5000余种非均匀表面侵蚀工况，建立了误差查找表。结果表明当壁厚仅为长度和外径的10%时，建模误差可达-11%至+31%。该研究为可降解聚合物器件（如PHA吸管）的设计测试提供了关键理论工具，解决了"多薄才算足够薄"的核心设计难题。

　　

在材料降解研究领域，科学家们常常使用薄片状样品来简化质量损失的数学建模。理想情况下，所有的质量损失都归因于样品厚度这一个维度的减少。然而在现实中，当样品的边缘也开始腐蚀时，数学模型就会变得复杂，涉及表面侵蚀量的三次方程。对于空心圆柱体这类常见几何形状（从薄壁垫圈到吸管），情况则更为复杂。尽管出于简便，研究人员可能会采用一维降解模型，但该模型仅在壁厚趋于无限薄时才完全精确。一个长期困扰研究人员的关键问题是：对于具有实际尺寸的圆柱体，应用一维模型会引入多大误差？或者说，圆柱体的壁厚需要多薄，才能放心地使用简单的一维模型？

更普遍地说，由于一维模型将所有质量损失完全归因于物体单个维度的减少，而忽略了其他暴露表面的收缩贡献，因此一维模型总是会高估表面侵蚀量以及材料降解速率。相应地，基于一维模型预测的聚合物降解物体寿命总是短于基于全面三维模型的评估。当涉及到空心圆柱体的数学建模时，类似的问题也存在，例如圆柱体端部的侵蚀对总质量损失的贡献有多大？特别是在非均匀表面侵蚀的情况下（例如，圆柱体内壁、外壁和端面的侵蚀速率可能不同），准确建模变得极具挑战性。尽管2024年已经开发了针对三维空心圆柱体所有表面均匀侵蚀的模型，并能通过广义理论处理不同表面的差异侵蚀，但直接求解控制方程（一个三次方程）在代数上十分繁琐。如果测试样品的几何形状允许使用一维模型处理数据，那么分析和解释质量损失测量结果将大大简化。

因此，确定空心圆柱体的理想尺寸，使得简单的一维降解模型足够精确，成为一个重要的研究目标。这项研究旨在填补这一空白，为研究人员规划和进行空心圆柱体降解实验提供实用指导。本研究的主要目标是识别适用于一维降解模型的空心圆柱体理想尺寸。如果降解实验使用具有这些尺寸的测试样品，那么质量损失测量的分析和解释将大大简化，因为一维模型易于实施。此外，本研究还通过数值方法定义了任意空心圆柱体几何形状下控制三次方程的关键根，并将这些精确结果用一维近似的表面侵蚀量表示出来，提供了便于研究人员使用的查找表，避免了为特定圆柱体尺寸求解控制三次方程的挑战。

该论文发表于《Results in Surfaces and Interfaces》。

为开展研究，作者采用了关键的技术方法。首先，将控制三维空心圆柱体非均匀侵蚀的方程转化为无量纲形式。其次，明确了一维近似模型，该模型仅在壁厚与长度之比（W/L）和壁厚与外径之比（W/D_o）都趋近于零时精确。研究的核心是使用牛顿-拉夫森数值方法求解精确的三维模型控制方程，并以此为基础，系统量化一维模型相对于精确解的百分比误差。该方法涵盖了超过5000种模拟工况，参数范围包括不同的W/L和W/D_o比值，并考虑了侵蚀仅发生在内表面、外表面或均匀发生在所有表面的不同情况。最终，研究结果以误差查找表的形式呈现，并辅以线性曲线拟合公式，便于研究人员快速评估特定几何形状下一维模型的适用性及其误差范围。

4.1. 控制方程的无量纲形式

研究将描述三维空心圆柱体降解的控制方程转化为无量纲形式，其系数由圆柱体的初始尺寸（内径、外径、长度）定义。方程的解必须满足在降解完成时壁厚或长度为零的物理约束。本研究主要关注壁厚先于长度降为零的“管状”圆柱体。

4.2. 一维近似

当圆柱体的壁厚与半径和长度相比非常小时（即W/L和W/D_o都趋近于零），三维模型简化为简单的一维模型。在此模型下，表面侵蚀量与质量损失呈线性关系，计算简便，无需求解复杂方程。

5.1. 数值收敛性

通过牛顿-拉夫森法求解控制方程，显示出良好的收敛性。例如，在特定参数下，经过数次迭代，解便能快速收敛到精确值，证明了数值方法的鲁棒性。

5.2. 参数化研究

参数化研究是本文的核心。通过系统改变W/L和W/D_o，并考虑不同的表面侵蚀场景（仅内壁侵蚀、仅外壁侵蚀、所有表面均匀侵蚀），研究量化了一维模型的误差。

• •
  当侵蚀仅发生在圆柱体外表面和端面时，一维模型的误差为正值（即一维模型高估了侵蚀量），且随着W/L或W/D_o的增大而增加。
• •
  当侵蚀仅发生在圆柱体内表面时，一维模型的误差为负值（即一维模型低估了侵蚀量），且误差大小与W/D_o有关，与W/L无关。
• •
  当所有表面均匀侵蚀时，误差行为介于上述两者之间。
  研究还发现，在降解过程开始和结束时，一维模型是精确的。误差在降解初期达到最大。研究提供了详细的误差查找表，展示了在不同W/L和W/D_o组合下，一维模型可能产生的最大正误差、最大负误差以及总误差范围。

5.3. 一维模型假设引入的误差

本节明确了“一维模型何时适用”的判据。研究表明，没有绝对的答案，取决于研究人员可容忍的误差范围。本研究首次量化了这种误差，为选择测试样品形状提供了依据。误差百分比定义为（一维模型近似值 - 精确解）/ 精确解 × 100%。通过分析，研究指出当W/L和W/D_o都小于或等于0.1时，一维模型在大多数情况下的误差可能处于可接受范围内。例如，当W/L = W/D_o= 0.1时，最大误差范围在-11%到+31%之间。研究还提供了误差的线性拟合公式，便于快速估算。

5.4. 厚壁圆柱体

虽然重点在薄壁圆柱体，但研究也指出提供的工具和查找表同样适用于厚壁圆柱体（W/D_o较大）。误差等高线图显示，随着圆柱体壁厚增加或长径比减小，一维模型的误差显著增大，强调了在分析厚壁圆柱体时使用精确三维模型的重要性。

本研究为进行空心圆柱体材料降解研究的研究人员提供了实用的指导。通过量化一维模型在不同几何尺寸和侵蚀条件下的误差，并提供易于使用的误差查找表及拟合公式，使研究人员能够根据所需的精度水平，明智地选择测试样品的几何形状，或在使用一维模型时了解其误差范围。这对于可生物降解聚合物制品（如PHA吸管）的生命周期评估、生物医学领域可降解空心圆柱体器件（如药物输送系统、神经导尿管、血管支架）的性能预测，以及其他涉及圆柱形结构侵蚀的领域（如石油工业中的砂粒产生、管道侵蚀等）具有重要价值。

该研究首次提供了涵盖所有实际环境条件的、用户可访问的全面数据集，可用于空心圆柱形部件的设计、建模和测试阶段。尽管生命周期评估本身是一个重要领域，但本文所采用的数学模型是目前对任意尺寸空心圆柱体降解的唯一可靠表征，填补了三维形状生命周期评估分析中的一个重要空白。研究也指出，未来需要更多关于圆柱体内外表面相对降解速率的研究，以进一步完善模型在实际环境中的应用。