重力應(yīng)激下心血管循環(huán)系統(tǒng)建模仿真研究進(jìn)展

王漢青，趙金萍，陳娟

（湖南工業(yè)大學(xué)土木工程學(xué)院，湖南株洲412007）

重力應(yīng)激下心血管循環(huán)系統(tǒng)建模仿真研究進(jìn)展

王漢青，趙金萍，陳娟

（湖南工業(yè)大學(xué)土木工程學(xué)院，湖南株洲412007）

重力應(yīng)激通過改變血液流體靜壓對人體心血管循環(huán)系統(tǒng)產(chǎn)生不利影響。了解心血管循環(huán)系統(tǒng)在不同重力應(yīng)激下的工作機(jī)制，有利于建立合理有效的防護(hù)措施?？偨Y(jié)了人體心血管循環(huán)系統(tǒng)的建模方法，回顧了適用于重力應(yīng)激下的心血管循環(huán)系統(tǒng)模型的發(fā)展進(jìn)程，闡述了重力變化對心血管循環(huán)系統(tǒng)影響的研究及應(yīng)用現(xiàn)狀，介紹了作者的研究成果，并指明了未來的研究發(fā)展方向。

心血管循環(huán)系統(tǒng)；數(shù)學(xué)模型；重力應(yīng)激；建模仿真；熱舒適

0 引言

人體心血管循環(huán)系統(tǒng)負(fù)責(zé)向組織、器官輸送用于新陳代謝的O2和其他營養(yǎng)物質(zhì)，是維持人體各項生理機(jī)能正常運行的重要保障[1]。由于人類長期在地球上生活，心血管循環(huán)系統(tǒng)已經(jīng)適應(yīng)了1 G的重力環(huán)境，所以當(dāng)重力改變時，血液的流體靜壓將會隨之改變甚至消失，進(jìn)而對心血管循環(huán)系統(tǒng)產(chǎn)生不利影響。隨著人類進(jìn)入宇宙空間機(jī)會的增多以及停留時間的延長，失重環(huán)境下人體心血管循環(huán)系統(tǒng)響應(yīng)越來越受到國內(nèi)外科學(xué)工作者的重視。失重條件下血液的重新分布會影響人體的熱調(diào)節(jié)系統(tǒng)，改變?nèi)梭w溫度分布，使航天員處于不舒適狀態(tài)，進(jìn)而降低工作效率、危害航天員的身體健康。因此，研究失重/微重力環(huán)境對人體心血管循環(huán)系統(tǒng)各項生理參數(shù)的影響以及各參數(shù)隨重力改變的變化規(guī)律，對于預(yù)防重力變化對人體心血管循環(huán)系統(tǒng)的不良影響，提高航天飛行員的健康舒適以及保證正常工作能力具有重要意義。國內(nèi)外的研究表明，對心血管循環(huán)系統(tǒng)進(jìn)行建模仿真是研究心血管系統(tǒng)受內(nèi)外環(huán)境變化以及各種生理、病理狀態(tài)影響的有效手段。

心血管循環(huán)系統(tǒng)由心臟、動脈、毛細(xì)血管網(wǎng)絡(luò)及靜脈構(gòu)成。血液循環(huán)分為體循環(huán)和肺循環(huán)，機(jī)體通過神經(jīng)和體液機(jī)制對心血管循環(huán)系統(tǒng)的活動進(jìn)行調(diào)節(jié)，協(xié)調(diào)各器官間的血液分配，保障系統(tǒng)的正常運行[2]，重力對于心血管循環(huán)系統(tǒng)的影響通過改變血液的流體靜壓實現(xiàn)[3]。近年來，不同重力應(yīng)激對心血管循環(huán)系統(tǒng)模型的影響成為國內(nèi)外學(xué)者關(guān)注的熱點。本文總結(jié)了人體心血管循環(huán)系統(tǒng)的建模方法，回顧了適用于重力應(yīng)激下的心血管循環(huán)系統(tǒng)模型的發(fā)展進(jìn)程，闡述了重力變化對心血管循環(huán)系統(tǒng)影響的研究及應(yīng)用現(xiàn)狀，指出了存在的問題與解決方法以及未來的發(fā)展方向。

