国
立
大
学
法
人
室
蘭
工
業
大
学
航
空
宇
宙
機
シ
ス
テ
ム
研
究
セ
ン
タ
ー
年
次
報
告
書
２
０
１
３

Muroran Institute of Technology
Aerospace Plane Research Center

Annual Report 2013

年次報告書２０１３

２０１４年８月
国立大学法人 室蘭工業大学

航空宇宙機システム研究センター

巻頭言
超音速飛行に向けた基盤研究促進と活溌な共同研究の実施
センター長 東野和幸
平成２６年度より，棚次亘弘前センター長から本職務を引き継ぎしましたのでよろしくお願いします．
平成２５年度は特別経費（プロジェクト分）から一般経費へ組替えられた２年目の年度です．一定期間（最低３年間）は，当初計画どおり事業が進捗しているかを把握するため，通常の特別経費と同様に，事業の進捗状況の報告が求められています．
本学は研究活動の更なる向上を図るため，研究の実施体制や研究成果等について平成２４年
度に自己点検・評価を行い，さらに学外有識者からの評価を受けました．評価の概要は年次報告書に記載しました．航空宇宙システム工学分野は本学のミッションの再定義にあげられ，さらに平成２７年度までの大学の第２期中期目標期間においても重点研究分野になっています．
航空宇宙は高度なシステム工学の象徴であり，主要な構成要素である機体，エンジン，誘導制
御（データ伝送を含む）そして飛行力学の間でシステム整合性を図る必要があります．また，この高度なシステムを安全に効率よく実験するための運用や関係する法規についても検討中です．
機体については，本格的な超音速飛行試験が可能な「オオワシ２」の概念設計を行い，実物大のモックアップを製作しました．これを用いて，搭載機器の配置やメンテナンス性等の検討し，離着陸時の荷重や衝撃に対応する脚についても検討をしています．
超音速飛行を行うためのエンジンとして，小型で大推力を発生するエアーターボラムジェットエン
ジン（GG-ATR）のコンポーネント設計製造を進め，ファンやタービン等いわゆる回転系について部
品製造を完了し，動バランスをとり，組立を完了しています．また，ラム燃焼器やＧＧの高温部分の設計検討を継続中です．
２０１１年度に実施しました「オオワシ１」の飛行試験の結果，低速飛行時の操縦の難しさを再認識しましたので，オンボードコンピューターによる全自動操縦を行うための誘導制御の研究を進行中です．飛行力学の観点からは飛行に必要な空力制御について風洞試験や解析により制御能力を高める工夫を進めています．
大型試験設備の高速走行軌道試験設備については，民間企業との共同研究で本格的な運用段階に入り，航空宇宙機に搭載する機器の高耐 G 試験や高速空気力学実験ができます．また，他大学とは共同研究で先端エンジンによるサブサイズ軌道による走行実験も継続しており推進系の立証装置としての役割が増しています．このことは飛行試験を行わないで，地上で繰り返し，安全に試験ができ，開発コストの低減や開発期間の短縮に繋がります．超音速風洞においてもインテークの基礎実験をＪＡＸＡ殿や大学と共同研究を実施しました．今後，需要がさらに増加する見込みです．

以上のようなプロジェクト研究には，クリーン推進燃料，飛行制御，空力制御，構造解析等に関する基盤技術の研究も含めて進めています．推進燃料に関する研究では，アルミニウムと水を触媒を介して高圧水素をを発生させられることを実験によって確認し，小型衛星用推進装置への適用の観点からＪＡＸＡ殿との共同研究を進めています．また，学外設備の白老エンジン実験場にて民間企業と新規ロケットエンジン基礎燃焼実験等の共同研究を継続して進めています．これら基盤技術研究の詳細についても本報告書の各項を参照してください．本研究センターの研究開発の進捗状況や組織および試験設備等の詳細については，本学ホームページの「航空宇宙機システ
ム研究センター」の項を参照ください． （http://www.muroran-it.ac.jp/aprec/）目次
巻頭言～超音速飛行に向けた基盤研究促進と活溌な共同研究の実施
連携および共同研究 ............................................................................................................... 1
講演等の啓蒙活動の概要および見学者 .................................................................................. 3

研究概要[機体関連]
第二世代小型超音速飛行実験機の操舵空力 ........................................................................... 5
第二世代小型超音速飛行実験機のエリアルールに基づく抗力低減 .................................... 10
小型超音速飛行実験機の飛行性能の予測 ............................................................................ 14
小型超音速飛行実験機の予備的飛行試験のための縮小機体の設計・製作………………….19
小型超音速飛行実験機のエアインテークの空力特性 .......................................................... 24
小翼列によるアクティブ制御を搭載したラジコン機風洞試験 ........................................... 29
小翼列によるアクティブ制御を搭載したラジコン機実証飛行試験 .................................... 33
層流後退翼の境界層遷移制御に関する風洞実験 ................................................................. 37
超音速風洞における小型トラバース装置 ............................................................................ 41
小型無人超音速機の応力解析と剛性解析 ............................................................................ 43
小型無人超音速機の翼胴一体型構造による機器搭載性向上対策と１／２スケールモックアップの製作 ............................................................................................................................ 48[エンジン関連]
GG-ATR エンジンの回転系試験とガスジェネレータ・ラム燃焼器の燃焼、冷却特性について .......................................................................................................................................... 53
GG-ATR エンジン用超音速インテークの基礎研究 ............................................................. 56
アルミー水反応の衛星推進系への適用 ................................................................................ 58
バイオエタノールの熱分解吸熱反応と触媒効果について ................................................... 63
高圧，高温環境下におけるバイオエタノールの熱分解吸熱反応，サルファアタック・コー
キングの実験装置について .................................................................................................. 70
反転軸流ファン試験装置 ...................................................................................................... 73

[誘導制御関連]
小型無人航空機向け遠隔監視制御系用無線通信装置の構築………………………………….78
小型無人超音速機向け誘導制御システムの構築とラジコン機による飛行検証 .................. 81
小型無人超音速機の縦系着陸制御系設計と外乱に対する性能評価 .................................... 85
無人航空機向けシステム同定の研究（-連続系伝達関数の開ループ同定-） ....................... 90[設備関連]
高速走行軌道実験設備に関する基盤技術実証 ..................................................................... 93

発表論文 ........................................................................................................................... 100

連携および共同研究○東野 和幸（航空宇宙機システム研究センター 教授）
中田 大将（航空宇宙機システム研究センター 特任助教）

１ ＪＡＸＡ宇宙輸送系システム技術研究開発センターとの共同研究：「将来輸送系リファレ
ンスシステムの推進系に関する研究」
平成２４年度に引き続き，バイオエタノール/液体酸素を用いた推力 21ton 級，250ton 級の各ロケットエンジンシステムについて成立性，性能比較さらに技術課題についての検討を実施した．昨年度の結果を踏まえ，GG サイクルに焦点を絞りタービン効率・ポンプ効率が与える影響について精査した．また，タンクについての技術課題についても列挙した．
２ ＪＡＸＡ研究開発本部推進系グループとの共同研究：「衛星推進系信頼性向上のための基盤研究（その３）」
アルミニウム粉末と水との反応で水素が生成されるプロセスを衛星推進系に利用するための基礎研究を実施した．水素生成量に及ぼすアルミ粉末粒径の影響や，生成率を高めるようなアルミニウム合金について検討した．
３ 川崎重工業株式会社との共同研究：「高加速度環境下における高速走行軌道実験の発展」
白老エンジン実験場に設置したフルサイズ高速軌道試験設備を用いて高加速度環境の実験を実施した．川崎重工業株式会社は，高加速度環境下において航空宇宙機器の健全性を確認した．室蘭工業大学は高速走行する走行台車の走行プロファイルと水制動の特性を取得した．走行中の最大加減速度は約 13 Ｇ，最高速度は約 350 ｋｍ/h であった．
４ 名古屋大学との共同研究：「高速走行軌道装置における Rotating Detonation Engine 滑走試験」
白老エンジン実験場に設置したサブスケール高速軌道試験設備を用い名古屋大学が研究中の Rotating Detonation Engine の滑走試験を実施した．目的は滑走により推力の発生とその規模を確認することである．実験は 2013 年 12 月と 2014 年 3 月に行われたが，点火系の不具合等により滑走には至らなかった．
5 早稲田大学・JAXA との共同研究:「エアインテーク性能取得風洞試験」
早稲田大学・JAXA と共同で超音速インテークの風洞試験を実施した．JAXA では小型予冷
ターボジェットエンジンの飛行試験を計画（kBOV 計画）しており，今回風洞試験を実施したのは，同飛行試験で用いられるフライトモデル用インテークである．この風洞試験では，インテーク出口のオリフィスを変えることで流量の調節を行った他，インテークのスロート高さ

1

も変化させて，全圧回復率と流量捕獲率の計測を行った．同時に室蘭工大で検討されているオオワシ 2 号機用超音速インテークの，予備的風洞試験も実施して，インテーク設計に必要なデータを計測した．

図１：インテーク風洞試験にて得られたシュリーレン写真
６ 大阪府立大学との共同研究：「小型超音速飛行実験機の空力特性の改良と評価」
小型超音速飛行実験機の空力形状を改良し，亜音速ならびに超音速風洞試験によって空力特性を評価した．操舵および姿勢変化による空力特性の変化や，推進器ながれと外部ながれの干渉による空力特性の変化を評価した．
７ 東京都市大学との共同研究：「大学における教育用小型ロケットに関する研究」
東京都市大の液体窒素・水ロケットは燃焼を伴わない安全なロケットであり，高校生などの教育用に適している．H23 年度に白老エンジン実験場において水平走行実験を実施し，H24 年度には打ち上げシステムの検討を行った．
室蘭工大ではハイブリッドロケットと呼ばれる亜酸化窒素＝プラスチック燃料のロケットを運用しているが，このシステムと比べた時の運用性・コストなどについて整理すると共に液体窒素・水ロケットの打ち上げ実験を北海道で行う際の予備検討など
も行った．

2

講演等の啓蒙活動の概要および見学者○ 東野 和幸（航空宇宙機システム研究センター 教授）
中田 大将（航空宇宙機システム研究センター 特任助教）

航空宇宙機システム研究センターの専任教員および併任教員は，依頼講演，高校訪問，オー
プンキャンパス，大学開放推進事業に協力しています．平成２５年度に行いましたこれらの活動の概要を表１にまとめました．ここでは主に社会に対する啓蒙や本学の広報に関する活動を示しました．学会やシンポジウム等の研究に関する活動は除いています．
表１ 航空宇宙機システム研究センターが協力した啓蒙・広報活動の概要
団体

開催名

日本航空

第 55 回構造

宇宙学会

強度に関す

構造部門

る講演会

北海道職

第１１回北海業能力開道ポリテック

発大学校

ビジョン

日時

場所

平成 25 年 8 月 8 日

室蘭工業

15:20～15：50

講演者

企画講演
棚次亘弘

大学

平成 26 年 2 月 21 日
10:30～12:00

講演名「大学でできる航空宇宙
分野の実践研究」

北海道職

記念講演

業能力開

棚次亘弘

発大学校

「宇宙旅行に出かけるためのものづくり」

航空宇宙機システム研究センターには，報道機関の取材，国外の大学関係者，中学・高校の教諭が見学されます．見学の対象は主に超音速風洞設備，オオワシ２号機モックアップ，反転ファン
試験設備，フライトシミュレーター，高速走行軌道実験設備、白老エンジン実験場です．平成２５年度に訪問された学外の見学者を表２に示します．
表２ 航空宇宙機システム研究センターを訪問された見学者

見学の目的

日時

参加人数

研究センターの活動説明、風洞装置とシミュレーター見学

平成 25 年 4 月 16 日（火）アルミ／水反応に関する情報

14：00～18：30

東京大学名誉教授：荒川義博
3

洞装置とシミュレーター見学
IHI の新規事業について

平成 25 年 5 月 30 日（木）13：30～17：00
研究センターの活動説明、風

平成 25 年 6 月 26 日（水）洞装置等 所内設備の見学

13：50～14：30

3静岡大学 助教：松井 信
日本エクスクロン 技術顧問：白形弘文

交換研究センターの活動説明、風参加者（敬称略）

1

17（株）IHI 理事、新事業推進部長
山本 洋一

昭和４５年電気工学科卒業生

研究センターの活動説明、風

平成 25 年 6 月 26 日（水）洞装置等 所内設備の見学

15：30～16：30

研究センターの活動説明、風

平成 25 年 7 月 5 日（金）洞装置等 所内設備の見学

12：40～13：10

研究センターの活動説明、風

平成 25 年 8 月 9 日

洞装置等 所内設備の見学

15：20～16:30

研究センターの活動説明、白平成 25 年 9 月 6 日

老実験場案内10：30-12：00

研究センターの活動説明、風

平成 25 年 9 月 7 日

洞装置等 所内設備の見学

13：30～14：30

研究センターの活動説明、風

平成 25 年 10 月 15 日

洞装置等 所内設備の見学

10：00～12：00

5

風洞装置、オオワシモックア

平成 25 年 10 月 29 日
11：00～12:30

学科教員 1 名，国際交流センター3 名
文部科学省科学技術・学術政策局長
1

土屋定之
学長，事務局長，総務課長

10

JSASS 構造強度に関する講演会参加者

16

ISAS ハイブリッドロケットＷＧ有志

2

三菱重工業 課長 恩河忠興、金子敬郎
独立行政法人高齢・障害・求職者雇用支援

1

機構、北海道職業能力開発大学校
校長 前田康二

GG-ATR エンジンのファンとタービン

フィンランド関係者 5 名，情報電子工学系

独立行政法人高齢・障害・求職者雇用支援
2

機構、北海道職業能力開発大学校
山口力、川畑雅司

ップ、シミュレーター見学
研究センターの活動説明、風

平成 25 年 10 月 28 日

洞装置等 所内設備の見学

14：35～14：50

研究センターの活動説明、風平成 26 年 3 月 14 日

洞装置等 所内設備の見学

14：50～16：00

48

児玉隆夫、土屋 俊、西永、渡邊一衛、学位授与機構職員 4 名
文部科学省官房参事官 佐伯 浩治
2

文部科学省官房人事課 総務班 総務係長
熊谷 達也

第二世代小型超音速飛行実験機の操舵空力○鈴木

祥弘 （航空宇宙システム工学専攻 博士前期 1 年）

上村

隆太 （機械航空創造系学科 4 年）

溝端

一秀 （もの創造系領域 准教授）

１．はじめに
第二世代小型超音速飛行実験機（オオワシ）の M2011 空力形状について，エレベータ操舵によるピッチトリム性能，エルロン操舵によるロールモーメント，およびエルロンのみによるロール
制御性能劣化すなわち AAD（Aileron Alone Departure）が，昨年度の亜音速風洞試験で評価されている[1]．そこで本年度は，引き続き，遷音速域でのピッチトリム性能，亜音速域でのラダー
性能，翼面スポイラーによるロール制御性能，エルロンとラダーの併用によるロール制御性能，および外翼フラッペロンのフラップとしての性能を，風洞試験によって評価する．
２．小型超音速飛行実験機の舵面設定
第一世代および第二世代小型超音速飛行実験機に当面設定されている蛇面を図１および表１に示す．第二世代機体に固有の舵面設定について以下に述べる．
２－１．内翼フラップの廃止と外翼フラッペロンの追加
第二世代の空力形状 M2011 においては，胴体と主翼の結合部強化と Area rule に基づく抗力低減のために翼胴フィレットの搭載が予定されている．第一世代の空力形状 M2006-prototype に採用された内翼フラップは，この翼胴フィレットと空間的に干渉する可能性がある．そこで，当面，
空力形状 M2011 では内翼フラップを廃し，代わりの高揚力装置として外翼エルロンにフラップ機能を追加した外翼フラッペロンを採用している．
２－２．内翼スポイラーの追加
フラッペロンをフラップモードで操舵している場合，さらにエルロンとして操舵するのは，機構上の舵角の限界がある．そこで，ロール制御の補助として内翼スポイラーを採用する．これは，主翼片翼上面のスポイラーを立ち上げることによってその片翼の揚力を減じ，ローリングモーメ
ントを発生させるものである．同時にその片翼の抗力が増すため，アドバースヨーの抑制，ひいてはロール制御性能劣化（LCD：Lateral Control Departure）の緩和が期待される．また，両翼のスポイラーを立てることによってエアブレーキとして機能させることも期待される．

(a) 第一世代（M2006prototype 形状）(b) 第二世代（M2011 形状）
図１．小型超音速飛行実験機の舵面

5

表１．小型超音速飛行実験機の舵面
Function

M2006prototype

Roll Control

Outboard Aileron,
ElevonPitch Control
Yaw Control
High Lift Device
Air Brake

ElevonRudderInboard Flap
NA

M2011
Outboard Flaperon,Inboard Spoiler，ElevonElevonRudderOntboard Flaperon
Inboard Spoiler

３．風洞試験
３－１．風洞試験装置
風洞試験装置として，JAXA/ISAS の吹出し式遷音速風洞および大阪府立大の回流式亜音速風洞を用いる．通風条件は，遷音速風洞でマッハ数 0.5~1.3, 迎角範囲±12°，亜音速風洞で流速約30m/sec，迎角範囲±20°ないしは±30°である. 空力測定には 6 分力内挿天秤を用いる．
３－２．風洞試験模型
２つの空力形状 M2006prototype および M2011 は，主翼・尾翼の形状および位置関係，ならびに翼幅と胴体径の関係は相似である．従って，風試模型のうち主翼部品，尾翼部品，および中胴・後胴部品は共用できる．M2006prototype と M2011 の寸法比が 1.5 倍であることから，風試模型の縮小比は M2006prototype に対して 7/40 倍，M2011 に対して 7/60 倍としている．
舵面については，図２の通り外翼エルロン付き主翼部品（エルロン舵角 10 度，20 度）およびエレボン部品（エレベータ舵角 10 度，15 度，エルロン舵角 10 度，15 度）を昨年度製作済みである[1]．さらに，図３の通り外翼フラップ付き主翼部品（フラップ舵角 20 度・30 度），ラダー付き垂直尾翼（ラダー舵角 10 度・20 度），スポイラーマウント付き主翼部品，およびスポイラー
部品（舵角 0，10 度，20 度，30 度）を今年度製作した．

(a) 外翼エルロン付き主翼部品(b) エレボン部品

図２．昨年度製作された風試模型部品(a) 外翼フラップ付き主翼部品(b) ラダー付き垂直尾翼部品

6

(c) スポイラーマウント付き主翼部品とスポイラー部品
図３．今年度製作された風試模型部品
３－３．風試結果と考察
天秤計測上の空力モーメントの回転中心は天秤中心であるが，機体重心周りに換算して空力モーメントを評価する．第二世代実験機の重心位置は構造設計・艤装計画に依存してやや流動的であるため，ここでは，第一世代実験機と同じく重心位置を平均空力翼弦の 25%に仮定する．
３－３－１．遷音速域でのピッチトリム性能
空力形状 M2011 について，種々のマッハ数におけるエレベータ操舵によるピッチングモーメン
ト係数を図４に示す．0.5～1.3 のマッハ数全域に渡ってピッチング静安定を維持しており，エレ
ベータ舵角±10°の範囲では舵角にほぼ比例してピッチングモーメントが増減している．一方，舵角が 10°を超えると，それ以上舵角を増やしてもモーメント増分が得られなくなる．また，
M0.5 の場合に比べて M0.9～1.3 では，エレベータ操舵によるピッチングモーメント増分が半分程度に減じており，これに対応して，ピッチング静安定を取ることのできる迎角範囲も， M0.5
では-15°~12°であったものが M1.3 では-10°~5°となっている．オオワシの高亜音速～超音速飛行において所要迎角は 1°程度以下であることから，このエレベータ性能によって十分にピッ
チトリムを維持できると判定される．また，これらピッチトリム性能は重心位置によって変わる．
３－３－２．亜音速域でのラダー性能
M2011 形状について，マッハ 0.5 におけるラダー操舵によるヨーイングモーメント係数を図５に示す．ラダー操舵によるヨーイングモーメント生成は良好と判定できる．

図４．遷音速域でのピッチングモーメント 図５．ラダー操舵によるヨーイングモーメント
３－３－３．低亜音速域および遷音速域でのエルロンおよびスポイラー性能
種々のマッハ数における外翼エルロン操舵及びスポイラー操舵によるローリングモーメント係数を図６に示す．全てのマッハ数でエルロン操舵によって十分なローリングモーメント増分が得られており，エルロン性能は良好であると判断できる．マッハ数が大きいほどローリングモーメ
ント係数が小さくなる傾向が見られるが，オオワシの飛行経路においては飛行マッハ数が大きい
7

ほど飛行動圧が高くなるため，ロール制御力に問題は無い．
一方，スポイラーに関しては，エルロンと比較してモーメント増分が極端に小さいことがわかる．今後スポイラーの位置および面積を変えて，さらに詳細な評価・検証を行う必要がある．

(a) 低亜音速(b) 遷音速

図６．外翼エルロンおよび翼面スポイラーによるローリングモーメント
３－３－４．横制御発散
小型超音速飛行実験機オオワシのように大きな後退角および高翼形態を有する機体は上反角効
果が大きく，離着陸時などの大迎角条件においてアドバースヨーの増加および方向安定の劣化と相俟 って，エル ロン操舵に よる ロール 制御性が劣 化 し，更に は逆転す る Lateral Control
Departure 現象が発生する可能性がある．その発生可能性を評価するために，以下の横制御発散係数(Lateral Control Departure Parameter : LCDP)が用いられる：
LCDP = 𝑪𝒏𝜷 −ここで，𝑘 =

𝑪𝒏𝜹 +𝒌𝑪𝒏𝜹
𝒂𝒓𝑪𝒍𝜹 +𝒌𝑪𝒍𝜹
𝒂

𝒓

𝑪𝒍𝜷

（1）．

𝛿𝑟
⁄𝛿 であり，ラダーゲインと呼ばれる．LCDP が負となるときにエルロンの効き
𝑎が逆転し，機体は操舵方向とは逆方向にロールする（ロール・リバーサル）．この現象は，ラダー操舵によるヨーイングモーメントによって緩和される．ラダー操舵を伴わない場合（ラダーゲイ
ンゼロ）の LCDP を特に区別して Aileron Alone Departure Parameter (AADP)と呼ぶ．LCDPと AADP の値を比較することによって，ロール制御に対するラダーの効果を評価できる．ラダー
ゲイン k を 1.0 としたときの LCD および AADP を図７に示す．M2006prototype 形状に比べてM2011 形状は AADP が正となる迎角範囲が狭いことがわかる．その迎角範囲は，ラダーを操舵無しの場合-10°~12°であり，ラダー操舵有りの場合は-10°~17°へと拡大している．
３－３－５．亜音速域での外翼フラッペロン性能
M2011 形状について，亜音速域において外翼フラッペロンをフラップとして操舵することによる揚力係数増分を図８に示す．離陸時の迎角を 10 度とし，全備質量 350，400，450kg，フラップ舵角 0 度，30 度に対する離陸速度を推算すると表２の通りとなる．フラップ舵角 30 度の揚力増分によって離陸速度が 10 [m/sec]程度（14%程度）低減される．

8

図７．亜音速域における LCDP

図８．外翼フラップによる揚力増分

表２．離陸速度の推算
Flap deflection
[deg]
30

0

Take-off mass
[kg]
350
400
450
350
400
450

Take-off airspeed[m/sec]
64.8
69.2
73.4
75.0
80.2
85.1

Lift coefficient

AOA
[deg]

0.59
10
0.44

４．まとめ
小型超音速飛行実験機の操舵空力特性の評価のために風洞試験を行ったところ，以下の結果を得た．
１） エレベータの効きは良好であり，重心位置が 25%MAC の場合，ピッチングトリムを取れる
迎角範囲は低亜音速では±10°，M0.5 では-15°~12°，M1.3 では-10°~5°である．
２）エルロンの効きは良好であるが，それに比較して内翼スポイラーの効果は非常に小さい．ス
ポイラーの搭載位置や面積の改良が必要である．
３）ラダーの効きは良好であり，ロール・リバーサルをラダー操舵によって緩和できる．
４）外翼フラッペロンをフラップとして操舵することにより，離陸速度を 10 [m/sec]（14%）
程度低減できる．
今後，引き続き，スポイラーの改良と空力特性評価を進めると共に，エレボンによるロール制
御性能，内翼フラップによる離陸速度低減効果，等を調べる計画である．

9

第二世代小型超音速飛行実験機のエリアルールに基づく抗力低減○大石

栄 （航空宇宙システム工学専攻 博士前期 2 年）

行澤

陽 （機械航空創造系学科 4 年）

溝端

一秀 （もの創造系領域 准教授）

１．はじめに
第二世代小型超音速飛行実験機の空力形状として，別途開発中のガスジェネレータサイクル・
エアターボラムジェット（GG-ATR）エンジン１基を搭載する M2011 形状が提案されている[1]．
これまでの風洞試験結果と飛行経路解析の結果，GG-ATR エンジンと M2011 空力形状の組合せ
ではマッハ数 1.1～1.6 の遷音速域で推力余裕が厳しくなり，目標とするマッハ数 2.0 まで加速させるには 8%程度の抗力低減が必要であることが判明した[2]．そこで，本研究では遷音速・超音
速域特有の造波抗力に着目し，機体の抗力低減手法の一つであるエリアルール[3]に基づいてM2011 空力形状を改良する．
２．エリアルールを適用した空力形状
エリアルールは，機軸上の任意点から発生するマッハコーンで機体を切断し，その円錐断面の
底 面 積 の 機 軸 方 向 分 布 が 滑 ら か に な る よ う に 機 体 形 状 を 修 正 す る 方 法 で あ る [4] ． 元 のM2011Nose-C 空力形状の風試模型と，それをエリアルールに従って修正した結果の断面積分布を図１に示す．

(a) 元の M2011Nose-C 空力形状．

(b) エリアルールによって修正された空力形状．

図１．元の M2011Nose-C 空力形状とエリアルール適用形状の断面積分布．

10

AR-NoseC

BulgeA

FilletBulgeB

図２．エリアルール適用形状におけるノーズ，フィレット，およびバルジの風試模型部品．
ここでは，主翼・尾翼の断面積が造波抗力最小となる Sears-Haack 曲線に極力収まるように，
ノーズコーンの底面を後退させ（エリアルールノーズ C），主翼・尾翼は相対的位置関係を変えずにまとめて前方へ移動させている．また，主翼尾翼間のくびれをなくすために胴体にバルジ（膨らみ）を設けている．フィレットは翼胴結合部での空気流を整えるために設ける．それぞれの模型部品の三面図および立体図を図２に示す．
３．風洞試験と造波抗力解析
エリアルールによる形状修正の効果を確かめるために，JAXA/ISAS の遷音速風洞を用いて風洞試験を実施した．併せて，NASA で開発された造波抗力計算プログラム WAVEDRAG を使用
し造波抗力を推算した．
造波抗力解析における模型形態条件は表１の 15 通りである．JAXA/ISASにおける遷音速風試の模型形態条件は，各種制約から Case No. 5 のみである．
Case No. 5 について風洞試験と造波抗力解析の結果を図３に示す．風試結果によれば，マッハ1.08～1.32 の全域でエリアルール適用の効果が確認でき，5%～20%程度抗力が低減されている．
造波抗力解析は細長物体理論すなわち非粘性微小擾乱ポテンシャル理論に基づいており，そこには粘性抗力が勘定されないため，その推算値は一般的には風試による抗力値と一致しない．ここでは，マッハ数依存性の概略を比較するために，グラフ上で造波抗力解析結果を上方へ一定量平行移動している．その結果，マッハ数依存性は，風試結果と造波抗力解析結果で，概ねよく一致していることが確認できる．
次に各模型形態について造波抗力解析結果を図４に示す．ノーズコーン底面を後退させることでマッハ 1.0～2.0 の全域で造波抗力低減が可能と推定される．また，バルジ B よりバルジ Aのほうが抗力低減効果は高いと推算されており，バルジの効きには限度があるものと推定される．
更に，マッハ数 1.0～1.1 においては主翼・尾翼を前方へ移動させることで造波抗力を一層低減できると考えられる．

11

表１．風試および造波抗力推算の模型形態．
Case No.
1
2
3
4

エリアルール Nose
○
☓
○
☓

バルジ
☓
☓
☓
A

フィレット
☓
○
○
☓

翼前方移動
☓
☓
☓
☓

5/8

○

A

☓

☓/○

6
7/9
10
11/14
12
13/15

☓
○
☓
○
☓
○

A
A
B
B
B
B

○
○
☓
☓
○
○

☓
☓/○
☓
☓/○
☓
☓/○

備考

JAXA/ISAS 遷音速
風試は Case No. 5 のみ(a) 全機抗力係数(b) 抗力低減率
図３．風洞試験および造波抗力解析の結果（Case No. 5：エリアルールノーズ C，バルジ A，迎角ゼロ）

12

50%

=

低減率

エリアルール用−エリアルール
エリアルール
用

40%

Case No.1
Case No.2
Case No.3
Case No.4
Case No.5
Case No.6
Case No.7
Case No.8
Case No.9
Case No.10
Case No.11
Case No.12
Case No.13
Case No.14
Case No.1530%造波抗力低減率

用

20%

10%
0%

-10%
-20%
-30%

0.9

1.0

1.1

1.2

1.3

1.4 1.5 1.6
マッハ数 M

1.7

1.8

1.9

2.0

2.1

図４．種々の模型形態についての造波抗力推算結果
４．まとめ
室蘭工大で研究開発している第二世代小型超音速飛行実験機の空力形状について，遷音速抗力低減のためにエリアルールに従って形状修正を試みた．マッハコーンによる断面の底面積分布がSears-Haack 曲線に極力沿うようノーズコーン底面を後退させ，主翼・尾翼の間にバルジを設けた形態について風洞試験と造波抗力解析を実施したところ，遷音速域において 5%~20%の抗力低減が推定された．更に主翼・尾翼を前方へ移動させることによって一層抗力を低減できるものと推測された．今後，エリアルール適用形状の種々の形態について遷音速・超音速風試を実施し，それらの抗力低減効果を総合的に評価する計画である．

参考文献[1] 溝端一秀，上村威揮，大石栄，「第二世代小型超音速機飛行実験機の空力設計と空力評価」，
国立大学法人室蘭工業大学 航空宇宙機システム研究センター 年次報告書 2011，pp.17-20，2012
年 7 月．
[2] 近藤賢，溝端一秀，「第二世代小型超音速機飛行実験機の飛行性能予測」，国立大学法人室蘭工業大学 航空宇宙機システム研究センター 年次報告書 2012，pp.19-22，2013 年 7 月．
[3] Jones, R. T., “Theory of Wing-Body Drag at Supersonic Speeds,” NACA Report 1284, Jan.
1956.
[4] Nikolic, V. R., and Jumper, E. J., “Zero-Lift Wave Drag Calculations Using SupersonicArea Rule and Its Modifications,” AIAA Paper 2004-217, Jan. 2004．

13

小型超音速飛行実験機の飛行性能の予測○近藤賢 （航空宇宙システム工学専攻 博士前期１年）

溝端 一秀 （もの創造系領域 准教授）

１．はじめにエアターボラムジェット・ガスジェネレータサイクル（ATR-GG）エンジンを１基搭載しマッハ２程度の超音速までの一連の飛行が可能な第二世代超音速飛行実験機（第二世代オオワシ）の設計が進められており，その空力形状 M2011 が提案されている[1]．ミッションによる推進剤搭載量の違いに対応させて，図１のように胴体長の異なる３形状が提案されており，最も小規模な M2011 Nose-A 形状の機体システム構築が当面の目標とされている．
昨年，M2011 Nose-C 形状機体の直線飛行について３自由度飛行解析が実施され[2]，マッハ 2.0 に到達するためにはエンジン回転数の 5%程度の増強または機体抗力の 8%程度の低減が必要であることが示された．本研究では，一層実際的な飛行解析ツールとして６自由度飛行解析システムを構築する．そして，当面のシステム構築目標である M2011 Nose-A
形状機体の周回・帰還飛行性能を３自由度および６自由度で解析する．

