JP2008045651A - Shift simulation device, shift simulation program, and automobile - Google Patents
Shift simulation device, shift simulation program, and automobile Download PDFInfo
- Publication number
- JP2008045651A JP2008045651A JP2006221312A JP2006221312A JP2008045651A JP 2008045651 A JP2008045651 A JP 2008045651A JP 2006221312 A JP2006221312 A JP 2006221312A JP 2006221312 A JP2006221312 A JP 2006221312A JP 2008045651 A JP2008045651 A JP 2008045651A
- Authority
- JP
- Japan
- Prior art keywords
- command value
- input command
- value
- shift
- simulation
- 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
Images
Landscapes
- Control Of Transmission Device (AREA)
- Feedback Control In General (AREA)
Abstract
Description
本発明は、モデルを利用して自動車用の自動変速機システムの変速動作のシミュレーションを行う変速シミュレーションに関する。 The present invention relates to a shift simulation for simulating a shift operation of an automatic transmission system for an automobile using a model.
従来より、多くの車両に自動変速機が搭載されており、この自動変速機における変速動作は、出力軸トルクおよびタービン回転数が目標通りに推移するように制御される。そして、この制御指令値(入力指令値)を得るためには、各種条件による実車走行試験を行う必要がある。 Conventionally, an automatic transmission is mounted on many vehicles, and the speed change operation in the automatic transmission is controlled so that the output shaft torque and the turbine rotational speed change as intended. In order to obtain this control command value (input command value), it is necessary to perform an actual vehicle running test under various conditions.
しかし、このような走行試験は、時間と作業量が膨大となるため、これに代わるものとして自動変速機のモデルを利用したシミュレーションなども行われている。 However, since such a running test requires a large amount of time and work, a simulation using a model of an automatic transmission is also performed as an alternative.
例えば、特許文献1では、クラッチ・トゥ・クラッチ変速において、入力トルクTin、係合側摩擦要素のトルク容量T1、解放側摩擦要素のトルク容量T2、TinとT1とのトルク比を表すA、Tinと変速前の変速段において釣り合うT2とのトルク比を表すB、補正項Cにより設定される運動方程式Tin=AT1+BT2+Cに従い、解放側・係合側摩擦要素のどちらか一方の油圧指令を求めることが示されている。
For example,
なお、変速制御についてファジィ推論を用いるものとして、特許文献2、特許文献3等の提案がある。特許文献4には、ファジィ推論を用いた自動車エンジンの空燃比制御について記載されている。また、品質工学において、設計諸元を決定するに当たり、品質をばらつかせる誤差因子(ノイズ)を減衰させる水準を見つけ出し、ノイズに強い安定設計を行う手法があり、非特許文献1等に示されており、特許文献5,6にはこの品質工学の考えを利用したロボット制御について記載されている。
In addition, there exist proposals of
前記特許文献1では、変速時の解放側または係合側摩擦要素のどちらか一方の油圧指令は、予め適合作業によって求めておく必要がある。また、性能バラツキへの対応は、オンボードでのフィードバック制御や学習制御(これらは補正項に含まれる)によって実行されており、フィードフォワード制御によって性能バラツキに対応する試みはなされていない。
In
また、特許文献2〜4では、ファジィ推論を用いてはいるが、対象となる装置のパラメータについての適合作業は従来通り必要であり、抜本的な工数削減にはつながらない。
Further, in
さらに、特許文献5,6は、フィードフォワード制御ではなく、自動変速機制御には適用しにくく、さらに非特許文献1は対象とする製品の機構・形状・材質などの設計諸元の決定に関するものであり、動的なシステムを対象としていない。
Furthermore,
本発明は、モデルを利用して自動車用の自動変速機システムの変速動作のシミュレーションを行う変速シミュレーション装置であって、前記システムにおいて変速動作を行う際の入力指令値もしくは変速動作を制御する装置の構成要素の特性値について、変動を付与してシミュレーションし、このシミュレーションの結果であるシステムの出力値を、前記変動に対する影響の大きさに基づいて評価し、この評価結果から、前記変動の出力値への影響が要求仕様を満たす目標出力値を決定し、決定された目標出力値に一致する出力を得るための入力指令値を学習計算によって求めることで、要求仕様を満たす入力指令値を求めることを特徴とする。 The present invention relates to a shift simulation device that simulates a shift operation of an automatic transmission system for an automobile using a model, and is an apparatus that controls an input command value or a shift operation when performing a shift operation in the system. The characteristic value of the component is simulated with a variation, and the output value of the system as a result of the simulation is evaluated based on the magnitude of the influence on the variation. From the evaluation result, the output value of the variation is evaluated. The target output value that satisfies the required specification is determined by learning calculation to determine the target output value that satisfies the required specification and obtains the input command value to obtain the output that matches the determined target output value. It is characterized by.
また、本発明は、モデルを利用して自動車用の自動変速機システムの変速動作のシミュレーションを行う変速シミュレーション装置であって、
(a)前記システムにおいて変速動作を行う際の入力指令値もしくは変速動作を制御する装置の構成要素の特性値について、変動を付与してシミュレーションし、
(b)このシミュレーションの結果であるシステムの出力値を、前記変動に対する影響の大きさに基づいて評価し、
(c)評価結果に基づき、ファジィ推論によって、目標出力値を修正し、
(d)修正した目標出力値を用いて、シミュレーションを行い、その学習計算により新たな入力指令値を算出し、
(e)算出した入力指令値により前記入力指令値を更新し、前記(a)に戻り、前記(a)〜(d)の処理を(b)の評価が要求仕様を満たすまで繰り返し、要求仕様を満たす入力指令値を求めることを特徴とする。
Further, the present invention is a shift simulation device that simulates a shift operation of an automatic transmission system for an automobile using a model,
(A) A simulation is performed by giving a variation to an input command value or a characteristic value of a component of a device that controls the shift operation when performing a shift operation in the system;
(B) evaluating the output value of the system, which is the result of this simulation, based on the magnitude of the effect on the variation;
(C) Based on the evaluation result, the target output value is corrected by fuzzy inference,
(D) A simulation is performed using the corrected target output value, a new input command value is calculated by the learning calculation,
(E) Update the input command value with the calculated input command value, return to (a), and repeat the processes (a) to (d) until the evaluation in (b) satisfies the required specification. An input command value satisfying the above is obtained.
また、本発明は、モデルを利用して自動車用の自動変速機システムの変速動作のシミュレーションを行う変速シミュレーション装置であって、
(a)目標出力値における複数の制御因子についてそれぞれ複数の水準を設定し、
(b)設定された複数の制御因子についての複数の水準の組み合わせで決定される複数の目標出力値を設定し、
(c)設定された複数の目標出力値を用い、前記システムにおいて変速動作を行う際の入力指令値もしくは変速動作を制御する装置の構成要素の特性値について変動を付与してそれぞれシミュレーションし、
(d)このシミュレーションの結果であるシステムの複数の出力値を、前記変動に対する影響の大きさに基づいて評価し、得られた評価結果が要求仕様を満たす目標出力値を決定し、
(e)決定された目標出力値を用いて、学習計算によって要求仕様を満たす入力指令値を求める、ことを特徴とする。
Further, the present invention is a shift simulation device that simulates a shift operation of an automatic transmission system for an automobile using a model,
(A) A plurality of levels are set for a plurality of control factors at the target output value,
(B) setting a plurality of target output values determined by a combination of a plurality of levels for a plurality of set control factors;
(C) Using a plurality of set target output values, simulating each of the input command value or the characteristic value of the component of the device that controls the shift operation by performing a shift operation in the system, respectively,
(D) evaluating a plurality of output values of the system as a result of the simulation based on the magnitude of the influence on the fluctuation, and determining a target output value that the obtained evaluation result satisfies the required specifications;
(E) Using the determined target output value, an input command value satisfying the required specifications is obtained by learning calculation.
また、前記(b)において、複数の制御因子の複数の水準の組み合わせは、直交表を用いて、決定し、前記(c)は決定された全ての組み合わせについて行われることが好適である。 In (b), it is preferable that combinations of a plurality of levels of a plurality of control factors are determined using an orthogonal table, and (c) is performed for all the determined combinations.
さらに、本発明は、上述のような変速シミュレーション装置をコンピュータを用いて実現するための変速プログラムに関する。 Furthermore, the present invention relates to a shift program for realizing the shift simulation apparatus as described above using a computer.
本発明によれば、入力指令値もしくは変速動作を制御する装置の構成要素の特性値を変動させて自動変速機システムの変速動作のシミュレーションを行い、対応する出力値がどう変動するかを検出する。そして、入力指令値もしくは変速動作を制御する装置の構成要素の特性値の変動に対し頑強な目標出力値を算出し、得られた目標出力値に基づいて入力指令値を得る。このため、変動に強い入力指令値を得ることができ、適切な自動変速機の制御値を得ることができる。 According to the present invention, the input command value or the characteristic value of the component of the device that controls the shift operation is varied to simulate the shift operation of the automatic transmission system and detect how the corresponding output value varies. . Then, a target output value that is robust against fluctuations in the input command value or the characteristic value of the component of the device that controls the shift operation is calculated, and the input command value is obtained based on the obtained target output value. Therefore, an input command value that is resistant to fluctuations can be obtained, and an appropriate control value for the automatic transmission can be obtained.
以下、本発明の実施形態について、図面に基づいて説明する。 Hereinafter, embodiments of the present invention will be described with reference to the drawings.
本実施形態は、複数の摩擦係合装置の係合・解放状態に応じて変速比が設定される自動変速機の入力指令値についてシミュレーションにより適合化するものである。特に、目標出力値を適切に設定することで、入力指令値もしくは変速動作を制御する制御装置の構成要素の特性値のバラツキに強い入力指令値を得る。 In the present embodiment, an input command value of an automatic transmission in which a gear ratio is set according to the engagement / release states of a plurality of friction engagement devices is adapted by simulation. In particular, by appropriately setting the target output value, an input command value that is resistant to variations in the input command value or the characteristic values of the components of the control device that controls the shift operation is obtained.
この適合対象としている自動変速機は、複数の遊星歯車と、複数の油圧アクチュエータによって摩擦係合される油圧式摩擦係合装置によって構成されており、この油圧アクチュエータは、車載コンピュータによって電子制御される。なお、シミュレーションは、コンピュータを用いて行われ、その結果で得られた入力指令値が車載メモリ等に記憶される。そして、車載コンピュータが車載メモリに記憶されている入力指令値を用いて変速時の動作を制御する。 The automatic transmission to be adapted is configured by a plurality of planetary gears and a hydraulic friction engagement device that is frictionally engaged by a plurality of hydraulic actuators, and the hydraulic actuators are electronically controlled by an in-vehicle computer. . The simulation is performed using a computer, and the input command value obtained as a result is stored in an in-vehicle memory or the like. Then, the in-vehicle computer controls the operation at the time of shifting using the input command value stored in the in-vehicle memory.
「モデルによるシミュレーション」
図8に、CAE(computer aided engineering)モデル上での適合作業の概要を示す。CAEモデルは、エンジン・トルクコンバータモデル10と、自動変速機の遊星歯車によって構成されるギヤトレーンモデル12と、解放側ブレーキモデル14、係合側ブレーキモデル16とを含み、これらモデルへ入力される入力指令値の修正量が修正量計算部18において算出される。この修正量計算部18は、実際の出力値と目標出力値の比較に基づく学習計算によって修正量を算出する。
"Simulation by model"
FIG. 8 shows an outline of the adaptation work on a CAE (computer aided engineering) model. The CAE model includes an engine /
エンジン遅角量についての初期値は加算器20に供給される。ここには、修正量計算部18からエンジン遅角量についての修正量が供給されており、両者が加算されることでエンジン遅角量が修正される。修正されたエンジン遅角量はエンジン・トルクコンバータモデル10に供給される。従って、供給されたエンジン遅角量に応じてエンジン出力が制御されトルクコンバータからの出力タービントルクが決定される。
An initial value for the engine retard amount is supplied to the
また、解放側油圧指令の初期値が加算器22に供給され、ここには修正量計算部18からの解放側油圧指令の修正量が供給されている。従って、加算器22において修正された解放側油圧指令が得られ、これが解放側ブレーキモデル14に供給され、この解放側ブレーキモデル14から対応する解放側トルク容量T1が出力される。
The initial value of the release side hydraulic pressure command is supplied to the
さらに、係合側油圧指令の初期値が加算器24に供給され、ここには修正量計算部18からの係合側油圧指令の修正量が供給されている。従って、加算器24において修正された係合側油圧指令が得られ、これが係合側ブレーキモデル16に供給され、係合側ブレーキモデル16から対応する係合側トルク容量T2が出力される。
Further, the initial value of the engagement side hydraulic pressure command is supplied to the adder 24, and here, the correction amount of the engagement side hydraulic pressure command from the correction
エンジン・トルクコンバータモデル10から出力されたタービントルク、解放側ブレーキモデル14から出力された解放側トルク容量T1、および係合側ブレーキモデル16から対応する係合側トルク容量T2がギヤトレーンモデル12に供給される。これによって、このギヤトレーンモデル12から入力に応じた出力軸トルクToおよびタービン角加速度(d/dt)ωtが出力される。
The turbine torque output from the engine
そして、出力軸トルクToは加算器26に負の値として供給され、ここで目標値から減算され、得られた誤差が修正量計算部18に供給される。また、タービン角加速度(d/dt)ωtは加算器28に負の値として供給され、ここで目標値Tod、(d/dt)ωtdからそれぞれ減算され、得られた誤差が修正量計算部18に供給される。
Then, the output shaft torque T o is supplied as a negative value to the
修正量計算部18は、供給される誤差に基づいて、エンジン遅角量、解放側油圧指令、係合側油圧指令についての修正量を演算算出し、これらを加算器20,22,24にそれぞれ供給する。
Based on the supplied error, the correction
「ファジィ推論を用いる実施形態」
本実施形態では、自動車に搭載する自動変速機を対象とした数値計算シミュレーションを通じて、理想的な変速動作を実現するためのフィードフォワード指令値を導き出す。自動変速機の変速段については、多段化が進んであり、その各段間の変速時において、目標出力軸トルク波形および目標タービン角加速度波形を最適なものとして、その目標波形を実現するフィードフォワード指令値(解放側および係合側ブレーキ油圧・エンジン遅角量)を学習計算により導き出す。
"Embodiment using fuzzy reasoning"
In the present embodiment, a feedforward command value for realizing an ideal shift operation is derived through a numerical calculation simulation for an automatic transmission mounted on an automobile. The shift stages of the automatic transmission are becoming multi-stage, and the feed-forward that realizes the target waveform by optimizing the target output shaft torque waveform and the target turbine angular acceleration waveform at the time of shifting between the stages. The command value (release side and engagement side brake hydraulic pressure and engine retard amount) is derived by learning calculation.
ここで、実際の自動変速機においては、部品の個体バラツキや、走行状況などの変化によって、自動変速機の機械的な特性には変動が生じる。この特性変動を考慮せず、自動変速機が設計時の定数に従う(この場合をノミナル状態と呼ぶ)と仮定し、入力指令値を導出した場合には、出力結果が変動し、その変動幅が許容範囲を大きく逸脱する可能性もある。 Here, in an actual automatic transmission, the mechanical characteristics of the automatic transmission vary due to variations in individual parts, changes in traveling conditions, and the like. Without considering this characteristic variation, assuming that the automatic transmission follows the design constants (this case is called the nominal state) and deriving the input command value, the output result varies, and the fluctuation range is There is also a possibility that the tolerance is greatly deviated.
そこで、本実施形態では、入力指令値を学習計算によって導き出す際に用いる目標出力値(目標波形)について、入力指令値もしくは変速動作を制御する装置の構成要素の特性値の変動に強いものにする。すなわち、入力指令値もしくは変速動作を制御する装置の構成要素の特性値について変動があった場合にも出力結果の変動が要求される許容範囲内となる目標出力値を決定する。 Therefore, in this embodiment, the target output value (target waveform) used when the input command value is derived by learning calculation is made resistant to fluctuations in the input command value or the characteristic value of the component of the device that controls the shift operation. . That is, a target output value that falls within an allowable range in which a change in the output result is required even when there is a change in the input command value or the characteristic value of the component of the device that controls the shift operation is determined.
このために、本実施形態では、コンピュータ上での数値シミュレーションにより、目標波形(本実施形態の場合には目標出力軸トルク波形)を設定し、さらにこの目標波形を実現するための入力指令値を導出する。特に、目標波形の決定には、ファジィ推論を活用する。また、これらの処理はコンピュータを用い、自動的に行う。 Therefore, in this embodiment, a target waveform (target output shaft torque waveform in this embodiment) is set by numerical simulation on a computer, and an input command value for realizing this target waveform is set. To derive. In particular, fuzzy reasoning is used to determine the target waveform. These processes are automatically performed using a computer.
本実施形態における、自動変速機の入力指令値を求める適合作業の流れを図1に示す。まず、ステップS1では、自動変速機に入力する入力指令値の初期値を設定する。この初期指令値は、単に所定の変速段の切り換え作業を行うことができればよい。このため、図2の(i)に示すようなランプ入力などのごく簡単な入力波形を与えればよい。図2(i)では、ブレーキ掴み換えの入力指令値(解放側油圧指令および係合側油圧指令)についての初期値の例を示してある。なお、上述の図8によれば、入力指令値は、係合側油圧指令、解放側油圧指令、遅角量の3つである。 FIG. 1 shows a flow of adapting work for obtaining the input command value of the automatic transmission in this embodiment. First, in step S1, an initial value of an input command value input to the automatic transmission is set. The initial command value only needs to be able to perform a predetermined shift stage switching operation. Therefore, a very simple input waveform such as a ramp input as shown in FIG. FIG. 2 (i) shows an example of initial values for the input command values (release side hydraulic pressure command and engagement side hydraulic pressure command) for brake grip replacement. According to FIG. 8 described above, there are three input command values: an engagement side hydraulic pressure command, a release side hydraulic pressure command, and a retard amount.
次に、ステップS2では、入力指令値もしくは変速動作を制御する装置の構成要素の特性値に意図的に変動を加えたものを入力として、変速シミュレーションを実行する。図2(ii)の例では、係合側の油圧指令値について、+バラツキのものと、−バラツキのものを加えた3入力とし、それぞれについてシミュレーションする。また、図3では、図2(ii)に示したように、係合側油圧波形に±のバラツキを加えた状態で変速シミュレーションを行ったときの3つの結果(出力軸トルク)を示している。この図3のように、入力指令値の変動は、出力結果に大きな差異をもたらす。 Next, in step S2, a shift simulation is executed with an input command value or a characteristic value of a component of a device that controls the shift operation as an input. In the example of FIG. 2 (ii), the hydraulic pressure command value on the engagement side is assumed to be three inputs including a plus variation and a minus variation, and each is simulated. Further, FIG. 3 shows three results (output shaft torque) when a shift simulation is performed in a state where ± variation is added to the engagement-side hydraulic waveform as shown in FIG. 2 (ii). . As shown in FIG. 3, the change in the input command value causes a large difference in the output result.
ステップS3では、ステップS2の結果から、入力指令値もしくは変速動作を制御する装置の構成要素の特性値に変動が加えられた場合においても、変速シミュレーションの出力結果について、予め定められている複数の評価指標が基準値(要求仕様)を満足しているかどうかを判断する。評価指標には、図2の(ii)、(iii)に示してあるように、変速動作中の出力軸トルクの増減、出力軸トルクの最大値・最小値、変速動作の開始から変速が完了するまでの変速時間、変速動作中に摩擦材が回転体間の相対回転を止める際に発生する熱の吸収量などが挙げられる。また、この他にも、瞬間的なショックを避けるという意図において、出力軸トルクの勾配、摩擦材の熱吸収率などが評価指標として採用される。 In step S3, a plurality of predetermined output results of the shift simulation are determined even when the input command value or the characteristic value of the component of the device that controls the shift operation is changed from the result of step S2. It is determined whether the evaluation index satisfies a standard value (required specification). As shown in (ii) and (iii) of FIG. 2, the evaluation index includes increase / decrease in output shaft torque during the shifting operation, maximum / minimum values of the output shaft torque, and shifting is completed from the start of shifting operation. And the amount of heat absorbed when the friction material stops relative rotation between the rotating bodies during the shifting operation. In addition, with the intention of avoiding an instantaneous shock, the gradient of the output shaft torque, the heat absorption rate of the friction material, and the like are employed as evaluation indexes.
これらの評価指標は、ファジィロジックの前件部のメンバシップ関数へ入力され、それぞれの指標がどのファジィ集合に、どの度合い(グレード)で属するかが判定される。前件部のメンバシップ関数の例を図4に示す。図4(a)のトルク増減の場合では、トルク増減が小さい(S)、中間(M)、大きい(B)の3つの集合に分けられる。また、図4(b)の変速時間に関しては、短い(S)、中間(M)、長い(L)の3つの集合に分けられる。 These evaluation indices are input to the membership function of the antecedent part of the fuzzy logic, and it is determined to which fuzzy set each index belongs to which degree (grade). An example of the membership function of the antecedent part is shown in FIG. In the case of torque increase / decrease in FIG. 4A, the torque increase / decrease is divided into three sets of small (S), middle (M), and large (B). 4B is divided into three sets of short (S), middle (M), and long (L).
そして、各評価指標が、望ましいファジィ集合に所属しているかによって、ステップS3における判定を行う。 Then, the determination in step S3 is performed depending on whether each evaluation index belongs to a desirable fuzzy set.
ステップS3の判定において、Noであれば、判定がYesになるような修正を行う。このために、ステップS4では、ファジィ推論により目標波形を修正する。 If the determination in step S3 is No, the correction is made so that the determination is Yes. For this purpose, in step S4, the target waveform is corrected by fuzzy inference.
すなわち、各評価指標が望ましいファジィ集合に所属していなかった場合、図2の(iii)に示すように、表1に示すようなIF−THEN形式で記述される複数の規則に従い、目標出力軸トルク波形の修正量を決定する。 That is, when each evaluation index does not belong to the desired fuzzy set, as shown in (iii) of FIG. 2, the target output axis is set according to a plurality of rules described in the IF-THEN format as shown in Table 1. Determine the amount of torque waveform correction.
後件部に関しても、図5に示すようなメンバシップ関数が用意されており、各規則の前件部の度合いから、後件部のファジィ集合の度合いがファジィ推論により求められる。そして、各規則の推論結果が集められて、最終的な推論結果を得る。この推論結果は、”良い”、”悪い”といったファジィ変数ではなく、非ファジィ化、つまり数値化されたものであり、これが目標出力軸トルク波形の修正量となる。なお、図8の例であれば、目標波形として、タービン角速度波形があり、これについても同様にして修正量が求められる。 Regarding the consequent part, membership functions as shown in FIG. 5 are prepared, and the degree of the fuzzy set of the consequent part is obtained by fuzzy inference from the degree of the antecedent part of each rule. Then, the inference results of each rule are collected to obtain a final inference result. This inference result is not a fuzzy variable such as “good” or “bad”, but is defuzzified, that is, a numerical value, and this is a correction amount of the target output shaft torque waveform. In the example of FIG. 8, there is a turbine angular velocity waveform as the target waveform, and the correction amount is obtained in the same manner.
このようにして、目標波形の修正量が求められた場合には、ステップS5において、得られた目標波形を用いて、図8のモデルにより適合作業を行い、学習計算によって修正された入力指令値を得る。このステップS5における適合作業については、後述する。 In this way, when the correction amount of the target waveform is obtained, the input command value corrected by the learning calculation is performed in step S5 by using the target waveform obtained and performing the adaptation work using the model of FIG. Get. The fitting operation in step S5 will be described later.
ステップS5において、修正された目標波形を用いて得た入力指令値が得られた場合には、新しい入力指令値を用いてステップS2からの処理を繰り返す。そして、シミュレーションを繰り返し実行しながら、ステップS3においてYesとなる修正された目標波形の修正量を求める。出力波形である出力軸トルクの目標波形についての修正の状態の一例を図6に示す。このように、初期波形について、1回目、2回目と徐々に修正が行われる。 In step S5, when the input command value obtained using the corrected target waveform is obtained, the processing from step S2 is repeated using the new input command value. Then, while repeatedly executing the simulation, the correction amount of the corrected target waveform that becomes Yes in step S3 is obtained. An example of the state of correction for the target waveform of the output shaft torque, which is an output waveform, is shown in FIG. As described above, the initial waveform is gradually corrected first time and second time.
このような処理を繰り返すことで、図7に示すように、ステップS3の判定結果がYesとなる特性変動が発生した場合においても、許容範囲内に収まる目標波形を作り出すことができる。そして、この目標波形に対応する入力指令値を得ることができる。従って、このようにして得られた入力指令値を用いて、自動変速機を制御することで、入力指令値もしくは変速動作を制御する装置の構成要素の特性値についての変動があっても、出力への影響を抑制した変速動作が達成される。 By repeating such a process, as shown in FIG. 7, even when a characteristic variation occurs in which the determination result in step S3 is Yes, a target waveform that falls within the allowable range can be created. An input command value corresponding to the target waveform can be obtained. Therefore, by controlling the automatic transmission using the input command value obtained in this way, even if there is a change in the input command value or the characteristic value of the component of the device that controls the shift operation, the output is performed. A speed change operation that suppresses the influence on is achieved.
なお、ファジィ集合や規則群の数を増やせば、問題は複雑になるが、その一方で、よりきめ細やかな目標波形の調整が可能となる。 If the number of fuzzy sets and rule groups is increased, the problem becomes complicated, but on the other hand, the target waveform can be adjusted more finely.
上述した方法により、入力指令値もしくは変速動作を制御する装置の構成要素の特性値の変動に対して頑強な目標波形、およびそれを実現するためのフィードフォワード入力指令値を、実際にシミュレーションにて導き出した結果を以下に示す。 By the above-described method, the target waveform that is robust against fluctuations in the input command value or the characteristic value of the component of the device that controls the shift operation, and the feedforward input command value for realizing the target waveform are actually simulated. The derived results are shown below.
シミュレーションに使用したCAEモデルを図8に示したものであり、この変速機の2速-3速アップシフトを対象とした。 The CAE model used for the simulation is shown in FIG. 8, and the 2nd-3rd speed upshift of this transmission was targeted.
図9に、初期入力指令値に変動を加えた状態で、変速シミュレーションを行った結果を示す(図1におけるステップS2の結果)。 FIG. 9 shows the result of the shift simulation performed with the initial input command value varied (result of step S2 in FIG. 1).
図9(c)は、係合側・解放側の摩擦係合要素へ供給される油圧波形の初期状態を表す。このように、両ブレーキ油圧とも、比較的簡単なランプ波形としている。図9(d)は、エンジンの点火時期の遅角量を表す。このように、初期状態では遅角量は常に零とした。 FIG. 9C shows the initial state of the hydraulic waveform supplied to the friction engagement elements on the engagement side / release side. Thus, both brake hydraulic pressures have a relatively simple ramp waveform. FIG. 9D shows the retard amount of the ignition timing of the engine. Thus, in the initial state, the retardation amount is always zero.
図9(a),(b)は、変速シミュレーションの結果として出力されたエンジン回転数、変速機の出力軸トルクをそれぞれ表す。図9(a),(b)の実線は、図9(c),(d)に表されている各入力指令値波形を、シミュレーション実行時にそのまま車両への入力とした場合(これをノミナル状態と呼ぶ)のシミュレーション結果を表す。 FIGS. 9A and 9B respectively show the engine speed and the output shaft torque of the transmission that are output as a result of the shift simulation. The solid lines in FIGS. 9 (a) and 9 (b) show the cases where the input command value waveforms shown in FIGS. 9 (c) and 9 (d) are directly input to the vehicle when the simulation is executed (this is the nominal state). The simulation results).
一方、係合側油圧波形に意図的に±10[kPa]のオフセット(バラツキ)を加えて、シミュレーションを実行した結果を図9(a),(b)に破線で示す。この実線と破線との違いから、入力指令値の変動により、変速ショックが大きくなったり、変速時間が間延びしてしまうことが分かる。図9(a),(b)には、許容範囲についての記載は省略してあるが、この場合図1のステップS3の判定はNoになることにする。 On the other hand, the result of executing simulation by intentionally adding an offset (variation) of ± 10 [kPa] to the engagement-side hydraulic waveform is shown by broken lines in FIGS. 9 (a) and 9 (b). From the difference between the solid line and the broken line, it can be seen that the shift shock increases or the shift time is prolonged due to the fluctuation of the input command value. In FIG. 9A and FIG. 9B, description of the allowable range is omitted, but in this case, the determination in step S3 of FIG. 1 is No.
そこで、図9で示したシミュレーション結果をファジィ推論器に入力し、目標出力軸トルク波形を修正する(図1のステップS4)。その結果を図10に示す。なお、この修正は上述のように、ファジィ推論の後件部の処理による。 Therefore, the simulation result shown in FIG. 9 is input to the fuzzy inference device to correct the target output shaft torque waveform (step S4 in FIG. 1). The result is shown in FIG. Note that this correction is based on the processing of the consequent part of fuzzy inference as described above.
入力指令値の変動の影響を抑えるために、各時間での出力軸トルクの勾配や、高さが変更されていることが見て取れる。なお、図10には、目標出力軸トルク波形しか示していないが、出力軸トルクの波形が変更になった場合には、物理現象として矛盾が生じないように、同時に目標タービン回転数の波形の修正も行い、変速動作において、ギヤトレーン内部の各回転体の回転数が同期した瞬間に、出力軸トルクも変速が終了した状態になるように調整を行う。 It can be seen that the gradient and height of the output shaft torque are changed at each time in order to suppress the influence of fluctuations in the input command value. FIG. 10 shows only the target output shaft torque waveform. However, when the output shaft torque waveform is changed, the waveform of the target turbine speed is simultaneously changed so that no contradiction occurs as a physical phenomenon. Correction is also performed, and in the speed change operation, the output shaft torque is also adjusted so that the speed change is completed at the moment when the rotational speeds of the rotating bodies in the gear train are synchronized.
このようにして、目標とする出力軸トルク波形が設定された場合には、図1におけるS5に進み、この目標波形に近い出力結果を得るための入力指令値を、学習計算に基づく繰り返し計算手法を利用して算出する。この学習計算手法は、システムの運動方程式を用いて、適切な指令値の修正量を逐次求めることで、目標値との偏差を徐々に減少させることを特徴としている。図11に、この学習計算の流れを示す。 When the target output shaft torque waveform is set in this way, the process proceeds to S5 in FIG. 1, and an input command value for obtaining an output result close to the target waveform is calculated repeatedly based on learning calculation. Calculate using. This learning calculation method is characterized in that the deviation from the target value is gradually reduced by sequentially obtaining an appropriate correction value of the command value using the equation of motion of the system. FIG. 11 shows the flow of this learning calculation.
ここでは、図8におけるブレーキB1,B2の掴み換え動作を例に取り、学習計算手法の具体的な説明を行う。自動変速機の変速時に関する運動方程式は、下式に示すように記述される。 Here, a specific description of the learning calculation method will be given, taking as an example the gripping and changing operation of the brakes B1 and B2 in FIG. The equation of motion relating to the shift of the automatic transmission is described as shown in the following equation.
ここで、TOは出力軸トルク、ωt・(=dωt/dt)はタービン角加速度、TB1は解放側ブレーキB1のトルク容量、TB2は係合側ブレーキB2のトルク容量、Ttはタービントルク、Twは走行負荷トルクを表す。そして、上式内の各係数 A11,A12,・・・,A24は,遊星歯車によって構成されるギヤトレーン系のギヤ比や各部の慣性モーメントによって決定される定数を表す。 Here, T O is the output shaft torque, ω t · (= dω t / dt) is the turbine angular acceleration, T B1 is the torque capacity of the disengagement side brake B1, T B2 is the torque capacity of the engagement side brake B2, and T t Is the turbine torque, and Tw is the running load torque. Each coefficient A 11 , A 12 ,..., A 24 in the above expression represents a constant determined by the gear ratio of the gear train system constituted by the planetary gears and the moment of inertia of each part.
行列表現を用いて、前述の式を下式に示すように記述し直す。 Using the matrix representation, the above equation is rewritten as shown in the following equation.
そして、図1のS1と同様に、図11のS11において油圧波形および遅角量の初期値を決定して、ブレーキトルク容量・タービントルクを入力として、変速シミュレーションを実行する(図11ステップS12)。 Then, similarly to S1 in FIG. 1, the initial values of the hydraulic pressure waveform and the retard amount are determined in S11 in FIG. 11, and the shift simulation is executed with the brake torque capacity and the turbine torque as inputs (step S12 in FIG. 11). .
図11のステップS13において、出力値波形が目標波形に十分近いなどの学習の終了条件を満たしているか、すなわち出力軸トルク・タービン角加速度の目標値との間の誤差ΔTO,Δωtが所定値以下かを判定し、Noの場合には図11のステップS14に進み、Aの擬似逆行列A+及び目標値への収束性を調整するゲイン行列Kを用いてブレーキトルク容量・タービントルクの修正量を次式を用いて求める。 In step S13 in FIG. 11, whether the output value waveform is sufficiently close to the target waveform or the like, ie, the errors ΔT O and Δω t between the output shaft torque and the target value of the turbine angular acceleration are predetermined. If it is No, the process proceeds to step S14 in FIG. 11 and the pseudo inverse matrix A + of A and the gain matrix K for adjusting the convergence to the target value are used to determine the brake torque capacity / turbine torque. The correction amount is obtained using the following equation.
ブレーキトルク容量に関しては、この修正量をこれまでのブレーキトルク容量に加算することで、新たなブレーキトルク容量とする。 With respect to the brake torque capacity, a new brake torque capacity is obtained by adding this correction amount to the previous brake torque capacity.
図11ステップS15においては,ステップS14の結果を受けてフィードフォワード指令値の修正を行う。係合側ブレーキに関しては、トルク相・イナーシャ相を通じて、以下に示す関係式より、トルク容量を発生させるブレーキ油圧へ変換される。 In step S15 of FIG. 11, the feedforward command value is corrected based on the result of step S14. The engaging brake is converted to a brake hydraulic pressure that generates torque capacity through the torque phase and inertia phase according to the following relational expression.
ここで、PB2はブレーキ油圧、μB2は摩擦係数、NB2は摩擦材枚数、RB2は摩擦材有効半径、SB2はピストン受圧面積、FB2はリターンスプリング荷重を表す。 Here, P B2 is the brake hydraulic pressure, μ B2 is the friction coefficient, N B2 is the number of friction materials, R B2 is the friction material effective radius, S B2 is the piston pressure receiving area, and F B2 is the return spring load.
また、トルク相中においては、係合側ブレーキのトルク容量TB2の修正に対して,この係合側ブレーキのトルク容量とタービントルクとの関係から、タイアップ及びエンジン回転数の吹き現象が発生しない様に、解放側ブレーキの油圧波形が修正される。 In the torque phase, to the correction of the torque capacity T B2 of the engaging brake, the relationship between the torque capacity and the turbine torque of the engagement side brake, is blown phenomenon of tie-up and the engine speed occurs The hydraulic waveform of the release side brake is corrected so that it does not.
さらに、イナーシャ相中においては、出力軸トルクのピーク値を抑えるために、エンジンの遅角制御が行われる。タービントルクの修正量から、次式に従いエンジン遅角量の修正量を求める。 Further, during the inertia phase, the retard control of the engine is performed in order to suppress the peak value of the output shaft torque. From the turbine torque correction amount, the engine retardation amount correction amount is obtained according to the following equation.
ここで、αは遅角量、Δαは遅角量修正量、tはエンジントルクとタービントルクのトルク比、kは定数を表す。 Here, α represents a retard amount, Δα represents a retard amount correction amount, t represents a torque ratio between the engine torque and the turbine torque, and k represents a constant.
各フィードフォワード指令値の修正後は、図11ステップS12に戻り,修正された指令値による変速シミュレーションを再度行う。 After the correction of each feedforward command value, the process returns to step S12 in FIG. 11, and the shift simulation using the corrected command value is performed again.
この指令値の修正作業は、図11ステップS13において、学習計算の終了条件が満たされるまで繰り返し実行される。学習の終了条件の例には、学習計算の繰り返し回数が所定の回数を超えた場合や、目標値との誤差が十分に小さい値になった場合などが挙げられる。 This command value correction operation is repeatedly executed until the learning calculation end condition is satisfied in step S13 of FIG. Examples of the learning end condition include a case where the number of repetitions of the learning calculation exceeds a predetermined number or a case where an error from the target value becomes a sufficiently small value.
このような学習計算を用いて、目標波形を実現するための入力指令値を求めた結果を、図12に示す。図12の(a)〜(e)において、破線は初期状態での波形を、実線は学習後の波形を示す。 FIG. 12 shows the result of obtaining the input command value for realizing the target waveform using such learning calculation. In (a) to (e) of FIG. 12, a broken line indicates a waveform in an initial state, and a solid line indicates a waveform after learning.
図12(f)には、出力軸トルクの誤差の履歴を示す。学習計算を繰り返し行い、入力指令値を繰り返し修正していくことで、適切なフィードフォワード入力指令値が求められている。 FIG. 12F shows a history of output shaft torque errors. By repeating the learning calculation and repeatedly correcting the input command value, an appropriate feedforward input command value is obtained.
上述した手法に従い、目標波形の修正とそれに対応する入力指令値の導出を交互に繰り返した結果、最終的に得られた入力指令値での変速シミュレーション結果を図13に示す。 FIG. 13 shows a shift simulation result with the input command value finally obtained as a result of alternately repeating the correction of the target waveform and the derivation of the input command value corresponding to the target waveform according to the above-described method.
図13(c)は係合側・解放側の摩擦係合要素へ供給される油圧波形、図13(d)はエンジンの点火時期の遅角量を表す。また、図13(a),(b)は、変速シミュレーションの結果として出力されたエンジン回転数、変速機の出力軸トルクをそれぞれ表す。図13(a),(b)の実線は、図13(c),(d)に表されている各入力波形を、シミュレーション実行時にそのまま車両への入力とした場合のシミュレーション結果を表す。 FIG. 13C shows a hydraulic pressure waveform supplied to the engagement side / release side friction engagement element, and FIG. 13D shows a retard amount of the ignition timing of the engine. FIGS. 13A and 13B respectively show the engine speed and the output shaft torque of the transmission that are output as a result of the shift simulation. The solid lines in FIGS. 13A and 13B represent the simulation results when the input waveforms shown in FIGS. 13C and 13D are directly input to the vehicle when the simulation is executed.
一方、係合側油圧波形に意図的に±10[kPa]のオフセットを加えて、シミュレーションを実行した結果を図13(a),(b)の破線にて示す。図9に示した初期の入力指令値に対して同様の変動を加えたときの結果と比較して、全体的に変速中の出力軸トルク波形のピークが抑えられており、変速ショックが小さく抑えられていることが分かる。 On the other hand, the result of executing simulation by intentionally adding an offset of ± 10 [kPa] to the engagement-side hydraulic waveform is shown by broken lines in FIGS. Compared with the result when the same fluctuation is applied to the initial input command value shown in FIG. 9, the peak of the output shaft torque waveform during the shift is suppressed as a whole, and the shift shock is reduced. You can see that
ここで、上述の説明では、出力に影響する変動として入力指令値の変動を取り上げた。しかし変動要因は、入力指令値だけではなく、係合摩擦材の摩擦係数の変動や油圧回路内の経路・絞り径(オリフィス)の寸法誤差、エンジン遅角制御によるエンジントルクの落ち込みタイミングのずれなど、変速動作を制御する装置内のあらゆる特性変動が考えられ、これら変動を対象とすることが可能である。 Here, in the above description, the variation of the input command value is taken up as the variation that affects the output. However, the fluctuation factors are not only the input command value, but also the fluctuation of the friction coefficient of the engagement friction material, the dimensional error of the path / throttle diameter (orifice) in the hydraulic circuit, the deviation of the engine torque drop timing due to the engine retarding control, etc. Any characteristic variation in the device for controlling the shifting operation can be considered, and these variations can be targeted.
このように、本実施形態の手法を用いれば、上記のように入力指令値もしくは変速動作を制御する装置の構成要素の特性値の変動に対して頑強な目標波形、およびそれに対応する各入力指令値を、ベンチでの実機試験に必要な工数と比較して、ごく短時間で計算機シミュレーションにより求めることができる。 As described above, by using the method of the present embodiment, the target waveform that is robust against fluctuations in the input command value or the characteristic value of the component of the device that controls the shift operation as described above, and each input command corresponding thereto. The value can be obtained by computer simulation in a very short time compared with the man-hours required for the actual machine test on the bench.
「感度解析を利用する実施形態」
上記実施形態においては、ファジィ推論を用いて、入力指令値もしくは変速動作を制御する装置の構成要素の特性値の変動に強い目標波形を求め、その目標波形を用いて入力指令値を決定した。
"Embodiment using sensitivity analysis"
In the above embodiment, a fuzzy inference is used to obtain a target waveform that is resistant to fluctuations in the input command value or the characteristic value of the component of the device that controls the shift operation, and the input command value is determined using the target waveform.
本実施形態では、時系列のフィードフォワード指令値を設計するために品質工学の手法である感度解析を応用する。 In the present embodiment, sensitivity analysis, which is a quality engineering technique, is applied to design time-series feedforward command values.
本実施形態での適合作業の流れを図14に示す。なお、以下の説明では、図8に示した自動変速機を搭載した自動車が低速段で走行中に、一度スロットル開度を零まで戻して中速段を飛び越え高速段へアップシフトした後、再度スロットル開度を上げてそのまま高速段で走行するパワオフ・パワオン多重変速を対象とする。 FIG. 14 shows the flow of adaptation work in this embodiment. In the following description, while a vehicle equipped with the automatic transmission shown in FIG. 8 is traveling at a low speed, the throttle opening is once returned to zero, jumped over the medium speed and upshifted to a high speed, and then again. Power-off / power-on multiple shifts that run at high speeds with the throttle opening increased are the targets.
ステップS21では、まず自動変速機に入力する指令値の初期値を設定する。この初期指令値は、単に所定の変速段の切り換え作業を行うことができればよいため、図15の左側に示すような簡単な形状の入力指令値を与えればよい。 In step S21, first, an initial value of a command value input to the automatic transmission is set. Since the initial command value only needs to be able to perform a predetermined gear shift operation, an input command value having a simple shape as shown on the left side of FIG. 15 may be given.
次に,図14のステップS22では、図15の右側に示すように、ステップS21での出力結果を参考にして、出力軸トルク波形の特徴量に関して3つの水準を設定する。ここでは、特徴量として、スロットル開度が立ち上がった瞬間の出力軸トルクの値To1、変速終了時に出力軸トルクが立ち下がる直前の出力軸トルクTo2、変速時間tscの3つを採用し、これを目標値における制御因子とする。 Next, in step S22 of FIG. 14, as shown on the right side of FIG. 15, with reference to the output result in step S21, three levels are set for the feature quantity of the output shaft torque waveform. Here, three characteristic values are adopted: the output shaft torque value To1 at the moment when the throttle opening rises, the output shaft torque To2 immediately before the output shaft torque falls at the end of the shift, and the shift time tsc. The control factor at the target value.
これ以後では、この制御因子の組み合わせに関して検討するが,図16に示すように制御因子の組み合わせは無数に存在する。このため、制御因子の数が増えると設計のために実行しなければならないシミュレーション回数が膨大となる。そこで、本実施形態においては、これらの制御因子を図17に示す直交表に割り付けることにより、水準の組み合わせに関する検討を効率よく実施する。以後、図17の直交表の制御因子AにTo1、BにTo2、Cにtscが割り付けられているとする。 In the following, this combination of control factors will be discussed, but there are innumerable combinations of control factors as shown in FIG. For this reason, when the number of control factors increases, the number of simulations that must be executed for design becomes enormous. Therefore, in this embodiment, these control factors are allocated to the orthogonal table shown in FIG. Hereinafter, it is assumed that To1 is assigned to control factor A in the orthogonal table of FIG. 17, To2 is assigned to B, and tsc is assigned to C.
図14のステップS23においては、図17の直交表に従って設定された目標出力軸トルク波形を実現するフィードフォワード入力指令値を学習計算にて求める。この学習計算による入力指令値の算出は、上述の実施形態の場合と同様にして行われる。 In step S23 of FIG. 14, a feedforward input command value for realizing the target output shaft torque waveform set according to the orthogonal table of FIG. 17 is obtained by learning calculation. The calculation of the input command value by this learning calculation is performed in the same manner as in the above-described embodiment.
図14ステップS24では、学習計算手法により導き出されたブレーキ油圧のフィードフォワード入力指令値に±βのオフセットを与えて、これを誤差(変動)因子として図17の直交表の外側に割り付ける。なお,誤差(変動)因子βには、ブレーキ油圧のオフセット(バラツキ)だけではなく、係合摩擦材の摩擦係数の変動や油圧回路内の経路・絞り径(オリフィス)の寸法誤差、エンジン遅角制御によるエンジントルクの落ち込みタイミングのずれなど、変速動作を制御する装置内のあらゆる特性変動を与えることが可能である。変動として、入力指令値だけでなく変速動作を制御する装置の構成要素の特性値の変動を与えることができる。 In step S24 of FIG. 14, an offset of ± β is given to the feed pressure input command value of the brake hydraulic pressure derived by the learning calculation method, and this is assigned outside the orthogonal table of FIG. 17 as an error (variation) factor. The error (variation) factor β includes not only the brake hydraulic pressure offset (variation) but also the variation of the friction coefficient of the engagement friction material, the path / throttle diameter (orifice) dimensional error in the hydraulic circuit, and the engine delay angle. It is possible to give any characteristic variation in the device that controls the shift operation, such as a shift in engine torque drop timing due to control. As the variation, not only the input command value but also the variation of the characteristic value of the component of the device that controls the shift operation can be given.
図14ステップS25では、この外乱が加わった入力指令値にて変速シミュレーションを実行する。そして、得られた結果から変速動作が円滑に終了したかを判断するために定められた各評価指標のSN比・平均値を求める。 In step S25 of FIG. 14, a shift simulation is executed with the input command value to which this disturbance has been added. Then, the S / N ratio / average value of each evaluation index determined to determine whether the speed change operation has been completed smoothly is obtained from the obtained result.
図18に評価指標の例を示す。この例では、評価指標として図中に表されているように出力軸トルクの増減から計算する変速ショック、変速時間、変速時のエンジンの吹き上がり回転数などを採用する。また、変速時には変速機内部の歯車の回転同期をとるために発生する熱の処理が問題となるため、摩擦係合要素が吸収しなければならない熱吸収量、熱吸収率も評価指標とする。なお、この他にも、変速中の車両の加速度やジャーク、動力をタイヤに伝えるドライブシャフトの捻り振動なども評価指標として採用することができる。 FIG. 18 shows an example of the evaluation index. In this example, as shown in the figure, a shift shock, shift time, engine speed at the time of shift, and the like calculated from increase / decrease in output shaft torque are employed as evaluation indexes. In addition, since the processing of the heat generated to synchronize the rotation of the gears inside the transmission becomes a problem at the time of shifting, the heat absorption amount and the heat absorption rate that the friction engagement element must absorb are also used as evaluation indexes. In addition, the acceleration, jerk, and torsional vibration of the drive shaft that transmits the power to the tire can be employed as an evaluation index.
また、図17の直交表のn番目の制御因子の組み合わせに関するSN比、平均値は、以下に示す式により求められる。 Moreover, the SN ratio and the average value regarding the combination of the nth control factor in the orthogonal table of FIG. 17 are obtained by the following equations.
ここで、yn,1,yn,2,yn,3は,各誤差(変動)因子を加えた場合での変速シミュレーション結果から求められた評価指標の値を表し、ηはそのSN比、μは平均値を表す。 Here, y n, 1 , y n, 2 , y n, 3 represent the values of evaluation indices obtained from the shift simulation results when each error (variation) factor is added, and η is the SN ratio , Μ represents an average value.
上記の図14ステップS23からステップS25までの作業を、図17の直交表で定められた全ての因子の組み合わせについて実行する。すなわち、ステップS26において、全ての組み合わせが終了が終了したかを判定し、終了していなかった場合にはステップS23に戻り処理を繰り返す。 The above-described operations from step S23 to step S25 in FIG. 14 are executed for all factor combinations defined in the orthogonal table in FIG. That is, in step S26, it is determined whether or not the end of all combinations has been completed. If not completed, the process returns to step S23 to repeat the process.
図14ステップS27では、全ての変速シミュレーション結果から得られたデータを基に、要因効果図を作成し、各評価指標のバラツキを抑えつつ、平均値を設計要求を満たす水準まで移動させることができる制御因子の組み合わせについて検討を行う。 In step S27 of FIG. 14, a factorial effect diagram is created based on data obtained from all shift simulation results, and the average value can be moved to a level that satisfies the design requirements while suppressing variations in each evaluation index. Consider combinations of control factors.
図19から図23に、変速ショック、変速時間、エンジンの吹き上がり回転数、摩擦材の熱吸収量、熱吸収率に関する要因効果図をそれぞれ示す。これらの効果図において、SN比は値が大きいほうが、平均値は値が小さいほうが望ましい傾向を示している。 FIG. 19 to FIG. 23 show factor effect diagrams relating to shift shock, shift time, engine speed of rotation, friction material heat absorption, and heat absorption rate, respectively. In these effect diagrams, it is desirable that the SN ratio has a larger value and the average value has a smaller value.
これらの要因効果図から各制御因子の値をどのように設定すれば、誤差(変動)因子の影響を減衰させるのかについて検討を行うが、図19から図23を見比べると分かるように、図19のSN比においてはC3の水準が望ましいが、それ以外ではC1の水準が望ましいといった具合に相反する傾向を持つ場合がある。 It is examined how to set the value of each control factor from these factor effect diagrams to attenuate the influence of the error (fluctuation) factor. As can be seen by comparing FIG. 19 to FIG. 23, FIG. In the S / N ratio, the C3 level is desirable, but the C1 level is desirable in other cases.
そこで、このような矛盾点を回避しつつ、最適な解を求めるため非線形最適化の手法を用いる。まず、要因効果図内の特性を2次関数で近似し、以下に示すような評価関数を作る。 Therefore, a nonlinear optimization technique is used to obtain an optimal solution while avoiding such contradictions. First, the characteristics in the factor-effect diagram are approximated by a quadratic function to create an evaluation function as shown below.
ここで、変数xは制御因子が水準1から水準3までの間のどの値を持つのかを表し、ai,bi,ci(i=1,2,・・・,10)は2次関数の係数を表す。なお、上記の式において、本来SN比は大きければ大きいほど良い特性を表すが、符号を反対にすることにより、関数ベクトルF(x)の最適化をする際に、最小値を求めるように変更する。
Here, the variable x represents what value the control factor has between
上記のF(x)のような非線形の多目的関数ベクトルを最小化する際には、ファジィ推論、遺伝的アルゴリズム、ニューロコンピューティングなどの手法を利用することも考えられるが、ここでは数理計画法の一種である多目的ゴール到達法と呼ばれる解法を利用する。多目的ゴール到達法というのは、以下の式で示すように、重み関数wiと設計ゴール f* iとで表される制約条件のもとで、変数γを最小化する手法である。 When minimizing a nonlinear multi-objective function vector such as F (x) above, it is possible to use techniques such as fuzzy inference, genetic algorithm, and neurocomputing. A so-called multi-objective goal attainment method is used. The multi-objective goal attainment method is a method of minimizing the variable γ under the constraint condition represented by the weight function w i and the design goal f * i as shown by the following equation.
この手法を適用して最適化を行った結果、表2に示す各水準が最適値として得られた。なお、実際の計算には The MathWorks,Inc.製のソフトウェアであるMATLAB(登録商標)を使用した。なお、各水準についての最適値を得る手法としては他の最適化手法を採用してもよい。 As a result of optimization by applying this method, the levels shown in Table 2 were obtained as optimum values. For the actual calculation, The MathWorks, Inc. The software manufactured by MATLAB (registered trademark) was used. It should be noted that other optimization methods may be employed as a method for obtaining the optimum value for each level.
図14ステップS28では、上記の最適化された制御因子から設計された目標出力軸トルク波形を実現する各フィードフォワード入力指令値を、前述の学習計算手法により導出する。 In FIG. 14 step S28, each feedforward input command value which implement | achieves the target output shaft torque waveform designed from the said optimized control factor is derived | led-out by the above-mentioned learning calculation method.
図24には初期入力指令値による変速シミュレーション結果、図25には最終的に得られた入力指令値によるシミュレーション結果を示す。各図の(a)は出力軸トルク、(b)はエンジン回転数、(c)はエンジントルク、(d)は係合側・解放側の摩擦係合要素へ供給されるブレーキ油圧、(e)はスロットル開度、(f)はエンジンの点火時期の遅角量をそれぞれ表す。 FIG. 24 shows a shift simulation result based on the initial input command value, and FIG. 25 shows a simulation result based on the finally obtained input command value. In each figure, (a) is the output shaft torque, (b) is the engine speed, (c) is the engine torque, (d) is the brake hydraulic pressure supplied to the engagement / release side frictional engagement elements, (e ) Represents the throttle opening, and (f) represents the retard amount of the ignition timing of the engine.
図24では、ブレーキ油圧を急激に変化させているため、図24(a)の矢印で指摘している箇所で大きな変速ショックが発生しているが、最適化された指令値を用いた図25では、ブレーキ油圧を緩やかに変化させると同時に、エンジンの点火時期を制御することでエンジンから出力されるトルクを減少させたことにより、出力軸トルク波形に現れる変速ショックを軽減させている。 In FIG. 24, since the brake hydraulic pressure is changed abruptly, a large shift shock is generated at the position indicated by the arrow in FIG. 24A, but FIG. 25 using an optimized command value is used. In this case, the brake hydraulic pressure is gradually changed and at the same time the torque output from the engine is reduced by controlling the ignition timing of the engine, thereby reducing the shift shock appearing in the output shaft torque waveform.
また、図26には、最適化された係合側ブレーキへの入力指令値に、外乱として−10[kPa]のオフセットを加えた場合のシミュレーション結果を示す。予めこのような外乱の影響を軽減することを念頭におき、目標出力軸トルク波形、およびフィードフォワード入力指令値を設計したため、図25と図26の出力結果の差異は、最小限に抑えられている。 FIG. 26 shows a simulation result when an offset of −10 [kPa] is added as a disturbance to the optimized input command value to the engagement-side brake. Since the target output shaft torque waveform and the feedforward input command value are designed in advance in consideration of reducing the influence of such disturbances, the difference between the output results in FIG. 25 and FIG. 26 is minimized. Yes.
このように、本実施形態に係る手法を用いれば、上記のようにシステム内部の特性変動(入力指令値もしくは変速動作を制御する装置の構成要素の特性値)に対して頑強な目標出力波形、およびそれに対応する各フィードフォワード入力指令値を、ベンチでの実機試験に必要な工数と比較して、ごく短時間で計算機シミュレーションにより求めることができる。また、設計の際に、外乱・ノイズの影響を定量的に求めることにより、そのシステムの特性を熟知した熟練者でなくとも、容易に最適設計を行うことが可能となる。 As described above, when the method according to the present embodiment is used, the target output waveform that is robust against the characteristic fluctuation (input command value or characteristic value of the component of the device that controls the shift operation) inside the system as described above, And each feedforward input command value corresponding to it can be calculated | required by computer simulation in a very short time compared with the man-hour required for the real machine test on a bench. In addition, by quantitatively determining the influence of disturbance and noise at the time of designing, it is possible to easily carry out an optimum design even if it is not an expert who is familiar with the characteristics of the system.
「ファジィ推論を用いる場合と感度解析を用いる場合の比較」
ファジィ推論を用いた実施形態と、感度解析を用いた実施形態との共通点・相違点について、下に示す。
"Comparison of using fuzzy reasoning and sensitivity analysis"
Common points and differences between the embodiment using fuzzy inference and the embodiment using sensitivity analysis are shown below.
[共通点]
・自動車用自動変速機からの出力値が,システム内部に特性変動が発生した場合においても、常にある一定の要求仕様を満たすことを目的とする。
・上記の目的を実現するための制御指令値を自動的に求める手法である。
・上記の要求仕様は、複数の評価指標(変速ショックを表す官能指数・摩擦材の熱吸収量など)の値によって定められる。
・上記で述べた要求仕様を満たすために、特性変動によって出力値に悪影響が出た場合の評価指標の値を参考にして、特性変動に対して頑強な目標出力値を導き出す。
・続いて、実際に目標出力値に一致する出力結果を得るための制御入力を学習計算手法にて求める。
[Common point]
・ The output value from the automatic transmission for automobiles is intended to always satisfy a certain required specification even when characteristic fluctuations occur in the system.
A method for automatically obtaining a control command value for realizing the above object.
The above-mentioned required specifications are determined by the values of a plurality of evaluation indexes (sensory index representing shift shock, heat absorption amount of friction material, etc.).
In order to satisfy the required specifications described above, a target output value that is robust against characteristic fluctuation is derived with reference to the value of the evaluation index when the output value is adversely affected by characteristic fluctuation.
Subsequently, a control input for obtaining an output result that actually matches the target output value is obtained by a learning calculation method.
[相違点]
「ファジィ推論を利用する場合」
・メンバシップ関数・修正ルール群の作成に、熟練適合者の経験と、試行錯誤のために工数が必要である。
・上記の関数・ルール群が適切に設定されていれば、正しい修正量を推論結果として得ることができる。
[Difference]
"When fuzzy reasoning is used"
・ Establishing membership functions / correction rules requires man-hours for the experience of skilled adaptors and trial and error.
If the above function / rule group is set appropriately, the correct correction amount can be obtained as an inference result.
「感度解析を利用する場合」
・入力の変動値による出力値の変動幅を定量的に評価する。
・混合系直交表を利用することで、最小限の試行回数で,効率よく感度解析を実行することが可能となる。
・定量的な解析結果を基に、そこから最適な解を得る手法が別途必要になる。
"When using sensitivity analysis"
・ Quantitatively evaluate the fluctuation range of the output value due to the fluctuation value of the input.
-By using a mixed system orthogonal table, sensitivity analysis can be performed efficiently with a minimum number of trials.
・ A separate method for obtaining the optimal solution based on the quantitative analysis results is required.
10 トルクコンバータモデル、12 ギヤトレーンモデル、14 解放側ブレーキモデル、16 係合側ブレーキモデル、18 修正量計算部、20,22,24,26,28 加算器。 10 Torque converter model, 12 Gear train model, 14 Release side brake model, 16 Engagement side brake model, 18 Correction amount calculation unit, 20, 22, 24, 26, 28 Adder.
Claims (8)
前記システムにおいて変速動作を行う際の入力指令値もしくは変速動作を制御する装置の構成要素の特性値について、変動を付与してシミュレーションし、
このシミュレーションの結果であるシステムの出力値を、前記変動に対する影響の大きさに基づいて評価し、
この評価結果から、前記変動の出力値への影響が要求仕様を満たす目標出力値を決定し、
決定された目標出力値に一致する出力を得るための入力指令値を学習計算によって求めることで、要求仕様を満たす入力指令値を求めることを特徴とする変速シミュレーション装置。 A shift simulation device that simulates a shift operation of an automatic transmission system for an automobile using a model,
For the input command value when performing a shift operation in the system or the characteristic value of the component of the device that controls the shift operation, a simulation is performed by giving a variation,
The output value of the system as a result of this simulation is evaluated based on the magnitude of the influence on the fluctuation,
From this evaluation result, a target output value that satisfies the required specifications is determined for the influence of the variation on the output value,
A shift simulation device characterized by obtaining an input command value that satisfies a required specification by obtaining an input command value for obtaining an output that matches the determined target output value by learning calculation.
(a)前記システムにおいて変速動作を行う際の入力指令値もしくは変速動作を制御する装置の構成要素の特性値について、変動を付与してシミュレーションし、
(b)このシミュレーションの結果であるシステムの出力値を、前記変動に対する影響の大きさに基づいて評価し、
(c)評価結果に基づき、ファジィ推論によって、目標出力値を修正し、
(d)修正した目標出力値を用いて、シミュレーションを行い、その学習計算により新たな入力指令値を算出し、
(e)算出した入力指令値により前記入力指令値を更新し、前記(a)に戻り、前記(a)〜(d)の処理を(b)の評価が要求仕様を満たすまで繰り返し、
要求仕様を満たす入力指令値を求めることを特徴とする変速シミュレーション装置。 A shift simulation device that simulates a shift operation of an automatic transmission system for an automobile using a model,
(A) A simulation is performed by giving a variation to an input command value or a characteristic value of a component of a device that controls the shift operation when performing a shift operation in the system;
(B) evaluating the output value of the system, which is the result of this simulation, based on the magnitude of the effect on the variation;
(C) Based on the evaluation result, the target output value is corrected by fuzzy inference,
(D) A simulation is performed using the corrected target output value, a new input command value is calculated by the learning calculation,
(E) Update the input command value with the calculated input command value, return to (a), repeat the processes (a) to (d) until the evaluation in (b) satisfies the required specifications,
A speed change simulation apparatus characterized by obtaining an input command value satisfying a required specification.
(a)目標出力値における複数の制御因子についてそれぞれ複数の水準を設定し、
(b)設定された複数の制御因子についての複数の水準の組み合わせで決定される複数の目標出力値を設定し、
(c)設定された複数の目標出力値を用い、前記システムにおいて変速動作を行う際の入力指令値もしくは変速動作を制御する装置の構成要素の特性値について変動を付与してそれぞれシミュレーションし、
(d)このシミュレーションの結果であるシステムの複数の出力値を、前記変動に対する影響の大きさに基づいて評価し、得られた評価結果が要求仕様を満たす目標出力値を決定し、
(e)決定された目標出力値を用いて、学習計算によって要求仕様を満たす入力指令値を求める、
ことを特徴とする変速シミュレーション装置。 A shift simulation device that simulates a shift operation of an automatic transmission system for an automobile using a model,
(A) A plurality of levels are set for a plurality of control factors at the target output value,
(B) setting a plurality of target output values determined by a combination of a plurality of levels for a plurality of set control factors;
(C) Using a plurality of set target output values, simulating each of the input command value or the characteristic value of the component of the device that controls the shift operation by performing a shift operation in the system, respectively,
(D) evaluating a plurality of output values of the system as a result of the simulation based on the magnitude of the influence on the fluctuation, and determining a target output value that the obtained evaluation result satisfies the required specifications;
(E) Using the determined target output value, find an input command value that satisfies the required specifications by learning calculation;
A shift simulation apparatus characterized by the above.
前記(b)において、複数の制御因子の複数の水準の組み合わせは、直交表を用いて、決定し、前記(c)は決定された全ての組み合わせについて行われることを特徴とする変速シミュレーション装置。 In the shift simulation device according to claim 3,
In (b), a combination of a plurality of levels of a plurality of control factors is determined using an orthogonal table, and (c) is performed for all the determined combinations.
コンピュータに、
前記システムにおいて変速動作を行う際の入力指令値もしくは変速動作を制御する装置の構成要素の特性値について、変動を付与してシミュレーションさせ、
このシミュレーションの結果であるシステムの出力値を、前記変動に対する影響の大きさに基づいて評価させ、
この評価結果から、前記変動の出力値への影響が要求仕様を満たす目標出力値を決定させ、
決定された目標出力値に一致する出力を得るための入力指令値を学習計算によって求めさせることで、要求仕様を満たす入力指令値を求めさせることを特徴とする変速シミュレーションプログラム。 A shift simulation program for simulating a shift operation of an automatic transmission system for an automobile using a model,
On the computer,
For the input command value or the characteristic value of the component of the device that controls the shift operation when performing the shift operation in the system, the simulation is performed by giving a variation,
The output value of the system as a result of this simulation is evaluated based on the magnitude of the influence on the fluctuation,
From this evaluation result, the target output value that satisfies the required specifications is determined by the influence of the variation on the output value,
A shift simulation program characterized in that an input command value for obtaining an output that matches a determined target output value is obtained by learning calculation to obtain an input command value that satisfies a required specification.
コンピュータに、
(a)前記システムにおいて変速動作を行う際の入力指令値もしくは変速動作を制御する装置の構成要素の特性値について、変動を付与してシミュレーションさせ、
(b)このシミュレーションの結果であるシステムの出力値を、前記変動に対する影響の大きさに基づいて評価させ、
(c)評価結果に基づき、ファジィ推論によって、目標出力値を修正させ、
(d)修正した目標出力値を用いて、シミュレーションを行い、その学習計算により新たな入力指令値を算出させ、
(e)算出した入力指令値により前記入力指令値を更新し、前記(a)に戻り、前記(a)〜(d)の処理を(b)の評価が要求仕様を満たすまで繰り返させ、
要求仕様を満たす入力指令値を求めることを特徴とする変速シミュレーションプログラム。 A shift simulation program for simulating a shift operation of an automatic transmission system for an automobile using a model,
On the computer,
(A) The input command value when performing a shift operation in the system or the characteristic value of the component of the device that controls the shift operation is simulated by giving a variation,
(B) let the system output value, which is the result of this simulation, be evaluated based on the magnitude of the effect on the variation;
(C) Based on the evaluation result, the target output value is corrected by fuzzy inference,
(D) A simulation is performed using the corrected target output value, and a new input command value is calculated by the learning calculation.
(E) Update the input command value with the calculated input command value, return to (a), repeat the processes (a) to (d) until the evaluation in (b) satisfies the required specifications,
A shift simulation program characterized by obtaining an input command value that satisfies a required specification.
コンピュータに、
(a)目標出力値における複数の制御因子についてそれぞれ複数の水準を設定させ、
(b)設定された複数の制御因子についての複数の水準の組み合わせで決定される複数の目標出力値を設定させ、
(c)設定された複数の目標出力値を用い、前記システムにおいて変速動作を行う際の入力指令値もしくは変速動作を制御する装置の構成要素の特性値について変動を付与してそれぞれシミュレーションさせ、
(d)このシミュレーションの結果であるシステムの複数の出力値を、前記変動に対する影響の大きさに基づいて評価し、得られた評価結果が要求仕様を満たす目標出力値を決定させ、
(e)決定された目標出力値を用いて、学習計算によって要求仕様を満たす入力指令値を求めさせる、
ことを特徴とする変速シミュレーションプログラム。 A shift simulation program for simulating a shift operation of an automatic transmission system for an automobile using a model,
On the computer,
(A) A plurality of levels are set for a plurality of control factors in the target output value,
(B) setting a plurality of target output values determined by a combination of a plurality of levels for a plurality of set control factors;
(C) Using a plurality of set target output values, the input command value when performing a shift operation in the system or a characteristic value of a component of a device that controls the shift operation is given a variation, and each is simulated.
(D) Evaluating a plurality of output values of the system as a result of the simulation based on the magnitude of the influence on the variation, and determining a target output value that the obtained evaluation result satisfies the required specifications;
(E) Using the determined target output value, an input command value that satisfies the required specifications is determined by learning calculation.
A shift simulation program characterized by the above.
A memory for storing a shift control command value based on an input command value obtained by using any one shift simulation device according to any one of claims 1 to 4 or any one shift simulation program according to any one of claims 5 to 7 is mounted. An automobile for controlling an automatic transmission mounted based on a shift control command value stored in the memory.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006221312A JP5134220B2 (en) | 2006-08-14 | 2006-08-14 | Shift simulation device, shift simulation program, and automobile |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006221312A JP5134220B2 (en) | 2006-08-14 | 2006-08-14 | Shift simulation device, shift simulation program, and automobile |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2008045651A true JP2008045651A (en) | 2008-02-28 |
JP5134220B2 JP5134220B2 (en) | 2013-01-30 |
Family
ID=39179587
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2006221312A Expired - Fee Related JP5134220B2 (en) | 2006-08-14 | 2006-08-14 | Shift simulation device, shift simulation program, and automobile |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5134220B2 (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010276120A (en) * | 2009-05-28 | 2010-12-09 | Toyota Central R&D Labs Inc | Automatic transmission parameter identification device and identification program |
WO2013168226A1 (en) * | 2012-05-08 | 2013-11-14 | トヨタ自動車株式会社 | Shift control device for vehicle |
JPWO2013168225A1 (en) * | 2012-05-08 | 2015-12-24 | トヨタ自動車株式会社 | Vehicle shift control device |
JP2018200539A (en) * | 2017-05-26 | 2018-12-20 | オムロン株式会社 | Learning device, learning control method, and its program |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0348303A (en) * | 1989-04-27 | 1991-03-01 | Omron Corp | Method of evaluating in fuzzy inference device |
JPH0979362A (en) * | 1995-09-12 | 1997-03-25 | Mitsubishi Motors Corp | Gear shift control device for automatic transmission for vehicle |
JP2001337716A (en) * | 2000-05-26 | 2001-12-07 | Isuzu Motors Ltd | Attribute identification method and its system |
JP2002122223A (en) * | 2000-08-11 | 2002-04-26 | Honda Motor Co Ltd | Development supporting device for controller of automatic transmission for vehicle |
-
2006
- 2006-08-14 JP JP2006221312A patent/JP5134220B2/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0348303A (en) * | 1989-04-27 | 1991-03-01 | Omron Corp | Method of evaluating in fuzzy inference device |
JPH0979362A (en) * | 1995-09-12 | 1997-03-25 | Mitsubishi Motors Corp | Gear shift control device for automatic transmission for vehicle |
JP2001337716A (en) * | 2000-05-26 | 2001-12-07 | Isuzu Motors Ltd | Attribute identification method and its system |
JP2002122223A (en) * | 2000-08-11 | 2002-04-26 | Honda Motor Co Ltd | Development supporting device for controller of automatic transmission for vehicle |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010276120A (en) * | 2009-05-28 | 2010-12-09 | Toyota Central R&D Labs Inc | Automatic transmission parameter identification device and identification program |
WO2013168226A1 (en) * | 2012-05-08 | 2013-11-14 | トヨタ自動車株式会社 | Shift control device for vehicle |
CN104285083A (en) * | 2012-05-08 | 2015-01-14 | 丰田自动车株式会社 | Shift control device for vehicle |
JPWO2013168225A1 (en) * | 2012-05-08 | 2015-12-24 | トヨタ自動車株式会社 | Vehicle shift control device |
JPWO2013168226A1 (en) * | 2012-05-08 | 2015-12-24 | トヨタ自動車株式会社 | Vehicle shift control device |
JP2018200539A (en) * | 2017-05-26 | 2018-12-20 | オムロン株式会社 | Learning device, learning control method, and its program |
Also Published As
Publication number | Publication date |
---|---|
JP5134220B2 (en) | 2013-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gao et al. | Design of clutch-slip controller for automatic transmission using backstepping | |
CN108454607B (en) | Model Predictive Control (MPC) based torque protection for powertrain control | |
Hu et al. | Nonlinear gearshifts control of dual-clutch transmissions during inertia phase | |
Gao et al. | Design of nonlinear shaft torque observer for trucks with Automated Manual Transmission | |
Chen et al. | Nonlinear estimation and control of automotive drivetrains | |
Lu et al. | Model predictive control of AMT clutch during start-up process | |
Čorić et al. | Optimization of shift control trajectories for step gear automatic transmissions | |
JP5134220B2 (en) | Shift simulation device, shift simulation program, and automobile | |
Lee et al. | A systematic approach for dynamic analysis of vehicles with eight or more speed automatic transmission | |
Ranogajec et al. | Multi-objective parameter optimization of automatic transmission shift control profiles | |
JP4769111B2 (en) | Automatic transmission simulation apparatus and program thereof | |
Reddy et al. | Real-time predictive clunk control using a reference governor | |
JP5237589B2 (en) | In-vehicle computer manufacturing method, in-vehicle computer data generation apparatus, and in-vehicle computer data generation program | |
Cvok et al. | Model predictive control for automatic transmission upshift inertia phase | |
Cvok et al. | An LQR approach of automatic transmission upshift control including use of off-going clutch within inertia phase | |
Jiang et al. | Drivability Improvement Control for Vehicle Start‐Up Applied to an Automated Manual Transmission | |
EP3891419B1 (en) | Method for automated calibration and adaptation of automatic transmission controllers | |
Canale et al. | A Model Predictive Control approach to micro-slip control of Automated Manual Transmission systems | |
Dong et al. | Gearshift overlap for multistep downshift of automatic transmissions based on iterative learning control | |
Hao et al. | Double‐Target Switching Control of Vehicle Longitudinal Low‐Frequency Vibration Based on Fuzzy Logic | |
Soldo et al. | Automatic transmission upshift control using a linearized reduced-order model-based LQR approach | |
Lee et al. | Model based automated calibration for shift control of automatic transmission | |
JP7487709B2 (en) | Vehicle drive unit control device and control method | |
Taratorkin et al. | Improving the quality of transmission shifting transients due to controlling torque redistribution | |
Dong et al. | Adaptive control for clutch engagement on starting up of vehicle |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20090805 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20111013 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120306 |
|
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: 20121106 |
|
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: 20121109 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20151116 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |