徐淑波，趙晨浩，劉建營(yíng)，張森，薛現(xiàn)猛，柳婷婷，林少宇

高孔隙度可再生骨支架仿真與實(shí)驗(yàn)研究

徐淑波1,2，趙晨浩1，劉建營(yíng)1，張森1，薛現(xiàn)猛1，柳婷婷1，林少宇1

（1.山東建筑大學(xué) 材料科學(xué)與工程學(xué)院，濟(jì)南 250101；2.華南理工大學(xué) 國(guó)家金屬材料近凈成形工程技術(shù)研究中心，廣州 510640）

確定單元體與股骨的最佳孔隙度骨支架結(jié)構(gòu)。通過(guò)掃描電鏡分析選區(qū)激光熔化（Selective laser melting，SLM）成形試樣的微觀結(jié)構(gòu)；通過(guò)靜力學(xué)模擬與實(shí)驗(yàn)分析不同孔隙度下標(biāo)準(zhǔn)結(jié)構(gòu)與Voronoi多孔結(jié)構(gòu)的壓縮變形規(guī)律；通過(guò)生物力學(xué)仿真實(shí)驗(yàn)分析步態(tài)周期下標(biāo)準(zhǔn)結(jié)構(gòu)與Voronoi多孔結(jié)構(gòu)的應(yīng)力分布情況。在選區(qū)激光熔化成形的316L不銹鋼微觀組織中，均勻分布著細(xì)小的近六邊形、伸長(zhǎng)六邊形的胞狀結(jié)構(gòu)和條柱狀亞結(jié)構(gòu)，受壓時(shí)有利于分散應(yīng)力，提高整體結(jié)構(gòu)的穩(wěn)定性；在壓縮變形時(shí)，標(biāo)準(zhǔn)結(jié)構(gòu)應(yīng)力集中于垂直棱柱，易導(dǎo)致棱柱斷裂引起試樣傾斜；Voronoi結(jié)構(gòu)連接桿的不均勻分布有利于分散應(yīng)力，使Voronoi結(jié)構(gòu)的最大等效應(yīng)力（250.34 MPa）遠(yuǎn)低于標(biāo)準(zhǔn)結(jié)構(gòu)（738.07 MPa），保證了整體受力均勻與結(jié)構(gòu)穩(wěn)定；在步態(tài)周期下，2種骨支架結(jié)構(gòu)的應(yīng)力隨孔隙度的增加而增加，75%孔隙度的Voronoi結(jié)構(gòu)具有更優(yōu)異的承壓能力與緩解應(yīng)力屏蔽的作用。通過(guò)模擬與試驗(yàn)分析，確定了單元體與股骨的最佳孔隙度及Voronoi結(jié)構(gòu)優(yōu)異的力學(xué)性能，驗(yàn)證了在步態(tài)周期下高孔隙度Voronoi骨支架結(jié)構(gòu)的可靠性，為股骨置換手術(shù)提供了理論依據(jù)。

選區(qū)激光熔化；靜力學(xué)仿真模擬；生物力學(xué)仿真實(shí)驗(yàn)；Voronoi多孔結(jié)構(gòu)；高孔隙度；步態(tài)周期

隨著激光選區(qū)熔化技術(shù)的發(fā)展，其在醫(yī)療上的應(yīng)用也越來(lái)越廣泛，不僅可根據(jù)患者的不同情況進(jìn)行精準(zhǔn)的個(gè)性化定制[1-2]，也可對(duì)孔的類型、孔徑尺寸、孔壁厚度及孔隙度進(jìn)行設(shè)計(jì)，完成特有的多孔結(jié)構(gòu)與胞格結(jié)構(gòu)制造[3]，在骨生物醫(yī)療領(lǐng)域有著重要的應(yīng)用前景[4-5]。

近年來(lái)，國(guó)內(nèi)外專家學(xué)者對(duì)SLM打印多孔結(jié)構(gòu)進(jìn)行了大量理論與實(shí)驗(yàn)研究[6-9]。崔新建等[10]采用SLM技術(shù)打印了BCC、FCC和FCCZ等3種多孔結(jié)構(gòu)，驗(yàn)證了該技術(shù)成形復(fù)雜多孔結(jié)構(gòu)件的可行性。高芮寧等[11]采用有限元分析與實(shí)驗(yàn)對(duì)比分析了多孔支架的力學(xué)性能，驗(yàn)證了多孔支架有限元分析的可行性。李姝博等[12]研究了單元體結(jié)構(gòu)對(duì)假體多孔結(jié)構(gòu)的微觀組織與力學(xué)性能的影響，并模擬了在不同多孔結(jié)構(gòu)下人體站立時(shí)髖關(guān)節(jié)假體載荷的傳遞與應(yīng)力分布規(guī)律。張?zhí)m等[13]建立了一種具有高孔隙率的骨小梁結(jié)構(gòu)，并采用有限元軟件分析了骨小梁的拉伸性能。

