許李悅,張文光,賀雨欣,周旭暉
(上海交通大學(xué)機(jī)械系統(tǒng)與振動(dòng)國(guó)家重點(diǎn)實(shí)驗(yàn)室,200240,上海)
腦部植入式神經(jīng)電極是對(duì)大腦進(jìn)行研究和神經(jīng)疾病治療的重要手段[1]。傳統(tǒng)的腦部植入式神經(jīng)電極多采用硅等剛性材料制成[2]。由于大腦存在微動(dòng),剛性材料與大腦之間的力學(xué)失配會(huì)導(dǎo)致植入電極周圍腦組織炎癥反應(yīng)增加、神經(jīng)元密度下降,并產(chǎn)生膠質(zhì)瘢痕對(duì)電極形成包裹,使電極的電信號(hào)記錄能力下降[3-4]。
使用柔性材料替代剛性材料作為電極基底已被證明可以減少這些不良反應(yīng)[5],但由于柔性電極的基底材料楊氏模量較小,且形狀細(xì)長(zhǎng),極易在植入大腦的過(guò)程中發(fā)生屈曲[6],因此一般都采用輔助手段臨時(shí)增加電極的強(qiáng)度來(lái)實(shí)現(xiàn)柔性電極的植入。一種常用的方法是采用剛性的輔助工具暫時(shí)與電極相連,在植入后再將工具回撤[7-8],但這種方法容易導(dǎo)致電極植入后定位的不準(zhǔn)確[9]。
另一種方法是使用生物可溶性材料如聚乙二醇(PEG)[10]、絲蛋白[11]、葡聚糖[12]等作為涂層,臨時(shí)增強(qiáng)柔性電極的剛度,植入到指定位置后,生物可溶性涂層能在一定時(shí)間內(nèi)在腦組織和體液中溶解并降解,使電極恢復(fù)柔性[13]。這種方法不存在回撤的問(wèn)題,且制作加工較為方便靈活。常用的涂層制作方法有浸涂法[14]和模具注塑法[11],其中浸涂法無(wú)法較好地控制涂層尺寸,且表面精度較差,結(jié)構(gòu)粗糙[15],容易對(duì)腦組織產(chǎn)生更大的損傷。本文采用改進(jìn)的模具注塑法制造涂層,可以較好地控制涂層的形狀、尺寸及表面精度。
生物可溶性涂層的截面積大小決定了植入時(shí)對(duì)腦組織的損傷范圍[16],是電極生物相容性的重要評(píng)判指標(biāo)。目前尚未有對(duì)涂層的截面形狀和尺寸進(jìn)行優(yōu)化設(shè)計(jì)的研究。對(duì)此,本文提出一種T形截面涂層,優(yōu)化了涂層的截面積,并基于植入力、加工工藝等約束,通過(guò)遺傳算法設(shè)計(jì)了適用于不同寬度電極的最佳尺寸涂層。
目前對(duì)于使用涂層輔助柔性神經(jīng)電極植入的植入力暫未有系統(tǒng)的研究,已有研究中測(cè)得腦部植入式神經(jīng)電極的植入力范圍較大,一般在幾到幾百毫牛[11-12,17],無(wú)法為涂層的優(yōu)化設(shè)計(jì)提供有力的力學(xué)支撐。本文針對(duì)不同尺寸的涂層進(jìn)行了植入力的實(shí)驗(yàn)測(cè)量,并分析了植入力與涂層的尺寸關(guān)系,為T(mén)形截面涂層的尺寸參數(shù)設(shè)計(jì)提供了參考。
植入力是指神經(jīng)電極植入腦組織時(shí)所受到的力?,F(xiàn)有研究表明腦部植入式神經(jīng)電極在植入過(guò)程中所受的最大植入力與植入物的形狀、尺寸、材料和植入速度有關(guān)[16,18-19]。但目前尚無(wú)針對(duì)涂層輔助柔性電極植入的植入力情況的研究,本文主要針對(duì)涂層的尺寸對(duì)植入力的影響進(jìn)行了研究,具體方法為控制矩形截面涂層的長(zhǎng)度和厚度不變,僅改變涂層寬度來(lái)得到不同尺寸的涂層,同時(shí)涂層材料、植入速度和深度也保持不變。
植入實(shí)驗(yàn)所使用的平臺(tái)如圖1所示,包括數(shù)碼顯微鏡、NMB力傳感器(UT-100GR,量程980.7 mN)、數(shù)據(jù)采集系統(tǒng)(PCD-30A)、微驅(qū)動(dòng)器、Z軸升降平臺(tái)。為簡(jiǎn)化實(shí)驗(yàn)并使實(shí)驗(yàn)具有較好的重復(fù)性,植入對(duì)象為0.6%瓊脂糖凝膠制作的模擬腦組織,其力學(xué)性能已被證明與腦組織一致[18]。
圖1 電極植入力測(cè)量系統(tǒng)
植入時(shí)將帶涂層的電極固定在微驅(qū)動(dòng)器上,力傳感器與微驅(qū)動(dòng)器相連,所測(cè)得的電極受力即為植入力??刂泼看沃踩雽?shí)驗(yàn)的植入深度為2 mm,植入速度為1 mm/s,力傳感器的數(shù)據(jù)采樣頻率設(shè)為500 Hz。數(shù)碼顯微鏡可觀察電極與模擬腦組織的對(duì)準(zhǔn)情況以及電極在植入過(guò)程中的形態(tài)。模擬腦組織固定于Z軸升降平臺(tái)上,可調(diào)整與電極的相對(duì)高度,改變植入點(diǎn)的位置。
圖2給出了柔性電極與涂層示意圖。實(shí)驗(yàn)所用柔性電極尺寸如圖2a所示,長(zhǎng)5 mm,寬300 μm,厚15 μm,材料為光敏聚酰亞胺(PSPI),通過(guò)光刻工藝加工,為簡(jiǎn)化實(shí)驗(yàn),電極不具有電學(xué)功能,僅用于驗(yàn)證其機(jī)械性能。
(a)柔性電極尺寸 (b)電極與涂層組裝示意圖
涂層材料為生物可溶性材料聚乙二醇(PEG,相對(duì)分子質(zhì)量20 000)和抗炎藥物地塞米松(DEX)的混合物,其中PEG與DEX配制的質(zhì)量比例為7∶1,混合物的彈性模量為(58.07±1.13) MPa。通過(guò)微模具注塑法在PSPI電極的一側(cè)制作涂層,組裝后的電極如圖2b所示。
常見(jiàn)輔助柔性電極植入的涂層寬度一般在350~600 μm[11-12,16],因此,實(shí)驗(yàn)中使用的涂層的寬度分別設(shè)計(jì)為350、400、450、500、550、600 μm,厚度均為400 μm。對(duì)每個(gè)寬度的帶涂層電極進(jìn)行12次植入實(shí)驗(yàn),去掉未成功植入的數(shù)據(jù)。
實(shí)驗(yàn)測(cè)得每次植入力均有如圖3所示趨勢(shì):從植入開(kāi)始,植入力一直呈上升趨勢(shì),在電極到達(dá)目標(biāo)位置后停止植入,此時(shí)植入力達(dá)到最大值,之后下降,并不降為0,說(shuō)明植入完成后,腦組織仍持續(xù)受到電極的植入壓力作用。
圖3 電極植入力曲線
將每次植入實(shí)驗(yàn)中所測(cè)得的最大植入力作為評(píng)判植入力的指標(biāo),也可作為后續(xù)T形截面涂層的設(shè)計(jì)依據(jù)。記錄下每次實(shí)驗(yàn)中植入力的最大值,進(jìn)行后續(xù)的實(shí)驗(yàn)數(shù)據(jù)處理。
圖4 不同涂層寬度的最大植入力
由圖4可以看出,同一涂層寬度下所測(cè)得的最大植入力的分布較廣,具體體現(xiàn)為每個(gè)寬度下植入力的標(biāo)準(zhǔn)差較大。因此,為研究涂層寬度確實(shí)對(duì)植入力造成了影響,而不是環(huán)境等其他因素造成的不同涂層寬度植入力數(shù)據(jù)的差異,對(duì)所測(cè)最大植入力數(shù)據(jù)進(jìn)行單因素方差分析。
假設(shè)除涂層寬度外的其他因素,如環(huán)境和人為操作等對(duì)每一寬度下的植入力測(cè)量的影響相同,則可將每一寬度下測(cè)得的最大植入力分別視為獨(dú)立服從同方差的正態(tài)分布總體,即Fi~N(μi,σ2),其中μ為總體均值,σ2為總體方差。
假設(shè)涂層寬度對(duì)最大植入力不存在影響,即
H0:μ1=μ2=μ3=μ4=μ5=μ6
(1)
對(duì)不同寬度下測(cè)得的最大植入力做單因素方差分析,計(jì)算得檢驗(yàn)統(tǒng)計(jì)量F=3.8,顯著性水平在0.01的情況下,可拒絕原假設(shè),說(shuō)明涂層寬度對(duì)最大植入力存在影響。
F=0.014 6W+0.557 3
(2)
再利用F檢驗(yàn)判斷線性關(guān)系是否顯著,取顯著性水平為0.01時(shí),F0.01(1,8)=11.26,計(jì)算得檢驗(yàn)統(tǒng)計(jì)量F=18.95>F0.01(1,8),同時(shí)相關(guān)系數(shù)r=0.908 7,由此可以認(rèn)為植入過(guò)程中的最大植入力與涂層寬度具有顯著的線性關(guān)系,表明對(duì)于涂層輔助植入的電極,存在涂層尺寸越大,植入力越大的關(guān)系。
對(duì)于每一寬度下最大植入力的標(biāo)準(zhǔn)差較大的原因進(jìn)行分析,從實(shí)驗(yàn)中來(lái)看,主要是由于瓊脂糖凝膠模擬腦組織不均勻所造成的,即多次在同一瓊脂糖凝膠中植入后,瓊脂糖會(huì)發(fā)生變形、表面產(chǎn)生褶皺等現(xiàn)象,而實(shí)際腦組織也是不均勻的,因此該現(xiàn)象可以接受。腦組織不均勻的因素很難在實(shí)際應(yīng)用中徹底排除,因此后續(xù)設(shè)計(jì)將以所有實(shí)驗(yàn)中測(cè)得的最大植入力作為涂層臨界屈曲力的設(shè)計(jì)指標(biāo),對(duì)于厚度為400 μm,寬度在350~600 μm的矩形涂層,實(shí)驗(yàn)中測(cè)得的F的最大值Fmax為14.2 mN。
基于上述實(shí)驗(yàn)所測(cè)得的植入力的約束,進(jìn)行涂層截面的優(yōu)化設(shè)計(jì)。由于縮小植入物的截面積能減少植入時(shí)對(duì)腦組織的損傷[16],從而減少腦組織的炎癥反應(yīng)及對(duì)電極的膠質(zhì)瘢痕包裹,因此從電極的生物相容性和長(zhǎng)期使用的角度考慮,本研究所設(shè)計(jì)的涂層采用一種T形截面結(jié)構(gòu),并以最小化涂層截面積為設(shè)計(jì)目標(biāo)。
涂層輔助植入的電極在植入過(guò)程中可簡(jiǎn)化為一端固定、一端鉸支的細(xì)長(zhǎng)壓桿:其尖端與腦組織接觸,不可移動(dòng)但可發(fā)生轉(zhuǎn)動(dòng),表現(xiàn)為鉸支;其末端與植入平臺(tái)固定連接,并受到徑向力的作用。
輔助柔性電極植入涂層的一個(gè)重要設(shè)計(jì)指標(biāo)是其臨界屈曲力,這是其能否成功植入腦組織的關(guān)鍵。若電極發(fā)生屈曲,即產(chǎn)生不可逆的變形,則無(wú)法繼續(xù)植入,且電極也可能產(chǎn)生損傷甚至破裂,因此需使涂層輔助植入的電極整體結(jié)構(gòu)的臨界屈曲力大于植入中可能受到的最大植入力。由于電極本體的體積在整體結(jié)構(gòu)中所占的比例很小,因此在后續(xù)計(jì)算中忽略電極,直接計(jì)算涂層的屈曲力情況。理論計(jì)算所得到的結(jié)果為線性屈曲,實(shí)際工程中一般取3~5倍的安全系數(shù)[20]。
由此可得涂層的臨界屈曲力Fb的力學(xué)約束為
Fb≥kFmax
(3)
式中:k為安全系數(shù),為保證結(jié)構(gòu)的可靠性,取為5。
涂層的屈曲力Fb可通過(guò)歐拉公式計(jì)算
(4)
式中:E為涂層材料的彈性模量;I為橫截面積的慣性矩;L為涂層長(zhǎng)度;u為有效長(zhǎng)度系數(shù),在該種受力情況下取0.7。
本文所提出的T形截面涂層結(jié)構(gòu)如圖5所示。涂層長(zhǎng)度L為5 mm,尖端角度為60°,尖端有R為0.1 mm的圓角以方便模具的加工和涂層的制造。橫截面的形狀為T(mén)字形,待優(yōu)化的設(shè)計(jì)參數(shù)W、H、w、h如圖5所示。
圖5 T形截面涂層結(jié)構(gòu)示意圖
涂層的設(shè)計(jì)目標(biāo)是使其截面積在滿足約束的條件下達(dá)到最小,以減少對(duì)大腦的損傷,對(duì)此優(yōu)化模型的目標(biāo)函數(shù)可用涂層截面積表示
S=WH+wh
(5)
基于原厚度為400 μm的矩形涂層進(jìn)行截面的優(yōu)化設(shè)計(jì),限制涂層的總厚度H+h不超過(guò)400 μm,涂層寬度W不超過(guò)600 μm,以便使用1.1節(jié)實(shí)驗(yàn)中測(cè)得的植入力對(duì)涂層力學(xué)性能進(jìn)行約束。涂層的力學(xué)約束條件為:其臨界屈曲力應(yīng)大于植入所需的最大植入力(5倍安全系數(shù)情況下)。
涂層的寬度還由電極的寬度約束,為了便于組裝,要求涂層寬度W至少比電極寬100 μm,即W≥B+100 μm(B為PSPI電極寬度),并且保證T形結(jié)構(gòu)中w小于W。
涂層模具的加工采用精密數(shù)控銑的方法,其要求寬度最小尺寸大于200 μm,厚度最小尺寸大于100 μm,因此有如下加工工藝約束:w≥200 μm,H≥100 μm,h≥100 μm。
分別對(duì)不考慮加工工藝約束和考慮加工工藝約束這兩種情況進(jìn)行涂層的尺寸優(yōu)化設(shè)計(jì)??傮w設(shè)計(jì)過(guò)程如圖6所示。
圖6 T形截面涂層設(shè)計(jì)流程圖
由此抽象出優(yōu)化模型的方程
(6)
由于該優(yōu)化問(wèn)題為非凸規(guī)劃,用Matlab遺傳算法進(jìn)行全局最優(yōu)解的近似求解。初始種群隨機(jī)產(chǎn)生,種群大小設(shè)為100,最大迭代次數(shù)為300次,最大變異概率為0.08,交叉概率為0.5,得到優(yōu)化結(jié)果(取近似值)。此參數(shù)條件下,得到的結(jié)果具有較好的收斂性,遺傳算法優(yōu)化T形截面涂層設(shè)計(jì)參數(shù)的收斂曲線如圖7所示,某次遺傳算法結(jié)果的極值和平均值在250代時(shí)已收斂。
圖7 遺傳算法優(yōu)化T形截面涂層設(shè)計(jì)參數(shù)的收斂曲線
涂層寬度W在400~600 μm范圍內(nèi),分別通過(guò)遺傳算法獲得無(wú)加工工藝約束和有加工工藝約束下的最優(yōu)w、H、h值,將有加工工藝約束時(shí)的最優(yōu)設(shè)計(jì)參數(shù)見(jiàn)表1,并計(jì)算涂層橫截面積S,不同寬度下T形截面涂層最優(yōu)截面積曲線如圖8所示,由于寬度為400 μm時(shí),無(wú)法得到滿足約束結(jié)果的解,因此以下均只討論寬度大于等于425 μm的情況。
表1 不同寬度T形截面涂層最優(yōu)設(shè)計(jì)參數(shù)
圖8 不同寬度下T形截面涂層最優(yōu)截面積曲線
由圖8可得,涂層寬度在425~600 μm的范圍內(nèi),在無(wú)加工工藝約束的情況下,理論上可得當(dāng)涂層寬度W越大時(shí),橫截面積S越小,因此當(dāng)電極寬度B一定時(shí),設(shè)計(jì)涂層時(shí)可以在合理的情況下盡可能地取最大的涂層寬度。
但在本文所選擇的加工工藝條件下,涂層橫截面積在寬度為525 μm時(shí)取得最小,寬度越遠(yuǎn)離525 μm,橫截面積越大。由此可得出當(dāng)電極寬度B小于425 μm時(shí),都可以采用寬度為525 μm的T形截面涂層,當(dāng)B大于425 μm時(shí),則應(yīng)當(dāng)選擇適合電極寬度的最小涂層寬度。
由圖8可以看出,加工工藝的約束使各寬度下最優(yōu)橫截面積相比無(wú)加工工藝約束下均有增加,平均增加幅度為12.7%??梢赃x擇用精度更高的加工工藝如光刻等方法來(lái)進(jìn)行涂層模具的加工,但這也會(huì)增加加工成本。在設(shè)計(jì)時(shí)需綜合考慮成本和實(shí)際使用要求。
針對(duì)加工工藝約束條件下的T形截面涂層與原矩形截面涂層進(jìn)行性能比較。在相同寬度下,T形截面涂層的橫截面積均有減少,平均面積僅為原矩形涂層的58.24%,表明T形截面形狀對(duì)涂層橫截面積有非常好的優(yōu)化效果。
相比于原矩形截面涂層,若要達(dá)到與T形截面涂層相同的屈曲力,均需要設(shè)計(jì)更大的橫截面積,T形截面涂層與矩形涂層橫截面積對(duì)比如圖9所示,T形截面涂層的平均截面積僅為同屈曲力的矩形涂層的70.32%??梢?jiàn)T形截面涂層相比于傳統(tǒng)的矩形涂層具有很大的優(yōu)越性。
圖9 相同屈曲力T形截面涂層與矩形涂層橫截面積對(duì)比
針對(duì)寬度為300 μm的PSPI柔性電極,選取寬度為525 μm的最優(yōu)尺寸T形截面涂層作為輔助柔性電極植入的手段。通過(guò)Ansys有限元仿真涂層植入時(shí)的受力情況,并加工制造出涂層進(jìn)行實(shí)物的植入實(shí)驗(yàn)。
在實(shí)物實(shí)驗(yàn)前先對(duì)所設(shè)計(jì)的涂層進(jìn)行有限元植入力學(xué)仿真。選擇表1中W為525 μm時(shí)的設(shè)計(jì)參數(shù)進(jìn)行建模(僅考慮涂層植入部分),涂層末端為固定約束,涂層尖端為鉸支。在Ansys Workbench軟件中進(jìn)行靜力學(xué)仿真,網(wǎng)格劃分采用四面體單元,網(wǎng)格尺寸為0.08 mm。1.1節(jié)中植入實(shí)驗(yàn)測(cè)得涂層寬度為525 μm左右時(shí)的植入力范圍在3.55~12.26 mN之間,在此范圍內(nèi)取多個(gè)植入力進(jìn)行仿真,得到不同植入力下的變形量曲線如圖10所示,最大變形量在2~12 μm之間,變形量很小,且均在彈性限度范圍內(nèi),可以保證電極的順利植入。
圖10 不同植入力下涂層變形量
再聯(lián)合靜力學(xué)和線性屈曲模塊進(jìn)行屈曲仿真,得涂層的臨界屈曲力為72.713 mN,符合設(shè)計(jì)要求,發(fā)生屈曲時(shí)的最大位移為1.58 mm,T形截面涂層的屈曲仿真變形云圖如圖11所示。
圖11 T形截面涂層的屈曲仿真變形云圖
為使涂層便于制造,模具的制造分為剛性的凸模和柔性的凹模。凸模材料為有機(jī)玻璃(PMMA),采用精密數(shù)控銑加工。凹模材料為聚二甲基硅氧烷(PDMS),將PDMS與固化劑(SYLGARD-184B)以10∶1的比例混合均勻后沉積在PMMA凸模上,在80 ℃的熱板上加熱固化4 h,固化完成后用鑷子分離取下,即可得到柔性凹模,如圖12所示。PDMS制作的凹模柔軟可變形,方便將制作好的涂層從凹模中完整地取出。
圖12 PDMS凹模
涂層材料依舊采用PEG與DEX的混合物。本文采用一種適合于T形截面涂層的改進(jìn)模具注塑方式,具體流程如圖13所示:在PSPI柔性電極上注塑涂層時(shí),先將PDMS模具放置在90 ℃的熱板上,將混合涂層材料填入凹模內(nèi),待完全熔化后,將鑷子深入凹模通道內(nèi),將多余涂層材料刮除;從熱板上取下模具,待涂層材料冷卻凝固后,用鑷子將PSPI電極水平放置在凹模通道中央;再將模具放回90 ℃熱板上加熱,并再次用涂層材料填滿模具,待涂層材料熔化后用鑷子刮平,冷卻后將涂層和電極整體取出即可。此方法制作的電極兩側(cè)均被涂層包裹覆蓋,在植入過(guò)程中不易產(chǎn)生電極與涂層的分層現(xiàn)象。
圖13 T形截面涂層制作過(guò)程
制作完涂層的電極如圖13所示,涂層結(jié)構(gòu)表面光滑,與電極連接緊密。制造10個(gè)如圖14所示的帶涂層電極,在圖1所示的實(shí)驗(yàn)平臺(tái)上將電極植入瓊脂糖凝膠模擬腦組織,植入實(shí)驗(yàn)參數(shù)依舊控制為植入深度2 mm,植入速度1 mm/s,在力傳感器的數(shù)據(jù)采樣頻率為500 Hz的情況下,獲得實(shí)驗(yàn)中典型的植入力曲線,與矩形截面涂層的植入力曲線對(duì)比如圖15所示。實(shí)驗(yàn)結(jié)果是10個(gè)電極均成功植入,測(cè)得植入力為(4.84±1.04) mN,植入力較未優(yōu)化的矩形涂層大大減小,僅為原矩形涂層平均植入力的64.62%。
(a)T形截面涂層外觀
圖15 兩種涂層的植入力曲線對(duì)比
本文針對(duì)輔助柔性神經(jīng)電極植入的涂層進(jìn)行了植入力的研究,并提出一種T形截面的新型涂層,對(duì)涂層的截面積進(jìn)行了優(yōu)化,大幅減少了植入過(guò)程中對(duì)大腦造成的損傷,主要結(jié)論如下。
(1)通過(guò)實(shí)驗(yàn)測(cè)量不同寬度的帶涂層電極的植入力,得出涂層厚度為400 μm、寬度在350~600 μm的范圍內(nèi)時(shí),植入力與涂層寬度存在顯著的線性關(guān)系,其線性模型可表示為F=0.014 6W+0.557 3。同時(shí)實(shí)驗(yàn)過(guò)程中測(cè)得最大植入力為14.2 mN,該值可作為涂層參數(shù)設(shè)計(jì)時(shí)的屈曲力參考指標(biāo)。
(2)提出一種T形截面的涂層,以最小化橫截面積為設(shè)計(jì)目標(biāo)建立了優(yōu)化模型,并通過(guò)遺傳算法分別對(duì)有、無(wú)加工工藝兩種約束條件下的T形截面涂層進(jìn)行了參數(shù)優(yōu)化設(shè)計(jì),得出加工工藝約束對(duì)涂層設(shè)計(jì)存在較大的制約,并提出針對(duì)不同寬度電極的最佳涂層參數(shù)選擇策略。同時(shí)將T形截面涂層與傳統(tǒng)矩形截面涂層進(jìn)行對(duì)比,得出其平均橫截面積僅為同寬度矩形涂層的58.24%、同屈曲力矩形涂層的70.32%,大幅減小了橫截面積,也即是減少了對(duì)腦組織的損傷面積。
(3)通過(guò)Ansys有限元植入仿真和實(shí)物實(shí)驗(yàn)證明T形截面涂層輔助柔性電極植入的可行性及可靠性;同時(shí)T形截面還大幅優(yōu)化了植入時(shí)的植入力,僅為原矩形截面涂層的64.62%,進(jìn)一步說(shuō)明其能減輕對(duì)大腦的植入損傷。