数値サーマルマネキンの開発
体温調節機能（発汗調節、血流調節、ふるえ産熱）

呼吸による熱放散 Qres④

皮膚温度、発汗量

室内空気流動

体温調節モデル蒸発による熱放散 Esk

対流熱放散 Qcv③①

人体内部核層代謝熱生成(M-W)
皮膚表面との内部熱伝達(体温調節モデル)⑤②
壁面との放射による熱放散 Qr

人体－環境系の熱交換CFD

放射伝熱解析

対流放熱量

放射放熱量

数値サーマルマネキン

数値サ－マルマネキンは、人間形状の人体モデルを室内空間に取り込み、その環境下の人体における温熱生理応答、及び人体起源の室内温熱環境への影響を数値的に解明する総合数値解析手法である。数値サーマルマネキンを用いて、実験せずに室内温熱快適性状を予測・評価できるため、経済的な温熱環境設計ツールともなる。

1/4

加藤研究室・大岡研究室・菊本研究室
Kato Lab., Ooka Lab., and Kikumoto Lab.

被験者実験と数値解析の概要
0.67m

吸込口yes人 体？

放射パネル

no1.83m

C F D

放射パネル

9009001830

放射計算ルーチン

対流伝熱量 QCi

Uin：0.05m/s

QR i = å C jisT j -Aie isTi

実験室の概要実験概要

44q = ‐q i

4j =1

風速0.05m/s以下
0.915mQRi = åC jisTj -AieisTi
全熱流束
q i= (QRi +QCi )/Ai

QRi = QR ,in - QR ,out
4n

j =1

QRi + QCi + QWi + QGi = 0n

放射伝熱量0.915m

固体面温度分布

数値解析のモデル空間

Ti (i=1 ～ n)

体温調節モデ皮膚温度（Fanger中立）
ModelルtSmith’s= 36.4-0.054 qs

連成解析の流れ

健康な大学生男女6名ずつ計12名を被験者とした。実験時間は60分、50～60分の平均値を皮膚温の分析対象とする。

数値解析概要Cij ：全交換面積
T ： 固体面温度
σ ： ｽﾃﾌｧﾝ・ﾎﾞﾙﾂﾏﾝ定数
ε ： 放射率n ： 固体面分割数
QR ：正味放射伝熱量
Qc ： 対流伝熱量Qw ： 熱伝導伝熱量
QG ： 熱発生量
QR,in ： 吸収放射伝熱量
QR,out ： 射出放射伝熱量①②

空間の総メッシュ数は268,421、表面総メッシュ数は14,038。
乱流モデルは低Re数型k-εモデルを用い、差分スキームは速度項に対しMARSスキーム、温度項に対しで一次風上スキームを採用する。
③ Smith Modelの計算の設定に関し、代謝量を1.2Met、気温と吸気の温度を共に28℃に設定する。

検討ケース

Case 気温(℃) 相対湿度(％) 前パネル温度(℃) 後パネル温度(℃)
1502135250352132850282844014425404214

人体生理モデル

マッピング

Smith’s Model

2/4

人体形状
加藤研究室・大岡研究室・菊本研究室
Kato Lab., Ooka Lab., and Kikumoto Lab.

被験者実験と数値解析による皮膚温分布の比較
男子の解析結果

男子の実験結果腰(前)左大腿(上)
右大腿(上)(℃) 顔面
胸 38
36

背3432

左下腿(前)腰(後)左大腿(下)
右大腿(下)腰(前)左大腿(上)
右大腿(上)

左下腿(後)30
28
26

右下腿(前)

女子の実験結果

左下腿(前)

右下腿(後)

右下腿(前)

女子の解析結果(℃) 顔面
胸 38
36
34
32
30
28
26

背腰(後)左大腿(下)右大腿(下)
左下腿(後)

右下腿(後)

左上腕(前)

左上腕(後)

左上腕(前)

左上腕(後)

右上腕(前)

右上腕(後)

右上腕(前)

右上腕(後)

左前腕(前)
右前腕(前)
左足背
右足背

平均

右手背

左前腕(後)
右前腕(後)
左手背

左前腕(前)
右前腕(前)
左足背
右足背

Case4
腰(前)
左大腿(上)
右大腿(上)
左下腿(前)
右下腿(前)(℃) 顔面
胸 3836343230
2826

Case5
背

腰(後)左大腿(下)右大腿(下)
左下腿(後)
右下腿(後)

左上腕(前)

左上腕(後)

右上腕(前)

右上腕(後)

左前腕(前)
右前腕(前)
左足背
右足背

平均

平均

左前腕(後)
右前腕
左手背(後)
右手背

左前腕(後)
右前腕
左手背(後)
右手背差異が生じる原因の考察①

実験による各部位の皮膚温が該当部位のある点の温度により代表され、数値解析による皮膚温は該当部位の各表面メッシュの温度の面積重みつけ平均値となる。不均一環境では、各部位でも皮膚温が一定の分布となっているため、部位の代表点の温度と平均温度は必ずしも一致しない。
② 数値解析手法自体は、連成手法・Smith Model・着衣の扱い方などに様々な問題点がある。

新たな人体熱モデルの開発が必要

Case3

3/4

加藤研究室・大岡研究室・菊本研究室
Kato Lab., Ooka Lab., and Kikumoto Lab.

新体温調節モデルの開発 （Sakoiモデル）
概要
部位周辺の外界条件、部位内の熱産生、皮下脂肪による断熱、血流による耐熱の再配分、発汗、ふるえに着目して、Smith Modelに基づき開発された。3次元体温調節モデルの形状、組成にはSmith Modelを用いた。ただし、温度分布を反映した熱移動計算を行えるように、より微細な要素分割とする。
特徴✓全身の血液に関する質量保存を改良する。
血液体積が変化する血管拡張・収縮を行わない。血管体積一定の場合、血流量設定値のみを変化する。
✓体内での3次元温度分布を再現できる。
✓断熱性の高い皮下脂肪などの伝導熱量を、組成自体の熱伝道率と3次元形状、組成内温度分布自体に基づき計算する。
✓高温条件ほど表在静脈を介して戻る血液の割合が大きくなって血流の1方向性が強くなる傾向を表現できる。
皮膚組織から戻る静脈血のみは表在静脈を介して搬出され、それ以外の組織から戻る血液は深部の静脈を介して搬出される。また、断熱性の高い皮下脂肪を通過する皮膚血流のみを体温に応じて変動し、高温条件では皮膚血流が増加し、低温条件では皮膚血流が減少する。
℃[℃]
36.5

本モデル：T a = 25.5℃
Smithモデル：T a = 25.5℃
実測値：T a = 25.5℃35.9

手前腕35.3

上腕首34.7

脳
肺
内臓
骨（頭蓋）
骨（椎、脇、腰）骨（四肢）筋
脂肪皮膚足下腿37.0

局所皮膚温 Tsk i [℃]

頭大腿モデルの組成

本モデル：T a = 28.0℃
Smithモデル：T a = 28.0℃
実測値：T a = 28.0℃

34.1

35.0
33.0
31.0
29.0

椅座, 静穏,T a = T r (均一),
RH50%,トランクスのみ

27.0
頭33.5

気温 = 放射温 = 28.0℃, 椅座, トランクスのみ, RH50％, 静穏
図-17a

背腰前腰
後上腕前腕

手大腿

下腿

足

モデルによる計算結果と実測値の比較

気温 28.0℃均一における計算皮膚温モデルの計算例（皮膚温分布）
分布（椅座、気温=放射温、RH50％、胸●新体温調節モデルは産業技術総合研究所との共同研究として行われた。

静穏、トランクスのみ）

4/4

加藤研究室・大岡研究室・菊本研究室
Kato Lab., Ooka Lab., and Kikumoto Lab.