綜上所述，研究人員對(duì)多孔結(jié)構(gòu)件進(jìn)行了大量的理論與實(shí)驗(yàn)研究，但仍未考慮人體在實(shí)際步態(tài)周期下應(yīng)力分布的變化規(guī)律，且人體步態(tài)周期股骨的受力主要以壓應(yīng)力為主，因此有必要從壓縮變形方面探索不同孔隙度多孔結(jié)構(gòu)的變化規(guī)律。采用有限元分析與實(shí)驗(yàn)研究相結(jié)合，系統(tǒng)研究了不同孔隙度、不同骨支架結(jié)構(gòu)的等效應(yīng)力分布，并從生物載荷的角度驗(yàn)證了骨支架多孔結(jié)構(gòu)的可靠性，可為股骨置換手術(shù)提供理論依據(jù)。

1 研究對(duì)象與方法

1.1 單元體模型設(shè)計(jì)

1.1.1 單元體多孔結(jié)構(gòu)設(shè)計(jì)

為適應(yīng)天然骨小梁微觀結(jié)構(gòu)，設(shè)計(jì)了一種基于種子點(diǎn)的空間（平面）剖分結(jié)構(gòu)的不規(guī)則Voronoi多孔結(jié)構(gòu)。同時(shí)，為便于突出Voronoi多孔結(jié)構(gòu)的優(yōu)點(diǎn)，研究孔隙度對(duì)股骨力學(xué)性能的影響，設(shè)計(jì)了四棱柱標(biāo)準(zhǔn)多孔結(jié)構(gòu)作為對(duì)照組。通過(guò)大量研究發(fā)現(xiàn)[14-15]，骨組織支架結(jié)構(gòu)的理想多孔尺寸為300～900mm，孔隙度為60%～95%。將多孔支架結(jié)構(gòu)的孔尺寸設(shè)為300mm，孔隙度設(shè)為75%～85%，相關(guān)設(shè)計(jì)參數(shù)見(jiàn)表1。

表1 Voronoi孔型多孔結(jié)構(gòu)及對(duì)照組設(shè)計(jì)參數(shù)

Tab.1 Design parameters of Voronoi porous structure and control group

1.1.2 單元多孔模型建立

在Rhino中，利用插件Grasshopper設(shè)計(jì)隨機(jī)多孔結(jié)構(gòu)的建模流程見(jiàn)圖1。首先，在20 mm×20 mm區(qū)域隨機(jī)生成所需的種子點(diǎn)，通過(guò)相應(yīng)的概率球半徑生成二維Voronoi圖；然后，利用Voronoi 3D模塊生成含隨機(jī)單元體的三維Voronoi結(jié)構(gòu)（圖1b），并提取三維結(jié)構(gòu)內(nèi)部的不規(guī)則種子點(diǎn)，將相鄰的2點(diǎn)連接，構(gòu)建棱柱為1 mm×1 mm截面的隨機(jī)多孔模型（圖1c）；最后，用T-splines插件對(duì)棱柱光滑處理，生成棱柱直徑為1 mm截面的Voronoi多孔結(jié)構(gòu)模型(圖1d)。標(biāo)準(zhǔn)多孔結(jié)構(gòu)的建模需在二維建模區(qū)域生成規(guī)則點(diǎn)陣，才會(huì)得到相應(yīng)的三維結(jié)構(gòu)。

Voronoi多孔結(jié)構(gòu)建模流程的電池組圖見(jiàn)圖2?？赏ㄟ^(guò)設(shè)置體積為20 mm×20 mm×20 mm的正方體區(qū)域、種子點(diǎn)數(shù)量與棱柱直徑，創(chuàng)建出孔徑為1 mm、孔隙率為75%～85%的多孔三維結(jié)構(gòu)。將Voronoi多孔結(jié)構(gòu)以.xt格式導(dǎo)入有限元軟件Ansys，完成單元體多孔結(jié)構(gòu)的建模。

1.1.3 材料屬性設(shè)置

有限元模擬與實(shí)驗(yàn)的材料均選用316L不銹鋼，其具體應(yīng)力–應(yīng)變曲線見(jiàn)圖3，實(shí)驗(yàn)階段選用SLM技術(shù)制備多孔單元體。

1.1.4 邊界和約束條件

