JP2021189909A - Autonomous travel system for work vehicle - Google Patents
Autonomous travel system for work vehicle Download PDFInfo
- Publication number
- JP2021189909A JP2021189909A JP2020096400A JP2020096400A JP2021189909A JP 2021189909 A JP2021189909 A JP 2021189909A JP 2020096400 A JP2020096400 A JP 2020096400A JP 2020096400 A JP2020096400 A JP 2020096400A JP 2021189909 A JP2021189909 A JP 2021189909A
- Authority
- JP
- Japan
- Prior art keywords
- work vehicle
- steering
- control
- travel
- equation
- 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
- 238000001514 detection method Methods 0.000 claims abstract description 24
- 238000011156 evaluation Methods 0.000 claims abstract description 13
- 238000005259 measurement Methods 0.000 abstract description 9
- 238000004891 communication Methods 0.000 description 25
- 238000000034 method Methods 0.000 description 20
- 241000196324 Embryophyta Species 0.000 description 13
- 239000011159 matrix material Substances 0.000 description 13
- 238000010586 diagram Methods 0.000 description 11
- 238000012937 correction Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 10
- 238000012545 processing Methods 0.000 description 10
- 230000033001 locomotion Effects 0.000 description 9
- 230000008859 change Effects 0.000 description 8
- 238000005070 sampling Methods 0.000 description 8
- 230000004044 response Effects 0.000 description 4
- 230000003111 delayed effect Effects 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- ZZUFCTLCJUWOSV-UHFFFAOYSA-N furosemide Chemical compound C1=C(Cl)C(S(=O)(=O)N)=CC(C(O)=O)=C1NCC1=CC=CO1 ZZUFCTLCJUWOSV-UHFFFAOYSA-N 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000009347 mechanical transmission Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
Images
Abstract
Description
本発明は、作業車両の自律走行システムに関し、特に、前輪が操舵輪である作業車両を自律走行させる作業車両の自律走行システムに関するものである。 The present invention relates to an autonomous traveling system of a work vehicle, and more particularly to an autonomous traveling system of a work vehicle in which a work vehicle whose front wheel is a steering wheel is autonomously traveled.
例えば圃場や枕地で予め定められた走行経路に沿って前進直進と前進旋回とを繰り返す等といった具合に、予め定められ走行パターンにて作業車両を走行させる場合には、オペレータの負担を軽減するべく、作業車両を自律走行させることが従来から行われている。 For example, when the work vehicle is driven in a predetermined driving pattern, such as repeating forward straight movement and forward turning along a predetermined traveling route in a field or a headland, the burden on the operator is reduced. Therefore, it has been conventionally practiced to autonomously drive a work vehicle.
例えば特許文献1には、圃場内に予め設定された走行経路上の前方目標点と作業車両とを結ぶ線分と、走行経路と、がなす前方目標角度を算出するとともに、作業車両に配置された角速度センサから取得したヨーレートを実現する演算タイヤ角度を算出し、前方目標角度と演算タイヤ角度との比較結果に基づいて舵角を制御する、前方注視方式を採用した自律操舵装置が開示されている。
For example, in
ところで、特許文献1のもののように、前方注視方式を採用した自律操舵装置では、作業車両に設けられたGNSSアンテナが走行経路上に位置するようにフィードバック補正を行うことが、具体的には、操舵制御後のGNSSアンテナの位置を検出し、検出結果に基づいて舵角の補正を行うのが一般的である。
By the way, in the autonomous steering device adopting the forward gaze method as in
しかしながら、前輪が操舵輪である作業車両の後進旋回時に、かかる前方注視方式を採用した場合には、以下のような問題を生ずるおそれがある。 However, when the forward gaze method is adopted when the front wheel is a steering wheel for a reverse turn of a work vehicle, the following problems may occur.
すなわち、旋回の向きと前輪の操舵の向きとが一致する前進旋回とは異なり、後進旋回の場合には、旋回の向きとは反対の向きに前輪を操舵する必要がある。それ故、後進旋回を開始すると、前輪側が例えば円弧状の走行経路の外側に膨らむとともに、GNSSアンテナ位置が走行経路の内側に入ることになるため、旋回の向きと同じ向きに前輪を操舵する補正が入ることになる。すると、今度は、GNSSアンテナ位置が走行経路の外側に出ることになるため、旋回の向きとは反対の向きに前輪を操舵する補正が入ることになる。このような操舵制御とフィードバック補正とが繰り返されることで、操舵制御が発散し、作業車両が蛇行するという問題が生じることになる。 That is, unlike forward turning in which the turning direction and the steering direction of the front wheels match, in the case of reverse turning, it is necessary to steer the front wheels in the direction opposite to the turning direction. Therefore, when the reverse turn is started, the front wheel side bulges to the outside of the arcuate traveling path, and the GNSS antenna position enters the inside of the traveling path. Will enter. Then, since the GNSS antenna position goes out of the traveling path this time, the correction for steering the front wheels in the direction opposite to the turning direction is entered. By repeating such steering control and feedback correction, the steering control diverges, causing a problem that the work vehicle meanders.
本発明はかかる点に鑑みてなされたものであり、その目的とするところは、前輪が操舵輪である作業車両を、自律走行で後進旋回させる場合に、作業車両が蛇行するのを抑える技術を提供することにある。 The present invention has been made in view of this point, and an object of the present invention is to provide a technique for suppressing meandering of a work vehicle whose front wheels are steered wheels when the work vehicle is autonomously turned backward. To provide.
前記目的を達成するため、本発明に係る作業車両の自律走行システムでは、予測モデルを用いたモデル予測制御により、所定の制御周期毎に最適な制御量を算出するようにしている。 In order to achieve the above object, in the autonomous traveling system of the work vehicle according to the present invention, the optimum control amount is calculated for each predetermined control cycle by model predictive control using a predictive model.
具体的には、本発明は、前輪が操舵輪である作業車両を自律走行させる作業車両の自律走行システムを対象としている。 Specifically, the present invention is intended for an autonomous travel system of a work vehicle in which a work vehicle whose front wheel is a steering wheel is autonomously traveled.
そして、この作業車両の自律走行システムは、上記作業車両が走行する目標となる走行経路を生成する経路生成手段と、上記作業車両に設定された制御中心の位置を検出可能な位置検出手段と、上記前輪の操舵角の制御量を演算して出力する演算手段と、上記演算手段によって出力された制御量に基づいて前輪を制御する操舵手段と、を備え、上記演算手段は、後進旋回時には、上記走行経路からの上記制御中心の偏差と上記操舵角の制御量との相関を状態方程式で表した予測モデルに基づき、所定の制御周期毎に、所定の評価関数を用いて最適な制御量を所定のステップ先まで算出し、第1番目のステップの制御量のみを出力するとともに、当該制御量に対応する上記制御中心の予測位置と、上記位置検出手段によって検出された当該制御中心の実測位置とに基づいて、当該予測モデルを補正するように構成されていることを特徴とするものである。 The autonomous traveling system of the work vehicle includes a route generation means for generating a target travel route on which the work vehicle travels, a position detection means capable of detecting the position of the control center set in the work vehicle, and a position detection means. The calculation means includes a calculation means for calculating and outputting a control amount of the steering angle of the front wheels and a steering means for controlling the front wheels based on the control amount output by the calculation means. Based on a prediction model that expresses the correlation between the deviation of the control center from the travel path and the control amount of the steering angle by a state equation, the optimum control amount is obtained by using a predetermined evaluation function for each predetermined control cycle. It calculates up to a predetermined step destination, outputs only the control amount of the first step, predicts the predicted position of the control center corresponding to the control amount, and actually measures the control center detected by the position detecting means. It is characterized in that it is configured to correct the prediction model based on the above.
なお、本発明において「自律走行」とは、少なくとも操舵を自律的に行って、オペレータの運転操作なしで作業車両を走行させることを意味する。 In the present invention, "autonomous driving" means that the work vehicle is driven by at least steering autonomously without any driving operation of the operator.
この構成によれば、予測モデルに基づき所定の制御周期毎に、最適な制御量を所定のステップ先まで算出することから、換言すると、数ステップ先を予測した上で最適な制御量を算出することから、制御中心の予測位置と実測位置とに基づいて予測モデルをフィードバック補正する際に、操舵制御が発散するのを抑えることができる。したがって、前輪が操舵輪である作業車両を自律走行で後進旋回させる場合にも、作業車両が蛇行するのを抑えることができる。 According to this configuration, the optimum control amount is calculated up to a predetermined step destination for each predetermined control cycle based on the prediction model. In other words, the optimum control amount is calculated after predicting several steps ahead. Therefore, it is possible to suppress the divergence of steering control when the prediction model is feedback-corrected based on the predicted position of the control center and the actually measured position. Therefore, even when the work vehicle whose front wheel is the steering wheel is autonomously turned backward, it is possible to suppress the work vehicle from meandering.
また、上記作業車両の自律走行システムでは、上記制御中心は、上記作業車両の後輪の車軸における車幅方向の中心位置に設定されており、上記位置検出手段は、上記作業車両に設けられたGNSSアンテナの現在位置に基づいて、上記制御中心の現在位置を検出するように構成されていることが好ましい。 Further, in the autonomous traveling system of the work vehicle, the control center is set at the center position in the vehicle width direction on the axle of the rear wheel of the work vehicle, and the position detecting means is provided in the work vehicle. It is preferable that the control center is configured to detect the current position of the control center based on the current position of the GNSS antenna.
作業車両を自律走行させる場合には、通常、GNSSアンテナの位置で作業車両を制御するところ、この構成によれば、制御中心をGNSSアンテナの位置ではなく、作業車両の後輪の車軸における車幅方向の中心位置に設定することから、換言すると、仮想のGNSSアンテナの位置(制御中心)を、後進旋回時における車両中心に設定することから、蛇行の影響をより一層抑えることができる。 When the work vehicle is autonomously driven, the work vehicle is usually controlled at the position of the GNSS antenna. According to this configuration, the control center is not the position of the GNSS antenna but the width of the rear wheel axle of the work vehicle. Since it is set at the center position in the direction, in other words, the position (control center) of the virtual GNSS antenna is set at the center of the vehicle at the time of reverse turning, so that the influence of meandering can be further suppressed.
ところで、前進直進してきた作業車両を、後進旋回させて直角に方向転換させるような場合には、作業車両を前進直進からいきなり後進旋回させる訳ではなく、短い後進直進を経て後進旋回を開始するのが一般的である。このように後進直進から後進旋回へ切り替える場合には、曲率が0から1/曲率半径へ急激に変化するが、操舵手段にはダイナミクス(動的特性)があるため、不可避的に応答遅れが生じてしまい、旋回のタイミングが遅れることになる。かかる応答遅れを抑えるために、制御量を大きくすると、実際の操舵角が大きくオーバーシュートしてしまうおそれがある。なお、このような問題は、後進旋回から前進直進または後進直進へ切り替える場合にも生じ得る。 By the way, when a work vehicle that has traveled forward and straight is turned backward and turned at a right angle, the work vehicle is not suddenly turned backward from forward straight, but starts a reverse turn after a short reverse straight. Is common. When switching from reverse straight to reverse turning in this way, the curvature suddenly changes from 0 to 1 / radius of curvature, but the steering means has dynamics (dynamic characteristics), which inevitably causes a response delay. Therefore, the timing of turning will be delayed. If the control amount is increased in order to suppress such a response delay, the actual steering angle may be significantly overshooted. It should be noted that such a problem may also occur when switching from a reverse turn to a forward straight or a reverse straight.
そこで、上記作業車両の自律走行システムでは、上記演算手段は、後進直進から後進旋回へ切り替える場合、または、後進旋回から直進へ切り替える場合には、上記走行経路における切替え予定位置に到達する前に、切替え後の制御量を出力するように構成されていることが好ましい。 Therefore, in the autonomous traveling system of the work vehicle, when the calculation means is switched from the reverse straight direction to the reverse turn, or when the reverse turn is switched to the straight direction, the calculation means is used before reaching the planned switching position in the travel path. It is preferable that the control amount after switching is output.
この構成によれば、走行経路における切替え予定位置に到達する前に、先を見越して演算手段が早めに切替え後の制御量を出力することで、オーバーシュートの発生を抑えつつ、操舵手段のダイナミクスに起因する遅れを抑えることができる。 According to this configuration, the calculation means proactively outputs the control amount after the switching before reaching the planned switching position on the traveling route, thereby suppressing the occurrence of overshoot and the dynamics of the steering means. It is possible to suppress the delay caused by.
以上説明したように、本発明に係る作業車両の自律走行システムによれば、前輪が操舵輪である作業車両を、自律走行で後進旋回させる場合に、作業車両が蛇行するのを抑えることができる。 As described above, according to the autonomous traveling system of the work vehicle according to the present invention, it is possible to prevent the work vehicle from meandering when the work vehicle whose front wheel is the steering wheel is turned backward by autonomous travel. ..
以下、本発明を実施するための形態を図面に基づいて説明する。 Hereinafter, embodiments for carrying out the present invention will be described with reference to the drawings.
−自律走行システム−
図1は、本実施形態に係る自律走行システム1の概要を模式的に示す図である。この自律走行システム1は、GNSS衛星65を利用してロボットトラクタ10の位置を把握しながら、無線通信端末50に予め設定された走行経路に沿って、圃場内や枕地等においてロボットトラクタ10を自律走行させるものである。なお、「自律走行」とは、少なくとも操舵を自律的に行って、オペレータ(図示せず)の運転操作なしでロボットトラクタ10を走行させることを意味する。
-Autonomous driving system-
FIG. 1 is a diagram schematically showing an outline of the
ロボットトラクタ(作業車両)10は、図1に示すように、車体部11と、左右1対の前輪13と、左右一対の後輪15と、作業機17と、を備えている。車体部11の前部には、エンジン19が配置されており、エンジン19で発生した動力がHMT(Hydro Mechanical Transmission)21(図2参照)を介して前輪13および後輪15に伝達されるようになっている。一方、車体部11の後部には、オペレータが搭乗するキャビン23が設けられている。キャビン23の内部には、操舵輪である前輪13をオペレータが操舵するための操舵ハンドル25が配置されている。また、キャビン23の上部には、GNSS衛星65からの搬送波位相を検出することが可能なGNSSアンテナ27や、無線通信端末50と無線通信を行うための通信アンテナ29や、後述する慣性計測装置45等が設けられている。作業機17は、車体部11の後端部に設けられたリンク機構11aに装着されている。
As shown in FIG. 1, the robot tractor (working vehicle) 10 includes a
このロボットトラクタ10は、オペレータがキャビン23内に搭乗して操舵ハンドル25等を実際に操作することで走行する他、オペレータによって設定された無線通信端末50からの指示等に基づいて自律的に走行することが可能に構成されている。なお、ロボットトラクタ10は、無線通信端末50を操作するオペレータの乗車の有無に関係なく、換言すると、オペレータがキャビン23内に搭乗していてもしていなくても、自律走行を行うことが可能となっている。以下、このような自律走行を可能とする本実施形態の自律走行システム1について詳細に説明する。
The
図2は、自律走行システム1を模式的に示すブロック図である。本実施形態の自律走行システム1は、図2に示すように、ロボットトラクタ10に搭載された走行制御装置30と、オペレータによって操作される無線通信端末50と、基準局60と、を中心に構成されている。
FIG. 2 is a block diagram schematically showing the
〈走行制御装置〉
ロボットトラクタ10の走行制御装置30は、公知のコンピュータからなり、図示は省略するが、CPU(Central Processing Unit)、ROM(Read Only Memory)、RAM(Random Access Memory)、入力インターフェース、出力インターフェース等を含んでいる。CPUは、ROMに記憶された各種制御プログラムやマップに基づいて演算処理を実行する。ROMには、各種制御プログラムや、それら各種制御プログラムを実行する際に参照されるマップ等が記憶されている。RAMは、CPUによる演算結果や各センサの検出結果等を一時的に記憶するメモリである。この走行制御装置30には、図2に示すように、車速センサ41と、舵角センサ43と、慣性計測装置45と、GNSSアンテナ27と、通信アンテナ29と、エンジンECU31と、HMTECU33と、ステアリングECU35と、が接続されている。
<Driving control device>
The
車速センサ41は、例えば前輪13の車軸に配置されていて、前輪13の車軸の回転に応じたパルス信号を発生するように構成されている。車速センサ41による検出結果(パルス信号)は、走行制御装置30に入力される。
The
舵角センサ43は、ステアリングシャフト(図示せず)に配置されていて、ステアリング角度(ステアリングシャフトの回転角度)を検出するように構成されている。なお、舵角センサ43は、タイヤ角度(ロボットトラクタ10の前後方向軸に対する前輪13の車輪の傾斜角)を検出するように構成されていてもよい。舵角センサ43による検出結果は、走行制御装置30に入力される。
The
慣性計測装置45は、3つのジャイロセンサ(角速度センサ)と3つの加速度センサを備えている。これにより、慣性計測装置45は、上下方向を回転軸としたロボットトラクタ10の角速度であるヨーレートの他、ロールや、ピッチを検出することが可能となっている。慣性計測装置45が検出するロボットトラクタ10の角速度および加速度は、走行制御装置30に入力される。
The
GNSSアンテナ27は、GNSS衛星65からの搬送波位相を検出して、検出結果を走行制御装置30へ出力する。GNSSアンテナ27は、ロボットトラクタ10における車幅方向の中央部に設けられている。走行制御装置30は、後述するように、基準局60からの補正情報を受信した上で、RTKを利用して測位を行うことで、ロボットトラクタ10の位置を例えば緯度および経度の情報として算出する。
The
通信アンテナ29は、無線通信端末50の通信アンテナ53や、基準局60の無線通信装置63との間で各種の信号を送受信可能に構成されている。これにより、走行制御装置30は、適宜の方式で変調処理または復調処理を行うことで、無線通信端末50や基準局60との間で情報の遣り取りを行うことが可能となっている。
The
エンジンECU31は走行制御装置30の指令に基づいてエンジン19の制御を行うとともに、エンジン19の状態に関する情報を走行制御装置30へ出力するように構成されている。また、HMTECU33は走行制御装置30の指令に基づいてHMT21の制御を行うとともに、HMT21の状態に関する情報を走行制御装置30へ出力するように構成されている。走行制御装置30は、車速センサ41の検出結果により得られた、ロボットトラクタ10の現在の車速Vが目標の車速に近づくように、エンジン19の回転速度および/またはHMT21の変速比を変更することで、ロボットトラクタ10の車速Vを自律的に変更する車速制御を行うようになっている。
The
ステアリングECU35は、走行制御装置30との間で指令や情報の遣り取りを行うことが可能になっており、走行制御装置30の指令に基づいて、前輪13の電動パワーステアリング37の制御を行うように構成されている。走行制御装置30は、舵角センサ43の検出結果により得られた現在の舵角に基づいて、例えば電動パワーステアリング37の操舵トルクを制御するアシストモータ(図示せず)を駆動させることで、ロボットトラクタ10の舵角を自律的に変更する操舵制御を行うことが可能となっている。
The
〈無線通信端末〉
無線通信端末50は、タブレット型のコンピュータであり、図2に示すように、制御部51と、通信アンテナ53と、ディスプレイ55と、操作部57と、を備えている。
<Wireless communication terminal>
The
制御部51は、公知のコンピュータとして構成されており、図示は省略するが、CPU、ROM、RAM、入力インターフェース、出力インターフェース等を含んでいる。そして、制御部51は、走行経路生成部51aと、記憶部51bと、を備えている。
The
走行経路生成部51aは、オペレータの操作に基づいて、ロボットトラクタ10が走行する目標となる走行経路を生成する処理を行う。走行経路には、直線状走行経路71(図3参照)、前進旋回用走行経路73(図5参照)および後進旋回用走行経路75(図6参照)、並びに、これらを組み合わせた経路等がある。記憶部51bには、ロボットトラクタ10の設定に関する情報、オペレータが作成した走行経路71,73,75、および、走行経路71,73,75を生成するための各種条件等が記憶されている。
The travel
なお、請求項との関係では、本実施形態の無線通信端末50の走行経路生成部51aが、本発明でいうところの「作業車両が走行する目標となる走行経路を生成する経路生成手段」に相当する。
In relation to the claims, the travel
通信アンテナ53は、ロボットトラクタ10と無線通信を行うための近距離通信用のアンテナで構成されている。制御部51は、適宜の方式で変調処理または復調処理を行うことで、ロボットトラクタ10の走行制御装置30や基準局60との間で情報を遣り取りすることが可能となっている。
The
ディスプレイ55は、液晶ディスプレイや有機ELディスプレイ等であり、画像を表示可能に構成されている。ディスプレイ55は、制御部51の制御によって、例えば、自律走行に関する情報、ロボットトラクタ10の設定に関する情報、各種センサの検出結果および警告情報等を表示することが可能となっている。操作部57は、例えば、ディスプレイ55に重ねて配置されたタッチパネルであり、ディスプレイ55上におけるオペレータの指等による操作を検出するように構成されている。
The
〈基準局〉
本実施形態では、衛星測位システムを用いた測位方法としてRTK(Real Time Kinematic)を採用しており、基準点からの補正情報と、側位点の衛星測位情報と、を用いて、ロボットトラクタ10の現在位置を求めるようにしている。
<Standards Bureau>
In this embodiment, RTK (Real Time Kinematic) is adopted as a positioning method using a satellite positioning system, and the
具体的には、ロボットトラクタ10の走行の邪魔にならない位置に基準点となる基準局60を設置している。なお、基準局60を設置した箇所は、ロボットトラクタ10の走行制御装置30や無線通信端末50によって記憶される。基準局60は、GNSSアンテナ61と、ロボットトラクタ10の通信アンテナ29との間で各種の信号を送受信可能な無線通信装置63と、を備えている。これにより、基準点でも、GNSS衛星65からの搬送波位相を検出することができるとともに、基準局60とロボットトラクタ10との間で各種の情報が送受信可能になっている。
Specifically, the
そうして、基準点に設置された基準局60のGNSSアンテナ61と、側位点となるロボットトラクタ10のGNSSアンテナ27との両方でGNSS衛星65からの搬送波位相を検出する。基準局60では、GNSS衛星65から搬送波位相を検出する毎に、検出した搬送波位相と基準点の位置情報等を含む補正情報を生成して、無線通信装置63からロボットトラクタ10の通信アンテナ29に補正情報を送信する。これに対し、ロボットトラクタ10の走行制御装置30では、GNSSアンテナ27にて検出した搬送波位相と、基準局60から送信された補正情報とを用いて、上述の如く、ロボットトラクタ10の現在位置情報を緯度および経度の形式で算出する。
Then, the carrier phase from the
−自律操舵制御(前進時)−
〈直線状走行経路〉
次に、図3および図4を参照して、ロボットトラクタ10を直線状走行経路71に沿って自律走行させる場合の操舵制御について説明する。
-Autonomous steering control (when moving forward)-
<Linear travel route>
Next, with reference to FIGS. 3 and 4, steering control when the
図3は、直線状走行経路71に沿ってロボットトラクタ10を自律走行させる手法を模式的に説明する図であり、図4は、ヨーレートから演算タイヤ角度βを算出する手法を模式的に説明する図である。図3に示す直線状走行経路71から離れた位置を走行しているロボットトラクタ10を、当該直線状走行経路71に沿って自律走行させる場合、走行制御装置30は、以下に説明するように、前方注視方式にて目標ステアリング角速度を算出するように構成されている。この場合、操舵制御の入力は、前方目標角度φおよび演算タイヤ角度βである。
FIG. 3 is a diagram schematically explaining a method of autonomously traveling the
前方目標角度φは、ロボットトラクタ10を直線状走行経路71に沿って走行させる際に用いられる角度であり、前方目標点P1に基づいて、走行制御装置30によって算出される。ロボットトラクタ10には予め前方目標距離が設定されており、走行制御装置30は、ロボットトラクタ10のGNSSアンテナ27から直線状走行経路71までの距離が前方目標距離となる点(前方目標点P1)を特定する。また、走行制御装置30は、ロボットトラクタ10のGNSSアンテナ27の位置と前方目標点P1とを接続する線分(線L1)を引き、線L1と直線状走行経路71とがなす角を前方目標角度φとして算出する。
The forward target angle φ is an angle used when the
また、走行制御装置30は、慣性計測装置45が検出したヨーレートに基づいて、演算タイヤ角度βを算出する。演算タイヤ角度βは、検出したヨーレートでロボットトラクタ10を旋回させるための理想的な(理論上の)タイヤ角度、換言すれば、検出されたヨーレートをタイヤ角度に変換した値である。演算タイヤ角度βは、4輪車であるロボットトラクタ10を、図4に示すように、前輪13および後輪15が左右一対ではなく中央に1つずつあると仮定した2輪車モデルに変換して算出される。
Further, the
図4において、ωはヨーレートであり、Vは車速であり、Lはホイールベースであり、Rは2輪車モデルにおける旋回半径である。なお、車速Vは車速センサ41によって検出され、また、ホイールベースLは予め走行制御装置30に記憶されている。そうして、2輪車モデルにおける旋回半径Rは、下記の(式1)で表される。
In FIG. 4, ω is the yaw rate, V is the vehicle speed, L is the wheelbase, and R is the turning radius in the two-wheeled vehicle model. The vehicle speed V is detected by the
また、ヨーレート(角速度)ωと車速Vと旋回半径Rの関係は、(式1)を用いて下記の(式2)で表すことができる。 Further, the relationship between the yaw rate (angular velocity) ω, the vehicle speed V, and the turning radius R can be expressed by the following (Equation 2) using (Equation 1).
さらに、演算タイヤ角度βの値は、(式2)を変形して、下記の(式3)表すことができる。 Further, the value of the calculated tire angle β can be expressed by the following (Equation 3) by modifying (Equation 2).
このような処理を行うことで、慣性計測装置45が検出したヨーレートωから演算タイヤ角度βを算出することが可能となる。
By performing such processing, it becomes possible to calculate the calculated tire angle β from the yaw rate ω detected by the
そうして、走行制御装置30は、前方目標角度φから演算タイヤ角度βを減算する。これにより、タイヤ角度と進行方向との差異が含まれないように、換言すると、横滑りの影響がなくなるように、前方目標角度φが補正されることになる。
Then, the
走行制御装置30は、前方目標角度φから演算タイヤ角度βを減算した処理結果に基づいて目標ステアリング角速度を算出し、ステアリングECU35に出力する。ステアリングECU35が、目標ステアリング角速度に基づいて、前輪13の電動パワーステアリング37の制御を行うことで、ロボットトラクタ10の舵角が自律的に変更される。これにより、前方目標角度φおよび演算タイヤ角度βの値が変化するため、前方目標角度φから演算タイヤ角度βを減算した値も変化し、新たな目標ステアリング角速度が出力される。以上の処理を継続して行うことで、ロボットトラクタ10を直線状走行経路71に沿って走行させることができる。
The
〈前進旋回用走行経路〉
次に、図5を参照して、ロボットトラクタ10を前進旋回用走行経路73に沿って自律走行させる場合の操舵制御について説明する。
<Traveling route for forward turning>
Next, with reference to FIG. 5, steering control when the
図5は、前進旋回用走行経路73に沿ってロボットトラクタ10を自律走行させる手法を模式的に説明する図である。前進旋回用走行経路73に沿ってロボットトラクタ10を自律走行させる場合には、走行制御装置30は、直線状走行経路71の場合と同様に、前方注視方式にて目標ステアリング角速度を算出するように構成されている。この場合、操舵制御の入力は、前方目標角度φ、演算タイヤ角度βおよび旋回演算タイヤ角度である。
FIG. 5 is a diagram schematically illustrating a method of autonomously traveling the
前方目標角度φは、ロボットトラクタ10を前進旋回用走行経路73に沿って走行させる際に用いられる角度であり、前方目標点P1に基づいて、走行制御装置30によって算出される。走行制御装置30は、ロボットトラクタ10から前進旋回用走行経路73までの距離が前方目標距離となる点(前方目標点P1)を特定する。また、走行制御装置30は、ロボットトラクタ10のGNSSアンテナ27の位置と前方目標点P1とを接続する線分(線L1)を引く。
The forward target angle φ is an angle used when the
次に、走行制御装置30は、ロボットトラクタ10のGNSSアンテナ27の位置と、前進旋回用走行経路73の旋回中心P2と、を結ぶ線分である線L2を引く。そうして、走行制御装置30は、線L2に垂直であってGNSSアンテナ27の位置を通る線L3を引き、線L1と線L3がなす角を前方目標角度φとして算出する。ここで、線L3は、線L2と前進旋回用走行経路73との交点における前進旋回用走行経路73の接線であるため、現在のロボットトラクタ10位置における前進旋回用走行経路73の方向を示している。
Next, the traveling
また、走行制御装置30は、直線状走行経路71の場合と同様に、慣性計測装置45が検出したヨーレートωに基づいて、演算タイヤ角度βを算出する。そうして、走行制御装置30は、前方目標角度φから演算タイヤ角度βを減算する。
Further, the
さらに、走行制御装置30は、前進旋回用走行経路73に沿って旋回するためのヨーレートである旋回ヨーレートを、前進旋回用走行経路73の曲率半径等に基づいて推測し、推測した旋回ヨーレートを図4に示す手法で旋回演算タイヤ角度に変換する。そうして、走行制御装置30は、減算後の前方目標角度φに対して、旋回演算タイヤ角度を加算する。後の処理は、直線状走行経路71の場合と同様なので、説明を省略する。
Further, the traveling
−自律操舵制御(後進旋回時)−
次に、ロボットトラクタ10を後進旋回用走行経路75に沿って自律走行させる場合の操舵制御について説明する。
-Autonomous steering control (when turning backward)-
Next, steering control when the
図6は、枕地5において、ロボットトラクタ10を進行方向と直角に方向転換させる態様の一例を模式的に説明する図である。本実施形態では、圃場3の外側に設定された枕地5において、ロボットトラクタ10を進行方向と直角に方向転換させる態様として、図6に示すように、αターンを採用している。
FIG. 6 is a diagram schematically illustrating an example of a mode in which the
ここで、「αターン」とは、図6の白抜き矢印で示すように、直線状走行経路71aに沿って前進直進してきたロボットトラクタ10を、図6の黒塗り矢印で示すように、後進旋回用走行経路75に沿って後進旋回させて直角に方向転換させた後、図6のハッチング矢印で示すように、方向転換後のロボットトラクタ10を直線状走行経路71bに沿って前進直進させる方向転換態様を指す。
Here, the “α-turn” means that the
なお、後進旋回用走行経路75には、図6に示すように、円弧経路75bの両側に短い直線経路75aおよび直線経路75cが含まれている。それ故、より正確には、αターンとは、直線状走行経路71aに沿って前進直進してきたロボットトラクタ10を、一旦直線経路75aに沿って短く後進させてから、円弧経路75bに沿って後進旋回させて直角に方向転換させた後、直線経路75cに沿って短く後進させてから、方向転換後のロボットトラクタ10を直線状走行経路71bに沿って前進直進させる方向転換態様を指す。
As shown in FIG. 6, the reverse turning
ここで、ロボットトラクタ10を、後進旋回用走行経路75に沿って後進旋回させる場合にも、前進の場合と同様に、前方注視方式を採用することが考えられる。しかしながら、前輪13が操舵輪であるロボットトラクタ10の後進旋回時に、前方注視方式を採用した場合には、以下のような問題を生ずるおそれがある。
Here, when the
すなわち、旋回の向きと前輪13の操舵の向きとが一致する前進旋回とは異なり、後進旋回の場合には、旋回の向きとは反対の向きに前輪13を操舵する必要がある。それ故、後進旋回を開始すると、前輪13側が後進旋回用走行経路75の外側に膨らむとともに、GNSSアンテナ27の位置が後進旋回用走行経路75の内側に入ることになるため、旋回の向きと同じ向きに前輪13を操舵する補正が入ることになる。すると、今度は、GNSSアンテナ27の位置が後進旋回用走行経路75の外側に出ることになるため、旋回の向きとは反対の向きに前輪13を操舵する補正が入ることになる。このような操舵制御とフィードバック補正とが繰り返されることで、操舵制御が発散し、ロボットトラクタ10が蛇行するという問題が生じることになる。
That is, unlike forward turning in which the turning direction and the steering direction of the
そこで、本実施形態の自律走行システム1では、予測モデルを用いたモデル予測制御により、所定の制御周期毎に最適な制御量を算出するようにしている。
Therefore, in the
具体的には、後進旋回時には、後進旋回用走行経路75からの制御中心CCの偏差と操舵角の制御量との相関を状態方程式で表した予測モデルに基づき、所定の制御周期毎に、所定の評価関数Jを用いて最適な制御量を所定のステップ先まで算出し、第1番目のステップの制御量のみを出力するとともに、制御量に対応する制御中心CCの予測位置と、GNSSアンテナ27等によって検出される制御中心CCの実測位置とに基づいて、予測モデルを補正するように、走行制御装置30を構成している。以下、ロボットトラクタ10を後進旋回用走行経路75に沿って自律走行させる場合のモデル予測制御について詳細に説明する。
Specifically, at the time of reverse turning, the correlation between the deviation of the control center CC from the reverse turning
〈状態方程式〉
本実施形態では、状態方程式(差分方程式)を下記の(式4)のように定義する。
<Equation of state>
In this embodiment, the equation of state (difference equation) is defined as shown in (Equation 4) below.
ここで、k(0以上の整数)は予測ステップ数であり、Aは状態係数ベクトルであり、x[k]は第kステップにおけるモデルの状態ベクトルであり、b1は入力係数ベクトルであり、u[k]は第kステップにおける入力(制御量)であり、b2は外乱係数ベクトルであり、d[k]は第kステップにおける既知外乱である。 Here, k (an integer of 0 or more) is the number of predicted steps, A is the state coefficient vector, x [k] is the state vector of the model in the kth step, and b 1 is the input coefficient vector. u [k] is the input (control amount) in the k-th step, b 2 is the disturbance coefficient vector, and d [k] is the known disturbance in the k-th step.
そうして、上記(式4)から、以下のような状態時系列を考える。 Then, from the above (Equation 4), consider the following state time series.
以上の状態系列をまとめて、行列で表現すると、下記の(式5)のように表すことができる。なお、(式5)におけるEは単位行列である。 When the above state series are collectively expressed by a matrix, it can be expressed as (Equation 5) below. Note that E in (Equation 5) is an identity matrix.
ここで、状態系列ベクトルX、入力系列ベクトルU、外乱系列ベクトルD、行列H1、行列H2および行列H3を、下記の(式6)〜(式11)のように定義する。 Here, the state sequence vector X, the input sequence vector U, the disturbance sequence vector D, the matrix H 1 , the matrix H 2, and the matrix H 3 are defined as the following (Equation 6) to (Equation 11).
すると、上記(式5)は、下記の(式12)のように表すことができ、この(式12)で表される状態方程式が、下記の評価関数Jによる評価対象となる。 Then, the above (Equation 5) can be expressed as the following (Equation 12), and the equation of state represented by this (Equation 12) becomes an evaluation target by the following evaluation function J.
〈評価関数〉
本実施形態では、評価関数Jを下記の(式13)のように定義する。なお、W1およびW2は重み行列である。
<Evaluation function>
In this embodiment, the evaluation function J is defined as shown in (Equation 13) below. Note that W 1 and W 2 are weight matrices.
上記(式13)における1/2は、後の計算に出てくる2を打ち消すために与えられた単なる係数である。それ故、上記(式13)の評価関数Jは、状態系列ベクトルXの2乗と入力系列ベクトルUの2乗との和を重み付けして評価しているといえる。 1/2 in the above (Equation 13) is just a coefficient given to cancel the 2 that appears in the later calculation. Therefore, it can be said that the evaluation function J of the above (Equation 13) weights and evaluates the sum of the square of the state series vector X and the square of the input series vector U.
上記(式13)に上記(式12)の状態方程式を代入して整理すると、評価関数Jを下記の(式14)のように表すことができる。なお、(式12)を整理する際に出てくる、UH2 TW1H1x0、UTH2 TW1H3DおよびDTH1 TW1H1x0は、最終的にはスカラーになることから、それぞれx0 TH1 TW1H2U、DTH3 TW1H2Uおよびx0 TH1 TW1H1Dに変換可能である。
By substituting the equation of state of the above (Equation 12) into the above (Equation 13) and rearranging it, the evaluation function J can be expressed as the following (Equation 14). In addition, UH 2 T W 1 H 1 x 0 , U T H 2 T W 1 H 3 D and D T H 1 T W 1 H 1 x 0 , which appear when (Equation 12) is arranged, are the final. thereof include from becoming scalar can be converted into x 0 T H 1 T W 1 H 2 U,
上記(式14)に示す評価関数Jを、入力系列ベクトルUで微分すると、下記の(式15)が得られる。 When the evaluation function J shown in the above (Equation 14) is differentiated by the input series vector U, the following (Equation 15) is obtained.
ここで、∂(J(U))/∂U=0となる入力系列ベクトルUで、評価関数Jが最小となることから、∂(J(U))/∂U=0を解くことで、下記の(式16)で示す、最適な入力系列ベクトルUが得られる。 Here, since the evaluation function J is the minimum in the input sequence vector U such that ∂ (J (U)) / ∂U = 0, by solving ∂ (J (U)) / ∂U = 0, The optimum input sequence vector U shown by the following (Equation 16) can be obtained.
そうして、走行制御装置30は、例えば、予測サンプリング時間dt=0.1(秒)毎に、上記(式16)に基づき最適な入力(制御量)を5ステップ先まで算出し、第1番目のステップの制御量のみを出力する。それ故、請求項との関係では、予測サンプリング時間dtが、本発明でいうところの「所定の制御周期」に相当し、5ステップが、本発明でいうところの「所定のステップ」に相当する。
Then, for example, the
〈プラントモデル〉
図7は、時計回りの後進旋回時のプラントモデルを模式的に説明する図であり、図8は、反時計回りの後進旋回時のプラントモデルを模式的に説明する図である。本実施形態では、図7および図8に示すように、後進旋回時のプラントモデルとして、前進時の操舵制御と同様に、2輪車モデルを採用している。
<Plant model>
FIG. 7 is a diagram schematically explaining a plant model in a clockwise reverse turn, and FIG. 8 is a diagram schematically explaining a plant model in a counterclockwise reverse turn. In this embodiment, as shown in FIGS. 7 and 8, a two-wheeled vehicle model is adopted as a plant model at the time of reverse turning, as in the case of steering control at the time of forward movement.
もっとも、後進旋回時のプラントモデルでは、GNSSアンテナ27の位置を基準として制御を行う前進時の操舵制御と異なり、ロボットトラクタ10の制御中心CCを、後輪15の車軸における車幅方向の中心位置に設定している。GNSSアンテナ27と後輪15の車軸の中心位置との位置関係は予め走行制御装置30に記憶されている。走行制御装置30は、上述の如く、GNSSアンテナ27、基準局60およびGNSS衛星65等を用いて、GNSSアンテナ27の現在位置を算出し、かかるGNSSアンテナ27の現在位置に基づいて、制御中心CCの現在位置を検出するように構成されている。
However, in the plant model during reverse turning, unlike the steering control during forward control, which is controlled based on the position of the
このように、制御中心CCをロボットトラクタ10の後輪15の車軸における車幅方向の中心位置に設定することから、換言すると、GNSSアンテナ27の仮想位置(制御中心CC)を、後進旋回時の車両中心に設定することから、蛇行の影響をより一層抑えることができる。
In this way, since the control center CC is set at the center position in the vehicle width direction on the axle of the
なお、請求項との関係では、走行制御装置30、GNSSアンテナ27、基準局60およびGNSS衛星65等が、本発明でいうところの「作業車両に設定された制御中心を検出可能な位置検出手段」に相当する。なお、以下では、これらをまとめて「位置検出システム」とも称する。
In relation to the claims, the traveling
そうして、本実施形態では、上記(式4)の状態方程式の状態ベクトルxを、下記の(式17)のように定義している。 Then, in this embodiment, the state vector x of the state equation of the above (Equation 4) is defined as the following (Equation 17).
なお、図7および図8に示すように、yは、後進旋回用走行経路75からの制御中心CCの横方向偏差であり、θは、後進旋回用走行経路75の接線75’に対する後輪15のタイヤ角度偏差であり、δは、前輪13のタイヤ角度である。
As shown in FIGS. 7 and 8, y is the lateral deviation of the control center CC from the reverse turning
ここで、時計回りの後進旋回時における、第kステップでの横方向偏差y[k]と第k+1ステップでの横方向偏差y[k+1]との関係は、図7を参照して、車速V、第kステップでのタイヤ角度偏差θ[k]、および、予測サンプリング時間dtを用いて、下記の(式18)のように表すことができる。 Here, the relationship between the lateral deviation y [k] in the kth step and the lateral deviation y [k + 1] in the k + 1th step during the clockwise reverse turn is the vehicle speed V with reference to FIG. 7. , The tire angle deviation θ [k] in the kth step, and the predicted sampling time dt can be expressed as the following (Equation 18).
また、時計回りの後進旋回時における、第kステップでのタイヤ角度偏差θ[k]と第k+1ステップでのタイヤ角度偏差θ[k+1]との関係は、図7を参照して、車速V、ホイールベースL、第kステップでの前輪13のタイヤ角度δ[k]、第kステップでの曲率c[k]、および、予測サンプリング時間dtを用いて、下記の(式19)のように表すことができる。
Further, the relationship between the tire angle deviation θ [k] in the kth step and the tire angle deviation θ [k + 1] in the k + 1th step during the clockwise reverse turn is as follows. Using the wheel base L, the tire angle δ [k] of the
さらに、ステアリングダイナミクスの時定数をτと規定し、第kステップでの前輪13のタイヤ角度の制御量(ステアリング角度)をu[k]とすると、第kステップでのタイヤ角度δ[k]と第k+1ステップでのタイヤ角度δ[k+1]との関係は、図7を参照して、時定数τ、第kステップでの制御量u[k]、および、予測サンプリング時間dtを用いて、下記の(式20)のように表すことができる。
Further, assuming that the time constant of steering dynamics is defined as τ and the control amount (steering angle) of the tire angle of the
以上の(式18)〜(式20)をまとめて、時計回りの後進旋回時におけるプラントモデルの運動方程式は、下記の(式21)のように定義される。 Summarizing the above (Equation 18) to (Equation 20), the equation of motion of the plant model at the time of clockwise reverse turning is defined as the following (Equation 21).
一方、反時計回りの後進旋回時においても、第kステップでの横方向偏差y[k]と第k+1ステップでの横方向偏差y[k+1]との関係、および、第kステップでのタイヤ角度δ[k]と第k+1ステップでのタイヤ角度δ[k+1]との関係は、上記(式18)および(式20)と同じになる。 On the other hand, even during a counterclockwise reverse turn, the relationship between the lateral deviation y [k] in the kth step and the lateral deviation y [k + 1] in the k + 1 step, and the tire angle in the kth step. The relationship between δ [k] and the tire angle δ [k + 1] in the k + 1th step is the same as in (Equation 18) and (Equation 20) above.
これに対し、反時計回りの後進旋回時における、第kステップでのタイヤ角度偏差θ[k]と第k+1ステップでのタイヤ角度偏差θ[k+1]との関係は、図8を参照して、下記の(式22)のように表すことができる。 On the other hand, the relationship between the tire angle deviation θ [k] in the kth step and the tire angle deviation θ [k + 1] in the k + 1th step during the counterclockwise reverse turn is described with reference to FIG. It can be expressed as (Equation 22) below.
よって、(式18)、(式22)および(式20)をまとめて、反時計回りの後進旋回時におけるプラントモデルの運動方程式は、下記の(式23)のように定義される。 Therefore, (Equation 18), (Equation 22) and (Equation 20) are put together, and the equation of motion of the plant model at the time of reverse turning counterclockwise is defined as the following (Equation 23).
これら(式21)および(式23)で表されるプラントモデルに上記(式4)の状態方程式を適用する場合には、上記(式4)の状態方程式と上記(式21)のプラントモデルの運動方程式とを見比べれば分かるように、状態方程式における状態係数ベクトルA、入力係数ベクトルb1および外乱係数ベクトルb2は、それぞれ下記の(式24)、(式25)および(式26)で規定される。 When the state equation of the above (Equation 4) is applied to the plant models represented by these (Equation 21) and (Equation 23), the state equation of the above (Equation 4) and the plant model of the above (Equation 21) are used. As can be seen by comparing with the equation of motion, the state coefficient vector A, the input coefficient vector b 1 and the disturbance coefficient vector b 2 in the state equation are expressed in the following (Equation 24), (Equation 25) and (Equation 26), respectively. It is stipulated.
これら状態係数ベクトルA、入力係数ベクトルb1および外乱係数ベクトルb2を、上記(式9)〜(式11)に代入することで、これらを要素とする行列H1、行列H2および行列H3を算出することができる。そうして、これら行列H1、行列H2および行列H3を上記(式12)に代入したものが、本実施形態の予測モデルとなる。 By substituting the state coefficient vector A, the input coefficient vector b 1 and the disturbance coefficient vector b 2 into the above (Equation 9) to (Equation 11), the matrix H 1 , the matrix H 2 and the matrix H having these as elements are substituted. 3 can be calculated. Then, the matrix H 1 , the matrix H 2, and the matrix H 3 are substituted into the above (Equation 12), which is the prediction model of the present embodiment.
また、重み行列W1および重み行列W2は、適宜設定されることから、上記(式16)に基づいて、最適な制御量を算出することが可能となる。 Further, since the weight matrix W 1 and the weight matrix W 2 are appropriately set, it is possible to calculate the optimum control amount based on the above (Equation 16).
〈モデル予測制御〉
次に、時計回りの後進旋回を例として、上記予測モデルおよび評価関数Jを用いたモデル予測制御について説明する。なお、図9は、モデル予測制御の概念を模式的に説明する図である。また、各制御周期は、予測サンプリング時間dt=0.1である。
<Model prediction control>
Next, the model prediction control using the prediction model and the evaluation function J will be described by taking a clockwise reverse turn as an example. Note that FIG. 9 is a diagram schematically explaining the concept of model prediction control. Further, each control cycle has a predicted sampling time dt = 0.1.
≪第1制御周期≫
先ず、本実施形態のモデル予測制御は、図6に示すように、直線状走行経路71aに沿って前進直進してきたロボットトラクタ10の制御中心CCが、圃場3のコーナー近傍の枕地5に設定された制御開始位置SPに到達したことが、位置検出システムによって検出されたときに開始される。これにより、ロボットトラクタ10が直線経路75aに沿って後進し始める。
<< 1st control cycle >>
First, in the model prediction control of the present embodiment, as shown in FIG. 6, the control center CC of the
そうして、本実施形態のモデル予測制御では、図9の(A)に示すように、位置検出システム等によって検出されたプラントモデルの毎ステップの観測値x0を、上記(式12)の予測モデルに代入して、上記(式16)から最適な入力系列を算出する。 Then, in the model prediction control of the present embodiment, as shown in FIG. 9A, the observed value x 0 of each step of the plant model detected by the position detection system or the like is set to the above (Equation 12). Substituting into the prediction model, the optimum input sequence is calculated from the above (Equation 16).
具体的には、後進旋回開始時(k=0)における、制御中心CCの横方向偏差y0は、位置検出システムによって検出されたロボットトラクタ10の現在位置情報から求めることができる。また、後進旋回開始時(k=0)における、タイヤ角度偏差θ0は、ロボットトラクタ10の位置変化情報から求めることができる。なお、タイヤ角度偏差θ0は、方位センサ(図示せず)で求めてもよい。さらに、後進旋回開始時(k=0)における、前輪13のタイヤ角度δ0は、舵角センサ43の検出結果から求めることができる。よって、プラントモデルの毎ステップの観測値であるx0=[y0 θ0 δ0]Tを算出することができる。
Specifically, the reverse turning start at (k = 0), lateral deviation y 0 of the control center CC may be determined from the current position information of the
また、外乱系列ベクトルD=[d[1] d[1] d[2]…d[N−1]]Tは、後進旋回用走行経路75が円弧であり、曲率半径(旋回半径R)が一定であれば、c[k]=1/Rなので一定となる。
Further, in the disturbance series vector D = [d [1] d [1] d [2] ... d [N-1]] T , the traveling
これらにより、上記(式16)から、5ステップ分の(0.5秒先までの)入力系列ベクトルU=[u[1] u[2] u[3] u[4] u[5]]Tを算出することができる。もっとも、走行制御装置30は、算出した第1〜第5ステップのうち第1番目のステップの制御量u[1](ステアリング角度(操舵角))のみをステアリングECU35へ出力し、ステアリングECU35は制御量u[1]に基づいて前輪13を制御する。それ故、請求項との関係では、走行制御装置30が、本発明でいうところの「前輪の操舵角の制御量を演算して出力する演算手段」に相当し、ステアリングECU35が、本発明でいうところの「演算手段によって出力された制御量に基づいて前輪を制御する操舵手段」に相当する。
As a result, from the above (Equation 16), the input sequence vector U = [u [1] u [2] u [3] u [4] u [5] for 5 steps (up to 0.5 seconds ahead)]. T can be calculated. However, the
なお、位置検出システムによって検出された制御中心CCの現在位置を、図9の(A)に示すような毎ステップの観測値として用いるのは、第1制御周期だけである。 It should be noted that the current position of the control center CC detected by the position detection system is used only in the first control cycle as an observation value at each step as shown in FIG. 9A.
≪第2制御周期以降≫
第2制御周期では、走行制御装置30は、制御量に対応する制御中心CCの予測位置と、位置検出システムによって検出された制御中心CCの実測位置とに基づいて、予測モデルを補正する。より詳しくは、走行制御装置30は、制御量u[1]に基づいて前輪13が操舵された後に、位置検出システムによって検出された制御中心CCの現在位置と、制御量u[1]に対応する制御中心CCの予測位置との誤差を、図9の(B)に示すように、外乱系列ベクトルD=[d[1] d[2] d[3] d[4] d[5]]Tとして予測モデルに反映させて、5ステップ分の入力系列ベクトルU=[u[2] u[3] u[4] u[5] u[6]]Tを算出する。そうして、走行制御装置30は、算出した第2〜第5ステップのうち第1番目のステップの制御量u[2]のみをステアリングECU35へ出力し、ステアリングECU35は制御量u[2]に基づいて前輪13を制御する。
≪After the second control cycle≫
In the second control cycle, the
同様に、第3制御周期では、外乱系列ベクトルD=[d[2] d[3] d[4] u[5] u[6]]Tとして予測モデルに反映させて、5ステップ分の入力系列ベクトルU=[u[3] u[4] u[5] u[6] u[7]]Tを算出し、制御量u[3]のみをステアリングECU35へ出力する。以後、このような制御を継続して行うことで、ロボットトラクタ10を後進旋回用走行経路75に沿って走行させる。
Similarly, in the third control cycle, the disturbance sequence vector D = [d [2] d [3] d [4] u [5] u [6]] T is reflected in the prediction model and input for 5 steps. The series vector U = [u [3] u [4] u [5] u [6] u [7]] T is calculated, and only the control amount u [3] is output to the
このように、本実施形態によれば、予測モデルに基づき0.1(秒)毎に、最適な制御量を5ステップ先まで算出することから、換言すると、5ステップ先を予測した上で最適な制御量を算出することから、制御中心CCの予測位置と実測位置とに基づいて、予測モデルをフィードバック補正する際に、操舵制御が発散するのを抑えることができる。つまり、ロボットトラクタ10の状態を0.5秒先まで予測した上で、最適な制御量を算出することから、制御中心CCの実測位置とのずれが予測の範囲内であれば、蛇行を生じるような制御量が出力されるのを抑えることができる。したがって、前輪13が操舵輪であるロボットトラクタ10を自律走行で後進旋回させる場合にも、ロボットトラクタ10が蛇行するのを抑えることができる。
As described above, according to the present embodiment, the optimum control amount is calculated up to 5 steps ahead every 0.1 (seconds) based on the prediction model. In other words, it is optimal after predicting 5 steps ahead. Since a large amount of control is calculated, it is possible to suppress divergence of steering control when feedback-correcting the prediction model based on the predicted position and the actually measured position of the control center CC. That is, since the optimum control amount is calculated after predicting the state of the
〈制御タイミング〉
図10は、後進直進から後進旋回に切り替える場合、および、後進旋回から後進直進に切り替える場合における、曲率cの変化を模式的に説明する図である。上述の如く、前進直進してきたロボットトラクタ10をαターンさせる場合には、ロボットトラクタ10を前進直進からいきなり後進旋回させる訳ではなく、直線経路75aに沿った後進直進を経てから、円弧経路75bに沿った後進旋回を開始する。このように後進直進から後進旋回へ切り替える場合には、図10に示すように、曲率cが0から1/Rへ急激に変化するが、電動パワーステアリング37にはダイナミクス(動的特性)があるため、不可避的に応答遅れが生じてしまい、旋回のタイミングが遅れることになる。
<Control timing>
FIG. 10 is a diagram schematically illustrating a change in curvature c when switching from reverse straight to reverse turn and when switching from reverse turn to reverse straight. As described above, when the
かかる応答遅れを抑えるために、制御量u[k]を大きくすると、実際の操舵角が大きくオーバーシュートしてしまうおそれがある。なお、このような問題は、円弧経路75bに沿った後進旋回から前進直進または直線経路75cに沿った後進直進へ切り替える場合にも生じ得る。
If the control amount u [k] is increased in order to suppress such a response delay, the actual steering angle may be significantly overshooted. It should be noted that such a problem may also occur when switching from a reverse turn along the
そこで、本実施形態の自律走行システム1では、後進直進から後進旋回へ切り替える場合、または、後進旋回から直進へ切り替える場合には、後進旋回用走行経路75における切替え予定位置に到達する前に、切替え後の制御量を出力するように走行制御装置30を構成している。以下、このような制御について詳細に説明する。
Therefore, in the
≪後進直進から後進旋回へ切り替える場合≫
先ず、直線経路75aの終点E1(切替え予定位置)では、曲率c=1/Rとなることが好ましい。また、仮にステップ状に曲率cを増加させても電動パワーステアリング37は追従することができないため、終点E1に到達するまでに曲率cを滑らかに増加させるのが好ましい。これらを考慮すると、下記の(式27)で表される、図11に示すようなグラフを描くことができる。
≪When switching from reverse straight to reverse turn≫
First, at the end point E1 (planned switching position) of the
なお、(式27)において、AEDは、直線経路75a上での終点E1(=0)までの距離であり、また、XP1は、直線経路75a上において曲率cの増加を開始する位置であり、どこから曲率cを増加させるかに応じて適宜設定される値である。
In (Equation 27), AED is the distance to the end point E1 (= 0) on the
ここで、速度の横方向成分は微小であると考えて、d(AED)/dt=車速Vとすると、(式27)における傾きの時間変化は、下記の(式28)で表される。 Here, assuming that the lateral component of the speed is minute and d (AED) / dt = vehicle speed V, the time change of the inclination in (Equation 27) is expressed by the following (Equation 28).
これを予測サンプリング時間dtで離散化(オイラー近似)すると、第k+1ステップでの曲率c[k+1]と、第kステップでの曲率c[k]との関係は、下記の(式29)で表される。 When this is discretized by the predicted sampling time dt (Euler approximation), the relationship between the curvature c [k + 1] in the k + 1 step and the curvature c [k] in the kth step is shown in the following (Equation 29). Will be done.
そうして、初期値c[0]を下記の(式30)のように定め、上記(式29)を用いて予測区間を求める。ただし、曲率c[k]は1/Rで飽和させる。 Then, the initial value c [0] is set as shown in (Equation 30) below, and the prediction interval is obtained using the above (Equation 29). However, the curvature c [k] is saturated at 1 / R.
≪後進旋回から直進へ切り替える場合≫
円弧経路75bに沿った後進旋回が終わるまで同じ曲率1/Rで旋回し続けると、円弧経路75bの終点E2から惰性で進むときにロボットトラクタ10が切れ込み過ぎてしまうおそれがある。
≪When switching from reverse turning to straight ahead≫
If the robot continues to turn with the
それ故、後進旋回から直進へ切り替える場合には、円弧経路75bの終点E2(切替え予定位置)では、曲率c=0となることが好ましい。また、仮にステップ状に曲率cを減少させても電動パワーステアリング37は追従することができないため、終点E2に到達するまでに曲率cを滑らかに減少させるのが好ましい。これらを考慮すると、下記の(式31)で表される、図12に示すようなグラフを描くことができる。
Therefore, when switching from reverse turning to straight traveling, it is preferable that the curvature c = 0 at the end point E2 (planned switching position) of the
なお、(式31)において、AEDは、円弧経路75b上での終点E2(=0)までの距離であり、また、XP2は、円弧経路75b上において曲率cの減少を開始する位置であり、どこから曲率cを減少させるかに応じて適宜設定される値である。
In (Equation 31), AED is the distance to the end point E2 (= 0) on the
ここで、d(AED)/dt=車速Vとすると、(式31)における傾きの時間変化は、下記の(式32)で表される。 Here, assuming that d (AED) / dt = vehicle speed V, the time change of the inclination in (Equation 31) is expressed by the following (Equation 32).
これを予測サンプリング時間dtで離散化すると、第k+1ステップでの曲率c[k+1]と、第kステップでの曲率c[k]との関係は、下記の(式33)で表される。 When this is discretized by the predicted sampling time dt, the relationship between the curvature c [k + 1] in the k + 1th step and the curvature c [k] in the kth step is expressed by the following (Equation 33).
そうして、初期値c[0]を下記の(式34)のように定め、上記(式33)を用いて予測区間を求める。ただし、曲率c[k]は0で飽和させる。 Then, the initial value c [0] is set as shown in (Equation 34) below, and the prediction interval is obtained using the above (Equation 33). However, the curvature c [k] is saturated at 0.
以上の操作を施すことで、図12に示すように、後進直進から後進旋回へ切り替える場合には、終点E1よりも手前で、また、後進旋回から直進へ切り替える場合には、終点E2よりも手前で、走行制御装置30が、先を見越して早めに切替え後の制御量を出力することで、オーバーシュートの発生を抑えつつ、操舵手段のダイナミクスに起因する遅れを抑えることができる。
By performing the above operation, as shown in FIG. 12, when switching from reverse straight to reverse turn, it is before the end point E1, and when switching from reverse turn to straight, it is before end point E2. Therefore, the
(その他の実施形態)
本発明は、実施形態に限定されず、その精神又は主要な特徴から逸脱することなく他の色々な形で実施することができる。
(Other embodiments)
The present invention is not limited to embodiments and can be practiced in various other forms without departing from its spirit or key features.
上記実施形態では、トラクタに本発明を適用したが、これに限らず、例えば、田植機、コンバイン、土木・建築作業装置、除雪車等に本発明を適用してもよい。 In the above embodiment, the present invention is applied to a tractor, but the present invention is not limited to this, and the present invention may be applied to, for example, a rice transplanter, a combine, a civil engineering / building work device, a snowplow, and the like.
また、上記実施形態では、衛星測位システムを用いた測位方法としてRTKを採用したが、これに限らず、例えば、DGPS(Differential GPS)や、基準点を用いない単独測位を採用してもよい。 Further, in the above embodiment, RTK is adopted as a positioning method using a satellite positioning system, but the present invention is not limited to this, and for example, DGPS (Differential GPS) or independent positioning without using a reference point may be adopted.
さらに、上記実施形態では、前進については前方注視方式を採用したが、これに限らず、前進についてもモデル予測制御を採用するようにしてもよい。 Further, in the above embodiment, the forward gaze method is adopted for the forward movement, but the present invention is not limited to this, and the model prediction control may be adopted for the forward movement.
このように、上述の実施形態はあらゆる点で単なる例示に過ぎず、限定的に解釈してはならない。さらに、特許請求の範囲の均等範囲に属する変形や変更は、全て本発明の範囲内のものである。 Thus, the above embodiments are merely exemplary in all respects and should not be construed in a limited way. Further, all modifications and modifications belonging to the equivalent scope of the claims are within the scope of the present invention.
本発明によると、前輪が操舵輪である作業車両を、自律走行で後進旋回させる場合に、作業車両が蛇行するのを抑えることができるので、作業車両の自律走行システムに適用して極めて有益である。 According to the present invention, when a work vehicle whose front wheel is a steering wheel is turned backward by autonomous driving, it is possible to suppress meandering of the work vehicle, which is extremely useful in application to an autonomous driving system of a work vehicle. be.
1 自律走行システム
10 ロボットトラクタ(作業車両)
13 前輪
15 後輪
27 GNSSアンテナ(位置検出手段)
30 走行制御装置(演算手段)(位置検出手段)
35 ステアリングECU(操舵手段)
51a 走行経路生成部(経路生成手段)
60 基準局(位置検出手段)
65 GNSS衛星(位置検出手段)
75 後進旋回用走行経路
CC 制御中心
E2 終点(切替え予定位置)
E1 終点(切替え予定位置)
u[k] ステアリング角度(操舵角)
1
13
30 Travel control device (calculation means) (position detection means)
35 Steering ECU (steering means)
51a Travel route generation unit (route generation means)
60 Reference station (position detection means)
65 GNSS satellite (position detection means)
75 Travel path for reverse turning CC Control center E2 End point (planned switching position)
E1 end point (planned switching position)
u [k] Steering angle (steering angle)
Claims (3)
上記作業車両が走行する目標となる走行経路を生成する経路生成手段と、
上記作業車両に設定された制御中心の位置を検出可能な位置検出手段と、
上記前輪の操舵角の制御量を演算して出力する演算手段と、
上記演算手段によって出力された制御量に基づいて前輪を制御する操舵手段と、を備え、
上記演算手段は、後進旋回時には、上記走行経路からの上記制御中心の偏差と上記操舵角の制御量との相関を状態方程式で表した予測モデルに基づき、所定の制御周期毎に、所定の評価関数を用いて最適な制御量を所定のステップ先まで算出し、第1番目のステップの制御量のみを出力するとともに、当該制御量に対応する上記制御中心の予測位置と、上記位置検出手段によって検出された当該制御中心の実測位置とに基づいて、当該予測モデルを補正するように構成されていることを特徴とする作業車両の自律走行システム。 It is an autonomous driving system of a work vehicle that autonomously travels a work vehicle whose front wheel is a steering wheel.
A route generation means for generating a target travel route on which the work vehicle travels, and
A position detecting means capable of detecting the position of the control center set on the work vehicle, and
A calculation means that calculates and outputs the control amount of the steering angle of the front wheels, and
A steering means for controlling the front wheels based on the control amount output by the above calculation means is provided.
The calculation means performs a predetermined evaluation for each predetermined control cycle based on a prediction model in which the correlation between the deviation of the control center from the travel path and the control amount of the steering angle is expressed by a state equation at the time of reverse turning. The optimum control amount is calculated up to a predetermined step destination using a function, and only the control amount of the first step is output, and the predicted position of the control center corresponding to the control amount and the position detection means are used. An autonomous travel system for a work vehicle, characterized in that it is configured to correct the prediction model based on the detected measured position of the control center.
上記制御中心は、上記作業車両の後輪の車軸における車幅方向の中心位置に設定されており、
上記位置検出手段は、上記作業車両に設けられたGNSSアンテナの現在位置に基づいて、上記制御中心の現在位置を検出するように構成されていることを特徴とする作業車両の自律走行システム。 In the autonomous traveling system of the work vehicle according to claim 1,
The control center is set at the center position in the vehicle width direction on the axle of the rear wheel of the work vehicle.
The position detecting means is an autonomous traveling system for a work vehicle, characterized in that it is configured to detect the current position of the control center based on the current position of the GNSS antenna provided on the work vehicle.
上記演算手段は、後進直進から後進旋回へ切り替える場合、または、後進旋回から直進へ切り替える場合には、上記走行経路における切替え予定位置に到達する前に、切替え後の制御量を出力するように構成されていることを特徴とする作業車両の自律走行システム。 In the autonomous traveling system of the work vehicle according to claim 1 or 2,
The calculation means is configured to output the control amount after switching before reaching the planned switching position in the travel route when switching from reverse straight to reverse turn or from reverse turn to straight. An autonomous driving system for work vehicles, characterized by being
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020096400A JP7386757B2 (en) | 2020-06-02 | 2020-06-02 | Autonomous driving system for work vehicles |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020096400A JP7386757B2 (en) | 2020-06-02 | 2020-06-02 | Autonomous driving system for work vehicles |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2021189909A true JP2021189909A (en) | 2021-12-13 |
JP7386757B2 JP7386757B2 (en) | 2023-11-27 |
Family
ID=78849553
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020096400A Active JP7386757B2 (en) | 2020-06-02 | 2020-06-02 | Autonomous driving system for work vehicles |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7386757B2 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117075485A (en) * | 2023-10-18 | 2023-11-17 | 舜泰汽车有限公司 | Agricultural machinery autopilot navigation equipment |
JP7388580B1 (en) | 2023-02-15 | 2023-11-29 | いすゞ自動車株式会社 | Information output device and information output method |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5174636B2 (en) | 2008-11-28 | 2013-04-03 | ヤマハ発動機株式会社 | Remote control system and remote control device |
JP6422912B2 (en) | 2016-04-06 | 2018-11-14 | 株式会社クボタ | POSITIONING DETECTION DEVICE AND WORKING MACHINE HAVING POSITIONING DETECTION DEVICE |
-
2020
- 2020-06-02 JP JP2020096400A patent/JP7386757B2/en active Active
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7388580B1 (en) | 2023-02-15 | 2023-11-29 | いすゞ自動車株式会社 | Information output device and information output method |
CN117075485A (en) * | 2023-10-18 | 2023-11-17 | 舜泰汽车有限公司 | Agricultural machinery autopilot navigation equipment |
Also Published As
Publication number | Publication date |
---|---|
JP7386757B2 (en) | 2023-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8392104B2 (en) | Avoidance maneuver calculation device, avoidance control device, vehicle having each device, avoidance maneuver calculating method, and avoidance controlling method | |
US9090285B2 (en) | Method for providing a lanekeeping assistance based on modifying mechanical sources of steering torques | |
EP3467606B1 (en) | Autonomous work vehicle | |
CN108614551B (en) | Remote operation carrier and carrier control device and control method thereof | |
JP7072496B2 (en) | Control device for self-driving work vehicles | |
Chebly et al. | Coupled longitudinal/lateral controllers for autonomous vehicles navigation, with experimental validation | |
US20230211786A1 (en) | Path-controlling module, associated path-controlling device and associated method | |
JP2000302055A (en) | Traffic lane followup control device | |
JP2021189909A (en) | Autonomous travel system for work vehicle | |
JP2007022169A (en) | Vehicle controller and cant state deciding method | |
JP2011008385A (en) | Vehicle running support apparatus | |
KR20170124455A (en) | Travelling apparatus, control method of travelling apparatus, and control program of travelling apparatus | |
JP6850759B2 (en) | Autonomous steering device | |
JP6216545B2 (en) | Autonomous traveling vehicle, autonomous vehicle traveling system, and autonomous traveling vehicle control method | |
WO2020183988A1 (en) | Work vehicle | |
JP2019053470A (en) | Work vehicle-purpose autonomous travelling system | |
JPH0522926B2 (en) | ||
JP2022180938A (en) | Apparatus for setting travel route | |
Schein | Optimal Vehicle Control of Four-Wheel Steering | |
JP2020080743A (en) | Vehicular posture estimation device | |
US20230161355A1 (en) | Control system of dump truck | |
US20230311988A1 (en) | Systems and methods for machine steering | |
JP7115816B2 (en) | Autonomous vehicles and systems | |
JP7459733B2 (en) | Self-location estimation device | |
JP2023137764A (en) | Automatic travel vehicle |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230220 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20231018 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20231024 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20231114 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7386757 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |