許長(zhǎng)興，張　華，胡文超

(重慶大學(xué) 航空航天學(xué)院，重慶　400044)

?

腦出血中豆紋動(dòng)脈出血易發(fā)點(diǎn)的流固耦合分析

許長(zhǎng)興，張華，胡文超

(重慶大學(xué) 航空航天學(xué)院，重慶400044)

摘要：建立了大腦中動(dòng)脈M1段與豆紋動(dòng)脈(LSA)管壁/血液耦合模型。采用有限元法，利用ANSYS WORKBENCH中Transient Structural 和CFX 相互結(jié)合的專用流固耦合算法Fluid Solid Interface 實(shí)現(xiàn)結(jié)構(gòu)分析和流體分析的雙向耦合計(jì)算，其中假設(shè)血管壁為各項(xiàng)同性線彈性體，血液為不可壓縮牛頓流體，用有限元方法對(duì)大腦中動(dòng)脈段M1段與豆紋動(dòng)脈血管內(nèi)的血流情況進(jìn)行了數(shù)值模擬，在給定進(jìn)出口血壓條件下模擬出豆紋動(dòng)脈彈性管壁壓強(qiáng)、應(yīng)力、應(yīng)變分布與血液流場(chǎng)分布，并由此構(gòu)建了豆紋動(dòng)脈近似等強(qiáng)度模型。臨床顯示顱內(nèi)壓增高與血壓增高可導(dǎo)致腦部丘腦部位供血不足，即豆紋動(dòng)脈末端流速減慢。在近似等強(qiáng)度模型下，通過對(duì)比分析發(fā)現(xiàn)：隨著豆紋動(dòng)脈末端流速的減慢，豆紋動(dòng)脈內(nèi)血壓增高，且距離末端越近，血壓增高幅度越大，導(dǎo)致末段附近等效應(yīng)力大幅度提高，形成管壁破裂易發(fā)點(diǎn)，即為豆紋動(dòng)脈出血危險(xiǎn)部位,與臨床統(tǒng)計(jì)的腦動(dòng)脈出血高發(fā)部位相吻合。

關(guān)鍵詞：豆紋動(dòng)脈；血管模型；有限元分析；流固耦合；等效應(yīng)力

在中國(guó)，腦卒中是僅次于腫瘤和心臟病的第三大致死性疾病(衛(wèi)生部《2011中國(guó)衛(wèi)生統(tǒng)計(jì)年鑒》)。腦卒中包括缺血性腦卒中和出血性腦卒中，其中腦出血防治仍未取得明顯的進(jìn)展，成為腦血管病中致殘和致死率最高的疾病。目前，中國(guó)高血壓腦出血年發(fā)病人數(shù)為200萬，累計(jì)每年約150萬人死于腦出血[1]。近年來國(guó)際上已開始重視對(duì)腦出血的研究，投入了巨大的人力和物力，腦出血的發(fā)病率和死亡率開始下降。而中國(guó)在腦出血領(lǐng)域的研究卻相對(duì)滯后，研究的投入和支持力度有待加強(qiáng)。結(jié)合我國(guó)高血壓腦出血的發(fā)病特點(diǎn)和規(guī)律，深入探索其發(fā)生、發(fā)展機(jī)制，尋找有效預(yù)警標(biāo)志與防治新策略，降低發(fā)生率和死殘率，對(duì)保障國(guó)民身體健康具有重大意義。

高血壓腦出血高發(fā)部位在豆紋動(dòng)脈供血的丘腦基底節(jié)區(qū)，具有獨(dú)特的解剖生理特點(diǎn)：① 獨(dú)特的血流動(dòng)力學(xué)特點(diǎn)。豆紋動(dòng)脈從大腦中動(dòng)脈以直角發(fā)出供應(yīng)基底節(jié)、尾狀核和內(nèi)囊區(qū)，豆紋動(dòng)脈易破裂節(jié)段承受了上級(jí)分支與下級(jí)分支之間約有6.7 kPa 的脈壓差[2-3]。② 特殊的力學(xué)邊界條件。高血壓腦出血的責(zé)任血管以直徑為150～660 μm的小動(dòng)脈為主[4]，該區(qū)域血管周圍主要由疏松的白質(zhì)纖維束構(gòu)成，靠近腦室，血管內(nèi)外壓力差明顯，具備血管易破的力學(xué)邊界條件。因此，明確豆紋動(dòng)脈出血易發(fā)部位的結(jié)構(gòu)特征、流場(chǎng)特征、血管力學(xué)特性和受力特征，為豆紋動(dòng)脈出血預(yù)警提供參照依據(jù)是十分必要的。

近年來，流固耦合分析研究和應(yīng)用取得了飛速的發(fā)展，在生物力學(xué)(如血液的流動(dòng))、新型工程材料(如磁流變液)特性的研究[5]等方面得到了廣泛應(yīng)用。

1豆紋動(dòng)脈

動(dòng)脈起自心臟，從心臟發(fā)出后不斷分支，最后發(fā)展為毛細(xì)血管，其管徑隨分支逐漸減小。在不斷分支與遠(yuǎn)離心臟的過程中，動(dòng)脈血壓逐漸降低，血液流速相對(duì)變緩，管壁隨著血壓的降低與血流作用的減弱逐漸變薄。在人體的生長(zhǎng)進(jìn)化過程中血管最終成為等強(qiáng)度模型，符合自然界不斷進(jìn)化的變化規(guī)律。由于豆紋動(dòng)脈過于細(xì)小，只能測(cè)量出其管徑大小而無法通過檢測(cè)手段得知其管臂厚度及變化，所以通過有限元模擬建立一種在人體正常生理狀態(tài)下的一種血管近似等強(qiáng)度模型是十分必要的。該模型的確定為模擬豆紋動(dòng)脈在顱內(nèi)壓增高、不同等級(jí)的高血壓等病理狀態(tài)下的血流、管壁應(yīng)力、應(yīng)變分析以及出血易發(fā)點(diǎn)預(yù)測(cè)提供了分析基礎(chǔ)。

豆紋動(dòng)脈起自打包中動(dòng)脈M1段，是供應(yīng)紋狀體和內(nèi)囊的主要?jiǎng)用}，易破裂出血，故又名“出血?jiǎng)用}”。胡福廣等[6]測(cè)得大腦中動(dòng)脈的 M1 段長(zhǎng)度為 18.58±4.51 mm，外徑約為3.6 mm，豆紋動(dòng)脈外徑約為5 mm。人體生長(zhǎng)發(fā)育過程中血管分支有兩種方式：絕大部分血管分支是以與主干偏折較小角度的方式來發(fā)出分支血管，其目的是避免血液流向改變過大造成對(duì)管壁的沖擊損傷；另一種分支方式遵循的原則為最短路徑原則[7]。大腦中動(dòng)脈分支的方式即遵循最短路徑原則，也就導(dǎo)致了較多的大腦血管相關(guān)性疾病的發(fā)生[8]。大腦中動(dòng)脈發(fā)出的豆紋動(dòng)脈多以直角或近似直角的方式發(fā)出[9]，豆紋動(dòng)脈幾乎存在于所有的大腦半球，主要起自 M1 段外側(cè)部，形成“S”形彎曲后進(jìn)入前穿質(zhì)。這種曲折的行程可盡可能減緩血液的流速，起到調(diào)節(jié)血壓的作用[10]。

在生理學(xué)上血液為非牛頓流體，血管壁為黏彈性體，在血液與管壁的流固耦合中所有的研究者都做了不同程度的簡(jiǎn)化假設(shè),引入了各種計(jì)算模型?；镜暮?jiǎn)化假設(shè)有兩個(gè)：其一，將血管壁看做剛性或線彈性管壁；其二,把血液看成牛頓流體。Aenis等[11]的研究結(jié)果表明：采用牛頓流體與非牛頓流體模型時(shí),流場(chǎng)的差別不大。體外模擬實(shí)驗(yàn)也得出相同結(jié)果[12]。本文主要分析豆紋動(dòng)脈破裂危險(xiǎn)部位，即出血易發(fā)點(diǎn)，更側(cè)重于血管壁部分的受力變形，所以本文血液采用不可壓縮牛頓流體模型，而血管壁采用線彈性體，相對(duì)于剛性管壁更接近人體真實(shí)環(huán)境。

2腦動(dòng)脈的流固耦合動(dòng)力學(xué)模型

1) 血流動(dòng)力學(xué)方程

假定血液是層流且黏性的、不可壓縮的牛頓流體[13-14]?？刂屏鲃?dòng)的基本方程是不可壓縮Navier-Stokes方程。

(1)

(2)

式中um是流體網(wǎng)格速度。

(3)

(4)

(5)

初始條件：速度壓力均為0。其中：ρ為密度，數(shù)值為1 050 kg/m3；P為壓強(qiáng)；u為速度；μ為運(yùn)動(dòng)黏度系數(shù)，數(shù)值為0.003 5 Pa/s。

2) 血管壁的動(dòng)力學(xué)方程

壁面條件設(shè)定：假定血管為不可壓縮的各項(xiàng)同性材料[15-17],密度ρ=1.06×103kg/m3，彈性模量E=2.7 MPa，泊松比為0.45(接近不可壓縮材料)

3) 流固耦合面條件

① 無滑移；

③ df=ds；

④ σf·nf=σs·ns；

d,σ,n分別標(biāo)示位移、應(yīng)力張量、邊界法向；下標(biāo)f,s分別代表流體和固體。

3模型建立、劃分網(wǎng)格與邊界條件

本文所建模型為M1段直徑為3.6mm、長(zhǎng)度為20mm的血管。豆紋動(dòng)脈從距離M1段起始部位為10mm的位置以直角形式發(fā)出，發(fā)出后豆紋動(dòng)脈在近端有“s”型走向，長(zhǎng)度約25mm。本模型豆紋動(dòng)脈部分血管內(nèi)徑依照血液橫截面積保持不變?yōu)樵瓌t，內(nèi)徑為 0.4mm，分支后內(nèi)徑為為0.283mm，起始部位外徑5mm，隨距離起始部位長(zhǎng)度變大而逐漸減小，保證豆紋動(dòng)脈血管為近似等強(qiáng)度血管。血管三維模型如圖1所示。該模型豆紋動(dòng)脈形態(tài)尺寸與人體實(shí)際構(gòu)造近似，且本文所建立模型為近似等強(qiáng)度模型，使豆紋動(dòng)脈在正常血壓狀態(tài)下各部位等效應(yīng)力近似相等。與實(shí)際情況的主要差別在于，本模型只有一根豆紋動(dòng)脈從大腦動(dòng)脈分支出來，而實(shí)際豆紋動(dòng)脈為多支。本文主要分析豆紋動(dòng)脈的出血易發(fā)點(diǎn)，所以對(duì)模型進(jìn)行簡(jiǎn)化，只進(jìn)行單根豆紋動(dòng)脈的流固耦合分析。

在FSI(fluid solid interface)流固耦合分析技術(shù)中，固液的網(wǎng)格劃分采取不同的方式來進(jìn)行。管壁運(yùn)用Ansys Workbench里的Meshing來進(jìn)行掃掠劃分，并在分支處做局部細(xì)化，共有287 083個(gè)節(jié)點(diǎn)，145 446個(gè)單元，如圖2所示。流體采用ICEM CFD來劃分并在細(xì)節(jié)部位提高精度，共 72 375 個(gè)節(jié)點(diǎn)、348 937個(gè)單元，如圖3所示。在初始條件中設(shè)定壁面無滑移。在邊界條件中，由于豆紋動(dòng)脈易破裂節(jié)段承受了上級(jí)分支與下級(jí)分支之間約6.7 kPa的脈壓差[2-3]，所以假設(shè)進(jìn)出口管壁面的自由度為0，大腦中動(dòng)脈出口壓力為10 kPa，進(jìn)口處壓力為10.5 kPa，豆紋動(dòng)脈出口處壓力設(shè)為4 kPa。

圖2　血管網(wǎng)格劃分

圖3　血液網(wǎng)格劃分

4計(jì)算結(jié)果與分析

4.1模型尺寸確定與參考點(diǎn)設(shè)定