由于該多孔結(jié)構(gòu)應(yīng)用于人體骨骼，且人體的骨骼受力以壓應(yīng)力為主[16]。因此，應(yīng)對(duì)多孔結(jié)構(gòu)試樣進(jìn)行壓縮測(cè)試實(shí)驗(yàn)。通過(guò)Ansys模擬軟件對(duì)多孔模型上下表面加入厚度1 mm的實(shí)心板，實(shí)心板材料設(shè)為工程鋼材料，接觸面設(shè)為摩擦接觸。整體模型采用四面體網(wǎng)格劃分，網(wǎng)格尺寸設(shè)為0.3 mm。對(duì)底板實(shí)施完全固定約束，對(duì)上板施加100 N的載荷，載荷方向見(jiàn)圖4。為確保有限元分析的準(zhǔn)確性，保證多孔結(jié)構(gòu)植入人體后不產(chǎn)生塑性形變，其有限元力學(xué)分析僅限于線彈性形變。

圖2 隨機(jī)多孔結(jié)構(gòu)的Grasshopper主程序電池組圖

圖3 SLM技術(shù)制備316L不銹鋼壓縮應(yīng)力–應(yīng)變曲線

圖4 有限元壓縮模擬受力方向

1.2 股骨模型設(shè)計(jì)

1.2.1 股骨三維有限元模型建立

通過(guò)螺旋CT掃描機(jī)，實(shí)現(xiàn)人體股骨長(zhǎng)軸方向橫斷面的螺旋掃描。將股骨CT掃描數(shù)據(jù)以Dicom 3.0格式導(dǎo)入Mimics軟件，閾值設(shè)為662～161 3，并通過(guò)Calculate 3D界面創(chuàng)建股骨的三維模型（圖5a）。但此時(shí)的三維股骨模型由若干個(gè)曲面片組成（圖5b），故需將三維模型以STL格式保存并導(dǎo)入曲面重構(gòu)軟件Geomagic，通過(guò)表面裝填、縫合、封閉等操作將曲面片擬合成相連曲面，進(jìn)而生成NURBS曲面模型，實(shí)現(xiàn)股骨三維模型的光滑處理。以step的格式導(dǎo)入有限元軟件ABAQUS，完成股骨三維模型建模。

圖5 股骨實(shí)體模型

1.2.2 股骨與假體裝配

為設(shè)計(jì)出最符合骨小梁的多孔骨支架，將股骨中間部分用支架代替，在植入過(guò)程中保證假體中心線與股骨髓腔中心線吻合。假體與骨骼的接觸設(shè)置為“Bonded”，并在Mesh模塊中對(duì)股骨-假體各部分進(jìn)行自由網(wǎng)格劃分，單元設(shè)為四面體形狀。根據(jù)不同孔隙度的骨支架單元體結(jié)構(gòu)建立股骨-假體模型（圖6）。

圖6 股骨–假體模型

1.2.3 材料屬性設(shè)置

在模擬階段假定骨結(jié)構(gòu)、軟骨為連續(xù)、均勻、各向同性的彈性材料[17-20]，股骨各部分模型具體材料屬性見(jiàn)表2。

1.2.4 股骨約束

采用anybody軟件模擬體質(zhì)量為70 kg的成人完成整個(gè)行走周期[21-22]，忽略其余因素影響，導(dǎo)出完整步態(tài)周期下股骨反作用力（HMF）、髖關(guān)節(jié)后側(cè)力（HAPH）和髖關(guān)節(jié)近端力（HPDF）的載荷數(shù)據(jù)，具體的方向及大小見(jiàn)圖7。通過(guò)ABAQUS載荷模塊建立3組數(shù)據(jù)的載荷模型，同時(shí)將遠(yuǎn)端股骨節(jié)點(diǎn)固定約束，模擬步態(tài)周期下股骨受力對(duì)假體的影響。

表2 模型材料屬性

Tab.2 Model material properties

圖7 股骨近端各載荷的角度與大小

2 結(jié)果分析

2.1 單元體有限元結(jié)果分析

不同孔隙度多孔結(jié)構(gòu)壓縮試樣的等效應(yīng)力圖見(jiàn)圖8。由圖8可知，標(biāo)準(zhǔn)結(jié)構(gòu)的最大等效應(yīng)力分別為738.07、1 014.4、1 445.2 MPa，而Voronoi結(jié)構(gòu)的最大等效應(yīng)力分別為250.34、425.58、876.2 MPa?？梢?jiàn)，在相同加載力下，Voronoi結(jié)構(gòu)的最大等效應(yīng)力遠(yuǎn)低于標(biāo)準(zhǔn)結(jié)構(gòu)。對(duì)試樣各部分的應(yīng)力分析可知，隨著孔隙度增大，有限元模型的等效應(yīng)力逐漸增大，應(yīng)力對(duì)上下板的影響亦逐漸增大，這將導(dǎo)致骨支架與骨骼接觸面的應(yīng)力分布不均勻，從而對(duì)骨組織造成更為嚴(yán)重的應(yīng)力屏蔽。對(duì)單元體結(jié)構(gòu)的分析可知，標(biāo)準(zhǔn)結(jié)構(gòu)的應(yīng)力主要分布在垂直固定面的棱柱上，這說(shuō)明試樣破壞時(shí)垂直棱柱會(huì)先行斷裂，且由于棱柱的不均勻受力，試樣易產(chǎn)生傾斜而對(duì)人體造成二次損傷。Voronoi結(jié)構(gòu)的應(yīng)力分布在每個(gè)桿件連接處，試樣整體的應(yīng)力分布更為均勻，因而試樣在壓縮時(shí)并無(wú)傾斜現(xiàn)象發(fā)生。

2.2 實(shí)驗(yàn)結(jié)果分析

為分析多孔結(jié)構(gòu)在打印過(guò)程中的組織演變過(guò)程，驗(yàn)證有限元分析的可靠性，對(duì)標(biāo)準(zhǔn)結(jié)構(gòu)與Voronoi結(jié)構(gòu)進(jìn)行SEM分析與力學(xué)性能測(cè)試。需完成不同孔隙度單元體多孔結(jié)構(gòu)的打印，試驗(yàn)材料選用實(shí)驗(yàn)室316L不銹鋼粉末，釆用FS271M金屬激光3D打印機(jī)快速打印成形，打印試樣見(jiàn)圖9。

2.2.1 多孔316L不銹鋼微觀組織分析

SLM成形試樣的SEM圖見(jiàn)圖10。由圖10可知，晶粒內(nèi)部均勻分布著近六邊形、伸長(zhǎng)六邊形的胞狀結(jié)構(gòu)和條柱狀亞結(jié)構(gòu)，且相比于傳統(tǒng)316L不銹鋼，組織中并未發(fā)現(xiàn)孿晶。柱狀亞結(jié)構(gòu)條帶之間的寬度為400～500 nm（圖 10a），胞狀組織結(jié)構(gòu)的平均尺寸為1 μm左右（圖 10b）。這些細(xì)小亞結(jié)構(gòu)在承受壓縮應(yīng)力時(shí)可由每個(gè)部位均勻分?jǐn)?，能夠承受更高的壓縮應(yīng)力。通過(guò)對(duì)不銹鋼材料的大量調(diào)研發(fā)現(xiàn)[23-26]，316L不銹鋼常被用作人體牙齒與骨骼的植入，能有效促進(jìn)骨組織再生且不產(chǎn)生排異反應(yīng)，具有良好的生物相容性。因此，采用3D打印成形的316L不銹鋼結(jié)構(gòu)作為人骨的實(shí)驗(yàn)材料具有較優(yōu)越的組織優(yōu)勢(shì)。

圖8 不同孔隙度多孔結(jié)構(gòu)有限元模擬等效應(yīng)力圖

圖9 SLM技術(shù)打印試樣

圖10 316L不銹鋼試樣SEM圖像

2.2.2 多孔316L不銹鋼壓縮性能分析

分別對(duì)變形量在2%、10%、30%時(shí)的標(biāo)準(zhǔn)結(jié)構(gòu)與Voronoi結(jié)構(gòu)進(jìn)行取樣分析（圖11），2種多孔結(jié)構(gòu)與實(shí)體金屬類似，在壓縮過(guò)程中均變成鼓形，體現(xiàn)出良好的塑性變形能力。但由于標(biāo)準(zhǔn)結(jié)構(gòu)在壓縮過(guò)程中連接上、下層的棱柱斷裂，導(dǎo)致整體試樣出現(xiàn)偏移現(xiàn)象；Voronoi結(jié)構(gòu)的連接桿呈不均勻分布，可用于分散應(yīng)力，保證試樣受力均勻，不易引起試樣傾斜?？梢?jiàn)，有限元仿真與實(shí)驗(yàn)結(jié)果相同，Voronoi結(jié)構(gòu)可有效反應(yīng)試樣的實(shí)際受力情況。

圖11 不同變形量多孔結(jié)構(gòu)試樣壓縮實(shí)驗(yàn)

多孔結(jié)構(gòu)壓縮應(yīng)力–應(yīng)變曲線見(jiàn)圖12。與圖10所示的壓縮過(guò)程相似，多孔結(jié)構(gòu)的壓縮變形分為彈性、塑性與壓縮密實(shí)等3個(gè)階段。其中，彈性階段的應(yīng)力–應(yīng)變呈線性增長(zhǎng)；塑性階段由于多孔結(jié)構(gòu)壓縮后晶粒的不均勻性，試樣局部發(fā)生塑性變形，應(yīng)力–應(yīng)變曲線呈現(xiàn)平臺(tái)狀；壓縮密實(shí)階段由于變形和作用力方向的不同，會(huì)產(chǎn)生一定的剪切應(yīng)力，最終在孔壁尖角形成應(yīng)力集中，從而造成裂紋源。隨著應(yīng)力的逐漸增加，裂紋逐漸擴(kuò)展，致使整體結(jié)構(gòu)完全斷裂。整個(gè)變形過(guò)程無(wú)明顯的屈服階段。

圖12 多孔結(jié)構(gòu)壓縮應(yīng)力–應(yīng)變曲線

由圖12還可知，標(biāo)準(zhǔn)件在壓縮密實(shí)階段應(yīng)力出現(xiàn)明顯的偏移現(xiàn)象，這說(shuō)明應(yīng)力集中于垂直連接柱處，若植入到人體中極易產(chǎn)生應(yīng)力集中現(xiàn)象，對(duì)人體造成二次損傷；而Voronoi結(jié)構(gòu)壓縮曲線在經(jīng)過(guò)屈服點(diǎn)后依舊呈線性增長(zhǎng)趨勢(shì)，且增長(zhǎng)率非常小，壓縮過(guò)程并未產(chǎn)生偏移，無(wú)明顯應(yīng)力集中。另外，壓縮應(yīng)力隨孔隙度的減小而不斷增加，孔隙度為75%的Voronoi結(jié)構(gòu)具有更好的承壓能力，因而Voronoi結(jié)構(gòu)比標(biāo)準(zhǔn)結(jié)構(gòu)更適合作為人體骨骼植入物。

2.3 多孔316L不銹鋼生物力學(xué)仿真分析

股骨不同孔隙度、不同結(jié)構(gòu)的等效應(yīng)力云圖見(jiàn)圖13。由圖13可知，在步態(tài)周期下主要應(yīng)力為股骨反作用力，這說(shuō)明股骨在步態(tài)周期主要受壓縮力，驗(yàn)證了壓縮實(shí)驗(yàn)的必要性。股骨生物力學(xué)實(shí)驗(yàn)數(shù)據(jù)與單元體壓縮變形相似，其最大等效應(yīng)力隨孔隙的增大逐漸增加，且最大應(yīng)力點(diǎn)位置逐漸下移；同時(shí)，Voronoi結(jié)構(gòu)的最大等效應(yīng)力遠(yuǎn)低于標(biāo)準(zhǔn)結(jié)構(gòu)，最大應(yīng)力點(diǎn)位置較高，且試樣應(yīng)力分布更為均勻，這表明Voronoi結(jié)構(gòu)具有更優(yōu)異的抗壓縮變形能力。

圖13 股骨不同孔隙度不同結(jié)構(gòu)的Mises應(yīng)力云圖

對(duì)比圖13a—c或d—f可知，孔隙度越高，金屬假體對(duì)骨組織造成的應(yīng)力屏蔽越嚴(yán)重，假體與骨組織接觸面的應(yīng)力也逐漸增大，易發(fā)生假體與骨組織間的摩擦作用，從而減少植入物的使用壽命。因此，75%孔隙度的Voronoi結(jié)構(gòu)最適合作為人骨植入物。

3 結(jié)論

通過(guò)有限元模擬與試驗(yàn)研究，獲得了不同孔隙度、不同結(jié)構(gòu)下單元體多孔結(jié)構(gòu)與股骨-假體多孔結(jié)構(gòu)的壓縮變形規(guī)律。

1）SLM打印的316L不銹鋼均由1 μm以下的細(xì)小胞狀結(jié)構(gòu)、柱狀亞結(jié)構(gòu)構(gòu)成，在壓縮變形時(shí)可均勻分?jǐn)倯?yīng)力，保證整體結(jié)構(gòu)穩(wěn)定。

2）多孔結(jié)構(gòu)壓縮大致可分為彈性階段、塑性平臺(tái)階段和壓縮密實(shí)階段；在壓縮過(guò)程中，標(biāo)準(zhǔn)多孔結(jié)構(gòu)應(yīng)力集中在縱向棱柱處，易造成桿件斷裂而引起試樣傾斜；Voronoi多孔結(jié)構(gòu)的應(yīng)力集中在多桿交點(diǎn)處，試樣受力均勻且無(wú)傾斜現(xiàn)象。

3）對(duì)6種不同骨支架多孔結(jié)構(gòu)的生物力學(xué)仿真實(shí)驗(yàn)分析可知，等效應(yīng)力隨孔隙度的減小而逐漸減弱，Voronoi多孔結(jié)構(gòu)能有效降低人體步態(tài)周期下的應(yīng)力屏蔽效應(yīng)，減少多孔支架與骨組織的摩擦，75%孔隙度的Voronoi多孔結(jié)構(gòu)植入人骨具有更優(yōu)異的性能。

