JP6868530B2

JP6868530B2 - シミュレーション方法

Info

Publication number: JP6868530B2
Application number: JP2017191939A
Authority: JP
Inventors: 山崎　勝; 勝山崎; 宮崎　泰三; 泰三宮崎
Original assignee: Hitachi Astemo Ltd
Current assignee: Hitachi Astemo Ltd
Priority date: 2017-09-29
Filing date: 2017-09-29
Publication date: 2021-05-12
Anticipated expiration: 2037-09-29
Also published as: JP2019068621A; WO2019064832A1

Description

本発明は、コンピューターを用いたシミュレーション方法に関する。

本技術分野の背景技術として、下記特許文献１が知られている。特許文献１に記載の技術では、モデル選択部が、条件入力部から設定された選択条件に基づいて、シミュレーションモデルを選択し、当該シミュレーションモデルをモデルデータベースから読み出す。シミュレーション計算部が、この読み出されたシミュレーションモデルを用い、条件入力部に設定された初期状態及びシミュレーション条件に基づいて、シミュレーション計算を行う。これにより、モデルの選択条件に基づいて、詳細度の異なる各シミュレーションモデルを切り替えてシミュレーション計算を行う。例えば、重要な部分は詳細度の高いモデルを使って高い精度のシミュレーションを行い、あまり重要でない部分は詳細度の低いモデルを使って短時間でシミュレーションを行う。

特開２００２−２５９８８８号公報

特許文献１には、計算途中にシミュレーションモデルを切り替える際の手法が記載されている。しかしながら、特許文献１の切り替え手法では、シミュレーションに用いられる変数値を切り替えの前後で変換するためのプログラムを別途作成する必要がある。このプログラムの作成には、未知の変数値を推定する必要があるなどの理由から、時間がかかることがある。そのため、シミュレーター作成から計算結果の評価までのシミュレーション全体の時間を長大化させてしまい、結果として、シミュレーションを利用したシステム開発効率が低下する場合がある。

本発明によるシミュレーション方法は、コンピューターを用いてシミュレーション対象に関する計算を行うものであって、前記シミュレーション対象は、第１要素および第２要素と、前記第１要素と前記第２要素の間に介在する中間要素と、を含み、前記コンピューターを、前記第１要素をモデル化した第１モデルと、前記第２要素をモデル化した第２モデルと、同一の前記中間要素を異なる計算モデルでそれぞれモデル化した第１中間モデルおよび第２中間モデルと、としてそれぞれ機能させ、前記第１モデルと前記第２モデルが前記第１中間モデルを介して情報を相互に送受信する第１モードと、前記第１モデルと前記第２モデルが前記第２中間モデルを介して情報を相互に送受信する第２モードと、のいずれかにより、前記コンピューターに前記計算を実行させ、前記第１モードでは、前記第１モデルから出力される第１パラメーターを前記第１中間モデルに入力し、入力された前記第１パラメーターに応じて前記第１中間モデルから出力される第２パラメーターを前記第２モデルに入力すると共に、前記第２モデルから出力される第３パラメーターを前記第１中間モデルに入力し、入力された前記第３パラメーターに応じて前記第１中間モデルから出力される第４パラメーターを前記第１モデルに入力するように、前記コンピューターを動作させ、前記第２モードでは、前記第１モデルから出力される第１パラメーターを前記第２中間モデルに入力し、入力された前記第１パラメーターに応じて前記第２中間モデルから出力される第２パラメーターを前記第２モデルに入力すると共に、前記第２モデルから出力される第３パラメーターを前記第２中間モデルに入力し、入力された前記第３パラメーターに応じて前記第２中間モデルから出力される第４パラメーターを前記第１モデルに入力するように、前記コンピューターを動作させ、前記第１モードから前記第２モードに切り替える際に、前記第１パラメーターまたは前記第３パラメーターを、前記第１中間モデルおよび前記第２中間モデルに入力し、入力された前記第１パラメーターまたは前記第３パラメーターに応じて前記第２中間モデルから出力される前記第２パラメーターまたは前記第４パラメーターと、入力された前記第１パラメーターまたは前記第３パラメーターに応じて前記第１中間モデルから出力される前記第２パラメーターまたは前記第４パラメーターと、の比較を行い、前記比較の結果が所定の基準を満たしたときに、前記第１モデルおよび前記第２モデルと前記第２中間モデルとを接続すると共に、前記第１モデルおよび前記第２モデルと前記第１中間モデルとを切り離すことで、前記第１モードから前記第２モードへの切り替えを実施するように、前記コンピューターを動作させる。

本発明によれば、システム開発効率の高いシミュレーション方法を提供できる。

本発明の第１の実施形態に係るシミュレーターの構成例を示す図である。本発明の第１の実施形態に関する従来のシミュレーターの構成例を示す図である。シミュレーターに用いられる各物理ドメインでのパラメーターの例を示した図である。本発明の第１の実施形態に係るシミュレーターを用いて計算されたローター回転速度の例を示す図である。本発明の第１の実施形態に係るシミュレーターにおいてモードを切り替える際の制御手順例を示したフロー図である。第１モードから第２モードへの切り替えを行ったときのローター回転速度Ｖ、Ｖｂの変化例を示す図である。本発明の第２の実施形態に係るシミュレーターの構成例を示す図である。本発明の第２の実施形態に係るシミュレーターにおいてモードを切り替える際の制御手順例を示したフロー図である。本発明の第３の実施形態に係るシミュレーターの構成例を示す図である。本発明の第４の実施形態に関する従来のシミュレーターの構成例を示す図である。本発明の第４の実施形態に係るシミュレーターの構成例を示す図である。

（第１の実施形態）
本実施形態では、電気自動車の駆動装置のシステムシミュレーターへの適用例を説明する。

図１は、本発明の第１の実施形態に係るシミュレーターの構成例を示す図である。図１に示すシミュレーターは、電気自動車の駆動装置をシミュレーション対象としたシステムシミュレーターであり、たとえばコンピューターを用いて実現される。このシミュレーターは、モーターモデル１の発生する電磁トルク２１を高速ギヤモデル３を通して車軸トルク２３に変換し、車両モデル２に伝達する。また、車両モデル２の運動状態に応じた車軸回転速度２４を高速ギヤモデル３を通してローター回転速度２２に変換し、モーターモデル１に入力する。このときのモーターモデル１、高速ギヤモデル３および車両モデル２の運動をそれぞれ計算することで、電気自動車の駆動装置の動きをシミュレーションする。すなわち、本実施形態のシミュレーターでは、シミュレーション対象である電気自動車の駆動装置を構成する各要素（モーター、ギヤ、車両）を、モーターモデル１、高速ギヤモデル３、車両モデル２としてそれぞれモデル化し、これらのモデル間で相互に影響するパラメーターを互いに受け渡すことで、各要素の連成運動を計算することを可能としている。これにより、シミュレーション対象である電気自動車の駆動装置に関する計算を行うシミュレーターを実現している。

ここで図１に示すように、本実施形態のシミュレーターでは、さらに高精度ギヤモデル３ｂが高速ギヤモデル３と並列に配置されている。高精度ギヤモデル３ｂは、高速ギヤモデル３と同様に、シミュレーション対象である電気自動車の駆動装置を構成する要素の一つであるギヤをモデル化したものである。なお、高速ギヤモデル３は、ギヤの減速比によるトルクと回転速度の変換を解く計算内容をモデル化したものであり、計算負荷は軽いが軸振動の現象などを解くことが出来ないモデルである。一方、高精度ギヤモデル３ｂは、ギヤ内部の軸ねじれを考慮した減速比によるトルクと回転速度の変換を解く計算内容をモデル化したものであり、高速ギヤモデル３の計算内容に加えて、内部のギヤの挙動やドライブシャフトの捩れなどの計算も含んでいる。すなわち高精度ギヤモデル３ｂは、計算負荷は重いが軸振動の現象などを解くことができるモデルである。

なお、図１では、高速ギヤモデル３に入出力される各パラメーターと区別するため、高精度ギヤモデル３ｂに入出力される各パラメーターの符号の末尾に「ｂ」を付している。すなわち、高精度ギヤモデル３ｂは、モーターモデル１から入力された電磁トルク２１ｂを車軸トルク２３ｂに変換し、車両モデル２に伝達する。また、車両モデル２から入力された車軸回転速度２４ｂをローター回転速度２２ｂに変換し、モーターモデル１に伝達する。

本実施形態のシミュレーターにおいて、高速ギヤモデル３と高精度ギヤモデル３ｂのどちらを計算に用いるかは、各モデル間に配置されたスイッチ３１，３２，３３，３４，３１ｂ，３２ｂ，３３ｂ，３４ｂをそれぞれ切り替えることで決定される。図１では、モーターモデル１および車両モデル２と高速ギヤモデル３の間にそれぞれ配置されたスイッチ３１，３２，３３，３４が接続状態である一方で、モーターモデル１および車両モデル２と高精度ギヤモデル３ｂの間にそれぞれ配置されたスイッチ３１ｂ，３２ｂ，３３ｂ，３４ｂが切断状態である。すなわち、図１では、モーターモデル１と車両モデル２が高速ギヤモデル３を介して接続されており、高速ギヤモデル３を用いた計算が実施される状態となっている。以下では、高速ギヤモデル３を用いて計算を行う場合を「第１モード」と称する。一方、各スイッチの状態を入れ換え、モーターモデル１と車両モデル２が高精度ギヤモデル３ｂを介して接続されるようにすると、高精度ギヤモデル３ｂを用いた計算が実施される状態となる。以下では、高精度ギヤモデル３ｂを用いて計算を行う場合を「第２モード」と称する。

上記各スイッチの制御は、切り替え制御器４で実施される。切り替え制御器４と各スイッチとは不図示の配線により接続されており、切り替え制御器４から各スイッチを操作できるようになっている。なお、切り替え制御器４と各スイッチ間の配線は、物理的な電線ではなく、コンピューター上で仮想的に実現されるものである。

切り替え制御器４には、高速ギヤモデル３から出力されるローター回転速度２２と、高精度ギヤモデル３ｂから出力されるローター回転速度２２ｂとが入力される。また、モーターモデル１から出力される電磁トルク２１、２１ｂも切り替え制御器４に入力される。切り替え制御器４は、これらのパラメーターに基づいて、各スイッチの操作タイミングを決定する。なお、各スイッチの操作タイミングの具体的な決定方法については、後で説明する。

本実施形態のシミュレーターでは、以上説明したような各構成がコンピューター上で実現されている。すなわち、本実施形態のシミュレーターでは、コンピューターを、モーターをモデル化したモーターモデル１と、車両をモデル化した車両モデル２と、モーターと車両の間に介在するギヤをそれぞれモデル化した高速ギヤモデル３および高精度ギヤモデル３ｂと、としてそれぞれ機能させている。そして、モーターモデル１と車両モデル２が高速ギヤモデル３を介して情報を送受信する第１モードと、モーターモデル１と車両モデル２が高精度ギヤモデル３ｂを介して情報を送受信する第２モードと、のいずれかにより、コンピューターに計算を実行させる。具体的には、第１モードでは、モーターモデル１から出力される電磁トルク２１を高速ギヤモデル３に入力し、入力された電磁トルク２１に応じて高速ギヤモデル３から出力される車軸トルク２３を車両モデル２に入力すると共に、車両モデル２から出力される車軸回転速度２４を高速ギヤモデル３に入力し、入力された車軸回転速度２４に応じて高速ギヤモデル３から出力されるローター回転速度２２をモーターモデル１に入力するように、コンピューターを動作させる。また、第２モードでは、モーターモデル１から出力される電磁トルク２１ｂを高精度ギヤモデル３ｂに入力し、入力された電磁トルク２１ｂに応じて高精度ギヤモデル３ｂから出力される車軸トルク２３ｂを車両モデル２に入力すると共に、車両モデル２から出力される車軸回転速度２４ｂを高精度ギヤモデル３ｂに入力し、入力された車軸回転速度２４ｂに応じて高精度ギヤモデル３ｂから出力されるローター回転速度２２ｂをモーターモデル１に入力するように、コンピューターを動作させる。

次に、本実施形態のシミュレーターにおいて各モデル間で受け渡されるパラメーターについて説明する。なお、以下では理解を簡単にするため、図２に示す従来のシミュレーターを用いて、各パラメーターについて解説する。

図２は、本発明の第１の実施形態に関する従来のシミュレーターの構成例を示す図である。図２に示す従来のシミュレーターでは、モーターモデル１と車両モデル２が高速ギヤモデル３のみを介して接続されており、図１の高精度ギヤモデル３ｂが存在していない。そのため、図１で示した高速ギヤモデル３と高精度ギヤモデル３ｂとを切り替えるための各スイッチは、図２には設けられていない。

図２のシミュレーターでは、モーターモデル１と高速ギヤモデル３の間で入出力される電磁トルク２１およびローター回転速度２２と、高速ギヤモデル３と車両モデル２の間で入出力される車軸トルク２３および車軸回転速度２４とが、各モデル間で受け渡されるパラメーターに該当する。なお、図２の矢印の向きは、各パラメーターによる情報の伝達方向を示している。このように、従来のシミュレーターでは、各モデル間での情報伝達をパラメーターの受け渡しによって実現している。この点は、図１に示した本実施形態のシミュレーターについても同様である。

図３は、シミュレーターに用いられる各物理ドメインでのパラメーターの例を示した図である。シミュレーターで用いられる各種パラメーターは、アクロス変数とスルー変数の２種類に分類される。スルー変数とは、物理ドメインが電気である場合のキルヒホッフの第一法則や、他の物理ドメインに対して成り立つ同様の法則、すなわち一点に流れ込む量の和がゼロとなる法則を満たす変数である。一方、アクロス変数とは，物理ドメインが電気である場合のキルヒホッフの第二法則や、他の物理ドメインに対して成り立つ同様の法則、すなわちモデル内でループを形成する要素の積算値がゼロとなる法則を満たす変数である。これらの変数は、同一の物理ドメインでそれぞれ対応関係にあるものが対にして用いられる。

具体的には、図３に示すように、たとえば物理ドメインが機械回転の場合には、アクロス変数である角速度（回転速度）と、スルー変数であるトルクとが、シミュレーターにおいて対になるパラメーターとして利用される。同様に、物理ドメインが機械並進の場合には、アクロス変数である速度と、スルー変数である力とが、シミュレーターにおいて対になるパラメーターとして利用される。また、物理ドメインが電気の場合には、アクロス変数である電圧と、スルー変数である電流とが、シミュレーターにおいて対になるパラメーターとして利用され、物理ドメインが油や空圧の場合には、アクロス変数である圧力と、スルー変数である流量とが、シミュレーターにおいて対になるパラメーターとして利用される。これにより、物理ドメインに応じてシミュレーター内に設定された各モデル同士を連結し、シミュレーション対象を構成する各要素の連成運動の計算を行うことができるようにしている。なお、図３に示した物理ドメインと対応するアクロス変数およびスルー変数とは一例であり、本実施形態のシミュレーターで利用されるパラメーターはこれに限定されない。

次に、本実施形態のシミュレーターにおいて計算途中にモードを切り替える場合の動作について解説する。以下の説明では、電気自動車を停止している状態から加速させ、所定の速度に達したらブレーキをかけて減速し停止する走行パターンについて、電力の消費量と減速時に発生する軸振動とを計算する場合での、本実施形態のシミュレーターの動作例を説明する。

図４は、本発明の第１の実施形態に係るシミュレーターを用いて計算されたローター回転速度の例を示している。図４のグラフ４１は、電気自動車の加速を開始してからブレーキ開始直後までの前半部分を、高速ギヤモデル３を用いた第１モードにより計算を行い、ブレーキによる減速を開始してから電気自動車が停止するまでの後半部分を、高精度ギヤモデル３ｂを用いた第２モードにより計算を行った場合での、ローター回転速度の時間変化の様子を示している。なお、図４のグラフ４１では、第１モードから第２モードに切り替える際の時間経過を０としているが、実際には後述するように、高精度ギヤモデル３ｂの状態量が切り替え可能な状態へと変化するまでにある程度の時間を要する。したがって、この時間を加味してローター回転速度の時間変化を求めてもよい。

上記の走行パターンにおけるローター回転速度を計算する場合、前半の加速部分では、モーターの発生する電磁トルクがどれだけ車両を加速させるかを解くことが課題である。そのため、高速ギヤモデル３を用いた計算で十分である。一方、後半の減速部分では、軸振動が発生するため、この挙動を算出することも課題となる。そのため、高速ギヤモデル３を用いた計算では不十分であり、高精度ギヤモデル３ｂを用いた計算を行うことが必須となる。

なお、前半部分を含めた全区間で高精度ギヤモデル３ｂを用いた計算を行っても、所望の結果を得ることができる。しかし、高精度ギヤモデル３ｂを用いた計算は時間がかかるため、計算コストを含めて考えた場合、全区間で高精度ギヤモデル３ｂを用いた計算を行うのはあまり好ましくない。そこで、本実施形態のシミュレーターのように、計算途中で第１モードから第２モードに切り替えて計算を行うことにより、効率的なシミュレーションを実現できる。

図２に示した従来のシミュレーターにおいて、第１モードから第２モードへの切り替えを実現する場合には、切り替え前後での変数値の変換を行った上で、高速ギヤモデル３を高精度ギヤモデル３ｂに置き換える必要がある。具体的には、前半の第１モードでの計算が終了した後に、高速ギヤモデル３内の各状態量を表す変数値、たとえばローターの回転速度、ドライブシャフトの回転速度、内部の伝達トルクなどを、高精度ギヤモデル３ｂ内の各状態量を表す変数値に変換する。そして、変換されたこれらの変数値を高精度ギヤモデル３ｂに設定して、後半の第２モードでの計算を継続すればよい。

しかし、上記のような変数値の変換計算を実現するためには、高速ギヤモデル３と高精度ギヤモデル３ｂそれぞれの物理特性や拘束条件を加味した変換プログラムを、シミュレーターとは別に作成する必要がある。また、高速ギヤモデル３から高精度ギヤモデル３ｂへの切り替えでは、内部の状態量の数が増加するため、単純な変換式では高精度ギヤモデル３ｂ内の各状態量を表す変数値を求めることができない。そのため、過去の計算履歴を参照して状態量を推定したり、最小二乗法等の手法によりカーブフィッティングを行って状態量を推定したりする作業が必要となる。さらに、高速ギヤモデル３や高精度ギヤモデル３ｂの内容が変更された場合には、それに応じて変換プログラムも作り直す必要が生じる。これらの原因から、変換プログラムの作成に時間がかかってしまい、その結果、シミュレーションを利用したシステム開発効率が低下する場合がある。

一方、本実施形態のシミュレーターでは、図１のように高精度ギヤモデル３ｂを高速ギヤモデル３と並列に配置して、各モデル間に配置されたスイッチを適切な操作タイミングで切り替える。これにより、変数値の変換を必要とせずに、計算途中でのモード切り替えを実現している。

図５は、本発明の第１の実施形態に係るシミュレーターにおいてモードを切り替える際の制御手順例を示したフロー図である。以下では、図５のフロー図に従って、本実施形態のシミュレーターの動作を説明する。なお、図１のシミュレーターにおいて、切り替え制御器４は、予め設定されたモード切り替えのタイミングになると、図５のフロー図に示す処理を開始する。

まずステップＳ５１において、切り替え制御器４は、スイッチ３４ｂを切断状態から接続状態に切り替えて、高精度ギヤモデル３ｂへ車軸回転速度２４ｂの入力を接続する。すなわち、それまでは高速ギヤモデル３のみに入力されていた車両モデル２からの車軸回転速度２４を、高速ギヤモデル３と高精度ギヤモデル３ｂに車軸回転速度２４、２４ｂとしてそれぞれ入力されるようにする。これにより、高速ギヤモデル３での計算結果を用いて、高精度ギヤモデル３ｂにおける内部の状態量の変動が開始される。

前述の図３を参照すると、車軸回転速度２４、２４ｂは、物理ドメインが機械回転の場合におけるアクロス変数に相当する。すなわちステップＳ５１では、アクロス変数である車軸回転速度２４、２４ｂを、モードを切り替える際に最初に高速ギヤモデル３と高精度ギヤモデル３ｂに接続するパラメーターとすることで、高速ギヤモデル３内部のアクロス変数と高精度ギヤモデル３ｂ内部のアクロス変数とを合致させるようにしている。ここで、アクロス変数とは、前述のように角速度、速度、電圧、圧力などであり、これらは内部エネルギーと密接な関係がある。すなわち、上記のようにステップＳ５１で高速ギヤモデル３と高精度ギヤモデル３ｂの内部のアクロス変数を合致させることは、これらのモデルの内部エネルギーを合致させることと近似である。したがって、高精度ギヤモデル３ｂ内部の状態量を切り替え可能な状態へと迅速に変化させることができる。