1 心血管循環(huán)系統(tǒng)的建模方法

目前，主要通過臨床試驗、動物實驗以及建模仿真等方法研究人體心血管循環(huán)系統(tǒng)。臨床試驗方法直接對人體進(jìn)行測量和試驗，盡管此種方法最準(zhǔn)確和直觀，但是不可避免地受到了倫理道德和設(shè)備技術(shù)手段上的限制；某些生理參數(shù)不可能在人體活體上直接測量，而間接測得的許多生理參數(shù)不能保證其準(zhǔn)確性，這些成為臨床試驗方法無法規(guī)避的弊端。由于人體與動物在結(jié)構(gòu)、生理上存在差異，因此動物實驗測得的數(shù)據(jù)可信度較低，不能用于分析人體相應(yīng)參數(shù)的變化規(guī)律。而仿真模型不僅可以突破上述兩種方法的限制，而且能模擬完成極端條件下的實驗，實現(xiàn)對未來事件的預(yù)測[4]。

心血管循環(huán)系統(tǒng)的建模仿真方法主要包括電網(wǎng)絡(luò)模型法、鍵合圖法、有限元分析法以及系統(tǒng)辨識法等[5]。電網(wǎng)絡(luò)模型法的理論基礎(chǔ)是利用電氣網(wǎng)絡(luò)的相關(guān)特性類比于流體網(wǎng)絡(luò)，用電學(xué)量表示血液動力學(xué)參數(shù)，用電路理論中的基爾霍夫定律描述血液循環(huán)系統(tǒng)各血管間的串/并聯(lián)關(guān)系，建立血液循環(huán)系統(tǒng)的參數(shù)模型；鍵合圖法以圖形方法來描述系統(tǒng)的結(jié)構(gòu)，是對流體力學(xué)進(jìn)行動態(tài)數(shù)字仿真時有效的建模方法；有限元分析法主要用于非線性系統(tǒng)，要求對系統(tǒng)的有關(guān)生理、解剖及物理參數(shù)有詳盡的了解，可用于準(zhǔn)確模擬整個心血管循環(huán)系統(tǒng)；系統(tǒng)辨識法的理論基礎(chǔ)是系統(tǒng)辨識技術(shù)，建模過程中把整個系統(tǒng)看作一個非線性黑箱，利用參數(shù)辨識方法分析估計系統(tǒng)的輸入輸出關(guān)系，建立相應(yīng)的數(shù)學(xué)模型[5]。利用建模仿真方法研究人體心血管循環(huán)系統(tǒng)的運行機(jī)制及其動態(tài)變化規(guī)律，已成為近年來心血管系統(tǒng)研究的發(fā)展趨勢。

2 國內(nèi)外研究進(jìn)展及現(xiàn)狀

2.1 國外研究進(jìn)展

為模擬重力應(yīng)激對人體心血管循環(huán)系統(tǒng)的影響，需要在原有心血管循環(huán)系統(tǒng)數(shù)學(xué)模型中引入重力對血流靜壓的影響。David G. Boyers等人[6]通過改進(jìn)A. C. Guyton等人[7]提出的心血管系統(tǒng)模型，建立了能用于模擬體位變化對人體心血管系統(tǒng)生理參數(shù)影響的集總參數(shù)模型。模型由心血管受控系統(tǒng)和控制系統(tǒng)構(gòu)成，受控系統(tǒng)由手臂模塊、內(nèi)臟模塊、腿部模塊、外周組織模塊、靜脈腔以及右心模塊、肺循環(huán)系統(tǒng)以及左心模塊及動脈系統(tǒng)模塊7部分構(gòu)成，模型中血液分布受自主神經(jīng)系統(tǒng)調(diào)節(jié)，通過接收動脈壓及肺部血容量的變化信號，調(diào)節(jié)心輸出量及外周血管阻力及順應(yīng)性。該模型能夠準(zhǔn)確預(yù)測心輸出量、心率、每搏輸出量、血容量及動脈壓力等參數(shù)在不同體位時的穩(wěn)態(tài)值。

Francois M. Melchior等人[8-9]總結(jié)了模擬立位應(yīng)激需要考慮的生理機(jī)制及已有模型的解決方法。心血管循環(huán)系統(tǒng)模型分為血液動力學(xué)模型即受控系統(tǒng)，反射控制系統(tǒng)。血液動力學(xué)模型主要考慮血管的壓力-血流量關(guān)系及壓力-容積關(guān)系、心臟泵血作用、血管床局部循環(huán)的自動調(diào)節(jié)作用、毛細(xì)血管壁滲透性以及血管外壓力對血管半徑的影響；反射控制模型大多為經(jīng)驗式，主要有S曲線模型、積分器模型及線性-非線性函數(shù)模型。S曲線模型將外周阻力、心臟收縮性及心率的變化集總為一個壓力信號，頸動脈壓力感受器通過反饋壓力信號直接調(diào)節(jié)系統(tǒng)動脈壓力，該信號與系統(tǒng)動脈壓力的關(guān)系為一條S型曲線；積分器模型將中樞神經(jīng)系統(tǒng)（central nervous system，CNS）視作積分器，根據(jù)CNS激發(fā)水平將輸入壓力信號積分為單位函數(shù)（unique function）作用到受控位置；線性-非線性函數(shù)模型分別考慮心率、心臟收縮性、動脈阻力及靜脈緊張度對系統(tǒng)動脈壓力的反射控制作用、建立線性/非線性函數(shù)。由于S曲線模型中各影響因素集總為一個壓力信號，模型的物理意義不明顯，而積分器模型的建立理論依據(jù)不充分，且作者沒有進(jìn)行實驗驗證模型的正確性，所以線性-非線性函數(shù)模型因其以實驗數(shù)據(jù)為基礎(chǔ)、控制機(jī)理清晰而得到廣泛應(yīng)用。根據(jù)以上理論基礎(chǔ)及建模方法建立了用于0～40 mmHg下體負(fù)壓（lower body negative pressure，LBNP）條件下心血管系統(tǒng)短期響應(yīng)的數(shù)學(xué)模型，模擬結(jié)果與實驗值一致。未來模型的發(fā)展應(yīng)致力于擴(kuò)大適用范圍，使同一模型能夠用于預(yù)測不同模擬失重條件下心血管循環(huán)系統(tǒng)的生理響應(yīng)，有助于找出系統(tǒng)在各種立位應(yīng)激下生理參數(shù)變化的區(qū)別。

為了模擬心血管循環(huán)系統(tǒng)在不同立位應(yīng)激下的短期瞬態(tài)響應(yīng)行為，Thomas Heldt等人[10-12]建立了十二腔集總參數(shù)模型代表血液動力學(xué)系統(tǒng)，控制系統(tǒng)考慮動脈壓力反射作用和心肺反射作用，分別仿真頭上位傾斜（head up tilt，HUT）、臺架試驗（stand test）及LBNP條件下的血液動力學(xué)瞬態(tài)響應(yīng)。模擬結(jié)果與試驗結(jié)果對比表明，該模型在上述立位應(yīng)激下均有較高的準(zhǔn)確性，因而可用于重力應(yīng)激防護(hù)措施的研究開發(fā)。

考慮到呼吸時血液中氧氣與二氧化碳含量的變化均會影響局部血管的阻力，F(xiàn)ranz Kappel等人[13]將心血管系統(tǒng)與呼吸系統(tǒng)結(jié)合，模擬了LBNP條件下系統(tǒng)阻力、非壓力靜脈容積、靜脈順應(yīng)性、心率以及心臟收縮性等參數(shù)的變化規(guī)律及其相互作用，并且與實驗值進(jìn)行了對比，得知該模型能夠準(zhǔn)確地預(yù)測各參數(shù)變化趨勢以及變化范圍，得出心血管循環(huán)系統(tǒng)主要通過非壓力容積即血管順應(yīng)性的變化消除LBNP影響的結(jié)論。

K. Van Heusden等人[14]在心血管循環(huán)系統(tǒng)模型中引入了生物力學(xué)的觀點，考慮靜脈血管的粘彈性應(yīng)力松弛行為的基礎(chǔ)，準(zhǔn)確預(yù)測了HUT初期心血管循環(huán)系統(tǒng)的響應(yīng)情況。當(dāng)人體頭上位傾斜某一角度時，靜脈由于所受重力作用減小，因此血壓及血管容量降低，應(yīng)力松弛作用能夠減緩靜脈血管壁在血壓降低初期的收縮作用，在一定程度上降低血管容量下降速度。模擬得到在頭上位傾斜前30 s，調(diào)節(jié)動脈血壓的主要因素是下體靜脈淤積（venous pooling）作用，而非壓力反射作用。了解體循環(huán)系統(tǒng)中靜脈系統(tǒng)的生理機(jī)制有助于研究HUT時動脈壓力的瞬態(tài)響應(yīng)規(guī)律。

2.2 國內(nèi)研究進(jìn)展

國內(nèi)對重力應(yīng)激下心血管循環(huán)系統(tǒng)的研究起步較晚，其中最具代表性的是清華大學(xué)的白凈。1988年，Jaron和白凈建立了用于加速應(yīng)力下的心血管系統(tǒng)模型，此后在該模型基礎(chǔ)上不斷改進(jìn)。利用改進(jìn)模型仿真了人體在LBNP時血壓和心率的變化，結(jié)果與同步進(jìn)行的人體LBNP實驗結(jié)果吻合；仿真研究了立位應(yīng)激時血容量降低對血壓、心率和休克指數(shù)的影響，證明血容量降低是引發(fā)人體失重后立位耐力不良的主要因素之一[15]。馮岱雅等人[16]依據(jù)標(biāo)準(zhǔn)人體生理參數(shù)，建立了適用于各種重力環(huán)境下的人體心血管系統(tǒng)仿真模型，仿真模擬了失重條件下人體腦血流量與血壓的變化，分析血容量降低對人體心血管系統(tǒng)的影響，為研究失重對人體心血管功能的影響及其機(jī)制提供了新的思路和方法。劉洋等人[17]建立了分布式心血管系統(tǒng)數(shù)學(xué)模型，通過調(diào)整重力沿血管方向的分量模擬人體不同的傾斜角度，研究傾斜角度變化對人體心血管系統(tǒng)運行的影響。結(jié)果表明隨著傾斜角度的增加，血管的平均血壓近似呈現(xiàn)指數(shù)變化趨勢，仿真結(jié)果與相關(guān)的人體實驗數(shù)據(jù)基本一致。該模型在模擬仿真人體心血管循環(huán)系統(tǒng)對于傾斜測試中角度變化的瞬態(tài)響應(yīng)方面具有較好的適用性，但是由于人們越來越需要長時間停留在失重或模擬失重環(huán)境中，需要考慮長期失重使心血管系統(tǒng)產(chǎn)生的適應(yīng)性變化，所以需要將心臟與血管功能和結(jié)構(gòu)的慢性適應(yīng)變化、神經(jīng)-體液調(diào)節(jié)以及血容量調(diào)節(jié)等機(jī)制耦合進(jìn)來。

2.3 應(yīng)用現(xiàn)狀

目前，心血管循環(huán)系統(tǒng)模型已經(jīng)被廣泛地應(yīng)用于醫(yī)學(xué)、軍事、航空航天等多個領(lǐng)域。醫(yī)學(xué)研究中采用傾斜測試（tilt test）檢測血液動力學(xué)參數(shù)的響應(yīng)行為，用于心血管系統(tǒng)風(fēng)險等級評價以及病癥診斷。Virginie Le Rolle等人[18]通過耦合心室模型、循環(huán)系統(tǒng)模型以及自主神經(jīng)系統(tǒng)模型，并且利用進(jìn)化算法辨識具體患者的生理參數(shù)，實現(xiàn)了對糖尿病患者在傾斜測試中的心率變異（heart rate variability，HRV）以及迷走神經(jīng)、交感神經(jīng)活動水平的預(yù)測。該模型不僅能夠準(zhǔn)確地模擬不同患者的HRV，進(jìn)而評估患者自主神經(jīng)系統(tǒng)損傷程度，而且可以用于疾病治療手段有效性的預(yù)測。Einly Lim等人[19]基于13位健康者的傾斜測試數(shù)據(jù)建立了心血管循環(huán)系統(tǒng)模型，模擬健康者與充血性心力衰竭患者在不同角度傾斜測試中心血管系統(tǒng)的血液動力學(xué)響應(yīng)。為了能將模型用于充血性心力衰竭患者，原有模型中降低了心臟收縮性，增大了肺循環(huán)及體循環(huán)的基本阻力，增大了總血容量，同時降低了壓力反射系統(tǒng)敏感度，因此保證了平均動脈壓、平均左心房壓力、平均肺動脈壓力、平均中央靜脈壓力、心率、每搏量及心輸出量等參數(shù)的預(yù)測值在試驗值范圍內(nèi)。結(jié)果表明，與健康者相比，充血性心力衰竭患者進(jìn)行傾斜測試時，心率及體血管阻力的增加幅度較小，而左心室舒張末期血容量、心輸出量及每搏量下降幅度均減小，并且各參數(shù)響應(yīng)明顯滯后。

心血管循環(huán)系統(tǒng)的正常運行在軍事活動中尤為重要。軍用機(jī)飛行員或特技飛行員在飛行中受加速度變化影響時，心血管循環(huán)系統(tǒng)將受到極大考驗，嚴(yán)重時將會阻礙飛行員的安全操作。預(yù)測飛行員在特定飛行演習(xí)中血壓、心率等心血管系統(tǒng)響應(yīng)的變化規(guī)律，有助于合理設(shè)計演習(xí)方案。Liu Yang等人[20]建立了分布式心血管系統(tǒng)模型，并模擬進(jìn)行推拉動作（push-pull maneuver，PPM）時頸動脈壓力反射作用對+Gz耐力的影響，以及推拉效應(yīng)（push-pull effect，PPE）引發(fā)的交感神經(jīng)響應(yīng)行為。比較PPM飛行及對比飛行（control run）下的模擬結(jié)果，可知這兩種情況下，縮短頸動脈壓力反射的延遲時間及增大增益系數(shù)對+Gz耐力的影響基本一致，即PPE不會影響頸動脈壓力反射功能；PPM模式主要由預(yù)設(shè)G值、“推”階段初始速度、“推”階段G值及時長、“拉”階段初始速度以及“拉”階段G值6個參數(shù)描述，預(yù)測比較不同模式下PPM中交感傳出神經(jīng)脈動頻率發(fā)現(xiàn)，“推”“拉”階段的G值是輸出脈動頻率的重要影響因素。

在航空航天領(lǐng)域，飛行后立位耐力（postfligt orthostatic intolerance, POI）限制了宇航員著陸時處理突發(fā)狀況的能力，并且減緩了飛行后工作能力的恢復(fù)，成為亟待解決的問題。Justin Broskey和M. Keith Sharp[21]建立了包含心臟/肺、體動脈、頭靜脈、中央靜脈及尾靜脈的五腔模型，研究了立位耐力對心血管系統(tǒng)的影響，得到毛細(xì)血管濾過作用（capillary filtration）是引發(fā)POI的主要因素。通過降低毛細(xì)血管滲透性及激發(fā)血容量代償機(jī)制可以減輕POI現(xiàn)象，而后者的作用明顯優(yōu)于前者，為防護(hù)措施的開發(fā)提供了理論依據(jù)。

3 作者研究成果

筆者以Mauro Ursino[22]建立的心血管循環(huán)系統(tǒng)數(shù)學(xué)模型為基礎(chǔ)，引入重力變化對心血管靜脈壓的影響，利用體位變化來模擬重力變化對靜脈壓的影響，利用Simulink/Matlab建立包含左心模塊、右心模塊、體循環(huán)模塊、肺循環(huán)模塊及神經(jīng)控制模塊的數(shù)字仿真模型。分別模擬正常重力環(huán)境、失重環(huán)境及微重力環(huán)境下人體心血管基本生理參數(shù)的響應(yīng)，為研究重力改變時人體熱舒適程度的變化奠定了基礎(chǔ)。

