JP3428232B2 - 電磁界強度算出装置 - Google Patents
電磁界強度算出装置Info
- Publication number
- JP3428232B2 JP3428232B2 JP14994795A JP14994795A JP3428232B2 JP 3428232 B2 JP3428232 B2 JP 3428232B2 JP 14994795 A JP14994795 A JP 14994795A JP 14994795 A JP14994795 A JP 14994795A JP 3428232 B2 JP3428232 B2 JP 3428232B2
- Authority
- JP
- Japan
- Prior art keywords
- current
- electromagnetic field
- processing means
- display processing
- calculated
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Electromagnetism (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
【０００１】
【産業上の利用分野】本発明は，モーメント法に基づい
て電気回路装置の放射する電磁界強度を算出する電磁界
強度算出装置に関し，特に，電磁界強度の算出において
計算した解析対象としている電気回路装置に流れる電流
を可視化して，わかりやすく表示することを可能とした
電磁界強度算出装置に関する。
て電気回路装置の放射する電磁界強度を算出する電磁界
強度算出装置に関し，特に，電磁界強度の算出において
計算した解析対象としている電気回路装置に流れる電流
を可視化して，わかりやすく表示することを可能とした
電磁界強度算出装置に関する。
【０００２】電気回路装置から不要放射される電波は，
テレビやラジオ等の他の電波と干渉することから，最
近，各国で厳しく規制されるようになってきた。このよ
うな規制の規格として，日本ではＶＣＣＩ規格があり，
米国ではＦＣＣ規格があり，ドイツではＶＤＥ規格があ
る。
テレビやラジオ等の他の電波と干渉することから，最
近，各国で厳しく規制されるようになってきた。このよ
うな規制の規格として，日本ではＶＣＣＩ規格があり，
米国ではＦＣＣ規格があり，ドイツではＶＤＥ規格があ
る。
【０００３】このような電波規制を満足させるために
は，シールド技術やフィルタ技術等のような種々の対策
技術を使う必要があり，これらの対策技術がどの程度電
波を減少させるかを定量的にシミュレートする技術が必
要となる。このような電磁波解析のシミュレーションで
は，電気回路装置の放射する電磁界強度を算出するにあ
たって，電気回路装置の各要素に流れる電流を計算する
が，この電流の流れる様子を適切な方法で可視化できれ
ば，計算の不具合や電磁波の発生源となっている位置の
把握が可能となり，さらに放射メカニズムの解明にも極
めて有用であると考えられる。
は，シールド技術やフィルタ技術等のような種々の対策
技術を使う必要があり，これらの対策技術がどの程度電
波を減少させるかを定量的にシミュレートする技術が必
要となる。このような電磁波解析のシミュレーションで
は，電気回路装置の放射する電磁界強度を算出するにあ
たって，電気回路装置の各要素に流れる電流を計算する
が，この電流の流れる様子を適切な方法で可視化できれ
ば，計算の不具合や電磁波の発生源となっている位置の
把握が可能となり，さらに放射メカニズムの解明にも極
めて有用であると考えられる。
【０００４】
【従来の技術】任意形状の物体から放射される電磁界強
度は，物体各部に流れる電流が分かれば，公知の理論式
を用いて容易に計算することができる。この電流値は，
理論的には，マックスウェルの電磁波動方程式を与えら
れた境界条件の下で解くことで得られるが，任意形状物
体を対象とした複雑な境界条件の下での直接的な数式解
は現在知られていない。
度は，物体各部に流れる電流が分かれば，公知の理論式
を用いて容易に計算することができる。この電流値は，
理論的には，マックスウェルの電磁波動方程式を与えら
れた境界条件の下で解くことで得られるが，任意形状物
体を対象とした複雑な境界条件の下での直接的な数式解
は現在知られていない。
【０００５】したがって，現在の電磁界強度算出装置で
用いられている電流を求める解法は，難易の程度はある
にせよ，すべて近似的なものである。現在，この近似的
な解法の代表的なものとして，微小ループアンテナ近似
法と，分布定数線路近似法と，モーメント法という３種
類のものが知られている。
用いられている電流を求める解法は，難易の程度はある
にせよ，すべて近似的なものである。現在，この近似的
な解法の代表的なものとして，微小ループアンテナ近似
法と，分布定数線路近似法と，モーメント法という３種
類のものが知られている。
【０００６】微小ループアンテナ近似法は，波源回路と
負荷回路との間を接続する配線をループアンテナとして
扱い，ループ上の電流は平坦なものと仮定して，これを
集中定数回路の計算手法で求める方法である。この微小
ループアンテナ近似法による計算は最も簡単であるが，
ループの寸法が電磁波の波長に比べて無視できない条件
では精度が極めて低下することから，現実にはほとんど
用いられていない。
負荷回路との間を接続する配線をループアンテナとして
扱い，ループ上の電流は平坦なものと仮定して，これを
集中定数回路の計算手法で求める方法である。この微小
ループアンテナ近似法による計算は最も簡単であるが，
ループの寸法が電磁波の波長に比べて無視できない条件
では精度が極めて低下することから，現実にはほとんど
用いられていない。
【０００７】また，分布定数線路近似法は，１次元の構
造物として近似できる物体に対して，分布定数線路の方
程式を適用することで電流を求める方法である。計算は
比較的簡単であり，計算時間および記憶容量も解析要素
数にほぼ比例して増加するだけで，線路の反射や共振等
の現象も含めた解析ができることから，１次元の近似が
成立する物体については，高速・高精度の解析ができ
る。この分布定数線路近似法による計算は，１次元の構
造物として近似できる物体については，高速・高精度の
解析ができるものの，近似できない物体については解析
できないという問題点がある。
造物として近似できる物体に対して，分布定数線路の方
程式を適用することで電流を求める方法である。計算は
比較的簡単であり，計算時間および記憶容量も解析要素
数にほぼ比例して増加するだけで，線路の反射や共振等
の現象も含めた解析ができることから，１次元の近似が
成立する物体については，高速・高精度の解析ができ
る。この分布定数線路近似法による計算は，１次元の構
造物として近似できる物体については，高速・高精度の
解析ができるものの，近似できない物体については解析
できないという問題点がある。
【０００８】一方，モーメント法は，マックスウェルの
電磁波動方程式から導かれる積分方程式の解法の１つで
あり，３次元の任意形状物体を扱うことができる。具体
的には，物体を小さな要素に分割して電流の計算を行う
ものである。このように，モーメント法が３次元の任意
形状物体を扱えることから，電磁界強度算出装置では，
モーメント法を使って，電気回路装置の放射する電磁界
強度を算出するという構成が有力である。
電磁波動方程式から導かれる積分方程式の解法の１つで
あり，３次元の任意形状物体を扱うことができる。具体
的には，物体を小さな要素に分割して電流の計算を行う
ものである。このように，モーメント法が３次元の任意
形状物体を扱えることから，電磁界強度算出装置では，
モーメント法を使って，電気回路装置の放射する電磁界
強度を算出するという構成が有力である。
【０００９】このモーメント法を用いる方式では，金属
対象物を扱うときには，金属部分を解析対象としてメッ
シュ化し，分割した金属間の相互インピーダンスＺijを
求め，この相互インピーダンスＺijと，波源Ｖi と，分
割した金属に流れる電流Ｉiとの間に成立するモーメン
ト法の連立方程式 〔Ｚij〕〔Ｉi 〕＝〔Ｖi 〕 を解いて電流Ｉi を求め，この結果から放射される電磁
界強度を算出するという方法を採る。ここで，「〔
〕」はマトリクスを表している。
対象物を扱うときには，金属部分を解析対象としてメッ
シュ化し，分割した金属間の相互インピーダンスＺijを
求め，この相互インピーダンスＺijと，波源Ｖi と，分
割した金属に流れる電流Ｉiとの間に成立するモーメン
ト法の連立方程式 〔Ｚij〕〔Ｉi 〕＝〔Ｖi 〕 を解いて電流Ｉi を求め，この結果から放射される電磁
界強度を算出するという方法を採る。ここで，「〔
〕」はマトリクスを表している。
【００１０】なお，モーメント法についての参考文献と
しては，以下のものがある。 〔参考文献１〕H.N.Wang, J.H.Richmond and M.C.Gilre
ath ：“Sinusoidal reaction formulation for radiat
ion and scattering from conducting surface”IEEE T
RANSACTIONS ON ANTENNAS AND PROPAGATION vol.AP-23,
1975． 以上のように電磁界強度解析装置において，モーメント
法により算出した電流は，公知の式から電磁界強度を算
出するためのものである。従来，このような電流は，電
磁界強度を求めるための途中結果であり，これを可視化
して表示することは考えられていなかった。また，モー
メント法では，電流はメッシュ化した要素（これをパッ
チという）の辺ごとに計算されるため，要素が粗かった
り，または要素の並びや大きさがバラバラであると，算
出した電流を可視化して表示しても活用することができ
ないという問題があった。
しては，以下のものがある。 〔参考文献１〕H.N.Wang, J.H.Richmond and M.C.Gilre
ath ：“Sinusoidal reaction formulation for radiat
ion and scattering from conducting surface”IEEE T
RANSACTIONS ON ANTENNAS AND PROPAGATION vol.AP-23,
1975． 以上のように電磁界強度解析装置において，モーメント
法により算出した電流は，公知の式から電磁界強度を算
出するためのものである。従来，このような電流は，電
磁界強度を求めるための途中結果であり，これを可視化
して表示することは考えられていなかった。また，モー
メント法では，電流はメッシュ化した要素（これをパッ
チという）の辺ごとに計算されるため，要素が粗かった
り，または要素の並びや大きさがバラバラであると，算
出した電流を可視化して表示しても活用することができ
ないという問題があった。
【００１１】
【発明が解決しようとする課題】電気回路装置から放射
される電磁界強度の精度を上げるためには，その電磁界
発生の原因となる電流を把握できるようにすることが重
要である。また，これらの電流を詳しく表示すること
は，電磁波の放射メカニズムを把握する上で，非常に重
要である。さらに，電流の流れる様子を適切に表示すれ
ば，一度解析した電気回路装置またはこれから解析する
似たような形状の電気回路装置において，メッシュ分割
を行う際に，電流の強い箇所を密に，電流の弱い箇所を
粗に調整するための有力な情報となる。
される電磁界強度の精度を上げるためには，その電磁界
発生の原因となる電流を把握できるようにすることが重
要である。また，これらの電流を詳しく表示すること
は，電磁波の放射メカニズムを把握する上で，非常に重
要である。さらに，電流の流れる様子を適切に表示すれ
ば，一度解析した電気回路装置またはこれから解析する
似たような形状の電気回路装置において，メッシュ分割
を行う際に，電流の強い箇所を密に，電流の弱い箇所を
粗に調整するための有力な情報となる。
【００１２】本発明は上記問題点の解決を図り，電気回
路装置の要素上の任意の点における電流を求め，電流の
流れる様子を，要素の形状によらずに見やすく，わかり
やすい態様で表示する手段を提供することを目的とす
る。
路装置の要素上の任意の点における電流を求め，電流の
流れる様子を，要素の形状によらずに見やすく，わかり
やすい態様で表示する手段を提供することを目的とす
る。
【００１３】
【課題を解決するための手段】図１は本発明の原理構成
図である。図中，１０はＣＰＵおよびメモリ等からなる
処理装置により，モーメント法もしくは分布定数線路近
似法に基づいて解析対象となる電気回路装置の放射する
電磁界強度を算出する電磁界強度算出装置である。
図である。図中，１０はＣＰＵおよびメモリ等からなる
処理装置により，モーメント法もしくは分布定数線路近
似法に基づいて解析対象となる電気回路装置の放射する
電磁界強度を算出する電磁界強度算出装置である。
【００１４】データ入力部１１は，解析対象となる電気
回路装置の構造データ等の電気回路装置データ２２を入
力する入力手段である。電磁界強度計算部１２は，電磁
界強度の解析対象を微小要素であるパッチに分割し，モ
ーメント法に基づいて各要素に流れる電流を算出して電
磁界強度を計算する処理手段である。なお，電気回路装
置において分布定数線路近似法を適用できる部分には，
分布定数線路近似法を用いて電流を算出することもあ
る。
回路装置の構造データ等の電気回路装置データ２２を入
力する入力手段である。電磁界強度計算部１２は，電磁
界強度の解析対象を微小要素であるパッチに分割し，モ
ーメント法に基づいて各要素に流れる電流を算出して電
磁界強度を計算する処理手段である。なお，電気回路装
置において分布定数線路近似法を適用できる部分には，
分布定数線路近似法を用いて電流を算出することもあ
る。
【００１５】解析情報出力部１３は，電磁界強度計算部
１２で算出した電磁界強度を指定された形式で図化して
表示する手段である。本発明は，特に電磁界強度だけで
はなく，電磁界強度計算部１２で計算した電気回路装置
を流れる電流を表示する。そのため，以下の手段を有す
る。
１２で算出した電磁界強度を指定された形式で図化して
表示する手段である。本発明は，特に電磁界強度だけで
はなく，電磁界強度計算部１２で計算した電気回路装置
を流れる電流を表示する。そのため，以下の手段を有す
る。
【００１６】本発明は，第１の特徴として，解析対象の
電気回路装置モデルにおける平面を，指定された間隔ま
たは所定の間隔の格子状の線で区分する格子設定処理手
段１４と，モーメント法により算出された電気回路装置
の要素に流れる電流に基づいて，前記格子状の線の交点
または格子内の点の電流を計算する格子点電流計算手段
１５と，この格子点電流計算手段１５によって計算した
電流を可視化して表示する電流表示処理手段１６とを備
える。特に，格子設定処理手段１４が，電気回路装置モ
デルにおける任意平面を，少なくともその任意平面を包
含する矩形エリアで囲み，その矩形エリア内を前記モー
メント法での計算単位となる要素の大きさとは別の指定
された間隔または所定の間隔の格子状の線で区分するこ
とにより，任意平面を格子状の線で区分された領域の集
合体として表すように構成されていることを特徴とす
る。
電気回路装置モデルにおける平面を，指定された間隔ま
たは所定の間隔の格子状の線で区分する格子設定処理手
段１４と，モーメント法により算出された電気回路装置
の要素に流れる電流に基づいて，前記格子状の線の交点
または格子内の点の電流を計算する格子点電流計算手段
１５と，この格子点電流計算手段１５によって計算した
電流を可視化して表示する電流表示処理手段１６とを備
える。特に，格子設定処理手段１４が，電気回路装置モ
デルにおける任意平面を，少なくともその任意平面を包
含する矩形エリアで囲み，その矩形エリア内を前記モー
メント法での計算単位となる要素の大きさとは別の指定
された間隔または所定の間隔の格子状の線で区分するこ
とにより，任意平面を格子状の線で区分された領域の集
合体として表すように構成されていることを特徴とす
る。
【００１７】また，第１の特徴において，電流表示処理
手段１６が格子点電流計算手段１５により計算した電流
を直交座標系にベクトル分解し，その結果を立体的に分
布図として隠面消去した表示面で表示する電流分布図表
示処理手段１７を持つことを，第２の特徴とする。電流
分布図は，ある平面に対して，ＸＹＺの３次元直交座標
系を考えて，そのＺ成分を電流の大きさとして，３次元
グラフにしたものである。
手段１６が格子点電流計算手段１５により計算した電流
を直交座標系にベクトル分解し，その結果を立体的に分
布図として隠面消去した表示面で表示する電流分布図表
示処理手段１７を持つことを，第２の特徴とする。電流
分布図は，ある平面に対して，ＸＹＺの３次元直交座標
系を考えて，そのＺ成分を電流の大きさとして，３次元
グラフにしたものである。
【００１８】また，第１の特徴において，電流表示処理
手段１６が格子点電流計算手段１５により計算した電流
を，電気回路装置のモデル図における該当する面上に，
電流の方向を矢印の向き，電流の大きさを矢印の長さで
表した矢印図形で，前記各格子に対応して表示する電流
ベクトル図表示処理手段１８を持つことを，第３の特徴
とする。電流ベクトル図は，モデルのすべての面に対し
て計算した格子点の電流の大きさと向きをモデルの表面
上に表したものである。
手段１６が格子点電流計算手段１５により計算した電流
を，電気回路装置のモデル図における該当する面上に，
電流の方向を矢印の向き，電流の大きさを矢印の長さで
表した矢印図形で，前記各格子に対応して表示する電流
ベクトル図表示処理手段１８を持つことを，第３の特徴
とする。電流ベクトル図は，モデルのすべての面に対し
て計算した格子点の電流の大きさと向きをモデルの表面
上に表したものである。
【００１９】
【００２０】さらに，第１の特徴において，格子設定処
理手段１４が電気回路装置モデルにおいて三角形または
四角形で構成される任意平面を，格子状の線で区分した
領域のみで表すように構成されていることを，第４の特
徴とする。
理手段１４が電気回路装置モデルにおいて三角形または
四角形で構成される任意平面を，格子状の線で区分した
領域のみで表すように構成されていることを，第４の特
徴とする。
【００２１】また，第１の特徴において，電気回路装置
のモデル図内の１点または１要素を指示する入力を受け
付け，その指示された１点または１要素を含む平面を，
電流表示処理手段１６における電流の表示対象平面とす
る平面選択処理手段２０を備えたことを，第５の特徴と
する。
のモデル図内の１点または１要素を指示する入力を受け
付け，その指示された１点または１要素を含む平面を，
電流表示処理手段１６における電流の表示対象平面とす
る平面選択処理手段２０を備えたことを，第５の特徴と
する。
【００２２】また，第１の特徴において，電流表示処理
手段１６が格子点電流計算手段１５により計算した電流
の結果における実数部と虚数部とを，外部からの選択指
示によって別に表示する手段を持つことを，第６の特徴
とする。
手段１６が格子点電流計算手段１５により計算した電流
の結果における実数部と虚数部とを，外部からの選択指
示によって別に表示する手段を持つことを，第６の特徴
とする。
【００２３】また，第１の特徴において，電流表示処理
手段１６が格子点電流計算手段１５により計算した電流
の結果を絶対値と位相とで表示する手段を持つことを，
第７の特徴とする。
手段１６が格子点電流計算手段１５により計算した電流
の結果を絶対値と位相とで表示する手段を持つことを，
第７の特徴とする。
【００２４】さらに，第２の特徴において，電流表示処
理手段１６が電流の表示対象としている平面が，電気回
路装置の３次元モデル図におけるどの平面であるかを示
すモデル図を表示するモデル図表示処理手段２１を備え
ることを，第８の特徴とする。
理手段１６が電流の表示対象としている平面が，電気回
路装置の３次元モデル図におけるどの平面であるかを示
すモデル図を表示するモデル図表示処理手段２１を備え
ることを，第８の特徴とする。
【００２５】さらに，第３の特徴において，電流表示処
理手段１６の電流ベクトル図表示処理手段１８は，電流
の大きさに応じて矢印の表示態様を変えることにより，
複数のスケールを使い分けるように構成されていること
を，第９の特徴とする。
理手段１６の電流ベクトル図表示処理手段１８は，電流
の大きさに応じて矢印の表示態様を変えることにより，
複数のスケールを使い分けるように構成されていること
を，第９の特徴とする。
【００２６】また，第１の特徴において，電流表示処理
手段１６が格子点電流計算手段１５により計算した電流
の結果について，表示対象として指定された複数個の各
点の位置を横軸とし，縦軸を各点における電流値とした
ＸＹグラフ形式で電流を表示する電流ＸＹグラフ図表示
処理手段１９を持つことを，第１０の特徴とする。電流
ＸＹグラフ図は，任意の点をいくつか選択し，横軸に置
き，その点での電流値をＹ軸で表示した図である。特
に，ケーブル上の電流値を比較する場合などに，有用で
ある。
手段１６が格子点電流計算手段１５により計算した電流
の結果について，表示対象として指定された複数個の各
点の位置を横軸とし，縦軸を各点における電流値とした
ＸＹグラフ形式で電流を表示する電流ＸＹグラフ図表示
処理手段１９を持つことを，第１０の特徴とする。電流
ＸＹグラフ図は，任意の点をいくつか選択し，横軸に置
き，その点での電流値をＹ軸で表示した図である。特
に，ケーブル上の電流値を比較する場合などに，有用で
ある。
【００２７】さらに，第１０の特徴において，電流ＸＹ
グラフ図表示処理手段１９によりＸＹグラフ形式で表示
する点列の範囲を，その範囲の始点と終点とを指定する
入力装置２３からの入力で定めるように構成されている
ことを，第１１の特徴とする。
グラフ図表示処理手段１９によりＸＹグラフ形式で表示
する点列の範囲を，その範囲の始点と終点とを指定する
入力装置２３からの入力で定めるように構成されている
ことを，第１１の特徴とする。
【００２８】
【００２９】
【作用】図２は，本発明における格子設定処理説明図で
ある。解析対象の電気回路装置について，モーメント法
により電流を計算する場合には，解析対象を微小要素で
あるパッチに分割し，各パッチ（要素）間の相互インピ
ーダンスを計算して，それらの相互インピーダンスをも
とに連立方程式を解くことにより，各要素に流れる電流
を算出する。
ある。解析対象の電気回路装置について，モーメント法
により電流を計算する場合には，解析対象を微小要素で
あるパッチに分割し，各パッチ（要素）間の相互インピ
ーダンスを計算して，それらの相互インピーダンスをも
とに連立方程式を解くことにより，各要素に流れる電流
を算出する。
【００３０】図２（Ａ）に示すＥ１〜Ｅ９は，電気回路
装置の任意平面においてモーメント法での計算単位とな
る要素を表す。電流は，これらの要素Ｅ１〜Ｅ９の辺ご
とに計算される。そのため，これらの要素Ｅ１〜Ｅ９が
粗かったり，サイズがバラバラであったりすると，計算
された電流を表示しても利用者にとって活用することは
難しい。
装置の任意平面においてモーメント法での計算単位とな
る要素を表す。電流は，これらの要素Ｅ１〜Ｅ９の辺ご
とに計算される。そのため，これらの要素Ｅ１〜Ｅ９が
粗かったり，サイズがバラバラであったりすると，計算
された電流を表示しても利用者にとって活用することは
難しい。
【００３１】本発明の第１の特徴によれば，格子設定処
理手段１４により，図２（Ａ）に示す電気回路装置の平
面上に，図２（Ｂ）に示すように，指定された間隔また
は所定の間隔の格子を設定する。そして，格子点電流計
算手段１５により，先に電磁界強度計算部１２において
モーメント法により算出した要素Ｅ１〜Ｅ９の辺に流れ
る電流に基づいて，格子点Ｐ１〜Ｐ１１の電流を計算
し，電流表示処理手段１６により，この電流を表示装置
２４に可視化して表示する。
理手段１４により，図２（Ａ）に示す電気回路装置の平
面上に，図２（Ｂ）に示すように，指定された間隔また
は所定の間隔の格子を設定する。そして，格子点電流計
算手段１５により，先に電磁界強度計算部１２において
モーメント法により算出した要素Ｅ１〜Ｅ９の辺に流れ
る電流に基づいて，格子点Ｐ１〜Ｐ１１の電流を計算
し，電流表示処理手段１６により，この電流を表示装置
２４に可視化して表示する。
【００３２】これにより，装置の各部に流れる電流が同
じスケールで統一して可視化されることになり，電磁界
強度に影響する電流の把握が容易になる。したがって，
電磁波対策をする上で，重要な手掛かりを得ることがで
きる。また，この第１の特徴によれば，電気回路装置モ
デルにおいて電流を表示する平面が不規則な形状であっ
ても，図２（Ｂ）に示すように，それを包含する矩形エ
リアで囲み，その矩形エリア内を格子状の線で区分する
ので，規則的に電流を可視化することができる。
じスケールで統一して可視化されることになり，電磁界
強度に影響する電流の把握が容易になる。したがって，
電磁波対策をする上で，重要な手掛かりを得ることがで
きる。また，この第１の特徴によれば，電気回路装置モ
デルにおいて電流を表示する平面が不規則な形状であっ
ても，図２（Ｂ）に示すように，それを包含する矩形エ
リアで囲み，その矩形エリア内を格子状の線で区分する
ので，規則的に電流を可視化することができる。
【００３３】図３は，本発明による電流分布図の表示例
を示す。本発明の第２の特徴によれば，格子点電流計算
手段１５によって計算した電流を直交座標系にベクトル
分解し，その結果のＸ方向またはＹ方向の電流の大きさ
を，例えば図３に示すように立体的に電流分布図として
表示する。したがって，ベクトル分解した方向における
電磁界強度に対する電流の影響を容易に把握できるよう
になる。
を示す。本発明の第２の特徴によれば，格子点電流計算
手段１５によって計算した電流を直交座標系にベクトル
分解し，その結果のＸ方向またはＹ方向の電流の大きさ
を，例えば図３に示すように立体的に電流分布図として
表示する。したがって，ベクトル分解した方向における
電磁界強度に対する電流の影響を容易に把握できるよう
になる。
【００３４】図４は，本発明による電流ベクトル図の表
示例を示す。図４（Ａ）では，前面の電流ベクトル図の
みを表示し，他の面は省略している。（Ｂ）は部分拡大
図である。太い実線は電気回路装置のモデル形状，細い
実線はモーメント法におけるメッシュ化された要素，点
線は格子設定処理手段１４によって設定された格子を表
している。
示例を示す。図４（Ａ）では，前面の電流ベクトル図の
みを表示し，他の面は省略している。（Ｂ）は部分拡大
図である。太い実線は電気回路装置のモデル形状，細い
実線はモーメント法におけるメッシュ化された要素，点
線は格子設定処理手段１４によって設定された格子を表
している。
【００３５】本発明の第３の特徴によれば，図４に示す
ように，格子点電流計算手段１５によって計算された電
流の結果が，モデル図の面上に矢印の長さが大きさ，矢
印の向きが流れる方向となるようにベクトル表示される
ので，装置のどの部分に電流が集中しているかを一目で
把握することができるようになり，電磁波対策をする上
で重要な手掛かりを得ることが可能になる。
ように，格子点電流計算手段１５によって計算された電
流の結果が，モデル図の面上に矢印の長さが大きさ，矢
印の向きが流れる方向となるようにベクトル表示される
ので，装置のどの部分に電流が集中しているかを一目で
把握することができるようになり，電磁波対策をする上
で重要な手掛かりを得ることが可能になる。
【００３６】
【００３７】また，本発明の第４の特徴によれば，電流
を表示する任意平面が，例えば三角形であっても，格子
状の線で区分した領域のみで表すことにより，規則的に
電流を可視化することができる。
を表示する任意平面が，例えば三角形であっても，格子
状の線で区分した領域のみで表すことにより，規則的に
電流を可視化することができる。
【００３８】図５は，本発明によるモデル図の表示例を
示す。本発明の第５の特徴によれば，例えば，図５
（Ａ）に示すような，前面に開口部がある電気回路装置
の筐体に流れる電流を表示する場合，図５（Ｂ）に示す
ようなモデル図の表示において，マウスカーソルＣ等に
より１点または１要素を指定し，電流を表示させたい平
面を選択することができる。また，平面上のパッチ番号
を指定することによって，そのパッチ番号の要素が含ま
れる平面を選択することもできる。したがって，平面が
入り組んで複雑な場合にも，任意の平面だけについて，
簡単に電流の流れる様子を表示させることができる。
示す。本発明の第５の特徴によれば，例えば，図５
（Ａ）に示すような，前面に開口部がある電気回路装置
の筐体に流れる電流を表示する場合，図５（Ｂ）に示す
ようなモデル図の表示において，マウスカーソルＣ等に
より１点または１要素を指定し，電流を表示させたい平
面を選択することができる。また，平面上のパッチ番号
を指定することによって，そのパッチ番号の要素が含ま
れる平面を選択することもできる。したがって，平面が
入り組んで複雑な場合にも，任意の平面だけについて，
簡単に電流の流れる様子を表示させることができる。
【００３９】本発明の第６の特徴によれば，格子点電流
計算手段１５により計算した電流の結果について，実数
部と虚数部とを分けて可視化することができるので，電
磁波対策に有用な情報を得ることができる。
計算手段１５により計算した電流の結果について，実数
部と虚数部とを分けて可視化することができるので，電
磁波対策に有用な情報を得ることができる。
【００４０】同様に，本発明の第７の特徴によれば，格
子点電流計算手段１５により計算した電流の結果から，
例えば電流ＸＹグラフ図上に電流の大きさの絶対値と位
相とを重ねて表示する。これにより電磁波対策に有用な
情報を得ることができる。
子点電流計算手段１５により計算した電流の結果から，
例えば電流ＸＹグラフ図上に電流の大きさの絶対値と位
相とを重ねて表示する。これにより電磁波対策に有用な
情報を得ることができる。
【００４１】本発明の第８の特徴によれば，電流分布図
表示処理手段１７により表示している平面が，３次元モ
デル図におけるどの平面であるかを，モデル図表示処理
手段２１によって表示するので，装置のどの部分が現在
表示されているかを一目で把握することができる。な
お，この面の表示は，例えば図５（Ｂ）に示すようなモ
デル図において，該当する面を色分けや輝度を強調する
ことによって行う。
表示処理手段１７により表示している平面が，３次元モ
デル図におけるどの平面であるかを，モデル図表示処理
手段２１によって表示するので，装置のどの部分が現在
表示されているかを一目で把握することができる。な
お，この面の表示は，例えば図５（Ｂ）に示すようなモ
デル図において，該当する面を色分けや輝度を強調する
ことによって行う。
【００４２】図６は，本発明による矢印の色分け表示の
例を説明する図である。例えば図４に示すような電流ベ
クトル図表示において，矢印の長さを電流の大きさに比
例させて表示すると，大きな電流の場合には矢印が長く
なり過ぎて見にくくなり，また小さな電流の場合には矢
印が短くなり過ぎて見えなくなるということが生じる。
例を説明する図である。例えば図４に示すような電流ベ
クトル図表示において，矢印の長さを電流の大きさに比
例させて表示すると，大きな電流の場合には矢印が長く
なり過ぎて見にくくなり，また小さな電流の場合には矢
印が短くなり過ぎて見えなくなるということが生じる。
【００４３】本発明の第９の特徴によれば，例えば図６
に示すように，電流値が基準のスケールに対して，１０
０倍以上の場合には矢印（ベクトル）の色を赤色，１０
倍以上の場合にはオレンジ色，１／１０〜１０倍の場合
には黄色，１／１０の場合には緑色，１／１００以下の
場合には青色というように，色を変えて表示するので，
電流ベクトル図を小さな電流から大きな電流まで一度に
見やすく表示することができるようになる。色分けに代
えて，矢印の太さを変えるなどして区別できるようにし
てもよい。
に示すように，電流値が基準のスケールに対して，１０
０倍以上の場合には矢印（ベクトル）の色を赤色，１０
倍以上の場合にはオレンジ色，１／１０〜１０倍の場合
には黄色，１／１０の場合には緑色，１／１００以下の
場合には青色というように，色を変えて表示するので，
電流ベクトル図を小さな電流から大きな電流まで一度に
見やすく表示することができるようになる。色分けに代
えて，矢印の太さを変えるなどして区別できるようにし
てもよい。
【００４４】図７は，本発明による電流ＸＹグラフ図の
表示例を示す。本発明の第１０の特徴によれば，例えば
図７に示すように，電流ＸＹグラフ形式で，電流値を実
数・虚数または絶対値・位相で重ね書き表示することが
できるので，特定の線上における電流の変化を一目で把
握することができるようになり，電磁波対策の有力な手
掛かりを得ることができる。
表示例を示す。本発明の第１０の特徴によれば，例えば
図７に示すように，電流ＸＹグラフ形式で，電流値を実
数・虚数または絶対値・位相で重ね書き表示することが
できるので，特定の線上における電流の変化を一目で把
握することができるようになり，電磁波対策の有力な手
掛かりを得ることができる。
【００４５】さらに，本発明の第１１の特徴によれば，
例えば図７に示すような電流ＸＹグラフ図を表示するに
あたって，表示対象とする格子点の始点と終点とを，格
子点番号の数値入力またはモデル図上のマウス等による
２点入力することにより，表示範囲を指定することがで
きるので，簡易で迅速な入力が可能である。
例えば図７に示すような電流ＸＹグラフ図を表示するに
あたって，表示対象とする格子点の始点と終点とを，格
子点番号の数値入力またはモデル図上のマウス等による
２点入力することにより，表示範囲を指定することがで
きるので，簡易で迅速な入力が可能である。
【００４６】
【００４７】
【実施例】図８は本発明の実施例の処理フローの概要図
である。まず，ステップＳ１０により，ＣＡＤ等によっ
て生成された解析対象である電気回路装置データ２２を
入力し，モーメント法を適用する部分について，解析対
象をワイヤーやサーフェイスパッチに分割し，電流の計
算に必要なデータを構造体，配列として設定する。
である。まず，ステップＳ１０により，ＣＡＤ等によっ
て生成された解析対象である電気回路装置データ２２を
入力し，モーメント法を適用する部分について，解析対
象をワイヤーやサーフェイスパッチに分割し，電流の計
算に必要なデータを構造体，配列として設定する。
【００４８】次に，あらかじめ定義された周波数につい
て以下の処理を繰り返す。そのため，まずステップＳ１
１で，処理済の周波数をカウントすることで，登録され
ている全周波数についての処理が終了したか否かを判断
する。終了した場合にはステップＳ１９へ移り，未処理
の周波数がある場合には，次に処理する周波数を未処理
の中から１つ選択して，続くステップＳ１２で，この選
択した周波数の波長λを計算する。
て以下の処理を繰り返す。そのため，まずステップＳ１
１で，処理済の周波数をカウントすることで，登録され
ている全周波数についての処理が終了したか否かを判断
する。終了した場合にはステップＳ１９へ移り，未処理
の周波数がある場合には，次に処理する周波数を未処理
の中から１つ選択して，続くステップＳ１２で，この選
択した周波数の波長λを計算する。
【００４９】続いて，順次，モーメント法の連立方程式
で用いる相互インピーダンスＺijを算出するために，ス
テップＳ１３によりｍ個ある金属要素の中から要素ｉ
（ｉ＝１〜ｍ）と要素ｊ（ｊ＝１〜ｍ）の組み合わせを
生成する。そして，次のステップＳ１４で，組み合わせ
た要素間の相互インピーダンスＺijを計算する。すべて
の組み合わせの計算が終了するまで，ステップＳ１３，
Ｓ１４を繰り返し，ステップＳ１５により，すべて要素
の組み合わせについて相互インピーダンスの計算が終了
したことを確認したならば，次のステップＳ１６へ進
む。
で用いる相互インピーダンスＺijを算出するために，ス
テップＳ１３によりｍ個ある金属要素の中から要素ｉ
（ｉ＝１〜ｍ）と要素ｊ（ｊ＝１〜ｍ）の組み合わせを
生成する。そして，次のステップＳ１４で，組み合わせ
た要素間の相互インピーダンスＺijを計算する。すべて
の組み合わせの計算が終了するまで，ステップＳ１３，
Ｓ１４を繰り返し，ステップＳ１５により，すべて要素
の組み合わせについて相互インピーダンスの計算が終了
したことを確認したならば，次のステップＳ１６へ進
む。
【００５０】ステップＳ１６では，算出した相互インピ
ーダンスＺijと，入力データとして与えられた波源Ｖi
とを使って，金属要素に流れる電流Ｉi を未知数とする
モーメント法の連立方程式 〔Ｚij〕〔Ｉi 〕＝〔Ｖi 〕 を導出し，これを解くことで，金属要素に流れる電流Ｉ
i を求める。求めた電流Ｉi は，電流ファイル４０に保
存する。
ーダンスＺijと，入力データとして与えられた波源Ｖi
とを使って，金属要素に流れる電流Ｉi を未知数とする
モーメント法の連立方程式 〔Ｚij〕〔Ｉi 〕＝〔Ｖi 〕 を導出し，これを解くことで，金属要素に流れる電流Ｉ
i を求める。求めた電流Ｉi は，電流ファイル４０に保
存する。
【００５１】続くステップＳ１７で，処理済の観測点を
カウントすることで，登録されている全観測点について
の処理が終了したか否かを判断する。終了していない場
合には，ステップＳ１８に進んで，算出した電流Ｉi が
観測点にもたらす電磁界強度を規定の算出式に従って算
出する。その算出結果を電界・磁界ファイル４１に格納
してからステップＳ１７に戻り，すべての観測点につい
て電界・磁界を計算する。全観測点についての処理が終
了した場合にはステップＳ１１に戻り，次の周波数につ
いて同様に処理を繰り返す。
カウントすることで，登録されている全観測点について
の処理が終了したか否かを判断する。終了していない場
合には，ステップＳ１８に進んで，算出した電流Ｉi が
観測点にもたらす電磁界強度を規定の算出式に従って算
出する。その算出結果を電界・磁界ファイル４１に格納
してからステップＳ１７に戻り，すべての観測点につい
て電界・磁界を計算する。全観測点についての処理が終
了した場合にはステップＳ１１に戻り，次の周波数につ
いて同様に処理を繰り返す。
【００５２】上記ステップＳ１４による相互インピーダ
ンスＺijの計算，ステップＳ１６による電流計算および
ステップＳ１８による電界，磁界の計算法については，
従来から知られている方法を用いることができるので，
ここでの詳細な説明は省略する。
ンスＺijの計算，ステップＳ１６による電流計算および
ステップＳ１８による電界，磁界の計算法については，
従来から知られている方法を用いることができるので，
ここでの詳細な説明は省略する。
【００５３】なお，説明を簡単にするために，電磁波の
解析についてモーメント法を用いる場合の例だけを説明
したが，本発明は，電気回路装置のすべての部分にモー
メント法を用いる場合だけではなく，一部にのみモーメ
ント法を用いる場合にも適用できることは言うまでもな
い。例えば１次元の構造物として近似できる部分に対し
ては，分布定数線路近似法を適用し，分布定数線路の方
程式によって各部の電流等を算出することもできる。
解析についてモーメント法を用いる場合の例だけを説明
したが，本発明は，電気回路装置のすべての部分にモー
メント法を用いる場合だけではなく，一部にのみモーメ
ント法を用いる場合にも適用できることは言うまでもな
い。例えば１次元の構造物として近似できる部分に対し
ては，分布定数線路近似法を適用し，分布定数線路の方
程式によって各部の電流等を算出することもできる。
【００５４】ステップＳ１１の判定により，全周波数に
ついての計算が終了したことを確認したならば，次にス
テップＳ１９により，電流の表示要求があるかどうかを
判定する。電流の表示要求がない場合には，その他の要
求された処理を実行する。その他の処理として，例えば
電界・磁界ファイル４１に格納された計算結果を図化し
て表示する処理などがあるが，本発明の要旨には直接関
係しないため，ここでの詳しい説明は省略する。
ついての計算が終了したことを確認したならば，次にス
テップＳ１９により，電流の表示要求があるかどうかを
判定する。電流の表示要求がない場合には，その他の要
求された処理を実行する。その他の処理として，例えば
電界・磁界ファイル４１に格納された計算結果を図化し
て表示する処理などがあるが，本発明の要旨には直接関
係しないため，ここでの詳しい説明は省略する。
【００５５】電流表示要求があった場合には，ステップ
Ｓ２０によって，金属要素ごとに求められた電流値が格
納された電流ファイル４０を参照し，本発明に係る電流
表示処理を行う。その後，ステップＳ２１によって終了
要求があったかどうかを判定し，終了要求があったなら
ば，処理を終了する。終了要求がないならば，ステップ
Ｓ１９に戻り，同様に処理を繰り返す。
Ｓ２０によって，金属要素ごとに求められた電流値が格
納された電流ファイル４０を参照し，本発明に係る電流
表示処理を行う。その後，ステップＳ２１によって終了
要求があったかどうかを判定し，終了要求があったなら
ば，処理を終了する。終了要求がないならば，ステップ
Ｓ１９に戻り，同様に処理を繰り返す。
【００５６】図９は，本発明の実施例における電流表示
処理手段１６の処理概要説明図である。前述した図８の
ステップＳ２０では，電流表示処理手段１６により図９
に示す処理を行う。
処理手段１６の処理概要説明図である。前述した図８の
ステップＳ２０では，電流表示処理手段１６により図９
に示す処理を行う。
【００５７】まず，図９のステップＳ３０により，電流
ファイル４０から各金属要素に流れる電流値を入力し，
指定された間隔または所定の間隔で設定された各格子の
交点または格子内の点の電流を計算する。次に，ステッ
プＳ３１によって，電流の表示方式を利用者からの指示
によって選択する。
ファイル４０から各金属要素に流れる電流値を入力し，
指定された間隔または所定の間隔で設定された各格子の
交点または格子内の点の電流を計算する。次に，ステッ
プＳ３１によって，電流の表示方式を利用者からの指示
によって選択する。
【００５８】電流ベクトル図の表示を要求されたなら
ば，ステップＳ３２によって表示すべきスケールを選択
し，図４に示すような電流ベクトル図を表示する。電流
分布図の表示を要求されたならば，ステップＳ３３によ
って表示すべきスケールを選択し，また，ステップＳ３
４により，電流を絶対値によって表示するか，実数また
は虚数によって表示するかを選択し，さらにステップＳ
３５により，表示平面を選択して，図３に示すような電
流分布図を表示する。
ば，ステップＳ３２によって表示すべきスケールを選択
し，図４に示すような電流ベクトル図を表示する。電流
分布図の表示を要求されたならば，ステップＳ３３によ
って表示すべきスケールを選択し，また，ステップＳ３
４により，電流を絶対値によって表示するか，実数また
は虚数によって表示するかを選択し，さらにステップＳ
３５により，表示平面を選択して，図３に示すような電
流分布図を表示する。
【００５９】電流ＸＹグラフ図の表示を要求されたなら
ば，ステップＳ３６によって表示すべきスケールを選択
し，また，ステップＳ３７により，電流を絶対値によっ
て表示するか，実数または虚数によって表示するかを選
択し，さらにステップＳ３８により，表示平面または表
示すべきワイヤを選択して，図７に示すような電流ＸＹ
グラフ図を表示する。
ば，ステップＳ３６によって表示すべきスケールを選択
し，また，ステップＳ３７により，電流を絶対値によっ
て表示するか，実数または虚数によって表示するかを選
択し，さらにステップＳ３８により，表示平面または表
示すべきワイヤを選択して，図７に示すような電流ＸＹ
グラフ図を表示する。
【００６０】図１０は，本発明の実施例における電流表
示用パラメータ入力画面の例を示す。解析対象の電気回
路装置を流れる電流の表示が要求された場合，図１０に
示すような電流表示用パラメータ入力画面をポップアッ
プ表示する。この画面において，周波数入力領域５０で
は，キーボードやマウス等を用いて，表示する周波数を
選択することができる。ドロウタイプ入力領域５１にお
いて，「Ｖｅｃｔｏｒ」をクリックすると，電流ベクト
ル図の表示要求となり，「Ｄｉｓｔｒｉｂｕｔｉｏｎ−
ｐｌａｎｅ」をクリックすると，電流分布図の表示要求
となり，「Ｄｉｓｔｒｉｂｕｔｉｏｎ−ＸＹ」をクリッ
クすると，電流ＸＹグラフ図の表示要求となる。
示用パラメータ入力画面の例を示す。解析対象の電気回
路装置を流れる電流の表示が要求された場合，図１０に
示すような電流表示用パラメータ入力画面をポップアッ
プ表示する。この画面において，周波数入力領域５０で
は，キーボードやマウス等を用いて，表示する周波数を
選択することができる。ドロウタイプ入力領域５１にお
いて，「Ｖｅｃｔｏｒ」をクリックすると，電流ベクト
ル図の表示要求となり，「Ｄｉｓｔｒｉｂｕｔｉｏｎ−
ｐｌａｎｅ」をクリックすると，電流分布図の表示要求
となり，「Ｄｉｓｔｒｉｂｕｔｉｏｎ−ＸＹ」をクリッ
クすると，電流ＸＹグラフ図の表示要求となる。
【００６１】実行ボタン５２のクリックによって，入力
されたパラメータに基づく電流の表示が行われる。終了
ボタン５３がクリックされると，電流の表示を実行しな
いで，処理を終了する。
されたパラメータに基づく電流の表示が行われる。終了
ボタン５３がクリックされると，電流の表示を実行しな
いで，処理を終了する。
【００６２】図１１は，本発明の実施例における電流ベ
クトル図表示画面の例を示す。図１０に示すドロウタイ
プ入力領域５１において「Ｖｅｃｔｏｒ」が選択され，
実行ボタン５２がクリックされると，例えば図１１に示
すような電流ベクトル図が表示される。この画面におけ
るモード指定領域６０において，「Ａｂｓｏｌｕｔ
ｅ」，「Ｒｅａｌ」，「Ｉｍａｇｉｎａｒｙ」の選択が
可能であり，「Ａｂｓｏｌｕｔｅ」が選択されると，電
流値を絶対値で表示する。「Ｒｅａｌ」が選択される
と，複素数で計算された電流値の実数部分の大きさに従
って電流を表示する。また，「Ｉｍａｇｉｎａｒｙ」が
選択されると，複素数で計算された電流値の虚数部分の
大きさに基づいて電流を表示する。
クトル図表示画面の例を示す。図１０に示すドロウタイ
プ入力領域５１において「Ｖｅｃｔｏｒ」が選択され，
実行ボタン５２がクリックされると，例えば図１１に示
すような電流ベクトル図が表示される。この画面におけ
るモード指定領域６０において，「Ａｂｓｏｌｕｔ
ｅ」，「Ｒｅａｌ」，「Ｉｍａｇｉｎａｒｙ」の選択が
可能であり，「Ａｂｓｏｌｕｔｅ」が選択されると，電
流値を絶対値で表示する。「Ｒｅａｌ」が選択される
と，複素数で計算された電流値の実数部分の大きさに従
って電流を表示する。また，「Ｉｍａｇｉｎａｒｙ」が
選択されると，複素数で計算された電流値の虚数部分の
大きさに基づいて電流を表示する。
【００６３】スケール入力領域６１では，電流値の最大
スケールの設定を行うことができる。スケール入力領域
６１におけるボタンをマウスの右でクリックすると，ス
ケールのメニューが表示される。
スケールの設定を行うことができる。スケール入力領域
６１におけるボタンをマウスの右でクリックすると，ス
ケールのメニューが表示される。
【００６４】実行ボタン６２のクリックによって，新た
に設定されたパラメータによる表示を実行する。ズーム
ボタン６３のクリックによって，拡大表示を行う。この
ボタンをクリックした後，矩形領域をマウスによってド
ラッグすると，その領域が拡大して表示される。リセッ
トボタン６４のクリックによって，画面がリセットされ
る。
に設定されたパラメータによる表示を実行する。ズーム
ボタン６３のクリックによって，拡大表示を行う。この
ボタンをクリックした後，矩形領域をマウスによってド
ラッグすると，その領域が拡大して表示される。リセッ
トボタン６４のクリックによって，画面がリセットされ
る。
【００６５】図１２は，本発明の実施例における電流分
布図表示画面の例を示す。図１０に示すドロウタイプ入
力領域５１において「Ｄｉｓｔｒｉｂｕｔｉｏｎ−ｐｌ
ａｎｅ」が選択され，実行ボタン５２がクリックされる
と，例えば図１２に示すような電流分布図が表示され
る。この画面におけるモード指定領域７０において，
「Ａｂｓｏｌｕｔｅ」，「Ｒｅａｌ」，「Ｉｍａｇｉｎ
ａｒｙ」の選択が可能であり，これらのいずれかの選択
によって，それぞれ電流の絶対値による表示，実数値に
よる表示，虚数値による表示が行われる。
布図表示画面の例を示す。図１０に示すドロウタイプ入
力領域５１において「Ｄｉｓｔｒｉｂｕｔｉｏｎ−ｐｌ
ａｎｅ」が選択され，実行ボタン５２がクリックされる
と，例えば図１２に示すような電流分布図が表示され
る。この画面におけるモード指定領域７０において，
「Ａｂｓｏｌｕｔｅ」，「Ｒｅａｌ」，「Ｉｍａｇｉｎ
ａｒｙ」の選択が可能であり，これらのいずれかの選択
によって，それぞれ電流の絶対値による表示，実数値に
よる表示，虚数値による表示が行われる。
【００６６】パッチ番号入力領域７１では，表示する平
面上の１つのパッチ番号を指定することができる。ここ
でパッチ番号を指定することにより，そのパッチ番号の
要素を含む平面が，電流分布図の表示対象平面となる。
面上の１つのパッチ番号を指定することができる。ここ
でパッチ番号を指定することにより，そのパッチ番号の
要素を含む平面が，電流分布図の表示対象平面となる。
【００６７】方向入力領域７２によって，Ｘ軸方向にベ
クトル分解した電流値成分の分布を表示するか，Ｙ軸方
向にベクトル分解した電流値成分の分布を表示するかを
選択することができる。
クトル分解した電流値成分の分布を表示するか，Ｙ軸方
向にベクトル分解した電流値成分の分布を表示するかを
選択することができる。
【００６８】回転指示領域７３では，グラフィック領域
における電流分布図の回転を指示することができる。回
転角度は９０度単位であり，「Ｆｒｏｎｔ」，「Ｌｅｆ
ｔ」，「Ｒｉｇｈｔ」，「Ｂａｃｋ」のいずれかが選択
されると，それぞれ前方向，左方向，右方向，後方向に
電流分布図を回転させる。
における電流分布図の回転を指示することができる。回
転角度は９０度単位であり，「Ｆｒｏｎｔ」，「Ｌｅｆ
ｔ」，「Ｒｉｇｈｔ」，「Ｂａｃｋ」のいずれかが選択
されると，それぞれ前方向，左方向，右方向，後方向に
電流分布図を回転させる。
【００６９】スケール入力領域７４では，電流値の分布
図の最大目盛りを選択することができる。実行ボタン７
５のクリックによって，新たに設定されたパラメータに
よる表示を実行する。
図の最大目盛りを選択することができる。実行ボタン７
５のクリックによって，新たに設定されたパラメータに
よる表示を実行する。
【００７０】モデルボタン７６のクリックによって，グ
ラフィック領域に表示されている電流分布図の平面が電
気回路装置のどの部分の平面かを確認させるためのモデ
ル図を表示する。図１３は，モデルボタン７６のクリッ
クによって表示されたモデル図表示画面の例を示す。こ
のモデル図の表示では，電流分布図の表示対象となって
いる平面を特定の色で塗りつぶして表示する。また，使
用する際の目安になるように，基準線は太線で表示す
る。
ラフィック領域に表示されている電流分布図の平面が電
気回路装置のどの部分の平面かを確認させるためのモデ
ル図を表示する。図１３は，モデルボタン７６のクリッ
クによって表示されたモデル図表示画面の例を示す。こ
のモデル図の表示では，電流分布図の表示対象となって
いる平面を特定の色で塗りつぶして表示する。また，使
用する際の目安になるように，基準線は太線で表示す
る。
【００７１】図１４は，本発明の実施例における電流Ｘ
Ｙグラフ図表示画面の例を示す。図１０に示すドロウタ
イプ入力領域５１において「Ｄｉｓｔｒｉｂｕｔｉｏｎ
−ＸＹ」が選択され，実行ボタン５２がクリックされる
と，例えば図１４に示すような電流ＸＹグラフ図が表示
される。この画面におけるモード指定領域８０におい
て，「Ｒｅａｌ／Ｉｍａｇｉｎａｒｙ」，「Ａｂｓｏｌ
ｕｔｅ／Ｐｈａｓｅ」の選択ボタンがあり，これによっ
て，電流の実数値と虚数値の重ね書きグラフ表示，また
は絶対値と位相の重ね書きグラフ表示のいずれかが選択
可能になっている。
Ｙグラフ図表示画面の例を示す。図１０に示すドロウタ
イプ入力領域５１において「Ｄｉｓｔｒｉｂｕｔｉｏｎ
−ＸＹ」が選択され，実行ボタン５２がクリックされる
と，例えば図１４に示すような電流ＸＹグラフ図が表示
される。この画面におけるモード指定領域８０におい
て，「Ｒｅａｌ／Ｉｍａｇｉｎａｒｙ」，「Ａｂｓｏｌ
ｕｔｅ／Ｐｈａｓｅ」の選択ボタンがあり，これによっ
て，電流の実数値と虚数値の重ね書きグラフ表示，また
は絶対値と位相の重ね書きグラフ表示のいずれかが選択
可能になっている。
【００７２】方向入力領域８１によって，Ｘ軸方向の電
流またはＹ軸方向の電流のいずれかを選択することがで
きる。スケール入力領域８２では，グラフの電流値の最
大目盛りを選択することができる。
流またはＹ軸方向の電流のいずれかを選択することがで
きる。スケール入力領域８２では，グラフの電流値の最
大目盛りを選択することができる。
【００７３】ワイヤ番号入力領域８３では，選択した軸
のワイヤ番号を指定することができる。ここで，ワイヤ
番号が連続している場合には，『始点番号−終点番号』
というように，ハイフンを使用して，入力を効率化させ
ることもできるようになっている。なお，ワイヤではな
く，平面上に流れる電流を表示する場合には，次のパッ
チ番号入力領域８４によって，パッチ番号の範囲を入力
する。ワイヤ番号入力領域８３とパッチ番号入力領域８
４との同時入力は不可である。
のワイヤ番号を指定することができる。ここで，ワイヤ
番号が連続している場合には，『始点番号−終点番号』
というように，ハイフンを使用して，入力を効率化させ
ることもできるようになっている。なお，ワイヤではな
く，平面上に流れる電流を表示する場合には，次のパッ
チ番号入力領域８４によって，パッチ番号の範囲を入力
する。ワイヤ番号入力領域８３とパッチ番号入力領域８
４との同時入力は不可である。
【００７４】実行ボタン８５のクリックによって，新た
に設定されたパラメータによる表示を実行する。次に，
モーメント法で算出した各要素の辺に流れる電流から，
格子設定処理手段１４により設定した格子の交点に流れ
る電流を計算する方法について，図１５ないし図１７に
示す電流計算説明図に従って説明する。
に設定されたパラメータによる表示を実行する。次に，
モーメント法で算出した各要素の辺に流れる電流から，
格子設定処理手段１４により設定した格子の交点に流れ
る電流を計算する方法について，図１５ないし図１７に
示す電流計算説明図に従って説明する。
【００７５】最初に，パッチが四角形要素である場合に
ついて説明する。図１５（Ａ）において，Ｍはパッチ境
界であって，モーメント法の要素の辺を表す。モーメン
ト法では，よく知られているように各要素の辺を流れる
電流値ｉi 〜ｉ7 が求められる。これらの要素（パッ
チ）は，大きさや形状が均一でないため，本発明では，
例えば図１５（Ｂ）に示すように，所定の間隔または利
用者により指定された間隔で格子Ｇを設定し，この格子
ごとに電流値を算出する。本実施例では，格子の交点Ｐ
の電流値を求めるものとする。格子の交点Ｐではなく，
中央の点の電流を求めるようにしても全く同様である。
格子の交点Ｐを流れる電流は，図１５（Ｂ）に斜線で示
した範囲を流れる電流値を代表したものとして計算す
る。
ついて説明する。図１５（Ａ）において，Ｍはパッチ境
界であって，モーメント法の要素の辺を表す。モーメン
ト法では，よく知られているように各要素の辺を流れる
電流値ｉi 〜ｉ7 が求められる。これらの要素（パッ
チ）は，大きさや形状が均一でないため，本発明では，
例えば図１５（Ｂ）に示すように，所定の間隔または利
用者により指定された間隔で格子Ｇを設定し，この格子
ごとに電流値を算出する。本実施例では，格子の交点Ｐ
の電流値を求めるものとする。格子の交点Ｐではなく，
中央の点の電流を求めるようにしても全く同様である。
格子の交点Ｐを流れる電流は，図１５（Ｂ）に斜線で示
した範囲を流れる電流値を代表したものとして計算す
る。
【００７６】格子の交点Ｐでの電流は，水平方向（Ｘ方
向），垂直方向（Ｙ方向）に分けて求める。例えば，図
１５（Ａ）に示す電流値ｉ1 〜ｉ7 から，図１５（Ｂ）
に斜線で示す範囲の電流値を求める場合には，６個のパ
ッチについて電流値を考慮する必要がある。以下，水平
方向の電流の計算について説明するが，垂直方向の計算
も同様に行うことができる。大まかな，処理の流れは以
下のとおりである。
向），垂直方向（Ｙ方向）に分けて求める。例えば，図
１５（Ａ）に示す電流値ｉ1 〜ｉ7 から，図１５（Ｂ）
に斜線で示す範囲の電流値を求める場合には，６個のパ
ッチについて電流値を考慮する必要がある。以下，水平
方向の電流の計算について説明するが，垂直方向の計算
も同様に行うことができる。大まかな，処理の流れは以
下のとおりである。
【００７７】(a) 電流値を求めようとする任意の点の範
囲がどのパッチに含まれるかを求め，その辺（エッジ）
を流れる電流を得る。 (b) 対辺の電流から電流密度を求める。
囲がどのパッチに含まれるかを求め，その辺（エッジ）
を流れる電流を得る。 (b) 対辺の電流から電流密度を求める。
【００７８】(c) 得られた電流密度を格子の方向に分
解，合成する。 (d) この電流密度から，その点が範囲とする電流を電流
密度として求める。 対辺の電流から電流密度を求める場合について，図１６
に従って説明する。今，水平方向の対辺の電流ｉ1 ，ｉ
2 から，図１６（Ａ）の点Ｑの電流密度を求めるものと
する。与えられるパラメータは，周波数ｆ，辺の長さｌ
1 ，ｌ2 ，辺に流れる電流ｉ1 ，ｉ2 ，要素長ｘ0 およ
び点Ｑの辺からの距離ｘである。
解，合成する。 (d) この電流密度から，その点が範囲とする電流を電流
密度として求める。 対辺の電流から電流密度を求める場合について，図１６
に従って説明する。今，水平方向の対辺の電流ｉ1 ，ｉ
2 から，図１６（Ａ）の点Ｑの電流密度を求めるものと
する。与えられるパラメータは，周波数ｆ，辺の長さｌ
1 ，ｌ2 ，辺に流れる電流ｉ1 ，ｉ2 ，要素長ｘ0 およ
び点Ｑの辺からの距離ｘである。
【００７９】この２つのモードに流れる電流ｉ1 ，ｉ2
の辺における電流密度Ｉ1 ，Ｉ2 は，それぞれＩ1 ＝ｉ
1 ／ｌ1 ，Ｉ2 ＝ｉ2 ／ｌ2 である。このパッチに流れ
る電流は，正弦（ｓｉｎ）カーブであるから，対辺に向
かって電流の大きさを図化すると，図１６（Ｂ）に示す
ようになる。
の辺における電流密度Ｉ1 ，Ｉ2 は，それぞれＩ1 ＝ｉ
1 ／ｌ1 ，Ｉ2 ＝ｉ2 ／ｌ2 である。このパッチに流れ
る電流は，正弦（ｓｉｎ）カーブであるから，対辺に向
かって電流の大きさを図化すると，図１６（Ｂ）に示す
ようになる。
【００８０】そこで，図１６（Ｃ）に示すように，電流
Ｉ1 の点Ｑにおける大きさをＩA として，ＩA を次のよ
うに求める。 ＩA ＝Ａ sin（２πλ）＝Ａ sin（２πｘ＊ｆ／ｃ） ここで，ｃは光速であり，３０万Ｋｍ／ｓｅｃである。
ｘの単位をメール，ｆの単位をメガヘルツとして， ＩA ＝Ａ sin（２πｘ＊ｆ／３００） となる。ただし，＊は乗算を表し，Ａは以下の値であ
る。
Ｉ1 の点Ｑにおける大きさをＩA として，ＩA を次のよ
うに求める。 ＩA ＝Ａ sin（２πλ）＝Ａ sin（２πｘ＊ｆ／ｃ） ここで，ｃは光速であり，３０万Ｋｍ／ｓｅｃである。
ｘの単位をメール，ｆの単位をメガヘルツとして， ＩA ＝Ａ sin（２πｘ＊ｆ／３００） となる。ただし，＊は乗算を表し，Ａは以下の値であ
る。
【００８１】
Ａ＝Ｉ1 ／｛ sin（２πｘ0 ＊ｆ／３００）
同様に，電流Ｉ2 の点Ｑにおける大きさをＩB とする
と，ＩB は次のようになる。
と，ＩB は次のようになる。
【００８２】
ＩB ＝Ｂ sin（２π（ｘ0 −ｘ）＊ｆ／３００）
となる。ただし，
Ｂ＝Ｉ2 ／｛ sin（２πｘ0 ＊ｆ／３００）
である。以上から，点Ｑにおける水平方向の電流密度Ｉ
は， Ｉ＝ＩA ＋ＩB として求まる。
は， Ｉ＝ＩA ＋ＩB として求まる。
【００８３】以上のようにして算出した任意の点の電流
密度から，その点が範囲とする電流を電流密度として求
める。この計算に必要となるパラメータは，パッチ境界
によって区分される格子の対象範囲についての計算済み
の電流密度（水平方向および垂直方向）と，格子の１辺
の長さａである。例えば，図１５（Ｂ）に示す斜線部の
範囲の電流密度を求める場合は，以下のようになる。
密度から，その点が範囲とする電流を電流密度として求
める。この計算に必要となるパラメータは，パッチ境界
によって区分される格子の対象範囲についての計算済み
の電流密度（水平方向および垂直方向）と，格子の１辺
の長さａである。例えば，図１５（Ｂ）に示す斜線部の
範囲の電流密度を求める場合は，以下のようになる。
【００８４】格子点Ｐを含む斜線部の範囲の電流密度を
求める場合，図１５（Ｂ）において，水平方向について
は点Ｐからそれぞれ１／（２ａ）の範囲の点Ａ〜点Ｄ間
を流れる電流の電流密度を求める。点Ａ〜点Ｂ間に流れ
る電流は，点Ａを含むパッチ（第１パッチという）に関
して求めた電流密度を使用し，点Ｂ〜点Ｃ間に流れる電
流は，点Ｐを含むパッチ（第２パッチという）に関して
求めた電流密度を使用し，点Ｃ〜点Ｄ間に流れる電流
は，点Ｄを含むパッチ（第３パッチという）に関して求
めた電流密度を使用する。
求める場合，図１５（Ｂ）において，水平方向について
は点Ｐからそれぞれ１／（２ａ）の範囲の点Ａ〜点Ｄ間
を流れる電流の電流密度を求める。点Ａ〜点Ｂ間に流れ
る電流は，点Ａを含むパッチ（第１パッチという）に関
して求めた電流密度を使用し，点Ｂ〜点Ｃ間に流れる電
流は，点Ｐを含むパッチ（第２パッチという）に関して
求めた電流密度を使用し，点Ｃ〜点Ｄ間に流れる電流
は，点Ｄを含むパッチ（第３パッチという）に関して求
めた電流密度を使用する。
【００８５】第１パッチの水平方向の電流密度をＩx1，
第２パッチの水平方向の電流密度をＩx2，第３パッチの
水平方向の電流密度をＩx3とする。点Ａ〜点Ｂ間に流れ
る電流ｉABは， ｉAB＝Ｉx1＊ＡＢ（ＡＢは線分） 点Ｂ〜点Ｃ間に流れる電流ｉBCは， ｉBC＝Ｉx2＊ＢＣ（ＢＣは線分） 点Ｃ〜点Ｄ間に流れる電流ｉCDは， ｉCD＝Ｉx3＊ＣＤ（ＣＤは線分） 以上から，図１５（Ｂ）に示す斜線部分の範囲に流れる
電流の電流密度Ｉx0は， Ｉx0＝（ｉAB＋ｉBC＋ｉCD）／ａ として求まる。垂直方向（Ｙ方向）についても同様に求
まる。
第２パッチの水平方向の電流密度をＩx2，第３パッチの
水平方向の電流密度をＩx3とする。点Ａ〜点Ｂ間に流れ
る電流ｉABは， ｉAB＝Ｉx1＊ＡＢ（ＡＢは線分） 点Ｂ〜点Ｃ間に流れる電流ｉBCは， ｉBC＝Ｉx2＊ＢＣ（ＢＣは線分） 点Ｃ〜点Ｄ間に流れる電流ｉCDは， ｉCD＝Ｉx3＊ＣＤ（ＣＤは線分） 以上から，図１５（Ｂ）に示す斜線部分の範囲に流れる
電流の電流密度Ｉx0は， Ｉx0＝（ｉAB＋ｉBC＋ｉCD）／ａ として求まる。垂直方向（Ｙ方向）についても同様に求
まる。
【００８６】以上は四角形要素についての電流密度の計
算であるが，モーメント法でのパッチが三角形要素であ
る場合には，図１７（Ａ）に示すように，３方向の電流
密度を考えて，３方向についてそれぞれ電流密度を求め
る。エッジに流れる電流値をｉ1 とした場合に，四角形
の対辺に相当するものは三角形の頂点であり，その部分
の電流値ｉ2 は，ｉ2 ＝０とする。入力として，次のも
のが与えられる。
算であるが，モーメント法でのパッチが三角形要素であ
る場合には，図１７（Ａ）に示すように，３方向の電流
密度を考えて，３方向についてそれぞれ電流密度を求め
る。エッジに流れる電流値をｉ1 とした場合に，四角形
の対辺に相当するものは三角形の頂点であり，その部分
の電流値ｉ2 は，ｉ2 ＝０とする。入力として，次のも
のが与えられる。
【００８７】ｉ1 ＝エッジに流れる電流値
ｉ2 ＝０
ｘ ＝頂点からの距離
ｌ1 ＝辺の長さｌ
これらの値から，図１６で説明した四角形の場合と同様
な計算により，図１７（Ｂ）に示すように，任意の点Ｐ
に対して，１つの頂点に向かう方向に対する電流密度Ｉ
を算出することができる。これを３辺について計算し，
合成することにより，水平および垂直方向の電流を算出
することができる。
な計算により，図１７（Ｂ）に示すように，任意の点Ｐ
に対して，１つの頂点に向かう方向に対する電流密度Ｉ
を算出することができる。これを３辺について計算し，
合成することにより，水平および垂直方向の電流を算出
することができる。
【００８８】また，パッチ要素が図１７（Ｃ）に示すよ
うな一般の四角形である場合，２方向の電流密度Ｉx ，
Ｉy を計算する。この場合，任意の点Ｐを通り，対辺を
ｍ対ｎに分ける点を通る直線を考える。この直線を利用
してｘ，ｘ0 を求め，Ｉ1 およびＩ2 からＩx を，Ｉ3
およびＩ4 からＩy を計算する。
うな一般の四角形である場合，２方向の電流密度Ｉx ，
Ｉy を計算する。この場合，任意の点Ｐを通り，対辺を
ｍ対ｎに分ける点を通る直線を考える。この直線を利用
してｘ，ｘ0 を求め，Ｉ1 およびＩ2 からＩx を，Ｉ3
およびＩ4 からＩy を計算する。
【００８９】図１７（Ｃ）で得られたＩx ，Ｉy を，格
子の座標系の電流密度Ｉx0，Ｉy0に変換する場合には，
図１７（Ｄ）に示すような角度φ，θの関係から， Ｉx0＝ cosφ＊Ｉy − cosθ＊Ｉx Ｉy0＝ sinφ＊Ｉy ＋ sinθ＊Ｉx の式により計算する。
子の座標系の電流密度Ｉx0，Ｉy0に変換する場合には，
図１７（Ｄ）に示すような角度φ，θの関係から， Ｉx0＝ cosφ＊Ｉy − cosθ＊Ｉx Ｉy0＝ sinφ＊Ｉy ＋ sinθ＊Ｉx の式により計算する。
【００９０】図１８は，電流ベクトル図の計算説明図で
ある。前述した計算により，ある格子点における水平方
向の電流ｘ（Ａ）と，垂直方向の電流ｙ（Ａ）とが求ま
ったとする。四角形要素の場合，図１８（Ａ）に示すよ
うに，電流の方向と大きさを求める。電流の大きさは， ｜Ｉ｜＝（ｘ2 ＋ｙ2 ）1/2 である。ただし，ｘ，ｙは複素数の電流値で，Ｉは求め
る電流である。
ある。前述した計算により，ある格子点における水平方
向の電流ｘ（Ａ）と，垂直方向の電流ｙ（Ａ）とが求ま
ったとする。四角形要素の場合，図１８（Ａ）に示すよ
うに，電流の方向と大きさを求める。電流の大きさは， ｜Ｉ｜＝（ｘ2 ＋ｙ2 ）1/2 である。ただし，ｘ，ｙは複素数の電流値で，Ｉは求め
る電流である。
【００９１】方向θを図１８（Ａ）に示す方向とする
と， θ＝ tan-1（ｙ／ｘ） である。ただし，ｘ＝０なら，θ＝π／２である。
と， θ＝ tan-1（ｙ／ｘ） である。ただし，ｘ＝０なら，θ＝π／２である。
【００９２】三角形要素の場合にも，図１８（Ｂ）に示
すように，四角形要素の場合と同様に電流の大きさと方
向を計算して，結果を表示する。 ｜Ｉ｜＝（ｘ2 ＋ｙ2 ）1/2 θ＝ tan-1（ｙ／ｘ） ただし，ｘ＝０なら，θ＝π／２である。
すように，四角形要素の場合と同様に電流の大きさと方
向を計算して，結果を表示する。 ｜Ｉ｜＝（ｘ2 ＋ｙ2 ）1/2 θ＝ tan-1（ｙ／ｘ） ただし，ｘ＝０なら，θ＝π／２である。
【００９３】電流分布図を表示するための電流の絶対
値，実数値，虚数値等も，以上の計算から同様に算出が
可能である。また，電流ＸＹグラフ図についても表示に
必要な値（絶対値と位相）を算出することができる。
値，実数値，虚数値等も，以上の計算から同様に算出が
可能である。また，電流ＸＹグラフ図についても表示に
必要な値（絶対値と位相）を算出することができる。
【００９４】なお，本実施例では，さらに解析対象の電
気回路装置モデルにおけるケーブル部分については，モ
ーメント法または分布定数線路近似法によりワイヤ要素
の電流分布を計算し，計算した電流の結果について，同
様に電流ＸＹグラフ図により電流表示を行うようにして
いる。
気回路装置モデルにおけるケーブル部分については，モ
ーメント法または分布定数線路近似法によりワイヤ要素
の電流分布を計算し，計算した電流の結果について，同
様に電流ＸＹグラフ図により電流表示を行うようにして
いる。
【００９５】図１９は，本発明を適用するシステムの構
成例を示す。部品配置・配線ＣＡＤデータファイル９０
は，電子系ＣＡＤシステムで作成したプリント板外形，
部品ドライバ／レシーバのピン情報，配線ライン，ビア
等のデータが格納されたファイルである。筐体構造ＣＡ
Ｄデータファイル９１は，構造系ＣＡＤシステムで作成
された解析分野での標準的なインタフェースであるＮＡ
ＳＴＲＡＮ形式の筐体構造データを保持するファイルで
ある。
成例を示す。部品配置・配線ＣＡＤデータファイル９０
は，電子系ＣＡＤシステムで作成したプリント板外形，
部品ドライバ／レシーバのピン情報，配線ライン，ビア
等のデータが格納されたファイルである。筐体構造ＣＡ
Ｄデータファイル９１は，構造系ＣＡＤシステムで作成
された解析分野での標準的なインタフェースであるＮＡ
ＳＴＲＡＮ形式の筐体構造データを保持するファイルで
ある。
【００９６】電磁界解析システム１００は，図１に示す
電磁界強度算出装置１０に相当するものであり，プレ処
理部１１０，ソルバ入力ファイル１２０，ライブラリ１
３０，ソルバ部１４０，補助機能部１５０，解析結果の
出力ファイル１６０，ポスト処理部１７０等からなる。
電磁界強度算出装置１０に相当するものであり，プレ処
理部１１０，ソルバ入力ファイル１２０，ライブラリ１
３０，ソルバ部１４０，補助機能部１５０，解析結果の
出力ファイル１６０，ポスト処理部１７０等からなる。
【００９７】プレ処理部１１０は，部品配置・配線ＣＡ
Ｄデータファイル９０からプリント板に関するデータを
入力して，ソルバ入力データに変換し，ソルバ入力ファ
イル１２０に出力する。また，筐体構造ＣＡＤデータフ
ァイル９１から筐体レベルの解析に必要となる筐体構造
データを入力して，ソルバ入力データに変換し，ソルバ
入力ファイル１２０に出力する。ソルバ入力データへの
変換では，プリント板間あるいは筐体外部へ出る各種ケ
ーブル経路を作成する。また，一度取り込んだプリント
板配線情報や筐体構造の変更を行う。また，プレ処理部
１１０は，筐体内部でのプリント板実装位置を指定する
機能などを持つ。また，パッチ番号付与部１１１を有す
る。
Ｄデータファイル９０からプリント板に関するデータを
入力して，ソルバ入力データに変換し，ソルバ入力ファ
イル１２０に出力する。また，筐体構造ＣＡＤデータフ
ァイル９１から筐体レベルの解析に必要となる筐体構造
データを入力して，ソルバ入力データに変換し，ソルバ
入力ファイル１２０に出力する。ソルバ入力データへの
変換では，プリント板間あるいは筐体外部へ出る各種ケ
ーブル経路を作成する。また，一度取り込んだプリント
板配線情報や筐体構造の変更を行う。また，プレ処理部
１１０は，筐体内部でのプリント板実装位置を指定する
機能などを持つ。また，パッチ番号付与部１１１を有す
る。
【００９８】ソルバ部１４０は，電磁界強度を解析する
モジュールで，ソルバ入力ファイル１２０からの入力デ
ータおよびライブラリ１３０から入力したプリント板で
使用されている部品のライブラリ情報をもとに，解析対
象種類毎にモーメント法または分布定数線路近似法によ
り解析を行い，その結果を解析結果の出力ファイル１６
０に出力する。
モジュールで，ソルバ入力ファイル１２０からの入力デ
ータおよびライブラリ１３０から入力したプリント板で
使用されている部品のライブラリ情報をもとに，解析対
象種類毎にモーメント法または分布定数線路近似法によ
り解析を行い，その結果を解析結果の出力ファイル１６
０に出力する。
【００９９】ポスト処理部１７０は，解析結果の出力フ
ァイル１６０から解析結果データを入力し，ディスプレ
イ１８０に指定された形式で図化して表示する。例え
ば，次のような形式でディスプレイ１８０に表示するこ
とができる。
ァイル１６０から解析結果データを入力し，ディスプレ
イ１８０に指定された形式で図化して表示する。例え
ば，次のような形式でディスプレイ１８０に表示するこ
とができる。
【０１００】 周波数スペクトル図
放射パターン図
電磁界マップ図
インピーダンススペクトル図
電流分布図
電流ベクトル図
電流ＸＹグラフ図
ワーストネットパターン図
補助機能部１５０は，ソルバ部１４０を処理途中で一旦
終了し，途中結果を確認し，再度処理を続行するソルバ
中断・リスタート機能を提供する。また，ソルバ部１４
０の処理がどの程度進んでいるかの処理段階を表示する
状況表示機能を提供する。
終了し，途中結果を確認し，再度処理を続行するソルバ
中断・リスタート機能を提供する。また，ソルバ部１４
０の処理がどの程度進んでいるかの処理段階を表示する
状況表示機能を提供する。
【０１０１】この電磁界解析システム１００によって，
プリント板，ケーブル，筐体を対象とした実用的な解析
を実現し，プリント板とケーブルのカップリング効果，
金属板によるシールド効果，物体表面の電流分布など，
電磁界放射の様子や放射メカニズムの把握に有力な情報
を得ることができる。
プリント板，ケーブル，筐体を対象とした実用的な解析
を実現し，プリント板とケーブルのカップリング効果，
金属板によるシールド効果，物体表面の電流分布など，
電磁界放射の様子や放射メカニズムの把握に有力な情報
を得ることができる。
【０１０２】
【発明の効果】以上説明したように，本発明によれば，
電気回路装置モデルのメッシュ状の微小要素に依存する
ことなく，均一的に分散する格子状の点が求められ，そ
れらの点を流れる電流が各種態様で可視化されて表示さ
れるので，電磁波の放射メカニズムを解明する有力な情
報を得ることが可能になる。
電気回路装置モデルのメッシュ状の微小要素に依存する
ことなく，均一的に分散する格子状の点が求められ，そ
れらの点を流れる電流が各種態様で可視化されて表示さ
れるので，電磁波の放射メカニズムを解明する有力な情
報を得ることが可能になる。
【図１】本発明の原理構成図である。
【図２】本発明における格子設定処理説明図である。
【図３】本発明による電流分布図の表示例を示す図であ
る。
る。
【図４】本発明による電流ベクトル図の表示例を示す図
である。
である。
【図５】本発明によるモデル図の表示例を示す図であ
る。
る。
【図６】本発明による矢印の色分け表示の例を説明する
図である。
図である。
【図７】本発明による電流ＸＹグラフ図の表示例を示す
図である。
図である。
【図８】本発明の実施例の処理フローの概要図である。
【図９】本発明の実施例における電流表示処理手段の処
理概要説明図である。
理概要説明図である。
【図１０】本発明の実施例の電流表示用パラメータ入力
画面の例を示す図である。
画面の例を示す図である。
【図１１】本発明の実施例における電流ベクトル図表示
画面の例を示す図である。
画面の例を示す図である。
【図１２】本発明の実施例における電流分布図表示画面
の例を示す図である。
の例を示す図である。
【図１３】本発明の実施例におけるモデル図表示画面の
例を示す図である。
例を示す図である。
【図１４】本発明の実施例における電流ＸＹグラフ図表
示画面の例を示す図である。
示画面の例を示す図である。
【図１５】本発明の実施例における電流計算説明図であ
る。
る。
【図１６】本発明の実施例における電流計算説明図であ
る。
る。
【図１７】本発明の実施例における電流計算説明図であ
る。
る。
【図１８】本発明の実施例における電流ベクトル図の計
算説明図である。
算説明図である。
【図１９】本発明を適用するシステムの構成例を示す図
である。
である。
１０ 電磁界強度算出装置
１１ データ入力部
１２ 電磁界強度計算部
１３ 解析情報出力部
１４ 格子設定処理手段
１５ 格子点電流計算手段
１６ 電流表示処理手段
１７ 電流分布図表示処理手段
１８ 電流ベクトル図表示処理手段
１９ 電流ＸＹグラフ図表示処理手段
２０ 平面選択処理手段
２１ モデル図表示処理手段
２２ 電気回路装置データ
２３ 入力装置
２４ 表示装置
─────────────────────────────────────────────────────
フロントページの続き
(56)参考文献 特開 平６−148247（ＪＰ，Ａ)
特開 平５−26927（ＪＰ，Ａ)
特開 平５−26930（ＪＰ，Ａ)
特開 平６−249899（ＪＰ，Ａ)
特開 平７−121579（ＪＰ，Ａ)
(58)調査した分野(Int.Cl.7，ＤＢ名)
G01R 29/08
G06F 17/50
Claims (11)
- 【請求項１】 解析対象となる電気回路装置の要素に流
れる電流をモーメント法により算出し，それらの算出し
た電流に基づいて電気回路装置が放射する電磁界強度を
計算して表示する電磁界強度算出装置において， 解析対象の電気回路装置モデルにおける平面を，少なく
ともその平面を包含する矩形エリアで囲み，その矩形エ
リア内を前記モーメント法での計算単位となる要素の大
きさとは別の指定された間隔または所定の間隔の格子状
の線で区分する格子設定処理手段と， モーメント法により算出された電気回路装置の要素に流
れる電流に基づいて，前記格子状の線の交点または格子
内の点の電流を計算する格子点電流計算手段と， この格子点電流計算手段によって計算した電流を可視化
して表示する電流表示処理手段とを備えたことを特徴と
する電磁界強度算出装置。 - 【請求項２】 請求項１記載の電磁界強度算出装置にお
いて， 前記電流表示処理手段は，前記格子点電流計算手段によ
り計算した電流を直交座標系にベクトル分解し，その結
果を立体的に分布図として隠面消去した表示面で表示す
る電流分布図表示処理手段を持つことを特徴とする電磁
界強度算出装置。 - 【請求項３】 請求項１記載の電磁界強度算出装置にお
いて， 前記電流表示処理手段は，前記格子点電流計算手段によ
り計算した電流を，電気回路装置のモデル図における該
当する面上に，電流の方向を矢印の向き，電流の大きさ
を矢印の長さで表した矢印図形で，前記各格子に対応し
て表示する電流ベクトル図表示処理手段を持つことを特
徴とする電磁界強度算出装置。 - 【請求項４】 請求項１記載の電磁界強度算出装置にお
いて， 前記格子設定処理手段は，前記電気回路装置モデルにお
いて三角形または四角形で構成される任意平面を，格子
状の線で区分した領域のみで表すようにしたことを特徴
とする電磁界強度算出装置。 - 【請求項５】 請求項１記載の電磁界強度算出装置にお
いて， 前記電気回路装置のモデル図内の１点または１要素を指
示する入力を受け付け，その指示された１点または１要
素を含む平面を，前記電流表示処理手段における電流の
表示対象平面とする平面選択処理手段を備えたことを特
徴とする電磁界強度算出装置。 - 【請求項６】 請求項１記載の電磁界強度算出装置にお
いて， 前記電流表示処理手段は，前記格子点電流計算手段によ
り計算した電流の結果における実数部と虚数部とを，外
部からの選択指示によって別に表示する手段を持つこと
を特徴とする電磁界強度算出装置。 - 【請求項７】 請求項１記載の電磁界強度算出装置にお
いて， 前記電流表示処理手段は，前記格子点電流計算手段によ
り計算した電流の結果を絶対値と位相とで表示する手段
を持つことを特徴とする電磁界強度算出装置。 - 【請求項８】 請求項２記載の電磁界強度算出装置にお
いて， 前記電流表示処理手段において電流の表示対象としてい
る平面が，電気回路装置の３次元モデル図におけるどの
平面であるかを示すモデル図を表示するモデル図表示処
理手段を備えたことを特徴とする電磁界強度算出装置。 - 【請求項９】 請求項３記載の電磁界強度算出装置にお
いて， 前記電流ベクトル図表示処理手段は，電流の大きさに応
じて矢印の表示態様を変えることにより，複数のスケー
ルを使い分けるようにしたことを特徴とする電磁界強度
算出装置。 - 【請求項１０】 請求項１記載の電磁界強度算出装置に
おいて， 前記電流表示処理手段は，前記格子点電流計算手段によ
り計算した電流の結果について，表示対象として指定さ
れた複数個の各点の位置を横軸とし，縦軸を各点におけ
る電流値としたＸＹグラフ形式で電流を表示する電流Ｘ
Ｙグラフ図表示処理手段を持つことを特徴とする電磁界
強度算出装置。 - 【請求項１１】 請求項１０記載の電磁界強度算出装置
において， 前記電流ＸＹグラフ図表示処理手段によりＸＹグラフ形
式で表示する点列の範囲を，その範囲の始点と終点とを
指定する入力で定めるようにしたことを特徴とする電磁
界強度算出装置。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP14994795A JP3428232B2 (ja) | 1995-06-16 | 1995-06-16 | 電磁界強度算出装置 |
| US08/661,916 US5812434A (en) | 1995-06-16 | 1996-06-12 | Electromagnetic field strength calculator having function of displaying currents to be analyzed |
| DE19623688A DE19623688B4 (de) | 1995-06-16 | 1996-06-14 | Berechnungsgerät und Berechnungsverfahren für eine elektromagnetische Feldstärke |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP14994795A JP3428232B2 (ja) | 1995-06-16 | 1995-06-16 | 電磁界強度算出装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH095375A JPH095375A (ja) | 1997-01-10 |
| JP3428232B2 true JP3428232B2 (ja) | 2003-07-22 |
Family
ID=15486065
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP14994795A Expired - Fee Related JP3428232B2 (ja) | 1995-06-16 | 1995-06-16 | 電磁界強度算出装置 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US5812434A (ja) |
| JP (1) | JP3428232B2 (ja) |
| DE (1) | DE19623688B4 (ja) |
Families Citing this family (48)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE19616772A1 (de) * | 1995-04-28 | 1996-10-31 | Fujitsu Ltd | Berechnungsgerät für eine elektromagnetische Feldintensität |
| JP2858559B2 (ja) * | 1996-04-25 | 1999-02-17 | 日本電気株式会社 | デバイスシミュレーション方法 |
| JP3825133B2 (ja) * | 1997-04-30 | 2006-09-20 | 富士通株式会社 | 電磁界強度算出装置及び方法並びにプログラム記録媒体 |
| JPH1115814A (ja) * | 1997-06-26 | 1999-01-22 | Fujitsu Ltd | モーメント法を用いたシミュレーション装置及び方法並びにプログラム記憶媒体 |
| JP3405905B2 (ja) * | 1997-08-12 | 2003-05-12 | 富士通株式会社 | 電磁界強度算出装置及び方法並びにプログラム記憶媒体 |
| JPH1194889A (ja) * | 1997-09-19 | 1999-04-09 | Fujitsu Ltd | 多層基板からの放射電磁波解析装置 |
| JPH11161690A (ja) * | 1997-11-26 | 1999-06-18 | Fujitsu Ltd | 電磁界強度算出装置、電磁界強度算出方法、及び電磁界強度算出手段を有するプログラムを格納した記録媒体 |
| DE19818290A1 (de) * | 1998-04-23 | 1999-10-28 | Bsh Bosch Siemens Hausgeraete | Meßeinrichtung für eine elektromagnetische Verträglichkeitsprüfung und -bestimmung |
| US6665849B2 (en) * | 1999-06-09 | 2003-12-16 | Interuniversitair Microelektronica Centrum Vzw | Method and apparatus for simulating physical fields |
| US6253202B1 (en) | 1998-09-18 | 2001-06-26 | Tacit Knowledge Systems, Inc. | Method, system and apparatus for authorizing access by a first user to a knowledge profile of a second user responsive to an access request from the first user |
| US6115709A (en) | 1998-09-18 | 2000-09-05 | Tacit Knowledge Systems, Inc. | Method and system for constructing a knowledge profile of a user having unrestricted and restricted access portions according to respective levels of confidence of content of the portions |
| WO2000017727A2 (en) | 1998-09-18 | 2000-03-30 | Tacit Knowledge Systems | Method and apparatus for querying a user knowledge profile |
| US6377949B1 (en) | 1998-09-18 | 2002-04-23 | Tacit Knowledge Systems, Inc. | Method and apparatus for assigning a confidence level to a term within a user knowledge profile |
| US8380875B1 (en) | 1998-09-18 | 2013-02-19 | Oracle International Corporation | Method and system for addressing a communication document for transmission over a network based on the content thereof |
| JPWO2000025145A1 (ja) * | 1998-10-26 | 2004-04-02 | 株式会社日立製作所 | 回転機の電磁場解析方法および電磁場解析装置 |
| JP3363808B2 (ja) * | 1998-11-24 | 2003-01-08 | 富士通株式会社 | シミュレーション装置及び方法並びにプログラム記録媒体 |
| JP3878345B2 (ja) * | 1998-12-01 | 2007-02-07 | 富士通株式会社 | シミュレーション装置及び方法並びにプログラム記録媒体 |
| US7016794B2 (en) * | 1999-03-16 | 2006-03-21 | Lsi Logic Corporation | Floor plan development electromigration and voltage drop analysis tool |
| US6381729B1 (en) * | 1999-03-31 | 2002-04-30 | International Business Machines Corporation | Data processing system and method for efficient determination of return current spread |
| JP3652168B2 (ja) * | 1999-04-28 | 2005-05-25 | 富士通株式会社 | 電磁界強度算出方法及びプログラム記録媒体 |
| JP3348709B2 (ja) * | 1999-11-24 | 2002-11-20 | 日本電気株式会社 | プリント回路基板設計支援装置及び制御プログラム記録媒体 |
| FR2810136B1 (fr) * | 2000-06-09 | 2002-10-18 | Thomson Csf | Algorithme de simulation electromagnetique, notamment des performances d'une antenne |
| US6711511B2 (en) * | 2000-07-06 | 2004-03-23 | Fujitsu Limited | Electromagnetic wave analyzing apparatus and computer readable storage medium |
| JP2002041598A (ja) * | 2000-07-27 | 2002-02-08 | Dainippon Printing Co Ltd | 電波障害対策システム |
| US6668251B1 (en) | 2000-11-01 | 2003-12-23 | Tacit Knowledge Systems, Inc. | Rendering discriminator members from an initial set of result data |
| WO2002039387A1 (en) * | 2000-11-07 | 2002-05-16 | Holographic Imaging Llc | Improved three dimensional display |
| US20020082812A1 (en) * | 2000-12-26 | 2002-06-27 | Fujitsu Limited | Apparatus and method for simulating the receiving characteristic of radio waves |
| US20040111386A1 (en) * | 2001-01-08 | 2004-06-10 | Goldberg Jonathan M. | Knowledge neighborhoods |
| JP2003157376A (ja) * | 2001-11-21 | 2003-05-30 | Ricoh Co Ltd | ネットワークシステム、識別情報管理方法、サーバ装置、プログラム、および記録媒体 |
| JP2003296395A (ja) * | 2002-04-01 | 2003-10-17 | Fujitsu Ltd | 相互インミッタンス計算プログラム、相互インミッタンス計算装置、相互インミッタンス計算方法および電磁界強度計算プログラム |
| JP4096611B2 (ja) * | 2002-05-07 | 2008-06-04 | 富士ゼロックス株式会社 | 設計支援装置 |
| JP2004020403A (ja) * | 2002-06-18 | 2004-01-22 | Hitachi Ltd | 電子機器表示装置、電子機器、及びその製造方法 |
| US9805373B1 (en) | 2002-11-19 | 2017-10-31 | Oracle International Corporation | Expertise services platform |
| JP4136673B2 (ja) * | 2003-01-10 | 2008-08-20 | 株式会社リコー | 機器安全性検証装置、機器安全性検証方法、および記憶媒体 |
| US7043705B2 (en) * | 2003-07-08 | 2006-05-09 | Texas Instruments Incorporated | Estimating current density parameters on signal leads of an integrated circuit |
| US7222033B1 (en) | 2003-08-18 | 2007-05-22 | Steven Lynn Newson | Electromagnetic emissions and susceptibility calculating method and apparatus |
| US7248263B2 (en) * | 2004-05-10 | 2007-07-24 | International Business Machines Corporation | User interface for representing multiple data sets on a common graph |
| US7200825B2 (en) * | 2004-08-27 | 2007-04-03 | International Business Machines Corporation | Methodology of quantification of transmission probability for minority carrier collection in a semiconductor chip |
| JP4499781B2 (ja) | 2005-03-28 | 2010-07-07 | 富士通株式会社 | 電磁界強度算出方法、電磁界強度算出装置、制御プログラム |
| JP4855853B2 (ja) * | 2006-07-05 | 2012-01-18 | 富士通株式会社 | 解析装置、コンピュータの制御方法およびモデル作成プログラム |
| JP2008059274A (ja) * | 2006-08-31 | 2008-03-13 | Toyota Motor Corp | 振動解析システム |
| JP2011008524A (ja) * | 2009-06-25 | 2011-01-13 | Toshiba Corp | 設計支援装置および設計支援方法 |
| JP5293501B2 (ja) * | 2009-08-26 | 2013-09-18 | 富士通株式会社 | ３次元データ表示装置およびプログラム |
| US9443046B2 (en) * | 2012-10-16 | 2016-09-13 | Toyota Motor Engineering & Manufacturing North America, Inc. | Determining the electromagnetic field in a computer aided design environment |
| JP6458373B2 (ja) * | 2014-06-30 | 2019-01-30 | キヤノンマーケティングジャパン株式会社 | 情報処理装置、情報処理システム、その制御方法及びプログラム |
| JP2016143082A (ja) * | 2015-01-29 | 2016-08-08 | 富士通株式会社 | 算出プログラム、算出装置および算出方法 |
| JP6586637B2 (ja) * | 2015-12-25 | 2019-10-09 | 富士通株式会社 | 物理量分布計算プログラム、物理量分布計算方法、および情報処理装置 |
| JP2021076942A (ja) * | 2019-11-05 | 2021-05-20 | 株式会社Ｎｔｔファシリティーズ | 表示制御装置及び表示制御方法 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2635617B2 (ja) * | 1987-09-29 | 1997-07-30 | 株式会社東芝 | 半導体素子特性評価用の直交格子点の発生方法 |
| US5005147A (en) * | 1988-12-30 | 1991-04-02 | The United States Of America As Represented By The Administrator, The National Aeronautics And Space Administration | Method and apparatus for sensor fusion |
| DE4013017A1 (de) * | 1990-04-24 | 1991-10-31 | Siemens Ag | Verfahren zur ueberpruefung von in einem koerper verdeckt verlaufenden strompfaden |
| DE69427417T2 (de) * | 1993-03-08 | 2002-05-29 | Koninklijke Philips Electronics N.V., Eindhoven | PCB-Simulation auf der Basis von reduzierten Ersatzschaltungen |
| US5537329A (en) * | 1994-06-30 | 1996-07-16 | At&T Corp. | Apparatus and method for analyzing circuits |
| US5566083A (en) * | 1994-10-18 | 1996-10-15 | The Research Foundation Of State University Of New York | Method for analyzing voltage fluctuations in multilayered electronic packaging structures |
| DE19616772A1 (de) * | 1995-04-28 | 1996-10-31 | Fujitsu Ltd | Berechnungsgerät für eine elektromagnetische Feldintensität |
-
1995
- 1995-06-16 JP JP14994795A patent/JP3428232B2/ja not_active Expired - Fee Related
-
1996
- 1996-06-12 US US08/661,916 patent/US5812434A/en not_active Expired - Lifetime
- 1996-06-14 DE DE19623688A patent/DE19623688B4/de not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| US5812434A (en) | 1998-09-22 |
| JPH095375A (ja) | 1997-01-10 |
| DE19623688B4 (de) | 2006-06-14 |
| DE19623688A1 (de) | 1996-12-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3428232B2 (ja) | 電磁界強度算出装置 | |
| JP2768900B2 (ja) | 電磁界強度算出装置 | |
| US20030204343A1 (en) | Electromagnetic field analysis method based on FDTD method, medium representation method in electromagnetic field analysis, simulation device, and storage medium | |
| Chang et al. | A surface formulation for characteristic modes of material bodies | |
| US5650935A (en) | Apparatus for calculating intensity of an electromagnetic field | |
| US5903477A (en) | Simulation apparatus and simulation method for electromagnetic field intensity using moment method | |
| Tirkas et al. | Finite-difference time-domain method for electromagnetic radiation, interference, and interaction with complex structures | |
| DE3706509A1 (de) | Elektronische geometrieprozessorschaltung fuer ein graphikanzeigesystem | |
| JP3488629B2 (ja) | 放射電磁界耐性算出装置及び方法並びにプログラム記録媒体 | |
| US20020140435A1 (en) | Method and device generating data input to an electromagnetic field intensity calculating device | |
| JPH1115814A (ja) | モーメント法を用いたシミュレーション装置及び方法並びにプログラム記憶媒体 | |
| CN117494447A (zh) | 一种基于波谱展开和计算电磁算法的辐射近场源构建方法 | |
| CN106709150B (zh) | 一种电流分布和近场电磁分布三维可视化方法 | |
| US6555998B1 (en) | Apparatus and method for simulation of electromagnetic field and storage medium storing programs therfor | |
| JP3505062B2 (ja) | モーメント法を用いたシミュレーション装置及びシミュレーション方法並びにプログラム記憶媒体 | |
| Park et al. | Rapid calculation of the Green's function in the shielded planar structures | |
| JP3443162B2 (ja) | 電磁界強度算出装置 | |
| Mosig | Scattering by arbitrarily-shaped slots in thick conducting screens: an approximate solution | |
| Park et al. | Rapid summation of the Green's function for the rectangular waveguide | |
| Ramesan et al. | Modeling of radiation source using an equivalent dipole moment model | |
| JPH0915281A (ja) | 電磁場解析方法及び解析装置 | |
| JP3363657B2 (ja) | 電磁界強度算出装置 | |
| Öing et al. | Calculation of radiated electromagnetic fields from electronic systems | |
| Smith et al. | Goal-orientated computational steering | |
| Wnuk et al. | The analysis of microstrip antennas using the FDTD method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20030114 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20030422 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080516 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090516 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090516 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100516 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100516 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110516 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120516 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130516 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140516 Year of fee payment: 11 |
|
| LAPS | Cancellation because of no payment of annual fees |