JPH07234890A - 電磁界強度算出装置 - Google Patents
電磁界強度算出装置Info
- Publication number
- JPH07234890A JPH07234890A JP2710994A JP2710994A JPH07234890A JP H07234890 A JPH07234890 A JP H07234890A JP 2710994 A JP2710994 A JP 2710994A JP 2710994 A JP2710994 A JP 2710994A JP H07234890 A JPH07234890 A JP H07234890A
- Authority
- JP
- Japan
- Prior art keywords
- electromagnetic field
- calculation
- calculating
- distributed constant
- constant line
- 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
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
(57)【要約】
【目的】本発明は、電気回路装置の放射する電磁界強度
を算出する電磁界強度算出装置に関し、高い精度で、高
速に、しかも少ない記憶資源でもって電磁界強度を算出
できるようにすることを目的とする。 【構成】電気回路装置の各装置部分を、分布定数線路近
似法の適用可能な装置部分と、適用不可能な装置部分と
に区分けする区分手段１１と、適用可能と区分けされた
装置部分に対して分布定数線路近似法を施すことで、そ
の装置部分の電流分布を算出する第１の算出手段１２
と、第１の算出手段１２の算出する電流分布を用いつ
つ、適用不可能と区分けされた装置部分に対してモーメ
ント法を施すことで、その装置部分の電流分布を算出す
る第２の算出手段１３と、第１及び第２の算出手段１
２，１３の算出する電流分布から、電気回路装置の放射
する電磁界強度を算出する第３の算出手段１４とを備え
るように構成する。
を算出する電磁界強度算出装置に関し、高い精度で、高
速に、しかも少ない記憶資源でもって電磁界強度を算出
できるようにすることを目的とする。 【構成】電気回路装置の各装置部分を、分布定数線路近
似法の適用可能な装置部分と、適用不可能な装置部分と
に区分けする区分手段１１と、適用可能と区分けされた
装置部分に対して分布定数線路近似法を施すことで、そ
の装置部分の電流分布を算出する第１の算出手段１２
と、第１の算出手段１２の算出する電流分布を用いつ
つ、適用不可能と区分けされた装置部分に対してモーメ
ント法を施すことで、その装置部分の電流分布を算出す
る第２の算出手段１３と、第１及び第２の算出手段１
２，１３の算出する電流分布から、電気回路装置の放射
する電磁界強度を算出する第３の算出手段１４とを備え
るように構成する。
Description
【０００１】
【産業上の利用分野】本発明は、電気回路装置の放射す
る電磁界強度を算出する電磁界強度算出装置に関し、特
に、高い精度で、高速に、しかも少ない記憶資源でもっ
て電磁界強度を算出できる電磁界強度算出装置に関す
る。
る電磁界強度を算出する電磁界強度算出装置に関し、特
に、高い精度で、高速に、しかも少ない記憶資源でもっ
て電磁界強度を算出できる電磁界強度算出装置に関す
る。
【０００２】電気回路装置を設計していくにあたって、
電気回路装置の放射する電磁界強度を算出することが要
求されることがある。このような電磁界強度を算出する
電磁界強度算出装置は、高い精度で、高速に、しかも少
ない記憶資源でもって電磁界強度を算出できることが好
ましい。
電気回路装置の放射する電磁界強度を算出することが要
求されることがある。このような電磁界強度を算出する
電磁界強度算出装置は、高い精度で、高速に、しかも少
ない記憶資源でもって電磁界強度を算出できることが好
ましい。
【０００３】
【従来の技術】任意形状の物体から放射される電磁波の
電磁界強度は、物体各部に流れる電流が分かれば、公知
の理論式を用いて容易に計算することができる。この電
流値は、理論的には、マックスウェルの電磁波動方程式
を与えられた境界条件の下で解くことで得られるもので
あるが、任意形状物体を対象とした複雑な境界条件の下
での理論的な解法は現在知られていない。
電磁界強度は、物体各部に流れる電流が分かれば、公知
の理論式を用いて容易に計算することができる。この電
流値は、理論的には、マックスウェルの電磁波動方程式
を与えられた境界条件の下で解くことで得られるもので
あるが、任意形状物体を対象とした複雑な境界条件の下
での理論的な解法は現在知られていない。
【０００４】従って、現在の電磁界強度算出装置で用い
られている電流を求める解法は、難易の程度はあるにせ
よ、すべて近似的なものである。現在、この近似的な解
法として、微小ループアンテナ近似法と、分布定数線路
近似法と、モーメント法という３種類のものが知られて
いる。
られている電流を求める解法は、難易の程度はあるにせ
よ、すべて近似的なものである。現在、この近似的な解
法として、微小ループアンテナ近似法と、分布定数線路
近似法と、モーメント法という３種類のものが知られて
いる。
【０００５】微小ループアンテナ近似法は、波源回路と
負荷回路との間を接続する配線をループアンテナとして
扱い、ループ上の電流は平坦なものと仮定して、これを
集中定数回路の計算手法で求める方法である。図１８
に、この微小ループアンテナ近似法の構成を図示する。
この微小ループアンテナ近似法による計算は最も簡単で
あるが、ループの寸法が電磁波の波長に比べて無視でき
ない条件では精度が極めて低下することから、現実には
殆ど用いられていない。
負荷回路との間を接続する配線をループアンテナとして
扱い、ループ上の電流は平坦なものと仮定して、これを
集中定数回路の計算手法で求める方法である。図１８
に、この微小ループアンテナ近似法の構成を図示する。
この微小ループアンテナ近似法による計算は最も簡単で
あるが、ループの寸法が電磁波の波長に比べて無視でき
ない条件では精度が極めて低下することから、現実には
殆ど用いられていない。
【０００６】一方、分布定数線路近似法は、１次元の構
造物として近似できる物体に対して、分布定数線路の方
程式を適用することで電流を求める方法である。計算は
比較的簡単であり、計算時間及び記憶容量も解析要素数
にほぼ比例して増加するに止まり、線路の反射や共振等
の現象も含めた解析ができることから、１次元の近似が
成立する物体に対しては、高速・高精度の解析ができ
る。図１９に、この分布定数線路近似法の構成を図示す
る。
造物として近似できる物体に対して、分布定数線路の方
程式を適用することで電流を求める方法である。計算は
比較的簡単であり、計算時間及び記憶容量も解析要素数
にほぼ比例して増加するに止まり、線路の反射や共振等
の現象も含めた解析ができることから、１次元の近似が
成立する物体に対しては、高速・高精度の解析ができ
る。図１９に、この分布定数線路近似法の構成を図示す
る。
【０００７】一方、モーメント法は、マックスウェルの
電磁波動方程式から導かれる積分方程式の解法の１つで
あり、図２０に示すような３次元の任意形状物体を扱う
ことができる。具体的には、物体を小さな要素に分割し
て電流の計算を行うものであり、計算時間及び記憶容量
は、解析要素数のほぼ２乗に比例して増大することにな
る。
電磁波動方程式から導かれる積分方程式の解法の１つで
あり、図２０に示すような３次元の任意形状物体を扱う
ことができる。具体的には、物体を小さな要素に分割し
て電流の計算を行うものであり、計算時間及び記憶容量
は、解析要素数のほぼ２乗に比例して増大することにな
る。
【０００８】
【発明が解決しようとする課題】確かに、解析対象の電
気回路装置が１次元の構造物として近似できる場合に
は、分布定数線路近似法を用いることで、物体を流れる
電流を高速・高精度に求めることができ、電気回路装置
の放射する電磁界強度を高速・高精度に求めることがで
きる。
気回路装置が１次元の構造物として近似できる場合に
は、分布定数線路近似法を用いることで、物体を流れる
電流を高速・高精度に求めることができ、電気回路装置
の放射する電磁界強度を高速・高精度に求めることがで
きる。
【０００９】しかしながら、この分布定数線路近似法で
は、解析対象の電気回路装置が１次元の構造物として近
似できない場合には、原理的に電流を求めることができ
ないという問題点がある。
は、解析対象の電気回路装置が１次元の構造物として近
似できない場合には、原理的に電流を求めることができ
ないという問題点がある。
【００１０】すなわち、この分布定数線路近似法では、
接地層を有する多層プリント板や同軸ケーブル等の電流
については、高速・高精度に求めることができるもの
の、接地層を持たないことで周囲への電磁界の放射が多
いとともに、周囲の電磁界による影響を受ける特性を持
つ片面・両面プリント板や、筐体等の任意形状物体につ
いては電流を求めることができないことから、プリント
板から筐体まで含めた装置全体の解析は原理的に不可能
であるという問題点がある。
接地層を有する多層プリント板や同軸ケーブル等の電流
については、高速・高精度に求めることができるもの
の、接地層を持たないことで周囲への電磁界の放射が多
いとともに、周囲の電磁界による影響を受ける特性を持
つ片面・両面プリント板や、筐体等の任意形状物体につ
いては電流を求めることができないことから、プリント
板から筐体まで含めた装置全体の解析は原理的に不可能
であるという問題点がある。
【００１１】一方、確かに、モーメント法を用いると、
プリント板から筐体まで含めた装置全体の解析が可能で
あるものの、解析対象となる電気回路装置の規模が大き
くなると、計算量が膨大なものとなるために、現在のコ
ンピュータの能力では、実用的な時間内に解析を行うこ
とができないという問題点がある。そして、記憶資源が
膨大なものとなるという問題点もある。
プリント板から筐体まで含めた装置全体の解析が可能で
あるものの、解析対象となる電気回路装置の規模が大き
くなると、計算量が膨大なものとなるために、現在のコ
ンピュータの能力では、実用的な時間内に解析を行うこ
とができないという問題点がある。そして、記憶資源が
膨大なものとなるという問題点もある。
【００１２】本発明はかかる事情に鑑みてなされたもの
であって、高い精度で、高速に、しかも少ない記憶資源
でもって、複雑な構造体の電気回路装置の放射する電磁
界強度を算出できるようにする新たな電磁界強度算出装
置の提供を目的とする。
であって、高い精度で、高速に、しかも少ない記憶資源
でもって、複雑な構造体の電気回路装置の放射する電磁
界強度を算出できるようにする新たな電磁界強度算出装
置の提供を目的とする。
【００１３】
【課題を解決するための手段】図１に本発明の原理構成
を図示する。図中、１は本発明により構成される電磁界
強度算出装置であって、本発明に従って、解析対象とな
る電気回路装置の放射する電磁界強度を算出する。
を図示する。図中、１は本発明により構成される電磁界
強度算出装置であって、本発明に従って、解析対象とな
る電気回路装置の放射する電磁界強度を算出する。
【００１４】この電磁界強度算出装置１は、入力手段１
０と、区分手段１１と、第１の算出手段１２と、第２の
算出手段１３と、第３の算出手段１４と、出力手段１５
とを備える。
０と、区分手段１１と、第１の算出手段１２と、第２の
算出手段１３と、第３の算出手段１４と、出力手段１５
とを備える。
【００１５】入力手段１０は、解析対象となる電気回路
装置の構造情報を外部から受け取る。区分手段１１は、
解析対象となる電気回路装置の各装置部分を、分布定数
線路近似法の適用可能な装置部分と、分布定数線路近似
法の適用不可能な装置部分とに区分けする。
装置の構造情報を外部から受け取る。区分手段１１は、
解析対象となる電気回路装置の各装置部分を、分布定数
線路近似法の適用可能な装置部分と、分布定数線路近似
法の適用不可能な装置部分とに区分けする。
【００１６】第１の算出手段１２は、区分手段１１が適
用可能と区分けする装置部分に対して分布定数線路近似
法を施すことで、その装置部分の電流分布を算出する。
第２の算出手段１３は、相互インピーダンス算出手段１
３０と、誘起波源算出手段１３１と、電流分布算出手段
１３２とから構成されて、第１の算出手段１２の算出す
る電流分布を用いつつ、区分手段１１が適用不可能と区
分けする装置部分に対してモーメント法を施すことで、
その装置部分の電流分布を算出する。
用可能と区分けする装置部分に対して分布定数線路近似
法を施すことで、その装置部分の電流分布を算出する。
第２の算出手段１３は、相互インピーダンス算出手段１
３０と、誘起波源算出手段１３１と、電流分布算出手段
１３２とから構成されて、第１の算出手段１２の算出す
る電流分布を用いつつ、区分手段１１が適用不可能と区
分けする装置部分に対してモーメント法を施すことで、
その装置部分の電流分布を算出する。
【００１７】第３の算出手段１４は、第１の算出手段１
２の算出する電流分布と、第２の算出手段１３の算出す
る電流分布とから、電気回路装置の放射する電磁界強度
を算出する。
２の算出する電流分布と、第２の算出手段１３の算出す
る電流分布とから、電気回路装置の放射する電磁界強度
を算出する。
【００１８】出力手段１５は、第３の算出手段１４の算
出する電磁界強度を外部に出力する。
出する電磁界強度を外部に出力する。
【００１９】
【作用】本発明では、入力手段１０が、外部から解析対
象となる電気回路装置の持つ形状や寸法や材料特性や電
気特性等の構造情報を受け取って区分手段１１に与える
と、区分手段１１は、その構造情報を分析することで、
解析対象となる電気回路装置の各装置部分を、分布定数
線路近似法の適用可能な装置部分と、分布定数線路近似
法の適用不可能な装置部分とに区分けする。
象となる電気回路装置の持つ形状や寸法や材料特性や電
気特性等の構造情報を受け取って区分手段１１に与える
と、区分手段１１は、その構造情報を分析することで、
解析対象となる電気回路装置の各装置部分を、分布定数
線路近似法の適用可能な装置部分と、分布定数線路近似
法の適用不可能な装置部分とに区分けする。
【００２０】例えば、信号配線層と接地配線層とを有す
るプリント板部分や、同軸構造を有するケーブル部分
や、平行配線構造部分や、平行平板構造部分や、平板と
配線とが平行に配設される構造部分といった１次元の構
造物として近似できる装置部分については分布定数線路
近似法の適用可能な装置部分と判断し、それ以外につい
ては、分布定数線路近似法の適用不可能な装置部分と判
断して区分けしていくのである。なお、このとき、例え
ば、同軸構造を有するケーブルの外導体部分については
分布定数線路近似法では十分精度がでないことに対応さ
せて、このような原則的には分布定数線路近似法の適用
可能な装置部分であっても、分布定数線路近似法の適用
不可能な装置部分と判断していくことがある。
るプリント板部分や、同軸構造を有するケーブル部分
や、平行配線構造部分や、平行平板構造部分や、平板と
配線とが平行に配設される構造部分といった１次元の構
造物として近似できる装置部分については分布定数線路
近似法の適用可能な装置部分と判断し、それ以外につい
ては、分布定数線路近似法の適用不可能な装置部分と判
断して区分けしていくのである。なお、このとき、例え
ば、同軸構造を有するケーブルの外導体部分については
分布定数線路近似法では十分精度がでないことに対応さ
せて、このような原則的には分布定数線路近似法の適用
可能な装置部分であっても、分布定数線路近似法の適用
不可能な装置部分と判断していくことがある。
【００２１】そして、この区分け処理を実行すると、区
分手段１１は、第１の算出手段１２に対して、分布定数
線路近似法の適用可能な装置部分の構造情報を通知する
とともに、第２の算出手段１３に対して、分布定数線路
近似法の適用可能な装置部分の構造情報と、分布定数線
路近似法の適用不可能な装置部分の構造情報とを区別し
つつ通知する。
分手段１１は、第１の算出手段１２に対して、分布定数
線路近似法の適用可能な装置部分の構造情報を通知する
とともに、第２の算出手段１３に対して、分布定数線路
近似法の適用可能な装置部分の構造情報と、分布定数線
路近似法の適用不可能な装置部分の構造情報とを区別し
つつ通知する。
【００２２】区分手段１１から、分布定数線路近似法の
適用可能な装置部分の構造情報を受け取ると、第１の算
出手段１２は、その装置部分に含まれる波源を用いて、
その装置部分に対して分布定数線路近似法を施すこと
で、その装置部分の電流分布を算出する。
適用可能な装置部分の構造情報を受け取ると、第１の算
出手段１２は、その装置部分に含まれる波源を用いて、
その装置部分に対して分布定数線路近似法を施すこと
で、その装置部分の電流分布を算出する。
【００２３】一方、区分手段１１から、分布定数線路近
似法の適用可能な装置部分の構造情報と、分布定数線路
近似法の適用不可能な装置部分の構造情報とを区別しつ
つ受け取ると、第２の算出手段１３は、最初に、相互イ
ンピーダンス算出手段１３０の算出処理に従って、分布
定数線路近似法の適用不可能な装置部分間の相互インピ
ーダンスを算出するとともに、分布定数線路近似法の適
用不可能な装置部分と、分布定数線路近似法の適用可能
な装置部分との間の相互インピーダンスを算出する。
似法の適用可能な装置部分の構造情報と、分布定数線路
近似法の適用不可能な装置部分の構造情報とを区別しつ
つ受け取ると、第２の算出手段１３は、最初に、相互イ
ンピーダンス算出手段１３０の算出処理に従って、分布
定数線路近似法の適用不可能な装置部分間の相互インピ
ーダンスを算出するとともに、分布定数線路近似法の適
用不可能な装置部分と、分布定数線路近似法の適用可能
な装置部分との間の相互インピーダンスを算出する。
【００２４】続いて、誘起波源算出手段１３１の算出処
理に従って、第１の算出手段１２の算出した電流分布
と、相互インピーダンス算出手段１３０の算出した異種
間の相互インピーダンスとのマトリクス演算から、第１
の算出手段１２の算出した電流分布が分布定数線路近似
法の適用不可能な装置部分に誘起する波源を算出する。
理に従って、第１の算出手段１２の算出した電流分布
と、相互インピーダンス算出手段１３０の算出した異種
間の相互インピーダンスとのマトリクス演算から、第１
の算出手段１２の算出した電流分布が分布定数線路近似
法の適用不可能な装置部分に誘起する波源を算出する。
【００２５】続いて、電流分布算出手段１３２の算出処
理に従って、誘起波源算出手段１３１の算出した波源
と、分布定数線路近似法の適用不可能な装置部分の持つ
本来の波源との合成値を波源として使用して、相互イン
ピーダンス算出手段１３０の算出した同種間の相互イン
ピーダンスを係数行列とするモーメント法の連立一次方
程式を解くことで、分布定数線路近似法の適用不可能な
装置部分の電流分布を算出する。
理に従って、誘起波源算出手段１３１の算出した波源
と、分布定数線路近似法の適用不可能な装置部分の持つ
本来の波源との合成値を波源として使用して、相互イン
ピーダンス算出手段１３０の算出した同種間の相互イン
ピーダンスを係数行列とするモーメント法の連立一次方
程式を解くことで、分布定数線路近似法の適用不可能な
装置部分の電流分布を算出する。
【００２６】このようにして、第１の算出手段１２が、
分布定数線路近似法の適用可能な装置部分の電流分布を
分布定数線路近似法に従って算出し、第２の算出手段１
３が、分布定数線路近似法の適用不可能な装置部分の電
流分布をモーメント法に従って算出すると、第３の算出
手段１４は、これらの電流分布から、電気回路装置の放
射する電磁界強度を算出し、この算出結果を受けて、出
力手段１５は、算出された電磁界強度を外部に出力す
る。
分布定数線路近似法の適用可能な装置部分の電流分布を
分布定数線路近似法に従って算出し、第２の算出手段１
３が、分布定数線路近似法の適用不可能な装置部分の電
流分布をモーメント法に従って算出すると、第３の算出
手段１４は、これらの電流分布から、電気回路装置の放
射する電磁界強度を算出し、この算出結果を受けて、出
力手段１５は、算出された電磁界強度を外部に出力す
る。
【００２７】このように、本発明では、分布定数線路近
似法の適用可能な装置部分については、適用不可能な装
置部分の持つ電流値が小さいことに対応して、その電流
値の誘起する波源を考慮することなく、計算の簡単な分
布定数線路近似法に従って電流分布を求め、一方、分布
定数線路近似法の適用不可能な装置部分については、適
用可能な装置部分の持つ電流値が大きいことに対応し
て、その電流値の誘起する波源を考慮しつつ、高精度の
モーメント法に従って電流分布を求める構成を採るもの
であることから、高い精度で、高速に、しかも少ない記
憶資源でもって、複雑な構造体の電気回路装置の放射す
る電磁界強度を算出できるようになる。
似法の適用可能な装置部分については、適用不可能な装
置部分の持つ電流値が小さいことに対応して、その電流
値の誘起する波源を考慮することなく、計算の簡単な分
布定数線路近似法に従って電流分布を求め、一方、分布
定数線路近似法の適用不可能な装置部分については、適
用可能な装置部分の持つ電流値が大きいことに対応し
て、その電流値の誘起する波源を考慮しつつ、高精度の
モーメント法に従って電流分布を求める構成を採るもの
であることから、高い精度で、高速に、しかも少ない記
憶資源でもって、複雑な構造体の電気回路装置の放射す
る電磁界強度を算出できるようになる。
【００２８】
【実施例】以下、実施例に従って本発明を詳細に説明す
る。図２に、図１で説明した本発明の電磁界強度算出装
置１の実行する処理フローの一実施例を図示する。
る。図２に、図１で説明した本発明の電磁界強度算出装
置１の実行する処理フローの一実施例を図示する。
【００２９】ここで、図中、２０は入力データファイル
であって、解析対象となる電気回路装置の構造体情報を
管理するもの、２１は出力データファイルであって、算
出された電磁界強度を格納するものである。
であって、解析対象となる電気回路装置の構造体情報を
管理するもの、２１は出力データファイルであって、算
出された電磁界強度を格納するものである。
【００３０】図３に、入力データファイル２０の管理す
るデータのデータ構造を図示する。この管理データは、
手入力やＣＡＤシステムから入力されるものであって、
この図に示すように、解析ステートメントブロックと、
モデル記述ステートメントブロックと、出力制御ステー
トメントブロックとから構成される。
るデータのデータ構造を図示する。この管理データは、
手入力やＣＡＤシステムから入力されるものであって、
この図に示すように、解析ステートメントブロックと、
モデル記述ステートメントブロックと、出力制御ステー
トメントブロックとから構成される。
【００３１】この解析ステートメントブロックは、図４
に示すように、周波数や収束判定値といった各種の解析
条件を管理するものであり、出力制御ステートメントブ
ロックは、図５に示すように、電磁放射パターン図の出
力可否や、電磁等値面図の出力可否等の各種の出力条件
を管理するものである。
に示すように、周波数や収束判定値といった各種の解析
条件を管理するものであり、出力制御ステートメントブ
ロックは、図５に示すように、電磁放射パターン図の出
力可否や、電磁等値面図の出力可否等の各種の出力条件
を管理するものである。
【００３２】一方、モデル記述ステートメントブロック
は、解析対象となる電気回路装置の構造体情報を管理す
るものである。この構造体情報は、分布定数線路近似法
でモデル化できる装置部分については、分布定数線路近
似法でモデル化され、分布定数線路近似法でモデル化で
きない装置部分については、モーメント法でモデル化さ
れることになる。図６に、プリント板の配線をセグメン
トに分割していくことでモデル化する分布定数線路近似
法のモデル化例、図７に、構造体をパッチに分割してい
くことでモデル化するモーメント法のモデル化例を図示
する。
は、解析対象となる電気回路装置の構造体情報を管理す
るものである。この構造体情報は、分布定数線路近似法
でモデル化できる装置部分については、分布定数線路近
似法でモデル化され、分布定数線路近似法でモデル化で
きない装置部分については、モーメント法でモデル化さ
れることになる。図６に、プリント板の配線をセグメン
トに分割していくことでモデル化する分布定数線路近似
法のモデル化例、図７に、構造体をパッチに分割してい
くことでモデル化するモーメント法のモデル化例を図示
する。
【００３３】図８に、このモデル記述ステートメントブ
ロックのデータ構造の詳細構成を図示する。この図に示
すように、モデル記述ステートメントブロックは、具体
的には、「＄POINT 」で、電気回路装置上の各ポイント
についてのポイント番号／座標値を管理し、「＄WIRE」
で、電気回路装置に実装される各ワイヤについてのワイ
ヤ番号／材料番号／配置ポイント番号を管理し、「＄PA
TCH 」で、電気回路装置上の各パッチについてのパッチ
番号／材料番号／配置ポイント番号を管理し、「＄MATE
RIAL」で、電気回路装置に実装される各材料の材料番号
／種別／電気特性／形状特性を管理し、「＄BOARD 」
で、電気回路装置に実装される各プリント板のボード番
号／材料番号／配置ポイント番号を管理し、「＄SEGMEN
T 」で、プリント板上の各セグメントについてのボード
番号／セグメント番号／配置座標値／プリント板の層番
号を管理し、「＄PIN 」で、プリント板上の各ピンにつ
いてのボード番号／ピン番号／配置座標値／プリント板
の層番号／回路番号／接続される波源番号を管理し、
「＄SOURCE」で、電気回路装置に実装される各波源の波
源番号／接続されるポイント番号／波源の電圧値を管理
し、「＄LOAD」で、電気回路装置に実装される各負荷の
負荷番号／接続されるポイント番号／負荷となるインピ
ーダンスを管理するというように、解析対象となる電気
回路装置の構造体情報を管理するのである。
ロックのデータ構造の詳細構成を図示する。この図に示
すように、モデル記述ステートメントブロックは、具体
的には、「＄POINT 」で、電気回路装置上の各ポイント
についてのポイント番号／座標値を管理し、「＄WIRE」
で、電気回路装置に実装される各ワイヤについてのワイ
ヤ番号／材料番号／配置ポイント番号を管理し、「＄PA
TCH 」で、電気回路装置上の各パッチについてのパッチ
番号／材料番号／配置ポイント番号を管理し、「＄MATE
RIAL」で、電気回路装置に実装される各材料の材料番号
／種別／電気特性／形状特性を管理し、「＄BOARD 」
で、電気回路装置に実装される各プリント板のボード番
号／材料番号／配置ポイント番号を管理し、「＄SEGMEN
T 」で、プリント板上の各セグメントについてのボード
番号／セグメント番号／配置座標値／プリント板の層番
号を管理し、「＄PIN 」で、プリント板上の各ピンにつ
いてのボード番号／ピン番号／配置座標値／プリント板
の層番号／回路番号／接続される波源番号を管理し、
「＄SOURCE」で、電気回路装置に実装される各波源の波
源番号／接続されるポイント番号／波源の電圧値を管理
し、「＄LOAD」で、電気回路装置に実装される各負荷の
負荷番号／接続されるポイント番号／負荷となるインピ
ーダンスを管理するというように、解析対象となる電気
回路装置の構造体情報を管理するのである。
【００３４】次に、図２に示す処理フローに従って、本
発明の電磁界強度算出装置１の実行する電磁界強度算出
処理について詳細に説明する。本発明の電磁界強度算出
装置１は、起動されると、図２の処理フローに示すよう
に、先ず最初に、ステップ１で、入力データファイル２
０から解析対象となる電気回路装置の構造体情報を読み
込んで初期設定を施す。
発明の電磁界強度算出装置１の実行する電磁界強度算出
処理について詳細に説明する。本発明の電磁界強度算出
装置１は、起動されると、図２の処理フローに示すよう
に、先ず最初に、ステップ１で、入力データファイル２
０から解析対象となる電気回路装置の構造体情報を読み
込んで初期設定を施す。
【００３５】次に、ステップ２で、処理済の周波数をカ
ウントすることで、解析制御ステートメントブロックに
登録されている全周波数についての処理を終了したのか
否かを判断して、終了したことを判断するときには全処
理を終了し、終了していないことを判断するときには、
以下の処理を行う周波数を未処理の中から１つ選択す
る。
ウントすることで、解析制御ステートメントブロックに
登録されている全周波数についての処理を終了したのか
否かを判断して、終了したことを判断するときには全処
理を終了し、終了していないことを判断するときには、
以下の処理を行う周波数を未処理の中から１つ選択す
る。
【００３６】続いて、ステップ３で、読み込んだ構造体
情報から、分布定数線路近似法を適用する装置部分であ
るのか、モーメント法を適用する装置部分であるのかを
判断する。この判断処理は、読み込んだ構造体情報の
内、上述の区分標識「＄WIRE」及び「＄PATCH 」の指す
装置部分についてはモーメント法を適用する装置部分で
あると判断し、一方、読み込んだ構造体情報の内、上述
の区分標識「＄BOARD 」、「＄SEGMENT 」及び「＄PIN
」の指す装置部分については分布定数線路近似法を適
用する装置部分であると判断していくことで行う。
情報から、分布定数線路近似法を適用する装置部分であ
るのか、モーメント法を適用する装置部分であるのかを
判断する。この判断処理は、読み込んだ構造体情報の
内、上述の区分標識「＄WIRE」及び「＄PATCH 」の指す
装置部分についてはモーメント法を適用する装置部分で
あると判断し、一方、読み込んだ構造体情報の内、上述
の区分標識「＄BOARD 」、「＄SEGMENT 」及び「＄PIN
」の指す装置部分については分布定数線路近似法を適
用する装置部分であると判断していくことで行う。
【００３７】ここで、このような区分標識による判断を
用いるのではなくて、構造体情報の形状や寸法等から直
接判断していく方法を採ることも可能である。すなわ
ち、信号配線層と接地配線層とを有するプリント板部分
ついては、分布定数線路近似法の適用可能な装置部分と
判断するとともに、図９に示すように、一定の径で、か
つ一定の間隔を持つ平行配線構造部分や、一定の幅で、
かつ一定の間隔を持つ平行平板構造部分や、平板と一定
の径の線分とが一定の間隔を持つ構造部分や、内部線分
が一定の径で、かつ一定の間隔を持つ外部円筒又は外部
多角柱からなる同軸構造部分については、分布定数線路
近似法の適用可能な装置部分と判断する。一方、筐体の
構造部分や、筐体の配線構造部分や、筐体から大地への
接地配線構造部分や、筐体に取り付けられたシールド構
造部分や、筐体の内部に実装される電気回路ユニットの
導体面構造部分といったような、それ以外の構造部分に
ついては、モーメント法を適用する装置部分と判断する
のである。なお、このとき、精度を優先させて、原則的
には分布定数線路近似法の適用可能な装置部分であって
も、モーメント法を適用する装置部分と判断していくこ
とがある。
用いるのではなくて、構造体情報の形状や寸法等から直
接判断していく方法を採ることも可能である。すなわ
ち、信号配線層と接地配線層とを有するプリント板部分
ついては、分布定数線路近似法の適用可能な装置部分と
判断するとともに、図９に示すように、一定の径で、か
つ一定の間隔を持つ平行配線構造部分や、一定の幅で、
かつ一定の間隔を持つ平行平板構造部分や、平板と一定
の径の線分とが一定の間隔を持つ構造部分や、内部線分
が一定の径で、かつ一定の間隔を持つ外部円筒又は外部
多角柱からなる同軸構造部分については、分布定数線路
近似法の適用可能な装置部分と判断する。一方、筐体の
構造部分や、筐体の配線構造部分や、筐体から大地への
接地配線構造部分や、筐体に取り付けられたシールド構
造部分や、筐体の内部に実装される電気回路ユニットの
導体面構造部分といったような、それ以外の構造部分に
ついては、モーメント法を適用する装置部分と判断する
のである。なお、このとき、精度を優先させて、原則的
には分布定数線路近似法の適用可能な装置部分であって
も、モーメント法を適用する装置部分と判断していくこ
とがある。
【００３８】ステップ３の処理を終了すると、分布定数
線路近似法を適用する装置部分の電流分布を算出すべ
く、ステップ４に進んで、ステップ３で分布定数線路近
似法を適用する装置部分であると判断されたｎ個の線路
に対して、最下流の線路の持つ負荷インピーダンス（上
述の「＄LOAD」の管理するもの）を使い、下流側の線路
から順番に各線路の入力インピーダンスＺk （ｋ＝ｎ〜
１）を計算する。
線路近似法を適用する装置部分の電流分布を算出すべ
く、ステップ４に進んで、ステップ３で分布定数線路近
似法を適用する装置部分であると判断されたｎ個の線路
に対して、最下流の線路の持つ負荷インピーダンス（上
述の「＄LOAD」の管理するもの）を使い、下流側の線路
から順番に各線路の入力インピーダンスＺk （ｋ＝ｎ〜
１）を計算する。
【００３９】この入力インピーダンスＺk （≡Ｚin）の
計算は、図１０に示すように、線路の特性インピーダン
スをＺ0 、伝播定数をγ、線路長をｌ、下流側の負荷イ
ンピーダンスをＺL で表すならば、分布定数線路近似法
に従い、
計算は、図１０に示すように、線路の特性インピーダン
スをＺ0 、伝播定数をγ、線路長をｌ、下流側の負荷イ
ンピーダンスをＺL で表すならば、分布定数線路近似法
に従い、
【００４０】
【数１】
【００４１】により算出され、そして、この線路の上流
に位置する線路の入力インピーダンスＺinは、算出され
た下流側線路の入力インピーダンスＺinを新たにＺL と
して用いることで算出されていくことになる。なお、図
１１に示すように線路が分岐する場合には、オームの法
則に従い、
に位置する線路の入力インピーダンスＺinは、算出され
た下流側線路の入力インピーダンスＺinを新たにＺL と
して用いることで算出されていくことになる。なお、図
１１に示すように線路が分岐する場合には、オームの法
則に従い、
【００４２】
【数２】
【００４３】により算出されていく。続いて、ステップ
５で、ステップ３で分布定数線路近似法を適用する装置
部分であると判断されたｎ個の線路に対して、最上流の
線路の持つ波源（上述の「＄SOURCE」の管理するもの）
と、算出された入力インピーダンスＺk とを用いて、上
流側から順番に各線路の電流Ｉk （ｋ＝１〜ｎ）を計算
する。
５で、ステップ３で分布定数線路近似法を適用する装置
部分であると判断されたｎ個の線路に対して、最上流の
線路の持つ波源（上述の「＄SOURCE」の管理するもの）
と、算出された入力インピーダンスＺk とを用いて、上
流側から順番に各線路の電流Ｉk （ｋ＝１〜ｎ）を計算
する。
【００４４】この電流Ｉk の計算は、図１２に示すよう
に、線路の特性インピーダンスをＺ 0 、伝播定数をγ、
線路長をｌ、線路の上流側インピーダンスをＺinA 、線
路の下流側インピーダンスをＺinB 、線路の上流側電流
をＩA 、線路の下流側電流をＩB 、線路の上流側波源を
ＶA （＝ＺinA ×ＩA ）で表すならば、分布定数線路近
似法に従い、
に、線路の特性インピーダンスをＺ 0 、伝播定数をγ、
線路長をｌ、線路の上流側インピーダンスをＺinA 、線
路の下流側インピーダンスをＺinB 、線路の上流側電流
をＩA 、線路の下流側電流をＩB 、線路の上流側波源を
ＶA （＝ＺinA ×ＩA ）で表すならば、分布定数線路近
似法に従い、
【００４５】
【数３】
【００４６】により算出されていくことになる。一方、
ステップ３の処理を終了すると、モーメント法を適用す
る装置部分の電流分布を算出すべく、ステップ６に進ん
で、グリーン関数を使って、ステップ３でモーメント法
を適用する装置部分であると判断されたｍ個の装置部分
間の相互インピーダンスＺij（ｉ＝１〜ｍ、ｊ＝ｉ〜
ｍ）を計算するとともに、ステップ３でモーメント法を
適用する装置部分であると判断されたｍ個の装置部分
と、ステップ３で分布定数線路近似法を適用する装置部
分であると判断されたｎ個の線路との間の相互インピー
ダンスＺik（ｉ＝１〜ｍ、ｋ＝１〜ｎ）を計算する。す
なわち、図１３に示すような相互インピーダンスＺij，
Ｚikとを計算するのである。
ステップ３の処理を終了すると、モーメント法を適用す
る装置部分の電流分布を算出すべく、ステップ６に進ん
で、グリーン関数を使って、ステップ３でモーメント法
を適用する装置部分であると判断されたｍ個の装置部分
間の相互インピーダンスＺij（ｉ＝１〜ｍ、ｊ＝ｉ〜
ｍ）を計算するとともに、ステップ３でモーメント法を
適用する装置部分であると判断されたｍ個の装置部分
と、ステップ３で分布定数線路近似法を適用する装置部
分であると判断されたｎ個の線路との間の相互インピー
ダンスＺik（ｉ＝１〜ｍ、ｋ＝１〜ｎ）を計算する。す
なわち、図１３に示すような相互インピーダンスＺij，
Ｚikとを計算するのである。
【００４７】続いて、ステップ７で、ステップ６で算出
した相互インピーダンスＺij，Ｚikと、上述の「＄SOUR
CE」で指定されるモーメント法の波源である電圧値Ｖi
（ｉ＝１〜ｍ）と、ステップ５で算出した分布定数線路
近似法を適用する装置部分の電流Ｉk とを使い、図１４
に示すモーメント法の積分方程式を解くことで、モーメ
ント法を適用する装置部分の電流Ｉi （ｉ＝１〜ｍ）を
算出する。
した相互インピーダンスＺij，Ｚikと、上述の「＄SOUR
CE」で指定されるモーメント法の波源である電圧値Ｖi
（ｉ＝１〜ｍ）と、ステップ５で算出した分布定数線路
近似法を適用する装置部分の電流Ｉk とを使い、図１４
に示すモーメント法の積分方程式を解くことで、モーメ
ント法を適用する装置部分の電流Ｉi （ｉ＝１〜ｍ）を
算出する。
【００４８】ここで、この算出処理にあって、相互イン
ピーダンスＺikは、ステップ５で算出した電流Ｉk との
マトリクス演算のみに使用されて、電流Ｉk がモーメン
ト法の適用装置部分に誘起する波源を算出するためにの
み用いられる。すなわち、従来のように、解析対象とな
る電気回路装置全体にモーメント法を適用する場合に
は、図１５（ａ）に示すように、「（ｎ＋ｍ）×（ｎ＋
ｍ）」個の相互インピーダンスを用いてモーメント法の
積分方程式を解くのに対して、本発明では、図１５
（ｂ）に示すように、「ｍ×ｍ」の相互インピーダンス
Ｚijのみを用いてモーメント法の積分方程式を解くこと
になるのである。従って、高速に積分方程式を解くこと
ができる。なお、相互インピーダンスは対称行列である
ことから、図１５の破線に示すように実際には計算量は
半分となる。
ピーダンスＺikは、ステップ５で算出した電流Ｉk との
マトリクス演算のみに使用されて、電流Ｉk がモーメン
ト法の適用装置部分に誘起する波源を算出するためにの
み用いられる。すなわち、従来のように、解析対象とな
る電気回路装置全体にモーメント法を適用する場合に
は、図１５（ａ）に示すように、「（ｎ＋ｍ）×（ｎ＋
ｍ）」個の相互インピーダンスを用いてモーメント法の
積分方程式を解くのに対して、本発明では、図１５
（ｂ）に示すように、「ｍ×ｍ」の相互インピーダンス
Ｚijのみを用いてモーメント法の積分方程式を解くこと
になるのである。従って、高速に積分方程式を解くこと
ができる。なお、相互インピーダンスは対称行列である
ことから、図１５の破線に示すように実際には計算量は
半分となる。
【００４９】このようにして、ステップ５で分布定数線
路近似法を適用する装置部分の電流Ｉk が算出され、ス
テップ７でモーメント法を適用する装置部分の電流Ｉi
が算出されると、続いて、ステップ８で、処理済の観測
点をカウントすることで、出力制御ステートメントブロ
ックに登録されている全観測点についての処理を終了し
たのか否かを判断して、終了したことを判断するときに
はステップ２に戻り、終了していないことを判断すると
きには、以下の処理を行う観測点を未処理の中から１つ
選択する。
路近似法を適用する装置部分の電流Ｉk が算出され、ス
テップ７でモーメント法を適用する装置部分の電流Ｉi
が算出されると、続いて、ステップ８で、処理済の観測
点をカウントすることで、出力制御ステートメントブロ
ックに登録されている全観測点についての処理を終了し
たのか否かを判断して、終了したことを判断するときに
はステップ２に戻り、終了していないことを判断すると
きには、以下の処理を行う観測点を未処理の中から１つ
選択する。
【００５０】続いて、ステップ９で、ステップ５で算出
された電流Ｉk と、ステップ７で算出された電流Ｉi と
を使って、それらの電流Ｉk,Ｉi が観測点にもたらす電
磁界強度を算出して、その算出結果を出力データファイ
ル２１に格納してからステップ８に戻っていく。
された電流Ｉk と、ステップ７で算出された電流Ｉi と
を使って、それらの電流Ｉk,Ｉi が観測点にもたらす電
磁界強度を算出して、その算出結果を出力データファイ
ル２１に格納してからステップ８に戻っていく。
【００５１】このようにして、本発明では、図１６に示
すように、解析対象となる電気回路装置が、モーメント
法を適用する装置部分である２個のパッチと、分布定数
線路近似法を適用する装置部分である３個のセグメント
とからなる場合には、図１７に示すように、パッチとセ
グメントとの間の相互インピーダンスと、ステップ５の
処理で求まるセグメントの電流値とのマトリクス演算に
より、それらの電流値が誘起する波源を算出して、それ
をパッチの波源に加算しつつ、積分方程式を解くことで
パッチの電流値を求める構成を採るものである。
すように、解析対象となる電気回路装置が、モーメント
法を適用する装置部分である２個のパッチと、分布定数
線路近似法を適用する装置部分である３個のセグメント
とからなる場合には、図１７に示すように、パッチとセ
グメントとの間の相互インピーダンスと、ステップ５の
処理で求まるセグメントの電流値とのマトリクス演算に
より、それらの電流値が誘起する波源を算出して、それ
をパッチの波源に加算しつつ、積分方程式を解くことで
パッチの電流値を求める構成を採るものである。
【００５２】
【発明の効果】以上説明したように、本発明の電磁界強
度算出装置は、分布定数線路近似法の適用可能な装置部
分については、適用不可能な装置部分の持つ電流値が小
さいことに対応して、その電流値の誘起する波源を考慮
することなく、計算の簡単な分布定数線路近似法に従っ
て電流分布を求め、一方、分布定数線路近似法の適用不
可能な装置部分については、適用可能な装置部分の持つ
電流値が大きいことに対応して、その電流値の誘起する
波源を考慮しつつ、高精度のモーメント法に従って電流
分布を求める構成を採るものであることから、高い精度
で、高速に、しかも少ない記憶資源でもって、複雑な構
造体の電気回路装置の放射する電磁界強度を算出できる
ようになる。
度算出装置は、分布定数線路近似法の適用可能な装置部
分については、適用不可能な装置部分の持つ電流値が小
さいことに対応して、その電流値の誘起する波源を考慮
することなく、計算の簡単な分布定数線路近似法に従っ
て電流分布を求め、一方、分布定数線路近似法の適用不
可能な装置部分については、適用可能な装置部分の持つ
電流値が大きいことに対応して、その電流値の誘起する
波源を考慮しつつ、高精度のモーメント法に従って電流
分布を求める構成を採るものであることから、高い精度
で、高速に、しかも少ない記憶資源でもって、複雑な構
造体の電気回路装置の放射する電磁界強度を算出できる
ようになる。
【図１】本発明の原理構成図である。
【図２】本発明の実行する処理フローの一実施例であ
る。
る。
【図３】入力データファイルの管理データのデータ構造
例である。
例である。
【図４】解析制御ステートメントブロックのデータ構造
例である。
例である。
【図５】出力制御ステートメントブロックのデータ構造
例である。
例である。
【図６】分布定数線路近似法のモデル化例である。
【図７】モーメント法のモデル化例である。
【図８】モデル記述ステートメントブロックのデータ構
造例である。
造例である。
【図９】分布定数線路近似法の適用可能な構造例であ
る。
る。
【図１０】入力インピーダンスの算出処理の説明図であ
る。
る。
【図１１】入力インピーダンスの算出処理の説明図であ
る。
る。
【図１２】電流算出処理の説明図である。
【図１３】相互インピーダンスの説明図である。
【図１４】本発明の積分方程式の説明図である。
【図１５】積分方程式に用いられる相互インピーダンス
の説明図である。
の説明図である。
【図１６】電気回路装置の構造例である。
【図１７】本発明の処理説明図である。
【図１８】微小ループアンテナ近似法の説明図である。
【図１９】分布定数線路近似法の説明図である。
【図２０】モーメント法が適用される構造体の一例であ
る。
る。
１ 電磁界強度算出装置 １０ 入力手段 １１ 区分手段 １２ 第１の算出手段 １３ 第２の算出手段 １４ 第３の算出手段 １５ 出力手段 １３０ 相互インピーダンス算出手段 １３１ 誘起波源算出手段 １３２ 電流分布算出手段
───────────────────────────────────────────────────── フロントページの続き (72)発明者 長瀬 健二 神奈川県川崎市中原区上小田中1015番地 富士通株式会社内
Claims (13)
- 【請求項１】 電気回路装置の放射する電磁界強度を算
出する電磁界強度算出装置において、 解析対象となる電気回路装置の各装置部分を、分布定数
線路近似法の適用可能な装置部分と、適用不可能な装置
部分とに区分けする区分手段(11)と、 上記区分手段(11)が適用可能と区分けする装置部分に対
して分布定数線路近似法を施すことで、該装置部分の電
流分布を算出する第１の算出手段(12)と、 上記第１の算出手段(12)の算出する電流分布を用いつ
つ、上記区分手段(11)が適用不可能と区分けする装置部
分に対してモーメント法を施すことで、該装置部分の電
流分布を算出する第２の算出手段(13)と、 上記第１の算出手段(12)の算出する電流分布と、上記第
２の算出手段(13)の算出する電流分布とから、電気回路
装置の放射する電磁界強度を算出する第３の算出手段(1
4)とを備えることを、 特徴とする電磁界強度算出装置。 - 【請求項２】 請求項１記載の電磁界強度算出装置にお
いて、 第２の算出手段(13)は、 区分手段(11)が適用不可能と判断する装置部分間の相互
インピーダンスを算出するとともに、該装置部分と、区
分手段(11)が適用可能と判断する装置部分との間の相互
インピーダンスを算出する相互インピーダンス算出手段
(130) と、 第１の算出手段(12)の算出する電流分布と、上記相互イ
ンピーダンス算出手段(130) の算出する異種間の相互イ
ンピーダンスとのマトリクス演算から、該電流分布が適
用不可能装置部分に誘起する波源を算出する誘起波源算
出手段(131) と、 上記誘起波源算出手段(131) の算出する波源を考慮しつ
つ、上記相互インピーダンス算出手段(130) の算出する
同種間の相互インピーダンスを使い、モーメント法に従
って区分手段(11)が適用不可能と判断する装置部分の電
流分布を算出する電流分布算出手段(132) とから構成さ
れることを、 特徴とする電磁界強度算出装置。 - 【請求項３】 請求項２記載の電磁界強度算出装置にお
いて、 電流分布算出手段(132) は、誘起波源算出手段(131) の
算出する波源と、区分手段(11)が適用不可能と判断する
装置部分の持つ本来の波源との合成値を、モーメント法
に用いる波源として使用していくよう処理することを、 特徴とする電磁界強度算出装置。 - 【請求項４】 請求項１ないし３記載の電磁界強度算出
装置において、 区分手段(11)は、電気回路装置の構造情報を分析するこ
とで、分布定数線路近似法の適用可能な装置部分と、適
用不可能な装置部分とを区分けしていくよう処理するこ
とを、 特徴とする電磁界強度算出装置。 - 【請求項５】 請求項１ないし３記載の電磁界強度算出
装置において、 電気回路装置の構造情報が区分標識を使って表現され
て、区分手段(11)に与えられるよう構成されることを、 特徴とする電磁界強度算出装置。 - 【請求項６】 請求項５記載の電磁界強度算出装置にお
いて、 区分手段(11)は、区分標識に従って、分布定数線路近似
法の適用可能な装置部分と、適用不可能な装置部分とを
区分けしていくよう処理することを、 特徴とする電磁界強度算出装置。 - 【請求項７】 請求項１ないし６記載の電磁界強度算出
装置において、 区分手段(11)は、１次元の構造物として近似できる装置
部分を、分布定数線路近似法の適用可能な装置部分と判
断していくよう処理することを、 特徴とする電磁界強度算出装置。 - 【請求項８】 請求項１ないし７記載の電磁界強度算出
装置において、 区分手段(11)は、信号配線層と接地配線層とを有するプ
リント板部分を、分布定数線路近似法の適用可能な装置
部分と判断していくよう処理することを、 特徴とする電磁界強度算出装置。 - 【請求項９】 請求項１ないし７記載の電磁界強度算出
装置において、 区分手段(11)は、同軸構造を有するケーブル部分を、分
布定数線路近似法の適用可能な装置部分と判断していく
よう処理することを、 特徴とする電磁界強度算出装置。 - 【請求項１０】 請求項１ないし７記載の電磁界強度算
出装置において、 区分手段(11)は、平行配線構造部分を、分布定数線路近
似法の適用可能な装置部分と判断していくよう処理する
ことを、 特徴とする電磁界強度算出装置。 - 【請求項１１】 請求項１ないし７記載の電磁界強度算
出装置において、 区分手段(11)は、平行平板構造部分を、分布定数線路近
似法の適用可能な装置部分と判断していくよう処理する
ことを、 特徴とする電磁界強度算出装置。 - 【請求項１２】 請求項１ないし７記載の電磁界強度算
出装置において、 区分手段(11)は、平板と配線とが平行に配設される構造
部分を、分布定数線路近似法の適用可能な装置部分と判
断していくよう処理することを、 特徴とする電磁界強度算出装置。 - 【請求項１３】 請求項１ないし１２記載の電磁界強度
算出装置において、 区分手段(11)は、原則的には分布定数線路近似法の適用
可能な装置部分であっても、分布定数線路近似法の適用
不可能な装置部分と判断していくよう処理することを、 特徴とする電磁界強度算出装置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2710994A JP3091815B2 (ja) | 1994-02-25 | 1994-02-25 | 電磁界強度算出装置 |
| US08/351,213 US5650935A (en) | 1994-02-25 | 1994-11-30 | Apparatus for calculating intensity of an electromagnetic field |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2710994A JP3091815B2 (ja) | 1994-02-25 | 1994-02-25 | 電磁界強度算出装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH07234890A true JPH07234890A (ja) | 1995-09-05 |
| JP3091815B2 JP3091815B2 (ja) | 2000-09-25 |
Family
ID=12211920
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2710994A Expired - Lifetime JP3091815B2 (ja) | 1994-02-25 | 1994-02-25 | 電磁界強度算出装置 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US5650935A (ja) |
| JP (1) | JP3091815B2 (ja) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0778533A2 (en) | 1995-11-16 | 1997-06-11 | Fujitsu Limited | Electromagnetic fed intensity calculation apparatus and method |
| WO1999012110A1 (fr) * | 1997-09-02 | 1999-03-11 | Hitachi Chemical Company, Ltd. | Procede et systeme d'approximation d'une ligne a constantes reparties |
| US6285973B1 (en) | 1997-11-26 | 2001-09-04 | Fujitsu Limited | Apparatus for, and method of, calculating electromagnetic field strength and medium for recording program to achieve the method |
| US6766261B2 (en) | 2001-03-27 | 2004-07-20 | Fujitsu Limited | Method and device generating data input to an electromagnetic field intensity calculating device |
| DE19710687B4 (de) * | 1996-04-10 | 2006-02-16 | Fujitsu Ltd., Kawasaki | Simulationsvorrichtung und Simulationsverfahren für die Intensität eines elektromagnetischen Feldes mit Hilfe der Momentenmethode |
| CN105929253A (zh) * | 2016-06-15 | 2016-09-07 | 北京航空航天大学 | 采用偶极矩模型对pcb电路电磁场进行反演的方法 |
Families Citing this family (20)
| 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 |
| EP0897731A4 (en) * | 1996-11-26 | 2003-07-30 | Mitsubishi Electric Corp | METHOD FOR CHANGING THE POWER DISTRIBUTION |
| 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 | 多層基板からの放射電磁波解析装置 |
| JPH11296499A (ja) | 1998-04-07 | 1999-10-29 | Fujitsu Ltd | モーメント法を用いたシミュレーション装置及び方法並びにプログラム記録媒体 |
| US6263288B1 (en) | 1998-07-22 | 2001-07-17 | Eastman Kodak Company | Method and apparatus for indicating proximity of film scanner to CRT display monitor |
| JP3394202B2 (ja) * | 1999-01-13 | 2003-04-07 | 太陽誘電株式会社 | 電磁界強度の測定方法及び装置並びに電流電圧分布の測定方法及び装置 |
| JP2000214191A (ja) | 1999-01-27 | 2000-08-04 | Fujitsu Ltd | シミュレ―ション装置及び方法並びにプログラム記録媒体 |
| JP2002092414A (ja) * | 2000-09-18 | 2002-03-29 | Saga Sanyo Industries Co Ltd | 電子部品の機種選択方法および電子部品メーカ・サーバ |
| JP3943326B2 (ja) * | 2000-11-27 | 2007-07-11 | 松下電器産業株式会社 | 不要輻射解析方法および不要輻射解析装置 |
| US20020082812A1 (en) * | 2000-12-26 | 2002-06-27 | Fujitsu Limited | Apparatus and method for simulating the receiving characteristic of radio waves |
| JP4003403B2 (ja) * | 2001-03-16 | 2007-11-07 | 株式会社日立製作所 | 電磁波発生源探査方法 |
| JP2003296395A (ja) * | 2002-04-01 | 2003-10-17 | Fujitsu Ltd | 相互インミッタンス計算プログラム、相互インミッタンス計算装置、相互インミッタンス計算方法および電磁界強度計算プログラム |
| US7222033B1 (en) | 2003-08-18 | 2007-05-22 | Steven Lynn Newson | Electromagnetic emissions and susceptibility calculating method and apparatus |
| WO2006103741A1 (ja) * | 2005-03-28 | 2006-10-05 | Fujitsu Limited | 電磁界強度算出方法、電磁界強度算出装置、制御プログラム |
| US8026731B1 (en) * | 2006-04-20 | 2011-09-27 | Jack Scott Emery | Method and apparatus for constructing images from measurements of impedance |
| CN110705172A (zh) * | 2019-10-23 | 2020-01-17 | 西安交通大学 | 雷电冲击电压下变压器瞬时电场分布的计算方法 |
| JP7424133B2 (ja) * | 2020-03-17 | 2024-01-30 | 富士通株式会社 | 放射電磁波推定プログラム，情報処理装置および放射電磁波推定方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
-
1994
- 1994-02-25 JP JP2710994A patent/JP3091815B2/ja not_active Expired - Lifetime
- 1994-11-30 US US08/351,213 patent/US5650935A/en not_active Expired - Lifetime
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0778533A2 (en) | 1995-11-16 | 1997-06-11 | Fujitsu Limited | Electromagnetic fed intensity calculation apparatus and method |
| DE19710687B4 (de) * | 1996-04-10 | 2006-02-16 | Fujitsu Ltd., Kawasaki | Simulationsvorrichtung und Simulationsverfahren für die Intensität eines elektromagnetischen Feldes mit Hilfe der Momentenmethode |
| WO1999012110A1 (fr) * | 1997-09-02 | 1999-03-11 | Hitachi Chemical Company, Ltd. | Procede et systeme d'approximation d'une ligne a constantes reparties |
| US6285973B1 (en) | 1997-11-26 | 2001-09-04 | Fujitsu Limited | Apparatus for, and method of, calculating electromagnetic field strength and medium for recording program to achieve the method |
| DE19826236B4 (de) * | 1997-11-26 | 2007-09-20 | Fujitsu Ltd., Kawasaki | Vorrichtung und Verfahren zum Berechnen der Elektromagnetischen Feldstärke |
| US6766261B2 (en) | 2001-03-27 | 2004-07-20 | Fujitsu Limited | Method and device generating data input to an electromagnetic field intensity calculating device |
| CN105929253A (zh) * | 2016-06-15 | 2016-09-07 | 北京航空航天大学 | 采用偶极矩模型对pcb电路电磁场进行反演的方法 |
| CN105929253B (zh) * | 2016-06-15 | 2018-10-02 | 北京航空航天大学 | 采用偶极矩模型对pcb电路电磁场进行反演的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| US5650935A (en) | 1997-07-22 |
| JP3091815B2 (ja) | 2000-09-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH07234890A (ja) | 電磁界強度算出装置 | |
| Hockanson et al. | Quantifying EMI resulting from finite-impedance reference planes | |
| Luebbers et al. | A simple feed model that reduces time steps needed for FDTD antenna and microstrip calculations | |
| JP2768900B2 (ja) | 電磁界強度算出装置 | |
| US5812434A (en) | Electromagnetic field strength calculator having function of displaying currents to be analyzed | |
| US6185517B1 (en) | Electromagnetic field intensity computing device | |
| Gay‐Balmaz et al. | Three‐dimensional planar radiating structures in stratified media | |
| US6879942B1 (en) | Apparatus for calculating immunity from radiated electromagnetic field, method for achieving calculation, and storage medium storing programs therefor | |
| US6129459A (en) | Apparatus for analyzing radiating electromagnetic wave from multilayer substrate | |
| US6766261B2 (en) | Method and device generating data input to an electromagnetic field intensity calculating device | |
| Li et al. | FD-TD analysis of electromagnetic radiation from modules-on-backplane configurations | |
| JP3980811B2 (ja) | 電磁界強度算出装置、電磁界強度算出方法およびその方法をコンピュータに実行させるプログラムを記録したコンピュータ読み取り可能な記録媒体 | |
| JP3825133B2 (ja) | 電磁界強度算出装置及び方法並びにプログラム記録媒体 | |
| JP3443162B2 (ja) | 電磁界強度算出装置 | |
| US11378605B2 (en) | Method for high-intensity radiated field (HIRF) and electromagnetic pulse (EMP) analysis of a vehicle | |
| Toscano et al. | A new efficient moment method formulation for the design of microstrip antennas over a chiral grounded slab | |
| JP3703812B2 (ja) | Ｆｄｔｄ法を用いた電磁界解析方法、電磁界解析における媒質表現方法、シミュレーション装置、及びプログラム | |
| Paul et al. | Fast finite difference time domain method for the analysis of planar microstrip circuits | |
| JP3887365B2 (ja) | 放射電磁界耐性算出装置及びプログラム記録媒体 | |
| JP3363657B2 (ja) | 電磁界強度算出装置 | |
| JPH08304494A (ja) | 電磁界強度算出装置 | |
| Singer et al. | State of the art in the method of moments | |
| Mocker et al. | Comparison of electromagnetic solvers for antennas mounted on vehicles | |
| Liu et al. | Sphere-TEM mode in rectangular top horn offset coaxial transmission line | |
| Zhao et al. | Efficient solution of electromagnetic problems based on integral equation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20000711 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080721 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090721 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100721 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100721 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110721 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110721 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120721 Year of fee payment: 12 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120721 Year of fee payment: 12 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130721 Year of fee payment: 13 |
|
| EXPY | Cancellation because of completion of term |