體動脈壓力、左右心室壓力、左心室容積及左右心室血流量在失重與微重力環(huán)境下的波形與正常重力環(huán)境下基本一致（見圖1至圖6）。其中人體的收縮壓與舒張壓在失重環(huán)境下分別為112.62 mmHg及66.06 mmHg，在微重力環(huán)境下分別為112.64 mmHg及66.09 mmHg。與正常重力環(huán)境下的模擬結(jié)果對比，失重環(huán)境下人體舒張壓降低14.86 mmHg，收縮壓降低4.85 mmHg，馮岱雅[16]通過建模仿真得到失重條件下人體舒張壓與收縮壓相比正常重力環(huán)境下分別降低17 mmHg及3 mmHg，與本文結(jié)果基本一致。另外圖1顯示了人體心動周期隨重力減小而增大，即正常重力環(huán)境下人體心率較失重環(huán)境下升高。

圖1 不同重力環(huán)境下體動脈壓力-時間仿真曲線Fig.1Emulational curves of systemic arterial pressure-time under different gravitional conditions

圖2 不同重力環(huán)境下左心室壓力-時間仿真曲線Fig.2Emulational curves of left ventricle pressure-time under different gravitional conditions

圖3 不同重力環(huán)境下右心室壓力-時間仿真曲線Fig.3Emulational curves of right ventricle pressure-time under different gravitional conditions

圖4 不同重力環(huán)境下左心室容積-時間仿真曲線Fig.4Emulational curves of left ventricle volume-time under different gravitional conditions

圖5 不同重力環(huán)境下左心室血流量仿真曲線Fig.5Emulational curves of left ventricle blood flow rate under different gravitional conditions

圖6 不同重力環(huán)境下右心室血流量仿真曲線Fig.6Emulational curves of right ventricle blood flow rate under different gravitional conditions

左心室壓力的變化范圍在失重及微重力環(huán)境下分別為4.04 ～ 137.62 mmHg及4.03～137.61 mmHg（見圖2）；右心室壓力在失重條件下為2.42 ～28.22 mmHg。

微重力環(huán)境下為2.42～ 28.20 mmHg（見圖3）。與正常重力環(huán)境下相比較，二者的變化范圍均會隨著重力的減小而增大，即失重條件下人體左右心室壓力均升高。

由圖4可知，左心室容積隨著重力減小而增大。當(dāng)人體處于失重環(huán)境中時，左心室容積為54.25～136.37 mL，微重力環(huán)境下為54.24～ 136.36 mL，相比于正常重力環(huán)境下的50.65～ 111.15 mL，有了大幅度升高。由于左心室容積的最大值與最小值之差為心臟每搏輸出量，因此由圖4亦可知重力減小將引起每搏輸出量增大。

與左心室容積的變化趨勢相同，左右心室血流量在重力減小時同樣呈現(xiàn)出升高的趨勢。由圖5及圖6可知，失重條件下左心室血流量最大值為717.45 mL/ min，微重力環(huán)境下為716.18 mL/min，較正常重力環(huán)境下的579.55 mL/min分別升高了137.90 mL/min及136.63 mL/min；右心室血流量在失重與微重力環(huán)境下的最大值分別為486.71 mL/min及486.65 mL/min，相對正常重力條件下分別升高59.20 mL/min及59.14 mL/ min，相比左心室血流量升高幅度較小。

4 結(jié)論與展望

通過對國內(nèi)外適用于重力應(yīng)激下的心血管循環(huán)系統(tǒng)模型的總結(jié)分析，得出如下結(jié)論：

受試驗條件所限，特定疾病患者的生理參數(shù)無法獲得，阻礙了心血管循環(huán)系統(tǒng)模型在病理分析及疾病防護(hù)措施研發(fā)等方面的發(fā)展。采用模糊辨識等方法獲得具體生理參數(shù)，建立適用于特定患者的心血管循環(huán)系統(tǒng)模型，將是這一問題的有效解決辦法。

心血管模型正朝著功能多樣化的方向發(fā)展，不僅能用于研究心血管系統(tǒng)的運行機(jī)制以及重力應(yīng)激下的防護(hù)措施，也能用于預(yù)測患者心血管疾病的發(fā)展過程，還可以用于測試醫(yī)療儀器的臨床應(yīng)用效果。

目前少有重力變化對人體熱舒適影響的研究，因此，模擬失重/微重力環(huán)境下人體各部位血液分布，有助于深化對失重條件下人體熱舒適機(jī)理的認(rèn)識，建立失重/微重力環(huán)境下人體熱調(diào)節(jié)系統(tǒng)數(shù)學(xué)模型，進(jìn)而提高航天飛行過程中人體的熱舒適程度并建立相應(yīng)的熱舒適評價指標(biāo)。

[1]Shim E B，Sah J Y，Youn C H. Mathematical Modeling of Cardiovascular System Dynamics Using a Lumped Parameter Method[J].The Japanese Journal of Physiology，2004，54(6)，545-553.

[2]何為，余傳祥.心血管動力學(xué)參數(shù)測量原理和臨床應(yīng)用[M].北京：科學(xué)出版社，2010：276-278. He Wei，Yu Chuanxiang. Cardiovascular Kinetic ParametersMeasurement Principle and Clinical Applications[M]. Beijing：Science Press，2010：276-278.

[3]陳杰.長期模擬失重下心臟適應(yīng)性變化的實驗研究[D].西安：第四軍醫(yī)大學(xué)，1992. Chen Jie.Cardiac Adaptation to Long-Term Simulated Weightlessness in Rat[D]. Xi’an The Fourth Military Medical University，1992.

[4]Liang F Y，Takagi S， Himeno R，et al.Biomechanical Characterization of Ventricular-Arterial Coupling During Aging：A Multi-Scale Model Study[J]. Journal of Biomechanics，2009，42(6)：692-704.

[5]代開勇.心血管系統(tǒng)鍵合圖模型仿真研究[D].杭州：浙江大學(xué)，2006. Dai Kaiyong. Study on Bound Graph Modeling and Simulation for Cardiovascular System[D].Hangzhou：Zhejiang University，2006.

[6]Boyers D G，Cuthbertson J G，Luetscher J A. Simulation of the Human Cardiovascular System：a Model with Normal Responses to Change of Posture，Blood Loss，Transfusion，and Autonomic Blockade[J]. Simulation，1972，18(6)：197-206.

[7]Guyton A C，Coleman T G. Quantitative Analysis of the Pathophysiology of Hypertension[J]. Circulation Research，1969，24(S5)：1-19.

[8]Melchior F M，Srinivasan R S，Charles J B. Mathematical Modeling of Human Cardiovascular System for Simulation of Orthostatic Response[J]. American Journal of Physiology-Heart and Circulatory Physiology，1992，262(6)：H1920-H1933.

[9]Melchior F M，Srinivasan R S，Thullier P H，et al. Simulation of Cardiovascular Response to Lower Body Negative Pressure from 0 to 40 mmHg[J]. Journal of Applied Physiology，1994，77(2)：630-640.

[10]Heldt T，Shim E B，Kamm R D，et al. Computational Modeling of Cardiovascular Response to Orthostatic Stress [J]. Journal of Applied Physiology，2002，92(3)：1239-1254.

[11]Heldt Thomas. Computational Models of Cardiovascular Response to Orthostatic Stress[D]. Cambridge：Massachusetts Institute of Technology，2004.

[12]Heldt T，Mark R G. Understanding Post-Spaceflight Orthostatic Intolerance：a Simulation Study[C]//In Computers in Cardiology 2005. Lyon：IEEE, 2005：631-634.

[13]Kappel F，F(xiàn)ink M，Batzel J J. Aspects of Control of the Cardiovascular-Respiratory System During Orthostatic Stress Induced by Lower Body Negative Pressure[J]. Mathematical Biosciences，2007，206(2)：273-308.

[14]Van Heusden K，Gisolf J，Stok W J，et al.Mathematical Modeling of Gravitational Effects on the Circulation：Importance of the Time Course of Venous Pooling and Blood Volume Changes in the Lungs[J]. American Journal of Physiology-Heart and Circulatory Physiology，2006，291 (5)：H2152-H2165.

[15]郝衛(wèi)亞，白凈，張立藩，等.失重后血量減少致立位應(yīng)激時心血管反應(yīng)改變的仿真研究[J].生物醫(yī)學(xué)工程學(xué)雜志，2002，19(1)：48-52. Hao Weiya，Bai Jing，Zhang Lifan，et al. The Changes of Cardiovascular Response to Orthostatic Stress Caused by Hypovolemia Induced by Weightlessness：A Simulation Study[J]. Journal of Biomedical Engineering，2002，19(1)：48-52.

[16]馮岱雅. 失重條件下人體心血管系統(tǒng)反應(yīng)的仿真研究[D].西安：第四軍醫(yī)大學(xué)，2005. Feng Daiya. A Simulation Study on the Responses of Cardiovascular System to the Weightlessness[D]. Xi’an：The Fourth Military Medical University，2005.

[17]劉洋，焦純，張國鵬，等.不同角度傾斜對人體心血管系統(tǒng)影響的仿真研究[J].航天醫(yī)學(xué)與醫(yī)學(xué)工程，2009，22(2)：89-93. Liu Yang，Jiao Chun，Zhang Guopeng，et al. Simulation Study on Cardiovascular Effect Induced with Tilting At Different Angles[J]. Space Medicine & Medical Engineering，2009，22(2)：89-93.

[18]Le Rolle V，Herná Ndez A I，Richard P Y，et al.An Autonomic Nervous System Model Applied to the Analysis of Orthostatic Tests[J]. Modelling and Simulation in Engineering, 2008 (2008)：2.

[19]Lim Einly，Chan G S H，Dokos S， et al. A Cardiovascular Mathematical Model of Graded Head-Up Tilt[J]. Plos One，2013，8(10).

[20]Liu Y，Zhang L F，Zhang K L，et al. Role of Carotid Baroreflex and Sympathetic Responses in the Push-Pull Effect：A Simulation Study[J]. Aviation，Space，and Environmental Medicine，2012，83(9)：841-849.

[21]Broskey J，Sharp M K. Evaluation of Mechanisms of Postflight Orthostatic Intolerance with a Simple Cardiovascular System Model[J]. Annals of Biomedical Engineering，2007，35(10)：1800-1811.

[22]Ursino M.Interaction Between Carotid Baroregulation and the Pulsating Heart：a Mathematical Model[J]. American Journal of Physiology-Heart and Circulatory Physiology，1998，275(5)：1733-1747.

（責(zé)任編輯：申劍）

Modeling and Simulation of Cardiovascular System under Gravitational Stress

Wang Hanqing，Zhao Jinping，Chen Juan
（School of civil Engineering, Hunan University of Technology，Zhuzhou Hunan 412007，China）

Cardiovascular system is affected by blood hydrostatic pressure change due to gravitational stress. To understand the mechanism of cardiovascular system under gravitational stress is helpful for countermeasures developing. Summarizes the cardiovascular system modeling methods, reviews the model development progress, expounds the research and application status of effects of gravity changes on cardiovascular system, introduces the authors' research achievements, and finally points out the research trend of cardiovascular circulation system modeling and simulation.

cardiovascular system；mathematical model；gravitational stress；modeling and simulation；thermal comfort

Q463；TU83

A

1673-9833(2014)03-0099-06

10.3969/j.issn.1673-9833.2014.03.020

2014-02-13

國家自然科學(xué)基金資助項目（51276057），湖南省研究生創(chuàng)新基金資助項目（CX2013B411）

王漢青（1963-），男，湖南益陽人，湖南工業(yè)大學(xué)教授，博士，主要從事室內(nèi)空氣品質(zhì)計算機(jī)仿真及建筑環(huán)境與設(shè)備節(jié)能技術(shù)方面的研究，E-mail：hqwang2011@126.com