本文在Transient Structrual和CFX兩個(gè)模塊的基礎(chǔ)上，通過專用的流固耦合算法實(shí)現(xiàn)結(jié)構(gòu)分析和流體分析的雙向耦合計(jì)算，在Transient Structrual和CFX中設(shè)置相同的耦合時(shí)間和耦合時(shí)間步，分別為1 s和0.1 s，以實(shí)現(xiàn)流場(chǎng)與結(jié)構(gòu)上的耦合計(jì)算。

計(jì)算結(jié)果在CFX-Post中顯現(xiàn)。為了直觀地分析，本文在豆紋動(dòng)脈內(nèi)壁上選取6個(gè)參考點(diǎn)測(cè)量其等效應(yīng)力值，參考點(diǎn)位置如圖4所示。通過對(duì)模擬結(jié)果的分析發(fā)現(xiàn)：大腦中動(dòng)脈中的流速在0.85～0.90 m/s，比較符合正常人大腦中動(dòng)脈流速，此時(shí)豆紋動(dòng)脈出口處速度約為1.25 m/s。流場(chǎng)如圖5所示，此時(shí)通過不斷調(diào)整豆紋動(dòng)脈外徑使得豆紋動(dòng)脈血管內(nèi)壁各處等效應(yīng)力在42～45 kPa 之間，等效應(yīng)力近似相等，使其基本符合設(shè)計(jì)要求的等強(qiáng)度模型。在模型中豆紋動(dòng)脈起始端外徑為0.5 mm，分叉前外徑為0.45 mm，分叉后外徑為0.318 mm，豆紋動(dòng)脈末端外徑為0.31 mm。等效應(yīng)力如圖6所示。

圖4　參考點(diǎn)設(shè)置

圖5　血液流場(chǎng)

圖6　等效應(yīng)力

4.2出血易發(fā)點(diǎn)結(jié)果分析

顱內(nèi)壓簡(jiǎn)稱ICP，顱內(nèi)壓增高會(huì)導(dǎo)致腦供血不足。Barcroft[18]在1943 年研究發(fā)現(xiàn)：ICP增高至一定值時(shí)腦血流量開始有明顯下降。Mcintosh[19]在動(dòng)物實(shí)驗(yàn)中發(fā)現(xiàn)：當(dāng)動(dòng)物顱內(nèi)壓增高后，動(dòng)物大腦皮質(zhì)平均局部腦血流量降低。同樣臨床醫(yī)學(xué)研究結(jié)果表明：高血壓會(huì)引起腦血管痙攣以及血管調(diào)節(jié)功能減退,導(dǎo)致嚴(yán)重的腦血流量降低,出現(xiàn)不同程度的腦動(dòng)脈供血不足癥狀。上述顱內(nèi)壓與血壓的升高，都導(dǎo)致丘腦和基底節(jié)區(qū)供血的豆紋動(dòng)脈流速降低，所以在本文模擬過程中，將模型中豆紋動(dòng)脈末端的血液流速逐漸降低，分析其流場(chǎng)與管壁的受力變化，找到此狀態(tài)下可能導(dǎo)致豆紋動(dòng)脈出血的危險(xiǎn)點(diǎn)，為臨床醫(yī)學(xué)提供參考。

當(dāng)人體血壓增高、大腦顱內(nèi)壓增高時(shí)，豆紋動(dòng)脈末端血液流速減慢，對(duì)應(yīng)此種人體狀態(tài)，有限元模擬將邊界條件的豆紋動(dòng)脈出口壓強(qiáng)4 kPa改為速度條件，速度線性變化從1.25 m/s降至0.25 m/s，如圖7所示。在此邊界條件下6個(gè)參考點(diǎn)的等效應(yīng)力值隨豆紋動(dòng)脈末端出口速度的變化如圖8所示，從中可以看出：參考點(diǎn)1,2在末端附近，等效應(yīng)力值明顯增加，達(dá)到90 kPa以上；而其余4個(gè)參考點(diǎn)幅值變化不大，且越接近大腦中動(dòng)脈變化幅值越小。由此結(jié)果可以看出：隨著大腦顱內(nèi)壓與血壓的增高，豆紋動(dòng)脈末端血液流速降低，豆紋動(dòng)脈遠(yuǎn)端等效應(yīng)力值顯著增加，相對(duì)于其他部位，更易導(dǎo)致血管破裂出血，形成豆紋動(dòng)脈出血易發(fā)點(diǎn)。