なお、モードを切り替える際に最初に接続するパラメーターは、必ずしもアクロス変数に限定されるものではない。すなわち、モードを切り替える際に最初に接続するパラメーターをアクロス変数ではなくスルー変数としてもよいし、アクロス変数とスルー変数とを混在させてもよい。

次にステップＳ５２において、切り替え制御器４は、入力された車軸回転速度２４に応じて高速ギヤモデル３から出力されるローター回転速度２２と、入力された車軸回転速度２４ｂに応じて高精度ギヤモデル３ｂから出力されるローター回転速度２２ｂとを検出する。なお、以下の説明では、高速ギヤモデル３からのローター回転速度２２の検出値をローター回転速度Ｖとし、高精度ギヤモデル３ｂからのローター回転速度２２ｂの検出値をローター回転速度Ｖｂとする。

次にステップＳ５３において、切り替え制御器４は、ステップＳ５２で検出したローター回転速度Ｖとローター回転速度Ｖｂとを比較する。ここでは、ローター回転速度Ｖｂとローター回転速度Ｖとの差分δｖを計算することで、これらの比較を行う。

次にステップＳ５４において、切り替え制御器４は、ステップＳ５３で算出した差分δｖに基づいて、高精度ギヤモデル３ｂに対する制御量Ｔａを計算する。ここでは、モーターモデル１から出力された電磁トルク２１の値を検出し、その電磁トルク２１の値に応じて、差分δｖの絶対値が縮小するように制御量Ｔａを決定する。このとき、たとえば一般的なＰＩＤ制御を用いて制御量Ｔａを演算してもよいし、予め設定されたマップ情報などを用いて制御量Ｔａを決定してもよい。少なくとも差分δｖの絶対値を縮小できるものであれば、どのような方法で制御量Ｔａを決定してもよい。

次にステップＳ５５において、切り替え制御器４は、ステップＳ５４で計算した制御量Ｔａを、高精度ギヤモデル３ｂの電磁トルク入力へ供給する。すなわち、制御量Ｔａに応じた電磁トルク２１ｂを、切り替え制御器４からスイッチ３１ｂを介して高精度ギヤモデル３ｂに入力する。なお、スイッチ３１ｂは、モーターモデル１と高精度ギヤモデル３ｂとの間が切断状態であるときには、切り替え制御器４と高精度ギヤモデル３ｂとの間が接続状態になっている。これにより、車軸回転速度２４ｂに加えて、差分δｖに応じて調整した電磁トルク２１ｂをさらに高精度ギヤモデル３ｂに入力して、差分δｖの絶対値が小さくなるように高精度ギヤモデル３ｂを動作させる。

次にステップＳ５６において、切り替え制御器４は、ステップＳ５２で検出したローター回転速度Ｖとローター回転速度Ｖｂとがおおよそ等しいか否かを判定する。ここでは、たとえばステップＳ５３で計算した差分δｖの絶対値が所定の基準値以下であるか否かを判断する。その結果、差分δｖの絶対値が基準値以下であれば、ローター回転速度Ｖとローター回転速度Ｖｂとがおおよそ等しいと判定し、ステップＳ５７に進む。一方、差分δｖの絶対値が基準値を超えている場合は、ローター回転速度Ｖとローター回転速度Ｖｂとがおおよそ等しくないと判定し、ステップＳ５２に戻ってステップＳ５２以降の処理を再度実行する。

ステップＳ５７において、切り替え制御器４は、高速ギヤモデル３の全端子、すなわちスイッチ３１，３２，３３，３４を接続状態から切断状態に切り替えると共に、高精度ギヤモデル３ｂの全端子、すなわちステップＳ５１で接続済みのスイッチ３４ｂを除いた残りのスイッチ３１ｂ，３２ｂ，３３ｂを切断状態から接続状態に切り替える。これにより、モーターモデル１および車両モデル２と高精度ギヤモデル３ｂとを接続すると共に、モーターモデル１および車両モデル２と高速ギヤモデル３とを切り離して、第１モードから第２モードへの切り替えが実施されるようにする。ステップＳ５７の処理を実行したら、切り替え制御器４は図５のフロー図を終了する。

図６は、図５のフロー図に従って第１モードから第２モードへの切り替えを行ったときのローター回転速度Ｖ、Ｖｂの変化例を示す図である。図６において、グラフ６１はローター回転速度Ｖの変化の様子を示し、グラフ６２はローター回転速度Ｖｂの変化の様子を示している。符号６３に示すタイミングで図５のフロー図によりモード切り替えを開始すると、グラフ６２に示すローター回転速度Ｖｂの値が、次第にグラフ６１に示すローター回転速度Ｖの値に近づいていく。そして、ある一定の計算時間が経過した符号６４に示すタイミングで、ローター回転速度Ｖとローター回転速度Ｖｂとがほぼ同一の値になると、図５のステップＳ５６が肯定判定され、第１モードから第２モードへの切り替えが終了する。

以上説明した図５のフロー図に示す処理を切り替え制御器４が行うことで、第１モードから第２モードへの切り替えが実施される。第２モードへの切り替え後は、モーターモデル１からの電磁トルク２１ｂが高精度ギヤモデル３ｂに入力されて車軸トルク２３ｂに変換され、車両モデル２に伝達される。また、車両モデル２からの車軸回転速度２４ｂが高精度ギヤモデル３ｂに入力されてローター回転速度２２ｂに変換され、モーターモデル１に伝達される。したがって、前述のような変数値の変換を行わずに、計算途中でのモード切り替えを実現できる。

以上説明した本発明の第１の実施形態によれば、以下の作用効果を奏する。

（１）本実施形態のシミュレーターは、コンピューターを用いて、シミュレーション対象である電気自動車の駆動装置に関する計算を行う。このシミュレーション対象は、モーターおよび車両と、モーターと車両の間に介在するギヤとを含む。本実施形態のシミュレーターでは、コンピューターを、モーターをモデル化したモーターモデル１と、車両をモデル化した車両モデル２と、ギヤをそれぞれモデル化した高速ギヤモデル３および高精度ギヤモデル３ｂと、としてそれぞれ機能させる。そして、モーターモデル１と車両モデル２が高速ギヤモデル３を介して情報を送受信する第１モードと、モーターモデル１と車両モデル２が高精度ギヤモデル３ｂを介して情報を送受信する第２モードと、のいずれかにより、コンピューターに計算を実行させる。第１モードでは、モーターモデル１から出力される電磁トルク２１を高速ギヤモデル３に入力し、入力された電磁トルク２１に応じて高速ギヤモデル３から出力される車軸トルク２３を車両モデル２に入力すると共に、車両モデル２から出力される車軸回転速度２４を高速ギヤモデル３に入力し、入力された車軸回転速度２４に応じて高速ギヤモデル３から出力されるローター回転速度２２をモーターモデル１に入力するように、コンピューターを動作させる。第２モードでは、モーターモデル１から出力される電磁トルク２１ｂを高精度ギヤモデル３ｂに入力し、入力された電磁トルク２１ｂに応じて高精度ギヤモデル３ｂから出力される車軸トルク２３ｂを車両モデル２に入力すると共に、車両モデル２から出力される車軸回転速度２４ｂを高精度ギヤモデル３ｂに入力し、入力された車軸回転速度２４ｂに応じて高精度ギヤモデル３ｂから出力されるローター回転速度２２ｂをモーターモデル１に入力するように、コンピューターを動作させる。さらに、第１モードから第２モードに切り替える際に、車軸回転速度２４、２４ｂを、高速ギヤモデル３および高精度ギヤモデル３ｂにそれぞれ入力し（ステップＳ５１）、入力された車軸回転速度２４ｂに応じて高精度ギヤモデル３ｂから出力されるローター回転速度２２ｂと、入力された車軸回転速度２４に応じて高速ギヤモデル３から出力されるローター回転速度２２と、の比較を行い（ステップＳ５３）、その比較の結果が所定の基準を満たしたときに（ステップＳ５６：Ｙｅｓ）、モーターモデル１および車両モデル２と高精度ギヤモデル３ｂとを接続すると共に、モーターモデル１および車両モデル２と高速ギヤモデル３とを切り離すことで、第１モードから第２モードへの切り替えを実施する（ステップＳ５７）ように、コンピューターを動作させる。このようにしたので、変数値の変換を行わずに計算途中でのモード切り替えを実現でき、その結果、システム開発効率の高いシミュレーション方法を提供できる。

（２）ステップＳ５３の比較では、高精度ギヤモデル３ｂから出力されるローター回転速度２２ｂの値Ｖｂと、高速ギヤモデル３から出力されるローター回転速度２２の値Ｖとの差分δｖを算出する。そして、算出した差分δｖの絶対値が所定の基準値以下となったときに、ステップＳ５６を肯定判定してステップＳ５７を実行し、モーターモデル１および車両モデル２と高精度ギヤモデル３ｂとを接続すると共に、モーターモデル１および車両モデル２と高速ギヤモデル３とを切り離すことで、第１モードから第２モードへの切り替えを実施するように、コンピューターを動作させる。このようにしたので、適切なタイミングで第１モードから第２モードへの切り替えを行うことができる。

（３）本実施形態のシミュレーターでは、第１モードから第２モードへの切り替えの際に、高精度ギヤモデル３ｂに入力される電磁トルク２１ｂの値を調整する（ステップＳ５４、Ｓ５５）ことで差分δｖが小さくなるように、コンピューターを動作させる。具体的には、高精度ギヤモデル３ｂから出力されるローター回転速度２２ｂの値Ｖｂと高速ギヤモデル３から出力されるローター回転速度２２の値Ｖとの差分δｖに応じて調整した制御量ＴａをステップＳ５４で計算し、この制御量Ｔａに応じた電磁トルク２１ｂをステップＳ５５でさらに高精度ギヤモデル３ｂに入力することで、差分δｖが小さくなるようにコンピューターを動作させる。このようにしたので、第１モードから第２モードへの切り替えをスムーズに行うことができる。

（４）なお、第１モードから第２モードに切り替える際に、ステップＳ５１で高速ギヤモデル３および高精度ギヤモデル３ｂにそれぞれ入力される車軸回転速度２４、２４ｂは、アクロス変数である。したがって、高速ギヤモデル３と高精度ギヤモデル３ｂとの内部エネルギーを合致させ、高精度ギヤモデル３ｂ内部の状態量を切り替え可能な状態へと迅速に変化させることができる。

なお、本実施形態では、以下の変形例も可能である。

（変形例１）
図５のステップＳ５５で制御量Ｔａに応じた電磁トルク２１ｂを高精度ギヤモデル３ｂに入力する際に、その値を徐々に変化させてもよい。具体的には、電磁トルク２１ｂの単位時間当たり、たとえば演算周期当たりの入力変化量を所定値以下に制限することで、高精度ギヤモデル３ｂに入力される電磁トルク２１ｂの値が急激に変化しないようにする。このようにすれば、高精度ギヤモデル３ｂ内部の状態量が急変するのを防いで、第１モードから第２モードへの切り替えをより適切に行うことができる。

（変形例２）
本実施形態では、第１モードから第２モードへの切り替えを行う前に、図５のステップＳ５３〜Ｓ５５の処理を行うことで、高精度ギヤモデル３ｂから出力されるローター回転速度２２ｂの値Ｖｂと、高速ギヤモデル３から出力されるローター回転速度２２の値Ｖとが迅速に近づくようにしたが、これに限定するものではない。たとえば、ステップＳ５３〜Ｓ５５の処理を実施しなくてもよい。

（変形例３）
本実施形態では、図５のステップＳ５６において、ローター回転速度Ｖとローター回転速度Ｖｂとがおおよそ等しいか否かを判定することで、第１モードから第２モードへの切り替えが可能であるか否かを判定したが、判定方法はこれに限定するものではない。たとえば、所定の計算時間が経過することで判定を行ってもよいし、複数の方法を併用して判定を行ってもよい。

なお、本実施形態では、第１モードから第２モードへの切り替えを行う際の手順を説明したが、第２モードから第１モードへの切り替えについても、高速ギヤモデル３と高精度ギヤモデル３ｂとを入れ替えることで、同様の手順により行うことが可能である。すなわち、本実施形態で説明したモード切り替えの手法は、第１モードから第２モードへの切り替え、第２モードから第１モードへの切り替えのどちらにも適用可能である。

（第２の実施形態）
次に本発明の第２の実施形態について説明する。本実施形態では、アクロス変数である車軸回転速度２４、２４ｂではなく、スルー変数である電磁トルク２１、２１ｂを最初に高速ギヤモデル３と高精度ギヤモデル３ｂに接続して、モード切り替えを行う例を説明する。

図７は、本発明の第２の実施形態に係るシミュレーターの構成例を示す図である。図７に示すシミュレーターは、第１の実施形態で説明した図１のシミュレーターと同様に、電気自動車の駆動装置をシミュレーション対象としたシステムシミュレーターであり、たとえばコンピューターを用いて実現される。本実施形態のシミュレーターは、第１の実施形態のシミュレーターと比較して、切り替え制御器４に入出力されるパラメーターの種類が異なっている。具体的には、第１の実施形態では、高速ギヤモデル３から出力されるローター回転速度２２と、高精度ギヤモデル３ｂから出力されるローター回転速度２２ｂと、モーターモデル１から出力される電磁トルク２１とが、切り替え制御器４に入力され、切り替え制御器４から高速ギヤモデル３と高精度ギヤモデル３ｂに対して電磁トルク２１、２１ｂがそれぞれ出力されていた。一方、本実施形態では、高速ギヤモデル３から出力される車軸トルク２３と、高精度ギヤモデル３ｂから出力される車軸トルク２３ｂと、車両モデル２から出力される車軸回転速度２４とが、切り替え制御器４に入力され、切り替え制御器４から高速ギヤモデル３と高精度ギヤモデル３ｂに対して車軸回転速度２４、２４ｂがそれぞれ出力されている。これ以外の点では、本実施形態のシミュレーターは、第１の実施形態のシミュレーターと同様である。

図８は、本発明の第２の実施形態に係るシミュレーターにおいてモードを切り替える際の制御手順例を示したフロー図である。以下では、図８のフロー図に従って、本実施形態のシミュレーターの動作を説明する。なお、図７のシミュレーターにおいて、切り替え制御器４は、予め設定されたモード切り替えのタイミングになると、図８のフロー図に示す処理を開始する。

まずステップＳ７１において、切り替え制御器４は、スイッチ３１ｂを切断状態から接続状態に切り替えて、高精度ギヤモデル３ｂへモータートルク、すなわち電磁トルク２１ｂの入力を接続する。すなわち、それまでは高速ギヤモデル３のみに入力されていたモーターモデル１からの電磁トルク２１を、高速ギヤモデル３と高精度ギヤモデル３ｂに電磁トルク２１、２１ｂとしてそれぞれ入力されるようにする。これにより、高速ギヤモデル３での計算結果を用いて、高精度ギヤモデル３ｂにおける内部の状態量の変動が開始される。

次にステップＳ７２において、切り替え制御器４は、入力された電磁トルク２１に応じて高速ギヤモデル３から出力される車軸トルク２３と、入力された電磁トルク２１ｂに応じて高精度ギヤモデル３ｂから出力される車軸トルク２３ｂとを検出する。なお、以下の説明では、高速ギヤモデル３からの車軸トルク２３の検出値を車軸トルクＴとし、高精度ギヤモデル３ｂからの車軸トルク２３ｂの検出値を車軸トルクＴｂとする。

次にステップＳ７３において、切り替え制御器４は、ステップＳ７２で検出した車軸トルクＴと車軸トルクＴｂとを比較する。ここでは、車軸トルクＴｂと車軸トルクＴとの差分δｔを計算することで、これらの比較を行う。

次にステップＳ７４において、切り替え制御器４は、ステップＳ７３で算出した差分δｔに基づいて、高精度ギヤモデル３ｂに対する制御量Ｖａを計算する。ここでは、車両モデル２から出力された車軸回転速度２４の値を検出し、その車軸回転速度２４の値に応じて、差分δｔの絶対値が縮小するように制御量Ｖａを決定する。このとき、たとえば一般的なＰＩＤ制御を用いて制御量Ｖａを演算してもよいし、予め設定されたマップ情報などを用いて制御量Ｖａを決定してもよい。少なくとも差分δｔの絶対値を縮小できるものであれば、どのような方法で制御量Ｖａを決定してもよい。

次にステップＳ７５において、切り替え制御器４は、ステップＳ７４で計算した制御量Ｖａを、高精度ギヤモデル３ｂの車軸回転速度入力へ供給する。すなわち、制御量Ｖａに応じた車軸回転速度２４ｂを、切り替え制御器４からスイッチ３４ｂを介して高精度ギヤモデル３ｂに入力する。なお、スイッチ３４ｂは、車両モデル２と高精度ギヤモデル３ｂとの間が切断状態であるときには、切り替え制御器４と高精度ギヤモデル３ｂとの間が接続状態になっている。これにより、電磁トルク２１ｂに加えて、差分δｔに応じて調整した車軸回転速度２４ｂをさらに高精度ギヤモデル３ｂに入力して、差分δｔの絶対値が小さくなるように高精度ギヤモデル３ｂを動作させる。

次にステップＳ７６において、切り替え制御器４は、ステップＳ７２で検出した車軸トルクＴと車軸トルクＴｂとがおおよそ等しいか否かを判定する。ここでは、たとえばステップＳ７３で計算した差分δｔの絶対値が所定の基準値以下であるか否かを判断する。その結果、差分δｔの絶対値が基準値以下であれば、車軸トルクＴと車軸トルクＴｂとがおおよそ等しいと判定し、ステップＳ７７に進む。一方、差分δｔの絶対値が基準値を超えている場合は、車軸トルクＴと車軸トルクＴｂとがおおよそ等しくないと判定し、ステップＳ７２に戻ってステップＳ７２以降の処理を再度実行する。

ステップＳ７７において、切り替え制御器４は、高速ギヤモデル３の全端子、すなわちスイッチ３１，３２，３３，３４を接続状態から切断状態に切り替えると共に、高精度ギヤモデル３ｂの全端子、すなわちステップＳ７１で接続済みのスイッチ３１ｂを除いた残りのスイッチ３２ｂ，３３ｂ，３４ｂを切断状態から接続状態に切り替える。これにより、モーターモデル１および車両モデル２と高精度ギヤモデル３ｂとを接続すると共に、モーターモデル１および車両モデル２と高速ギヤモデル３とを切り離して、第１モードから第２モードへの切り替えが実施されるようにする。ステップＳ７７の処理を実行したら、切り替え制御器４は図８のフロー図を終了する。

以上説明した本発明の第２の実施形態によれば、以下の作用効果を奏する。

（１）本実施形態のシミュレーターでは、第１モードと第２モードのそれぞれにおいて、第１の実施形態と同様の動作をコンピューターに行わせる。さらに、第１モードから第２モードに切り替える際に、電磁トルク２１、２１ｂを、高速ギヤモデル３および高精度ギヤモデル３ｂにそれぞれ入力し（ステップＳ７１）、入力された電磁トルク２１ｂに応じて高精度ギヤモデル３ｂから出力される車軸トルク２３ｂと、入力された電磁トルク２１に応じて高速ギヤモデル３から出力される車軸トルク２３と、の比較を行い（ステップＳ７３）、その比較の結果が所定の基準を満たしたときに（ステップＳ７６：Ｙｅｓ）、モーターモデル１および車両モデル２と高精度ギヤモデル３ｂとを接続すると共に、モーターモデル１および車両モデル２と高速ギヤモデル３とを切り離すことで、第１モードから第２モードへの切り替えを実施する（ステップＳ７７）ように、コンピューターを動作させる。このようにしたので、第１の実施形態で説明したのと同様に、変数値の変換を行わずに計算途中でのモード切り替えを実現でき、その結果、システム開発効率の高いシミュレーション方法を提供できる。

（２）ステップＳ７３の比較では、高精度ギヤモデル３ｂから出力される車軸トルク２３ｂの値Ｔｂと、高速ギヤモデル３から出力される車軸トルク２３の値Ｔとの差分δｔを算出する。そして、算出した差分δｔの絶対値が所定の基準値以下となったときに、ステップＳ７６を肯定判定してステップＳ７７を実行し、モーターモデル１および車両モデル２と高精度ギヤモデル３ｂとを接続すると共に、モーターモデル１および車両モデル２と高速ギヤモデル３とを切り離すことで、第１モードから第２モードへの切り替えを実施するように、コンピューターを動作させる。このようにしたので、第１の実施形態で説明したのと同様に、適切なタイミングで第１モードから第２モードへの切り替えを行うことができる。

（３）本実施形態のシミュレーターでは、第１モードから第２モードへの切り替えの際に、高精度ギヤモデル３ｂに入力される車軸回転速度２４ｂの値を調整する（ステップＳ７４、Ｓ７５）ことで差分δｔが小さくなるように、コンピューターを動作させる。具体的には、高精度ギヤモデル３ｂから出力される車軸トルク２３ｂの値Ｔｂと高速ギヤモデル３から出力される車軸トルク２３の値Ｔとの差分δｔに応じて調整した制御量ＶａをステップＳ７４で計算し、この制御量Ｖａに応じた車軸回転速度２４ｂをステップＳ７５でさらに高精度ギヤモデル３ｂに入力することで、差分δｔが小さくなるようにコンピューターを動作させる。このようにしたので、第１の実施形態で説明したのと同様に、第１モードから第２モードへの切り替えをスムーズに行うことができる。

なお、本実施形態においても、第１の実施形態で説明した各変形例が可能である。また、第２モードから第１モードへの切り替えについても、高速ギヤモデル３と高精度ギヤモデル３ｂとを入れ替えることで、同様の手順により行うことが可能である。

（第３の実施形態）
以上説明した第１および第２の実施形態では、単独のコンピューター上で実行されるシミュレーターの例を説明したが、本発明はこれに限定されず、複数のコンピューターに分散させてシミュレーターを実現してもよい。本実施形態では、それぞれ部分的に計算を実行する複数のコンピューターを用いたシミュレーターの例を説明する。

図９は、本発明の第３の実施形態に係るシミュレーターの構成例を示す図である。図９に示すシミュレーターは、第１の実施形態で説明した図１のシミュレーターと比較して、高速ギヤモデル３と高精度ギヤモデル３ｂが配置されていた部分にメイン通信Ｉ／Ｆ８３が配置されており、このメイン通信Ｉ／Ｆ８３と、モーターモデル１、車両モデル２および切り替え制御器４とが、メインコンピューター８１上で実現されている。また、高速ギヤモデル３と高精度ギヤモデル３ｂとは、それぞれ別個のサブコンピューター８２、８２ｂ上で実現されている。サブコンピューター８２では、高速ギヤモデル３がサブ通信Ｉ／Ｆ８４に接続されており、サブコンピューター８２ｂでは、高精度ギヤモデル３ｂがサブ通信Ｉ／Ｆ８４ｂに接続されている。メイン通信Ｉ／Ｆ８３、サブ通信Ｉ／Ｆ８４およびサブ通信Ｉ／Ｆ８４ｂは、不図示のネットワーク等を介して相互に接続されており、互いに通信を行って各パラメーターの受け渡しを実施できるように構成されている。

以上説明した本発明の第３の実施形態によれば、複数のコンピューターを用いて、第１の実施形態と同様の計算結果を得ることができる。また、複数のコンピューターに計算負荷を分散させることにより、少ないハードウェア資源でも高速に計算結果を得ることができる。さらに、各コンピューターの設置場所には制約がないため、物理的な専有面積や電力消費を分散させることも可能である。そのため、シミュレーション実行効率のさらなる向上を図ることができる。

なお、本実施形態では、第１の実施形態で説明したシミュレーターを複数のコンピューターに分散させた例を説明したが、第２の実施形態で説明したシミュレーターを複数のコンピューターに分散させることも可能である。また、本実施形態においても、第１の実施形態で説明した各変形例が可能である。

（第４の実施形態）
以上説明した第１〜第３の実施形態では、電気自動車の駆動装置をシミュレーション対象としたシステムシミュレーターの例を説明したが、本発明はこれに限定されず、様々な種類のシミュレーターに適用可能である。本実施形態では、電気モーターを用いた電動装置をシミュレーション対象としたシステムシミュレーターの例を説明する。

図１０は、本発明の第４の実施形態に関する従来のシミュレーターの構成例を示す図である。図１０に示す従来のシミュレーターでは、バッテリーモデル５と機構モデル６が高精度モーターモデル１ｂを介して接続されている。バッテリーモデル５、高精度モーターモデル１ｂおよび機構モデル６は、シミュレーション対象である電動装置を構成する各要素（バッテリー、モーター、駆動機構）をそれぞれモデル化したものである。なお、モーターが交流モーターである場合、高精度モーターモデル１ｂはインバーターを含んでモデル化することが好ましい。

図１０のシミュレーターでは、各モデル間の情報伝達用のパラメーターとして、電磁トルク２１、ローター回転速度２２、直流電圧２５、直流電流２６が用いられる。図１０の矢印の向きは、各パラメーターによる情報の伝達方向を示している。すなわち、図１０のシミュレーターは、バッテリーモデル５の発生する直流電圧２５を高精度モーターモデル１ｂを通して電磁トルク２１に変換し、機構モデル６に伝達する。また、機構モデル６の状態に応じたローター回転速度２２を高精度モーターモデル１ｂを通して直流電流２６に変換し、バッテリーモデル５に入力する。このときのバッテリーモデル５、高精度モーターモデル１ｂおよび機構モデル６の状態をそれぞれ計算することで、電動装置の動きをシミュレーションする。

図１１は、本発明の第４の実施形態に係るシミュレーターの構成例を示す図である。図１１に示すシミュレーターは、図１０に示した従来のシミュレーターと比較して、さらに高速モーターモデル１ｃが高精度モーターモデル１ｂと並列に配置されている。高速モーターモデル１ｃは、高精度モーターモデル１ｂと同様に、シミュレーション対象である電動装置を構成する要素の一つであるモーターをモデル化したものである。

図１１のシミュレーターは、第１の実施形態で説明した図１のシミュレーターと同様に、切り替え制御器４を有している。この切り替え制御器４の制御により、高精度モーターモデル１ｂを用いて計算を行うモードと、高速モーターモデル１ｃを用いて計算を行うモードとを、第１の実施形態と同様の手法で切り替えることができる。

以上説明した本発明の第４の実施形態によれば、電気モーターを用いた電動装置をシミュレーション対象としたシステムシミュレーターにおいても、第１の実施形態で説明したのと同様の効果を得ることができる。

なお、本実施形態では、第１の実施形態と同様の手法でモード切り替えを行う電動装置のシミュレーターの例を説明したが、第２の実施形態で説明したのと同様の手法でモード切り替えを行う電動装置のシミュレーターにも適用可能である。また、本実施形態においても、第１の実施形態で説明した各変形例が可能である。さらに、単独のコンピューター上で実現されるシミュレーターとしてもよいし、あるいは第３の実施形態で説明したように、複数のコンピューターに分散されたシミュレーターとしてもよい。

なお、本発明は上記した各実施形態や各種変形例に限定されるものではなく、それ以外にも様々な変形が可能である。たとえば、上記の各実施形態は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されない。また、ある実施形態の構成の一部を他の実施形態の構成に置き換えたり、ある実施形態の構成に他の実施形態の構成を加えたりすることも可能である。すなわち、各実施形態の構成の一部について、他の構成の追加・削除・置換をすることが可能である。

また、上記の各構成、機能、処理部、処理手段等は、それらの一部又は全部を、例えば集積回路のハードウェアで実現してもよい。あるいは、上記の各構成、機能等をソフトウェアで実現し、プロセッサがそれぞれの機能を実現するプログラムを解釈して実行することとしてもよい。また、上記の各実施形態において、制御線や情報線は説明上必要と考えられるものを示しており、製品上必ずしも全ての制御線や情報線を示しているとは限らない。実際には殆ど全ての構成が相互に接続されていると考えてもよい。

以上説明した各実施形態や各種変形例はあくまで一例であり、発明の特徴が損なわれない限り、本発明はこれらの内容に限定されるものではない。また、上記では種々の実施形態や変形例を説明したが、本発明はこれらの内容に限定されるものではない。本発明の技術的思想の範囲内で考えられるその他の態様も本発明の範囲内に含まれる。

１モーターモデル
１ｂ高精度モーターモデル
１ｃ高速モーターモデル
２車両モデル
３高速ギヤモデル
３ｂ高精度ギヤモデル
４切り替え制御器
５バッテリーモデル
６機構モデル
２１，２１ｂ電磁トルク
２２，２２ｂローター回転速度
２３，２３ｂ車軸トルク
２４，２４ｂ車軸回転速度
３１，３２，３３，３４，３１ｂ，３２ｂ，３３ｂ，３４ｂスイッチ

Claims

コンピューターを用いてシミュレーション対象に関する計算を行うシミュレーション方法であって、
前記シミュレーション対象は、第１要素および第２要素と、前記第１要素と前記第２要素の間に介在する中間要素と、を含み、
前記コンピューターを、前記第１要素をモデル化した第１モデルと、前記第２要素をモデル化した第２モデルと、同一の前記中間要素を異なる計算モデルでそれぞれモデル化した第１中間モデルおよび第２中間モデルと、としてそれぞれ機能させ、
前記第１モデルと前記第２モデルが前記第１中間モデルを介して情報を相互に送受信する第１モードと、前記第１モデルと前記第２モデルが前記第２中間モデルを介して情報を相互に送受信する第２モードと、のいずれかにより、前記コンピューターに前記計算を実行させ、
前記第１モードでは、前記第１モデルから出力される第１パラメーターを前記第１中間モデルに入力し、入力された前記第１パラメーターに応じて前記第１中間モデルから出力される第２パラメーターを前記第２モデルに入力すると共に、前記第２モデルから出力される第３パラメーターを前記第１中間モデルに入力し、入力された前記第３パラメーターに応じて前記第１中間モデルから出力される第４パラメーターを前記第１モデルに入力するように、前記コンピューターを動作させ、
前記第２モードでは、前記第１モデルから出力される第１パラメーターを前記第２中間モデルに入力し、入力された前記第１パラメーターに応じて前記第２中間モデルから出力される第２パラメーターを前記第２モデルに入力すると共に、前記第２モデルから出力される第３パラメーターを前記第２中間モデルに入力し、入力された前記第３パラメーターに応じて前記第２中間モデルから出力される第４パラメーターを前記第１モデルに入力するように、前記コンピューターを動作させ、
前記第１モードから前記第２モードに切り替える際に、
前記第１パラメーターまたは前記第３パラメーターを、前記第１中間モデルおよび前記第２中間モデルに入力し、
入力された前記第１パラメーターまたは前記第３パラメーターに応じて前記第２中間モデルから出力される前記第２パラメーターまたは前記第４パラメーターと、入力された前記第１パラメーターまたは前記第３パラメーターに応じて前記第１中間モデルから出力される前記第２パラメーターまたは前記第４パラメーターと、の比較を行い、
前記比較の結果が所定の基準を満たしたときに、前記第１モデルおよび前記第２モデルと前記第２中間モデルとを接続すると共に、前記第１モデルおよび前記第２モデルと前記第１中間モデルとを切り離すことで、前記第１モードから前記第２モードへの切り替えを実施するように、前記コンピューターを動作させるシミュレーション方法。
請求項１に記載のシミュレーション方法において、
前記比較では、前記第２中間モデルから出力される前記第２パラメーターまたは前記第４パラメーターと、前記第１中間モデルから出力される前記第２パラメーターまたは前記第４パラメーターとの差分を算出し、
算出した前記差分の絶対値が所定の基準値以下となったときに、前記第１モデルおよび前記第２モデルと前記第２中間モデルとを接続すると共に、前記第１モデルおよび前記第２モデルと前記第１中間モデルとを切り離すことで、前記第１モードから前記第２モードへの切り替えを実施するように、前記コンピューターを動作させるシミュレーション方法。
請求項２に記載のシミュレーション方法において、
前記第２中間モデルに入力される前記第１パラメーターおよび前記第３パラメーターの少なくとも一方の値を調整することで前記差分が小さくなるように、前記コンピューターを動作させるシミュレーション方法。
請求項３に記載のシミュレーション方法において、
前記第１モードから前記第２モードに切り替える際に、前記第３パラメーターを前記第１中間モデルおよび前記第２中間モデルに入力する場合には、前記第２中間モデルから出力される前記第４パラメーターと前記第１中間モデルから出力される前記第４パラメーターとの差分に応じて調整した前記第１パラメーターをさらに前記第２中間モデルに入力することで、前記差分が小さくなるように前記コンピューターを動作させ、
前記第１モードから前記第２モードに切り替える際に、前記第１パラメーターを前記第１中間モデルおよび前記第２中間モデルに入力する場合には、前記第２中間モデルから出力される前記第２パラメーターと前記第１中間モデルから出力される前記第２パラメーターとの差分に応じて調整した前記第３パラメーターをさらに前記第２中間モデルに入力することで、前記差分が小さくなるように前記コンピューターを動作させるシミュレーション方法。
請求項４に記載のシミュレーション方法において、
調整した前記第１パラメーターまたは前記第３パラメーターをさらに前記第２中間モデルに入力する際には、単位時間当たりの入力変化量を所定値以下に制限するシミュレーション方法。
請求項１から請求項５までのいずれか一項に記載のシミュレーション方法において、
前記第１モードから前記第２モードに切り替える際に前記第１中間モデルおよび前記第２中間モデルに入力される前記第１パラメーターまたは前記第３パラメーターは、アクロス変数であるシミュレーション方法。