楊永青，劉曦雨，顏澤棟，譚芷芬，羅二平，景 達

（空軍軍醫(yī)大學(xué)軍事生物醫(yī)學(xué)工程學(xué)系，西安 710032）

0 引言

骨質(zhì)疏松癥是一種全身性骨骼疾病，多見于絕經(jīng)后婦女和老年男性，其特征是骨量低和骨組織微結(jié)構(gòu)惡化，導(dǎo)致骨脆性增加和易發(fā)生骨折[1-2]。隨著人口老齡化的加劇，老年性骨質(zhì)疏松癥越來越普遍，由此帶來的經(jīng)濟和社會負擔(dān)也逐年加重[3]。骨骼在日?；顒又幸恢背惺苤鞣N形式的應(yīng)力，Wolff定律表明骨骼具有功能適應(yīng)性，能夠根據(jù)應(yīng)力環(huán)境的變化改變自身結(jié)構(gòu)和功能，進而影響自身代謝和健康狀態(tài)[4]。Frost的力學(xué)穩(wěn)定性理論進一步定義了骨塑建和骨重塑的2個基本過程，并將應(yīng)變閾值范圍劃分為廢用區(qū)（小于 200 με）、適應(yīng)區(qū)（200～1000 με）、過載荷區(qū)（1000～25000 με）和骨折區(qū)（大于 25000 με），骨骼在不同閾值范圍內(nèi)會發(fā)生相應(yīng)的變化[5-6]，因此對應(yīng)力集中區(qū)組織應(yīng)變分布的了解有助于理解應(yīng)力集中區(qū)失效的機制[7]。相關(guān)研究表明，隨著骨量下降，老年骨骼力學(xué)響應(yīng)衰退，承受負載能力下降，骨骼易發(fā)生變形和骨折[8]。目前缺少對于老年骨骼和青年骨骼力學(xué)響應(yīng)的系統(tǒng)比較和對其生物學(xué)變化的研究，因此，研究不同年齡骨骼受載的力學(xué)響應(yīng)的差異對于揭示老年骨骼力學(xué)響應(yīng)衰退發(fā)生機制、闡明骨質(zhì)疏松機制、提高骨質(zhì)量具有重要意義，有助于對老年骨質(zhì)疏松性骨折起始、預(yù)防和治療機制的更好理解。

課題組前期已開發(fā)出大鼠尺骨在體軸向周期性壓縮加載裝置對骨骼受載的力學(xué)響應(yīng)進行研究[9]。由于不同骨骼形狀大小不同，維持恒定的應(yīng)變變量對于研究不同骨骼的力學(xué)響應(yīng)至關(guān)重要，因此首先應(yīng)確定老年和青年大鼠尺骨表面相同位置達到相同的壓縮應(yīng)變需要的載荷大小。目前對于骨骼受載的應(yīng)變分布情況的研究主要是通過有限元分析和應(yīng)變片測量進行[10]。有限元分析是一種計算機數(shù)值分析程序，可根據(jù)需要建立生物模型并進行相應(yīng)實驗條件下的仿真，模擬不同實驗條件下任意部位的變形（如壓縮、拉伸、彎曲和扭轉(zhuǎn)等），進而分析其應(yīng)力應(yīng)變分布[11]。應(yīng)變片測量多采用電阻式應(yīng)變測量方法，即在被測骨骼表面粘貼應(yīng)變片，應(yīng)變片通過連接信號調(diào)理器接入采集電路中，在施加載荷后可測出應(yīng)變片粘貼位置骨骼表面相應(yīng)的應(yīng)變數(shù)據(jù)[12]。

本研究旨在通過基于Micro-CT三維重建的有限元分析和應(yīng)變片測量結(jié)合的方式，分析老年和青年大鼠的尺骨表面應(yīng)變分布情況，獲得后續(xù)實驗研究所需的尺骨表面最大壓縮應(yīng)變?yōu)?500 με時需施加的軸向壓縮載荷大小。本研究從不同層面、不同尺度對大鼠不同年齡骨骼力學(xué)響應(yīng)進行研究，以為下一步對老年性骨質(zhì)疏松的研究提供理論依據(jù)。

1 材料與設(shè)備

實驗動物：7只青年（3月齡，300～350 g）和 7只老年（22月齡，700～800 g）的 SPF級雄性 SD 大鼠[許可證號：SCXK（川）2020-030]，購自成都達碩生物科技有限公司。

Micro-CT：采用西安交通大學(xué)機械制造系統(tǒng)工程實驗室的微米X射線三維成像系統(tǒng)。