當(dāng)豆紋動(dòng)脈末端血液流速分別為v=1.25 m/s(左)和v=0.25 m/s(右)時(shí)，6個(gè)參考點(diǎn)附近局部等效應(yīng)力如圖9所示，從中對(duì)比可以看出：參考點(diǎn)1與參考點(diǎn)2附近的等效應(yīng)力明顯增大，由45 kPa增加至90 kPa以上，約為原來的2倍；而參考點(diǎn)3附近的等效應(yīng)力增加至65 kPa左右，有一定幅度的提高；參考點(diǎn)4附近增加至60 kPa左右；參考點(diǎn)5、6兩點(diǎn)等效應(yīng)力變?yōu)?0 kPa左右，變化幅度較小。參考點(diǎn)1附近隨豆紋動(dòng)脈末端血液流速減慢對(duì)應(yīng)的局部等效應(yīng)力如圖10所示。從圖9、10可以發(fā)現(xiàn)：隨豆紋動(dòng)脈末端流速的減慢，豆紋動(dòng)脈等效應(yīng)力升高，且距離末端越近，增高幅度越大。

圖8　各參考點(diǎn)等效應(yīng)力值隨末端速度變化的曲線

圖9　豆紋動(dòng)脈末端血液流速分別為v=1.25 m/s(左)，v=0.25 m/s(右)時(shí)6個(gè)參考點(diǎn)附近局部等效應(yīng)力

圖10　參考點(diǎn)1附近隨豆紋動(dòng)脈末端流速血液減慢

5討論

目前大多數(shù)血管生物力學(xué)主要針對(duì)冠狀動(dòng)脈、頸動(dòng)脈與動(dòng)脈瘤的血管血流分析，針對(duì)豆紋動(dòng)脈的生物力學(xué)分析相對(duì)較少。大多數(shù)生物力學(xué)主要單獨(dú)分析血管壁的力學(xué)屬性或人體血液流動(dòng)情況，而對(duì)血流-血管耦聯(lián)的流固耦合的有限元分析相對(duì)較少。本文血管采用彈性材料，用有限元方法對(duì)豆紋動(dòng)脈血管內(nèi)的血流情況進(jìn)行了數(shù)值模擬，在給定正常生理血壓條件下模擬出豆紋動(dòng)脈彈性管壁應(yīng)力、應(yīng)變分布與血液流場(chǎng)分布，為臨床醫(yī)學(xué)提供參考。

本文模擬結(jié)果顯示：當(dāng)人體顱內(nèi)壓增高、血壓增高導(dǎo)致豆紋動(dòng)脈末端血液流速減慢時(shí)，豆紋動(dòng)脈內(nèi)等效應(yīng)力增高，且距離末端越近，應(yīng)力增高幅度越大，形成管壁破裂易發(fā)點(diǎn)，即豆紋動(dòng)脈出血危險(xiǎn)部位；豆紋動(dòng)脈出血易發(fā)點(diǎn)位于豆紋動(dòng)脈遠(yuǎn)端。根據(jù)國(guó)際疾病分類(International Classification of Diseases-10,ICD-10) 的分級(jí)方法[20]，將患者的出血部位分成深部(基底核區(qū)和丘腦)、腦葉(單個(gè)或多腦葉)、腦干、小腦、腦室(排除腦實(shí)質(zhì)出血破入腦室)、多個(gè)部位(同時(shí)存在2個(gè)部位的出血)。Wang和Ruiz-Sandoval等[21-22]的臨床調(diào)查研究結(jié)果表明：在腦出血的各出血部位中，以深部(基底核區(qū)、丘腦)出血居多，超過總量的60%。本文的模擬結(jié)果顯示：豆紋動(dòng)脈出血易發(fā)點(diǎn)位于豆紋動(dòng)脈遠(yuǎn)端，即豆紋動(dòng)脈供血的丘腦、基底節(jié)區(qū)，與臨床調(diào)查結(jié)果相符。

