徐 新,鄧 斌,劉寧波
(1.西南交通大學先進驅動節能技術教育部工程研究中心,四川 成都 610031;2.海軍航空大學信息融合研究所,山東 煙臺 264000)
如今,艾灸[1]已普及到百姓人家,遍布大街小巷。人人均可在理療店甚至在家進行艾灸。但是在艾灸裝備領域始終存在自動化程度不高、艾灸溫度控制不精確等問題。因此人們目前更傾向于人工艾灸,可以自己根據皮膚的溫度承受度手動調節艾灸的距離和時間提高艾灸的效果及舒適度。人工艾灸雖效果好,但無法解放雙手,艾灸穴位多時間長,長時間的手持艾條,耗時耗力。故適應快速增長的生活高質量的意愿,那么迫切的需要提高艾灸裝置的適用性和實用性。
目前,艾灸裝置發展有兩個方向:第一個方向是向著大型化、集成化、全身各部位同時艾灸的方向發展,文獻[2]研制的智能灸療床,還存在諸多問題,一般在專業人員的指導下使用,不能普及;第二個方向是向著小型化、便攜化、單穴位艾灸方向發展。這類艾灸裝置多種多樣,文獻[3]設計的電子艾灸盒,通過滾珠絲杠帶動點燃的艾條進行艾灸,但是采用的是DS18B20溫度傳感器,只能采集到艾灸時的環境溫度,溫度控制不精準。文獻[4]設計了一種新型智能艾灸儀,分析了艾灸的溫度場,但通過送風來控制溫度,不可控。目前艾灸裝置普遍存在不能自動清灰、艾灸溫度控制不精確的問題。所以設計一種新的艾灸裝置具有一定的時代意義和現實意義。
為了使用方便和成本低減,設計的艾灸裝置在保證各方面功能的同時需做到輕量化、保證力學性能不受影響。這就需要對艾灸裝置中對質量影響最大的部件艾灸支撐座進行優化設計。傳統的優化設計是將各個參數變成設計變量進行優化,但是將每個參數都進行考慮則計算工作量大且效率低下,文獻[5-6]通過改變結構進行優化,但沒有系統的優化設計方法。文獻[7]將影響立式加工中心的幾個關鍵尺寸作為設計變量,對其進行了輕量化優化設計,得到了較好的效果。文獻[8]在曲軸結構優化中利用響應面法結合多島遺傳算法得到Pareto 解集,具有很好的工程意義。
針對上述情況,這里欲設計一種既能提高艾灸效率又能降低自身質量的艾灸裝置,并以關鍵部件艾灸支撐座為對象進行多目標的優化。擬定設計變量為艾灸支撐座的三個結構尺寸,利用OSF 方法作為設計變量的試驗點的選取方法,再構建神經網絡代理模型,最后檢驗該模型的精度。以艾灸支撐座的質量和最大變形量為多目標優化函數,以最大變形量和最大應力為約束條件,利用多島遺傳算法選取最優解,為其他結構的優化設計提供了參考。
設計艾灸裝置的目的是提高艾灸的自動化程度和效率,因此對艾灸裝置的設計要求有:
(1)覆蓋人體各個穴位;
(2)實現艾條自動推送;
(3)具有清灰功能;
(4)盡可能質量輕、結構小巧。
基于以上要求設計出的艾灸裝置,如圖1所示。艾灸裝置的工作原理是通過金屬定型軟管和旋轉關節將艾條豎直對準人體將要艾灸的區域,通過推送部件調整艾條與艾灸點皮膚之間的距離,通過模糊PID控制算法改變艾條與艾灸點皮膚之間的距離實現艾灸點皮膚溫度的恒溫控制,并設計專門的清灰裝置在需要清理艾灰時清理艾灰,保證艾灸效果。整個艾灸裝置具有結構輕巧、能自動清灰、精準控制艾灸溫度等優點。
圖1 艾灸裝置結構示意圖Fig.1 Schematic Diagram of Moxibustion Device Structure
艾灸裝置包括艾灸裝置本體、清灰裝置以及集灰裝置。清灰裝置和集灰裝置從上到下依次設置在艾灸支撐座上,艾灸支撐座上還設有推送部件,艾灸裝置本體、清灰裝置、集灰裝置以及推送部件分別與控制裝置電連接。
艾灸裝置本體包括:底座、金屬定型軟管、旋轉關節、艾灸支撐座以及艾條推送組件。艾條推送組件的另一端設有艾條。旋轉關節為舵機驅動,通過舵機安裝座和軸承的連接實現艾灸支撐座的旋轉。推送部件是由步進電機驅動,通過同步帶傳動帶動滾珠絲杠運動,艾條通過“Z”字形連桿固定在絲杠螺母上實現直線運動。艾灸裝置本體為艾灸的主體功能,通過金屬定型軟管的粗調整和旋轉關節精確調整將艾條豎直對準人體將要艾灸的區域,工作范圍可覆蓋人體各個穴位,可滿足設計要求(1)。推送部件實現艾灸的推送,可滿足設計要求(2)。
清灰裝置結構,如圖2所示。
圖2 清灰裝置結構示意圖Fig.2 Schematic Diagram of the Structure of the Cleaning Device
清灰裝置用于清理艾灸燃燒端的艾草灰,因此該清灰裝置整體設置在外殼下端,艾條在推送部件的帶動下可進入清灰裝置中進行清灰,在旋轉電機的驅動下旋轉殼不斷轉動,刷體組件在旋轉的同時還可以在導向盤的波浪形導向面的導向作用不斷上下移動,縱向導向組件帶動刷頭上下移動,橫向導向組件帶動刷頭橫向移動,從而對艾條完成多方位清灰,在旋轉殼中沿艾條的周向設置多組清灰刷體組件,每組清灰刷體組件采用不同形狀的刷頭,加強清灰效果;集灰裝置用于收集清灰裝置清理下來的艾草灰,集灰裝置中的滑蓋可以在清灰時將集灰裝置外殼上供艾條穿出的通孔封閉,避免清灰時,艾草灰從通孔掉落到人體皮膚上,集灰裝置內部設有溫度傳感器,當溫度傳感器檢測到的溫度達到設定值時,可判斷艾條需要清灰,因此觸發控制裝置進行清灰操作。
清灰功能只有一個電機驅動,通過各個機械結構的設計,完成艾灸條的清灰操作,滿足功能(3)的需求,解決了以往艾灸時需要手動清灰的難題,并且提高了艾灸時的溫度控制精度。
艾灸裝置裝有紅外溫度傳感器、編碼器來控制艾灸裝置的正常工作。紅外傳感器1負責檢測艾灸點皮膚表面溫度,紅外傳感器2負責檢測艾灸燃燒端溫度。設定艾灸時的溫度,利用傳感器1作為反饋,滾珠絲杠為執行機構,在模糊PID算法的控制下調整艾條與艾灸點皮膚之間的距離,達到恒溫控制效果。當傳感器2檢測到艾條燃燒端溫度高于設定值,起動清灰功能完成艾灰的清理。因此,通過機械硬件和控制算法的結合,實現了艾灸時的高精度溫度控制。
艾灸支撐座為艾灸推送組件和艾灸清灰集灰組件提供支撐。由于艾灸支撐座為艾灸裝置主要的受力部件以及質量最大的部件,因此以艾灸支撐座為研究對象,進行有限元仿真分析,在保證力學性能的基礎上實現輕量化以滿足上述功能(4)的需求。
利用三維建模軟件建立艾灸支撐座實體模型,將凹槽、倒角等結構忽略以方便計算,將模型的各個結構尺寸進行參數化,將具有參數化結構尺寸的模型導入有限元軟件進行分析,艾灸支撐座Workbench模型,如圖3所示。
圖3 艾灸支撐座Workbench模型Fig.3 Workbench Model of Moxibustion Support Seat
艾灸支撐座以Nylon610為材料,彈性模量為8300MPa,泊松比為0.28,密度為1400kg/m3,屈服強度為68.9MPa。此時,艾灸支撐座的質量0.13643kg。選擇通用性較強的四面體網格,網格最小單元尺寸1mm,單元數為1135572個,節點數為706669。
根據使用時極端工況,艾灸支撐座主要有幾個力與約束如下:
(1)艾灸清灰上外殼通過螺栓連接在艾灸支撐座的正面,受艾灸清灰組件重力;
(2)步進電機安裝座與艾灸支撐座通過螺栓連接,受步進電機及其安裝座重力;
(3)艾條推送組件重力;
(4)艾灸清灰時產生的彎矩。
將以上力分別加在艾灸支撐座上,圓柱端固定,對其進行靜力學分析,應力圖和總變形圖,如圖4、圖5所示??梢?,艾灸支撐座在應力出現在背面支撐點的附近,但應力較小,小于標準的10%,滿足設計及使用要求;最大變形出現在前端的底部,小于標準值,且考慮實際使用情況,不會對艾灸裝置的工作產生較大影響,也滿足要求。
圖4 艾灸支撐座應力圖Fig.4 Stress Diagram of Moxibustion Support Seat
圖5 艾灸支撐座總變形圖Fig.5 The Total Deformation of the Moxibustion Support Seat
該部件的應力應變結果,如表1所示。
表1 艾灸支撐座的最大變形和最大應力Tab.1 Maximum Deformation and Maximum Stress of Moxibustion Support
為了進一步使艾灸裝置使用方便和成本低減,需要對影響艾灸裝置整體結構性能最大的艾灸支撐座做進一步分析。
響應曲面法(Response Surface Methodology)[9-10]是一種結合數學方法和統計方法的優化方法。為了使艾灸裝置輕量化且滿足使用需求,在不改變艾灸支撐座的基本結構的前提下,以最大程度減輕質量和減小最大變形量為目標,用響應面法對艾灸支撐座進行優化設計。
對于優化設計來說,一般響應與設計變量之間的關系式是未知的,因此利用響應曲面法得到它們之間的真實函數關系。由于艾灸支撐座三個變量之間對響應不是線性關系,因此使用二次多項式建立擬合函數:
式中:y—目標函數;
xi,xj—設計變量;
β,ε,k—常量。
試驗點的選取對響應面的精度至關重要。本次設計采取Optimal Space-Filling Design(簡稱OSF)的試驗設計方法選取試驗點。OSF方法是根據拉丁超立方采樣算法進一步優化得到的,有著更為均勻的點空間分布。
相對于一般的CCD(Central Composite Design)算法,OSF 可以以最小的數量獲得對設計點的最大洞察,相應的需要更多設計點以免影響響應預測質量。對艾灸支撐座影響較大的3個參數選為設計變量:長方形板厚記錄為d1,初始值為2mm;長方形寬度記錄為d2,初始值為25mm;上方減重孔寬度記錄為d3,初始值為30mm。根據設計經驗和實際能利用的空間確定了設計變量的取值范圍,如表2 所示。利用OSF 方法計算得到艾灸支撐座設計試驗取值,如表3所示。
表2 設計變量的取值范圍(mm)Tab.2 Value Range of Design Variables(mm)
表3 設計變量的設計試驗取值(mm)Tab.3 Design Test Values of Design Variables(mm)
以表2所示的艾灸支撐座d1、d2、d3三個參數為影響因素,對艾灸支撐座進行靈敏度分析,如圖6所示。
圖6 艾灸支撐座靈敏度分析Fig.6 Sensitivity Analysis of Moxibustion Support Base
由圖6可知,三個設計變量d1、d2、d3對艾灸支撐座最大變形量、最大位移量、質量的靈敏度分別為10、-58、-6,7、-68、5,22、75、-3。說明三個設計變量對艾灸支撐座的最大變形量、最大位移量、質量均有影響,這三個設計變量的選取具有實際意義。其中,設計變量d2的靈敏度值最大,說明對艾灸支撐座的影響最大。
徑向基函數是形如(zx)=Υ(‖x-u‖)函數,其中,‖x-u‖是向量x與中心μ之間的歐氏距離。徑向基神經網絡[11]具有結構簡潔,很好的解決非線性問題的能力,并且學習速度快等優點。因此,構建了艾灸支撐座結構尺寸與靜態特性的徑向基神經網絡代理模型,如圖7所示。
圖7 徑向基神經網絡模型Fig.7 Radial Basis Function Neural Network Model
輸出層第k個節點的輸出的計算方法如下:
式中:wki—隱含層第i個節點到輸出層第k個節點的連接權重;
wk0—輸出層第k個節點的偏差項。
徑向基神經網絡的誤差函數為:
式中:yd,k—節點k的期望輸出。
徑向基神經網絡的權值參數可用最小二乘法求解,可將式(2)寫成矩陣:
通過式(2)~式(4)即可求得艾灸支撐座質量、最大形變量,最大應力關于結構尺寸的徑向基神經網絡代理模型。
對于神經網絡響應面精度的好壞,需要進行檢驗。一般采用RMSE(均方根差),R(2復相關系數)(修正的復相關系數)來檢驗。其公式如下:
式中:m—模型驗證的樣本量;
yi—真實響應值;
—由回歸方程得到的觀測值;
—真實響應值的均值。
三個評價標準,如表4所示??傻庙憫娴木葷M足要求。
表4 評價標準擬合度值Tab.4 Evaluation Standard Fitness Value
優化設計的目標是以艾灸支撐座的質量和最大變形量為多目標優化函數,約束條件為最大變形量和最大應力,設計變量為艾灸支撐座的三個結構尺寸,因此構建如下多目標優化的數學模型:
式中:X—決策矢量,X∈Ω,Ω—決策空間;
dmin—設計尺寸變量下限值;
dmax—設計尺寸變量上限值;
m—整體質量;
λ—允許變形量;
λ(X)—最大變形量;
[σ]—許用應力;
σ(X)—最大應力。
在多目標優化過程中無法實現多個目標同時最優,所以通過多目標優化計算求解出多目標優化解集,即Pareto最優解集[12],再通過要求和經驗選取最優解。有很多種方法可以選取Pareto最優解,這里采用具有多島效應的遺傳算法對進行優化,該算法可以同時處理多個目標,魯棒性強。
各個子目標的權重都設為1,初始種群數量設為50,島嶼個數10,最大遺傳代數20,交叉概率1.0,變異概率0.01,遷徙率0.5。Pareto解集在經過2000次迭代后得到的結果,如圖8所示。得到的最終優化結果,如表5、表6所示。結果表明,利用多島遺傳算法對徑向基神經網絡代理模型進行優化,質量、最大變形量、最大應力的誤差分別為0.026%、5.53%、10.30%,優化精度較為理想。與初始值相比,質量從0.13643kg 下降到0.11595kg,下降了15.01%。最大變形量從2.3395mm 變為2.3613mm 小于3mm,僅僅增加了0.0218mm,對整個艾灸裝置的運行沒有影響。最大應力從5.9446MPa 上升到6.2052MPa,由于遠小于68.9MPa。優化結果滿足要求。
表5 設計變量優化前、后和工程化取值(mm)Tab.5 Design Variables Before,After and Engineering Value of Optimization(mm)
表6 優化前、后結果對比Tab.6 Comparison of Results Before and After Optimization
圖8 Pareto最優解集Fig.8 Pareto Optimal Solution Set
(1)設計了一種結構輕巧、能自動清灰且能精準控制艾灸溫度的艾灸裝置。具有實用性,能夠有效代替人工手動艾灸,省時省力。
(2)對艾灸支撐座的設計優化和分析表明,基于徑向基神經網絡的艾灸支撐座響應面模型和利用多島遺傳算法的多目標優化設計,質量從0.13643kg 下降到0.11595kg,下降了15.01%,有效降低了支撐座的質量,實現了艾灸裝置的輕量化目標,同時保證了力學性能。
(3)該艾灸裝置結構的優化成功,為其他艾灸設備的機械結構設計優化提供了理論依據和設計參考。