本文属于类型a，即单项原始研究报告。以下是对该研究的学术报告。

该研究由波兰比亚韦斯托克技术大学（Bialystok University of Technology）机械工程学院的Slawomir Grycuk和Piotr Mrozek完成，其中Slawomir Grycuk来自该校博士学院，Piotr Mrozek来自生物医学工程研究所。论文发表于《Applied Sciences》期刊，2022年4月12日出版，文章编号为3876。

在学术背景方面，该研究属于骨科矫形器生物力学与有限元分析领域。脊柱侧弯（scoliosis）是一种严重的三维脊柱畸形，波士顿矫形器（Boston brace）是常用的非手术治疗手段。此类矫形器通常较为坚硬且沉重，患者每天需佩戴超过23小时。矫形器通过三点压力系统（three-point pressure system）产生矫正力：在脊柱弯曲顶部施加力，并在对侧两个位置提供支撑。然而，目前关于矫形器本身机械特性的研究较少，大多数文献关注的是矫正力的分布，而忽视了矫形器结构如何承担和传递这些力。该研究的核心目标是建立一个经过实验验证的波士顿矫形器有限元模型（finite element method, FEM），并基于该模型分析其内部应力分布和力流线（force flow lines），从而为矫形器的结构优化提供理论依据。研究假设为：基于一个在实践中已被证明有效的真实矫形器建立可靠的数值模型，是分析矫形器机械性能的最佳起点。

研究流程包括三个主要程序。第一个程序是建立矫形器的几何模型。研究团队使用ATOS Core 200三维光学扫描仪对真实的波士顿矫形器进行逆向工程（reverse engineering）扫描。扫描前对矫形器进行了准备处理：用带子固定矫形器以增加刚度，表面喷涂消光材料，并粘贴28个参考标记点以进行多次扫描的空间配准。共采集35次曝光，生成包含3,547,641个三角形单元的STL网格。随后使用Geomagic Wrap 2017软件清理网格，去除扣件、修复伪影和填补缺损，并在处理后的外表面基础上通过NURBS样条生成了厚度为4毫米的壳体模型。在脊柱弯曲顶点水平处，还建模了内侧衬垫和软性矫正环。第二个程序是建立有限元模型并施加边界条件。矫形器外壳材料为聚丙烯（polypropylene），弹性模量E=1000 MPa，泊松比v=0.2；内衬为软泡沫，E=100 MPa，v=0.45。采用10节点四面体单元Solid187进行网格划分，并在几何复杂的关键区域进行了网格加密。网格独立性验证以von Mises应力的变化不超过3%为标准，最终得到199,187个节点和104,742个单元。边界条件按照三点力系统设置：两个固定支撑B和C，以及大小为36.25 N的预加载力A，方向与躯干对矫形器的作用一致；同时考虑了重力作用。在完成第一步计算后，在A力方向附近施加了一个0.1 N的附加力D，以对应实验中干涉仪的测量范围。由于干涉仪测量的是沿相机镜头轴线方向的位移，因此数值计算的位移结果也沿该方向输出，并取施加D力前后的位移差作为最终模拟结果。第三个程序是实验验证。研究团队搭建了基于电子散斑干涉法（electronic speckle pattern interferometry, ESPI）的测试台。ESPI系统使用功率为50 mW、波长为532 nm的DPSS激光器，光束经分束器分为参考光和物光，分别经扩束后照射参考漫射面和被测矫形器表面，两束散斑图在CMOS相机成像面上干涉。矫形器固定在专用支架上，通过带子施加36.25 N的预紧力，由压力传感器测量。随后通过滑轮系统施加0.1 N的附加载荷，记录附加载荷引起的位移增量。根据干涉条纹级数n和波长λ，利用公式Δz = n×λ/2计算离面位移（out-of-plane displacement）。

研究结果分为两个部分。第一部分是模型验证结果。有限元数值分析得到的z方向最大位移为1.606 μm，ESPI实验测量得到的最大位移为1.596 μm。实验干涉图显示了0至6级干涉条纹。两种方法在矫形器前壁中心区域的相对位移差异约为0.6%。位移等值线分布在FEM图像与ESPI干涉条纹图案之间具有明显的相似性。由于实验测量覆盖了矫形器几乎一半的前表面，并且位移符合良好，因此可以认为该验证方法足以确认模型的可信度。基于线弹性关系和胡克定律σ=Eε，位移场的符合意味着应变场和应力场也具有对应关系，从而表明该有限元模型能够以较高的保真度反映真实矫形器的力学特性。第二部分是利用已验证的数值模型进行的分析。在30 N的生理载荷水平下，von Mises应力分布显示，除施力点和支撑点等局部小区域外，矫形器大部分区域的应力低于4.4 MPa，处于聚丙烯应力-应变曲线的线性范围内，验证了模型采用线弹性假设的合理性。von Mises应力图本身难以清晰展示矫形器的整体受力工作方式。进一步的主应力（principal stress）矢量分析显示，在矫形器前壁区域，外表面主要承受压应力，而内表面承受拉应力，内表面应力的绝对值约为外表面的两倍。这种应力分布特征源于拉伸与弯曲两种载荷状态的叠加效应。基于最大拉应力σ1矢量绘制了主应力轨迹线，这些轨迹线沿应力矢量方向以光滑样条的形式连接各外力施加点。轨迹线清楚展示了力流线的分布，表明矫形器以相对较短的路径在矫正力施加点之间传递载荷。根据这些力流线，可以区分出承担矫正功能的核心结构区域和主要起稳定作用的辅助区域。

在结论方面，该研究成功建立了一个经过实验验证的波士顿矫形器有限元模型，并通过主应力矢量轨迹分析定义了矫形器的力学工作模式。研究明确了矫形器中承担矫正载荷的关键区域和几乎不参与力学工作的区域，为结构优化提供了直接依据。从科学价值看，该研究首次将主应力轨迹和力流线分析系统地应用于矫形器结构研究，弥补了现有文献中对矫形器本身力学特性分析不足的空白。从应用价值看，研究结果支持在关键承载区域通过添加碳纤维或玻璃纤维增强复合材料薄层来增加刚度，同时在其他区域减少材料以减轻重量，从而在不改变矫正力分布的前提下优化矫形器的功能性和佩戴舒适性。此外，该分析方法具有通用性，可推广至其他类型的矫形器。

研究的亮点包括：第一，采用ESPI激光干涉测量与有限元数值模拟相结合的方法，实现了矫形器离面位移的高精度验证，相对差异仅为0.6%；第二，通过绘制主应力轨迹来定义力流线，将矫形器的受力传递路径可视化，为结构优化提供了直观且可靠的依据；第三，研究对象为真实临床使用的波士顿矫形器，而非简化模型，几何保真度高；第四，研究还提出了具体的结构优化方案，包括在关键区域添加复合材料加固层、在辅助区域开孔去除材料，以及与3D打印技术结合的开孔式轻量化设计思路，具有直接的工程转化前景。