軟件：Mimics Research 19．0軟件（Materialise公司），Geomagic Studio 2014軟件（Geomagic公司），SolidWorks 2019軟件（SolidWorks公司），Abaqus 6．14-4軟件（達索SIMULIA公司）。

2 實驗方法

過量麻醉處死14只大鼠后，取下右前肢（青年大鼠7只、老年大鼠7只）置于生理鹽水中備用。各隨機選取1只青年大鼠和老年大鼠的右前肢，將組織剝離干凈后行Micro-CT掃描（如圖1所示），另外各6只用作應(yīng)變片測量。

圖1 雄性老年（上）和青年（下）SD大鼠的尺骨和橈骨

2.1 三維重建及有限元分析

對大鼠前肢行Micro-CT掃描，掃描基本參數(shù)設(shè)置為：管電壓為80 kV、管電流為 80 μA，曝光時間為 2．96 s，總旋轉(zhuǎn)角為 210°，旋轉(zhuǎn)角增量為 0．4°。沿長軸方向施加連續(xù)斷層掃描，獲得960張381×439像素的DICOM格式圖像數(shù)據(jù)。

通過對Micro-CT掃描所得的數(shù)據(jù)重建獲得三維有限元模型，并對此模型進行有限元分析計算，仿真分析流程圖如圖2所示。

圖2 三維重建及有限元分析流程圖

2.1.1 三維重建

將960張DICOM數(shù)據(jù)文件導(dǎo)入到醫(yī)學(xué)影像分割處理軟件Mimics Research 19．0中，經(jīng)過閾值分割（閾值范圍為26866～65535 HU）和分離蒙版功能除去多余組織和分離尺骨、橈骨，使用區(qū)域增長、蒙版編輯和空洞注水工具消除相應(yīng)蒙版中骨骼模型內(nèi)部空洞，最后進行3D模型計算，獲得三維重建的尺骨和橈骨曲面模型，在光滑處理后，分別以STL格式導(dǎo)出。將STL格式文件導(dǎo)入到逆向工程軟件Geomagic Studio 2014中進行曲面修補，通過偏移工具對尺骨和橈骨處理后獲得髓腔和松質(zhì)骨的結(jié)構(gòu)，經(jīng)過網(wǎng)格醫(yī)生、松弛、刪除釘狀物和快速光順等命令處理后，消除骨骼模型表面分離出的和不光滑的部分，重劃網(wǎng)格使網(wǎng)格均勻，經(jīng)網(wǎng)格醫(yī)生檢查無誤后對模型曲面化處理，依此步驟分別獲得尺骨、橈骨和髓腔、松質(zhì)骨結(jié)構(gòu)，完成后以STP格式導(dǎo)出。將STP格式文件導(dǎo)入到三維CAD軟件SolidWorks 2019中進行尺骨和橈骨的裝配，并通過布爾運算完成皮質(zhì)骨和髓腔及松質(zhì)骨結(jié)構(gòu)的裝配，生成實體后保存成X_T格式以進行后續(xù)的有限元分析。

2.1.2 有限元分析

使用有限元軟件Abaqus 6.14-4對實體模型進行有限元分析計算。首先將X_T文件導(dǎo)入到Abaqus 6.14-4中進行前處理操作，對模型分別設(shè)置材料屬性，青年大鼠尺骨和橈骨的皮質(zhì)骨彈性模量分別設(shè)置為17和15 GPa，老年大鼠尺骨和橈骨的皮質(zhì)骨彈性模量分別設(shè)置為15和13 GPa，松質(zhì)骨彈性模量均設(shè)置為300 MPa，泊松比均設(shè)置為0.3[13-15]，設(shè)置完成后建立模型裝配體。然后進行求解操作，依次定義靜力學(xué)分析步并設(shè)置相互作用，設(shè)置載荷及邊界條件，在尺骨鷹嘴端施加沿尺骨骨干方向的軸向壓力，另一端設(shè)置為6個自由度完全固定的邊界條件。設(shè)置完成后劃分四面體體網(wǎng)格，老年大鼠裝配體模型的體網(wǎng)格數(shù)和節(jié)點數(shù)分別為115235和190074，青年大鼠裝配體模型的體網(wǎng)格數(shù)和節(jié)點數(shù)分別為144923和248110。

一切設(shè)置完成后提交運算，并對求解結(jié)果進行后處理操作，查看應(yīng)變云圖并對應(yīng)變數(shù)據(jù)進行處理分析。課題組前期研究已證實在軸向壓縮載荷下峰值應(yīng)變發(fā)生在大鼠尺骨最大彎曲區(qū)域即骨干中部遠端區(qū)域，因此在尺骨遠端1/3處進行應(yīng)變評估[16]。

