三维有限元模型力学分析可预测视乳头的形状变化

doi:10.3969/j.issn.2095-4344.2013.50.015

中国组织工程研究

• 组织构建基础实验 basic experiments in tissue construction • 上一篇下一篇

三维有限元模型力学分析可预测视乳头的形状变化

祁昕征，魏超，杨佳燕，穆晶，张昆亚，刘志成，钱秀清

首都医科大学生物医学工程学院，北京市 100069

修回日期:2013-11-09 出版日期:2013-12-10 发布日期:2013-12-10
通讯作者: 钱秀清，博士，副教授，首都医科大学生物医学工程学院，北京市 100069 qianxq@ccmu.edu.cn
作者简介:祁昕征，女，1992年生，山东省聊城市人，汉族，2013年首都医科大学毕业，主要从事生物力学方面的研究。 Qi.xinzheng@163.com
基金资助:
国家自然科学基金(11102123, 31070840)**；北京市教育委员会科技计划面上项目(KM201110025009)*；北京市属市管高校人才强教计划资助项目(PHR201008398, PHR201110506)**

Shape variation of optic nerve head by mechanical analysis using three-dimensional finite element model

Qi Xin-zheng, Wei Chao, Yang Jia-yan, Mu Jing, Zhang Kun-ya, Liu Zhi-cheng, Qian Xiu-qing

School of Biomedical Engineering, Capital Medical University, Beijing 100069, China

Revised:2013-11-09 Online:2013-12-10 Published:2013-12-10
Contact: Qian Xiu-qing, Ph.D., Associate professor, School of Biomedical Engineering, Capital Medical University, Beijing 100069, China qianxq@ccmu.edu.cne
About author:Qi Xin-zheng, School of Biomedical Engineering, Capital Medical University, Beijing 100069, China Qi.xinzheng@163.com
Supported by:
the National Natural Science Foundation of China, No. 11102123*, 31070840*; Scientific General Program of Beijing Municipal Education Commission, No. KM201110025009*; Funding Program for Academic Human Resources Development in Institutions of Higher Learning Under the Jurisdiction of Beijing Municipality of China, No. PHR201008398*, PHR201110506*

摘要/Abstract

摘要：

背景：视乳头形状变化的准确测量及预测是青光眼的早期诊断及病程预测的关键，因此高眼压下视乳头形状的测量与分析具有重要意义。
目的：建立视乳头三维有限元模型，并分析急性高眼压下视乳头形状及视网膜厚度的变化。
方法：利用前房灌注的方法建立急性高眼压动物模型，利用光相干断层扫描仪获得猫在正常眼压及5 320， 7 980，10 640，13 300，15 960 Pa眼压时的断层图片，测量视网膜特征点处厚度的变化；基于正常眼压猫眼视乳头的断层数据，利用MIMICS软件获得视网膜及脉络膜的三维结构，并进行组装获得视乳头三维模型；利用有限元软件ABAQUS分析高眼压状态下视乳头的形状变化及视网膜厚度变化，并与实验测试结果进行比较。
结果与结论：随着眼压的升高，发现视网膜厚度逐渐变薄，视乳头深度及视乳头宽度均逐渐增加，视乳头宽度与深度的比值也逐渐增加。提示急性高眼压导致视网膜变薄，视乳头宽度增大，视乳头深度增大。结果证明，利用光学相干断层扫描仪进行活体视乳头建模和分析是可行的，力学分析可预测视乳头的形状变化，对进一步确定青光眼的病理进程有一定的指导意义。

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

全文链接：

关键词: 组织构建, 组织构建基础实验, 急性高眼压, 三维重建, 有限元分析, 视乳头, 视网膜, 国家自然科学基金

Abstract:

BACKGROUND: The measure and prediction of the shape of optic nerve head are the key issues for early diagnosis and duration prediction of glaucoma. Therefore, it is significant to research the shape variation of optic nerve head under high intraocular pressure.
OBJECTIVE: To establish three-dimensional finite element model of optic nerve head and to analyze the deformation of optic nerve head and the change in retina thickness after acute high intraocular pressure.
METHODS: Animal model of acute ocular hypertension was established by methods of anterior chamber perfusion. The tomographic images of the optic nerve head of a cat were obtained at the normal intraocular pressure and 5 320 Pa, 7 980 Pa, 10 640 Pa, 13 300 Pa, 15 960 Pa intraocular pressures using optical coherence tomography. Then, we measured the variation of retinal thickness at typical location. Basing on the tomographic images of optic nerve head of a cat at normal intraocular pressure, we obtained three-dimensional structure of retina and choroid using software MIMICS. Then three-dimensional model of optic nerve head was established by assembling the retina and choroid. The deformation of optic nerve head and the change in retinal thickness under high intraocular pressure were observed using software ABAQUS. The effectiveness of the model was verified by the experimental result.
RESULTS AND CONCLUSION: With the increase of intraocular pressure, the retinal thickness is thinner, the optic nerve head depth, width and the ratio of the depth to width are gradually increased. It suggests that the acute high intraocular pressure causes retinal thinning, optic nerve head widening and deepening. It is feasible to establish optic nerve head modeling in vivo by using optical coherence tomography, mechanical analysis can be applied to predict the shape variation of optic nerve head, which is significant to further deduce the pathological process of glaucoma.

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

全文链接：

Key words: high intraocular pressure, glaucoma, imaging, three-dimensional, finite element analysis, retina

中图分类号:

R318

祁昕征，魏超，杨佳燕，穆晶，张昆亚，刘志成，钱秀清. 三维有限元模型力学分析可预测视乳头的形状变化[J]. 中国组织工程研究, doi: 10.3969/j.issn.2095-4344.2013.50.015.

Qi Xin-zheng, Wei Chao, Yang Jia-yan, Mu Jing, Zhang Kun-ya, Liu Zhi-cheng, Qian Xiu-qing. Shape variation of optic nerve head by mechanical analysis using three-dimensional finite element model[J]. Chinese Journal of Tissue Engineering Research, doi: 10.3969/j.issn.2095-4344.2013.50.015.

图/表（结果） 1

参考文献

[1] Quigley HA, Addicks EM, Green WR, et al. Optic nerve damage in human glaucoma. II: The site of injury and susceptibility to damage. Arch Ophthalmol.1981;99(4): 635-649.

[2] 季宝玲, 游逸安, 方爱武, 等. OCT测量视网膜神经纤维层和黄斑厚度在青光眼早期诊断中的意义[J]. 眼视光学杂志,2008, 10(1): 54-58.

[3] 胡赛静, 邱翎. 正常眼与原发性开角型青光眼的OCT视盘参数分析的对比研究[J]. 临床眼科杂志, 2007, 15(1):26-28.

[4] 戴惟葭, 边俊杰, 杨惠青, 等. 急性闭角型青光眼视网膜神经纤维层改变的一年动态观察[J]. 眼科, 2010, 19(5): 331-335.

[5] Sigal IA, Flanagan JG, Tertinegg I, et al. Modeling individual-specific human optic nerve head biomechanics. Part II: influence of material properties. Biomech Model Mechanobiol. 2009;8:99-109.

[6] Yang H, Downs JC, Sigal IA, et al. Deformation of the normal monkey optic nerve head connective tissue after acute IOP elevation within 3-D histomorphometric reconstructions. Invest Ophthalmol Vis Sci. 2009;50:5785-5799.

[7] 杨翎, 袁志兰, 占慧琴. FD-OCT测量视盘周围视网膜神经纤维层和视神经乳头参数在青光眼早期诊断中的意义[J]. 南京医科大学学报:自然科学版, 2010,30(7): 1026-1030.

[8] 杨新光, 刘钊, 于敬妮, 等. 原发性青光眼早期视盘参数与视野缺损的关系[J]. 临床眼科杂志, 2009, 17(4):295-299.

[9] 徐亮, 夏翠然, 杨桦, 等. 正常人不同类型视乳头及早期青光眼患者视乳头形态学研究[J]. 中华眼科杂志, 2002, 38(6): 325-328.

[10] Quigley HA, Addicks EM. Quantitative studies of retinal nerve fiber layer defects. Arch Ophthalmol. 1982;100: 807-814.

[11] Quigley HA, Addicks EM. Regional differences in the stucture of the laminal cribeosa and their relation to glaucomatous optic nerve damage. Arch Ophthalmol. 1981;99: 137-143.

[12] He DQ, Ren ZQ. A biomathematical model for pressure- dependent lamina cribrosa behavior. J Biomech. 1999; 32: 579-584.

[13] Belleza AJ, Hart RT, Burgoyne CF. The optic nerve head as a biomechanical structure: initial finite element modeling. Invest OphthalmolVis Sci. 2000;41(10): 2991-3000.

[14] Yang H, Downs JC, Bellezza A, et al. 3-D histomorphometry of the normal and early glaucomatous monkey optic nerve head: prelaminar neural tissues and cupping. Invest Ophthalmol Vis Sci. 2007;48:5068-5084.

