JPS59220496A - 船舶の自動操舵装置 - Google Patents
船舶の自動操舵装置Info
- Publication number
- JPS59220496A JPS59220496A JP58092088A JP9208883A JPS59220496A JP S59220496 A JPS59220496 A JP S59220496A JP 58092088 A JP58092088 A JP 58092088A JP 9208883 A JP9208883 A JP 9208883A JP S59220496 A JPS59220496 A JP S59220496A
- Authority
- JP
- Japan
- Prior art keywords
- angle
- gain
- disturbance
- controller
- ship
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004088 simulation Methods 0.000 claims abstract description 17
- 238000011156 evaluation Methods 0.000 claims description 32
- 238000004364 calculation method Methods 0.000 claims description 26
- 239000000446 fuel Substances 0.000 abstract description 16
- 230000003595 spectral effect Effects 0.000 description 14
- 238000000034 method Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 3
- 239000003795 chemical substances by application Substances 0.000 description 2
- 230000009194 climbing Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 210000003127 knee Anatomy 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000002747 voluntary effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
(イ)産業上の利用分野
本発明は船舶の自動操舵装置に関し、詳しくは、船舶が
波・風などの外乱下で航行する際、その推進エネルギの
直進航行時の推進エネルギからの増加率を評価関数とし
、この評価関数値が最小となるように操舵させる省エネ
ルギ型の船舶の自動操舵装置に関する。
波・風などの外乱下で航行する際、その推進エネルギの
直進航行時の推進エネルギからの増加率を評価関数とし
、この評価関数値が最小となるように操舵させる省エネ
ルギ型の船舶の自動操舵装置に関する。
(ロ)従来技術
従来の船舶の自動操舵装置は、周知の通り、設定された
方位に船舶を直進航行させること、すなわち、方位偏差
をできるだ′け小さくすることにある。このために、ジ
ャイロコンパスなどにより検出された実方位と設定方位
との差に、例えば、比例、微分、積分などの演算を施し
て舵を取らせるフィードバック制御が行なわれている。
方位に船舶を直進航行させること、すなわち、方位偏差
をできるだ′け小さくすることにある。このために、ジ
ャイロコンパスなどにより検出された実方位と設定方位
との差に、例えば、比例、微分、積分などの演算を施し
て舵を取らせるフィードバック制御が行なわれている。
さらに最近においては、設定された方位に航行し、しか
も、省エネルギを目指すように制御している自動操舵装
置が出現している。そのような装置で制御を行なう際、
ある評価関数を最小にするようそのゲインを調整してい
る。そして、現在までに提案されなお、δ2は舵角の自
乗平均値、 ψ2は方位角の自乗平均値、 rは方位角速度の自乗平均値、 λ1、λ2、λ3は重み係数、 であって、最適レギュレータ理論などを用いて最適制御
則を計算によって求めている。
も、省エネルギを目指すように制御している自動操舵装
置が出現している。そのような装置で制御を行なう際、
ある評価関数を最小にするようそのゲインを調整してい
る。そして、現在までに提案されなお、δ2は舵角の自
乗平均値、 ψ2は方位角の自乗平均値、 rは方位角速度の自乗平均値、 λ1、λ2、λ3は重み係数、 であって、最適レギュレータ理論などを用いて最適制御
則を計算によって求めている。
このような制御において、航行中の船体の運動状態を諸
検出量より同定するものや、その状態によって重み係数
を可変にする方法などが採られているが、いずれにして
もジャイロコンパスで検出された方位角に基づいて上記
評価関数値を最小にするものである。このような評価関
数月は上述したような自乗項のみで表わされており、後
述する斜杭角を加味する場合にも斜杭角を舵角と方位角
速度とを用いて近似的に求めている。このような評価関
数が用いられているのは、それが自乗項のみからなって
いてその値を最小にする制御が、旧来から存在する近代
制御理論を用いて極めて容易に行ない得ること、および
、斜杭角βを高精度で検出するのが容易でないという理
由に基づいている。
検出量より同定するものや、その状態によって重み係数
を可変にする方法などが採られているが、いずれにして
もジャイロコンパスで検出された方位角に基づいて上記
評価関数値を最小にするものである。このような評価関
数月は上述したような自乗項のみで表わされており、後
述する斜杭角を加味する場合にも斜杭角を舵角と方位角
速度とを用いて近似的に求めている。このような評価関
数が用いられているのは、それが自乗項のみからなって
いてその値を最小にする制御が、旧来から存在する近代
制御理論を用いて極めて容易に行ない得ること、および
、斜杭角βを高精度で検出するのが容易でないという理
由に基づいている。
しかし、このような評価関数を用いた制御は、次のよう
な欠点がある。上述の評価関数月それ自体できるだけ小
さな操作量によりできるだけ小さい偏差に留めようとい
う意味、船舶の制御に適用して言い換えると、船舶が直
進状態からずれたとき、それを直進状態に戻すのに最小
の月で戻すことができるという意味しか持っていない。
な欠点がある。上述の評価関数月それ自体できるだけ小
さな操作量によりできるだけ小さい偏差に留めようとい
う意味、船舶の制御に適用して言い換えると、船舶が直
進状態からずれたとき、それを直進状態に戻すのに最小
の月で戻すことができるという意味しか持っていない。
したがって、後で詳しく述べるような外乱を受けて航行
する船舶における真の意味での省エネルギ機能を有して
いない。すなわち、船体は舵を取ったときのみならず外
乱を横または斜めから受けたときも斜杭するのであり、
斜杭角βを舵角δまたは方位角速度φのみで近似すると
、操舵に無関係な斜杭が生じた場合上記の返信は全く不
合理となり、燃料消II!、量が最小となる適切な操舵
制御がなされないことになる。
する船舶における真の意味での省エネルギ機能を有して
いない。すなわち、船体は舵を取ったときのみならず外
乱を横または斜めから受けたときも斜杭するのであり、
斜杭角βを舵角δまたは方位角速度φのみで近似すると
、操舵に無関係な斜杭が生じた場合上記の返信は全く不
合理となり、燃料消II!、量が最小となる適切な操舵
制御がなされないことになる。
(ハ)発明の目的
本発明は上述の欠点を解消するためになされたもので、
コントローラにより舵角を調整し目標針路に船舶を航行
させる場合、上述した外乱をも考慮し、航行中ある目標
に向かう航程において燃料消費量を最小に維持する省エ
ネルギ航行を可能にした船舶の自動操舵装置を提供する
ことを目的とする。
コントローラにより舵角を調整し目標針路に船舶を航行
させる場合、上述した外乱をも考慮し、航行中ある目標
に向かう航程において燃料消費量を最小に維持する省エ
ネルギ航行を可能にした船舶の自動操舵装置を提供する
ことを目的とする。
なお、本発明においては、省エネルギを可能にする新し
い評価関数を見mすと共に、斜杭角の検出を精度良く行
なうことができるようになったので、実現できたもので
あるが、後者に関しては本、発明の要旨とするところで
ないので説明を省く。
い評価関数を見mすと共に、斜杭角の検出を精度良く行
なうことができるようになったので、実現できたもので
あるが、後者に関しては本、発明の要旨とするところで
ないので説明を省く。
なお、前者の新しい評価関数の詳細については後述する
が、概略は以下の通りである。従来の評価関数月は前述
したように方位角を小さくすることを目的としているが
、船舶がある目標までたとえ方位角が大きい状態で航行
しても燃料消費量が最小になれば何ら差し支えないとい
うことを意味している。本発明でtよそれに採用される
評価関数を最小にすることによって、推進エネルギの増
加を最小に留めることができるという知見に基づくもの
である。
が、概略は以下の通りである。従来の評価関数月は前述
したように方位角を小さくすることを目的としているが
、船舶がある目標までたとえ方位角が大きい状態で航行
しても燃料消費量が最小になれば何ら差し支えないとい
うことを意味している。本発明でtよそれに採用される
評価関数を最小にすることによって、推進エネルギの増
加を最小に留めることができるという知見に基づくもの
である。
(ニ)発明の構成
本発明の構成を第1図を参照して説明すると、コントロ
ーラ4により舵角を調整し目標針路に船舶を航行させる
場合、航行中のある状態における船体2の斜杭角を加味
して波・風などによる外乱の強さを推定し、この外乱の
強さから燃料消費量に関する評価関数値を最小とすると
きのコントローラ4のゲインを推定し、このゲインを有
するコントローラ4により船舶を目標に向かって航行さ
せる自動操舵装置としたことで、航行中に検出された舵
角、方位角、斜読角からそのときの船体の動特性を推定
する動特性演算手段5と、検出された方位角、斜読角か
ら外乱の強さを推定する外乱演算手段6と、この外乱演
算手段6および動特性演算手段5から出力される外乱の
強さと船体2の動特性を用いてコントローラ4で制御さ
れる船舶の航行状態をシミュレーションする航行シミュ
レーション部10と、この航行シミュレーション部10
から出力される舵角、方位角、方位角速度、斜読角に基
づいて評価関数値を計算する評価関数演算部11と、こ
の評価関数値が最小となる航行シミュレーション部11
のコントローラ13のゲインを探索する最適利得探索部
12とを有すると共に、コントローラ4にその最適ゲイ
ンを出力するゲイン出力手段7と、を具備する船舶の自
動操舵装置1である。
ーラ4により舵角を調整し目標針路に船舶を航行させる
場合、航行中のある状態における船体2の斜杭角を加味
して波・風などによる外乱の強さを推定し、この外乱の
強さから燃料消費量に関する評価関数値を最小とすると
きのコントローラ4のゲインを推定し、このゲインを有
するコントローラ4により船舶を目標に向かって航行さ
せる自動操舵装置としたことで、航行中に検出された舵
角、方位角、斜読角からそのときの船体の動特性を推定
する動特性演算手段5と、検出された方位角、斜読角か
ら外乱の強さを推定する外乱演算手段6と、この外乱演
算手段6および動特性演算手段5から出力される外乱の
強さと船体2の動特性を用いてコントローラ4で制御さ
れる船舶の航行状態をシミュレーションする航行シミュ
レーション部10と、この航行シミュレーション部10
から出力される舵角、方位角、方位角速度、斜読角に基
づいて評価関数値を計算する評価関数演算部11と、こ
の評価関数値が最小となる航行シミュレーション部11
のコントローラ13のゲインを探索する最適利得探索部
12とを有すると共に、コントローラ4にその最適ゲイ
ンを出力するゲイン出力手段7と、を具備する船舶の自
動操舵装置1である。
(ホ)実施例
以下、本発明をその実施例に基づいて詳細に説明する。
第1図は本発明の実施例である船舶の自動操舵装置1の
全体系統図で、これは、船体2の舵3に舵角を指令する
コントローラ4と、船体2の動特性を表わす運動微係数
を推定する動特性演算手段5と、外乱の強さを推定する
外乱演算手段6と、最適ゲインを出力するゲイン出力手
段7とを、その主たる構成としている。
全体系統図で、これは、船体2の舵3に舵角を指令する
コントローラ4と、船体2の動特性を表わす運動微係数
を推定する動特性演算手段5と、外乱の強さを推定する
外乱演算手段6と、最適ゲインを出力するゲイン出力手
段7とを、その主たる構成としている。
詳述すると、前記動特性演算手段5は、船舶が航行中の
ある状態下で船速を変更した時などからある一定の時間
間隔で検出された舵角δ、方位角ψ、斜読角βと、方位
角ψより求められた方位角速度Φとから、そのときの船
体2の船速、積荷や燃料残存状態により定まる動特性を
推定するもので、船体2の運動方程式に上記値を与える
ことにより、その時点における運動微係数を演算し、こ
の係数をゲイン出力手段7の後述する航行シミュレーシ
ョン部10に出力するものである。
ある状態下で船速を変更した時などからある一定の時間
間隔で検出された舵角δ、方位角ψ、斜読角βと、方位
角ψより求められた方位角速度Φとから、そのときの船
体2の船速、積荷や燃料残存状態により定まる動特性を
推定するもので、船体2の運動方程式に上記値を与える
ことにより、その時点における運動微係数を演算し、こ
の係数をゲイン出力手段7の後述する航行シミュレーシ
ョン部10に出力するものである。
前記外乱演算手段6は、上述と同様の方位角ψと斜読角
βおよび求められた方位角速度ψとから、そのときの波
・風などによる外乱の強さを演算し、その外乱の強さを
航行シミュレーション部10に出力するものである。な
お、この外乱の強さとは、波・風により船体2に生じる
旋回モーメントと横力(サイドフォース)などであり、
これらを周波数領域でスペクトル処理するものであって
も、また、時間領域で処理するものであってもよいが、
図では前者によるものを示している。この場合、船体運
動スペクトル密度計算部8と外乱スペクトル密度計算部
9が設けられ、前者においては、方位角ψ、斜読角β、
方位角速度φより、周波数領域でそれぞれの運動スペク
トル密度が計算され、これらのスペクトル密度を受けて
後者では周波数領域で外乱スペクトル密度が計算される
ようになっている。そして、これらの外乱スペクトル密
度が航行シミュレーション部10に出力され、前述した
運動微係数と共に幾種類もの仮想操舵における航行状態
のシミュレーションに供される。
βおよび求められた方位角速度ψとから、そのときの波
・風などによる外乱の強さを演算し、その外乱の強さを
航行シミュレーション部10に出力するものである。な
お、この外乱の強さとは、波・風により船体2に生じる
旋回モーメントと横力(サイドフォース)などであり、
これらを周波数領域でスペクトル処理するものであって
も、また、時間領域で処理するものであってもよいが、
図では前者によるものを示している。この場合、船体運
動スペクトル密度計算部8と外乱スペクトル密度計算部
9が設けられ、前者においては、方位角ψ、斜読角β、
方位角速度φより、周波数領域でそれぞれの運動スペク
トル密度が計算され、これらのスペクトル密度を受けて
後者では周波数領域で外乱スペクトル密度が計算される
ようになっている。そして、これらの外乱スペクトル密
度が航行シミュレーション部10に出力され、前述した
運動微係数と共に幾種類もの仮想操舵における航行状態
のシミュレーションに供される。
前記ゲイン出力手段7は、航行状態をシミュレーション
する航行シミュレーション部10と、評価関数演算部1
1と、最適利得探索部12とを有し、コントローラ4に
その最適ゲインを出力するものである。その航行シミュ
レーション部10は、数学モデルのコントローラ13、
舵14、船体15とを用いて演算するものであって、前
述した船体2の運動微係数と外乱スペクトル密度を用い
て、コントローラ13による幾種類もの舵角調整信号、
すなわち、各ゲインにおける異なる値でもってそれぞれ
に対応する幾通りものシミュレーションを行なうもので
ある。そして、この各シミュレーションにおける舵角δ
S、方位角ψS、斜読角βSおよび方位角ψSより求め
られる方位角速度tfis、または、それぞれのスペク
トル密度および相互スペクトル密度を評価関数演算部1
1に出力する。
する航行シミュレーション部10と、評価関数演算部1
1と、最適利得探索部12とを有し、コントローラ4に
その最適ゲインを出力するものである。その航行シミュ
レーション部10は、数学モデルのコントローラ13、
舵14、船体15とを用いて演算するものであって、前
述した船体2の運動微係数と外乱スペクトル密度を用い
て、コントローラ13による幾種類もの舵角調整信号、
すなわち、各ゲインにおける異なる値でもってそれぞれ
に対応する幾通りものシミュレーションを行なうもので
ある。そして、この各シミュレーションにおける舵角δ
S、方位角ψS、斜読角βSおよび方位角ψSより求め
られる方位角速度tfis、または、それぞれのスペク
トル密度および相互スペクトル密度を評価関数演算部1
1に出力する。
ここで、本発明で採用される新しい評価関数について説
明する。新しい評価関数J2は、外乱による斜読角を加
味した推進エネルギの増加率を示す下記の式(6)で与
えられるものである。この評価関数J2は次のような論
理により導かれる。
明する。新しい評価関数J2は、外乱による斜読角を加
味した推進エネルギの増加率を示す下記の式(6)で与
えられるものである。この評価関数J2は次のような論
理により導かれる。
第2図に示すように船体2の座標原点をその重心Gに採
ると、外乱下で蛇行したり針路からのず0 れを起こしなから一航海を行なうに必要な推進工と表わ
される。なお、Tは推力、Uは船体2の前進方向速度す
なわち船速Uの斜読角β余弦、τは一航海に要する時間
である。ところで、船舶が舵と表わせ、船体2が蛇行や
斜航することなく、舵角を零に保ったま\船首方向に理
想的に航行する場合の推進エネルギである。なお、τ0
は一般にτより小さい。これらの式(2)および(3)
から推進エネルギの増加率J2は、 J 2−(E−Eo ) /Eo ・−−−−−
−−−(4)と表わされ、自動操舵により波・風・潮流
などの外乱下を航行するときに必要な推進エネルギと、
平水状態を直進するときに要するそれとを比較した場合
の増加率である。この式(4)に船体2の前進J2−μ
m (ψ−β)”十μ2 (β+Tφ)φ十μ3 β2
+μ4 δ” −−−−−−(51表わされるも
のである。なお、μm、μ2、μ3、μ4は定数、Tは
距離を表わす定数である。加えて、上記式(5)の各項
の物理的意味は次の通りである。第1項の(ψ−β)は
、船舶が蛇行することによって航路が直進時より長くな
ることによる推進エネルギの損失、第2項の(β十rf
i)ψは、船体が遠心力すなわち慣性抵抗を受けること
によr っで生ずる推進エネルギの損失、第3項のβは、船体が
斜航することによって受ける抵抗による推] 進エネルギの損失、第4項のδは、操舵操作による推進
エネルギの損失を表わしている。
ると、外乱下で蛇行したり針路からのず0 れを起こしなから一航海を行なうに必要な推進工と表わ
される。なお、Tは推力、Uは船体2の前進方向速度す
なわち船速Uの斜読角β余弦、τは一航海に要する時間
である。ところで、船舶が舵と表わせ、船体2が蛇行や
斜航することなく、舵角を零に保ったま\船首方向に理
想的に航行する場合の推進エネルギである。なお、τ0
は一般にτより小さい。これらの式(2)および(3)
から推進エネルギの増加率J2は、 J 2−(E−Eo ) /Eo ・−−−−−
−−−(4)と表わされ、自動操舵により波・風・潮流
などの外乱下を航行するときに必要な推進エネルギと、
平水状態を直進するときに要するそれとを比較した場合
の増加率である。この式(4)に船体2の前進J2−μ
m (ψ−β)”十μ2 (β+Tφ)φ十μ3 β2
+μ4 δ” −−−−−−(51表わされるも
のである。なお、μm、μ2、μ3、μ4は定数、Tは
距離を表わす定数である。加えて、上記式(5)の各項
の物理的意味は次の通りである。第1項の(ψ−β)は
、船舶が蛇行することによって航路が直進時より長くな
ることによる推進エネルギの損失、第2項の(β十rf
i)ψは、船体が遠心力すなわち慣性抵抗を受けること
によr っで生ずる推進エネルギの損失、第3項のβは、船体が
斜航することによって受ける抵抗による推] 進エネルギの損失、第4項のδは、操舵操作による推進
エネルギの損失を表わしている。
上記の式(5)をさらに整理すると、評価関数J2は、
次式で表わされる。
次式で表わされる。
J2−λ17(λ2yトλ30−+λ4y了
+λ5ψ+λ6βψ 〜−−−−−−−−T6
)なお、δは舵角、βは斜読角、ψは方位角、↓は方位
角速度、λ1〜λ6は定数である。
)なお、δは舵角、βは斜読角、ψは方位角、↓は方位
角速度、λ1〜λ6は定数である。
1
再び、前述したゲイン出力手段7の残りの構成を述べる
。最適利得探索部12は、前記評価関数演算部11で演
算される評価関数値を最小にするように、例えば、山登
り法などの手法を用いてコントローラ13のゲインを探
索するもので、そのゲインを最適ゲインとして船体2を
操舵調整するコントローラ4に出力するものである。な
お、コントローラ4へのゲイン入力のために、最適制御
ゲイン自動設定部16が必要に応じてゲイン出力手段7
の内部または外部に設けられる。
。最適利得探索部12は、前記評価関数演算部11で演
算される評価関数値を最小にするように、例えば、山登
り法などの手法を用いてコントローラ13のゲインを探
索するもので、そのゲインを最適ゲインとして船体2を
操舵調整するコントローラ4に出力するものである。な
お、コントローラ4へのゲイン入力のために、最適制御
ゲイン自動設定部16が必要に応じてゲイン出力手段7
の内部または外部に設けられる。
このような実施例によれば、次のようにして燃料消費量
が最小となるよう船舶を自動操舵させることができる。
が最小となるよう船舶を自動操舵させることができる。
第1図において、まず、目標に向かう針路ψ0がコント
ローラ4に入力される。いま、コントローラ4において
例えばPID制御がなされるとすると、そこに予め入力
されているPIDの各ゲインKpO+ Kio、 Kd
oにより、舵3が舵角δ0となるよう操舵指令がなされ
る。その結果、船体2は舵角δ0にしたがって目標針路
φ0を採って航行3 2 する。ところが、このような航行中に波・風などの外乱
が船体2に作用すると、第2図に示すように船体2が針
路ψ0より方位角ψ外れた針路ψbを採ったり、また、
場合によっては横滑りを起こして左右方向速度Vが生じ
、これによる斜読角βが発生したりする。ところで、船
体2に作用する外乱は急変することもあるが、通常はあ
る一定時間持続するのて゛、その外乱が持続する間にお
ける航行中設定されている船速Uの下で、ある目標に向
かう間の燃料消費量を最小にする必要がある。
ローラ4に入力される。いま、コントローラ4において
例えばPID制御がなされるとすると、そこに予め入力
されているPIDの各ゲインKpO+ Kio、 Kd
oにより、舵3が舵角δ0となるよう操舵指令がなされ
る。その結果、船体2は舵角δ0にしたがって目標針路
φ0を採って航行3 2 する。ところが、このような航行中に波・風などの外乱
が船体2に作用すると、第2図に示すように船体2が針
路ψ0より方位角ψ外れた針路ψbを採ったり、また、
場合によっては横滑りを起こして左右方向速度Vが生じ
、これによる斜読角βが発生したりする。ところで、船
体2に作用する外乱は急変することもあるが、通常はあ
る一定時間持続するのて゛、その外乱が持続する間にお
ける航行中設定されている船速Uの下で、ある目標に向
かう間の燃料消費量を最小にする必要がある。
そのための操舵量を計算するために、上記方位角ψと斜
読角βが図示しない計測機器により検出され、その検出
量と前記舵角δ0が動特性演算手段5に入力される。な
お、船体2のそのときの動特性は前述したように、船速
、重心位置などにより異なるので、航行中に異なった船
速を設定したり、排水量や重心位置が異なる事情が生じ
たときなどは、その都度上述の検出が行なわれてその状
態における動特性が演算される。
読角βが図示しない計測機器により検出され、その検出
量と前記舵角δ0が動特性演算手段5に入力される。な
お、船体2のそのときの動特性は前述したように、船速
、重心位置などにより異なるので、航行中に異なった船
速を設定したり、排水量や重心位置が異なる事情が生じ
たときなどは、その都度上述の検出が行なわれてその状
態における動特性が演算される。
上述した方位角ψと斜読角βおよびそのときの4
舵角δ0が検出され、かつ、その方位角ψより方位角速
度φが計算され、これらを用いて動特性演算手段5にお
いてその外乱下における船体2の動特性を指標する運動
微係数が演算され、この微係数がゲイン出力手段7の航
行シミュレーション部10に出力されて、船体15の動
特性として扱われる。
度φが計算され、これらを用いて動特性演算手段5にお
いてその外乱下における船体2の動特性を指標する運動
微係数が演算され、この微係数がゲイン出力手段7の航
行シミュレーション部10に出力されて、船体15の動
特性として扱われる。
一方、外乱演算手段6にも同様に、前記方位角ψと斜読
角βが例えば数分ないし十数骨の一定時間ごとに入力さ
れ、かつ、方位角速度ξが計算され、これらから外乱の
強さが演算される。その過程は、まず、船体運動スペク
トル密度計算部8において、第3図(a)〜(C)に示
すように方位角ψ、斜読角β、方位角速度かの周波数に
対する運動スペクトル密度ΦPI’、Φf訃Φ忰令が計
算される。次に、このスペクトル密度を受けて外乱スペ
クトル密度計算部9では、周波数に対する外乱の強さ、
すなわち、旋回モーメントや横力などのスペクトル密度
Φyy、Φnn、Φynが第4図(a)、(b)に示す
ように推定計算される。
角βが例えば数分ないし十数骨の一定時間ごとに入力さ
れ、かつ、方位角速度ξが計算され、これらから外乱の
強さが演算される。その過程は、まず、船体運動スペク
トル密度計算部8において、第3図(a)〜(C)に示
すように方位角ψ、斜読角β、方位角速度かの周波数に
対する運動スペクトル密度ΦPI’、Φf訃Φ忰令が計
算される。次に、このスペクトル密度を受けて外乱スペ
クトル密度計算部9では、周波数に対する外乱の強さ、
すなわち、旋回モーメントや横力などのスペクトル密度
Φyy、Φnn、Φynが第4図(a)、(b)に示す
ように推定計算される。
5
この外乱のスペクトル密度φyy、Φnn、Φynがゲ
イン出力手段7に入力されると共に、すてに動特性演算
手段5から船体15に入力されている運動微係数に基づ
いて、船体2の航行状態のシミュレーションが、数学モ
デルである航行シミュレーション部10において行なわ
れる。この際、最適利得探索部12は、前記評価関数演
算部11で演算される評価関数値を最小とするように例
えば、山登り法などの手法を用いてコントローラ13の
ゲインを探索する。そして、このゲインが上述した状態
にある船体2の航行における最小燃料消費量となる最適
ゲインとして、最適制御ゲイン自動設定部16を介して
実船体2を操舵調整するコントローラ4に入力される。
イン出力手段7に入力されると共に、すてに動特性演算
手段5から船体15に入力されている運動微係数に基づ
いて、船体2の航行状態のシミュレーションが、数学モ
デルである航行シミュレーション部10において行なわ
れる。この際、最適利得探索部12は、前記評価関数演
算部11で演算される評価関数値を最小とするように例
えば、山登り法などの手法を用いてコントローラ13の
ゲインを探索する。そして、このゲインが上述した状態
にある船体2の航行における最小燃料消費量となる最適
ゲインとして、最適制御ゲイン自動設定部16を介して
実船体2を操舵調整するコントローラ4に入力される。
すなわち、コントローラ4に記憶されている前述のゲイ
ンKPo+ K to。
ンKPo+ K to。
Kdoが、最適ゲインであるKpl、 Kil、 Kd
lに置き換えられる。したがって、コントローラ4は新
しいゲインに基づいて操舵信号を発信する。船体2は上
述した外乱の強さや船体2の速度などが変わらない状態
で所望の目標に向けて航行する間、6 実質的に燃料消費量が最小となる航海がなされる。
lに置き換えられる。したがって、コントローラ4は新
しいゲインに基づいて操舵信号を発信する。船体2は上
述した外乱の強さや船体2の速度などが変わらない状態
で所望の目標に向けて航行する間、6 実質的に燃料消費量が最小となる航海がなされる。
なお、外乱の強さや船体2の状態が変れば、その都度上
述した一連の検出、演算などの処理が行なわれることは
いうまでもない。したがって、船舶が最終目的地に到達
するまでの間に、設定船速を変えてもその状態おける目
標までの航行中の燃料消費量を最小にすることができ、
総じて、−航海における燃料消費量を少なくできる。
述した一連の検出、演算などの処理が行なわれることは
いうまでもない。したがって、船舶が最終目的地に到達
するまでの間に、設定船速を変えてもその状態おける目
標までの航行中の燃料消費量を最小にすることができ、
総じて、−航海における燃料消費量を少なくできる。
以上を簡潔にまとめると、コントローラにより舵角を調
整し、目標針路に船舶を航行させる場合1、航行中のあ
る状態における船体の斜読角を加味して波・風などによ
る外乱゛の強さを推定し、この外乱の強さから燃料消費
量に関する評価関数値を最小とするときのゲインを推定
し、このゲインによりコントローラを稼働させ、船舶を
目標に向かって最小の燃料消費量で航行させることが゛
できる。
整し、目標針路に船舶を航行させる場合1、航行中のあ
る状態における船体の斜読角を加味して波・風などによ
る外乱゛の強さを推定し、この外乱の強さから燃料消費
量に関する評価関数値を最小とするときのゲインを推定
し、このゲインによりコントローラを稼働させ、船舶を
目標に向かって最小の燃料消費量で航行させることが゛
できる。
(f)発明の効果
本発明は以上の実施例の説明から判るように、コントロ
ーラにより舵角を調整し目標針路に船舶を航行させる場
合、設定された船速における航行7 中船体に作用する外乱をも考慮し、そのような外乱や船
体の状態が維持される限り、燃料消費量を最小に維持し
て船舶の真の省エネルギ航行を行なわせることができる
。したがって、外乱などに変化があれば、その都度、そ
れに応じた省エネルギ航行が行なわれ、最終目的地まで
の一航海を通じて、最も少ない燃料消費量でもって船舶
を航行さ、せることができる。
ーラにより舵角を調整し目標針路に船舶を航行させる場
合、設定された船速における航行7 中船体に作用する外乱をも考慮し、そのような外乱や船
体の状態が維持される限り、燃料消費量を最小に維持し
て船舶の真の省エネルギ航行を行なわせることができる
。したがって、外乱などに変化があれば、その都度、そ
れに応じた省エネルギ航行が行なわれ、最終目的地まで
の一航海を通じて、最も少ない燃料消費量でもって船舶
を航行さ、せることができる。
、4、図面の簡単な説明
第1図は本発明の実施例である船舶の自動操舵装置の全
体系統図、警2図は外乱を受けている船体の力学モデル
図、第3TI!l1(a)〜(C)は運動スペクトル図
、第4図(a)、(b)は外乱のスペクトル図である。
体系統図、警2図は外乱を受けている船体の力学モデル
図、第3TI!l1(a)〜(C)は運動スペクトル図
、第4図(a)、(b)は外乱のスペクトル図である。
1−自動操舵装置、2−船体、4−コントロ“−ラ、5
−動特性演算手段、6−外乱演算手段、7−ゲイン出力
手段、1〇−航行シミュレーション部、11−評価関数
演算部、12−最適利得探索部、13−数学モデルのコ
ントローラ、68 0、δa−舵角、β−斜斜角角ψ−−位角特許出願人
川崎重工業株式会社 代理人 弁理士 吉村勝俊(はが1名)第1頁の続き 0発 明 者 松本反相 神戸市中央区東用崎町3丁目1 番1号川崎重工業株式会社神戸 工場内 0発 明 者 西義和 明石市川崎町1番1号川崎重工 業株式会社技術研究所内 手続補正書(自発) 1、事件の表示 昭和58年特許願第092088号 2、発明の名称 船舶の自動操舵装置 3、?!正をする者 事件との関係 特許出願人 代表者長谷用 謙浩 4、代理人 明細書の「発明の詳細な説明」の欄 6、補正の内容 (1〕 発明の詳細な説明の欄 と訂正する。
−動特性演算手段、6−外乱演算手段、7−ゲイン出力
手段、1〇−航行シミュレーション部、11−評価関数
演算部、12−最適利得探索部、13−数学モデルのコ
ントローラ、68 0、δa−舵角、β−斜斜角角ψ−−位角特許出願人
川崎重工業株式会社 代理人 弁理士 吉村勝俊(はが1名)第1頁の続き 0発 明 者 松本反相 神戸市中央区東用崎町3丁目1 番1号川崎重工業株式会社神戸 工場内 0発 明 者 西義和 明石市川崎町1番1号川崎重工 業株式会社技術研究所内 手続補正書(自発) 1、事件の表示 昭和58年特許願第092088号 2、発明の名称 船舶の自動操舵装置 3、?!正をする者 事件との関係 特許出願人 代表者長谷用 謙浩 4、代理人 明細書の「発明の詳細な説明」の欄 6、補正の内容 (1〕 発明の詳細な説明の欄 と訂正する。
■ 同第12W第1行目
] 1
「十μ3β十μ4δ −−〜−−−−−−(5)」と
あるのを、 「 十 μ 3 β + メ! 4 δe −
−・−−−(F+) jと訂正する。
あるのを、 「 十 μ 3 β + メ! 4 δe −
−・−−−(F+) jと訂正する。
■ 同第12頁第5行目「定数である。」の後に、「ま
た、δeは有効舵角であり、a、bを定数としてδe−
δ−aβ−bjと表わすことができる。
た、δeは有効舵角であり、a、bを定数としてδe−
δ−aβ−bjと表わすことができる。
」を挿入する。
■ 同第12頁第13行目。
] −1
「δ]とあるのを、「δe」と訂正する。
■ 同第12頁第18j〒目
「+λ5ψ十λ6βψ −−−−−−(61J−・−
ニー+61 J と訂正する。
ニー+61 J と訂正する。
■ 同第12頁第20行目
「λ1〜λ6」とあるのを、「λ1〜λ8」と訂正する
。
。
Claims (1)
- (1) 目標方位と実方位の差に応じて舵を取らせる閉
ループを構成する自動操舵装置において、航行中に検出
された舵角、方位角、斜杭角からそのときの船体の動特
性を推定する動特性演算手段と、 検出された方位角、斜杭角から外乱の強さを推定する外
乱演算手段と、 この外乱演算手段および前記動特性演算手段から出力さ
れる外乱の強さと船体の動特性を用いて、前記閉ループ
の船舶の航行状態をシミュレーションする航行シミュレ
ーション部と、この航行シミュレーション部から出力さ
れる舵角、方位角、方位角速度、斜杭角の関数として示
される推進エネルギ増加率を表わす評価関数値を計算す
る評価関数演算部と、この評価関数値を最小とする航行
シミュレーション部のコントローラのゲインを探索する
最適利f!を探索部とを有すると共に、前記閉ループを
構成するコントローラにその@適ゲインを出力するゲイ
ン出力手段と、 を具備することを特徴とする船舶の自動操舵装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP58092088A JPS59220496A (ja) | 1983-05-24 | 1983-05-24 | 船舶の自動操舵装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP58092088A JPS59220496A (ja) | 1983-05-24 | 1983-05-24 | 船舶の自動操舵装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPS59220496A true JPS59220496A (ja) | 1984-12-11 |
JPH0330557B2 JPH0330557B2 (ja) | 1991-04-30 |
Family
ID=14044682
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP58092088A Granted JPS59220496A (ja) | 1983-05-24 | 1983-05-24 | 船舶の自動操舵装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPS59220496A (ja) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61139594A (ja) * | 1984-12-12 | 1986-06-26 | Tokyo Keiki Co Ltd | 舵の自動天候調整装置 |
JPH10332416A (ja) * | 1997-06-02 | 1998-12-18 | Kansai Electric Power Co Inc:The | 自動移動体、自動運航船およびダム堆砂自動測量船 |
JP2006188095A (ja) * | 2004-12-28 | 2006-07-20 | Tokimec Inc | 船舶用自動操舵装置 |
JP2011042284A (ja) * | 2009-08-21 | 2011-03-03 | Yokogawa Denshikiki Co Ltd | 自動操舵装置及び方法 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5067541B2 (ja) | 2007-03-30 | 2012-11-07 | Tdk株式会社 | 誘電体磁器組成物、複合電子部品および積層セラミックコンデンサ |
JP4618383B2 (ja) | 2008-05-12 | 2011-01-26 | Tdk株式会社 | 誘電体磁器組成物、積層複合電子部品、積層コモンモードフィルタ、積層セラミックコイルおよび積層セラミックコンデンサ |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5511968A (en) * | 1978-07-14 | 1980-01-28 | Mitsubishi Heavy Ind Ltd | Navigation control system |
-
1983
- 1983-05-24 JP JP58092088A patent/JPS59220496A/ja active Granted
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5511968A (en) * | 1978-07-14 | 1980-01-28 | Mitsubishi Heavy Ind Ltd | Navigation control system |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61139594A (ja) * | 1984-12-12 | 1986-06-26 | Tokyo Keiki Co Ltd | 舵の自動天候調整装置 |
JPH0341398B2 (ja) * | 1984-12-12 | 1991-06-21 | ||
JPH10332416A (ja) * | 1997-06-02 | 1998-12-18 | Kansai Electric Power Co Inc:The | 自動移動体、自動運航船およびダム堆砂自動測量船 |
JP2006188095A (ja) * | 2004-12-28 | 2006-07-20 | Tokimec Inc | 船舶用自動操舵装置 |
JP4669283B2 (ja) * | 2004-12-28 | 2011-04-13 | 東京計器株式会社 | 船舶用自動操舵装置 |
JP2011042284A (ja) * | 2009-08-21 | 2011-03-03 | Yokogawa Denshikiki Co Ltd | 自動操舵装置及び方法 |
Also Published As
Publication number | Publication date |
---|---|
JPH0330557B2 (ja) | 1991-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Fossen et al. | Direct and indirect adaptive integral line‐of‐sight path‐following controllers for marine craft exposed to ocean currents | |
Sarda et al. | Station-keeping control of an unmanned surface vehicle exposed to current and wind disturbances | |
Xiang et al. | Smooth transition of AUV motion control: From fully-actuated to under-actuated configuration | |
Lekkas et al. | Line-of-sight guidance for path following of marine vehicles | |
Fossen et al. | Nonlinear output feedback control of dynamically positioned ships using vectorial observer backstepping | |
Borhaug et al. | Cross-track control for underactuated autonomous vehicles | |
Xu et al. | L1 adaptive backstepping control for path-following of underactuated marine surface ships | |
CN106444806B (zh) | 基于生物速度调节的欠驱动auv三维轨迹跟踪控制方法 | |
CN115113524B (zh) | 一种基于干预lvs制导的asv多端口事件触发路径跟踪控制方法 | |
JP3038209B1 (ja) | 自動方位設定方法とその装置 | |
Liu | Improved ELOS based path following control for underactuated surface vessels with roll constraint | |
Li et al. | Nonlinear heading control of an autonomous underwater vehicle with internal actuators | |
CN116736709A (zh) | 一种海洋机器人的动态补偿型自抗扰艏向控制方法 | |
CN115014355A (zh) | 一种双体无人船的定点返航调控方法及装置 | |
JPS59220496A (ja) | 船舶の自動操舵装置 | |
Liu et al. | Hierarchical robust path following control of fully submerged hydrofoil vessels | |
Liu et al. | Ship heading control with speed keeping via a nonlinear disturbance observer | |
JP3033571B1 (ja) | 潜水船の深度保持制御方式 | |
Huang et al. | SHSA-based adaptive roll-safety 3D tracking control of a X-Rudder AUV with actuator dynamics | |
CN111506090A (zh) | 一种水下机器人深度区间控制方法和系统 | |
Han et al. | A Ship Path Tracking Control Method Using a Fuzzy Control Integrated Line-of-Sight Guidance Law | |
JP2008213682A (ja) | 船舶用自動操舵装置 | |
JP6487365B2 (ja) | 船舶用自動操舵装置 | |
Unar et al. | A neural network based heading and position control system of a ship | |
CN117519136B (zh) | 一种考虑大曲率转弯的无人艇路径跟踪方法 |