JPH07302277A - 電磁界強度算出装置 - Google Patents

電磁界強度算出装置

Info

Publication number
JPH07302277A
JPH07302277A JP9536294A JP9536294A JPH07302277A JP H07302277 A JPH07302277 A JP H07302277A JP 9536294 A JP9536294 A JP 9536294A JP 9536294 A JP9536294 A JP 9536294A JP H07302277 A JPH07302277 A JP H07302277A
Authority
JP
Japan
Prior art keywords
electromagnetic field
calculation
mutual
field strength
calculating
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.)
Pending
Application number
JP9536294A
Other languages
English (en)
Inventor
Kenji Nagase
健二 長瀬
Makoto Mukai
誠 向井
Shinichi Otsu
信一 大津
Kanji Nishino
関司 西野
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP9536294A priority Critical patent/JPH07302277A/ja
Publication of JPH07302277A publication Critical patent/JPH07302277A/ja
Pending legal-status Critical Current

Links

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30—Computing systems specially adapted for manufacturing

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

(57)【要約】 【目的】本発明は、モーメント法に基づいて電気回路装
置の放射する電磁界強度を算出する電磁界強度算出装置
に関し、高速かつ高精度に電磁界強度を算出できるよう
にすることを目的とする。 【構成】モーメント法の連立方程式を構成する相互イン
ピーダンスの算出にあたって、相互インピーダンスの算
出対象となる要素の電気長と、要素間距離の電気長とを
評価する評価手段１１２と、評価手段１１２の評価結果
に応じて、相互インピーダンスの算出手法を決定する決
定手段１１３とを備え、決定手段１１３の決定する算出
手法に従って相互インピーダンスを算出することで、モ
ーメント法の連立方程式を導出していくように構成す
る。

Description

