(57)【要約】【課題】人体熱モデルに周期的変動風の風速を入力して人の生理状態を予測することによって、周期的変動風全体の総合的な温冷感を適切に反映することができるようにする。【解決手段】周期的変動風の風速を人体熱モデルに入力して体感温度を算出する工程と、前記体感温度から平均温冷感予測を算出する工程と、前記平均温冷感予測から不満足率予測を算出する工程と、前記不満足率予測から空調強度又は周期的変動風発生器の強度の補正量を決定する工程と、を含む。【選択図】図１
2018075913.tif 000002
【特許請求の範囲】
【請求項１】
周期的変動風の風速を人体熱モデルに入力して体感温度を算出する工程と、
前記体感温度から平均温冷感予測を算出する工程と、
前記平均温冷感予測から不満足率予測を算出する工程と、
前記不満足率予測から空調強度又は周期的変動風発生器の強度の補正量を決定する工程と、
を含むことを特徴とする周期的変動風を利用する車内空調方法。
【請求項２】
前記人体熱モデルに入力される周期的変動風の風速は、周期的変動風体感時と予測平均皮膚温が同等となる定常風の風速に換算された定常風等価換算風速である請求項１に記載の周期的変動風を利用する車内空調方法。
【請求項３】
前記定常風等価換算風速は、周期的変動風体感時の予測平均皮膚温と定常風体感時の予測平均皮膚温との誤差が極小となる定常風の風速である請求項２に記載の周期的変動風を利用する車内空調方法。
【請求項４】
前記周期的変動風発生器の強度が所定の閾値より小さい場合、前記人体熱モデルに入力される周期的変動風の風速は、実測された風速である請求項１に記載の周期的変動風を利用する車内空調方法。
【請求項５】
前記不満足率予測から空調強度又は周期的変動風発生器の強度の補正量を決定する工程では、できる限り空調強度を低く抑えるように補正量を決定する請求項１〜４のいずれか１項に記載の周期的変動風を利用する車内空調方法。
【請求項６】
前記周期的変動風発生器は、往復しながら送風を行う横流ファンである請求項１〜５のいずれか１項に記載の周期的変動風を利用する車内空調方法。
【請求項７】
外気温及び乗客が発生する熱量から、温冷感変動を抑制する空調強度の他の補正量を計算する工程と、
前記空調強度の他の補正量と、前記不満足率予測から決定された空調強度又は周期的変動風発生器の強度の補正量とを加算する工程とを更に含む請求項１〜６のいずれか１項に記載の周期的変動風を利用する車内空調方法。
【請求項８】
車両の客室内の空調を行う空調装置と、前記客室内に周期的変動風を発生させる周期的変動風発生器と、前記空調装置及び周期的変動風発生器の制御を行う制御器とを有する周期的変動風を利用する車内空調システムであって、
周期的変動風の風速が入力される人体熱モデルを使用して体感温度を算出する人体熱モデル部と、
前記体感温度から平均温冷感予測を算出する平均温冷感予測部と、
前記平均温冷感予測から不満足率予測を算出する熱的不満足率予測部と、
前記不満足率予測から空調強度又は周期的変動風発生器の強度の補正量を決定する補正量決定部と、
を備えることを特徴とする周期的変動風を利用する車内空調システム。
【請求項９】
周期的変動風の風速を、周期的変動風体感時と予測平均皮膚温が同等となる定常風の風速である定常風等価換算風速に換算して前記人体熱モデルに入力する定常風換算部を更に備える請求項８に記載の周期的変動風を利用する車内空調システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、周期的変動風を利用する車内空調方法及びシステムに関するものである。
【背景技術】
【０００２】
従来、複数人が滞在する空間では、冷房の補助として、扇風機や送風機による送風が利用されている。また、より多くの人に風が当たるように、扇風機や送風機の首振りや回転が行われるので、各人が受ける風は、首振りや回転の周期に対応する変動風、すなわち、周期的変動風となる（例えば、特許文献１参照。）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開平９−２７３８００号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、前記従来の技術では、室内全体の空気流を変動させて人の涼感を高めるだけであって、周期的変動風を体感した際の人の生理状態を十分に考慮していない。
【０００５】
ここでは、前記従来の技術の問題点を解決して、人体熱モデルに周期的変動風の風速を入力して人の生理状態を予測することによって、周期的変動風全体の総合的な温冷感を適切に反映することができる周期的変動風を利用する車内空調方法及びシステムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
そのために、周期的変動風を利用する車内空調方法においては、周期的変動風の風速を人体熱モデルに入力して体感温度を算出する工程と、前記体感温度から平均温冷感予測を算出する工程と、前記平均温冷感予測から不満足率予測を算出する工程と、前記不満足率予測から空調強度又は周期的変動風発生器の強度の補正量を決定する工程と、を含む。
【０００７】
他の周期的変動風を利用する車内空調方法においては、さらに、前記人体熱モデルに入力される周期的変動風の風速は、周期的変動風体感時と予測平均皮膚温が同等となる定常風の風速に換算された定常風等価換算風速である。
【０００８】
更に他の周期的変動風を利用する車内空調方法においては、さらに、前記定常風等価換算風速は、周期的変動風体感時の予測平均皮膚温と定常風体感時の予測平均皮膚温との誤差が極小となる定常風の風速である。
【０００９】
更に他の周期的変動風を利用する車内空調方法においては、さらに、前記周期的変動風発生器の強度が所定の閾値より小さい場合、前記人体熱モデルに入力される周期的変動風の風速は、実測された風速である。
【００１０】
更に他の周期的変動風を利用する車内空調方法においては、さらに、前記不満足率予測から空調強度又は周期的変動風発生器の強度の補正量を決定する工程では、できる限り空調強度を低く抑えるように補正量を決定する。
【００１１】
更に他の周期的変動風を利用する車内空調方法においては、さらに、前記周期的変動風発生器は、往復しながら送風を行う横流ファンである。
【００１２】
更に他の周期的変動風を利用する車内空調方法においては、さらに、外気温及び乗客が発生する熱量から、温冷感変動を抑制する空調強度の他の補正量を計算する工程と、前記空調強度の他の補正量と、前記不満足率予測から決定された空調強度又は周期的変動風発生器の強度の補正量とを加算する工程とを更に含む。
【００１３】
周期的変動風を利用する車内空調システムにおいては、車両の客室内の空調を行う空調装置と、前記客室内に周期的変動風を発生させる周期的変動風発生器と、前記空調装置及び周期的変動風発生器の制御を行う制御器とを有する周期的変動風を利用する車内空調システムであって、周期的変動風の風速が入力される人体熱モデルを使用して体感温度を算出する人体熱モデル部と、前記体感温度から平均温冷感予測を算出する平均温冷感予測部と、前記平均温冷感予測から不満足率予測を算出する熱的不満足率予測部と、前記不満足率予測から空調強度又は周期的変動風発生器の強度の補正量を決定する補正量決定部と、を備える。
【００１４】
他の周期的変動風を利用する車内空調システムにおいては、さらに、周期的変動風の風速を、周期的変動風体感時と予測平均皮膚温が同等となる定常風の風速である定常風等価換算風速に換算して前記人体熱モデルに入力する定常風換算部を更に備える。
【発明の効果】
【００１５】
本開示によれば、人体熱モデルに周期的変動風の風速を入力して人の生理状態を予測する。これにより、周期的変動風全体の総合的な温冷感を適切に反映することができる。
【図面の簡単な説明】
【００１６】
【図１】本実施の形態における車内空調システムが搭載された車両を示す模式断面図である。
【図２】本実施の形態における車内空調システムの要部の機能構成を示すブロック図である。
【図３】本実施の形態における車内空調システムの第１補正量を得るための動作を説明する図である。
【図４】本実施の形態における定常風換算部の動作を説明する図である。
【図５】本実施の形態における定常風等価換算風速の計算例を示す図である。
【図６】本実施の形態における周期的変動風の風速の実測例を示す図である。
【図７】本実施の形態におけるモデルパラメータの例を示す図である。
【図８】本実施の形態における熱的不満足率を説明する図である。
【図９】本実施の形態における補正対象決定部の動作を説明する図である。
【図１０】本実施の形態におけるファジールールの例を説明する図である。
【発明を実施するための形態】
【００１７】
以下、本実施の形態について図面を参照しながら詳細に説明する。
【００１８】
図１は本実施の形態における車内空調システムが搭載された車両を示す模式断面図である。なお、図において、（ａ）は側断面図、（ｂ）は横断面図である。
【００１９】
図において、１１は、本実施の形態における車内空調システムが搭載された車両であって、例えば、通勤列車等の列車、乗り合いバス、モノレール等の車両であって、いかなる種類の車両であってもよいが、ここでは、説明の都合上、通勤列車である場合についてのみ説明する。また、２１は、本実施の形態における車内空調システムの空調装置であって、例えば、冷媒圧縮機、蒸発器、凝縮器、ファン等を備え、冷房、暖房、除湿、通風等の空調を行う装置であるが、ここでは、説明の都合上、冷房を行う場合についてのみ説明する。さらに、前記空調装置２１は、必ずしも車両１１の頂部に搭載される必要はなく、床下等のいかなる箇所に搭載されていてもよいが、ここでは、説明の都合上、車両１１の頂部に搭載されている場合についてのみ説明する。
【００２０】
なお、本実施の形態において、車両１１、空調装置２１及びその他の装置の各部の構成及び動作を説明するために使用される上、下、左、右、前、後等の方向を示す表現は、絶対的なものでなく相対的なものであり、前記車両１１、空調装置２１及びその他の装置の各部が図に示される姿勢である場合に適切であるが、その姿勢が変化した場合には姿勢の変化に応じて変更して解釈されるべきものである。
【００２１】
前記車両１１は、図示されない乗客が乗車する客室１２と、該客室１２の天井１６の上方に位置する天井内空間１３とを備える。そして、空調装置２１は、矢印で示されるように、空調風としての冷却風２７を天井内空間１３内に送出する。すると、前記冷却風２７は、天井１６に形成された複数の吹出開口１７を通って、客室１２内に流入する。なお、該客室１２内の空気は、天井１６に形成されたリターン開口１８を通って天井内空間１３内に流入し、空調装置２１に吸引される。
【００２２】
また、２３は、本実施の形態における車内空調システムの横流ファンであり、天井１６に複数個配設され、矢印２８で示される角度範囲内で往復しながら、送風を行う。前記横流ファン２３は、別名、クロスファン、ラインフローファン、タンジェンシャルファン、又は、貫流ファンとも呼ばれ、客室１２内の空気を攪拌するための送風機であって、客室１２内に周期的変動風を発生させる周期的変動風発生器として機能する。なお、該周期的変動風発生器は、客室１２内に周期的変動風を発生させる機器であればいかなるタイプの機器であってもよく、例えば、首振り可能乃至回転可能な扇風機であってもよく、必ずしも天井１６に配設される必要はなく、横流ファン２３である必要もなく、複数でなく単数であってもよいが、ここでは、説明の都合上、天井１６に複数個配設された横流ファン２３である場合についてのみ説明する。
【００２３】
さらに、２２は、本実施の形態における車内空調システムの制御器であり、例えば、ＣＰＵ等の演算装置、半導体メモリ、ハードディスク等の記憶装置、有線又は無線の通信装置等を備える一種のコンピュータであって、前記記憶装置にインストールされたアプリケーションソフトウェア等のプログラムに従って動作し、前記空調装置２１及び横流ファン２３の制御を行う。なお、図に示される例において、制御器２２は、空調装置２１に一体的に取り付けられているが、必ずしも、空調装置２１に一体的に取り付けられる必要はなく、空調装置２１から離れた箇所に配設されていてもよい。
【００２４】
図に示される例において、前記客室１２内には、床１５の近傍の温度、湿度、風速、放射等を計測する床部センサ２５ａや、壁１４の近傍の温度、湿度、風速、放射等を計測する壁部センサ２５ｂが配設され、天井内空間１３内には、リターン開口１８の近傍の温度、湿度、風速、放射等を計測する天井部センサ２５ｃが配設されている。なお、前記床部センサ２５ａ、壁部センサ２５ｂ及び天井部センサ２５ｃを統合的に説明する場合には、センサ２５として説明する。本実施の形態における車内空調システムでは、前記センサ２５の数、位置等は、図に示される例に限定されるものでなく、いくつであってもよいし、いかなる位置に配設されていてもよい。
【００２５】
次に、前記構成の車内空調システムの構成について説明する。
【００２６】
図２は本実施の形態における車内空調システムの要部の機能構成を示すブロック図である。
【００２７】
図において、３１は、前記制御器２２が備える情報取得部であり、センサ２５及び設定情報発信部２６から、客室１２内外の温度、湿度等の情報を含む環境情報や、乗客の人数、性別、服装等の情報を含む乗客情報を取得する。なお、前記設定情報発信部２６は、乗客情報に加え、事前の設定値やその他の情報を前記情報取得部３１に送信する。３２は、前記制御器２２が備える温熱環境取得部であり、後述される定常風換算部３２ａ及び登録風速データベース３２ｂを含み、前記情報取得部３１が取得した情報の中から、客室１２内の温熱環境である車内温熱環境、乗客の平均代謝量、乗客の平均着衣量、乗客の男女比、風速等の情報を取得する。また、３３は、前記制御器２２が備える平均温冷感計算部であり、後述される人体熱モデル部３３ａ及び平均温冷感予測部３３ｂを含み、前記温熱環境取得部３２が取得した情報に基づき、平均温冷感予測（ＤｙＰＭＶ）を算出する。さらに、３４は、前記制御器２２が備える第１補正量計算部であり、後述される熱的不満足率予測部３４ａ及び第１補正量決定部３４ｂを含み、前記平均温冷感計算部３３から、平均温冷感予測値を取得すると、乗客の熱的不満足率を推定し、第１補正量を算出する。なお、前記温熱環境取得部３２、平均温冷感計算部３３及び第１補正量計算部３４が行う動作は、現在の乗客の平均温冷感を最適値に近付けるための補正である。
【００２８】
３５は、前記制御器２２が備える外気温及び乗車人数取得部であり、前記情報取得部３１が取得した情報の中から、客室１２外の温熱環境である外気温、乗客の人数等の情報を取得する。また、３６は、前記制御器２２が備える第２補正量計算部であり、前記外気温及び乗車人数取得部３５が取得した情報に基づき、第２補正量を算出する。なお、前記外気温及び乗車人数取得部３５並びに第２補正量計算部３６が行う動作は、外気温及び乗客が発生する熱量から温冷感変動を抑制するための補正であって、今後の温冷感の変動を抑制するための補正である。
【００２９】
そして、３７は、前記制御器２２が備える補正量取得部であり、前記第１補正量計算部３４が算出した第１補正量と、前記第２補正量計算部３６が算出した第２補正量とを加算した値を、空調装置２１の制御量を補正するための補正量として取得する。
【００３０】
次に、前記構成の車内空調システムの動作について説明する。ここでは、前述のように、空調として冷房を行う場合についてのみ説明する。
【００３１】
図３は本実施の形態における車内空調システムの第１補正量を得るための動作を説明する図、図４は本実施の形態における定常風換算部の動作を説明する図、図５は本実施の形態における定常風等価換算風速の計算例を示す図、図６は本実施の形態における周期的変動風の風速の実測例を示す図、図７は本実施の形態におけるモデルパラメータの例を示す図、図８は本実施の形態における熱的不満足率を説明する図、図９は本実施の形態における補正対象決定部の動作を説明する図、図１０は本実施の形態におけるファジールールの例を説明する図である。なお、図８において、（ａ）は季節調整パラメータの変化を示す図、（ｂ）は不満足率の変化を示す図である。
【００３２】
本実施の形態における制御器２２は、制御変数を乗客の平均温冷感とし、該平均温冷感が乗客の熱的不満足率が最小となる最適平均温冷感になるように、冷房強度及び／又は横流ファン２３の強度である横流ファン強度を補正する。なお、最適平均温冷感は、車内の乗客の男女比や季節によって異なる。そこで、前記温熱環境取得部３２、平均温冷感計算部３３及び第１補正量計算部３４が行う現在の乗客の平均温冷感を最適値に近付けるための補正では、乗客の性差や季節差を考慮した上で、現在の平均温冷感を最適平均温冷感に近付けるための冷房強度及び／又は横流ファン強度の補正量を計算する。また、外気温が車内温度より高い場合、熱が客室１２内に貫流する。さらに、乗客も熱を発生する。そこで、前記外気温及び乗車人数取得部３５並びに第２補正量計算部３６が行う今後の温冷感の変動を抑制するための補正では、これらの熱負荷による温冷感変動を抑制するための冷房強度補正量を計算する。
【００３３】
既存の温熱指標として、ＩＳＯ７７３０で採用されている予測平均温冷感（ＰＭＶ：Ｐｒｅｄｉｃｔｅｄ  Ｍｅａｎ  Ｖｏｔｅ）、及び、予測不満足率（ＰＰＤ：Ｐｒｅｄｉｃｔｅｄ  Ｐｅｒｃｅｎｔａｇｅ  ｏｆ  Ｄｉｓｓａｔｉｓｆｉｅｄ）がある。ＰＭＶは、周囲温熱環境による人体の熱負荷量を熱平衡式によって計算され、かかる熱負荷量と被験者実験での温冷感申告値とを結び付けた回帰式により、計算される。ＰＰＤは、あるＰＭＶに対して何割の人が不満足と感じるかを予測するもので、ＰＭＶと被験者実験での不満足者の割合とを結び付けた回帰式によって、計算される。ただし、ＰＭＶ及びＰＰＤの適用範囲は定常環境であり、ＰＭＶとＰＰＤとの関係には温熱快適性の季節差は考慮されていない。
【００３４】
図３には、本実施の形態における現在の乗客の平均温冷感を最適値に近付けるための補正に用いられる快適性指標の計算モデルが、制御器２２の平均温冷感計算部３３及び第１補正量計算部３４の機能構成とともに示されている。前記計算モデルでは、前述のＰＭＶ及びＰＰＤの計算方法を踏襲して、体感温度（ＤｙＳＥＴ* ）から平均温冷感予測（ＤｙＰＭＶ）を算出する計算と、平均温冷感予測（ＤｙＰＭＶ）から不満足率予測（ＤｙＰＰＤ）を算出する計算との２段階の計算を行う。本実施の形態においては、体感温度（ＤｙＳＥＴ* ）から平均温冷感予測（ＤｙＰＭＶ）を算出する計算に人体熱モデルを使用することによって、適用範囲を非定常温熱環境に拡張し（例えば、特許文献２参照。）、さらに、平均温冷感予測（ＤｙＰＭＶ）から不満足率予測（ＤｙＰＰＤ）を算出する計算に、各季節で実施した被験者実験データを反映させることによって、温熱快適性の季節特性を反映させる。ちなみに、特許文献２では、夏季にのみ実施された被験者実験データが用いられている。
【特許文献２】特開２０１４−２２８１７２号公報
【００３５】
現在の乗客の平均温冷感を最適値に近付けるための補正では、図３に示されるように、平均温冷感計算部３３は、人体熱モデル（例えば、非特許文献１参照。）を使用して、体感温度（ＤｙＳＥＴ* ）を算出する人体熱モデル部３３ａと、平均温冷感予測式に従って平均温冷感予測（ＤｙＰＭＶ）を算出する平均温冷感予測部３３ｂとを備える。
【非特許文献１】田辺他、「温熱環境評価のための６５分割体温調整モデルに関する研究」、日本建築学会計画系論文集、Ｎｏ．５４１、ｐｐ．９−１６、２００１
【００３６】
前記人体熱モデル部３３ａは、非定常温熱環境下での体温調節系を有する人体熱モデルを使用して、環境−人体間及び人体内の熱移動計算を行う。なお、前記人体熱モデルには、温度（Ｔａ（ｔ））、湿度（ＲＨ（ｔ））、風速、放射（Ｔｒ（ｔ））、平均着衣量（Ｉｃｌｏ）、及び、平均代謝量（Ｍ（ｔ））が入力される。該平均代謝量（Ｍ（ｔ））は、重量及び乗車人数に基づいて算出される。
【００３７】
ＳＥＴ* は、「温熱生理状態（皮膚温、発汗）及び放熱量が実在環境におけるものと同等になるような相対湿度５０〔％〕の標準環境の気温」（単位は〔℃〕）として定義される温熱指標である。本実施の形態における体感温度（ＤｙＳＥＴ* ）は、温熱環境の非定常性を考慮したもので、非特許文献１でも紹介されている公知の指標である。
【００３８】
平均着衣量（Ｉｃｌｏ）は、例えば、次の式（１）によって算出することができる（例えば、非特許文献２参照。）。
【００３９】
【数１】
2018075913.tif 000003
【００４０】
ここで、平均年齢は、線区又は時間帯による利用者の平均年齢の情報を利用することができる場合には該情報に基づいて決定される。未知の場合には、例えば、男性は４０歳、女性は３０歳と仮定する。また、日平均気温には、過去数年の同日の気温の平均値や当日の気象予報による気温などを用いる。
【非特許文献２】大井元、石井昭夫、斎藤基之、城野修、「屋外における着衣に関するアンケート調査  その３：着衣量推定式の構築」、日本建築学会大会学術講演梗概集、ｐｐ．４２５−４２６、２０００
【００４１】
また、平均代謝量（Ｍ（ｔ））は、例えば、次の式（２）によって算出することができる。
【００４２】
【数２】
2018075913.tif 000004
【００４３】
ここで、乗車人数が座席数以下の場合には、全員が座っていると仮定して、座位安静時に相当する代謝量である１．０〔ｍｅｔ〕とする。また、乗車人数が座席数より多い場合には、立位乗客（代謝量は１．３〔ｍｅｔ〕と仮定）を想定し、各々の人数で重付けした乗客の代謝量の平均値を算出する。
【００４４】
ところで、温冷感は、風速変動に対応して時々刻々と変動するので、周期的変動風全体の総合的な温冷感（１周期以上体感した場合の総合的な温冷感）を一意に決定することはできない。しかし、周期的変動風の人体への影響については、ある時刻の一瞬の温冷感というよりも、周期的変動風全体の総合的な温冷感を評価乃至予測することが有用な場合が多いと考えられる。そこで、本実施の形態においては、周期的変動風全体の総合的な温冷感を評価乃至予測するようになっている。
【００４５】
そのために、温熱環境取得部３２は、図３に示されるような定常風換算部３２ａ及び登録風速データベース３２ｂを含んでいる。該登録風速データベース３２ｂには、あらかじめ、横流ファン２３の送風の強度設定（横流ファン強度（Ｎ））毎に、周期的変動風の風速の時系列データが格納されて登録されている。例えば、事前の風速の測定結果等に基づいて、実測データ又はそれを周期関数等で近似した模擬波形データを登録風速データとして登録する。当該時系列データは時間長さＮｔ〔ｓ〕の時系列データである。各部位の重み値としては、その表面積、温冷感への寄与度等の適当な値を用いることができる。なお、定常風換算部３２ａは、図４に示されるように、等価換算風速計算部３２ｃを含み、該等価換算風速計算部３２ｃが定常風等価換算風速を算出して、平均温冷感計算部３３に送信する。
【００４６】
定常風換算部３２ａは、次のステップＳ１〜Ｓ４に示されるような動作を行って、定常風等価換算風速を出力する。なお、横流ファン強度（Ｎ）が低い場合には、周期的変動風がほとんど感じられないことが想定され、定常風と近似することができるので、定常風等価換算風速の算出を実行する横流ファン強度（Ｎ）の閾値Ｎmin があらかじめ設定されている。すなわち、図３に示されるように、横流ファン強度（Ｎ）が所定の閾値Ｎmin より小さい場合、定常風換算部３２ａは定常風等価換算風速を算出せず、実測された風速がそのまま平均温冷感計算部３３の人体熱モデル部３３ａに入力されて使用される。なお、横流ファン強度（Ｎ）には、１〜Ｎmax 段階の設定があるものとする。そして、横流ファン強度（Ｎ）が閾値Ｎmin 段階以上の場合に、定常風等価換算風速の算出が実行される。
【００４７】
ステップＳ１：定常風換算部３２ａは、横流ファン強度（Ｎ）に対応した登録風速データ（登録風速データベース３２ｂに登録されている風速のデータ）を人体熱モデル部３３ａの人体熱モデルに入力する。前記登録風速データは、前述のように、実測データ又はそれを周期関数等で近似した模擬波形データであって、時間長さはＮｔ〔ｓ〕である。また、風速以外の温熱環境値、すなわち、風速以外の温熱要素のデータである温度（Ｔａ（ｔ））、湿度（ＲＨ（ｔ））、放射（Ｔｒ（ｔ））、平均着衣量（Ｉｃｌｏ）、及び、平均代謝量（Ｍ（ｔ））が、人体熱モデルに入力される。これらのデータは、計算開始時点の測定値（又は、現時点から過去Ｔｐ〔ｓ〕間の平均値）が、Ｎｔ〔ｓ〕間継続するものと仮定した時系列データである。すると、人体熱モデルは、予測平均皮膚温（Ｔｓｋ＿ｖ）を出力する。
【００４８】
ステップＳ２：定常風換算部３２ａは、定常風風速初期値（Ｖｃ（０））を人体熱モデル部３３ａの人体熱モデルに入力する。前記定常風風速初期値（Ｖｃ（０））は、Ｖｃ（０）がＮｔ〔ｓ〕間継続するものと仮定した時系列データである。例えば、Ｖｃ（０）＝０．０５〔ｍ／ｓ〕である。また、前記ステップＳ１で説明したものと同様の風速以外の温熱要素のデータが、人体熱モデルに入力される。すると、該人体熱モデルは、予測平均皮膚温（Ｔｓｋ＿ｃ（０））を出力する。そして、定常風換算部３２ａは、変動風体感時の予測平均皮膚温と定常風体感時の予測平均皮膚温との平均二乗誤差として、Ｔｓｋ＿ｃ（０）とＴｓｋ＿ｖとの平均二乗誤差（ＲＭＳ（０））を、次の式（３）によって算出する。
【００４９】
【数３】
2018075913.tif 000005
【００５０】
ステップＳ３：定常風換算部３２ａの等価換算風速計算部３２ｃは、定常風風速をｄＶ〔ｍ／ｓ〕（例えば、０．０５〔ｍ／ｓ〕）刻みで上昇させていき、Ｔｓｋ＿ｃ（ｎ）とＴｓｋ＿ｖとの平均二乗誤差（ＲＭＳ（ｎ））を、次の式（４）によって算出する。
【００５１】
【数４】
2018075913.tif 000006
【００５２】
ここで、ｎは定常風風速初期値（Ｖｃ（０））からの定常風風速の上昇回数である。ｎ回目の定常風風速（Ｖｃ（ｎ））は、Ｖｃ（ｎ）＝Ｖｃ（０）＋ｎ×ｄＶである。そして、等価換算風速計算部３２ｃは、ｎ回目の平均二乗誤差（ＲＭＳ（ｎ））とｎ−１回目の平均二乗誤差（ＲＭＳ（ｎ−１））との差（ｄＲＭＳ（ｎ））が０以上となったら（平均二乗誤差が減少から増加に転じたら）、計算を終了し、平均二乗誤差が減少から増加に転じる直前の定常風風速（Ｖｃ（ｎ−１））を定常風等価換算風速とする。つまり、周期的変動風体感時の予測平均皮膚温と定常風体感時の予測平均皮膚温との誤差が極小となる定常風の風速を定常風等価換算風速とする。これにより、周期的変動風体感時と予測平均皮膚温が同等となる定常風風速を決定することができる。
【００５３】
ステップＳ４：定常風換算部３２ａの等価換算風速計算部３２ｃは、前記ステップＳ３で決定した定常風等価換算風速（Ｖｃ（ｎ−１））を平均温冷感計算部３３の人体熱モデル部３３ａに入力する。これにより、周期的変動風全体の総合的な温冷感の予測値を得ることができる。
【００５４】
図５には、定常風換算部３２ａの等価換算風速計算部３２ｃによって算出された定常風等価換算風速の例が示されている。これは、気温３０〔℃〕、湿度７０〔％ｒｈ〕に制御された人工気候室内に６０分間滞在した後に、図６に示されるような風速の実測波形を有する周期的変動風（周期２０秒、風速ピーク値２．０〔ｍ／ｓ〕）を体感した場合の例である。図５には、平均皮膚温の実測値（被験者８名の平均値）、周期的変動風風速を人体熱モデルに入力した場合に出力される予測平均皮膚温（Ｔｓｋ＿ｖ）の例、前記ステップＳ２で説明した定常風風速を人体熱モデルに入力した場合に出力される予測平均皮膚温（Ｔｓｋ＿ｃ）の例、及び、最終的に得られた定常風等価換算風速１．１５〔ｍ／ｓ〕の場合の予測平均皮膚温（Ｔｓｋ＿ｃ）の例が示されている。周期的変動風に対する平均皮膚温の計算では、図６に示されるような風速の実測波形に基づき、最大値２．０〔ｍ／ｓ〕、最小値０．２〔ｍ／ｓ〕、周期２０秒となるｓｉｎ関数で近似した模擬波形データが、人体熱モデルの入力風速として使用された。
【００５５】
ここで、定常風等価換算風速とは、周期的変動風体感時と平均皮膚温が同等となる定常風風速として定義されたものであり、周期的変動風体感時の温冷感を、瞬時的なレベルではなく、複数周期体感した場合の総合的な温冷感を予測するための仮想風速である。周期的変動風の総合的な温冷感を予測するための風速換算方法として、変動風風速の時系列に対する統計値（平均値、標準偏差等）を採用する方法も考えられるが、それらは生理的根拠に乏しい。本実施の形態による方法は、人体熱モデルの生理状態予測に基づく風速換算方法である。
【００５６】
そして、平均温冷感予測部３３ｂは、体感温度（ＤｙＳＥＴ* ）から、次の式（５）のような平均温冷感予測式に従って、平均温冷感予測（ＤｙＰＭＶ）を算出する。
【００５７】
【数５】
2018075913.tif 000007
【００５８】
ここで、ｐ0 、ｐ1 及びｐ2 は、被験者実験データによって導出されるモデルパラメータである。また、前記式（５）は、男女別に用意され、男女別の平均温冷感予測（ＤｙＰＭＶ）が算出される。
【００５９】
算出された平均温冷感予測（ＤｙＰＭＶ）の値は、多数の人が感じる温冷感の平均値であって、−４：寒い、−３：やや寒い、−２：涼しい、−１：やや涼しい、０：中立、１：やや暖かい、２：暖かい、３：やや暑い、４：暑い、に、それぞれ、相当する。
【００６０】
そして、第１補正量計算部３４は、熱的不満足率予測式に従って不満足率（ＤｙＰＰＤ）を算出する熱的不満足率予測部３４ａと、ファジールールに従って空調強度の補正量としての第１補正量を決定する第１補正量決定部３４ｂとを備える。該第１補正量決定部３４ｂは、平均温冷感最適補正を冷房強度の補正で行うか、横流ファン強度の補正で行うかを決定する補正対象決定部３４ｃと、補正量決定関数を適用して第１補正量を決定する補正量決定関数適用部３４ｄとを含んでいる。
【００６１】
前記熱的不満足率予測部３４ａは、乗客の性差及び季節差を考慮した上で、乗客の熱的不満足率を推定する。具体的には、前記熱的不満足率予測部３４ａは、平均温冷感計算部３３が算出した平均温冷感予測（ＤｙＰＭＶ）を使用し、次の式（６）〜（８）のような熱的不満足率予測式に従って、寒い不満足率（ＤｙＰＰＤ＿ｃ又はＤｙＰＰＤ＿ｃｏｌｄ）、暑い不満足率（ＤｙＰＰＤ＿ｈ又はＤｙＰＰＤ＿ｈｏｔ）、及び、合計の不満足率（ＤｙＰＰＤ）を算出する。
【００６２】
【数６】
2018075913.tif 000008
【００６３】
ここで、寒い不満足率（ＤｙＰＰＤ＿ｃｏｌｄ）は乗客全体のうちで寒くて不満足と感じる人の割合であり、暑い不満足率（ＤｙＰＰＤ＿ｈｏｔ）は乗客全体のうちで暑くて不満足と感じる人の割合である。また、ａh 、ａc 、ｂh 、ｂc は、被験者実験データによって導出されるモデルパラメータである。
【００６４】
さらに、ｂc-season及びｂh-seasonは、季節差を考慮可能とするために、乗客の温熱的な生理及び心理状態の季節変化を総合的に表現する季節調整パラメータであって、日付を変数とする周期関数として、次の式（９）〜（１２）で定義される。
【００６５】
【数７】
2018075913.tif 000009
【００６６】
ここで、ｄａｙは、４月１日を第１日（ｄａｙ  ｏｎｅ）とした場合の日付であって、ｄａｙ＝１〜３６５である。そして、Ｚ0 は、４月１日から始まる１年の前期（前半）においては１となり、後期（後半）においては０となるパラメータであり、Ｚ1 は、その逆の値となるパラメータである。
【００６７】
図７には、前記式（５）〜（１２）における各種モデルパラメータの男女別の数値例が示されている。図７に示される例は、着衣及び姿勢を統一した優等車両実験のデータを使用し、全季節の実験結果に最もよく適合するものを選定して同定した結果である。なお、線形モデルは最小二乗法、ロジスティックモデルは最尤法によって同定した。
【００６８】
人の体温調節系は、春から夏にかけては、体表面からの熱放散を促すように血管拡張が進み、発汗量も増加する傾向にある。一方、秋から冬にかけては、体表面からの熱放散を抑制するように血管収縮が進み、発汗量も減少する傾向にある。このような体温調節系の季節変化が温熱快適性に影響すると考えられ（例えば、非特許文献３〜６参照。）、前記季節調整パラメータは、その影響の度合いを表現するために導入されている。また、１年の前期と後期とでは体温調節系の季節変化が温熱快適性に及ぼす影響が異なる可能性があるので、前期と後期とではｓｉｎ関数の振幅は異なるものと仮定した（ｂc-seasonの前期及び後期におけるｓｉｎ関数の振幅はＣｏｌｄ0 及びＣｏｌｄ1 に相当し、ｂh-seasonの前期及び後期におけるｓｉｎ関数の振幅はＨｏｔ0 及びＨｏｔ1 に相当する。）。なお、季節調整パラメータにおけるｓｉｎ関数表現は、日本の平均外気温の変動傾向にも概ね対応する。
【非特許文献３】中山昭雄編、「温熱生理学」、理工学社、１９８１
【非特許文献４】深井他、「標準新有効温度（ＳＥＴ* ）と日本人の温熱感覚に関する実験的研究」、空気調和・衛生工学論文集、Ｎｏ．５１、ｐｐ．１３９−１４７、１９９２
【非特許文献５】久保他、「温冷感と快適感の季節差に関する実験的研究」、人間と生活環境、１（１）、５１−５７、１９９４
【非特許文献６】安岡他、「住空間における生理心理反応からみた温熱的快適範囲の季節差に関する研究」、日本建築学会環境系論文集、第７６巻、第６６３号、ｐｐ．４７９−４８４、２０１１
【００６９】
前記季節調整パラメータによる不満足率の季節変化は、例えば、図８のように表される。図８は、着衣及び姿勢を統一した優等車両実験のデータを使用し、全季節の実験結果に最もよく適合するものを選定して同定した結果であって、乗客が女性の場合を示している。なお、図において、（ａ）は季節調整パラメータの変化を示し、（ｂ）は不満足率（ＤｙＰＰＤ）の変化を示している。
【００７０】
図８に示される例は、１年の前期のうち、春から夏にかけては、寒い不快がより生じにくくなるが、暑い不快が生じやすくなることを示している。また、１年の後期のうち、秋から冬にかけては、寒い不快がより生じやすくなるが、秋から春にかけて、暑い不快は変化がない、すなわち、暑い不快に季節変化がないと近似することができることを示している。
【００７１】
なお、一部の鉄道車両では、空調制御において季節補正を行っているが、これは季節による乗客の着衣量の変化を推定して設定温度を補正するものであって、温熱快適性の季節特性を十分に考慮しているものではない（例えば、特許文献３参照。）。過去の実験研究において、同じ温熱環境（温度、湿度、風速、放射）、着衣量及び代謝量であっても、季節によって、人の感じる快適性が異なることが報告されている。したがって、年間を通して快適な車内空調の環境を実現するためには、本実施の形態において説明したような季節調整パラメータを使用し、温熱快適性の季節特性を考慮した車内空調を行う必要がある。
【特許文献３】特開２００１−１９９３３２号公報
【００７２】
そして、第１補正量決定部３４ｂの補正対象決定部３４ｃは、図９に示されるように、乗客男女比の情報に基づく男女比による重付け加算を行って、算出された寒い不満足率及び暑い不満足率に基づいて、平均温冷感最適補正を冷房強度の補正で行うか、又は、横流ファン強度の補正で行うか、すなわち、冷房強度又は横流ファン強度のいずれに対して補正量決定関数を適用するかを決定する。なお、図９に示される方法は、できる限り冷房強度を低く抑えて補正するための方法である。
【００７３】
乗客男女比の情報は、例えば、女性専用車両として使用される時間帯には女性が１で男性が０であるとする。また、例えば、線区又は時間帯毎の男女比に関する乗客情報を利用することができる場合には該乗客情報に基づいて設定される。さらに、例えば、全く未知である場合には、デフォルトとして男女が半々（０．５）であるものとする。
【００７４】
補正対象決定部３４ｃは、男女比による重付け加算を行った後、暑い不満足率（ＤｙＰＰＤ＿ｈ）が寒い不満足率（ＤｙＰＰＤ＿ｃ）以上であるか否かを判断する。
【００７５】
そして、暑い不満足率が寒い不満足率以上である場合、すなわち、暑い不満足率が大きい場合、補正対象決定部３４ｃは、横流ファン強度（Ｎ）がＮmax であるか否かを判断する。前述のように、横流ファン強度には、１〜Ｎmax 段階の設定があるものとする。そして、横流ファン強度がＮmax である場合、補正対象決定部３４ｃは、冷房強度に対して補正量決定関数を適用することに決定する。また、横流ファン強度がＮmax でない場合、補正対象決定部３４ｃは、できる限り横流ファン強度の補正で快適性を改善するように、横流ファン強度に対して補正量決定関数を適用することに決定する。
【００７６】
一方、暑い不満足率が寒い不満足率以上でない場合、すなわち、寒い不満足率が大きい場合、補正対象決定部３４ｃは、冷房強度が１であるか否かを判断する。そして、冷房強度が１である場合、補正対象決定部３４ｃは、横流ファン強度に対して補正量決定関数を適用することに決定する。また、冷房強度が１でない場合、補正対象決定部３４ｃは、できる限り冷房強度を下げて快適性を改善するように、冷房強度に対して補正量決定関数を適用することに決定する。
【００７７】
そして、第１補正量決定部３４ｂの補正量決定関数適用部３４ｄは、補正対象決定部３４ｃの決定に従って、冷房強度又は横流ファン強度に対して補正量決定関数を適用し、第１補正量を決定する。なお、前記補正量決定関数は、例えば、以下の式（１３）で表されるファジールールである。
【００７８】
【数８】
2018075913.tif 000010
【００７９】
なお、該式（１３）では、冷房強度及び横流ファン強度が数段階（例えば、５段階）に設定されていることを想定している。
【００８０】
例えば、前記式（１３）における最初の行（１）は、暑い不満足率が寒い不満足率よりも大きく、かつ、前時刻よりも暑い不満足率が増大している、という状況なので、冷房強度又は横流ファン強度をより大きくする、という補正を行う旨を意味している。行（１）の意味する状況は、図１０に示されている。
【００８１】
また、例えば、前記式（１３）における最後の行（９）は、寒い不満足率が暑い不満足率よりも大きく、かつ、前時刻よりも寒い不満足率が増大している、という状況なので、冷房強度又は横流ファン強度をより小さくする、という補正を行う旨を意味している。
【００８２】
このようにして、第１補正量を決定することによって、乗客の寒い不満足率及び暑い不満足率の合計の不満足率が最小となるような最適平均温冷感を実現することができる。なお、暖房を行う場合にも、同様の考え方によって、第１補正量を決定することができる。そして、決定された空調強度の補正量である第１補正量は、図２に示されるように、第２補正量計算部３６が算出した他の空調強度の補正量である第２補正量に加算され、補正量取得部３７に送信される。
【００８３】
本実施の形態において、外気温及び乗車人数取得部３５並びに第２補正量計算部３６が行う動作である今後の温冷感の変動を抑制するための補正の動作についての説明は、既知の先行技術（例えば、非特許文献７参照。）と同様であるので、その詳細な説明を省略する。なお、前記先行技術は、例えば、外気温と乗車率とを変数とするファジー推論による制御であり、冷房を想定した場合には、「ＩＦ  外気温が高い  ＡＮＤ  乗車率が高い  ＴＨＥＮ  冷房強度をより大きくする」というような制御である。
【非特許文献７】白石他、「車内環境向上を目指した空調システム」、三菱電機技報、１８（１０）、２７−３７、２００７
【００８４】
このように、本実施の形態において、周期的変動風を利用する車内空調方法は、周期的変動風の風速を人体熱モデルに入力して体感温度を算出する工程と、体感温度から平均温冷感予測を算出する工程と、平均温冷感予測から不満足率予測を算出する工程と、不満足率予測から空調強度又は周期的変動風発生器としての横流ファン２３の強度の補正量を決定する工程と、を含んでいる。これにより、周期的変動風全体の総合的な温冷感を適切に反映することができる。したがって、快適な室内空調環境を乗客に対して提供することができる。
【００８５】
また、人体熱モデルに入力される周期的変動風の風速は、周期的変動風体感時と予測平均皮膚温が同等となる定常風の風速に換算された定常風等価換算風速である。さらに、定常風等価換算風速は、周期的変動風体感時の予測平均皮膚温と定常風体感時の予測平均皮膚温との誤差が極小となる定常風の風速である。さらに、横流ファン２３の強度が所定の閾値より小さい場合、人体熱モデルに入力される周期的変動風の風速は、実測された風速である。さらに、不満足率予測から空調強度又は横流ファン２３の強度の補正量を決定する工程では、できる限り空調強度を低く抑えるように補正量を決定する。さらに、周期的変動風を利用する車内空調方法は、外気温及び乗客が発生する熱量から、温冷感変動を抑制する空調強度の他の補正量を計算する工程と、空調強度の他の補正量と、不満足率予測から決定された空調強度又は横流ファン２３の強度の補正量とを加算する工程とを更に含んでいる。
【００８６】
なお、本明細書の開示は、好適で例示的な実施の形態に関する特徴を述べたものである。ここに添付された特許請求の範囲内及びその趣旨内における種々の他の実施の形態、修正及び変形は、当業者であれば、本明細書の開示を総覧することにより、当然に考え付くことである。
【産業上の利用可能性】
【００８７】
本発明は、周期的変動風を利用する車内空調方法及びシステムに適用することができる。
【符号の説明】
【００８８】
１１    車両
１２    客室
２１    空調装置
２２    制御器
２３    横流ファン
３２ａ    定常風換算部
３３ａ    人体熱モデル部
３３ｂ    平均温冷感予測部
３４ａ    熱的不満足率予測部
【図１】
2018075913.tif 000011
【図２】
2018075913.tif 000012
【図３】
2018075913.tif 000013
【図４】
2018075913.tif 000014
【図５】
2018075913.tif 000015
【図６】
2018075913.tif 000016
【図７】
2018075913.tif 000017
【図８】
2018075913.tif 000018
【図９】
2018075913.tif 000019
【図１０】
2018075913.tif 000020