[1] 吳小寶, 馮坤, 田大為. 3D打印個(gè)性化導(dǎo)航截骨導(dǎo)板在膝關(guān)節(jié)置換中的應(yīng)用研究進(jìn)展[J]. 中國(guó)醫(yī)藥導(dǎo)報(bào), 2022, 19(4): 52-56.

WU Xiao-bao, FENG Kun, TIAN Da-wei. Research Progress on the Application of 3D Printing Personalized Navigation Osteotomy Guide in Knee Replacement[J]. China Medical Herald, 2022, 19(4): 52-56.

[2] 魏志遠(yuǎn), 李浩, 熊旭，等. 雙峰3D打印個(gè)性化接骨板內(nèi)固定治療嚴(yán)重脛骨平臺(tái)骨折[J]. 中國(guó)骨與關(guān)節(jié)損傷雜志, 2021, 36(10): 1087-1089.

WEI Zhi-yuan, LI Hao, XIONG Xu, et al. Internal Fixation of Bimodal 3D Printing Personalized Bone Plate for Severe Tibial Plateau Fractures[J]. Chinese Journal of Bone and Joint Injury, 2021, 36(10): 1087-1089.

[3] 宋穎, 鄧久鵬, 張煒. 骨組織工程多孔支架制備方法的研究進(jìn)展[J]. 生物醫(yī)學(xué)工程與臨床, 2021, 25(2): 246-250.

SONG Ying, DENG Jiu-peng, ZHANG Wei. Research Progress in the Preparation of Porous Scaffolds for Bone Tissue Engineering[J]. Biomedical Engineering and Clinical, 2021, 25(2): 246-250.

[4] FANTINI M, CURTO M, CRESCENZIO D F. A Method to Design Biomimetic Scaffolds for Bone Tissue Engineering Based on Voronoi Lattices[J]. Virtual Physi Prototyping, 2016, 11(2): 77-90.

[5] MACHADO B I, LOPEZ M I, BRACKE J, et al. Next-generation Biomedical Implants Using Additive Manufacturing of Complex, Cellular and Functional Mesh Arrays[J]. Philosophical Transactions of the Royal Society, 2010, 368(1917): 1999-2032.

[6] Gomez S, Vlad M D, Lopez J, et al. Design and Properties of 3D Scaffolds for Bone Tissue Engineering[J]. Acta Biomaterialia, 2016, 42: 341-350.

[7] Tan X P, Tan Y J, Chow C S L, et al. Metallic Powder-bed Based 3D Printing of Cellular Scaffolds for Orthopaedic Implants: A State-of-the-art Review on Manufacturing, Topological Design, Mechanical Properties and Biocompatibility[J]. Materials Science & Engineering, C. Materials for Biogical Applications, 2017, 76: 1328-1343.

[8] Fantini M, Curto M, De Crescenzio F. A Method to Design Biomimetic Scaffolds for Bone Tissue engineering Based on Voronoi Lattices[J]. Virtual & Physical Prototyping, 2016, 11(2): 77-90.

[9] Fantini M, Curto M. Interactive Design and Manufacturing of a Voronoi-based Biomimetic Done Scaffold for Morphological Characterization[J]. International Journal on Interactive Design & Manufacturing, 2018, 12(2): 585-596.

[10] 崔新建, 劉玉德, 石文天，等. 選區(qū)激光熔化成形316L不銹鋼多孔結(jié)構(gòu)及力學(xué)性能[J]. 特種鑄造及有色合金, 2020, 40(1): 75-80.

CUI Xin-jian, LIU Yu-de, SHI Wen-tian, et al. Porous Structure and Mechanical Properties of 316L Stainless Steel Formed by Selective Laser Melting[J]. Special Casting and Nonferrous Alloys, 2020, 40(1): 75-80.

[11] 高芮寧, 李祥. 徑向梯度多孔支架設(shè)計(jì)與力學(xué)性能分析[J]. 機(jī)械工程學(xué)報(bào), 2021, 57(3):220-226.

GAO Rui-ning, LI Xiang. Design and Mechanical Properties Analysis of Radial Gradient Porous Scaffolds[J]. Chinese Journal of Mechanical Engineering, 2021, 57(3):220-226.

[12] 李姝博. 3D打印髖關(guān)節(jié)假體多孔結(jié)構(gòu)力學(xué)性能的數(shù)值仿真和實(shí)驗(yàn)研究[D]. 吉林: 吉林大學(xué), 2020: 51-57.

LI Shu-bo. Numerical Simulation and Experimental Study on Mechanical Properties of 3D Printed Hip Prosthesis Porous Structure[D]. Ji Lin: Jilin University, 2020: 51-57.

[13] 張?zhí)m, 王翔, 劉軍，等. 3D打印鈦合金骨小梁多孔結(jié)構(gòu)的拉伸性能[J]. 中國(guó)組織工程研究, 2020, 24(22): 3498-3503.

ZHANG Lan, WANG Xiang, LIU Jun, et al. Tensile Properties of 3D Printed Titanium Alloy Trabecular Bone Porous Structures[J]. Chinese Journal of Tissue Engineering Research, 2020, 24(22): 3498-3503.

[14] 馮雨欣, 王玥涵, 曲佳欣，等. 基于CT的股骨精確建模與三維有限元分析[J]. 軟件, 2020, 41(1): 46-49.

FENG Yu-xin, WANG Yue-han, QU Jia-xin, et al. Accurate Modeling and 3D Finite Element Analysis of Femur Based on CT[J]. Software, 2020, 41(1): 46-49.

[15] 孫海波, 徐淑波, 張森，等. SLM成形不同孔隙結(jié)構(gòu)骨支架的仿真與實(shí)驗(yàn)研究[J]. 精密成形工程, 2022, 14(2): 123-128.

Sun Hai-bo, Xu Shu-bo, Zhang Sen, et al. Simulation and Experimental Study on SLM Forming of Bone Scaffolds with Different Pore Structures[J]. Precision Forming Engineering, 2022, 14(2): 123-128.

[16] 任啟超. 正常重力環(huán)境下人體下肢力學(xué)特性研究[D]. 太原理工大學(xué), 2017: 23-45.

REN Qi-chao. Research on mechanical properties of human lower limbs in normal gravity environment[D]. Taiyuan University of Technology, 2017: 23-45.

[17] 石志良, 黃琛, 盧小龍，等. Sphere梯度孔結(jié)構(gòu)力學(xué)性能有限元分析[J]. 中國(guó)生物醫(yī)學(xué)工程學(xué)報(bào), 2019, 38(3): 339-347.

SHI Zhi-liang, HUANG Chen, Lu Xiao-long, et al. Finite Element Analysis of Mechanical Properties of Sphere Gradient Holes[J]. Chinese Journal of Biomedical Engineering, 2019, 38(3): 339-347.

[18] 王學(xué)斌, 龐清江, 余霄. 人工股骨頭置換術(shù)中置入不同直徑假體球頭的有限元分析[J]. 中國(guó)骨傷, 2020, 3(6): 558-563.

WANG Xue-bin, PANG Qing-jiang, YU Xiao. Finite Element Analysis of Different Diameter Ball Head Implant in Artificial Femoral Head Replacement[J]. China Journal of Bone Trauma, 2020, 3(6): 558-563.

[19] 張升. 醫(yī)用合金粉末激光選區(qū)熔化成形工藝與性能研究[D]. 武漢: 華中科技大學(xué), 2014: 29-57．

Zhang Sheng. Laser Selective Melting Forming of Med-ical Alloy Powder and Its Properties[D]. Wuhan: Huazhong University of Science and Technology, 2014: 29-57.

[20] 顏繼英. 不同材料賦值下股骨靜力學(xué)有限元模型的力學(xué)仿真分析[J]. 中國(guó)組織工程研究, 2020, 24(9): 1390-1394.

YAN Ji-ying. Mechanical Simulation Analysis of Femur Static Finite Element Model with Different Material Assignment[J]. Chinese Journal of Tissue Engineering Research, 2020, 24(9): 1390-1394.

[21] 徐歡歡, 何育民, 孫朝陽(yáng)，等. AnyBody環(huán)境下人體步態(tài)的逆向動(dòng)力學(xué)研究[J]. 機(jī)械科學(xué)與技術(shù), 2019, 38(12): 1819-1824.

XU Huan-huan, HE Yu-min, SUN Chao-yang, et al. Study on Inverse Dynamics of Human Gait in the Environment of AnyBody[J]. Mechanical Science and Technology, 2019, 38(12): 1819-1824.

[22] 羅林聰, 馬立敏, 林澤，等. 基于AnyBody骨骼肌肉多體動(dòng)力學(xué)分析的有限元仿真[J]. 醫(yī)用生物力學(xué), 2019, 34(3): 237-242.

LUO Lin-cong, MA Li-min, Lin Ze, et al. Finite Element Simulation of Multibody Dynamics Analysis of Skeletal Muscle of AnyBody[J]. Medical Biomechanics, 2019, 34(3): 237-242.

[23] LAWRENCE J, CHEW H R, CHONG C K, et al. Laser Modification of the Wettability Characteristics of a 316L Stainless Steel Bio-metal and the Effects Thereof on Human Fibroblast Cell Response[J]. Lasers in Engineering, 2005, 15: 75-90.

[24] 吳啟銳. 選區(qū)激光熔化成形316L不銹鋼多孔結(jié)構(gòu)的力學(xué)性能及生物相容性研究[D]. 福建: 福建工程學(xué)院, 2021: 56-62.

WU Qi-rui. Study on Mechanical Properties and Biocompatibility of 316L Stainless Steel Porous Structure Formed by Selective Laser Melting[D]. Fu Jian: Fujian Institute of Technology, 2021: 56-62.

[25] 李超. 激光增材制造成形316L不銹鋼的研究進(jìn)展[J]. 機(jī)械工程材料, 2022, 46(8): 1-7.

LI Chao. Progress of laser additive manufacturing for forming 316L stainless steel [J]. Mechanical Engineering Materials, 2022, 46(8): 1-7.

[26] 謝焯俊, 方輝, 沈顯峰，等. 面向基于Voronoi結(jié)構(gòu)的隨機(jī)多孔骨植入材料的小孔結(jié)構(gòu)參數(shù)化建模方法研究[J]. 生物骨科材料與臨床研究, 2020, 17(4): 1-6.

XIE Chao-jun, FANG Hui, SHEN Xian-feng, et al. A parametric modeling method for small pore structures of random porous bone implant materials based on Voronoi structure[J]. Bio-Orthopaedic Materials and Clinical Research, 2020, 17(4): 1-6.

Simulation and experiment of high porosity renewable bone scaffold

XU Shu-bo1,2, ZHAO Chen-hao1, LIU Jian-ying1, ZHANG Sen1, XUE Xian-meng1, LIU Ting-ting1, LIN Shao-yu1

(1. School of Materials Science and Engineering, Shandong Jianzhu University, Jinan 250101, China; 2. National Engineering Technology Research Center for Material Forming, South China University of Technology, Guangzhou 510640, China)

The work aims to determine the optimal porosity bone scaffold structure between the unit body and the femur.The microstructure of the sample formed by selective laser melting (SLM) was analyzed by scanning electron microscope. The compression deformation law of standard structure and Voronoi porous structure under different porosity was analyzed by static simulation and experiment. The stress distribution of the standard structure and the Voronoi porous structure under the gait cycle was analyzed by biomechanical simulation experiment.The microstructure of 316L stainless steel formed by selective laser melting was evenly distributed with fine cell structures of near hexagon, elongated hexagon and strip columnar substructure, which was conducive to the stress dispersion and improved the stability of the whole structure. During the compression deformation, the stress of the standard structure was concentrated in the vertical prism, which easily led to the fracture of the prism and caused the sample tilt. The non-uniform distribution of Voronoi structure was conducive to the stress dispersion, so that the maximum equivalent stress of Voronoi structure (250.34 MPa) was much lower than that of the standard structure (738.07 MPa), ensuring the uniform stress and stability of the whole structure. Under the gait cycle, the stress of the two bone scaffold structures increased with the increase of porosity, and the Voronoi structure with 75% porosity had better bearing capacity and relieved the effect of stress shielding. Through simulation and experimental analysis, the optimal porosity of the unit body and the femur and the excellent mechanical properties of the Voronoi structure are determined and the reliability of the high porosity Voronoi bone scaffold structure under the gait cycle is verified, which provides a theoretical basis for femoral replacement surgery.

selective laser melting; static simulation; biomechanical simulation experiment; Voronoi porous structure ; high porosity; gait cycle

10.3969/j.issn.1674-6457.2023.01.016

R318.08

A

1674-6457(2023)01-0120-08

2022–04–16

2022-04-16

國(guó)家自然科學(xué)基金（41305124）；山東省科技型中小企業(yè)創(chuàng)新能力提升工程（2022TSGC2108，2022TSGC2402）；山東省自然科學(xué)基金面上項(xiàng)目（ZR2021ME182）；大學(xué)生國(guó)創(chuàng)計(jì)劃（S20211043001，202210430010，202210430008）

Natural Science Foundation of Shandong Province (41305124); Science and Technology Enterprise Innovation Program of Shandong Province (2022TSGC2108, 2022TSGC2402); Natural Science Foundation of Shandong Province (ZR2021ME182); National College Student Innovation and Entrepreneurship Program (S20211043001, 202210430010, 202210430008)

徐淑波（1977—），男，博士，教授，主要研究方向?yàn)樗苄约庸ず驮霾闹圃臁?/p>

XU Shu-bo (1977-), Male, Doctor, Professor, Research focus: plastic forming and additive manufacturing.

徐淑波, 趙晨浩, 劉建營(yíng), 等. 高孔隙度可再生骨支架仿真與實(shí)驗(yàn)研究[J]. 精密成形工程, 2023, 15(1): 120-127.

XU Shu-bo, ZHAO Chen-hao, LIU Jian-ying, et al. Simulation and experiment of high porosity renewable bone scaffold[J]. Journal of Netshape Forming Engineering, 2023, 15(1): 120-127.