JPH09311617A - 地下空間における空気調和・換気シミュレータ - Google Patents
地下空間における空気調和・換気シミュレータInfo
- Publication number
- JPH09311617A JPH09311617A JP8127982A JP12798296A JPH09311617A JP H09311617 A JPH09311617 A JP H09311617A JP 8127982 A JP8127982 A JP 8127982A JP 12798296 A JP12798296 A JP 12798296A JP H09311617 A JPH09311617 A JP H09311617A
- Authority
- JP
- Japan
- Prior art keywords
- data
- temperature
- branch
- processing unit
- airflow
- 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
Links
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
- F24F11/00—Control or safety arrangements
- F24F11/30—Control or safety arrangements for purposes related to the operation of the system, e.g. for safety or monitoring
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
- F24F11/00—Control or safety arrangements
- F24F11/62—Control or safety arrangements characterised by the type of control or by internal processing, e.g. using fuzzy logic, adaptive control or estimation of values
- F24F11/63—Electronic processing
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D27/00—Simultaneous control of variables covered by two or more of main groups G05D1/00 - G05D25/00
- G05D27/02—Simultaneous control of variables covered by two or more of main groups G05D1/00 - G05D25/00 characterised by the use of electric means
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Ventilation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
(57)【要約】 (修正有)
【課題】 地下鉄道等の地下空間を構成する各ブランチ
内の気流速度、圧力、温度、湿度、CO2 濃度、粉じん
濃度を高い精度で予測する。 【解決手段】 トンネル、階段、通路、換気口等の各ブ
ランチの気流速度を求める第1の処理部1、各ブランチ
の壁温度、土中温度の短期変動量及び各ブランチ内の温
度、水蒸気密度、CO2 濃度を演算する第2の処理部
2、各ブランチの壁、土中の温度、気流温度の演算処理
を行う第3の処理部3並びに各ブランチの気流の温度、
水蒸気密度、CO2 濃度の時間変化データを演算する気
流速度、温湿度、CO2 濃度の詳細な演算処理を行う第
4の処理部4を備える。
内の気流速度、圧力、温度、湿度、CO2 濃度、粉じん
濃度を高い精度で予測する。 【解決手段】 トンネル、階段、通路、換気口等の各ブ
ランチの気流速度を求める第1の処理部1、各ブランチ
の壁温度、土中温度の短期変動量及び各ブランチ内の温
度、水蒸気密度、CO2 濃度を演算する第2の処理部
2、各ブランチの壁、土中の温度、気流温度の演算処理
を行う第3の処理部3並びに各ブランチの気流の温度、
水蒸気密度、CO2 濃度の時間変化データを演算する気
流速度、温湿度、CO2 濃度の詳細な演算処理を行う第
4の処理部4を備える。
Description
【0001】
【発明が属する技術分野】本発明は、地下空間における
空気調和・換気シミュレータに係り、特に、地下鉄道、
地下街等の地下空間における風、温度、湿度、紛塵濃
度、二酸化炭素(以下、CO2 という)濃度等の環境を
シミュレーションし、冷凍機を含む空気調和機の設備容
量等の設定のための情報を作成する空気調和・換気シミ
ュレータに関する。
空気調和・換気シミュレータに係り、特に、地下鉄道、
地下街等の地下空間における風、温度、湿度、紛塵濃
度、二酸化炭素(以下、CO2 という)濃度等の環境を
シミュレーションし、冷凍機を含む空気調和機の設備容
量等の設定のための情報を作成する空気調和・換気シミ
ュレータに関する。
【0002】
【従来の技術】地下空間における空気調和・換気シミュ
レータに関する従来技術として、例えば、U.S.Departme
nt of Transportation,FEBRURY発行のレポート“HANDBO
OK VALIDATION THROUGH FIELD TESTING AND SUBWAY ENV
IRONMENTAL DESIGN”に記載されたSUBWAY ENVIRONMENTA
L SIMULATION PROGRAMが知られている。
レータに関する従来技術として、例えば、U.S.Departme
nt of Transportation,FEBRURY発行のレポート“HANDBO
OK VALIDATION THROUGH FIELD TESTING AND SUBWAY ENV
IRONMENTAL DESIGN”に記載されたSUBWAY ENVIRONMENTA
L SIMULATION PROGRAMが知られている。
【0003】この従来技術は、地下鉄道におけるトンネ
ルを含む空間に対する空気調和・換気シミュレータに関
するもので、トンネル内の空気流を非圧縮性の1次元の
流れとして捕らえ、トンネル内の空気がいたるところで
一様な速度で流れるものとして扱われている。また、こ
の従来技術において、トンネル内に列車が走行する場合
に列車が受ける空気抵抗は、流体としての空気の粘性と
は無関係に、列車と気流との相対速度の2乗に比例する
ものとして、空気抵抗係数と列車の断面積とから決定さ
れており、また、空気抵抗係数は、トンネル内気流速
度、列車速度、トンネル壁の摩擦抵抗係数、及び、トン
ネルと列車との間の幾何学的条件によって決定されてい
る。
ルを含む空間に対する空気調和・換気シミュレータに関
するもので、トンネル内の空気流を非圧縮性の1次元の
流れとして捕らえ、トンネル内の空気がいたるところで
一様な速度で流れるものとして扱われている。また、こ
の従来技術において、トンネル内に列車が走行する場合
に列車が受ける空気抵抗は、流体としての空気の粘性と
は無関係に、列車と気流との相対速度の2乗に比例する
ものとして、空気抵抗係数と列車の断面積とから決定さ
れており、また、空気抵抗係数は、トンネル内気流速
度、列車速度、トンネル壁の摩擦抵抗係数、及び、トン
ネルと列車との間の幾何学的条件によって決定されてい
る。
【0004】しかし、トンネル断面積に対する列車断面
積の割合(以下、閉塞率という)が高く、この閉塞率が
0.5程度を越えるトンネル内では、気体の圧縮性の効
果が顕著になり、空気抵抗は、流体の粘性と無関係に列
車と気流との相対速度の2乗に比例するものではなくな
ってしまう。また、前述の閉塞率が低い値の場合にも、
トンネル内を走行する列車が、短い距離で、加速、高速
走行、減速を繰り返す走行パターンで運行される場合、
空気抵抗係数の時間依存性が顕著となり、空気抵抗係数
を一定値として取り扱うことができなくなってしまう。
積の割合(以下、閉塞率という)が高く、この閉塞率が
0.5程度を越えるトンネル内では、気体の圧縮性の効
果が顕著になり、空気抵抗は、流体の粘性と無関係に列
車と気流との相対速度の2乗に比例するものではなくな
ってしまう。また、前述の閉塞率が低い値の場合にも、
トンネル内を走行する列車が、短い距離で、加速、高速
走行、減速を繰り返す走行パターンで運行される場合、
空気抵抗係数の時間依存性が顕著となり、空気抵抗係数
を一定値として取り扱うことができなくなってしまう。
【0005】
【発明が解決しようとする課題】前述したように、従来
技術による空気調和・換気シミュレータは、トンネルと
列車との間の幾何学的条件を必要とするため、これらの
条件によるデータを求めるための膨大な実験を必要と
し、また、空気抵抗が、流体としての空気の粘性とは無
関係であるという実際とは異なる条件の下にシミュレー
ションを行っているため、地下空間における気流の状態
を正確にシミュレートすることができないという問題点
を有している。従って、この従来技術により得られたデ
ータを用いても、空気調和機の設備容量等を正確に設定
することができなかった。
技術による空気調和・換気シミュレータは、トンネルと
列車との間の幾何学的条件を必要とするため、これらの
条件によるデータを求めるための膨大な実験を必要と
し、また、空気抵抗が、流体としての空気の粘性とは無
関係であるという実際とは異なる条件の下にシミュレー
ションを行っているため、地下空間における気流の状態
を正確にシミュレートすることができないという問題点
を有している。従って、この従来技術により得られたデ
ータを用いても、空気調和機の設備容量等を正確に設定
することができなかった。
【0006】本発明の目的は、前記従来技術の問題点を
解決し、地下空間における気流の状態を正確にシミュレ
ートし、地下空間における温度、湿度、CO2 濃度、粉
じん濃度を正確に予測することができる地下空間におけ
る空気調和・換気シミュレータを提供することにある。
解決し、地下空間における気流の状態を正確にシミュレ
ートし、地下空間における温度、湿度、CO2 濃度、粉
じん濃度を正確に予測することができる地下空間におけ
る空気調和・換気シミュレータを提供することにある。
【0007】
【課題を解決するための手段】本発明によれば前記目的
は、地下空間における環境状態の予測を行う空気調和・
換気シミュレータにおいて、地下空間を構成するトンネ
ル、階段、通路、換気孔等のブランチによる地下空間構
造モデルデータ、前記地下空間内で空気を移動させる運
動物体の運動データを入力データとして、前記各ブラン
チの気流速度を求める第1の処理部と、前記空間構造モ
デルデータ、前記第1の処理部からの解析結果データ、
運動物体の熱特性に関するデータ、外気温湿度、外気C
O2 濃度データ、前記空間へ出入りする客データを入力
として、各ブランチの壁温度、土中温度の初期値と短期
変動量及び各ブランチ内の温度、水蒸気密度、CO2 濃
度の初期値と境界値を演算する第2の処理部と、前記空
間構造モデルデータ、第1の処理部からの気流解析結果
データ、第2の処理部からの温熱解析結果データを入力
として、各ブランチの壁、土中の温度、各ブランチの気
流温度の演算処理を行う第3の処理部と、前記空間構造
モデルデータ、第1の処理部からの気流解析結果デー
タ、第2の処理部からの温熱解析結果データ、第3の処
理部からの壁温度、土中温度データを入力として、各ブ
ランチの気流の温度、水蒸気密度、CO2 濃度の時間変
化データを演算する気流速度、温湿度、CO2 濃度の詳
細な演算処理を行う第4の処理部とを備えることにより
達成される。
は、地下空間における環境状態の予測を行う空気調和・
換気シミュレータにおいて、地下空間を構成するトンネ
ル、階段、通路、換気孔等のブランチによる地下空間構
造モデルデータ、前記地下空間内で空気を移動させる運
動物体の運動データを入力データとして、前記各ブラン
チの気流速度を求める第1の処理部と、前記空間構造モ
デルデータ、前記第1の処理部からの解析結果データ、
運動物体の熱特性に関するデータ、外気温湿度、外気C
O2 濃度データ、前記空間へ出入りする客データを入力
として、各ブランチの壁温度、土中温度の初期値と短期
変動量及び各ブランチ内の温度、水蒸気密度、CO2 濃
度の初期値と境界値を演算する第2の処理部と、前記空
間構造モデルデータ、第1の処理部からの気流解析結果
データ、第2の処理部からの温熱解析結果データを入力
として、各ブランチの壁、土中の温度、各ブランチの気
流温度の演算処理を行う第3の処理部と、前記空間構造
モデルデータ、第1の処理部からの気流解析結果デー
タ、第2の処理部からの温熱解析結果データ、第3の処
理部からの壁温度、土中温度データを入力として、各ブ
ランチの気流の温度、水蒸気密度、CO2 濃度の時間変
化データを演算する気流速度、温湿度、CO2 濃度の詳
細な演算処理を行う第4の処理部とを備えることにより
達成される。
【0008】
【発明の実施の形態】以下、本発明による空気調和・換
気シミュレータの一実施形態を図面により詳細に説明す
る。
気シミュレータの一実施形態を図面により詳細に説明す
る。
【0009】図1は本発明の一実施形態による空気調和
・換気シミュレータの機能構成を示すブロック図、図2
は第1の処理部である気流・圧力解析処理部の処理動作
を説明するフローチャート、図3は第2の処理部である
温熱解析処理部の処理動作を説明するフローチャート、
図4は第3の処理部である土中温度・壁温度演算処理部
の処理動作を説明するフローチャート、図5は第4の処
理部である気流速度・温湿度・CO2 濃度詳細演算処理
部の処理動作を説明するフローチャート(その1)、図
6は第4の処理部である気流速度・温湿度・CO2 濃度
詳細演算処理部の処理動作を説明するフローチャート
(その2)、図7は本発明の実施形態でシミュレーショ
ンの対象とする地下空間の構造例を示す図である。図
1、図7において、1は気流・圧力解析処理部、2は温
熱解析処理部、3は土中温度・壁温度演算処理部、4は
気流速度・温湿度・CO2 濃度詳細演算処理部、11は
排気ファン、12は給気ファン、13はホーム、14は
コンコース、15は空気調和機、16は冷凍機、17は
制御装置、18a、18bは温度センサ、19はトンネ
ル、20は排気口、21は給気口、22は階段である。
・換気シミュレータの機能構成を示すブロック図、図2
は第1の処理部である気流・圧力解析処理部の処理動作
を説明するフローチャート、図3は第2の処理部である
温熱解析処理部の処理動作を説明するフローチャート、
図4は第3の処理部である土中温度・壁温度演算処理部
の処理動作を説明するフローチャート、図5は第4の処
理部である気流速度・温湿度・CO2 濃度詳細演算処理
部の処理動作を説明するフローチャート(その1)、図
6は第4の処理部である気流速度・温湿度・CO2 濃度
詳細演算処理部の処理動作を説明するフローチャート
(その2)、図7は本発明の実施形態でシミュレーショ
ンの対象とする地下空間の構造例を示す図である。図
1、図7において、1は気流・圧力解析処理部、2は温
熱解析処理部、3は土中温度・壁温度演算処理部、4は
気流速度・温湿度・CO2 濃度詳細演算処理部、11は
排気ファン、12は給気ファン、13はホーム、14は
コンコース、15は空気調和機、16は冷凍機、17は
制御装置、18a、18bは温度センサ、19はトンネ
ル、20は排気口、21は給気口、22は階段である。
【0010】以下に説明する本発明の一実施形態は、本
発明を地下鉄道システムに適用したものとして説明す
る。
発明を地下鉄道システムに適用したものとして説明す
る。
【0011】本発明の一実施形態による空気調和・換気
シミュレータは、図1に示すように、気流・圧力解析処
理部(以下、第1の処理部という)1と、温熱解析処理
部(以下、第2の処理部という)2と、土中温度・壁温
度演算処理部(以下、第3の処理部という)3と、気流
速度・温湿度・CO2 濃度詳細演算処理部(以下、第4
の処理部という)4とを備えて構成される。
シミュレータは、図1に示すように、気流・圧力解析処
理部(以下、第1の処理部という)1と、温熱解析処理
部(以下、第2の処理部という)2と、土中温度・壁温
度演算処理部(以下、第3の処理部という)3と、気流
速度・温湿度・CO2 濃度詳細演算処理部(以下、第4
の処理部という)4とを備えて構成される。
【0012】本発明において、前述の各処理部は、公知
のワークステーション等の電子計算機の中に構築され、
処理に必要なデータ等は、外部メモリ等に構築されるフ
ァイルから入力され、処理結果のデータは、必要に応じ
て他の他の処理部に渡され、また、ディスプレイに表示
され、必要によりプリントアウトされるものとする。従
って、ここではワークステーション等の電子計算機のハ
ードウェア構成についての説明を省略する。
のワークステーション等の電子計算機の中に構築され、
処理に必要なデータ等は、外部メモリ等に構築されるフ
ァイルから入力され、処理結果のデータは、必要に応じ
て他の他の処理部に渡され、また、ディスプレイに表示
され、必要によりプリントアウトされるものとする。従
って、ここではワークステーション等の電子計算機のハ
ードウェア構成についての説明を省略する。
【0013】まず、図7を参照して、本発明の実施形態
でシミュレーションの対象とする地下空間の構造例であ
る地下鉄の駅構造の例について説明する。
でシミュレーションの対象とする地下空間の構造例であ
る地下鉄の駅構造の例について説明する。
【0014】図7は地下鉄の駅構造を極めて簡単に示す
もので、駅構造は、隣接する駅に接続され列車が運行す
るトンネル19と、駅ホーム13と、地上に通じる階段
22及びコンコース14と、外部との間で空気の交換を
行うための排気ファン11を有する排気口20及び給気
ファン12を有する給気口21によりなり、駅構内の空
気調和のための冷凍機16、温度センサ18a、18b
からの情報を受ける制御装置17により運転される空気
調和装置15を備えている。
もので、駅構造は、隣接する駅に接続され列車が運行す
るトンネル19と、駅ホーム13と、地上に通じる階段
22及びコンコース14と、外部との間で空気の交換を
行うための排気ファン11を有する排気口20及び給気
ファン12を有する給気口21によりなり、駅構内の空
気調和のための冷凍機16、温度センサ18a、18b
からの情報を受ける制御装置17により運転される空気
調和装置15を備えている。
【0015】本発明の実施形態は、前述した例のように
構成される、駅構造のトンネル19、階段22、コンコ
ース14、排気口11、給気口21等の気流の速度、温
度、湿度、CO2 濃度を予測し、また、これらの構造物
の壁面温度、土中温度等を予測しようとするものであ
る。
構成される、駅構造のトンネル19、階段22、コンコ
ース14、排気口11、給気口21等の気流の速度、温
度、湿度、CO2 濃度を予測し、また、これらの構造物
の壁面温度、土中温度等を予測しようとするものであ
る。
【0016】次に、図1を参照して、各処理部の機能を
簡単に説明する。
簡単に説明する。
【0017】気流、圧力の解析を行う第1の処理部1
は、後に詳述するが、駅構造モデルデータ、列車ダイヤ
データ、列車特性データ等を入力データとして、列車が
存在するブランチ内の列車側方領域毎の列車圧力差を演
算し、さらに、トンネルブランチ、ホームブランチ、一
般ブランチ、給気ファンブランチ、排気ファンブラン
チ、及び、空気調和機ブランチの気流速度を算出する。
これらの各ブランチの気流速度のデータは、時間変化デ
ータとして出力され、10秒毎のデータが24時間分出
力され、第2〜第4の処理部2〜4に与えられる。
は、後に詳述するが、駅構造モデルデータ、列車ダイヤ
データ、列車特性データ等を入力データとして、列車が
存在するブランチ内の列車側方領域毎の列車圧力差を演
算し、さらに、トンネルブランチ、ホームブランチ、一
般ブランチ、給気ファンブランチ、排気ファンブラン
チ、及び、空気調和機ブランチの気流速度を算出する。
これらの各ブランチの気流速度のデータは、時間変化デ
ータとして出力され、10秒毎のデータが24時間分出
力され、第2〜第4の処理部2〜4に与えられる。
【0018】前述した列車圧力差による気流生成量の演
算は、従来技術の場合のように列車の空気抵抗係数を用
いることなく行われるものであり、例えば、「ながれ」
第12巻第4号(平成5年12月30日、日本流体学会
発行)等に記載された緒理論を用いて行うことができ
る。
算は、従来技術の場合のように列車の空気抵抗係数を用
いることなく行われるものであり、例えば、「ながれ」
第12巻第4号(平成5年12月30日、日本流体学会
発行)等に記載された緒理論を用いて行うことができ
る。
【0019】温熱解析を行う第2の処理部は、駅構造モ
デルデータ、外気温データ、外気湿度データ、外気CO
2 濃度データ、駅乗降客データ、列車の熱特性に関する
データ、第1の処理部1による気流解析データ等を入力
として、各ブランチ、駅、トンネル関連の温度と水蒸気
密度とCO2 濃度を演算し、第3の処理部3で用いるト
ンネルの壁と土中温度の初期値と短期変動量、第4の処
理部4で用いる列車、ブランチ及び駅、トンネルのそれ
ぞれの温度と水蒸気密度とCO2 濃度との初期値及び境
界値を演算する。これらの各データは、所定の時間間隔
毎に365日分が出力される。
デルデータ、外気温データ、外気湿度データ、外気CO
2 濃度データ、駅乗降客データ、列車の熱特性に関する
データ、第1の処理部1による気流解析データ等を入力
として、各ブランチ、駅、トンネル関連の温度と水蒸気
密度とCO2 濃度を演算し、第3の処理部3で用いるト
ンネルの壁と土中温度の初期値と短期変動量、第4の処
理部4で用いる列車、ブランチ及び駅、トンネルのそれ
ぞれの温度と水蒸気密度とCO2 濃度との初期値及び境
界値を演算する。これらの各データは、所定の時間間隔
毎に365日分が出力される。
【0020】土中温度、壁温度の演算処理を行う第3の
処理部3は、駅構造モデルデータ、第1の処理部1から
の気流解析結果データ、第2の処理部2からの温熱解析
結果データを入力として、第4の処理部4で用いる駅及
びトンネルの壁の温度と土中の温度、各ブランチの気流
温度の6時間毎の時間変化の20年間にわたるデータを
演算する。
処理部3は、駅構造モデルデータ、第1の処理部1から
の気流解析結果データ、第2の処理部2からの温熱解析
結果データを入力として、第4の処理部4で用いる駅及
びトンネルの壁の温度と土中の温度、各ブランチの気流
温度の6時間毎の時間変化の20年間にわたるデータを
演算する。
【0021】気流速度・温湿度・CO2 濃度の詳細な演
算処理を行う第4の処理部4は、駅構造モデルデータ、
第1の処理部1からの気流解析結果データ、第2の処理
部2からの温熱解析結果データ、第3の処理部からの土
中温度、壁温度の演算結果データを入力として、各ブラ
ンチの気流、列車、ブランチ、駅、トンネルのそれぞれ
の温度と水蒸気密度とCO2 濃度との時間変化データ、
統計解析(平均値、最大値、最小値)データを演算して
出力する。
算処理を行う第4の処理部4は、駅構造モデルデータ、
第1の処理部1からの気流解析結果データ、第2の処理
部2からの温熱解析結果データ、第3の処理部からの土
中温度、壁温度の演算結果データを入力として、各ブラ
ンチの気流、列車、ブランチ、駅、トンネルのそれぞれ
の温度と水蒸気密度とCO2 濃度との時間変化データ、
統計解析(平均値、最大値、最小値)データを演算して
出力する。
【0022】次に、前述で説明した本発明の一実施形態
を構成する各処理部の動作をフローを参照して説明す
る。
を構成する各処理部の動作をフローを参照して説明す
る。
【0023】第1の処理部1は、詳細には、実行条件デ
ータ、ブランチ詳細構造データ、構造構築データ、ブラ
ンチ依存の管理データ、続行用の初期データ、及び、定
数データを入力データとする。
ータ、ブランチ詳細構造データ、構造構築データ、ブラ
ンチ依存の管理データ、続行用の初期データ、及び、定
数データを入力データとする。
【0024】実行条件データは、シミュレーションの開
始時刻、終了時刻、演算に用いる内部ファイル(構造モ
デルデータ、ダイヤデータ)名称、続行用初期値データ
の使用の有無、処理結果のプリントアウトの有無を与え
るデータである。
始時刻、終了時刻、演算に用いる内部ファイル(構造モ
デルデータ、ダイヤデータ)名称、続行用初期値データ
の使用の有無、処理結果のプリントアウトの有無を与え
るデータである。
【0025】ブランチ詳細構造データは各ブランチ(ト
ンネル、ホーム、階段、通路、出入口、緩和縦坑、ダク
ト、給気ファン、排気ファン、空気調和機)の物理特性
データ、壁面及びブランチ端の圧力損失係数、各種発熱
量、乗降客の滞在時間、及び、各種の係数を与えるデー
タである。
ンネル、ホーム、階段、通路、出入口、緩和縦坑、ダク
ト、給気ファン、排気ファン、空気調和機)の物理特性
データ、壁面及びブランチ端の圧力損失係数、各種発熱
量、乗降客の滞在時間、及び、各種の係数を与えるデー
タである。
【0026】構造構築データは、各ブランチを接続する
内部ポイントに接続されるブランチの総数、各内部ポイ
ントに接続されるブランチの総称種別と総称種別内イン
デックス、及び、接続される相手側ポイントの総称種別
と総称種別内インデックスを与えるデータであり、図示
しない地下鉄系構造構築処理部から与えられる。
内部ポイントに接続されるブランチの総数、各内部ポイ
ントに接続されるブランチの総称種別と総称種別内イン
デックス、及び、接続される相手側ポイントの総称種別
と総称種別内インデックスを与えるデータであり、図示
しない地下鉄系構造構築処理部から与えられる。
【0027】ブランチ依存の管理データは、列車が存在
するブランチ内の演算すべき列車数、列車前後端の位
置、列車の速度、列車の加速度、列車の長さを与えるデ
ータであり、図示しない列車運行管理処理部から与えら
れる。
するブランチ内の演算すべき列車数、列車前後端の位
置、列車の速度、列車の加速度、列車の長さを与えるデ
ータであり、図示しない列車運行管理処理部から与えら
れる。
【0028】続行用の初期データは、第1の処理部にお
けるシミュレーション開始時の各ブランチの気流流量
値、前述した「ながれ」等に記載された理論による列車
近傍の気流速度である列車ゆらぎ比速度を与えるデータ
である。
けるシミュレーション開始時の各ブランチの気流流量
値、前述した「ながれ」等に記載された理論による列車
近傍の気流速度である列車ゆらぎ比速度を与えるデータ
である。
【0029】第1の処理部1は、前述したデータを受
け、各ブランチの気流速度と、各ポイントの圧力とを演
算し、解析結果データ、プリントアウトデータ、及び、
他の処理部のための続行用初期値データを出力する。解
析結果データは、第2の処理部2から第4の処理部4で
用いる各ブランチの気流速度の時間変化データとして出
力され、10秒毎のデータが24時間分出力される。
け、各ブランチの気流速度と、各ポイントの圧力とを演
算し、解析結果データ、プリントアウトデータ、及び、
他の処理部のための続行用初期値データを出力する。解
析結果データは、第2の処理部2から第4の処理部4で
用いる各ブランチの気流速度の時間変化データとして出
力され、10秒毎のデータが24時間分出力される。
【0030】本発明とは直接的には関係しないが、前述
した構造構築データを与える地下鉄系構造構築処理部
と、ブランチ依存の管理データを与える列車運行管理処
理部について、以下に、簡単に説明する。
した構造構築データを与える地下鉄系構造構築処理部
と、ブランチ依存の管理データを与える列車運行管理処
理部について、以下に、簡単に説明する。
【0031】地下鉄系構造構築処理部は、地下鉄道の全
体を構成する駅総数、全ポイント総数、各種別毎のブラ
ンチ総数と全ブランチ総数、系内全体の各駅の構造をも
含む全体構造一覧データと全体構造ブランチ一覧デー
タ、及び、列車が運行する路線構成データを入力データ
とし、その出力データの一部として、前述したような本
発明に必要な構造構築データを出力するものである。
体を構成する駅総数、全ポイント総数、各種別毎のブラ
ンチ総数と全ブランチ総数、系内全体の各駅の構造をも
含む全体構造一覧データと全体構造ブランチ一覧デー
タ、及び、列車が運行する路線構成データを入力データ
とし、その出力データの一部として、前述したような本
発明に必要な構造構築データを出力するものである。
【0032】また、列車運行管理処理部は、運用列車総
数、車両運用データとダイヤパターンデータとを含むダ
イヤデータ、列車の各駅間での走行特性、列車及び軌道
プレートの発熱量を与えるランカーブデータ、線路構築
データを入力データとし、その出力データの一部とし
て、前述したような本発明に必要なブランチ依存の管理
データを出力するものである。
数、車両運用データとダイヤパターンデータとを含むダ
イヤデータ、列車の各駅間での走行特性、列車及び軌道
プレートの発熱量を与えるランカーブデータ、線路構築
データを入力データとし、その出力データの一部とし
て、前述したような本発明に必要なブランチ依存の管理
データを出力するものである。
【0033】次に、図2に示すフローを参照して、第1
の処理部の動作を説明する。
の処理部の動作を説明する。
【0034】(1)まず、図示しない実行条件データフ
ァイルより、この処理部での処理開始時刻と終了時刻、
構造モデルデータ内部ファイルとダイヤデータ内部ファ
イルの名称、続行用初期値データ使用の有無、これを使
用する場合の続行初期値データファイル名称を読み込む
(ステップ201)。
ァイルより、この処理部での処理開始時刻と終了時刻、
構造モデルデータ内部ファイルとダイヤデータ内部ファ
イルの名称、続行用初期値データ使用の有無、これを使
用する場合の続行初期値データファイル名称を読み込む
(ステップ201)。
【0035】(2)次に、構造モデルデータ内部ファイ
ルを読み込み、各データを内部変数にとして格納し、ダ
イヤデータ内部ファイルを読み込み、各データを内部変
数として格納する(ステップ202、203)。
ルを読み込み、各データを内部変数にとして格納し、ダ
イヤデータ内部ファイルを読み込み、各データを内部変
数として格納する(ステップ202、203)。
【0036】(3)初期設定として、定数データファイ
ルを読み込み、各定数データを内部変数として格納し、
前述で説明した図示しない地下鉄系構造構築処理部から
与えられる構造構築データを格納する。また、地下鉄系
構造構築処理部から与えられる圧力計算のための連立一
次方程式解くために用いる気流流量係数マトリクスを設
定する(ステップ204〜206)。
ルを読み込み、各定数データを内部変数として格納し、
前述で説明した図示しない地下鉄系構造構築処理部から
与えられる構造構築データを格納する。また、地下鉄系
構造構築処理部から与えられる圧力計算のための連立一
次方程式解くために用いる気流流量係数マトリクスを設
定する(ステップ204〜206)。
【0037】(4)各ブランチの気流速度に標準値を設
定し、続行用積分計算初期値データとして各ブランチの
気流速度値を設定すると共に、列車の初期状態を設定す
る(ステップ207)。
定し、続行用積分計算初期値データとして各ブランチの
気流速度値を設定すると共に、列車の初期状態を設定す
る(ステップ207)。
【0038】(5)前述までの処理でデータファイルの
入力処理と初期設定の処理とを終了し、次に、各時刻毎
(1秒毎)の列車登録、列車情報設定、ブランチ情報の
設定を行う。これらのデータは、図示しない列車運行管
理処理部から得られるデータであり、各1秒毎に、列車
登録処理として、運用列車が系内に出現中であるか、出
現前あるいは消滅後であるかを登録し、列車情報とし
て、存在列車について、運行中の線路指標、基準位置か
らの距離、列車の速度と加速度、及び、列車と軌道プレ
ート(リニアモータのための2次導体)との発熱量を設
定し、さらに、ブランチ情報設定として、列車の存在す
るブランチ番号を設定する(ステップ209〜21
2)。
入力処理と初期設定の処理とを終了し、次に、各時刻毎
(1秒毎)の列車登録、列車情報設定、ブランチ情報の
設定を行う。これらのデータは、図示しない列車運行管
理処理部から得られるデータであり、各1秒毎に、列車
登録処理として、運用列車が系内に出現中であるか、出
現前あるいは消滅後であるかを登録し、列車情報とし
て、存在列車について、運行中の線路指標、基準位置か
らの距離、列車の速度と加速度、及び、列車と軌道プレ
ート(リニアモータのための2次導体)との発熱量を設
定し、さらに、ブランチ情報設定として、列車の存在す
るブランチ番号を設定する(ステップ209〜21
2)。
【0039】(6)次に、時間依存データの設定とし
て、給気ファン、排気ファン、空気調和機等のON/O
FF時刻データを設定する(ステップ213)。
て、給気ファン、排気ファン、空気調和機等のON/O
FF時刻データを設定する(ステップ213)。
【0040】(7)前述までに設定された各データを使
用し、ルンゲークッタ法による積分計算のコントロール
を行い、各段の微係数演算と気流積分演算を行うことに
より、列車が存在するブランチ内の列車側方領域毎の列
車ゆらぎ比速度係数を演算すると共に、列車が存在する
ブランチ内の列車側方領域毎の列車圧力差を演算し、さ
らに、トンネルブランチ、ホームブランチ、一般ブラン
チ、給気ファンブランチ、排気ファンブランチ、及び、
空気調和機ブランチの気流流量微係数を演算する(ステ
ップ214〜217)。
用し、ルンゲークッタ法による積分計算のコントロール
を行い、各段の微係数演算と気流積分演算を行うことに
より、列車が存在するブランチ内の列車側方領域毎の列
車ゆらぎ比速度係数を演算すると共に、列車が存在する
ブランチ内の列車側方領域毎の列車圧力差を演算し、さ
らに、トンネルブランチ、ホームブランチ、一般ブラン
チ、給気ファンブランチ、排気ファンブランチ、及び、
空気調和機ブランチの気流流量微係数を演算する(ステ
ップ214〜217)。
【0041】(8)次に、気流流量係数を計算し、連立
一次方程式によるアクティブカラム法により各ポイント
の圧力を演算すると共に、この結果を基に各ブランチの
気流流量微係数を演算する(ステップ218、21
9)。
一次方程式によるアクティブカラム法により各ポイント
の圧力を演算すると共に、この結果を基に各ブランチの
気流流量微係数を演算する(ステップ218、21
9)。
【0042】(9)単純統計解析処理を行い、平均値計
算により各ブランチの1時間毎のプラス側及びマイナス
側の平均値と出現時刻回数を計算し、最小値最大計算に
より各ブランチの1時間毎の最小値及び最大値とそれぞ
れの発生時刻、出現時刻回数を計算する(ステップ22
0)。
算により各ブランチの1時間毎のプラス側及びマイナス
側の平均値と出現時刻回数を計算し、最小値最大計算に
より各ブランチの1時間毎の最小値及び最大値とそれぞ
れの発生時刻、出現時刻回数を計算する(ステップ22
0)。
【0043】(10)最後に、アウトプット用データファ
イルとして、続行用初期値データファイル、解析結果デ
ータファイルを出力する(ステップ221)。
イルとして、続行用初期値データファイル、解析結果デ
ータファイルを出力する(ステップ221)。
【0044】前述したステップ214〜217における
気流速度と圧力の演算は、列車の運動によるトンネル内
空気の運動と換気装置による換気能力とからトンネル、
駅の各区間の接続点の圧力を求め、接続点間の圧力差と
接続点間の圧力損失とから接続点間すなわちブランチの
風速を求めるもので、以下に式(1)として示すバラン
ス式を解いて求めることができる。
気流速度と圧力の演算は、列車の運動によるトンネル内
空気の運動と換気装置による換気能力とからトンネル、
駅の各区間の接続点の圧力を求め、接続点間の圧力差と
接続点間の圧力損失とから接続点間すなわちブランチの
風速を求めるもので、以下に式(1)として示すバラン
ス式を解いて求めることができる。
【0045】 気流流量変化量(dQt/dt) =入口端圧力による気流引き込み量((At/ρ(Lt−lg))Pin) +出口端圧力による気流引き込み量((At/ρ(Lt−lg))Pex) +列車圧力差による気流生成量((At/ρ(Lt−lg))δP) +区間壁面流体抵抗による気流流量減少量(−λ/2AtDt)|Dt|Qt +区間端分岐合流抵抗による気流流量減少量((-ζ/2At(Lt-lg))|Dt|Qt ……(1) 式(1)において、Qt:気流流量(m3/s)、t:時
間(s)、At:流路段面積(m3)、ρ:気流密度(K
g/m3)、Lt:区間長(m)、lg:列車長(m)、
Pin:入口端圧力(Pa)、Pex:出口端圧力(p
a)、δP:列車圧力差(Pa)、λ:壁面抵抗係数
(−)、Dt:トンネル等価直径(m)、ζ:入口端抵
抗係数(−)である。
間(s)、At:流路段面積(m3)、ρ:気流密度(K
g/m3)、Lt:区間長(m)、lg:列車長(m)、
Pin:入口端圧力(Pa)、Pex:出口端圧力(p
a)、δP:列車圧力差(Pa)、λ:壁面抵抗係数
(−)、Dt:トンネル等価直径(m)、ζ:入口端抵
抗係数(−)である。
【0046】第2の処理部2は、詳細には、実行条件デ
ータ、処理部1からの気流解析結果データ、ブランチ詳
細構造データ、構造構築データ、列車依存の管理デー
タ、ブランチ依存の管理データ、駅乗降客データ、駅間
乗客データ、続行用の初期データ、及び、定数データを
入力データとする。
ータ、処理部1からの気流解析結果データ、ブランチ詳
細構造データ、構造構築データ、列車依存の管理デー
タ、ブランチ依存の管理データ、駅乗降客データ、駅間
乗客データ、続行用の初期データ、及び、定数データを
入力データとする。
【0047】実行条件データは、処理部1からの気流解
析結果データファイル名称、内部データ(駅乗降客デー
タ、駅間乗客データ)ファイル名称、続行用初期値デー
タの使用の有無、処理結果のプリントアウトの有無を与
えるデータである。
析結果データファイル名称、内部データ(駅乗降客デー
タ、駅間乗客データ)ファイル名称、続行用初期値デー
タの使用の有無、処理結果のプリントアウトの有無を与
えるデータである。
【0048】続行用の初期値データは、第2の処理部2
の処理開始時の各ブランチに関連する温度と水蒸気密度
とCO2 濃度値、列車に関連する温度と水蒸気密度、及
び、駅、トンネルに関連する温度値を与えるデータであ
る。
の処理開始時の各ブランチに関連する温度と水蒸気密度
とCO2 濃度値、列車に関連する温度と水蒸気密度、及
び、駅、トンネルに関連する温度値を与えるデータであ
る。
【0049】第2の処理部2は、前述したデータを受
け、各ブランチ関連、駅、トンネル関連の温度と水蒸気
密度とCO2 濃度とを演算し、解析結果データ、プリン
トアウトデータ、及び、他の処理部のための続行用初期
値データを出力する。解析結果データは、第3の処理部
3で用いる駅及びトンネルの壁温度と土中温度との初期
値、短期変動量を与えるデータ、及び、第4の処理部4
で用いる列車、ブランチ、駅、及びトンネルのそれぞれ
の温度と水蒸気密度とCO2 濃度との初期値と境界値と
を与えるデータであり、365日分が出力される。
け、各ブランチ関連、駅、トンネル関連の温度と水蒸気
密度とCO2 濃度とを演算し、解析結果データ、プリン
トアウトデータ、及び、他の処理部のための続行用初期
値データを出力する。解析結果データは、第3の処理部
3で用いる駅及びトンネルの壁温度と土中温度との初期
値、短期変動量を与えるデータ、及び、第4の処理部4
で用いる列車、ブランチ、駅、及びトンネルのそれぞれ
の温度と水蒸気密度とCO2 濃度との初期値と境界値と
を与えるデータであり、365日分が出力される。
【0050】次に、図3に示すフローを参照して、第2
の処理部の動作を説明する。
の処理部の動作を説明する。
【0051】(1)まず、図示しない実行条件データフ
ァイルより、第1の処理部1による気流解析結果データ
ファイル名称、駅乗降客データ内部ファイル名称、駅間
乗客データ内部ファイル名称、続行用初期値データ使用
の有無、結果プリントアウトの有無を読み込む(ステッ
プ301)。
ァイルより、第1の処理部1による気流解析結果データ
ファイル名称、駅乗降客データ内部ファイル名称、駅間
乗客データ内部ファイル名称、続行用初期値データ使用
の有無、結果プリントアウトの有無を読み込む(ステッ
プ301)。
【0052】(2)第1の処理部1の気流解析結果デー
タ内部ファイルを読み込み、各データを内部変数として
格納し、構造モデルデータ内部ファイルを読み込み、各
データ内部変数として格納し、さらに、ダイヤデータ内
部ファイルを読み込み、各データを内部変数として格納
する(ステップ302〜304)。
タ内部ファイルを読み込み、各データを内部変数として
格納し、構造モデルデータ内部ファイルを読み込み、各
データ内部変数として格納し、さらに、ダイヤデータ内
部ファイルを読み込み、各データを内部変数として格納
する(ステップ302〜304)。
【0053】(3)駅乗降客データ内部ファイル、駅間
乗客データ内部ファイルを読み込み、各データを内部変
数として格納し、定数データファイル読み込み、各デー
タを内部変数として設定格納する(ステップ305〜3
07)。
乗客データ内部ファイルを読み込み、各データを内部変
数として格納し、定数データファイル読み込み、各デー
タを内部変数として設定格納する(ステップ305〜3
07)。
【0054】(4)前述で説明した図示しない地下鉄系
構造構築処理部から与えられる構造構築データを格納
し、積分計算初期値設定のために、各ブランチに関連す
る温度と水蒸気密度とCO2濃度 、列車に関連する温度
と水蒸気密度、及び、駅、トンネルに関連する温度のそ
れぞれに標準値を設定し、さらに、続行用初期値データ
ファイルを読み込み各ブランチに関連する温度と水蒸気
密度とCO2 濃度値、列車に関連する温度と水蒸気密度
値、及び、駅、トンネルに関連する温度値を設定する
(ステップ308〜310)。
構造構築処理部から与えられる構造構築データを格納
し、積分計算初期値設定のために、各ブランチに関連す
る温度と水蒸気密度とCO2濃度 、列車に関連する温度
と水蒸気密度、及び、駅、トンネルに関連する温度のそ
れぞれに標準値を設定し、さらに、続行用初期値データ
ファイルを読み込み各ブランチに関連する温度と水蒸気
密度とCO2 濃度値、列車に関連する温度と水蒸気密度
値、及び、駅、トンネルに関連する温度値を設定する
(ステップ308〜310)。
【0055】(5)前述までの処理でデータファイルの
入力処理と初期設定の処理とを終了し、次に、各時刻毎
(1秒毎)の列車登録、列車実施形態設定、ブランチ情
報の設定を行う。これらのデータは、図示しない列車運
行管理処理部から得られるデータであり、各10秒毎
に、列車登録処理として、運用列車が系内に出現中であ
るか、出現前あるいは消滅後であるかを登録し、列車情
報として、存在列車について、運行中の線路指標、基準
位置からの距離、列車の速度と加速度、及び、列車と軌
道プレート(リニアモータのための2次導体)との発熱
量を設定し、さらに、ブランチ情報設定として、列車の
存在するブランチ番号を設定する(ステップ312〜3
14)。
入力処理と初期設定の処理とを終了し、次に、各時刻毎
(1秒毎)の列車登録、列車実施形態設定、ブランチ情
報の設定を行う。これらのデータは、図示しない列車運
行管理処理部から得られるデータであり、各10秒毎
に、列車登録処理として、運用列車が系内に出現中であ
るか、出現前あるいは消滅後であるかを登録し、列車情
報として、存在列車について、運行中の線路指標、基準
位置からの距離、列車の速度と加速度、及び、列車と軌
道プレート(リニアモータのための2次導体)との発熱
量を設定し、さらに、ブランチ情報設定として、列車の
存在するブランチ番号を設定する(ステップ312〜3
14)。
【0056】(6)次に、時間依存データ設定として、
各ブランチの気流速度、外気の温度と水蒸気密度、駅に
存在する乗降客数と駅間を走行する列車に存在する乗客
数、各ブランチ壁面とチャンバー壁面、列車表面、及
び、軌道プレート表面の熱伝達率を設定する(ステップ
315〜318)。
各ブランチの気流速度、外気の温度と水蒸気密度、駅に
存在する乗降客数と駅間を走行する列車に存在する乗客
数、各ブランチ壁面とチャンバー壁面、列車表面、及
び、軌道プレート表面の熱伝達率を設定する(ステップ
315〜318)。
【0057】(7)次に、給気ファン、排気ファン、空
気調和機等に対する1日の営業時間に依存したON/O
FFデータと1年の季節に依存したON/OFFデータ
とを設定する(ステップ319)。
気調和機等に対する1日の営業時間に依存したON/O
FFデータと1年の季節に依存したON/OFFデータ
とを設定する(ステップ319)。
【0058】(8)前述までに設定された各データを使
用し、ルンゲークッタ法による積分計算のコントロール
を行い、各段の微係数演算と温熱積分演算を行うことに
より、系内に出現中の列車内気流温度と水蒸気密度、列
車温度の微係数を求め、系内に出現中の列車下部にあた
る発熱軌道プレートと列車の走行していない放熱軌道プ
レートとの温度微係数を求める。また、系内に出現中の
列車の側方領域の気流温度と水蒸気密度微係数、列車の
走行していないトンネル、ホームブランチ部の気流温度
と水蒸気密度とCO2 濃度とのそれぞれの微係数を求め
る。
用し、ルンゲークッタ法による積分計算のコントロール
を行い、各段の微係数演算と温熱積分演算を行うことに
より、系内に出現中の列車内気流温度と水蒸気密度、列
車温度の微係数を求め、系内に出現中の列車下部にあた
る発熱軌道プレートと列車の走行していない放熱軌道プ
レートとの温度微係数を求める。また、系内に出現中の
列車の側方領域の気流温度と水蒸気密度微係数、列車の
走行していないトンネル、ホームブランチ部の気流温度
と水蒸気密度とCO2 濃度とのそれぞれの微係数を求め
る。
【0059】さらに、一般ブランチの気流温度、水蒸気
密度、CO2 濃度のそれぞれの微係数、給気ファン及び
排気ファンブランチの気流温度、水蒸気密度、CO2 濃
度のそれぞれの微係数、空気調和機ブランチの気流温
度、水蒸気密度、CO2 濃度のそれぞれの微係数、チャ
ンバー内の気流温度、水蒸気密度、CO2 濃度のそれぞ
れの微係数、駅、トンネルブランチの壁温度微係数及び
土中温度微係数をを求める(ステップ320〜32
3)。
密度、CO2 濃度のそれぞれの微係数、給気ファン及び
排気ファンブランチの気流温度、水蒸気密度、CO2 濃
度のそれぞれの微係数、空気調和機ブランチの気流温
度、水蒸気密度、CO2 濃度のそれぞれの微係数、チャ
ンバー内の気流温度、水蒸気密度、CO2 濃度のそれぞ
れの微係数、駅、トンネルブランチの壁温度微係数及び
土中温度微係数をを求める(ステップ320〜32
3)。
【0060】(9)単純統計解析処理を行い、平均値計
算により各ブランチの気流温度と水蒸気密度及びCO2
濃度の季節毎1日のプラス側及びマイナス側の平均値を
求め、また、各ブランチの気流温度と水蒸気密度および
CO2 濃度の季節毎1日の最小値及び最大値とそれぞれ
の発生時刻を求め、さらに、熱収支に関する計算を行う
(ステップ324、325)。
算により各ブランチの気流温度と水蒸気密度及びCO2
濃度の季節毎1日のプラス側及びマイナス側の平均値を
求め、また、各ブランチの気流温度と水蒸気密度および
CO2 濃度の季節毎1日の最小値及び最大値とそれぞれ
の発生時刻を求め、さらに、熱収支に関する計算を行う
(ステップ324、325)。
【0061】(10)最後に、アウトプット用データファ
イルとして、続行用初期値データファイル、解析結果デ
ータファイルを出力する(ステップ326)。
イルとして、続行用初期値データファイル、解析結果デ
ータファイルを出力する(ステップ326)。
【0062】前述したステップ320〜323の処理に
よる熱状態解析モデルによる温度、湿度状態の演算は、
トンネル内風速と換気機による換気能力と、トンネル、
駅構内の車両の発熱量と、外気、土中への放熱量とに基
づいて、気流熱量と列車熱量とについて、それぞれ、式
(2)、式(3)の以下のバランス式を解いて、トンネ
ル、駅構内を分割した各区間ごとの発熱量から求めるこ
とができる。
よる熱状態解析モデルによる温度、湿度状態の演算は、
トンネル内風速と換気機による換気能力と、トンネル、
駅構内の車両の発熱量と、外気、土中への放熱量とに基
づいて、気流熱量と列車熱量とについて、それぞれ、式
(2)、式(3)の以下のバランス式を解いて、トンネ
ル、駅構内を分割した各区間ごとの発熱量から求めるこ
とができる。
【0063】 気流熱量の時間変化率(ρcΩ(dT/dt)) =隣接区間から流入する熱量(ρcASV(T1−T)) +列車からの顕熱量(agSg(U−V)(Tg−T) +列車からの放射熱量(εgBSg(Tg−T) +壁面からの顕熱量(aWSW(TW−T)) +壁面からの放射熱量(εWBSW(TW 4−T4)) +機器発熱量(QE) +乗降客からの発熱量(nCQC) +列車からの換気熱量(ncQc) +列車からの冷房放熱量(ρcqR(TR−T)) ……(2) 式(2)において、ρ:気流密度(kg/m3)、c:
気流比熱(J/kgK)、Ω:気流体積(m3)、T:
気流温度(K)、t:時間(秒)、AS:気流流路断面
積(m2)、V:気流速度(m/s)、TI:隣接区間気
流温度(K)、ag:列車表面熱伝達率(W/m2K)、
Sg:列車表面積(m2)、Tg:列車温度(K)、U:
列車速度(m/s)、aW:壁面熱伝達率(W/m
2K)、TW:壁面温度(K)、εg:列車放射率
(−)、εW:壁面車放射率(−)、B:ボルツマン定
数、QE:機器発熱量(W)、nC:乗降客数(人)、Q
C:乗降客発熱量(W)、qR:列車換気量(m3)、
TR:列車内気流温度(K)である。
気流比熱(J/kgK)、Ω:気流体積(m3)、T:
気流温度(K)、t:時間(秒)、AS:気流流路断面
積(m2)、V:気流速度(m/s)、TI:隣接区間気
流温度(K)、ag:列車表面熱伝達率(W/m2K)、
Sg:列車表面積(m2)、Tg:列車温度(K)、U:
列車速度(m/s)、aW:壁面熱伝達率(W/m
2K)、TW:壁面温度(K)、εg:列車放射率
(−)、εW:壁面車放射率(−)、B:ボルツマン定
数、QE:機器発熱量(W)、nC:乗降客数(人)、Q
C:乗降客発熱量(W)、qR:列車換気量(m3)、
TR:列車内気流温度(K)である。
【0064】 列車熱量の時間変化(Ω(dH/dt)) =隣接区間より流入する水蒸気量(AT|U−V|(HT−H)) +列車からの換気放出水蒸気量(qR(HR−H)) +湧水からの水蒸気量(LdWW2) +壁面での蒸発(凝縮)量(LHSWW1) +乗降客供給水蒸気量(nCWC) +列車冷房による放出水蒸気量(WR) ……(3) 式(3)において、H:水蒸気密度(kg/m3)、
AT:気流断面積(m2)、HT:隣接区間水蒸気密度
(kg/m3)、qR:列車換気量(m3)、HR:列車内
水蒸気密度(kg/m3)、L:区間長(m)、dW:湧
水側溝幅(m)、W2:湧水水蒸気量(kg/m3)、L
H:蒸発潜熱(W)、SW:壁表面積(m2)、W1:壁面
蒸発(凝縮)量(kg/m3)、WC:乗降客水蒸気放出
量(kg/m3)、WR:列車冷房放出水蒸気量(kg/
m3)である。
AT:気流断面積(m2)、HT:隣接区間水蒸気密度
(kg/m3)、qR:列車換気量(m3)、HR:列車内
水蒸気密度(kg/m3)、L:区間長(m)、dW:湧
水側溝幅(m)、W2:湧水水蒸気量(kg/m3)、L
H:蒸発潜熱(W)、SW:壁表面積(m2)、W1:壁面
蒸発(凝縮)量(kg/m3)、WC:乗降客水蒸気放出
量(kg/m3)、WR:列車冷房放出水蒸気量(kg/
m3)である。
【0065】第3の処理部3は、詳細には、実行条件デ
ータ、第2の処理部2からの温熱解析結果データ、第1
の処理部1からの気流解析結果データ、ブランチ詳細構
造データ、構造構築データ、続行用の初期値データ、及
び、定数データを入力データとする。
ータ、第2の処理部2からの温熱解析結果データ、第1
の処理部1からの気流解析結果データ、ブランチ詳細構
造データ、構造構築データ、続行用の初期値データ、及
び、定数データを入力データとする。
【0066】実行条件データは、この処理部における処
理期間、第2の処理部からの温熱解析結果データファイ
ル名称、第1処理部1からの気流解析結果データファイ
ル名称、処理結果のプリントアウトの有無を与えるデー
タである。
理期間、第2の処理部からの温熱解析結果データファイ
ル名称、第1処理部1からの気流解析結果データファイ
ル名称、処理結果のプリントアウトの有無を与えるデー
タである。
【0067】第3の処理部3は、前述したデータを受
け、駅、トンネルブランチ毎の壁及び土中の温度、各ブ
ランチ内の気流温度を求め、解析結果データ、プリント
アウトデータとして出力する。解析結果データは、第4
の処理部4で用いる駅、トンネルの壁及び土中の温度、
各ブランチ内の気流温度の6時間毎の時間変化の20年
間にわたるデータである。
け、駅、トンネルブランチ毎の壁及び土中の温度、各ブ
ランチ内の気流温度を求め、解析結果データ、プリント
アウトデータとして出力する。解析結果データは、第4
の処理部4で用いる駅、トンネルの壁及び土中の温度、
各ブランチ内の気流温度の6時間毎の時間変化の20年
間にわたるデータである。
【0068】次に、図4に示すフローを参照して、第3
の処理部の動作を説明する。
の処理部の動作を説明する。
【0069】(1)まず、図示しない実行条件データフ
ァイルより、処理期間、第2の処理部2による温熱解析
結果データファイル名称を読み込む(ステップ40
1)。
ァイルより、処理期間、第2の処理部2による温熱解析
結果データファイル名称を読み込む(ステップ40
1)。
【0070】(2)第2の処理部2の温熱解析結果デー
タ内部ファイルを読み込み、第1の処理部1気流解析結
果データファイル名称と各データを内部変数として格納
し、第1の処理部の気流解析結果データ内部ファイルを
読み込み、構造モデルデータファイル名称と各データを
内部変数に格納する。さらに、構造モデルデータ内部フ
ァイルを読み込み、各データ内部変数に格納し、定数デ
ータファイルを読み込み、各データを内部変数に格納す
る(ステップ402〜404)。
タ内部ファイルを読み込み、第1の処理部1気流解析結
果データファイル名称と各データを内部変数として格納
し、第1の処理部の気流解析結果データ内部ファイルを
読み込み、構造モデルデータファイル名称と各データを
内部変数に格納する。さらに、構造モデルデータ内部フ
ァイルを読み込み、各データ内部変数に格納し、定数デ
ータファイルを読み込み、各データを内部変数に格納す
る(ステップ402〜404)。
【0071】(3)前述で説明した図示しない地下鉄系
構造構築処理部から与えられる構造構築データを格納
し、積分計算初期値設定のために、第2の処理部による
温熱解析結果の最終時刻における各ブランチの壁温度、
土中温度として標準値を設定する(ステップ405〜4
06)。
構造構築処理部から与えられる構造構築データを格納
し、積分計算初期値設定のために、第2の処理部による
温熱解析結果の最終時刻における各ブランチの壁温度、
土中温度として標準値を設定する(ステップ405〜4
06)。
【0072】(4)次に、駅を構成する各ブランチ及び
トンネルブランチの気流速度と前年度の壁温度微係数を
設定し、さらに、駅を構成する各ブランチ、チャンバ
ー、トンネルブランチの熱伝達率を設定する(ステップ
407〜409)。
トンネルブランチの気流速度と前年度の壁温度微係数を
設定し、さらに、駅を構成する各ブランチ、チャンバ
ー、トンネルブランチの熱伝達率を設定する(ステップ
407〜409)。
【0073】(5)前述までに設定された各データを使
用し、ルンゲークッタ法による積分計算のコントロール
を行い、各段の微係数演算と温熱積分演算を行うことに
より、駅、トンネルブランチの壁温度微係数を計算し、
また、駅、トンネルブランチの土中温度微係数を計算す
る(ステップ410〜413)。
用し、ルンゲークッタ法による積分計算のコントロール
を行い、各段の微係数演算と温熱積分演算を行うことに
より、駅、トンネルブランチの壁温度微係数を計算し、
また、駅、トンネルブランチの土中温度微係数を計算す
る(ステップ410〜413)。
【0074】(6)単純統計解析の処理により、駅及び
駅を構成する各ブランチの気流温度の平均値、最大値、
最小値を計算する(ステップ414)。
駅を構成する各ブランチの気流温度の平均値、最大値、
最小値を計算する(ステップ414)。
【0075】(7)最後に、アウトプット用データファ
イルとして、続行用初期値データファイル、解析結果デ
ータファイルを出力する(ステップ415)。
イルとして、続行用初期値データファイル、解析結果デ
ータファイルを出力する(ステップ415)。
【0076】第4の処理部4は、詳細には、実行条件デ
ータ、第1の処理部1からの気流解析結果データ、第2
の処理部2からの温熱解析結果データ、第3の処理部3
からの解析結果データ、構造構築データ、ブランチ詳細
構造データ、列車依存の管理データ、ブランチ依存の管
理データ、駅乗降客データ、駅間乗客データ、外気デー
タ、続行用の初期値データ、及び、定数データを入力デ
ータとする。
ータ、第1の処理部1からの気流解析結果データ、第2
の処理部2からの温熱解析結果データ、第3の処理部3
からの解析結果データ、構造構築データ、ブランチ詳細
構造データ、列車依存の管理データ、ブランチ依存の管
理データ、駅乗降客データ、駅間乗客データ、外気デー
タ、続行用の初期値データ、及び、定数データを入力デ
ータとする。
【0077】実行条件データは、第1の処理部1からの
気流解析結果データファイル名称、第2の処理部2から
の温熱解析結果データファイル名称、第3の処理部3か
らの解析結果データファイル名称、内部ファイル(構造
モデル、外気データ、駅乗降客、駅間乗客データの各フ
ァイル)名称、続行用初期値データ使用の有無、結果プ
リントアウトの有無を与えるデータである。
気流解析結果データファイル名称、第2の処理部2から
の温熱解析結果データファイル名称、第3の処理部3か
らの解析結果データファイル名称、内部ファイル(構造
モデル、外気データ、駅乗降客、駅間乗客データの各フ
ァイル)名称、続行用初期値データ使用の有無、結果プ
リントアウトの有無を与えるデータである。
【0078】続行用の初期値データは、第4の処理部に
よる処理開始時の各ブランチに関連する気流速度と温度
と水蒸気密度とCO2 濃度値、列車に関する温度と水蒸
気密度値、及び、駅、トンネルに関連する温度値を与え
るデータである。なお、他のデータは、第1〜第3の処
理部に対する入力データとして説明したものと同一であ
る。
よる処理開始時の各ブランチに関連する気流速度と温度
と水蒸気密度とCO2 濃度値、列車に関する温度と水蒸
気密度値、及び、駅、トンネルに関連する温度値を与え
るデータである。なお、他のデータは、第1〜第3の処
理部に対する入力データとして説明したものと同一であ
る。
【0079】第4の処理部4は、前述のデータを読み込
み、各ブランチの気流と温熱データとを演算し、プリン
トアウトデータ及び続行用初期値データを出力する。プ
リントアウトデータは、第4の処理部4による演算結果
のグラフ、あるいは、帳票出力のためのブランチ気流、
列車とブランチ及び駅、トンネルの温度と水蒸気密度と
CO2 濃度の時間変化、統計解析(平均値、最大値、最
小値)を与える。
み、各ブランチの気流と温熱データとを演算し、プリン
トアウトデータ及び続行用初期値データを出力する。プ
リントアウトデータは、第4の処理部4による演算結果
のグラフ、あるいは、帳票出力のためのブランチ気流、
列車とブランチ及び駅、トンネルの温度と水蒸気密度と
CO2 濃度の時間変化、統計解析(平均値、最大値、最
小値)を与える。
【0080】次に、図5、図6に示すフローを参照し
て、第4の処理部の動作を説明する。
て、第4の処理部の動作を説明する。
【0081】(1)まず、図示しない実行条件データフ
ァイルより、第1の処理部1による気流解析結果データ
ファイル名称、第2の処理部2による温熱解析結果デー
タファイル名称、第3の処理部3による解析結果データ
ファイル名称、構造モデルデータ内部ファイル名称、外
気データ内部ファイル使用の有無と名称、駅乗降客デー
タ内部ファイル使用の有無と名称、駅間乗客データ内部
ファイル使用の有無と名称、続行用初期値データ使用の
有無、結果プリントアウトの有無を読み込む(ステップ
501)。
ァイルより、第1の処理部1による気流解析結果データ
ファイル名称、第2の処理部2による温熱解析結果デー
タファイル名称、第3の処理部3による解析結果データ
ファイル名称、構造モデルデータ内部ファイル名称、外
気データ内部ファイル使用の有無と名称、駅乗降客デー
タ内部ファイル使用の有無と名称、駅間乗客データ内部
ファイル使用の有無と名称、続行用初期値データ使用の
有無、結果プリントアウトの有無を読み込む(ステップ
501)。
【0082】(2)第1の処理部1による気流解析結果
データ内部ファイル、第2の処理部2による温熱解析結
果データ内部ファイル、及び、第3の処理部3による解
析結果データ内部ファイルを読み込み、各データを内部
変数に格納する。また、構造モデルデータ内部ファイ
ル、ダイヤデータ内部ファイル、駅乗降客データ内部フ
ァイル、及び、駅間乗客データファイルを読み込み、各
データを内部変数に格納する。さらに、駅構内等に設置
される給気ファン、排気ファン、空気調和機等のスケジ
ュール運転用のデータファイルを読み込み、スケジュー
ル運転データを生成する(ステップ502〜509)。
データ内部ファイル、第2の処理部2による温熱解析結
果データ内部ファイル、及び、第3の処理部3による解
析結果データ内部ファイルを読み込み、各データを内部
変数に格納する。また、構造モデルデータ内部ファイ
ル、ダイヤデータ内部ファイル、駅乗降客データ内部フ
ァイル、及び、駅間乗客データファイルを読み込み、各
データを内部変数に格納する。さらに、駅構内等に設置
される給気ファン、排気ファン、空気調和機等のスケジ
ュール運転用のデータファイルを読み込み、スケジュー
ル運転データを生成する(ステップ502〜509)。
【0083】(3)次に、定数データファイル読み込
み、各データを内部変数に格納する(ステップ51
0)。
み、各データを内部変数に格納する(ステップ51
0)。
【0084】(4)前述で説明した図示しない地下鉄系
構造構築処理部から与えられる構造構築データを格納す
る。また、地下鉄系構造構築処理部から与えられる圧力
計算のための連立一次方程式解くために用いる気流流量
係数マトリクスを設定する(ステップ511、51
2)。
構造構築処理部から与えられる構造構築データを格納す
る。また、地下鉄系構造構築処理部から与えられる圧力
計算のための連立一次方程式解くために用いる気流流量
係数マトリクスを設定する(ステップ511、51
2)。
【0085】(5)各ブランチに関連する気流速度と温
度と水蒸気密度とCO2 濃度、列車に関連する温度と水
蒸気密度、及び、駅、トンネルに関連する温度のそれぞ
れに標準値を設定し、また、続行用初期値データファイ
ルを読み込み各ブランチに関連する気流速度と温度と水
蒸気密度とCO2 濃度、列車に関連する温度と水蒸気密
度、及び、駅、トンネルに関連する温度値を設定する
(ステップ513、514)。
度と水蒸気密度とCO2 濃度、列車に関連する温度と水
蒸気密度、及び、駅、トンネルに関連する温度のそれぞ
れに標準値を設定し、また、続行用初期値データファイ
ルを読み込み各ブランチに関連する気流速度と温度と水
蒸気密度とCO2 濃度、列車に関連する温度と水蒸気密
度、及び、駅、トンネルに関連する温度値を設定する
(ステップ513、514)。
【0086】(6)前述までの処理でデータファイルの
入力処理と初期設定の処理とを終了し、次に、各時刻毎
(1秒毎)の列車登録、列車実施形態設定、ブランチ情
報の設定を行う。これらのデータは、図示しない列車運
行管理処理部から得られるデータであり、各1秒毎に、
列車登録処理として、運用列車が系内に出現中である
か、出現前あるいは消滅後であるかを登録し、列車情報
として、存在列車について、運行中の線路指標、基準位
置からの距離、列車の速度と加速度、及び、列車と軌道
プレートとの発熱量を設定し、さらに、ブランチ情報設
定として、列車の存在するブランチ番号を設定する(ス
テップ516〜518)。
入力処理と初期設定の処理とを終了し、次に、各時刻毎
(1秒毎)の列車登録、列車実施形態設定、ブランチ情
報の設定を行う。これらのデータは、図示しない列車運
行管理処理部から得られるデータであり、各1秒毎に、
列車登録処理として、運用列車が系内に出現中である
か、出現前あるいは消滅後であるかを登録し、列車情報
として、存在列車について、運行中の線路指標、基準位
置からの距離、列車の速度と加速度、及び、列車と軌道
プレートとの発熱量を設定し、さらに、ブランチ情報設
定として、列車の存在するブランチ番号を設定する(ス
テップ516〜518)。
【0087】(7)次に、境界データ設定として、系内
接続点の圧力と系内接続点を通過する列車関連データを
設定し、さらに、外気の温度と水蒸気密度、駅に存在す
る乗降客数と駅間を走行する列車に存在する乗客数、各
ブランチ壁面とチャンバー壁面、列車表面、及び、軌道
プレート表面の熱伝達率を設定する(ステップ519〜
522)。
接続点の圧力と系内接続点を通過する列車関連データを
設定し、さらに、外気の温度と水蒸気密度、駅に存在す
る乗降客数と駅間を走行する列車に存在する乗客数、各
ブランチ壁面とチャンバー壁面、列車表面、及び、軌道
プレート表面の熱伝達率を設定する(ステップ519〜
522)。
【0088】(8)次に、給気ファン、排気ファン、空
気調和機等に対する1日の営業時間に依存したON/O
FFデータ、1年の季節に依存したON/OFFデー
タ、及び、空気調和機に対する補正データとして、風量
割合の設定及び熱交換割合の設定を行う(ステップ52
3〜526)。
気調和機等に対する1日の営業時間に依存したON/O
FFデータ、1年の季節に依存したON/OFFデー
タ、及び、空気調和機に対する補正データとして、風量
割合の設定及び熱交換割合の設定を行う(ステップ52
3〜526)。
【0089】(9)前述までに設定された各データを使
用し、ルンゲークッタ法による積分計算のコントロール
を行い、各段の微係数演算と気流、温熱積分演算を行う
ことにより、列車が存在するブランチ内の列車側方領域
毎の列車ゆらぎ比速度微係数を演算すると共に、列車が
存在するブランチ内の列車側方領域毎の列車圧力差を演
算し、さらに、トンネル・ホームブランチ、一般ブラン
チ、給気ファンブランチ、排気ファンブランチ、及び、
空気調和機ブランチの気流流量微係数を演算する(ステ
ップ527〜530)。
用し、ルンゲークッタ法による積分計算のコントロール
を行い、各段の微係数演算と気流、温熱積分演算を行う
ことにより、列車が存在するブランチ内の列車側方領域
毎の列車ゆらぎ比速度微係数を演算すると共に、列車が
存在するブランチ内の列車側方領域毎の列車圧力差を演
算し、さらに、トンネル・ホームブランチ、一般ブラン
チ、給気ファンブランチ、排気ファンブランチ、及び、
空気調和機ブランチの気流流量微係数を演算する(ステ
ップ527〜530)。
【0090】(10)次に、気流流量係数を計算し、連立
一次方程式によるアクティブカラム法により各ポイント
の圧力を演算する(ステップ531)。
一次方程式によるアクティブカラム法により各ポイント
の圧力を演算する(ステップ531)。
【0091】(11)ステップ531による演算結果圧力
に基づいて、各ブランチの気流流量微係数、系内に出現
中の列車の列車内気流温度と水蒸気密度との微係数、列
車温度の微係数、系内に出現中の列車下部にあたる発熱
軌道プレートの温度微係数と、列車の走行していない放
熱軌道プレートの温度微係数、系内に出現中の列車の側
方領域の気流温度の微係数と水蒸気密度微係数、列車の
走行していないトンネル、ホームブランチ部の気流温
度、水蒸気密度及びCO2 濃度のそれぞれの微係数、一
般ブランチの気流温度と水蒸気密度及びCO2 濃度微係
数、給気ファン及び排気ファンブランチの気流温度、水
蒸気密度及びCO2 濃度のそれぞれの微係数、空気調和
機ブランチの気流温度、水蒸気密度及びCO2 濃度のそ
れぞれの微係数、チャンバー内の気流温度、水蒸気密度
及びCO2 濃度のそれぞれの微係数、駅、トンネルブラ
ンチの壁温度微係数及び土中温度微係数を求める(ステ
ップ532)。
に基づいて、各ブランチの気流流量微係数、系内に出現
中の列車の列車内気流温度と水蒸気密度との微係数、列
車温度の微係数、系内に出現中の列車下部にあたる発熱
軌道プレートの温度微係数と、列車の走行していない放
熱軌道プレートの温度微係数、系内に出現中の列車の側
方領域の気流温度の微係数と水蒸気密度微係数、列車の
走行していないトンネル、ホームブランチ部の気流温
度、水蒸気密度及びCO2 濃度のそれぞれの微係数、一
般ブランチの気流温度と水蒸気密度及びCO2 濃度微係
数、給気ファン及び排気ファンブランチの気流温度、水
蒸気密度及びCO2 濃度のそれぞれの微係数、空気調和
機ブランチの気流温度、水蒸気密度及びCO2 濃度のそ
れぞれの微係数、チャンバー内の気流温度、水蒸気密度
及びCO2 濃度のそれぞれの微係数、駅、トンネルブラ
ンチの壁温度微係数及び土中温度微係数を求める(ステ
ップ532)。
【0092】(12)単純統計解析処理を行い、平均値計
算により各ブランチの気流温度と水蒸気密度及びCO2
濃度の季節毎1日のプラス側及びマイナス側の平均値を
求め、また、各ブランチの気流温度と水蒸気密度および
CO2 濃度の季節毎の1日の最小値及び最大値とそれぞ
れの発生時刻を求め、さらに、熱収支に関する計算を行
う(ステップ533、534)。
算により各ブランチの気流温度と水蒸気密度及びCO2
濃度の季節毎1日のプラス側及びマイナス側の平均値を
求め、また、各ブランチの気流温度と水蒸気密度および
CO2 濃度の季節毎の1日の最小値及び最大値とそれぞ
れの発生時刻を求め、さらに、熱収支に関する計算を行
う(ステップ533、534)。
【0093】(13)最後に、アウトプット用データファ
イルとして、続行用初期値データファイル、解析結果デ
ータファイルを出力する(ステップ535)。
イルとして、続行用初期値データファイル、解析結果デ
ータファイルを出力する(ステップ535)。
【0094】前述した本発明による空気調和・換気シミ
ュレータによれば、入力データとして、トンネル、駅構
内の各ブランチの構造、伝熱、発熱に関する条件、空気
調和設備、換気設備の設備データ及び時系列の運転条
件、列車の構造、伝熱、運動、時系列な発熱に関する条
件、トンネル、駅構内の各ブランチの複数の列車の時系
列な運転状況、外気の時系列な温度、湿度の条件、土中
の土質、地下水の流速等のデータを与えることにより、
1日の路線全体のトンネル、駅構内の各ブランチの気流
速度、圧力を予測することができる。
ュレータによれば、入力データとして、トンネル、駅構
内の各ブランチの構造、伝熱、発熱に関する条件、空気
調和設備、換気設備の設備データ及び時系列の運転条
件、列車の構造、伝熱、運動、時系列な発熱に関する条
件、トンネル、駅構内の各ブランチの複数の列車の時系
列な運転状況、外気の時系列な温度、湿度の条件、土中
の土質、地下水の流速等のデータを与えることにより、
1日の路線全体のトンネル、駅構内の各ブランチの気流
速度、圧力を予測することができる。
【0095】また、1日の路線全体のトンネル、駅構内
の各ブランチの気流速度の解析結果を元に、1年間の路
線全体のトンネル、駅構内の各ブランチの温度、湿度を
予測でき、さらに、1日の路線全体のトンネル、駅構内
の各ブランチの気流速度の解析結果、1年間の路線全体
のトンネル、駅構内の各ブランチの温度、湿度の解析結
果に基づいて、20年間にわたる駅構造物及び土中への
蓄熱を予測することができる。
の各ブランチの気流速度の解析結果を元に、1年間の路
線全体のトンネル、駅構内の各ブランチの温度、湿度を
予測でき、さらに、1日の路線全体のトンネル、駅構内
の各ブランチの気流速度の解析結果、1年間の路線全体
のトンネル、駅構内の各ブランチの温度、湿度の解析結
果に基づいて、20年間にわたる駅構造物及び土中への
蓄熱を予測することができる。
【0096】さらに、駅構内の各ブランチの構造、伝
熱、発熱に関する条件及び空気調和設備、換気設備の設
備データとを前記解析結果と等価になるように詳細に細
分し、設定した解析対象日の前記解析結果を入力条件と
して駅構内各ブランチの1日の間の気流速度、圧力、温
度、湿度、CO2濃度を予測することができる。また、
この気流速度に基づき、乗降する客が持ってくる紛じん
量、その大きさ等を仮定することにより、粉じん濃度を
も予測することができる。
熱、発熱に関する条件及び空気調和設備、換気設備の設
備データとを前記解析結果と等価になるように詳細に細
分し、設定した解析対象日の前記解析結果を入力条件と
して駅構内各ブランチの1日の間の気流速度、圧力、温
度、湿度、CO2濃度を予測することができる。また、
この気流速度に基づき、乗降する客が持ってくる紛じん
量、その大きさ等を仮定することにより、粉じん濃度を
も予測することができる。
【0097】前述した本発明の実施形態は、本発明を地
下鉄道システムに適用したものとして説明したが、本発
明は、地下街等の地下空間、あるいは、極めて密閉度の
高い空間における環境シミュレーションに対しても適用
することができる。
下鉄道システムに適用したものとして説明したが、本発
明は、地下街等の地下空間、あるいは、極めて密閉度の
高い空間における環境シミュレーションに対しても適用
することができる。
【0098】
【発明の効果】以上説明したように本発明によれば、地
下空間、特に、地下鉄道におけるトンネル、駅構造を構
成する各ブランチ内の気流速度、圧力、温度、湿度、C
O2 濃度、粉じん濃度を高い精度で予測することがで
き、これを利用することにより、列車風抑制のためのト
ンネル、駅構造上の設計指針を得ることができると共
に、熱負荷、CO2 濃度、粉じん濃度の予測結果に基づ
いて、空気調和設備、換気設備等の適切な設備容量の設
定を行うことができる。
下空間、特に、地下鉄道におけるトンネル、駅構造を構
成する各ブランチ内の気流速度、圧力、温度、湿度、C
O2 濃度、粉じん濃度を高い精度で予測することがで
き、これを利用することにより、列車風抑制のためのト
ンネル、駅構造上の設計指針を得ることができると共
に、熱負荷、CO2 濃度、粉じん濃度の予測結果に基づ
いて、空気調和設備、換気設備等の適切な設備容量の設
定を行うことができる。
【図1】本発明の一実施形態による空気調和・換気シミ
ュレータの機能構成を示すブロック図である。
ュレータの機能構成を示すブロック図である。
【図2】第1の処理部である気流・圧力解析処理部の処
理動作を説明するフローチャートである。
理動作を説明するフローチャートである。
【図3】第2の処理部である温熱解析処理部の処理動作
を説明するフローチャートである。
を説明するフローチャートである。
【図4】第3の処理部である土中温度・壁温度演算処理
部の処理動作を説明するフローチャートである。
部の処理動作を説明するフローチャートである。
【図5】第4の処理部である気流速度・温湿度・CO2
濃度詳細演算処理部の処理動作を説明するフローチャー
ト(その1)である。
濃度詳細演算処理部の処理動作を説明するフローチャー
ト(その1)である。
【図6】第4の処理部である気流速度・温湿度・CO2
濃度詳細演算処理部の処理動作を説明するフローチャー
ト(その2)である。
濃度詳細演算処理部の処理動作を説明するフローチャー
ト(その2)である。
【図7】本発明の実施形態でシミュレーションの対象と
する地下空間の構造例を示す図である。
する地下空間の構造例を示す図である。
1 気流・圧力解析処理部 2 温熱解析処理部 3 土中温度・壁温度演算処理部 4 気流速度・温湿度・CO2 濃度詳細演算処理部 11 排気ファン 12 給気ファン 13 ホーム 14 コンコース 15 空気調和機 16 冷凍機 17 制御装置 18a、18b 温度センサ 19 トンネル 20 排気口 21 給気口 22 階段
Claims (5)
- 【請求項1】 地下空間における環境状態の予測を行う
空気調和・換気シミュレータにおいて、地下空間を構成
するトンネル、階段、通路、換気孔等のブランチによる
地下空間構造モデルデータ、前記地下空間内で空気を移
動させる運動物体の運動データを入力データとして、前
記各ブランチの気流速度を求める第1の処理部と、前記
空間構造モデルデータ、前記第1の処理部からの解析結
果データ、運動物体の熱特性に関するデータ、外気温、
外気湿度、外気二酸化炭素濃度データ、前記空間へ出入
りする客データを入力として、各ブランチの壁温度、土
中温度の初期値と短期変動量及び各ブランチ内の温度、
水蒸気密度、二酸化炭素濃度の初期値と境界値を演算す
る第2の処理部と、前記空間構造モデルデータ、第1の
処理部からの気流解析結果データ、第2の処理部からの
温熱解析結果データを入力として、各ブランチの壁、土
中の温度、各ブランチの気流温度の演算処理を行う第3
の処理部と、前記空間構造モデルデータ、第1の処理部
からの気流解析結果データ、第2の処理部からの温熱解
析結果データ、第3の処理部からの壁温度、土中温度デ
ータを入力として、各ブランチの気流の温度、水蒸気密
度、二酸化炭素濃度の時間変化データを演算する気流速
度、温湿度、二酸化炭素濃度の詳細な演算処理を行う第
4の処理部とを備えることを特徴とする地下空間におけ
る空気調和・換気シミュレータ。 - 【請求項2】 前記第1の処理部は、前記各ブランチの
気流速度を1日分の時間変化情報として出力することを
特徴とする請求項1記載の地下空間における空気調和・
換気シミュレータ。 - 【請求項3】 前記第2の処理部は、1年間の季節の経
過に応じた各ブランチ内の気流温度、気流湿度、ブラン
チ壁温度、土中温度を時間変化情報として出力すること
を特徴とする請求項1または2記載の地下空間における
空気調和・換気シミュレータ。 - 【請求項4】 前記第3の処理部は、20年間にわたる
ブランチ壁温度、土中温度を時間変化情報として出力す
ることを特徴とする請求項1、2または3記載の地下空
間における空気調和・換気シミュレータ。 - 【請求項5】 前記地下空間が地下鉄道を構成する地下
空間であり、前記地下空間構造モデルデータが駅構造モ
デルデータであり、前記運動物体がトンネル内を走行す
る車両であることを特徴とする請求項1ないし4のうち
1記載の地下空間における空気調和・換気シミュレー
タ。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP8127982A JP2765630B2 (ja) | 1996-05-23 | 1996-05-23 | 地下空間における空気調和・換気シミュレータ |
US08/771,163 US5822227A (en) | 1996-05-23 | 1996-12-20 | Air-conditioning and ventilation simulator in subsurface space |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP8127982A JP2765630B2 (ja) | 1996-05-23 | 1996-05-23 | 地下空間における空気調和・換気シミュレータ |
US08/771,163 US5822227A (en) | 1996-05-23 | 1996-12-20 | Air-conditioning and ventilation simulator in subsurface space |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH09311617A true JPH09311617A (ja) | 1997-12-02 |
JP2765630B2 JP2765630B2 (ja) | 1998-06-18 |
Family
ID=26463785
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP8127982A Expired - Fee Related JP2765630B2 (ja) | 1996-05-23 | 1996-05-23 | 地下空間における空気調和・換気シミュレータ |
Country Status (2)
Country | Link |
---|---|
US (1) | US5822227A (ja) |
JP (1) | JP2765630B2 (ja) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010112576A (ja) * | 2008-11-04 | 2010-05-20 | Takasago Thermal Eng Co Ltd | 地下空間に対する連絡通路 |
JP2013210112A (ja) * | 2012-03-30 | 2013-10-10 | East Japan Railway Co | 地下駅の空調システム |
CN114484815A (zh) * | 2021-12-20 | 2022-05-13 | 彩虹显示器件股份有限公司 | 一种tft-lcd光电玻璃成型间环境控制系统及方法 |
CN114660937A (zh) * | 2022-03-29 | 2022-06-24 | 西南交通大学 | 二氧化碳约束下的车内压力波动迭代学习控制方法及系统 |
CN115727506A (zh) * | 2022-11-21 | 2023-03-03 | 珠海格力电器股份有限公司 | 空调器自适应控制方法、空调器以及计算机可读存储介质 |
JP2023511478A (ja) * | 2020-12-30 | 2023-03-20 | 中南大学 | 列車車両用空気調和方法、装置、記憶媒体及びプログラム製品 |
CN116951690A (zh) * | 2023-07-24 | 2023-10-27 | 中建安装集团有限公司 | 一种地下室微风道防潮系统控制方法、介质及系统 |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6866203B2 (en) | 2003-02-26 | 2005-03-15 | Unico, Inc. | Method and apparatus for sizing an environmental control system |
FR2866945B1 (fr) * | 2004-03-01 | 2006-05-19 | Florence Tantot | Systeme et procede de commande d'equipements de conditionnement d'ambiance dans une enceinte |
AT7319U3 (de) * | 2004-08-19 | 2005-12-15 | Roman Wagner & Partner Gmbh | Verfahren zur bewertung der raumqualität |
KR100851520B1 (ko) * | 2005-12-28 | 2008-08-11 | 한국건설기술연구원 | 통계적 분석기법을 이용한 지하 공간의 대기 환경 예측장치 및 방법 |
CN104864528B (zh) * | 2015-05-07 | 2016-12-28 | 广东申菱环境系统股份有限公司 | 一种地铁站竖井型蒸发冷凝直膨冷风型通风空调系统 |
CN104864527B (zh) * | 2015-05-07 | 2016-10-12 | 广东申菱环境系统股份有限公司 | 一种地铁站坑道型蒸发冷凝直膨冷风型通风空调系统 |
CN104864529B (zh) * | 2015-05-07 | 2016-10-12 | 广东申菱环境系统股份有限公司 | 一种地铁站机房型蒸发冷凝直膨冷风型通风空调系统 |
CN106338127B (zh) * | 2016-09-20 | 2018-06-22 | 珠海格力电器股份有限公司 | 用于地铁暖通空调系统的负荷预测和控制系统及其方法 |
CN108444076B (zh) * | 2018-03-08 | 2020-03-06 | 广东美的制冷设备有限公司 | 推送方法、空气调节设备、移动终端及存储介质 |
CN108488987B (zh) * | 2018-03-08 | 2020-04-07 | 广东美的制冷设备有限公司 | 空气调节设备的控制方法、存储介质和设备 |
CN108469855A (zh) * | 2018-05-16 | 2018-08-31 | 贵州省原态上品农业科技有限责任公司 | 一种袖珍菇菇房环境信息采集及控制系统 |
CN109615966A (zh) * | 2018-12-29 | 2019-04-12 | 珠海格力电器股份有限公司 | 基于虚拟现实技术的空调器安装与维修实训方法和系统 |
CN109799858B (zh) * | 2019-01-28 | 2021-01-26 | 浙江理工大学 | 一种夯土建筑智能监测系统的监测参数训练方法 |
CN110595667B (zh) * | 2019-09-06 | 2020-05-12 | 江苏航天大为科技股份有限公司 | 一种检测轨道交通通风系统风压的方法 |
KR102230856B1 (ko) * | 2019-11-01 | 2021-03-23 | (주)한국필터시험원 | 지하철 승강장 공기정화시스템 |
CN111322735B (zh) * | 2020-02-24 | 2021-11-05 | 上海理工大学 | 过渡季节有效增大屏蔽门系统地铁车站新风量的方法 |
ES2850377A1 (es) * | 2020-02-27 | 2021-08-27 | Sener Ing & Sist | Sistema de control de ventilación para un subterráneo |
CN111930165B (zh) * | 2020-08-14 | 2022-01-25 | 珠海格力电器股份有限公司 | 环境模拟系统、环境模拟系统的控制方法及装置 |
CN114861423B (zh) * | 2022-04-22 | 2024-06-21 | 西北核技术研究所 | 一种热辐射作用下地介质动态响应数值模拟方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3982330A (en) * | 1975-08-06 | 1976-09-28 | Bituminous Coal Research, Inc. | Method and apparatus for simulating in a model the disbursement of gas and dust in a mine |
US4552059A (en) * | 1984-09-18 | 1985-11-12 | Cambridge Engineering, Inc. | Flow measurement for exhaust-type canopy and ventilating hood |
EP0179625B1 (en) * | 1984-10-19 | 1990-03-07 | Nippondenso Co., Ltd. | Air conditioner for automobiles |
US4591093A (en) * | 1985-07-02 | 1986-05-27 | E-Zee Company | Calibration apparatus for air flow controllers |
US5042470A (en) * | 1989-05-30 | 1991-08-27 | Nozomi Kanesaka | Ventilating system for respiration of a patient |
US5226818A (en) * | 1991-11-25 | 1993-07-13 | Aai Corporation | Firefighter training system with thermal agent detection |
FR2702035B1 (fr) * | 1993-02-24 | 1995-05-12 | Gec Alsthom Transport Sa | Dispositif de régulation de climatisation d'un local. |
US5590062A (en) * | 1993-07-02 | 1996-12-31 | Matsushita Electric Industrial Co., Ltd. | Simulator for producing various living environments mainly for visual perception |
KR0160410B1 (ko) * | 1994-03-31 | 1999-01-15 | 김광호 | 공기조화기및 그 제어방법 |
-
1996
- 1996-05-23 JP JP8127982A patent/JP2765630B2/ja not_active Expired - Fee Related
- 1996-12-20 US US08/771,163 patent/US5822227A/en not_active Expired - Fee Related
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010112576A (ja) * | 2008-11-04 | 2010-05-20 | Takasago Thermal Eng Co Ltd | 地下空間に対する連絡通路 |
JP2013210112A (ja) * | 2012-03-30 | 2013-10-10 | East Japan Railway Co | 地下駅の空調システム |
JP2023511478A (ja) * | 2020-12-30 | 2023-03-20 | 中南大学 | 列車車両用空気調和方法、装置、記憶媒体及びプログラム製品 |
CN114484815A (zh) * | 2021-12-20 | 2022-05-13 | 彩虹显示器件股份有限公司 | 一种tft-lcd光电玻璃成型间环境控制系统及方法 |
CN114484815B (zh) * | 2021-12-20 | 2023-12-26 | 彩虹显示器件股份有限公司 | 一种tft-lcd光电玻璃成型间环境控制系统及方法 |
CN114660937A (zh) * | 2022-03-29 | 2022-06-24 | 西南交通大学 | 二氧化碳约束下的车内压力波动迭代学习控制方法及系统 |
CN114660937B (zh) * | 2022-03-29 | 2023-04-14 | 西南交通大学 | 二氧化碳约束下的车内压力波动迭代学习控制方法及系统 |
CN115727506A (zh) * | 2022-11-21 | 2023-03-03 | 珠海格力电器股份有限公司 | 空调器自适应控制方法、空调器以及计算机可读存储介质 |
CN116951690A (zh) * | 2023-07-24 | 2023-10-27 | 中建安装集团有限公司 | 一种地下室微风道防潮系统控制方法、介质及系统 |
Also Published As
Publication number | Publication date |
---|---|
US5822227A (en) | 1998-10-13 |
JP2765630B2 (ja) | 1998-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2765630B2 (ja) | 地下空間における空気調和・換気シミュレータ | |
Yuan et al. | CFD simulation and optimization of the ventilation for subway side-platform | |
Riederer et al. | Room thermal modelling adapted to the test of HVAC control systems | |
Axley | Multizone airflow modeling in buildings: History and theory | |
Yang et al. | An innovative environmental control system of subway | |
Cejudo et al. | Physical and neural network models of a silica-gel desiccant wheel | |
Liu et al. | CFD prediction of surface condensation on walls and its experimental validation | |
Kosonen et al. | The effect of typical buoyant flow elements and heat load combinations on room air temperature profile with displacement ventilation | |
Liu et al. | A network model for natural ventilation simulation in deep buried underground structures | |
Papakonstantinou et al. | Computational analysis of thermal comfort: the case of the archaeological museum of Athens | |
Jokisalo et al. | A comparison of measured and simulated air pressure conditions of a detached house in a cold climate | |
Krishnan et al. | Development of a small-scale test facility for effectiveness evaluation of fixed-bed regenerators | |
Chew et al. | Improving thermal model predictions for naturally ventilated buildings using large eddy simulations | |
JP4456007B2 (ja) | 地下空間の空調システム | |
Coydon et al. | Energy performance of façade integrated decentralised ventilation systems | |
Fruhwirt et al. | Development of a methodology for studying tunnel climate in long railway tunnels and for optimizing the design process of cross-passage cooling systems | |
Howell et al. | On the natural displacement flow through a full-scale enclosure, and the importance of the radiative participation of the water vapour content of the ambient air | |
Sambolek | Model testing of road tunnel ventilation in normal traffic conditions | |
Karlsson | Air infiltration through building entrances | |
Kanam et al. | A review on underground mine ventilation system | |
Wang et al. | Comparison of Simulation Programs for Airflow and Smoke Movement during High-Rise Fires. | |
Desta et al. | Combining CFD and data-based mechanistic (DBM) modelling approaches | |
Papakonstantinou et al. | Numerical simulation of CO2 dispersion in an auditorium | |
Xu et al. | Simulation of mixed convection flow in a room with a two-layer turbulence model. | |
Papakonstantinou et al. | Mathematical modeling of environmental conditions inside historical buildings. The case of the archaeological museum of Athens |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |