(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】6856930
(24)【登録日】2021年3月23日
(45)【発行日】2021年4月14日
(54)【発明の名称】ロータ、ドローン及びヘリコプタ
(51)【国際特許分類】
B64C 27/467 20060101AFI20210405BHJP
B64C 11/18 20060101ALN20210405BHJP
【FI】
B64C27/467
!B64C11/18
【請求項の数】4
【全頁数】11
(21)【出願番号】特願2017-25535(P2017-25535)
(22)【出願日】2017年2月15日
(65)【公開番号】特開2018-131038(P2018-131038A)
(43)【公開日】2018年8月23日
【審査請求日】2019年12月19日
(73)【特許権者】
【識別番号】503361400
【氏名又は名称】国立研究開発法人宇宙航空研究開発機構
(74)【代理人】
【識別番号】100104215
【弁理士】
【氏名又は名称】大森 純一
(74)【代理人】
【識別番号】100196575
【弁理士】
【氏名又は名称】高橋 満
(74)【代理人】
【識別番号】100168181
【弁理士】
【氏名又は名称】中村 哲平
(74)【代理人】
【識別番号】100117330
【弁理士】
【氏名又は名称】折居 章
(74)【代理人】
【識別番号】100160989
【弁理士】
【氏名又は名称】関根 正好
(74)【代理人】
【識別番号】100168745
【弁理士】
【氏名又は名称】金子 彩子
(74)【代理人】
【識別番号】100176131
【弁理士】
【氏名又は名称】金山 慎太郎
(74)【代理人】
【識別番号】100197398
【弁理士】
【氏名又は名称】千葉 絢子
(74)【代理人】
【識別番号】100197619
【弁理士】
【氏名又は名称】白鹿 智久
(72)【発明者】
【氏名】原田 正志
【審査官】
志水 裕司
(56)【参考文献】
【文献】
特表2017−512695(JP,A)
【文献】
特表2016−524075(JP,A)
【文献】
特許第6778440(JP,B2)
【文献】
赤坂剛史,外2名,小型無人飛行体に搭載するプロペラの高効率化について〜長時間飛行の実現に向けて〜,川田技報,川田工業株式会社,2010年,Vol. 29,p. 1-6,URL,https://www.kawada.co.jp/technology/gihou/vol_29.html
(58)【調査した分野】(Int.Cl.,DB名)
B64C 27/467
B64C 11/18
(57)【特許請求の範囲】
【請求項1】
ソリディティが10%以下となるロータであり、
2枚のブレードを有し、
前記各ブレードの形状が、
半径10%から35%の領域に最大翼弦長をもち、
前記最大翼弦長の位置から先端方向にかけて翼弦長が減少し、
前記翼弦長の減少が所定の位置で最大になり、その位置から先端方向にかけて前記翼弦長の減少がゆるやかになり、
半径70%から95%の領域に変曲点があり、
先端が楕円形であり、
翼弦に2点で接する接線が翼端で作る仮想の翼端の翼弦長に対し、最大翼弦長は3倍以上の長さであり、
半径中程のくびれ量は仮想の翼端の30%以下であり、
取付け角は前記翼弦長が最大になる付近で最大となり、前記翼弦長が最大に位置から先端にかけて減少し、
付け根の前記取付け角は当該取り付け角が最大値を取る位置から付け根に行くほど減少している
ロータ。
【請求項2】
請求項1に記載のロータであって、
前記各ブレードの形状が、
翼弦長は付根から先端にかけて1階微分が連続な曲線である
ロータ。
【請求項3】
請求項1又は2に記載のロータを有するドローン。
【請求項4】
請求項1又は2に記載のロータを有するヘリコプタ。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、ドローンやヘリコプタなどに用いられるロータ、そのようなロータを有するドローン及びヘリコプタに関する。
【背景技術】
【0002】
近年、4発以上のロータを持つドローンと呼ばれる飛行体が撮影、観測目的で盛んに利用されている。ドローンは通常プラスチック製の2枚羽の固定ピッチロータを使用しており、このロータの性能によりペイロードの重さ及び飛行時間が左右される。ロータ形状は様々であり、単純テーパから古典プロペラ理論(非特許文献1)を応用したと思われる、付け根の翼弦長が広く、そこから先端に向けて急激に翼弦長を減らし、ある小さなテーパ比に移行し、先端を丸めた形状のものなどが見られる。また有人小型ヘリコプタとしてGEN CORPORATIONが開発したGEN H?4(非特許文献2)は同様に固定ピッチロータを使用しており、ロータに単純テーパを採用している。
【0003】
ヘリコプタの分野では英国のWestrand LinxがBERP(非特許文献2,特許文献1)と呼ばれる先端の翼弦長を拡大し、後退角を持たせたロータを使用し、世界速度記録を達成している。本ロータは先端の翼弦を拡大した不連続な前縁を飛行機で使用されるドッグツースとして機能させ、縦渦を発生させて上面の剥離を抑えることで後退側での揚力係数を稼ぐとともに、先端の後退角によって前進側での造波抵抗を減少させる事を目的としている。ヘリコプタは単純運動量理論とロータ回転面積と入力パワから求められる理論最大揚力で実際の揚力を除したフィギュアオブメリット(非特許文献3)と最大速度をトレードオフしなければならないため、ロータ形状は矩形か先端に弱いテーパがついた形状のものが多く使用されている。これに対してドローンはホバリングする時間が長いためフィギュアオブメリットの高いロータを使用し、飛行時間及びペイロードを増加させようとする要求が高い。
【先行技術文献】
【特許文献】
【0004】
【特許文献1】米国特許第5174721号明細書
【非特許文献】
【0005】
【非特許文献1】Adkins, C. N., Liebeck, R. H., "Design of Optimum Propellers," AIAA−83−0190, Jan. 1983
【非特許文献2】柳沢 源内、二重反転一人乗りヘリコフターGEN H?4 開発と将来展望、Heli Japan 2002、 11月 2002
【非特許文献3】https://en.wikipedia.org/wiki/BERP_rotor
【非特許文献4】R. W. Prouty,"Helicopter Performance, Stability, and Control,"1990
【非特許文献5】原田正志、"低レイノルズ数プロペラの設計法、"JAXA-RR-06?032、2007
【発明の概要】
【発明が解決しようとする課題】
【0006】
フィギュアオブメリットの大きなロータを設計する際に有効な手法としてAdkinsらの方法(非特許文献4)がある。しかしAdkinsらの方法は実際のロータの後流に見られる縮流(筒状の流れの半径の縮小)や、渦の巻き込み(ロールアップ)を反映出来なかった。そのためフィギアオブメリットは最大とはなっていなかった。またAdkinsらの方法では付根付近で翼弦長が大きくなりすぎるため、ブレード付根付近の翼弦長を一定値以下に修正する処理が必要となるが、この際に取付け角が修正されていないため、翼弦長を制限した際の最適な取付け角となっていなかった。ドローンに用いられるロータは一般に小型であるため、翼型のレイノルズ数が低い。この低いレイノルズ数で大きな揚力係数を発生させるため、薄く大きなキャンバーを持った翼型が使用される。大きなキャンバーを持つためモーメント係数Cmが負の大きな値をとる。このためブレードを捩じりさげるモーメントが働くが、Adkinsらの方法で得られるブレードはブレード半径の中程がくびれている上、前述のように薄い翼型を用いているため、捻れに耐える剛性が不足し、捩じれる問題がある。
【0007】
以上のような事情に鑑み、本発明の目的は、高いフィギュアオブメリットを得つつ、高い捻れ剛性を確保出来るロータ、そのようなロータを有するドローン及びヘリコプタを提供する事にある。
【課題を解決するための手段】
【0008】
上記目的を達成するため、本発明の一形態に係るロータは、
ソリディティが10%以下となるロータであり、
2枚のブレードを有し、
前記各ブレードの形状が、
半径10%から35%の領域に最大翼弦長をもち、
前記最大翼弦長の位置から先端方向にかけて翼弦長が減少し、
前記翼弦長の減少が所定の位置で最大になり、その位置から先端方向にかけて前記翼弦長の減少がゆるやかになり、
半径70%から95%の領域に変曲点があり、
先端が楕円形であり、
翼弦に2点で接する接線が翼端で作る仮想の翼端の翼弦長に対し、最大翼弦長は3倍以上の長さであり、
半径中程のくびれ量は仮想の翼端の30%以下であり、
取付け角は前記翼弦長が最大になる付近で最大となり、前記翼弦長が最大に位置から先端にかけて減少し、
付け根の前記取付け角は当該取り付け角が最大値を取る位置から付け根に行くほど減少している。
【0009】
前記各ブレードの形状は、翼弦長は付根から先端にかけて1階微分が連続な曲線であることが好ましい形態である。
【0010】
本発明の別形態に係るドローンは、上記構成のロータを有する。
【0011】
本発明のさらに別形態に係るヘリコプタは、上記構成のロータを有する。
【発明の効果】
【0012】
本発明により、ロータのフィギュアオブメリットが高く、ドローン、ヘリコプタに使用すれば、飛行時間、ペイロード重量を増す事が出来る。またAdkinsらの方法ではブレード付け根付近の翼弦長が大きくなり過ぎ、これを一定の値以下に制限すると取付け角が最適ではなくなっていた。本発明によれば翼弦長を制限しつつも取付け角が最適化される。Adkinsらの方法では半径の中程の翼弦長にくびれが現れるが、本発明で適度に翼弦長を制限する事で半径中程の翼弦長を太らす事が出来、またテーパ比を大きくし、捻れに強い平面形状とする事が出来る。
【図面の簡単な説明】
【0013】
【
図1】本発明の一実施形態に係るロータの上面図である。
【
図2】ロータの翼端渦の軌跡を示している。a)はAdkinsらの方法によるロータの翼端渦の軌跡であり、b)は本発明者である原田の渦法によるロータの翼端渦の軌跡である。
【
図3】a)は実際の翼端渦のロールアップであり、b)はモデル化した翼端渦のロールアップである。
【
図4】従来方法の最適な翼弦長分布を示すグラフである。
【
図5】従来方法の最適な取付角を示すグラフである。
【
図6】翼弦長を制限した従来方法の最適な翼弦長分布を示すグラフである。
【
図7】翼弦長を制限した従来方法の最適な取付角を示すグラフである。
【
図8】本発明の翼弦長分布の特徴を示すグラフである。
【
図10】原田の渦法を説明するための図であり、プロペラの座標系と放出渦を示している。
【
図11】原田の渦法を説明するための図であり、
図10における一枚目のブレードの拡大図を示している。
【
図12】原田の渦法を説明するための図であり、ブレードの断面と流入速度を示している。
【
図13】原田の渦法を説明するための図であり、ブレード翼素に働く力を示している。
【発明を実施するための形態】
【0014】
以下、図面を参照しながら、本発明の実施形態を説明する。
【0015】
図1は本発明の一実施形態に係るロータの上面図である。
【0016】
図1に示すように、ロータ6は、2枚のブレード3と、シャフト穴5を有するハブ4とを有する。
【0017】
ここで、Adkinsらの方法では、
図2のa)に示す様に、一定ピッチかつ一定半径の後流しか扱えない。これに対して、本発明者である原田の渦法(非特許文献5参照。後述する。)を用いることで、さらに後流の形状を次の様に定めることが出来る。
1)
図2のb)に示すように翼端渦半径rが次式
r/R=A+(1−A)e
−λψ (1)
に従い収縮する。ここでRはロータ半径、Aは実験から得られる定数であり0.76から0.8、好ましくは0.78である。またλも実験により得られる定数であり、ψはブレード方位角である。
2)
図3のa)に示す様に翼端から螺旋渦のピッチ程度の範囲にある放出渦は、直ちにロールアップして翼端渦1となり一群となって移動する。
3)翼端渦1は翼端渦1より内部の渦の移動速度のおよそ半分の速度で移動する。
4)ロータ回転面で内部の渦は単純運動量理論から導かれる流入速度vで下方に移動する。上記の式(1)に従い放出渦の半径が収縮するにつれ速度が速くなり、次式で与えられるv'の速度で移動する。
【0018】
v'=v(R/r)
2 (2)
5)翼端渦1の模擬は困難であるため、
図3のb)で示す様に翼端渦として一群になる渦を半径方向に移動させず等速度で移動させるモデル2を用いる。
【0019】
このような後流渦を用い、Adkinsらの方法では大きくなりすぎる付根付近の翼弦長を、原田の渦法から得られる評価関数を翼弦長の制約条件の下で最小化する事で、適度な翼弦長を持つロータ6を得る事が出来る。
【0020】
この実施形態では、以下の形状のロータ6を得る。
1.ソリディティが0.1以下であるロータ6でかつ
2.ブレード3の枚数が2枚であり、
3.翼弦長は付根から先端にかけて1階微分が連続な曲線であり、
4.半径10%から35%に最大翼弦長をもち、
5.そこから先端方向にかけて翼弦長が減少し
6.減少が一旦最大になった後に、減少がゆるやかになり、半径70%から95%に変曲点があり、
7.先端が楕円形であり、
8.翼弦に2点で接する接線が翼端で作る仮想の翼端の翼弦長に対し、最大翼弦長は3倍以上の長さであり、
9.半径中程のくびれ量は仮想の翼端の30%以下であり、
10.取付け角は翼弦長が最大になる付近で最大となり、そこから先端にかけて減少し
11.付け根の取付け角は最大値を取る位置から付け根に行くほど減少している。
【0021】
この形状のロータ6を使用する事で高いフィギュアオブメリットを得つつ、高い捻れ剛性を確保出来る。
【0022】
図4に原田の渦法を
図2のa)の渦モデルを用いて最適化した際の翼弦長分布を、また
図5に取付け角分布を示す。ここで計算条件を表1に示す。
【0024】
本形状は
図4に示すように付け根の翼弦長が12cmに達しているため、実際的ではない。市販のドローンの中にはこのブレード3の付け根の翼弦長を単に実用的な長さに切り落として使用している例がある。
【0025】
原田の渦法は評価関数を最小化する事で解を得るので、翼弦長が4cmを超えると値が大きくなるペナルティ関数を評価関数に加えることで翼弦長の制限を課す事が出来る。原田の渦法において
図2のa)の渦モデルを用い、さらに翼弦長を制限して最適化した際の翼弦長分布を
図6に、また取付け角分布を
図7に示す。このときのフィギュアオブメリットは89.4%である。
【0026】
得られた翼弦分布に二点で接する線分ABの翼端Cにおける値、仮想の翼端翼弦長Dは1.36cm程度である。一方、半径Eにおいて翼弦長Fは最大値4.79cmをとり、仮想の翼端翼弦長Dの3.5倍である。また翼端付近で翼弦長はなだらかに減少するが、翼端Cにおいて切り落とした形状になっている。Mで示した半径中程のくびれは0.32cmであり、仮想の翼端の翼弦長の23%である。
【0027】
また得られた取付け角は翼弦長が最大値を取る半径Eに近い半径Iにおいて最大値9.32度を取る。付け根付近の取付け角分布への接線KLは正の勾配を持つ。
【0028】
原田の渦法において
図3のb)のモデルを用い、さらに翼弦長を制限して最適化した際の翼弦長分布を
図8に、また取付け角分布を
図9に示す。このときのフィギュアオブメリットは89.5%である。
【0029】
得られた翼弦分布に二点で接する線分ABが翼端Cにおける値、仮想の翼端翼弦長Dは1.25cm程度である。一方、半径Eにおいて翼弦長Fは最大値4.94cmをとり、仮想の翼端翼弦長Dの4.0倍である。また翼端付近で翼弦長はなだらかに減少するが、翼端C近傍で急に減少し、翼端が楕円形状になっている。Mで示した半径の中程で現れるくびれ量は、0.14cmと仮想の翼端の翼弦長の11%に抑えられる。
【0030】
また得られた取付け角は翼弦長が最大値を取る半径Eに近い半径Iにおいて最大値9.0度を取る。付け根付近の取付け角分布への接線KLは正の勾配を持つ。
【0031】
図9で示した翼弦長を用いて設計したロータ6が
図1に示したものである。
【0032】
本発明に係るロータ6は、ドローンのロータブレード、固定ピッチ同軸二反転ロータ式有人ヘリコプタのロータブレード、固定ピッチ同軸二重反転ロータ式無人ヘリコプタのロータブレード等に用いる事が出来る。
【0033】
本発明は、上記の実施形態には限定されず、その技術思想の範囲で様々な変形が可能であり、その変形の範囲も本発明の技術的範囲に属するものである。
【0034】
次に、上述した原田の渦法を説明する。
【0035】
図10にプロペラの座標系と放出渦を示す。プロペラは回転しながらx軸方向に移動し、移動した軌跡に放出渦が残されると考える。
【0036】
図11に一枚目のブレードの拡大図を示す。回転軸からr
i離れたブレード上にi番目の代表点があり、j番目の放出渦は図中白丸で示された様に短い線分に離散化する。j番目の単位強度の放出渦がi番目の代表点に引き起こすx方向の誘導速度をX
ij、z方向の誘導速度をZ
ijとすると、i番目の代表点に引き起こされるx方向の誘導速度u
i、z方向の誘導速度w
iはそれぞれ
u
i=ΣX
ijг
j (3)
w
i=ΣZ
ijг
j (4)
で与えられる。ここで、X
ij、Z
ijはビオサバールの法則から得られる定数であり、г
jはj番目の放出渦である。
【0037】
図12にブレードの断面と流入速度を示す。ブレードに流入する空気の接線方向成分をU
TとするとU
Tは
U
T=r
iΩ−w
i (5)
で与えられる。ここでΩはプロペラの回転角速度である。またブレードに流入する空気の軸方向成分をU
PとするとU
Pは
U
P=U−u
i (6)
で与えられる。ここでUはプロペラの前進速度である。流入角φ
i及び流入速度V
iはそれぞれ次式で与えられる。
【0038】
φ
i=tan
−1(U
P/U
T) (7)
V
i=√(U
P2+U
T2) (8)
図13にブレード翼素に働く力を示す。i番目の翼素に働く局所揚力dL
iはクッタ=ジューコフスキーの定理により
dL
i=ρV
iг
idb (9)
で与えられる。ここでρは空気密度、dbは翼素の幅である。または、揚力係数C
Lを用いて
dL
i=1/2・ρV
i2C
Lc
idb (10)
で表される。ここで、c
iはi番目の翼素の翼弦長である。(9)、(10)式よりc
iは次式で与えられる。
【0039】
c
i=2г
i/C
LV
i (11)
i番目の翼素に働く局所抵抗dD
iは抵抗係数C
Dを用いて次式で与えられる。
【0040】
dD
i=1/2・ρV
i2C
Dc
idb (12)
C
Dはレイノルズ数とC
Lの関数であるが、C
Lに定数を用い、C
Dがレイノルズ数に対しては鈍感であるとして定数としてよい。
【0041】
局所揚力dL
iと局所抵抗dD
iの合力の軸方向分力は局所推力dT
iとなり、接線方向分力をdN
iとするとそれぞれ、
dT
i=dL
icosφ
i−dD
isinφ
i (13)
dN
i=dL
isinφ
i+dD
icosφ
i (14)
となる。局所吸収パワdP
iはdN
ir
iΩであるから次式で与えられる。
【0042】
dP
i=(dL
isinφ
i+dD
icosφ
i)r
iΩ (15)
結局、推力と吸収パワは(13)式と(15)式よりそれぞれ
T=BΣdT
i (16)
P=BΣdP
i (17)
で与えられる。ここでBはブレード枚数である。
【0043】
プロペラの最適設計問題はC
L、C
D、Ω、設計パワP
0を設定し、гiを未知数として
【0045】
[目的関数]:−T
を最小化する
最適化問題に帰着する。
【符号の説明】
【0046】
1 翼端渦
2 翼端渦のモデル
3 ブレード
4 ハブ
5 シャフト穴
6 ロータ