CT影像数据建立吞咽肌骨三维模型以及生物力学特征分析

doi:10.12307/2025.796

中国组织工程研究 ›› 2025, Vol. 29 ›› Issue (29): 6167-6163.doi: 10.12307/2025.796

• 骨组织构建 bone tissue construction • 上一篇下一篇

CT影像数据建立吞咽肌骨三维模型以及生物力学特征分析

方丹丹1，马睿杰2，黄怿3，何克林2，吴磊2

1浙江中医药大学附属第三医院浙北分院(湖州市吴兴区中医院)，浙江省湖州市 313002；2浙江中医药大学附属第三医院，浙江省杭州市 310005；3浙江中医药大学第三临床医学院，浙江省杭州市 310053

收稿日期:2024-09-10 接受日期:2024-11-12 出版日期:2025-10-18 发布日期:2025-03-01
通讯作者: 吴磊，博士，副主任中医师，浙江中医药大学附属第三医院，浙江省杭州市 310005
作者简介:方丹丹，女，1996年生，安徽省铜陵市人，汉族，2023年浙江中医药大学毕业，硕士，中医师，主要从事针灸防治神经系统疾病的研究。
基金资助:
浙江省中医药科技计划项目(2022ZB198)，项目负责人：吴磊；浙江省中医药现代化专项项目(2020ZX011)，项目负责人：马睿杰；浙江省中医药科技计划项目(2025ZX232)，项目负责人：方丹丹

Three-dimensional model of swallowing musculoskeletal system based on CT image data and biomechanical characteristics analysis

Fang Dandan1, Ma Ruijie2, Huang Yi3, He Kelin2, Wu Lei2

1Zhebei Branch of Affiliated Third Hospital of Zhejiang Chinese Medical University (Wuxing District Hospital of Traditional Chinese Medicine), Huzhou 313002, Zhejiang Province, China; 2Affiliated Third Hospital of Zhejiang Chinese Medical University, Hangzhou 310005, Zhejiang Province, China; 3Third School of Clinical Medicine, Zhejiang Chinese Medical University, Hangzhou 310053, Zhejiang Province, China

Received:2024-09-10 Accepted:2024-11-12 Online:2025-10-18 Published:2025-03-01
Contact: Wu Lei, MD, Associate chief physician, Third Affiliated Hospital of Zhejiang Chinese Medical University, Hangzhou 310005, Zhejiang Province, China
About author:Fang Dandan, MS, Physician, Zhebei Branch of Affiliated Third Hospital of Zhejiang Chinese Medical University (Wuxing District Hospital of Traditional Chinese Medicine), Huzhou 313002, Zhejiang Province, China
Supported by:
Zhejiang Province Traditional Chinese Medicine Science and Technology Plan Project, No. 2022ZB198 (to WL); Special Traditional Chinese Medicine Modernization Project of Zhejiang Province, No. 2020ZX011 (to MRJ); Traditional Chinese Medicine Science and Technology Plan Project of Zhejiang Province, No. 2025ZX232 (to FDD)

摘要/Abstract

摘要：

文题释义：
吞咽障碍：是指口腔、咽喉、食管结构和/或功能受损，导致机体将食物及唾液送至胃内的过程出现了异常，包括解剖学改变、感觉控制功能及动力学障碍等。
有限元分析：对研究对象进行数学化模拟，将一个整体分割成小单元的组合体，其实质是利用有限数量的模拟量去模仿无限单元的真实结构，将各形质构造数学化进行计算。

背景：目前对于吞咽障碍的临床研究较多，而吞咽障碍状态对骨骼及舌骨肌群应力影响的研究较少。
目的：建立吞咽肌骨三维模型，探讨吞咽运动中舌骨及舌骨上肌群的应力变化。
方法：选取1名健康成年男性并拍摄头颈部CT，通过影像数据提取法构建三维模型。根据舌骨的位移情况，设定舌骨的最低位置及最高位置，依次加载舌骨上肌群的肌力，观察生理状态下舌骨上肌群对舌骨的应力影响。通过改变肌肉的材料属性，模拟不同程度吞咽障碍，观察吞咽异常状态下舌骨上肌群的应力变化。
结果与结论：①建立了包含C2-C4椎体、部分颞骨、下颌骨、舌骨及下颌舌骨肌-颏舌骨肌的有限元模型，可反映患者吞咽时舌骨位移情况及应力变化。②生理状态下，舌骨上肌群中下颌舌骨肌对舌骨的应力影响最大，其次是二腹肌、颏舌骨肌及茎突舌骨肌。③通过模拟不同程度的吞咽障碍，发现随着吞咽障碍严重程度增加，舌骨上肌群的应力也逐渐增加。④结果说明，舌骨上肌群中下颌舌骨肌及二腹肌对舌骨应力影响较大，提示特定肌肉刺激的重要性。吞咽障碍的加重与舌骨上肌群的应力改变呈正相关，该结果为吞咽障碍的机制研究及临床诊疗提供了生物力学参考。
https://orcid.org/0000-0001-7055-7865(吴磊)

中国组织工程研究杂志出版内容重点：组织构建；骨细胞；软骨细胞；细胞培养；成纤维细胞；血管内皮细胞；骨质疏松；组织工程

关键词: 舌骨, 舌骨上肌群, 吞咽障碍, 生物力学, 有限元, 工程化骨组织

Abstract: BACKGROUND: At present, there are many clinical studies on dysphagia, but little is known about the stress of bones and hyoid muscles caused by dysphagia.
OBJECTIVE: To establish a three-dimensional model of swallowing musculoskeletal system to investigate the stress changes of hyoid and suprahyoid muscles during swallowing.
METHODS: A healthy adult male was selected for head and neck CT scanning. A three-dimensional model was constructed by image extraction method. According to hyoid motion, the lowest and highest positions of hyoid were set, and the muscle strength of suprahyoid muscles was successively loaded to observe the stress influence of suprahyoid muscles on hyoid under physiological state. By changing the material properties of muscles, we simulated different degrees of dysphagia and observed the stress changes of suprahyoid muscle group under the condition of dysphagia.
RESULTS AND CONCLUSION: (1) The finite element models including C2-C4 vertebrae, partial temporal bone, mandible, hyoid, and mylohyoid-geniohyoid muscle were established, which could reflect the displacement of hyoid and stress changes during swallowing. (2) Under physiological conditions, among suprahyoid muscles, mylohyoid muscle had the greatest influence on hyoid stress, followed by digastric muscle, geniohyoid muscle, and stylohyoid muscle. (3) By simulating different degrees of dysphagia, it was found that with the increase of the severity of dysphagia, the stress of suprahyoid muscles gradually increased. (4) These findings conclude that among suprahyoid muscles, mylohyoid muscle and digastric muscle have the greatest influence on hyoid stress, which suggests the importance of specific muscle stimulation. The aggravation of dysphagia was positively correlated with the stress change of suprahyoid muscles, which could provide a biomechanical reference for the mechanism study and clinical diagnosis and treatment of dysphagia.

Key words: hyoid, suprahyoid muscles, dysphagia, biomechanics, finite element, engineered bone tissue

中图分类号:

方丹丹, 马睿杰, 黄怿, 何克林, 吴磊. CT影像数据建立吞咽肌骨三维模型以及生物力学特征分析[J]. 中国组织工程研究, 2025, 29(29): 6167-6163.

Fang Dandan, Ma Ruijie, Huang Yi, He Kelin, Wu Lei. Three-dimensional model of swallowing musculoskeletal system based on CT image data and biomechanical characteristics analysis[J]. Chinese Journal of Tissue Engineering Research, 2025, 29(29): 6167-6163.

图/表（结果） 3

参考文献

[1] 刘悦文, 郭琪, 于莹. 老年人口咽吞咽障碍的筛查工具及康复治疗方法的研究进展[J].中国康复医学杂志,2020,35(3):361-365.
[2] CHRISTMAS C, ROGUS-PULIA N. Swallowing Disorders in the Older Population. J Am Geriatr Soc. 2019;67(12):2643-2649.
[3] 刘雅鑫, 蒋运兰, 黄孝星, 等. 中国老年人吞咽障碍患病率的Meta分析[J].中国全科医学,2023,26(12):1496-1502,1512.
[4] 何君芳, 任亚子, 白佳佳, 等.吞咽障碍治疗相关研究计量学分析[J].西部中医药,2018,31(6):69-72.
[5] 袁盈, 蔡向红, 陈枫, 等.天突深刺治疗中风后吞咽障碍临床疗效观察[J]. 针刺研究,2019,44(1):47-50.
[6] 李眺, 陈卓铭, 金小千, 等.基于超声可视化测量舌骨-下颌骨运动评估针刺治疗脑梗死后吞咽障碍疗效的临床研究[J].中国康复医学杂志,2023, 38(9):1247-1251.
[7] 刘潇, 张彬, 任力娟, 等.女性盆底功能障碍性疾病中肛提肌的有限元应力分析[J].中华临床医师杂志(电子版),2017,11(9):1609-1612.
[8] 陈博韬. 足踝三维有限元模型建立与分歧韧带损伤的有限元研究[D].广州:南方医科大学,2022.
[9] 曹子君, 王芳, 何耀广, 等.糖尿病患者足底压力和鞋垫减压结构的有限元分析[J].中国康复理论与实践,2021,27(7):852-858.
[10] 詹民民. 基于粘弹理论的咀嚼吞咽过程生物力学研究[D].无锡:江南大学, 2018.
[11] 王兴林.吞咽障碍的生物力学变化及电刺激治疗机制[J].中华物理医学与康复杂志,2013,35(12):938-940.
[12] Huang Q, Zhang C, Bai H, et al. Biomechanical evaluation of two modified intramedullary fixation system for treating unstable femoral neck fractures: A finite element analysis.Front Bioeng Biotechnol. 2023;11:1116976.
[13] Warren JL, Yoo JE, Meyer CA, et al. Automated finite element approach to generate anatomical patient-specific biomechanical models of atherosclerotic arteries from virtual histology-intravascular ultrasound. Front Med Technol. 2022;4:1008540.
[14] Kulczyk T, Rychlik M, Lorkiewicz-Muszyńska D, et al. Computed Tomography versus Optical Scanning: A Comparison of Different Methods of 3D Data Acquisition for Tooth Replication. Biomed Res Int. 2019;2019:4985121.
[15] 周宁, 杨帆, 时亚松, 等.基于ANSYS的真实软组织物理建模与穿刺仿真[J]. 河北大学学报(自然科学版),2015,35(2):188-192.
[16] 李鹏, 梁瑞, 申龙朵, 等. 应用Mimics和ANSYS软件建立下颌骨三维有限元模型的方法学研究 [J]. 实用口腔医学杂志,2012,28(6):709-713.
[17] 赵雪岩, 黄任含, 楼航迪, 等. 阻塞性睡眠呼吸暂停综合征的生物力学研究[J].北京大学学报(自然科学版),2009,45(5):737-742.
[18] TSOU L.The effects of muscle aging on hyoid motion during swallowing: a study using a 3D biomechanical model. University of British Columbia. 2012.
[19] 李宁, 李进让, 孙建军, 等. 健康成年人咽部吞咽功能客观参数的测量[J]. 中华耳鼻咽喉头颈外科杂志,2012,47(11):884-883.
[20] 杨晓琼,向莹媛,廖梅红,等.分期针刺治疗脑卒中后吞咽障碍的临床观察[J].江西中医药,2024,55(6):54-57.
[21] 邵天祥, 金欣, 陈晓锋, 等.舌骨喉复合体运动与咽期吞咽障碍关系的研究进展[J]. 临床医学研究与实践,2023,8(7):178-181.
[22] 刘春.不同舌骨上肌群加强训练方法在脑卒中吞咽障碍患者中的应用效果[J].中国医学创新,2024,21(30):82-86.
[23] ALVES MRM, OLIVEIRA NETO IC, ZICA GM, et al. Hyoid displacement patterns in healthy swallowing. Einstein (Sao Paulo). Einstein (Sao Paulo). 2022;20:eAO6268.
[24] 房芳芳,王孝文,郭振丽.基于VFSS检查结果观察脑卒中后咽期吞咽障碍患者吞咽咽期舌骨和喉向前、向上移动距离及速率变化[J].山东医药,2024, 64(21):57-59.
[25] 白林.针刺舌骨上肌群对脑卒中后口腔期吞咽障碍的临床疗效观察[D].哈尔滨:黑龙江中医药大学,2023.
[26] 常丽.超声引导下针刺环咽肌及舌骨上肌群对脑卒中后咽期吞咽障碍患者的临床研究[J].现代医学与健康研究电子杂志,2023,7(13):7-9.
[27] 贺凯,邢文华,李峰,等.有限元法在脊柱胸腰段骨折生物力学分析中的应用及发展方向[J].中国组织工程研究,2025,29(15):3244-3252.
[28] 卜寒梅,王旭,冯天笑,等.不同条件下旋提手法对颈椎关节突关节软骨应力的影响[J].中国中医骨伤科杂志,2024,32(3):35-39+44.
[29] 周宗昊,罗思阳,陈佳文,等.生物可吸收板与微型钛板在不同骨质下颌骨骨折固定中的有限元分析[J].中国组织工程研究,2025,29(4):818-826.
[30] PELTERET JP, REDDY BD. Computational model of soft tissues in the human upper airway. Int J Numer Method Biomed Eng. 2012;28(1):111-132.
[31] 张岳,乔海,李幸芳,等.颌骨的逆向建模及咀嚼肌牵动的有限元分析[J].科技与创新,2020(22):1-3.
[32] HANNAM AG, STAVNESS I, LLOYD JE, et al. A dynamic model of jaw and hyoid biomechanics during chewing. J Biomech. 2008;41(5):1069-1076.
[33] 叶萍, 李兆申, 许国铭, 等.反流性食管炎患者食管动力学临床研究[J].第二军医大学学报,2001,22(3):204-206.
[34] 任攀.食管的各向异性本构模型研究及模拟分析[D].成都:西南交通大学, 2021.
[35] 李振亚,王光明,徐元顺.吞咽造影录像检查定量评价脑卒中患者吞咽功能改变[J].中国介入影像与治疗学,2022,19(6):361-364.
[36] 唐安丽.脑卒中后吞咽障碍患者颏下肌群功能和舌骨运动的特征研究[D].广州:南方医科大学,2022.
[37] 王丽琦,田力,邓苗苗.超声定量评估脑卒中后吞咽障碍患者颏舌骨肌[J].中国医学影像技术,2023,39(3):425-428.
[38] 谢亚青,毛忠南,张晓凌,等.基于Drp1、Opa1蛋白与线粒体联系探讨针刺二腹肌前腹治疗卒中后吞咽障碍的作用机制[J].中医药信息,2022,39(3):85-88.
[39] 陈春花,卓春和,梁沛君.颏下超声测量舌肌厚度变化舌骨位移幅度与口咽期吞咽障碍脑性瘫痪患儿吞咽功能病情程度及肌肉纤维化的相关性分析[J].中国妇幼保健,2020,35(22):4263-4266.
[40] 陈兆聪,曹君妍,喻勇,等.鼻咽癌放疗后吞咽困难患者肌肉纤维化与舌骨位移的相关性研究[J]. 中华物理医学与康复杂志,2017,39(12):903-907.
[41] 徐雷鸣. 二腹肌后腹影像解剖研究及其意义[D]. 杭州:浙江大学,2005.
[42] 郝义勇, 徐浩铜, 向轲, 等. 利用可视人项目与CT研究左膈下腹膜外间隙空间关系[J]. 西南国防医药,2015,25(9):938-941.

引言

随着人口老龄化及疾病谱的变化，吞咽障碍患病率逐年增加，严重影响患者的生活质量，严重者可出现营养不良、吸入性肺炎等并发症[1]。研究显示一些疾病有加重吞咽障碍的风险[2-3]，其中脑卒中、颅脑损伤、帕金森等疾病与吞咽障碍密切相关[4]。目前不乏吞咽障碍相关临床试验[5-6]，然而关于其机制研究仍缺乏较好的手段，需要通过生物力学分析及客观的指标进行评价。有限元法是生物力学研究常用方法之一，适用于对三维人体结构建立的力学模型进行分析，具有可视化和无创性的优势，避免了尸源紧缺、实体研究不便的弊端[7-9]。有限元分析是将研究对象进行数学化模拟，将一个整体分割成小单元的组合体，其实质是利用有限数量的模拟量去模仿无限单元的真实结构，将各形质构造数学化进行计算[10]。有限元建模需要采集和提取原始数据，常用影像资料和DICOM数据建模法，肌骨模型的构建常用有限元软件进行处理，如Mimics、Geomagic、Ansys等。
由于人体吞咽系统的复杂性和体内数据测量的困难性，吞咽的生物力学研究受到了极大的限制，目前尚未有构建吞咽肌骨三维有限元模型的相关文献报道，吞咽运动中各结构的力学变化尚不清晰。舌骨运动过程复杂，涉及众多细小肌群的协同作用。其中，舌骨上肌群在舌骨向上向前运动中起到了非常重要的作用[11]。鉴于舌骨运动相关肌群众多，以及建模的复杂性，此次研究重点关注舌骨上肌群，拟利用舌骨的位移参数及头颈部CT构建吞咽模型。通过建立包括C2-C4椎体、下颌骨、舌骨、部分颞骨及舌骨上肌群的三维模型进行有限元分析，旨在观察舌骨运动状态下舌骨上肌群对舌骨的应力作用、异常吞咽下舌骨上肌群的应力改变及分布，为吞咽障碍的力学机制研究及临床诊疗提供思路。

中国组织工程研究杂志出版内容重点：组织构建；骨细胞；软骨细胞；细胞培养；成纤维细胞；血管内皮细胞；骨质疏松；组织工程

材料方法

1 对象和方法 Subjects and methods
1.1 设计基于有限元分析法，通过志愿者的CT图像构建三维模型，结合既往研究中健康人吞咽造影中舌骨位移参数，模拟舌骨运动的力学环境，进一步分析舌骨及舌骨上肌群的应力改变。
1.2 时间及地点实验于2022年8月至2023年1月在浙江中医药大学附属第三医院完成。
1.3 对象选取1名成年健康志愿者并采集头颈CT影像数据。健康志愿者年龄56岁，身高171.50 cm，体质量65.30 kg；要求符合以下条件：无头部创伤史，面型协调对称；无脑卒中及吞咽障碍病史；面部无创伤史，面型协调对称，牙列完整；无风湿、类风湿病史及其他系统性疾病。健康志愿者已签署知情同意书，且取得浙江中医药大学附属第三医院的伦理审核(伦理编号：ZSLL-KY-2022-027-01)。
1.4 设备及软件
1.4.1 CT扫描仪采用西门子螺旋扫描CT(SOMATOM EMOTION 16)对志愿者头颈部进行扫描。取仰卧位，头颈部两侧对称无偏斜，牙列自然对合，摄片中勿吞咽和咀嚼。层厚0.6 mm，电压120 kV，电流335 mAs，矩阵512×512，共获得断层影像387张，以DICOM数据格式导出。
1.4.2 计算机配置 ①处理器：AMD Ryzen 7 4800U；②CPU：1.80 GHz；③内存：16 GB；④系统软件：操作系统Windows 10×64位。
1.4.3 图像处理与有限元软件Mimics 20.0：具有强大的交互功能及多种医学图像处理功能，可根据各种二维扫描数据建立3D模型，并可对建立的3D模型进行距离及角度测量。该软件除了构建骨骼模型，还能进行心脏、肿瘤、血管等器官和软组织的模型重建，常用于手术规划及生物力学分析[12-13]。
Geomagic Studio 2021：一款较为常用的逆向工程软件，可以把任何产品或零件通过点云计算逆向生成精确的数字模型，在计算机上进行实物还原，并提供点、多边形及曲面模型等多种建模格式。该软件可用于病变器官及部位的医学重建，尤其在口腔科领域应用普遍[14]。
Ansys 2022 R1：该软件主要包括前处理模块、分析计算模块以及后处理模块3个部分。其中前处理模块具有强大的实体建模和网格划分功能，可进行复杂生物组织的数值模拟和仿真[15-16]。
1.5 方法
1.5.1 吞咽肌骨模型的构建下颌舌骨肌的起始位置骨性标志明确，颏舌骨肌位于其深面，二者解剖结构上十分贴近，形态相对易于观察和提取，故在研究肌肉时将下颌舌骨肌与颏舌骨肌共面处理，视为整体，即下颌舌骨肌-颏舌骨肌。余舌骨上肌群由于形态及走向复杂等原因，构建难度较大，故肌肉部分仅根据CT提取下颌舌骨肌-颏舌骨肌的三维模型。C2-C4椎体作为构建舌骨位移参照坐标系的骨性标志，部分颞骨用于测量二腹肌及茎突舌骨肌的肌力方向。故研究构建了包含下颌骨、下颌舌骨肌-颏舌骨肌、舌骨、C2-C4椎体及部分颞骨的吞咽模型。
1.5.2 Mimics构建三维模型首先，将DICOM格式的CT图像文件导入Mimics，设置图层的位置关系，软件根据文件自动生成冠状位、矢状位及横断面的断层图像，见图1。选择阈值设定“Bone (CT)”并适当调整阈值集，直到将骨骼与其他软组织明显区分且结构完整，见图2。不同组织选用不同颜色的蒙面，建立骨骼及软组织的蒙版，见图3。由于椎体、下颌骨及部分颞骨存在骨质相连的情况，使用Chop、Edit Mask、Region Growing等手动编辑功能，分别在冠状面、矢状面和水平面3个层面将不同组织分离，初步提取下颌骨、舌骨、椎体、下颌舌骨肌-颏舌骨肌及部分颞骨。经过软件计算后并初步光滑处理，得到粗糙的包含下颌骨、舌骨、椎体、下颌舌骨肌-颏舌骨肌及部分颞骨的三维模型，见图4。将重建好的三维模型分别以STL格式导出。

1.5.3 Geomagic平滑处理及构建曲面片将STL文件导入Geomagic中，设置单位为mm。此处以舌骨为例，根据模型表面情况选择多边形下松弛、降噪、磨砂、删除、填充等处理，强度设置为1/4，避免由于强度过度导致模型失真。光滑处理后，探测抽取轮廓，手动编辑轮廓线，

构建曲面片。对于问题部分，如相交，手动编辑调整，直到获得合理的曲面片。在曲面片的基础上进行格栅构造，最后进行曲面拟合，见图5。以同样的方法处理下颌骨、椎体、下颌舌骨肌-颏舌骨肌及部分颞骨，最终文件以SAT格式导出，见图6。

1.5.4 Ansys有限元分析将下颌骨、舌骨、下颌舌骨肌-颏舌骨肌和椎体分别以SAT格式导入ANSYS软件，设定单位长度(mm)、质量(kg)、力(N)、时间(s)。设定模型中各组织为各向同性、均质连续的线弹性材料，根据赵雪岩等[17]研究，定义材料属性：骨组织的弹性模量为1 137×1010 Pa，泊松比为0.30；软组织的弹性模量为1 100×104 Pa，泊松比为0.45。对不同组织划分网格，单元类型采用10节点的Solid187四面体单元。对舌骨、下颌骨、下颌舌骨肌-颏舌骨肌及部分颞骨进行网格划分，调至合适尺寸，此处以舌骨为例，单元尺寸调整为1，见图7。下颌骨、舌骨、下颌舌骨肌-颏舌骨肌及部分颞骨的单元数和节点数见表1。舌骨上肌群包含下颌舌骨肌、颏舌骨肌、茎突舌骨肌及二腹肌，肌肉收缩力则参考TSOU[18]研究中Fmax=PCSA×40 N/cm2，其中PCSA为生理横截面积，见表2。

1.5.5 舌骨最低及最高位置加载舌骨上肌群肌力将各组舌骨上肌群的肌力全部加载于舌骨的接触位置，约束舌骨中未与舌骨上肌群的接触部分，观察舌骨的应力分布及变化。首先要根据解剖骨性标志测量该肌肉与骨骼的相对位置，肌肉的解剖起止标志见表3。基于三维重建模型，根据解剖位置选取相应节点，测量每组肌肉的各方向分量，将各方向分量结合肌力，转换为X、Y、Z 轴3个方向的肌力值，舌骨最低位置时舌骨上肌群各肌肉的肌力分量见表 4。由于目前所得肌力为生理状态舌骨肌群的最大收缩力，同样将健康人的舌骨位移参数(前移10.50 mm，上抬14.40 mm)赋值给模型[19]，观察舌骨运动至最低位置时的应力分布及变化，并与最高位置的舌骨做比较。对于舌骨位移，以C2椎体前下端与C4椎体前下端连线为纵轴，穿过舌骨对称中心与纵轴的垂直线为横轴建立坐标系，进而使模型的舌骨产生预定的位移量。用同样的方法手动测得舌骨最高点时舌骨上肌群各方向肌力，见表5。

1.5.6 加载不同舌骨上肌群肌力根据之前测量的数据，将不同舌骨上肌群肌力分别加载于舌骨相应接触部位，每次加载一组肌群则约束其他未与该组肌群接触的舌骨部分，观察不同舌骨上肌群对舌骨的应力影响，舌骨上肌群肌力方向见图8。

1.5.7 改变材料属性以模拟不同程度的吞咽障碍设定下颌骨、舌骨与下颌舌骨肌-颏舌骨肌的接触关系为绑定，分别形成2个接触对。根据吞咽时舌骨运动情况，分别对下颌骨和部分颞骨施加全约束。由于给定材料属性视为肌肉的生理状态，通过改变肌肉材料的弹性模量，模拟不同程度吞咽障碍时肌肉的僵硬程度，将肌肉的弹性模量分别增加25%、50%及75%，分别对应轻度、中度及重度吞咽障碍。对舌骨加载的位移参数仍不变，观察不同吞咽障碍程度下舌骨及下颌舌骨肌-颏舌骨肌的应力变化。
1.6 主要观察指标 ①舌骨上肌群对舌骨的应力影响；②吞咽障碍严重程度对下颌舌骨肌-颏舌骨肌的应力影响。

讨论

吞咽是一个复杂的动态过程，一般可大致分为口腔期、咽期及食管期[20]。其中，发生在咽期的吞咽障碍会导致食物对气道的侵入，引起呛咳、窒息等症，并导致吸入性肺炎[21]。因而，咽期也是吞咽动作的最关键时期，当舌将食团塞进咽腔后，将软腭向上提起，同时保持舌骨向前和向上的动作，会厌返折将喉盖封闭后预防误吸，方便食管内食物顺利进入[22]。大量的研究表明，舌-喉复合体的上提和前移使得气道关闭和食管上括约肌开放，有助于气道保护，防止误吸[23-24]。吞咽反射接受多种肌肉的控制，促进舌骨以及喉向上、向前的运动，而舌骨肌群的有力收缩对于吞咽舌骨的运动十分重要。许多临床研究着重于加强舌骨上肌群的功能，从而改善吞咽功能[25-26]，但其中的力学机制并不明晰，可见对于舌骨运动及吞咽相关肌群的研究意义重大。
有限元法是生物力学研究常用分析方法，广泛应用于骨科、脊柱及口腔领域，适用于对三维人体结构建立的力学模型进行生物力学分析[27]。有限元法具有实验条件要求简单、可重复性高、结果形象直观等优势，通过简单而又相互作用的元素(单元)，用有限数量的未知量去逼近无限未知量的真实系统，能够计算模型整体及局部的应力大小和位移变化[28-29]。目前国内吞咽方面的有限元研究较少，鲜有吞咽肌骨三维模型的构建，且相关的生物力学分析甚少。
目前国内外关于吞咽建模的现状不容乐观，尚未有完整吞咽肌骨有限元模型的构建，且对于吞咽的生物力学研究多集中在口腔期及食管期，对于舌骨运动的研究则更为有限。PELTERET等[30]建立了包含舌体、舌骨及舌骨上肌群的有限元模型，然而研究的重点在于舌体并非舌骨运动。张岳等[31]构建了包含上颌骨、下颌骨及牙列的有限元模型，分析颌骨在咀嚼肌牵动下的应力分布。HANNAM等[32]仅建立了咀嚼过程中包含下颌和舌骨的模型。TSOU[18]构建了舌骨喉复合体的三维模型，发现肌肉老化会降低舌骨运动幅度，进而影响正常的吞咽过程。叶萍等[33]观察了吞咽过程中食管蠕动输送变化、吞咽反流情况和下食管括约肌的生理状态变化。任攀[34]进行了食管组织的有限元模拟，侧重于食管生理狭窄区域的力学分析。以上研究更多关注咀嚼动作或食管的力学分析，虽各有侧重，均未能反映吞咽运动中舌骨及相关肌群的应力应变。因而，对于吞咽肌骨模型的构建及舌骨运动的研究有助于阐明吞咽的生物力学机制，也可为吞咽障碍的临床诊疗提供参考。此次研究则首次基于CT影像数据，完善了吞咽咽期舌骨运动所包含的舌骨、下颌骨、舌骨上肌群等肌骨模型的构建并分析其生物力学特征，扩展了有限元在吞咽相关研究中的应用。对于生理状态下吞咽舌骨运动力学环境的研究，不仅有利于了解正常吞咽的应力变化，从而为异常吞咽模式的临床治疗提供思路。目前临床上多运用吞咽造影评估患者的吞咽功能，可客观、动态地观察舌骨运动并记录其位移等参数，被认为是吞咽障碍诊断的金标准[35]，但仅局限于正侧位的定量分析，不能够营造三维的力学环境。作者以吞咽造影中舌骨运动的轨迹为基础，选定了吞咽启动时的舌骨最低位置和舌骨最高位置两个关键点进行建模并受力分析，直观地展示了舌骨上肌群对舌骨的应力及分布情况。舌骨在最低位置时，吞咽初始启动，需要更大的肌肉收缩力使舌骨前移及上抬。当舌骨越接近最高位置，由位移导致肌肉力的方向也发生变化，舌骨运动开始出现下移的趋势，最高位置的应力值降低。结果显示舌骨在最低位置受到的应力值更高，符合吞咽的生理状态。
舌骨上肌群附着于下颌骨及颞骨，由4组肌肉组成，分别是下颌舌骨肌、颏舌骨肌、茎突舌骨肌及二腹肌，在舌骨向上向前移位中起到了非常重要的作用。唐安丽[36]发现健康人的颏下肌群(下颌舌骨肌、颏舌骨肌和二腹肌前部)肌电活动时间与舌骨运动时间、颏下肌群肌电峰值与舌骨最大运动距离均存在显著相关性，认为舌骨向上、向前运动的动力主要来源于颏下肌群收缩牵拉。王丽琦等[37]研究提示吞咽障碍患者康复治疗过程中颏舌骨肌运动幅度恢复较快。有研究提出针刺二腹肌前腹可以调控负责供能的线粒体蛋白的表达，进而提升吞咽质量[38]。以上研究通过表面肌电、超声等方式进行评估，然而由于测量的局限性，未见有不同肌群作用的比较。此次研究通过在模型上加载不同舌骨上肌群肌力，发现下颌舌骨肌对舌骨应力影响最大，其次是二腹肌、颏舌骨肌及茎突舌骨肌，这与上述颏下肌群肌电活动及舌骨运动特征也基本一致。此次研究从另一方面也提示关键肌肉的功能对于舌骨运动具有重要的作用，可为吞咽康复治疗提供基础，对揭示吞咽障碍的发病机制有一定的指导意义。对于吞咽障碍者，重点训练下颌舌骨肌、二腹肌也许是促进舌骨上抬和前移的有效手段，可进一步研发特定肌肉刺激仪或更为精准的治疗方法。陈春花等[39]研究发现肌肉纤维化水平增加是影响吞咽障碍患者舌骨位移幅度的独立危险因素。陈兆聪等[40]认为肌肉弹性模量与肌肉纤维化密切相关，可作为评价患者吞咽功能的辅助监测手段。此次研究以改变肌肉弹性模量的方式模拟不同程度吞咽障碍，发现应力峰值与吞咽障碍程度呈正相关，且应力多集中于舌骨接触的近端。在异常吞咽的模拟中，发现下颌舌骨肌-颏舌骨肌的应力集中区域主要分布在该肌肉与舌骨接触的近端，其中应力峰值位于舌骨体与舌骨大角的交界处。这提示改善舌骨上肌群应力分布、降低应力峰值可能是改善吞咽障碍的生物力学机制之一，具体需结合临床试验进一步扩大样本量，收集试验数据并分析舌骨及肌群应力变化，从而提高可信度。
此项研究还存在一定局限性：①由于肌肉建模的复杂性，研究仅提取了对舌骨运动影响较大的舌骨上肌群部分，其他细小肌群可能也会对结果产生一定影响，未来可优化建模方法，完善吞咽细小肌群的仿真建模[41-42]；②样本年龄单一，未考虑现实中结构退化的力学影响，与真实人体存在一定的误差，且未考虑性别因素，后续可扩大年龄范围、纳入女性志愿者行深入对比研究，以提供更全面的参考意见；③以肌肉材料参数的赋值模拟吞咽障碍程度，并未考虑病程的影响，未来可纳入病程因素进行对比研究，以提供更有力的数据支撑；④受技术限制，未能分析整个舌骨动态轨迹，只选择了舌骨运动两个瞬时状态进行静力分析。
综上所述，此项研究结合舌骨的位移参数，模拟了吞咽时舌骨最低及最高位置两个瞬时状态的生物力学环境，基本能够反映吞咽时舌骨运动情况下舌骨及舌骨上肌群的应力变化。其一，舌骨上肌群中，下颌舌骨肌及二腹肌对舌骨应力影响较大，提示特定肌肉靶向刺激的重要性，为基于特定肌肉刺激的吞咽障碍治疗方法提供指导。其二，模拟了不同程度吞咽障碍下应力改变情况，随着吞咽障碍程度加重，舌骨肌群的应力峰值增高，未来可进一步探索吞咽障碍的诊断与应力应变的相关性。其三，将有限元法引入吞咽肌骨建模及分析中，为吞咽障碍在医工交叉领域的深入研究打开思路。

