肛瘘封堵器的材料设计与生物力学研究.pdf

文章编号：1007 6735(2023)05 0488 07DOI:10.13255/ki.jusst.20220119004肛瘘封堵器的材料设计与生物力学研究刘颖，崔海坡，宋成利，李来福，郎景成，姜赫宣（上海理工大学教育部微创医疗器械工程中心，上海200093）摘要：针对肛瘘的手术治疗，设计一款可降解的肛瘘封堵器，主要利用高温熔融法制备不同分子量、不同基团的左旋聚乳酸试件，通过拉伸试验获得材料的应力应变曲线，确定了封堵器的最佳制作材料。利用有限元仿真分析软件 ABAQUS，系统地研究了封堵器在植入过程中的应力应变分布以及径向支撑力等参数的变化。结果表明：分子量为 110 万、基团为酯封端的左旋聚乳酸弹性模量值最高；封堵器在封堵过程中受到的最大应力没有超过其许用应力（4.55kPa），初始接触阶段应力仅集中在封堵器周围区域，并且距离接触中心越远，应力值越大；随着封堵器的不断侵入，接触应力以环形波的方式从接触中心向四周扩散；封堵器基本保持了初始形状，没有发生明显的挤压变形，在封堵过程中具有良好的封堵效果且不易被破坏。关键词：左旋聚乳酸；有限元分析；封堵器；肛瘘中图分类号：TB333文献标志码：AMaterial design and biomechanical study onsealing device for anal fistulaLIU Ying，CUI Haipo，SONG Chengli，LI Laifu，LANG Jingcheng，JIANG Hexuan(Shanghai Institute for Minimally Invasive Therapy,University of Shanghai for Science and Technology,Shanghai 200093,China)Abstract:Aimingatthesurgicaltreatmentofanalfistula,adegradableanalfistulasealingdevicewasdesigned.L-polylactic acid specimens with different molecular weights and different groups werepreparedbyhigh-temperaturemeltingmethod.Thestress-straincurveofthematerialwasobtainedbytensile test,and the best material for the sealing device was determined.Finite element simulationanalysissoftwareABAQUSwasusedtosystematicallystudythestress-straindistributionandradialsupportforceduringthesealingprocess.TheresultsshowedthattheelasticmodulusofL-polylacticacidisthehighestwiththemolecularweightof1.1millionandthegroupisester-sealed.Themaximumstressofthesealingdeviceintheprocessofsealingdidnotexceeditsallowablestress4.55kPa.Intheinitialcontactstage,thestresswasonlyconcentratedintheareaaroundthesealingdevice,andthefurtherawayfromthecontactcenter,thegreaterthestressvaluewas.Thecontactstressdiffusedfromthecontactcentertoallsidesintheformofringwavewiththecontinuousintrusionofthesealing上 海 理 工 大 学 学 报第45卷第5期J.UniversityofShanghaiforScienceandTechnologyVol.45No.52023收稿日期：20220119基金项目：科技部国家重点研发计划（2019YFC0120402）；上海市自然科学基金资助项目（19ZR1435300）第一作者：刘颖（1997），女，硕士研究生研究方向：材料设计与有限元仿真E-mail：通信作者：崔海坡（1978），男，教授研究方向：材料设计与应用E-mail：h_b_device.Thesealingdevicebasicallymaintainedtheinitialshapeanddidnothaveobviousextrusiondeformation.Itwouldnotbecrushedinthesealingprocessandhadagoodsealingeffect.Keywords:L-polylactic acid;finite element analysis;sealing device;anal fistula肛瘘是临床上较为常见和棘手的肛肠科疾病之一，发生于肛门直肠周围皮肤之间的病理性管道。大约有 33%的肛瘘由肛周脓肿继发而形成，其临床症状、治疗方式等相关情况都具有其特殊性，不能自愈，必须手术治疗。我国肛瘘发病率为 1.67%3.60%，中青年好发，男性是女性的26 倍，其中复杂性肛瘘占 5%15%，为一种难治的良性疾病，治疗不彻底有恶变倾向1-2。肛瘘治疗的关键在于正确找到内口，尽可能减少肛门括约肌的损伤，最大限度保护及维持肛门的功能，减少并发症的发生，降低复发率3。对于肛瘘的治疗多以手术为主，包括肛瘘切开术、肛瘘挂线术、括约肌保留术、肛瘘栓、纤维蛋白胶封堵术等。这些手术都会对肛门括约肌造成不可修复的损伤，且术后容易复发，复发风险在 10%57%。现阶段肛瘘治疗已经逐渐从创伤较大的开放手术转变为微创手术，从切开性手术转变为恢复性手术4-7。但是，对于复杂性肛瘘治疗、术后复发、愈合缓慢等临床难题，仍然没有很好的解决方案和治疗器械。现有的组织闭合器大多采用不可被人体吸收的材料所制作，如在做肛瘘内口闭合时采用不锈钢钉使切口相吻合，或采用钛合金闭合夹对肛瘘内口进行闭合。钢钉或闭合夹均需二次手术取出，增加了患者的痛苦和手术风险，且对病人的生活质量造成较大影响。左旋聚乳酸（PLLA）具有良好的生物降解性和生物相容性，且易于加工成型，已被美国食品和药物管理局批准用于生物医学的各个领域8-9。左旋聚乳酸作为内固定材料最大的优势在于刚性小，无应力遮挡作用，植入人体能够很好地被吸收，不需要二次手术取出，实现体内无异物的治疗效果，近年来已成为研究的热点10-11。基于上述分析，本文以左旋聚乳酸为原材料，设计一款新型可降解肛瘘封堵器。首先，通过实验确定封堵器最佳制作材料，然后利用有限元仿真分析技术，系统研究封堵器在植入过程中的应力应变分布以及径向支撑力等参数的变化，为肛瘘内口封堵器的设计和优化提供参考，从而为临床上肛瘘疾病的治疗提供新型的手术方案和相应理论依据。1 材料与方法 1.1 材料与实验本文设计的封堵器采用左旋聚乳酸材料制作，原材料购自山东省济南市岱罡生物聚乳酸有限公司。为选择最佳力学性能的制作材料，并为后续有限元仿真分析输入材料参数，首先针对3 种不同分子量以及 3 种不同基团的左旋聚乳酸，开展了拉伸实验。3 种型号左旋聚乳酸 DG-L200、DG-L400、DG-L600 的分子量分别为 22 万、55 万、110 万，3 种基团分别为酯封端、端羟基、端羧基。在分析分子量对材料力学性能影响规律时，固定基团为酯封端；在分析基团对材料力学性能影响规律时，固定分子量为 110 万。通过高温熔融法制备哑铃状试件，其尺寸如图 1 所示。根据国家标准高分子材料拉伸性能试验方法：GB/T1040.12006开展单向拉伸实验，拉伸速度为10mm/min。53.010.0R7.05.01.0单位：mm63.7100.0图 1 试件尺寸图Fig.1 Dimension of specimen 1.2 结果分析经拉伸实验测得的左旋聚乳酸系列材料的应力应变曲线如图 2 所示，弹性模量如表 1 所示，弹性模量低意味着在较小的外力作用下会发生较大的变形12。对于瘘道封堵器而言，其制作材料的弹性模量越高，则封堵器在使用过程中的变形越小，结构的稳定性就越好，从而能够支撑瘘道的挤压和变形，达到封堵的效果。根据表 1 所知，随着材料分子量的增加，其弹性模量呈现增大趋势，即分子量为 110 万的左旋聚乳酸具有较高的弹性模量。而对于酯封端、端羟基、端羧基第5期刘颖，等：肛瘘封堵器的材料设计与生物力学研究4893 种基团而言，酯封端左旋聚乳酸刚度更大，具有更佳的抗变形能力。根据实验结果，选择分子量为 110 万、基团为酯封端的左旋聚乳酸作为瘘道封堵器的制作材料。3045DG-I 200DG-L400DG-L60015应力/MPa0153045607590酯封端端羟基端羧基应力/MPa00.010.020.03应变(a)酯封端不同分子量试件(b)不同基团 DG-L600 试件0.010.02应变0.030.040.05图 2 左旋聚乳酸系列材料的应力应变曲线Fig.2 Stress-strain curves of L-polylactic acid seriesmaterials 表 1 左旋聚乳酸系列试件弹性模量Tab.1 Elastic modulus of L-polylactic acid series specimens样品酯封端DG-L200酯封端DG-L400酯封端DG-L600端羟基DG-L600端羧基DG-L600弹性模量E/MPa16811822203018991700 2 有限元模型 2.1 几何模型2.1.1封堵器模型肛瘘内口封堵器以线圈状植入物为主体，结构简单，无额外扭矩，与瘘道接触时容易定位且能够分散瘘道内压的冲击。将这种可吸收的线圈状植入物包裹在肛瘘的内口周围，可抵抗瘘道内口的周向压缩，在保证瘘道愈合的同时保留了括约肌功能，通过肛瘘内口封堵器的推送装置将植入物旋进正确的深度，自动调节患者之间的差异，最终达到治愈的目的。图 3 为封堵器的尺寸图，单位为 mm。1.257.001.203.001.80R2.901.100.521.983.024.485.006.254.25图 3 封堵器模型Fig.3 Model of sealing device2.1.2肌肉组织和瘘道模型肌肉组织由于没有特定的形状，因此可以忽略肌肉组织形状对于封堵过程的影响13，在建模时将其设计为一个规则的圆柱体。考虑到对称性，将肌肉组织模型简化为四分之一圆弧。由于只研究封堵器在局部瘘道进行封堵过程中的封堵效果，不考虑复杂瘘道以及多瘘道的封堵，故将瘘道简化为一个长圆柱体。2.2 材料参数在 ABAUQS2016 中进行有限元静力学分析时，肛瘘内口封堵器材料为上述试验确定的酯封端 DG-L600（分子量为 110 万）的左旋聚乳酸，将材料视为各向同性、连续均匀的线弹性材料。封堵器及肌肉组织的材料参数如表 2 所示。表 2 封堵器与肌肉组织的材料参数Tab.2 Material parameters of the sealing device and tissue部件材料模型弹性模量/MPa泊松比密度/（g.cm3）封堵器线弹性20300.301.24肌肉组织超弹性2.530.491.30对于肌肉组织的材料模型，综合考虑几类常见的超弹性模型，如Neo-Hooken14，Ogden，Mooney-Rivlin 和 Yeho15等，选择 Ogden 模型，其应变能函数模型采用幂级数表达形式，更易于拟合各向同性的非线性材料。而且肌肉组织材料应用三阶Ogden 超弹性模型，最大的特殊性在于理论算法上采用了主拉伸作为基准变量，而不是常见的应变张量不变量16。Ogden 模型中包括 2 个参数，分别为剪切模量 与应变硬化指数，这两个参数能够较好地反映外应力作用下肌肉组织硬化的力学性能。Ogden 模型的应变能函数如式（1）所示：490上海理工大学学报2023年第45卷W=Ni=1ii(i1+i2+i33)+Nk=11Dk(j1)2k式中：N 为模型的阶数，通常取 13 之间的数；和 为材料参数，代表压力，为无纲量参数；D 为不可压缩参数，表示体积变化。为了达到所需要的精度要求，保证在应变很高时还能很好地拟合数据，选用 N=3 的超弹性Ogden 模型应用于肌肉组织材料17。和 的取值如表 3 所示18。表 3 三阶 Ogden 模型参数Tab.3 Third order Ogden model parameters1/MPa12/MPa23/MPa30.00837.6250.000213.8750.00627.625 2.3 网格划分静力学仿真分析中封堵器采用二阶四面体网格，网格类型为 C3D10，单元数为 7504，节点数为12580。显式动力学仿真分析中采用线性四面体网格，网格类型为 C3D4，封堵器共划分为 175005个单元数，节点数为34604，肌肉组织共划分为89173个单元，节点数为 18348。2.4 载荷与边界条件2.4.1支撑性能分析进行支撑性能测试时，为了防止刚体位移，在封堵器主体上约束柱坐标下进行 Z 向位移，同时在周向对称面上施加压力载荷。根据周旺伟19的研究，手术前的瘘道静息压为（9.821.04）kPa，手术后的瘘道静息压为（11.132.20）kPa，理论上当封堵器的径向支撑力超过 4.55kPa 时视为被压坏，故本文仅分析封堵器植入人体后的力学性能，不考虑植入过程时血液和其他体液对封堵器的影响。对封堵器线圈状主体表面施加一个 4.55kPa的径向压力模拟瘘道外载荷的作用，通过计算封堵器在该载荷下的径向位移来分析封堵器的径向支撑性能。为了模拟封堵器植入人体的旋转过程，在封堵器的 4 个四边形截面上施加 1Nm 的力矩，根据右手定则为逆时针方向。同时，在旋进的一端施加一个 0.0025N 的载荷（f=ps，其中 p=9.82kPa，s=0.25mm2）来模拟封堵器植入过程的阻力。2.4.2封堵过程分析肛瘘内口封堵是一个大变形的过程，需考虑封堵器和组织的接触。根据蒋志航13的研究，当两个物体表面相互接触时，选择通用的罚摩擦参数进行滑移摩擦行为的判定，可以顺利完成整个封堵过程中的仿真分析。故本文在参数设定时选择罚接触的相互作用，封堵器与肌肉组织的摩擦系数为 0.15。在封堵过程中，封堵器与肌肉组织之间会产生相对滑移，因此选择动态接触方式的面面接触，将封堵器线圈状主体的外周设置为主接触面，组织接触的面设置为副接触面。对于边界条件的设置，约束封堵器全部 6 个自由度作为边界条件。在封堵器旋转轴上选择一个参考点，将其与封堵器整体建立耦合约束（distributioncoupling）。基于初始边界条件，对肌肉组织四分之一圆弧的外周施加固定约束。3 结果分析 3.1 静力学分析由于肛瘘封堵器大多用于肛瘘瘘道的愈合封堵，不仅要承受瘘道的变形，还会受到外力的挤压。封堵器仿真分析结果如图 4 所示，图中，表示应力，表示应变。从图中可以看出，封堵器在封堵过程中最大应力为 6.398109kPa，远小于其许用应力 4.55kPa，且总体变形也较小。因此，封堵器在封堵过程中能够承受瘘道的挤压，满足瘘道的支撑要求，从而达到封堵的效果。3.2 动力学分析图 5 为封堵器与组织接触过程中不同时刻的(a)应力云图(b)应变云图/MPa6.398e125.687e124.976e124.265e123.554e122.843e122.133e121.422e127.109e134.333e29/mm4.474e173.977e173.480e172.983e172.486e171.989e171.491e179.943e184.971e180图 4 封堵器仿真分析结果Fig.4 Simulation analysis results of sealing device第5期刘颖，等：肛瘘封堵器的材料设计与生物力学研究491等效变形云图。其中，图 5(a)接触时间为 0s，即封堵器未与组织接触，相应变形也为 0；图 5(b)(f)为封堵器与组织接触过程中，接触时间分别为0.05，0.15，0.25，0.35，0.45s 时，封堵器的变形图。从图 5(b)(f)中可以看出，当封堵器一端接触到组织边界时，模型开始发生变形。在初始接触阶段应力仅集中在封堵器三角状圆弧凸起区域，距离接触中心越远，应力值越大，且随着封堵器逐渐侵入组织，封堵器的变形也越来越大。接触时间为 0.15s 时，变形达到了 8.71mm；接触时间超过 0.15s 后，应力逐渐向封堵器中心传递，封堵器的线圈状主体承受一部分压力，开始发生变形。当接触时间为 0.35s 时，封堵器线圈状主体承受的应力达到了最大值，此时封堵器变形为7.54mm；接触时间为 0.45s 时，封堵器变形为3.27mm。封堵器变形有所减小是因为随着封堵器的不断侵入，接触应力以环形波的方式从接触中心向四周扩散开，分散了封堵器所承受的压力。图 5 显示，在封堵过程中，封堵器基本保持了初始形状，没有发生明显的挤压变形，表明封堵器在封堵过程中有良好的封堵效果且不至被破坏。图 6 为组织与封堵器接触过程中不同时刻组织的等效变形云图。其中，图 6(a)接触时间为 0s，即组织未与封堵器接触，相应变形也为 0；图 6(b)(f)为组织与封堵器接触过程中，接触时间分别为0.05，0.15，0.25，0.35，0.45s 时，组织的变形云图。从图 6(b)(f)中可以看出，在封堵器逐渐与组织接触过程中，组织模型从接触中心向四周以环形波的形式发生变形。随着接触时间的增加，组织的变形逐渐增大；当接触时间小于 0.15s时，组织的轮廓并没有发生明显的变形，此时应力主要集中于与封堵器接触的部分；当接触时间为 0.45s 时，组织的变形达到最大值，为 8.99mm，此时组织轮廓发生了肉眼可见的破坏变形。图 7 为封堵器与组织接触过程中不同时刻的等效应力云图。其中：图 7(a)接触时间为 0s，即封堵器未与组织接触，相应应力为 0；图 7(b)(f)为封堵器与组织接触过程中，接触时间分别为0.05，0.15，0.25，0.35，0.45s 时，封堵器的应力云图。由图 7(b)(e)可以看出，随着接触时间的增加，封堵器受到的应力逐渐增大，当接触时间为 0.35s 时，封堵器主体承受的应力达到了最大(a)0 s(b)0.05 s(c)0.15 s(d)0.25 s(e)0.35 s(f)0.45 s/mm0000000000/mm4.0583.6073.1562.7052.2541.8041.3530.9020.4510/mm7.0086.2295.4514.6723.8933.1152.3361.5570.7790/mm8.7147.7466.7785.8094.8413.8732.9051.9360.9680/mm3.2702.9052.5422.1791.8161.4531.0890.7260.3630/mm7.5416.7035.8655.0274.1903.3522.5141.6760.8380图 5 与组织接触不同时刻封堵器变形云图Fig.5 Deformation diagram of sealing device at different time of contact with tissue492上海理工大学学报2023年第45卷(a)0 s(b)0.05 s(c)0.15 s(d)0.25 s(e)0.35 s(f)0.45 s/mm0000000000/mm1.0000.8890.7780.6670.5560.4450.3330.2220.1110/mm5.0004.4443.8893.3332.7782.2221.6671.1110.5560/mm3.0002.6672.3332.0001.6671.3331.0000.6670.3330/mm8.9887.9896.9905.9924.9933.9942.9961.9970.9990/mm6.9896.2135.4364.6603.8833.1062.3301.5530.7770图 6 与封堵器接触不同时刻组织变形云图Fig.6 Deformation diagram of tissue at different time of contact with sealing device/kPa0000000000/kPa0.2930.2820.2610.1960.1630.1300.0980.0650.0330/kPa1.8191.6171.4151.2131.0110.8090.6060.4040.2020/kPa1.1791.0480.9170.7860.6550.5240.3930.2620.1310/kPa1.1411.0140.8870.7600.6340.5070.3800.0250.0130/kPa2.2281.9811.7331.4861.2380.9900.7430.4950.2480(a)0 s(b)0.05 s(c)0.15 s(d)0.25 s(e)0.35 s(f)0.45 s图 7 与组织接触不同时刻封堵器应力云图Fig.7 Stress distribution of sealing device at different time of contact with tissue第5期刘颖，等：肛瘘封堵器的材料设计与生物力学研究493值，为 2.228kPa，仍小于材料的许用应力 4.55kPa；随后应力值有所减小，当接触时间为 0.45s 时，最大应力值为 1.141kPa。应力值有所减小是因为随着封堵器的不断侵入，接触应力以环形波的方式从接触中心向四周扩散开，分散了封堵器所承受的压力。图中显示，在封堵过程中，封堵器受到的最大应力小于其许用应力，表明封堵器在封堵过程中不至被破坏。4 讨论肛瘘一旦形成，极难自愈，现今手术仍为肛瘘治疗最为有效的方式20。由于肛瘘的管道通过外括约肌和肛提肌，有时管道在内括约肌之间走行，传统手术中会切开管道，必然要切开有关的括约肌。特别是通过肛门直肠环的瘘管21，切开肛门直肠环，术后会引起肛门失禁并且极易感染，对患者日后的正常生活造成极大影响。本文设计了一款新型肛瘘内口封堵器，可采用微创的理念完成手术植入。封堵器选用左旋聚乳酸，其刚性小，无应力遮挡作用，植入人体后能够很好地被吸收，不需要二次手术取出，达到体内无异物的治疗效果，大大缩短手术时间，且术后护理安全简单，能够较大程度地改善患者的生活质量。本文利用有限元数值仿真分析技术，模拟了封堵器封堵肛瘘的过程。结果表明，封堵器在封堵过程中不会被压坏，具备良好的封堵效果。但是在封堵过程中，组织发生了肉眼可见的破坏变形。组织被破坏后会极大地影响患者的术后恢复，是一个比较棘手的问题，但是相较于传统手术的直接切开，本文所设计的肛瘘封堵器对括约肌造成的损伤是可以接受的。由于本文只研究封堵器在局部瘘道进行封堵过程中的封堵效果，没有考虑复杂瘘道以及多瘘道的封堵，所以瘘道模型还不能解决复杂瘘道问题。另外，理论上当封堵器的径向支撑力超过 4.55kPa 时视为被压坏，故本文仅模拟了封堵器植入人体后的力学性能，没有考虑植入过程时血液和其他体液对封堵器的影响，因此，后续还需要针对植入过程中人体生理环境对封堵效果的影响作进一步深入研究。参考文献：安阿玥.肛肠病学 M.北京:人民教育出版社,2005.1谷云飞,陈红锦,史仁杰,等.保留括约肌挂线法治疗复2杂性肛瘘的临床研究 J.南京中医药大学学报,2007,23(1):2023.GILMORE B,JACKSON K L,MIGALY J.NewinnovationsinanalfistulasurgeryJ.SeminarsinColonandRectalSurgery,2019,30(4):100707.3何耀彬,方镇国,季英,等.LIFT-plug 微创技术在肛瘘治疗中的临床研究 J.现代诊断与治疗,2019,30(1):2122.4YAOYB,XIAOCF,WANGQT,etal.VAAFTplusFiLaCTM:a combined procedure for complex analfistulaJ.Techniques in Coloproctology,2021,25(8):977979.5ANDREOUC,ZEINDLERJ,OERTLID,etal.Longtermoutcome of anal fistula a retrospective studyJ.ScientificReports,2020,10(1):6483.6SHIMOJIMAY,MATSUSHIMAM,MATSUSHIMAS,etal.Anovelsurgicaltechniqueforanalfistulasurgerydesigned to preserve the anal sphincter function andanodermJ.JournaloftheAnus,RectumandColon,2021,5(1):2533.7KOHS,LEES,LEED,etal.Mechanicalpropertiesandbioactivity of poly(lactic acid)composites containingpoly(glycolic acid)fiber and hydroxyapatite particlesJ.Nanomaterials,2021,11(1):249.8PITJAMIT S,THUNSIRI K,NAKKIEW W,et al.Thepossibility of interlocking nail fabrication from FFF 3Dprinting PLA/PCL/HA composites coated by local silkfibroin for canine bone fracture treatmentJ.Materials,2020,13(7):1564.9LICM,GUOCC,FITZPATRICKV,etal.Designofbiodegradable,implantabledevicestowardsclinicaltranslationJ.Nature Reviews Materials,2019,5(1):6181.10HUXL,MIS,LUJL,etal.In vitrodegradationbehaviorof shape memory PLLA-TMC random copolymersJ.ColloidsandSurfacesA:PhysicochemicalandEngineeringAspects,2021,615:126220.11谢兆丰,王树水,张智伟,等.新型生物可吸收房间隔缺损封堵器的生物相容性研究 J.北京生物医学工程,2016,35(6):582587.12蒋永航.肌肉组织穿刺过程的有限元仿真与实验研究 D.济南:山东大学,2020:2632.13RAD M H G,SHAHABIAN F,HOSSEINI S M.Geometricallynonlinearelastodynamicanalysisofhyper-elasticNeo-HookenFGcylindersubjectedtoshockloadingusing MLPG methodJ.Engineering Analysis withBoundaryElements,2015,50:8396.14 下转第 502 页 494上海理工大学学报2023年第45卷structuralidentificationwithoutput-onlymeasurementsJ.Measurement,2018,122:597610.王晓燕,黄维平,李华军.地震动反演及结构参数识别的EKF 算法 J.工程力学,2005,22(4):2023,111.12LEIY,JIANGYQ,XUZQ.Structuraldamagedetectionwith limited input and output measurement signalsJ.Mechanical Systems and Signal Processing,2012,28:229243.13雷鹰,周欢.有限观测下的结构损伤实时在线诊断 J.振动与冲击,2014,33(17):161166.14Khodarahmi M,Maihami V.A review on kalman filtermodelsJ.Archives of Computational Methods inEngineering,2023,30(1):727747.15KALMAN R E.A new approach to linear filtering and16predictionproblemsJ.JournalofBasicEngineering,1960,82(D):3545.欧进萍.结构振动控制主动、半主动和智能控制M.北京:科学出版社,2003.17KALMAN R E.A new approach to linear filtering andpredictionproblemsJ.TransactionsoftheASME-JournalofBasicEngineering,1960,82:3545.18CHUICK,CHENGR.Kalmanfiltering:withreal-timeapplicationsM.3rded.Berlin:Springer,1999:5558.19彭斌,汪澜涯,王冬冬,等.砌体墙抗震性能伪静力试验研究 J.上海理工大学学报,2016,38(4):402408.20王冬冬.砌体结构性能的多尺度数值模拟 D.上海:上海理工大学,2017.21（编辑：丁红艺）上接第 494 页 张琦,时剑文,索双富,等.基于 Mooney-Rivlin 模型和Yeoh 模型的橡胶材料有限元分析 J.合成橡胶工业,2020,43(6):468471.15AYUOL P O M.Investigation on muscle mechanicalproperties under compression loading for human bodymodel improvementD.Changsha:Hunan University,2020:3536.16田良.心脏运动过程超弹性材料建模与仿真研究 D.西安:西安工业大学,2019:2223.17TRNTN,NOVCEKV,TOLBAR,etal.ExperimentalandcomputationalapproachtostudycolorectalanastomosisC/ProceedingsoftheXXIIICongressofthe18InternationalSocietyofBiomechanics.Brussels,Belgium,2011:12.周旺伟.肛肠手术前后肛管直肠内压变化的观察 J.广西中医药,1997,20(3):78.19GIAMUNDO P,DE ANGELIS M.Treatment of analfistulawithFiLaC:resultsofa10-yearexperiencewith175 patientsJ.Techniques in Coloproctology,2021,25(8):941948.20LIMURAE,GIORDANOP.ModernmanagementofanalfistulaJ.WorldJournalofGastroenterology,2015,21(1):1220.21（编辑：丁红艺）502上海理工大学学报2023年第45卷