図１．第二世代小型超音速飛行実験機の空力形状

14

２．空力特性M2011 形状の総合的な空力特性を評価するために JAXA/ISAS 遷音速風洞および大阪府
立大・回流式低速風洞を用いて風洞試験が実施されている[1, 3]．その結果，マッハ数 0.3～2.0 の全域において揚力は迎角に対して十分線形であり，遷音速域の抗力ピーク（いわゆる音速の壁）は十分小さいこと，良好な空力静安定を持つこと，エルロン・エレベータ・
ラダーの操舵空力が良好であることが確認されている．これら風試データから空力係数の近似多項式を求めて飛行経路解析に用いる．また，６自由度飛行経路解析に必要なロール
ダンピング係数等の空力微係数は, 翼型の空力特性を翼幅方向に積算する手法の理論解析[4]によって推算する．これによって得られた空力微係数を表１に示す．
表１．理論解析で得られた空力微係数
CoefficientValue
0.6062114-0.2187821
0.0819851-5.0347575-0.48056994
0.0116804

CYr
Clp
Clr
Cmq
Cnr
Cnp
３．６自由度飛行経路解析システムの構築

地球の中心を原点とする三次元極座標系において，解析対象の機体を剛体として扱い，
重心の並進運動（３自由度）と重心周りの回転運動（３自由度）の計６自由度の運動方程式を解く．重心周りの回転運動の方程式は，機体固定座標系において以下の通り記述され
る[5]：
L




Mb  M   Ib  b  ( Ib) N 

ここで， Mb  L, M , N T は外力モーメントベクトル，I は慣性モーメント行列， b は角速度ベクトルである．外力モーメントを推算しつつ，この回転運動方程式を時間積分して角速度 b を 求 め ， さ ら に 時 間 積 分 し て 姿 勢 角 （ オ イ ラ ー 角 ） を 求 め る ． こ の 手 順 を
MATLAB/SIMULINK でプログラムする[6]．
風洞試験および理論解析で得られた各種空力係数，構造設計によって推算された機体重量，重心，および慣性モーメント行列，および熱サイクル解析によって推算された ATR-GG
エンジンの推力・比推力データを用い，図２のマンマシン・インターフェースによってス
ロットルおよび舵角をリアルタイムで調整しながら，上述のプログラムを実行することによって，飛行経路の解が得られる．

15

図２．６自由度飛行解析のマンマシン・インターフェイス
４．M2011 Nose-A 形状機体の周回・帰還飛行の解析
M2011 Nose-A 形状の機体（推進剤搭載量 70kg）について，大樹町滑走路から離陸して加速上昇し周回したのち帰還する飛行経路を，３自由度と６自由度で解析した．その結果を図３に示す．エンジンシステムのスロットリングは，推進剤燃え尽き前は 100%，燃え尽
き後は 0%としている．
３自由度解析に比べて６自由度解析の結果は，飛行距離および飛行時間が短く，到達マ
ッハ数も若干低い．これは，６自由度解析ではピッチトリム・ピッチ制御のためのエレベータ操舵による揚力低減やフラッペロン操舵による抗力増加を勘定に入れており，３自由度解析に比べて全体的に全機揚力が低下し全機抗力が増加しているためと考えられる．これらのことから，６自由度解析に比べて３自由度解析は楽観側の解を与えることがわかる．
６自由度解によれば，M2011 Nose-A 形状の機体は，大樹町の滑走路を東北東 80°方向へ離陸した後，時刻 147 秒に高度 11km において推進剤燃え尽きの瞬間に最大マッハ数 1.2
弱に到達する．その後，時刻 150～400 秒にかけて右旋回・左旋回を経て帰還経路に入り，滑走路近傍の-100°方向に帰還する．ダウンレンジ 50km 程度，ロール角 50 度程度，旋回半径 5km 程度である．延べ飛行時間は 550 秒程度，軸加速度（ax：機首方向）は離陸滑走時に最大 1G 程度，横加速度（ay：右翼方向）は旋回時に最大 1G 弱，垂直加速度（az：腹方向）は±2G 弱である．最大迎角は離陸引き起こしの際に 8 度程度である．総合して，概ね良好かつ実現可能な飛行経路である．

16

16km
22km

(a) 飛行経路(b) 飛行経路拡大図(d) 飛行マッハ数の履歴(c) 飛行高度の履歴(f) ６自由度解析による加速度の履歴(e) ６自由度解析による迎角と姿勢角の履歴

図２ ３自由度および６自由度解析による周回・帰還飛行解析結果
５．まとめ
第二世代小型超音速飛行実験機について，風洞試験および理論解析に基づく空力特性データ，ATR-GG エンジン性能解析データ，および構造設計に基づく重量・重心・慣性モー
メント推算データを用いて３自由度および６自由度飛行経路解析を行い，M2011 Nose-A 形状の機体の周回・帰還飛行性能を予測した．６自由度解析によれば，大樹町の滑走路を東
北東方向へ離陸した後，時刻 147 秒に高度 11km において最大マッハ数 1.2 弱に到達し，旋回ののち，滑走路近傍まで帰還できることが予測された．ダウンレンジは 50km 程度，
延べ飛行時間は 550 秒程度である．
今後は，タッチダウン時の対地速度を低減できる着陸接近飛行の方法を，３自由度および６自由度飛行解析によって予測する計画である．

17

参考文献
[1] 溝端，上村，大石，「第二世代小型超音速飛行実験機の空力設計と空力評価」，国立大学
法人室蘭工業大学 航空宇宙機システム研究センター 年次報告書 2011，pp.17-20，2012 年
7 月．
[2] 近藤，溝端，「第二世代小型超音速飛行実験機の飛行性能予測」，国立大学法人室蘭工業
大学 航空宇宙機システム研究センター 年次報告書 2012，pp.19-22，2013 年 7 月．
[3] 鈴木，近藤，大石，溝端，「第一世代および第二世代小型超音速飛行実験機の操舵空力の評価」，国立大学法人室蘭工業大学 航空宇宙機システム研究センター 年次報告書 2012,
pp.11-14, 2013 年 7 月．
[4] 加藤寛一郎，大屋昭男，柄沢研治，「航空機力学入門」，東大出版会，pp.51-110，1982．
[5] McFarland, Richard E., “A Standard Kinematic Model for Flight Simulation atNASA-Ames,” NASA CR-2497, 1975.
[6] Jackson, E. B., and Cruz, C. L., “Preliminary Subsonic Aerodynamic Model for
Simulation Studies of the HL-20 Lifting Body,” NASA TM4302, 1992.

18

小型超音速飛行実験機の予備的飛行試験のための縮小機体の設計・製作○渡口

翼 （航空宇宙システム工学専攻 博士前期 1 年）

山中 大樹（機械航空創造系学科 4 年）
梅村 諒磨（機械航空創造系学科 4 年）
溝端 一秀 （もの創造系領域 准教授）

１．はじめに風洞試験は，風洞壁や模型スティングの存在，流れの乱れ度やレイノルズ数，など実飛行環境と異なる制約条件を有するため，その取得データの信頼性に制約がある．したがって，有翼飛行
体の実飛行環境における空力特性・飛行特性を検証するためには，風洞試験だけでは不十分であり，飛行試験が欠かせない．そこで，小型超音速飛行実験機「オオワシ」について繰り返し安全に飛行試験を実施し，質・量ともに十分なデータを取得することを狙って，製作・飛行の比較的簡便な縮小機体を準備している．昨年度[1]に引き続き第一世代オオワシ（M2006prototype 形状）の 1/2 スケール機体の設計・製作を進めている．推進器としては電動ダクテッドファンユニット
を用いる．
２．機体構造の設計改良
昨年度設計・試作した縮小機体[1]において，推進用モーターマウントの構造強度の不足，胴体・尾翼接続箇所の整備性の不良，ラダー操舵系の機体表面への露出，等の課題が残っていた．そこで，推進用モーターマウントを木質角材からアルミアングル材に変え，胴体・尾翼接続構造を改良し，さらにラダー操舵系を後胴部に内蔵することとした．その概略を図１に示す．
また，後述の地上静止推力の計測の結果，エンジンナセル内面の構造部材が推進器の空気流を阻害しているものと推定されたため，推進器搭載位置をナセル前端に変更し，さらにナセル内面に滑らかな円筒壁を設けた．その概略を図２に示す．
これらの構造改良を施し，推進器，推進用バッテリー，無線操縦受信機，等の機器を搭載した
ところ，全備質量 4.38kg となり，公称推力 4.2kgf（41.2N）に対する推重比は 0.96 となった．

(b) 着脱の容易な尾翼・胴体接続構造と内蔵されたラダー操舵系(a) アルミ製のモーターマウント

19

(c) 着脱の容易な水平尾翼接続構造
図１．縮小機体の構造の改良

図２．推進器搭載方法の改良
３．地上静止推力の計測と推進器搭載方法の改良
飛行中に推力を計測することは困難であることから，縮小機体の推力を地上で計測した．そのための機材を図３に示す．水平な机の上に滑らかに動く台車を置き，それにオオワシ縮小機体を載せて結わえる．無線操縦によって機体搭載推進器を作動させると，推力によって機体・台車は前方に進もうとして，機首前方に設置されたロードセルを押す．このロードセルの出力信号と，無線操縦のスロットル制御信号を収録し，スロットル制御信号と推力の関係を得る．なお，スロ
ットル制御信号の収録のために，オオワシ 1 号機の飛行試験に用いたのと同一のオンボードアビ
オニクス機器（Sylphide）を用いる．また，予め 1～5kg の分銅を用いてロードセルを較正している．
得られたスロットル制御信号と推力の関係を図４に示す．推進器搭載方法改良前（Before
revision）において最大推力は 1.40kgf（13.7N）となっており，後述の定常水平飛行のための必要最低推力 1.37kgf（13.4N）をわずかに上回程度であり，加速上昇が困難である．ダクテッドファンユニットの公称推力 4.2kgf（41.2N）に比べてこのように低推力になった原因としては，エンジンナセル内に構造部材や降着装置が突き出ており，推進器の空気流を阻害しているものと推定される．また，グラフが非線形となっているのは，ラジコン機器の操縦信号とスピードコント
ローラ（操縦信号に基づいて推進モーターの回転数を制御する）の整合性が取れていなかったものと推定される．
そこで，推進器搭載方法を改良するとともに，ラジコン機器およびスピードコントローラを調整したところ（After revision），図４に示されたとおり，最大 2.4kgf の推力が得られるに至った．

20

(a) 台車，縮小機体，およびロードセル(b) スロットル制御信号収録器（Sylphide）

図３．推力計測のための機材

図４．計測されたスロットル信号と推力の関係
４．飛行性能の予測
M2006prototype 形状については亜音速風試によって空力特性データが得られている．その揚力係数・抗力係数データを利用して，縮小機体の海面上での定常水平飛行速度と必要推力の関係
を 推 算 す る と 図 ５ の 通 り で あ る ． 公 称 推 力 4.2kgf （ 41.2N ） で は 定 常 水 平 飛 行 速 度 は
61.0m/s(219.7km/h)となる．定常水平飛行のための必要最低推力は 1.37kgf(13.4N)であり，その時の定常水平飛行速度は 26.7[m/s](96.1[km/h])と推算される．推進器搭載方法改良後（After

revision）の定常水平飛行速度は 44.3[m/s]（159.5[km/h]）と推算される．
また，飛行速度と上昇率の関係は図６の通りとなる．最大上昇率は，公称推力 4.2kgf（41.2N）
では 20.08[m/s]，
推進器搭載方法改良前（Before revision）では 0.06[m/s]，改良後（After revision）
では 6.15[m/s]と推算される．

図５．飛行速度と必要推力

図６．飛行速度と上昇率
21

５．慣性モーメント・慣性乗積の計測・推算
捻り振り子の原理を用い，捻り振動周期を計測することによって，縮小機体の慣性モーメントを推算した．振動周期計測の様子を図７に示す．また，慣性乗積の測定は不可能であることから，公表されている実機の慣性乗積を元に，質量比例・寸法２乗比例の相似則を仮定して推定した．
これらの結果を表１に示す．慣性モーメントの推算結果は，質量・代表寸法 2 乗で無次元化した
相似パラメタがプロトタイプ機と概ね同等となっており，良好な結果といえる．慣性乗積の推定においては，この相似パラメタの値を機種によらず一定と仮定している．

(b) Y 軸周り(c) Z 軸周り

(a) X 軸周り
図７．慣性モーメントの計測

表１ 慣性モーメントの計測値および慣性乗積の推算値
機体全備質量 m [kg]
全長 L [m]全幅 b [m]全高（脚除く）h [m]慣性能率 Ixx [kgm2]
相似パラメタα
慣性能率 Iyy [kgm2]
相似パラメタβ
慣性能率 Izz [kgm2]
相似パラメタγ
慣性乗積 Ixz [kgm2]
相似パラメタδ

F16
9299
15.03
9.45
3.57
12870
0.015575670
0.036085550
0.02921331
0.00267

プロトタイプ機
27.3
3.2
1.61
0.665
0.94269
0.01338.30130.0297
11.2361
0.03210.15497
0.00267

縮小機体
4.38
1.59
0.8
0.324
0.029730.0106
0.38150.03450.48620.03510.00601
0.00267

６．タキシングおよび滑走試験
離陸滑走の安定性や加速性を検証するために，白老滑空場の滑走路にて，タキシングおよび滑走試験を行った．その様子を図８に示す．直進安定性の不足が確認されたため，降着装置の剛性やアラインメントの調整が課題となる．

22

図８．タキシング・滑走試験の様子
７．まとめと今後の展望
室蘭工大の小型超音速飛行実験機プロトタイプについて，飛行試験を繰り返し実施するため，
1/2 スケールの縮小機体の設計・製作を進め，強度ならびに整備性の観点から機体構造を改良した．
推力改善のために推進器搭載方法を改良し，2.4kgf の推力を得た．この推力による飛行性能は概ね良好と推測された．さらに慣性モーメント・慣性乗積を計測・推算した．今後，走行安定のために降着装置を改良し，本縮小機体を予備的飛行試験に供することによって， プロトタイプ形状（M2006prototype）の低速飛行特性，各舵面による空力性能，動的飛行特性，等の解明を進める計画である．
参考文献[1] 渡口翼，福士誠，溝端一秀，「小型超音速飛行実験機の縮小機体の設計・試作」，国立大学法人室蘭工業大学 航空宇宙機システム研究センター 年次報告書 2012，pp.23-26，2013 年 7 月．

23

小型超音速飛行実験機のエアインテークの空力特性○溝端

一秀 （もの創造系領域 准教授）

塩田光平 （機械航空創造系学科 4 年）

１．はじめに
第二世代小型超音速飛行実験機（M2011 空力形状）では，高速飛行に適合した GG-ATR エン
ジン一基が胴体後端に搭載され，胴体腹面にインテーク（空気取り入れ口）が設置される．そのエンジン性能はインテークの空力性能に大きく依存し，また飛行性能を見積もる上でインテーク
抗力は非常に重要である．そこで，本研究は第二世代実験機に搭載予定のエアインテークの空力特性を評価することを目的とし，JAXA/ISAS 遷音速風洞および室蘭工大超音速風洞において風洞試験を実施する．
２．インテークの設計
本研究で取り扱うインテークは，第二世代実験機の胴体腹面に設置される矩形インテークであり，GG-ATR エンジン所要の流量・圧力の空気を供給できるように，超音速ディフューザ・亜音
速ディフューザを組み合わせ，マッハ 2.0 で流量捕獲率 100%になるように設計されている．その形状・寸法を図１に示す．また，インテークを搭載した第二世代実験機の予想図を図２に示す．

図１．インテークの形状・寸法

図２．インテークを搭載した第二世代実験機

３．インテークの空力特性指標
インテークの空力性能を評価する指標として，インテーク抗力，全圧回復率，および流量捕獲率の３つがある．風洞試験において，インテーク抗力は天秤を用いて計測される．全圧回復率と流量捕獲率の推算手法を以下に述べる．
３－１． 全圧回復率
全圧回復率は一様流全圧 P0，一様流静圧 P∞ ，およびインテーク出口全圧 P0i を用いて，
p 

p0i p0i P

P0
P P0

・・・・・・・・・・(1)ここで，等エントロピー関係式
P0    1 2   1
 1 
M 
P 
2


24

・・・・・・・・・・(2)

が成り立つため，全圧回復率は p  1 

 1



 1 p0i
M2
2
P


・・・・・・・・・・(3)と表される．したがって，風試において一様流マッハ数 M，一様流静圧 P∞ ，およびインテーク
出口全圧 P0i を計測することで，全圧回復率を算出できる．
３－２．流量捕獲率
ランプから発生する斜め衝撃波によって曲げられた流れの偏角はランプの傾斜角に等しいため，
シュリーレン画像から図２のようにカウルリップにぶつかる流線を作図し，インテークに流入す
る流管の断面積を求め，インテーク入口前面投影面積で割ることで流量捕獲率を推算する．

図２．流量捕獲率推算法
４．風洞試験
４－１．風洞試験模型
インテーク入口およびダクト内の流れの計測・観察を容易にするために，風洞閉塞率が過大にならない範囲で模型をできるだけ大きくする必要から，縮尺比 1/3 の風試模型を設計製作した．
その形状・寸法を図３に，製作された風試模型を図４に示す．実機の胴体を模擬する胴体模型を設け，インテークダクト出口は胴体模型の対称軸上に置いている．実機のインテークダクト出口
は円筒形状であるが，この風試模型では製作の容易さの観点から同等流路面積を持つ矩形出口としている．
エンジンの空気吸い込み流量を模擬するために，インテークダクト出口にオリフィスを設け，その開口率（面積比）を 100%，87.5%，75%の３通りとしている．また，オリフィスにピトーレ
イクを共締めし，オリフィス直近上流の全圧分布および静圧を計測する．オリフィスとピトーレ
イクの寸法・配置を図５に示す．
インテーク模型と胴体模型にはたらく抗力を計測するために，インテーク模型のカウル部品を内装天秤に装着できるようにしている．今回は，胴体模型をインテーク模型に接続し，胴体模型
およびインテーク模型の総体の抗力を計測する．今後，胴体模型をスティング上面に直接接続し，胴体模型とインテーク模型を遊離させることによって，インテーク模型のみの抗力を計測することも可能である．
インテーク模型の側面は，当面アルミ板（t5）としているが，その一部を透明アクリル板とすることによって内部流れの光学観察も可能である．

25

出口にオリフィス（穴開きアクリル板 t5）を設置．
ピトーレイクにより全圧・静圧を計測．

側面は t5 アルミ板．
部分的に t5 アクリル板（内部流の可視化）．

図３. インテーク模型の設計

図４. 製作されたインテーク模型

図５．オリフィスとピトーレイクの寸法・配置
26

４－２．試験条件
通風条件は，ISAS 遷音速風洞においては，マッハ 1.3～0.7 のマッハスイープおよび迎角ゼロである．また，模型にはたらく抗力が天秤秤量を超えないように，一様流全圧として遷音速風洞の運転可能最小値 1.5kgf/cm2 に設定する．室蘭工大超音速風洞においては，マッハ 2.0 ノズルブ
ロックを用い，一様流全圧 100kPa（大気圧）である．実際は，風洞のラバルノズルの流路拡大部分にインテーク入口が配置されるため，インテーク入口マッハ数は 1.6 程度である．また，整流金網の圧損によって一様流全圧は 80kPa 程度である．
４－３．風洞試験の結果及び考察
３種類の開口率のオリフィスを用いて，マッハ 1.3～0.7 のマッハスイープ試験によって計測した抗力値を図６に示す．オリフィス開口率の違いによる抗力の差はほとんど見られない．また，
インテークダクト出口における全圧分布を図７に示す．各通風マッハ数において，３種類のオリフィスによるデータが１本の曲線に乗っていることから，オリフィス開口率による全圧分布の差は小さいものと判定される．さらに，紙面の都合から掲載していないが，インテーク入口近傍のシュリーレン画像も，３種類のオリフィスによってほとんど差が見られない．付け加えて，図７の全圧分布において上部 z = 0~20mm と下部 z = 45~60mm で全圧が小さいことがわかる．これらのことから，インテークダクトの上下壁面近くで流れが剥離しており，開口比 75%のオリフィ
スよりも実質的な内部流路面積が狭くなっていると考えられる．
実機の飛行解析から得られる飛行経路において，エンジン所要の空気流量・流量捕獲率の推算
値と，今回の風試で計測された流量捕獲率を，図８に示す．風試結果は，インタークダクト内流れが大きく剥離しているものの，所要流量捕獲率推算値を上回る結果となっている．このことから，内部流れの剥離を低減できればインテーク入口サイズを小さくでき，インテーク抗力の低減につながるものと期待される．

図６．３種類のオリフィスによる抗力の比較

27

図７．インテークダクト出口の全圧分布

図８．飛行経路上でのエンジン所要空気流量・所要流量捕獲率の推算値，および風試結果
５．まとめ
第二世代小型超音速飛行実験機（M2011 形状）の胴体腹面に設置するインテークの空力特性を評価するため，M2.0 を設計点とする 1/3 スケール・インテーク模型を設計製作し風洞試験を行った．エンジンの吸い込み空気流量を模擬するために，インテークダクト出口にオリフィスを設置
し，開口率を 100％，87.5％，および 75％の３通りとした．また，オリフィス直近の上流側でピトーレイクを用いて全圧分布および静圧を計測した．
その結果，オリフィス開口率の違いによって抗力，入口シュリーレン画像，出口全圧分布に差はほとんど見られなかった．また，上下壁面近くに全圧の低い領域があることが確認された．これは，インテークのダクト形状の曲げ勾配が大きいために，ダクト内で流れが剥離しているためと考えられる．インテークダクト形状の曲がりを緩やかに改善することで空気流の剥離を防ぎ，所定の圧力・速度の空気をエンジンに供給できると思われる．その場合，機体胴体下部に露出するインテーク部分は長くなるが，2011 年度の風洞試験[1]より，インテーク長さが M2011 形状の抗力または空力特性に与える影響は小さいことが判明している．また，その場合，所要空気流量を得るためのインテーク入口サイズを小さくできるため，インテーク抗力の低減に寄与できると考えられる．
今後，インテーク模型のみの抗力の計測，内部流れの観察，等を進め，インテークの詳細な空力特性を把握する計画である．

28

小翼列によるアクティブ制御を搭載したラジコン機風洞試験○ 髙木

正平（航空宇宙機システム研究センター

教授）

田中

清隆（航空宇宙システム工学専攻

博士前期 1 年）

上田祐士（航空宇宙システム工学専攻

博士前期 2 年）

１．研究背景および目的
現行の小型無人機は予期せぬ突発的な外乱に対しての安全飛行の対策が十分とは言い難い．汎用ラジコン機の主翼前縁近傍の境界層は層流であり，対気流角の急変により翼上面の流れが大規模に剥離しやすく，失速して墜落の危険に晒される．従来の剥離制御技術は強靭性，小型化，省電力化など克服すべき問題を抱える．ロバスト性を有し，必要な時のみ流れを制御できるような
Active 制御機構が求められている．
本研究の最終目的は，ラジコン飛行機の失速を検知し，必要に応じて境界層を乱流化させ失速回復させる Active 制御手法を考案し，それらを実装して飛行実証することである．今年度は昨年の風洞試験によって得られたボルテックスジェネレータ（VG と略記）による Active 制御の効果を実機に搭載し飛行実証を行った．ここでは，飛行実証の前段階として実施した風洞試験について報告する．

２．ボルテックスジェネレータを用いた Active 制御
平成 25 年度に実施したセンター所有の小型低乱風洞試験では，VG を用いた Active 制御についてその効果を確認した．その実施形態は，翼弦長の 30%位置に 40mm 間隔で左右 8 個ずつ取り付けた VG を旗のように主流に対して±30°往復運動させることによって強い縦渦を導入する
ことで剥離を制御する．VG の往復振動にはサーボモーターを使用し，振動周波数は最大 3.5Hzまで振動させることができる．図 1 に VG による Acitve 制御を行ったときの六分力天秤で評価した迎角に対する揚力と抗力係数特性を示す．

図1

失速特性

図 1 から，VG を 3.5Hz で振動させたときと振動させなかったときを比較すると，揚力係数は約8%の増加し，失速迎角は 3°上昇していることがわかる．今年度は，風洞試験で明らかにした
VG による Active 制御をラジコン機に搭載し，その効果を検証した．
29

３．VG を搭載したラジコン機の風洞試験
2 月 23 日から 28 日の期間，首都大学東京の大気開放型の回流式低速風洞を借用し，VG を搭載したラジコン機の風洞試験を行った．図 2 に示す機体は市販されているカルマートアルファ 40
トレーナーEP/GP（京商）を使用し，主翼は VG と制御機器を搭載するために独自に設計・製作
したものを用いた．表 1 に機体諸元を示す．
図 3 および図 4 に搭載したリン青銅板の VG の寸法と外観を示す．VG は 65mm 間隔で配置し，片翼 7 個の合計 14 個取り付けた．VG を往復運動させる動力源はラジコン用サーボモーターを使用し，ラジコン受信機に接続することで遠隔操作が可能となっている．なお，往復周波数は 3.5Hzで振動させた．
この風洞試験では 6 分力天秤による揚力や抗力などの力計測，オイルフロー法による剥離流の可視化を行い，実証飛行を行う前の基本特性を把握することが狙いである．
表1

機体諸元

全長

1300 mm

翼幅1600 mm

翼弦長

270 mm

翼型

オリジナル準対称

全備重量

約 2650 g

図２ 実験使用機体

30

図3

VG 寸法

図4

VG の概観

31

４．実験結果
図 5 に VG を作動させたときとさせなかったときの迎角に対する失速特性を示す．この特性は，ラジコン機のプロペラを停止した状態で，流速を 7.5m/s のときの結果である．この図から，VGの往復運動の有無にかかわらず揚力係数および抗力係数はほぼ一致していることがわかる．また，
図 6 および図 7 にスモークワイヤー法による可視化の結果を示す．図 5 から失速迎角は 18°であることが確認できたため，この迎角における VG の効果によって翼上面の流れを観察した．図6 および図 7 から VG の動作に関係なく流れはほぼ同じであることがわかる．以上から，今回の実験では，昨年度実施した風洞試験から得られた効果を確認することができなかった．

図5

図6

迎角 18°,

ラジコン機失速特性

VG 制御なし

図7

迎角 18°,VG 制御あり５．考察
今回の首都大学東京における風洞試験からは期待していた効果が得られなかった．この原因としては現在も考察中であるが，考えられる原因としては，昨年度実施した風洞試験翼模型の大きさと VG の取り付け位置を相似形のままラジコン機主翼を設計したことが挙げられる．これによって，主翼上面を流れる流体の剥離位置が上流側に移動してしまい，既に流れが剥離している位
置で VG を動作させていたことになる．そのため，揚力係数だけでなく，抗力係数にも影響を与えない結果となったと考えられる．以上の考察をふまえ，今後は風洞試験においてさらなる検証を行っていく予定である．
32

小翼列によるアクティブ制御を搭載したラジコン機実証飛行試験○ 髙木

正平（航空宇宙機システム研究センター

教授）

田中

清隆（航空宇宙システム工学専攻

博士前期 1 年）

上田祐士（航空宇宙システム工学専攻

博士前期 2 年）

１．目的
昨年度の本学風洞において行ったボルテックス・ジェネレータ（VG）による Active 制御の効
果は翼上面流れの剥離を遅延・防止させるために有効な手法であることが確認された．そこで本
年度は実機においても VG による剥離制御効果が有効であることを確認するために，市販ラジコン飛行機とボルテックスジェネレータ（VG）を搭載し独自に設計・製作した主翼を用いて実証飛行試験を行った．この実証飛行試験の目的は，VG の効果を検証することを主眼とし，実機を用いた実証飛行は初となるため，電波到達範囲の把握や自作した主翼の強度と舵面の利き，慣性航法装置，データロガー，五孔ピトー管のデータ収録，VG 動作基準となる逆流検出器など，飛行実証に必要な基本的な動作確認を行うことも主眼に準じた目的とした．

２．実施概要実証飛行試験は 3 月 3 日と 4 日の 2 日間で実施し，図 1 に示す静岡県富士川滑空場にて実証飛行を行った．図 2 に示すように，実証飛行に用いた機体は本年次報告書「小翼列によるアクティ
ブ制御を搭載したラジコン機風洞試験」で用いた機体と同様のものである．飛行試験は 1 日目に
2 回，2 日目に 5 回の計 7 回実施し，Active 制御を行ったときと行わなかったときの機体の姿勢変化や飛行経路，飛行速度などの計測を行った．

図2
図1

実証飛行機体富士川滑空場

３．計測機器
ラジコン機に搭載する計測装置は軽量かつ小型の装置が要求され，遠隔で操作・監視できる必要がある．そのため，図 3 に示すようなノートパソコンによってリアルタイムで飛行状態をモニタできる長田電機社製の TinyFeather を慣性航法装置として使用した．図 4 と図 5 に示す五孔ピトー管や逆流検出器の出力信号を記録するデータロガーとしては，図 6 に示す MUSIN 社製の
33

Data Logger を使用した．

図3

Tiny Feather 送受信回路(左)とリアルタイムモニタ用パソコン(右)図4

図5左翼端に取り付けた 5 孔ピトー管

左翼表面に実装した逆流センサーと検出器

34

図6

MUSIN データロガー

４．実験結果
ラジコン機を低速で高進入角飛行させることで主翼上面に大規模な剥離を発生させて不安定な飛行条件を作り，VG を動作させて剥離抑制ができれば安定飛行が期待される．実際このような飛行形態を実現し，VG の効果を調べた．飛行観察と収録データを処理した範囲ではその効果を十分に確認できなかった．この結果は首都大学東京における風洞試験とも整合している．
次に，各種計測装置の動作と機能確認を行った．自作した主翼の飛行特性と構造強度は問題な
く，購入した主翼に比べて軽量であるにもかかわらず，構造強度は十分であった．図 7 に飛行中の写真を示す．
慣性航法装置 TinyFeather は緯度・経度・高度・速度・ピッチ角・ヨー角・ロール角を計測することができ，XBee と組み合わせることで前記データをリアルタイムでパソコン画面で監視できる．図 8 は Tiny Feather で計測した飛行経路の一例である．また，図 9 にリアルタイムパソ
コン画面を示す．今回の動作確認で慣性航法装置である TinyFeather に関しては良好にデータ取得できることを確認した．
一方で，5 孔ピトー管と逆流検出器のデータを取得するための MUSIN Data Logger の電波到達距離はカタログ値に比べて小さく，リアルタイムモニタ及びデータ収録を行うことができなか
った．

図7

図8

飛行写真

35

飛行経路

図9

リアルタイムパソコン画面

５．まとめ
VG による Active 制御をラジコン機に搭載した飛行実験では，直前に実施した風洞試験と同様に明白な効果を確認することができなかった．しかし，自作した主翼の強度確認や各種計測装置の動作確認は完了し，その結果，データロガーの不具合が明らかとなった．今後は飛行試験範囲内の無線操作で安定して動作するデータロガーの選定を行う必要がある．
また主眼としていた VG による Active 制御に関しては，実機レイノルズ数の風洞試験で再検証
し，VG の形状や位置についてさらに追求していく必要がある．

36

