This document is provided by JAXA.
This document is provided by JAXA.
宇宙航空研究開発機構特別資料
JAXA Special Publication
第 78 回 風洞研究会議論文集
Proceedings of the Wind Tunnel Technology Association 78
thmeeting
2007 年 8 月
August 2007
宇宙航空研究開発機構
Japan Aerospace Exploration Agency
JAXA-SP-07-005
This document is provided by JAXA.
This document is provided by JAXA.
目 次
第 78 回風洞研究会議(平成 19 年 5 月 17、18 日)
1.小型超音速実験機開発における風洞実験の紹介(資料のみ) ‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥ 1 大貫 武(JAXA)
2.旅客機開発と音響風洞試験 ‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥ 12 柴田 眞(JAXA)
3.JAXA6.5m× 5.5m低速風洞における
PSP
計測システムの研究開発 ‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥ 19 満尾 和徳、栗田 充、口石 茂、藤井 啓介、渡辺 重哉(JAXA)4.JAXA風洞技術開発センター傾斜計校正装置と内装式傾斜計の応用 ‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥ 26 青木 良尚、細江 信幸、星野 秀雄、伊藤 健(JAXA)
5.0.5m極超音速風洞におけるマッハ 5AGARD-B標準模型 6 分力試験 ‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥ 33 津田 尚一、小山 忠勇、平林 則明、渡利 實(JAXA)
関根 英夫(JAST)、木伏 淳子(スペースサービス)、中村 晃祥(JAST)
6.JAXA高エンタルピ風洞における光学的気流計測技術の現状 ‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥‥ 40 水野 雅仁、伊藤 健、藤田 和央、石田 清道(JAXA)
長井 遵正(AES)、藤井 啓介(JAXA)
This document is provided by JAXA.
This document is provided by JAXA.
第 78 回 風洞研究会議論文集 1
This document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 2
This document is provided by JAXA.
第 78 回 風洞研究会議論文集 3
This document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 4
This document is provided by JAXA.
第 78 回 風洞研究会議論文集 5
This document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 6
This document is provided by JAXA.
第 78 回 風洞研究会議論文集 7
This document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 8
This document is provided by JAXA.
第 78 回 風洞研究会議論文集 9
This document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 10
This document is provided by JAXA.
第 78 回 風洞研究会議論文集 11
This document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 12
This document is provided by JAXA.
第 78 回 風洞研究会議論文集 13
This document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 14
This document is provided by JAXA.
第 78 回 風洞研究会議論文集 15
This document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 16
This document is provided by JAXA.
第 78 回 風洞研究会議論文集 17
This document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 18
This document is provided by JAXA.
1. はじめに
燐光物質の酸素消光現象を利用した感圧塗料計測
(Pressure-Sensitive Paint: PSP)と呼ばれる表面圧力場計 測技術が注目をされている.従来の静圧孔を用いた計測 では限られた離散的なデータしか得られないが,PSPを 用いると模型全面の詳細なデータが取得できる.また,
塗装するだけで良いので手間がかからず,安価に計測す ることができる利点がある.すでに,超音速/遷音速域 の高速流れ場における
PSP
計測は多くの実績があり,実用化されつつある[1], [2].JAXAの
PSP
技術は国産小型 旅客機開発に利用され,設計用圧力データを提供してい る[3].JAXA
では,低速における航空機性能向上を目的とし て,航空機の離着陸特性の改善に大きな効果がある高揚 力装置に注目して研究を行なっている.その研究の一環 として,JAXA6.5 m× 5.5m低速風洞(LWT1)におい て,感圧塗料(PSP)による高揚力半裁装置模型の表面 圧力場計測を実施した[4]-[6].JAXA大型低速風洞における初めての試みである.PSP計測は,模型全体の詳細な 圧力場(流れ場)情報を得ることができるため,先進風 洞計測ツールとして注目されている.
低速における計測では
PSP
の発光が弱く,また圧力の 変化が小さいため計測が困難であり,まだ実用レベルに は達していない.さらに,PSPは圧力感度を有するだけ ではなく,温度によっても発光強度が変化するため,圧 力変化の小さい低速では温度の影響を強く受ける[7]-[9].そこで,本研究では大型低速風洞に
PSP
計測システ ムを適用するため,高出力LED
励起照明を開発し,さ らに温度感度の低いPSP
を使用することにより,計測 精度向上を図った.PSP発光強度の圧力値への変換は,静圧孔データを参照した
In-situ
法を用いた.本文では,まず,PSP塗装環境,励起
LED
光源,カ メラシステムなど大型低速風洞対応型PSP
システムに ついて説明する.次に,PSPの計測精度および模型表面 圧力場の迎角および速度依存性に関する結果を紹介す る.JAXA 6.5m × 5.5m 低速風洞における PSP 計測システムの研究開発
満尾 和徳、栗田 充、口石 茂、藤井 啓介、渡辺 重哉
(宇宙航空研究開発機構 総合技術研究本部)
Research and Development of PSP Measurement System at JAXA 6.5m x 5.5m Low-Speed Wind Tunnel
Kazunori M ITSUO , Mitsuru K URITA , Shigeru K UCHI-ISHI , Keisuke F UJII and Shigeya W ATANABE
( Japan Aerospace Exploration Agency / Institute of Aerospace Technology )
Abstract
Pressure-Sensitive Paint (PSP) measurement system was developed in the JAXA 6.5m x 5.5m low-speed wind tunnel. Pressure images on a high-lift-device (HLD) model, which was built for R&D for next generation of civil transport aircraft, were measured by the low-speed PSP system. The dependencies of pressure patterns on flow speed and angle of attack were investigated. Pressure distribution peculiar to the HLD model were clearly visualized. The measurement accuracy of the present low-speed PSP system was approximately 0.16 -
0.2in Cp at 60m/s. These results indicated that this PSP system was a practical measurement tool to acquire pressure images on an aerodynamic model at low-speed.
第 78 回 風洞研究会議論文集 19
This document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 20
2. 低速 PSP 計測システム 2.1
PSP 塗装系PSP
は圧力感度だけではなく同時に温度感度を有する ため,温度による計測誤差を生じる.特に,圧力変化の 小さい低速ではその影響を強く受ける.そこで,本実験 では温度感度の低いISSI
社製のFIB-PSP
を使用し,温 度による計測誤差を軽減した.図 1 に特性を示す.圧力 感度の線形性が良く,温度感度は 1%/℃未満である.LWT1 で使用する模型(図 2)は巨大であるため,模
型塗装時も大掛かりな装置が必要となる.通常 1m程度 の模型であれば,塗装専用のブース内で塗装を行うが,2mを超えるサイズの模型の場合,専用ブースで塗装で きないため風洞カート内で塗装しなければならない.図 3 に示すように,模型周りを透明なカーテンで囲い,塗 料が飛散しないように養生して塗装した.PSP塗装は,
完全防備スーツを装着した作業員 2 人で行った.塗装作 業以外に,マスキングや
PSP
計測用マーカー取り付け などを含めて準備に 2 日を要した(塗装乾燥の週末 2 日 間を除く).2.2
光学系高迎角をなす高揚力装置模型の上面全体を計測するた め,カートの上流と下流に
CCD
カメラを 2 台設置し,PSP
計測を行なった.その 2 台を用いて主翼全体を計測 するグローバル計測と,重要な箇所を局所的に計測する クローズアップ計測の 2 種類の計測を行なった.図 4 の 模式図が示すように,16bit-CCDカメラ(HAMAMATSUPHOTONICS, ORCA-II-BT104)2 台と高出力 LED
励起 照明 2 台をカートに配置した.本実験で使用したPSP
0.7 0.75 0.8 0.85 0.9 0.95 1 1.05 1.1 1.15 1.2
0.70 0.75 0.80 0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 P/Pref (Pref=100kPa)
Iref/I
0℃
5℃
10℃
15℃
20℃
25℃
30℃
(a) 圧力感度特性
0.9 0.95 1 1.05 1.1 1.15 1.2
-0.2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
T/Tref (Tref=20℃)
I/Iref
@ 100kPa
(b) 温度感度特性
図
1
PSP の圧力・温度感度特性Single-slotted flap Slat
Engine nacelle
Inboard Slat
4.9 m 2.3 m
Double-slotted flap FTF Single-slotted flap Slat
Engine nacelle
Inboard Slat
4.9 m 2.3 m
Double-slotted flap FTF
図
2
JAXA 高揚力装置半裁標準模型の鳥瞰図図
3
PSP 塗装作業及び,PSP 塗装模型の写真(a) PSP
塗装(b)
マーカー取り付け(c) PSP
塗装された模型This document is provided by JAXA.
第 78 回 風洞研究会議論文集 21
の励起帯は 380-530nmにあり,発光ピークは 650nmに ある.PSPの発光のみを捉えるため
CCD
カメラの前面 に 650 ± 20nmバンドパスフィルタを設置して計測し た.大型風洞で光学系を設置する場合,手間と人手がかか るため,予め光学系の配置を検討しておくと作業効率が 良くなる.PSP計測では図 5 に示すような
CAD
データ を活用した電子モックアップを使用して光学系の配置を 事前に決定した.通風時と無風時における取得画像枚数はそれぞれ 64 枚で,計測に要した時間は各々 8 分程度であった.画像 処理で画質を評価した結果,16 枚程度で
S/N
が十分高 い処理画像が得られることを確認したので,処理に使用 する枚数は 16 枚とした.PSP
の発光を強めるため高出力LED
励起照明を作製 した(図 6,7,8).405nmにピークをもつ 1024 個のLED
素子を並べて大型化したものである.これまで使 用していた 300WのXe
光源よりも十分明るい.LEDが 発する熱による発光特性の変化を抑えるため,ファンに よりLED
基板の背面を空冷している.さらに,LED発 光強度を安定に作動させるため,フォトダイオード(PD)によるフィードバック制御機能を備えている.
2.3
計測方法PSP
画像から圧力画像を得るためには,通風中のPSP
画像と無風時のPSP
画像が必要である.先にも述べた ように,PSPは温度感度を有するため,通風時と無風時 の間に温度差がない方が良い.そこで,本試験では通風 画像取得後,直ちに風洞を停止し,通風直後のPSP
画 像を計測することでPSP
の温度依存性を軽減した[10].CC D CCD
LED Illuminator LED Illuminator
Optical Window Flow Direction
CC D CC
D CCD
CCD
LED Illuminator LED Illuminator
Optical Window Flow Direction
図
4
PSP 光学系図
5
LWT1 ストラットカートおよび模型の 電子モックアップ図
6
LED 励起照明により照らされた模型Photo-Diode
図
7
LED 励起照明図
8
LED 励起照明の内部構造Optical filter for cutting IR emission Cooling fan LEDs board
Air intake Air intake
Fin for cooling 1024 LEDs
Exhaust air Photo Diode
Optical filter for cutting IR emission Cooling fan LEDs board
Air intake Air intake
Fin for cooling 1024 LEDs
Exhaust air Photo Diode
This document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 22
試験の手順は以下の通りである.試験条件(風速,迎 角)設定後,模型温度の安定を待ってから通風中の
PSP
画像計測と静圧孔データを取得した.模型には温度セン サーが内蔵されており,模型温度が安定してから計測を 開始した.通風中の計測が終了して風洞停止後,無風時 におけるPSP
画像を計測する.1 ケース毎にこのシーケ ンスを繰り返して実験を行なった.2.4
画像処理画像処理は,これまで遷音速,超音速風洞試験におい て実績のあるプログラムを基本とし[11],画像フィルタ など,低速用に一部改良したプログラムを使用した.
PSP
画像を圧力に変換する手法は,圧力孔を参照するIn-situ
法を用いた.3. 実験結果
風速
U=60m/sec,迎角α=15deg
の実験結果を図 9 に 示す.カラーバーの赤色が高圧,青色が低圧を表してい る.赤褐色の部位は胴体などデータのない領域を表して いる.スラット後方の母翼前縁に強い負圧領域がみられ,さらに外舷フラップ外側にも低圧領域がみられる.また,
翼端では渦による低圧パターンが確認できる.
PSP
画像から圧力への変換は,静圧孔を参照したIn- situ
法を用いた.そのIn-situ
曲線を図 10 に示す.フィ ッティングカーブは 2 次式を用いた.フィッティングの 精度(RMS:Root-Mean-Square)は Cp
換算で約 0.16 であった.風速 60m/sにおける他の迎角のケースも約 0.16 〜 0.2 であった.データのバラつきは主としてPSP
の温度依存性によるものである.HLD模型の主翼は複 数の部位から構成されており,構造が複雑であるため温 度分布の不均一さを生じやすい.そのため,In-situ曲線 に主翼全体の静圧孔を参照すると,温度による計測誤差 の影響を受ける.しかし,翼面上の圧力分布はCp
換算 で約 5 の範囲で分布するので,0.2 の分解能があれば十 分に圧力分布を可視化することができる.図 11 に
PSP
と静圧孔データの比較を示す.S21 〜 S101 はそれぞれ静圧孔列の位置を表している.主翼全
体に亘ってPSP
データは静圧孔データと良く一致して いるのがわかる.しかしながら,局所的に少し差が大き くなっている箇所がある.たとえば,S61 圧力孔列の後 縁部分をみると,PSP
と静圧孔の差が大きくなっている.これは,翼下面にフラップ収納スペースがあり,断面形 状がこの部分で急激に変化するため温度の影響を強く受 けるからである.このような温度による誤差を解消し,
計測精度を高めるため,JAXAでは圧力と温度を同時に 計測できる複合
PSP
の開発を進めている[12].圧力分布の迎角依存性を調べるため,流速
U=60m/s,
迎角α=5, 10, 15degの
PSP
計測結果を図 12 に示す.迎 角が大きくなるに従い,スラット上面に負圧の領域が現 れはじめ,またスラット後方の母翼前縁に生じるサクシ ョンピークが強くなる.さらに,翼端では渦による低圧 領域がみられ,迎角により形成される低圧領域は変化し ている.図 13 にPSP
と静圧孔の比較を示す.迎角が変 わってもPSP
データは静圧孔データとよく合っている のがわかる.S21 line
S22 line S40 line
S61 line
S62 line
S101 line
図
9
主翼上の静圧孔位置(U=60m/s, α=15deg)図
10
In-situ 較正カーブ(U=60m/s, α=15deg)RMS=0.16 in Cp
This document is provided by JAXA.
第 78 回 風洞研究会議論文集 23
次に,圧力分布の速度依存性を図 14,15 に示す.速 度が遅くなるにしたがい画像の
S/N
が悪くなっている.30m/sの動圧は 60m/sの 4 分の 1 であり,PSP発光強度 変化が 4 分の 1 になるためノイズの影響が大きくなる.
U=30m/s
における計測精度はRMS
値で約 0.4(in Cp)であった.風速 30m/sではノイズが大きく計測精度は期 待できないが,可視化としては十分使用できる.
4. まとめ
高揚力形態旅客機模型の主翼表面圧力場を計測するた
め,JAXA6.5m× 5.5m低速風洞用
PSP
システムを開発 し,その実用性を評価した.PSP塗装系,光源/カメラ を含む光学系は適正に機能し,画像処理手法も有効であ ることが確認できた.本PSP
システム用いることによ り,風速 60m/sにおいてRMS
値約 0.2(Cp換算)の精 度で圧力場を計測することができる.また,高揚力装置 模型特有の圧力場を鮮明に可視化することができ,本PSP
計測システムは低速試験において有効な圧力場計測 ツールであることが実証された.x x
x x
x x
Cp Cp
Cp Cp
Cp Cp
図
11
PSP データと静圧孔データの比較(U=60m/s, α=15deg)S21 S22
S40 S61
S62 S101
This document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 24
(a) α=15deg (b) α=10deg (c) α=5deg
1 [Cp]
0
− 1
− 2
− 3
− 4
− 5
1 [Cp]
0
− 1
− 2
− 3
− 4
− 5
1 [Cp]
0
− 1
− 2
− 3
− 4
− 5
図
12
圧力分布の迎角依存性(U=60m/s)(a) 60m/s (b) 40m/s (c)30m/s
図
14
圧力分布の速度依存性(α=10 deg)図
13
PSP データと静圧孔データの比較(S21 line , U=60m/s)x x x
(a) α=15deg (b) α=10deg (c) α=5deg
Cp Cp Cp
1 [Cp]
0
− 1
− 2
− 3
− 4
− 5
1 [Cp]
0
− 1
− 2
− 3
− 4
− 5
1 [Cp]
0
− 1
− 2
− 3
− 4
− 5
This document is provided by JAXA.
第 78 回 風洞研究会議論文集 25
参考文献
[1] Bell, J.H, Schairer, E. T., Hand, L. A and Mehta, R. D.,
“Surface Pressure Measurements Using Luminescent Coatings,” Annu. Rev. Fluid Mech., 33(2001), pp.155-
206.[2] Liu, T., Campbell, B. T., Burns, S. P. and Sullivan, J. P.,
“Temperature- and Pressure-Sensitive Luminescent Paints in Aerodynamics,” Appl. Mech. Rev., 50-4, pp.227-246,
1997.[3] Nakakita, K., Kurita, M., Mitsuo, K. and Watanabe, S.,
“Practical pressure-sensitive paint measurement system for industrial wind tunnels at JAXA,” Meas. Sci.
Technol.
17No.
2, February, 2006, pp.359-366.
[4] Ito, T., Yokokawa, Y., Ura, H, Kato, H., Mitsuo, K. and Yamamoto, K., “Height-Lift Device Testing in JAXA
6.5m x 5.5m Low-speed Wind Tunnel,” 25th AIAAAerodynamic Measurement Technology and Ground Testing Conference, AIAA 2006-3643, San Francisco, California,.
5-8June
2006.[5]
満尾和徳,栗田充,口石茂,藤井啓介,渡辺重哉,伊藤正剛:
JAXA
高揚力形態旅客機模型の低速風 洞試験(低速PSP
計測システムの開発),第 44 回 飛行機シンポジウム,大宮ソニックシティー,2006 年 10 月.[6] Mitsuo, K., Kurita, M., Kuchi-Ishi, S., Fujii, K., Ito, T.
and Watanabe, S., “PSP Measurement of a High-Lift- Device Model in JAXA 6.5m
× 5.5m Low-Speed WindTunnel,” AIAA-2007-1065, Reno, Nevada, 2007.
[7]
坂上博隆,満尾和徳,中北和之:感圧塗料技術の近 年の動向について,可視化情報学会誌-特集記事,Vol.24, No.95,pp.218-223,
2004 年.[8] Mitsuo, K., Nakakita, K. and Kurita, M., “Application of Pressure-Sensitive Paint to Low-Speed Wind Tunnel Testing at Japan Aerospace Exploration Agency,” 24th International Congress of the Aeronautical Sciences, ICAS
2004-3.2.3, Yokohama, Japan, 2004.[9] Sant, Y. Le, Bouvier, F., Merienne, M. C. and Peron, J.
L., “Low Speed Tests using PSP at ONERA, ” 39th AIAA Aerospace Sciences Meeting & Exhibit, AIAA
2001-0555, Reno, Nevada, 2001.[10] Bell, J. H, “Applications of Pressure sensitive Paint to Testing at Very Low Flow Speeds,” 42nd AIAA Aerospace Sciences Meeting & Exhibit, Reno, Nevada, AIAA-2004-0878,
2004.[11]
栗田充,満尾和徳,口石茂,中北和之,藤井啓介,渡辺重哉:
JAXA
実用風洞におけるPSP
計測法処 理技術,第 34 回可視化情報シンポジウム,工学院 大学,2006 年 7 月 25 日.[12] Mitsuo, K., Kurita, M., Nakakita, K. Fujii, K. and Watanabe, S., “Development of Bi-Luminophore Pressure-Sensitive Paint Systems,” 22nd International Congress on Instrumentation in Aerospace Simulation Facilities, ICIASF'07,
2007.図
15
PSP データと静圧孔データの比較(S21 line , α=10 deg)x x x
(a) 60m/s (b) 40m/s (c) 30m/s
Cp Cp Cp
This document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 26
1. はじめに
風洞試験における重要な計測項目の一つに,模型の姿 勢角がある.模型の姿勢角は一般的に重力を基準とした 角度と一様流の偏流角から求められる.サーボ内装式傾 斜計は,重力を利用して傾斜角を計測するセンサであり 多くの産業分野で使用されているが,風洞試験における 重力を基準としたピッチ角とロール角の計測にも使用さ れているセンサである[1].JAXA6.5m× 5.5m低速風洞 でもサーボ内装式傾斜計
Sensorex41400 を購入し,風洞
試験に使用してきた.2003 年 8 月に完成したJAXA
風 洞技術開発センター傾斜計校正装置の概要と共に,課題 であったサーボ内装式傾斜計のロール角を含めた 2 軸校 正法とサーボ内装式傾斜計の設置誤差の補正法につい て,標準模型試験におけるサーボ内装式傾斜計の応用例 と合わせて報告する.2. JAXA 風洞技術開発センター傾斜計校正装置[2]
JAXA
風洞技術開発センター傾斜計校正装置は,風洞 試験における模型姿勢角計測精度向上を目的として導入 されたサーボ内装式傾斜計を正確に校正するために製作 が進められ,2003 年 8 月に完成した.装置の概観を図 1 に,主な仕様を表 1 に示す.この装置はピッチ角とロール角の 2 自由度変角機能を備え,装置の角度検出器と被 校正傾斜計の出力を比較する事による静止状態における 角度校正を行うことが可能である.恒温槽も備えており,
温度ドリフト確認機能を持つと共に,データの取得は全 自動で行うことが可能である.被校正傾斜計は
Sensorex
タイプと
Q-Flex
タイプを想定しているが,傾斜計取り付け台座を介して取り付ける為,この台座を製作すれば どのタイプの傾斜計でも校正することが出来る.また,
ISO
に対応するために,全ての計測機器のトレーサビリ ティは確保されている.実際の校正手順は以下の通りとなる.
1.準備操作
・ 校正装置と傾斜計設置面を水平に調整
・ 傾斜計を設置し,計測系の配線を接続 2.信号設定
・ 取得するデータの収集条件を設定 3.試験条件設定
・ 校正を行うピッチ角,ロール角を設定 4.試験
・ 設定した条件で校正試験を実行 5.校正記録
・ 取得したデータを保存
JAXA 風洞技術開発センター傾斜計校正装置と内装式傾斜計の応用
青木 良尚、細江 信幸、星野 秀雄、伊藤 健(JAXA)
The Introduction to the JAXA Wind Tunnel Technology Center Inclinometer Calibration Facility and Application of the
Servo Inclinometer to Wind Tunnel Testing
Aoki yoshihisa, Hosoe nobuyuki, Hoshino hideo, Ito takeshi (JAXA)
Abstract
The Sensorex and Q-Flex type inclinometer calibration facility in JAXA Wintec was developed in August, 2003.
The Inclinometer Calibration Facility is introduced and 2-axis servo inclinometer calibration formula based on servo inclinometer model is derived in this paper. With an example for a wind tunnel testing using the Sensorex41400 servo inclinometer conducted in JAXA 6.5- by 5.5- m Low Speed Wind Tunnel, we show the application of the inclinometer to the wind tunnel testing.
This document is provided by JAXA.
第 78 回 風洞研究会議論文集 27
取得したデータをオフラインで処理して傾斜計校正係 数を算出することにより,傾斜計の校正が完了する.
3. サーボ傾斜計の 2 軸較正
3. 1
多項式フィッティングによる校正式と課題 ピッチ角のみを計測可能な 1 軸サーボ内装式傾斜計の 校正式は,ピッチ角をθ,ピッチ角出力電圧をS
θ,校正 係数をA,B
とすると,・・・(1)
となる.この式から 2 軸サーボ内装式傾斜計の出力もピ ッチ角とロール角の三角関数の線形和と仮定し,それぞ れの角度の範囲を± 90 °と仮定して 3 次多項式近似を行 うと,2 軸校正式は,ロール角をφ,ロール角電圧出力 を
S
φ,校正係数をa
i, b
iとすると,・・・(2)
となる.
この校正式を検証する為に,JAXA風洞技術開発セン ター傾斜計校正装置でサーボ傾斜計
Sensorex41400 の校
正データを取得した.校正データは図 2 に示すような 2 種類を取得して結果を比較した.校正データ 1 は図 2 の 黒丸で示されるピッチ角 0 °上の点とロール角 0 °上の 点,ピッチ角とロール角の絶対値が等しくなる点から 25 点抜粋した点であり,校正データ 2 は図 2 の×で示 されるピッチ角・ロール角共に 5 °間隔で取得した 225 点である.それぞれの校正データから校正係数を算出し,他の校正データを使用して角度を算出することによる検 証結果を表 2 に示す.
表 2 より,校正点が細かく校正範囲全体に分布してい る場合には,校正点以外の角度計測精度は十分確保され るが,校正点が荒い場合には校正点以外の計測精度が落 ちることが判る.また,式(2)で示される校正式の校正 項はピッチ角・ロール角合わせて 32 項あり,複雑であ る.したがって,式(2)の校正式には,計測精度を高め る為には多数の校正点が必要,校正項の数が 32 項とな り,1 軸校正式と比較して非常に複雑であるという課題 があることが判明した.
測定項目 項目 精度
ピッチ角 測定範囲 ± 40 °
設定精度 0.005 °以下 検出精度 0.002 °以下
ロール角 測定範囲 ± 180 °
設定精度 0.005 °以下 検出精度 0.002 °以下 高温槽 温度設定範囲 5 〜 85 ℃
温度制御精度 ± 0.3 ℃以内 校正条件 自動校正点数 最大 500 点 表
1
JAXA 風洞技術開発センター傾斜計校正装置の主な仕様図
1
JAXA 風洞技術開発センター傾斜計校正装置概観誤差平均 誤差標準偏差 誤差平均 誤差標準偏差
θ[deg] φ[deg]
校正データ 確認データ
2 1
1
2 -0.0008 0.0000
0.3717 0.0019
0.0038 0.0009
0.1031 0.0103
表2
多項式フィッティング 2 軸校正式検証結果図
2
多項式フィッティング 2 軸校正式検証データ-4 0 -3 0 -2 0 -1 0 0 10 20 30 40
-4 0 - 30 -2 0 - 10 0 10 20 30 40
φ [deg]
θ[deg] 校正データ1
校正データ2
This document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 28
3. 2
サーボ傾斜計のモデル化と校正式の導出前節の多項式フィッティングによる校正式の課題を解 決する為に,サーボ傾斜計のモデルに基づいて校正式を 導出する.
図 3 に示す 1 軸サーボ傾斜計モデルでは,傾斜計が傾 いても傾斜計に固定された座標系における重りの位置が 変わらないように重りに力を加える動作を行い,この重 りに加える力に比例した信号が出力される.従って,サ ーボ傾斜計モデルに基づく 2 軸校正式を導出するために は,2 軸を変角した時に重力に反して傾斜計内部の一定 の位置に重りを固定する為に必要な力を求めればよい.
まず,図 4 に示すように傾斜計に固定された座標系を 取る.ここで,x軸はロール回転軸,y軸はピッチ回転 軸と一致する方向とする.また,z軸正方向を重力加速 度方向と一致する位置をピッチ角・ロール角 0 °とし,
これらの角度に関する回転行列
R
θ,Rφを,重りの固定 位置のピッチ角・ロール角オフセット値をそれぞれ αi,βi(i=
θ,
φ)として下記の通りに定義する.・・・(3)
・・・(4)
ピ ッ チ 角 ・ ロ ー ル 角 方 向 の サ ー ボ 入 力 軸 を
(εθx
,
εθy,
εθz),(εφx,
εφy,
εφz),重力によって重りに 加わる荷重をf
g,オフセット電圧をΔθ,Δφ,サーボ入 力を電圧に変換する比例定数をk
θ,kφとすると,2 軸変 角による 2 軸信号出力S
θ,Sφは下式となる.・・・(5)
式(5)は,
という一次独立な 7 個の基底関数の線形結合である為,
サーボ傾斜計モデルに基づく 2 軸校正式はこれらの基底 関数から下式となる.
・・・(6)
以上より,サーボ傾斜計モデルに基づく 2 軸校正式(式
(6))が導かれた.傾斜計校正データから最小二乗法に
よって式(6)の 14 項の校正係数を求める事が出来る.ピ ッチ角・ロール角は,求めた校正係数と信号出力から,Newton-Raphson
法のような一般的な非線形方程式解法を用いて式(6)を解く事で計測することが可能となる.
図
4
傾斜計内部に固定された座標系の定義 図3
1 軸サーボ傾斜計モデルThis document is provided by JAXA.
第 78 回 風洞研究会議論文集 29
3. 3
サーボ傾斜計 Sensorex41400 によるサーボ傾斜計モデルに基づく校正式の検証
前節で求めたサーボ傾斜計の 2 軸校正式の妥当性を検 証する為に,多項式フィッティング校正式の検証に使用 した校正データに加えて,傾斜計を再設置して計測した 校正データ 1 と同様の 25 点校正データ 3,校正データ 1
〜校正データ 3 の 1 年前に計測したピッチ角が−10 °,0 °, 10 °におけるロール角−35 °〜 35 °まで 5 °間隔,ロール 角が−10 °,0 °,10 °におけるピッチ角−35 °〜 35 °まで 5 °間隔の 175 点の校正データ 4 を追加した.校正式の 検証は 25 点の校正データ 1 から計算した校正係数を使 用して,校正データ 1 〜校正データ 4 の計測値から角度 を計算して角度設定値と計測値の差分を評価して行っ た.表 3 にサーボ傾斜計モデルに基づく 2 軸校正式の検 証結果を,図 5 にピッチ角計測誤差のヒストグラムを,
図 6 にロール角計測誤算のヒストグラムを示す.
校正データ 1,校正データ 2 の誤差は全て± 0.1 °以内 となり,多項式フィッティングによる 2 軸校正式と比較 して,少ない校正点数で校正点以外の計測点でも良い計 測結果が得られた.再設置時データでは,ピッチ角・ロ ール角共に角度のオフセットが発生しているが,大きな 計測誤差は発生していない.再設置時データのロール角 誤差標準偏差は 0.02 °程度と他の誤差標準偏差と比較し て大きくなる.この原因は不明であるが,センサの特性 である可能性が考えられる.以上より,サーボ傾斜計モ デルに基づく校正式(式(6))は妥当である事が確認で きた.
4. 設置誤差補正法の検討
4. 1
傾斜計設置誤差先に導出したサーボ傾斜計モデルに基づく傾斜計の 2 軸校正式は妥当である事を確認したが,表 3 より,傾斜 計再設置後には校正データ 3 でピッチ角 0.013 °,ロー ル角で 0.035 °,校正データ 4 でピッチ角 0.017 °の微小 な計測誤差平均値のずれが発生している.これらの計測 誤差平均値のずれは傾斜計の誤差範囲の可能性もある が,傾斜計は傾斜計校正装置で校正された後に一旦取り 外して模型内部に再設置されるので,微小な設置誤差が 発生した場合に,設置誤差が傾斜計のピッチ角・ロール 角計測値に与える影響を補正する必要がある.次節で,
微小な設置誤差が傾斜計の計測値に与える影響の補正法 を検討する.
4. 2
微小な傾斜計設置誤差補正法の検討微小な設置角度のオフセットをオイラー角で定義し,
微小な傾斜計設置誤差を持つ傾斜計が計測するピッチ角 とロール角から,実際の設定ピッチ角とロール角を求め る補正式を算出する.
微小な設置角度のオフセットのピッチ角θe・ロール角 φe・ヨー角ψeを定義する回転行列をそれぞれ
R
θe,R
φe,R
ψeとすると,・・・(7)
・・・(8)
誤差平均 誤差標準偏差 誤差平均 誤差標準偏差
θ[deg] φ[deg]
校正データ1 校正データ2 校正データ3 校正データ4
0.0000 0.0004 0.0134 0.0174
0.0014 0.0022 0.0044 0.0030
0.0000 - 0.0005 0.0347 - 0.0038
0.0019 0.0024 0.0234 0.0192
0 10 20 30 40 50 60 70
-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05 ピッチ角校正誤差 [deg]
頻度 [%]
校正データ1 校正データ2 校正データ3 校正データ4
図
5
ピッチ角校正誤差ヒストグラム 表3
サーボ傾斜計モデルに基づく 2 軸校正式検証結果0 10 20 30 40 50 60 70
-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 φ較正値との差 [deg]
頻度 [%]
校正データ1 校正データ2 校正データ3 校正データ4
図
6
ロール角校正誤差ヒストグラムThis document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 30
・・・(9)
となる.また,実際の設定ピッチ角θと設定ロール角φ を表現する回転行列
R'
θ,R'φは,式(3)と式(4)におけるピ ッチ角・ロール角オフセットを 0 °とすると,・・・(10)
・・・(11)
となる.式(7)〜式(9)で表される設置誤差を含む傾斜計 ピッチ角θ^とロール角φ^の回転行列
R"
θ,R"φは,・・・(12)
・・・(13)
となる.以上の式(7)〜式(13)より,傾斜計で計測される ピッチ角・ロール角と,実際に設定されたピッチ角・ロ ール角の関係は,以下の式で表される.
・・・(14)
傾斜計で計測されたピッチ角・ロール角から実際に設定 されたピッチ角・ロール角を算出する為に,右辺の設置 角度のオフセット回転行列の逆行列を両辺に掛けて整理 すると,ピッチ角が± 90 °の範囲内では,以下の関係式 が得られる.
・・・(15)
式(15)を簡略化するために,設置角度のオフセット角 0 °周りでテーラー展開し,傾斜計で計測されたピッチ 角・ロール角から実際に設定されたピッチ角・ロール角 の一次近似式を求めると,下式となる.
・・・(16)
式(16)の補正式から,供試体の基準面上で計測した姿勢 角 3 点以上から,設置角度のオフセット角を求める事が 出来る.
式(16)の補正式を使ってオフセット補正をした補正後 のピッチ角校正誤差を表 4 に,補正後のロール角校正誤 差を表 5 に,補正後のピッチ角校正誤差ヒストグラムを 図 7 に,補正後のロール角校正誤差ヒストグラムを図 8 に示す.
補正後のピッチ角では,誤差の平均値は 0 °に近づき,
誤差の標準偏差は校正データ 3 では半分程度に,校正デ ータ 4 では 1.4 倍となった.校正データ 3 は校正に使用 した校正データ 1 を取得直後に再設置して計測したデー タなので,傾斜計自体の特性の変化の影響が小さいため に設置角度補正の効果があったと考えられる.
補正後のロール角では,誤差の平均値は 0 °に近づく が,誤差の標準偏差はほとんど変化がない.したがって,
校正データ 3 と校正データ 4 で見られるロール角の誤差 の標準偏差が他の誤差の標準偏差と比較して大きいの は,設置角度の影響ではなく,センサの特性など他に原 因があると考えられる.
以上より,導出した設置角度補正近似式(式(16))は 誤差のオフセットを補正するのに妥当であることがわか った.
This document is provided by JAXA.
第 78 回 風洞研究会議論文集 31
5. JAXA6.5m × 5.5m 低速風洞における内装式 傾斜計の応用
5. 1
内装式傾斜計の使用法J A X A6.
5m× 5.
5m低 速 風 洞 で は , サ ー ボ 傾 斜 計Sensorex41400 を供試体の内部に設置面を製作して搭載
し,ピッチ角の計測に使用している.この風洞の支持装 置はピッチ&ヨー変角機構であるので,内装式傾斜計の出力で使用するのはピッチ角のみとしている.供試体内 部でピッチ角を計測するので,たわみを含めたピッチ角 を計測することが出来る為,たわみ補正は行わない.
5. 2
ピッチ角の傾斜計出力とエンコーダ出力の比較JAXA6.5m
× 5.5m低速風洞で行われたONERA M
5 2.5 倍相似模型による標準模型試験で取得したデータに より,ピッチ角の傾斜計出力とエンコーダ出力を比較し た結果を図 9 に示す.たわみ補正を行っていないが,ピ ッチ角の傾斜計とエンコーダ出力はほぼ 0.05 °範囲で一 致している.傾斜計はアナログ出力である為,エンコー ダよりも高精度・高分解能でピッチ角を計測可能とな る.6. まとめ
・
JAXA
風洞技術開発センター傾斜計校正装置の概 要を示した.・ サーボ式傾斜計のモデルを作成し,校正データに よってその妥当性を示した.
・ サーボ式傾斜計
Sensorex41400 の校正を行い,導出
した校正式が妥当であることを確認した.・ 校正条件とは異なった取り付け角を持った場合の 補正方法を検討し,妥当性を確認した.
・
JAXA6.5m
× 5.5m低速風洞における内装式傾斜計の応用を紹介した.
7. 謝辞
標準模型試験を行う為に,JAXA6.5m× 5.5m低速風 洞の皆様には大きな助力を頂きました.また,超音速風 洞の永井伸治セクションリーダーには,適切な助言を頂 きました.ここに感謝の意を表します.
0 10 20 30 40 50 60 70
-0.2 -0.15 -0.1 -0.05 0 0.05 0.1 0.15 0.2 傾斜計とエンコーダピッチ角差分 [deg]
頻度 [%]
傾斜計とエンコーダ出力差分
図
9
標準模型試験におけるピッチ角の傾斜計出力とエンコー ダ出力比較結果誤差平均 誤差標準偏差 誤差平均 誤差標準偏差 θe ψe φe
補正前φ[deg] 補正後φ[deg]
− 校正データ1
校正データ2 校正データ3 校正データ4
0.0000 -0.0005 0.0347 -0.0038
0.0019 0.0024 0.0234 0.0192
0.0000 0.0000 0.0000 0.0000
0.0019 0.0022 0.0232 0.0186
0.0003 0.0138 0.0179
-0.0016 -0.0125 0.0100
-0.0005 0.0347 -0.0038
表
5
補正後ロール角校正誤差誤差平均 誤差標準偏差 誤差平均 誤差標準偏差 θe ψe φe
補正前θ[deg] 補正後θ[deg] オフセット角[deg]
校正データ1 − 校正データ2 校正データ3 校正データ4
0.0000 0.0004 0.0134 0.0174
0.0014 0.0022 0.0044 0.0030
0.0000 0.0001 0.0001 0.0000
0.0014 0.0022 0.0022 0.0044
0.0003 0.0138 0.0179
-0.0016 -0.0125 0.0100
-0.0005 0.0347 -0.0038
表
4
補正後ピッチ角校正誤差0 10 20 30 40 50 60 70
-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05 補正後ピッチ角校正誤差 [deg]
頻度 [%]
補正後校正データ2 補正後校正データ3 補正後校正データ4
図
7
補正後ピッチ角校正誤差ヒストグラム0 10 20 30 40 50 60 70
-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05
補正後ロール角校正誤差 [deg]
頻度 [%]
補正後校正データ2 補正後校正データ3 補正後校正データ4
図
8
補正後ロール角校正誤差ヒストグラムThis document is provided by JAXA.
宇宙航空研究開発機構特別資料 JAXA-SP-07-005 32
参考文献
[1] T. Finley and P. Tcheng, “Model Attitude Measurements at NASA Langley Research Center”, AIAA 92-0763 [2]
川崎重工業株式会社, “模型姿勢角センサ・角度較正装置 完成図書”,2003
This document is provided by JAXA.
第 78 回 風洞研究会議論文集 33
1. はじめに
JAXA0.5m
極超音速風洞(以下,HWT1)は 1965 年に完成し,ノズル出口直径が 0.5m,通風時間が最大 120 秒,ノズルを交換することによりマッハ数 5,7 及び 9 の試験が可能である[1].1995 年には 1.27m極超音速風 洞(以下,HWT2)が既設の
HWT1 と設備を一部共用す
る形で増設された.HWT2 はノズル出口直径が 1.27m,通風時間が最大 60 秒,固定ノズルによるマッハ数 10 の 世界最大規模の極超音速風洞である[2].HWT1/HWT2 の 全体図を図 1.1 に,
HWT1 を図 1.
2 にそれぞれ示す.HWT1 のマッハ数 5 及び 7 ノズルについては近年,それ
ぞれ詳細なマッハ数校正試験が行われた[3].気流校正試験結果のみからでは得られない空力特性へ の影響は,標準的な単純形態の模型を用いた標準模型試 験 に よ り 詳 し く 調 べ る こ と が 重 要 で あ る .
H W T1 ,
HWT2 ではこれまで HB-1,HB-2,AGARD-E
等の標準模型を使い風洞の検証が行われてきたが,有翼形状の標 準模型試験は行われていない.近年,有翼形態の極超音 速飛行体に関する研究開発が予想されるため,有翼形態 における空力特性への影響を確実に把握する必要性が高 まってきた.AGARD-B模型は有翼形状の標準模型のひ
とつとして,古くから多くの国内外の超音速及び極超音 速風洞で試験が行われている.今回,HWT1 において
AGARD-B
標準模型によるマッハ数 5 の 6 分力試験を行い,空力特性データを取得した.本試験結果を他風洞の 実験結果及び異なるマッハ数での結果と比較し,検討す る.
2. 風洞試験
2.1
AGARD-B 標準模型とベース圧配管本試験に用いた
AGARD-B
模型の諸元を表 2.1,模型 概略を図 2.1 にそれぞれ示す.測定室内の模型支持装置(スティング)に取り付けた同模型を図 2.2 に示す.本 模型の寸法は後部胴体直径(48mm)を基準に定められ ている.両翼は後退角 60 度のデルタ翼,翼スパンは 192mm,全長は 408mmである.模型材質は
SUS304 で
ある.天秤の取り付け部には断熱材としてジルコニアを 使用し,天秤温度ドリフトを防ぐよう配慮した.両翼は 取り外しができるように設計し,翼無し形態の試験も行 える.ベ ー ス 圧 力 を 測 定 す る 圧 力 配 管 (SUS管 内 径 : 2.0mm,長さ:約 35cm)を模型後端部の上下対称位置 に配し,スティングの貫通孔を通して圧力センサと接続
0.5m 極超音速風洞におけるマッハ 5 AGARD-B 標準模型 6 分力試験
津田 尚一、小山 忠勇、平林 則明、渡利 實(JAXA)、関根 英夫(JAST)、 木伏 淳子(㈱スペースサービス)、中村 晃祥(JAST)