用于人-航天服仿真的人體熱模型發展研究

方明元，王 晶，李西園，畢研強

(1. 北京衛星環境工程研究所， 北京100094； 2. 可靠性與環境工程技術重點實驗室， 北京100094)

1 引言

人-航天服熱模型是一種將人體熱模型和航天服熱模型相耦合的數值仿真模型，能夠預測各種工況下航天服內溫度、濕度、CO2濃度等參數的變化情況，用于評估航天服的性能，預測航天員在穿著航天服進行艙外活動時的熱舒適度等指標，為航天服系統的性能分析提供便利，同時為航天服及其生命保障系統的研制提供參考。

航天服熱模型一般由液冷/通風服熱模型和便攜式生保系統熱模型組成，雖然實際結構復雜，但其傳熱機制相對簡單。 與此相對，人體作為生物體具有復雜的結構和生理調節機制，其建模難度要大于航天服。 同時人體作為一個熱源又與航天服內傳熱和傳質情況相耦合，因此是否準確建立人體熱模型將直接決定人-航天服模型的準確度。

自上世紀60 年代起，各航天機構均開展了面向航天服設計、優化的人體熱模型研究，為航天服、環控生保系統的研制提供支撐。 與建筑空調、載具、醫療等領域的人體熱模型相比，航天服用人體熱模型在人體的幾何構型、血液建模、個體差異影響等方面均存在明顯不同，并且其他領域的人體熱模型往往欠缺對航天環境低壓、失重因素的考慮。 本文整理近50 年來各機構研究的用于航天領域的人體熱模型，對比其他領域典型的人體熱模型，從幾何構型、血液建模、低壓失重、個體差異影響4 個方面給出未來航天服用人體熱模型的發展方向，為未來航天服的研制提供參考。

2 人體熱模型幾何精細度

自實施阿波羅計劃以來，美國開展了大量人-航天服的建模仿真研究。 從上世紀90 年代起，中國也開始相關研究。 早期大部分人體熱模型為一維10 節段模型，進入21 世紀后，人體熱模型的幾何精細度得到提升，主要體現在模型節段數量的增加和模型維數的提高。

對于模型節段的劃分，早期模型多基于一維10 節段的Stolwijk 模型[1]。 該模型將人體劃分成頭部、軀干、手臂、手、腿、腳等10 段，如圖1(a)所示。 隨著計算機性能的提升，相關軟件的發展以及人體生理學、解剖學數據的豐富，開始使用更細致的節段劃分人體熱模型。 如應用較為廣泛的基于Wissler 模型[2]和Fiala 模型[3]的15 節段人體熱模型。 相比于早期的10 節段模型，Wissler 模型對人的四肢和軀干進行了細致劃分(圖1(b))，將手臂分為上臂和下臂，腿部分為大腿和小腿，軀干分為胸部和腹部。 這樣的劃分不僅詳細地考慮了節段之間組織結構的差異，還考慮了人在做不同運動時不同部位肌肉產熱的差異。 而Fiala 模型則增加了面部、肩部、頸部和腰部節段(圖1(c))。 面部的組織分布以及對溫度的敏感程度與頭部其他地方有較明顯的區別，因此把面部從頭部節段分離出來。 其他增加的節段同理。 由此可以看出越細致的節段劃分使得人體熱模型所考慮的人體細節更加豐富，模型更能反映真實情況。

對于模型維數的選擇，早期模型多為一維模型，每一個節段上各層組織的物理參數沒有周向和軸向的差別(圖2(a))。 進入21 世紀之后，二維和三維模型開始發展起來，如Fiala 模型是典型的二維模型，該模型對每一個節段進行了徑向劃分，使得模型節段能夠反映周向溫度的變化(圖2(b))。 三維模型則一般通過對節段進行網格劃分來建立(圖2(c))，通過有限元或有限差分的數學方法求解，三維模型能夠比二維和一維模型獲得更詳細的溫度場分布，適應性更廣，尤其適用于溫度梯度較大的情況。 另外在某些特殊情況下，如航天服液冷服管路的分布設計需要三維的人體熱模型來反映人體各個節段周向和軸向的溫度分布差異。

圖1 人體熱模型節段劃分示意[1-3]Fig.1 Schematic diagram of segmentation of human thermal model[1-3]

圖2 不同維數的人體熱模型節段對比Fig.2 Comparison of segments with different dimensions

綜上所述，采用更多的節段劃分和三維構型是未來航天服熱仿真研究中的人體熱模型的發展趨勢。 需要注意的是，節段不是劃分得越多越好，一方面節段的增加會導致模型的結構變得復雜，增加建模難度和仿真時間；另一方面節段劃分越多，人體劃分越細致，需要測量的生理參數就越多，相鄰兩節段之間的生理差異也越小，整個模型預測精度的提升也越小。 此外，在工程實際應用中，還需要根據具體的熱邊界條件，如液冷服的分布等有針對性地對人體熱模型節段的劃分作出調整。 目前用于航天服仿真的人體熱模型普遍在15節段及15 節段以下，節段劃分標準也不統一。 綜合各種模型的劃分優點，可將比較流行的Wissler的模型、Fiala 的模型和Berkeley 模型[4]的劃分方法相結合，即在Fiala 模型的基礎上考慮上臂和下臂、大腿和小腿，并將胸部節段分成前胸和后背兩部分，構成20 節段的人體熱模型(圖1(d))。 對于手臂和腿部的分割，由于它們位于活動關節的兩側，一側為近心端，一側為遠心端，生理參數會有差異，組織溫度也不同，比如大腿溫度一般比小腿溫度要低。 對于航天服來說，液冷服管路中的液體和通風服中的氣體溫度，遠心端也會比近心端要低一些；對于前胸和后背的分割，除了考慮這兩部分對人體熱感覺的敏感程度不同，艙外航天服背部還安裝有便攜式生保系統，這也是液冷服液體的出口位置，因此艙外航天服內側前胸和后背的邊界條件也會略有不同。 綜上差異，將人體熱模型節段增加到20 節段可進一步提升模型精度，且不至于把模型變得過于復雜和冗余。

目前用于航天服仿真的三維人體熱模型較少，新模型的開發可參考其他領域的三維人體熱模型，如Shoji[5]開發的用于預測著裝人體生理反應的模型、趙陽[6]開發的用于暖通空調領域的模型、Ferreira[7]開發的通用人體模型、Sun[8]開發的用于非均勻熱環境的模型、Tang[9]開發的用于功能性服裝設計和舒適度評估的模型等。 將其他領域三維人體熱模型用于航天服仿真的難點在于人體熱模型與其邊界條件即航天服中液冷/通風服熱模型的耦合。 對于各節段需進行徑向扇面劃分的人體熱模型，液冷/通風服熱模型需要進行同樣的徑向劃分，使每一個服裝扇面與人體節段扇面相吻合；通過有限元方法建立的三維人體熱模型，液冷/通風服熱模型也要使用相同的方法進行三維建模。 李杰[10]建立了三維人-液冷/通風服熱模型來研究人與液冷/通風服的換熱情況，但模型中沒有考慮人體血液調節作用和人體與手套、頭盔的熱交換，三維人-航天服熱模型距離工程化應用還需要進一步完善。 另外，對于與液冷/通風服耦合的便攜生保系統模型，其輸入和輸出量一般為冷卻液流量、通風量等一維參數，因此其建模的方法和維數可以和人體熱模型、液冷/通風服熱模型不同，仍可保持其常用的一維集總參數結構。表1 列出了在航天服研究中使用的人體熱模型幾何結構。

3 血液循環系統換熱模型

人體的血液循環系統結構復雜，對其進行建模時需要合理簡化。 血液循環系統對人體內部的傳熱及熱調節過程有較大影響，有研究表明人體的血液循環轉移了人體50%～80%的熱量[31]，建模時不恰當的簡化會導致人體熱模型的瞬態響應結果偏離實際。 對于曾廣泛應用于航天領域的10 節段人體熱模型，研究人員為了簡化模型，對血液循環系統模型做了許多與解剖學不同的假設，例如假設心臟直接向所有節段供血，使得所有節段輸入的動脈血液溫度相同，與血液溫度沿肢端逐漸降低的事實相悖；假設血流速率和肌肉代謝率成正比[32]，使得人體熱模型在模擬寒冷環境出現寒顫時，血流速率會突然升高；假設動脈血液直接流向各個節段，使得模型各節段溫度會快速升高，甚至觸發出汗機制，這與實際結果存在較大偏差。 因此血液循環系統對模型的瞬態響應影響顯著。

近年來其他領域的人體熱模型已經使用了較為復雜的血液循環模型，大大提高了模型對瞬態環境的響應能力。 如UC Berkeley 模型、Fiala 模型、Salloum 模型[33]、JOS-2 模型[34]等，這些模型均模擬了動靜脈血液從心臟處依次流到各個人體節段的循環過程中的傳熱，并考慮了動靜脈逆流換熱、毛細血管與各組織的換熱(圖3)。 這些精細化考慮，保證了人體熱模型在瞬態條件下溫度預測不會明顯偏離實際。 在此基礎上，引入如JOS-2 模型中對人體熱模型的手、腳部構建的淺表靜脈和動靜脈吻合結構，使得模型中手部和腳部的溫度預測更加準確。 另外，可以考慮類似Salloum 模型中用Avolio[35]精確的生理數據構建真實尺寸的動脈樹結構，見圖4。 通過動脈樹可以計算每條動脈的特性阻抗、傳播常數、波反射系數和輸入阻抗，進而較為精確地預測每條動脈和靜脈中的血流速率，從而使得血液和組織內的傳熱計算更加準確。 但這種動脈樹模型仍存在一些問題，如各節段的血流量無法根據人體的熱狀態進行調整，動脈血液進入皮膚、肌肉和直接進入靜脈的血流量的比例還無法確定[31]。 對于研究人員來說，一方面可以對現有血液循環系統換熱模型所存在的問題進行改進，另一方面可將近年來生理學和解剖學領域涌現出的優秀血液循環系統換熱模型引入航天領域人體熱模型中，來提高模型的預測精度[31]。

表1 國內外人-航天服熱仿真模型結構Table 1 The structure of human thermal model used in space suit development

圖3 血液換熱示意 Fig.3 Schematic diagram of blood heat transfer

4 低壓和失重對人體的影響

圖4 Salloum 模型中的動脈樹結構Fig.4 Structure of arterial tree in Salloum model

目前絕大多數人體熱模型只適用于地面正常壓力和重力情況，人-航天服熱模型需要考慮低壓、失重對人體的影響。

4.1 低壓

中國“飛天” 系列航天服內的氣壓約為40 kPa，氧氣含量不低于95%[36]。 在該低壓環境下，液體的沸點降低，人體蒸發散熱量增加，進而使得人體皮膚的散熱及呼吸散熱增加；由于空氣密度改變，人體與環境之間的換熱系數也發生變化。 崔代秀等[37]的實驗表明，23 ℃環境下壓力從1 atm 降到0.47 atm，人體的汗液蒸發量增加了17%。 低氣壓下人體血液氧分壓降低，會引起呼吸系統和心血管系統的活動增加，使得人體代謝率增加。 童力[38]的研究表明，當大氣壓力從1 atm 降到0.8 atm 時，靜坐的人體代謝率增加22.2%。

低壓環境下的人體蒸發量可以通過改變模型中呼吸散熱和皮膚汗液蒸發散熱計算中的環境壓力項來計算；對于人體與環境的各項換熱系數的改變，也可以從修改模型中的相關公式入手。 對于低壓對人體呼吸和代謝率的影響，可通過人體實驗獲得。 低壓環境的人體實驗一般在低壓艙中進行，低壓艙的空調系統調節艙內的溫度和濕度，受試者在低壓艙中處于靜坐或運動狀態，通過皮膚上的溫度傳感器記錄受試者的皮膚溫度變化情況，通過心肺功能測試系統測量人的呼吸速率、呼吸氣體中的氧分壓和二氧化碳分壓等多項指標，獲得人體的呼吸速率、代謝率等指標[38]。

4.2 失重

失重環境下，人的生理情況會發生多方面的變化：①心臟供血能力下降，血液流速變慢，血管內徑變小，血液管壁收縮能力下降，和平號空間站在軌實驗表明6 個月的空間飛行使航天員的心率下降了20%[39]，SLS-2 火箭在軌數據表明航天員的血液循環總量在36 h 內降低了20%[40]，導致人體的熱調節能力下降；②血流量在四肢的分配發生變化，下肢血壓降低，上肢、頭部血壓升高[39]，導致人體上肢和頭部的皮膚層和核心層的溫度升高，下肢的皮膚層和核心層溫度降低，禮炮6 號數據表明在軌航天員的胸部溫度上升3 ℃，小腿溫度下降3 ℃[41]；③在軌飛行中，多種因素使得航天員的攝入能量減少，人體長時間生活在狹小的航天器內，肌肉活動受到限制，失重環境使得支撐人體所需的力量減少，人體的代謝率會降低[40]；④由于下肢不再需要力量來支撐人體，人體運動模式發生了顯著改變，人體在工作時的代謝率分布會發生變化；⑤由于代謝率的降低，進一步影響人體對環境的散熱量，人的蒸發散熱、呼吸散熱均會下降；⑥失重影響人的心理狀態，人在失重環境中容易出現抑郁、焦慮、緊張等負面情緒，這將降低人體的熱舒適程度。 綜上所述，人體熱模型建模時需要考慮失重的影響。

失重對人體影響的定量分析需要通過人體實驗來實現。 由于在軌實驗費用昂貴，失重的生理學研究一般采用大量地面實驗和少量在軌實驗相結合的方法。 在地面實驗中，一般采用臥床方式來模擬人體在失重情況下的生理變化情況[42]，當臥床角度為頭低位-6°時，人體的心血管、肌肉等組織的工作情況與失重時最為相似。 另外也可使用大鼠尾吊法模擬微重力，通過測量大鼠生理指標變化探究失重對生物體的影響機理，從而推測失重環境下人體的生理指標變化情況[36]。

目前中國在失重條件下人體生理學研究還不能滿足相關工程的需求，失重對人體生理指標的影響尚未完全量化，失重環境下人在各種工作狀態下的代謝率、血流速率等數據較少，需要通過實驗來獲取相關數據。 考慮到航天員常常同時處于低壓和失重環境，2 種環境對人的影響是互相耦合的，因此需要將低壓和失重環境實驗相結合來獲得更準確的數據。 在地面模擬實驗中，可以在低壓艙中營造低壓純氧環境，充分做好防火措施，受試者采用頭低位臥床姿勢模擬失重環境，空調系統控制不同的溫度環境，通過貼在受試者皮膚上的溫度傳感器獲取人體溫度變化數據，通過心肺測試系統獲取人體呼吸及代謝率相關數據。 未來也可以通過在軌實驗獲得更為真實可信的數據以修正地面模擬實驗。

5 個體差異

一般來說，人體的年齡、性別、健康狀況、脂肪含量、BMI(Body Mass Index)指數、BMR(Basal Metabolic Rate)指數都會顯著影響人體的生理參數，進而影響人體的溫度分布和熱舒適度[43]。 研究表明個體差異對人體熱模型平均皮膚溫度的影響可達1.2 ℃[44]。 在現階段，航天員一般來自于層層選拔、經過訓練的空軍飛行員，體質優于普通人。 未來會有更多專業人員執行航天任務，如載荷專家、工程師等，他們的體質更接近于一般人，因此在未來航天服設計中，應當考慮針對種群、年齡、性別等引起的個體差異，通過重要特征參數化及參數特征化的迭代與回歸形成適用于不同群體的人體熱模型，以提升預測精度。

人體個性化生理參數的獲得方法可分為2類：間接估算和直接測量。 在有實驗條件的情況下，對航天員的生理數據進行直接測量能夠獲得更加準確的結果；但某些參數，如熱導率、血液灌注率和組織代謝率等測量比較困難[45]，在條件限制的情況下，可采用間接估算的方法，通過測量人體一些容易獲得的參數來估算那些不易測量的參數。 許多人體熱模型都通過估算的方法來考慮個體生理參數差異影響，如Berkeley 模型、Fiala 模型[46]、ThermoSEM 模型[47]等，但這幾種模型考慮的參數都不夠全面，Davoodi 的人體模型[48]較為全面地考慮了各種個性化參數，使用身高、體重、性別、年齡等容易獲得的參數來估算其他不易獲得的參數，如皮膚面積、各層組織的厚度、組織的熱容和熱阻、皮膚血流量、代謝率以及人體在熱調節過程中的血流變化、出汗率和寒顫產熱量等，計算公式由實驗數據回歸分析獲得。 這種通過簡單易測的數據來估算難測的生理參數的方法值得借鑒。 不過由于Davoodi 人體模型使用的是2 節點模型，并且各參數的估算方法基于外國人的體質，無法直接用于中國的多節點人體熱模型。 Zhou等[44]基于Fiala 模型建立了適用中國人體質的人體模型，模型在皮膚面積、脂肪含量等生理參數上按中國人的體質進行了調整，其數據和方法可作為參考。 由于缺少基于中國人體質的熱調節系統的相關生理數據，模型中有關人體熱調節系統的參數如寒顫產熱量、出汗率以及人體對環境溫度變化的響應特征等并未進行個性化修正。

綜上，為中國航天員建立覆蓋更多人群的個性化人體熱模型關鍵在于獲取他們的人體生理數據。 對于各節段的皮膚面積、脂肪含量等生理解剖方面的數據可通過直接測量與估算相結合的方式獲取；對于人體對環境變化做出的響應差異，如出汗臨界溫度、寒顫產熱量等需要通過實驗來測量；對于人體主動熱調節系統中的各項控制系數可以通過對皮膚溫度在穩態和瞬態溫度環境中的響應數據進行回歸分析來反解算獲得。 隨著中國載人航天活動逐漸頻繁，航天員在軌時間逐漸增加，研究人員可考慮通過在軌實驗采集航天員在低壓失重環境下的個性化生理數據，在軌實驗數據相比地面模擬實驗數據更加真實可信，并通過地面模擬實驗數據與在軌數據對比，提出地面實驗方法及數據處理的修正方案，以此來完善地面模擬實驗。

6 結論

1)用于人-航天服仿真系統中的人體熱模型的幾何結構應向更加精細的方向發展，以三維20節段的構型為宜。 三維結構使得模型的應用范圍更廣，20 節段的劃分形式使得模型在溫度場計算精度與計算復雜度之間獲得平衡。

2)引入生理學、解剖學領域的先進血液循環系統模型可以精確模擬血管與組織間的換熱過程，改善人體熱調節模型的瞬態響應。

3)低壓失重環境對人體影響的數據較為缺乏，需要進一步開展相應的生理學實驗獲取數據來提高人體熱模型在低壓失重環境下的仿真精度。

4)考慮到個體差異的影響，應為中國航天員建立覆蓋更多人群的個性化生理數據庫，經過個體參數優化的人體熱模型，預測結果才能更貼近航天員的實際生理情況。