[15] Sigal IA, Flinagan JG, Ethier CR. Finite element modeling of optic nerve head. Invest Ophthalmol Vis Sci. 2004;45(12): 4378- 4387.

[16] Sigal IA, Flinagan JG, Tertinegg I, et al. Reconstruction of human nerve heads for finite element modeling. Tech Health Care. 2005;13: 313 -329.

[17] Sigal IA, Flanagan JG, Tertinegg I, et al. 3D morphometry of the human optic nerve head. Exp Eye Res. 2010;90:70-80.

[18] Roberts MD, Grau V, Grimm J, et al. Remodeling of the connective tissue microarchitecture of the lamina cribrosa in early experimental glaucoma. Invest Ophthalmol Vis Sci. 2009; 50(2): 681-690.

[19] Roberts MD, Liang Y, Sigal IA, et al. Correlation between Local Stress and Strain and Lamina Cribrosa Connective Tissue Volume Fraction in Normal Monkey Eyes. Invest Ophthalmol Vis Sci. 2010;51(1):295-307.

[20] Sigal IA. An applet to estimate the IOP-induced stress and strain within the optic nerve head. Invest Ophthalmol Vis Sci. 2011;52(8): 5497-5506.

[21] Sigal IA, Yang H, Roberts MD, et al. IOP-induced lamina cribrosa displacement and scleral canal expansion: an analysis of factor interactions using parameterized eye-specific models. Invest Ophthalmol Vis Sci. 2011; 52(3):1896-1907.

[22] Norman RE, Flanagan JG, Sigal IA, et al. Finite element modeling of the human sclera: Influence on optic nerve head biomechanics and connections with glaucoma. Exp Eye Res. 2011;93: 4-12.

[23] Chen TC, Cense B, Miller JW, et al. Histologic correlation of in vivo optical coherence tomography images of the human retina. Am J Ophthalmol. 2006;141:1165-1168.

[24] Toth CA, Birngruber R, Boppart SA, et al. Argon laser retinal lesions evaluated in vivo by optical coherence tomography. Am J Ophthalmol. 1997;123:188-198.

[25] Toth CA, Narayan DG, Boppart SA, et al. A comparison of retinal morphology viewed by optical coherence tomography and by light microscopy. Arch Ophthalmol. 1997;115: 1425-1428.

[26] Strouthidis NG, Fortune B, Yang H, et al. longitudinal change detected by spectral domain optical coherence tomography in the optic nerve head and peripapillary retina in experimental glaucoma. Invest Ophthalmol Vis Sci. 2011;52(3): 1206-1219.

[27] Yang H, Qi J, Hardin C, et al. Spectral-domain optical coherence tomography enhanced depth imaging of the normal and glaucomatous nonhuman primate optic nerve head. Invest Ophthalmol Vis Sci. 2012;53: 394-405.

[28] Reis AS, Sharpe GP, Yang H, et al. Optic disc margin anatomy in patients with glaucoma and normal controls with spectral domain optical coherence tomography. Ophthalmology. 2012;119:738-747.

[29] Inoue R, Hangai M, Kotera Y, et al. Three-dimensional high-speed optical coherence tomography imaging of lamina cribrosa in glaucoma. Ophthalmology. 2009;116: 214-222.

[30] Ruggeri M, Wehbe H, Jiao S, et al. In vivo three-dimensional high-resolution imaging of rodent retina with spectral-domain optical coherence tomography. Invest Ophthalmol Vis Sci. 2007;48:1808-1814.

[31] Abbott CJ, McBrien NA, Grunert U, et al. Relationship of the optical coherence tomography signal to underlying retinal histology in the tree shrew (Tupaia belangeri). Invest Ophthalmol Vis Sci. 2009;50:414-423

[32] Qiu JF, Qian XQ, Cui QQ, et al. Three-Dimensional Reconstruction and Finite Element Modeling Analysis of the Rabbit Optic Nerve Head in Acute High Intraocular Pressure. Jap J Appl Phys. 2012;51(6): 067001.

[33] Chen K, Rowley AP, Weiland JD. Elastic properties of porcine ocular posterior soft tissues. J Biomed Mater Res A. 2010; 93(2):634-645.

[34] Del Priore LV. Stiffness of retinal and choroidal tissue: a surface wrinkling analysis of epiretinal membranes and choroidal folds. Am J Ophthalmol. 2006;142(3):435-440.

[35] 钱勇, 刘菊, 段文静, 等. 频域OCT检测青光眼视乳头及黄斑区容积改变的初步研究[J]. 临床眼科杂志, 2011,19(3): 200-203.

[36] Maruko I, Iida T, Sugano Y, et al. Subfoveal retinal and choroidal thickness after verteporfin photodynamic therapy for polypoidal choroidal vasculopathy. Am J Ophthalmol. 2011; 151(4):594-603.

[37] Park HY, Jeon SH, Park CK. Enhanced depth imaging detects lamina cribrosa thickness differences in normal tension glaucoma and primary open-angle glaucoma. Ophthalmology. 2012;119: 10-20.

[38] Park SC, De Moraes CG, Teng CC, et al. Enhanced depth imaging optical coherence tomography of deep optic nerve complex structures in glaucoma. Ophthalmology. 2012;119: 3-9.

[39] 张昆亚, 殷莹, 张伟, 等. 利用连通器原理设计体内加压及眼球体积监测装置观测兔眼球蠕变特性[J]. 中国组织工程研究与临床康复, 2009, 43(13): 8499- 8502.

[40] 王兰, 梁远波, 王宁利, 等.原发性开角型和闭角型青光眼眼压降低前后视盘改变和筛板顺应性的对比研究[J]. 眼科, 2009, 18(4): 264-269.

引言

青光眼是全球第一大不可修复致盲眼病，临床上以视神经乳头凹陷进行性扩大及视野缺损进行性加重为特点。虽然目前标准化视野检查仍被认为是评价青光眼各种诊断指标的“金标准”，但众多研究证实，当检查出视野缺损时，视网膜神经节细胞丢失已达30%，在青光眼视野发生改变以前已有视乳头和视网膜神经纤维层的改变^[1-4]；同时多数研究认为青光眼中视网膜创伤的主要原因是高眼压下的视乳头变形所导致^[5-6]；因此对视乳头形状变化的准确测量及预测是青光眼早期诊断及病程预测的关键^[7-9]。

利用力学分析可实现对视乳头形状变化的预测。自从Quigley等^[1,10-11]研究发现青光眼视神经损害的原发部位在视乳头周围组织，人们就开始用力学手段对视乳头的变形进行分析，但由于当时的测试手段和计算技术发展水平有限，因此仅仅限于视乳头几何形态的测量或利用解析法对模型进行简化分析^[12]。Belleza等^[13]第一次建立了三维几何模型，利用有限元方法分析某一水平眼压及眼压有关的压力对巩膜管的大小和形状及视乳头周围巩膜厚度的影响。Deves眼科研究所的Yang 等^[6,14]重建了慢性青光眼的猴眼视神经乳头和筛板模型，并进行了有限元计算；猴眼模型采用激光电凝小梁网，造成持续数周的慢性青光眼；并对正常猴眼和青光眼猴眼模型做离体的石蜡组织切片，对图像配准后，定量建模。Sigal等^[15-17]把人尸眼经过一系列组织学处理后，对典型切片进行照相，建立了有限元模型，此模型在几何学、变形上与真实眼有很好的一致性。Roberts等^[18-19]则利用冷冻切片及图像处理技术对三只猴眼的筛板及其连接组织进行了三维重建，得到了很漂亮的筛板的三维结构，并进行了分析。Sigal^[20]建立了参数化模型，基于有限元分析结果编制了软件，可方便地分析视乳头各部分几何尺寸及材料特性变化时筛板的应力及应变变化；同时，他们基于猴视乳头的几何尺寸，建立了多参数有限元模型，分析了各参数的敏感性及参数间的交叉作用^[21]。他们最近的研究工作则将巩膜三维模型与视乳头模型组装分析^[22]，得到了更合理的结果。

目前视乳头的三维重建与有限元分析多是基于解剖数据或组织切片，组织切片染色的不均匀使切片的选取也具有主观性；材料特性的不确定性等使所建模型有丢失重要几何信息的危险^[16]。因此，依据在体的视乳头结构进行力学分析会得到更加合理的结果。20世纪末出现的光学相干成像技术，给眼睛后节的观察，提供了新的途径。早期的时域光学相干断层扫描仪利用可见光光源，可观察到眼底结构，获得较大范围的、能够同时捕捉黄斑和视神经乳头的眼底像。因此，临床实践中，光学相干断层扫描仪被广泛应用于眼底视神经纤维层的观察，用于青光眼诊断。 Chen^[23]和Toth等^[24-25]使用光学相干断层扫描仪获得了灵长类动物和人眼的视网膜断层影像，并进行了形态学的相关测量，但是时域光学相干断层扫描仪的光源穿透能力差，无法获取组织深部信息，采集速度也比较慢。

随着光学相干断层扫描技术的发展，应用频域光学相干断层扫描仪可得到清晰的眼底结构图像^[26-28]，分析正常人和青光眼患者视乳头形状的改变，研究发现青光眼患者的视乳头形状与正常人明显不同。最重要的是，频域光学相干断层扫描仪可以帮助我们看到青光眼试验研究中的实时的视乳头图像。对于高速频域光学相干断层扫描仪的连续平行扫描功能来说，序列图像能以超高的成像速度和轴向分辨率来显示视乳头的三维形态，并且可以得到活体视乳头三维模型用于有限元分析，这种模型能够克服离体切片建模的缺点^[29-31]。

课题组基于频域光学相干断层扫描仪的体内实验结果，得到了兔视乳头的三维结构，并进行了有限元分析，但由于图像质量限制，没有区分视网膜与脉络膜^[32]，根据前人的研究结果，视网膜与脉络膜的力学特性相差很远^[33-34]，应当加以区分。最新的研究结果发现，利用加强脉络膜扫描模式的光学相干断层扫描仪则可获得脉络膜的清晰图像^[35-38]，因此文章拟基于光学相干断层扫描仪测量活体猫眼视网膜、脉络膜的断层数据，利用三维重建技术获得视乳头的三维结构，最后利用有限元方法分析高眼压状态下视乳头变形及视网膜厚度变化，并与实验结果进行比较。

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

全文链接：

材料方法

1 视乳头断层图像的获取 Obtaining tomographic images of optic nerve head

1.1 实验动物选用健康、无眼病家猫，年龄一至两岁，体质量1.5-2.5 kg，由首都医科大学实验动物部提供，许可证号：SYXK(京)2009-0020。饲养条件为室温20-25 ℃，光照12 h，自由饮食。实验过程中对动物的处理方法符合动物伦理学要求。

1.2 实验设备光学相干断层扫描仪(3D OCT-1000 MarkⅡ，拓扑康公司，日本)，轴向分辨率最低可以达到10 µm，纵向分辨率已达到5 µm，扫描模式有三维扫描、辐射状扫描、线扫等，压力传感器。

1.3 实验方法见图1。

测量实验动物正常眼压，利用光学相干断层扫描仪采取以视乳头为中心、6 mm×6 mm×2.3 mm的三维扫描模式进行扫描，获得视乳头断层图像，如图1所示。

然后利用本实验室自行搭建的眼球加压装置^[39]，将眼压分别保持5 320，7 980，10 640，13 300， 15 960 Pa，利用光学相干断层扫描仪进行扫描，获取不同眼压时视乳头断层图像，中间凹陷部分为筛板区域。

2 视乳头有限元分析 Finite element analysis of optic nerve head

2.1 视乳头三维重建获取成年家猫正常眼压时光学相干断层扫描仪图像作为建模对象，由于乳头区上层为视网膜，视网膜下面为脉络膜，研究可知，视网膜与脉络膜的力学特性相差很远^[23-24]，筛板与视网膜力学特性相近，为了真实地模拟视乳头的变形，作者将区分视网膜与脉络膜，将获得的断层图像在Mimics软件中分别建立脉络膜和视网膜(包含筛板)的三维模型，然后在Solidworks中组装视乳头三维模型，如图2所示，导入大型有限元软件Abaqus进行分析。

2.2 有限元分析建立有限元模型，划分单元271 012个，节点数58 690。由于视乳头结构和材料性质的复杂性，考虑到软组织为不可压缩材料，采用非线性Ogden模型，Ogden模型广泛应用于橡胶的力学特性描述，其应变能函数为：

式中：λ(_)_i=J λ_i，J为变形前后体积之比。α_i，μ_i和N为材料参数。考虑到软组织的不可压缩性，不考虑上式中第2项。选取Ogden模型的主要原因是这种模型已广泛存在于各有限元软件，如ANSYS、ABAQUS、ADINA等，方便进行数值模拟。

材料参数取为一阶Ogden模型，其中脉络膜 μ₁= 90 MPa，α=6 600；视网膜μ₁=0.145 MPa，α=0.01。为了模拟巩膜对视乳头的支撑作用，在脉络膜、下表面加弹簧约束，其刚度值为6 MPa/节点数；脉络膜与视网膜之间设为绑定；视网膜上表面加均布载荷模拟高眼压的作用。为了验证模型的可靠性，将实验结果与有限元分析结果进行比较。由于Abaqus只能直接输出节点位置的变化，因此选择视乳头3个典型位置的变形进行比较，3个位置如图3所示，结果比较见图4-6。由图示可知，数值模拟结果与实验结果比较接近，认为模型可信。

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

全文链接：

三维有限元模型力学分析可预测视乳头的形状变化

Shape variation of optic nerve head by mechanical analysis using three-dimensional finite element model

PDF

可视化

被引次数

摘要/Abstract

引用本文

使用本文

图/表（结果） 1

参考文献

相关文章 15

引言

材料方法

讨论

文章快阅

延伸阅读

编辑推荐

Metrics

本文评价

[1]	刘杰, 范若寻, 高甲子, 曾晟, 刘军. 有限元法分析车辆低频振动时不同退变程度人体腰椎的力学响应[J]. 中国组织工程研究, 2020, 24(9): 1371-1377.
[2]	王祎娅, 张涵, 兰海. 椎体后凸成形骨水泥注射治疗的有限元模型数字化评价[J]. 中国组织工程研究, 2020, 24(9): 1378-1383.
[3]	刘炎, 葛鸿庆, 管华, 陈文治. 内侧柱缺失型肱骨近端骨折不同固定方式的有限元分析[J]. 中国组织工程研究, 2020, 24(9): 1384-1389.
[4]	赵春涛, 卿明松, 余浪波, 彭笳宸. 运动学对线和机械力学对线指导全膝关节置换效果的Meta分析[J]. 中国组织工程研究, 2020, 24(9): 1435-1442.
[5]	涂鹏程, 郭杨, 马勇, 潘娅岚, 郑苏阳, 王礼宁, 吴承杰, 过俊杰. 威灵仙提取物可促进体外牵张应力环境下软骨细胞表型的维持[J]. 中国组织工程研究, 2020, 24(8): 1182-1187.
[6]	张聪, 赵岩, 杜小宇, 杜欣瑞, 庞廷娟, 付怡凝, 张浩, 张步洲, 李筱贺, 王利东. 青少年特发性脊柱侧凸腰主弯患者腰椎-骨盆的生物力学分析[J]. 中国组织工程研究, 2020, 24(8): 1155-1161.
[7]	许国峰, 李学斌, 唐一钒, 赵寅, 周盛源, 陈雄生, 贾连顺. 人黄韧带细胞骨化发生过程中的自噬[J]. 中国组织工程研究, 2020, 24(8): 1174-1181.
[8]	周玉, 龙小安, 李宁, 王纯. 有限元法分析髌腱炎状态时的生物力学变化[J]. 中国组织工程研究, 2020, 24(8): 1280-1286.
[9]	林明, 潘金勇, 张惠荣. 敲除NIPBL基因可下调小鼠骨髓间充质干细胞增殖及成骨分化能力[J]. 中国组织工程研究, 2020, 24(7): 1002-1008.
[10]	岑妍慧, 夏猛, 贾微, 罗伟生, 林江, 陈松林, 陈炜, 刘鹏, 黎明星, 李景云, 李曼莉, 艾丁丁, 蒋云霞. 黄芩素能下调诱捕受体3表达抑制肝癌干细胞的生物学行为[J]. 中国组织工程研究, 2020, 24(7): 1023-1029.
[11]	刘群, 孙东东, 高丽兰, 何志江, 孙明林. 经皮椎体后凸成形后再发骨折相关因素的Meta分析[J]. 中国组织工程研究, 2020, 24(6): 976-984.
[12]	王宏, 伍权, 汤耿, 张巩. 个性化膝关节抗生素骨水泥占位器的设计及临床应用[J]. 中国组织工程研究, 2020, 24(6): 821-826.
[13]	王晓辉, 王坤, 胡志勇, 田红亮. 假肢接受腔设计及界面应力的有限元分析[J]. 中国组织工程研究, 2020, 24(6): 862-868.
[14]	和雨洁, 王海燕, 李志军, 李筱贺, 蔡永强, 戴丽娜, 许阳阳, 王一丹, 徐雪彬. 学龄期儿童胸椎经“椎弓根-肋骨”单元内固定的数字化测量[J]. 中国组织工程研究, 2020, 24(6): 869-876.
[15]	严舒, 卢岩, 欧阳昭连. 中国组织工程领域2013至2018年国家自然科学基金资助项目分析[J]. 中国组织工程研究, 2020, 24(5): 731-735.