腦卒中問題現(xiàn)主要處在醫(yī)療治療方面，針對(duì)其發(fā)生機(jī)理與預(yù)防的研究相對(duì)較少。本文主要從出血性腦卒中的發(fā)病機(jī)制出發(fā)，通過力學(xué)手段分析人體的醫(yī)學(xué)問題，通過CAD軟件對(duì)人體大腦血管進(jìn)行三維建模，建立大腦中動(dòng)脈與豆紋動(dòng)脈模型，通過血流血管流固耦合的有限元分析研究大腦中動(dòng)脈與豆紋動(dòng)脈的血液流動(dòng)情況和血管壁的應(yīng)力應(yīng)變狀況，得出高血壓腦出血高發(fā)部位位于基底節(jié)區(qū),從而為臨床醫(yī)學(xué)提供力學(xué)參數(shù)參考。

參考文獻(xiàn)：

[1]FEIGIN V L,LAWES C M,BENNETT D A,et al. Worldwide stroke incidence and early case fatality reported in 56 population-based studies:a systematic review[J].Lancet Neurol,2009,8(4):355-369.

[2]RHODIN J A.The ultrastructure of mammalian arterioles and precapillary sphincters[J].J Ultrastruct Res,1967,18(1):181-223.

[3]WIEDERHIELM C A,WOODBURY J W,KIRK S,et al. Pulsatile Pressures in the Microcirculation of Frog’s Mesentery[J].Am J Physiol,1964,207:173-176.

[4]TAKAKEBAYASHI S,KANEKO M.Electron microscopic studies of ruptured arteries in hypertensive intracerebral hemorrhage[J].Stroke,1983,14(1):28-36.

[5]CHEN S,HUANG J,SHU H Y,et al.Analysis and Testing of Chain Characteristics and Rheological Properties for Magnetorheological Fluid[J].Advances in Materials Science and Engineering,2013,29:691.

[6]胡福廣,張皓峰,岳立峰.腦外側(cè)裂區(qū)顯微外科解剖學(xué)研究[J].中國(guó)神經(jīng)精神疾病雜志,2010(5):272-275.

[7]秦建軍,邊進(jìn)才,劉俊文，等.大腦中動(dòng)脈的顯微解剖[J].河北職工醫(yī)學(xué)院學(xué)報(bào),1995(4):36-41.

[8]毛增榮,張鋆.腦底諸動(dòng)脈及其分枝[J].解剖學(xué)報(bào)， 1958(2):37-56.

[9]王克臻,李巖,馬威,等.正常成人腦動(dòng)脈的 MRA 測(cè)量和應(yīng)用解剖研究[J].中國(guó)臨床解剖學(xué)雜志,2008(4):392-395.

[10]韓 卉,張為龍.人大腦中動(dòng)脈中央支腦內(nèi)段的研究[J].解剖學(xué)報(bào)，1993,24(1):12-16.

[11]AENIS M,STANCAMPIANO A P.Modeling of flow in a straight stented and nonstented side wall aneurysm model[J].Journal of Biomechanical Engineering,1997,119(2):206-212.

[12]ERNEMANN U U,GRONEWALLER E,et al.Influence of geometric and hemodynamic parameters on aneurysm visualization during three-dimensional rotational angiography:An in vitro study[J].AMERICAN JOURNAL OF NEURORADIOLOGY,2003,24(4):597-603.

[13]MILNER J S,MOORE J A,RUTT B K,et al.Hemodynamics of human carotid artery bifurcations:computational studies with models reconstructed from magnetic resonance imaging of normal subjects[J].Journal of vascular surgery,1998,28(1):143-156.

[14]WOOD N B,WESTON S J.KILNER P J,et al.Combined MR imaging and CFD Simulation of flow in the human descending aorta[J].Journal of Magnetic Resonance Imaging,2001,13(5):699-713.

[15]Di MARTINO E S,GUADAGNI G,FUMERO A,et al.Fluid-structure interaction within realistic three-dimensional models of the aneurysmatic aorta as a guidance to assess the risk of rupture of the aneurysm[J].Med Eng Phys,2001,23(9):647-655.

[16]GAO F,GUO Z H,SAKAMOTO M,et al.Fluid-structure Interaction within a Layered Aortic Arch Model[J].Journal of Biological Physics,2006,32(5):435-454.

[17]KELLY S,O’ROURKE M.Fluid,solid and fluid-structure interaction simulations on patient-based abdominal aortic aneurysm models[J].Proc Inst Mech Eng H,2012,226(4):288-304.

[18]BARCROFT H,EDHOLM O G.The effect of temperature on blood flow and deep temperature in the human forearm[J].J Physiol,1943,10(2):5-20.

[19]MCINTOSH T K,NOBEL L,ANDEREW B,et al.Traumatic brain injury in the rat:characterization of a midline fluid-percussion model[J].Cent Nerv Syst Trauma,1987,4:119-134.

[20]Van DRIMMELEN-KRABBE J J,BRADLEY W G,ORGOGOZO J M,et al.The application of the international statistical classification of diseases to neurology:ICD-10 Na[J].J Neurol Sci,1998,161(1):2-9.

[21]WANG W J,LU J J,WANG Y J,et al.Clinical Characteristics,Management,and Functional Outcomes in Chinese Patients Within the First Year After Intracerebral Hemorrhage:Analysis from China National Stroke Registry[J].CNS Neurosci Ther,2012,18(9):773-780.

[22]RUIZ-SANDOVAL J L,CHIQUETE E,GARATE-CARRILLO A,et al.Spontaneous intracerebral hemorrhage in Mexico:results from a Multicenter Nationwide Hospital-based Registry on Cerebrovascular Disease (RENAMEVASC)[J].Rev Neurol,2011,53(12):705-712.

(責(zé)任編輯何杰玲)

Fluid-Solid Interactive Analysis of Bleeding-Prone Points of LSA in Cerebral Hemorrhage

XU Chang-xing, ZHANG Hua, HU Wen-chao

(College of Aerospace Engineering, Chongqing University, Chongqing 400044, China)

Abstract:The wall-blood interactive model of M1 segment of the middle cerebral artery and the lenticulostriate artery (LSA) of brain was established, and then two-way interactive calculation of structural analysis and fluid analysis were achieved by finite element method (FEM) and the special Fluid-Solid interactive algorithm which combinesd Transient Structural and CFX. In experiment, the vascular wall was assumed as the isotropic linear and elastomeric, and the blood was incompressible Newtonian fluid. The numerical simulation was carried out in terms of the blood flow of M1 segment of the middle cerebral artery segment and the LSA by FEM. Under the given conditions of entrance and outlet of vessels’ blood pressure, we simulated the pressure, stresses, strains distribution and flow field distribution and built lenticulostriate artery approximately equal strength model. Under the model, the results of experiment and simulation show that the blood pressure within the LSA increases when the flowing velocity of the LSA is slow. Especially, the magnitude of the blood pressure is greater when it close the end. It leads to the substantial increasing of the equivalent stress near the end, and forms the rupture-prone points of the wall which are dangerous parts of the LSA bleeding. These results of experiment and simulation have a high match with the clinical statistics on the cerebral arterial bleeding site.

Key words:lenticulostriate artery; vessel model; finite element analysis; fluid-solid interaction; equivalent stress

文章編號(hào)：1674-8425(2016)04-0066-07

中圖分類號(hào)：R54

文獻(xiàn)標(biāo)識(shí)碼：A

doi:10.3969/j.issn.1674-8425(z).2016.04.012

作者簡(jiǎn)介：許長(zhǎng)興(1991—)，男，碩士，主要從事流固耦合分析和結(jié)構(gòu)優(yōu)化設(shè)計(jì)研究。

基金項(xiàng)目：國(guó)家重點(diǎn)基礎(chǔ)研究計(jì)劃(973計(jì)劃)資助項(xiàng)目“出血性腦卒中發(fā)生發(fā)展機(jī)制及其防治的基礎(chǔ)研究”(2014CB541600)

收稿日期：2016-01-21

引用格式：許長(zhǎng)興，張華，胡文超.腦出血中豆紋動(dòng)脈出血易發(fā)點(diǎn)的流固耦合分析[J].重慶理工大學(xué)學(xué)報(bào)(自然科學(xué))，2016(4):66-72.

Citation format：XU Chang-xing, ZHANG Hua, HU Wen-chao.Fluid-Solid Interactive Analysis of Bleeding-Prone Points of LSA in Cerebral Hemorrhage[J].Journal of Chongqing University of Technology(Natural Science)，2016(4):66-72.