參考文獻(References）

[ 1] Stolwijk J A J. A mathematical model of physiological temperature regulation[R]. NASA-CR-1855， 1971.

[ 2] Wissler E. Mathematical simulation of human thermal behavior using whole body models[J]. Heat Transfer in Medicine and Biology， 1985， 1(13)： 325-373.

[ 3] Fiala D， Lomas K J， Stohrer M. A computer model of human thermoregulation for a wide range of environmental conditions：the passive system[J]. Journal of Applied Physiology， 1999，87(5)： 72-1957.

[ 4] Huizenga C， Hui Z， Arens E. A model of human physiology and comfort for assessing complex thermal environments[J].Building and Environment， 2001， 36(6)： 691-699.

[ 5] Shoji Y， Takemori T， Matsunami H， et al. A three-dimensional，clothed human thermal model[J]. Heat Transfer-Japanese Research， 1997， 26(8)： 554-566.

[ 6] 趙陽. 基于Smith 模型的三維有限元人體熱調節模型[D].大連： 大連理工大學， 2009.

Zhao Y. A Finite Element Based Transient Three-dimensional Model of the Human Thermal System Evaluated from Smith Model[D]. Dalian： Dalian University of Technology， 2009.(in Chinese)

[ 7] Ferreira M S， Yanagihara J I. A transient three-dimensional heat transfer model of the human body[J]. International Communications in Heat and Mass Transfer， 2009， 36(7)： 718-724.

[ 8] Sun X， Eckels S， Zheng Z C. An improved thermal model of the human body[J]. HVAC ＆ R Research， 2012， 18(3)：323-338.

[ 9] Tang Y， He Y， Shao H， et al. Assessment of comfortable clothing thermal resistance using a multi-scale human thermoregulatory model[J]. International Journal of Heat and Mass Transfer， 2016， 98： 568-583.

[10] 李杰. 艙外航天服—人體熱耦合數值模擬[D]. 南京： 南京航空航天大學， 2017.

Li J. Numerical Simulation of Coupling Heat Transfer between/in the Human Body and the Spacesuit for Extravehicular Activity[D]. Nanjing： Nanjing University of Aeronautics and Astronautics， 2017. (in Chinese)

[11] Kuznetz L H. Automatic control of human thermal comfort by a liquid-cooled garment[J]. Journal of Biomechanical Engineering， 1980， 102(2)： 155-161.

[12] Kuznetz L H. Control of Thermal Balance by a Liquid Circulating Garment Based on a Mathematical Representation of the Human Thermoregulatory System.[D]. Berkeley： University of California， 1976.

[13] Wissler E H. Simulation of fluid-cooled or heated garments that allow man to function in hostile environments[J]. Chemical Engineering Science， 1986， 41(6)： 1689-1698.

[14] 袁修干，沙斌，王俊清. 液冷式個體熱調節的二維體溫調節數學模型[J]. 航空學報， 1989， 10(5)： 242-248.

Yuan X G， Sha B， Wang J Q. A two-dimensional mathematical model of human thermoregulation for personal thermal conditioning with water cooling[J]. Journal of Aeronautics，1989， 10(5)： 242-248. (in Chinese)

[15] Bue G C， Conger B C， Iovine J V， et al. Asda-advanced suit design analyzer computer program[J]. SAE Transactions，1992， 101： 1398-1407.

[16] Zaratuichenko I，Kuznetz E，Filipenkov S. Human thermoregulation model for space suit： mathematical model for human thermal regulation in a suited mode with ventilation and liquid cooling capabilities provided[J]. SAE Transactions， 1997，106(1996)： 567-578.

[17] Campbell A B， French J D， Nair S S， et al. Thermal analysis and design of an advanced space suit[J]. Journal of Thermophysics and Heat Transfer， 2000， 14(2)： 151-160.

[18] Campbell A B， French J D， Nair S S， et al. Dynamic modeling of the minimum consumables plss[R]. SAE 1999-01-1999， 1999.

[19] Mays D C， French J， Nair S S， et al. Design of a transient thermal model of the cryogenic plss[R]. SAE 1999-01-2000，1999.

[20] 邱義芬，袁修干，梅志光，等. 艙外航天服熱設計技術[J].北京航空航天大學學報， 2003， 29(7)： 616-619.

Qiu Y F，Yuan X G，Mei Z G，et al. Thermal design technology of space suit used for extravehicular activity[J]. Journal of Beijing University of Aeronautics and Astronautics， 2003， 29(7)： 616-619. (in Chinese)

[21] 邱義芬，袁修干，梅志光，等. 艙外航天液冷服傳熱分析[J]. 航天醫學與醫學工程， 2001， 14(5)： 364-367.

大賽的題目都是面對當前的實際問題，結合最前沿的技術和主流技術的發展趨勢，對軟件進行設計和開發，學生在實際參賽過程中，會通過各種手段去學習最新最前沿的技術方法，對問題進行解決，從而進一步提升了學生學習新技術的能力[4]。

Qiu Y F， Yuan X G， Mei Z G， et al， Heat transfer analysis of liquid cooling garment used for extravehicular activity[J].Space Medicine ＆ Medical Engineering， 2001，14(5)：364-367. (in Chinese)

[22] 邱義芬，陳擁華，袁修干. 基于人體模型的艙外航天服熱系統分析與計算[J]. 北京航空航天大學學報， 2006， 32(2)： 167-170.

Qiu Y F，Chen Y H，Yuan X G. Thermal system analysis and calculation for extravehicular spacesuit based on human body model[J]. Journal of Beijing University of Aeronautics and Astronautics， 2006， 32(2)： 167-170. (in Chinese)

[23] Matjushev T V， Khabarovsky N N， Kurmazenko E A， et al.Simulation model of the human body thermoregulation system and its applications for design of air/space autonomous life support systems[R]. SAE 2001-01-2097， 2001.

[24] Nyberg K L， Diller K R， Wissler E H. Model of human/liquid cooling garment interaction for space suit automatic thermal control[J]. Journal of Biomechanical Engineering，2001，123(1)： 114-120.

[25] Thornton S， Xu W， Nair S. Thermal comfort control of an advanced space suit[R]. SAE 2001-01-2268， 2001.

[26] 余敏賢，林貴平，王亮. 艙外航天服熱控系統仿真[J]. 空間科學學報， 2004， 24(5)： 373-381.

Yu M X， Lin G P， Wang L. A simulation of EVA space suit thermal control system[J]. Chinese Journal of Space Science，2004， 24(5)： 373-381. (in Chinese)

[28] 王晶，袁衛星，袁修干. 艙外航天服便攜式生命保障系統研究進展[J]. 航天醫學與醫學工程， 2009， 22(1)： 67-71.

Wang J，Yuan W X，Yuan X G. Research progress of portable life support system for extravehicular activity space suit[J].Space Medicine ＆ Medical Engineering， 2009， 22(1)： 67-71. (in Chinese)

[29] Westin J K， Kapat J S， Chow L C. An improved thermoregulatory model for automatic cooling control development in liquid cooling garment systems[J]. Journal of Thermal Science and Engineering Applications， 2010， 2(3)： 3-14.

[30] G?ser J， Olthoff C T. Results of a dynamic liquid cooling garment simulation in v-suit[C]/ /44thInternational Conference on Environmental Systems， Tucson， Arizona： Research gate，2014： 213-234.

[31] 趙陽，端木琳，李祥立，等. 數值人體模型的建立方法及其研究發展綜述[J]. 建筑熱能通風空調，2010，29(1)：20-25.

Zhao Y， Duanmu L， LI X L， et al. Literature review and development of the numerical human thermal model[J]. Building Energy ＆ Environment， 2010， 29(1)： 20-25. (in Chinese)

[32] Smith L， Nair S， Miles J， et al. Human thermal model evaluation for portable life support system control development[C]/ /Space Programs and Technologies Conference and Exhibit， Washington D.C.： AIAA， 1994： 4621.

[33] Salloum M， Ghaddar N， Ghali K. A new transient bioheat model of the human body and its integration to clothing models[J]. International Journal of Thermal Sciences， 2007， 46(4)： 371-384.

[34] Kobayashi Y，Tanabe S. Development of JOS-2 human thermoregulation model with detailed vascular system[J]. Building and Environment， 2013， 66： 1-10.

[35] Avolio A. Multi-branched model of the human arterial system[J]. Medical and Biological Engineering and Computing，1980， 18(6)： 709-718.

[36] 李潭秋. “飛天”艙外航天服的研制[J]. 載人航天， 2008，14(4)： 8-18.

Li T Q. Development of Feitian EVA space suit[J]. Manned Spaceflight， 2008， 14(4)： 8-18. (in Chinese)

[37] 崔代秀，祁章年，王憲章，等. 座艙溫、濕度及壓力改變對人體同環境熱交換的影響[J]. 航天醫學與醫學工程，1996， 9(1)： 42-48.

Cui D X， Qi Z N， Wang X Z， et al. Effect of temperature，humidity， pressure in cabin on heat exchange of human with environment[J]. Space Medicine ＆ Medical Engineering，1996，9(1)：42-48.(in Chinese)

[38] 童力. 基于局部換熱影響的低壓人體熱調節模型研究[D]. 青島： 青島理工大學， 2014.

Tong L. Study on Human Thermoregulation Model based on Local Heat Exchange in Low Pressure Environment [D].Qingdao： Qingdao Technological University， 2014. (in Chinese)

[39] Zhu H， Wang H， Liu Z. Effects of real and simulated weightlessness on the cardiac and peripheral vascular functions of humans： a review[J]. International Journal of Occupational Medicine ＆ Environmental Health， 2015， 28(5)： 793-802.

[40] 段波. 微重力下熱濕環境對人體熱舒適的影響研究[D].重慶： 重慶大學， 2017.

Duan B，Study on the Effect of Microgravity Thermal Environment on Human Thermal Comfort[D]. Chongqing： Chongqing University， 2017. (in Chinese)

[41] Novak L，Ulicny B，Wilhelm Z，et al. The simulation of thermal microclimate in the garment similar to those observed in the weightlessness[J]. Physiologist， 1988， 31(1)： 38-39.

[42] Greenleaf J E. Physiological responses to prolonged bed rest and fluid immersion in humans[J]. Journal of Applied Physiology， 1984， 57(3)： 619-633.

[43] van Ooijen A M J， van Marken Lichtenbelt W D， Westerterp K R. Individual differences in body temperature and the relation to energy expenditure： The influence of mild cold[J].Journal of Thermal Biology， 2001， 26(4)： 455-459.

[44] Zhou X， Lian Z， Lan L. An individualized human thermoregulation model for Chinese adults[J]. Building and Environment， 2013， 70： 257-265.

[45] 李世剛. 人體參數差異及其對人體熱舒適模型影響的研究[D]. 上海： 上海交通大學， 2012.

Li S G. Research for Differences of Human Body Parameters and Effect on Human Thermal Comfort Model[D]. Shanghai：Shanghai Jiao Tong University， 2012. (in Chinese)

[46] Fiala D，Havenith G. The Mechanobiology and Mechanophysiology of Military-related Injuries[M]. Springer， 2015： 265-302.

[47] van Marken Lichtenbelt W D， Frijns A J H， Fiala D， et al.Effect of individual characteristics on a mathematical model of human thermoregulation[J]. Journal of Thermal Biology，2004， 29(7)： 577-581.

[48] Davoodi F， Hassanzadeh H， Zolfaghari S A， et al. A new individualized thermoregulatory bio-heat model for evaluating the effects of personal characteristics on human body thermal response[J]. Building and Environment， 2018， 136： 62-76.