【発明の詳細な説明】
【０００１】
【産業上の利用分野】本発明は、モーメント法に基づい
て電気回路装置の放射する電磁界強度を算出する電磁界
強度算出装置に関し、特に、高速かつ高精度に電磁界強
度を算出できるようにする電磁界強度算出装置に関す
る。
【０００２】電気回路装置から不要放射される電波は、
テレビやラジオ等の他の電波と干渉することから、最
近、各国で厳しく規制されるようになってきた。このよ
うな規制の規格として、日本ではＶＣＣＩ規格があり、
米国ではＦＣＣ規格があり、ドイツではＶＤＥ規格があ
る。
【０００３】このような電波規制を満足させるために
は、シールド技術やフィルタ技術等のような種々の対策
技術を使う必要があり、これらの対策技術がどの程度電
波を減少させるかを定量的にシミュレートしていく必要
がある。これから、高い精度でもって電気回路装置の放
射する電磁界強度をシミュレートできるようにする電磁
界強度算出装置の構築が叫ばれている。
【０００４】
【従来の技術】任意形状の物体から放射される電磁界強
度は、物体各部に流れる電流が分かれば、公知の理論式
を用いて容易に計算することができる。この電流値は、
理論的には、マックスウェルの電磁波動方程式を与えら
れた境界条件の下で解くことで得られるものであるが、
任意形状物体を対象とした複雑な境界条件の下での直接
的な数式解は現在知られていない。
【０００５】従って、現在の電磁界強度算出装置で用い
られている電流を求める解法は、難易の程度はあるにせ
よ、すべて近似的なものである。現在、この近似的な解
法として、微小ループアンテナ近似法と、分布定数線路
近似法と、モーメント法という３種類のものが知られて
いる。
【０００６】微小ループアンテナ近似法は、波源回路と
負荷回路との間を接続する配線をループアンテナとして
扱い、ループ上の電流は平坦なものと仮定して、これを
集中定数回路の計算手法で求める方法である。図１１
に、この微小ループアンテナ近似法の構成を図示する。
【０００７】この微小ループアンテナ近似法による計算
は、最も簡単であるが、ループの寸法が電磁波の波長に
比べて無視できない条件では精度が極めて低下すること
から、現実には殆ど用いられていない。
【０００８】一方、分布定数線路近似法は、１次元の構
造物として近似できる物体に対して、分布定数線路の方
程式を適用することで電流を求める方法である。計算は
比較的簡単であり、計算時間及び記憶容量も解析要素数
にほぼ比例して増加するに止まり、線路の反射や共振等
の現象も含めた解析ができることから、１次元の近似が
成立する物体については、高速・高精度の解析ができ
る。図１２に、この分布定数線路近似法の構成を図示す
る。
【０００９】この分布定数線路近似法による計算は、１
次元の構造物として近似できる物体については、高速・
高精度の解析ができるものの、近似できない物体につい
ては解析できないという問題点がある。
【００１０】一方、モーメント法は、マックスウェルの
電磁波動方程式から導かれる積分方程式の解法の１つで
あり、３次元の任意形状物体を扱うことができる。具体
的には、物体を小さな要素に分割して電流の計算を行う
ものである。
【００１１】このように、モーメント法が３次元の任意
形状物体を扱えることから、電磁界強度算出装置では、
モーメント法を使って、電気回路装置の放射する電磁界
強度を算出するという構成が有力である。
【００１２】このモーメント法では、金属対象物を扱う
ときには、金属部分を解析対象としてメッシュ化し、分
割した金属間の相互インピーダンスＺijを求め、この相
互インピーダンスＺijと、波源Ｖi と、分割した金属に
流れる電流Ｉi との間に成立するモーメント法の連立方
程式、 〔Ｚij〕〔Ｉi 〕＝〔Ｖi 〕 を解いて電流Ｉi を求めていくことで、放射される電磁
界強度を算出するという方法を採っている。ここで、
「〔 〕」はマトリクスを表している。
【００１３】このようなモーメント法を用いるときにあ
って、従来の電磁界強度算出装置では、高速演算を実現
するために、この相互インピーダンスＺijを倍精度実数
の算出手法を使って算出していくという構成を採ってい
る。
【００１４】
【発明が解決しようとする課題】しかしながら、倍精度
実数の算出手法では、積演算のときに下位の有効数字が
失われることから、倍精度実数の算出手法を使って相互
インピーダンスＺijを算出していくと、金属メッシュの
電気長（放射される電磁波の波長を尺度する長さ）が短
くなるときに桁落ちが発生することで、相互インピーダ
ンスＺijを正確に算出できないという問題点があった。
【００１５】これから、従来の電磁界強度算出装置で
は、高い精度でもって電磁界強度を算出することができ
ないという問題点があった。本発明はかかる事情に鑑み
てなされたものであって、電気回路装置の放射する電磁
界強度を算出する構成を採るときにあって、高速かつ高
精度に電磁界強度を算出できるようにする新たな電磁界
強度算出装置の提供を目的とする。
【００１６】
【課題を解決するための手段】図１に本発明の原理構成
を図示する。図中、１は本発明により構成される電磁界
強度算出装置であって、モーメント法に基づいて解析対
象となる電気回路装置の放射する電磁界強度を算出す
る。
【００１７】この電磁界強度算出装置１は、解析対象と
なる電気回路装置の構造体を入力する入力手段１０と、
モーメント法に基づいて解析対象となる電気回路装置の
放射する電磁界強度を算出する電磁界強度算出手段１１
と、電磁界強度算出手段１１の算出結果を出力する出力
手段１２とを備える。
【００１８】電磁界強度算出手段１１は、モーメント法
に基づいて解析対象となる電気回路装置の放射する電磁
界強度を算出するために、分割手段１１０と、導出手段
１１１と、評価手段１１２と、決定手段１１３と、算出
手段１１４と、計算手段１１５とを備える。
【００１９】この分割手段１１０は、入力された電気回
路装置の構造体を要素（メッシュ）に分割する。導出手
段１１１は、金属のみを解析対象とする場合には、要素
間の相互インピーダンスのみを相互係数として算出し、
金属と誘電体とを解析対象とする場合には、要素間の相
互インピーダンス，相互アドミッタンス，相互リアクシ
ョンを相互係数として算出して、これらの相互係数を使
ってモーメント法の連立方程式を導出する。
【００２０】評価手段１１２は、導出手段１１１に展開
されて、相互係数の算出対象となる要素の電気長と、要
素間距離の電気長とを評価する。決定手段１１３は、導
出手段１１１に展開されて、評価手段１１２の評価結果
に応じて、相互係数の算出手法を決定する。
【００２１】算出手段１１４は、導出されたモーメント
法の連立方程式を解くことで、金属のみを解析対象とす
る場合には、金属に流れる電流を算出し、金属と誘電体
とを解析対象とする場合には、金属に流れる電流と、誘
電体の表面に流れる等価電流及び等価磁流とを算出す
る。計算手段１１５は、算出手段１１４の算出する算出
値から、電気回路装置の放射する電磁界強度を計算す
る。
【００２２】
【作用】本発明では、入力手段１０により解析対象とな
る電気回路装置の構造体が入力されると、先ず最初に、
電磁界強度算出手段１１の分割手段１１０は、その電気
回路装置を要素に分割する。
【００２３】続いて、評価手段１１２は、モーメント法
の連立方程式を構成する相互係数の算出対象となる要素
の電気長と、要素間距離の電気長とを評価し、この評価
結果を受けて、決定手段１１３は、要素の電気長が規定
値よりも短くなるときにあって、その電気長と比較して
要素間距離の電気長が長いことが評価されるときには、
高精度の算出手法（例えば、多倍精度実数や多倍長整数
の算出手法）を用いることを決定し、そうでないことが
評価されるときには、高速の算出手法（例えば、倍精度
実数の算出手法）を用いることを決定する。
【００２４】この決定手段１１３の決定結果を受けて、
導出手段１１１は、決定された算出手法を使って、金属
のみを解析対象とする場合には、要素間の相互係数とし
て、要素間の相互インピーダンスを算出し、金属と誘電
体とを解析対象とする場合には、要素間の相互係数とし
て、要素間の相互インピーダンス，相互アドミッタン
ス，相互リアクションを算出する。
【００２５】続いて、導出手段１１１は、算出した相互
係数を用いて、金属のみを解析対象とする場合には、金
属に流れる電流を未知数とするモーメント法の連立方程
式を導出し、金属と誘電体とを解析対象とする場合に
は、金属に流れる電流と、誘電体の表面に流れる等価電
流，等価磁流とを未知数とするモーメント法の連立方程
式を導出する。
【００２６】導出手段１１１によりモーメント法の連立
方程式が導出されると、算出手段１１４は、この連立方
程式を解くことで、金属のみを解析対象とする場合に
は、金属に流れる電流を算出し、金属と誘電体とを解析
対象とする場合には、金属に流れる電流と、誘電体の表
面に流れる等価電流，等価磁流とを算出する。そして、
この算出結果を受けて、計算手段１１５は、この算出さ
れた算出値から、電気回路装置の放射する電磁界強度を
計算し、この計算結果を受けて、出力手段１２は、計算
された電磁界強度を例えば図式化して出力する。
【００２７】このように、本発明では、モーメント法に
基づいて電気回路装置の放射する電磁界強度を算出する
ときにあって、モーメント法の連立方程式を構成する相
互係数の算出にあたって、通常は高速の算出手法を用い
つつ、桁落ちが発生する可能性があるときには桁落ちし
ない高精度の算出手法を用いることで、相互係数の正確
な値を算出していく構成を採るものであることから、電
気回路装置の放射する電磁界強度を高速かつ高精度に算
出できるようになる。
【００２８】
【実施例】以下、金属のみを解析対象とするモーメント
法に適用した実施例に従って本発明を詳細に説明する。
【００２９】図２及び図３に、本発明の電磁界強度算出
装置１が実行する処理フローの一実施例を図示する。こ
こで、図中、２０は入力データファイルであって、解析
対象となるメッシュ化された電気回路装置の構造体情報
を管理するもの、２１は判断テーブルであって、相互イ
ンピーダンスの算出手法の選択に用いる判断データを管
理するもの、２２は出力データファイルであって、算出
された電磁界強度を格納するものである。
【００３０】図４に、判断テーブル２１の管理する判断
データの一実施例を図示する。この図に示すように、判
断テーブル２１は、相互インピーダンスの算出対象とな
る２つの金属要素の電気長の代表値と、２つの金属要素
の距離の電気長とを検索キーにして、２つの金属要素の
電気長の代表値が規定値よりも短くなるときにあって、
その電気長と比較して２つの金属要素の距離の電気長が
長い領域を領域Ｂ、それ以外の領域を領域Ａとして管理
する構成を採るものである。ここで、２つの金属要素の
電気長の代表値としては、２つの電気長の小さい方や、
２つの電気長の平均値や、２つの電気長の掛け算値の平
方根等が用いられる。また、電気長とは、上述したよう
に、電磁波の波長λを尺度する長さである。
【００３１】次に、図２及び図３に示す処理フローに従
って、本発明の電磁界強度算出装置１の実行する電磁界
強度算出処理について詳細に説明する。本発明の電磁界
強度算出装置１は、起動されると、図２及び図３の処理
フローに示すように、先ず最初に、ステップ１で、入力
データファイル２０からメッシュ化された電気回路装置
の構造体情報を読み込んで、金属要素や他のデータ（周
波数等）を構造体，配列として設定する。
【００３２】次に、ステップ２で、処理済の周波数をカ
ウントすることで、登録されている全周波数についての
処理が終了したのか否かを判断して、終了したことを判
断するときには全ての処理を終了し、終了していないこ
とを判断するときには、以下の処理を行う周波数を未処
理の中から１つ選択して、続くステップ３で、この選択
した周波数の波長λを計算する。
【００３３】続いて、ステップ４で、相互インピーダン
スＺijを算出するために、ｍ個ある金属要素の中から未
処理のものを１つ選択する。このステップ４で、未処理
の金属要素を選択できたことを判断すると、続いて、ス
テップ５で、相互インピーダンスＺijの算出対象となる
もう一方の金属要素をｍ個の中から１つ選択する。
【００３４】このステップ５で、全ての金属要素を選択
したことを判断するときには、ステップ４に戻ってい
く。一方、このステップ５で、未処理の金属要素を選択
できたことを判断するときには、ステップ６に進んで、
相互インピーダンスＺijの算出対象として選択した２つ
の金属要素の距離の電気長を特定するとともに、２つの
金属要素の電気長の代表値を算出して、この２つの電気
長，代表値を検索キーにして判断テーブル２１を検索す
ることで、相互インピーダンスＺijの算出対象として選
択した２つの金属要素が、領域Ａか領域Ｂのいずれに属
するのかを判断する。
【００３５】このステップ６で、領域Ａに属することを
判断するときには、ステップ７に進んで、演算速度を優
先して、高速演算処理を実現する倍精度実数の算出手法
に従って、選択した２つの金属要素の相互インピーダン
スＺijを算出してステップ５に戻っていく。一方、この
ステップ６で、領域Ｂに属することを判断するときに
は、ステップ８に進んで、演算精度を優先して、高精度
演算処理を実現する多倍精度実数や、固定小数点数（多
倍長整数）の算出手法に従って、選択した２つの金属要
素の相互インピーダンスＺijを算出してステップ５に戻
っていく。
【００３６】一方、ステップ４で、全ての金属要素を選
択したことを判断するとき、すなわち、全ての相互イン
ピーダンスＺijの算出完了を判断するときには、ステッ
プ９に進んで、算出した相互インピーダンスＺijと、入
力データファイル２０から読み込んだ波源Ｖi とを使っ
て、金属要素に流れる電流Ｉi を未知数とするモーメン
ト法の連立方程式 〔Ｚij〕〔Ｉi 〕＝〔Ｖi 〕 を導出し、これを解くことで、金属要素に流れる電流Ｉ
i を求める。
【００３７】そして、続くステップ１０で、処理済の観
測点をカウントすることで、登録されている全観測点に
ついての処理が終了したのか否かを判断して、終了した
ことを判断するときにはステップ２に戻り、終了してい
ないことを判断するときには、ステップ１１に進んで、
算出した電流Ｉi が観測点にもたらす電磁界強度を規定
の算出式に従って算出して、その算出結果を出力データ
ファイル２２に格納してからステップ１０に戻ってい
く。
【００３８】このようにして、本発明では、モーメント
法に基づいて電気回路装置の放射する電磁界強度を算出
するときにあって、相互インピーダンスＺijの算出にあ
たって、倍精度実数の算出手法では桁落ちが発生する可
能性があるときには、大きな作業域が必要とされ、プロ
グラムも複雑になり、計算時間も大幅に増大するもの
の、桁落ちせずに高精度演算を実現する多倍精度実数や
多倍長整数の算出手法を用いていく構成を採ることか
ら、電気回路装置の放射する電磁界強度を正確に算出で
きるようになる。
【００３９】次に、図２及び図３の処理フローで用いた
判断テーブル２１の管理する判断データについて具体的
に説明する。本発明の電磁界強度算出装置１は、金属の
ワイヤ部分について線分割することで、図５に示すよう
にワイヤのダイポールを生成する。
【００４０】電磁界強度算出装置１は、空気中であるの
か誘電体内であるのかを表す識別子をｑ（以下、空気中
であるときにはｑ＝０、誘電体内であるときにはｑ＝ｄ
とする）、２πを波長λで割った値をｋで表すならば、
このようにして生成した金属のワイヤのダイポールにつ
いては、
【００４１】
【数１】
【００４２】という展開関数Ｊ（式中に示す矢印の付加
されたＺは、ｚ方向の単位ベクトルであり、ａはワイヤ
の半径である）を想定して、金属のワイヤ部分に流れる
電流をＪs で表すならば、この展開関数を使って、電流
Ｊs を
【００４３】
【数２】
【００４４】のようにＮc 個の展開関数のセットに展開
する。ここで、この係数Ｉn が電流の大きさを表し、最
終的にはモーメント法の連立方程式の未知数となる。こ
のようなことを前提にして、図６に示すような２つの金
属ワイヤダイポール間の相互インピーダンスの算出処理
について考えることにする。ここで、この金属ワイヤダ
イポール間の相互インピーダンスＺmnは、図中にも示す
ように、モノポール間の相互インピーダンスＺ00，
Ｚ01，Ｚ10，Ｚ11の和で表されるものである。
【００４５】金属メッシュｕのｍ番目の展開関数の電流
が金属メッシュｖに与える電界をＥ q （Ｊu,ｍ）で表す
ならば、金属メッシュｕのｍ番目の展開関数の電流と、
金属メッシュｖのｎ番目の展開関数の電流との間の相互
インピーダンスは、
【００４６】
【数３】
【００４７】で表される。ここで、図７に示すような波
源と観測点との関係を想定するとともに、波源の電流分
布として図中に示すＪ+ を想定すると、電流Ｊ+ による
電界のｚ方向成分Ｅz + と、ρ方向成分Ｅρ+ は、図８
に示す算出式に従って算出されることになる。
【００４８】これから、２つのモノポール間の座標系と
して、図９に示す座標系を想定するとともに、図中に示
すワイヤ１の電流Ｊ1 + 、ワイヤ２の電流Ｊ2 + を想定
すると、ワイヤ１の電流Ｊ1 + によりワイヤ２の点に生
ずる電界Ｅz + ，Ｅρ+ の内、ワイヤ２の接線方向の成
分はＥz + のみであり、また、ワイヤ２の電流Ｊ2 +に
よりワイヤ２の点に生ずる電界Ｅz + ，Ｅρ+ の内、ワ
イヤ１の接線方向の成分はＥz + のみであるという点を
考慮すると、図６に示すモノポール間の相互インピーダ
ンスＺ00は、
【００４９】
【数４】
【００５０】のように算出される。ここで、式中のＸ
は、
【００５１】
【数５】
【００５２】で定義され、Ｙは、
【００５３】
【数６】
【００５４】で定義され、また、
【００５５】
【数７】
【００５６】という定義がある。同様にして、図６に示
すモノポール間の相互インピーダンスＺ01は、
【００５７】
【数８】
【００５８】のように算出され、図６に示すモノポール
間の相互インピーダンスＺ10は、
【００５９】
【数９】
【００６０】のように算出され、図６に示すモノポール
間の相互インピーダンスＺ11は、
【００６１】
【数１０】
【００６２】のように算出される。このようにして、モ
ノポール間の相互インピーダンスＺ00，Ｚ01，Ｚ10，Ｚ
11が算出されると、これらの和を算出することで、図６
に示すような２つの金属ワイヤダイポール間の相互イン
ピーダンスＺmnが算出されることになる。
【００６３】図１０に、この相互インピーダンスＺmnの
数値計算例を図示する。ここで、この数値計算は、２つ
の金属ワイヤダイポールの長さを同一のものとして、ダ
イポールの長さと、ダイポール間の距離とを電気長で表
現して行った。
【００６４】図中の左側部分が、高速演算処理を可能に
する倍精度実数の算出手法を用いた場合の算出値であ
り、右側部分が、高精度演算処理を可能にする固定小数
点数（多倍長整数）の算出手法を用いた場合の算出値で
ある。
【００６５】この２つの算出値を比較すれば分かるよう
に、図中の破線で囲った倍精度実数の算出値の精度が落
ちることが分かる。これは、この領域が桁落ちの発生す
る領域だからである。これから、判断テーブル２１は、
この破線で囲った領域を領域Ｂとして管理し、それ以外
の領域を領域Ａとして管理することで、相互インピーダ
ンスが高速かつ高精度に算出される構成を実現するので
ある。
【００６６】以上に説明した実施例では、金属のみを解
析対象とするモーメント法に適用した実施例に従って本
発明を開示したが、本発明は、金属と誘電体とを解析対
象とするモーメント法にもそのまま適用できるものであ
る。
【００６７】すなわち、金属と誘電体とを解析対象とす
るモーメント法では、分割した金属同士の相互インピー
ダンスＺ0 CCと、分割した金属と誘電体との相互インピ
ーダンスＺ0 Cd，Ｚ0 dcと、分割した誘電体同士の相互
インピーダンスＺ0 dd，Ｚd ddと、分割した誘電体同士
の相互アドミッタンスＹ0 dd，Ｙd ddと、分割した金属
と誘電体との相互リアクションＢ0 Cd，Ｂ0 dcと、分割
した誘電体同士の相互リアクションＢ0 dd，Ｂd ddとを
算出する。
【００６８】ここで、肩付き文字“０”は空気中での計
算値を表し、肩付き文字“ｄ”は誘電体中での計算値を
表し、添え字“ｃ”は金属を表し、添え字“ｄ”は誘電
体を表し、添え字“ｃｃ”は、金属から金属への相互関
係を表し、添え字“ｄｄ”は、誘電体から誘電体への相
互関係を表し、添え字“ｃｄ”は、誘電体から金属への
相互関係を表し、添え字“ｄｃ”は、金属から誘電体へ
の相互関係を表している。
【００６９】そして、相互インピーダンス，相互アドミ
ッタンス，相互リアクションを算出すると、波源をＶで
表すならば、この波源Ｖと、分割した金属に流れる電流
Ｉcと、分割した誘電体の表面に流れる等価電流Ｉ
d と、分割した誘電体の表面に流れる等価磁流Ｍとを使
い、金属の表面電界値がゼロを示すという境界条件から
導かれるモーメント法の連立方程式 〔Ｚ0 CC〕〔Ｉc 〕＋〔Ｚ0 Cd〕〔Ｉd 〕＋〔Ｂ0 Cd〕
〔Ｍ〕＝〔Ｖ〕 と、誘電体の境界面の両側で電界の接線成分が等しいと
いう境界条件から導かれるモーメント法の連立方程式 〔Ｚ0 dc〕〔Ｉc 〕＋〔Ｚ0 dd＋Ｚd dd〕〔Ｉd 〕＋
〔Ｂ0 dd＋Ｂd dd〕〔Ｍ〕＝
〔０〕 と、誘電体の境界面の両側で磁界の接線成分が等しいと
いう境界条件から導かれるモーメント法の連立方程式 〔Ｂ0 dc〕〔Ｉc 〕＋〔Ｂ0 dd＋Ｂd dd〕〔Ｉd 〕＋
〔−Ｙ0 dd−Ｙd dd〕〔Ｍ〕＝
〔０〕 とを導出する。
【００７０】そして、このモーメント法の連立方程式を
解くことで、金属に流れる電流Ｉcと、分割した誘電体
の表面に流れる等価電流Ｉd と、分割した誘電体の表面
に流れる等価磁流Ｍとを求めて、これらの値から、電気
回路装置の放射する電磁界強度を算出していくことにな
る。
【００７１】このように、金属と誘電体とを解析対象と
するモーメント法では、相互インピーダンス，相互アド
ミッタンス，相互リアクションを算出していく構成を採
るので、この相互係数の算出にあたって本発明を用いれ
ば、電気回路装置の放射する電磁界強度を高速かつ高精
度に算出できるようになる。
【００７２】
【発明の効果】以上説明したように、本発明によれば、
モーメント法に基づいて電気回路装置の放射する電磁界
強度を算出するときにあって、モーメント法の連立方程
式を構成する相互係数の算出にあたって、通常は高速の
算出手法を用いつつ、桁落ちが発生する可能性があると
きには桁落ちしない高精度の算出手法を用いることで、
相互係数の正確な値を算出していく構成を採るものであ
ることから、電気回路装置の放射する電磁界強度を高速
かつ高精度に算出できるようになる。
【図面の簡単な説明】
【図１】本発明の原理構成図である。
【図２】本発明の実行する処理フローの一実施例であ
る。
【図３】本発明の実行する処理フローの一実施例であ
る。
【図４】判断テーブルの管理する判断データの一実施例
である。
【図５】ワイヤダイポールの説明図である。
【図６】２つのダイポール間の相互インピーダンスの説
明図である。
【図７】波源と観測点との関係図である。
【図８】電界の算出式の説明図である。
【図９】２つのモノポール間の相互インピーダンスの説
明図である。
【図１０】相互インピーダンスの数値計算例である。
【図１１】微小ループアンテナ近似法の説明図である。
【図１２】分布定数線路近似法の説明図である。
【符号の説明】
１ 電磁界強度算出装置 １０ 入力手段 １１ 電磁界強度算出手段 １２ 出力手段 １１０ 分割手段 １１１ 導出手段 １１２ 評価手段 １１３ 決定手段 １１４ 算出手段 １１５ 計算手段
───────────────────────────────────────────────────── フロントページの続き (72)発明者 西野 関司 神奈川県川崎市中原区上小田中1015番地 富士通株式会社内

Claims (6)

    【特許請求の範囲】
  1. 【請求項１】 金属対象物を扱うモーメント法に基づい
    て電気回路装置の放射する電磁界強度を算出する電磁界
    強度算出装置において、 モーメント法の連立方程式を構成する相互インピーダン
    スの算出にあたって、該相互インピーダンスの算出対象
    となる要素の電気長と、要素間距離の電気長とを評価す
    る評価手段(112) と、 上記評価手段(112) の評価結果に応じて、相互インピー
    ダンスの算出手法を決定する決定手段(113) とを備え、 上記決定手段(113) の決定する算出手法に従って相互イ
    ンピーダンスを算出することで、モーメント法の連立方
    程式を導出していくよう構成されることを、 特徴とする電磁界強度算出装置。
  2. 【請求項２】 金属対象物及び誘電体対象物を扱うモー
    メント法に基づいて電気回路装置の放射する電磁界強度
    を算出する電磁界強度算出装置において、 モーメント法の連立方程式を構成する相互インピーダン
    ス、相互アドミッタンス及び相互リアクションの内の一
    部又は全ての相互係数の算出にあたって、該相互係数の
    算出対象となる要素の電気長と、要素間距離の電気長と
    を評価する評価手段(112) と、 上記評価手段(112) の評価結果に応じて、評価対象とな
    る相互係数の算出手法を決定する決定手段(113) とを備
    え、 上記決定手段(113) の決定する算出手法に従って相互係
    数を算出することで、モーメント法の連立方程式を導出
    していくよう構成されることを、 特徴とする電磁界強度算出装置。
  3. 【請求項３】 請求項１又は２記載の電磁界強度算出装
    置において、 評価手段(112) は、評価対象となる２つの要素の電気長
    に規定の演算を施すことで、評価処理に用いる要素の電
    気長を特定していくよう処理することを、 特徴とする電磁界強度算出装置。
  4. 【請求項４】 請求項１、２又は３記載の電磁界強度算
    出装置において、 決定手段(113) は、要素の電気長が規定値よりも短くな
    るときにあって、該電気長と比較して要素間距離の電気
    長が長いことが評価されるときには、高精度の算出手法
    を用いることを決定し、そうでないことが評価されると
    きには、高速の算出手法を用いることを決定していくよ
    う処理することを、 特徴とする電磁界強度算出装置。
  5. 【請求項５】 請求項４記載の電磁界強度算出装置にお
    いて、 決定手段(113) は、高精度の算出手法として、多倍精度
    実数の算出手法を用いることを決定していくよう処理す
    ることを、 特徴とする電磁界強度算出装置。
  6. 【請求項６】 請求項４記載の電磁界強度算出装置にお
    いて、 決定手段(113) は、高精度の算出手法として、多倍長整
    数の算出手法を用いることを決定していくよう処理する
    ことを、 特徴とする電磁界強度算出装置。
JP9536294A 1994-05-10 1994-05-10 電磁界強度算出装置 Pending JPH07302277A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP9536294A JPH07302277A (ja) 1994-05-10 1994-05-10 電磁界強度算出装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP9536294A JPH07302277A (ja) 1994-05-10 1994-05-10 電磁界強度算出装置

Publications (1)

Publication Number Publication Date
JPH07302277A true JPH07302277A (ja) 1995-11-14

Family

ID=14135527

Family Applications (1)

Application Number Title Priority Date Filing Date
JP9536294A Pending JPH07302277A (ja) 1994-05-10 1994-05-10 電磁界強度算出装置

Country Status (1)

Country Link
JP (1) JPH07302277A (ja)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09237291A (ja) * 1995-12-28 1997-09-09 Fujitsu Ltd 電磁界強度算出装置及び電磁界強度算出方法
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
US6456949B1 (en) 1999-10-20 2002-09-24 Fujitsu Limited Method and apparatus for calculating electromagnetic field intensity, and a computer-readable recording medium
US6604225B1 (en) 1995-11-16 2003-08-05 Fujitsu Limited Calculation of electromagnetic field intensity by moment method
US7266794B2 (en) 2004-04-28 2007-09-04 Fujitsu Limited Apparatus for and method of analyzing transmission characteristics of a circuit apparatus
CN101794355A (zh) * 2010-03-26 2010-08-04 中国人民解放军空军工程大学 电大物体电磁辐射和散射的计算机核外并行计算方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6604225B1 (en) 1995-11-16 2003-08-05 Fujitsu Limited Calculation of electromagnetic field intensity by moment method
JPH09237291A (ja) * 1995-12-28 1997-09-09 Fujitsu Ltd 電磁界強度算出装置及び電磁界強度算出方法
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
US6456949B1 (en) 1999-10-20 2002-09-24 Fujitsu Limited Method and apparatus for calculating electromagnetic field intensity, and a computer-readable recording medium
US7266794B2 (en) 2004-04-28 2007-09-04 Fujitsu Limited Apparatus for and method of analyzing transmission characteristics of a circuit apparatus
CN101794355A (zh) * 2010-03-26 2010-08-04 中国人民解放军空军工程大学 电大物体电磁辐射和散射的计算机核外并行计算方法

Similar Documents

Publication Publication Date Title
Rao Time domain electromagnetics
EP0778533B1 (en) Electromagnetic field intensity calculation apparatus and method
JP2768900B2 (ja) 電磁界強度算出装置
Pinchuk et al. Error estimation for automatic adaptive finite element mesh generation
Obelleiro et al. Hybrid moment‐method physical‐optics formulation for modeling the electromagnetic behavior of on‐board antennas
Sun et al. Analysis of singly and doubly periodic absorbers by frequency-domain finite-difference method
Tanaka et al. Improved primary characteristic basis function method for monostatic radar cross section analysis of specific coordinate plane
JPH07302277A (ja) 電磁界強度算出装置
US6456949B1 (en) Method and apparatus for calculating electromagnetic field intensity, and a computer-readable recording medium
Awadalla et al. Monopole antenna at center of circular ground plane: Input impedance and radiation pattern
CN116070562A (zh) 一种复合平面波激励下的传输线高频耦合方法及系统
JP3443162B2 (ja) 電磁界強度算出装置
JP3714499B2 (ja) 電磁界強度算出装置及び電磁界強度算出方法
Chumachenko Domain-product technique solution for the problem of electromagnetic scattering from multiangular composite cylinders
JPH08304491A (ja) 電磁界強度算出装置
Zhu et al. Numerical computation of incomplete Lipschitz-Hankel integrals of the Hankel type for complex-valued arguments
US6324493B1 (en) Method and apparatus for modeling electromagnetic interactions in electrical circuit metalizations to simulate their electrical characteristics
JPH08304494A (ja) 電磁界強度算出装置
Öing et al. Calculation of radiated electromagnetic fields from electronic systems
JPH08304493A (ja) 電磁界強度算出装置
Helmers et al. Analyzing electromagnetic pulse coupling by combining TLT, MoM, and GTD/UTD
JPH09237291A (ja) 電磁界強度算出装置及び電磁界強度算出方法
Chen Application of Eigenmode Theory to the Analysis of Low-frequency Problems
Altun et al. A hybrid numerical model for long-range electromagnetic wave propagation
Wang et al. A study of the reactions between rooftop test and expansion functions in inhomogeneously-filled stripline structures

Legal Events

Date Code Title Description
A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20030610