2.2 應(yīng)變片測量

使用課題組前期研制的小動物長骨軸向壓縮應(yīng)力載荷加載系統(tǒng)的應(yīng)變測量模塊分別對老年和青年大鼠的尺骨表面（各6個）進行應(yīng)變測量[17]。將取下的大鼠右前肢內(nèi)表面軟組織充分剝離后露出尺骨內(nèi)表面，使用M-Bond 200黏接劑將應(yīng)變片粘貼到尺骨內(nèi)表面遠端1/3處[18]，如圖3（a）所示。將大鼠前肢固定后，以0.2 mm/s的速度施加軸向壓縮載荷，如圖3（b）所示，應(yīng)變片電阻變化信息通過信號調(diào)理器處理后被采集為應(yīng)力應(yīng)變數(shù)據(jù)。

圖3 應(yīng)變片測量示意圖

3 結(jié)果

3.1 有限元分析結(jié)果

老年和青年大鼠模型的應(yīng)變云圖如圖4（a）所示，與課題組前期研究結(jié)果一致，應(yīng)變最大處均位于尺骨遠端骨干[18]。由于老年大鼠與青年大鼠的尺骨形態(tài)大小不一致，相同位置及壓縮應(yīng)變發(fā)生需要的載荷不相同，根據(jù)研究需要，僅對最大壓縮應(yīng)變發(fā)生位置（選取尺骨遠端1/3處）進行分析。應(yīng)變云圖結(jié)果顯示，老年和青年大鼠尺骨表面應(yīng)變在軸向分布情況上相同，而青年大鼠尺骨表面最大壓縮應(yīng)變發(fā)生在內(nèi)表面尺側(cè)，老年大鼠尺骨表面最大壓縮應(yīng)變發(fā)生在內(nèi)表面橈側(cè)。在尺骨內(nèi)表面達到相同的最大壓縮應(yīng)變老年大鼠模型需要的載荷比青年大鼠大，為使尺骨中下段相同位置內(nèi)表面最大壓縮應(yīng)變達到3500 με，應(yīng)分別對老年和青年大鼠施加約20和38 N 的軸向壓縮載荷[如圖4（b）所示]。

圖4 有限元分析結(jié)果

3.2 應(yīng)變片測量結(jié)果

與有限元分析結(jié)果一致，應(yīng)變片測量結(jié)果顯示，達到相同的最大壓縮應(yīng)變老年大鼠需要的載荷比青年大鼠大，為使尺骨遠端1/3處表面應(yīng)變達到-3500 με，應(yīng)分別對老年和青年大鼠施加約19和35 N的軸向壓縮載荷。軸向壓縮載荷分別為10、20、30和40 N時，老年和青年大鼠尺骨遠端1/3處應(yīng)變片測得的應(yīng)變大小和有限元分析獲得數(shù)據(jù)對比如圖5所示，應(yīng)變片測得數(shù)據(jù)與有限元分析數(shù)據(jù)偏差均小于10%。

圖5 不同加載強度下應(yīng)變片測量和有限元分析數(shù)據(jù)對比（相對差異均小于10%）

4 討論