層流後退翼の境界層遷移制御に関する風洞実験○ 髙木 正平（航空宇宙機システム研究センター 教授）
伏木田 裕介（航空宇宙システム工学専攻 博士前期 1 年）
正木 陽（航空宇宙システム工学専攻 博士前期 1 年）

１． 研究目的
後退翼面上に発達する三次元境界層は，横流れ不安定によって二次元境界層と比較して早期に乱流化することは周知である．翼が後退すると圧力勾配の方向と主流方向が異なることから外部
流線と直角方向に横流れ速度成分が誘導され，その分布は変曲点を持つために非粘性型不安定が生じることに起因している．この不安定は横流れ不安定と呼ばれ，この不安定に起因して成長する定在型の縦渦は，翼表面粗さに極めて敏感であることも良く知られている．その不安定性を利用し，本研究では後退翼前縁近傍に微小な粗度列を添付することで縦渦の成長を制御し，乱流遷移の促進，あるいは遅延技術を獲得することを目的とする．すなわち，自然遷移で見られる縦渦の波長に粗度間隔を選べば遷移が促進され，波長の 60～75%程度に選べば遅延されることが期待できる．

２． 実験装置
２.１ 低乱熱伝達風洞
本実験では，東北大学流体科学研究所所有の低乱熱伝達風洞を用いた．開口部断面は正八角形，
開口面積は 0.544[m2]であり，また風速は 28～35[m/s]で実験を行った．
２.２ 層流翼
本実験では翼弦長 c=600[mm]，翼幅 S=910[mm]の層流翼 NLF(2)-0415 を用いた．模型設定条
件は迎角α=-4[deg]，後退角Λ=45[deg]とした．設置した模型側面写真を図１に示す．

図.1 設置した模型側面写真 (NLF(2)-0415)
２.３ 赤外線カメラ
本実験では日本アビオニクス社製の赤外線カメラである Advanced Thermo TVS-500 を使用した．主な仕様は動作環境-10[℃]～50[℃]，有効表示画素数 320(H)×240(V)，最小温度分解能 0.0637

[℃]以下，空間分解能 1.68[mrad]である．

３． 実験結果
３.１ 翼表面静圧分布計測

主流速度 U=28[m/s]において 16 点の静圧孔から得られた静圧分布を図 2 に示す．図 2 より70[%]翼弦長までは順圧力勾配であり，翼の前縁近傍は特に圧力勾配が大きいことから，この領域では強い横流れ速度成分が誘導される．従って境界層遷移を制御するための粗度の貼付には効
果的な位置と考えられる．

図 2 NLF(2)-0415 の静圧分布
３.２ ナフタレン昇華法による可視化実験
ナフタレン昇華法による模型表面流れの縦渦及び乱流遷移点の可視化実験を行った．図 3 に粗度のない滑面状態，すなわち自然遷移状態において，ナフタレン昇華法による翼表面の可視化結果を示す．縦渦の挙動に伴って生じる筋状の痕跡から，縦渦の平均波長はλ=7[mm]であった．
その結果より，縦渦を成長促進させるには粗度間隔をλ=7[mm]，逆に成長を抑制するにはλ=4.7～5.2[mm]程度に決定すればいいことが分かった．そこで厚さ t=50[μm]，幅 1.2mm，長さ 10mmの粗度を水平線に対して角度 15 度で粗度の上流先端が x/c=0.023C の位置に貼り付け，粗度間隔をλ=7[mm]とλ=5.2 [mm]，また貼付しないケースを同一模型内で，主流速度 U=35[m/s]で実施した．その結果を図 4 に示す．ナフタレンの昇華が早いほど乱流遷移が早いと診断されることから，粗度間隔がλ=5.2[mm]，粗度なし，粗度間隔λ=7[mm]の順に乱流遷移が進行することが分かった．この結果より，粗度の間隔を自然遷移で観察される波長の 75%にすれば乱流遷移を人工的に遅延できることが実証できた．

38

図 3 ナフタレン昇華法による滑面の可視化(a) 通風前のナフタレン塗布状況(b) 通風後のナフタレン塗布状況

図 4 ナフタレン昇華法による粗さ貼付した表面の可視化
３.３ 赤外線カメラを用いた可視化実験
図 5 には，赤外線カメラによる模型表面の温度分布画像を示す．模型表面に断熱フィルム(50μm)を貼付し，通風前の予熱した温度分布と通風後の温度分布との温度差をマッピングした画像であるが，模型表面を断熱フィルムで覆ったけれども，模型が熱伝導性の良いアルミ製であるため模型表面温度の均一化が予想以上に早く，縦渦の存在を示す筋状の痕跡は確認することができなかった．

図 5 赤外線カメラによる模型表面温度分布
３.４ 熱線風速計による縦渦の計測
主流速度を U=28[m/s]として図 6 には c=50%において前縁に平行なスパン方向に熱線をトラ39

バースして計測した平均速度分布を示す．粗度を貼付すると縦渦の存在を表す正弦波構造の振幅をはっきりと捉えており，ナフタレン昇華法の結果と整合している．ただ，滑面の自然遷移過程
では，振幅は非常に小さい．c=60%，70%位置でも同様の計測を行ったところ，下流方向に振幅の減衰傾向を示した．

図 6 熱線風速計による前縁平行方向の 50％翼弦位置における速度分布 u/U(U=28[m/s])

４.まとめ
横流れ不安定が遷移過程を支配している三次元境界層において，横流れ不安定の強い模型前縁近傍に粗度間隔を選択することで，遷移を遅延，ないしは促進できることを実証できた．自然遷
移で観察される縦渦間隔に粗度を添付すると遷移が促進され，75%程度の間隔で粗度を貼付すると遷移が遅延されるという結果であった．また，熱線計測により縦渦の存在を示す平均速度の分布を取得でき，可視化実験から得られた結果と整合していることを確認した．今後はより鮮明な可視化写真を得るためにナフタレン昇華法の習熟が求められる．粗度の形状，厚み等パラメータの影響も詳しく調査する必要があり，さらに明確な確証試験(熱線計測等)が望まれる．

40

超音速風洞における小型トラバース装置
○ 正木
髙木

陽

（航空宇宙システム工学専攻 博士前期 1 年）

正平 （航空宇宙機システム研究センター 教授）

1.はじめに航空宇宙機システム研究センターでは，大気中を高速で飛行するための基盤技術の創出を目的としてプロジェクト研究が進められている．
大気中を高速で飛行するための推進機である ATR-GG 開発の一環として，要求する推力を満たすた
めのインテークの捕獲流量の数値解析と超音速風洞実験が行われている．しかし，固定された位置における総圧・静圧計測から正確な捕獲流量を推定することは難しいことから，圧力等のプローブを通風中にインテーク断面内を移動させるための小型トラバース装置を新規に製作した．この装置の概略と使用方法を説明する．
２.計測装置２.１小型トラバース装置の三面図27製作した小型トラバース装置の三面図と外観写真を図１に示す．

40
台座

側面図

135

正面図

台座
平面図

概観図

図 1 小型トラバース装置の三面図と外観写真．寸法数値単位 mm
41流れと直角方向の厚みは可能限り薄くなるよう設計した．また，装置の剛性を確保するために，材料はス
テンレスを用いた．高さも低姿勢となるよう配慮したため，モーターは横置きとした．プローブを搭載する
移動部は超音速流中でも滑らかに摺動するよう，２本の超薄型ガイドレールを用いた．

２．２ 装置概要
装置のモーターは，オリエンタルモーター社の PK513PB を使用している．移動長さは，最大で約 60mm
であり，最大移動速度は約 4mm/s である．性能については以下の表としてまとめる．
表 1 装置性能
五相モーター

PK513PB（オリエンタルモーター社）

装置の大きさ

108×27×135 [mm]プローブ固定台座

44.4×40 [mm]移動長さ

約 60mm

最大移動速度

約 4mm/s

２．３ 使用方法
製作した小型トラバース装置は，既存のスティングと固定できるよう台座は加工されているが，それ以外の固定については新規に台座を製作する必要がある．
トラバース装置の制御システム構成を図２に示す．装置は超音速風洞計測室に用意された PC にインストールされた MEXP01 を用いて遠隔制御する．MEXP01 から送られた制御信号は，測定部近傍に設置されたコントローラを経たのち，ステッピングモータを駆動するためのドライバーに送られる．
なお，使用に当たっては最大通風時間で制約されるので，プローブの移動範囲が制限を受けるので注意が必要である．

モーター

ドライバー

DC24V 電源コントローラDC24V 電源

超音速風洞測定部

PC
超音速風洞計測室

図2

トラバース装置の制御システム構成図
42

小型無人超音速機の応力解析と剛性解析○勝又 暢久 （もの創造系領域 助教）
樋口 健（もの創造系領域 教授）

石田 貴大 （航空宇宙システム工学専攻 博士前期 1 年）
小杉 義和 （航空宇宙システム工学専攻 博士前期 1 年）
田中 啓太 （機械航空創造系学科 4 年）
丹治 陸 （機械航空創造系学科 4 年）

１．はじめに室蘭工業大学・航空宇宙機システム研究センターでは，小型無人超音速飛行実験機の研究開発が進められている．2010 年度夏期および 2011 年度夏期には，プロトタイプ機であるオオワシ 1号機を用いて亜音速飛行試験を行った．この飛行試験結果に基づいて，現在はオオワシ 2 号機の開発，設計と研究が進められている．オオワシ 2 号機は目標飛行速度をマッハ 2 程度に設定して
おり，1 号機に比べて全長は約 2 倍，全幅は約 1.5 倍の大きさで設計を進めている．そのため機体構造だけでも重量は 10 倍以上増加する．さらに超音速飛行に適したエアターボラムジェット・
ガスジェネレーターサイクル（ATR-GG）エンジンやその燃料となるバイオエタノール，LOX な
どの搭載機器を考慮すると燃料満載時で約 400kg 程度の機体重量になることが予想されている．
そこで本稿では，搭載機器質量を加味した機体構造の静荷重解析と固有振動数解析を行い，オ
オワシ 2 号機の胴体部構造の基本構造設計について検討した結果を以下に示す．
２．胴体部構造のモデル化
本研究では，3 次元 CAD ソフト「CATIA」を用いてオオワシ 2 号機胴体部のソリッドモデルを作成した．また作成した 3 次元 CAD モデルを汎用有限要素法(FEM)解析ソフトウェア「Patran/ Nastran」を用いてソリッド要素により有限要素モデル化し，胴体部構造の静荷重解析を行った．

図１．3D CAD モデル
表１ 材料物性値と胴体部重量
Material E [ kgf/mm2 ] ν [ - ] G [ kgf/mm2 ] t [ mm ] ρ [ kgf･s 2 /mm4 ] Weiht[kg]Frame A2024-T3
7390
0.3
2818
2.824×10-10
Model A
42.9
Cover
CFRP
7143
0.3
2747
3
1.531×10-10
Frame‐
Model B
A2024-T3
7390
0.3
2818
66.1
2.824×10-10
Cover
3

43

解析モデルにおいては，胴体部構造に用いる材料の違いにより 2 種類もモデルを作成した．Model
A においては，フレーム材に「A2024-T3」を適用し，胴体表面カバー材に「CFRP」を適用した．
Mdel B においては，フレーム材とカバー材の両方に「A2024-T3」を適用した．また搭載機器の有無についても比較を行うため，Model A,B それぞれにおいて搭載機器がある場合とない場合についてもモデル化を行った．つまり計 4 種類の構造モデルを用いて静荷重解析と固有振動数解析を行った．本解析に用いたオオワシ 2 号機胴体部構造の 3D CAD モデルを図１に示す．また Mdel
A と Model A に適用した構造材の材料物性値とその構造材を用いた場合の胴体部重量を表１に示
す．
また搭載機器の有無においては，種搭載機器をそれぞれモデル化し搭載機器モデルを実際のように胴体内に配置するのではなく，搭載機器が配置されるリングフレーム部の重量を搭載機器分増加（搭載部リングフレームの密度を増加）することで，仮想的に搭載機器の重量を模擬してモデル化した．表２に各種搭載機器の重量を，また図２に搭載機器の配置イメージを示す．
表２ 各種搭載機器の重量と搭載機器総重量

Component
GG-ATR Engine・Intake
LOX Tank
Ethamol Tank
GN2 Tank No.1
GN2 Tank No.2
GN2 Tank No.3
GN2 Tank No.4
Nose Landing Gear
Main Landing Gear

Total Weight[kg] Component
41.8 Inertial navigation equipment
25.21 GNC Circuit
57.12 TTC Circuit
5.21 Wireless Communication Device
5.21 Battery
5.21 Valve・Plumbing 1
5.21 Valve・Plumbing 2
10.54 Valve・Plumbing 3
34.79 Components' Total Weight

図２ 搭載機器の配置イメージ（側面図）

44

Total Weight[kg]
0.5
0.5
1
1
0.3
5
5
30
233.6

３．静荷重解析結果と固有振動数解析結果
まず静荷重解析の結果を示す．静荷重解析は，機体に 1G が作用した場合について計算した．
また弾性変形範囲内と仮定し，1G 作用時の結果を 9 倍することで 9G の結果を求めた．拘束条件
は，胴体構造部と翼構造部を結合する 8 か所のボルト穴が完全拘束されているという条件で解析した．図３に搭載機器がない場合で 1G の荷重が作用した場合の Model A の変位を示し，図４に搭載機器がある場合で 1G の荷重が作用した場合の Model A の変位を示す．図中で水色の矢印で示された 8 点が拘束点である．
図３，４の比較からも搭載機器があることにより変位が増加していることが分かるが，それぞれの解析より得られた主ひずみ，ミーゼスひずみ，主応力と主ひずみを表３，４にまとめる．表
３に搭載機器がない場合，表４に搭載機器がある場合を示す．それぞれの場合において，応力，
ひずみが最大となった場所は拘束点まわり，つまり翼胴結合部で生じる結果となった．

[mm]

図３ 1G 静荷重作用時の Model A の変位（搭載機器なしの場合）[mm]

図４ 1G 静荷重作用時の Model A の変位（搭載機器ありの場合）

45

表３ 主応力，主ひずみ，ミーゼス応力，ミーゼスひずみ値（搭載機器なしの場合）
Gravitational Acceleration

Major Strain
εMa

1G

1.90×10-4Model A

9G
1G

Model B

9G

Mises Strain Major Stress Mises Stress
εMi
σMa
σ Mi
1.87×10

-4

1.43

1.58

1.71×10-3

1.68×10

-3

12.9

14.22

2.96×10-4

2.89×10

-4

2.23

2.47

2.66×10-3

2.60×10

-3

20.1

22.23

表４ 主応力，主ひずみ，ミーゼス応力，ミーゼスひずみ値（搭載機器ありの場合）
Gravitational Acceleration

Major Strain
εMa

1G

3.81×10-3

3.43×10-2

4.14×10-3

3.73×10-2

Model A
Model B

9G
1G
9G

Mises Strain Major Stress Mises Stress
εMi
σMa
σ Mi
3.39×10

-3

29.9

28.9

3.05×10

-2

269.1

260.1

3.68×10

-3

32.5

31.4

3.31×10

-3

292.5

292.5

搭載機器なしの場合においては，Model A, B ともに使用構造材の限界ひずみを超えていないた
め，破壊が起こらないと考えられる．しかし，搭載機器ありの場合においては，Model A, B ともに翼胴結合穴で構造材の限界ひずみを超える結果となった．以上の解析結果より，今後の構造設計において翼胴結合部の板厚の増加や翼胴取り付け箇所の増加など，応力緩和の対策が必要であることが明らかになった．
次に固有振動数解析の結果を示す．図５に剛体モード以外の最低次モードとして得られた機体
ピッチ軸方向の曲げ振動モードを示す．また表４に搭載機器なしの場合の振動モードと固有振動数，表５に搭載機器ありの場合の振動モードと固有振動数を示す．

図５ ピッチ軸曲げ振動モード46

表 4 振動モードと固有振動数（搭載機器なし） 表 5 振動モードと固有振動数（搭載機器あり）
Mode
1
2
3
4
5
6
7
8
9
10
11
12

Natural Frequency [Hz]Model A
9.1330E-04
7.5243E-04
6.2235E-04
5.3718E-04
3.7544E-04
1.8273E-04
54.87961.727119.12165.86167.56237.26Mode

Mode Shape

Model B
8.9067E-04
6.2585E-04
4.3271E-04
Rigid Body Mode
3.0016E-04
2.3843E-04
4.4120E-04
44.519Airclaft Pitch Axis Bending
50.908
Airclaft Yaw Axis Bending
92.017 Vertical Tail Yaw Axis Bending
134.05
Airclaft Pitch Axis Bending
137.52
Airclaft Yaw Axis Bending
191.87
Airclaft Twist

1
2
3
4
5
6
7
8
9
10
11
12

Natural Frequency [Hz]Model A
3.5319E-04
3.2389E-04
2.5239E-04
1.2003E-04
7.2090E-05
2.0431E-04
23.41624.32264.88671.42686.650112.12Mode Shape

Model B
4.0027E-04
2.4424E-04
2.2622E-04
Rigid Body Mode
5.9475E-05
6.1811E-05
2.1115E-04
22.603
Airclaft Pitch Axis Bending
23.598
Airclaft Yaw Axis Bending
63.187
Airclaft Yaw Axis Bending
68.849Airclaft Pitch Axis Bending
72.000 Vertical Tail Yaw Axis Bending
102.43 Vertical Tail Yaw Axis Bending

搭載機器なしおよび搭載機器ありの Model A, B いずれも，制御系との連成を防ぐための基準として設定した 20 [Hz]を最低次モードで超えているため共振は問題にならないと考えられる．
４．まとめ
胴体部構造の静荷重解析と固有振動数解析において，搭載機器重量の有無によって比較を行った．強度面では翼胴結合部の高強度化が，また振動面では未計上である脚，各種取付金具や艤装の質量を加味した場合の固有振動数低下に備え，固有振動数を高くするための軽量化，高剛性化が今後の構造設計に必要である．

47

小型無人超音速機の翼胴一体型構造による機器搭載性向上対策と１／２スケー
ルモックアップの製作○樋口 健 （もの創造系領域 教授）
勝又 暢久 （もの創造系領域 助教）
小杉 義和 （航空宇宙システム工学専攻 博士前期１年）

１．はじめに室蘭工業大学航空宇宙機システム研究センターで開発中のオオワシ２号機の基本構造概念を再検討した．今後の搭載機器配置の自由度を確保するために，主翼直下の胴体内径が小さくなっている構造を見直し，同一直径の搭載空間を作ることにした．胴体の上に主翼全体を載せる構造様式ではなく，主翼中央部と胴体中央部を一体型とし，その両端に左右主翼をつなげる構造様式とした．併せて，主翼と胴体の結合部で鋭角になっている空間にフィレットを施すことにより空力的な損失を減じるとともに，同一直径の搭載空間を作るための強度上の困難をフィレット部分で克服する考え方を採ることとした．このモックアップを実スケールで製作することは費用の問題があったので，樹脂３Ｄプリンティングによる１／２スケールモデルを製作した．
２．ベースライン形状
図１にベースライン形状を示す．このベースライン形状にしたがって実機大モックアップが昨年度に製作された．

図１．オオワシ２号機ベースライン形状
48

オオワシ２号機の機体構造は，このベースライン形状にしたがって設計される．昨年度までの構造設計では，図２，３に示すように主翼直下が胴体内径に円筒が収められない領域ができる．今後の設計変更や搭載機器の配置調整には大きな制約となるため，胴体の前から後ろまで同一内径の空間を確保できるように考慮する必要がある．また，主翼と胴体の結合部には鋭角に凹んだ流路が機体の左右にできるため，高速飛行時に空力抵抗が大きくなる要因となる懸念がある．

図２．胴体の上に主翼を載せる構造様式

図３．リングフレーム形状（左から胴体前方，右へ胴体後方）

49

３．翼胴一体型構造
翼胴結合部にフィレットを施し，併せて胴体内部搭載容積を拡大する改修を検討した．
フレットを付けることにより，主翼と胴体の結合部の鋭角に凹んだ流路を廃することができるだけでなく，主翼直下の胴体内部に円筒タンクを積めるように改修する際の主翼桁高さを減ずることによる主翼曲げ強度不足を補い，図３のようなリングフレーム内径の凹部をなくする形状を検討した．フィレット形状は，できるだけ単純化するために図５に示す
ように曲面の曲率を１方向とする形状案とした．これにより図６に示す前面形状となる．
なお，ストリンガー（縦通材）を胴体前部から胴体後部まで直線で通せなくなるため，ス
トリンガーの形状とストリンガー同士のつなぎ方には工夫を要する．フィレット内部は，今後の艤装設計により，配管・配線経路として利用したり，電池を分散して配置するための空間として，あるいは舵面駆動のために空間として利用したりすることが可能となる．

(a) 現状構造概念(b) 翼胴一体型構造概念

図４．翼胴一体型構造におけるリングフレームと主翼主桁案

図５．提案したフィレットと翼胴一体型構造（鳥瞰図）50

図６．提案したフィレットと翼胴一体型構造（前立面図）(a) 上面(b) 下面(c) 中央胴体内部

図７．翼胴一体型構造部の 1/2 スケールモックアップ（樹脂３Ｄプリンタ製）
翼胴結合部構造を他サブシステムから見てインターフェースが理解しやすい縮小三次元模型（いわゆる樹脂製３Ｄプリンタによる造形）を製作した．図７(c)に示すように，胴体中央部を貫通するように円筒タンク搭載部が確保でき，機器配置の自由度が大幅に向上した．これにより，フィレットにより高速飛行時の空力抵抗を減らす工夫も考慮されるとともに，図８の機器配置例に示すように，今後の搭載機器配置を見直す際の大きな制約が取り除かれた．

51

図８．オオワシ２号機の搭載機器配置案
４．まとめ
以上の検討の結果，図９に示すようにオオワシ２号機構造概念が固まったので，他サブシステムとのインターフェース調整により次年度は一部実機部品の製作に取り掛かることとする．

図９．オオワシ２号機構造概念

52

GG-ATR エンジンの回転系試験とガスジェネレータ・ラム燃焼器の燃焼、冷却特性について○湊 亮二郎 （もの創造系領域 助教）
渡邉義昭（機械航空創造系学科 4 年）
渡辺翔平（機械航空創造系学科 4 年）
中田大将（航空宇宙機システム研究センター特任助教）
東野和幸（航空宇宙機システム研究センター教授）

１．はじめに室蘭工業大学・航空宇宙機システム研究センターでは大気中を高速度で飛行するための革新的基盤技術に関する研究開発を推進しており、その基盤技術を実際の高速飛行環境で飛行実証するための実験機（フライング・テストベッド FTB）オオワシ 2 号機の研究開発を進めている。同実験機の推進エンジンは，ガスジェネレータサイクル・エアターボラム
ジェットエンジン( Gas-Generator Cycle Air Turbo Ramjet Engine, GG-ATR )を採用している．
そこで本報では，GN2 による GG-ATR エンジンのターボ系要素の回転試験とガスジェネ
レータ(GG)及びラム燃焼器の燃焼・冷却特性について報告する．
２．GG-ATR エンジンの製作と GN2 による回転系冷走試験
2013 年度に GG-ATR エンジンの回転系要素の全部品の製作発注が終わり，エンジンの組
み立てが完成した．このエンジンにはまだ GG とラム燃焼器は取り付けられていないが，まずは GN2 によるターボ系要素の回転冷走試験を実施することを計画している．図 1 に完成した GG-ATR エンジンと GN2 ガス冷走試験系統図を示す．この試験では，回転軸系の振動変位や振動加速度を計測し，安定に定格回転数まで運転できることを実証するのが目
的である．

53

図 1 製作した GG-ATR エンジンと試験系統図
この冷走試験を実施するのに必要なエンジン推力架台と試験設備の整備も順次進めている．
３．GG-ATR エンジン用ガスジェネレータに関する研究
ガスジェネレータの燃焼器とインジェクターの概念設計を行った．図 3 に GG 燃焼器の概観とインジェクター部の断面図を示す．GG の設計点条件は燃焼圧 1.35MPa，燃焼温度1100K である．また GG 燃焼器の体積は，GG 燃焼ガスが安定に燃焼できるように 4.0msecの滞留時間が確保できるように決定し，燃焼器チャンバー構造は GG 燃焼圧力と熱応力に耐えられるように設計した．この GG 燃焼器にはエタノール燃料によるフィルム冷却を行うことを想定してあり，
フィルム冷却を行った場合の GG 燃焼器の内壁壁面温度は 850K 程
度になるものと予想した．

図 3 設計したインジェクター形状案(左)とインジェクター部の断面図（右）54図 4 ラム燃焼器ライナ冷却概念図

図 5 ラム燃焼器ライナ表面の冷却解析

４．GG-ATR エンジン用ラム燃焼器に関する研究
GG-ATR エンジンでは，圧縮機で圧縮された空気とタービンを駆動した GG 燃焼ガスがラム燃焼器で混合・燃焼し，推力を発生する仕組みになっている．このラム燃焼器の燃焼温
度は 2300K 近くまで上昇するので，ラム燃焼器のライナには効果的な冷却システムが必要になる．
GG-ATR エンジンのライナはアフターバーナのような二重ライナ構造にして，二つのライ
ナの間に冷却空気を流して，内側ライナの燃焼ガス側表面をフィルム冷却することを検討している．図 4 はラム燃焼器ライナ壁の冷却システムの概念図を示しており，図 5 はライ
ナ表面をフィルム冷却した場合の，壁面温度の冷却解析した結果である．
これらの冷却解析した結果，ラム燃焼器の内側ライナの材質には，SiC-SiC を用いると耐熱性を確保できると考えられている．
５．まとめ
小型無人超音速機用 GG-ATR エンジンの製作を進め，ターボ回転系要素組立を完成させた．また燃焼器要素（GG 燃焼器，ラム燃焼器）の設計と冷却効果の解析を行った．今後はGN2 によるターボ系要素の冷走試験を行い，定格回転数まで安定に運転できることを実証することを目指す．更に，今回設計検討した燃焼器要素を製作し，タービンをエタノールと LOX の燃焼ガスによる熱走試験を実施し，GG-ATR エンジンの推力や Isp などが所定の特性を確保していることを確認していく．

55

GG-ATR エンジン用超音速インテークの基礎研究○湊 亮二郎 （もの創造系領域 助教）

１．はじめに室蘭工業大学・航空宇宙機システム研究センターでは大気中を高速度で飛行するための革新的基盤技術に関する研究開発を推進しており、その基盤技術を実際の高速飛行環境で飛行実証するための実験機（フライング・テストベッド FTB）オオワシ 2 号機の研究開発を進めている。同実験機の推進エンジンはガスジェネレータサイクル・エアターボラムジェットエンジン( Gas Generator Cycle Air Turbo Ramjet Engine )を採用しているが，同エンジンに取り付けられたインテークで，効率良く空気を取り込めるかが，この実験機の超音速飛行の成否の鍵を握っている．
そこで本研究では，超音速飛行を想定してオオワシ 2 号機搭載用インテーク・ダクトの形状設計とその解析を行った．また超音速風洞でのインテーク・ダクト性能試験を行うことを想定し，その予備的風洞試験も実施したので，それを報告する．

図 1 オオワシ 2 号機用インテークの概念設計案
２．オオワシ 2 号機搭載用インテークの設計
オオワシ 2 号機搭載用の超音速インテーク・ダクトの設計と解析を行った．図 1 に現在，検討中の超音速インテーク・ダクトの形状を示す．このインテークの設計点マッハ数は 1.3である． 図 2 にはこのインテーク・ダクトモデルの CFD 解析結果を示した．CFD 結果では，ダクト内の圧力回復係数が 95%以上を達成できた．今後は機体構造との兼ね合いから，更にインテーク・ダクト形状の改善を図り，超音速風洞試験によりインテーク・ダクトの形状を決定させる．

56

３．インテーク側壁に関する超音速風洞空力試験
前節で設計したインテーク・ダクトモデルは超音速風洞試験にて，空力性能を評価することになる．その風洞試験の前段階として，インテーク側壁形状の違いよる空力性能評価
を本学中型超音速風洞で実施した．この風洞試験の目的は,超音速インテークの側壁が図 3に示されているように，Classical タイプと Inverted タイプについて，性能にどのような違いが現れるか，風洞試験で評価するものである．風洞試験の結果は，オオワシ 2 号機搭載用インテークに反映させる予定である．

図 2 オオワシ 2 号機用インテークの CFD 解析（マッハ数分布）

図 3 Classical タイプ（左）と Inverted タイプ（右）の側壁形状を持つ超音速インテーク
４．まとめ
小型無人超音速機用 GG-ATR エンジンの超音速インテークダクトの設計と，その CFD
解析及び風洞試験を行った．現時点では，インテーク・ダクトの空力性能は M=1.3 条件で圧力回復率 95%以上を達成できた．今後も CFD 解析と風洞試験の両面から，空力性能を評価して行き，オオワシ 2 号機搭載用インテークの形状を決定する予定である．

57

アルミー水反応の衛星推進系への適用

東野 和幸（航空宇宙機システム研究センター 教授）

杉岡 正敏（航空宇宙機システム研究センター 特任教授）○小野寺 英之 （機械航空創造系学科 4 年）
増田 井出男 （JAXA 研究開発本部）

１．はじめに
現在，宇宙機推進システムの燃料としてヒドラジンが使用されている．ヒドラジンは高比推力で貯蔵性に優れた燃料である．しかし，ヒドラジンは人体に極めて有毒である．そこで本学ではAl と水の反応から得られる水素を宇宙機推進システムに適用する研究を行っている．水素は化学
推進薬の中でも最も比推力が高く，人体に無害である．しかし，従来の水素製造法である（1），（2）式の水蒸気改質法では，化石燃料が必要であり，地球温暖化の原因である CO2 を副生する．
CnHm + n H2O → n CO + (m/2 + n) H2（1）

CO + H2O → CO2 + H2（2）

また，水素は金属中に吸収され金属タンクの結晶格子内で金属水素化物になり，金属材料の強度を低下させる（水素脆性）．さらに水素分子は非常に小さく，水素漏れにより貯蔵性が低い．これらのため，長期ミッションが目的である人工衛星への適応性は低い．しかし，(3)式に示す Alと水による水素製造方法では以下の利点が存在する．
2Al+6H2O→2Al(OH)3+3H2（3）（3）式で得られる生成物は水素と水酸化アルミニウムのみで，どちらも人体に無害であり，環境への負担も小さい．また人工衛星での噴射に必要な水素のみを随時製造することができるため，
長期間水素を貯蔵する必要がない．
これまで本学では Al 合金である Al-40%Sn-10%Bi を用いることで短時間で高圧水素を製造することに成功している 1)．しかし，重量検討の結果，ヒドラジン系推進システムと比較して 1.2
倍重量が増加することが判明した．さらに反応温度を 90[℃]に設定する必要があり，電力を多く必要とする．また，宇宙機のエンジンは何度も始動する必要があるが，Al/水反応を停止した後，再度水素製造が可能か確認していない．
そこで，本研究では，軽量で反応温度が低く，反応の制御が可能である Al と Al 合金粉末の開発，反応条件の解明，およびアルミ－水系の比推力の測定を目的とした実験を行った．また，ア
ルミ－水系と他の推薬との比較を行った．
２．実験概要２－１．実験装置
本研究での水素製造は金属製高圧反応容器（オートクレーブ）を用いた．ヒータで実験温度まで加熱し，水素製造量を水上置換で測定した．Al 合金の作製にはセラミック製管状
電気炉を用いた．金属試料をアルミナ製タンマン管に入れ，(a)オートクレーブ (b)捕集管
図１ 実験装置58

電気炉を使用して溶解させて Al 合金を作製した．
表１ 実験条件

２－２．実験方法
表１に①加熱実験，②粒径変化
実験，③Al-Sn 合金と Al-Zn 合金のAl-Sn,Al-Z 合金の実験項目

粒径変化実験

加熱実験
水素製造実験

水素製造実験の条件を示す．①加熱実験では水と純 Al を実験温度ま

試料粉末重量 [g]5で上昇させ，水素を発生させる．

水量 [ml]80

30 分間水素を発生させた後ヒータを止め，30[℃]まで自然冷却させ1350
撹拌速度 [rpm]

1350

1350
なし

る．温度が 30[℃]まで低下後，再び実験温度まで上昇させる．これ

30,50,60
温度 [℃]

60

60

を 4 回繰り返す．90[℃]の場合で70,80,90は，無撹拌の実験も行う．②粒径変化実験では，水と純 Al を実験温開始，停止を
実験時間 [min]

360

360
4 回繰り返す度まで上昇させる．実験温度まで上昇した後に，温度を一定に保ち，実験時間に達するまで水素製造量を計測する．③Al-Sn 合金と Al-Zn

粒径 [μm]

500～0

試料の組成

150～0

100～0

Al-5～40%Sn
純 Aｌ[質量パーセント]純 Al
Al-5～40%Zn

合金の水素製造実験では，過去に優れた結果を残した Sn1）と，水素製造量を増加させた Zn を選定した．実験では，水と Al 合金
を実験温度まで上昇させる．実験温度まで上昇した後に，温度を一定に保ち，実験時間に達す
るまで水素製造量を計測する．
３． 実験結果と考察
３－１ 加熱実験
図 2 の (a) ～ (e) に 撹 拌 速 度
1350[rpm]加熱実験の結果を示す．温
度が高いほど制御性が向上する傾向があるが水素製造量は減少する傾向があった．
以上より，反応させる温度が高い場
合は，人工衛星での噴射に必要な水素

図 230 分ごとに温度を変化させた時の水素製造

量を製造する際に適していることが判明した．しかし反応温度が低い場合には，水素製造量をコ
ントロールすることは困難なので人工衛星のエンジン操作に適していないと考えられる．高温で制御が可能になる理由として，高温の水中で生成される水和酸化皮膜(Al2O3・H2O，β-Al2O3・
3H2O 等)が反応を妨げているためであると考えられる．
図 2 の(f)に撹拌速度 0[rpm]，90[℃]の場合を示す．無撹拌の場合においても，制御性があることが判明した．しかし，撹拌がある場合より水素製造量が大幅に減少することが判明した．
水素製造量が減少した原因としては， Al 粒子が堆積し，隣接する Al 粒子どうしが密着するた59

め，水との接触面積が小さくなったためであると考えられる．
３－２ 粒径変化実験の結果
図 3 に粒径変化実験の結果を示
す．粒径が細かいほど水素製造量は大幅に増加することが判明した．実験開始から 360[min]後における Al(粒径：20～0[μm])の水素製造効率(計測値を理論値で割った値)は 77.8%である．また，Al(粒径：20～0[μm])は，実験開始から 2270[min]後に 6650[ml]の水

図 3 種々の粒径の Al 粉末による水素製造(60[℃])

素を計測した．水素製造効率は 99.5%であった．これは，粒径が小さくなるに従い，水との接触面積が大きくなり，水素製造量が増加するためと考えられる．
３－３ Al-Sn 合金と Al-Zn 合金による水素製造実験
図 4 に Al-Sn 合金を示し，図 5 に Al-Zn 合金の水素製造結果を示す．どちらの合金においても添加金属が多いほど水素製造量が多くなる傾向がある．しかし，Al-5%Sn では添加金属が少ないにも関わらず水素製造量が大きいことが判明した．また，全ての Al-Sn 合金において，常温で水と接触するだけで水素が発生することを目視観察で確認した．

図 4 Al-Sn 合金による水素製造

図 5 Al-Zn 合金による水素製造

一方，Al-Zn 合金は水と常温で接触しただけでは水素は発生しなかったが，常温より高い温度では水素が発生する可能性がある．
４. 衛星等への適用検討
４－1

純アルミ－水反応方式と他の方式との比較

質量 50kg 程度の小型衛星の推進系に純アルミ－水反応を利用した際と他の方式（窒素ガ
スジェット，ヒドラジンガスジェット等）を適用した場合のシステムトレードオフ検討を行った．純アルミ－水反応を利用した方式は推進系として搭載質量的に成立する可能性があり，他の方式に比べ「完全無毒」や「高圧ガス非該当」等，安全上有利であることが判明した．
60

４－２ 噴射試験
図６に噴射試験の結果，図７に実験装置の詳細を

表２ 噴射試験の結果示す．水素ガスの噴射試験を行い，発生水素ガスのみで理論上に近い比推力が得られる可能性が示された．また，水素ガスに発生液を加えて噴射試験を実施したところ，ノズル（スロート径 0.5[mm]）は氷等で閉塞せず，比推力が増加（水素＋水では 639.4[s]，水素＋30[％]実液では 523.7[s]）することが判明した．
発生液が推力増加に寄与する可能性がある．
図 6 噴射試験の結果

図 7 実験装置
５．結言
本研究により明らかになった結果を以下に示す．
(1) 加熱実験：純 Al では，反応させる温度が 70～90[℃]の場合は水素製造の制御は容易であり，宇宙機の推進系に適していると考えられる．
(2) 粒径変化実験：粒径が 20～0[μm]の場合，水素製造量は大幅に増えることが判明した．
(3) Al-Sn 合金と Al-Zn 合金の水素製造実験：Al-Sn 合金は常温で水と反応して水素を発生させるが，Al-Zn 合金においても常温以上では水素が発生する可能性がある．
(4)小型衛星に本反応を利用した際と他推薬（窒素，ヒドラジン等）と比較を行い，成立の可能性があること及び他の方式に比べて有利なことを示した．
(5)水素ガスに発生液を追加すると性能が向上することが試験で判明した．残物質がデッドマスにならない可能性が示された．
今後の課題を以下に示す．
①Al 合金推進系の制御方法を確立．
本実験では純 Al の制御方法を確立したが，Al 合金の制御方法は確立していない．水素製造を制御できる Al 合金を開発する．
②Al-Sn 合金以外の常温でも水と反応する Al 合金粉末の開発．
Al-Sn 合金のような添加金属が少なくても常温で水と反応する Al 合金を開発し，宇宙機推進系に適用する．
③水素製造能力の高い Al-5[%]Zn-5[%]X，Al-5[%]Sn-5[%]X (X=In，Ga，Sb，Ni 等)を開発．
Al-5[%]Zn 合金，Al-5[%]Sn 合金に他の合金を加え，さらに水素製造能力の高い合金開発を行う．

61

④アルミ－水反応の現象の詳細を把握及び最適条件の調査を行う．
⑤アルミ－水反応の生成物の物性評価及び分析を行う．
⑥微小重力下での固気液分離装置を考案し，発生した水素を固液から分離する．
⑦本推進系のシステム検討を継続する．
参考文献
1) 近藤光輝，東野和幸，杉岡正敏：宇宙機推進システムとしての Al と水との反応を利用した常圧および高圧水素製造に関する研究，第 53 回航空原動機・宇宙推進講演会要旨集（2013）62

バイオエタノールの熱分解吸熱反応と触媒効果について
東野 和幸（航空宇宙機システム研究センター 教授）

杉岡 正敏（航空宇宙機システム研究センター 特任教授）

塚野 徹（航空宇宙システム工学専攻 博士後期 3 年）

山本 康平

（航空宇宙システム工学専攻 博士前期 2 年）

○ 飯島 明日香 （航空宇宙システム工学専攻 博士前期 1 年）
笹木 康平

（航空宇宙システム工学専攻 博士前期 1 年）

1． はじめに宇宙開発が重要視される中，今後宇宙輸送機の打ち上げ回数の増加が見込まれ，宇宙輸送機が環境に与える影響が無視できなくなる．そこで，再生可能な自然エネルギーである燃料としてバイオエタノールが注目されている．しかし，バイオエタノールはロケ
ットエンジンの燃料として用いられた経験がないため，実用化のためには燃焼特性や材
料適合性，冷却特性などの基礎特性を解明する必要がある．このうち冷却特性について，
バイオエタノールはアルコール燃料であることから，高温環境下で吸熱を伴う熱分解反応を生じる可能性がある．この化学的な吸熱反応をロケットエンジンの冷却に利用でき
れば，バイオエタノールの冷却能力向上に寄与する．そのため，熱分解による吸熱を含めた冷却能力を把握することはバイオエタノールロケットエンジンを開発する上で重要である．
昨年度までに行った研究より，準静的環境下でバイオエタノールを流通させた場合，適切な金属触媒を用いることで熱分解開始温度が低下し熱分解反応が促進されることが判明した．また，触媒により特定の反応を促進していることを確認し，触媒の有用性を確認した．本研究では流量を増加させた動的環境下での加熱流通実験を行い，バイオエ
タノールの吸熱効果について評価した．また，触媒を使用する際の設計指針を把握するため，触媒付近での滞留時間と流れの状態が熱分解吸熱反応や冷却能力に及ぼす影響を評価した．
2． 実験概要
動的環境下の加熱実験装置の概要を図 1 に示す．本実験ではバイオエタノールを自己加圧により流動させる．触媒を使用しない場合は電気炉，触媒を使用する場合は触媒リアクタを加熱しバイオエタノールを熱分解させる．表 1 に動的環境下の加熱実験の実験条件を示す．バイオエタノールの流量は約 1.0~1.5 g/s，実験圧力は 0.3 MPaA とし，実験
時間は加熱部のバイオエタノールの温度変化が安定的になるまでの約 5 min とした．

63

表 1 実験条件
実験目的
実験種別
実験温度

触媒を使用しない場合
等温～約950 K

触媒種類

使用しない

反応容器径

-

触媒を使用する場合
等温～約950 Kα -アルミナ
γ -アルミナ
白金アルミナφ 81

触媒反応容器の形状を変化させた場合
等温～約950 Kα -アルミナ
白金アルミナ
φ 81，φ 58
V7

触媒リアクタサンプリングボトル

PCRIN

V3

PSHIN

TCRH

TSHIN

TCR-1
THTOUT
PHTOUT

V9

V8

TCR-3

差圧計 PSHOUT

THTIN
PHTIN
PCROUT

V5

V6

オリフィス

R1

V2

TCR-2

加熱管

TCROUT

POROUTPSHOUT

TOROUTTSHOUT

TIW-1~9
TOW-1~9

ドレンポート

気密試験用プラグ

POROUT
オリフィス

TLIQ

R3

V4

TVAP

PTANK

V1

R2

供試液
シリコンオイル

電気炉
マイクロケーブルエアヒーター

投げ込みヒーター
GN2ボンベ

図 1 バイオエタノール加熱流通実験装置
3． 実験結果
3-1．触媒を使用しない場合の吸熱量
図 2 に触媒を使用せずバイオエタノールを熱分解させた実験結果を示す．GN2 流通時は投入熱量の増加に従って電気炉出口のエンタルピーは増加したが，バイオエタノール流
通時は熱分解吸熱反応によりエンタルピー上昇が抑えられた．
3-2．触媒を使用した場合の吸熱量
図 3 に触媒を使用してバイオエタノールを熱分解させた実験結果を示す．図 3 では，実験の設定温度として実験開始直後の触媒付近ガス温度の最大値を横軸にとり，縦軸は触媒付近のガス温度が安定した後とのエンタルピー差を示している．γ-アルミナおよび
白金アルミナを使用した場合は，触媒効果の無いα-アルミナと比較して，エンタルピー
低下が大きくなり，触媒効果による熱分解吸熱反応の促進が確認された．
3-3．触媒反応容器の形状を変化させた場合の吸熱量
図 4 に触媒反応容器の形状を変化させてバイオエタノールを熱分解させた実験結果を示す．実験の設定温度として実験開始直後の触媒付近ガス温度の最大値を横軸にとり，縦軸は単位時間あたりに低下した触媒付近ガス温度の最大値を示している．滞留時間は冷却能力へ与える影響が小さく，熱伝達率が大きいほど冷却能力が向上することが確認
された．

64

2000
触媒付近ガス温度TCR-3の最大値と安定化のエンタルピー差[kJ/kg]

300

BEh(Tout)-h(Tin)[kJ/kg]

1500

GN2

1000
807K
515K 740K
680K

500

929K

932K

640K

656K
928K
825K

374K
0

567K

864K

720K

-500
0

500

1000

1500

2000

100
-100
-300
-500
α-アルミナ300g-700

γ-アルミナ300g-900

白金アルミナ300g

-1100
300

500

700

900

1100

触媒付近ガス温度TCR-3の最大値[K]⊿qinput[kJ/kg]
単位質量当たりの加熱量と比エンタルピー上昇の関係
図 2 電気炉出入口のガスエンタルピー変化

図 3 触媒付近のガスエンタルピー変化

図 4 触媒付近ガス最大温度低下の変化
4． 結言
本研究では，バイオエタノールの吸熱効果について評価すること，触媒を使用する際の設計指針を把握することを目的としてバイオエタノールの加熱流通実験を行った．実験の結果，バイオエタノールの熱分解時にエンタルピー低下が起きる事を初めて確認した．また，触媒を使用した場合，比較的低温におけるエンタルピー低下が大きくなる事
や，触媒付近での熱伝達率が大きいほど冷却能力の向上が見られることを確認した．今後は，高圧条件（10 MPaA 程度）でのバイオエタノールの熱分解による化学的な吸熱量の定量評価を行う予定である．

65

ＡＴＲ－ＧＧ推進剤供給系の検討
東野 和幸 （航空宇宙システム研究センター 教授）
中田 大将 （航空宇宙システム研究センター 特任助教）○林 祐一郎 （機械航空創造系学科

4 年）

１． はじめに
室蘭工業大学航空宇宙機システム研究センター(APReC)では，小型無人超音速実験機オオワシ 2号機の研究・開発を通じて次世代クリーン宇宙輸送システムの革新的基盤技術の確立を目指している．本機にはバイオエタノール(BE)／液体酸素(LOX)を組み合わせた推進システム使用が検討されている．BE は表 1-1 で分かるように液体水素(LH2)と比べて密度比推力が優れている．そのため，燃料タンクを小型にでき，機体重量の削減が期待できる．また，燃料自体のコストが安く，今後のさらなる値段の低下も期待出来る．単位質量あたりの炭素含有量がケロシンと比較して少ないので煤の発生量も少なく整備コストの低減に寄与することが出来ると考えられる．また，常温で液体である BE は機器への負担も小さく部品コストも抑えることが出来る．
表 1-1 推進薬の製造コストの比較
比推力

密度比推力
3

整備性

3

コスト

安全性

(sec)

×10 (kg・s/m )(円/kg)

LH2

455

145

○

4000

△

BE

328

328

○

215

○

ケロシン

368

368

△

450

○

ヒドラジン

397

397

△

8800

△

本機体の推進剤供給方式としては，部品点数を減らせることにより故障の可能性を減らすことを狙い加圧供給方式とした．

図 1-1 GG 燃焼圧力と推進剤流量に関する変化

66

本稿では，オオワシ 2 号機に搭載されるエアターボラムジェット・ガスジェネレータサイクル(ATR-GG)エンジンにおける推進剤供給系において軽量で耐腐食性のある推進剤タンクについて検討する．また，本機は推進剤流量を変化させることで図 1-1 のようにガスジェネレータ(GG)の燃焼圧力を変化させ，推力のスロットリングを行う．連続的な流量調整弁に代わって，複数のバルブの組み合わせで推進剤流量を変化させる機構を検討しているが，その選定候補についても述べる．
２． 推進剤供給系検討
本機体の推進剤供給方式としては，部品点数を減らせことにより故障の原因を減らせるため，加圧ガスを用いた加圧供給方式を採用した．また，今回は加圧ガスにガス窒素(GN2)を用いる．
推進剤タンクに関しては，通常ではアルミニウム合金を使用することが多いが，本機の推進剤である BE はアルミニウム合金に対して強い腐食性を有しているため，過去の BE に対する材料適合性の研究をもとに炭素繊維強化プラスチック(CFRP)とニッケルを組み合わせることで BE に対する耐腐食性を有するだけでなく，アルミニウム合金でタンクを製作するよりも軽量に仕上げることを狙い検討対象とした (1)．
内圧 2.6MPaA，安全率 2.0 とした場合の燃料タンク重量を比較した結果が表 2-1 である(2)．以下の結果より CFRP と Ni メッキを組み合わせたものが最も軽量になることを確認した．
表 2-1

タンク重量比較

アルミニウム合金

21.3kgステンレス合金

61.7kg

CFRP+Ni メッキ

11.5kg３． 加圧ガス混入対策
オオワシ 2 号機には，推進剤を約 130kg 搭載する予定であるので，それに伴い推進剤タンクも約 3m と非常に長くなることが予想されている．これらの事により推進剤減少時に GN2 も供給されてエンジンに影響を与えることが懸念される．そのため，下記の 3 つの案を提案しトレードオ
フを行なった．
上記のトレードオフより 3 つの対策案の中では案(3)が適している．

図 3-1

案(1)

図 3-2

案(2)

67

図 3-3

案(3)

表 3-2

利点
対策案(1)

欠点
利点

対策案(2)対策案(3)

欠点
利点
欠点

各案のトレードオフ

軽量なディバイス．
混入防止効果が小さい．
推進剤の無駄が多い．
推進剤減少時の加速度の変化に対応出来る．
柔軟性のある配管の製造が難しい．
重りがタンク内部を損傷させる危険性がある．
タンク内部に入れることが難しい．
液体減少時でも液面を高く保つことが出来る．
スロッシング抑制効果も期待できる．
タンク内部に構造物を入れることが難しい．

４． バルブ選定
本機は推力を制御することで離陸から着陸まで行うため飛行状態にあわせて推進剤流量を変化させる機構として流量調節バルブの選定を行なった．表 4-1 の選定条件をもとにマッハ数変動
と共に変動するバルブ容量係数(CV)を図 4-1 に示す．

表 4-1 選定を使用した検討諸元

流体の種類
種類

BE

LOX

密度 (kg/m3)
温度(℃)
バルブ前後差圧
O/F
バルブ接続口径 (in)

789

1140
-183

図 4-1

15
30
0.454
1/2

マッハ数変動による CV 変化68本年度の検討諸元を用いて選定をおこなった結果としては，LOX 側に関しては，適合する既存のバルブを確認することが出来たが BE 側に関してはさらなる検討が必要であることが判明した．
５． まとめ推進剤タンクの材質を CFRP にすることにより従来のタンクより軽量化が可能であるという結果が得られた．
バルブ選定を行った結果，BE 側に関しては継続して検討を行う．
参考文献(1) 寺田利幸：バイオエタノールロケットエンジンの高温高圧での材料適合性に関する実験的研究，平成 23 年度修士論文，(2012)．
(2) JAXA：宇宙用高圧ガス機器技術基準，(2009)．

69

高圧，高温環境下におけるバイオエタノールの熱分解吸熱反応，サルファアタック・
コーキングの実験装置について
東野 和幸（航空宇宙機システム研究センター 教授）

杉岡 正敏（航空宇宙機システム研究センター 特任教授）

塚野 徹（航空宇宙システム工学専攻 博士後期 3 年）

山本 康平

（航空宇宙システム工学専攻 博士前期 2 年）

○ 笹木 康平（航空宇宙システム工学専攻 博士前期 1 年）

飯島 明日香 （航空宇宙システム工学専攻 博士前期 1 年）

1．はじめに
宇宙利用の活発化に向けて，環境適合性や再使用性に優れた新しい燃料としてバイオエタノール（BE）が注目されている．バイオエタノールは常温で液体であり，取り扱いが容易かつ毒性も有しておらず管理および運用コストを低減させられる可能性がある．また，植物等のバイオマスを原料として製造されているためカーボンニュートラルであることから，環境へ配慮した燃料である．
しかし，バイオエタノールを燃料としたエンジンシステムの開発実績は少なく，実用化のためには
燃焼特性，冷却特性，材料適合性等の基礎特性を解明する必要がある．
このうち冷却特性について，バイオエタノールはアルコール燃料である事から，高温にさらされた際，熱分解吸熱反応が進行し，その吸熱を機体やエンジンの冷却に利用することで冷却能力を向上できる可能性がある．さらに，適切な触媒により，吸熱量及び反応温度を制御できる可能性がある．
しかし，BE を冷却材として用いた場合，ロケットの高性能化（高燃焼圧化）に伴う冷却材の高温化によって，再生冷却溝表面へのコーキング（熱分解に伴う煤の析出）が生じ，冷却材と再生冷
却溝表面間の熱伝達率の低下や再生冷却流路の下流に位置する推進剤噴射器流路の閉塞を招く恐れがある．そこで，BE を対象とし，実機環境を模した熱分解吸熱反応の詳細を解明することが重要である．
一方で，BE は植物由来であるために 1[ppm]程度の有機硫黄化合物を含有する．そのため，含有硫黄成分によりロケットエンジンの再生冷却溝材料に対する腐食（サルファアタック）が生じる恐れがある．再生冷却溝にサルファアタック（含有硫黄成分による腐食）が発生した場合，その金属材料の構造強度を低下させる．これらが原因となって再使用性（耐久性）や信頼性が低下する．そこで，サルファアタックについても実機環境を模した実験により詳細を解明することが重要である．

図 1 コーキング進行メカニズム

図 2 サルファアタック進行メカニズム70

2．実験装置，実験内容
実機環境を模した環境において，熱分解吸熱反応・コーキング・サルファアタックの特性を詳細に解明するために図３のような装置を検討，製作している．
・熱分解吸熱反応
無触媒･各種触媒を使用した BE の加熱流通試験を行い，供試流体の温度変化の取得および試験後の BE のサンプリングを行う．そして，高温高圧の動的環境における，反応温度に対する吸熱量などの熱分解吸熱特性を検証する．
・サルファアタック，コーキング
実機における使用実績のある金属材料を配管とする供試体を製作する．また，サルファアタ
ック・コーキング形成，剥離の重要なパラメータとしては壁面温度，せん断力が挙げられる．しかし，サルファアタックの原理から，せん断力が効くのは堆積した硫化金属，炭素の剥離に対してだと考えるため，表面に硫化金属，炭素が形成されるかの確認から行う．
表１ サルファアタック・コーキング分析項目

図 3 高温・高圧実験装置 イメージ図

713．実験装置製作の進捗状況
現在，実験の実施に向けて装置を製作中である．図 4 は ＢＥタンクの製作途中における写真である．なお，本装置は JAXA 殿との共同研究としての位置づけである．

図4

水タンク・BE タンク

72

反転軸流ファン試験装置○ 中田 大将（航空宇宙機システム研究センター 特任助教）
棚次 亘弘（航空宇宙機システム研究センター 特任教授）

１． 背景
本学では超音速無人実験機オオワシにおける推進系の一候補として反転ファン式超音速用ジェットエンジン（図 1）の採用を想定し CFD 解析などを主体とした設計を行ってきた[1]．１段動翼＝静翼の軸流ファンを用いる場合と比べファン効率の改善が期待できるが２枚のローターの相対速度は超音速となるため通常の軸流ファンに比べ低い回転数で衝撃波による効率の低下が発生し長所を発揮するための適切な設計は難しい[2]．重要なパラメタとしては 1)前段＝後段の回転数比，2)ファン間距離などが挙げられる
が，これらの影響について実験装置を構築し，基礎実験を行う体制を整えた．

図1

反転軸流ファン式超音速ジェットエンジン断面図

図2

反転ファン設計点

73

２．実験装置２－１

羽根主要設計諸元

図２に設計した反転ファンの 3 次元形状を示す．前/後段の羽根数はそれぞれ 15/11 枚であり，ファン間距離は 54 mm でファン径に対する比は 0.35，前段ブレードの翼弦長に対する比は 0.89
である．一般にファン間距離が近すぎる場合，相互の干渉が大きくなり，逆にファン間距離が遠すぎると翼列間での全圧ロスにつながる[3]．今後，ファン間距離を変化させた試験が重要となるが，現段階では固定したファン間距離にて試験する．
２－２

動力装置

動力装置としては外径が 60 mm の電動 RC 機用 11kW 直流ブラシレスモーターを使用した（図3）．真空槽内で使用するため，全体をモールドし，かつモーター周囲を不凍液で水冷している．
電源としては，4000mAh・セル数 14 の LiPo バッテリを４並列で使用している．これにより，定格回転数での運転を連続で 6 分間程度行うことが出来る．
２－３

計測系の確認図４に示す通り，静圧ポート 15 箇所（１段ファン上流，１段動翼上，１－２段間，２段動翼
上，２段出口の 5 か所×周方向 3 箇所），3 孔ピトー管を４か所（１段ファン上流，１－２段間，２段出口，オリフィス直前）に配置した．ピトー管は r 方向およびθ方向にステッピングモータ
ーにより駆動でき，左右の孔に差圧計を取付け，中央の孔に全圧計を取り付けて用いる．差圧計の出力がゼロとなる方向にピトー管を向け，その地点での角度と全圧から流速ベクトルを推算できる．

図3

ファン駆動用モーターとバッテリーボックス

74

図4

計測点一覧

図 5 に示すようにピトー管挿入位置にはシース熱電対も組み込まれているが，温度についてはピ
トー管の熱容量が大きく，静定まで長い時間を要する．これまでの試験では 3 分程度経っても静定する様子が見られない．バッテリーの連続駆動時間は定格回転数で 6 分程度のため，素線がむ
き出しとなっているような熱容量の小さい熱電対を用いることで静定までの時間を短くしたい．

図5

3 孔ピトー管概要75

図6

反転ファンリグ試験機セッティング

図 6 に反転ファン試験機概要を示す．現在までに大気圧環境下で 18000rpm までの回転数で試験を行った．流量は２段ファン出口に径の異なるオリフィスプレートを取り付けて変化させる．図7 はファン１のみを 9600rpm で回転した際の１段上流のピトー管で測定された流速ベクトルの例である．チップとハブの中間に急激な流速ベクトルの変化が認められる．上流側を整流せず，狭い真空容器に入れたことによる影響などが考えられ，今後ベルマウスまたはハニカム板で整流するなどの対策を施して再試験する予定である．
また，各段の回転数を上昇させる過程でインペラの空力負荷が変化することにより，ブラシレス
モータードライバに大きな負荷がかかり，故障する事例が２度発生した．これを防止するには回76

転数上昇をよりゆるやかに行う，パワー変動を吸収するよう電力ラインへのコンデンサを増設するなどの対策が必要と考えられる．

図７ １段ファン上流の流れ場の様子（角度情報のみ）．Hub=Tip の中間に流速ベクトルの急変する箇所がある．
３．まとめ
超音速機エンジンへの適用を目指した反転ファンリグ試験機を整備した．今後，真空槽内で定格回転数までの試験を行い，ファン効率を低下させないようなファン間距離，回転数比などについて知見を得る予定である．

参考文献
1.

Ryojiro Minato, et. al., “Development of Counter Rotating Axial Fan Turbojet Engine for Supersonic Unmanned Plane,”
AIAA2007-5023

2.

Xiao-He Yang, et. al., “Design of Two Counter-Rotating Fan Types and CFD Investigation of Their Aerodynamic
Characteristics” ASME 2011 Turbo Expo: Turbine Technical Conference and Exposition, 20113.

岡部，福富，重光，二重反転型小型軸流ファンの前後段翼列間距離の影響，日本機械学会流体工学部門講演
会論文集, 201077

小型無人航空機向け遠隔監視制御系用無線通信装置の構築○

上羽 正純 （もの創造系領域 教授）
野口 慎太郎（機械航空創造系学科

4 年）

竹内 僚太郎（航空宇宙システム専攻 博士前期 1 年）

1．はじめに室蘭工業大学・航空宇宙機システム研

小型無人実験機究センターでは大気中を高速度で飛行（送受信装置搭載）

するための革新的な基盤技術を創出す
アンテナる研究開発を推進しており，地上で研究

テレメトリー回線（高速）された基盤技術を小規模ながらも実際の高速飛行環境で飛行実証するためのリアルタイム表示

コマンド回線（低速）

無線通信

小型無人超音速実験機（フライング・テ
ストベッド FTB）の構築を進めている．

地上局

本実験機を長距離飛行させるには，実図１ 遠隔監視制御系構成

験機に搭載したセンサにより位置，速度，
姿勢等の飛行状態を計測し，それらデータをリアルタイムで無線伝送する，あるいは地上からのコマンドを実験機に伝える遠隔監視制御系を構成（図１）することが必須である．無線伝送するための装置は，最終的な飛行距離である 100km の伝送できる能力が必要である．本報告では，本遠隔監視制御系実現のために，実験規模に応じて使用可能な無線通信装置を市販製品あるいは特注品により構築を進めている状況を報告する．
2．無線通信装置の構築
最終的な実験機は 100km 飛行するためにそれだけの距離を無線伝送可能な装置が必要である
が，最終的な実験機が完成するまでの途中段階の誘導制御系技術等の実証の観点からは，プロペ
ララジコン機，ジェットエンジンラジコン機等の小規模，中規模サイズの実験機が使用される．
この場合，これら中間段階の実験機では最終的な実験機に比して搭載スペースが小さく，飛行距離も短いため，規模に応じて搭載可能な市販の無線通信機器を使用しつつ，最終的な無線通信装
置を段階的に準備することが現実的である．
飛行実証のための距離を勘案して，1km，10km，100km の 3 段階の無線伝送距離に区分して，
前 2 つの伝送距離の場合には市販の無線通信を，100km の伝送距離においては特注品を製作することとした．これら準備した無線通信装置の概要及び性能を以下に示す．
2.1

1km 用無線通信装置

1km 用無線通信装置として XBee-PRO S1，XBee-PRO S2B， XBee Wi-Fi の 3 種類の市販の無線通信モジュール（図２）を選択し，それらの通信試験を行った（図３）．その結果，最も受信率の高く，500m 以上の伝送距離を有する XBee-PRO S1 を 0～1km 用無線通信装置として選定
した．

78

XBee-PRO S1

XBee-PRO S2B

XBee Wi-Fi

距離(参考値)

750m

1500m

120m

伝送速度

250kbps

250kbps

65Mbps

図 2 対象市販無線通信モジュールと性能

2.2

図３ 無線通信試験結果

10km 用無線通信装置

0～10km 用無線通信装置についても調査の結果，図 4 に示すように市販品である東京コスモス
電機製の無線通信モジュール TWE-STRONG を選定した．本無線通信モジュールは特定小電力製品であり，無線従事者免許不要で操作可能である．本モジュールについて通信試験を行った結果を図５に示す．現時点では，屋内の距離 5m での動作確認の完了にとどまっているため，今後，
指向アンテナを使用することにより，10km 程度の伝送距離の通信試験を行い，性能を確認する．
100.0%
99.5%

外観

伝送距離

99.0%

内臓アンテナ

1.7km

98.5%

無指向性アンテナ

3.1km

八木宇田アンテナ

10km

受信率

使用アンテナ

98.0%

97.5%

250kbps
500kbps
667kbps

97.0%
96.5%

伝送距離

96.0%
1kbytes

図４ TWE-STRONG の外観と伝送性能

屋内試験結果（5m）

10kbytes

100kbytes

1Mbytes

送信データ量図５ TWE-STRONG 通信試験結果
2.3100km 用無線通信装置

100kｍの長距離伝送を可能とする汎用の表１ テレメトリー，コマンド回線性能諸元

無線通信装置は市販されていないため，特注品として製作することとした．この場合，使テレメトリー回線

コマンド回線用周波数を決定した上で，無線通信システム

使用周波数(帯)1281.5MHz*

400MHｚ帯

として送信装置と受信装置をセットで性能

送信装置
アンテナ利得

0dB

0dBi

送信電力

2Wmax

1Wmax

受信装置
アンテナ利得

10 dBi

10dBi

情報データ速度128kbps

10kbpsを決定して，設計・製作を進める必要がある．
遠隔監視制御系に必要な無線通信回線は，
テレメトリー回線とコマンド回線があるた
め，表１に示すようにテレメトリー及びコマ
ンドそれぞれの性能諸元を設計により決定*:中心周波数し，これに基づき送信機，受信機の性能を決79定した．平成 25 年度は，図６に示すような外観，性能を有するテレメトリー回線におけるテレメ
トリー受信装置を製作した．

項目

性能値

最小受信感度-111.3dBm

FEC

Reed Solomon符号化

復調器

MSKあるいはFSKに対応

電力変動対応幅-57 ～―111.3dBm

受信帯域幅

1.5MHz

復調データ出力

シリアルデジタルデータ

サイズ（㎜）

300（Ｌ）×100（Ｗ）×40（Ｈ）(b) 性能諸元(a) 外観図６ 製作テレメトリー受信装置

3．まとめ使用する実験機の飛行距離に応じた無線通信装置を準備し，10km 以下のものについては市販品を選定し，その性能を確認した．100km 用にはテレメトリー回線の受信装置を製作し，所望の性能を有することを確認した．今後は，10km 級無線通信装置については，通信試験により 10km程度までの無線伝送性能の確認を行うとともに，100ｋｍ用無線通信装置については，テレメトリ
ー回線用の送信装置の製作を行うとともにコマンド回線用無線通信装置の基本設計を進める．

80

小型無人超音速機向け誘導制御システムの構築とラジコン機による飛行検証
上羽 正純（もの創造系領域 教授）○ 竹内 僚太郎（航空宇宙システム工学専攻 博士前期 1 年）
山下 智也（航空宇宙システム工学専攻 博士前期 1 年）

1.

はじめに
超音速機をはじめとする大気中を高速・高々度まで飛行する飛翔体実現のための基盤技術を確立するための Flying Test Bed と称する小型無人超音速機を飛行させ検証を行うことが大変重要である．本小型無人超音速機は，完全自律の長距離飛行とならざるを得ない．このため，滑走，
離陸，上昇，旋回，水平飛行，下降，着陸の各種飛行モードを実現するための飛行制御則，セン
サ，制御回路，アクチュエータ，無線通信装置等のハードウエアを開発していく必要がある．
本報告では，飛行制御則の検証ならびに実験機への搭載を視野にいれつつ，ラジコン飛行機を転用した小型自律無人航空機を構築し，水平定常飛行制御則について検証した結果を報告する．
２.

誘導制御システム構成構築した誘導制御システム構成を図１に示す．模型飛行機であるラジコン機にセンサとして

INS/GPS（慣性航法装置），ADS（エアーデータセンサ），試作品の制御回路を搭載し，ラジコン機の電動モーター並びにサーボモータに制御信号を送り，プロペラ回転，舵面角度を駆動する．また，本制御回路は，自律飛行のみならず，ラジコン操縦用プロポからの信号により，手動操縦・自動操縦の切替機能を有しており，特定の飛行モードの制御則の検証が可能である．

サーボモータ(SG-5010)
ch1:エルロン
ch2:エレベータ
ch4:ラダー

無線通信モジュール（XBee PRO S1）

地上PC

パワーアンプ
INS/GPS（VSAS-T1G）

誘導制御回路（試作品）

ADS
５孔ピトー管

モータ

LiPoバッテリー

受信機
ch1:エルロン
ch2:エレベータ
ch3:モータ
ch4:ラダー
ch5:M/A切り替え

圧力センサボード（LBA S500 B E 6 S）

プロポ

図１ 誘導制御システム構成
３. 構成機器の性能
図１に示した構成機器のうち，自律飛行に必要な計測・制御機器について述べる．

81a) INS/GPS
機体の位置，姿勢等の情報を計測する INS/GPS として，東京計器 VSAS-T1G を使用している．このセンサにより取得される姿勢角データは局地水平固定座標系からのずれ角であるため，任意の方角を基準座標系に設定可能なよう，制御回路内で座標変換計算を行っている．この基準座標系からのずれとして姿勢角が出力され，自律飛行時の制御に用いられる．
b) Air Data Sensor (ADS)
圧力センサと５孔ピトー管より ADS を構築した．ADS により機体の対気速度，迎角，横滑り角を計測し，制御器への入力として使用される．
c) 無線通信モジュール
INS/GPS により計測された位置，姿勢データ，搭載ハードウエアの動作状況を送信する．小型で手軽に使用する可能な無線通信モジュールとして XBee PRO S1 を選定した．飛行試験においては誘導制御回路の処理速度の制約により，伝送速度 19.2kbps となっている．
d) 地上端末
機体搭載の無線通信モジュールから地上に送信されたデータをｃ）と同じ無線通信モジュール
により受信し，ソフトウエア（TeraTerm）にてそれらデータが表示・蓄積される．
e) 誘導制御回路
マイコンボードを用いた誘導制御回路では，センサより得られるデータを用いて，制御則計算，
その結果に基づく各舵面への PWM 出力，手動・自動の切り替えを行う．手動モードではプロポより受信機を介して送られた舵角コマンドを，回路をそのまま通過し，サーボモータへ伝える．自動モードではセンサデータを用いて制御則により

表1

物理特性

計算された舵面角をサーボモータへ伝える．制御

質量(kg)

m

2.2

周期は 10ms である．

翼面積(

S

0.364

Ix

0.132

Iy

0.114

Iz

0.219

Ixz

0.006

)

飛行中の Manual/Auto はプロポの特定のチャネ
ル（スイッチ）の操作により行う．誘導制御回路

慣性能率(

)

側でそのチャネルのパルスを常に監視し，プロポのスイッチを操作することによるパルス幅が設定

慣性乗積(

)

しきい値以上，以下になることで自動・手動の切替を行う．

表 2 縦系無次元空力微係数f)使用機体-0.18607

高翼式京商カルマートであり，ブラシレスモー
タにより推進力を得る電動機である．Solidworks
を用いて機体の 3DCAD データを作成し，質量特性より表１に示す物理特性を求め，得られた値から推算式を用いて表２に示す機体の空力微係数を推
算した．

0
̇

-0.75456-9.25670-15.968

0.009618-5.53557-5.30549-0.33601

5.324942-0.96873

４. 制御則設計
無人航空機の自律飛行において，ピッチ角の制御を行う縦系とロール角，ヨー角の制御を行う
横・方向系の制御を行う．しかし，現在の飛行試験では，INS/GPS の方位角出力が大きく変動し，
82

その値を用いた姿勢角の算出が安定しないことから，ヨー角制御を実装していない．
a) 縦系制御則設計

有次元微係数を使用し，

𝑋𝑢̇𝛼̇[ 𝑞̇ ] = 𝑍̅𝑀′𝜃̇[

入力はエルロンのみの（１）式に示す運動方程式を縦系の機体ダイナミクスとして用い

𝑋
𝑍̅
𝑀′

−𝑔
57.3

𝑀′

𝑢
𝛼
𝑍̅
[𝑞 ] + ′
𝑀
𝑀′ 𝜃
𝜃
[
]

1

𝛿𝑒
[ ] ( )
𝛿𝑇]る（𝑢：x 軸方向速度，𝛼：迎角，𝑞：ピッチ角速度，𝜃：ピッチ角）．
このダイナミクスを用いて図２に示すピッチ角
制御則を用いて設計を行った．水平定常飛行時の制御要素は PID 制御とし，ピッチ角𝜃をゼロで一定となるように制御を行う．
縦系の機体
ダイナミクス

図２

ピッチ角制御b) 横・方向系制御則設計

図３

ピッチ角制御系ボード線図

横・方向系の制御則は非常に複雑であるため，
ロール角制御とヨー角制御に分解し，慣性能率のみを用いた回転体の運動方程式として簡略化し，設計を行った．
（φ：ロール角，p：ロール角速度，ψ：ヨー角，r：ヨー角速度）
姿勢角の変化がゼロとなるよう各軸における角度，
角速度をフィードバックする PD 制御を行っ
ている．

図４

図６

ロール角制御図５ ヨー角制御

ロール角制御系ボード線図

図７83ヨー角制御系ボード線図

５. 飛行結果
前述の機体及び誘導
制御システム構成にて，
白老滑空場にて水平定
常飛行試験を行った．
この試験での飛行軌跡を図８に示す．図中白
線は手動操縦，赤線は誘導制御システムによる自律飛行である．ま100m

た，この時の飛行状態のうち，ピッチ姿勢角，
図８

GPS 高度，GPS 速度履歴

飛行軌跡を図９～１１に示す．
飛行軌跡よりは，機体が徐々に右方向に進んでいくこと，ピッチ姿勢角履歴より，
ほとんど水平に姿勢が保たれており，その間の高度変動も 2m 程度，速度変動も 2m/sであることが確認できた．

図９ 姿勢角図１０

GPS 高度

図１１

GPS 速度

６. まとめラジコン機を対象に誘導制御システムを構築し，水平定常飛行試験を実施し，概ね良好な結果を得た．今後は，直進性の高い飛行実現に向けて，問題の明確化と対策を行うとともに，今回組
み込まれていないヨー角制御の実装による問題のない水平定常飛行をさせるとともに，今回使用されていない ADS を用いた旋回飛行制御則を始め，各種飛行モードに対する飛行試験を行う．

84

小型無人超音速機の縦系着陸制御系設計と外乱に対する性能評価
上羽 正純（もの創造系領域 教授）○ 横田 滋弘（機械航空創造系学科 4 年）

1. 研究の背景と目的
小型無人超音速実験機は，燃料重量の占める割合が高くそのため，時間とともに質量と重心位置が大きく変化する．これにより，重心位置を考慮した制御則が必要となる．また，着陸時の制御は飛行モード中最も複雑で，高い姿勢制御精度及び小さな降下率が要求される．降下率小さくすることは，破損や脚の強度，重量において重要な要素である．
本研究では，全ての飛行状態からの着陸を可能にするため，重心位置等による変化を反映した運動方程式の導出し，制御系設計用ダイナミクスを明らかにする．そのダイナミクスを使用して小型無人超音速実験機を対象に着陸制御系設計を行い，シミュレーションにより突風外乱を考慮した降下率を評価する．この結果を接地時降下率，接地高度の決定に反映する．

2. 機体の運動方程式
2.1 導出
図 1 に示すように，機体固定座標系𝑥𝑦𝑧の原点を定め，回転ベクトルをω，機体座標原点から重心までの位置ベクトル，重心から機体微小重量Δmまでの位置ベクトルをそれぞれ𝑟𝑐𝑔 ，𝑟𝑔 ，慣性座標系𝑋𝐸 𝑌𝐸 𝑍𝐸 からみた機体
座標原点と機体微小重量Δmの位置ベクトルをそれぞれ𝑟𝑐 ，𝑟 とし，重心と機体座標原点が異なる場合の運動方程式(1)，(2)を求めた．

図１

慣性能率は重心回りで求めている．







重心と機体固定座標現原点が異なる機体







(1)

 
L  my w  pv  qu   z v  ru  pw I p  I  I qr  I q  r   I q  rp  I r  pq M  mz u  qw  rv  x w  pv  qu   I q  I  I rp  I r  p   I  p  qr   I r  pq (2)
 N  mx v  ru  pw  y u  qw  rv I r  I  I  pq  I  p  q   I  p  qr   I q  rpm u  qw  rv  q 2  r 2 xcg   pq  r  y cg   pr  q z cg  Fx22
m v  ru  pw   pq  r xcg  p  r y cg  qr  p z cg  Fym w  pv  qu   pr  q xcg  qr  p  y cg  q 2  p 2 z cg  Fz

2

cg

cg

cg

cg

cg

cg

x

z

y

xy

2

y

x

z

y

x

xy

xz

2

xz

2

z

2

yz

xy

yz

xz

yz

2

2.2 妥当性の確認
外力・外モーメントが作用しない状態であれば，角運動量保存の法則に基づき角運動量は一定となる．これを用いて変数𝑢, 𝑣, 𝑤, 𝑝, 𝑞, 𝑟およびオイラー角𝜓, 𝜃, 𝜙をルンゲクッタ法により時間変化85

を求めたのち，機体座標原点での変数を重心位置に変換し，慣性座標系からみた機体重心位置での角運動量を求める，時間シミュレーションの結果を図 2 に示す．これにより，角運度量が一定となり，導出した運動方程式の妥当性が確認できた．
図2

3. 縦系制御系設計

慣性座標系での機体角運動量3.1 対象航空機モデル
小型無人超音速実験機としては M2006，M2011noseA を対象とし，風洞試験により或いは推算により得られた空力微係数を使用する．昇降舵のアクチュエータダイナミクスは暫定的に固有振動数𝜔𝑎 =20 Hz，ζ =0.1 の 2 次遅れ系，エンジンの取り付け角を 0°，サーボ，エンジン特性はそれ
ぞれ時定数 0.1 と 10 の 1 次遅れ系とした．また昇降舵の最大舵角は±20°，最大迎え角を±10°とした．
3.2 飛行経路における着陸時制御系
自動着陸における飛行経路を図 3 に示す．高度 1050m を水平定常飛行していた機体が，グライドスロープとのずれ 500m の状態から着陸を開始するとした．この時のトランスミッタからの距離を 1000m とした．
グライドスロープ角 3°にて水平速度一定にて降下する．
ピッチ角，グライドスロープ，オートスロッ
図3トル，フレア制御系は，それぞれ図 4，図 5，図飛行経路6，図 7 に示すもの制御系を採用した．
アクチュエータ

縦系の機体
ダイナミクス

図 4 ピッチ角制御系
サーボ

補助器
+

+

エンジン特性
縦系の機体
ダイナミクス)

図 5 オートスロットル制御系

86

オートパイロット（ピッチ角制御）

カプラ

図6

グライドスロープ制御系
オートパイロット（ピッチ角制御）

カプラ

図7

フレア制御系

4. シミュレーション

表1

4.1 シミュレーション条件

τの値

対象機体毎に設定されている着陸開始時

速度V(m/s)の水平定常速度を用いて，開始高度とフレ

トランスミッタからフレア開開始高度(m)
始地点までの距離a (m)
300
15.7
150
7.85
100
5.23

ア制御の比例定数τを表１のように設定した．これは，開始高度によらず，降下率をM2011 M2006
87
59.83
τ (s)
3.45
1.72
1.15

5.01
2.51
1.67

一定の条件より導かれている．この設計条
件において外乱の影響を調べるためにグライドスロープ制御時に突風外乱を上下±5m/s，フレア制御時には突風外乱を下向き 5m/s，それぞれ 1 秒間発生するとした．また，突風の外乱発生のタイミングについては，着地時の沈下率を重視して，高度 2m および 1m 時とした．
4.2 シミレーション結果
グライドスロープ，オートスロットル同時制御のシミュレーションでの外乱が作用しない場合
を図 8，作用した場合を図 9 に示す．
図 7 から，M2006，M2011 ともにグライドスロープとのずれ角が 0.07°である．また図 8 から，最大でグライドスロープとのずれが M2006 で 1°，M2011noseA で 0.9°になる．
M2011noseA

M2006

図8

外乱作用しない場合のシミュレーション結果（グライドスロープ）87

M2006

図9

M2011noseA

外乱作用した場合のシミュレーション結果（グライドスロープ）

表 2 に外乱作用によるフレア時の降下率とフレア区間の影響を示す．このとき，外乱なしでは理想に近いフレアを行っているので，漸近的に降下率が 0 に近づく．そのため，ある高度で接地すると考える必要がある．これを接地高度とする．外乱なしの降下率，つまり，式ℎ̇ = ( /τ) hよ
り求められる目標の降下率 0.1m/s のときの高度で接地するとした．図 10 示すように外乱作用高度が 1m の場合，接地時の降下率が大きくなっている．また，外乱作用高度は，理想時の高度であるため，接地高度決めると外乱作用高度がその分だけ低くなる．
表2

外乱作用高度2m

降下率とフレア区間のまとめ

M2006
開始高度 外乱作用 接地降下率 接地高度 フレア区間
(m)
高度(m)
(m/s)
(m)
(m)
なし
0.1
0.50
1000
15.7
2m
0.1
0.50
1040
1m
0.69
0.50
840
なし
0.1
0.25
480
7.85
2m
0.1
0.25
530
1m
0.66
0.25
350
なし
0.1
0.17
310
5.23
2m
0.1
0.17
470
1m
0.69
0.17
220
M20011noseA
開始高度 外乱作用 接地降下率 接地高度 フレア区間
(m)
高度(m)
(m/s)
(m)
(m)
なし
0.1
0.34
1100
15.7
2m
0.1
0.34
1200
1m
1.07
0.34
850
なし
0.1
0.17
530
7.85
2m
0.1
0.17
670
1m
1.12
0.17
360
なし
0.1
0.11
380
5.23
2m
0.1
0.11
520
1m
1.13
0.11
220

目標のモデル
降下率の応答外乱作用高度1m

図 10 M2006 外乱作用時のフレアシミュレーション結果(開始高度 15.7m)

5. まとめ小型無人超音速実験機 M2006 及び M2011 を対象に着陸時縦系制御系設計を行い，突風外乱を組み入れたシミュレーションにより接地時の降下率を評価した．この結果，接地高度の誤差が概ね0.1～0.3m 程度を想定した場合，突風外乱が発生する高度が低いと降下率への影響が大きく，それ以上の高度での突風外乱に対して，1 ケタ高い 1m/s 程度となることが明らかとなった．
今後は，着陸時の高迎角にも耐えうる横・方向系の制御系設計を行い，総合的なシミュレーシ
ョンによる評価を行う．
88

参考文献
1) 片柳亮二：航空機の飛行力学と制御，森北出版，2011．
2) 片柳亮二：航空機の飛行制御の実際，森北出版，2011．

89

無人航空機向けシステム同定の研究（-連続系伝達関数の開ループ同定-）
上羽 正純（もの創造系領域 教授）○ 山下 智也（航空宇宙システム工学専攻 博士前期 1 年）

1．研究の背景と目的
大陸間輸送及び地球周回軌道への位置輸送システムの実現には，高速・高高度を飛行するための空力，構造，制御系の基盤技術の確立が不可欠である．超音速で飛行する輸送システムにおいては，燃料の大量消費によって質量，慣性能率，重心変化ならびそれらを介して固有振動，固有振動モード等の構造特性を含んだダイナミクスが変化する．ダイナミクスの変化は，要求姿勢制
御精度に影響を与えるため，飛行中であっても正確に知る必要がある．ダイナミクスの変化は地上においてある程度事前に解析可能であるが，要求姿勢制御精度によっては十分な精度で把握できない場合も想定される．また，すべての飛行モードに対して風洞試験により空力微係数を測定することは不可能であり，ダイナミクスにはある程度の不確定性が存在する．さらに無人航空機の場合は，航空機重量に対する燃料重量の占める割合が高く大きなダイナミクス変化が発生することが想定される．したがって本研究では飛行しながらダイナミクスを精度よくリアルタイムで同定することを目指す．
ダイナミクスの同定は，システム同定法として確立されている最小 2 乗法の手法を用いて同定する．空力微係数の推定には多々使用されている(1)が，飛行している航空機に対して上記の手法を適用し，制御系設計に必要な伝達関数，及びそれを推定できる加振条件，センサの種類及びノイズの条件を検討した報告はない．したがって本研究で，無人航空機への適用を念頭に，まず亜音速の無人航空機としてラジコン機を対象に，伝達関数を精度良く同定する条件を明確化する．
2．対象とする運動方程式と航空機
2.1 縦系微小擾乱運動方程式
航空機ダイナミクスの内，2 次元平面内の運動に限定して縦系のダイナミクスを対象とする．この場合の運動方程式(2)は(1)式を使用する．
𝑋𝑢̇𝛼̇[ ] = 𝑍̅
𝑞̇
𝑀′
𝜃̇
[

𝑋
𝑍̅
𝑀′

−𝑔

𝑢
𝛼
𝑍̅
[ ] + [ ′ ] 𝛿𝑒
𝑞
𝑀
𝑀′ 𝜃
𝜃
]
57.3

𝑀′

(1)

2.2 対象の航空機
将来の飛行実験による実証を想定して図 1 に示す

図 1 対象とするラジコン機

ようなラジコン機を使用する．本ラジコン機の機体形
状より，空力微係数を推算し，シミュレーションに使表 1 ラジコン機の物理特性用した．表Ⅰにラジコン機の物理特性を示す．

機体重量(kg)

m

1.90

)

S

0.364

翼面積(
慣性能率(𝑘

90

)

0.132

3．システム同定
3.1 連続時間と離散時間
現実に存在するシステムの多くは時間的に連続の入出力値を有する連続システムである．しかしながら，データ処理においては，離散的に入出力信号を取得するため，必然的に離散システ
ムとして扱う．そこで今回扱う連続時間システムを次のような差分方程式として扱う．
𝜃𝑡 = −[𝑎1 𝜃𝑡−1 + ⋯ + 𝑎𝑛 𝜃𝑡−𝑛 ] + 𝑏0 𝛿𝑒𝑡 + 𝑏1 𝛿𝑒𝑡−1 ⋯ + 𝑏𝑛 𝛿𝑒𝑡−𝑛(2)

3.2 逐次最小 2 乗法式（2）を変形すると式（3），（4）に示すようなアルゴリズムが得られる．これが入出力から伝達関数の係数を求めることができる逐次最小 2 乗法と呼ばれるアルゴリズム(3)である．これは入力と出力が得られるたびに推定ゲインに代入することで，推定値をリアルタイムで更新することができる．ここで𝝎𝑁 は推定値，𝑷𝑁 は共分散行列，𝒛𝑁 は推定ゲイン，𝜃𝑁 は最新の出力値である．
𝑷

𝒛

𝑁

𝑁−1 𝑁

𝑁−1 𝑁
(𝜃𝑁 − 𝒛𝑁 𝑇 𝝎
̂𝑁 = 𝝎
̂ 𝑁−1 +
̂ 𝑁−1 )
𝝎
1+𝒛 𝑇 𝑷𝒛𝑷𝒛 𝒛 𝑇 𝑷𝑁−1𝑁
𝑁−1 𝒛𝑁𝑁−1 𝑁 𝑁𝑷𝑁 = 𝑷𝑁−1 −1 − 1+𝒛
𝑇𝑷

(3)
(4)

また，この式から得られる推定値は，離散時間領域の伝達関数の係数であるため，双一次変換
により連続時間系の伝達関数にし，推定精度をオリジナルの伝達関数と比較する．

4．シミュレーション
観測ノイズ分散値と計測時間をパラメータとして求められる伝達関数が，オリジナルの伝達関数と比較してどの程度復元できるかシミュレーションを行い，確認した．
4.1 シミュレーション条件
以下に示す条件にてシミュレーションを行う．
・観測ノイズ分散値 0.001°～1.0°
・計測時間 10 秒～50 秒
・入力 M 系列入力（-1°～1°）
4.2 シミュレーション結果
上記条件を用いてシミュレーションを行った結果のうち，計測時間 20 秒，50 秒の結果をそれぞれ図 2，図 3 に示す．
図 2，図 3 より，1Hz 以上の高周波数領域のゲイン曲線は，オリジナルのゲイン曲線とよく一致している．またノイズ分散値 1.0°では低周波数領域の共振点が推定できていない．
計測時間 20 秒（図 2）において，観測ノイズ分散値 0.001°，0.01°のケースでオリジナルのゲイン曲線とよく一致していることがわかる．計測時間を 50 秒（図 3）とした場合，観測ノイズ
分散値 0.001°では推定できたが，20 秒時に推定できていた観測ノイズ分散値 0.01°について精度が下がった．
以上より，計測時間 20 秒，ノイズ分散値が 0.01°のときオリジナルのゲイン曲線と最もよく一致し，また，計測時間を増加させても推定精度の向上につながっていない結果となった．

91

80

分散値0.01° 計測時間20秒
60

分散値0.1°
分散値1.0°ゲイン[dB]40
200

分散値0.001°

-20

オリジナル

-40
0.01

0.1

1
周波数[Hz]10100

図 2 推定ゲイン曲線 (計測時間 20 秒)80

分散値0.01°60

分散値0.1°

分散値1.0°40

ゲイン[dB]

計測時間50秒

20

オリジナル

0

分散値0.001°

-20
-40

0.01

0.1

1
周波数[Hz]10100

図 3 推定ゲイン曲線 (計測時間 50 秒)
5．まとめ
無人航空機としてラジコン機を対象に，その縦系のダイナミクスを逐次最小 2 乗法で推定し，計測時間 20 秒，観測ノイズ分散値 0.01°の条件で，ある程度伝達関数が復元できることを確認した．また計測時間の増加に対する推定精度の向上について理論的検討を進め，その後，閉ループ系に対して推定を行う．最終的には機体に組み込み，実証を行う．
参考文献
1) 加藤，「航空機力学入門」，東京大学出版会，1982
2) 片柳，「航空機の飛行力学と制御」，森北出版株式会社，2007
3) 中溝，「信号解析とシステム同定」，コロナ社，1988

92

高速走行軌道実験設備に関する基盤技術実証○ 中田 大将（航空宇宙機システム研究センター 特任助教）
小倉 達也（もの創造系 学部４年）
笹尾 鎮矢（もの創造系 学部４年）
ムハマド ハフィズ（もの創造系 学部４年）
東野 和幸（航空宇宙機システム研究センター 教授）
棚次 亘弘（航空宇宙機システム研究センター長 特任教授）

6. 概要
室蘭工業大学航空宇宙機システム研究センターでは地上で繰り返し安全に高速度環境を作り出すことの出来る高速走行軌道実験設備の基盤研究を進めている[1-9]．2013 年度は 2012 年度に引き続きフルサイズ軌道設備にて民間企業との共同研究により高加速度・高速度環境試験を実施した[10]．また，名古屋大学と共にサブスケール高速軌道装置を用いたローテーティング・デトネーショ
ン・エンジンの実証試験を行った[11]．これらについては本年次報告書の共同研究の項で述べる．
学内研究としては加速度キャンセル型空力天秤の改良，水制動装置の振動特性の解明，ハイブ
リッドロケットに関する着火遅れ特性の基礎研究[5,6,9]等を行った．長期的な展望として，大樹町に3 km の軌道設備を設置する計画について対外的な発表を行った[2,3]．
２．加速度補償型空力測定天秤の改良[7]大きな加速度がかかるスレッド上での空気力学測定では加速度に耐える秤量のロードセルを用いなければならず，測定精度の悪化要因となる．このため，機上での加速度をカウンターウェイトによってキャン
セルし，空気力に見合うサイズのロードセルを使用できる加速度補償型空力測定天秤の研究開発を
2010 年度より進めている．2013 年度は 2012 年度までの単純シーソー型の設計を見直し，上皿天秤等で使用されるロバーバル型と呼ばれる機構を採用したタイプについて製作した（図１）．

図1

サブスケール高速軌道実験装置に搭載したロバーバル型空力天秤
93

上皿天秤のような平行四辺形のリンク機構を，約 15cm の高さをとって上下二層に配置し，模型取付によるピッチングモーメントに充分耐える構造設計としている．ロバーバル型となること
で，左右のアーム作用点の誤差による精度の悪化は避けられる．この天秤を上向きに設置し，1Gの重力によるキャリブレーションを行ったところ，左右のおもり重量差分をロードセルひょう量（F.S. 5 N）の 1％以下（0.05 N/G）の誤差で検出した．このことから，静的な環境下における天秤そのものの精度は極めて高いと考えられる．

図2

1G の重力によるキャリブレーション試験

図１に示すようなセッティングで垂直平板を供試体として速度 20 m/s 程度までの走行試験を行った．動圧を計算する際にはピトー管で実測した対気速度または加速度センサの積分値による対地速度を用いた．対地速度が 15-20 m/s の範囲で Cd は概ね収束し，結果としては表１の通りとなった．
表1

新型天秤で求められた平板の抗力係数

対地速度を用いて求められた抗力係数1.0

対気速度を用いて求められた抗力係数0.8

文献値における垂直平板の抗力係数1.17

結果として，対気速度を用いて求められた抗力係数は文献値から 2 割程度のズレがあった．このように，静的環境下では良い精度を示したにも関わらず走行試験では良好な結果が得られなかった理由として，風速の測定誤差が支配的な要因であると推察された．今後，風速の取得精度を向上することが重要である．とりわけ，スレッド試験では短時間に速度が変化するため，ピトー
管や差圧計の時定数について考慮する必要がある．対地速度を用いて求められた抗力係数については，自然風の影響を考慮するとエラーバーに収まる[1]．
94３．水制動装置の振動特性に関する研究[8]高速軌道実験設備では加速されたスレッドの減速において，制動板を水路に浸す水制動装置を利用している．制動板の形状としては，１）単位面積当たり出来るだけ大きな制動力を発揮するようなもの，２）
水位に比例して的確な制動力が予測可能なもの などが研究されてきたが，これらに加えて３）制動時の振動の少ない形状を新たに検討した．従来の平板バケット，ダウンフォースを発揮して上下振動を抑制
するダウンフォース・バケット，中心力を発揮して横振動を抑制するスプリット・バケットを図 3 に示す．

図 3

振動を抑制する水制動板の試作（左から 平板バケット，ダウンフォース・バケット，ス

プリット・バケット）
これらの走行試験を実施した結果，ダウンフォース・バケットは明らかに上下方向の振動を抑制したが，（平板バケットにおいて実効値 1.24G であった振動を 0.53G に抑制）スプリット・バケットによって横方向の振動は抑制されなかった．この理由としては，サブスケール高速軌道装置のバケットはスレッド先頭に取り付けられており，重心よりも前にあることから擾乱に対し静安定を得られないことが考えられる．重心よりも後方に取り付けすることで今後，横方向の振動を抑制できるのではないかと考える．

４．ハイブリッドロケット着火遅れ特性の解明[5,6,9]
推進装置としてクラスタ化されたハイブリッドロケットを用いているが，図 4 に示すようなスタート時の着火のばらつきは予測された走行プロファイルからのズレを招く要因となる．これを解消するため，着火のばらつきを低減するための着火特性解明のための基礎実験を進めている．図 5 に示すような実験装置を用い，
図 6 に示すように燃焼室圧がほぼ定常となるまでの時間を着火遅れ時間と定義して諸条件に対する依存性を調べた．実験結果からは，酸化剤流量を増やすにつれて着火遅れおよびそのばらつきは短くなる傾向があることや，グレーン全長の依存性は小さく，ある時間以降は加速度的に燃え広がることなどが示唆された．また実用的な解決策として，内面に引火点の低いグリスを塗布することで平均着火時間を約 0.1
秒にまで短縮できることを実証した．

図 4 加速用ハイブリッドロケット４本クラスタ時の着火ばらつきの様子（左端から右端まで 0.5 秒）
95

図 5 ハイブリッドロケット着火遅れ特性確認実験装置レイアウト

図 6 ハイブリッドロケット着火遅れフェーズの分類（t1, t2, t3）
表 2 酸化剤流量と平均着火遅れ時間の関係

96５．3 km 軌道実験設備への展望
室蘭工業大学航空宇宙機システム研究センターでは，白老実験場における 300 m 軌道で得られた基盤技術を元に，北海道大樹町に 3 km クラスの軌道建設を計画している．建設予定地は大樹町航空公園
滑走路に隣接する区域で，海岸に向かって開けており，飛翔体の射出補助を行うことが可能である．
216 トンの 2 段式スペースプレーン親機をスレッド射出した場合の重量削減効果について検討した[7]．この計算によれば，スレッドで与えられる加速による燃料削減効果が 4 トン，ランディング・ギヤを省くことによる効果が 7 トンであり，後者が大きい（この場合，着陸はスキッドにて行う）．また，離陸速度を充分に大きく取れるため，主翼面積も僅かに小さくなる．これらの重量減により巡航中の重量や空気抵抗も小さくなることから，さらなる軽量化のフィードバックも期待できる．
このようなスペースプレーンを加速・射出するには推力 1400 トンの推進装置により 2.7G の加速度で 10秒間の加速を行い，終端速度 276 m/s で射出する．射出は x=1400-2000m の任意の点で 2 秒間の間に行われる．スペースプレーンを射出後のスレッド台車は最大 10G の多段階水路で減速され，2669 m で安全に停止する．なお，1400 トンの推力のうちほとんどは地上での動圧に打ち勝つために使用される．
表3

スペースプレーンの射出補助としてロケットスレッドを用いた場合の重量減

Parameter
FuselageWing
Landing-gear
Tank
ABE
Fuel (LH2)
Gross

Original
18.28 ton
15.56 ton
19.24 ton
9.56 ton
62.32 ton
61.4 ton
216.38 ton

表4

Sled assist
18.28 ton
14.56 ton
12.34 ton
9.06 ton
62.32 ton
57.6 ton
204.18 ton

Reduction

スペースプレーンを射出するスレッドの仕様

Gross initial weight
Payload length
Payload weight
Sled weight
Motor thrust
Motor initial weight
Burn time
Acceleration G
Top speed
Deceleration G (max)Running Distance

456 ton
76 m
374 ton
26 ton
1400 ton
56 ton10.0 sec
2.7 G
276 m/s
10 G
2669 m97

0.0 ton-1.0 ton-6.9 ton-0.5 ton
0 ton-3.8 ton-12.2 ton

図 7 スペースプレーンを射出するスレッドの距離＝速度線図．スペースプレーンは x =1400-2000 m の地点でリリースされる．

図8

スペースプレーンを射出するスレッドの時間＝速度線図．スペースプレーンは t = 10-12 sの間にリリースされる．

既に白老実験場において行われている他大学・民間企業との共同利用やこのような事例も含め，国内の多数の需要を取り込み，有用な装置として用いられることが期待できる[8]．
参考文献
[1]

D. Nakata, K. Nishine, K. Tateoke, K. Higashino, N. Tanatsugu, “Aerodynamic Measurement on theHigh Speed Test Track,” Aerospace Technology Japan, Accepted on Aug. 2013.
[2]

棚次亘弘，東野和幸，樋口健，中田大将，姫野武洋，渡邊力夫、「高速走行軌道試験設備」、

日本航空宇宙学会 第 44 期年会講演会、東京大学 山上会館、2013 年 4 月 18 日－19 日．
98

[3] D. Nakata, N. Tanatsugu, K. Higashino, K. Higuchi, T. Tsuchiya T. Himeno, “Rocket Sled Track
Facility as a Test Platform and a Launch Assist of Space Planes” Ground Based Space Facility Symposium,Jun. 2013, Paris.
[4] 中田大将，樋口健，棚次亘弘，東野和幸：ロケットスレッドの構造設計について，日本航空
宇宙学会・日本機械学会共催，第 55 回構造強度に関する講演会，1A09, JSASS-2013-3012, pp.25-27[5] 中田大将、ムハマド ハフィズ、棚次亘弘、東野和幸、廣田光智,ハイブリッドロケットの着火遅れに関する基礎実験，
宇宙輸送シンポジウム，2014 年 1 月 16－17 日，相模原市，STCP-2013-044.
[6] モハメドハフィズ、中田大将、棚次亘弘、東野和幸、廣田光智「ハイブリッドロケットの着火遅れ」、第１１回 HASTIC 学術講演会、北大、2014 年 2 月 24 日．
[7] 小倉達也，” サブスケール高速軌道走行搭載用ロバーバル型天秤に関する研究”，平成 25 年
度室蘭工業大学 卒業論文
[8]

笹尾鎮矢，” 高速軌道装置における振動を抑制する制動板の研究”，平成 25 年度室蘭工業大学 卒業論文[9]

ムハマド

ハフィズ，”ハイブリットロケットの着火遅れに関する基礎実験と考察”， 平成

25 年度室蘭工業大学

卒業論文[10] 室工大 B-043，

高加速度環境下における高速走行軌道実験の発展[11]2014 年 3 月

側原圭太、加藤優一、松岡健、笠原次郎、松尾亜紀子、船木一幸、中田大将、東野和幸、

棚次亘弘、高速走行軌道装置における回転デトネーションエンジン滑走試験に関する研究,宇宙輸
送シンポジウム，2014 年 1 月 16－17 日，相模原市，STCP-2013-033.

99

発表論文
2013 年度
査読付き学術誌［1］ 髙木正平、山谷直広、伊藤信毅： 人工音響ループによる二次元噴流不安定波周波数の制御, 日本航空宇宙学会論文集、第 61 巻、3 号、pp.57-63, (2013)［2］ 石川温士、今井良二、田中貴博、SF6-水を用いた水平管群まわりの気液二相流計測 : 高
温高圧の PWR/SG 二次冷却水側の気液二相流を常圧で模擬、配管技術、56(1), P.28-34,
2014.
［3］ Jun KOYANAGI, Akihito WATANABE, Nobuyoshi KAWABATA, Tsuyoshi OZAKI, Ken
HIGUCHI, Kosei ISHIMURA and Yasumasa KASABA, "Long-term durability of tri-axial
woven CFRP tube structure extended along the spin axis of spinning platforms for the SCOPE
mission," Advanced Composite Materials, Vol.23, No.2, pp.115–128, 2014.
［4］ Higashino，K., Sugioka, M., Izumi, S., Terada, T., Azuma，N., Kobayashi, T., “Study on
Material Compatibility for Bioethanol Fueled Rocket Engine,”JSASS Aerospace Technology
Japan, Accepted, on Aug. 23, 2013.
［5］ Mizobata, K., Minato, R., Higuchi,K., Ueba, M., Takagi, S., Nakata, D., Higashino, K., andTanatsugu, N., “Development of a Small-scale Supersonic Flight Experiment Vehicle as a
Flying Test Bed for Future Space Transportation Research,” ISTS Special Issue: Selected
papers from the 29th International Symposium on Space Technology and Science,
Transactions of JSASS, Aerospace Technology Japan, Accepted on Sep. 5, 2013.
［6］ Nakata, D., Nishine, K., Tateoke, K., Higashino, K., Tanatsugu, N., “Aerodynamic
Measurement on the High Speed Test Track,” Aerospace Technology Japan, Accepted, on Aug.
23, 2013.
［7］ Ken HIGUCHI, Yasuyuki MIYAZAKI, Kosei ISHIMURA, Hiroshi FURUYA, Hiroaki
TSUNODA, Kei SENDA, Akihito WATANABE, Nobuyoshi KAWABATA, Takeshi
KURATOMI, and SIMPLE Project Team, "Initial Operation and Deployment Experiment of
Inflatable Extension Mast in SIMPLE on JEM Exposure Platform in ISS," Trans. JSASS
Aerospace Tech. Japan, Accepted on Jan. 24, 2014.
［8］ 塚野徹, 棚次亘弘, 杉岡正敏, 東野和幸, 白金触媒を用いた熱分解吸熱性燃料の吸熱量の測定と評価, 日本航空宇宙学会論文集，2014 年 2 月 掲載可100

国際会議発表論文
［1］

［2］

［3］

［4］

［5］

Atsushi Ishikawa,Ryoji Imai, Takahiro Tanaka, Measurement of Two-Phase Cross-Flow of SF6
Gas-Water in Horizontal Tube Bundle, The 15th International Topical Meeting on Nuclear Reactor
Thermal-Hydraulics, Pisa, Italy, May 2013.
Kazuyuki Higashino，Masatoshi Sugioka,Shuntaro Izumi,Toshiyuki Terada,Noboyuki Azuma，Teiu Kobayashi, “Study on Material Compatibility for Bioethanol Fueled Rocket Engine” a30 The
29th International Symposium on Space Technology and Science,Nagoya， Aichi, Japan, 2013
Kazuhide Mizobata, Ryojiro Minato，Ken Higuchi，Masazumi Ueba，Syohei Takagi，Daisuke
Nakata， Kazuyuki Higashino，Nobohiro Tanatsugu, “Development of a Small-scale Supersonic
Flight Experiment Vehicle as a Flying Test Bed for Future Space Transportation Research,”
2013-o-3-02, The 29th International Symposium on Space Technology and Science, Nagoya， Aichi,
Japan, June 2013
Masazumi Ueba and Kazuhide Mizobata, “Study on landing control for a small-scale unmanned
supersonic airplane,” 2013-g-01, The 29th International Symposium on Space Technology and
Science, Nagoya，Aichi, Japan, June 2013
Nobukatsu Okuizumi, Yuki Ito, M.C.Natori, Nobuhisa Katsumata and Hiroshi Yamakawa,“Numerical Simulation of Stepwise Deployment of Membrane Structure with Booms using
Multi-Particle Approximation Method” The 29th International Symposium on Space Technology
and Science, Nagoya， Aichi, Japan, 2013

［6］

Takahira AOKI, Ken HIGUCHI, Kazuki WATANABE, and SIMPLE Project Team, "Progress
Report of SIMPLE Space Experiment Project on ISS Japan Experiment Module," 2013-c-07, The
29th International Symposium on Space Technology and Science, Nagoya， Aichi, Japan, June
2013.

［7］

Ken HIGUCHI, Yasuyuki MIYAZAKI, Kosei ISHIMURA, Hiroshi FURUYA, Hiroaki TSUNODA,
Kei SENDA, Akihito WATANABE, Nobuyoshi KAWABATA, Takeshi KURATOMI, and
SIMPLE Project Team, "Initial Operation and Deployment Experiment of Inflatable Extension Mast
in SIMPLE on JEM Exposure Platform in ISS," 2013-c-09, The 29th International Symposium onSpace Technology and Science, Nagoya， Aichi, Japan, June 2013.

［8］

Daisuke NAKATA, Nobuhiro TANATSUGU, Kazuyuki HIGASHINO, Ken HIGUCHI, Takeshi
TSUCHIYA and Takehiro HIMENO“ Rocket Sled Track Facility as a Test Platform and Launch
Assist of Space Planes”，Ground Based Space Facility conference，Paris, 2013

［9］

Ryojiro Minato, Kazuhide Mizobata、Kazuyuki Higashino, Nobuhiro Tanatsugu、“Characteristics
of Bio-Ethanol Fueled GG-cycle Air Turbo Ramjet Engine for Supersonic UAV along Flight
Trajectory” , 49th Joint Propulsion Conference, San Jose, 2013[10]

Ryojiro Minato, Kazuyuki Higashino, Nobuhiro Tanatsugu, ”Design and Development of
Bio-Ethanol Fueled GG-cycle Air Turbo Ramjet Engine for Supersonic UAV” 21 st International
Society for Air Breathing Engines, Busan, Korea, 2013.

[11]

Takagi, S., and Ueda, Y.：Active control of flow separation on an airfoil wing with the use of a rowof vortex generators, The 14th Asian Congress of Fluid Mechanics (14ACFM), Hanoi and Halong,
Vietnam, (10/15-19, 2013)[12]

Nobuhisa Katsumata, Masaaki Kaimori, Ken Higuchi, M.C. Natori, Hiroshi Yamakawa,"Deployment Characteristics of Braided Coated Bi-Convex tape Booms partly used SMA for101

Deployable Membrane Structures," 24th International Conference on Adaptive Structures andTechnologies (ICAST2013), October 7-9, 2013, Aruba.
[13]

Kazuhide MIZOBATA, Ryojiro MINATO, Kazuyuki HIGASHINO, and Nobuhiro TANATSUGU,“Aerodynamic Characterization of a Small-scale Supersonic Flight Experiment Vehicle with a
Cranked-arrow Main Wing and Tails by Wind-tunnel and Flight Tests,” The 2013 Asia-PacificInternational Symposium on Aerospace Technology, Nov. 20-22, 2013, Takamatsu, JAPAN.

[14] Ohtateme, Y. and Takagi, S.: On frequency-selection mechanism of vortex street behind 2-D airfoil,
Joint Symposium on Mechanical - Industrial Engineering, and Robotics 2013(MIER2013), (Chiang
Mai, 11/15-16).

国内学会発表論文［1］ 棚次亘弘，東野和幸，樋口 健，中田大将(室工大)，姫野武洋(東大)，渡邊力夫(都市大)、「P05
高速走行軌道試験設備」、日本航空宇宙学会

第44期年会講演会、東京大学 山上会館、

2013年4月18日－19日．
［2］ 髙木 正平、横山慶輔：後流不安定と超音速縮流胴境界層の再遷移、第 53 回乱流遷移研究会、(東京, 9/24-25, 2013).
［3］ 山本康平、塚野徹（室蘭工業大学大学院）、湊 亮二郎、杉岡正敏、東野和幸（室蘭工業大学）「２Ａ１０ バイオエタノールの熱分解吸熱反応に関する研究～特に準静的環境における熱分解特性～」、第５７回宇宙科学技術連合講演会、鳥取県、2013年10月09－10月11日．
［4］ 宮森康臣、畠中和明、髙木正平：石鹸膜に形成される柱体後流渦列のシュリーレン法による可視化, 高速度イメージングとフォトニクスに関する総合シンポジウム 2013
(JCHSIP2013), (室蘭, 10/17-19, 2013).
［5］ 髙木正平、上田祐士：小翼列を用いた翼剥離流のActive制御、第51回飛行機シンポジウム、（高松、11/20-11/22）.
［6］ 中田大将，樋口 健，棚次亘弘，東野和幸：ロケットスレッドの構造設計について，日本航空宇宙学会・日本機械学会共催，第 55 回構造強度に関する講演会，1A09, JSASS-2013-3012,
pp.25-27, (2013.8.7-9), 室蘭.
［7］ 相原弘匡，岩井達也，似鳥 透，樋口 健，岩佐貴史，岸本直子，藤垣元治：格子投影法２
カメラ計測法による大型構造物の面形状計測と精度評価，日本航空宇宙学会・日本機械学会共催，第 55 回構造強度に関する講演会，2A06, JSASS-2013-3050, pp.136-138, (2013.8.7-9),
室蘭.
［8］ 岩佐貴史，小野静之，後垣雄太，岸本直子，相原弘匡，似鳥 透，樋口 健，小木曽望，藤垣元治，ISAS/JAXA 戦略開発費高精度計測チーム：点計測と面計測を統合した大型宇宙構造システムの表面形状計測法に関する検討実験，日本航空宇宙学会・日本機械学会共催，
第 55 回構造強度に関する講演会，2A07, JSASS-2013-3051, pp.139-141, (2013.8.7-9), 室蘭.
［9］ 近江信耶，蒔田愛道，樋口 健，宮崎康行，古谷 寛，青木隆平，川端信義，渡邊秋人，坂本信臣，SIMPLE-IEM 実験チーム：宇宙インフレータブル伸展マストの軌道上初期運用，
日本航空 宇宙学 会・日 本機械学 会共催 ，第 55 回構造 強度に 関する 講演会， 3A10,102

JSASS-2013-3080, pp.225-227, (2013.8.7-9), 室蘭.
［10］ 浅野等, 五明

泰作, 河南治, 廣川

良二, 川崎 春夫, 藤井 清澄, 栗本

智己, 新本 康久, 大田 治彦, 鈴木卓, 駒崎康一, 今井雅人,二相流体ループ ISS 実験装置のシステム特性,日本混相流学会混相流シンポジウム 2013,2013 年 8 月 9 日,長野市（信州大学）．
［11］ 和田大輔，今井良二，三原礼, 極低温流体による予冷時の沸騰熱伝達特性の計測, 日本混
相流学会混相流シンポジウム 2013,2013 年 8 月 9 日，長野市（信州大学）．
［12］ 岸本直子，小野静之，後垣雄太，岩佐貴史，樋口 健，藤垣元治：格子投影法による大型構造物の面形状計測と精度評価, 日本機械学会, 2013 年度年次大会, 岡山, (2013.9.8-11).
［13］ 岸本直子，岩佐貴史，樋口健：格子投影法を用いた３次元形状計測の高精度化実証実験について，第 57 回宇宙科学技術連合講演会，1K14，JSASS-2013-4191，米子（米子コンベンションセンターBigShip），(2013.10.9-11).
［14］ 樋口 健，青木隆平，宮崎康行，渡邊秋人，川端信義，堀 利行，伊藤裕明：SIMPLE 宇宙
実験に供されたインフレータブル伸展マストの開発試験，第 57 回宇宙科学技術連合講演会，3M12，JSASS-2013-4676，米子，(2013.10.9-11).
［15］ 栗本 卓, 松本 聡, 駒崎 雅人, 澤田 健一郎, 川崎 春夫, 大田 治彦, 新本 康久, 鈴木 康
一, 浅野 等, 河南 治, 今井 良二,「きぼう」における沸騰二相流実験装置の開発状況,日
本マイクログラビティ応用学会第 27 回学術講演会(JASMAC-27),2013 年 11 月 29 日．
［16］ 樋口 健，相原弘匡，似鳥 透，松尾剣都，岩佐貴史，岸本直子：格子投影法２カメラ法における繰り返し計測精度，第 29 回宇宙構造・材料シンポジウム，A11，(2013.12.3),
ISAS/JAXA 相模原．
［17］ 貝森政明，勝又暢久，樋口 健：形状記憶合金を用いた組紐被服ブームの展開挙動評価，
日本機械学会宇宙工学部門，
第 22 回スペース・エンジニアリング・コンファレンス(SEC'13)，(2013.12.20-21)，東京．
［18］ 藤原直也，勝又暢久，樋口 健：大型宇宙構造物の実現に向けた自動構築メカニズムの開発，日本機械学会宇宙工学部門，第 22 回スペース・エンジニアリング・コンファレンス(SEC'13)，(2013.12.20-21)，東京．
［19］ 岩佐貴史，牧 恒男，原田 卓，岸本直子，亀井宏貴，樋口 健，相原弘匡，藤垣元治，吉川隆章，小木曽望，田中宏明，石村康生，河野太郎，神谷友裕，古谷 寛，坂本 啓，佐藤泰貴，中篠恭一，荻 芳郎，泉田 啓，仙場敦彦：高精度大型宇宙構造システムの開発研究－高精度形状計測技術の開発研究－，ISAS/JAXA 第 14 回宇宙科学シンポジウム，相模原，(2014.1.9-10).
［20］ 溝端―秀、湊亮二郎、樋口 健、勝又暢久、上羽正純、中田大将、高木正平、東野和幸、
棚次亘弘(室蘭工大) ，室蘭工大の小型超音速飛行実験機(オオワシ)の研究開発の現状，宇宙輸送シンポジウム，2014 年 1 月 16－17 日，相模原市，STCP-2013-026.
［21］ 鈴木祥弘、大石栄(室蘭工大院)、上村隆太、塩田光平、行澤陽、溝端―秀、棚次亘弘(室蘭工大)、新井隆景(阪府大), 室蘭工大の小型超音速飛行実験機(オオワシ)の操舵空力特性,宇宙輸送シンポジウム，2014 年 1 月 16－17 日，相模原市，STCP-2013-027.
［22］ 側原圭太(筑波大院)、加藤優―(名大工院)、松岡健(名大)、笠原次郎(名大)、松尾亜紀子(慶大)、船木―幸(ISAS/)AXA)、中田大将(室工大)、東野和幸(室工大)、棚次亘弘(室工大)、高103

速走行軌道装置における回転デトネーションエンジン滑走試験に関する研究,宇宙輸送シンポジウム，2014 年 1 月 16－17 日，相模原市，STCP-2013-033.
［23］ 中田大将(室工大)、棚次亘弘、湊亮二郎、軸流反転ファンの基礎実験, 宇宙輸送シンポジ
ウム，2014 年 1 月 16－17 日，相模原市，STCP-2013-036.
［24］ 高木正平(室工大)、上田祐士、伏木田裕介，平田裕, 低速から超音速境界層の空力制御について, 宇宙輸送シンポジウム，2014 年 1 月 16－17 日，相模原市，STCP-2013-041.
［25］ 正木陽(室工大・院)、高木正平、棚次亘弘(室工大), ISAS 高速風洞の気流特性, 宇宙輸送シンポジウム，2014 年 1 月 16－17 日，相模原市，STCP-2013-042.
［26］ 中田大将[164]、ムハマド ハフィズ、棚次亘弘、東野和幸、廣田光智（室工大）,ハイブリ
ッドロケットの着火遅れに関する基礎実験，宇宙輸送シンポジウム，2014 年 1 月 16－17
日，相模原市，STCP-2013-044.
［27］ 杉岡正敏，東野和幸，神田康晴，板倉賢一（室蘭工業大学）UGC ガスの脱硫と水素製造における触媒開発、第 14 回北海道資源環境研究発表会、北大、2014 年 1 月 28 日．
［28］ 塚野徹、山本康平、杉岡正敏、棚次亘弘、東野和幸（室蘭工業大学）「 炭化水素燃料における熱分解吸熱反応」、第１１回 HASTIC 学術講演会、北大、2014 年 2 月 24 日．
［29］ モハメドハフィズ、中田大将、棚次亘弘、東野和幸、廣田光智（室蘭工業大学）「ハイブ
リッドロケットの着火遅れ」、第１１回 HASTIC 学術講演会、北大、2014 年 2 月 24 日．
［30］ 飯島 明日香、笹木 康平、山本 康平、塚野 徹(室蘭工業大学大学院),湊 亮二郎、杉岡 正
敏、棚次 亘弘、東野 和幸(室蘭工業大学

航空宇宙機システム研究センター)： バイオ

エタノールの吸熱特性に関する研究～特に白金アルミナ触媒を用いた場合、日本航空宇宙学会北部支部 2014 年講演会ならびに第 15 回再使用型推進系シンポジウム、東北大、2014
年 3 月 10－11 日．
［31］ 近藤賢，鈴木祥弘，ムハマッドハナン，溝端一秀，湊亮二郎，東野和幸，棚次亘弘，「室蘭工大の小型超音速飛行実験機の空力特性と飛行性能予測」，日本航空宇宙学会北部支部
2014 年講演会ならびに第 15 回再使用型推進系シンポジウム、東北大、2014 年 3 月 10－
11 日．
［32］ 渡口 翼，小野正貴，溝端一秀，湊亮二郎，東野和幸，棚次亘弘，福士誠，山中大樹，梅村諒磨、「室蘭工大の小型超音速飛行実験機の予備的飛行試験」，日本航空宇宙学会北部支部 2014 年講演会ならびに第 15 回再使用型推進系シンポジウム、東北大、2014 年 3 月 10
－11 日．
［33］ 上羽正純、溝端一秀（室蘭工業大学大学院）「1C10

小型無人超音速実験機の着陸制御の一検討」、第５１回飛行機シンポジウム、香川県、2013 年 11 月 20－11 月 22 日．
［34］ 山下智也、上羽正純（室蘭工業大学大学院）「１C１２

小型無人超音速実験機のリ

アルタイムダイナミクス同定法」、第５１回飛行機シンポジウム、香川県、2013 年 11 月 20
－11 月 22 日．

104

2008～2012 年度
査読付き学術誌［1］

棚次亘弘，「室蘭工業大学における航空宇宙分野の実践教育」，工学教育，第 56 巻，第 5 号，pp.79-83，2008．

［2］

湊亮二郎，溝端一秀，桑田耕明，「吸込み式超音速風洞におけるスターティングロードの計測と供試体の挙動」，日本航空宇宙学会論文集，Vol.56, No.659，pp.580-589, 2008．

［3］

吹場活佳，井上翔，佐藤哲也，大久保英敏，「予冷ターボジェットエンジンにおけるジェット噴射を用いた熱交換器の除霜」，日本航空宇宙学会論文集，Vol.56, No. 657, pp. 464-470, 2008．

［4］

渡邉真也，湊亮二郎，「多数非劣解集合からの設計支援手法の開発－ジェットエンジン最適化を通して」，人工知能学会論文誌，24 巻 1 号，2009．

［5］

Katsuyoshi Fukiba, Shou Inoue, Hidetoshi Ohkubo and Tetsuya Sato, “A New Defrosting Method
Using Jet Impingement for Precooled Turbojet Engines”, Journal of Thermophysics and Heat
Transfer, Volume 23, Number 3, pp. 533-542,2009．

［6］

東野和幸，杉岡正敏，小林隆夫，湊亮二郎，丸祐介，笹山容資，大塚雅也，牧野隆，坂口裕之，「ＬＮＧロケットエンジンのコーキング特性に関する基盤研究」，日本航空宇宙学会論文集，Vol. 57,
pp. 210-216, 2009.

［7］

Nishizawa, A., N. Tokugawa and Takagi, S.,”Critical Reynolds number off theattachment-line boundary layer on a yawed cylinder”, Fluid Dynamics Research, Vol.
41, No.3, June (2009) pp.385-392.

［8］

Atobe, T., Tuinstra, M., and Takagi, S., “A study on airfoil tonal noise generation inresonant environments”, Transactions of Japan Society of Aeronautical and Space
Sciences, Vo.52, No.176 (2009) pp.74-80.

［9］

Takagi, S. and Konishi, Y.,”Suppression of Trailing-Edge Noise Emitted by TwoDimensional Airfoils”, Transactions of Japan Society of Aeronautical and Space
Sciences, Vol.53, No.179(2010), pp.19-23.

［10］ Minato. R, Mizobata, K., and Kuwada, K., “Experimental Measurements of Starting
Loads and Model behaviors in the Indraft Supersonic Wind Tunnel”, Transaction of theJapan Society for Aeronautics and Space Science Vol.53, No.179(2010), pp.54-62．
［11］ 東野和幸，杉岡正敏、小林隆夫、境昌弘、湊亮二郎、「LNG ロケットエンジンのサルファアタック
及びコーキングに関する基礎的研究」、日本航空宇宙学会論文集、 vol57, No.670(2010),
pp.445-452.
［12］ 東野和幸、杉岡正敏、小林隆夫、湊亮二郎、大屋俊輔、笹山容資、「グラファイトコーティングを利用した再生冷却型 LNG ロケットエンジンのコーキング抑制」、日本航空宇宙学会論文集、vol58、No.676)2010)、 pp.138-145．
［13］ Takagi, S. and Konishi, Y., "On the Frequency Selection Mechanism of Airfoil
Trailing-Edge Noise", Journal of Aircraft, Vol. 47, No.4, (2010) pp.1111-1116.
［14］ 東野和幸、駒崎慎一（室蘭工大）、木元健一（ＩＨＩ）、坂口裕之（ＩＨＩエアロスペース）、石橋利幸（ＮＥＧＥＮＴ）、「ロケット燃焼室におけるニッケル電鋳特性向上に関する基礎研究」、 日本航空宇宙学会論文集、vol58, No.676(2010)、pp.146-151．
105

［15］ 吹場活佳、坪井伸幸、小林弘明、「ノーズコーンをもつ超音速機の壁面圧力を用いた機体姿勢およびマッハ数計測に関する数値解析」、日本航空宇宙学会論文集、Vol.58、No. 676(2010)、pp.
130-137.
［16］ .廃アルミニウムからの水素製造、神田康晴、上道芳夫、杉岡正敏、化学工学、74 巻、2 号、p.60-62(2010)［17］ Y. Kanda, C. Temma, K. Nakata, T. Kobayashi, M. Sugioka, Y. Uemichi,” Pｒeparationand performance of noble metal phosphides on silica as new hydrodesulfurization
catalysts” , Applied Catalysis A;General,Vol.386, pp.171-178，2010.
［18］ Tetsuya Sato, Hideyuki Taguchi, Hiroaki Kobayashi, Takayuki Kojima,Katsuyoshi
Fukiba, Daisaku Masaki Keiichi Okai, Kazuhisa Fujita,Motoyuki Hongo, Shujiro Sawai,“Development study of a precooled turbojet engine”, Acta Astronautica, Volume 66,
Issue 7-8s, April-May 2010, pp. 1169-1176.
［19］ Akira Iwakawa, Daisuke Nakata, and Hitoshi Kuninaka, “Experimental Study of a
Two-Dimensional Applied Field Magnetoplasmadynamic Thruster,” Aerospace
Technology Japan, Vol. 8 (2010) pp.Pb_13-Pb_17, 2010．
［20］ Y. Kasaba, A. Kumamoto, K. Ishisaka, H. Kojima, K. Higuchi, A. Watanabe, and K.
Watanabe, “Development of stiff and extensible electromagnetic sensors for space
missions,” Advances in Geosciences, Solar & Terrestrial Science, Vol.21, pp.447-459,(2010)［21］ 神田康晴、小林隆夫、上道芳夫、杉岡正敏、「廃アルミニウムによる新規水素製造法の開発―ア
ルミニウム廃材を用いた水素製造法のスケールアップについてー」、配管技術、Vol.53, No.1,
pp.7-11(2011).
［22］ Daisuke Nakata, Kyoichiro Toki, Ikkoh Funaki, Hitoshi Kuninaka、” Performance of
ThO2-W, Y2O3-W and La2O3-W cathodes in Quasi-Steady Magnetoplasmadynamic
Thrusters” , Journal of Propulsion and Power, Volume 27, Number 4, Jul. – Aug.
2011.pp. 912-915..
［23］ Asai, M., Inasawa, A., Konishi, Y., Hoshino, S. and Takagi, S.: "Experimental
investigation of the instability of wakes of axisymmetric streamline body", Journal of
Fluid Mechanics, 675 (2011), pp. 574-595.
［24］ 笹山容資、境昌宏、小林隆夫、杉岡正敏、東野和幸、東伸幸、青木賢司、小林完、沖田耕一、「LNGロケットエンジンにおけるサルファアタック防止に関する研究 -特に金メッキの特性評価-」，
日本航空宇宙学会論文集，59 No.688（201１）， pp.138-145.
［25］ Yasuharu KANDA, Shinji KONDO, Shunsuke OOYA, Takao KOBAYASHI, Yoshio UEMICHI,
Kazuyuki HIGASHINO and Masatoshi SUGIOKA, “Green Hydrogen Production by Mechanical
Mixing of Aluminum with Water”, Journal of Chemical Engineering of Japan, Vol.44, No.10,ｐｐ803-808，2011［26］ T. Ikeda, T. Atobe and S. Takagi: Direct simulations of trailing-edge noise generationfrom two-dimensional airfoils at low-Reynolds numbers, Journal of Sound andVibration, 331 (2011) pp. 556–574.
［27］ M.Ohira, T.Umaba, S.Kitazawa, H.Ban, M.Ueba, “Experimental Characterization of
Microwave Radio Propagation in ICT Equipment for Wireless Harness
Communications”, IEEE Transactions on Antennas and Propagation, Vol.50, No.12, Dec.
106

2011, pp.4757-4765[28]

東野和幸、杉岡正敏、棚次亘弘、湊亮二郎、笹山容資、磯田浩志、「熱分解吸熱反応燃料の吸熱量の測定と評価」、日本航空宇宙学会論文集，60 No.3（2012）， pp.115-120[29]

寺島 修、大西 一弘、酒井 康彦、長田 孝二、髙木 正平；定温度型熱線流速計の周波数応答特性を向上させる試み、日本機械学会論文集（B 編）, 78 巻, 788 号 (2012-4), pp. 862-866．

[30]

Daisuke Nakata, Ami Kozu, Jun Yajima, Kenji Nishine, Kazuyuki Higashino andNobuhiro Tanatsugu, “Predicted and Experimented Acceleration Profile of the Rocket
Sled,” Aerospace Technology Japan, Vol. 10, No. ists28, pp. Ta_1-Ta_5, 2012[31] Itoh, N., Takagi, S. and Ikeda T., “Instability and frequency selection of the wakebehind a flat plate”, Transactions of the Japan Society for Aeronautics and Space
Science, Vol.55, No.6(2012)356-363．
[32]

Y. Ogi, K. Higuchi and K. Ishimura: “Effect of Attachment Errors of Flexible
Appendages on the Spin Axis of a Rigid Body,” ISTS Special Issue, Transactions of theJapan Society for Aeronautical and Space Sciences, Aerospace Technology Japan, Vol.10,
No.ists28, 2012, pp. Pc_7-Pc_12.

[33] Kosei ISHIMURA, Tsuneo KII, Keiji KOMATSU, Ken GOTO, Ken HIGUCHI, Kazuro
MATSUMOTO, Shoichi IIKURA, Makoto YOSHIHARA, Masaharu TSUCHIYA,“Shape
Prediction of Large Deployable Antenna Structure on Orbit,” TRANSACTIONS OF
THE JAPAN SOCIETY FOR AERONAUTICAL AND SPACE SCIENCES,
AEROSPACE TECHNOLOGY JAPAN, Vol.10, No.ists28, pp.Pc_35-Pc_40, (2012).
[34] Hiroshi KAWAMURA, Takashi IWASA, Keita UNEMITSU, Ryo SUGIYAMA, Naoko
KISHIMOTO, Ken HIGUCHI, Motoharu FUJIGAKI and Takayuki SHIOKAWA,“Visualization on Formation Process of Wrinkling Phenomena on Thin Membranes
Using Grating Projection Method and Investigation on Their Nonlinear Mechanical
Properties,” TRANSACTIONS OF THE JAPAN SOCIETY FOR AERONAUTICAL AND
SPACE SCIENCES, AEROSPACE TECHNOLOGY JAPAN, Vol.10, No.ists28,
pp.Pc_91-Pc_97, (2012).
[35]

岩佐貴史，采光啓太，杉山 遼，岸本直子，樋口 健，藤垣元治，塩川貴之：格子投影法によるポリイミドフィルムの表面形状計測，航空宇宙技術，Vol.11, pp.69-78, (2012).

国際会議発表論文［1］

Minato, R., Arai, T., Himeno, T., Kobayashi, H., Mizobata, K. and Tanatsugu, N., “Off-design
Analysis of Turbojet Engine for Supersonic Unmanned Plane”, Asian Joint Conference on
Propulsion and Power, Kitakyushu, 2005.

［2］

Kaneko, H., Minato, R., Arai, T., Himeno, T., Kobayashi, H., Mizobata, K. and Tanatsugu, N.
“Performance Analysis of Turbojet Engine with After Burner for Supersonic Unmanned Vehicle”,
Asian Joint Conference on Propulsion and Power, Kitakyusu, 2005.

［3］

Minato, R., Niioka, T., Sugiyama, H. and Mizobata, K., “Numerical Analysis of Supersonic107Combustion by a Plasma Torch,” 13th AIAA/CIRA International Conference on Space Planes andHypersonic Systems and Technologies, AIAA Paper 2005-3424, Naples, Italy, 2005.
［4］

Minato, R., Arai, T., Himeno, T., Kobayashi, H., Mizobata, K. and Tanatsugu, N. “Preliminary
Analysis of Turbojet Engine at Off-design Condition for Subscale Supersonic Unmanned Plane,”
13th AIAA/CIRA International Conference on Space Planes and Hypersonic Systems andTechnologies, AIAA Paper 2005-3415, Naples, Italy, 2005.

［5］

Arai, T., Mizobata, K., Minato, R., Tanatsugu, N, Mori, Y. and Kudo, T., “Correlation between
Fluctuation of Mass Flux and Hydrogen Concentration in Supersonic Mixing,” 13th AIAA/CIRA
International Conference on Space Planes and Hypersonic Systems and Technologies, AIAA Paper
2005-3236, Naples, Italy, 2005.

［6］

Mizobata, K., Minato. R, Kimura, H., Arai, T., Sugiyama, H., and Tanatsugu, N., “Development
Study on a Small-scale Supersonic Flight Experiment Vehicle with Air-breathing Propulsion,” 13th
AIAA/CIRA International Conference on Space Planes and Hypersonic Systems and Technologies,
AIAA Paper 2005-3347, Naples, Italy, 2005.

［7］

Minato, R., Tokura, I., Tanatsugu, N., Himeno, T., Kojima, T., Kobayashi, H., “Preliminary Designof a Small-Scale Turbojet Engine for a Supersonic Flight Experimental Vehcle.” 25 th International
Symposium on Space Techmology and Science, ISTS 2006-g-22, Kanazawa, 2006.

［8］

Kazuhide Mizobata, Ryojiro Minato, Nobuhiro Tanatsugu, Hiroyuki Kkimura, Takehiro Himeno,
Hiroaki Kobayashi, Takayuki Kojima, Sigeru Aso, Yasuhiro Tani, Takakage Arai, “Development
Study on a Small-scale Supersonic Flight Experiment Vehicle with Jet Propulsion in Cooperation
between Universities”, 14th AIAA/AHI Space Plane and Hypersonic System and Technologies
Conference, AIAA-2006-7985, Canberra, Australia, 2006.

［9］

H. Kimura, S. Hada, K. Mizobata, R. Minato, N. Tanatsugu, “Flight Trajectory Analysis of a
Small-scale Supersonic Flight Experimental Vehicle”, 25th International Symposium on SpaceTechnology and Science, ISTS 2006-g-21, Kanazawa, 2006.

［10］ Minato, R., Tanatsugu, N., Mizobata, K., Kojima, T., Kobayashi, H., “Development of Turbojet
Engine for Supersonic Unmanned Plane”, 14th AIAA/AHI International Conference on Space
Planes and Hypersonic Systems and Technologies, AIAA Paper 2006-8145, Canberra, Australia,
2006.
［11］ Mizobata, K., Minato, R., Tanatsugu, N., Kimura, H., Hada, S., Kudo, M., et al., “Development
Study on a Small-scale Supersonic Flight Experiment Vehicle with Jet Propulsion in Cooperation
between Universities”, KSAS-JSASS Joint International Symposium on Aerospace Engineering,
Busan, Korea, 2006.
［12］ Minato, R., Ota, T., Fukutomi, K., Tanatsugu, N., Mizobata, K., Kojima, T., Kobayashi, H.
“Development of Turbojet Engine for Supersonic Unmanned Plane”, KSAS-JSASS Joint
International Symposium on Aerospace Engineering, Busan, Korea, 2006.
［13］ Minato, R., Ota, T., Fukutomi, K., Tanatsugu, N., Mizobata, K., Kojima, T., Kobayashi, H.,“Development of Counter Rotating Axial Fan Turbojet Engine for Supersonic Unmanned Plane”,
Joint Propulsion Conference, AIAA Paper 2007-5023, Cincinnati, America, 2007.
［14］ Minato, R., Himeno, T., Kojima, T., Kobayashi, H., Taguchi, H., Tetsuya Sato, T., Arai, T., Mizobata,K., Sugiyama, H., Tanatsugu, N. “Development of Counter Rotating Axial Fan Turbojet Engine forSupersonic Unmanned Plane at Muroran Institute of Technology”, International Gas Turbine
108

Congress , Tokyo, 2007.
［15］ K, Higashino. M, Sugioka. T, Kobayashi. R, Minato. Y, Maru. Y, Sasayama. M, Otsuka. T, Makino.
and H, Sakaguchi, “Fundamental Study on Coking Characteristics of LNG Rocket Engines,”
AIAA-2008-4753, AIAA/SME/ASME/ASEE 44th Joint Propulsion Conference & Exhibit, July
20-23, Hartford, 2008.
［16］ Watanabe, S. and Minato, R., “Development of a Design Support System that can Efficiently Utilize
Non-Dominated Solutions,” Joint 4th International Conference on Soft Computing and Intelligent
Systems and. 9th. International Symposium on advanced Intelligent Systems (SCIS & ISIS 2008),
CD-paper, 2008.
［17］ Minato, R., “Study of Subscale Turbojet Engine for Supersonic Unmanned Plane in MuroranInstitute of Technology”, International Symposium of Advanced Engineering, Pusan, 2009.
［18］ K, Higashino. M, Sugioka. T, Kobayashi. R, Minato. Y, Maru. Y, Sasayama. M, Otsuka. T, Makino.
and H, Sakaguchi, “Fundamental Study on Coking Characteristics of LNG Rocket Engines,”
AIAA-2008-4753, AIAA/SME/ASME/ASEE 44th Joint Propulsion Conference & Exhibit, July
20-23, 2008, Hartford.
［19］ Fukiba, K., Tsuboi, N, and Minato, R., 「Numerical Study on the Heat Transfer of theflow with Endothermic Chemical Reaction」, AIAA-2009-5461,AIAA Joint Propulsion
Conference 2009, Denver, 2009.
［20］ Fukiba, K., Tsuboi, N., Kobayashi, H. and Hongoh, M., 「Flow Angles and Mach
Number Measurement Using Surface Pressures on the Nose Cone 」 , 27th
International Symposium on Space Technology and Science, ISTS2009-g-15, Tsukuba,
Japan, July 2009.
［21］ Ryojiro Minato, Kazuyuki Higashino, Masatoshi Sugioka, Takao Kobayashi, Shunsuke,
Ooya, Yosuke, Sasayama (Muroran Instsitute of Technogy), “LNG Rocket Engine
with Coking Inhibited Regenerative Cooling System, ” 16th AIAA/DLR/DGLRInternational Space Planes and Hypersonic Systems and Technologies Conference，
AIAA-2009-7392(2009).
［22］ Yoshida, K., Sugiura, H., Ueda, Y., Ishikawa, H., Tokugawa, N., Atobe, A., Takagi, S., Arnal, D.,
Archambaud, Jean-Pierre and Seraude, A.:”Experimental and numerical research on boundary layertransition analysis at supersonic speed: JAXA-ONERA cooperative research project”, JAXA-RR08-007E, March 2009.
［23］ Y. Kanda, S.Ooya, T. Kobayashi, Y.Uemichi, K. Higashino, M. Sugioka, “Hydrogen production
reaction of aluminum with water using stainless steel reactor,” Joint Seminar on Environmental
Science and Disaster Mitigation Research 2010(JSED2010), Muroran Institute of Technology, 10
March 2010, Poster P-22B, Abstract book p.52-53.
［24］ Takashi Sakai and Shohei Takagi, “On the Strouhal number behind a series of symmetrical airfoils
at low Reynolds number.”, The Thailand-Japan International Symposium in Industrial Engineering,
Mechanical Engineering and Robotics 2010, Chiang Mai, Thailand (11/22-23, 2010).
［25］ Kazuhide MIZOBATA, Ryojiro MINATO, Katsuyoshi FUKIBA, Kazuyuki HIGASHINO and
Nobuhiro TANATSUGU, “Design of a Small-scale Supersonic Flight Experiment Vehicle as a
Flying Test Bed and Construction of its Prototype for Subsonic Flights,” 3rd Asia-PacificInternational Symposium on Aerospace Technology （ 14th Australian International Aerospace109

Congress）, Melbourne Convention Centre, Melbourne, Australia, Feb 28-March 3, 2011.
［26］ Y. Kanda, S.Ooya, T. Kobayashi, Y.Uemichi, K. Higashino, M. Sugioka, “Development of greenhydrogen production using aluminum and water,” Joint Seminar on Environmental Science andDisaster Mitigation Research 2011(JSED2011), Muroran Institute of Technology, 9 March
2011, Poster P-B-2, Abstract book p.62-63.
［27］ Kazuhide MIZOBATA, Ryojiro MINATO, Katsuyoshi FUKIBA, Kazuyuki HIGASHINO, andNobuhiro TANATSUGU, “Development of a Small-scale Supersonic Flight Experiment Vehicle as a
Flying Test Bed,” 17th AIAA International Space Planes and Hypersonic Systems and Technologies
Conference, San Francisco, CA, USA, 11-14 April 2011-2289［28］ Takagi S., Nishizawa A., and Itoh N., Dispersive Disturbances due to Cross-Flow and StreamlineCurvature Instabilities in 3-D Boundary Layers, Hawaii, June 4-7, 2011.(Invited)［29］ Yuichi Noguchi, Yasuhiro Ishikawa, Takashi Makino, Yosuke Sasayama, Kazuyuki
Higashino, “The Result of the Single Element Firing Tests of the Supercritical
LOX/methane Coaxial Injector,” ISTS, 2011［30］ Daisuke Nakata, Ami Kozu, Jun Yajima, Kenji Nishine, Kazuyuki Higashino andNobuhiro Tanatsugu, “Predicted and Experimented Acceleration Profile of the Rocket
Sled,” A-01 (oral) The 28th International Symposium on Space Technology and Science,
Okinawa, Japan, 2011［31］ Toru TSUKANO, “Thermal Energy Storage System for Stirling Cycle Electric
Generator in Orbit” q06s(oral) The 28th International Symposium on Space Technology
and Science, Okinawa, Japan, 2011［32］ Ken Goto, Tsuneo Kii, Keisuke Yoshihara, Kosei Ishimura, Keiji Komatsu, Ken Higuchi,
Taro Maruyama, Yuichi Ikeda, Kazuro Matsumoto, Makoto Yoshihara, Masaharu
Tsuchiya, “Long Term Durability of a Rib and Cable Tensioned Structure for HighAccuracy Large Deployable Antenna Reflector,” 28th International Symposium onSpace Technology and Science (ISTS), 2011-c-12, Ginowan,Okinawa, (2011.6.5-6.12).
［33］ Hiroshi Kawamura, Takashi Iwasa, Keita Unemitsu, Ryo Sugiyama, Naoko Kishimoto,Ken Higuchi, Motoharu Fujigaki, Takayuki Shiokawa, “Visualization on Formation
Process of Wrinkling Phenomena on Thin Membranes using Grating Projection Method
and Investigation on Their Nonlinear Mechanical Properties,” 28th International
Symposium on Space Technology and Science (ISTS), 2011-c-43, Ginowan, Okinawa,(2011.6.5-6.12).
［34］ Yoshiro Ogi, Ken Higuchi, Kosei Ishimura, “Effect of Attachment Errors of Flexible
Appendages on the Spin Axis of a Rigid Body,” 28th International Symposium on SpaceTechnology and Science (ISTS), 2011-c-43, Ginowan, Okinawa, (2011.6.5-6.12).
［35］ Kosei Ishimura, Tsuneo Kii, Keiji Komatsu, Ken Goto, Ken Higuchi, Kazuro Matsumoto,
Shoichi Iikura, Makoto Yoshihara, Masaharu Tsuchiya, “Shape Prediction of Large
Deployable Antenna Structure on Orbit,” 28th International Symposium on SpaceTechnology and Science (ISTS), 2011-c-45, Ginowan, Okinawa,(Okinawa Convention
Center) (2011.6.5-6.12).
［36］ Yasuharu Kanda, Yoshio Uemichi, Kazuyuki Higashino, Masatoshi Sugioka, “New hydrogenproduction by mechano-chemical reaction of aluminum with water,” The 8th Pacific Conference on
110Sustainable Energy & Environmental Technologies (APSCEET 2011), The University of Adelaide,
S.A. Australia（採択済み）, 10-13 July, 2011.
［37］ Minato, R, Kato, D. Higashino, K., Tanatsugu, N “Development Study on Counter
Rotating Fan Jet Engine for Supersonic Flight” ISABE 2011-1233, Goteburg, Sweden,
2011.
［38］ Minato, R, Higashino, K. and Tanatsugu, N. “Design and Performance Analysis of
Bio-Ethanol Fueled GG-cycle Air Turbo Ramjet Engine” AIAA Science Meeting 2012,
Nashville, Tennessee, USA 2012.
［39］ Daisuke Nakata, Jun Yajima, Kenji Nishine, Kazuyuki Higashino and Nobuhiro
Tanatsugu, “Research and Development of High Speed Test Track Facility in Japan,”
AIAA Science Meeting 2012, Nashville, Tennessee, USA 2012.
［40］ Masatoshi
Sugioka,Kazuyuki
Higashino,
YoshioUemichi,Yasuharu
Kanda, ”Production of Green Hydrogen by Mechanical Mixing of Aluminum with Water
using Stainless Steel Reactor” 14th Asia Pacific Confederation of Chemical Engineering
Congress (APCChE2012), 21-24 February, 2012 、Singapore［41］ Masatoshi Sugioka, Kazuyuki Higashino, Tadashi Momono,Yoshio Uemichi, andYasuharu Kanda, ”Development of Green Hydrogen Production in Muroran-IT”
JSED2012, 9 March, 2012，Muroran［42］ Ryojiro Akiba, Nobuhiro Tanatsugu, “A Vision of Future Space Transportation Systems”,
13th International Space Conference of Pacific-basin Societies(ISCOPS), Kyoto,(2012.5.15-18)［43］ K. Mizobata, R. Minato, K. Higashino, and N. Tanatsugu, “Development of a
Small-scale Supersonic Flight Experiment Vehicle as a Flying Test Bed,” 28 th Congressof the International Council of the Aeronautical Sciences, Brisbane, Australia,
September 24-28, 2012.
［44］ Masatoshi Sugioka, Kazuyuki Higashino, Yoshio Uemichi
and Yasuharu
Kanda, ”Production of Green Hydrogen by Reaction of Aluminum and Water” 19 thRegional Symposium of Chemical Engineering Congress (RSCE2012), 7-8 November,
2012，Bali，Indonesia．
［45］ Masatoshi Sugioka, Kazuyuki Higashino, Kohki Satoh and Yasuharu Kanda ， ”
Development of Green Hydrogen Production by Reaction of Aluminum and Water”
International Symposium on Chemical and Bioprocess Engineering 2013

国内学会発表論文［1］

棚次亘弘，「世界と我が国の宇宙輸送機開発の近況と，大学主導による超音速飛行実験計画」，
第 2 回 HASTIC 学術講演会，札幌，2005 年 3 月．

［2］

溝端 一秀，湊亮二郎，新井隆景，木村博幸，金子広樹，棚次亘弘，「超音速飛行実験のための小型ターボジェットエンジンと小型有翼機体の概念設計」，第 2 回 HASTIC 学術講演会，札幌，
2005 年 3 月．

［3］

棚次亘弘，溝端一秀，杉山 弘，「航空宇宙機システム研究センター」設立について」，日本航空
宇宙学会第 39 期年会講演会，調布，2005 年 4 月．
111

［4］

木村博幸，羽田尚太，溝端一秀，湊亮二郎，棚次亘弘，「小型超音速無人飛行機の飛行経路解
析」，第 49 回宇宙科学連合講演会，広島，2005 年 11 月．

［5］

溝端一秀，湊亮二郎，棚次亘弘，木村博幸，姫野武洋，小林弘明，小島孝之，麻生 茂，谷 泰
寛，新井隆景，「大学連携による小型超音速無人飛行機の開発構想」，第 49 回宇宙科学連合講演会，広島，2005 年 11 月．

［6］

谷 泰寛, 麻生 茂, 井上慶, 山口耕平, 伊藤睦高, 宮澤宏章，溝端一秀, 湊亮二郎, 棚次亘弘，
木村博幸, 羽田尚太, 本郷素行，「小型超音速無人飛行機の機体空力設計と風洞試験」，第 49回宇宙科学連合講演会，広島，2005 年 11 月．

［7］

溝端一秀，湊亮二郎，棚次亘弘，木村博幸，姫野武洋，小林弘明，小島孝之，麻生茂，谷泰寛，
新井隆景，「大学連携による小型超音速無人飛行機の開発構想」，第 3 回 HASTIC 学術講演会，札幌，2006 年 3 月．

［8］

湊亮二郎，戸倉郁夫，棚次亘弘，姫野武洋，小島孝之，小林弘明，「小型超音速無人飛行機のための小型ターボジェットエンジンの予備設計」，第 3 回 HASTIC 学術講演会，札幌，2006 年 3月．

［9］

木村博幸，羽田尚太，溝端一秀，湊亮二郎，棚次亘弘，谷泰寛，麻生茂，入門朋子，本郷素行，「小型超音速無人飛行機の機体空力設計，風洞試験，および飛行性能予測」，第 3 回 HASTIC
学術講演会，札幌，2006 年 3 月．

［10］ 羽田 尚太，「小型超音速飛行実験のための有翼機体の空力設計と風試結果」，平成 18 年度宇宙輸送シンポジウム，相模原，2006 年 1 月 18-19 日．
［11］ 棚次亘弘、「室蘭工業大学における航空宇宙工学の実践研究」、第５回 宇宙ミッションシンポジウ
ム、ＪＡＸＡ航空宇宙技術研究センター、平成１９年８月２８日［12］ 丸 祐介，工藤摩耶，笹山容資，桑田耕明，溝端一秀，坪井伸幸，「機体系および飛行システム技
術：小型超音速飛行実験の空力設計と軌道最適化」，第 51 回宇宙科学技術連合講演会，札幌，
2007 年 10 月 29-31 日．
［13］ 溝端 一秀，棚次 亘弘，東野 和幸，湊 亮二郎，「室蘭工大における航空宇宙機の基盤研究(3）
機体系および飛行システム技術)，第５１回宇宙科学技術連合講演会，札幌，2007 年 10 月 29-31
日．
［14］ 東野和幸，杉岡正敏，小林隆夫，湊亮二郎，丸祐介，笹山容資，大塚雅也，牧野隆，坂口裕之，「推進系基盤技術：LNGロケットエンジンのコーキング（炭素析出）特性に関する基盤研究」，第51回宇宙科学技術連合講演会，札幌，2007年10月29-31日．
［15］ 棚次亘弘，東野和幸，溝端一秀，湊亮二郎，「室蘭工大における航空宇宙機の基盤研究 （1）その意義と構想」，第51回宇宙科学技術連合講演会，札幌，2007年10月29-31日．
［16］ 東野和幸，棚次亘弘，湊亮二郎，「室工大における航空宇宙機の基盤研究（2）推進系基盤技術」，
第51回宇宙科学技術連合講演会，札幌，2007年10月29-31日．
［17］ 湊亮二郎，竹田広人，西村宗真，溝端一秀，東野和幸，棚次亘弘，太田豊彦，「小型超音速機用
反転軸流ファンの性能解析」，第51回宇宙科学技術連合講演会，札幌，2007年10月．
［18］ 丸 祐介，宮瀬 宗彦，安田 有佑，棚次 亘弘，溝端 一秀，「高速走行軌道装置構想と予備的走行実験」，平成 19 年度宇宙輸送シンポジウム，相模原，2007 年 1 月．
［19］ 工藤 摩耶，棚次 亘弘，溝端 一秀，丸 祐介，笹山 容資，桑田 耕明，新井 隆景，楠亀 拓也，
久保 良介，坪井 伸幸，「小型超音速飛行実験機の空力設計と空力性能評価」， 平成 19 年度
112

宇宙輸送シンポジウム，相模原，2007 年 1 月 28-29 日．
［20］ 棚次亘弘，「室蘭工業大学における航空宇宙工学の実践研究」，第 5 回 宇宙ミッションシンポジウ
ム，調布，2007 年 8 月．
［21］ 東野和幸、杉岡正敏、小林隆夫、湊亮二郎、丸 祐介、笹山容資、大塚雅也、「次世代型ロケット
用燃料としてのメタンの特性」、第８回北海道エネルギー資源環境研究発表会、札幌、2008 年 1
月.
［22］ 東野和幸、杉岡正敏、日比野政裕、小林隆夫、笹山容資、大屋俊輔、立川藍子、三輪峻之、「次世代型ロケット用燃料としての LNG の特性」、第９回北海道エネルギー資源環境研究発表会、
北海道大学 学術交流会館、平成２０年２月３日［23］ 溝端一秀，東野和幸，湊亮二郎，新井隆景，「FTB としての小型超音速飛行実験機の構想」，日
本航空宇宙学会第 39 期年会講演会，調布，2008 年 4 月．
［24］ 東野和幸，杉岡正敏，小林隆夫，境昌宏，湊亮二郎，笹山容資，大塚雅也，沖田耕一，青木賢司，
川島秀人，東伸，「LNG ロケットエンジンにおけるサルファアタック・コーキングに関する基礎研究」，
第 52 回 宇宙科学技術連合講演会，淡路島，2008 年 11 月 5-7 日．
［25］ 東野和幸、駒崎慎一、笹山容資（室工大）、木元健一(IHI)、坂口裕之(IHIエアロスペース)、石橋
利幸(NEGENT)、「ロケット燃焼室におけるニッケル電鋳特性向上に関する基礎研究」、第52回
宇宙科学技術連合講演会、淡路島、２００８年１１月５日～７日［26］ 湊亮二郎，溝端一秀，桑田耕明，「超音速風洞試験におけるスターティングロードの計測と推算”，
第 52 回宇宙科学技術連合講演，淡路，2008 年 11 月．
［27］ 吹場活佳，坪井伸幸，小林弘明，「円錐状ノーズコーンをもつ超音速機の壁面圧力計測による機
体姿勢の算出」，平成 20 年度宇宙航行の力学シンポジウム，相模原，2008 年 12 月．
［28］ 坪井伸幸，吹場活佳，“前処理法による非定常圧縮性低速流れの数値解析”，宇宙航行の力学シンポジウム，相模原，2008 年 12 月．
［29］ 坪井伸幸，吹場活佳，“超臨界燃焼流れにおける数値解析、「低速高レイノルズ数に対する解析とその課題」，平成 20 年度科学研究費空力班シンポジウム，嬉野，2009 年 1 月［30］ 溝端一秀，湊亮二郎，吹場活佳，東野和幸，棚次亘弘，「FTB としての小型超音速飛行実験機の構想と予備的な亜音速機の試作」，平成 20 年度宇宙輸送シンポジウム，相模原，2009 年 1 月．
［31］ 吹場活佳，前田大輔，磯田浩志，東野和幸，湊亮二郎，「吸熱反応を利用した再生冷却に関する基礎研究 －室蘭工業大学における取り組み－」，平成 20 年度宇宙輸送シンポジウム，相模原，2009 年１月．
［32］ 桑田耕明，笹山容資，飯村拓哉，渡辺侑也，棚次亘弘，溝端一秀，吹場活桂，坪井伸幸，「小型無人超音速有翼機の横および方向の姿勢安定と操縦性に関する風洞試験」，平成 20 年度宇宙輸送シンポジウム，相模原，2009 年１月．
［33］ 坪井伸幸、吹場活佳、「超臨界燃焼流れにおける数値解析：低速高レイノルズ数に対する解析とその課題」、平成 20 年度科学研究費空力班シンポジウム、嬉野、2009 年 1 月［34］ 東野和幸，杉岡正敏，日比野政裕，小林隆夫，笹山容資，大屋俊輔，立川藍子，三輪峻之，「次世代型ロケット用燃料としての LNG の特性」，第 9 回北海道エネルギー資源環境研究発表会，札幌，2009 年 2 月．
［35］ 坪井伸幸，吹場活佳，嶋田徹，「前処理法による非定常圧縮性低速流れの数値解析：ハイブリッド
113

ロケット燃焼器内流れの解析に向けて」，第 49 回航空原動機・宇宙推進講演会，長崎，2009 年 3月．
［36］ 安田有佑，溝端一秀，棚次亘弘，神津亜実，「サブスケール高速走行軌道装置を用いた予備的走行実験」，第 6 回 HASTIC 学術講演会，札幌，2009 年 3 月．
［37］ 吹場活佳，坪井伸幸，小林弘明，“気球を用いた微小重力実験機による超音速飛行実験と機体姿勢計測”，第 6 回 HASTIC 学術講演会，札幌，2009 年 3 月．
［38］ 髙木正平：「空力音と流れの不安定性」、第４１回流体力学講演会/航空宇宙数値シミュレーション
技術シンポジウム（基調講演）、6 月 18-19 日(2009), 宇宙航空研究開発機構.
［39］ 東野和幸、杉岡正敏、小林隆夫、湊 亮二郎、笹山容資、大屋俊輔（室蘭工大）、「次世代型ロケット燃料としてのＬＮＧの化学的特性」、北海道支部夏季研究発表会（2009-7）［40］ 笹山容資、杉岡正敏、小林隆夫、境昌宏、東野和幸（室蘭工大）、東伸幸、青木賢司、沖田耕一（ＪＡＸＡ）：ＬＮＧロケットエンジンにおけるサルファアタック防止に関する研究－特に金メッキの特性評価”，第５３回宇宙科学技術連合講演会（2009-9）［41］ 東野和幸、杉岡正敏、小林隆夫、湊 亮二郎、大屋 俊輔、笹山容資（室蘭工大）：グラファイトコー
ティングによる LNG ロケットエンジンの再生冷却管内のコーキング抑制”，第５３回宇宙科学技術
連合講演会（2009-9）［42］ 成尾芳博、徳留真一郎、稲谷芳文（ＪＡＸＡ）、東野和幸（室蘭工大）、「高度制御／着陸誘導可能な加圧供給方式エンジンの開発」， 第５３回宇宙科学技術連合講演会（2009-9）［43］ 塚野徹，吹場活佳，棚次亘弘（室蘭工大）、「宇宙用スターリングサイクル発電機に関する実験的研究」，第５３回宇宙科学技術連合講演会（2009-9）［44］ 髙木正平、「二次元物体後流の安定性と空力音」、平成２１年度第１回流体工学研究会・懇話会（依頼講演）（北大工学部）、（11/27/2009）．
［45］ 桑田耕明、飯村拓哉、酒井貴志、溝端一秀、中田大将、吹場活佳、高木正平、東野和幸、棚次亘
弘（室蘭工大）、丸 祐介（JAXA）、「小型無人超音速機において エンジン流量特性が全機空力特性に及ぼす影響」，平成 21 年度宇宙輸送シンポジウム（2010-1）［46］ 安田有佑、神津亜実、矢島淳、中田大将、棚次亘弘、吹場活佳（室蘭工大）、「高速走行軌道装
置の加減速に関する研究」，平成 21 年度宇宙輸送シンポジウム（2010-1）［47］ 石川昂紀、小林隆夫、神田康晴、桃野 正、東野和幸、杉岡正敏(室工大)、「アルミニウム合金と水との反応による水素製造」、第１０回北海道エネルギー資源環境研究発表会（2010-1）［48］ 笹山容資、湊 亮二郎、境 昌宏、小林隆夫、杉岡正敏、東野和幸（室蘭工大）、「LNG ロケットエ
ンジンにおけるサルファタック及びコーキングに関する研究」、第７回 HASTIC 学術講演会、北大、2010 年 3 月 4 日［49］ 前田大輔、笹山容資、湊 亮二郎、小林隆夫、杉岡正敏、東野和幸、棚次亘弘（室蘭工大）「ＥＦとしてのメチルシクロヘキサンの熱分解特性」、第７回 HASTIC 学術講演会、北大、2010 年 3 月 4日［50］ 高津武人、神津亜実、安田有佑、中田大将、溝端一秀、棚次亘弘（室蘭工大）「超音速小型実験機の研究開発と高速走行試験装置について」、第 7 回 HASTIC 学術講演会、北大、2010 年 3月4日［51］ 棚次亘弘、「大気中を高速度で飛行するための基盤技術の研究開発」、日本航空宇宙学会北部114

支部講演会、JSASS-2010-H053、東北大学、2010 年３月 1１日［52］ 片山範将，石村康生，宮崎康行，樋口 健、「マルチセルインフレータブルサンドイッチパネルの曲
げ特性に関する研究」, 第 52 回構造強度に関する講演会, 1A17, (JSASS-2010-3020), 鳥取,
pp.49-51, (2010.7).
［53］ 塩路義行, 樋口 健，岸本直子，藤垣元治，塩川貴之，森本吉春、「格子投影法によるメッシュ膜
面の面外形状の計測」, 第 52 回構造強度に関する講演会, 1B06, (JSASS-2010-3024), 鳥取,
pp.61-63, (2010.7).
［54］ 樋口 健，青木隆平，宮崎康行，古谷 寛，泉田 啓，岸本直子，角田博明，石村康生，石澤淳一
郎，酒井良次，渡邊秋人，川端信義，堀 利行，伊藤裕明，渡辺和樹，及川 祐，倉冨 剛、「インフレータブル方式伸展マスト宇宙実証のための「きぼう」曝露部搭載実験装置の開発」, 第 52 回
構造強度に関する講演会, 1A18, (JSASS-2010-3021), 鳥取, pp.52-54, (2010.7).
［55］ 荻 芳郎，樋口 健，石村康生：「剛体スピン軸への梁の取り付け誤差の影響」, 第 52 回構造強度に関する講演会, 3A04, 鳥取, pp.179-181, (2010.7).
［56］ 青木隆平, 古谷 寛, 樋口 健,宮崎康行, 石村康生, 石澤淳一郎, 泉田 啓, 岸本直子, 角田博明, 酒井良治, 渡邉秋人, 川端信義, 堀 利行, 伊藤裕明,渡辺和樹,及川 祐, 倉富 剛：「宇宙インフレータブル構造の宇宙実証－SIMPLE－」, S1903-1-6, 2010 年度日本機械学会年次大会, 名古屋, (2010.9).
［57］ 伊藤信毅、髙木正平、池田友明、「平板後流の不安定と振動数選択機構」、日本流体力学会年会
2010、(9/9-11、 2010).
［58］ 田中宏明, 森島駿一, 岡田侑樹, 池田忠繁, 坂本 啓, 秋田 剛, 荻 芳郎, 石村康生, 樋口 健、「スマート構造システムの低温低圧環境試験」, 第 54 回宇宙科学技術連合講演会, 1K06, 静岡,(2010.11).
［59］ 神津 亜実、中田 大将、矢島 淳、棚次 亘弘、「2L11 高速走行軌道実験装置における水ブレーキの最適化」、第５４回宇宙科学技術連合講演会、静岡県コンベンションアーツセンター、2010
年１１月１７－１９日［60］ 塚野 徹、吹場活佳、棚次亘弘、「3B09 軌道上での運用を目的としたスターリングサイクル発電機の研究」、第５４回宇宙科学技術連合講演会、静岡県コンベンションアーツセンター、2010 年１１
月１７－１９日［61］ 東伸幸、森谷信一、只野真、佐藤政裕、青木賢司、小林完、沖田耕一（宇宙航空研究開発機構）、
水野孝彦、丹生謙一（三菱重工）、笹山容資、東野和幸（室蘭工大）、「3H02 LOX/メタン再生冷却エンジンにおける金めっきによるサルファアタック影響抑制効果」、第５４回宇宙科学技術連合
講演会、静岡県コンベンションアーツセンター、2010 年１１月１７－１９日［62］ 東野和幸、棚次亘弘、杉岡正敏、湊 亮二郎、笹山容資、磯田浩志、「3H04 メチルシクロヘキサンの熱分解吸熱反応による冷却効果について」、第５４回宇宙科学技術連合講演会、静岡県コンベンションアーツセンター、2010 年１１月１７－１９日［63］ 石川昂紀、小林隆夫、神田康晴、杉岡正敏、東野和幸、「3H07 アルミ-水系反応による高圧水素の発生と宇宙機推進システムへの適用」、第５４回宇宙科学技術連合講演会、静岡県コンベンショ
ンアーツセンター、2010 年１１月１７－１９日［64］ 笹山容資、杉岡正敏、東野和幸（室工大）、東伸幸、平岩徹夫、沖田耕一（JAXA）、「3H03 バイ
オエタノールロケットエンジンの材料適合性に関する実験的研究」、第５４回宇宙科学技術連合講
演会、静岡県コンベンションアーツセンター、2010 年１１月１７－１９日
115

［65］ 溝端一秀、 湊亮二郎、 東野和幸、 棚次亘弘、「フライングテストベッドとしての小型超音速飛行
実験機およびそのプロトタイプの設計と試作、第 48 回飛行機シンポジウム、静岡県コンベンション
アーツセンター、2010 年１１月３０日－１２月２日［66］ 吹場活佳、 上村卓也、 溝端一秀、 湊亮二郎、高木正平、 東野和幸、 棚次亘弘、「室蘭工業大学吸い込み式超音速風洞の気流特性」、第 48 回飛行機シンポジウム、静岡県コンベンション
アーツセンター、2010 年１１月３０日－１２月２日［67］ 采光啓太, 岩佐貴史, 川村寛志, 樋口 健, 岸本直子, 藤垣元治, 塩川貴之、「格子投影法による皺の生じた膜面形状の載荷経路依存性に関する実験的検証」, 第 26 回宇宙構造・材料シンポ
ジウム, B19, 相模原, (2010.12).
［68］ 田中宏明, 森島駿一, 岡田侑樹, 池田忠繁, 坂本 啓, 秋田 剛, 荻 芳郎, 石村康生, 南部陽介, 樋口 健：「スマート構造システムの低温低圧環境における有効性実証試験」, 第 26 回宇宙構造・材料シンポジウム, 相模原, (2010.12).
［69］ 青木隆平, 樋口 健, 石澤淳一郎, 宮崎康行, 古谷 寛，石村康生, 角田博明, 岸本直子, 泉田
啓, 酒井良次, 渡邊秋人, 川端信義, 堀 利行, 伊藤裕明, 渡辺和樹, 及川 祐, 倉富 剛：「SIMPLE 膜構造宇宙実験の概要」, 第 26 回宇宙構造・材料シンポジウム, 相模原, (2010.12).
［70］ 小嶋浩嗣, 笠羽康正, 八木谷聡, 笠原禎也, 石坂圭吾, 熊本篤志, 樋口 健、「次期磁気圏観測衛星検討 WG SCOPE 計画 －プラズマ波動観測－」, 第 11 回宇宙科学シンポジウム,
P3-060,相模原, (2011.1).
［71］ 小 松 敬 治 , 紀 伊 恒 男 , 後 藤 健 , 樋 口 健 , 石 村 康 生 , 吉 原 眞 , 飯 倉 省 一 , 松 元 和 郎 ,ASTRO-G 技術実証チーム：「ASTRO-G LDR の EM 試験」, 第 11 回宇宙科学シンポジウム,
P2-043, 相模原, (2011.1).
［72］ 石村康生, 紀伊恒男, 小松敬治, 後藤 健, 樋口 健, 村田泰宏, 岩田隆浩, 朝木義晴, 土居明
宏, 萩野慎二, 土屋正治, 入門 寛, 隈下恭介, 吉原 眞, 松元和郎, ASTRO-G 技術実証チー
ム：「ASTRO-G LDR の総合評価」, 第 11 回宇宙科学シンポジウム, P2-042, 相模原, (2011.1).
［73］ 中田大将、神津亜実、矢島淳、西根賢治、東野和幸、棚次亘弘、「ハイブリッドロケット搭載高速走
行軌道実験設備」 、平成 22 年度宇宙輸送シンポジウム、相模原、2011 年 1 月 20-21 日［74］ 飯村拓哉、吹場活佳、上村卓也、「遷音速領域における ADS の計測特性に関する研究」、平成
22 年度宇宙輸送シンポジウム、相模原、2011 年 1 月 20-21 日［75］ 溝端一秀、湊亮二郎、吹場活佳、東野和幸、棚次亘弘、「フライングテストベッドとしての小型超音速飛行実験機の設計、およびプロトタイプの試作と試験飛行」、平成 22 年度宇宙輸送シンポジウ
ム、相模原、2011 年 1 月 20-21 日［76］ 前田大輔、笹山容資、湊亮二郎、小林隆夫、杉岡正敏、棚次亘弘、東野和幸、「熱分解吸熱反応
燃料としてのメチルシクロヘキサンの分解特性」、第 11 回北海道エネルギー資源環境研究発表会、札幌、2011 年 2 月.
［77］ 石川昂紀、小林隆夫、神田康晴、桃野 正、東野和幸、杉岡正敏、「アルミニウム合金と水との反応による高圧水素の製造」第１１回北海道エネルギー資源環境研究発表会、札幌、2011 年 2 月.
［78］ 東野和幸、湊亮二郎、中田大将、吹場活佳、笹山容資、牧野隆、石川康弘、野口裕一、「LOX/LNG 高圧ロケットエンジン燃焼器噴射エレメントに関する実験的研究」第５１回航空原動
機・宇宙推進講演会、広島、2011 年 3 月 3-4 日［79］ 中田 大将、神津 亜実、矢島 淳、西根 賢治、東野 和幸、棚次 亘弘、「ハイブリッドロケットで加速される高速走行軌道実験設備」第 8 回 HASTIC 学術講演会、札幌、2011 年 3 月 4 日
116

［80］ 湊亮二郎，溝端一秀，吹場活佳，東野和幸，棚次亘弘「室蘭工大の小型超音速飛行実験機の設計とプロトタイプ機体の飛行試験」、第8回HASTIC学術講演会、札幌、2011年3月4日［81］ 川村寛志，岩佐貴史，樋口 健，岸本直子，藤垣元治，塩川貴之、「格子投影法による皺の生じた膜面の全視野形状計測」，日本機械学会中国四国支部第 49 期総会講演会，岡山, 2011 年 3 月.
［82］ 髙木正平、山谷直広、「自励による噴流不安定波の周波数選択」、第 48 回「乱流遷移の解明と制御｣研究会、日本大学理工学部駿河台校舎、2011 年 3 月 29-30 日.
［83］ 加藤 大貴、髙木 正平、「低レイノルズ数における角柱からの渦放出特性」、日本航空宇宙学会
第４２期年会講演会、東京、2011 年 4 月 14-15 日．
［84］ 神田康晴、小林隆夫、上道芳夫、東野和幸、杉岡正敏、「アルミニウムと水との反応による水素製造とその応用、日本化学学会北海道支部、室蘭、２０１１年 7 月２３日．
［85］ 髙木正平、山谷直広、伊藤信毅：「人工音響ループによる二次元噴流不安定波周波数の選択」、
日本流体力学会年会 2011, (9/7-9/9/2011).
［86］ 横山慶輔, 髙木 正平、「縮流胴出口境界層の高周波変動の発生機構」、日本流体力学会年会
2011、(9/7-9/9, 2011).
［87］ 髙木正平：「さらに簡単化した自作 CTA 回路」、日本機械学会 2011 年度年次大会（依頼講演）（9/11-9/14/2011）．
［88］ 髙木正平、加藤大貴：「二次元翼後流における時間発展攪乱の観察と周波数選択機構」、第 49
回「乱流遷移の解明と制御」研究会, (9/22-9/23/2011).
［89］ 笹山容資（室蘭工業大学大学院）、東野和幸（室蘭工業大学）、牧野隆、石川康弘（株式会社 Ｉ
ＨＩ）、野口裕一（株式会社 ＩＨＩエアロスペース）「２Ｋ１１ バイオエタノールロケットエンジン燃焼器
エレメントに関する実験的研究」、第５５回宇宙科学技術連合講演会、愛媛県県民文化会館、
2011 年１１月 30－12 月 2 日．
［90］ 東野和幸（室蘭工業大学）、笹山容資（室蘭工業大学大学院）、杉岡正敏（室蘭工業大学）、笹山容資（室蘭工業大学大学院）「2K12 バイオエタノールロケットエンジンの高温高圧での材料適合性に関する実験的研究」、第５５回宇宙科学技術連合講演会、愛媛県県民文化会館、2011 年１１
月 30－12 月 2 日．
［91］ 前田大輔、笹山容資（室蘭工業大学大学院）、杉岡正敏、東野和幸（室蘭工業大学）「３Ｃ０９ 空
気吸込式エンジン冷却システムに用いる熱分解吸熱反応燃料に関する実験的研究（金属の触媒効果について）」、第５５回宇宙科学技術連合講演会、愛媛県県民文化会館、2011 年１１月 30－
12 月 2 日．
［92］ 石村康生, 紀伊恒男, 小松敬治, 後藤 健, 樋口 健, 村田泰宏, 岩田隆浩, 朝木義晴, 土居明
宏, 萩野慎二, 土屋正治, 入門 寛, 隈下恭介, 吉原 眞, 松元和郎, ASTRO-G 技術実証チー
ム：ASTRO-G LDR の総合評価, 第 11 回宇宙科学シンポジウム, P2-042, 相模原, (2011.1).
［93］ 川村寛志，岩佐貴史，樋口 健，岸本直子，藤垣元治，塩川貴之：格子投影法による皺の生じた膜面の全視野形状計測，日本機械学会中国四国支部第 49 期総会講演会，岡山, (2011.3.5).
［94］ 片山範将，石村康生，二橋勇気，杉山彩香，青木隆平，樋口 健，宮崎康行，岸本直子：マルチ
セルインフレータブルサンドイッチパネルの開発，第 53 回構造強度に関する講演会，1A12,(2011.7)，（秋田）［95］ 小野寺佑介，樋口 健，渡邊秋人：開断面一次元伸展構造物の力学特性取得実験，第 53 回構
造強度に関する講演会，1A17, JSASS-2011-53050, (2011.7)，（秋田）117［96］ 石村康生，梶川隆史，川端信義，渡邊秋人，尾崎毅志，樋口健，荻芳郎，渡辺和樹，笠羽康正，
小嶋浩嗣，斎藤義文，酒井良次：三軸織 STEM のスピン時における動的特性評価，日本機械学会 2011 年度年次大会，J191045，（東京），(2011.9).
［97］ 川村寛志，岩佐貴史，岸本直子，樋口 健，藤垣元治，塩川貴之：格子投影法を用いた膜面に生じる皺の減少過程の形状計測第 48 回日本航空宇宙学会関西・中部支部合同秋期大会，A33，(2011.11), （神戸）［98］ 三輪武史，樋口 健，藤垣元治，塩川貴之，岩井達也，似鳥透：格子投影法における 外挿法の
適 用 と 面 計 測 ， JSASS-2011-4441, 第 55 回 宇 宙 科 学 技 術 連 合 講 演 会 講 演 集 , 3A04,(2011.12),（松山）．
［99］ 川村寛志，岩佐貴史，岸本直子，樋口 健，藤垣元治，塩川貴之：格子投影法を用いた薄膜に生じる皺の挙動計測日本実験力学会(JSEM)，分科会合同ワークショップ 2011，(2011.12), （米子）．
［100］樋口 健，三輪武史，岩井達也，似島 透，岩佐貴史，岸本直子：格子投影法外挿法における計測精度向上の試み，第 27 回宇宙構造・材料シンポジウム，B4, （相模原），(2011.12).
［101］樋口 健，三輪武史，岩井達也，似島 透，岩佐貴史，岸本直子：小型凹面鏡の格子投影法による形状計測，第 27 回宇宙構造・材料シンポジウム，B5,（相模原），(2011.12).
［102］中田大将，矢島淳，西根賢治，森木嵩人，ムハマドファクルラー，東野和幸，棚次亘弘，”高速走行軌道実験設備の基盤技術について,” 平成 23 年度宇宙輸送シンポジウム，2012 年 1 月，相模原市
［103］矢島淳，中田大将，棚次亘弘，”高速走行軌道実験装置の摩擦低減に関する研究,” 平成 23 年
度宇宙輸送シンポジウム，2012 年 1 月，相模原市［104］友常雄太郎，姫野武洋，渡辺紀徳，中田大将，棚次亘弘，東野和幸，”高速走行軌道実験の水制動特性に関する数値解析,” 平成 23 年度宇宙輸送シンポジウム，2012 年 1 月，相模原市［105］北沢祥一、阿野進、伴弘司、上羽正純、“環境電波のハーベスティングに関する一検討”、電子情報通信学会マイクロウエーブ研究会、MW2011-83（2011-9）, pp.27-32［106］伴 弘 司 、中 本 成 洋 、岡 智 広 、上 羽 正 純 、 “ 狭 小 な ICT 機 器 内 空 間 に お ける 電 波 伝 搬 ”、Microwave Workshops & Exhibition 2011(MWE2011)、 2011.11.30-12.02、 横浜［107］牟田龍平， 東野和幸，杉岡正敏，笹山容資，バイオエタノールロケットエンジン冷却における熱分解特性に関する研究，平成 23 年度宇宙輸送シンポジウム，2012 年 1 月 19 日～20 日，相模原．
［108］溝端一秀，湊亮二郎，東野和幸，棚次亘弘，「フライングテストベッドとしての小型超音速飛行実験機の設計とプロトタイプ機の飛行試験」，第 49 回飛行機シンポジウム，2011 年 10 月 26 日～
28 日，金沢．
［109］溝端―秀，湊亮二郎，東野和幸，棚次亘弘，新井隆景，「羹蘭工大の小型超音速飛行実験機の空力設計と空力特性評価」，STCP-2011-069，平成 23 年度宇宙輸送シンポジウム，2012 年 1
月 19 日～20 日，相模原．
［110］渡邉 力夫（東京都市大）、中田 大将（室工大）、東野 和幸（室工大）、 棚次 亘弘（室工大）、室蘭工業大学における水／液体窒素ロケット軌道上稼働試験について、第 20 回スペース・エンジ
ニアリング・コンファレンス(SEC'11)、日本機械学会、2012 年 1 月 26-27 日［111］笹山容資（室蘭工大院），東野和幸(室蘭工大)，クリーンロケット燃料の燃焼特性と冷却特性に関118

する研究，日本航空宇宙学会北部支部 2011 年講演会，2012，室蘭．
［112］牟田龍平，笹山容資（室蘭工大院），東野和幸，杉岡正敏（室蘭工大），バイオエタノールロケット
エンジンの準静的環境下における熱分解特性と再生冷却システムに関する研究，日本航空宇宙学会北部支部 2011 年講演会，2012，室蘭．
［113］寺田利幸，泉俊太郎，笹山容資（室蘭工大院），杉岡正敏，東野和幸（室蘭工大），東伸幸，平岩哲夫（ＪＡＸＡ），バイオエタノールロケットエンジンに対する材料適合性に関する研究，日本航空宇宙学会北部支部 2011 年講演会，2012，室蘭．
［114］溝端一秀，湊亮二郎，東野和幸，棚次亘弘，フライングテストベッドとしての小型超音速飛行実験機の設計とプロトタイプ機の飛行試験，第 49 回飛行機シンポジウム，2011 年 10 月 26 日～ 28
日，金沢．
［115］溝端―秀，湊亮二郎，東野和幸，棚次亘弘，新井隆景，室蘭工大の小型超音速飛行実験機の空
力設計と空力特性評価，STCP-2011-069，平成 23 年度宇宙輸送シンポジウム，2012 年 1 月 19
日～20 日，相模原．
［116］江上 侑雄(HASTIC)、棚次 亘弘、秋葉 鐐二郎(HASTIC), 無重力実験装置 Ez-SpaceⅤの開発、第９回 HASTIC 学術講演会、2012 年３月 13 日、北大百年記念会館［117］中田大将(室蘭工大)，矢島淳,西根賢治（室蘭工大院），森木崇人，ムハマドﾞファクラー,東野和幸,
棚次亘弘(室蘭工大)，高速走行軌道実験設備の展望，第９回 HASTIC 学術講演会、2012 年３
月 13 日、北大百年記念会館［118］泉俊太郎，笹山容資，寺田利幸（室蘭工大院），牟田龍平，東野和幸，杉岡正敏（室蘭工大），東伸幸，平岩哲夫（ＪＡＸＡ），バイオエタノールロケットエンジンにおけるＣＦＲＰの材料適合性に関する評価，第９回 HASTIC 学術講演会、2012 年３月 13 日、北大百年記念会館［119］三輪武史，岩井達也（室蘭工大院），似島 透，樋口 健（室蘭工大），格子投影法における外挿法の提案と計測精度評価，第９回 HASTIC 学術講演会、2012 年３月 13 日、北大百年記念会館［120］岩井達也，三輪武史（室蘭工大院），似島 透，樋口 健（室蘭工大），大型宇宙構造物の面形状計測法としての格子投影法，第９回 HASTIC 学術講演会、2012 年３月 13 日、北大百年記念会
館［121］鷹取一哉（室蘭工大院），金谷良平，谷川 豊，石田貴大，樋口 健（室蘭工大），無人小型超音速実験機（オオワシ１号機）の着陸衝撃吸収脚，第９回 HASTIC 学術講演会、2012 年３月 13 日、
北大百年記念会館［122］塚野徹，佐藤紘之，棚次亘弘（室蘭工大）, 宇宙用フリーピストン形スターリングサイクル発電機の負荷変動による発電特性、第９回 HASTIC 学術講演会、2012 年３月 13 日、北大百年記念会館［123］近藤光輝，笹山容資（室蘭工大院），東野和幸，杉岡正敏（室蘭工大），宇宙機推進システムとしての Al/水反応を利用した高圧水素製造に関する研究，日本航空宇宙学会北部支部 2011 年講演会，2012，室蘭．
［124］山口研太，杉岡正敏，東野和幸，塚野徹(室蘭工大), 熱分解吸熱反応燃料の触媒脱水素反応
特性に関する研究，日本航空宇宙学会北部支部 2011 年講演会，2012，室蘭．
［125］佐藤紘之，塚野徹，棚次亘弘(室蘭工大)，電気的な負荷変動によるフリーピストン形スターリング
サイクル発電機の特性，日本航空宇宙学会北部支部 2011 年講演会，2012，室蘭．
［126］西根賢治，森木崇人，中田大将，東野和幸,棚次亘弘(室蘭工大)，ロケットスレッド搭載用標準模型の設計と天秤装置の検証，日本航空宇宙学会北部支部 2011 年講演会，2012，室蘭．
119

［127］中田大将，矢島淳,西根賢治，森木崇人，ムハマドﾞファクラー,東野和幸,棚次亘弘(室蘭工大)，高速走行軌道実験設備の基盤研究，日本航空宇宙学会北部支部 2011 年講演会，2012，室蘭．
［128］溝端一秀，湊亮二郎，東野和幸，棚次亘弘，小型超音速飛行実験機「オオワシ」の開発と予備的飛行実験，日本航空宇宙学会北部支部 2011 年講演会，2012，室蘭．
［129］湊亮二郎，東野和幸，棚次亘弘(室蘭工大)，バイオエタノール燃料を用いたガスジェネレータサイ
クル・ＡＴＲエンジンの開発，日本航空宇宙学会北部支部 2011 年講演会，2012，室蘭．
［130］髙木正平、坂上昇史：「室蘭工大における遷移研究の２，３の話題」、第 50 回「乱流遷移の解明と制御」研究会, (3/29-3/30/2012).
［131］岩井達也，似鳥 透，樋口 健，三輪武史：格子投影法による大型構造物の平面形状計測，日本
航空宇宙学会・日本機械学会共催，第 54 回構造強度に関する講演会，3A01,
JSASS-2012-3065, pp.185-187, (2012.8.1-3)，（熊本市国際交流会館）［132］崎海里，青田 一輝，岩佐貴史，樋口 健，岸本直子：格子投影法を用いた全視野形状計測時のノイズ特性につて，日本機械学会 2012 年次大会，J191032，金沢大学，2012.9.9－12.
［133］手塚亜聖（早大理工）、髙木正平；三角柱後流でみられる振動流の各位相に加えられた擾乱の安定解析(Stability analysis for disturbances added at each phase of oscillatory flow
observed in the wake of triangular cylinder)、第 50 回飛行機シンポジウム、新潟トキメッセ(11/5-7, 2012)［134］溝端一秀、湊亮二郎、東野和幸、棚次亘弘、髙木正平、樋口健、上羽正純、中田大将：室蘭工大の小型超音速飛行実験機（オオワシ）の開発の現状と展望，第 50 回飛行機シンポジウム，特殊航空機部門提案の企画「無人機の新技術と多面的応用」，新潟市朱鷺メッセ，2012.11.5-7.
［135］塚野 徹（室蘭工業大学大学院）、東野和幸、杉岡正敏（室蘭工業大学）「３Ｈ１８ 白金触媒を用いた熱分解吸熱反応燃料の吸熱量評価と再生冷却システムに関する研究」、第５６回宇宙科学技術連合講演会、大分県、２０１２年１１月２0－1１月２２日．
［136］泉俊太郎（室蘭工業大学大学院）、杉岡正敏、東野和幸（室蘭工業大学）「２Ｈ０４ バイオエタノ
ールロケットエンジンにおけるアルミニウムの材料適合性に関する研究」、第５６回宇宙科学技術連
合講演会、大分県、２０１２年１１月２0－1１月２２日．
［137］湊 亮二郎（室蘭工業大学）、東野和幸、棚次亘弘（室蘭工業大学航空宇宙機システム研究センター）、「小型無人超音速機用ガスジェネレータサイクルＡＴＲエンジンの開発」、第５６回宇宙科学技術連合講演会、大分県、２０１２年１１月２0－1１月２２日．
［138］渡辺豪、渡邊力夫（東京都市大学）、中田大将、東野和幸、棚次亘弘（室蘭工業大学）、「２Ｈ０９
高速走行軌道装置を用いた水／液体窒素ロケットの稼働試験」、第５６回宇宙科学技術連合講演会、大分県、２０１２年１１月２0－1１月２２日．
［139］上羽 正純、溝端一秀（室蘭工業大学）「１Ｈ１０ 小型無人超音速実験機の着陸制御の一検討」、
第５６回宇宙科学技術連合講演会、大分県、２０１２年１１月２0－1１月２２日［140］似鳥 透，岩井達也，樋口 健，相原弘匡，村瀬祥宏：格子投影法外挿法における面形状計測座
標の幾何算出法の適用，第 56 回宇宙科学技術連合講演会講演集 1O10，JSASS-2012-4286，
別府国際コンベンションセンター（ビーコンプラザ），2012.11.20～22.
［141］樋口 健，岩佐貴史，岸本直子，岩井達也，似鳥 透：２カメラを用いた格子投影法による表面形
状計測，第 28 回宇宙構造材料シンポジウム，ISAS/JAXA 相模原，2012.12.4.
［142］岩佐貴史，石村康生，田中宏明，樋口 健，岸本直子，古谷 寛，神谷友裕，小木曽望，目黒 在，120

泉田 啓，荻 芳郎，中篠恭一，坂本 啓，仙場淳彦，佐藤泰貴，松尾雄太郎：点計測と面計測を統合した大型宇宙構造システムの高精度形状計測法の提案，第 28 回宇宙構造材料シンポジウ
ム，ISAS/JAXA 相模原，2012.12.4.
［143］江上郁雄（HASTIC）、棚次亘弘、秋葉鐐二郎(HASTIC)、Ez-SpaceVI、第 10 回 HASTIC 学術
講演会、北海道大学-JAXA 連携宇宙環境利用分科会、北海道大学学術交流会館、2013 年 2月 22 日［144］塚野 徹、東野和幸、杉岡正敏（室蘭工業大学）「 白金触媒を用いたメチルシクロヘキサンの熱
分解吸熱反応特性」、第１０回ＨＡＳＴＩ学術講演会、北大、２０１３年２月２２日．
［145］西根賢治、立桶 薫、中田大将、棚次亘弘、東野和幸（室蘭工業大学）「 ロケットスレッドを用いた飛翔体の空力測定」、第１０回ＨＡＳＴＩ学術講演会、北大、２０１３年２月２２日．
［146］近藤光輝、杉岡正敏､東野和幸（室蘭工大），宇宙機推進システムとしての Al/水反応を利用した常圧および高圧水素製造に関する研究、日本航空宇宙学会 第５３回航空原動機・宇宙推進講
演会、倉敷市、2013 年 3 月 4－５日．
［147］山本庸平、湊 亮二郎、杉岡正敏、東野和幸（室蘭工業大学）「 白金アルミナ触媒を使用したバ
イオエタノールの熱分解吸熱反応に関する研究」、日本航空宇宙学会北部支部２０１３年講演会
ならびに第１４回再使用型推進系シンポジウム、東北大、２０１３年３月１４－１５日．
［148］泉俊太郎、杉岡正敏、東野和幸（室蘭工業大学）、東伸幸、小林悌宇「液体ロケットエンジンにおけるバイオエタノールの材料適合性に関する実験的研究」、日本航空宇宙学会北部支部２０１３年
講演会ならびに第１４回再使用型推進系シンポジウム、東北大、２０１３年３月１４－１５日．
［149］髙木正平、上田祐士（室蘭工業大学）「小翼列を用いた翼剥離流のＡctive 制御」、日本航空宇宙学会北部支部２０１３年講演会ならびに第１４回再使用型推進系シンポジウム、東北大、２０１３年３
月１４－１５日．
［150］髙木正平、宮森康臣、大立目浩幸（室蘭工業大学）「低レイノルズ数における二次元物体後流の渦構造の可視化」、日本航空宇宙学会北部支部２０１３年講演会ならびに第１４回再使用型推進系シンポジウム、東北大、２０１３年３月１４－１５日．
［151］上羽 正純、樋口 健、東野 和幸（室蘭工業大学）、「B-3-5 緊急時リアルタイム画像伝送高速小型無人航空機システムの提案」、電子情報通信学会 2013 年総合大会、岐阜大、2013 年 3 月 19
日～22 日［152］竹内僚太朗、上羽 正純（室蘭工業大学）、矢野一人、小林聖（ATR）「B-５-７７ 複数無人航空機
及び追尾アンテナ付地上局間 MIMO 空間相関行列特性の検討」、電子情報通信学会 2013 年
総合大会、岐阜大、2013 年 3 月 19 日～22 日

121

著書
［1］ Minato, R., Higashino, H, Sugioka, M. and Sasayama, Y. “Control of LNG pyrolysis andApplication to Regenerative Cooling Rocket Engine” , Heat Exchanger InTech (2011)

室蘭工業大学紀要[1] 室蘭工業大学紀要第５８号 航空宇宙機システム研究センター特集号(平成 21 年 2 月 ISSN 1344-2708)
1． 小型推進エンジンプロジェクトグループ(1)湊 亮二郎，竹田広人，西村宗真，「反転ファンターボジェットエンジンの研究」(2)東野和幸，杉岡正敏，小林隆夫，湊 亮二郎，丸 祐介，笹山容資，大塚雅也，牧野 隆，坂口
裕之，「LNG ロケットエンジンのメタンコーキング特性に関する基礎研究」(3) 東野和幸，杉岡正敏，小林隆夫，境 昌宏，湊 亮二郎，笹山容資，大塚雅也，沖田耕一，青木賢司，川島秀人，東 伸幸，「LNG ロケットエンジンにおけるサルファアタック・コーキングに関する基礎研究」(4) 東野和幸，駒崎慎一，笹山容資，木元健一，坂口裕之，石橋利幸，「ロケット燃焼室におけるニッケル電鋳特性向上に関する基礎研究」(5) 高橋将人，磯田浩志，棚次亘弘，東野和幸，湊 亮二郎，「熱分解吸熱反応燃料に関する研究」
2．小型無人実験機プロジェクトグループ(6) 溝端一秀，羽田尚太，工藤摩耶，笹山容資，桑田耕明，丸 祐介，湊 亮二郎，棚次亘弘，新井隆景，坪井伸幸，「小型超音速飛行実験機の空力特性」
3．地上・飛行試験設備支援グループ(7) 丸 祐介，宮瀬宗彦，安田有佑，棚次亘弘，溝端一秀，「サブスケール高速走行軌道装置の試
作と走行実験」(8) 湊 亮二郎，溝端一秀，桑田耕明，「吸込み式超音速風洞におけるスターティングロードの計測
と評価」

122

室蘭工業大学航空宇宙機システム研究センター 年次報告書 20132014 年 8 月発行
編集：国立大学法人室蘭工業大学 航空宇宙機システム研究センター
発行：国立大学法人室蘭工業大学 航空宇宙機システム研究センター〒050-8585 室蘭市水元町 27 番 1 号
TEL/FAX: 0143-46-5316
E-mail: higashi@mmm.muroran-it.ac.jp
URL: http://www.muroran-it.ac.jp/aprec/

123