平成
25 年度
修士論文
大型輸送車両の空気抵抗低減装置の
実験及び開発
指導 教員
筒井 康賢 教授
高知工科大学
知能機械システム工学科
筒井研究室
学生番号
1165045
矢野 佑樹
第
1 章 序論
1.1 研究背景
1.1.1 走行抵抗 1.1.2 境界層とレイノルズ数 Re1.2 研究目的
第
2 章 従来の研究
第
3 章 実験装置
3.1 風洞
3.2 計測機器
3.2.1 差圧計 DSA3217/16Px 3.2.2 トラバース装置 3.2.3 主流・後流計測用ピトー管3.3 風洞模型
3.4 空気抵抗低減装置
3.4.1 エアディフレクター 3.4.2 空気抵抗低減装置第
4 章 風洞実験
4.1 実験環境
4.2 実験手順
4.3 実験内容
4.3.1 計側部
4.3.2 実験
目次
第
5 章 実験結果及び考察
5.1 空気抵抗係数
5.2 考察
第
6 章 結論
謝辞
参考文献
1
第 1 章 序論
1.1 研究背景 日本における貨物輸送は図 1(1)に示すように現在大きな問題を抱えている。トラック類 (営業用トラック・自家用トラック)が日本の物流サービスの約 9 割を占めており、一般 的な生活必需品から、工業用製品・材料の供給には欠かせない存在となっている。(参考文 献 8 より輸送機関別輸送量を参考としている)また、近年のネットショッピングの増加に より今後の利用拡大なども考えられる。しかし、トラック運送事業は、図 2 に示すように 近年の軽油価格高騰で、2004 年~2013 年の期間に軽油価格が全国平均で約 50 円の値上が りとなっており、経営状況は悪化している。ここで燃料消費量の抑制を実現し、トラック 運送事業者の経営状況改善に貢献したいと考えた。 図1 トラック業界の現状(1) 図2 近年の軽油価格の推移2 1.1.1 走行抵抗 走行時の車両に働く抵抗には大別して次の4 種類がある。 走行抵抗=転がり抵抗 + 空気抵抗 + 勾配抵抗 + 加速抵抗 転がり抵抗とは走行時におけるタイヤの変形や、地面との摩擦などにより生じる抵抗で あり、一般的に以下の式で求められる。
𝑅
𝑟
= μ
r
Mg
𝑅𝑟 : 転がり抵抗[N] μr : 転がり抵抗係数[-] M : 車両質量[N] g : 重力加速度 空気抵抗とは車両が走行する際に周囲の空気の作用によって生じる抵抗のことで、大別 して圧力抵抗と摩擦抵抗に分類される。特に圧力抵抗が空気抵抗の大部分を占めており、 その中でも車両形状を起因とする形状抵抗が空気抵抗の主要素である。以下に空気抵抗を 求める際の一般的な式を示す。𝐷 =
1
2
ρAV
2
C
𝐷
𝐷 : 空気抵抗[N] ρ : 空気密度[kg/m3] A : 車両前面投影面積[m2] V : 車両走行速度[m/s] 𝐶𝐷 : 空気抵抗係数[-] 勾配抵抗とは図 3 のように、坂道など傾斜がある道路を走行するときに生じる抵抗のこ とである。以下に勾配抵抗を求める際の一般的な式を示す。𝑅
𝑠
= Mg sin θ
𝑅𝑠 : 転がり抵抗[N] θ : 路面傾角 式(1.1) 式(1.2) 式(1.3) )3 図 3 傾斜面における勾配抵抗 加速抵抗とは、運転手によって車両が加速される際に生じる抵抗のことである。車両の 特性というよりも、運転手の個人差によるところが大きい。 1.1.2 境界層とレイノルズ数 Re 1.1.2.1 境界層 主流速度U の流体中に設置された物体の表面では、流体の粘性によって流体が付着し、 速度が0 となる。他方、物体表面から十分離れた地点では、粘性の影響はほとんどみられ なくなり、主流速度U となる。以上から、物体表面の地点から、物体表面から十分離れた 地点までの間に速度勾配が生じるのは容易に理解できる。この速度勾配が顕著に現れる物 体表面近傍を境界層と呼ぶ。 境界層には層流境界層と乱流境界層が存在する。例として一様な主流速度U の流体中に 設置された平行平板を考える。図4 に示す様に、層流境界層では速度勾配が大きく、乱流 境界層では速度勾配が小さい。乱流境界層内部では流れの撹拌作用が大きく、主流速度U を有する流体を境界層内部に引き込むからである。 また、物体表面から境界層が剥がれる剥離と呼ばれる現象が生じることがある。剥離は 境界層内部の流体の速度が0 となり生じる。この点を剥離点と呼び、剥離点後方の剥離領 域では流れが逆流し、それを起因とする渦が生じ、非常に複雑な流れとなる。この領域の ことを死水領域とも呼ぶ。死水領域では主流圧力に比べ、圧力が低下していることが知ら れている。この圧力差が物体に作用する空気抵抗となって現れてくる。層流境界層で剥離 することを層流剥離、乱流境界層で剥離することを乱流剥離と呼ぶ。
4 1.1.2.2 レイノルズ数 Re 流れが層流(一様)であるか乱流(乱れている)であるかを評価する際に、無次元量であるレ イノルズ数Re を使用する。Reは以下の式で表される。
𝑅
𝑒
=
UL
ν
𝑅𝑒 : レイノルズ数[-] U : 主流相対速度[m/s] L : 代表長さ[m] ν : 流体の動粘度[m2/s] 平行平板の場合、代表長さは主流方向の平板の長さである。 流れが層流から乱流に遷移する Reを臨界レイノルズ数 Recと呼び、実際の車両周りの流 れは乱流であるので、風洞試験を行うときはRecを超える様に注意する必要がある。 図 4 平行平板における境界層内の流速分布(速度勾配)のイメージ 式(1.4) )5 1.2 研究目的 1.1.1 で示した走行抵抗の各成分は、車両重量と走行速度に大きく依存しているのがわか る。本研究では、走行速度に大きく依存している空気抵抗に注目する。走行抵抗は約80 km/h 未満の低速走行時では、走行抵抗の中では転がり抵抗の影響が大きいが、約80 km/h を超 える高速走行時においては、空気抵抗が車速の2 乗に比例するので影響が大きくなる。 研究対象となる大型輸送車両の多くは、高速道路など高速走行が必要とされる場所を走 行する場合が多い。以上のことを踏まえ、本研究の目的は既存の大型トラックに対して空 気抵抗低減装置を取り付け、装置の導入コストを限りなく抑制し、燃費向上を目指すこと である。また車両制限令(2)により原則として、車両の幅2.5 m、高さ 3.8 m、長さ 12 m 以 内に取り付けることを考慮する。
6
第 2 章 従来の研究
従来の研究としては図 5、6(3)に示すように、トラック荷台の後端にフラップや巨大な流
線型の空気抵抗低減装置を取り付けるなどその形状で空気の流れを変化させるような受動 的な空気抵抗低減装置が多い。この例としては2011 年にメルセデス・ベンツがコンセプト モデルとして示した"Mercedes Benz Aero Trailer Concept"などがあげられると考えられ る。この"Mercedes Benz Aero Trailer Concept"に関しては、空気抵抗を約 18%改善し、燃 費を約5%低減することが可能である。(図7(4)に"Mercedes Benz Aero Trailer Concept" を示す)
図5 Transtex Composites folding rear trailer deflector(3)
7
8 また、現在使用されている空気抵抗低減装置としては、運転台(キャブ)上部に設置さ れるエアディフレクタ・運転台(キャブ)と貨物部の間に設置されるギャップシールド・ シャーシ側面に設置されるサイドスカートといったものがあげられる。設置位置及び外観 を図8 にて示す。 図8 空気抵抗低減装置の設置位置及び外観(5)
9
第 3 章 実験装置
3.1 風洞 実験にて使用した風洞の仕様・特性を表1 に、外観を図 9 に示す。本実験に使用した風 洞は2012 年に高知工科大学に設置されたものである。 表 1 研究にて使用した風洞の仕様・特性 主要仕様・特性 全長 mm 10232 吸い込み口断面積 mm2 7680000 吹き出し口断面積 mm2 1000000 図9 研究に使用した風洞の外観10 3.2 計測機器 実験にて使用した差圧計、トラバース装置、ピトー管の主要仕様・特性を以下に示す。 3.2.1 差圧計 DSA3217/16Px 差圧計DSA3217/16Px(以下 DSA)の主要仕様・特性を表 2 に、外観を図 10(5)に示す。 表 2 差圧計 DSA3217/16Px の主要仕様・特性 主要仕様・特性 ピエゾ抵抗式圧力センサ[Ch] 16 圧力レンジ 差圧[kPa] ± 2.5、7.0、17、35、100、200、350、700、 1725、3500、5200 絶対圧[kPa] 100、200、350、700 精度 センサ圧力レンジ[kPa] 精度[% F.S.] ±2.5 ±0.20 ±7、±17 ±0.12 ±35~3500 ±0.05 ~±5200 ±0.08 分解能[bit] 16 スキャンレート[Hz/Ch] 500 接続 Ethernet 10baseT 接続プロトコル TCP/IP、UDP コネクタタイプ 入出力 RJ-45 電源 Bendix PTO6A-8-3S-SR、3 ピンメス トリガ Bendix PTO6A-8-6S-SR、6 ピンメス 電源 28Vdc nominal @ 400mA(20-36Vdc)
外部トリガ 6.5mA @ 9Vdc minimum edge sensing 許容圧力 2.5[kPa] = 13.79[kPa] 7[kPa] = 35[kPa]
1379[kPa] = 200[%] 3447[kPa] = 150[%] 4137[kPa] = 125[%] 5171[kPa] = 100[%] 最大リファレンス圧[kPa] 1724 動作温度[℃] 0~60 温度補正レンジ[℃] 0~60(標準) 温度補正レンジ内の温度影響[% F.S./℃] ±0.001 重量[kg] 2.9
11 図10 A3217/16Px の外観(5) 3.2.1.1 差圧計測概要 DSA の差圧計測の概要を説明する。図 11(5)に DSA 上面を示す。リファレンス孔に基準 圧力を与え、16Ch(チャンネル)ある計測孔に任意の圧力を与える。各チャンネルに入力さ れる圧力と基準圧の差が計測値となる。ピトー管を例に挙げると、ピトー管の静圧孔をリ ファレンス孔に接続し、全圧孔を計測孔に接続すると、計測値として動圧が得られる。 図 81 A3217/16Px の上面(5)
12
3.2.1.2 データの送受信概要
DSA は Ethernet、いわゆる LAN ケーブル接続でパソコンとデータの送受信を行う。接 続プロトコルはTCP/IP、UDP の 2 種類あるが、本研究では TCP/IP を使用した。送受信 されるデータの形式はASCII とバイナリの 2 種類あるが、ASCII 形式を使用した。
13 3.2.2 トラバース装置 トラバース装置の主要仕様・特性を表3 に、外観を図 12 に示す。使用したトラバース装 置は、計測機器を取り付け、上下・左右の 2 次元方向に任意に移動させること目的とした ものである。 表 3 トラバース装置の主要仕様・特性 主要仕様・特性 横軸 縦軸 ステッピングモータ AR46AA AR69MA モータ回転分解能[°/パルス] 0.36 0.36 軸送り分解能[mm] 0.01 0.01 ローター慣性モーメント[kg・m2] 58×10-7 870×10-7 励磁最大静止トルク[N・m] 0.3 2.0 許容スラスト荷重[N] 4.6 16.7 図12 トラバース装置の外観
14 3.2.3 主流・後流計測用ピトー管 主流計測用ピトー管及び後流計測用ピトー管の主要仕様・特性を表4 に、外観を図 13・ 14 に示す。 表 4 ピトー管の主要仕様・特性 主流計測用ピトー管 主要仕様・特性 ツクバリカセイキ(株)製 F-202 JIS 型ピトー管 管外径[mm] 6 全圧孔径[mm] 3 静圧孔径[mm] 1 後流計測用ピトー管 主要使用・特性 ツクバリカセイキ(株)製 管外形[mm] 3 全圧孔径[mm] 2 静圧孔径[mm] 1 図13 主流計測用ピトー管の外観 図14 後流計測用ピトー管の外観
15 3.3 風洞模型 本研究では風洞模型の形状決定のため、トラック製造会社が掲示する大型トラックの外 観寸法を参考にした。また模型のスケールに関してはトラック長手方向を代表長さとおき、 球や平板の臨界レイノルズ数Rec = 3.0 × 105、5.0 × 105を十分超えるようにレイノルズ数 Re ≧1.0×106を目標とした。本研究に使用した風洞模型に関しては、1.1.2 の境界層とレイ ノルズ数で示している式1.4 より計算した結果 Re=1.1×105となっている。表5 に主要外 観寸法を、図15 に風洞モデルの外観を示す。 表 5 風洞模型の主要外観寸法 主要仕様・特性 全長[mm] 795 全幅[mm] 166 全高[mm] 252 図15 風洞模型の外観
16 3.4 空気抵抗低減装置 ここでは本研究に使用した空気抵抗低減装置について説明する。 ・3.4.1 エアディフレクター 本研究では空気抵抗低減装置として、ギャップシールド付きのエアディフレクタを図 16 に示す。 図16 上:エアディフレクタ外観 中:装着時側面 下:装着時外観
17 3.4.2 空気抵抗低減装置 空気抵抗低減装置として、今回はNACA0024 翼型をベースとして前縁をなだらかにし たものを使用している。トラック後方部では周辺に比べ圧力が低下していると考えられ るので、トラック後方部分に流れを意図的に流すことにより、圧力を向上させることを 目的とする。 本研究で使用した空気抵抗低減装置として、以下の図17 に示す。 図17 空気抵抗低減装置
18
第 4 章 風洞実験
4.1 実験環境 第3 章で説明した実験機器を風洞計測部にて設置した概略図を図 18 に、表 6 に概略図内 の対応記号を示す。 図18 風洞計測部概略図 表 6 風洞計測部概略図対応記号 風洞計測部番号 ① 風洞吹出口 ⑤ テストベッド ② 主流計測用ピトー管 ⑥ 後流計測用ピトー管 ③ 主流計測用ピトー管用固定フレーム ⑦ トラバース装置 ④ 風洞モデル19
4.2 実験手順
本実験ではトラック後端と後流計測用ピトー管との間隔を100 mm とり、主流計測用ピ トー管の静圧Pa を基準圧力として、主流総圧 Pall、後流総圧 Pwall および後流静圧 Pwa との差圧を差圧計(ScanivalveDSA3217)で計測した。実験手順を図 19 に示す。
20 4.3 実験内容 4.3.1 計側部 本実験ではトラック後方での圧力変化を計測し考察を行っている。後流計測用ピトー管 をトラバース装置にマウントし、横軸に600 mm・縦軸に 300 mm を 10 mm 間隔で稼 働させることにより、稼働範囲内の計 1800 点を計測している。また、本実験に関して は後流計測用ピトー管を6 本使用しており 1 回の計測時間としては約 25 分となってい る。後流計測用ピトー管マウント及びマウント後の状態に関して図20 にて示す。 図20 後流計測用ピトー管マウント部及びマウント後の状態
21 ・4.3.2 実験 本実験では 3.4.2 で示した空気抵抗低減装置を使用して実験を行った。本実験に関 しては、①トラックのみ ②貨物部後方上面に設置した状態 ③②の状態に翼端板を装 着した状態 ④貨物部後方上面及び側面に設置した状態 ⑤④の状態の角を丸めた状態 の計5つの実験を行った。以上の実験に使用したモデル及び空気抵抗低減装置に関して 以下に示す。 図21 トラックのみ 図22 貨物部後方上面に設置した状態及び空気抵抗低減装置
22
図23 貨物部後方上面に空気抵抗低減装置及び翼端板設置
23
24
第 5 章 実験結果及び考察
5.1 空気抵抗係数 実験結果を評価する際に空気抵抗係数CDを使用する。式(1.2)を CDについて式変形した ものを以下に示す。 抗力係数は、𝐶
𝐷
=
1𝐷
2𝜌𝑆𝑈
2 ここですべての圧力を風洞出口の静圧を基準に計測しているので、1
2
𝜌𝑈
2
= 𝑃
1
,
1
2
𝜌𝑢
2
=𝑃
2
―
𝑃
3
となり、
𝐶
𝐷
=
1
𝑆
∫ {2 (√
𝑃
2𝑃
−𝑃
1−
𝑃
2−𝑃
𝑃
1) −
𝑃
𝑃
1} 𝑑𝑆
2
ρ : 空気密度[kg/m3] S : 車両前方投影面積[m2] Ua : 主流速度[m/s] 𝑃1 : 主流動圧[Pa] 𝑃2 : 後流計測点動圧[m/s] 𝑝3 : 後流計測点静圧[Pa] d𝑆2 : 計測点代表面積[m2] であり、実験結果から得られたトラック後流の圧力分布を積分して空気抵抗D を求める。 以上の計算より実際に求めたCDを表7 に示す。また、同じ条件で解析を行った場合の CD を表8 に示す。 表 7 空気抵抗低減装置の空気抵抗係数 Cdへの影響・比較 状態 Cd[-] 改善[%] トラックのみ 0.580 貨物部後方上面に設置 0.518 10 貨物部後方上面に空気抵抗低減装置及び翼端板設置 0.524 9 貨物部後方上面及び側面に設置 0.470 18 貨物部後方上面及び側面に設置し角を丸めたものを設置 0.469 19 式(5.1) ) 式(5.2) )25 表 8 空気抵抗低減装置の空気抵抗係数 CDへの影響(解析) 状態 Cd[-] 改善[%] トラックのみ 0.370 貨物部後方上面に設置 0.346 6.5 貨物部後方上面に空気抵抗低減装置及び翼端板設置 0.318 14 貨物部後方上面及び側面に設置 0.311 16 貨物部後方上面及び側面に設置し角を丸めたものを設置 0.322 13 次に計測面であるトラック後流の差圧(Pwall-Pa)を主流動圧(Pall-Pa)で無次元化し、コン ター図表記したものを以下の図に示す。横軸は計測開始点からのトラック幅方向の距離 X mm、縦軸はトラック高さ方向の距離 Y mm である。主流動圧に対して後流圧力の変動が 大きいと、その計測点において圧力が低いことを意味する。また、比較として同じ実験条 件で解析を行ったものを以下の図に示す。 実験結果 図 269 トラックのみでの後方圧力
26
図27 貨物部後方上面に設置時の後方圧力
27
図29 貨物部後方上面及び側面に設置
28 解析結果
図31 トラックのみでの解析結果
29
図33 貨物部後方上面に空気抵抗低減装置及び翼端板の解析結果
30 5.2 考察 本実験の結果において空気抵抗低減装置非装着時のトラック後方での圧力分布に比べ 空気抵抗低減装置を装着したすべての結果において、圧力分布の高さ方向 y=20~25・横 方向x=25~35 の範囲で圧力上昇が確認できた。また、これは図 35 に示すように、空気抵 抗低減装置の作用によりトラック後方での負圧となっている領域に流れができたために、 圧力向上が図られた結果だと考えられる。(図36:実際に実験を行った際の流れをタフト 法で確認した結果) 実験結果より求めた空気抵抗係数CDに関してはトラックのみの場合に比べすべての結 果で向上していることが確認できた。(表7)この結果、空気抵抗低減装置を使用するこ とにより、空気抵抗を改善できていることが確認できた。実験での空気抵抗低減装置を上 面のみ・上面及び翼端版を設置した2 つに関して、空気抵抗低減装置のみを使用した際に 翼端渦の影響が懸念されたため、翼端板を使用した結果、空気抵抗係数 CD値が翼端板を 使用したもののほうが大きくなってしまった。しかし解析結果での CD値は翼端板を使用 したほうが改善されていた。この結果より、空気抵抗低減装置自体に翼端渦が発生してい るため、翼端板を使用したほうが良いと考えられるが、本実験では翼端板に厚紙を使用し ていたため解析結果との違いが出たと考えられる。 図35 空気抵抗低減装置によるトラック後方への風の流れイメージ 上:非装着時 下: 装着時
31 図36 タフト法 本実験に関しては、空気抵抗低減装置を貨物部後方上面に設置したトラック後方の圧力 分布は高さ方向y=20~25 のみに変動が見て取れた。それに比べ、貨物部後方上面及び側面 に設置したものは、高さ方向 y=20~25 での変動並びに多少ではあるが高さ方向 y=100~ 200 mm の範囲で左右にも変動を確認することができた。ただし、各圧力分布の結果で分布 図高さ方向y=20~25 の左右に角のようなものが残った。これに関しては設置した空気抵抗 低減装置に対して翼端渦が発生しているためだと考えられた。(図38:赤で囲っている部分) そこで、空気抵抗低減装置及び翼端板を設置することで他の結果に比べ翼端渦の影響が小 さくなっていることが確認できた。また、翼端板を設置したうえで残っている影響に関し ては、翼端板自体の影響だと考えられる。上面及び側面に空気抵抗低減装置を設置した場 合、その角を埋めることにより翼端渦の発生を抑えられるような結果が出た。(図 37 青で 囲っている結果) しかし、上面及び側面に空気抵抗低減装置を設置した実験では、高さ方向 y=0~10・横 方向 x=20~25・35~40 の範囲で圧力が低下しており、側面に設置した空気抵抗低減装置 の地面側(図38:青で囲っている部分)に翼端渦が発生している可能性があると考えられ た。そのためにこの部分にも翼端板を設置するか、丸めるといった改善が必要だと考えら れる。比較部分を図37 に示す(赤で囲っているものが翼端板を使用した結果となっている)。
32
図37 上部の圧力分布の比較
33
第 6 章 結論
本研究では空気抵抗低減装置としてNACA0024 をベースとしたものを使用し、5 パター ンの実験を行った。トラックに対して空気抵抗低減装置を使用した結果、圧力向上を確認 することができた。また、実験結果より求めたCD値ではトラックのみに比べ空気抵抗低減 装置を装着したもので最大19%向上していることが分かった。(表7)この結果より空気抵 抗低減装置の優位性を見いだせたと考えられる。 空気抵抗低減装置自体の抵抗として翼端渦が発生して影響していることがわかった。そ こで、翼端板を使用したが、翼端渦による抵抗を改善することができなかった。また、ト ラック前方などの角にも翼端渦が発生している可能性があると考えられるので、今後の問 題であると考えられる。 結果として、本実験において空気抵抗低減装置を使用することにより圧力が向上する結 果が得られた。また、空気抵抗係数CD値に関しても大きな改善を図ることができた。ここ で、本研究より考えられた今後の問題及び対策を以下に示す。 ① 空気抵抗低減装置の寸法を改善する (本研究では効果を確認したかったので寸法を考慮していない。そのため、車両制限 令に収まっていないので、今後この寸法内に収める必要がある) ② 空気抵抗低減装置に働く抵抗への対策を行う (翼端板以外にも対策を考える) ③ 空気抵抗低減装置の最適な取付け位置を把握するためにトラック上面の流れを確認す る ④ 空気抵抗低減装置の最適形状を CFD(数値流体力学)の解析ソフトなどを使用して検 討を行う ⑤ 風洞試験と解析の結果を比較し、①、②、③、④の妥当性を評価する34
謝辞
本研究を行うに当たり、本研究全体にわたって指導していただきました高知工科大学 筒井康賢先生にこころより深く御礼申しあげます。 本研究を行うにあたり、共に研究活動を行った小笠原慧さん、廣江綾斗さん、野口秀幸 さんに深く感謝すると共に、実験設備の設置などに協力してくださった筒井研究室のメン バーに深く感謝します。35
参考文献
(1) 公益社団法人全日本トラック協会
http://www.jta.or.jp/coho/kigyobutsuryu/kigyobutsuryu.html (2012)
(2) e-Gov、車両制限令、第 3 条、http://law.e-gov.go.jp/htmldata/S36/S36SE265.html (2011) (3) Jason Leuschen、Kevin R. Cooper,Full-Scale Wind Tunnel Tests of Production and
Prototype, Second-Generation Aerodynamic Drag-Reducing Devices for
Tractor-Trailers,http://www.jta.or.jp/coho/kigyobutsuryu/img12/kigyobutsuryu12.pdf (2006)
(4) "Mercedes Benz Aero Trailer Concept"
http://blog.livedoor.jp/motersound/archives/51680826.html (5) ISUZU
http://blog.livedoor.jp/motersound/archives/51680826.html
(6) R. Spivey,R. Hewitt,H. Othman and T. Corke, Flow Separation Control on Trailing Edge Radiiusing Single Dielectric Barrier Discharge Plasma Actuators: An Application to VehicleDrag Control, Lecture Notes in Applied and Computational Mechanics Vol 41,pp.136-149 (2009)
(7) Scanivalve Inc.,DSA3217 Hardware Manual、
http://www.scanivalve.com/products/pressure-measurement/ethernet-intelligent-pr essure-scanners/dsa3217-pressure-scanner-gas-measurement/、(2010)
(8) 日本内航海運組合総連合会