JPH0518592A - 温熱感覚演算方法および装置、予測平均温感演算方法および装置 - Google Patents

温熱感覚演算方法および装置、予測平均温感演算方法および装置

Info

Publication number
JPH0518592A
JPH0518592A JP3254305A JP25430591A JPH0518592A JP H0518592 A JPH0518592 A JP H0518592A JP 3254305 A JP3254305 A JP 3254305A JP 25430591 A JP25430591 A JP 25430591A JP H0518592 A JPH0518592 A JP H0518592A
Authority
JP
Japan
Prior art keywords
temperature
teq
set temperature
temperature value
coefficient
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP3254305A
Other languages
English (en)
Other versions
JP2581625B2 (ja
Inventor
Akihiko Kon
昭彦 昆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Azbil Corp
Original Assignee
Azbil Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Azbil Corp filed Critical Azbil Corp
Publication of JPH0518592A publication Critical patent/JPH0518592A/ja
Application granted granted Critical
Publication of JP2581625B2 publication Critical patent/JP2581625B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/17—Catathermometers for measuring "cooling value" related either to weather conditions or to comfort of other human environment

Landscapes

  • Environmental & Geological Engineering (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Air Conditioning Control Device (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)
  • Radiation Pyrometers (AREA)

Abstract

(57)【要約】 【目的】 気流速度Ｖair が大きい場合でも、またセン
サ温度Ｔcrが設定温度値θ(th)よりも高い場合でも、等
価温度Ｔeq* の演算値を人体が感じる等価温度Ｔeqに高
精度に合致させる。 【構成】 気温Ｔａと衣服の熱抵抗Ｉｃｌとに基づき設
定温度値θ(th)を演算し、このθ(th)とセンサ温度Ｔcr
とが等しくなるように、Ｔcrを調整し得る加熱手段（モ
ジュール本体の内部に配置された加熱ヒータ）へ熱エネ
ルギＨθ(th)を供与し、Ｔcrがθ(th)よりも高い場合、
Ｔａ，θ(th)，Ｔcrとθ(th)との差ΔＴcr，Ｖair ，Ｉ
ｃｌを含んで表現される係数ｃ1 ，係数ｃ2 〜ｃ6 およ
び係数ｎに基づき、下記（ｂ）式によりＴeq* を演算す
る。 Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ4 ×（１
＋ｃ5 ×Ｖairn）×ΔＴcr−ｃ6 ×Ｈθ(th) ・・・
（ｂ）

Description

【発明の詳細な説明】
【０００１】
【産業上の利用分野】この発明は、人体が感じる温熱感
覚としての等価温度Ｔeq（Equivalent Temperature）を
演算する温熱感覚演算方法および装置、並びに、室内環
境の快適さを示す予測平均温感ＰＭＶ（Predicted Mean
Vote ）を演算する予測平均温感演算方法および装置に
関するものである。
【０００２】
【従来の技術】人体が感じる温熱環境の評価指数、すな
わち温熱感覚として、等価温度Ｔeqなるものがある。こ
の等価温度Ｔeqを求めるための第１の方法として、輻射
温度Ｔｒ，気温Ｔａ，気流速度Ｖair を個別のセンサで
計測し、所定の演算を施すことによって、人体の等価温
度Ｔeqにほゞ等しい等価温度Ｔeq* を求める方法が提案
されている。この方法は、演算処理が煩わしく、処理時
間が長引くという問題がある。そこで、最近、等価温度
Ｔeqを求めるための第２の方法として、モジュール本体
に加熱ヒータを組み込み、そのモジュール本体の温度
（センサ温度）Ｔcrを常に一定の値（例えば、３６.5
℃）に保つべく、加熱ヒータへの供給電力量を制御する
方法が提案されている。この方法によれば、加熱ヒータ
への供給電力量を計測し、所定の演算を施すことによっ
て、等価温度Ｔeq* を算出することができる。一方、室
内環境の快適さを示す予測平均温感ＰＭＶの演算式とし
て、従来より、ＩＳＯ７７３０に規定されたＰＭＶ方程
式、すなわち下記（１）式が採用されている。
【０００３】
【数１】
【０００４】すなわち、上記（１）式において、各種パ
ラメータを計測して代入すれば、自動的に予測平均温感
ＰＭＶを算出することができる。しかし、現実にビル室
内で各種パラメータの全てを計測することは難しく、ま
たＰＭＶの関係式それ自体が複雑すぎる。このため、本
出願人は、予測平均温感ＰＭＶに精度良く合致するもの
として、上述した等価温度Ｔeq* を使用し、簡単な演算
式を用いて、予測平均温感ＰＭＶ* を求める手法を先に
提案した。
【０００５】
【発明が解決しようとする課題】しかしながら、等価温
度Ｔeqを求めるための上述した第２の方法では、当然で
はあるが、モジュール本体と加熱ヒータとからなる環境
計測部が人体よりはるかに小さいので、その気流感度が
人体の気流感度よりも大きな値となる。この気流感度の
相異は、気流速度Ｖair が0.1 ｍ／ｓ前後の微風速領域
では一応正しく等価温度Ｔeqを計測できるが、より大き
い気流速度の場合は等価温度Ｔeqより等価温度Ｔeq* の
ずれが大きくなり、正確な計測ができなくなるという問
題が生じることを意味している。すなわち、気流速度Ｖ
air が少しでも大きくなると、気流感度が大きいため、
求めるべき等価温度Ｔeqより等価温度Ｔeq* が大幅に低
い値を示すものとなる。つまり、上述した第２の方法で
は、輻射温度Ｔｒ，気温Ｔａ，気流速度Ｖairを一体で
計測することができるという利点があるにも拘わらず、
すなわち演算処理が簡略化されるという利点があるにも
拘わらず、気流速度Ｖair が大きい場合には正確な計測
ができないという致命的な欠点がある。そして、この等
価温度Ｔeqと等価温度Ｔeq* とのずれは、予測平均温感
ＰＭＶ* の演算に際しても影響し、気流速度Ｖair が大
きい場合、正確な予測平均温感ＰＭＶ* を求めることが
できなという不都合を生じさせる。
【０００６】
【課題を解決するための手段】本発明はこのような課題
を解決するために提案されたもので、その第１発明（請
求項１に係る発明），第４発明（請求項４に係る発明）
では、気温Ｔａと衣服の熱抵抗Ｉｃｌとに基づき設定温
度値θ(th)を演算し、この設定温度値θ(th)とセンサ温
度Ｔcrとが等しくなるように、センサ温度Ｔcrを調整し
得る加熱手段へ熱エネルギＨθ(th)を供与し、センサ温
度Ｔcrが設定温度値θ(th)よりも高い場合、気温Ｔａ，
設定温度値θ(th)，センサ温度Ｔcrと設定温度値θ(th)
との差ΔＴcr，気流速度Ｖair ，衣服の熱抵抗Ｉｃｌを
含んで表現される係数ｃ1 ，係数ｃ2 〜ｃ6 および係数
ｎに基づき、下記（ｂ）式により等価温度Ｔeq* を演算
するようにしている。 Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ4 ×（１＋ｃ5 ×Ｖairn）×ΔＴ cr −ｃ6 ×Ｈθ(th) ・・・（ｂ） また、その第２発明（請求項２に係る発明），第５発明
（請求項５に係る発明）では、気温Ｔａと衣服の熱抵抗
Ｉｃｌとに基づき設定温度値θ(th)を演算し、この設定
温度値θ(th)とセンサ温度Ｔcrとが等しくなるように、
センサ温度Ｔcrを調整し得る加熱手段へ熱エネルギＨθ
(th)を供与し、センサ温度Ｔcrが設定温度値θ(th)より
も低い場合、気温Ｔａ，設定温度値θ(th)，熱エネルギ
Ｈθ(th)および係数ｂ1 〜ｂ4 に基づき、下記（ａ）式
により等価温度Ｔeq* を演算する一方、センサ温度Ｔcr
が設定温度値θ(th)よりも高い場合、気温Ｔａ，設定温
度値θ(th)，センサ温度Ｔcrと設定温度値θ(th)との差
ΔＴcr，気流速度Ｖair ，衣服の熱抵抗Ｉｃｌを含んで
表現される係数ｃ1 ，係数ｃ2 〜ｃ6 および係数ｎに基
づき、下記（ｂ）式により等価温度Ｔeq* を演算するよ
うにしている。 Ｔeq* ＝ｂ1 ＋ｂ2 ×θ(th)＋ｂ3 ×Ｔａ−ｂ4 ×Ｈθ(th) ・・・（ａ） Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ4 ×（１＋ｃ5 ×Ｖairn）×ΔＴ cr −ｃ6 ×Ｈθ(th) ・・・（ｂ） また、その第３発明（請求項３に係る発明），第７発明
（請求項７に係る発明）では、気温Ｔａと衣服の熱抵抗
Ｉｃｌとに基づき設定温度値θ(th)を演算し、この設定
温度値θ(th)とセンサ温度Ｔcrとが等しくなるように、
センサ温度Ｔcrを調整し得る加熱手段へ熱エネルギＨθ
(th)を供与し、センサ温度Ｔcrが設定温度値θ(th)より
も高い場合、気温Ｔａ，設定温度値θ(th)，センサ温度
Ｔcrと設定温度値θ(th)との差ΔＴcr，気流速度Ｖair
，衣服の熱抵抗Ｉｃｌを含んで表現される係数ｃ1 ，
係数ｃ2 〜ｃ6 および係数ｎに基づき、下記（ｂ）式に
より等価温度Ｔeq* を演算するものとし、これにより得
た等価温度Ｔeq* ，湿度ＲＨ，活動量Ｍet，皮膚温Ｔsk
（Ｔsk＝35.7−0.0275×Ｍet×58.15 ），気温Ｔａ，気
温Ｔａにおける絶対水蒸気圧Ｐ（Ｔａ），衣服の熱抵抗
Ｉｃｌ，係数Ｅ１〜Ｅ３に基づく熱抵抗の効果Ｆｃｌ
（Ｆｃｌ＝Ｅ１／（Ｅ２＋Ｅ３×Ｉｃｌ））および係数
Ｄ1 〜Ｄ7 に基づき、下記（ｃ）式により予測平均温感
ＰＭＶ* を演算するようにしている。 Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ4 ×（１＋ｃ5 ×Ｖairn）×ΔＴ cr −ｃ6 ×Ｈθ(th) ・・・（ｂ） ＰＭＶ* ＝α×ＬＯＡＤ* ・・・（ｃ） 但し、（ｃ）式において、 α＝0.303 ×ＥＸＰ（−0.036 ×Ｍet×58.15 ）＋0.028 Ｍet＜1.0 のとき、 ＬＯＡＤ* ＝Ｄ1 ＋Ｄ2 ×Ｍet−Ｄ3 ×Ｔsk＋Ｄ4 ×Ｍet×Ｔａ＋（Ｄ5 ＋Ｄ6 ×Ｍet）×ＲＨ×Ｐ（Ｔａ）−Ｆｃｌ×（Ｔsk−Ｔeq* ） Ｍet≧1.0 のとき、 ＬＯＡＤ* ＝Ｄ1 ＋Ｄ2 ×Ｍet−Ｄ3 ×Ｔsk＋Ｄ4 ×Ｍet×Ｔａ＋（Ｄ5 ＋Ｄ6 ×Ｍet）×ＲＨ×Ｐ（Ｔａ）−Ｄ7 ×（Ｍet−１）−Ｆｃｌ×（Ｔsk−Ｔeq* ） また、その第４発明（請求項４に係る発明），第８発明
（請求項８に係る発明）では、気温Ｔａと衣服の熱抵抗
Ｉｃｌとに基づき設定温度値θ(th)を演算し、この設定
温度値θ(th)とセンサ温度Ｔcrとが等しくなるように、
センサ温度Ｔcrを調整し得る加熱手段へ熱エネルギＨθ
(th)を供与し、センサ温度Ｔcrが設定温度値θ(th)より
も低い場合、気温Ｔａ，設定温度値θ(th)，熱エネルギ
Ｈθ(th)および係数ｂ1 〜ｂ4 に基づき、下記（ａ）式
により等価温度Ｔeq* を演算する一方、センサ温度Ｔcr
が設定温度値θ(th)よりも高い場合、気温Ｔａ，設定温
度値θ(th)，センサ温度Ｔcrと設定温度値θ(th)との差
ΔＴcr，気流速度Ｖair ，衣服の熱抵抗Ｉｃｌを含んで
表現される係数ｃ1 ，係数ｃ2 〜ｃ6 および係数ｎに基
づき、下記（ｂ）式により等価温度Ｔeq* を演算するも
のとし、これにより得た等価温度Ｔeq* ，湿度ＲＨ，活
動量Ｍet，皮膚温Ｔsk（Ｔsk＝35.7−0.0275×Ｍet×5
8.15 ），気温Ｔａ，気温Ｔａにおける絶対水蒸気圧Ｐ
（Ｔａ），衣服の熱抵抗Ｉｃｌ，係数Ｅ１〜Ｅ３に基づ
く熱抵抗の効果Ｆｃｌ（Ｆｃｌ＝Ｅ１／（Ｅ２＋Ｅ３×
Ｉｃｌ））および係数Ｄ1 〜Ｄ7 に基づき、下記（ｃ）
式により予測平均温感ＰＭＶ* を演算するようにしてい
る。 Ｔeq* ＝ｂ1 ＋ｂ2 ×θ(th)＋ｂ3 ×Ｔａ−ｂ4 ×Ｈθ(th) ・・・（ａ） Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ4 ×（１＋ｃ5 ×Ｖairn）×ΔＴ cr −ｃ6 ×Ｈθ(th) ・・・（ｂ） ＰＭＶ* ＝α×ＬＯＡＤ* ・・・（ｃ） 但し、（ｃ）式において、 α＝0.303 ×ＥＸＰ（−0.036 ×Ｍet×58.15 ）＋0.028 Ｍet＜1.0 のとき、 ＬＯＡＤ* ＝Ｄ1 ＋Ｄ2 ×Ｍet−Ｄ3 ×Ｔsk＋Ｄ4 ×Ｍet×Ｔａ＋（Ｄ5 ＋Ｄ6 ×Ｍet）×ＲＨ×Ｐ（Ｔａ）−Ｆｃｌ×（Ｔsk−Ｔeq* ） Ｍet≧1.0 のとき、 ＬＯＡＤ* ＝Ｄ1 ＋Ｄ2 ×Ｍet−Ｄ3 ×Ｔsk＋Ｄ4 ×Ｍet×Ｔａ＋（Ｄ5 ＋Ｄ6 ×Ｍet）×ＲＨ×Ｐ（Ｔａ）−Ｄ7 ×（Ｍet−１）−Ｆｃｌ×（Ｔsk−Ｔeq* ）
【０００７】
【作用】したがって、本願発明によれば、気温Ｔａと衣
服の熱抵抗Ｉｃｌとに基づき設定温度値θ(th)が定めら
れ、この設定温度値θ(th)にセンサ温度Ｔcrが合致する
ように、加熱手段へ熱エネルギＨθ(th)が与えられる。
すなわち、センサ温度Ｔcrが、一定の値ではなく、状況
に応じて変化する設定温度値θ(th)と等しくなるように
制御される。そして、本願の第１発明，第５発明によれ
ば、センサ温度Ｔcrが設定温度値θ(th)よりも高い場
合、すなわちセンサ温度Ｔcrが設定温度値θ(th)よりも
上昇している加熱状態である場合、気温Ｔａ，設定温度
値θ(th)，ΔＴcr，気流速度Ｖair ，係数ｃ1 〜ｃ6 お
よび係数ｎに基づき、等価温度Ｔeq* が演算される。さ
らに、本願の第２発明，第６発明によれば、第１発明，
第５発明の作用に加えて、センサ温度Ｔcrが設定温度値
θ(th)よりも低い場合、すなわち熱エネルギＨθ(th)を
与える必要がある冷却状態である場合、気温Ｔａ，設定
温度値θ(th)，熱エネルギＨθ(th)および係数ｂ1 〜ｂ
4 に基づき、等価温度Ｔeq* が演算される。また、本願
の第３発明，第７発明によれば、第１発明，第５発明に
より得られた等価温度Ｔeq* を使用し、簡単な演算式を
用いて、予測平均温感ＰＭＶ* が演算される。また、本
願の第４発明，第８発明によれば、第２発明，第６発明
により得られた等価温度Ｔeq* を使用し、簡単な演算式
を用いて、予測平均温感ＰＭＶ* が演算される。
【０００８】
【実施例】以下、本発明を詳細に説明する。
【０００９】図４は本発明に係る温熱感覚演算装置の一
実施例を示すブロック構成図であり、第１，第２，第
５，第６発明に関連する。
【００１０】この温熱感覚演算装置１は、入力部２，演
算部３および表示出力部４を備えている。入力部２に
は、衣服の熱抵抗Ｉｃｌが設定値として与えられている
一方、気流速度Ｖair ，気温Ｔａおよび環境計測部５か
らのセンサ温度Ｔcrが検出値として与えられている。そ
して、これらの設定値および検出値が入力部２より演算
部３へ送られ、演算部３にて、設定温度値θ(th)，熱エ
ネルギとしてのヒータパワーＨθ(th)，等価温度Ｔeq*
が求められる。
【００１１】すなわち、演算部３では、図１および図２
に示すように、その設定温度値演算部３-1にて、気温Ｔ
ａと衣服の熱抵抗Ｉｃｌとに基づき、下記（２）式によ
り設定温度値θ(th)を算出する。なお、下記（２）式に
おいて、ａ1 ，ａ2 については後述する。 θ(th)＝ａ1 ＋ａ2 ×Ｔａ ・・・（２）
【００１２】次に、ヒータ制御演算部３-2にて、設定温
度値θ(th)とセンサ温度Ｔcrとが等しくなるように、例
えばＰＩＤ制御又はファジー制御を行いつつ、この時の
ヒータパワーＨθ(th)を計測する。そして、この求めた
ヒータパワーＨθ(th)（Ｗ）を、入力部２を介して環境
計測部５へ与える。
【００１３】ここで、環境計測部５について、説明を加
える。環境計測部５は、例えば、図５（ａ）に示すよう
に、球形状のモジュール本体５-11 と、このモジュール
本体５-11 の内部に配置された加熱ヒータ５-12 と、モ
ジュール本体５-11 に付設された温度センサ５-13 とか
ら構成されている。モジュール本体５-11 は、銅やアル
ミのように熱伝導性の良い金属により形成されており、
その外表面は放射率が大きくされ、輻射温度Ｔｒに感じ
るものとされている。なお、環境計測部５でのモジュー
ル本体５-11 は、図５（ｂ），（ｃ），（ｄ）に示すよ
うに、半球形状，円筒形状，楕円形状などとしてもよ
い。
【００１４】上述したヒータパワーＨθ(th)は、環境計
測部５における加熱ヒータ５-12 へ与えられる。そし
て、この加熱ヒータ５-12 によってモジュール本体５-1
1 が加熱され、温度センサ５-13 の検出するセンサ温度
Ｔcrが設定温度値θ(th)と等しくなるように制御され
る。
【００１５】ここで、センサ温度Ｔcrが設定温度値θ
(th)よりも低い場合、すなわちヒータパワーＨθ(th)を
与える必要がある冷却状態である場合、図２に示した等
価温度演算部３-3にて、気温Ｔａ，衣服の熱抵抗Ｉｃ
ｌ，設定温度値θ(th)，ヒータパワーＨθ(th)に基づ
き、下記（３）式により、等価温度Ｔeq* を算出する。
なお、下記（３）式において、ｂ1 ，ｂ2 ，ｂ3 ，ｂ4
については後述する。 Ｔeq* ＝ｂ1 ＋ｂ2 ×θ(th)＋ｂ3 ×Ｔａ−ｂ4 ×Ｈθ(th) ・・・（３）
【００１６】これに対して、センサ温度Ｔcrが設定温度
値θ(th)よりも高い場合、すなわちヒータパワーＨθ
(th)が零でセンサ温度Ｔcrが設定温度値θ(th)よりも上
昇している加熱状態である場合、図１に示した等価温度
演算部３-3' にて、気温Ｔａ，設定温度値θ(th)，セン
サ温度Ｔcrと設定温度値θ(th)との差ΔＴcr，気流速度
Ｖair ，衣服の熱抵抗Ｉｃｌを含んで表現される係数ｃ
1 ，係数ｃ2 〜ｃ5 および係数ｎに基づき、下記（４）
式により、等価温度Ｔeq* を算出する。なお、下記
（４）式において、ｃ1 ，ｃ2 ，ｃ3 ，ｃ4 ，ｃ5 につ
いては後述する。 Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ4 ×（１＋ｃ5 ×Ｖairn）×ΔＴ cr ・・・（４）
【００１７】等価温度演算部３-3，３-3' にて得た等価
温度Ｔeq* は表示出力部４へ送り表示する。
【００１８】このように、本実施例による温熱感覚演算
装置１によれば、衣服の熱抵抗Ｉｃｌと気温Ｔａとに基
づき設定温度値θ(th)を求め、センサ温度Ｔcrを設定温
度値θ(th)に合致させるように、加熱ヒータ５-12 への
ヒータパワーＨθ(th)を制御し、センサ温度Ｔcrが設定
温度値θ(th)よりも低い場合には、気温Ｔａ，衣服の熱
抵抗Ｉｃｌ，設定温度値θ(th)，ヒータパワーＨθ(th)
に基づき等価温度Ｔeq* を算出するものとし、センサ温
度Ｔcrが設定温度値θ(th)よりも高い場合には、気温Ｔ
ａ，設定温度値θ(th)，センサ温度Ｔcrと設定温度値θ
(th)との差ΔＴcr，気流速度Ｖair ，係数ｃ1 〜ｃ5 お
よび係数ｎに基づき等価温度Ｔeq* を算出するものとし
ているため、気流速度Ｖair が大きい場合でも、またセ
ンサ温度Ｔcrが設定温度値θ(th)よりも高い場合でも、
等価温度Ｔeq* が人体の感じる等価温度Ｔeqに高精度で
合致するものとなり、正確な等価温度の計測が可能とな
る。
【００１９】次に、等価温度Ｔeq* と等価温度Ｔeqとが
高精度で合致する点について、どのようなメカニズムに
よるものかを補足して説明する。
【００２０】先ず、等価温度演算部３-3より得られる等
価温度Ｔeq* について、この等価温度Ｔeq* と等価温度
Ｔeqとが高精度で合致するメカニズムについて説明す
る。
【００２１】今、環境計測部５へ供給するヒータパワー
をＨとすると、このヒータパワーＨは、周囲環境（気温
Ｔａ，輻射温度Ｔｒ，気流速度Ｖair ）によりそのエネ
ルギが奪われる。このとき、モジュール本体５-11 の表
面の温度がＴ（℃）だとすると、 Ｈ＝Ｈｒ×（Ｔ−Ｔｒ）＋Ａ×Ｖairn×（Ｔ−Ｔａ） ・・・（５） という関係が成り立つ。Ｈｒ：輻射による熱伝達係数×
面積，Ａ×Ｖairn：対流による熱伝達係数×面積。
【００２２】（５）式を変形すると、下記（６）式が得
られる。 0.45×Ｔｒ−0.45×（Ａ／Ｈｒ）×Ｖairn×（Ｔ−Ｔａ）＝0.45×Ｔ−（0.45／ Ｈｒ）×Ｈ ・・・（６）
【００２３】一方、人体が感じる等価温度Ｔeqは下記7)
式により定義される。 Ｔeq＝0.45×Ｔｒ＋0.55×Ｔａ＋〔（0.24−0.75×Ｖair0.5）／（１＋Ｉｃｌ） 〕×（36.5−Ｔａ） ・・・（７）
【００２４】（６）式と（７）式とを比べると、（６）
式でＴ＝36.5とし、かつ、（６）式の両片に0.55×Ｔａ
＋〔0.24／（１＋Ｉｃｌ）〕×（36.5−Ｔａ）を加えれ
ば、等価温度Ｔeqにほゞ等しいものが得られるはずであ
る。すると、（６）式は、下記（８）式となる。 0.45×Ｔｒ＋0.55×Ｔａ＋｛〔0.24−0.45×Ａ×Ｖairn×（１＋Ｉｃｌ）／Ｈｒ 〕／（１＋Ｉｃｌ）｝×（36.5−Ｔａ）＝0.45×36.5＋0.55×Ｔａ＋〔0.24／（ １＋Ｉｃｌ）〕×（36.5−Ｔａ）−（0.45／Ｈｒ）×Ｈ ・・・（８）
【００２５】さて、（８）式と（７）式とを比べると、
その違いは、0.45×Ａ×〔（１＋Ｉｃｌ）／Ｈｒ〕×Ｖ
airnと0.75×Ｖair0.5の部分であることが分かる。
【００２６】さて、環境計測部５は小形に作っているの
で、その気流感度がどうしても人体の気流感度より大き
くなってしまう。すなわち、その表面効果により気流感
度部分が人体のそれよりも大きく、0.45×Ａ×〔（１＋
Ｉｃｌ）／Ｈｒ〕＞0.75となり、どうしてもこの不等号
が最後までつきまとい、気流速度Ｖair が大きいところ
では人体の感じる等価温度Ｔeqよりも大幅に低い値を示
し、正確な等価温度を求めることができなくなる。
【００２７】そこで、本実施例においては、気流感度を
一致させるために、気流感度の係数のみ注目するのでは
なく、上記（８）式および（７）式において、（36.5−
Ｔａ）の温度差まで含めて考え直した。これにより、0.
45×Ａ×〔（１＋Ｉｃｌ）／Ｈｒ〕×Ｖairn×（θ(th)
−Ｔａ）＝0.75×Ｖair0.5×（36.5−Ｔａ）となるよう
に、36.5℃ではなく、θ(th)℃という温度で環境計測部
５へのヒータパワーＨθ(th)を制御すればよいことに気
付いた。
【００２８】 すなわち、θ(th)＝｛0.75／〔0.45×Ａ×（１＋Ｉｃｌ）／Ｈｒ〕｝×（Ｖai r0.5／Ｖairn）×（36.5−Ｔａ）＋Ｔａ ・・・（９） と表され、0.1 〜1.0 ｍ／ｓの範囲では、ｎは0.35〜0.
6 の範囲の値をとるから、Ｖair0.5 ≒Ｖairnとみな
せ、 θ(th)＝1.667 ×｛Ｈｒ／〔Ａ×（１＋Ｉｃｌ）〕｝×（36.5−Ｔａ）＋Ｔａ ・・・（１０） と表される。
【００２９】したがって、衣服の熱抵抗Ｉｃｌと気温Ｔ
ａとに基づき設定温度値θ(th)を求め、センサ温度Ｔcr
を設定温度値θ(th)に合致させるように、加熱ヒータ５
-12へのヒータパワーＨθ(th)を制御し、気温Ｔａ，衣
服の熱抵抗Ｉｃｌ，設定温度値θ(th)，ヒータパワーＨ
θ(th)に基づき、 0.45×θ(th)＋0.55×Ｔａ＋〔0.24／（１＋Ｉｃｌ）〕×（36.5−Ｔａ）−（0. 45／Ｈｒ）×Ｈθ(th) ・・・（１１） を計算して等価温度Ｔeq* を求めるものとすれば、通常
のビルオフィスでは、気流速度Ｖair が大きい場合で
も、正しく等価温度を計測することができるものとな
る。
【００３０】なお、（１０）式において、ａ1 ＝1.667
×｛Ｈｒ／〔Ａ×（１＋Ｉｃｌ）〕｝×36.5、ａ2 ＝１
−1.667 ×｛Ｈｒ／〔Ａ×（１＋Ｉｃｌ）〕｝とする
と、θ(t h)＝ａ1 ＋ａ2 ×Ｔａとして表現される。
【００３１】また、（１１）式において、ｂ1 ＝〔0.24
／（１＋Ｉｃｌ）〕×36.5、ｂ2 ＝0.45、ｂ3 ＝0.55−
0.24／（１＋Ｉｃｌ）、ｂ4 ＝＋0.45／Ｈｒとすると、 Ｔeq* ＝ｂ1 ＋ｂ2 ×θ(th)＋ｂ3 ×Ｔａ−ｂ4 ×Ｈθ(th) と表現される。
【００３２】次に、等価温度演算部３-3' より得られる
等価温度Ｔeq* について、この等価温度Ｔeq* と等価温
度Ｔeqとが高精度で合致するメカニズムについて説明す
る。
【００３３】ヒータパワーＨθ(th)が零でセンサ温度Ｔ
crが設定温度値θ(th)よりも上昇している加熱状態と
は、何らかの形でモジュール本体５-11 が加熱されてい
ることを示している。ここで、再度、環境計測部５にお
ける関係式を調べてみる。Ｈθ(th)＝０なので、
【００３４】 ０＝Ｈｒ×（Ｔcr−Ｔｒ）＋Ａ×Ｖairn×（Ｔcr−Ｔａ） ・・・（１２） となっている。ところで、Ｔcr＞θ(th)なので、Ｔcr＝
θ(th)＋ΔＴcr（ΔＴcr≧０）と置くと、 ０＝Ｈｒ×〔（θ(th)＋ΔＴcr）−Ｔｒ）〕＋Ａ×Ｖairn×〔（θ(th)＋ΔＴcr ）−Ｔａ）〕 ・・・（１３） となり、 ０＝Ｈｒ×（θ(th)−Ｔｒ）＋Ａ×Ｖairn×（θ(th)−Ｔａ）＋（Ｈｒ＋Ａ×Ｖ airn）×ΔＴcr ・・・（１４） となる。
【００３５】式（１４）を再び変形すると、 0.45×Ｔｒ−0.45×（Ａ／Ｈｒ）×Ｖairn×（θ(th)−Ｔａ）＝0.45×θ(th)＋ 0.45×〔１＋（Ａ／Ｈｒ）×Ｖairn〕×ΔＴcr ・・・（１５） となる。
【００３６】前述した式（６）から式（８）への移行の
際と同様にして、式（１５）の両片に、0.55×Ｔａ＋
〔0.24／（１＋Ｉｃｌ）〕×（36.5−Ｔａ）を加える
と、 0.45×Ｔｒ＋0.55×Ｔａ＋〔0.24／（１＋Ｉｃｌ）〕×（36.5−Ｔａ）−0.45× （Ａ／Ｈｒ）×Ｖairn×（θ(th)−Ｔａ）＝0.45×θ(th)＋0.55×Ｔａ＋〔0.24 ／（１＋Ｉｃｌ）〕×（36.5−Ｔａ）＋0.45×（１＋（Ａ／Ｈｒ）×Ｖairn）× ΔＴcr ・・・（１６） となる。上式の左辺は、0.45×（Ａ／Ｈｒ）×Ｖairn×
（θ(th)−Ｔａ）＝〔0.75／（１＋Ｉｃｌ）〕×Ｖair
0.5×（36.5−Ｔａ）であるので、まさに等価温度Ｔeq
を表している。
【００３７】したがって、衣服の熱抵抗Ｉｃｌと気温Ｔ
ａとに基づき設定温度値θ(th)を求め、気温Ｔａ，設定
温度値θ(th)，センサ温度Ｔcrと設定温度値θ(th)との
差ΔＴcr，気流速度Ｖair ，衣服の熱抵抗Ｉｃｌに基づ
き、 0.45×θ(th)＋0.55×Ｔａ＋〔0.24／（１＋Ｉｃｌ）〕×（36.5−Ｔａ）＋0.45 ×（１＋（Ａ／Ｈｒ）×Ｖairn）×ΔＴcr ・・・（１７） を計算して等価温度Ｔeq* を求めるものとすれば、セン
サ温度Ｔcrが設定温度値θ(th)より高い場合であって
も、正しく等価温度Ｔeqを計測することができるものと
なる。
【００３８】なお、（１７）式において、ｃ1 ＝〔0.24
／（１＋Ｉｃｌ）〕×36.5、ｃ2 ＝0.45、ｃ3 ＝0.55−
0.24／（１＋Ｉｃｌ）、ｃ4 ＝0.45、ｃ5 ＝Ａ／Ｈｒと
すると、Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ
4 ×（１＋ｃ5 ×Ｖairn）×ΔＴcrと表現される。
【００３９】また、上述では、ヒータパワーＨθ(th)が
零でセンサ温度Ｔcrが設定温度値θ(th)よりも上昇して
いる加熱状態について述べたが、前記（４）式は下記
（４’）式においてヒータパワーＨθ(th)が零である場
合の一状態式に他ならない。 Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ4 ×（１＋ｃ5 ×Ｖairn）×ΔＴ cr −ｃ6 ×Ｈθ(th) ・・・（４’） なお、（４’）式において、係数ｃ6 はｃ6 ＝＋0.45／
Ｈｒである。
【００４０】また、上述の実施例においては、気流速度
Ｖair を求めるにあたり、気流速センサを使用すること
になるが、環境計測部５を室内に置いていることを予想
すると、また室内の気流速度は0.1 〜0.3 ｍ／ｓが多い
点を考え合わせると、気流速度Ｖair をある任意の値、
例えば0.2 ｍ／ｓの数値として定めたとしても、それ程
等価温度Ｔeqに対して誤差は生じない。すなわち、気流
速センサは必ずしも使用しなくてよく、これによって簡
素化を維持することができるものとなり、コストアップ
を来さないものともなる。もっとも、気流速度Ｖair が
強くなると、モジュール本体５-11 から熱エネルギが奪
われ、センサ温度Ｔcrもいずれ設定温度値θ(th)より低
下し、通常のヒータ加熱が再現するので、等価温度Ｔeq
は冷却状態に際しての方式でいずれ求まることになる。
【００４１】図６は本発明に係る予測平均温感演算装置
の一実施例を示すブロック構成図であり、第３，第４，
第７，第８発明に関連する。
【００４２】この予測平均温感演算装置６は、入力部
７，演算部８および表示出力部９を備えている。入力部
７には、人体の活動量Ｍet，衣服の熱抵抗Ｉｃｌが設定
値として与えられている一方、気温Ｔａ，湿度ＲＨ，等
価温度Ｔeq* が検出値として与えられている。等価温度
Ｔeq*は、図２に示した温熱感覚演算装置１を使用し、
その演算部３からの出力として与えられるものである。
そして、これらの設定値および検出値が入力部７より演
算部８へ送られ、演算部８にて、下記（１８）式により
予測平均温感ＰＭＶ* を算出する。
【００４３】 ＰＭＶ* ＝α×ＬＯＡＤ* ・・・（１８） 但し、（１８）式において、 α＝0.303 ×ＥＸＰ（−0.036 ×Ｍet×58.15 ）＋0.028 Ｍet＜1.0 のとき、 ＬＯＡＤ* ＝Ｄ1 ＋Ｄ2 ×Ｍet−Ｄ3 ×Ｔsk＋Ｄ4 ×Ｍet×Ｔａ＋（Ｄ5 ＋Ｄ6 ×Ｍet）×ＲＨ×Ｐ（Ｔａ）−Ｆｃｌ×（Ｔsk−Ｔeq* ） Ｍet≧1.0 のとき、 ＬＯＡＤ* ＝Ｄ1 ＋Ｄ2 ×Ｍet−Ｄ3 ×Ｔsk＋Ｄ4 ×Ｍet×Ｔａ＋（Ｄ5 ＋Ｄ6 ×Ｍet）×ＲＨ×Ｐ（Ｔａ）−Ｄ7 ×（Ｍet−１）−Ｆｃｌ×（Ｔsk−Ｔeq* ） Ｄ1 〜Ｄ7 の数値例として、Ｄ1 ＝10.2895 ，Ｄ2 ＝4
9.5152 ，Ｄ3 ＝0.7808，Ｄ4 ＝0.08142 ，Ｄ5 ＝0.004
067，Ｄ6 ＝0.00133338，Ｄ7 ＝24.423，Ｔsk：皮膚温
（Ｔsk＝35.7−0.0275×Ｍet×58.15 ），Ｐ（Ｔａ）：
気温Ｔａにおける絶対水蒸気圧〔Ｐ（Ｔａ）＝ＥＸＰ
（18.6686 −4030.183／（Ｔａ＋235 ）〕，Ｆｃｌ：熱
抵抗の効果（Ｆｃｌ＝Ｅ１／（Ｅ２＋Ｅ３×Ｉｃ
ｌ））。Ｅ１〜Ｅ３の数値例として、Ｅ１＝8.976，Ｅ
２＝１，Ｅ３＝1.466。
【００４４】なお、Ｐ（Ｔａ）について本実施例ではAn
toinneの関係式（飽和水蒸気圧）を使用しているが、別
の関係式でもかまわない。
【００４５】図３は演算部８での演算状況を示すブロッ
ク構成図である。すなわち、演算部８では、そのＦｃｌ
演算部８-1にて、衣服の熱抵抗Ｉｃｌを用い、熱抵抗の
効果Ｆｃｌを算出する。また、α・Ｔsk演算部８-2に
て、活動量Ｍetを用い、αおよび皮膚温Ｔskを算出す
る。また、ＰＭＶＬＯＡＤ演算部８-3にて、皮膚温Ｔs
k，活動量Ｍet，気温Ｔａ，相対湿度ＲＨ，絶対水蒸気
圧Ｐ（Ｔａ）を用い、ＰＭＶＬＯＡＤを算出する。ま
た、演算部８-4にて、熱抵抗の効果Ｆｃｌ，皮膚温Ｔs
k，等価温度Ｔeq* を用い、Ｆｃｌ×（Ｔsk−Ｔeq* ）
を算出する。そして、ＰＭＶ* 演算部８-5にて、演算部
８-3で得たＰＭＶＬＯＡＤより演算部８-4で得たＦｃｌ
×（Ｔsk−Ｔeq* ）を差し引いてＬＯＡＤ* を求め、こ
のＬＯＡＤ* に演算部８-2で得たαを乗じて、予測平均
温感ＰＭＶ* を算出する。
【００４６】ここで、上記（１８）式は、本出願人が従
来のＰＭＶ方程式（（１）式）を簡素化すべく、試行錯
誤して得た式であり、この式により得られるＰＭＶ*
は、等価温度Ｔeq* が精度良く得られることとの相乗効
果により、従来のＰＭＶ方程式で得られるＰＭＶに対し
てかなりの精度で合致する。表１に、有風，夏期，Ｍｅ
ｔ＝1.2 ，Ｉｃｌ＝0.5 ，ＲＨ＝50％，Ｔａ＝Ｔｒ＝30
℃，Ｖair ＝0.4 ｍ／ｓでの、ＰＭＶとＰＭＶ* との比
較を示す。また、表２に、有風，冬期，Ｍｅｔ＝1.2 ，
Ｉｃｌ＝1.0 ，ＲＨ＝50％，Ｔａ＝Ｔｒ＝18℃，Ｖair
＝0.4m/sでの、ＰＭＶとＰＭＶ* との比較を示す。これ
らの表からも、ＰＭＶとＰＭＶ* とは、微小誤差（0.1
前後）をもって合致することが分かる。
【００４７】
【表１】
【００４８】
【表２】
【００４９】なお、上表において、ＬＯＡＤとは、式
（１）をＰＭＶ＝（0.303 ｅ-0.036M＋0.028 ）×ＬＯ
ＡＤとした時、一つにまとめられる部分をさす。また、
Ｒ＋Ｃ＝3.96・10-8・ｆcl・｛（Ｔcl＋274 ）4 −（Ｔ
ｒ＋273 ）4 ｝＋ｆcl・ｈc ・（Ｔcl−Ｔa ）で、Ｒ*
＋Ｃ*はＦｃｌ×（Ｔsk−Ｔeq* ）を意味する。
【００５０】
【発明の効果】以上説明したことから明らかなように本
発明によると、センサ温度Ｔcrが、一定の値ではなく、
状況に応じて変化する設定温度値θ(th)と等しくなるよ
うに制御され、本願の第１発明，第５発明によれば、セ
ンサ温度Ｔcrが設定温度値θ(t h)よりも高い場合、気温
Ｔａ，設定温度値θ(th)，ΔＴcr，気流速度Ｖair ，係
数ｃ1 〜ｃ6 および係数ｎに基づき等価温度Ｔeq* が演
算され、本願の第２発明，第６発明によれば、さらに、
センサ温度Ｔcrが設定温度値θ(th)よりも低い場合、気
温Ｔａ，設定温度値θ(th)，熱エネルギＨθ(th)および
係数ｂ1 〜ｂ4 に基づき、等価温度Ｔeq* が演算される
ので、気流速度Ｖair が大きい場合でも、またセンサ温
度Ｔcrが設定温度値θ(th)よりも高い場合でも、その演
算して得た等価温度Ｔeq* が人体が感じる等価温度Ｔeq
に高精度に合致するものとなり、正確な等価温度Ｔeqの
計測が可能となる。また、本願の第３，第４，第７，第
８発明によれば、高精度で等価温度Ｔeqに合致する等価
温度Ｔeq*を使用するので、正確な予測平均温感ＰＭＶ*
を求めることができるようになる。
【図面の簡単な説明】
【図１】図４に示した温熱感覚演算装置においてその演
算部での演算状況を示すブロック構成図。
【図２】図４に示した温熱感覚演算装置においてその演
算部での演算状況を示すブロック構成図。
【図３】図６に示した予測平均温感演算装置においてそ
の演算部での演算状況を示すブロック構成図。
【図４】本発明に係る温熱感覚演算装置の一実施例を示
すブロック構成図。
【図５】環境計測部の一例を示す概略的な縦断面図。
【図６】本発明に係る予測平均温感演算装置の一実施例
を示すブロック構成図。
【符号の説明】
１ 温熱感覚演算装置 ３ 演算部 ３−１ 設定温度値演算部 ３−２ ヒータ制御演算部 ５ 環境計測部 ６ 予測平均温感演算装置 ８ 演算部 ８−５ ＰＭＶ*演算部 ３−３，３−３’ 等価温度演算部

Claims (8)

    【特許請求の範囲】
  1. 【請求項１】 気温Ｔａと衣服の熱抵抗Ｉｃｌとに基づ
    き設定温度値θ(th)を演算し、 この設定温度値θ(th)とセンサ温度Ｔcrとが等しくなる
    ように、前記センサ温度Ｔcrを調整し得る加熱手段へ熱
    エネルギＨθ(th)を供与し、 前記センサ温度Ｔcrが前記設定温度値θ(th)よりも高い
    場合、前記気温Ｔａ，前記設定温度値θ(th)，前記セン
    サ温度Ｔcrと前記設定温度値θ(th)との差ΔＴcr，気流
    速度Ｖair ，衣服の熱抵抗Ｉｃｌを含んで表現される係
    数ｃ1 ，係数ｃ2 〜ｃ6 および係数ｎに基づき、下記
    （ｂ）式により等価温度Ｔeq* を演算するようにしたこ
    とを特徴とする温熱感覚演算方法。 Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ4 ×（１＋ｃ5 ×Ｖairn）×ΔＴ cr −ｃ6 ×Ｈθ(th) ・・・（ｂ）
  2. 【請求項２】 気温Ｔａと衣服の熱抵抗Ｉｃｌとに基づ
    き設定温度値θ(th)を演算し、 この設定温度値θ(th)とセンサ温度Ｔcrとが等しくなる
    ように、前記センサ温度Ｔcrを調整し得る加熱手段へ熱
    エネルギＨθ(th)を供与し、 前記センサ温度Ｔcrが前記設定温度値θ(th)よりも低い
    場合、前記気温Ｔａ，前記設定温度値θ(th)，前記熱エ
    ネルギＨθ(th)および係数ｂ1 〜ｂ4 に基づき、下記
    （ａ）式により等価温度Ｔeq* を演算する一方、 前記センサ温度Ｔcrが前記設定温度値θ(th)よりも高い
    場合、前記気温Ｔａ，前記設定温度値θ(th)，前記セン
    サ温度Ｔcrと前記設定温度値θ(th)との差ΔＴcr，気流
    速度Ｖair ，衣服の熱抵抗Ｉｃｌを含んで表現される係
    数ｃ1 ，係数ｃ2 〜ｃ6 および係数ｎに基づき、下記
    （ｂ）式により等価温度Ｔeq* を演算するようにしたこ
    とを特徴とする温熱感覚演算方法。 Ｔeq* ＝ｂ1 ＋ｂ2 ×θ(th)＋ｂ3 ×Ｔａ−ｂ4 ×Ｈθ(th) ・・・（ａ） Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ4 ×（１＋ｃ5 ×Ｖairn）×ΔＴcr −ｃ6 ×Ｈθ(th) ・・・（ｂ）
  3. 【請求項３】 気温Ｔａと衣服の熱抵抗Ｉｃｌとに基づ
    き設定温度値θ(th)を演算し、 この設定温度値θ(th)とセンサ温度Ｔcrとが等しくなる
    ように、前記センサ温度Ｔcrを調整し得る加熱手段へ熱
    エネルギＨθ(th)を供与し、 前記センサ温度Ｔcrが前記設定温度値θ(th)よりも高い
    場合、前記気温Ｔａ，前記設定温度値θ(th)，前記セン
    サ温度Ｔcrと前記設定温度値θ(th)との差ΔＴcr，気流
    速度Ｖair ，衣服の熱抵抗Ｉｃｌを含んで表現される係
    数ｃ1 ，係数ｃ2 〜ｃ6 および係数ｎに基づき、下記
    （ｂ）式により等価温度Ｔeq* を演算するものとし、 この温熱感覚演算方法を適用して得た等価温度Ｔeq* ，
    湿度ＲＨ，活動量Ｍet，皮膚温Ｔsk（Ｔsk＝35.7−0.02
    75×Ｍet×58.15 ），気温Ｔａ，気温Ｔａにおける絶対
    水蒸気圧Ｐ（Ｔａ），衣服の熱抵抗Ｉｃｌ，係数Ｅ１〜
    Ｅ３に基づく熱抵抗の効果Ｆｃｌ（Ｆｃｌ＝Ｅ１／（Ｅ
    ２＋Ｅ３×Ｉｃｌ））および係数Ｄ1 〜Ｄ7 に基づき、
    下記（ｃ）式により予測平均温感ＰＭＶ* を演算するよ
    うにしたことを特徴とする予測平均温感演算方法。 Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ4 ×（１＋ｃ5 ×Ｖairn）×ΔＴ cr −ｃ6 ×Ｈθ(th) ・・・（ｂ） ＰＭＶ* ＝α×ＬＯＡＤ* ・・・（ｃ） 但し、（ｃ）式において、 α＝0.303 ×ＥＸＰ（−0.036 ×Ｍet×58.15 ）＋0.028 Ｍet＜1.0 のとき、 ＬＯＡＤ* ＝Ｄ1 ＋Ｄ2 ×Ｍet−Ｄ3 ×Ｔsk＋Ｄ4 ×Ｍet×Ｔａ＋（Ｄ5 ＋Ｄ6 ×Ｍet）×ＲＨ×Ｐ（Ｔａ）−Ｆｃｌ×（Ｔsk−Ｔeq* ） Ｍet≧1.0 のとき、 ＬＯＡＤ* ＝Ｄ1 ＋Ｄ2 ×Ｍet−Ｄ3 ×Ｔsk＋Ｄ4 ×Ｍet×Ｔａ＋（Ｄ5 ＋Ｄ6 ×Ｍet）×ＲＨ×Ｐ（Ｔａ）−Ｄ7 ×（Ｍet−１）−Ｆｃｌ×（Ｔsk−Ｔeq* ）
  4. 【請求項４】 気温Ｔａと衣服の熱抵抗Ｉｃｌとに基づ
    き設定温度値θ(th)を演算し、この設定温度値θ(th)と
    センサ温度Ｔcrとが等しくなるように、前記センサ温度
    Ｔcrを調整し得る加熱手段へ熱エネルギＨθ(th)を供与
    し、 前記センサ温度Ｔcrが前記設定温度値θ(th)よりも低い
    場合、前記気温Ｔａ，前記設定温度値θ(th)，前記熱エ
    ネルギＨθ(th)および係数ｂ1 〜ｂ4 に基づき、下記
    （ａ）式により等価温度Ｔeq* を演算する一方、 前記センサ温度Ｔcrが前記設定温度値θ(th)よりも高い
    場合、前記気温Ｔａ，前記設定温度値θ(th)，前記セン
    サ温度Ｔcrと前記設定温度値θ(th)との差ΔＴcr，気流
    速度Ｖair ，衣服の熱抵抗Ｉｃｌを含んで表現される係
    数ｃ1 ，係数ｃ2 〜ｃ6 および係数ｎに基づき、下記
    （ｂ）式により等価温度Ｔeq* を演算するものとし、 この温熱感覚演算方法を適用して得た等価温度Ｔeq* ，
    湿度ＲＨ，活動量Ｍet，皮膚温Ｔsk（Ｔsk＝35.7−0.02
    75×Ｍet×58.15 ），気温Ｔａ，気温Ｔａにおける絶対
    水蒸気圧Ｐ（Ｔａ），衣服の熱抵抗Ｉｃｌ，係数Ｅ１〜
    Ｅ３に基づく熱抵抗の効果Ｆｃｌ（Ｆｃｌ＝Ｅ１／（Ｅ
    ２＋Ｅ３×Ｉｃｌ））および係数Ｄ1 〜Ｄ7 に基づき、
    下記（ｃ）式により予測平均温感ＰＭＶ* を演算するよ
    うにしたことを特徴とする予測平均温感演算方法。 Ｔeq* ＝ｂ1 ＋ｂ2 ×θ(th)＋ｂ3 ×Ｔａ−ｂ4 ×Ｈθ(th) ・・・（ａ） Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ4 ×（１＋ｃ5 ×Ｖairn）×ΔＴ cr −ｃ6 ×Ｈθ(th) ・・・（ｂ） ＰＭＶ* ＝α×ＬＯＡＤ* ・・・（ｃ） 但し、（ｃ）式において、 α＝0.303 ×ＥＸＰ（−0.036 ×Ｍet×58.15 ）＋0.028 Ｍet＜1.0 のとき、 ＬＯＡＤ* ＝Ｄ1 ＋Ｄ2 ×Ｍet−Ｄ3 ×Ｔsk＋Ｄ4 ×Ｍet×Ｔａ＋（Ｄ5 ＋Ｄ6 ×Ｍet）×ＲＨ×Ｐ（Ｔａ）−Ｆｃｌ×（Ｔsk−Ｔeq* ） Ｍet≧1.0 のとき、 ＬＯＡＤ* ＝Ｄ1 ＋Ｄ2 ×Ｍet−Ｄ3 ×Ｔsk＋Ｄ4 ×Ｍet×Ｔａ＋（Ｄ5 ＋Ｄ6 ×Ｍet）×ＲＨ×Ｐ（Ｔａ）−Ｄ7 ×（Ｍet−１）−Ｆｃｌ×（Ｔsk−Ｔeq* ）
  5. 【請求項５】 センサ温度Ｔcrを調整し得る加熱手段を
    有する環境計測部に対して設けられる温熱感覚演算装置
    であって、 気温Ｔａと衣服の熱抵抗Ｉｃｌとに基づき設定温度値θ
    (th)を演算する設定温度値演算手段と、 前記設定温度値θ(th)と前記センサ温度Ｔcrとが等しく
    なるように前記加熱手段へ熱エネルギＨθ(th)を供与す
    る熱エネルギ供与手段と、 前記センサ温度Ｔcrが前記設定温度値θ(th)よりも高い
    場合、前記気温Ｔａ，前記設定温度値θ(th)，前記セン
    サ温度Ｔcrと前記設定温度値θ(th)との差ΔＴcr，気流
    速度Ｖair ，衣服の熱抵抗Ｉｃｌを含んで表現される係
    数ｃ1 ，係数ｃ2 〜ｃ6 および係数ｎに基づき、下記
    （ｂ）式により等価温度Ｔeq* を演算する温熱感覚演算
    手段とを備えたことを特徴とする温熱感覚演算装置。 Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ4×（１＋ｃ5 ×Ｖairn ）×ΔＴ cr −ｃ6 ×Ｈθ(th) ・・・（ｂ）
  6. 【請求項６】 センサ温度Ｔcrを調整し得る加熱手段を
    有する環境計測部に対して設けられる温熱感覚演算装置
    であって、 気温Ｔａと衣服の熱抵抗Ｉｃｌとに基づき設定温度値θ
    (th)を演算する設定温度値演算手段と、 前記設定温度値θ(th)と前記センサ温度Ｔcrとが等しく
    なるように前記加熱手段へ熱エネルギＨθ(th)を供与す
    る熱エネルギ供与手段と、 前記センサ温度Ｔcrが前記設定温度値θ(th)よりも低い
    場合、前記気温Ｔａ，前記設定温度値θ(th)，前記熱エ
    ネルギＨθ(th)および係数ｂ1 〜ｂ4 に基づき、下記
    （ａ）式により等価温度Ｔeq* を演算する一方、 前記センサ温度Ｔcrが前記設定温度値θ(th)よりも高い
    場合、前記気温Ｔａ，前記設定温度値θ(th)，前記セン
    サ温度Ｔcrと前記設定温度値θ(th)との差ΔＴcr，気流
    速度Ｖair ，衣服の熱抵抗Ｉｃｌを含んで表現される係
    数ｃ1 ，係数ｃ2 〜ｃ6 および係数ｎに基づき、下記
    （ｂ）式により等価温度Ｔeq* を演算する温熱感覚演算
    手段とを備えたことを特徴とする温熱感覚演算装置。 Ｔeq* ＝ｂ1 ＋ｂ2 ×θ(th)＋ｂ3 ×Ｔａ−ｂ4 ×Ｈθ(th) ・・・（ａ） Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ4 ×（１＋ｃ5 ×Ｖairn）×ΔＴ cr −ｃ6 ×Ｈθ(th) ・・・（ｂ）
  7. 【請求項７】 センサ温度Ｔcrを調整し得る加熱手段を
    有する環境計測部に対して設けられる温熱感覚演算装置
    であって、 気温Ｔａと衣服の熱抵抗Ｉｃｌとに基づき設定温度値θ
    (th)を演算する設定温度値演算手段と、 前記設定温度値θ(th)と前記センサ温度Ｔcrとが等しく
    なるように前記加熱手段へ熱エネルギＨθ(th)を供与す
    る熱エネルギ供与手段と、 前記センサ温度Ｔcrが前記設定温度値θ(th)よりも高い
    場合、前記気温Ｔａ，前記設定温度値θ(th)，前記セン
    サ温度Ｔcrと前記設定温度値θ(th)との差ΔＴcr，気流
    速度Ｖair ，衣服の熱抵抗Ｉｃｌを含んで表現される係
    数ｃ1 ，係数ｃ2 〜ｃ6 および係数ｎに基づき、下記
    （ｂ）式により等価温度Ｔeq* を演算する温熱感覚演算
    手段とを備えた温熱感覚演算装置と、 この温熱感覚演算装置を使用して求めた等価温度Ｔeq
    * ，湿度ＲＨ，活動量Ｍet，皮膚温Ｔsk（Ｔsk＝35.7−
    0.0275×Ｍet×58.15 ），気温Ｔａ，気温Ｔａにおける
    絶対水蒸気圧Ｐ（Ｔａ），衣服の熱抵抗Ｉｃｌ，係数Ｅ
    １〜Ｅ３に基づく熱抵抗の効果Ｆｃｌ（Ｆｃｌ＝Ｅ１／
    （Ｅ２＋Ｅ３×Ｉｃｌ））および係数Ｄ1〜Ｄ7 に基づ
    き、下記（ｃ）式により予測平均温感ＰＭＶ* を演算す
    る予測平均温感演算手段とを備えた予測平均温感演算装
    置。 Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ4 ×（１＋ｃ5 ×Ｖairn）×ΔＴ cr −ｃ6 ×Ｈθ(th) ・・・（ｂ） ＰＭＶ* ＝α×ＬＯＡＤ* ・・・（ｃ） 但し、（ｃ）式において、 α＝0.303 ×ＥＸＰ（−0.036 ×Ｍet×58.15 ）＋0.028 Ｍet＜1.0 のとき、 ＬＯＡＤ* ＝Ｄ1 ＋Ｄ2 ×Ｍet−Ｄ3 ×Ｔsk＋Ｄ4 ×Ｍet×Ｔａ＋（Ｄ5 ＋Ｄ6 ×Ｍet）×ＲＨ×Ｐ（Ｔａ）−Ｆｃｌ×（Ｔsk−Ｔeq* ） Ｍet≧1.0 のとき、 ＬＯＡＤ* ＝Ｄ1 ＋Ｄ2 ×Ｍet−Ｄ3 ×Ｔsk＋Ｄ4 ×Ｍet×Ｔａ＋（Ｄ5 ＋Ｄ6 ×Ｍet）×ＲＨ×Ｐ（Ｔａ）−Ｄ7 ×（Ｍet−１）−Ｆｃｌ×（Ｔsk−Ｔeq* ）
  8. 【請求項８】 センサ温度Ｔcrを調整し得る加熱手段を
    有する環境計測部に対して設けられる温熱感覚演算装置
    であって、 気温Ｔａと衣服の熱抵抗Ｉｃｌとに基づき設定温度値θ
    (th)を演算する設定温度値演算手段と、 前記設定温度値θ(th)と前記センサ温度Ｔcrとが等しく
    なるように前記加熱手段へ熱エネルギＨθ(th)を供与す
    る熱エネルギ供与手段と、 前記センサ温度Ｔcrが前記設定温度値θ(th)よりも低い
    場合、前記気温Ｔａ，前記設定温度値θ(th)，前記熱エ
    ネルギＨθ(th)および係数ｂ1 〜ｂ4 に基づき、下記
    （ａ）式により等価温度Ｔeq* を演算する一方、 前記センサ温度Ｔcrが前記設定温度値θ(th)よりも高い
    場合、前記気温Ｔａ，前記設定温度値θ(th)，前記セン
    サ温度Ｔcrと前記設定温度値θ(th)との差ΔＴcr，気流
    速度Ｖair ，衣服の熱抵抗Ｉｃｌを含んで表現される係
    数ｃ1 ，係数ｃ2 〜ｃ6 および係数ｎに基づき、下記
    （ｂ）式により等価温度Ｔeq* を演算する温熱感覚演算
    手段とを備えた温熱感覚演算装置と、 この温熱感覚演算装置を使用して求めた等価温度Ｔeq
    * ，湿度ＲＨ，活動量Ｍet，皮膚温Ｔsk（Ｔsk＝35.7−
    0.0275×Ｍet×58.15 ），気温Ｔａ，気温Ｔａにおける
    絶対水蒸気圧Ｐ（Ｔａ），衣服の熱抵抗Ｉｃｌ，係数Ｅ
    １〜Ｅ３に基づく熱抵抗の効果Ｆｃｌ（Ｆｃｌ＝Ｅ１／
    （Ｅ２＋Ｅ３×Ｉｃｌ））および係数Ｄ1〜Ｄ7 に基づ
    き、下記（ｃ）式により予測平均温感ＰＭＶ* を演算す
    る予測平均温感演算手段とを備えた予測平均温感演算装
    置。 Ｔeq* ＝ｂ1 ＋ｂ2 ×θ(th)＋ｂ3 ×Ｔａ−ｂ4 ×Ｈθ(th) ・・・（ａ） Ｔeq* ＝ｃ1 ＋ｃ2 ×θ(th)＋ｃ3 ×Ｔａ＋ｃ4 ×（１＋ｃ5 ×Ｖairn）×ΔＴ cr −ｃ6 ×Ｈθ(th) ・・・（ｂ） ＰＭＶ* ＝α×ＬＯＡＤ* ・・・（ｃ） 但し、（ｃ）式において、 α＝0.303 ×ＥＸＰ（−0.036 ×Ｍet×58.15 ）＋0.028 Ｍet＜1.0 のとき、 ＬＯＡＤ* ＝Ｄ1 ＋Ｄ2 ×Ｍet−Ｄ3 ×Ｔsk＋Ｄ4 ×Ｍet×Ｔａ＋（Ｄ5 ＋Ｄ6 ×Ｍet）×ＲＨ×Ｐ（Ｔａ）−Ｆｃｌ×（Ｔsk−Ｔeq* ） Ｍet≧1.0 のとき、 ＬＯＡＤ* ＝Ｄ1 ＋Ｄ2 ×Ｍet−Ｄ3 ×Ｔsk＋Ｄ4 ×Ｍet×Ｔａ＋（Ｄ5 ＋Ｄ6 ×Ｍet）×ＲＨ×Ｐ（Ｔａ）−Ｄ7 ×（Ｍet−１）−Ｆｃｌ×（Ｔsk−Ｔeq* ）
JP3254305A 1990-09-25 1991-09-06 温熱感覚演算方法および装置、予測平均温感演算方法および装置 Expired - Lifetime JP2581625B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2-256723 1990-09-25
JP25672390 1990-09-25

Publications (2)

Publication Number Publication Date
JPH0518592A true JPH0518592A (ja) 1993-01-26
JP2581625B2 JP2581625B2 (ja) 1997-02-12

Family

ID=17296556

Family Applications (1)

Application Number Title Priority Date Filing Date
JP3254305A Expired - Lifetime JP2581625B2 (ja) 1990-09-25 1991-09-06 温熱感覚演算方法および装置、予測平均温感演算方法および装置

Country Status (5)

Country Link
US (1) US5333953A (ja)
EP (1) EP0477796B1 (ja)
JP (1) JP2581625B2 (ja)
DE (1) DE69112162T2 (ja)
DK (1) DK0477796T3 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0618098A1 (en) * 1993-04-01 1994-10-05 FIAT AUTO S.p.A. A process and system for controlling environmental conditions
JPH08100941A (ja) * 1994-03-09 1996-04-16 Lg Electronics Inc 空調システムのｐｍｖ算出方法
JP2008089343A (ja) * 2006-09-29 2008-04-17 Toyota Motor Corp 接触温熱感計測用センサ装置および接触温熱感計測装置
CN109974244A (zh) * 2014-02-17 2019-07-05 松下电器产业株式会社 空调机的控制方法
CN113088359A (zh) * 2021-03-30 2021-07-09 重庆大学 一种工艺参数驱动的三甘醇脱水装置三甘醇损耗量在线预测方法

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5436852A (en) * 1990-09-21 1995-07-25 Yamatake-Honeywell Co., Ltd. Method and apparatus for calculating predicted mean thermal sensitivity
JP3203126B2 (ja) * 1994-04-19 2001-08-27 三洋電機株式会社 空気調和機の制御装置
US6145751A (en) * 1999-01-12 2000-11-14 Siemens Building Technologies, Inc. Method and apparatus for determining a thermal setpoint in a HVAC system
US6702456B2 (en) * 2002-04-04 2004-03-09 China Textile Institute Method and apparatus for evaluating the phase changing function of fabric
EP1989490B1 (en) * 2006-02-10 2011-12-21 Danfoss A/S Method and system for controlling the climate in a house
US8005655B2 (en) * 2008-02-26 2011-08-23 Kimberly-Clark Worldwide, Inc. Thermal comfort model having multiple fabric layers
US8577650B2 (en) * 2008-02-26 2013-11-05 Kimberly-Clark Worldwide, Inc. User interface for modeling thermal comfort
JP2013057476A (ja) * 2011-09-09 2013-03-28 Toshiba Corp Ｐｍｖ推定装置およびそのプログラム
IT201800020629A1 (it) * 2018-12-21 2020-06-21 Thinxs S R L Dispositivo per la misura della temperatura percepita

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4058254A (en) * 1968-01-26 1977-11-15 Danfoss A/S Temperature-monitoring and control system
DD105516A1 (ja) * 1973-05-23 1974-04-20
US4125012A (en) * 1975-06-25 1978-11-14 Madsen Thomas L Apparatus for measuring thermal discomfort originating from asymmetry in the thermal field or variations with time of the thermal influence on the skin
CA1035596A (en) * 1975-09-23 1978-08-01 Lloyd E. Machattie Still shade temperature measurement
JPS54139251A (en) * 1978-04-21 1979-10-29 Hitachi Ltd Method of controlling temperature and humidity
US4504157A (en) * 1982-09-30 1985-03-12 Austin Crabtree Chill temperature meter
WO1986004674A1 (en) * 1985-02-06 1986-08-14 Daikin Industries, Ltd. Heat detecting element
US4727500A (en) * 1985-05-01 1988-02-23 Sherwood Medical Company Electronic thermometer with fixed response time
US4702617A (en) * 1986-03-10 1987-10-27 Crabtree Austin A Chill temperature meter with selectable response
JPS6365318A (ja) * 1986-09-05 1988-03-23 Daikin Ind Ltd 温熱検知素子
US4799176A (en) * 1986-12-29 1989-01-17 Harper-Wyman Company Electronic digital thermostat
US4863279A (en) * 1988-02-22 1989-09-05 Morris L. Markel Operative temperature sensing system
DE8908474U1 (de) * 1989-07-12 1989-09-21 Katz, Philipp, Prof. Dr., 6300 Gießen Vorrichtung zur Messung der Behaglichkeit
JPH05332596A (ja) * 1990-07-30 1993-12-14 Yamatake Honeywell Co Ltd 等価温度演算方法および環境計測装置
JP2588792B2 (ja) * 1990-09-21 1997-03-12 山武ハネウエル株式会社 温熱感覚演算方法および装置

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0618098A1 (en) * 1993-04-01 1994-10-05 FIAT AUTO S.p.A. A process and system for controlling environmental conditions
JPH08100941A (ja) * 1994-03-09 1996-04-16 Lg Electronics Inc 空調システムのｐｍｖ算出方法
JP2008089343A (ja) * 2006-09-29 2008-04-17 Toyota Motor Corp 接触温熱感計測用センサ装置および接触温熱感計測装置
CN109974244A (zh) * 2014-02-17 2019-07-05 松下电器产业株式会社 空调机的控制方法
CN113088359A (zh) * 2021-03-30 2021-07-09 重庆大学 一种工艺参数驱动的三甘醇脱水装置三甘醇损耗量在线预测方法

Also Published As

Publication number Publication date
DK0477796T3 (da) 1996-01-02
US5333953A (en) 1994-08-02
JP2581625B2 (ja) 1997-02-12
DE69112162T2 (de) 1996-05-09
EP0477796B1 (en) 1995-08-16
EP0477796A3 (en) 1992-08-12
DE69112162D1 (de) 1995-09-21
EP0477796A2 (en) 1992-04-01

Similar Documents

Publication Publication Date Title
JP6505514B2 (ja) 空気調和機、センサシステムおよびその温冷感推定方法
US5436852A (en) Method and apparatus for calculating predicted mean thermal sensitivity
JPH0518592A (ja) 温熱感覚演算方法および装置、予測平均温感演算方法および装置
Mayer Objective criteria for thermal comfort
CN107102022A (zh) 基于暖体假人的热环境舒适性评价方法
US4433923A (en) Operative temperature sensing system
JP2588792B2 (ja) 温熱感覚演算方法および装置
JPH0577931B2 (ja)
JPH10103745A (ja) 温熱環境検出装置
WO1992002768A1 (en) Method and device for computing estimated mean temperature sensation
WO1992002767A1 (en) Method of computing equivalent temperature and instrument for environment measurement
KR101137662B1 (ko) 공조 시스템의 ｐｍｖ제어방법
JPH08100941A (ja) 空調システムのｐｍｖ算出方法
JP2734193B2 (ja) 空気調和機の制御装置
JPH0749141A (ja) 室内温熱環境の自動制御方法
Foda Evaluating local and overall thermal comfort in buildings using thermal manikins
JPH0713573B2 (ja) 予測平均温感演算方法および装置
JPS6012569B2 (ja) 熱環境測定器
JPH09280907A (ja) 温熱環境計測方法
US20240230429A1 (en) Thermal comfort measuring system
JP4309604B2 (ja) 空調設備の空調制御装置
JPH051842A (ja) 温熱感覚演算装置
CN114199603A (zh) 一种家用电器舒适性模拟自然风测试装置及测试方法
JPH0798142A (ja) 温感評価装置
Cena Environmental heat loss

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20071121

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20081121

Year of fee payment: 12

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20081121

Year of fee payment: 12

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091121

Year of fee payment: 13

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091121

Year of fee payment: 13

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101121

Year of fee payment: 14

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101121

Year of fee payment: 14

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111121

Year of fee payment: 15

EXPY Cancellation because of completion of term
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111121

Year of fee payment: 15