(57)【特許請求の範囲】
【請求項１】
回路用シミュレータが、直流電流が重畳された三角波の電流がインダクタに流れる場合のシミュレーションを前記インダクタのシミュレーションモデルを用いて実行する方法であって、
前記インダクタのシミュレーションモデルを、
直流抵抗と、交流損失に合わせた見掛けの交流抵抗とインダクタンスとの直列回路であるインダクタの等価回路であって、前記インダクタの直流抵抗を、前記三角波の振幅と前記三角波の直流重畳電流とに対する第１の関数、インダクタの見掛けの交流抵抗を、前記三角波の振幅と前記三角波の直流重畳電流と前記三角波の周波数とに対する第２の関数、インダクタのインダクタンスを、前記三角波の振幅と前記三角波の直流重畳電流と前記三角波の周波数とに対する第３の関数とした前記等価回路と、
前記直流抵抗の流れる電流を表す、前記直流抵抗の両端の電圧値を、前記第１の関数で得られた値で除した数式により制御される第１のビヘイビア電流源と、
前記見掛けの交流抵抗を流れる電流を表す、前記交流抵抗の両端の電圧値を、前記第２の関数で得られた値で除した数式により制御される第２のビヘイビア電流源と、
前記インダクタンスを流れる電流を表す、前記インダクタンスの両端の電圧の積分値を、前記第３の関数で得られた値で除した数式により制御される第３のビヘイビア電流源と、
前記第１のビヘイビア電流源に流れる電流値が入力される平滑化回路と、
前記第１のビヘイビア電流源に流れる電流値から前記平滑化回路の出力値を減じた値が入力されるピークホールド回路と、
前記第１のビヘイビア電流源に流れる電流値から前記平滑化回路の出力値を減じた値が入力される波形変換回路と、
前記波形変換回路の出力値が入力されるカウンタ回路と、
で構成し、
前記回路シミュレータが、
前記第１のビヘイビア電流源に流れる電流値を、平滑化回路で平滑化して前記三角波の直流重畳電流を求め、
前記第１のビヘイビア電流源に流れる電流値から前記三角波の直流重畳電流の値を減じた値を、ピークホールド回路でピークホールドして前記三角波の振幅を求め、
前記第１のビヘイビア電流源に流れる電流値から前記三角波の直流重畳電流の値を減じた値を、波形変換回路でパルス波に変換し、前記パルス波をカウンタ回路でカウントして前記三角波の周波数を求め、
前記求めた三角波の直流重畳電流、三角波の振幅及び周波数を前記第１から第３の関数に代入して、前記インダクタの等価回路の直流抵抗、見かけの交流抵抗及びインダクタンスを再計算する、
ことを特徴とするインダクタのシミュレーション方法。
【請求項２】
回路用シミュレータが、直流電流が重畳された三角波の電流がインダクタに流れる場合のシミュレーションを、前記インダクタのシミュレーションモデルを用いて実行する方法であって、
前記インダクタのシミュレーションモデルを、
直流抵抗と、交流損失に合わせた見掛けの交流抵抗とインダクタンスとの直列回路であるインダクタの等価回路であって、前記インダクタの直流抵抗を、前記三角波の振幅と前記三角波の直流重畳電流とに対する第１の関数、インダクタの見掛けの交流抵抗を、前記三角波の振幅と前記三角波の直流重畳電流と前記三角波の周波数とに対する第２の関数、
インダクタのインダクタンスを、前記三角波の振幅と前記三角波の直流重畳電流と前記三角波の周波数とに対する第３の関数とした前記等価回路と、
前記直流抵抗を流れる電流を表す、前記直流抵抗の両端の電圧値を、前記第１の関数で得られた値で除した数式により制御される第１のビヘイビア電流源と、
前記見掛けの交流抵抗を流れる電流を表す、前記交流抵抗の両端の電圧値を、前記第２の関数で得られた値で除した数式により制御される第２のビヘイビア電流源と、
前記インダクタンスを流れる電流を表す、前記インダクタンスの両端の電圧の積分値を、前記第３の関数で得られた値で除した数式により制御される第３のビヘイビア電流源と、
前記第１のビヘイビア電流源に流れる電流値が入力される平滑化回路と、
前記第１のビヘイビア電流源に流れる電流値が入力されるピークホールド回路と、
前記第１のビヘイビア電流源に流れる電流値が入力される波形変換回路と、
前記波形変換回路の出力値が入力されるカウンタ回路と、
で構成し、
前記回路シミュレータが、
前記第１のビヘイビア電流源に流れる電流値を、平滑化回路で平滑化して前記三角波の直流重畳電流を求め、
前記第１のビヘイビア電流源に流れる電流値を、ピークホールド回路でピークホールドし、前記三角波の直流重畳電流の値を減じて前記三角波の振幅を求め、
前記第１のビヘイビア電流源に流れる電流値を、波形変換回路でパルス波に変換し、前記パルス波をカウンタ回路でカウントして前記三角波の周波数を求め、
前記求めた三角波の直流重畳電流、三角波の振幅及び周波数を前記第１から第３の関数に代入して、前記インダクタの等価回路の直流抵抗、見かけの交流抵抗及びインダクタンスを再計算する、
ことを特徴とするインダクタのシミュレーション方法。
【請求項３】
前記等価回路は、
前記見掛けの交流抵抗とインダクタンスとの直列回路に並列に、さらに、インダクタの浮遊容量を備える、
請求項１乃至請求項２のいずれかに記載のインダクタのシミュレーション方法。
【請求項４】
回路用シミュレータが、直流電流が重畳された三角波の電流がインダクタに流れる場合のシミュレーションを、前記インダクタのシミュレーションモデルを用いて実行する方法であって、
前記インダクタのシミュレーションモデルを、
直流抵抗と見掛けの交流抵抗との和からなる抵抗と、インダクタンスとの直列回路であるインダクタの等価回路であって、前記抵抗を、前記三角波の振幅と前記三角波の直流重畳電流と前記三角波の周波数とに対する第１の関数、インダクタのインダクタンスを、前記三角波の振幅と前記三角波の直流重畳電流と前記三角波の周波数とに対する第２の関数とした前記等価回路と、
前記抵抗を流れる電流を表す、前記抵抗の両端の電圧値を、前記第１の関数により得られた値で除した数式により制御される第１のビヘイビア電流源と、
前記インダクタンスを流れる電流を表す、前記インダクタンスの両端の電圧の積分値を、前記第２の関数により得られた値で除した数式により制御される第２のビヘイビア電流源と、
前記第１のビヘイビア電流源に流れる電流値が入力される平滑化回路と、
前記第１のビヘイビア電流源に流れる電流値から前記平滑化回路の出力値を減じた値が入力されるピークホールド回路と、
前記第１のビヘイビア電流源に流れる電流値から前記平滑化回路の出力値を減じた値が入力される波形変換回路と、
前記波形変換回路の出力値が入力されるカウンタ回路と、
で構成し、
前記回路シミュレータが、
前記第１のビヘイビア電流源に流れる電流値を、平滑化回路で平滑化して前記三角波の直流重畳電流を求め、
前記第１のビヘイビア電流源に流れる電流値から前記三角波の直流重畳電流の値を減じた値を、ピークホールド回路でピークホールドして前記三角波の振幅を求め、
前記第１のビヘイビア電流源に流れる電流値から前記三角波の直流重畳電流の値を減じた値を、波形変換回路でパルス波に変換し、前記パルス波をカウンタ回路でカウントして前記三角波の周波数を求め、
前記求めた三角波の直流重畳電流、三角波の振幅及び周波数を前記第１及び第２の関数に代入して、前記インダクタの等価回路の直流抵抗、見かけの交流抵抗及びインダクタンスを再計算する、
ことを特徴とするインダクタのシミュレーション方法。
【請求項５】
前記等価回路は、
前記抵抗と前記インダクタンスとの直列回路とに並列に、インダクタの浮遊容量を備える、
請求項４に記載のインダクタのシミュレーション方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ＳＰＩＣＥ（Simulation Program with Integrated Circuit Emphasis）等の回路用シミュレータで用いられる、インダクタのシミュレーションモデルに関する。
【先行技術文献】
【非特許文献】
【０００２】
【非特許文献１】Christophe Basso, SPICE Analog Behavioral Modeling of Variable Passives, Power Electronics Technology, May 2005, p57
【発明の概要】
【発明が解決しようとする課題】
【０００３】
一般的なインダクタのシミュレーションモデルは、インピーダンスアナライザ等の小信号測定機器を用いて回路定数を測定し、その回路定数をパラメータとする等価回路を用いている。
【０００４】
このような従来のインダクタのシミュレーションモデルを、微小な電流しか流れない小電力のシミュレーションに対して用いた場合は特に問題は生じない。しかし、例えばＤＣ−ＤＣコンバータ用のパワーインダクタのような、インダクタに大振幅の交流電流が流れるシミュレーションに用いた場合は問題が生じる。
【０００５】
図１は、基本的なＤＣ−ＤＣコンバータにおいて、パワーインダクタに流れる電流を説明する図である。図１において、横軸は時間ｔ、縦軸は電流ｉを示す。
インダクタには、振幅ＩＡＣ、周期Ｔ（周波数Ｆ）の三角波の交流電流ｉＡＣに、直流電流ＩＤＣが重畳された電流ｉＬ（＝ｉＡＣ＋ＩＤＣ）が流れる。ＤＣ−ＤＣコンバータは、出力電流と出力電圧、インダクタのＬ値などに応じて、周波数Ｆ、振幅ＩＡＣ、直流電流ＩＤＣを変更する。
一般的に、インダクタＬのインピーダンスＺは、数１式のように表現することができる。
【数１】
0006657822.tif 000002
【０００６】
数１式において、実数部のＲＳはインダクタＬの抵抗成分、虚数部のＬＳはインダクタＬの等価インダクタンス成分、ωは角周波数である。
数１式から、インダクタＬは、電力損失に関連する抵抗成分ＲＳと、インダクタンス成分ＬＳの直列回路で表現できることがわかる。
【０００７】
直流電流ＩＤＣ＝０の場合、つまり、電流ｉＬが交流電流ｉＡＣのみの場合、抵抗成分ＲＳは、インダクタの導体抵抗に相当する直流抵抗ＲＤＣと、磁性体の渦損失抵抗などに相当する交流抵抗ＲＡＣとの和で表すことができる。したがって、インダクタの等価回路は、図２（ａ）のように表すことができる。ここで、直流抵抗ＲＤＣは４端子法で実測して、予め求めることができ、交流抵抗ＲＡＣとインダクタンスＬＳは、専用の測定システムで、実測して、予め求めることができる。
【０００８】
ただし、図１に示す電流ｉＬがインダクタに流れた場合は、図２（ａ）の等価回路が適用できない。なぜなら、直流抵抗ＲＤＣの損失はすべての電流ｉＬに依存するのに対して、交流抵抗ＲＡＣの損失は、電流ｉＬのうち交流成分ｉＡＣのみに依存する。しかし、図２（ａ）の等価回路は、直流抵抗ＲＤＣにも交流抵抗ＲＡＣにも同じ電流ｉＬが流れることを前提条件としている。その結果、交流損失は正しく計算できない。直流電流ＩＤＣがあっても、交流損失を正しく計算できるようにするためには、図２（ｂ）に示すように、実際の交流抵抗ＲＡＣを見かけの交流抵抗Ｒ’ＡＣに置換えた等価回路とすれば良い。
【０００９】
なお、従来のインダクタのシミュレーションモデルで、直流抵抗ＲＤＣ、交流抵抗ＲＡＣ、インダクタンスＬＳとも値が固定の受動素子が用いられる。しかし、パワーインダクタなどの大振幅の電流が流れる回路では、インダクタの電流の直流成分ＩＤＣ、交流成分ｉＡＣおよび周波数Ｆが変動すると、直流抵抗ＲＤＣ、交流抵抗ＲＡＣ、インダクタンスＬＳとも変動する。これが、シミュレーション結果と、実際の回路の動作とが大きく乖離する原因となっている。
【課題を解決するための手段】
【００１０】
本発明は、
直流電流が重畳された三角波の電流がインダクタに流れる場合のインダクタのシミュレーションモデルであって、
インダクタの等価回路を、直流抵抗と、交流損失に合わせた見掛けの交流抵抗とインダクタンスとの直列回路とし、
インダクタの直流抵抗を、前記三角波の振幅と前記三角波の直流重畳電流とに対する第１の関数、
インダクタの見掛けの交流抵抗を、前記三角波の振幅と前記三角波の直流重畳電流と前記三角波の周波数とに対する第２の関数、
インダクタのインダクタンスを、前記三角波の振幅と前記三角波の直流重畳電流と前記三角波の周波数とに対する第３の関数、
とし、
前記直流抵抗は、前記直流抵抗の両端の電圧値を、前記第１の関数で得られた値で除した数式により制御される第１のビヘイビア電流源とし、
前記見掛けの交流抵抗は、前記交流抵抗の両端の電圧値を、前記第２の関数で得られた値で除した数式により制御される第２のビヘイビア電流源とし、
前記インダクタンスは、前記インダクタンスの両端の電圧の積分値を、前記第３の関数で得られた値で除した数式により制御される第３のビヘイビア電流源とし、
前記三角波の直流重畳電流は、前記第１のビヘイビア電流源に流れる電流値を、平滑化回路で平滑化することにより求められ、
前記三角波の振幅は、前記第１のビヘイビア電流源に流れる電流値から前記三角波の重畳電流の値を減じた値を、ピークホールド回路でピークホールドすることにより求められ、
前記三角波の周波数は、前記第１のビヘイビア電流源に流れる電流値から前記三角波の直流重畳電流の値を減じた値を、波形変換回路でパルス波に変換し、前記パルス波をカウンタ回路でカウントすることにより求められることを特徴とする。
【発明の効果】
【００１１】
本発明のインダクタのシミュレーションモデルを回路シミュレータに用いることで、パワーインダクタのような、大振幅の電流が流れるインダクタのシミュレーションでも、良好なシミュレーション結果を得ることができる。
【図面の簡単な説明】
【００１２】
【図１】インダクタに流れる電流を示すグラフである。
【図２】インダクタの等価回路モデルを示す図である。
【図３】本発明のインダクタのシミュレーションモデルの第１の実施例である。
【図４】本発明のインダクタのシミュレーションモデルの第２の実施例である。
【図５】本発明のインダクタのシミュレーションモデルの第３の実施例である。
【図６】本発明のインダクタのシミュレーションモデルの第４の実施例である。
【図７】シミュレーション結果と実際の結果を比較するグラフである。
【図８】シミュレーション結果と実際の結果を比較するグラフである。
【発明を実施するための形態】
【実施例】
【００１３】
以下の実施例では、インダクタのシミュレーションモデルを検討するにあたり、図２（ｂ）に示した等価回路を用いることとする。
【００１４】
（第１の実施例）
振幅ＩＡＣ、周波数Ｆの三角波状の交流電流ｉＡＣに、直流電流ＩＤＣが重畳された電流ｉＬ（＝ｉＡＣ＋ＩＤＣ）がインダクタに流れた場合の、交流抵抗ＲＡＣおよびインダクタンスＬＳを、専用の測定システムを用いて測定する。この測定したデータを数式処理システム等を用いて、数２式、数３式、数４式を求める。
【数２】
0006657822.tif 000003
【数３】
0006657822.tif 000004
【数４】
0006657822.tif 000005
【００１５】
数２式は、直流抵抗ＲＤＣが、振幅ＩＡＣ、直流電流ＩＤＣの近似関数であることを示している。
数３式は、交流抵抗ＲＡＣが、周波数Ｆ、振幅ＩＡＣ、直流電流ＩＤＣの近似関数であることを示している。
数４式は、インダクタンスＬＳが、周波数Ｆ、振幅ＩＡＣ、直流電流ＩＤＣの近似関数であることを示している。
【００１６】
ここで、交流抵抗ＲＡＣと見掛けの交流抵抗Ｒ’ＡＣとの関係は次の通りである。
インダクタＬの損失ＰＬＯＳＳは、図２（ｂ）に示した抵抗成分Ｒｓの損失と同じなので、電流ｉＬの実効値をｉＬｒｍｓとすると、数５式で表される。
【数５】
0006657822.tif 000006
【００１７】
一方、インダクタＬの損失ＰＬＯＳＳは、直流抵抗ＲＤＣによる直流損失と交流抵抗ＲＡＣによる交流損失との和であるので、交流電流ｉＡＣの実効値ｉＡＣｒｍｓとすると、式６式で表される。
【数６】
0006657822.tif 000007
【００１８】
数５式と数６式とから、数７式が得られる。
【数７】
0006657822.tif 000008
【００１９】
数７式から、見掛けの交流抵抗Ｒ’ＡＣは、交流抵抗ＲＡＣを補正係数（ｉＡＣｒｍｓ／ｉＬｒｍｓ）２で補正した値であることがわかる。
【００２０】
ここで、周期関数が、φ＝α＋β（ｔ）で表されるとき、周期関数φの実効値φｒｍｓは数８式のようになる。ここで、αは直流成分、β（ｔ）は直流成分を含まない交流成分とする。
【数８】
0006657822.tif 000009
【００２１】
交流成分β（ｔ）の平均値は０なので、周期関数の平均値の定義式から、数９式が成り立つ。
【数９】
0006657822.tif 000010
【００２２】
また、β（ｔ）の実効値をβｒｍｓとすると、数１０式が成り立つ。
【数１０】
0006657822.tif 000011
【００２３】
数９式および数１０式を、数８式に代入して、数１１式が得られる。数１１は、任意の周期関数に適用可能な一般的な関係式である。
【数１１】
0006657822.tif 000012
【００２４】
数１１式から、電流ｉＬの実効値をｉＬｒｍｓ、交流電流ｉＡＣの実効値をｉＡＣｒｍｓとすると、電流ｉＬの実効値ｉＬｒｍｓ、交流電流ｉＡＣの実効値ｉＡＣｒｍｓ、直流電流ＩＤＣの間には、数１２式の関係が成り立つ。
【数１２】
0006657822.tif 000013
【００２５】
数１２式を数７式に代入して、数１３式が得られる。
【数１３】
0006657822.tif 000014
【００２６】
交流電流ｉＡＣは、三角波なので、数１４式が成り立つ。
【数１４】
0006657822.tif 000015
【００２７】
数１４式を数１３式に代入して、数１５式が得られる。
【数１５】
0006657822.tif 000016
【００２８】
数１５式から、三角波の場合には、交流抵抗ＲＡＣを補正係数１／（１＋３＊（ＩＤＣ／ＩＡＣ）２）で補正した値が、見掛けの交流抵抗Ｒ’ＡＣであると分かる。
【００２９】
数３式に示したように、交流抵抗ＲＡＣは、周波数Ｆ、振幅ＩＡＣ、直流電流ＩＤＣの関数であると定義したものなので、数１５式から、見かけの交流抵抗Ｒ’ＡＣもまた、周波数Ｆ、振幅ＩＡＣ、直流電流ＩＤＣの関数であることが分かる。すなわち、数１６式が成り立つ。
【数１６】
0006657822.tif 000017
【００３０】
したがって、インダクタに流れる電流ｉＬの、周波数Ｆ、振幅ＩＡＣ、直流電流ＩＤＣが分かれば、数２式、数４式、数１６式から、直流抵抗ＲＤＣ、見かけの交流抵抗Ｒ’ＡＣ、インダクタンスＬＳを計算することができる。
【００３１】
なお、数１５式において、ＩＤＣ＝０の時は、Ｒ’ＡＣ＝ＲＡＣとなるので、ＲＳ＝ＲＤＣ＋ＲＡＣが成り立つ。つまり、図２（ａ）の等価回路は、図２（ｂ）の等価回路において、ＩＤＣ＝０の場合の例を示している。
【００３２】
ＤＣ−ＤＣコンバータにおいて、ＩＤＣは負荷に流れた電流（負荷電流）と等しい。一般的に、ＤＣ−ＤＣコンバータの負荷電流が非常に小さいときは、Ｒ’ＡＣ≒ＲＡＣとなっている。しかし、ＩＤＣ＞＞ＩＡＣのときは、Ｒ’ＡＣ≒０となるので、ＲＳ≒ＲＤＣとなる。つまり、ＩＤＣ＞＞ＩＡＣの時、インダクタＬにおける損失の大部分は、直流抵抗ＲＤＣによるものである。なお、ＤＣ−ＤＣコンバータの負荷電流が大きいときは、ＲＳ≒ＲＤＣとなっている。
【００３３】
上記した、直流抵抗ＲＤＣ、見かけの交流抵抗Ｒ’ＡＣ、インダクタンスＬＳの値を用いてシミュレーションを行う場合は、周波数Ｆ、振幅ＩＡＣ、直流電流ＩＤＣの値から、逐次、直流抵抗ＲＤＣ、見かけの交流抵抗Ｒ’ＡＣとインダクタンスＬＳの値を計算し直す必要がある。しかし、ＳＰＩＣＥのような回路シミュレータは、抵抗素子やインダクタンス素子などの受動素子のパラメータを、解析中に変更できない。そこで、解析中でも再計算可能な数式で、多様な素子を定義できるビヘイビア電源を使って、受動素子を擬似的に表現する。
【００３４】
ここで、抵抗とインダクタをビヘイビア電流源で表現する方法について説明する。
【００３５】
抵抗Ｒに流れる電流をＩＲ、抵抗の両端の電圧をＶＲとすると、電流ＩＲ、電圧ＶＲ、抵抗値Ｒの関係は数１７式のように表される。
【数１７】
0006657822.tif 000018
【００３６】
したがって、抵抗Ｒは、ＶＲ／Ｒを用いたビヘイビア電流源で表現することができる。
【００３７】
インダクタＬに流れる電流ＩＬ、インダクタＬの両端の電圧ＶＬとすると、電流ＩＬと電圧ＶＬの関係は数１８式のように表される。
【数１８】
0006657822.tif 000019
【００３８】
数１２式を積分すると、数１９式が得られる。
【数１９】
0006657822.tif 000020
【００３９】
ここで、ＩＮＴは積分を意味する。
したがって、インダクタは、ＩＮＴ（ＶＬ）/Ｌを用いたビヘイビア電流源で表現することができる。
【００４０】
なお、数２式、数４式、数１６式のための、周波数Ｆ、振幅ＩＡＣ、直流電流ＩＤＣを求める必要がある。周波数Ｆ、振幅ＩＡＣ、直流電流ＩＤＣは、それぞれインダクタの電流ｉＬから求めることができる。
直流電流ＩＤＣは、電流ｉＬを平滑化することにより求めることができる。
振幅ＩＡＣは、電流ｉＬから直流電流ＩＤＣを減じた値の最大値により求めることができる。
周波数Ｆは、電流ｉＬから直流電流ＩＤＣを減じた値と、その値を平滑化した値とを比較して、結果から得られたパルス信号をカウンターで計測することにより求めることができる。
ここで、回路処理の都合により、電流信号は、ビヘイビア電圧源によって、等値の電圧信号に変換されて、周波数Ｆ、振幅ＩＡＣ、直流電流ＩＤＣの値が得られる。
【００４１】
上記の結果から、ビヘイビア電源を用いたインダクタのシミュレーションモデルは、図３のように表すことができる。
図３は、本発明のインダクタのシミュレーションモデルの第１の実施例を示す。
【００４２】
図３に示すように、インダクタＬのシミュレーションモデルは、
擬似的な受動素子を表現するためのビヘイビア電流源Ｂ＿ＲＤＣ、Ｂ＿Ｒ’ＡＣ、Ｂ＿ＬＳ、
電流信号を電圧信号に変換するためのビヘイビア電圧源Ｂ＿ｉＬ、Ｂ＿ｉＡＣ、
平滑回路ＡＶＥ、ピークホールド回路ＰＥＡＫ、波形変換回路ＣＯＮＶ、カウンタ回路ＣＯＵＮＴからなる。
【００４３】
ビヘイビア電流源Ｂ＿ＲＤＣは、図２（ｂ）に示した等価回路において、直流抵抗ＲＤＣに相当し、
ビヘイビア電流源Ｂ＿Ｒ’ＡＣは、図２（ｂ）に示した等価回路において、交流抵抗Ｒ’ＡＣに相当し、
ビヘイビア電流源Ｂ＿ＬＳは、図２（ｂ）に示した等価回路において、インダクタンスＬＳに相当する。
【００４４】
ビヘイビア電流源Ｂ＿ＲＤＣ、Ｂ＿Ｒ’ＡＣ、Ｂ＿ＬＳは、ノード１、２に直列に接続される。ここで、ビヘイビア電流源Ｂ＿ＲＤＣ、Ｂ＿Ｒ’ＡＣの間をノードｍ、ビヘイビア電流源Ｂ＿Ｒ’ＡＣ、Ｂ＿ＬＳの間をノードｎとする。
【００４５】
ビヘイビア電流源Ｂ＿ＲＤＣの電流Ｉ＿ＲＤＣは、数２０式により表現される。ここで、Ｖ（１，ｍ）は、ノード１とノードｍとの間の電圧、すなわち、ビヘイビア電流源Ｂ＿ＲＤＣの両端の電圧を示す。
【数２０】
0006657822.tif 000021
【００４６】
ビヘイビア電流源Ｂ＿Ｒ’ＡＣの電流Ｉ＿Ｒ’ＡＣは、数２１式により表現される。ここで、Ｖ（ｍ，ｎ）は、ノードｍとノードｎとの間の電圧、すなわち、ビヘイビア電流源Ｂ＿Ｒ’ＡＣの両端の電圧を示す。
【数２１】
0006657822.tif 000022
【００４７】
ビヘイビア電流源Ｂ＿ＬＳの電流Ｉ＿ＬＳは、数２２式により表現される。ここで、電圧Ｖ（ｎ，２）は、ノードｎとノード２との間の電圧、すなわち、ビヘイビア電流源Ｂ＿ＬＳの両端の電圧を示し、ＩＮＴは、積分を示す。
【数２２】
0006657822.tif 000023
【００４８】
次に、図３のインダクタのシミュレーションモデルの動作について説明する。
【００４９】
（ステップ１）
ビヘイビア電圧源Ｂ＿ｉＬのパラメータＶ１に、電流ｉＬ（つまり、ビヘイビア電流源Ｂ＿ＲＤＣに流れる電流ｉ（Ｂ＿ＲＤＣ））が与えられる。
【００５０】
ビヘイビア電圧源Ｂ＿ｉＬの出力は、ｉＬと等値の電圧信号が得られる。
【００５１】
（ステップ２）
ビヘイビア電圧源Ｂ＿ｉＬの出力電圧を、平滑回路ＡＶＥで平滑化することにより、直流電流ＩＤＣと等値の電圧値が得られる。
【００５２】
（ステップ３）
ビヘイビア電圧源Ｂ＿ｉＡＣのパラメータＶ２に、電流ｉＬと直流電流ＩＤＣとの差（つまり、（ｉＬ−ＩＤＣ＝ｉ（Ｂ＿ＲＤＣ）−ＩＤＣ））が与えられる。
これにより、ビヘイビア電圧源Ｂ＿ｉＡＣの出力は、ｉＬ−ＩＤＣ、つまり交流電流ｉＡＣと等値の電圧信号となる。
【００５３】
（ステップ４）
ビヘイビア電圧源Ｂ＿ｉＡＣの出力電圧を、ピークホールド回路ＰＥＡＫでピークホールドすることにより、振幅ＩＡＣと等値の電圧値が得られる。
【００５４】
（ステップ５）
ビヘイビア電圧源Ｂ＿ｉＡＣの出力電圧を、波形変換回路ＣＯＮＶで三角波から矩形波に変換し、さらに、その矩形波をカウンタ回路ＣＯＵＮＴで計測することにより、周波数Ｆが得られる。
【００５５】
（ステップ６）
ステップ１〜５で得られた周波数Ｆ、直流電流ＩＤＣ、振幅ＩＡＣを、数２式、数４式、数１６式に代入して、直流抵抗ＲＤＣ、見かけの交流抵抗Ｒ’ＡＣ、インダクタンスＬＳの値を再び計算することを、繰り返してシミュレーション結果を得ることができる。
【００５６】
なお、平滑回路ＡＶＥは、例えば、ＲＣ平滑回路などで構成することができる。
ピークホールド回路ＰＥＡＫは、例えば、コンデンサとダイオードからなるピークホールド回路、または、ＲＣ積分回路の出力を倍増する回路などで構成することができる。
波形変換回路ＣＯＮＶは、ビヘイビア電圧源Ｂ_ｉＡＣの出力電圧Ｖ２を用いた入力信号と、この入力信号の平均化した信号を比較するコンパレータなどで構成することができる。
波形変換回路ＣＯＮＶの入力信号として、ビヘイビア電圧源Ｂ_ｉＬの出力電圧Ｖ１を使っても同様な結果が得られる。
【００５７】
図３に示したインダクタのシミュレーションモデルは、一つのインダクタが複数の素子および複数の回路で表現されているため、取り扱いが面倒である。この複数の素子および複数の回路を、サブサーキットとしてモジュール化して、このサブサーキットにインダクタの回路図形を割り当てれば、一般的な受動部品と同様に、容易に取り扱うことができる。
【００５８】
なお、上記したインダクタのシミュレーションモデルにおいて、ビヘイビア電流源Ｂ＿ＲＤＣ、Ｂ＿Ｒ’ＡＣ、Ｂ＿ＬＳを主回路と呼び、ビヘイビア電圧源Ｂ＿ｉＬ、Ｂ＿ｉＡＣ、平滑回路ＡＶＥ、ピークホールド回路ＰＥＡＫ、波形変換回路ＣＯＮＶ、カウンタ回路ＣＯＵＮＴを補助回路と呼ぶ。この、主回路と補助回路との間は、回路で接続されているわけでなく、シミュレーション時に値のやり取りが同期して行われるだけなので、補助回路は主回路に対して電気的な影響を与えない。
【００５９】
ところで、インダクタに大電流が流れると、温度が上昇し、そのため直流抵抗ＲＤＣが多少変化する。
直流抵抗ＲＤＣを、導体の抵抗率の温度係数と温度上昇とで補正できれば、より正確な結果を得ることができるが、ＳＰＩＣＥのような回路シミュレータは、温度上昇を加味したシミュレーションすることができない。
このような場合は、温度上昇が電流ｉＬの実効値ｉＬｒｍｓに依存することを利用すれば、温度上昇を加味したシミュレーションすることができる。
【００６０】
数１２式によれば、電流ｉＬの実効値ｉＬｒｍｓは、直流電流ＩＤＣと交流電流ｉＡＣの実効値ｉＡＣｒｍｓとから求めることができる。つまり、直流抵抗ＲＤＣも、直流電流ＩＤＣと交流電流ｉＡＣの関数である数２式で表されるということが分かる。
よって、インダクタの温度上昇の、直流電流ＩＤＣと交流電流ｉＡＣの依存性を、予め実測して求めておけば、直流抵抗ＲＤＣの温度補正を加味したインダクタのシミュレーションモデルとすることができる。
【００６１】
（第２の実施例）
一般的に、ＤＣ−ＤＣコンバータ用のパワーインダクタが使用される周波数は、そのパワーインダクタの自己共振周波数より遥かに低い。例えば、ｕＨオーダーのチップ型のパワーインダクタの場合、自己共振周波数は数１０ＭＨｚ〜１００ＭＨｚであるが、ＤＣ−ＤＣコンバータで常用される周波数は数ＭＨｚである。したがってＤＣ−ＤＣコンバータ用のパワーインダクタのシミュレーションでは、共振を考慮しなくても問題ない。しかし、インダクタのシミュレーションモデルをより忠実に表現するため、シミュレーションモデルに浮遊容量を追加してもよい。
【００６２】
図４は、本発明のインダクタのシミュレーションモデルの第２の実施例を示す。
図４に示すように、第１の実施例のノードｍとノード２との間に、浮遊容量ＣＳＴＲＡＹを並列に接続している。
第２の実施例の図４において、第１の実施例の図３で用いたのと同一符号は同一の構成要素を指すものとし、詳しい説明は省略する。
【００６３】
常用される周波数範囲においては、浮遊容量ＣＳＴＲＡＹのインピーダンスは非常に高いため、インダクタの電流ｉＬから、浮遊容量ＣＳＴＲＡＹに分流される電流は極僅かである。したがって、電流ｉＬは、すべてインダクタＬＳに流れると見做して良い。なお、浮遊容量ＣＳＴＲＡＹは、自己共振周波数を実測して求めることができる。浮遊容量ＣＳＴＲＡＹは、ノード１とノード２との間に接続してもよい。
【００６４】
（第３の実施例）
図５は、本発明のインダクタのシミュレーションモデルの第３の実施例を示す。
インダクタの等価回路は、第１の実施と同じだが、補助回路にビヘイビア電圧源Ｂ＿ｉＡＣのみを用いて、周波数Ｆ、直流電流ＩＤＣ、振幅ＩＡＣを求めて点が異なる。
第３の実施例の図５において、第１の実施例１の図３で用いたのと同一符号は同一の構成要素を指すものとし、詳しい説明は省略する。
【００６５】
第１の実施例と同様に、ビヘイビア電圧源Ｂ＿ｉＬの出力電圧を、平滑回路ＡＶＥで平滑化することにより、直流電流ＩＤＣを得られる。
【００６６】
第３の実施例では、
電流ｉＬの最大値ｉＬｍａｘ（＝ＩＡＣ＋ＩＤＣ）は、ビヘイビア電圧源Ｂ＿ｉＬの出力電圧を、ピークホールド回路ＰＡＥＫでピークホールドして求めている。
また、振幅ＩＡＣは、ＩＡＣ＝ｉＬｍａｘ−ＩＤＣなので、最大値ｉＬｍａｘから直流電流ＩＤＣを減じて求めている。
また、周波数Ｆは、ビヘイビア電圧源Ｂ＿ｉＬの出力電圧Ｖ１を、波形変換回路ＣＯＮＶで三角波から矩形波に変換し、さらに、その矩形波をカウンタ回路ＣＯＵＮＴで計測することにより求めている。
【００６７】
（第４の実施例）
図６は、本発明のインダクタのシミュレーションモデルの第４の実施例である。
直流抵抗ＲＤＣと見掛けの交流抵抗Ｒ’ＡＣとで構成した要素を、まとめて一つの抵抗成分ＲＳ（＝ＲＤＣ＋Ｒ’ＡＣ）とし、ビヘイビア電流源Ｂ＿ＲＳで表現している点が第１の実施例と異なる。
【００６８】
図６において、ビヘイビア電流源Ｂ＿ＲＳに流れた電流ｉ（Ｂ＿ＲＳ）はインダクタの電流ｉＬなので、
ビヘイビア電圧源Ｂ＿ｉＬとビヘイビア電圧源Ｂ＿ｉＡＣのパラメータ中のｉ（Ｂ＿ＲＤＣ）が、ｉ（Ｂ＿ＲＳ）に置き換えられている。
第４の実施例の図６において、第１の実施例１の図３で用いたのと同一符号は同一の構成要素を指すものとし、詳しい説明は省略する。
【００６９】
ここで、抵抗成分ＲＳは数２３式で表される。
【数２３】
0006657822.tif 000024
【００７０】
ビヘイビア電流源Ｂ＿ＲＳの電流Ｉ＿ＲＳは、数２４式により表現される。
【数２４】
0006657822.tif 000025
【００７１】
図７および図８は、直流電流ＩＤＣを変化させたときの、交流抵抗ＲＡＣとインダクタンスＬＳを、シミュレーションモデルから計算したシミュレーション結果と、実際の結果を比較するグラフであり、図７では、Ｆ＝３ＭＨｚ、ＩＡＣ＝０．９Ａであり、図８では、Ｆ＝５ＭＨｚ、ＩＡＣ＝０．７Ａである。
図７と図８において、
◇は、実測したＬＳ、
△は、インダクタのシミュレーションモデルから得られた周波数Ｆ、直流電流ＩＤＣ、振幅ＩＡＣを用い、近似関数の数４式で計算したＬＳ、
□は、実測したＲａｃ、
×は、インダクタのシミュレーションモデルから得られた周波数Ｆ、直流電流ＩＤＣ、振幅ＩＡＣを用い、近似関数の数３式で計算したＲＡＣ、
を示している。
【００７２】
図７と図８の結果から、周波数Ｆ、直流電流ＩＤＣ、振幅ＩＡＣが変わっても、シミュレーション結果と実際の結果が良く一致していることが分かる。
したがって、本願によるインダクタのシミュレーションモデルは、大振幅の電流が流れるインダクタのシミュレーションであっても、精度の高い、良好なシミュレーション結果を得ることができることが分かる。
【符号の説明】
【００７３】
Ｂ＿ＲＤＣ、Ｂ＿Ｒ’ＡＣ、Ｂ＿ＬＳ  ビヘイビア電流源
Ｂ＿ｉＬ、Ｂ＿ｉＡＣ  ビヘイビア電圧源
ＡＶＥ  平滑回路
ＰＡＥＫ  ピークホールド回路
ＣＯＮＶ  波形変換回路
ＣＯＵＮＴ  カウンタ回路
【図１】
0006657822.tif 000026
【図２】
0006657822.tif 000027
【図３】
0006657822.tif 000028
【図４】
0006657822.tif 000029
【図５】
0006657822.tif 000030
【図６】
0006657822.tif 000031
【図７】
0006657822.tif 000032
【図８】
0006657822.tif 000033