








【0001】
【発明の属する技術分野】
本発明は、状態観測器を備えたトロイダル型無段変速機(以下、TCVTという。)に関するものである。
【0002】
【従来の技術】
状態観測器を用いてTCVTのトラニオン変位を推定する従来技術として、特開平8−270772号公報に記載の技術がある。この公報に記載の状態観測器は、変速アクチュエータの指令値(ステップモータのステップ数)uと傾転角度φとを入力とし、トラニオン軸方向変位yを推定する。従来技術では、推定対象であるTCVTを、次式のようにモデル化している。
【0003】
【数1】
式(1)
【0004】
【数2】
【0005】
ここで、a1、a2、bはTCVTの形状で決まる定数、f(φ、ωod)はφとωodとの非線形関数、gは変速制御弁のバルブゲインである。f(φ、ωod)は次式で表される。
【0006】
【数3】
式(2)
【0007】
ここで、η、θ、cfはTCVTの構造で決まる定数、ωodはTCVTの出力ディスク回転数である。このTCVTに対して、状態観測器を次のようにしている。
【0008】
【数4】
式(3)
【0009】
【数5】
【0010】
ここで、φtは傾転角度の推定値、ytはトラニオン変位の推定値、H0はオブザーバゲインである。オブザーバゲインは、次のように設定している。
【0011】
【数6】
式(4)
【0012】
ここで、ω0は、推定値fが真値xに収束する速さを決める正の定数である。
【0013】
そして、傾転角度φを目標値φ*に追従させる制御器は、次のような制御則を行い、目標値と傾転角度とトラニオンの推定値とに応じて、変速アクチュエータへの変位指令値を演算する。
【0014】
【数7】
【0015】
【数8】
式(5)
【0016】
【数9】
式(6)
【0017】
ここで、KPはPID制御器の比例ゲイン、KIはPID制御器の積分ゲイン、sはラプラス演算子、KDはPID制御器の微分ゲインである。
【0018】
【発明が解決しようとする課題】
式(1)から式(2)を両辺引いて、真値xと推定値xtとの推定誤差をe0=x−xtとすると、次式を得る。
【0019】
【数10】
式(7)
【0020】
この式(7)は、推定誤差の動特性を示す方程式である。この方程式の遷移行列(A0−H0C0)は、次式で表される。
【0021】
【数11】
式(8)
【0022】
式(8)と式(4)とから、(A0−H0C0)の固有値を求めると、(−ω0、−ω0)となる。
【0023】
しかしながら、この特開平8−270772号公報に記載の状態観測器において、式(8)に示す遷移行列は、傾転角度φとTCVT出力ディスク回転数ωodとに応じて変化する。これは、トラニオン変位推定値と実際のトラニオン変位との偏差の時間応答が、指数関数的に減衰しない非線形特性であることを示す。このため、偏差の時間応答に合わせて、理論的な裏付けがとれたコントローラを設計することが難しい。
【0024】
結果として、好ましい制御性能を達成するために、計算機シミュレーションや実験を繰り返してコントローラのパラメータ適合する必要があり、長い開発時間とコストを要する。つまり、短期間で理論的な裏付けがとれた制御装置を設計することが難しいという課題がある。
【0025】
さらに変速アクチュエータにステップモータを用いた場合、ステップモータが指令値どおりに動作せず、ステップ数指令値と実際のステップ数とのずれが定常的に残る脱調と呼ばれる現象が起きる場合がある。このずれを、状態観測器の入力uの定常的な外乱duと表す。式(1)にこの外乱を加えたTCVTモデルは、次式で表される。
【0026】
【数12】
式(9)
【0027】
式(9)から式(3)を引くと、次の誤差方程式を得る。
【0028】
【数13】
式(10)
【0029】
式(10)に示す誤差方程式は安定であると仮定する。このとき、推定誤差e0の定常値e00は、式(10)において微分値e´0=0として次式で表される。
【0030】
【数14】
【0031】
このように、脱調しているとき、推定値には定常偏差が残る。ステップモータのステップ数を定常的に指令値と一致させるには、センサ等を用いてステップモータが出力する変位を検出する方法もあるが、コストの増加を招いてしまう。適時ステップモータを基準位置へ移動して初期化を行う方法もあるが。走行中は初期化できないといった制約があるため、適切な時点で初期化できるとは限らない。つまり、安価にステップモータの実ステップ数を指令値に一致させることが難しいという課題がある。
【0032】
このような問題点を鑑み、本発明の目的は、制御誤差の小さい状態観測器を備えたTCVTを提供することである。
【0033】
【課題を解決するための手段】
第1の発明は、変速アクチュエータの指令値を入力とし、傾転角度または変速比を出力する無段変速機の状態方程式に基づいてトラニオン軸方向変位を推定する状態観測器と、このトラニオン軸方向変位推定値に基づいて変速アクチュエータの指令値を演算する指令値演算器とを備えるトロイダル型無段変速機において、前記状態観測器は、傾転角度を検出する傾転角度検出手段と、トラニオン軸方向変位に対する傾転角速度への感度のうち、傾転角度に依存する成分の逆数を現時刻の傾転角度まで積分する積分値算出手段とを有し、この算出された積分値を用いて前記トロイダル型無段変速機の状態方程式の遷移行列が定数行列となるように変換する。
【0034】
第2の発明は、第1の発明において、前記変速アクチュエータ指令値が、変速アクチュエータ変位であり、前記状態観測器は、トロイダル型無段変速機の回転数を検出する回転数検出手段と、トラニオン軸方向変位に対する傾転角速度への感度のうち、トロイダル型無段変速機の回転数に依存する成分f1を算出するf1算出手段と、前記積分値に比例し、トロイダル型無段変速機の回転数に依存する成分f1に反比例する新状態量を演算する新状態量演算手段と、この新状態量と変速アクチュエータの変位指令値とトロイダル型無段変速機の回転数とを入力とし、新状態量推定値とトラニオン変位推定値とを出力する状態観測用フィルタ手段とを有する。
【0035】
第3の発明は、第1の発明において、前記変速アクチュエータ指令値が、変速アクチュエータ駆動速度指令値であり、前記状態観測器は、トロイダル型無段変速機の回転数を検出または推定する回転数検出手段と、トラニオン軸方向変位に対する傾転角速度への感度のうち、トロイダル型無段変速機の回転数に依存する成分f1を算出するf1算出手段と、前記積分値に比例し、トロイダル型無段変速機の回転数に依存する成分f1に反比例する新状態量を演算する新状態量演算手段と、この新状態量と変速アクチュエータの駆動速度指令値とトロイダル型無段変速機の回転数とを入力とし、新状態量推定値とトラニオン変位推定値と変速アクチュエータ変位推定値とを出力する状態観測用フィルタ手段とを有する。
【0036】
第4の発明は、第1から3のいずれか一つの発明において、前記状態観測器が、トロイダル型無段変速機の回転数の時間微分値を演算する微分手段を有し、前記状態観測用フィルタ手段のフィルタ出力値をこの微分値で補正する。
【0037】
第5の発明は、第3または4の発明において、変速アクチュエータ駆動速度指令値の積分値と、変速アクチュエータ変位推定値とを比較して、変速アクチュエータの駆動状態を判断する駆動状態判断手段を有する。
【0038】
【発明の効果】
請求項1の発明によれば、状態観測器は、傾転角度検出手段と、トラニオン変位に対する傾転角速度への感度のうち、傾転角度に依存する成分を現時刻の傾転角度まで積分する積分値算出手段とを有し、この積分値を用いて、TCVTの状態方程式の遷移行列が定数行列となるように変換された状態方程式に基づいて構成するので、実際の値と推定値との誤差が指数関数的に減少することが理論的に証明されるため、制御装置を設計する際に行う計算機シミュレーションや実験が減り、開発時間とコストを削減できる。
【0039】
請求項2の発明によれば、前記変速アクチュエータ指令値は、変速アクチュエータ変位であり、前記状態観測器は、TCVTの回転数を検出または推定する回転数検出手段と、トラニオン軸方向変位に対する傾転角速度への感度のうち、TCVTの回転数に依存する成分f1を算出するf1算出手段と、前記積分値に比例し、TCVTの回転数に依存する成分f1に反比例する新状態量を演算する新状態量演算手段と、該新状態量と変速アクチュエータの変位指令値とTCVTの回転数とを入力とし、新状態量推定値とトラニオン変位推定値とを出力する状態観測用フィルタ手段とを有する構成としたので、目標値と傾転角度の偏差に応じて変速アクチュエータ位置を指令値として出力する従来技術の制御則の一部をそのまま使用できる。しかも、トラニオン変位の推定誤差が、一定の時定数のもとに指数関数的に減衰することが理論上裏付けられるため、制御装置の設計が一層容易になる。また、状態観測用フィルタは入力から出力への直達項を持たないことから、入力から出力までがローパスフィルタ特性となり、傾転角度に含まれる高周波数の検出ノイズを除去する性能を兼ね備える。
【0040】
請求項3の発明によれば、前記変速アクチュエータ指令値は、変速アクチュエータ駆動速度指令値であり、前記状態観測器は、TCVTの回転数を検出または推定する回転数検出手段と、トラニオン軸方向変位に対する傾転角速度への感度のうち、TCVTの回転数に依存する成分f1を算出するf1算出手段と、前記積分値に比例し、TCVTの回転数に依存する成分f1に反比例する新状態量を演算する新状態量演算手段と、該新状態量と変速アクチュエータの駆動速度指令値とTCVTの回転数とを入力とし、新状態量推定値とトラニオン変位推定値と変速アクチュエータ変位推定値とを出力する状態観測用フィルタ手段とを有する構成としたので、すべての推定値のずれは脱調の有無によらずゼロに収束するこれにより、ステップ数を検出するためのセンサを追加することによるコストの増加を防止できると共に、ステップ数の初期化を行う必要もなくなる。また、状態観測用フィルタも直達項がないため、傾転角度に含まれる高周波数の検出ノイズの除去する特長も有する。
【0041】
請求項4の発明によれば、前記状態観測器は、TCVTの回転数の時間微分値を演算する微分手段を有し、前記状態観測用フィルタ手段のフィルタ出力値を微分値で補正するので、変速の過渡状態における状態観測器の推定精度が向上する。
【0042】
請求項5の発明によれば、前記TCVTは、変速アクチュエータ駆動速度指令値の積分値と、変速アクチュエータ変位推定値とを比較して、変速アクチュエータの駆動状態を判断する駆動状態判断手段を有するので、ステップモータが駆動指令値どおりに駆動していないことが判断でき、モータ故障の判断等に利用できる。
【0043】
【発明の実施の形態】
以下、本発明の一実施形態を添付図面に基づいて説明する。
【0044】
図1はトロイダル型無段変速機10(以下、TCVTと示す。)の概略構成図であり、図2は、その断面、および、変速制御系の構成を示したものであり、図3はTCVTの変速比に対し一意に決まる傾転角度を目標値に従い制御するような変速指令値を決定するコントローラを含めた構成図であり、図4はパワーローラ軸方向変位を推定する状態観測器の構成図である。
【0045】
図1に示すTCVT10は、図中上方で連結される動力源としてのエンジン(図示せず)の回転が、トルクコンバータ12を介してTCVT10に入力されるようになっている。
【0046】
上記トルクコンバータ12は一般に良く知られているように、ポンプインペラ12a、タービンランナ12bおよびステータ12cを備え、特に該トルクコンバータ12ではロックアップクラッチ12dが設けられている。
【0047】
そして、TCVT10は、トルクコンバータ12の出力回転軸14と同軸上に配置されるトルク伝達軸16が設けられ、該トルク伝達軸16に第一トロイダル変速部18と第二トロイダル変速部20とがタンデム配置されている。
【0048】
上記トルク伝達軸16は中空に形成されると共に、ハウジング22に対し軸方向に若干の移動が可能に取付けられている。上記第一、第二トロイダル変速部18、20は、それぞれの対抗面がトロイダル曲面に形成される一対の第一入力ディスク18a、第一出力ディスク18bおよび第二入力ディスク20a、第二出力ディスク20bと、それぞれの対抗面間に摩擦接触されるパワーローラ18c、18d、および20c、20dとによって構成される。
【0049】
第一トロイダル変速部18は上記トルク伝達軸16の図中上方に配置されると共に、上記第二トロイダル変速部20は該トルク伝達軸16の図中下方に配置され、かつ、それぞれの第一入力ディスク18aおよび第二入力ディスク20aは互いに外側に配置されると共に、第一出力ディスク18bおよび第二出力ディスク20bは互いに内側に配置されている。そして、上記第一、第二入力ディスク18a、20aはボールスプライン24、26を介して上記トルク伝達軸16に、回転方向に係止されかつ軸方向の滑らかな移動が可能に取付けられている。
【0050】
一方、上記第一、第二出力ディスク18b、20bは、上記トルク伝達軸16と相対回転可能に軸支された出力ギア28にスプライン結合され、該第一、第二出力ディスク18b、20bに伝達された回転力は、出力ギア28と歯合する入力ギア30aを介してカウンターシャフト30に伝達され、さらに、回転力出力経路を介して図示しない出力軸に伝達される。
【0051】
ところで、上記第一入力ディスク18aの外側にはローディングカム装置34が設けられ、このローディングカム装置34には、回転力入力経路を介して伝達されるエンジン回転が入力され、この入力トルクに応じた押圧力がローディングカム装置34によって発生されるようになっている。なお、上記ローディングカム装置34のローディングカム34aは、上記トルク伝達軸16と相対回転可能に係合すると共に、スラストベアリング36を介して該トルク伝達軸16に係止される。
【0052】
また、上記第二入力ディスク20aと上記トルク伝達軸16の図中下方端部との間に皿バネ38が設けられている。従って、上記ローディングカム装置34で発生される押圧力は、第一入力ディスク18aに作用すると共に、上記トルク伝達軸16および上記皿バネ38を介して第二入力ディスク20aにも作用し、かつ、上記皿バネ38によって発生される予圧力は、第二入力ディスク20aに作用すると共に、上記トルク伝達軸16および上記ローディングカム装置34を介して第一入力ディスク18aにも作用するようになっている。
【0053】
ところで、上記ローディングカム装置34と上記トルクコンバータ12との間の回転力入力経路には、車両の前進時と後進時の回転方向を切り換える前後進切換装置40が設けられる。
【0054】
上記前後進切換装置40は、ダブルプラネタリー方式の遊星歯車機構42と、該遊星歯車機構42のキャリア42aを上記出力回転軸14に締結可能なフォワードクラッチ44と、該遊星歯車機構42のリングギア42bを上記ハウジング22に締結可能なリバースブレーキ46とによって構成される。
【0055】
そして、上記前後進切換装置40では、フォワードクラッチ44を締結すると共に、リバースブレーキ46を開放することにより、エンジン回転と同方向の回転が上記ローディングカム装置に入力され、かつ、フォワードクラッチ44を開放してリバースブレーキ46を締結することにより、逆方向の回転が入力されるようになっている。なお、上記遊星歯車機構42で、42cはサンギア、42d、42eは互いに噛み合いされるプラネタリギア、42bはリングギアである。
【0056】
ところで、上記第一トロイダル変速部18および第二トロイダル変速部20に設けられたパワーローラ18c、18dおよび20c、20dは、中心軸cに対して対称に配置され、それぞれのパワーローラは変速制御装置としての変速制御弁56および油圧サーボシリンダ50を介して、車両運転条件に応じて傾斜(傾転)され、もって前記第一、第二入力ディスク18a、20aの回転を無段階に変速して前記第一、第二出力ディスク18b、20bに伝達するようになっている。
【0057】
次に、図2はTCVTの変速を管理する油圧系の機械的構成図である。
【0058】
パワーローラ20c、20dはそれぞれトラニオン23、23で背面側から軸支されている。トラニオン23は油圧サーボシリンダ50のサーボピストン51と結合しており、油圧サーボシリンダ50内の油室50a内の油と油室50b内の油の差圧で変位する。油室50a、50bは、それぞれシフトコントロールバルブ(変速制御弁)56のHi側ポート56HiとLow側ポート56Lowに繋がっており、シフトコントロールバルブ56はバルブ内のスプール56Sが変位することにより、ライン圧の油をHi側ポート56HiまたはLow側ポート56Lowに流し、他方のポートからドレーン56Dへ油を流出させることで油圧サーボ内の差圧を変化させる。
【0059】
スプール56Sは、ステップモータ52(変速アクチュエータ)と後述するプリセスカム55とリンク構造で連結している。プリセスカム55は、4体のトラニオンのうちの1体に取り付けられており、パワーローラの上下方向変位とパワーローラの傾転角度をリンクの変位に変換し、スプール56Sへ伝達する。
【0060】
スプール56Sの変位はステップモータ52の変位(駆動位置)とプリセスカム55で伝えられる変位により決まる。TCVTの変速は、トラニオン23を平衡点(パワーローラの回転軸と入出力ディスクの回転軸が交差する状態で、図2の状態を示す)から上下に変位させることにより行い、この変位によりパワーローラ20cと両ディスク20aの回転方向ベクトルに差違が発生してパワーローラ20cは傾転する。
【0061】
変速の定常時には、パワーローラ20cおよびトラニオン23の軸方向変位yは平衡点に戻っており、スプール56Sの変位も中立点でバルブが閉じた状態であるので、このとき傾転角度は、プリセスカム比とリンク比で決まるステップモータ52の変位に対応した位置となる。
【0062】
プリセスカム55は、パワーローラ20cの傾転角度をスプール56Sの変位に負帰還して傾転角度の目標値とのずれを補償しながら、パワーローラ20cおよびトラニオン23の平衡点からの変位もスプール56Sの変位に負帰還して過渡状態においてダンピングの効果を与えて、変速のハンチングを抑える。
【0063】
すなわち、変速の到達点はステップモータ52の変位で決まり、一連の変速の過程を示すと、ステップモータ変位を変化させることでスプール56Sが変位してバルブが開き、サーボピストンの差圧が変化してパワーローラが平衡点から変位することで傾転し、傾転角度がステップモータ変位に対応した点でスプール56Sは中立点に戻り変速が終了する。
【0064】
図3は制御装置80を含んだTCVT10の構成図である。
【0065】
入力ディスク回転数センサ84は、入力ディスク18a、21aの何れか1つの回転に同期して発生するパルス信号を周期もしくは周波数として計測し、入力ディスク回転数ωidを検出する。出力ディスク回転数センサ83は、出力ディスク18b、21bの何れか1つの回転に同期して発生するパルス信号を周期もしくは周波数として計測し、出力ディスク回転数ωodを検出する。パワーローラ回転数センサ82は、パワーローラ18c、18d、20c、20dの何れか1つの回転に同期して発生するパルス信号を、周期もしくは周波数として計測してパワーローラ回転数ωprを検出する。
【0066】
傾転角度センサ85は、ロータリエンコーダ等を用いて傾転角度φを検出する。
【0067】
マイクロコンピュータを主体に構成された制御装置80は、入力ディスク回転数ωidと、出力ディスク回転数ωodと、パワーローラ回転数ωprと、傾転角度φを入力して、ステップモータ52の指令値を演算する。
【0068】
図4に、制御装置80を構成する本発明の状態観測器110の一例を示す。状態観測器110は、以下に説明する各手段からなり、状態推定値を指令値演算器200に出力する。
【0069】
傾転角度検出手段100は、例えば、前記傾転角度センサ85の出力を傾転角度φとする。もしくは、変速比icから、次式の変速比icと傾転角度φとの関係を用いて演算する。
【0070】
【数15】
式(11)
【0071】
ここで、η、θはTCVT10の形状で決まる定数である。式(11)に示す変速比icと傾転角度φとの関係を図示すると、図5となる。変速比icは、例えば、入力ディスク回転数センサの出力ωidと出力ディスク回転数センサの出力ωod0とから、次式を用いて演算するとよい。
【0072】
【数16】
式(12)
【0073】
また出力ディスク回転数検出値ωod0と入力ディスク回転数ωidとパワーローラ回転数ωprと傾転角度φとの間には、次のような関係がある。
【0074】
【数17】
【0075】
【数18】
【0076】
この関係を用いて、パワーローラ回転数ωprの検出値と傾転角度φの検出値とから、出力ディスク回転数ωod0と入力ディスク回転数ωidを算出し、式(12)に示す関係を用いて算出してもよい。
【0077】
回転数検出手段101では、例えば、出力ディスク回転数センサの出力を、TCVTの回転数検出値ωod0とする。
【0078】
微分手段102では、TCVT10の回転数検出値ωod0を入力として、次式で示すローパスフィルタの演算を行う。
【0079】
【数19】
式(13)
【0080】
式(13)で示すフィルタの出力ωod0をTCVT10の回転数推定値とし、中間変数のω´odをTCVT10の回転数微分値として出力する。
【0081】
積分値算出手段104では、傾転角度φを入力して、トラニオン軸方向変位yに対する傾転角速度への感度のうち、傾転角度φに依存する成分を現時刻の傾転角度まで積分する。トラニオン軸方向変位に対する傾転角速度への感度fは、前述したように式(2)で表される。この感度fは、次のように、TCVT10の回転数に依存する成分f1と傾転角度に依存する成分f2と定数部分cfとに分けられる。
【0082】
【数20】
式(14)
【0083】
【数21】
式(15)
【0084】
【数22】
式(16)
【0085】
現時刻の傾転角度φまでの、f2(φ)の逆数の傾転角度に関する積分値q(φ)は次式を用いて求める。
【0086】
【数23】
式(17)
【0087】
f1算出手段103では、TCVT10の回転数ωod0を入力して、トラニオン軸方向変位に対する傾転角速度への感度のうち、TCVT10の回転数ωodに依存する成分f1を、式(15)を用いて算出する。
【0088】
新状態量算出手段105では、前記積分値q(φ)とf1(ωod0)とを入力し、次式を用いて、新状態量ζ1を演算する。
【0089】
【数24】
式(18)
【0090】
状態観測用フィルタ手段106では、新状態量ζ1と、後述する指令値演算器200の出力であるステップモータ変位指令値uと傾転角度φとを入力とし、次式で表されるフィルタを用いて、新状態量ζ1とトラニオン変位yとの推定値ζ1t、ytを演算する。
【0091】
【数25】
式(19)
【0092】
【数26】
【0093】
Haはオブザーバゲインである。オブザーバゲインは、次のように設定する。
【0094】
【数27】
式(20)
【0095】
ここで、p10、p11はオブザーバの極に関係する定数であり、オブザーバ極は、
s=−p10、s=−p11
となる。
【0096】
以下に本発明の状態観測器110の導出と、特徴を示す。
【0097】
ステップモータのステップ数uを入力とし、傾転角度φとトラニオン変位yとを状態量として、TCVT10の動特性は式(1)で表される。式(18)を両辺時間微分すると、次式を得る。
【0098】
【数28】
式(21)
【0099】
式(21)と式(1)と式(14)とから、次式を得る。
【0100】
【数29】
式(22)
【0101】
このように、式(18)に示す状態変換を用いて、式(22)と式(1)とから、TCVTの動特性は、状態量を、
【0102】
【数30】
【0103】
として、次のように書き直すことができる。
【0104】
【数31】
式(23)
【0105】
以上の状態変換をまとめると、式(14)で表されるように、トラニオン変位yに対する傾転角速度の感度f(φ、ωod)を、式(16)で表される傾転角度に依存する成分f2(φ)とその他の部分と分け、新状態量を式(18)で表されるように設定する。このとき、この新状態量の微分値は式(21)で表され、傾転角速度φ´の係数は、式(14)で表されるf(φ、ωod0)の逆数となる。
【0106】
これにより、式(1)と式(21)とから、トラニオン変位yの係数の時変項がキャンセルされて、式(22)に示すように係数は一定値となる。そして、新状態量を、傾転角度の状態量の代わりとしてTCVTの動特性を考えると、式(23)で表されるように、遷移行列Aaが定数となる。
【0107】
この書き直されたTCVT10に対して、式(19)に示すオブザーバは設計した。
【0108】
ここで、状態量ζaの推定値誤差eaを、
【0109】
【数32】
式(24)
【0110】
として、式(23)から式(19)を引くと次式を得る。
【0111】
【数33】
式(25)
【0112】
式(25)は、誤差eaの減衰特性を表す誤差方程式である。Aa、Ha、Ca共に定数行列であるため、(Aa−HaCa)も定数行列である。これは、誤差eaが指数関数的に減衰することを示す。
【0113】
したがって、本発明の状態観測器110は、トラニオンの傾転角度を検出または推定する傾転角度検出手段100と、トラニオン変位に対する傾転角速度への感度のうち、傾転角度に依存する成分f2を現時刻の傾転角度まで積分する積分値算出手段104とを有し、この積分値を用いて、TCVT10の状態方程式の遷移行列が定数行列となるように変換された状態方程式に基づいて構成するので、実際の値と推定値との誤差が指数関数的に減少することが理論的に証明されるため、制御装置を設計する際に行う計算機シミュレーションや実験が減り、開発時間とコストを削減できる。
【0114】
さらに本発明の変速アクチュエータ指令値は、変速アクチュエータ変位であり、状態観測器110は、TCVT10の回転数を検出または推定する回転数検出手段101と、トラニオン軸方向変位に対する傾転角速度への感度のうち、TCVT10の回転数に依存する成分f1を算出するf1算出手段103と、前記積分値に比例し、回転数に依存する成分f1に反比例する新状態量を演算する新状態量演算手段105と、この新状態量と変速アクチュエータの変位指令値とTCVT10の回転数とを入力とし、新状態量推定値とトラニオン変位推定値とを出力する状態観測用フィルタ手段106とを有する構成としたので、目標値と傾転角度の偏差に応じて変速アクチュエータ位置を指令値として出力する従来技術の制御則の一部をそのまま使用できる。しかも、トラニオン変位の推定誤差が、一定の時定数のもとに指数関数的に減衰することが理論上裏付けられるため、制御装置の設計が一層容易になる。また、状態観測用フィルタ手段106は入力から出力への直達項を持たないことから、入力から出力までがローパスフィルタ特性となり、傾転角度に含まれる高周波数の検出ノイズを除去する性能を兼ね備える。
【0115】
次に、指令値演算器200では、状態観測用フィルタ手段106で演算したトラニオン変位推定値ytを入力して、式(5)、式(6)とで表される制御則を用いてステップモータ変位指令値uを出力する。
【0116】
以下において、制御装置で演算する変速制御演算を、図7に示すフローチャートを使い説明する。この変速制御演算は、ある所定の制御周期例えば10ms毎に実行される。
【0117】
ステップS100、S101では、TCVT入出力ディスク回転数ωid、ωco0を読み込む。
【0118】
ステップS102では、先ず、傾転角度検出手段100がTCVT入力ディスク回転数ωidとTCVT出力ディスク回転数検出値ωod0とから、式(12)を用いて変速比icを演算する。次に、図5に示す変速比icと傾転角度φとの関係を示すマップを用いて、変速比icから傾転角度φを算出する。
【0119】
ステップS103では、微分手段102がTCVT出力ディスク回転数検出値ωod0を入力し、式(13)で表されるローパスフィルタ演算を行い、TCVT出力ディスク回転数推定値ωodとTCVT出力ディスク回転数微分値ω´odとを算出する。
【0120】
ステップS104では、積分値算出手段104が、式(15)を用いて、TCVT出力ディスク回転数推定値ωodから回転数に依存する成分f1を算出する。
【0121】
ステップS105では、同じく積分値算出手段104が、傾転角度φを入力し、式(17)を用いて、積分値q(φ)を算出する。あるいは、予め式(17)を演算して求めた図6に示すマップを用いて、傾転角度φから積分値q(φ)を求める。
【0122】
ステップS106では、新状態量算出手段105が、式(18)を用いて、f1とq(φ)とから新状態量ζ1を算出する。
【0123】
ステップS107では、状態観測用フィルタ手段106が、傾転角度φとステップ数指令値uと新状態量ζ1とTCVT出力ディスク回転数推定値ωodとTCVT出力ディスク回転数微分値ω´odとを入力し、式(19)で表されるトラニオン変位推定値ytを算出する。
【0124】
ステップS108では、指令値演算器200が、トラニオン変位推定値ytと傾転角度φと目標値φ*とを入力し、式(5)と式(6)とで表される制御則を演算し、ステップ数指令値uを算出する。
【0125】
したがって、状態観測器110は、TCVT10の回転数の時間微分値ω´odを演算する微分手段102を有し、状態観測用フィルタ手段106のフィルタ出力値を微分値で補正するので、変速の過渡状態における状態観測器の推定精度が向上する。
【0126】
次に第2の実施形態について説明する。
【0127】
第2実施形態の状態観測器110は、第1の実施形態に比して状態観測用フィルタ手段106の演算方法が異なるのみであり、図1から3に示す構成及び図4の傾転角度検出手段100、回転数検出手段101、微分手段102、積分値算出手段103、f1算出手段104、新状態量算出手段105は、第1実施形態と同じであるので説明は省略する。
【0128】
第2実施形態の状態観測用フィルタ手段106では、新状態量ζ1と、後述する指令値演算器の出力であるステップモータ駆動速度指令値vと傾転角度φとを入力とし、次式で表されるフィルタを用いて、新状態量ζ1とトラニオン変位yとステップモータのステップ数uとを演算する。
【0129】
【数34】
式(26)
【0130】
【数35】
【0131】
ここで、Hbはオブザーバゲインである。オブザーバゲインは、次のように設定する。
【0132】
【数36】
式(27)
【0133】
ここで、pb0、pb1、pb2はオブザーバの極に関係する定数であり、オブザーバ極は、
s=−pb0、s=−pb1、s=−pb2
となる。
【0134】
以下に、この状態推定手段の導出と、特徴を示す。ステップモータステップ数uを入力とし、傾転角度φとトラニオン変位yとを状態量として、TCVTの動特性は式(1)で表される。また、ステップモータ駆動速度vとステップモータステップ数uとの関係は次式で表される。
【0135】
【数37】
式(28)
【0136】
実施例1と同様に、式(18)に示す状態変換を用いて、式(22)と式(1)と式(28)とから、TCVT10の動特性は、状態量
【0137】
【数38】
【0138】
として、次のように書き直すことができる。
【0139】
【数39】
式(29)
【0140】
このように第1実施形態と同じ状態変換により、遷移行列Abは定数となる。この書き直されたTCVT10に対して、式(26)に示すオブザーバは設計した。ζbの推定値誤差ebを、
【0141】
【数40】
式(30)
【0142】
として、式(29)から式(26)を引くと次式を得る。
【0143】
【数41】
式(31)
【0144】
式(26)は、誤差ebの減衰特性を表す誤差方程式である。Ab、Hb、Cb共に定数行列であるため、(Ab−HbCb)も定数行列である。これは、誤差ebが指数関数的に減衰することを示す。
【0145】
指令値演算器200では、状態観測用フィルタ手段で演算したトラニオン変位推定値とステップ数推定値とを入力して、ステップモータ駆動速度指令値を出力する。先ず、傾転角度φとトラニオン変位推定値ytとステップ数推定値utから、次に示す微分同相写像を用いて、傾転角速度推定値φ´tと傾転角加速度推定値φ´´tとを演算する。
【0146】
【数42】
式(32)
【0147】
【数43】
式(33)
【0148】
ここで、f´は次式で表される。
【0149】
【数44】
式(34)
【0150】
次に、次式に示す制御則を用いて、ステップモータ駆動速度指令値vを演算する。
【0151】
【数45】
式(35)
【0152】
【数46】
式(36)
【0153】
ここで、φ*は目標傾転角度、ωnは目標傾転角度に対する傾転角度の応答速度を決める定数、kはスイッチングゲイン、εはスイッチングゲインを制御誤差σのゼロ近傍で連続化する定数である。式(35)、式(36)はスライディングモード制御則である。式(35)により、制御誤差がゼロになる方向ヘステップモータは駆動される。式(36)において、σ=0とすると次式を得る。
【0154】
【数47】
式(37)
【0155】
このように、σ=0のとき、目標傾転角度に対して傾転角度は式(37)に示す動特性で応答する。
【0156】
また、ステップモータが指令値どおりに動作せず、脱調が起こっている場合を考える。式(29)に脱調による外乱duを加えたTCVTモデルは、次式で表される。
【0157】
【数48】
式(38)
【0158】
【数49】
【0159】
式(38)から式(26)を引くと、次の誤差方程式を得る。
【0160】
【数50】
式(39)
【0161】
式(39)に示す誤差方程式は安定であると仮定する。このとき、誤差ebの定常値eb0は、式(39)においてe´b=0として次式で表される。
【0162】
【数51】
式(40)
【0163】
このように、脱調が起こっても、傾転角度φの定常誤差とトラニオン変位yの定常誤差とはゼロである。さらに、式(40)から、ステップ数の推定値utは、
【0164】
【数52】
【0165】
となる。これは、外乱を含んだ実際のステップ数が推定できることを示しており、ステップ数の推定誤差はない。
【0166】
したがって本実施形態によれば、変速アクチュエータ指令値は、変速アクチュエータ駆動速度指令値であり、状態観測器110は、TCVT10の回転数を検出または推定する回転数検出手段101と、トラニオン軸方向変位に対する傾転角速度への感度のうち、TCVT10の回転数に依存する成分f1を算出するf1算出手段103と、前記積分値に比例し、f1に反比例する新状態量を演算する新状態量演算手段105と、この新状態量と変速アクチュエータの駆動速度指令値とTCVT10の回転数とを入力とし、新状態量推定値とトラニオン変位推定値と変速アクチュエータ変位推定値とを出力する状態観測用フィルタ手段106とを有する構成としたので、すべての推定値のずれは脱調の有無によらずゼロに収束するこれにより、ステップ数を検出するためのセンサを追加することによるコストの増加を防止できると共に、ステップ数の初期化を行う必要もなくなる。また、状態観測用フィルタも直達項がないため、傾転角度に含まれる高周波数の検出ノイズの除去できる。
【0167】
次に、図9のような駆動状態判断手段120を考える。駆動状態判断手段120は、ステップモータ駆動パルスの積算値をuItとし、uItと前記状態観測用フィルタ手段の出力utを比較する。そして、例えば、定常時にutとuItとが大きく異なるなら、脱調していると判断する。また例えば、utの変化とuItの変化とを比較して、ステップモータ異常(ステップモータの固着等)の判断を行う。
【0168】
以下において、制御装置80で演算する変速制御演算を、図7に示すフローチャートを使い説明する。この変速制御演算は、ある所定の制御周期、例えば10mS毎に実行される。
【0169】
ステップS100からステップS106は、第1実施形態と同じであるため、説明は省略する。
【0170】
ステップS107では、傾転角度φとステップモータ駆動速度指令値vとζ1とTCVT出力ディスク回転数推定値ωodとTCVT出力ディスク回転数微分値ω´odとを入力し、式(26)で表される状態観測器110を演算し、トラニオン変位推定値ytとステップモータステップ数推定値utとを算出する。
【0171】
ステップS108では、トラニオン変位推定値ytとステップモータステップ数推定値utと傾転角度φと目標値φ*とを入力し、式(32)から式(36)で表される制御則を演算し、ステップモータ駆動速度指令値vを算出する。
【0172】
図8に、この状態観測器110を用いて行った計算機シミュレーションの結果を示す。ステップ数に初期推定誤差duを与えた。傾転角度φ、トラニオン変位y、ステップ数uとも、推定誤差は安定に減衰し、指数関数的に減少することが確認できる。
【0173】
したがって、TCVT10は、変速アクチュエータ駆動速度指令値の積分値と、変速アクチュエータ変位推定値とを比較して、変速アクチュエータの駆動状態を判断する駆動状態判断手段を有するので、ステップモータが駆動指令値どおりに駆動していないことが判断でき、モータ故障の判断等に利用できる。
【0174】
本発明は、上記した実施形態に限定されるものではなく、本発明の技術的思想の範囲内でさまざまな変更がなしうることは明白である。
【図面の簡単な説明】
【図1】無段変速機の概略構成図である。
【図2】無段変速機の変速制御油圧系構成図である。
【図3】制御装置を含んだ制御系構成図である。
【図4】状態観測器の構成図である。
【図5】傾転角度と変速比との関係図である。
【図6】傾転角度とf2の逆数の積分値との関係図である。
【図7】制御装置の変速制御演算のフローチャートである。
【図8】計算機による本発明を用いたシミュレーション結果である。
【図9】駆動状態判断手段を含んだ制御系構成図である。
【符号の説明】
10 トロイダル型無段変速機
52 ステップモータ
56 変速制御弁
100 傾転角度検出装置
101 回転数検出手段
102 微分手段
103 f1算出手段
104 積分値算出手段
105 新状態量算出手段
106 状態観測用フィルタ手段
110 状態観測器
200 指令値演算器[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a toroidal continuously variable transmission (hereinafter referred to as TCVT) provided with a state observer.
[0002]
[Prior art]
As a conventional technique for estimating a TCVT trunnion displacement using a state observer, there is a technique described in Japanese Patent Laid-Open No. 8-270772. The state observer described in this publication estimates the trunnion axial displacement y using the command value (step motor step number) u of the speed change actuator and the tilt angle φ as inputs. In the prior art, the TCVT to be estimated is modeled as the following equation.
[0003]
[Expression 1]
Formula (1)
[0004]
[Expression 2]
[0005]
Where a1, A2, B are constants determined by the shape of TCVT, f (φ, ωod) Is φ and ωodAnd g is a valve gain of the shift control valve. f (φ, ωod) Is expressed by the following equation.
[0006]
[Equation 3]
Formula (2)
[0007]
Where η, θ, cfIs a constant determined by the structure of TCVT, ωodIs the TCVT output disk speed. For this TCVT, the state observer is as follows.
[0008]
[Expression 4]
Formula (3)
[0009]
[Equation 5]
[0010]
Where φtIs the estimated tilt angle, ytIs the estimated value of trunnion displacement, H0Is the observer gain. The observer gain is set as follows.
[0011]
[Formula 6]
Formula (4)
[0012]
Where ω0Is a positive constant that determines the speed at which the estimated value f converges to the true value x.
[0013]
Then, the tilt angle φ is set to the target value φ*The controller that follows the control performs the following control law, and calculates a displacement command value for the speed change actuator according to the target value, the tilt angle, and the estimated value of the trunnion.
[0014]
[Expression 7]
[0015]
[Equation 8]
Formula (5)
[0016]
[Equation 9]
Formula (6)
[0017]
Where KPIs the proportional gain of the PID controller, KIIs the integral gain of the PID controller, s is the Laplace operator, KDIs the differential gain of the PID controller.
[0018]
[Problems to be solved by the invention]
By subtracting both sides of equation (2) from equation (1), true value x and estimated value xtThe estimated error from0= XxtThen, the following equation is obtained.
[0019]
[Expression 10]
Formula (7)
[0020]
This equation (7) is an equation indicating the dynamic characteristics of the estimation error. Transition matrix (A0-H0C0) Is expressed by the following equation.
[0021]
## EQU11 ##
Formula (8)
[0022]
From Equation (8) and Equation (4), (A0-H0C0) Eigenvalue of (−ω0, -Ω0)
[0023]
However, in the state observer described in JP-A-8-270772, the transition matrix shown in the equation (8) has a tilt angle φ and a TCVT output disk rotational speed ω.odIt changes according to. This indicates that the time response of the deviation between the trunnion displacement estimated value and the actual trunnion displacement is a non-linear characteristic that does not decay exponentially. For this reason, it is difficult to design a controller that is theoretically supported in accordance with the time response of the deviation.
[0024]
As a result, in order to achieve the desired control performance, it is necessary to repeatedly perform computer simulations and experiments to adapt the parameters of the controller, which requires a long development time and cost. That is, there is a problem that it is difficult to design a control device that is theoretically supported in a short period of time.
[0025]
Further, when a step motor is used as the speed change actuator, there is a case where a phenomenon called step-out occurs in which the step motor does not operate according to the command value, and the deviation between the step number command value and the actual step number remains constant. This deviation is expressed as a steady disturbance d of the state observer input u.uIt expresses. A TCVT model obtained by adding this disturbance to the equation (1) is expressed by the following equation.
[0026]
[Expression 12]
Formula (9)
[0027]
Subtracting equation (3) from equation (9) yields the following error equation:
[0028]
[Formula 13]
Formula (10)
[0029]
It is assumed that the error equation shown in equation (10) is stable. At this time, the estimation error e0Steady value e00Is the differential value e ′ in equation (10).0= 0 and is represented by the following formula.
[0030]
[Expression 14]
[0031]
Thus, when stepping out, a steady-state deviation remains in the estimated value. Although there is a method of detecting the displacement output by the step motor using a sensor or the like in order to make the number of steps of the step motor consistent with the command value, the cost increases. There is also a method to perform initialization by moving the stepping motor to the reference position in a timely manner. Since there is a restriction that initialization cannot be performed during traveling, initialization cannot be performed at an appropriate time. That is, there is a problem that it is difficult to make the actual step number of the step motor coincide with the command value at a low cost.
[0032]
In view of such problems, an object of the present invention is to provide a TCVT including a state observer with a small control error.
[0033]
[Means for Solving the Problems]
According to a first aspect of the present invention, there is provided a state observer for estimating a trunnion axial displacement based on a state equation of a continuously variable transmission that receives a command value of a speed change actuator and outputs a tilt angle or a gear ratio, and the trunnion axial direction In a toroidal-type continuously variable transmission comprising a command value calculator for calculating a command value of a speed change actuator based on an estimated displacement value, the state observer includes a tilt angle detecting means for detecting a tilt angle, a trunnion shaft Of the sensitivity to the tilt angular velocity with respect to the directional displacement, the integrated value calculating means for integrating the reciprocal of the component depending on the tilt angle up to the tilt angle at the current time, and using the calculated integral value, The state equation of the toroidal continuously variable transmission is transformed so that the transition matrix becomes a constant matrix.
[0034]
According to a second invention, in the first invention, the shift actuator command value is a shift actuator displacement, and the state observer includes a rotation speed detection means for detecting the rotation speed of the toroidal type continuously variable transmission, a trunnion Of the sensitivity to the tilt angular velocity with respect to the axial displacement, the component f that depends on the rotational speed of the toroidal-type continuously variable transmission f1F to calculate1A calculation means and a component f proportional to the integral value and dependent on the rotational speed of the toroidal continuously variable transmission1New state quantity calculation means for calculating a new state quantity inversely proportional to the input, the new state quantity, the displacement command value of the shift actuator, and the rotational speed of the toroidal continuously variable transmission, and the new state quantity estimated value and trunnion displacement State observation filter means for outputting an estimated value.
[0035]
In a third aspect based on the first aspect, the shift actuator command value is a shift actuator drive speed command value, and the state observer detects or estimates the rotational speed of the toroidal continuously variable transmission. Of the sensitivity to the tilting angular velocity with respect to the detection means and the trunnion axial displacement, the component f depends on the rotational speed of the toroidal continuously variable transmission.1F to calculate1A calculation means and a component f proportional to the integral value and dependent on the rotational speed of the toroidal continuously variable transmission1New state quantity calculation means for calculating a new state quantity inversely proportional to the input, the new state quantity, the drive speed command value of the speed change actuator, and the rotational speed of the toroidal continuously variable transmission, and the new state quantity estimated value and trunnion State observation filter means for outputting the estimated displacement value and the estimated shift actuator displacement value.
[0036]
According to a fourth invention, in any one of the first to third inventions, the state observer has a differentiating means for calculating a time differential value of the rotational speed of the toroidal-type continuously variable transmission. The filter output value of the filter means is corrected with this differential value.
[0037]
According to a fifth invention, in the third or fourth invention, there is provided drive state determining means for comparing the integral value of the shift actuator drive speed command value with the estimated shift actuator displacement value to determine the drive state of the shift actuator. .
[0038]
【The invention's effect】
According to the first aspect of the present invention, the state observer integrates the tilt angle detection means and the component depending on the tilt angle of the sensitivity to the tilt angular velocity with respect to the trunnion displacement to the tilt angle at the current time. An integral value calculating means, and using this integral value, the transition matrix of the state equation of TCVT is configured based on the state equation converted into a constant matrix, so that the actual value and the estimated value Since it is theoretically proved that the error decreases exponentially, the number of computer simulations and experiments performed when designing the control device is reduced, and the development time and cost can be reduced.
[0039]
According to a second aspect of the present invention, the shift actuator command value is a shift actuator displacement, and the state observer detects a rotational speed of the TCVT, and a tilt relative to the trunnion axial displacement. Of the sensitivity to angular velocity, the component f that depends on the rotational speed of the TCVT1F to calculate1A calculation means and a component f proportional to the integral value and dependent on the rotational speed of the TCVT1New state quantity calculation means for calculating a new state quantity inversely proportional to the input, the new state quantity, a displacement actuator displacement command value and the TCVT rotation speed as inputs, and a new state quantity estimated value and a trunnion displacement estimated value are output. Therefore, a part of the control law of the prior art that outputs the shift actuator position as a command value according to the deviation between the target value and the tilt angle can be used as it is. Moreover, since it is theoretically supported that the estimation error of the trunnion displacement attenuates exponentially with a constant time constant, the design of the control device becomes easier. In addition, since the state observation filter does not have a direct term from input to output, it has low-pass filter characteristics from input to output, and has the ability to remove high-frequency detection noise included in the tilt angle.
[0040]
According to a third aspect of the present invention, the shift actuator command value is a shift actuator drive speed command value, and the state observer includes a rotation speed detection means for detecting or estimating the rotation speed of the TCVT, and trunnion axial displacement. Of the sensitivity to the tilt angular velocity for TCVT, the component f that depends on the rotational speed of TCVT1F to calculate1A calculation means and a component f proportional to the integral value and dependent on the rotational speed of the TCVT1New state quantity calculating means for calculating a new state quantity inversely proportional to the input, the new state quantity, the drive speed command value of the shift actuator, and the TCVT rotation speed, and the new state quantity estimated value, trunnion displacement estimated value, and shift Since the state observation filter means for outputting the actuator displacement estimated value is provided, the deviation of all estimated values converges to zero regardless of the presence or absence of step-out, thereby detecting the number of steps. It is possible to prevent an increase in cost due to the addition of and eliminate the need to initialize the number of steps. Further, since the state observation filter has no direct term, it has a feature of removing high-frequency detection noise included in the tilt angle.
[0041]
According to the invention of claim 4, the state observer has a differentiating means for calculating a time differential value of the rotational speed of the TCVT, and corrects the filter output value of the state observation filter means with the differential value. The estimation accuracy of the state observer in the transition state of the shift is improved.
[0042]
According to the fifth aspect of the present invention, the TCVT has drive state determination means for comparing the integral value of the shift actuator drive speed command value with the estimated shift actuator displacement value to determine the drive state of the shift actuator. Therefore, it can be determined that the step motor is not driven in accordance with the drive command value, which can be used for determining a motor failure or the like.
[0043]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, an embodiment of the present invention will be described with reference to the accompanying drawings.
[0044]
FIG. 1 is a schematic configuration diagram of a toroidal-type continuously variable transmission 10 (hereinafter referred to as TCVT), FIG. 2 shows its cross section and the configuration of a shift control system, and FIG. 3 shows TCVT. FIG. 4 is a block diagram including a controller that determines a shift command value that controls a tilt angle uniquely determined according to a target gear ratio according to a target value, and FIG. 4 is a configuration of a state observer that estimates a power roller axial displacement FIG.
[0045]
[0046]
As is generally known, the
[0047]
The
[0048]
The
[0049]
The first
[0050]
On the other hand, the first and
[0051]
By the way, a
[0052]
A
[0053]
By the way, a forward /
[0054]
The forward /
[0055]
In the forward /
[0056]
Incidentally, the
[0057]
Next, FIG. 2 is a mechanical configuration diagram of a hydraulic system that manages the shift of the TCVT.
[0058]
The
[0059]
The
[0060]
The displacement of the
[0061]
At the steady state of the speed change, the axial displacement y of the
[0062]
The
[0063]
That is, the reaching point of the speed change is determined by the displacement of the
[0064]
FIG. 3 is a configuration diagram of the
[0065]
The input disk
[0066]
The
[0067]
The
[0068]
FIG. 4 shows an example of the
[0069]
For example, the tilt
[0070]
[Expression 15]
Formula (11)
[0071]
Here, η and θ are constants determined by the shape of the
[0072]
[Expression 16]
Formula (12)
[0073]
Also output disk rotation speed detection value ωod0And input disk speed ωidAnd power roller rotation speed ωprAnd the tilt angle φ have the following relationship.
[0074]
[Expression 17]
[0075]
[Formula 18]
[0076]
Using this relationship, the power roller rotation speed ωprOutput disk rotational speed ω from the detected value of tilt and the detected value of tilt angle φod0And input disk speed ωidMay be calculated using the relationship shown in Equation (12).
[0077]
In the rotational speed detection means 101, for example, the output of the output disk rotational speed sensor is converted into the rotational speed detection value ω of the TCVT.od0And
[0078]
In the differentiating means 102, the rotational speed detection value ω of the
[0079]
[Equation 19]
Formula (13)
[0080]
The output ω of the filter shown in equation (13)od0Is the estimated rotational speed of TCVT10, and the intermediate variable ω ′odIs output as the rotational speed differential value of TCVT10.
[0081]
The integral value calculation means 104 inputs the tilt angle φ and integrates the component depending on the tilt angle φ of the sensitivity to the tilt angular velocity with respect to the trunnion axial displacement y to the tilt angle at the current time. The sensitivity f to the tilt angular velocity with respect to the trunnion axial displacement is expressed by the equation (2) as described above. This sensitivity f is a component f depending on the rotational speed of the
[0082]
[Expression 20]
Formula (14)
[0083]
[Expression 21]
Formula (15)
[0084]
[Expression 22]
Formula (16)
[0085]
F up to the tilt angle φ at the current time2The integral value q (φ) relating to the reciprocal tilt angle of (φ) is obtained using the following equation.
[0086]
[Expression 23]
Formula (17)
[0087]
f1In the calculation means 103, the rotational speed ω of the
[0088]
In the new state quantity calculation means 105, the integral values q (φ) and f1(Ωod0) And the new state quantity ζ1 is calculated using the following equation.
[0089]
[Expression 24]
Formula (18)
[0090]
In the state observation filter means 106, a new state quantity ζ1And a step motor displacement command value u, which is an output of a
[0091]
[Expression 25]
Formula (19)
[0092]
[Equation 26]
[0093]
HaIs the observer gain. The observer gain is set as follows.
[0094]
[Expression 27]
Formula (20)
[0095]
Where pTen, P11Is a constant related to the pole of the observer,
s = -pTen, S = -p11
It becomes.
[0096]
The derivation and characteristics of the
[0097]
The dynamic characteristic of the
[0098]
[Expression 28]
Formula (21)
[0099]
From the equations (21), (1), and (14), the following equation is obtained.
[0100]
[Expression 29]
Formula (22)
[0101]
Thus, using the state conversion shown in Equation (18), from Equation (22) and Equation (1), TCVT dynamics
[0102]
[30]
[0103]
Can be rewritten as:
[0104]
[31]
Formula (23)
[0105]
Summarizing the above state conversion, the sensitivity f (φ, ω) of the tilt angular velocity with respect to the trunnion displacement y is expressed by the equation (14).od) To the component f depending on the tilt angle represented by the equation (16)2(Φ) is separated from other parts, and the new state quantity is set as represented by the equation (18). At this time, the differential value of the new state quantity is expressed by Expression (21), and the coefficient of the tilt angular velocity φ ′ is expressed by f (φ, ω expressed by Expression (14).od0).
[0106]
As a result, the time-varying term of the coefficient of the trunnion displacement y is canceled from Expression (1) and Expression (21), and the coefficient becomes a constant value as shown in Expression (22). Then, considering the dynamic characteristics of the TCVT as a new state quantity instead of the state quantity of the tilt angle, the transition matrix A is expressed as shown in Expression (23).aIs a constant.
[0107]
The observer shown in Expression (19) was designed for the rewritten TCVT10.
[0108]
Where the state quantity ζaEstimated value error eaThe
[0109]
[Expression 32]
Formula (24)
[0110]
As a result, subtracting equation (19) from equation (23) yields the following equation:
[0111]
[Expression 33]
Formula (25)
[0112]
Equation (25) gives the error eaIt is an error equation showing the attenuation characteristic of Aa, Ha, CaSince both are constant matrices, (Aa-HaCa) Is also a constant matrix. This is the error eaIndicates that it decays exponentially.
[0113]
Therefore, the
[0114]
Furthermore, the speed change actuator command value of the present invention is the speed change actuator displacement, and the
[0115]
Next, the
[0116]
Hereinafter, the shift control calculation performed by the control device will be described using the flowchart shown in FIG. This shift control calculation is executed at a predetermined control cycle, for example, every 10 ms.
[0117]
In steps S100 and S101, the TCVT input / output disk rotational speed ωid, Ωco0Is read.
[0118]
In step S102, first, the tilt angle detecting means 100 performs the TCVT input disk rotational speed ω.idAnd TCVT output disk rotation speed detection value ωod0From the above, the gear ratio icIs calculated. Next, the gear ratio i shown in FIG.cAnd a map showing the relationship between the tilt angle φ and the gear ratio icFrom the above, the tilt angle φ is calculated.
[0119]
In step S103, the differentiating means 102 detects the TCVT output disk rotational speed detection value ω.od0Is input, the low-pass filter operation represented by Expression (13) is performed, and the TCVT output disk rotational speed estimated value ωodAnd TCVT output disc rotational speed differential value ω ′odAnd calculate.
[0120]
In step S104, the integral value calculation means 104 uses the equation (15) to calculate the TCVT output disk rotational speed estimated value ω.odTo the component f depending on the rotational speed1Is calculated.
[0121]
In step S105, similarly, the integral value calculation means 104 inputs the tilt angle φ and calculates the integral value q (φ) using the equation (17). Alternatively, the integral value q (φ) is obtained from the tilt angle φ by using the map shown in FIG. 6 obtained by calculating Equation (17) in advance.
[0122]
In step S106, the new state quantity calculation means 105 uses the equation (18) to calculate f1And q (φ), the new state quantity ζ1Is calculated.
[0123]
In step S107, the state observation filter means 106 performs the tilt angle φ, the step number command value u, and the new state quantity ζ.1And TCVT output disk rotational speed estimate ωodAnd TCVT output disc rotational speed differential value ω ′odAnd the trunnion displacement estimated value y expressed by the equation (19)tIs calculated.
[0124]
In step S108, the
[0125]
Therefore, the
[0126]
Next, a second embodiment will be described.
[0127]
The
[0128]
In the state observation filter means 106 of the second embodiment, the new state quantity ζ1, the step motor drive speed command value v and the tilt angle φ, which are the outputs of the command value calculator, which will be described later, are input, and the following equation is used. The new state quantity ζ1, the trunnion displacement y, and the step number u of the step motor are calculated using the filter.
[0129]
[Expression 34]
Formula (26)
[0130]
[Expression 35]
[0131]
Where HbIs the observer gain. The observer gain is set as follows.
[0132]
[Expression 36]
Formula (27)
[0133]
Where pb0, Pb1, Pb2Is a constant related to the pole of the observer,
s = -pb0, S = -pb1, S = -pb2
It becomes.
[0134]
The derivation and characteristics of this state estimation means are shown below. With the step motor step number u as input and the tilt angle φ and trunnion displacement y as state quantities, the dynamic characteristics of TCVT are expressed by equation (1). Further, the relationship between the step motor driving speed v and the step motor step number u is expressed by the following equation.
[0135]
[Expression 37]
Formula (28)
[0136]
Similar to the first embodiment, using the state conversion shown in equation (18), the dynamic characteristics of the
[0137]
[Formula 38]
[0138]
Can be rewritten as:
[0139]
[39]
Formula (29)
[0140]
Thus, the transition matrix A is obtained by the same state conversion as in the first embodiment.bIs a constant. The observer shown in Expression (26) was designed for the rewritten TCVT10. Estimated value error e of ζbbThe
[0141]
[Formula 40]
Formula (30)
[0142]
As a result, subtracting equation (26) from equation (29) yields the following equation:
[0143]
[Expression 41]
Formula (31)
[0144]
Equation (26) gives the error ebIt is an error equation showing the attenuation characteristic of Ab, Hb, CbSince both are constant matrices, (Ab-HbCb) Is also a constant matrix. This is the error ebIndicates that it decays exponentially.
[0145]
The
[0146]
[Expression 42]
Formula (32)
[0147]
[Expression 43]
Formula (33)
[0148]
Here, f ′ is expressed by the following equation.
[0149]
(44)
Formula (34)
[0150]
Next, the step motor drive speed command value v is calculated using the control law shown in the following equation.
[0151]
[Equation 45]
Formula (35)
[0152]
[Equation 46]
Formula (36)
[0153]
Where φ*Is the target tilt angle, ωnIs a constant that determines the response speed of the tilt angle with respect to the target tilt angle, k is a switching gain, and ε is a constant that makes the switching gain continuous in the vicinity of zero of the control error σ. Equations (35) and (36) are sliding mode control laws. According to the equation (35), the stepping motor is driven in the direction in which the control error becomes zero. In the equation (36), when σ = 0, the following equation is obtained.
[0154]
[Equation 47]
Formula (37)
[0155]
Thus, when σ = 0, the tilt angle responds with the dynamic characteristic shown in the equation (37) with respect to the target tilt angle.
[0156]
Further, consider a case where the stepping motor does not operate according to the command value and step-out occurs. Disturbance d caused by step-out in equation (29)uThe TCVT model to which is added is expressed by the following equation.
[0157]
[Formula 48]
Formula (38)
[0158]
[Equation 49]
[0159]
Subtracting equation (26) from equation (38) yields the following error equation:
[0160]
[Equation 50]
Formula (39)
[0161]
Assume that the error equation shown in equation (39) is stable. At this time, error ebSteady value eb0Is e ′ in equation (39).b= 0 and is represented by the following formula.
[0162]
[Equation 51]
Formula (40)
[0163]
Thus, even if a step-out occurs, the steady error of the tilt angle φ and the steady error of the trunnion displacement y are zero. Further, from equation (40), the estimated value u of the number of stepstIs
[0164]
[Formula 52]
[0165]
It becomes. This indicates that the actual number of steps including disturbance can be estimated, and there is no estimation error of the number of steps.
[0166]
Therefore, according to the present embodiment, the speed change actuator command value is the speed change actuator drive speed command value, and the
[0167]
Next, consider the drive state determination means 120 as shown in FIG. The drive state determination means 120 calculates the integrated value of the step motor drive pulse as u.ItAnd uItAnd the output u of the state observation filter meanstCompare And, for example, at steady state utAnd uItIf it is greatly different, it is determined that the step is out of step. For example, utChange and uItThe step motor abnormality (step motor sticking or the like) is determined by comparing with the change in the above.
[0168]
Hereinafter, the shift control calculation performed by the
[0169]
Steps S100 to S106 are the same as those in the first embodiment, and a description thereof will be omitted.
[0170]
In step S107, the tilt angle φ, the step motor drive speed command value v, and ζ1And TCVT output disk rotational speed estimate ωodAnd TCVT output disc rotational speed differential value ω ′odAnd the
[0171]
In step S108, the trunnion displacement estimated value ytAnd step motor step number estimate utAnd tilt angle φ and target value φ*And the control law expressed by the equations (32) to (36) is calculated to calculate the step motor drive speed command value v.
[0172]
FIG. 8 shows the result of a computer simulation performed using this
[0173]
Therefore, the
[0174]
The present invention is not limited to the above-described embodiment, and it is obvious that various modifications can be made within the scope of the technical idea of the present invention.
[Brief description of the drawings]
FIG. 1 is a schematic configuration diagram of a continuously variable transmission.
FIG. 2 is a configuration diagram of a shift control hydraulic system of a continuously variable transmission.
FIG. 3 is a configuration diagram of a control system including a control device.
FIG. 4 is a configuration diagram of a state observer.
FIG. 5 is a relationship diagram between a tilt angle and a gear ratio.
FIG. 6: Tilt angle and f2It is a relationship figure with the integral value of the reciprocal number of.
FIG. 7 is a flowchart of a shift control calculation of the control device.
FIG. 8 is a simulation result using the present invention by a computer.
FIG. 9 is a configuration diagram of a control system including drive state determination means.
[Explanation of symbols]
10 Toroidal type continuously variable transmission
52 step motor
56 Shift control valve
100 Tilt angle detector
101 Rotational speed detection means
102 Differentiation means
103 f1Calculation means
104 Integral value calculation means
105 New state quantity calculation means
106 Filter means for state observation
110 State observer
200 Command value calculator
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002168347AJP4214720B2 (en) | 2002-06-10 | 2002-06-10 | Toroidal continuously variable transmission |
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002168347AJP4214720B2 (en) | 2002-06-10 | 2002-06-10 | Toroidal continuously variable transmission |
| Publication Number | Publication Date |
|---|---|
| JP2004011834A JP2004011834A (en) | 2004-01-15 |
| JP4214720B2true JP4214720B2 (en) | 2009-01-28 |
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002168347AExpired - Fee RelatedJP4214720B2 (en) | 2002-06-10 | 2002-06-10 | Toroidal continuously variable transmission |
| Country | Link |
|---|---|
| JP (1) | JP4214720B2 (en) |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7011600B2 (en) | 2003-02-28 | 2006-03-14 | Fallbrook Technologies Inc. | Continuously variable transmission |
| JP4649907B2 (en)* | 2004-08-05 | 2011-03-16 | 株式会社豊田中央研究所 | Control device for continuously variable transmission |
| WO2006041718A2 (en) | 2004-10-05 | 2006-04-20 | Fallbrook Technologies, Inc. | Continuously variable transmission |
| WO2007070167A2 (en) | 2005-10-28 | 2007-06-21 | Fallbrook Technologies Inc. | Electromotive drives |
| PL1954959T3 (en) | 2005-11-22 | 2013-10-31 | Fallbrook Ip Co Llc | Continuously variable transmission |
| CN102221073B (en) | 2005-12-09 | 2013-03-27 | 福博科技术公司 | Continuously variable transmission |
| EP1811202A1 (en)* | 2005-12-30 | 2007-07-25 | Fallbrook Technologies, Inc. | A continuously variable gear transmission |
| CN102269055B (en) | 2006-06-26 | 2013-08-28 | 福博科技术公司 | Continuously variable transmission |
| EP2125469A2 (en) | 2007-02-01 | 2009-12-02 | Fallbrook Technologies Inc. | System and methods for control of transmission and/or prime mover |
| US20100093479A1 (en) | 2007-02-12 | 2010-04-15 | Fallbrook Technologies Inc. | Continuously variable transmissions and methods therefor |
| TWI461615B (en) | 2007-02-16 | 2014-11-21 | Fallbrook Ip Co Llc | Infinitely variable transmissions, continuously variable transmissions, methods, assemblies, subassemblies, and components therefor |
| EP2142826B1 (en) | 2007-04-24 | 2015-10-28 | Fallbrook Intellectual Property Company LLC | Electric traction drives |
| US8641577B2 (en) | 2007-06-11 | 2014-02-04 | Fallbrook Intellectual Property Company Llc | Continuously variable transmission |
| CN103697120B (en) | 2007-07-05 | 2017-04-12 | 福博科技术公司 | Continuously variable transmission |
| CN103939602B (en) | 2007-11-16 | 2016-12-07 | 福博科知识产权有限责任公司 | Controllers for variable speed drives |
| US8321097B2 (en) | 2007-12-21 | 2012-11-27 | Fallbrook Intellectual Property Company Llc | Automatic transmissions and methods therefor |
| US8313405B2 (en) | 2008-02-29 | 2012-11-20 | Fallbrook Intellectual Property Company Llc | Continuously and/or infinitely variable transmissions and methods therefor |
| US8317651B2 (en) | 2008-05-07 | 2012-11-27 | Fallbrook Intellectual Property Company Llc | Assemblies and methods for clamping force generation |
| CN102112778B (en) | 2008-06-06 | 2013-10-16 | 福博科技术公司 | Infinitely variable transmission, continuously variable transmission, methods, assemblies, subassemblies and components therefor |
| EP2304272B1 (en) | 2008-06-23 | 2017-03-08 | Fallbrook Intellectual Property Company LLC | Continuously variable transmission |
| CA2732668C (en) | 2008-08-05 | 2017-11-14 | Fallbrook Technologies Inc. | Methods for control of transmission and prime mover |
| US8469856B2 (en)* | 2008-08-26 | 2013-06-25 | Fallbrook Intellectual Property Company Llc | Continuously variable transmission |
| US8167759B2 (en) | 2008-10-14 | 2012-05-01 | Fallbrook Technologies Inc. | Continuously variable transmission |
| ES2439647T3 (en) | 2009-04-16 | 2014-01-24 | Fallbrook Intellectual Property Company Llc | Stator set and speed change mechanism for a continuously variable transmission |
| JP5346834B2 (en)* | 2010-02-08 | 2013-11-20 | 本田技研工業株式会社 | Control device for automatic transmission |
| US8512195B2 (en) | 2010-03-03 | 2013-08-20 | Fallbrook Intellectual Property Company Llc | Infinitely variable transmissions, continuously variable transmissions, methods, assemblies, subassemblies, and components therefor |
| US8888643B2 (en) | 2010-11-10 | 2014-11-18 | Fallbrook Intellectual Property Company Llc | Continuously variable transmission |
| CN104302949B (en) | 2012-01-23 | 2017-05-03 | 福博科知识产权有限责任公司 | Infinitely variable continuously variable transmission, continuously variable continuously variable transmission, method, assembly, subassembly, and parts thereof |
| KR102433297B1 (en) | 2013-04-19 | 2022-08-16 | 폴브룩 인텔렉츄얼 프로퍼티 컴퍼니 엘엘씨 | Continuously variable transmission |
| US10047861B2 (en) | 2016-01-15 | 2018-08-14 | Fallbrook Intellectual Property Company Llc | Systems and methods for controlling rollback in continuously variable transmissions |
| WO2017154037A1 (en)* | 2016-03-08 | 2017-09-14 | 川崎重工業株式会社 | Shift control device |
| KR102364407B1 (en) | 2016-03-18 | 2022-02-16 | 폴브룩 인텔렉츄얼 프로퍼티 컴퍼니 엘엘씨 | continuously variable transmission system and method |
| US10023266B2 (en) | 2016-05-11 | 2018-07-17 | Fallbrook Intellectual Property Company Llc | Systems and methods for automatic configuration and automatic calibration of continuously variable transmissions and bicycles having continuously variable transmissions |
| US11215268B2 (en) | 2018-11-06 | 2022-01-04 | Fallbrook Intellectual Property Company Llc | Continuously variable transmissions, synchronous shifting, twin countershafts and methods for control of same |
| WO2020176392A1 (en) | 2019-02-26 | 2020-09-03 | Fallbrook Intellectual Property Company Llc | Reversible variable drives and systems and methods for control in forward and reverse directions |
| Publication number | Publication date |
|---|---|
| JP2004011834A (en) | 2004-01-15 |
| Publication | Publication Date | Title |
|---|---|---|
| JP4214720B2 (en) | Toroidal continuously variable transmission | |
| JP3680739B2 (en) | Shift control device for continuously variable transmission | |
| JP4284905B2 (en) | Shift control device for continuously variable transmission | |
| US6931316B2 (en) | Toroidal continuously variable transmission control apparatus | |
| US6909953B2 (en) | Shift control of continuously-variable transmission | |
| JP3460547B2 (en) | Transmission control device for continuously variable transmission | |
| JP3161326B2 (en) | Control device for continuously variable automatic transmission | |
| JP2705383B2 (en) | Transmission control device for friction wheel type continuously variable transmission | |
| JP3531362B2 (en) | Fail safe device for toroidal type continuously variable transmission | |
| JP3794295B2 (en) | Shift control device for toroidal type continuously variable transmission | |
| JP3520593B2 (en) | Transmission control device for toroidal type continuously variable transmission | |
| JP3711619B2 (en) | Toroidal continuously variable transmission | |
| JP4304745B2 (en) | Powertrain control device | |
| JP3301365B2 (en) | Transmission control device for toroidal type continuously variable transmission | |
| JP2004068936A (en) | Gear ratio control device for continuously variable transmission | |
| JP3736471B2 (en) | Toroidal continuously variable transmission | |
| JP4253891B2 (en) | Powertrain control device | |
| JP3960165B2 (en) | Shift control device for continuously variable transmission | |
| JP4085679B2 (en) | Gear ratio control device for continuously variable transmission mechanism | |
| JP4007078B2 (en) | Shift control device for continuously variable transmission | |
| JPH09217806A (en) | Toroidal type continuously variable transmission | |
| JP2004011728A (en) | Transmission control device for continuously variable transmission | |
| JP4013690B2 (en) | Shift control device for continuously variable transmission | |
| JP2004324691A (en) | Transmission control device for continuously variable transmission | |
| JP2009063021A (en) | Continuously variable transmission |
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination | Free format text:JAPANESE INTERMEDIATE CODE: A621 Effective date:20050526 | |
| 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:20081014 | |
| A01 | Written decision to grant a patent or to grant a registration (utility model) | Free format text:JAPANESE INTERMEDIATE CODE: A01 | |
| A61 | First payment of annual fees (during grant procedure) | Free format text:JAPANESE INTERMEDIATE CODE: A61 Effective date:20081027 | |
| R150 | Certificate of patent or registration of utility model | Free format text:JAPANESE INTERMEDIATE CODE: R150 | |
| FPAY | Renewal fee payment (event date is renewal date of database) | Free format text:PAYMENT UNTIL: 20111114 Year of fee payment:3 | |
| LAPS | Cancellation because of no payment of annual fees |