"船体と、該船体の右側及び左側に少なくとも一つずつ備わる少なくとも一対の羽ばたき翼と、該一対の羽ばたき翼をそれぞれ駆動するための少なくとも一対のアクチュエータと、該アクチュエータのコントローラを備えた水中推進体において、前記羽ばたき翼は、前記アクチュエータによりストローク軸まわりに回転させることにより該羽ばたき翼を能動的にアップストローク及びダウンストロークさせるストローク運動を伝達するための前縁はりと、流体力を受ける翼板と、前記前縁はりと前記翼板を連結し、かつ該羽ばたき翼にねじり剛性を与えるバネ部材とを備え、前記ストローク運動によって流体力を翼板面に作用させることにより生じたバネ部材の変形によって、ストローク面に対して常に垂直な軸方向と、翼長方向に対して垂直な翼断面の翼弦方向との間のピッチ角に変位を生じさせて、ピッチ角変位の時間変化であるピッチング運動を受動的に生じさせることから成る水中推進体。前記翼板にねじり剛性を与えるバネ部材は、板バネ或いはコイルバネによって構成し、若しくは、翼板全体を板バネによって構成した請求項1に記載の水中推進体。前記翼板は、所定厚の板状部材によって形成した請求項1に記載の水中推進体。前記コントローラにより前記アクチュエータを制御して、前記ストローク運動のストローク角変位を周期関数に基づき与える請求項1に記載の水中推進体。前記コントローラにより前記アクチュエータを制御して、右側及び左側それぞれの羽ばたき翼のアップストローク運動における代表速さとダウンストローク運動における代表速さを等しく、かつ、右側及び左側それぞれの羽ばたき翼の代表速さを等しくすることで、右側及び左側の羽ばたき翼に同じ推進力を発生させて、船体を直進させ、かつ、右側及び左側それぞれの羽ばたき翼のアップストローク運動における代表速さとダウンストローク運動における代表速さを等しく、かつ、左側の羽ばたき翼の代表速さを右側よりも大きくすることで、左側の羽ばたき翼の推進力を右側よりも大きくして、船体を右へ旋回させ、或いは、その逆にして船体を左へ旋回させる請求項1に記載の水中推進体。右側及び左側それぞれの羽ばたき翼のアップストローク運動における代表速さよりも、ダウンストローク運動における代表速さを大きくして、船体を下降させ、或いは、その逆にして船体を上昇させる請求項5に記載の水中推進体。左側の羽ばたき翼をアップストロークさせると共に、右側の羽ばたき翼はダウンストロークさせるときの両翼の代表速さを、右側の羽ばたき翼をアップストロークさせると共に、左側の羽ばたき翼はダウンストロークさせるときの両翼の代表速さよりも大きくして、船体を反時計回りに回転させ、或いは、その逆にして船体を時計回りに回転させる請求項6に記載の水中推進体。例えば、養殖漁業において、魚病を防ぐための養殖魚の有効なモニタリング方法が無いため、生産効率が低いという問題を解決するために、潜水してモニタリングするロボットが求められている。この場合、養殖魚に刺激を与えるスクリュー推進よりも、生物型推進を用いるロボットが望ましい。この例に見られるように、様々な分野において、生物型推進による水中推進体が求められている。生物型推進のひとつとして、空中や水中の生物が広く利用している羽ばたき推進が挙げられる。従来の羽ばたき翼による水中推進体では、ストローク運動を与えるアクチュエータとは別に、ピッチング運動を能動的に与えるためのアクチュエータを必要としている。ゆえにその機構は複雑で大掛かりなものとなる。本発明は、係る問題点を解決するために、昆虫羽ばたき飛行の力学モデルを水中に適用することで、ピッチング運動を受動的に生じさせ、これによって、ピッチング運動を能動的に与えるためのアクチュエータを不要にして、機構制御系を単純化しつつ、推進効率が高く、多方向への推進が可能な水中推進体を提供することを目的としている。本発明の水中推進体は、船体と、該船体の右側及び左側に少なくとも一つずつ備わる少なくとも一対の羽ばたき翼と、該一対の羽ばたき翼をそれぞれ駆動するための少なくとも一対のアクチュエータと、該アクチュエータのコントローラを備えている。この羽ばたき翼は、アクチュエータによりストローク軸まわりに回転させることにより該羽ばたき翼を能動的にアップストローク及びダウンストロークさせるストローク運動を伝達するための前縁はりと、流体力を受ける翼板と、該前縁はりと該翼板を連結し、かつ該羽ばたき翼にねじり剛性を与えるバネ部材とを備えている。このストローク運動によって流体力を翼板面に作用させることにより生じたバネ部材の変形によって、ストローク面に対して常に垂直な軸方向と、翼長方向に対して垂直な翼断面の翼弦方向との間のピッチ角に変位を生じさせて、ピッチ角変位の時間変化であるピッチング運動を受動的に生じさせる。また、翼板にねじり剛性を与えるバネ部材は、板バネ或いはコイルバネによって構成し、若しくは、翼板面全体を板バネによって構成する。翼板は、所定厚の板状部材によって形成することができる。コントローラによるアクチュエータの制御の際、ストローク運動のストローク角変位を周期関数に基づき与える。コントローラによりアクチュエータを制御して、右側及び左側それぞれの羽ばたき翼のアップストローク運動における代表速さとダウンストローク運動における代表速さを等しく、かつ、右側及び左側それぞれの羽ばたき翼の代表速さを等しくすることで、右側及び左側の羽ばたき翼に同じ推進力を発生させて、船体を直進させる。右側及び左側それぞれの羽ばたき翼のアップストローク運動における代表速さとダウンストローク運動における代表速さを等しく、かつ、左側の羽ばたき翼の代表速さを右側よりも大きくすることで、左側の羽ばたき翼の推進力を右側よりも大きくして、船体を右へ旋回させ、或いは、その逆にして船体を左へ旋回させる。右側及び左側それぞれの羽ばたき翼のアップストローク運動における代表速さよりも、ダウンストローク運動における代表速さを大きくして、船体を下降させ、或いは、その逆にして船体を上昇させる。左側の羽ばたき翼をアップストロークさせると共に、右側の羽ばたき翼はダウンストロークさせるときの両翼の代表速さを、右側の羽ばたき翼をアップストロークさせると共に、左側の羽ばたき翼はダウンストロークさせるときの両翼の代表速さよりも大きくして、船体を反時計回りに回転させ、或いは、その逆にして船体を時計回りに回転させる。本発明ではピッチング運動が受動的に生じるので、ピッチング運動を能動的に与えるためのアクチュエータは全く不要である。翼の駆動制御すべき運動がストローク運動だけとなるので、従来よりも機構制御系を単純化できる。この結果、低コスト化につながる。しかも従来と同等以上の高い推進効率を発揮する。本発明により、例えば、羽ばたき推進型海洋ロボットを構成した際には、魚の成育に好ましくないスクリュー音がなく、かつ魚を刺激しない動きが実現できるので、養殖漁業に最適なものである。本発明に基づき構成した推進型海洋ロボットは、推進効率が高く、外部からのストローク運動だけのコントロールによって右転回、及び左転回させ、さらには、上昇、下降させることが容易に行える。図1は、羽ばたき翼による水中推進体の基本構造を例示する図である。同図に示すように、水中推進体は、潜水可能の船体、船体の右側及び左側にそれぞれ少なくとも1つ備わる少なくとも一対の羽ばたき翼、この羽ばたき翼を駆動するための少なくとも一対のアクチュエータ例えば、電動1軸アクチュエータを船体に内蔵、アクチュエータのコントローラ船体に内蔵からなる。図1には、一対の羽ばたき翼を例示したが、さらに、多く例えば、二対、三対などの羽ばたき翼を左右同数にして備えることができる。羽ばたき翼は、アクチュエータによるストローク運動を伝達するための前縁はり、流体力を受ける翼板、および、翼の柔軟部分であるバネ部材例えば、板バネから構成されている。バネ部材は、前縁はりと翼板の間で両者に結合されて、羽ばたき翼にねじり剛性を与える。前縁はりは水中から船体内部に貫通しており、水の船体内への流入を防ぐため、船体と前縁はりの間をシールする。同図に示すxyz座標系は、空間に対して固定された全体座標系であり、その原点を右側の羽ばたき翼の付け根に取り、x軸を船体の長軸と平行に取っている。またストローク軸はアクチュエータによる羽ばたき翼の回転の軸であり、x軸と平行である。またz軸と前縁はりの軸がなす角を、ストローク角変位φ図3参照という。このように構成した水中推進体は、図10~図16を参照して後述するように、コントローラとアクチュエータにより、羽ばたき翼のストローク運動を制御して、推進力と揺動力を変化させて、船体を推進するだけでなく、潜水可能の船体の直進、左右への旋回、下降、上昇、時計回り、反時計回りの操縦をすることができる。また、船体内には上下を感知するセンサーを搭載して、船体を、ほぼ上下対称に構成することができる。例えば養殖漁業用等として使用される場合に養殖されている魚と接触するなどして、転覆上下逆さまになるした場合、上下を感知するセンサーからの情報をコントローラに送り、上下を切り替えてアクチュエータのコントロールを行い、反転しなかった場合と同じ動きが出来る。その場合、監視カメラの上下の反転は機械的にすることもできるし、電子情報処理によって行い、受信側では上下の反転がない画像が提供できる構成にすることができる。従って、転覆しても、船体を自力で復元する必要がない。図2は、羽ばたき翼のストローク運動を説明する図である。羽ばたき翼のストローク運動を平面的に見た様子を図2に示している。前縁はりをアクチュエータによりストローク軸まわりに、所定の回転角の範囲内で上下に回転させることにより、羽ばたき翼を能動的に羽ばたかせる。以下、この羽ばたき運動をストローク運動と呼ぶ。図2において、空間に固定された全体座標系xyzの原点を右側の羽ばたき翼の付け根に取っている。ストローク運動する前縁はりが描く扇状の面をストローク面という。x軸はこのストローク面に垂直である。前縁はりの軸とz軸のなす角をストローク角変位φという。図2に示すように、ストローク角変位φを、Φ2≦φ≦Φ2Φ:ストローク角の範囲で周期的に変化させる。図1に示すコントローラにより、アクチュエータの制御を行って、ストローク角変位φの周期的変化の関数形、および、その制御パラメータであるストローク周波数fφとストローク角Φを与える。左右の羽ばたき翼の半周期分のストローク運動をアップストロークとダウンストロークといい、次のように定義される;アップストローク:右の羽ばたき翼の先端が点Aを出発し、点Gに至るまでのストローク運動、および、左の羽ばたき翼の先端が点aを出発し、点gに至るまでのストローク運動。ダウンストローク:右の羽ばたき翼の先端が点Gを出発し、点Aに至るまでのストローク運動、および、左の羽ばたき翼の先端が点gを出発し、点aに至るまでのストローク運動。次に、図3~図5を参照して、羽ばたき翼のピッチング運動を説明する。図3は、羽ばたき翼のストローク運動とピッチング運動を説明する図である。羽ばたき翼の運動を立体的に見た様子を図3に示している。同図に示すξηζ座標系は、羽ばたき翼に固定された局所座標系である。ξ軸は常にストローク面に対して垂直全体座標系xyzのx軸と常に平行である。またζ軸は前縁はりの軸と一致しており、前縁はりの軸とz軸のなすストローク角変位φが0のとき、z軸と重なり合う。またη軸は前縁はりのストローク運動の速度ベクトルの方向と平行である。羽ばたき翼のストローク運動によって、水の動圧が翼板面に作用し、板バネに曲げ変形が生じ、ピッチ角変位θが受動的に生じる。ピッチ角変位θは、ストローク面に対して常に垂直なξ軸方向と、翼弦方向の間の角度である。ここで翼弦とは、図1または図3に示すように、翼長方向ζ軸に対して垂直な翼断面において、その前縁と後縁を結んだ線のことである。船体内に内蔵したアクチュエータにより、前縁はりをアクチュエータの軸まわりに回転させて、羽ばたき翼を能動的にストローク運動させる。そのストローク角変位の時間変化を、例えば正弦波状に与える。羽ばたき翼のストローク運動によって、水の動圧が翼板面に作用し、板バネに曲げ変形が生じ、ピッチ角変位θが受動的に生じる。このピッチ角変位θの時間変化をピッチング運動という。本質的には、動圧による前縁はり軸周りのモーメントが作用したとき、羽ばたき翼が前縁はり軸周りにねじれることによって、ピッチ角が受動的に生じるバネ部材であれば、板バネに限定されるものではない。図4は、羽ばたき翼に作用する力を説明する図である。羽ばたき翼は、図1を参照して上述したように、アクチュエータによるストローク運動を伝達するための前縁はり、流体力を受ける翼板、および、羽ばたき翼の柔軟部分である板バネから構成されている。板バネは、その先端において、前縁はりに接合され、その後端において、翼板と接合されている。図4は、羽ばたき翼がη軸正方向にストローク運動ダウンストロークするときの羽ばたき翼の力学的状態を羽ばたき翼の断面において表したものである。即ち、翼板の先端は図2または図3の点Gから点Aに向かって動いているダウンストローク。ストローク運動によって生じる水の動圧が翼板面に垂直に作用し、板バネに曲げ変形が生じる。そして板バネの曲げ変形に伴う弾性反力と水の動圧がつりあうように、ピッチ角変位θストローク面に対して常に垂直なξ軸方向と、翼弦方向の間の角度が受動的に生じる。このピッチ角変位θの時間変化がピッチング運動である。図5は、羽ばたき翼の断面の動きを用いて、羽ばたき翼のピッチング運動を説明する図である。図中の点A~Gは図2または図3の点A~Gと対応している。上図aでは、羽ばたき翼が点Aから点Gに向かって、ストローク運動しアップストローク、下図bでは、羽ばたき翼が点Gから点Aに向かって、ストローク運動しているダウンストローク。ストローク運動は点Aと点Gの間での周期的運動なので、ストローク運動の切り返し点Aまたは点Gにおいて、ストローク角速度dφdtが0となり、ストローク運動の中盤φ0付近において、dφdtが最大となる。ここで動圧の大きさは、羽ばたき翼のストローク運動の速度の2乗に比例し、さらに図4に示すように、動圧はストローク運動に抵抗するように翼板面に作用する。ゆえに図5に示すように、ピッチ角変位θはストローク運動の切り返し点Aまたは点Gにおいて、その大きさがほぼ0となり、ストローク運動の中盤φ0付近において、その大きさが最大となり、しかもストローク運動の方向と逆方向に生じる。このように、本発明ではピッチング運動が受動的に生じるので、ピッチング運動を能動的に与えるためのアクチュエータは全く不要である。図4あるいは図5に示すように、翼弦にピッチ角変位θが生じると、動圧はξ方向x方向とη方向に成分を持つ。図5に示すように、動圧のξ方向成分は常に正である。そこで船体をξ方向x方向に推進するための力推進力として、動圧のξ方向成分を用いる。また図5に示すように、動圧のη方向成分はストローク運動の方向によって正負が変わり、船体をyz面内において揺動する力揺動力となる。船体の重量による慣性が十分でない場合、この揺動力およびそれによるモーメントにより、船体が大きく揺動する場合、船体の姿勢を安定化するための固定翼を取り付けても良い。図6は、翼板形状を例示する図である。上述したように、羽ばたき翼は、アクチュエータによるストローク運動を伝達するための前縁はり、流体力を受ける翼板、および、羽ばたき翼の柔軟部分である板バネから構成されている。前縁はりは、例えば、十分に剛なステンレス角棒で作成することができるが、必ずしも前縁はりが剛体的である必要は無い。但し、顕著な高次の変形モードが伴わないたわみの弾性変形の範囲であり、かつ、アクチュエータによるストローク運動の伝達において、多様な制御の下でも問題が生じない程度の剛性は必要である。翼板は、所定厚及び所定面積を有する板状部材で良いが、その平面形状は、例示のように、小型ガガンボ、大型ガガンボ、ハナアブ、クロバエ等の昆虫の羽形状と幾何学的相似にすることができる詳細は後述する。この翼板は、例えば、バネ部材に比べて、十分剛であるような板厚を有するプラスチック板によって構成する。なお翼断面形状は一般に流線型として表すことができる。そこで前出の図4,図5においては、翼板の断面形状を流線型として例示した。翼の断面形状は、流体の慣性力と粘性力の比であるレイノルズ数に応じて、有利な形状が変わる。航空機や鳥では通常の流線型であるのに対し、昆虫の翼はその極端な場合である薄板である。このことによって、図6で例示した翼板は平板となっている。但し、本発明は、単なる平板に限定されるものでなく、当然、流線型の翼断面形状を用いることもできる。本発明が水中推進体に関するものであることから、例えば、ペンギンなどの羽ばたき推進する水中生物の翼を参考にすることができる。本発明が昆虫飛行を背景に持つこととは関係無く、本発明の水中推進体は、流体の動圧に基づく翼のねじれをピッチング運動の発生受動的ピッチング運動に利用するものであれば、ペンギンなどの水中生物を模擬したものでも良いことは明らかである。なお、以上においては、板バネを用いる場合を想定してきたが、必ずしも板バネでなくてもよい。動圧によって翼板にζ軸図3参照周りのモーメントが作用したとき、翼板がζ軸周りにねじれることによって、ピッチ角が受動的に生じさえすればよいので、翼板にねじり剛性を与える構造としては、板バネに限らない。板バネの替わりとなるバネ部材の構造について、図7を参照して説明する。図7は、板バネのかわりにねじりバネコイルバネを用いた場合を例示する図である。簡単のため、翼板を長方形で表示している。板バネはその剛性を制御しやすいという利点があるものの、運転条件によっては、板バネに大きなたわみが集中することにより材料が破壊する可能性もある。これを回避するために、板バネの翼弦方向長さを大きくし、極端な場合、翼板全体を板バネにする構成や、板バネの替わりとなるバネ部材の構造として、図7に示すように、前縁はりの付け根付近に段付き軸と軸受の構造を設け、段付き軸と軸受の間をコイルバネで接続するなどの構成にすることもできる。羽ばたき翼がストローク軸周りに回転すると、翼板にζ軸周りのモーメントが作用し、コイルバネのねじりによって、翼板にピッチ角が生じる。[船体に作用する平均推進力]上述したように、推進力は動圧をξ方向x方向に分解した成分であるから、動圧の理論式に基づく次式を用いて、平均推進力を推定できる。平均推進力Fξは、FξCξ但し、Cξ:無次元化された推進力、ρf:流体質量密度、AW:翼面積、Vref:翼のストローク運動の代表速さであり、AWとVrefは次式で与えられる。c:平均翼弦長として、AWLWcVref2r2ΦLWfφ但し、LW:翼長長さ、r2:翼面積の2次モーメント、Φ:ストローク角、fφ:ストローク周波数であり、翼長長さLW、平均翼弦長c、と翼のアスペクト比rAの間には次の関係がある。rA2LWc但し、rA:翼のアスペクト比である。図8は、ピッチ角変位の最大θmaxと無次元化された推進力Cξの関係θmaxCξ特性曲線を模式的に示す図である。このθmaxCξ特性曲線は次の3つの性質を有する;θmaxが30~60°の範囲にあるときにCξひいては平均推進力Fξが最大となること、Cξの最大は約2であること、これらの性質が翼型に大きく左右されないこと。後述の[翼のねじり剛性]の項において、θmaxが30~60°となるように翼長方向の軸ζ軸周りの翼のねじり剛性GWを与える。ゆえに上の性質~から、式1のCξを2とすることができる。翼板形状を決めれば、翼のアスペクト比rAと翼面積の2次モーメントr2が決まる。ここでは、図6に例示した小型ガガンボrA107、r206、大型ガガンボrA113、r206、ハナアブrA72、r2054、クロバエrA59、r2052と幾何学的相似にする但し、小型ガガンボと大型ガガンボではrAとr2の値が近く、形状がほとんど同じある。但し、θmaxCξ特性曲線の性質から、本発明はこれらの翼型に限定されるものではない。図1のコントローラにより、アクチュエータを制御して、ストローク角変位φを周期関数に基づき与える。平均推進力Fξ式1参照は、周期関数の特性量のうち最も基本的なストローク角Φとストローク周波数fφでほぼ決まり、その具体的関数形への依存性が小さい。実際、例えば、ストローク角速度dφdtの時刻歴を台形波で与える場合、極端な場合である矩形波や三角波を含む様々な台形波に対しても、ストローク角Φとストローク周波数fφが同じであれば、平均推進力Fξはほとんど変わらないことを実験的に確認している。従って、本発明は、特定の関数形に限定されることなく、周期関数一般に適用される。上述の翼板形状に加えて、翼長LW、ストローク角Φ、ストローク周波数fφを決めれば、平均推進力Fξを推定するために必要な翼形状、翼寸法、および、ストローク運動の全てのパラメータが決まる。例えば、翼長を全ての翼型でLW0225mとし、ストローク角Φとストローク周波数fφをそれぞれ、小型ガガンボの翼型の場合:123°と0521Hz、大型ガガンボの翼型の場合:120°と0655Hz、ハナアブの翼型の場合:108°と0809Hz、クロバエの翼型の場合:150°と0391Hzと設定する。水の質量密度ρf1000kgm3とし、式1を用いて、平均推進力Fξを推定すると、小型ガガンボの翼型の場合:約09N実験値は約1N、大型ガガンボの翼型の場合:約12N実験値は約14N、ハナアブの翼型の場合:約19N実験値は約19N、クロバエの翼型の場合:約10N実験値は約12Nとなる。ここで括弧内の実験値は、水槽内での羽ばたき翼の性能試験により得たものである。但し、翼のねじり剛性を後述の方法により決めた。推定値が実験値と良く一致することから、式1は十分実用的である。この検証された式1を利用することで、本水中推進体を用途に応じて設計できる。すなわち用途で要求される目標推力を定め、それを達成できる翼型とストローク運動のパラメータを式1に基づき決める。すなわち本発明は、上で例として取り上げた具体的な翼形状、翼寸法、ストローク運動のパラメータに限定されるものではない。[翼のねじり剛性]図9は、コーシー数Ch流体動圧と弾性力の比を表す無次元数とピッチ角変位の最大θmaxとの関係を表す特性曲線Chθmax特性曲線の典型例を模式的に示す図である。Chθmax特性曲線は次の3つの性質を持つ;図9に示すように、特性曲線の左端において、x軸に漸近しながら原点に収束すること、図9に示すように、特性曲線の右端において、θ90°の水平線に漸近すること、これらの性質が翼型に依存しないこと。コーシー数Chは次式により与えられる。ChρfVrefc4fφGWここでGWは以下の式によって定義される翼長方向の軸ζ軸まわりの翼のねじり剛性である。MθGWθここでMθは翼に作用する翼長方向まわりのモーメントである。板バネを用いるとすると、このGWは板バネの曲げ剛性と一致するので、GWは次式によって与えられる。GWEpsIpscpsここでEps、Ips、cpsはそれぞれ板バネのヤング率、断面2次モーメント、板バネの翼弦方向長さである。Ipsは板バネの断面形状が長方形であることを仮定して次式で与えられる。Ipslpstps312ここでlps、tpsはそれぞれ板バネの翼長方向長さ、板厚である。以上に基づき、GWを次のように決める;まず図9に基づき、θmaxが30~60°となるChの値を求める。上述したように、このとき平均推進力Fξが最大となる。図9はあらかじめ実験的に求めておく。図9が利用できない場合は、次のようにChの値を決める。θmaxが30~60°となるChは、翼形状に依らず、おおよそCh001~01であることが、本発明者の実験により、明らかになっている。そこでこの範囲で当たりをつけて、Chの値を決める。次に上記のように決めたChの値とChの式5により、GWを求める。ここで式5の右辺の分子に含まれる全てのパラメータは、上述のようにして既に与えられていることに注意する。板バネの場合、上記のように求めたGWを式7に代入し、この式と式8を用いて、板バネの具体的な寸法を決めることができる。[船体の操縦]本推進体においては、翼のストローク運動を能動的に与えるだけで、船体を推進できる。さらに、それだけでなく、翼のストローク運動を能動的に与えるだけで、船体を操縦することも可能である。ゆえに本推進体は、通常の船舶と異なり、船体の操縦装置を別に備える必要が無いという利点を有する。しかも本推進体における船体の操縦は、全方位に方向転換できるものである。本推進体は、この全方位推進を、翼のストローク運動を能動的に与えるための2つのアクチュエータで実現することができる。全方位推進を通常の船舶で実現するのは容易でなく、多数のアクチュエータを含む複雑な操縦装置が必要である。平均推進力Fξの推定式は上述した式1に示す通りである。ここでは平均揺動力の推定式を示す。揺動力は推進力と同様に動圧に基づくので、平均推進力Fξと同様の次式により平均揺動力Fηを推定できる。平均揺動力Fηは、FηCη但し、Cηは無次元化された揺動力であり、その大きさはCξと同様に約2である。またρf:流体質量密度、AW:翼面積、Vref:翼のストローク運動の代表速さであり、AWとVrefはそれぞれ式2、式3の通りである。また揺動力の向きについては次の通りである。図2において、右の翼の先端が点Aから点Gにアップストロークの運動をするとき、図5aに示すように、揺動力の向きはη軸正であり、点Gから点Aにダウンストロークの運動をするとき、図5bに示すように、揺動力の向きはη軸負である。図2の左の翼についても、揺動力の発生の仕方は同じである。以下では、図1のコントローラとアクチュエータにより、翼のストローク運動を制御して、推進力と揺動力を変化させて、船体を操縦する例を示す。説明に用いる図10~図16において、左右の翼は、形状、大きさ、材質が同一であるとする。また図10~図16に示すxyz座標系において、x軸を船体の長軸と平行であるように取り、z軸を、翼のストローク角変位φが0のときに、前縁はりの軸と平行であるように取る。図10は、船体の直進船体の長軸正の方向に進むことの例を示す図である。は推進体をy方向から見た様子を示し、は推進体をx方向から見た様子を示している。右の翼において、点Aを出発し、点Gに至るまでのアップストロークの運動における代表速さと、点Gを出発し、点Aに至るまでのダウンストロークの運動における代表速さが等しいとする。ここでアップストロークの運動における代表速さとは、アップストロークの運動に関して、式1を適用して得られるものであり、アップストロークの運動の速さの平均に相当する量である。ダウンストロークの運動における代表速さも同様に定義する。このときアップストロークにおける平均揺動力とダウンストロークにおける平均揺動力は、式9より大きさが等しく、向きが逆となる。ゆえに右の翼のストローク運動全体アップストロークダウンストロークにおける平均揺動力は、ゼロとなる。左の翼においても、点aを出発し、点gに至るまでのアップストロークの運動における代表速さと、点gを出発し、点aに至るまでのダウンストロークの運動における代表速さが等しいとすると、ストローク運動全体アップストロークダウンストロークにおける平均揺動力はゼロとなる。また左右の翼の代表速さを同じにすることで、左右の翼に同じ推進力が発生し、それにより、船体が直進する。図11は、船体の左への旋回同図において、z軸負の方向に進むことの例を示す図である。は推進体をy方向から見た様子を示し、は推進体をx方向から見た様子を示している。右の翼において、アップストロークにおける代表速さとダウンストロークにおける代表速さが等しいとすると、上述した直進の例と同様に、ストローク運動全体における平均揺動力はゼロとなる。左の翼においても、アップストロークにおける代表速さとダウンストロークにおける代表速さが等しいとすると、ストローク運動全体における平均揺動力はゼロとなる。ここで上述の直進の場合と異なり、左の翼の代表速さを右の翼の代表速さよりも小さくすることにより、式1から、左の翼の平均推進力が右の翼の平均推進力よりも小さくなる。従って、図11aに示すように、船体が左へ旋回する。図12は、船体の右への旋回同図において、z軸正の方向に進むことの例を示す図である。は推進体をy方向から見た様子を示し、は推進体をx方向から見た様子を示している。上述の左への方向転換とは逆に、左の翼の代表速さを右の翼の代表速さよりも大きくすることにより、式1から、左の翼の平均推進力が右の翼の平均推進力よりも大きくなる。それ以外は左への方向転換と同じである。従って、図12に示すように、船体が右へ旋回する。図13は、船体の下降同図において、y軸負の方向に進むことの例を示す図である。は両翼ともアップストロークを示し、は両翼ともダウンストロークを示している。同図aの右に示すように、左の翼がアップストロークするとき、右の翼もアップストロークし、また同図bの右に示すように、左の翼がダウンストロークするとき、右の翼もダウンストロークする。このとき同図aにおける両翼の代表速さよりも、同図bにおける両翼の代表速さのほうが大きいので、同図aにおいて生じる平均揺動力よりも、同図bにおいて生じる平均揺動力のほうが大きい。ゆえにストローク全体としての平均揺動力は、y軸の負の方向に働く。従って、図13において、潜水可能の船体が下降する。図14は、船体の上昇同図において、y軸正の方向に進むことの例を示す図である。は両翼ともアップストロークを示し、は両翼ともダウンストロークを示している。上述した船体の下降の場合と同様に、同図aの右に示すように、左の翼がアップストロークするとき、右の翼もアップストロークし、また同図bの右に示すように、左の翼がダウンストロークするとき、右の翼もダウンストロークする。一方、船体の下降の場合とは逆に、同図bにおける両翼の代表速さよりも、同図aにおける両翼の代表速さのほうが大きいので、同図bにおいて生じる平均揺動力よりも、同図aにおいて生じる平均揺動力のほうが大きい。ゆえにストローク全体としての平均揺動力は、y軸の正の方向に働く。従って、図14において、船体が上昇する。図15は、船体の反時計回りの回転同図において、船体の長軸に関して反時計回りに回転することの例を示す図である。は左の翼がアップストローク,右の翼がダウンストロークを示し、は左の翼がダウンストローク,右の翼がアップストロークを示している。同図aの右に示すように、左の翼がアップストロークするとき、右の翼はダウンストロークし、同図bの右に示すように、左の翼がダウンストロークするとき、右の翼はアップストロークする。このとき同図aにおける両翼の代表速さのほうが、同図bにおける両翼の代表速さよりも大きいので、同図aにおいて生じる平均揺動力による回転モーメントのほうが、同図bにおいて生じる平均揺動力による回転モーメントよりも大きい。ゆえにストローク全体としての平均揺動力による回転モーメントは、船体の長軸に関して反時計回りに働く。従って、図15において、船体は反時計回りに回転する。図16は、船体の時計回りの回転同図において、船体の長軸に関して時計回りに回転することの例を示す図である。は左の翼がアップストローク,右の翼がダウンストロークを示し、は左の翼がダウンストローク,右の翼がアップストロークを示している。上述した反時計回りの回転の場合と同様に、同図aの右に示すように、左の翼がアップストロークするとき、右の翼はダウンストロークし、同図bの右に示すように、左の翼がダウンストロークするとき、右の翼はアップストロークする。一方、反時計回りの回転とは逆に、同図bにおける両翼の代表速さのほうが、同図aにおける両翼の代表速さよりも大きいので、同図bにおいて生じる平均揺動力による回転モーメントのほうが、同図aにおいて生じる平均揺動力による回転モーメントよりも大きい。ゆえにストローク全体としての平均揺動力による回転モーメントは、船体の長軸に関して時計回りに働く。従って、図16において、船体は時計回りに回転する。"