顧春玉 錢煒

摘 ?要：為探索皮質骨鉆削工藝參數(shù)對鉆削軸向力的影響，防止臨床手術過程中因鉆削力過大造成骨損傷，基于專業(yè)金屬切削軟件Advantage對麻花鉆鉆削皮質骨進行了有限元分析，并通過皮質骨鉆削實驗驗證了有限元分析的準確性。此外，在仿真數(shù)據(jù)的基礎上，通過多元線性回歸得到了皮質骨鉆削力的預測數(shù)學模型。結果表明：隨著鉆頭直徑和進給速度的增大，鉆削力增大;隨著轉速的增加，鉆削力減小。因此，在骨鉆孔過程中應選取較小的鉆頭直徑與進給速度，并適當提高轉速，以提高鉆骨手術的質量。

關鍵詞：Advantage;皮質骨;鉆削力;預測模型

中圖分類號：TP391.9 ? ? 文獻標識碼：A

Abstract： In order to explore the influence of cortical bone drilling process parameters on the axial force of drilling， and to prevent bone damage caused by excessive drilling force during clinical surgery， this paper proposes a finite element analysis of twist drill drilling cortical bone based on Advantage—professional metal cutting software. The accuracy of the finite element analysis is verified by the cortical bone drilling experiment. In addition， on the basis of simulation data， a predictive mathematical model of cortical bone drilling force is obtained through multiple linear regression. Results show that the drilling force increases with the increase of drill diameter and feed rate. With the increase of rotating speed， the drilling force decreases. Therefore， in the process of bone drilling， smaller diameter and feed speed should be selected， and the rotating speed should be increased appropriately to improve the quality drill of bone drilling surgery.

Keywords： Advantage; cortical bone; drilling force; prediction model

1 ? 引言（Introduction）

在傳統(tǒng)的骨科手術中，通常對骨材料進行鉆孔、銑削和鋸切等類似于工業(yè)制造應用的操作[1-2]。而在這些加工操作中，往往都要求骨材料在加工過程中具有低損傷、最小侵入時間和時間成本等要求[3]。其中，鉆孔在骨科手術過程中十分常見，鉆削力過高會對鉆孔周圍組織造成損傷，同時也會導致鉆削溫度升高，因此對骨骼鉆孔過程鉆削力的研究有著十分重要的意義。

很多學者對這方面進行了研究。楊毅欣等人[4]通過改變普通麻花鉆的頂角與螺旋角，降低了在相同工藝參數(shù)下的豬脛骨鉆削軸向力和鉆削溫度;白小帆等人[5]通過低頻軸向振動的輔助方式輔助骨鉆削加工，降低了傳統(tǒng)骨鉆削的鉆削軸向力和鉆削溫度;婁德鈺[6]通過對傳統(tǒng)克氏針進行分屑槽設計改變，設計了具有分屑槽和前角的新型克氏針，并研究了設計參數(shù)對克氏針鉆削皮質骨的力熱影響;ALAM等人[7]通過有限元分析仿真了皮質骨平面切割力，并進行了相應的實驗與數(shù)值分析，同時對骨材料切屑進行了觀測。

在這些研究分析中，他們大多著重優(yōu)化骨鉆孔的鉆削刀具參數(shù)和輔助加工方式，從而提高鉆骨質量，并沒有對工藝參數(shù)的影響進行深入研究。本文通過專業(yè)切削軟件Advantage進行皮質骨鉆削仿真，分析工藝參數(shù)對軸向力的影響，并通過實驗驗證模型準確性，在正交仿真實驗數(shù)據(jù)的基礎上，通過多元線性回歸模型建立了皮質骨鉆削力的預測數(shù)學模型，為骨鉆孔手術提供一定的理論指導。

2 ?Advantage三維皮質骨鉆削模型（3-D model of cortical bone drilling quality based on Advantage）

2.1 ? 鉆削三維模型的建立

Advantage切削仿真分為2D與3D兩大類，由于鉆削加工垂直方向上的進給運動，因此只能選用3D分析。在模型建立的開始，要求設定鉆削狀態(tài)，主要分為Entrance、Starting Depth、Exit，分別表示從外部鉆入工件、一定深度內的鉆孔和鉆出工件三種鉆削狀態(tài)。本文主要研究麻花鉆頭鉆入骨骼表層后，在鉆削穩(wěn)定狀態(tài)下的鉆削力和鉆削溫度狀況。由于鉆尖的頂角具有一定的長度，鉆頭需要在一定進給后才能完全鉆入工件，因此選用Starting Depth分析模式，使鉆頭在鉆削開始時即沒入工件，切削刃與工件完全接觸，如圖1（a）所示。由于人的股骨形狀并非規(guī)則圓柱，可假定從其內部取一個規(guī)則小塊進行分析，如圖1（b）所示，因此工件設定為規(guī)則長方體，長寬均為8 mm，高度為3 mm。

2.2 ? 皮質骨鉆削材料參數(shù)

皮質骨是一種天然復合材料，從橫觀來看它具有各向同性的特點，從正交方向上來看它又有各向異性的特點。本文將皮質骨作為各向同性材料進行研究，鉆頭采用通用醫(yī)療麻花鉆，其材料參數(shù)如表1所示。

材料的宏觀本質結構運動關系主要反映了一定微觀變形速度條件下的各種宏觀摩擦應力隨著微觀應變的速度變化而發(fā)生變化的運動規(guī)律，而這些變化規(guī)律又直接受到微觀應變、應變運動速率、溫度、時間、摩擦應力系數(shù)等各種因素的直接影響，因此選用合適的本構模型對于有限元仿真有著十分重要的意義[8]。Advantage軟件自帶有Power Law（冪指數(shù)材料模型）、Drucker Prager和用戶自定義三種材料本構模型設置。Johnson-Cook本構模型是一種描述材料在高溫度應變、大尺度變形等多種情況下溫度比較理想的溫度本體結構仿真模型，它不僅可以直接反映三維材料應變率硬化特性和應變硬化特性，也可以直接表現(xiàn)三維材料的應變溫度本體特點，在三維材料切削過程仿真中應用較廣，所以本文主要選用Johnson-Cook本構模型作為皮質骨骼類材料的溫度本體結構仿真模型。該模型的應力流動量和應力系數(shù)表達式結構如式（1）所示，皮質骨的Johnson-Cook模型參數(shù)如表2所示。

2.3 ? 鉆頭模型的建立

通用國際標準麻花鉆的鉆芯幾何形狀一般都是固定的，我國有關國家標準和國際有關工業(yè)刀具機械設計使用手冊一般都要給出六個幾何形狀描述性的參數(shù)：直徑、鉆芯厚度、頂角、螺旋角、外緣后角以及橫刃傾角[9]。Advantage軟件在鉆削仿真中可以通過設定標注麻花鉆參數(shù)自動生成麻花鉆頭。本文選用的鉆頭直徑為2—4 mm，螺旋角28°，橫刃傾角125°，頂角118°，如圖2（a）所示;自動生成的標準麻花鉆頭如圖2（b）所示。其余默認設置不變。

2.4 ? 鉆削工藝參數(shù)的選定

臨床手術中對骨骼的鉆削大多采用手持電鉆進行鉆孔，通常對骨骼進行鉆削的手持電鉆主軸轉速在空載狀態(tài)下不超過3，000 r/min[10]。對于鉆骨的進給速度，由于是手術操作者人為對手持電鉆施加壓力進行進給，過高的進給速度會增加手持電鉆進行手術的難度，且對骨骼造成較大的損傷，因此鉆骨手術中的進給速度不宜過高。綜合這些因素，參照其他學者對骨骼進行鉆削實驗研究過程中選用的切削速度和進給速度，在本文實驗過程中設定的切削參數(shù)如表3所示。在Advantage中，需要將進給速度轉變?yōu)槊魁X進給量，并設置初始深度與工件初始溫度。由于該分析是在一定深度下進行的，因此只需觀察鉆削的某一瞬間即可，不需要長時間鉆削，經過調試在旋轉90°之后軸向力趨于穩(wěn)定，并可在一定程度上節(jié)省分析時間。

3 ?皮質骨鉆削實驗（Experiments on cortical bone drilling）

3.1 ? 實驗試樣的制備

本文試樣選取牛股骨，因為牛骨的力學性能與人骨的性能最為相似。將牛股骨表面及內部軟組織清洗干凈，選取中間段圓柱部分，通過電鋸進行切段，確保斷面水平光滑，如圖3所示，以方便夾具夾持。

3.2 ? 實驗儀器及實驗設計

主要實驗儀器如表4所示。在多功能鉆床上通過壓板固定Kistler測力儀，在測力儀上安裝夾具，夾持皮質骨試樣。采用2 mm直徑鉆頭鉆通皮質骨，分別在定轉速為600 rpm與進給速度為60 mm/min的條件下進行測力實驗，每組重復進行三次，通過電腦采集鉆削力的變化情況。實驗設備和實驗流程如圖4所示。

4 ? 仿真與實驗結果分析（Analysis of simulation and experimental results）

圖5為皮質骨鉆削的三維有限元仿真圖，骨鉆在一定深度內按預定工藝參數(shù)鉆削，獲得皮質骨鉆削動態(tài)瞬時鉆削力，圖6、圖7為仿真結果。圖6表明在轉速不變的情況下，皮質骨鉆削軸向力隨著進給速度的增大而增大，大小為20—50 N;圖7表明在進給速度不變的情況下，皮質骨鉆削力隨著轉速的增大而減小，大小為20—40 N。對比皮質骨鉆削實驗數(shù)據(jù)，如表5所示，皮質骨鉆削力的仿真值與實際實驗值相差在5 N之內，仿真誤差在10%以內，具有可靠性。

5 ?皮質骨鉆削力預測模型的建立（Establishment of a prediction model for cortical bone drilling force）

5.1 ? 正交仿真實驗設計

正交實驗設計是研究多因素多水平的一種設計方法，它是根據(jù)正交性從全面實驗中挑選出部分有代表性的組合進行實驗，這些具有代表性的實驗具備“均勻分散，齊整可比”的特點。正交實驗設計是分式析因設計的主要方法，是一種高效率、快速、經濟的實驗設計方法[11]。通過數(shù)據(jù)分析軟件SPSS，根據(jù)上文中選取的工藝參數(shù)水平，設計正交實驗如表6所示，并根據(jù)實驗號進行仿真實驗，結果如下。

5.2 ? 多元線性回歸模型的建立

回歸分析用來評估變量之間關系的統(tǒng)計過程，解釋自變量X與因變量Y的關系，這里采用一般多元線性回歸模型[12]。將正交實驗的數(shù)據(jù)通過MATLAB進行多元線性回歸分析，得到軸向力預測模型為：

其中，F(xiàn)為軸向力，D為鉆頭直徑，f為鉆削進給速度，n為主軸轉速。表7為該預測數(shù)學模型的模型統(tǒng)計量表，其中包括R2統(tǒng)計量、F統(tǒng)計量及其P值，以及誤差方差的估計值。R2值為0.97，表明該模型可以擬合97%的實驗數(shù)據(jù)，且P值小于0.05，擬合效果顯著。

5.3 ? 結果分析

圖8為鉆削力與鉆頭直徑、鉆削進給速度以及鉆頭轉速之間的回歸曲面圖，從圖8（a）、圖8（b）、圖8（c）中可以看出，對鉆削力影響最大的參數(shù)是鉆頭直徑，其次是進給速度，最后是轉速，并且直徑與進給速度和轉速之間的交互影響同時影響著鉆削力的大小。這是因為鉆削力主要是由鉆頭的橫刃與切削刃產生的，鉆頭直徑的增大直接影響了橫刃與切削刃的大小，從而導致鉆削力的增大。而相對于轉速，進給速度的快慢在非高轉速的情況下對鉆削力的影響更大。