JPH11153634A - シミュレーション装置及びシミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体 - Google Patents
シミュレーション装置及びシミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体Info
- Publication number
- JPH11153634A JPH11153634A JP31789197A JP31789197A JPH11153634A JP H11153634 A JPH11153634 A JP H11153634A JP 31789197 A JP31789197 A JP 31789197A JP 31789197 A JP31789197 A JP 31789197A JP H11153634 A JPH11153634 A JP H11153634A
- Authority
- JP
- Japan
- Prior art keywords
- time
- simulation
- circuit
- value
- electric field
- 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
- 238000004088 simulation Methods 0.000 title claims abstract description 151
- 230000005684 electric field Effects 0.000 claims abstract description 97
- 238000004364 calculation method Methods 0.000 claims abstract description 31
- 230000001052 transient effect Effects 0.000 claims abstract description 12
- 238000000034 method Methods 0.000 claims description 37
- 238000013459 approach Methods 0.000 claims description 16
- 230000002250 progressing effect Effects 0.000 claims description 2
- 239000004020 conductor Substances 0.000 description 11
- 238000010586 diagram Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 6
- 230000006399 behavior Effects 0.000 description 4
- 230000004927 fusion Effects 0.000 description 3
- 230000005672 electromagnetic field Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 238000002940 Newton-Raphson method Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measuring Magnetic Variables (AREA)
- Tests Of Electronic Circuits (AREA)
Abstract
(57)【要約】
【課題】 電磁波解析と回路解析との融合シミュレーシ
ョンを少ない計算量で行い、かつ安定した解が得られる
ようにする。 【解決手段】 電磁波解析手段１は、時刻ｔem01で磁界
を計算し、時刻ｔem02で電界を計算することにより、過
渡的な電磁波解析を行っている。回路解析手段２が回路
方程式を解くべき時刻ｔcsは、状況に応じて値が変化す
る時間刻み幅Δｔ csで進行させられる。電流源値引き渡
し手段３は、時刻ｔem01と時刻ｔcsとが、時間差λ1 内
に接近した場合には、電磁波解析手段１により計算され
た磁界に基づいて回路の電流源値を算出し、算出した電
流源値を回路解析手段２へ引き渡す。電界値引き渡し手
段４は、時刻ｔem02と時刻ｔcsとが、時間差λ2 内に接
近した場合には、回路解析手段２により求められる回路
への電圧に基づいて、回路が存在している領域の電界値
を算出し、算出した電界値を電磁波解析手段１へ引き渡
す。
ョンを少ない計算量で行い、かつ安定した解が得られる
ようにする。 【解決手段】 電磁波解析手段１は、時刻ｔem01で磁界
を計算し、時刻ｔem02で電界を計算することにより、過
渡的な電磁波解析を行っている。回路解析手段２が回路
方程式を解くべき時刻ｔcsは、状況に応じて値が変化す
る時間刻み幅Δｔ csで進行させられる。電流源値引き渡
し手段３は、時刻ｔem01と時刻ｔcsとが、時間差λ1 内
に接近した場合には、電磁波解析手段１により計算され
た磁界に基づいて回路の電流源値を算出し、算出した電
流源値を回路解析手段２へ引き渡す。電界値引き渡し手
段４は、時刻ｔem02と時刻ｔcsとが、時間差λ2 内に接
近した場合には、回路解析手段２により求められる回路
への電圧に基づいて、回路が存在している領域の電界値
を算出し、算出した電界値を電磁波解析手段１へ引き渡
す。
Description
【０００１】
【発明の属する技術分野】本発明は電磁波解析と回路解
析との融合シミュレーションを行うシミュレーション装
置及びシミュレーションプログラムを記録したコンピュ
ータ読み取り可能な記録媒体に関し、特に回路解析にお
ける時間刻み幅が状況に応じて随時変化するような融合
シミュレーションを行うシミュレーション装置及びその
ようなシミュレーションのためのシミュレーションプロ
グラムを記録したコンピュータ読み取り可能な記録媒体
に関する。
析との融合シミュレーションを行うシミュレーション装
置及びシミュレーションプログラムを記録したコンピュ
ータ読み取り可能な記録媒体に関し、特に回路解析にお
ける時間刻み幅が状況に応じて随時変化するような融合
シミュレーションを行うシミュレーション装置及びその
ようなシミュレーションのためのシミュレーションプロ
グラムを記録したコンピュータ読み取り可能な記録媒体
に関する。
【０００２】
【従来の技術】電磁波の過渡的な挙動を電子計算機を使
った数値シミュレーションによって解析する方法の１つ
に、有限差分時間領域( ＦＤＴＤ：Finite Difference
Time Domain ）法がある。この方法は、マックスウェル
方程式を時間と空間において差分法で解く手法であり、
その適用範囲の広さ等から広く用いられている手法であ
る。
った数値シミュレーションによって解析する方法の１つ
に、有限差分時間領域( ＦＤＴＤ：Finite Difference
Time Domain ）法がある。この方法は、マックスウェル
方程式を時間と空間において差分法で解く手法であり、
その適用範囲の広さ等から広く用いられている手法であ
る。
【０００３】一方、回路の過渡的な挙動を電子計算機を
使った数値シミュレーションによって解析する方法があ
る。この方法は、一般的に、回路の節点電圧、節点電流
を時々刻々キルヒホッフの法則に従って求めるものであ
る。さらに、半導体素子などの非線形な電流電圧特性を
持つ素子などについては、その特性方程式を表現するプ
ログラムを装備し、非線形方程式を解く際にはニュート
ンラプソン法などを使用する。
使った数値シミュレーションによって解析する方法があ
る。この方法は、一般的に、回路の節点電圧、節点電流
を時々刻々キルヒホッフの法則に従って求めるものであ
る。さらに、半導体素子などの非線形な電流電圧特性を
持つ素子などについては、その特性方程式を表現するプ
ログラムを装備し、非線形方程式を解く際にはニュート
ンラプソン法などを使用する。
【０００４】また、電磁波解析と回路解析とを融合した
数値シミュレーション方法が提案されている。このシミ
ュレーション方法では、電磁波解析で定義される電界や
磁界と回路解析で定義される電圧や電流を関連付けなが
ら実行する。電磁波解析と回路解析とを融合した数値シ
ミュレーションは、回路素子の特性とその周囲の電磁界
現象を統一的に解析できるといった特徴を持っており、
回路中を伝搬する高周波信号の解析に非常に有用であ
る。
数値シミュレーション方法が提案されている。このシミ
ュレーション方法では、電磁波解析で定義される電界や
磁界と回路解析で定義される電圧や電流を関連付けなが
ら実行する。電磁波解析と回路解析とを融合した数値シ
ミュレーションは、回路素子の特性とその周囲の電磁界
現象を統一的に解析できるといった特徴を持っており、
回路中を伝搬する高周波信号の解析に非常に有用であ
る。
【０００５】ところで、電磁波の過渡的な挙動を有限時
間差分領域法を用いて数値シミュレーションする場合、
時間刻み幅Δｔemは次式
間差分領域法を用いて数値シミュレーションする場合、
時間刻み幅Δｔemは次式
【０００６】
【数１】
【０００７】の安定条件を満たす範囲でできるだけ大き
く設定することが望ましい。ここで、ｖは電磁波の速度
であり、Δｘmin 、Δｙmin 、Δｚmin 、はそれぞれＸ
ＹＺ方向の空間離散間隔の最小値である。このような時
間刻み幅を設定するのは、Δｔ emが式（１）の右辺の値
よりも大きい場合は数値シミュレーションが発散してし
まい、Δｔemが必要以上に小さいと、計算時間が余計に
かかるからである。そこで通常、次式
く設定することが望ましい。ここで、ｖは電磁波の速度
であり、Δｘmin 、Δｙmin 、Δｚmin 、はそれぞれＸ
ＹＺ方向の空間離散間隔の最小値である。このような時
間刻み幅を設定するのは、Δｔ emが式（１）の右辺の値
よりも大きい場合は数値シミュレーションが発散してし
まい、Δｔemが必要以上に小さいと、計算時間が余計に
かかるからである。そこで通常、次式
【０００８】
【数２】
【０００９】を満たすようなΔｔemが設定される。この
Δｔemはシミュレーションの開始から終了まで一定なの
が一般的である。ところが、回路の過渡的な挙動を数値
シミュレーションする場合、時間刻み幅Δｔcsは、状況
に応じて増減するのが一般的である。例えば、非線形方
程式を解く際に電位差の値が収束しない場合には、時間
刻み幅を小さくする必要がある。そこで、電磁波解析と
回路解析との間でのデータの受け渡しを行う場合には、
両者の時刻を確実かつ効率的に一致させ、その時刻でデ
ータの受け渡しをする。
Δｔemはシミュレーションの開始から終了まで一定なの
が一般的である。ところが、回路の過渡的な挙動を数値
シミュレーションする場合、時間刻み幅Δｔcsは、状況
に応じて増減するのが一般的である。例えば、非線形方
程式を解く際に電位差の値が収束しない場合には、時間
刻み幅を小さくする必要がある。そこで、電磁波解析と
回路解析との間でのデータの受け渡しを行う場合には、
両者の時刻を確実かつ効率的に一致させ、その時刻でデ
ータの受け渡しをする。
【００１０】以下に、電磁波解析と回路解析との間での
データの受け渡しのタイミングを説明する。図９は、従
来のシミュレーション装置によるデータの引き渡しタイ
ミングを示す図である。図中、上段に電磁波解析処理を
示しており、下段に回路解析処理を示している。なお、
以下に示す数式における各記号の右肩の数値は、ｋ×Δ
ｔemで表される時刻の「ｋ」の値を示したものである。 ［Ｓ３１］電磁波解析における時刻ｔem＝（ｎ＋１／
２）Δｔemにおいて、磁界を計算する。ここで、ｎは自
然数である。磁界Ｈの計算は、次の式によって行う。
データの受け渡しのタイミングを説明する。図９は、従
来のシミュレーション装置によるデータの引き渡しタイ
ミングを示す図である。図中、上段に電磁波解析処理を
示しており、下段に回路解析処理を示している。なお、
以下に示す数式における各記号の右肩の数値は、ｋ×Δ
ｔemで表される時刻の「ｋ」の値を示したものである。 ［Ｓ３１］電磁波解析における時刻ｔem＝（ｎ＋１／
２）Δｔemにおいて、磁界を計算する。ここで、ｎは自
然数である。磁界Ｈの計算は、次の式によって行う。
【００１１】
【数３】 ここで、μは透磁率である。 ［Ｓ３２］電磁波解析によって求め出された磁界Ｈに基
づいて電流源値Ｉを算出し、回路解析へ引き渡す。電流
源値Ｉの算出は、次の式によって行われる。
づいて電流源値Ｉを算出し、回路解析へ引き渡す。電流
源値Ｉの算出は、次の式によって行われる。
【００１２】
【数４】Ｉn+1/2 ＝ΣＨn+1/2 ・・・・（４） ［Ｓ３３］引き渡された電流源値Ｉを用い、回路解析に
おける時刻ｔcs＝（ｎ＋１／２）Δｔemにおいて回路解
析を行う。回路解析は、次の式によって行う。
おける時刻ｔcs＝（ｎ＋１／２）Δｔemにおいて回路解
析を行う。回路解析は、次の式によって行う。
【００１３】
【数５】Ｖn+1/2 ＝Ｚ-1Ｉn+1/2 ・・・・（５） ［Ｓ３４］回路解析によって求め出された電圧Ｖに基づ
いて、回路が存在する領域の電界値Ｅを算出し、電磁波
解析に引き渡す。電界値Ｅは、次の式によって求める。
いて、回路が存在する領域の電界値Ｅを算出し、電磁波
解析に引き渡す。電界値Ｅは、次の式によって求める。
【００１４】
【数６】Ｅn+1 ＝Ｖn+1/2 ／ｄ ・・・・（６） ここで、ｄは回路が存在する部分のセルの一辺の長さで
ある。 ［Ｓ３５］電磁波解析における時刻ｔem＝（ｎ＋１）Δ
ｔemにおいて、電界Ｅを計算する。電界Ｅの計算は、次
の式によって行う。
ある。 ［Ｓ３５］電磁波解析における時刻ｔem＝（ｎ＋１）Δ
ｔemにおいて、電界Ｅを計算する。電界Ｅの計算は、次
の式によって行う。
【００１５】
【数７】 ここで、εは誘電率である。ステップＳ３４で引き渡さ
れた電界値を、ステップＳ３５で計算した電界に反映さ
せることで、時刻ｔem＝（ｎ＋１）Δｔemにおける電界
が求められる。
れた電界値を、ステップＳ３５で計算した電界に反映さ
せることで、時刻ｔem＝（ｎ＋１）Δｔemにおける電界
が求められる。
【００１６】このように、電磁波解析と回路解析との間
で電流源値及び電界値との受け渡しを行うことで、電磁
波解析と回路解析とを結合したシミュレーションを行う
ことができる。
で電流源値及び電界値との受け渡しを行うことで、電磁
波解析と回路解析とを結合したシミュレーションを行う
ことができる。
【００１７】
【発明が解決しようとする課題】しかし、このような手
順でシミュレーションを行った場合、式（２）を満たす
ようなΔｔemを設定すると解が不安定になってしまうお
それがある。これは、電磁波解析の磁界の計算における
時刻と電界の計算における時刻とが「Δｔem/2」だけず
れていることが考慮されていないためであると考えられ
る。すなわち、式（６）において、時刻（ｎ＋１／２）
Δｔemの電圧を用いて、時刻（ｎ＋１）Δｔemの電界を
求めているため、この時刻の誤差が解を不安定にしてい
る。ここで、式（２）を満たすようなΔｔemをさらにそ
の１／２程度まで小さくすれば解の不安定化を防ぐこと
が可能であるが、そうすると計算時間の増加（約２倍）
を招くといった問題点があった。
順でシミュレーションを行った場合、式（２）を満たす
ようなΔｔemを設定すると解が不安定になってしまうお
それがある。これは、電磁波解析の磁界の計算における
時刻と電界の計算における時刻とが「Δｔem/2」だけず
れていることが考慮されていないためであると考えられ
る。すなわち、式（６）において、時刻（ｎ＋１／２）
Δｔemの電圧を用いて、時刻（ｎ＋１）Δｔemの電界を
求めているため、この時刻の誤差が解を不安定にしてい
る。ここで、式（２）を満たすようなΔｔemをさらにそ
の１／２程度まで小さくすれば解の不安定化を防ぐこと
が可能であるが、そうすると計算時間の増加（約２倍）
を招くといった問題点があった。
【００１８】また、ｔcsの値は回路解析の状況に応じて
随時変化するため、ｔem＝ｔcsとならずにｔcsだけが増
加して解析が破綻することがしばしば起きていた。本発
明はこのような点に鑑みてなされたものであり、電磁波
解析と回路解析との融合シミュレーションを少ない計算
量で行い、かつ安定した解を得ることのできるシミュレ
ーション装置を提供することを目的とする。
随時変化するため、ｔem＝ｔcsとならずにｔcsだけが増
加して解析が破綻することがしばしば起きていた。本発
明はこのような点に鑑みてなされたものであり、電磁波
解析と回路解析との融合シミュレーションを少ない計算
量で行い、かつ安定した解を得ることのできるシミュレ
ーション装置を提供することを目的とする。
【００１９】また、本発明の他の目的は、電磁波解析と
回路解析との融合シミュレーションを少ない計算量で行
い、かつ安定した解を得ることのできるようなシミュレ
ーションプログラムを記録したコンピュータ読み取り可
能な記録媒体を提供することである。
回路解析との融合シミュレーションを少ない計算量で行
い、かつ安定した解を得ることのできるようなシミュレ
ーションプログラムを記録したコンピュータ読み取り可
能な記録媒体を提供することである。
【００２０】
【課題を解決するための手段】本発明では上記課題を解
決するために、特定の空間内の過渡的な電磁波解析と前
記空間内に配置された回路の回路解析とを時間領域で連
携したシミュレーションを行うシミュレーション装置に
おいて、磁界を計算するたびに一定の第１の時間刻み幅
で進行させる第１の時刻と、前記第１の時刻から一定時
間ずらされており、電界を計算するたびに前記第１の時
間刻み幅で進行させる第２の時刻とを定義し、磁界の計
算と電界の計算とを交互に行うことにより電磁波解析を
行っており、電界を計算する際には、前記回路が存在し
ている領域の電界値が引き渡されるのを待ち、電界値が
引き渡されたら、引き渡された電界値を反映させて前記
空間内の電界を求める電磁波解析手段と、状況に応じて
値が変化する第２の時間刻み幅で進行させる第３の時刻
を定義し、前記第３の時刻における回路方程式を解くこ
とで前記回路にかかる電圧を求めており、前記第１の時
刻と前記第３の時刻とが、予め設定されている第１の時
間差内に接近した場合には、前記回路の電流源値が引き
渡されるのを待ち、前記回路の電流源値が引き渡された
ら、渡された値よって前記回路の電流源値を更新して回
路解析を続行する回路解析手段と、前記第１の時刻と前
記第３の時刻とが、前記第１の時間差内に接近した場合
には、前記電磁波解析手段により計算された磁界に基づ
いて前記回路の電流源値を算出し、算出した電流源値を
前記回路解析手段へ引き渡す電流源値引き渡し手段と、
前記第２の時刻と前記第３の時刻とが予め設定されてい
る第２の時間差内に接近した場合には、前記回路解析手
段により求められる前記回路にかかる電圧に基づいて、
前記回路が存在している領域の電界値を算出し、算出し
た電界値を前記電磁波解析手段へ引き渡す電界値引き渡
し手段と、を有することを特徴とするシミュレーション
装置が提供される。
決するために、特定の空間内の過渡的な電磁波解析と前
記空間内に配置された回路の回路解析とを時間領域で連
携したシミュレーションを行うシミュレーション装置に
おいて、磁界を計算するたびに一定の第１の時間刻み幅
で進行させる第１の時刻と、前記第１の時刻から一定時
間ずらされており、電界を計算するたびに前記第１の時
間刻み幅で進行させる第２の時刻とを定義し、磁界の計
算と電界の計算とを交互に行うことにより電磁波解析を
行っており、電界を計算する際には、前記回路が存在し
ている領域の電界値が引き渡されるのを待ち、電界値が
引き渡されたら、引き渡された電界値を反映させて前記
空間内の電界を求める電磁波解析手段と、状況に応じて
値が変化する第２の時間刻み幅で進行させる第３の時刻
を定義し、前記第３の時刻における回路方程式を解くこ
とで前記回路にかかる電圧を求めており、前記第１の時
刻と前記第３の時刻とが、予め設定されている第１の時
間差内に接近した場合には、前記回路の電流源値が引き
渡されるのを待ち、前記回路の電流源値が引き渡された
ら、渡された値よって前記回路の電流源値を更新して回
路解析を続行する回路解析手段と、前記第１の時刻と前
記第３の時刻とが、前記第１の時間差内に接近した場合
には、前記電磁波解析手段により計算された磁界に基づ
いて前記回路の電流源値を算出し、算出した電流源値を
前記回路解析手段へ引き渡す電流源値引き渡し手段と、
前記第２の時刻と前記第３の時刻とが予め設定されてい
る第２の時間差内に接近した場合には、前記回路解析手
段により求められる前記回路にかかる電圧に基づいて、
前記回路が存在している領域の電界値を算出し、算出し
た電界値を前記電磁波解析手段へ引き渡す電界値引き渡
し手段と、を有することを特徴とするシミュレーション
装置が提供される。
【００２１】このシミュレーション装置によりシミュレ
ーションを行えば、電磁波解析手段から回路解析手段へ
の電流源値の引き渡しは、電流源値引き渡し手段によ
り、電磁波解析において磁界が計算される第１の時刻と
回路解析の第３の時刻とが接近した際に行われる。同様
に、回路解析手段から電磁波解析手段への電界値の引き
渡しは、電界値引き渡し手段により、電磁波解析におい
て電界が計算される第２の時刻と回路解析の第３の時刻
とが接近した際に行われる。その結果、電磁波解析にお
ける磁界の計算と電界の計算との時刻のずれに合わせ
た、電流源値と電界値との引き渡しが行われ、安定した
解を得ることができる。
ーションを行えば、電磁波解析手段から回路解析手段へ
の電流源値の引き渡しは、電流源値引き渡し手段によ
り、電磁波解析において磁界が計算される第１の時刻と
回路解析の第３の時刻とが接近した際に行われる。同様
に、回路解析手段から電磁波解析手段への電界値の引き
渡しは、電界値引き渡し手段により、電磁波解析におい
て電界が計算される第２の時刻と回路解析の第３の時刻
とが接近した際に行われる。その結果、電磁波解析にお
ける磁界の計算と電界の計算との時刻のずれに合わせ
た、電流源値と電界値との引き渡しが行われ、安定した
解を得ることができる。
【００２２】また、時刻を一定の時間刻み幅で進めなが
ら行う第１シミュレーションと時刻を随時変化する時間
刻み幅で進めながら行う第２シミュレーションとを、互
いのシミュレーション結果を受け渡しながら行うシミュ
レーション装置において、第１シミュレーションの時刻
と第２シミュレーションの時刻との差が所定値内の場
合、第１シミュレーションにおける結果を第２シミュレ
ーションに引き渡す第１の引き渡し手段と、第１シミュ
レーションの結果を用いて、随時時間刻み幅を変化させ
ながら第２シミュレーションを行う解析手段と、第２シ
ミュレーションの時刻と前記一定の時間刻み幅で進行し
た第１シミュレーションの時刻との差が所定値内の場
合、第２シミュレーションにおける結果を第１シミュレ
ーションの結果に引き渡す第２の引き渡し手段と、を有
することを特徴とするシミュレーション装置が提供され
る。
ら行う第１シミュレーションと時刻を随時変化する時間
刻み幅で進めながら行う第２シミュレーションとを、互
いのシミュレーション結果を受け渡しながら行うシミュ
レーション装置において、第１シミュレーションの時刻
と第２シミュレーションの時刻との差が所定値内の場
合、第１シミュレーションにおける結果を第２シミュレ
ーションに引き渡す第１の引き渡し手段と、第１シミュ
レーションの結果を用いて、随時時間刻み幅を変化させ
ながら第２シミュレーションを行う解析手段と、第２シ
ミュレーションの時刻と前記一定の時間刻み幅で進行し
た第１シミュレーションの時刻との差が所定値内の場
合、第２シミュレーションにおける結果を第１シミュレ
ーションの結果に引き渡す第２の引き渡し手段と、を有
することを特徴とするシミュレーション装置が提供され
る。
【００２３】このようなシミュレーション装置によりシ
ミュレーションを行えば、第１シミュレーションの時刻
と第２シミュレーションの時刻との差が所定値内になる
と、第１の引き渡し手段により、第１シミュレーション
における結果が第２シミュレーションに引き渡される。
また、一定の時間刻み幅で進行した第１シミュレーショ
ンの時刻と第２シミュレーションの時刻との差が所定値
内になると、第２の引き渡し手段により、第２シミュレ
ーションにおける結果が第１シミュレーションに引き渡
される。
ミュレーションを行えば、第１シミュレーションの時刻
と第２シミュレーションの時刻との差が所定値内になる
と、第１の引き渡し手段により、第１シミュレーション
における結果が第２シミュレーションに引き渡される。
また、一定の時間刻み幅で進行した第１シミュレーショ
ンの時刻と第２シミュレーションの時刻との差が所定値
内になると、第２の引き渡し手段により、第２シミュレ
ーションにおける結果が第１シミュレーションに引き渡
される。
【００２４】また、特定の空間内の過渡的な電磁波解析
と前記空間内に配置された回路の回路解析とを時間領域
で連携したシミュレーションを行うためのシミュレーシ
ョンプログラムを記録したコンピュータ読み取り可能な
記録媒体において、磁界を計算するたびに一定の第１の
時間刻み幅で進行させる第１の時刻と、前記第１の時刻
から一定時間ずらされており、電界を計算するたびに前
記第１の時間刻み幅で進行させる第２の時刻とを定義
し、磁界の計算と電界の計算とを交互に行うことにより
電磁波解析を行っており、電界を計算する際には、前記
回路が存在している領域の電界値が引き渡されるのを待
ち、電界値が引き渡されたら、引き渡された電界値を反
映させて前記空間内の電界を求める電磁波解析手段、状
況に応じて値が変化する第２の時間刻み幅で進行させる
第３の時刻を定義し、前記第３の時刻における回路方程
式を解くことで前記回路にかかる電圧を求めており、前
記第１の時刻と前記第３の時刻とが、予め設定されてい
る第１の時間差内に接近した場合には、前記回路の電流
源値が引き渡されるのを待ち、前記回路の電流源値が引
き渡されたら、渡された値よって前記回路の電流源値を
更新して回路解析を続行する回路解析手段と、前記第１
の時刻と前記第３の時刻とが、前記第１の時間差内に接
近した場合には、前記電磁波解析手段により計算された
磁界に基づいて前記回路の電流源値を算出し、算出した
電流源値を前記回路解析手段へ引き渡す電流源値引き渡
し手段と、前記第２の時刻と前記第３の時刻とが予め設
定されている第２の時間差内に接近した場合には、前記
回路解析手段により求められる前記回路にかかる電圧に
基づいて、前記回路が存在している領域の電界値を算出
し、算出した電界値を前記電磁波解析手段へ引き渡す電
界値引き渡し手段、としてコンピュータを機能させるこ
とを特徴とするシミュレーションプログラムを記録した
コンピュータ読み取り可能な記録媒体が提供される。
と前記空間内に配置された回路の回路解析とを時間領域
で連携したシミュレーションを行うためのシミュレーシ
ョンプログラムを記録したコンピュータ読み取り可能な
記録媒体において、磁界を計算するたびに一定の第１の
時間刻み幅で進行させる第１の時刻と、前記第１の時刻
から一定時間ずらされており、電界を計算するたびに前
記第１の時間刻み幅で進行させる第２の時刻とを定義
し、磁界の計算と電界の計算とを交互に行うことにより
電磁波解析を行っており、電界を計算する際には、前記
回路が存在している領域の電界値が引き渡されるのを待
ち、電界値が引き渡されたら、引き渡された電界値を反
映させて前記空間内の電界を求める電磁波解析手段、状
況に応じて値が変化する第２の時間刻み幅で進行させる
第３の時刻を定義し、前記第３の時刻における回路方程
式を解くことで前記回路にかかる電圧を求めており、前
記第１の時刻と前記第３の時刻とが、予め設定されてい
る第１の時間差内に接近した場合には、前記回路の電流
源値が引き渡されるのを待ち、前記回路の電流源値が引
き渡されたら、渡された値よって前記回路の電流源値を
更新して回路解析を続行する回路解析手段と、前記第１
の時刻と前記第３の時刻とが、前記第１の時間差内に接
近した場合には、前記電磁波解析手段により計算された
磁界に基づいて前記回路の電流源値を算出し、算出した
電流源値を前記回路解析手段へ引き渡す電流源値引き渡
し手段と、前記第２の時刻と前記第３の時刻とが予め設
定されている第２の時間差内に接近した場合には、前記
回路解析手段により求められる前記回路にかかる電圧に
基づいて、前記回路が存在している領域の電界値を算出
し、算出した電界値を前記電磁波解析手段へ引き渡す電
界値引き渡し手段、としてコンピュータを機能させるこ
とを特徴とするシミュレーションプログラムを記録した
コンピュータ読み取り可能な記録媒体が提供される。
【００２５】このシミュレーションプログラムをコンピ
ュータで実行すれば、上記の電磁波解析と回路解析とを
連携させた本発明のシミュレーション装置に必要な各機
能がコンピュータにより実現される。
ュータで実行すれば、上記の電磁波解析と回路解析とを
連携させた本発明のシミュレーション装置に必要な各機
能がコンピュータにより実現される。
【００２６】また、時刻を一定の時間刻み幅で進めなが
ら行う第１シミュレーションと時刻を随時変化する時間
刻み幅で進めながら行う第２シミュレーションとを、互
いのシミュレーション結果を受け渡しながら行うための
シミュレーションプログラムを記録したコンピュータ読
み取り可能な記録媒体において、第１シミュレーション
の時刻と第２シミュレーションの時刻との差が所定値内
の場合、第１シミュレーションの時刻における結果を第
２シミュレーションに引き渡す第１の引き渡し手段、第
１シミュレーションの結果を用いて、随時時間刻み幅を
変化させながら第２シミュレーションを行う解析手段、
第２シミュレーションの時刻と前記一定の時間刻み幅で
進行した第１シミュレーションの時刻との差が所定値内
の場合、第２シミュレーションにおける結果を第１シミ
ュレーションの結果に引き渡す第２の引き渡し手段、と
してコンピュータを機能させることを特徴とするシミュ
レーションプログラムを記録したコンピュータ読み取り
可能な記録媒体が提供される。
ら行う第１シミュレーションと時刻を随時変化する時間
刻み幅で進めながら行う第２シミュレーションとを、互
いのシミュレーション結果を受け渡しながら行うための
シミュレーションプログラムを記録したコンピュータ読
み取り可能な記録媒体において、第１シミュレーション
の時刻と第２シミュレーションの時刻との差が所定値内
の場合、第１シミュレーションの時刻における結果を第
２シミュレーションに引き渡す第１の引き渡し手段、第
１シミュレーションの結果を用いて、随時時間刻み幅を
変化させながら第２シミュレーションを行う解析手段、
第２シミュレーションの時刻と前記一定の時間刻み幅で
進行した第１シミュレーションの時刻との差が所定値内
の場合、第２シミュレーションにおける結果を第１シミ
ュレーションの結果に引き渡す第２の引き渡し手段、と
してコンピュータを機能させることを特徴とするシミュ
レーションプログラムを記録したコンピュータ読み取り
可能な記録媒体が提供される。
【００２７】このシミュレーションプログラムをコンピ
ュータで実行すれば、上記の第１シミュレーションとだ
２シミュレーションとシミュレーション結果を受け渡し
を行う本発明のシミュレーション装置に必要な各機能が
コンピュータにより実現される。
ュータで実行すれば、上記の第１シミュレーションとだ
２シミュレーションとシミュレーション結果を受け渡し
を行う本発明のシミュレーション装置に必要な各機能が
コンピュータにより実現される。
【００２８】
【発明の実施の形態】以下、本発明の実施の形態を図面
を参照して説明する。図１は、本発明の原理構成図であ
る。本発明のシミュレーション装置は、ある特定の空間
の電磁波解析と、その空間内に存在する回路の回路解析
とを結合したものであり、電磁波解析手段１、回路解析
手段２、電流源値引き渡し手段３、及び電界値引き渡し
手段４からなる。なお、ここにいう「回路」には、単な
る１個の素子も含まれる。また、以下の説明に用いる
「時刻」とは、シミュレーションにおける時刻に関する
変数であり、現実の時刻とは異なる。
を参照して説明する。図１は、本発明の原理構成図であ
る。本発明のシミュレーション装置は、ある特定の空間
の電磁波解析と、その空間内に存在する回路の回路解析
とを結合したものであり、電磁波解析手段１、回路解析
手段２、電流源値引き渡し手段３、及び電界値引き渡し
手段４からなる。なお、ここにいう「回路」には、単な
る１個の素子も含まれる。また、以下の説明に用いる
「時刻」とは、シミュレーションにおける時刻に関する
変数であり、現実の時刻とは異なる。
【００２９】電磁波解析手段１は、特定の３次元空間を
多数のセルに分割し、ＦＤＴＤ法により過渡的な電磁波
解析を行う。すなわち、分割された各セルの辺における
電界Ｅと、その周囲の磁界Ｈとの過渡的な変化を計算す
る。そのために、電磁波解析手段１には、一定の時間刻
み幅Δｔem（第１の時間刻み幅）で進行させられる時刻
ｔem01（第１の時刻）と時刻ｔem02（第２の時刻）とが
定義されている。ここで、時刻ｔem02は、時刻ｔem01か
ら一定時間（Δｔem／２）ずらされている。
多数のセルに分割し、ＦＤＴＤ法により過渡的な電磁波
解析を行う。すなわち、分割された各セルの辺における
電界Ｅと、その周囲の磁界Ｈとの過渡的な変化を計算す
る。そのために、電磁波解析手段１には、一定の時間刻
み幅Δｔem（第１の時間刻み幅）で進行させられる時刻
ｔem01（第１の時刻）と時刻ｔem02（第２の時刻）とが
定義されている。ここで、時刻ｔem02は、時刻ｔem01か
ら一定時間（Δｔem／２）ずらされている。
【００３０】そして、電磁波解析手段１は、時刻ｔem01
で磁界を計算し、時刻ｔem02で電界を計算することによ
り、過渡的な電磁波解析を行っている。ただし、電界を
計算する際には、回路が存在している領域の電界値が引
き渡されるのを待ち、電界値が引き渡されたら、引き渡
された電界値を反映させて空間内の電界を求める。ここ
では、引き渡された電界値を反映させる方法として次の
ような方法を取る。すなわち、マックスウェル方程式を
時間と空間について解くことによって、回路が存在して
いる領域以外の空間の電界を求め、回路が存在している
領域の電界値として、引き渡された電界値を当てはめ
る。
で磁界を計算し、時刻ｔem02で電界を計算することによ
り、過渡的な電磁波解析を行っている。ただし、電界を
計算する際には、回路が存在している領域の電界値が引
き渡されるのを待ち、電界値が引き渡されたら、引き渡
された電界値を反映させて空間内の電界を求める。ここ
では、引き渡された電界値を反映させる方法として次の
ような方法を取る。すなわち、マックスウェル方程式を
時間と空間について解くことによって、回路が存在して
いる領域以外の空間の電界を求め、回路が存在している
領域の電界値として、引き渡された電界値を当てはめ
る。
【００３１】回路解析手段２では、電流源法を用いた回
路部の等価回路２ａに基づいて、例えば、ＳＰＩＣＥ
（Simulation Program with Integrated Circuit Empha
sis)による回路解析を行う。この回路解析では、電流源
から回路部の両極にかかる電圧Ｖを求める。等価回路２
ａは、回路部２ａａの両極にかかる電流源２ａｂと、電
流源２ａｂと並列に接続されたコンデンサＣとで表され
ている。また、回路解析手段２には、回路方程式を解く
べき時刻ｔcs（第３の時刻）が定義されている。この時
刻ｔcsは、時間刻み幅Δｔcs（第２の時間刻み幅）で進
行させられる。この時間刻み幅Δｔcsは、状況に応じて
値が変化する。そして、時刻ｔem01と時刻ｔcsとが予め
設定されている時間差λ1 （第１の時間差）内に接近し
た場合には、回路の電流源値が引き渡されるのを待ち、
回路の電流源値が引き渡されたら、渡された値によって
回路の電流源値を更新して回路解析を続行する。
路部の等価回路２ａに基づいて、例えば、ＳＰＩＣＥ
（Simulation Program with Integrated Circuit Empha
sis)による回路解析を行う。この回路解析では、電流源
から回路部の両極にかかる電圧Ｖを求める。等価回路２
ａは、回路部２ａａの両極にかかる電流源２ａｂと、電
流源２ａｂと並列に接続されたコンデンサＣとで表され
ている。また、回路解析手段２には、回路方程式を解く
べき時刻ｔcs（第３の時刻）が定義されている。この時
刻ｔcsは、時間刻み幅Δｔcs（第２の時間刻み幅）で進
行させられる。この時間刻み幅Δｔcsは、状況に応じて
値が変化する。そして、時刻ｔem01と時刻ｔcsとが予め
設定されている時間差λ1 （第１の時間差）内に接近し
た場合には、回路の電流源値が引き渡されるのを待ち、
回路の電流源値が引き渡されたら、渡された値によって
回路の電流源値を更新して回路解析を続行する。
【００３２】また、回路解析手段２は、時刻ｔcsを時間
刻み幅Δｔcsだけ進行させた際に、時刻ｔcsが時刻ｔem
01よりも進み、かつ時刻ｔem01と時刻ｔcsとが時間差λ
1 よりも離れている場合には、時刻ｔcs及び回路の状態
を進行前の状態に戻す。そして、時間刻み幅Δｔcsの値
を小さくした上で、時刻ｔcsを時間刻み幅Δｔcsだけ再
度進行させる。同様に、時刻ｔcsを時間刻み幅Δｔcsだ
け進行させた際に、時刻ｔcsが時刻ｔem02よりも進み、
かつ時刻ｔem02と時刻ｔcsとが時間差λ2 （第２の時間
差）よりも離れている場合には、時刻ｔcs及び回路の状
態を進行前の状態に戻す。そして、時間刻み幅Δｔcsの
値を小さくした上で、時刻ｔcsを時間刻み幅Δｔcsだけ
再度進行させる。ここで、時間差λ1 、λ2 は予め設定
された十分に小さな値であり、例えば、Δｔemの１００
分の１程度の値を用いる。
刻み幅Δｔcsだけ進行させた際に、時刻ｔcsが時刻ｔem
01よりも進み、かつ時刻ｔem01と時刻ｔcsとが時間差λ
1 よりも離れている場合には、時刻ｔcs及び回路の状態
を進行前の状態に戻す。そして、時間刻み幅Δｔcsの値
を小さくした上で、時刻ｔcsを時間刻み幅Δｔcsだけ再
度進行させる。同様に、時刻ｔcsを時間刻み幅Δｔcsだ
け進行させた際に、時刻ｔcsが時刻ｔem02よりも進み、
かつ時刻ｔem02と時刻ｔcsとが時間差λ2 （第２の時間
差）よりも離れている場合には、時刻ｔcs及び回路の状
態を進行前の状態に戻す。そして、時間刻み幅Δｔcsの
値を小さくした上で、時刻ｔcsを時間刻み幅Δｔcsだけ
再度進行させる。ここで、時間差λ1 、λ2 は予め設定
された十分に小さな値であり、例えば、Δｔemの１００
分の１程度の値を用いる。
【００３３】電流源値引き渡し手段３は、時刻ｔem01と
時刻ｔcsとが、時間差λ1 内に接近した場合には、電磁
波解析手段１により計算された磁界に基づいて回路の電
流源値を算出し、算出した電流源値を回路解析手段２へ
引き渡す。
時刻ｔcsとが、時間差λ1 内に接近した場合には、電磁
波解析手段１により計算された磁界に基づいて回路の電
流源値を算出し、算出した電流源値を回路解析手段２へ
引き渡す。
【００３４】電界値引き渡し手段４は、時刻ｔem02と時
刻ｔcsとが、予め設定されている時間差λ2 内に接近し
た場合には、回路解析手段２の回路解析により求められ
る回路への電圧に基づいて、回路が存在している領域の
電界値を算出する。そして、算出した電界値を電磁波解
析手段１へ引き渡す。
刻ｔcsとが、予め設定されている時間差λ2 内に接近し
た場合には、回路解析手段２の回路解析により求められ
る回路への電圧に基づいて、回路が存在している領域の
電界値を算出する。そして、算出した電界値を電磁波解
析手段１へ引き渡す。
【００３５】以下に、電磁波解析と回路解析との間での
データの受け渡し手順を説明する。図２は、電磁波解析
と回路解析との間でのデータの受け渡し手順を示す図で
ある。なお、以下に示す数式における各記号の右肩の数
値は、ｋ×Δｔemで表される時刻の「ｋ」の値を示した
ものである（図１中の式も同様）。 ［Ｓ１］電磁波解析手段１は、時刻ｔem＝（ｎ＋１／
２）Δｔemにおける磁界を計算する。ここで、ｎは自然
数である。磁界の計算は、次の式によって行う。
データの受け渡し手順を説明する。図２は、電磁波解析
と回路解析との間でのデータの受け渡し手順を示す図で
ある。なお、以下に示す数式における各記号の右肩の数
値は、ｋ×Δｔemで表される時刻の「ｋ」の値を示した
ものである（図１中の式も同様）。 ［Ｓ１］電磁波解析手段１は、時刻ｔem＝（ｎ＋１／
２）Δｔemにおける磁界を計算する。ここで、ｎは自然
数である。磁界の計算は、次の式によって行う。
【００３６】
【数８】 ［Ｓ２］回路解析手段２は、回路解析の間に｜ｔcs−ｔ
em01｜＜λ1 になると、解析処理を中断し、電流源値の
引き渡しを待つ。そして、｜ｔcs−ｔem01｜＜λ 1 にな
った際には、電流源値引き渡し手段３が、電磁波解析手
段１によって求め出された磁界に基づいて電流源値を算
出し、その電流源値を回路解析手段２へ引き渡す。電流
源値Ｉの算出は、次の式によって行われる。
em01｜＜λ1 になると、解析処理を中断し、電流源値の
引き渡しを待つ。そして、｜ｔcs−ｔem01｜＜λ 1 にな
った際には、電流源値引き渡し手段３が、電磁波解析手
段１によって求め出された磁界に基づいて電流源値を算
出し、その電流源値を回路解析手段２へ引き渡す。電流
源値Ｉの算出は、次の式によって行われる。
【００３７】
【数９】Ｉn+1/2 ＝ΣＨn+1/2 ・・・・（９） ［Ｓ３］回路解析手段２は、引き渡された電流源値を用
い、時刻ｔcs＝ｎ＋１／２において回路解析を行う。回
路解析は、次の式によって行う。
い、時刻ｔcs＝ｎ＋１／２において回路解析を行う。回
路解析は、次の式によって行う。
【００３８】
【数１０】Ｖn+1/2 ＝Ｚ-1Ｉn+1/2 ・・・・（１０） ここで、Ｚは、回路のインピーダンス行列である。回路
解析手段２は、ｔcsの値を時間刻み幅Δｔcsずつ進めな
がら回路解析を繰り返し行う。 ［Ｓ４］｜ｔcs−ｔem02｜＜λ2 となったら、電界値引
き渡し手段４が、回路解析によって求め出された電圧に
基づいて、回路が存在する領域の電界値を算出し、電磁
波解析手段１に引き渡す。電界値は、次の式によって求
める。
解析手段２は、ｔcsの値を時間刻み幅Δｔcsずつ進めな
がら回路解析を繰り返し行う。 ［Ｓ４］｜ｔcs−ｔem02｜＜λ2 となったら、電界値引
き渡し手段４が、回路解析によって求め出された電圧に
基づいて、回路が存在する領域の電界値を算出し、電磁
波解析手段１に引き渡す。電界値は、次の式によって求
める。
【００３９】
【数１１】Ｅn+1 ＝Ｖn+1 ／ｄ ・・・・（１１） ［Ｓ５］電磁波解析手段１は電界値引き渡されると、時
刻ｔem02＝（ｎ＋１）Δｔemにおける電界を計算する。
電界の計算は、次の式によって行う。
刻ｔem02＝（ｎ＋１）Δｔemにおける電界を計算する。
電界の計算は、次の式によって行う。
【００４０】
【数１２】 そして、電磁波解析手段１は、回路が存在する領域には
引き渡された電界値を設定し、その他の領域には式（１
２）で求められた電界を設定することで、所定の空間内
の電界を求める。
引き渡された電界値を設定し、その他の領域には式（１
２）で求められた電界を設定することで、所定の空間内
の電界を求める。
【００４１】このようなシミュレーション装置によれ
ば、電磁波解析手段１から回路解析手段２への電流源値
の引き渡しは、電磁波解析において磁界が計算される時
刻ｔem01と、回路解析の時刻ｔcsとが接近した際に行わ
れる。同様に、回路解析手段２から電磁波解析手段１へ
の電界値の引き渡しは、電磁波解析において電界が計算
される時刻ｔem02と、回路解析の時刻ｔcsとが接近した
際に行われる。その結果、電磁波解析における磁界の計
算と電界の計算との時刻のずれに合わせたデータとの引
き渡しが行われる。
ば、電磁波解析手段１から回路解析手段２への電流源値
の引き渡しは、電磁波解析において磁界が計算される時
刻ｔem01と、回路解析の時刻ｔcsとが接近した際に行わ
れる。同様に、回路解析手段２から電磁波解析手段１へ
の電界値の引き渡しは、電磁波解析において電界が計算
される時刻ｔem02と、回路解析の時刻ｔcsとが接近した
際に行われる。その結果、電磁波解析における磁界の計
算と電界の計算との時刻のずれに合わせたデータとの引
き渡しが行われる。
【００４２】なお、時刻ｔcsが時刻ｔemよりも大きく進
んでしまうことのないように制御することで、シミュレ
ーションの破綻が防止されている。以下に、そのような
制御をも含めたシミュレーションの手順を具体的に説明
する。
んでしまうことのないように制御することで、シミュレ
ーションの破綻が防止されている。以下に、そのような
制御をも含めたシミュレーションの手順を具体的に説明
する。
【００４３】図３は、本発明に係るシミュレーションの
処理手順を示す図である。なお、以下の処理は基本的に
図１に示した各構成要素が行う処理であるが、ステップ
Ｓ１１はシミュレーションの全体を統括的に制御する制
御部（図示せず）が行う処理である。 ［Ｓ１１］制御部が、データ係数等の主記憶領域を確保
し、変数等の初期化をする。初期化は、具体的には回路
解析における時刻を「ｔcs＝０．０」、電磁波解析にお
ける時刻ｔem01を「ｔem01＝０．０」、電磁波解析にお
ける時刻ｔem02を「ｔem02＝０．０」とする。 ［Ｓ１２］回路解析手段２が、時間変化しない係数等を
計算する。 ［Ｓ１３］電磁波解析手段１が、電磁波解析における時
刻ｔem01を「ｔem01＝ｔ em01＋Δｔem/2」、電磁波解析
における時刻ｔem02を「ｔem02＝ｔem02＋Δｔem」とす
る。 ［Ｓ１４］回路解析手段２が、回路解析における時刻を
「ｔcs＝ｔcs＋Δｔcs」とする。 ［Ｓ１５］回路解析手段２が、「｜ｔcs−ｔem01｜＜λ
1 」が満たされているか否かを判断する。この条件が満
たされていればステップＳ１６に進み、満たされていな
ければステップＳ１９に進む。なお、このときの判断条
件を「｜ｔcs−ｔ em01｜≦λ1 」としてもよい。 ［Ｓ１６］電磁波解析手段１が、磁界の計算を行う。 ［Ｓ１７］電流源値引き渡し手段３が、電磁界解析対応
の電流源値を回路解析手段２へ引き渡す。 ［Ｓ１８］電磁波解析手段１が、電磁波解析における時
刻ｔem01を「ｔem01＝ｔ em01＋Δｔem」とする。すなわ
ち、解析の時刻を、電磁波解析の時間刻み幅Δｔ emだけ
進める。そして、ステップＳ２０に進む。 ［Ｓ１９］回路解析手段２が、「ｔcs＞ｔem01」が満た
されているか否かを判断する。すなわち、回路解析にお
ける時刻ｔcsが電磁波解析における時刻ｔem01を超えて
いないことを確認する。回路解析における時刻ｔcsが電
磁波解析における時刻ｔem01を超えている場合にはステ
ップＳ２７に進み、そうでない場合にはステップＳ２０
に進む。 ［Ｓ２０］回路解析手段２が、回路方程式を解く。これ
により、時刻ｔcsにおいて、回路の両端にかかる電圧が
求められる。 ［Ｓ２１］回路解析手段２が、回路方程式の解が収束し
たか否かを判断する。収束したのであればステップＳ２
２に進み、収束していなければステップＳ２７に進む。 ［Ｓ２２］回路解析手段２が、「｜ｔcs−ｔem02｜＜λ
2 」が満たされているか否かを判断する。この条件が満
たされていればステップＳ２３に進み、満たされていな
ければステップＳ２６に進む。なお、このときの判断条
件を「｜ｔcs−ｔ em02｜≦λ2 」としてもよい。 ［Ｓ２３］電磁波解析手段１が、電界の計算を行う。こ
の電界の計算では、回路が存在している領域を除いた部
分の電界が計算される。 ［Ｓ２４］電界値引き渡し手段４が、回路解析対応の電
界値を電磁波解析手段１に引き渡す。そして、電磁波解
析手段１が、回路が存在している領域の電界として、引
き渡された電界値を設定する。 ［Ｓ２５］電磁波解析手段１が、電磁波解析における時
刻ｔem02を「ｔem02＝ｔ em02＋Δｔem」とする。すなわ
ち、電界を計算する時刻ｔem02を、電磁波解析の時間刻
み幅Δｔemだけ進める。そして、ステップＳ２８に進
む。 ［Ｓ２６］回路解析手段２が、「ｔcs＞ｔem02」が満た
されているか否かを判断する。すなわち、回路解析にお
ける時刻ｔcsが電磁波解析における時刻ｔem02を超えて
いないことを確認する。時刻ｔcsが時刻ｔem02を超えて
いる場合にはステップＳ２７に進み、そうでない場合に
はステップＳ２８に進む。 ［Ｓ２７］回路解析手段２が、回路解析における時刻を
「ｔcs＝ｔcs−Δｔcs」とし、回路解析における時間刻
み幅を「Δｔcs＝Δｔcs×χ」とする。ここで、χは、
予め設定された０＜χ＜１の実数である。この処理の
後、ステップＳ１４に戻る。 ［Ｓ２８］回路解析手段２が、「ｔcs≧ｔmax 」が満た
される否かを判断する。この条件が満たされていれば処
理を終了し、満たされていなければステップＳ２９に進
む。ここで、ｔmax は、シミュレーション最大時刻であ
り、予め設定されている。 ［Ｓ２９］回路解析手段２は、次のΔｔcsを決定し、ス
テップＳ１４に戻る。
処理手順を示す図である。なお、以下の処理は基本的に
図１に示した各構成要素が行う処理であるが、ステップ
Ｓ１１はシミュレーションの全体を統括的に制御する制
御部（図示せず）が行う処理である。 ［Ｓ１１］制御部が、データ係数等の主記憶領域を確保
し、変数等の初期化をする。初期化は、具体的には回路
解析における時刻を「ｔcs＝０．０」、電磁波解析にお
ける時刻ｔem01を「ｔem01＝０．０」、電磁波解析にお
ける時刻ｔem02を「ｔem02＝０．０」とする。 ［Ｓ１２］回路解析手段２が、時間変化しない係数等を
計算する。 ［Ｓ１３］電磁波解析手段１が、電磁波解析における時
刻ｔem01を「ｔem01＝ｔ em01＋Δｔem/2」、電磁波解析
における時刻ｔem02を「ｔem02＝ｔem02＋Δｔem」とす
る。 ［Ｓ１４］回路解析手段２が、回路解析における時刻を
「ｔcs＝ｔcs＋Δｔcs」とする。 ［Ｓ１５］回路解析手段２が、「｜ｔcs−ｔem01｜＜λ
1 」が満たされているか否かを判断する。この条件が満
たされていればステップＳ１６に進み、満たされていな
ければステップＳ１９に進む。なお、このときの判断条
件を「｜ｔcs−ｔ em01｜≦λ1 」としてもよい。 ［Ｓ１６］電磁波解析手段１が、磁界の計算を行う。 ［Ｓ１７］電流源値引き渡し手段３が、電磁界解析対応
の電流源値を回路解析手段２へ引き渡す。 ［Ｓ１８］電磁波解析手段１が、電磁波解析における時
刻ｔem01を「ｔem01＝ｔ em01＋Δｔem」とする。すなわ
ち、解析の時刻を、電磁波解析の時間刻み幅Δｔ emだけ
進める。そして、ステップＳ２０に進む。 ［Ｓ１９］回路解析手段２が、「ｔcs＞ｔem01」が満た
されているか否かを判断する。すなわち、回路解析にお
ける時刻ｔcsが電磁波解析における時刻ｔem01を超えて
いないことを確認する。回路解析における時刻ｔcsが電
磁波解析における時刻ｔem01を超えている場合にはステ
ップＳ２７に進み、そうでない場合にはステップＳ２０
に進む。 ［Ｓ２０］回路解析手段２が、回路方程式を解く。これ
により、時刻ｔcsにおいて、回路の両端にかかる電圧が
求められる。 ［Ｓ２１］回路解析手段２が、回路方程式の解が収束し
たか否かを判断する。収束したのであればステップＳ２
２に進み、収束していなければステップＳ２７に進む。 ［Ｓ２２］回路解析手段２が、「｜ｔcs−ｔem02｜＜λ
2 」が満たされているか否かを判断する。この条件が満
たされていればステップＳ２３に進み、満たされていな
ければステップＳ２６に進む。なお、このときの判断条
件を「｜ｔcs−ｔ em02｜≦λ2 」としてもよい。 ［Ｓ２３］電磁波解析手段１が、電界の計算を行う。こ
の電界の計算では、回路が存在している領域を除いた部
分の電界が計算される。 ［Ｓ２４］電界値引き渡し手段４が、回路解析対応の電
界値を電磁波解析手段１に引き渡す。そして、電磁波解
析手段１が、回路が存在している領域の電界として、引
き渡された電界値を設定する。 ［Ｓ２５］電磁波解析手段１が、電磁波解析における時
刻ｔem02を「ｔem02＝ｔ em02＋Δｔem」とする。すなわ
ち、電界を計算する時刻ｔem02を、電磁波解析の時間刻
み幅Δｔemだけ進める。そして、ステップＳ２８に進
む。 ［Ｓ２６］回路解析手段２が、「ｔcs＞ｔem02」が満た
されているか否かを判断する。すなわち、回路解析にお
ける時刻ｔcsが電磁波解析における時刻ｔem02を超えて
いないことを確認する。時刻ｔcsが時刻ｔem02を超えて
いる場合にはステップＳ２７に進み、そうでない場合に
はステップＳ２８に進む。 ［Ｓ２７］回路解析手段２が、回路解析における時刻を
「ｔcs＝ｔcs−Δｔcs」とし、回路解析における時間刻
み幅を「Δｔcs＝Δｔcs×χ」とする。ここで、χは、
予め設定された０＜χ＜１の実数である。この処理の
後、ステップＳ１４に戻る。 ［Ｓ２８］回路解析手段２が、「ｔcs≧ｔmax 」が満た
される否かを判断する。この条件が満たされていれば処
理を終了し、満たされていなければステップＳ２９に進
む。ここで、ｔmax は、シミュレーション最大時刻であ
り、予め設定されている。 ［Ｓ２９］回路解析手段２は、次のΔｔcsを決定し、ス
テップＳ１４に戻る。
【００４４】このようにして、電磁波解析と回路解析と
を融合したシミュレーションを行うことができる。この
シミュレーションによれば、回路解析における時刻ｔcs
が電磁波解析における時刻ｔemを超えないようにしてい
るため、シミュレーションが破綻することがない。
を融合したシミュレーションを行うことができる。この
シミュレーションによれば、回路解析における時刻ｔcs
が電磁波解析における時刻ｔemを超えないようにしてい
るため、シミュレーションが破綻することがない。
【００４５】また、電磁波解析の磁界の計算における時
刻と電界の計算における時刻とが「Δｔem／２」だけず
れていることを考慮して、電流源値や電界値の受け渡し
を行っているため、安定した解を得ることができる。し
かも、電磁波解析における時間刻み幅Δｔemを必要以上
に小さくせずにすむ（すなわち式（２）を満たすΔｔ em
を設定できる）ため、計算量が過大になることもない。
刻と電界の計算における時刻とが「Δｔem／２」だけず
れていることを考慮して、電流源値や電界値の受け渡し
を行っているため、安定した解を得ることができる。し
かも、電磁波解析における時間刻み幅Δｔemを必要以上
に小さくせずにすむ（すなわち式（２）を満たすΔｔ em
を設定できる）ため、計算量が過大になることもない。
【００４６】なお、上記の回路解析における時間刻み幅
Δｔcsは状況に応じて変化するが、ΔｔcsがΔｔem／２
を超えることはないような値が設定されるようにする。
これにより、時刻ｔcsの値が、時刻ｔem01や時刻ｔem02
より大きくなるような事態の発生を少なくすることがで
きる。すなわち、余分な処理の発生を少なくできる。
Δｔcsは状況に応じて変化するが、ΔｔcsがΔｔem／２
を超えることはないような値が設定されるようにする。
これにより、時刻ｔcsの値が、時刻ｔem01や時刻ｔem02
より大きくなるような事態の発生を少なくすることがで
きる。すなわち、余分な処理の発生を少なくできる。
【００４７】次に、本発明に係るシミュレーション装置
を用いて、マイクロストリップ線路のシミュレーション
を行った場合の例を説明する。図４は、マイクロストリ
ップ線路の側面図である。解析領域のセルの大きさは全
て「ｄｘ×ｄｙ×ｄｚ＝０．２６５×０．３８９×０．
４００ｍｍ3 」である。このマイクロストリップ線路
は、誘電率ε＝２．２、Ｘ軸方向の厚さ３セル分の基板
２１の両面に薄膜の導体板３１，３２が形成されてる。
図中、導体板３１の右端と、導体板３２との間にダイオ
ード３３が設けられている。ダイオード３３のアノード
が導線（Thin Wire)３４で導体板３１に接続され、ダイ
オード３３のカソードが導線(Thin Wire) ３５で導体板
３２に接続されている。また、２つの導体板３１，３２
の左端に交流電源３４が接続され、高周波の電源３６が
供給されている。
を用いて、マイクロストリップ線路のシミュレーション
を行った場合の例を説明する。図４は、マイクロストリ
ップ線路の側面図である。解析領域のセルの大きさは全
て「ｄｘ×ｄｙ×ｄｚ＝０．２６５×０．３８９×０．
４００ｍｍ3 」である。このマイクロストリップ線路
は、誘電率ε＝２．２、Ｘ軸方向の厚さ３セル分の基板
２１の両面に薄膜の導体板３１，３２が形成されてる。
図中、導体板３１の右端と、導体板３２との間にダイオ
ード３３が設けられている。ダイオード３３のアノード
が導線（Thin Wire)３４で導体板３１に接続され、ダイ
オード３３のカソードが導線(Thin Wire) ３５で導体板
３２に接続されている。また、２つの導体板３１，３２
の左端に交流電源３４が接続され、高周波の電源３６が
供給されている。
【００４８】基板２１上面の導体板３１は、Ｚ軸方向に
２５セル分の長さを有しており、基板の下面の導体板３
２は、Ｚ軸方向に３０セル分の長さを有している。これ
らの周囲の領域は、空気２２で満たされている。また、
図中の点線で示した外側の枠（Ｘ軸に１６セル、Ｚ軸に
３０セル）は、吸収境界４１を示している。
２５セル分の長さを有しており、基板の下面の導体板３
２は、Ｚ軸方向に３０セル分の長さを有している。これ
らの周囲の領域は、空気２２で満たされている。また、
図中の点線で示した外側の枠（Ｘ軸に１６セル、Ｚ軸に
３０セル）は、吸収境界４１を示している。
【００４９】図５は、マイクロストリップ線路の断面図
である。導体板３１はＹ軸方向に６セル分の幅を有して
おり、その中央にダイード３３が接続されてる。導体板
３２はＹ軸方向に２０セル分の幅を有しており、その中
央にダーオード３３が接続されている。また、吸収境界
４１のＹ軸方向の幅は２０セル分である。
である。導体板３１はＹ軸方向に６セル分の幅を有して
おり、その中央にダイード３３が接続されてる。導体板
３２はＹ軸方向に２０セル分の幅を有しており、その中
央にダーオード３３が接続されている。また、吸収境界
４１のＹ軸方向の幅は２０セル分である。
【００５０】このようなマイクロストリップ線路のシミ
ュレーションを、従来のシミュレーション装置で行った
場合と、本発明のシミュレーション装置で行った場合と
の結果を以下に示す。
ュレーションを、従来のシミュレーション装置で行った
場合と、本発明のシミュレーション装置で行った場合と
の結果を以下に示す。
【００５１】図６は、従来のシミュレーション装置によ
り、Δｔem＝０．３１５ｐｓでシミュレーションした結
果を示す図である。この図では、横軸が時間（ｐｓ）、
縦軸がダイオードにかかる電圧である（以下の図７、図
８も同様）。このように、Δｔem＝０．３１５ｐｓに設
定すれば、従来のシミュレーション装置を用いても安定
した結果を得ることができる。ところが、Δｔem＝０．
３１５ｐｓという時間刻み幅は、電磁波解析においては
非常に短い時間刻み幅であり、計算量が過大となってい
る。
り、Δｔem＝０．３１５ｐｓでシミュレーションした結
果を示す図である。この図では、横軸が時間（ｐｓ）、
縦軸がダイオードにかかる電圧である（以下の図７、図
８も同様）。このように、Δｔem＝０．３１５ｐｓに設
定すれば、従来のシミュレーション装置を用いても安定
した結果を得ることができる。ところが、Δｔem＝０．
３１５ｐｓという時間刻み幅は、電磁波解析においては
非常に短い時間刻み幅であり、計算量が過大となってい
る。
【００５２】そこで、時間刻み幅を２倍にして従来のシ
ミュレーション装置によりシミュレーションを行った。
図７は、従来のシミュレーション装置により、Δｔem＝
０．６３０ｐｓでシミュレーションした結果を示す図で
ある。このように、時間刻み幅を２倍にすると、解析の
時刻が進行するにしたがって電圧の振動が増大し、不安
定になってしまう。
ミュレーション装置によりシミュレーションを行った。
図７は、従来のシミュレーション装置により、Δｔem＝
０．６３０ｐｓでシミュレーションした結果を示す図で
ある。このように、時間刻み幅を２倍にすると、解析の
時刻が進行するにしたがって電圧の振動が増大し、不安
定になってしまう。
【００５３】そこで、本発明のシミュレーション装置を
用いて、図７の場合と同じ時間刻み幅でシミュレーショ
ンをした。図８は、本発明に係るミュレーション装置に
より、Δｔem＝０．６３０ｐｓでシミュレーションした
結果を示す図である。このように、本発明のシミュレー
ション装置を用いれば、電磁波解析における時間刻み幅
を大きくとっても安定した解析結果を得られる。すなわ
ち、回路解析と融合させるからといって、電磁波解析に
おける時間刻み幅を小さくする必要がない。その結果、
少ない計算量で安定したシミュレーションが可能となっ
ている。
用いて、図７の場合と同じ時間刻み幅でシミュレーショ
ンをした。図８は、本発明に係るミュレーション装置に
より、Δｔem＝０．６３０ｐｓでシミュレーションした
結果を示す図である。このように、本発明のシミュレー
ション装置を用いれば、電磁波解析における時間刻み幅
を大きくとっても安定した解析結果を得られる。すなわ
ち、回路解析と融合させるからといって、電磁波解析に
おける時間刻み幅を小さくする必要がない。その結果、
少ない計算量で安定したシミュレーションが可能となっ
ている。
【００５４】なお、上記の処理機能は、コンピュータに
よって実現することができる。その場合、シミュレーシ
ョン装置が有すべき機能の処理内容は、コンピュータで
読み取り可能な記録媒体に記録されたプログラムに記述
されており、このプログラムをコンピュータで実行する
ことにより、上記処理がコンピュータで実現される。コ
ンピュータで読み取り可能な記録媒体としては、磁気記
録装置や半導体メモリ等がある。市場を流通させる場合
には、ＣＤ−ＲＯＭ(Compact Disk Read OnlyMemory)
やフロッピーディスク等の可搬型記録媒体にプログラム
を格納して流通させたり、ネットワークを介して接続さ
れたコンピュータの記憶装置に格納しておき、ネットワ
ークを通じて他のコンピュータに転送することもでき
る。コンピュータで実行する際には、コンピュータ内の
ハードディスク装置等にプログラムを格納しておき、メ
インメモリにロードして実行する。
よって実現することができる。その場合、シミュレーシ
ョン装置が有すべき機能の処理内容は、コンピュータで
読み取り可能な記録媒体に記録されたプログラムに記述
されており、このプログラムをコンピュータで実行する
ことにより、上記処理がコンピュータで実現される。コ
ンピュータで読み取り可能な記録媒体としては、磁気記
録装置や半導体メモリ等がある。市場を流通させる場合
には、ＣＤ−ＲＯＭ(Compact Disk Read OnlyMemory)
やフロッピーディスク等の可搬型記録媒体にプログラム
を格納して流通させたり、ネットワークを介して接続さ
れたコンピュータの記憶装置に格納しておき、ネットワ
ークを通じて他のコンピュータに転送することもでき
る。コンピュータで実行する際には、コンピュータ内の
ハードディスク装置等にプログラムを格納しておき、メ
インメモリにロードして実行する。
【００５５】
【発明の効果】以上説明したように本発明の電磁波解析
と回路解析とを連携させたシミュレーション装置では、
回路解析の時刻が、電磁波解析において電界を求めるべ
き時刻に接近した際に、回路解析に基づく電界値を電磁
波解析へ引き渡すようにしたため、引き渡された電界値
を求めた時刻と、その電界値を反映させる電界の時刻と
の差が微少となり、安定した解析結果を得ることができ
る。
と回路解析とを連携させたシミュレーション装置では、
回路解析の時刻が、電磁波解析において電界を求めるべ
き時刻に接近した際に、回路解析に基づく電界値を電磁
波解析へ引き渡すようにしたため、引き渡された電界値
を求めた時刻と、その電界値を反映させる電界の時刻と
の差が微少となり、安定した解析結果を得ることができ
る。
【００５６】また、本発明の第１シミュレーションとだ
２シミュレーションとシミュレーション結果を受け渡し
を行うシミュレーション装置では、第１シミュレーショ
ンの時刻と第２シミュレーションの時刻との差が所定値
内の場合にのみシミュレーション結果の受け渡しを行う
ようにしたため、相手側のシミュレーション結果を利用
する際の時刻のずれが少なくなり、安定したシミュレー
ションが可能となる。
２シミュレーションとシミュレーション結果を受け渡し
を行うシミュレーション装置では、第１シミュレーショ
ンの時刻と第２シミュレーションの時刻との差が所定値
内の場合にのみシミュレーション結果の受け渡しを行う
ようにしたため、相手側のシミュレーション結果を利用
する際の時刻のずれが少なくなり、安定したシミュレー
ションが可能となる。
【００５７】また、本発明の電磁波解析と回路解析とを
連携させたシミュレーションプログラムを記録したコン
ピュータ読み取り可能な記録媒体では、記録されたシミ
ュレーションプログラムをコンピュータに実行させるこ
とで、電磁波解析と回路解析とを融合して、安定したシ
ミュレーションを行うことができる。
連携させたシミュレーションプログラムを記録したコン
ピュータ読み取り可能な記録媒体では、記録されたシミ
ュレーションプログラムをコンピュータに実行させるこ
とで、電磁波解析と回路解析とを融合して、安定したシ
ミュレーションを行うことができる。
【００５８】また、本発明の第１シミュレーションとだ
２シミュレーションとシミュレーション結果を受け渡し
を行うためのシミュレーションプログラムを記録したコ
ンピュータ読み取り可能な記録媒体では、記録されたシ
ミュレーションプログラムをコンピュータに実行させる
ことで、２つのシミュレーションを時間領域で連携さ
せ、安定したシミュレーションを行うことができる。
２シミュレーションとシミュレーション結果を受け渡し
を行うためのシミュレーションプログラムを記録したコ
ンピュータ読み取り可能な記録媒体では、記録されたシ
ミュレーションプログラムをコンピュータに実行させる
ことで、２つのシミュレーションを時間領域で連携さ
せ、安定したシミュレーションを行うことができる。
【図１】本発明の原理構成図である。
【図２】電磁波解析と回路解析との間でのデータの受け
渡し手順を示す図である。
渡し手順を示す図である。
【図３】本発明に係るシミュレーションの処理手順を示
す図である。
す図である。
【図４】マイクロストリップ線路の側面図である。
【図５】マイクロストリップ線路の断面図である。
【図６】従来のシミュレーション装置により、Δｔem＝
０．３１５ｐｓでシミュレーションした結果を示す図で
ある。
０．３１５ｐｓでシミュレーションした結果を示す図で
ある。
【図７】従来のシミュレーション装置により、Δｔem＝
０．６３０ｐｓでシミュレーションした結果を示す図で
ある。
０．６３０ｐｓでシミュレーションした結果を示す図で
ある。
【図８】本発明に係るミュレーション装置により、Δｔ
em＝０．６３０ｐｓでシミュレーションした結果を示す
図である。
em＝０．６３０ｐｓでシミュレーションした結果を示す
図である。
【図９】従来のシミュレーション装置によるデータの引
き渡しタイミングを示す図である。
き渡しタイミングを示す図である。
１ 電磁波解析手段 ２ 回路解析手段 ２ａ 等価回路 ３ 電流源値引き渡し手段 ４ 電界値引き渡し手段
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.6 識別記号 ＦＩ Ｇ０６Ｆ 15/60 ６６６Ｚ
Claims (6)
- 【請求項１】 特定の空間内の過渡的な電磁波解析と前
記空間内に配置された回路の回路解析とを時間領域で連
携したシミュレーションを行うシミュレーション装置に
おいて、 磁界を計算するたびに一定の第１の時間刻み幅で進行さ
せる第１の時刻と、前記第１の時刻から一定時間ずらさ
れており、電界を計算するたびに前記第１の時間刻み幅
で進行させる第２の時刻とを定義し、磁界の計算と電界
の計算とを交互に行うことにより電磁波解析を行ってお
り、電界を計算する際には、前記回路が存在している領
域の電界値が引き渡されるのを待ち、電界値が引き渡さ
れたら、引き渡された電界値を反映させて前記空間内の
電界を求める電磁波解析手段と、 状況に応じて値が変化する第２の時間刻み幅で進行させ
る第３の時刻を定義し、前記第３の時刻における回路方
程式を解くことで前記回路にかかる電圧を求めており、
前記第１の時刻と前記第３の時刻とが、予め設定されて
いる第１の時間差内に接近した場合には、前記回路の電
流源値が引き渡されるのを待ち、前記回路の電流源値が
引き渡されたら、渡された値によって前記回路の電流源
値を更新して回路解析を続行する回路解析手段と、 前記第１の時刻と前記第３の時刻とが前記第１の時間差
内に接近した場合には、前記電磁波解析手段により計算
された磁界に基づいて前記回路の電流源値を算出し、算
出した電流源値を前記回路解析手段へ引き渡す電流源値
引き渡し手段と、 前記第２の時刻と前記第３の時刻とが予め設定されてい
る第２の時間差内に接近した場合には、前記回路解析手
段により求められる前記回路にかかる電圧に基づいて、
前記回路が存在している領域の電界値を算出し、算出し
た電界値を前記電磁波解析手段へ引き渡す電界値引き渡
し手段と、 を有することを特徴とするシミュレーション装置。 - 【請求項２】 前記回路解析手段は、前記第３の時刻を
前記第２の時間刻み幅だけ進行させた際に、前記第３の
時刻が前記第１の時刻よりも進み、かつ前記第１の時刻
と前記第３の時刻とが前記第１の時間差よりも離れてい
る場合には、前記第３の時刻及び回路の状態を進行前の
状態に戻し、前記第２の時間刻み幅の値を小さくした上
で、前記第３の時刻を前記第２の時間刻み幅だけ再度進
行させるとともに、前記第３の時刻を前記第２の時間刻
み幅だけ進行させた際に、前記第３の時刻が前記第２の
時刻よりも進み、かつ前記第２の時刻と前記第３の時刻
とが前記第２の時間差よりも離れている場合には、前記
第３の時刻及び回路の状態を進行前の状態に戻し、前記
第２の時間刻み幅の値を小さくした上で、前記第３の時
刻を前記第２の時間刻み幅だけ再度進行させることを特
徴とする請求項１記載のシミュレーション装置。 - 【請求項３】 前記電磁波解析手段は、前記第１の時刻
と前記第２の時刻とを、前記第１の時間刻み幅の半分の
時間だけずらしており、 前記回路解析手段は、前記第２の時間刻み幅が、前記第
１の時間刻み幅の半分の時間を超えることのないよう
に、前記第２の時間刻み幅の値を決定していることを特
徴とする請求項１記載のシミュレーション装置。 - 【請求項４】 時刻を一定の時間刻み幅で進めながら行
う第１シミュレーションと時刻を随時変化する時間刻み
幅で進めながら行う第２シミュレーションとを、互いの
シミュレーション結果を受け渡しながら行うシミュレー
ション装置において、 第１シミュレーションの時刻と第２シミュレーションの
時刻との差が所定値内の場合、第１シミュレーションに
おける結果を第２シミュレーションに引き渡す第１の引
き渡し手段と、 第１シミュレーションの結果を用いて、随時時間刻み幅
を変化させながら第２シミュレーションを行う解析手段
と、 第２シミュレーションの時刻と前記一定の時間刻み幅で
進行した第１シミュレーションの時刻との差が所定値内
の場合、第２シミュレーションにおける結果を第１シミ
ュレーションの結果に引き渡す第２の引き渡し手段と、 を有することを特徴とするシミュレーション装置。 - 【請求項５】 特定の空間内の過渡的な電磁波解析と前
記空間内に配置された回路の回路解析とを時間領域で連
携したシミュレーションを行うためのシミュレーション
プログラムを記録したコンピュータ読み取り可能な記録
媒体において、 磁界を計算するたびに一定の第１の時間刻み幅で進行さ
せる第１の時刻と、前記第１の時刻から一定時間ずらさ
れており、電界を計算するたびに前記第１の時間刻み幅
で進行させる第２の時刻とを定義し、磁界の計算と電界
の計算とを交互に行うことにより電磁波解析を行ってお
り、電界を計算する際には、前記回路が存在している領
域の電界値が引き渡されるのを待ち、電界値が引き渡さ
れたら、引き渡された電界値を反映させて前記空間内の
電界を求める電磁波解析手段、 状況に応じて値が変化する第２の時間刻み幅で進行させ
る第３の時刻を定義し、前記第３の時刻における回路方
程式を解くことで前記回路にかかる電圧を求めており、
前記第１の時刻と前記第３の時刻とが予め設定されてい
る第１の時間差内に接近した場合には、前記回路の電流
源値が引き渡されるのを待ち、前記回路の電流源値が引
き渡されたら、渡された値によって前記回路の電流源値
を更新して回路解析を続行する回路解析手段、 前記第１の時刻と前記第３の時刻とが、前記第１の時間
差内に接近した場合には、前記電磁波解析手段により計
算された磁界に基づいて前記回路の電流源値を算出し、
算出した電流源値を前記回路解析手段へ引き渡す電流源
値引き渡し手段、 前記第２の時刻と前記第３の時刻とが予め設定されてい
る第２の時間差内に接近した場合には、前記回路解析手
段により求められる前記回路にかかる電圧に基づいて、
前記回路が存在している領域の電界値を算出し、算出し
た電界値を前記電磁波解析手段へ引き渡す電界値引き渡
し手段、 としてコンピュータを機能させることを特徴とするシミ
ュレーションプログラムを記録したコンピュータ読み取
り可能な記録媒体。 - 【請求項６】 時刻を一定の時間刻み幅で進めながら行
う第１シミュレーションと時刻を随時変化する時間刻み
幅で進めながら行う第２シミュレーションとを、互いの
シミュレーション結果を受け渡しながら行うためのシミ
ュレーションプログラムを記録したコンピュータ読み取
り可能な記録媒体において、 第１シミュレーションの時刻と第２シミュレーションの
時刻との差が所定値内の場合、第１シミュレーションの
時刻における結果を第２シミュレーションに引き渡す第
１の引き渡し手段、 第１シミュレーションの結果を用いて、随時時間刻み幅
を変化させながら第２シミュレーションを行う解析手
段、 第２シミュレーションの時刻と前記一定の時間刻み幅で
進行した第１シミュレーションの時刻との差が所定値内
の場合、第２シミュレーションにおける結果を第１シミ
ュレーションの結果に引き渡す第２の引き渡し手段、 としてコンピュータを機能させることを特徴とするシミ
ュレーションプログラムを記録したコンピュータ読み取
り可能な記録媒体。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP31789197A JP3633765B2 (ja) | 1997-11-19 | 1997-11-19 | シミュレーション装置及びシミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体 |
| US09/076,069 US6219629B1 (en) | 1997-11-19 | 1998-05-12 | Apparatus for combined simulation of electromagnetic wave analysis and circuit analysis, and computer-readable medium containing simulation program therefor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP31789197A JP3633765B2 (ja) | 1997-11-19 | 1997-11-19 | シミュレーション装置及びシミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH11153634A true JPH11153634A (ja) | 1999-06-08 |
| JP3633765B2 JP3633765B2 (ja) | 2005-03-30 |
Family
ID=18093214
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP31789197A Expired - Fee Related JP3633765B2 (ja) | 1997-11-19 | 1997-11-19 | シミュレーション装置及びシミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US6219629B1 (ja) |
| JP (1) | JP3633765B2 (ja) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100365822B1 (ko) * | 1998-10-26 | 2002-12-26 | 가부시끼가이샤 히다치 세이사꾸쇼 | 회전기의 전자장해석방법 및 전자장해석장치 |
| WO2004109560A1 (ja) * | 2003-06-05 | 2004-12-16 | Fujitsu Limited | 回路解析と電磁波解析との連携による電磁界強度算出装置および算出方法 |
| JP2008040756A (ja) * | 2006-08-04 | 2008-02-21 | Sharp Corp | 電磁界回路連携解析プログラム、電磁界回路連携解析プログラムを格納した記録媒体、電磁界回路連携解析装置および電磁界回路連携解析方法 |
| JP2008059153A (ja) * | 2006-08-30 | 2008-03-13 | Sharp Corp | モデル作成プログラム、モデル作成プログラムを格納したコンピュータ読み取り可能な記録媒体およびモデル作成装置 |
| WO2008081544A1 (ja) | 2006-12-28 | 2008-07-10 | Fujitsu Limited | シミュレーションシステム、シミュレーションプログラムおよび該プログラムを記録した記録媒体 |
| JP2008217327A (ja) * | 2007-03-02 | 2008-09-18 | Sharp Corp | 電磁界回路連携解析装置、電磁界回路連携解析プログラム、電磁界回路連携解析プログラムを格納した記録媒体および電磁界回路連携解析方法 |
| JP2009199574A (ja) * | 2008-01-23 | 2009-09-03 | Fujitsu Ltd | シミュレーションプログラム、該プログラムを記録した記録媒体、シミュレーション装置、およびシミュレーション方法 |
| US8019578B2 (en) | 2007-10-16 | 2011-09-13 | Sharp Kabushiki Kaisha | Simulation method of electromagnetic field and circuit analysis on first and second targets coupled together by a circuit element simulation apparatus, and computer-readable medium storing simulation program for performing the method |
| EP2378446A2 (en) | 2009-10-07 | 2011-10-19 | Fujitsu Limited | Coupled analysis simulation apparatus and coupled analysis simulation method |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3878345B2 (ja) * | 1998-12-01 | 2007-02-07 | 富士通株式会社 | シミュレーション装置及び方法並びにプログラム記録媒体 |
| IT1320393B1 (it) * | 2000-06-05 | 2003-11-26 | Cselt Centro Studi Lab Telecom | Procedimento per il modellamento elettromagnetico di componenti esistemi elettronici. |
| JP3943326B2 (ja) * | 2000-11-27 | 2007-07-11 | 松下電器産業株式会社 | 不要輻射解析方法および不要輻射解析装置 |
| JP3818874B2 (ja) | 2001-06-26 | 2006-09-06 | 富士通株式会社 | 電磁波解析装置および電磁波解析プログラム |
| JP2003216681A (ja) * | 2002-01-24 | 2003-07-31 | Fujitsu Ltd | 電磁界強度算出方法および算出装置 |
| DE10218404B4 (de) * | 2002-04-24 | 2007-09-13 | Infineon Technologies Ag | Verfahren zur numerischen Simulation einer elektrischen Schaltung und Trägermedium |
| US6931353B2 (en) * | 2002-08-16 | 2005-08-16 | Ispat Inland Inc. | Method of predicting motor performance based on CRML material properties |
| JP2005100067A (ja) * | 2003-09-24 | 2005-04-14 | Fujitsu Ltd | マイクロ磁化解析プログラムおよび解析装置 |
| GB2407660A (en) * | 2003-10-31 | 2005-05-04 | Hewlett Packard Development Co | Distributing components across resources by modelling physical properties |
| US20070113219A1 (en) * | 2005-11-17 | 2007-05-17 | Microsoft Corporation | Representing simulation values of variable in sharpley limited time and space |
| AU2010244981A1 (en) * | 2009-05-08 | 2011-12-15 | Hyperspace Pty Ltd | A useable electromagnetic blueprint of the structure of space |
| US10380292B1 (en) * | 2016-04-08 | 2019-08-13 | Cadence Design Systems, Inc. | Systems and methods for finite difference time domain simulation of an electronic design |
| CN108717500B (zh) * | 2018-05-24 | 2021-10-08 | 中国科学院电子学研究所 | 考虑杂散参数的瞬变电磁发射***建模分析方法 |
| CN110703028B (zh) * | 2019-11-01 | 2021-11-05 | 中国人民解放军军事科学院国防工程研究院工程防护研究所 | 一种磁场线圈对有界波模拟器空间场影响的等效电路分析模型 |
| JP2021110997A (ja) * | 2020-01-07 | 2021-08-02 | 住友重機械工業株式会社 | シミュレーション方法、シミュレーション装置、及びプログラム |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5995733A (en) * | 1997-01-27 | 1999-11-30 | Lucent Technologies, Inc. | Method and apparatus for efficient design and analysis of integrated circuits using multiple time scales |
| JPH1115814A (ja) * | 1997-06-26 | 1999-01-22 | Fujitsu Ltd | モーメント法を用いたシミュレーション装置及び方法並びにプログラム記憶媒体 |
-
1997
- 1997-11-19 JP JP31789197A patent/JP3633765B2/ja not_active Expired - Fee Related
-
1998
- 1998-05-12 US US09/076,069 patent/US6219629B1/en not_active Expired - Lifetime
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100365822B1 (ko) * | 1998-10-26 | 2002-12-26 | 가부시끼가이샤 히다치 세이사꾸쇼 | 회전기의 전자장해석방법 및 전자장해석장치 |
| WO2004109560A1 (ja) * | 2003-06-05 | 2004-12-16 | Fujitsu Limited | 回路解析と電磁波解析との連携による電磁界強度算出装置および算出方法 |
| US7027941B2 (en) | 2003-06-05 | 2006-04-11 | Fujitsu Limited | Device and method for calculating electro-magnetic field intensity by cooperation of circuit analyses and electro-magnetic wave analyses |
| JP2008040756A (ja) * | 2006-08-04 | 2008-02-21 | Sharp Corp | 電磁界回路連携解析プログラム、電磁界回路連携解析プログラムを格納した記録媒体、電磁界回路連携解析装置および電磁界回路連携解析方法 |
| JP2008059153A (ja) * | 2006-08-30 | 2008-03-13 | Sharp Corp | モデル作成プログラム、モデル作成プログラムを格納したコンピュータ読み取り可能な記録媒体およびモデル作成装置 |
| JP4684188B2 (ja) * | 2006-08-30 | 2011-05-18 | シャープ株式会社 | 回路モデル作成プログラム、回路モデル作成プログラムを格納したコンピュータ読み取り可能な記録媒体および回路モデル作成装置 |
| JP4644740B2 (ja) * | 2006-12-28 | 2011-03-02 | 富士通株式会社 | シミュレーションシステム、およびシミュレーションプログラム |
| JPWO2008081544A1 (ja) * | 2006-12-28 | 2010-04-30 | 富士通株式会社 | シミュレーションシステム、シミュレーションプログラムおよび該プログラムを記録した記録媒体 |
| WO2008081544A1 (ja) | 2006-12-28 | 2008-07-10 | Fujitsu Limited | シミュレーションシステム、シミュレーションプログラムおよび該プログラムを記録した記録媒体 |
| US8060313B2 (en) | 2006-12-28 | 2011-11-15 | Fujitsu Limited | Simulation system using electromagnetic field analysis and circuit analysis and recording medium storing simulation program |
| JP2008217327A (ja) * | 2007-03-02 | 2008-09-18 | Sharp Corp | 電磁界回路連携解析装置、電磁界回路連携解析プログラム、電磁界回路連携解析プログラムを格納した記録媒体および電磁界回路連携解析方法 |
| US8019578B2 (en) | 2007-10-16 | 2011-09-13 | Sharp Kabushiki Kaisha | Simulation method of electromagnetic field and circuit analysis on first and second targets coupled together by a circuit element simulation apparatus, and computer-readable medium storing simulation program for performing the method |
| JP2009199574A (ja) * | 2008-01-23 | 2009-09-03 | Fujitsu Ltd | シミュレーションプログラム、該プログラムを記録した記録媒体、シミュレーション装置、およびシミュレーション方法 |
| EP2378446A2 (en) | 2009-10-07 | 2011-10-19 | Fujitsu Limited | Coupled analysis simulation apparatus and coupled analysis simulation method |
| US8768677B2 (en) | 2009-10-07 | 2014-07-01 | Fujitsu Limited | Coupled analysis simulation apparatus and coupled analysis simulation method |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3633765B2 (ja) | 2005-03-30 |
| US6219629B1 (en) | 2001-04-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH11153634A (ja) | シミュレーション装置及びシミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体 | |
| Yefet et al. | A staggered fourth-order accurate explicit finite difference scheme for the time-domain Maxwell's equations | |
| Li et al. | A numerical study of electro-migration voiding by evolving level set functions on a fixed cartesian grid | |
| Haney et al. | Variational methods for studying tokamak stability in the presence of a resistive wall | |
| US6789241B2 (en) | Methodology for determining the placement of decoupling capacitors in a power distribution system | |
| Fang et al. | A locally conformed finite-difference time-domain algorithm of modeling arbitrary shape planar metal strips | |
| JP2002518763A (ja) | コンピュータ・システムを使用した電力分配システムのために所望の減結合構成要素を決定するシステムおよび方法 | |
| JP2002541531A (ja) | コンピュータ・システムを使用した電力分配システムのために所望の減結合構成要素を決定するシステムおよび方法 | |
| DE102005009503A1 (de) | Modellieren von Stromflüssen in dreidimensionalen leitfähigen und dielektrischen Körpern | |
| Ruehli et al. | Stable time domain solutions for EMC problems using PEEC circuit models | |
| US6850878B2 (en) | System and method for determining the required decoupling capacitors for a power distribution system using an improved capacitor model | |
| Geiss et al. | A regularization scheme for explicit level-set XFEM topology optimization | |
| Henshaw | Automatic grid generation | |
| AU2021101629A4 (en) | Boundary truncation layer method and apparatus for three-dimensional forward modelling of low-frequency | |
| CN115293356A (zh) | 实现双量子位门的方法、电路以及量子器件、量子芯片 | |
| Unfer et al. | Multi-scale gas discharge simulations using asynchronous adaptive mesh refinement | |
| JP4647886B2 (ja) | 電磁界強度算出装置、電磁界強度算出方法および電磁界強度算出プログラム | |
| Cao et al. | Top-layer inductance extraction for the pre-layout power integrity using the physics-based model size reduction (PMSR) method | |
| Sekine et al. | Stabilized mixed finite-element time-domain method for fast transient analysis of multiscale electromagnetic problems | |
| JP4459171B2 (ja) | 電磁界回路連携解析プログラム、記録媒体、および解析装置 | |
| Weeber et al. | Structural design optimization as a technology source for developments in the electromagnetics domain | |
| de Valcárcel et al. | Photorefractive four-wave mixing cannot be treated perturbatively. Revisiting the undepleted-pump approximation | |
| JPH10307160A (ja) | 電磁波解析装置及び電磁波解析プログラムを記録したコンピュータ読み取り可能な記録媒体 | |
| US10621387B2 (en) | On-die decoupling capacitor area optimization | |
| Naojima et al. | Stabilized Explicit Isogeometric Analysis (SE-IGA) for Efficient Modeling and Analysis of 2-D Curved Structures |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040811 |
|
| 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: 20041221 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20041221 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080107 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090107 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100107 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110107 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110107 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120107 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130107 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130107 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140107 Year of fee payment: 9 |
|
| LAPS | Cancellation because of no payment of annual fees |