近年來，隨著人口老齡化的加劇，老年性骨質(zhì)疏松癥已經(jīng)逐漸成為威脅人們生命健康的疾病，骨骼對力學(xué)響應(yīng)的能力下降是重要致病因素之一，因此，對于骨骼力學(xué)性能的研究已經(jīng)成為現(xiàn)代醫(yī)學(xué)研究的重要課題。骨骼生物力學(xué)測試是骨科基礎(chǔ)生物力學(xué)研究領(lǐng)域的重要組成部分，主要是利用工程學(xué)原理對骨骼和肌肉系統(tǒng)的力學(xué)性質(zhì)及功能進行檢驗和評價，直觀評估骨骼質(zhì)量，同時許多骨科器械和創(chuàng)新術(shù)式在應(yīng)用于臨床前都需要進行生物力學(xué)測試和仿真評估。骨骼承受的主要載荷是壓縮、拉伸、彎曲和扭轉(zhuǎn)，因此目前主要通過三點彎曲力學(xué)試驗、四點彎曲力學(xué)試驗、壓縮力學(xué)試驗和扭轉(zhuǎn)力學(xué)試驗等常規(guī)的骨骼生物力學(xué)測試獲取和分析壓縮、拉伸、彎曲和扭轉(zhuǎn)條件下的應(yīng)力應(yīng)變分布情況和彈性模量等力學(xué)指標[19]，其中壓縮測試一般用于負重骨的生物力學(xué)測試。本課題組針對雄性SD老年大鼠（22月齡）和青年大鼠（3月齡）尺骨在軸向壓縮載荷下的力學(xué)響應(yīng)進行研究，由于老年大鼠和青年大鼠尺骨形態(tài)大小不同，為控制恒定的應(yīng)變因素，應(yīng)首先明確達到相同應(yīng)變需施加的軸向壓縮載荷大小。目前應(yīng)變測量多使用有限元分析或應(yīng)變片測量的方式。有限元分析是一種準確、有效、低成本的方法，可對骨骼、肌肉和韌帶等組織以及植入物、輔助器械等醫(yī)療設(shè)備進行仿真，在骨骼生物力學(xué)領(lǐng)域被廣泛運用。該方法可模擬大鼠尺骨在軸向壓縮載荷下骨表面的應(yīng)變分布情況，有效彌補實驗所需嚴格的主客觀條件要求的不足。本研究使用Mimics Research、Geomagic Studio和SolidWorks軟件進行處理和三維逆向重建，使用Abaqus軟件對重建得到的三維有限元模型進行有限元分析。應(yīng)變片測量實驗可對骨骼表面應(yīng)變作實時動態(tài)測量，同時可對骨骼實體進行測量，較為準確，但影響因素較多，不易排除，會影響其精確度，而且僅能對應(yīng)變片粘貼位置的表面應(yīng)變進行研究，具有一定的局限性。目前還沒有更科學(xué)、更合理的技術(shù)提供應(yīng)變測量方法支撐，因此本研究采用2種方法相結(jié)合的方式探究軸向壓縮載荷下尺骨感興趣區(qū)域的表面應(yīng)變分布情況。

研究結(jié)果顯示，老年大鼠和青年大鼠尺骨表面應(yīng)變分布情況一致，最大壓縮應(yīng)變發(fā)生部位與前期研究一致[18]，均在約尺骨遠端1/3處，2種方法在相同應(yīng)變下測得所需的載荷大小的偏差在10%以內(nèi)。為使尺骨表面達到相同的壓縮微應(yīng)變，老年大鼠需要的載荷比青年大鼠大，且約為青年大鼠的2倍。雖然本研究結(jié)果與現(xiàn)實生活中隨著年齡的增長骨質(zhì)流失、骨骼強度下降，相同受力狀態(tài)下老年骨骼更易發(fā)生骨折的情況不相符，但本研究是以雄性SD大鼠為模型，與人類不同的是，SD大鼠成年后隨著年齡增長，體長仍會變長，骨骼結(jié)構(gòu)也會相應(yīng)發(fā)生改變，老年狀態(tài)下骨骼更長更粗，因此，雖然骨骼脆性增加，仍需要更大的載荷才能達到相同的應(yīng)變。

綜上所述，本研究方法可為應(yīng)力加載實驗提供指導(dǎo)，但仍存在一定的局限性。在有限元分析中，首先，本研究在理想狀態(tài)下將前肢簡化為尺骨和橈骨模型，考慮到生理狀態(tài)下骨骼表面有肌肉和韌帶等組織的覆蓋和與其他骨連接的情況，本模型的幾何結(jié)構(gòu)仍需進一步增強。其次，目前有限元分析中對材料屬性的賦值方法多是采用灰度賦值法和均一賦值法這2種方法?；叶荣x值法是在Mimics Research中根據(jù)CT圖像HU值不同通過經(jīng)驗公式計算并賦予彈性模量的值，這一方法雖較為簡單且更符合骨骼材質(zhì)分布特征，但需要通過參考文獻中對應(yīng)骨的經(jīng)驗公式[20]，且受CT圖像質(zhì)量的影響較大。本研究中采用的均一賦值法雖避免了灰度賦值法的缺陷，分別建立了皮質(zhì)骨、髓腔和松質(zhì)骨結(jié)構(gòu)的實體模型，并分別賦予對應(yīng)的材料屬性，但幾何結(jié)構(gòu)的建模存在較大的主觀性，且骨骼是各向異性材料，目前的有限元研究大都將骨骼簡化為分布均勻的各向同性材料[21]，因此會對結(jié)果的精確度產(chǎn)生一定的影響。在應(yīng)變片測量中，應(yīng)變片的粘貼也會影響測量數(shù)據(jù)的準確度，需要多次重復(fù)測量，增加了實驗的復(fù)雜性。后續(xù)研究會繼續(xù)優(yōu)化研究方法，以提高最終結(jié)果的精確度。