中国组织工程研究杂志出版内容重点：组织构建；骨细胞；软骨细胞；细胞培养；成纤维细胞；血管内皮细胞；骨质疏松；组织工程

CT影像数据建立吞咽肌骨三维模型以及生物力学特征分析

Three-dimensional model of swallowing musculoskeletal system based on CT image data and biomechanical characteristics analysis

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表（结果） 3

参考文献

相关文章 15

引言

材料方法

讨论

文章快阅

延伸阅读

编辑推荐

Metrics

本文评价

[1]	许浩, 丁璐, 李潇. 中央螺丝机械磨损对莫氏锥度种植系统力学影响的三维有限元分析[J]. 中国组织工程研究, 2025, 29(在线): 1-9.
[2]	孙晓君, 王华明, 张德宏, 宋学文, 黄晋, 张辰, 裴生太. 有限元法在儿童发育性髋关节发育不良及治疗中的作用[J]. 中国组织工程研究, 2025, 29(9): 1897-1904.
[3]	李良奎, 黄永灿, 王鹏, 于滨生. 颈椎前路椎体骨化物可控前移融合对后纵韧带骨化物和内植物影响的有限元分析[J]. 中国组织工程研究, 2025, 29(9): 1761-1767.
[4]	徐彪, 路坦, 姜亚琼, 阴玉娇. 有限元分析不同程度冈上肌断裂对肩关节应力的影响[J]. 中国组织工程研究, 2025, 29(9): 1768-1774.
[5]	周金海, 李江伟, 王序全, 庄颖, 赵瑛, 杨渝勇, 王嘉嘉, 杨阳, 周仕炼. 不同骨强度下全膝置换过程中发生股骨前皮质切迹的三维有限元分析[J]. 中国组织工程研究, 2025, 29(9): 1775-1782.
[6]	陈曦, 汤涛, 陈铜兵, 李青, 张文. 不同内固定系统治疗股骨转子间骨折的力学稳定性[J]. 中国组织工程研究, 2025, 29(9): 1783-1788.
[7]	付恩洪, 杨行, 梁成, 张小刚, 张亚丽, 靳忠民. 基于OpenSim预测青少年跖屈肌无力的下肢生物力学行为[J]. 中国组织工程研究, 2025, 29(9): 1789-1795.
[8]	芦劼明, 李亚静, 杜培洁, 徐冬青. 人工和天然草地对青年女性跳跃落地时下肢生物力学表现的影响 [J]. 中国组织工程研究, 2025, 29(6): 1101-1107.
[9]	张帅, 李子春, 徐亦豪, 谢晓峰, 郭忠圣, 赵清扬. 经颅磁声电刺激强度对小鼠前额叶皮质网络可塑性的影响[J]. 中国组织工程研究, 2025, 29(6): 1108-1117.
[10]	李帅, 刘桦, 商永慧, 刘义琮, 赵启航, 刘文. 安氏Ⅱ类错牙合畸形佩戴Twin-block矫治器时上颌骨的应力分布[J]. 中国组织工程研究, 2025, 29(5): 881-887.
[11]	周宗昊, 罗思阳, 陈佳文, 陈光能, 冯红超. 生物可吸收板与微型钛板在不同骨质下颌骨骨折固定中的有限元分析[J]. 中国组织工程研究, 2025, 29(4): 818-826.
[12]	陈倚珑, 张旭, 李红. 纤维桩联合不同冠部修复方式修复穿髓型非龋性颈部缺损牙体的力学分析[J]. 中国组织工程研究, 2025, 29(4): 866-871.
[13]	陈静, 张楠, 孟庆华, 鲍春雨. 不同年龄阶段膝关节有限元计算模型的材料特性分析[J]. 中国组织工程研究, 2025, 29(34): 7369-7375.
[14]	庄妍, 王心彧, 曹怡琳, 丁元欣, 王佳琦, 于淼, 栾春阳, 丁元圣. 3D打印上颌单根扭转牙个性化正畸装置的三维有限元分析[J]. 中国组织工程研究, 2025, 29(30): 6409-6415.
[15]	艾克丽亚•艾尼瓦尔, 娜菲莎•吾普尔, 白布加甫•叶力思, 孜拉来•居来提, 古丽再努•依布拉音, 尼加提•吐尔逊. 上颌窦内提升不植骨与植骨即刻负重的动态应力分析[J]. 中国组织工程研究, 2025, 29(30): 6416-6425.