JP3405905B2 - 電磁界強度算出装置及び方法並びにプログラム記憶媒体 - Google Patents
電磁界強度算出装置及び方法並びにプログラム記憶媒体Info
- Publication number
- JP3405905B2 JP3405905B2 JP21758397A JP21758397A JP3405905B2 JP 3405905 B2 JP3405905 B2 JP 3405905B2 JP 21758397 A JP21758397 A JP 21758397A JP 21758397 A JP21758397 A JP 21758397A JP 3405905 B2 JP3405905 B2 JP 3405905B2
- Authority
- JP
- Japan
- Prior art keywords
- cable
- electromagnetic field
- cable connection
- electric circuit
- circuit device
- 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
- 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
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/001—Measuring interference from external sources to, or emission from, the device under test, e.g. EMC, EMI, EMP or ESD testing
- G01R31/002—Measuring interference from external sources to, or emission from, the device under test, e.g. EMC, EMI, EMP or ESD testing where the device under test is an electronic circuit
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/16—Cables, cable trees or wire harnesses
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
【０００１】
【発明の属する技術分野】本発明は、モーメント法を使
って、電気回路装置の放射する電磁界強度を算出する電
磁界強度算出装置及び方法と、その電磁界強度算出装置
を実現するプログラムが記憶されるプログラム記憶媒体
とに関し、特に、電気回路装置の持つ２線ケーブルの放
射する電磁界強度をユーザフレンドリーに、しかも高精
度に算出できるようにする電磁界強度算出装置及び方法
と、その電磁界強度算出装置を実現するプログラムが記
憶されるプログラム記憶媒体とに関する。
って、電気回路装置の放射する電磁界強度を算出する電
磁界強度算出装置及び方法と、その電磁界強度算出装置
を実現するプログラムが記憶されるプログラム記憶媒体
とに関し、特に、電気回路装置の持つ２線ケーブルの放
射する電磁界強度をユーザフレンドリーに、しかも高精
度に算出できるようにする電磁界強度算出装置及び方法
と、その電磁界強度算出装置を実現するプログラムが記
憶されるプログラム記憶媒体とに関する。
【０００２】電気回路装置に対する社会的規制として、
一定のレベル以上の不要な電波やノイズを放射してはな
らないということがあり、各国の規格で厳しく規定され
るようになってきた。
一定のレベル以上の不要な電波やノイズを放射してはな
らないということがあり、各国の規格で厳しく規定され
るようになってきた。
【０００３】このような電波規格を満足させるために、
シールド技術やフィルタ技術等のような種々の対策技術
が用いられるが、これらの対策技術の採用に当たって、
それらがどの程度電波を減少できるのかを定量的に算出
できるようにするシミュレーション技術の開発が必要で
ある。
シールド技術やフィルタ技術等のような種々の対策技術
が用いられるが、これらの対策技術の採用に当たって、
それらがどの程度電波を減少できるのかを定量的に算出
できるようにするシミュレーション技術の開発が必要で
ある。
【０００４】このようなことを背景にして、本発明者ら
は、モーメント法を使って、電気回路装置の放射する電
磁界強度を算出するシミュレーション技術の発明を開示
してきた。このシミュレーション技術を確立していくに
は、電気回路装置の正確なモデルを構築していく必要が
ある。
は、モーメント法を使って、電気回路装置の放射する電
磁界強度を算出するシミュレーション技術の発明を開示
してきた。このシミュレーション技術を確立していくに
は、電気回路装置の正確なモデルを構築していく必要が
ある。
【０００５】
【従来の技術】物体の放射する電磁界強度は、物体各部
に流れる電流や磁流を求めて、それを公知の電磁波放射
の理論式に代入することでシミュレーションできる。こ
の物体各部に流れる電流や磁流は、理論的には、マック
スウェルの電磁波方程式を与えられた境界条件の下に解
くことで得られる。
に流れる電流や磁流を求めて、それを公知の電磁波放射
の理論式に代入することでシミュレーションできる。こ
の物体各部に流れる電流や磁流は、理論的には、マック
スウェルの電磁波方程式を与えられた境界条件の下に解
くことで得られる。
【０００６】これを解くものとしてモーメント法があ
る。モーメント法は、マックスウェルの電磁波方程式か
ら導かれる積分方程式の解法の１つで、物体を小さな要
素に分割して電流や磁流の計算を行う手法であり、３次
元の任意形状物体を扱うことができる。このモーメント
法についての参考文献としては、「H.N.Wang, J.H.Rich
mond and M.C.Gilreath:"Sinusoidal reaction formula
tion for radiation andscattering from cond-ucting
surface" IEEE TRANSACTIONS ANTENNAS PROPAGATION vo
l.AP-23 1975 」がある。
る。モーメント法は、マックスウェルの電磁波方程式か
ら導かれる積分方程式の解法の１つで、物体を小さな要
素に分割して電流や磁流の計算を行う手法であり、３次
元の任意形状物体を扱うことができる。このモーメント
法についての参考文献としては、「H.N.Wang, J.H.Rich
mond and M.C.Gilreath:"Sinusoidal reaction formula
tion for radiation andscattering from cond-ucting
surface" IEEE TRANSACTIONS ANTENNAS PROPAGATION vo
l.AP-23 1975 」がある。
【０００７】このモーメント法では、シミュレーション
対象となる電気回路装置の構造をメッシュ化し、処理対
象の周波数を選択すると、その周波数について、メッシ
ュ化された要素間の相互インピーダンスや相互アドミッ
タンスや相互リアクションを所定の計算処理によって求
めて、その求めた相互インピーダンス等と構造情報で指
定される波源とをモーメント法の連立方程式に代入し、
それを解くことで各要素に流れる電流や磁流を求めるこ
とになる。
対象となる電気回路装置の構造をメッシュ化し、処理対
象の周波数を選択すると、その周波数について、メッシ
ュ化された要素間の相互インピーダンスや相互アドミッ
タンスや相互リアクションを所定の計算処理によって求
めて、その求めた相互インピーダンス等と構造情報で指
定される波源とをモーメント法の連立方程式に代入し、
それを解くことで各要素に流れる電流や磁流を求めるこ
とになる。
【０００８】すなわち、金属対象物を扱うときには、金
属部分を解析対象としてメッシュ化し、メッシュ化した
金属要素間の相互インピーダンスＺijを求め、この相互
インピーダンスＺijと、波源Ｖi と、メッシュ化した金
属要素に流れる電流Ｉi との間に成立するモーメント法
の連立方程式 〔Ｖi 〕＝〔Ｚij〕〔Ｉi 〕 但し、〔 〕はマトリックス を解くことで電流Ｉi を求めて、この結果から電磁界強
度を算出する方法を採っている。
属部分を解析対象としてメッシュ化し、メッシュ化した
金属要素間の相互インピーダンスＺijを求め、この相互
インピーダンスＺijと、波源Ｖi と、メッシュ化した金
属要素に流れる電流Ｉi との間に成立するモーメント法
の連立方程式 〔Ｖi 〕＝〔Ｚij〕〔Ｉi 〕 但し、〔 〕はマトリックス を解くことで電流Ｉi を求めて、この結果から電磁界強
度を算出する方法を採っている。
【０００９】ここで、メッシュに、抵抗やコンデンサや
リアクタンスが存在するときには、それらは、そのメッ
シュの自己インピーダンス成分に加わることになる。と
ころで、電気回路装置では、ケーブルを使って、離れた
２点間を接続する構成を採ることが多い。
リアクタンスが存在するときには、それらは、そのメッ
シュの自己インピーダンス成分に加わることになる。と
ころで、電気回路装置では、ケーブルを使って、離れた
２点間を接続する構成を採ることが多い。
【００１０】しかるに、従来の電磁界強度算出装置で
は、電気回路装置がケーブルを持つ場合にあっても、特
別なユーザインタフェースを用意せずに、他の構造物と
同様の扱いに従って、ケーブルを固定的な構造物として
扱って、モーメント法のモデルを生成していくという方
法を採っていた。
は、電気回路装置がケーブルを持つ場合にあっても、特
別なユーザインタフェースを用意せずに、他の構造物と
同様の扱いに従って、ケーブルを固定的な構造物として
扱って、モーメント法のモデルを生成していくという方
法を採っていた。
【００１１】
【発明が解決しようとする課題】しかしながら、ケーブ
ルの配置は自由度が大きく、２点間を接続させるため
に、電気回路装置内で、折れ曲げられたり、捩じられた
りする。しかも、ケーブルの放射する電磁界強度は大き
なものである。
ルの配置は自由度が大きく、２点間を接続させるため
に、電気回路装置内で、折れ曲げられたり、捩じられた
りする。しかも、ケーブルの放射する電磁界強度は大き
なものである。
【００１２】これから、従来の電磁界強度算出装置のよ
うに、ケーブルを固定的な構造物として扱っていくとい
う構成を採っていると、実際のものに合ったシミュレー
ションを実行できないという問題点があるとともに、ケ
ーブルの配置を色々と変えてみて、その効果を見るとい
うようなことが簡単にできないという問題点があった。
うに、ケーブルを固定的な構造物として扱っていくとい
う構成を採っていると、実際のものに合ったシミュレー
ションを実行できないという問題点があるとともに、ケ
ーブルの配置を色々と変えてみて、その効果を見るとい
うようなことが簡単にできないという問題点があった。
【００１３】更に、このケーブルには、平行線ケーブル
や撚り線ケーブルといったように様々な種類があり、同
一種類でも様々な規格のものがある。しかるに、従来の
電磁界強度算出装置では、ケーブルに対して統一的な観
点からモデル化を行っておらず、これから、ケーブルの
種類が変わると、どのようにモデル化するのかをその都
度考えていかなくてはならないという問題点があった。
特に、撚り線ケーブルのモデル化は、実質的に構築され
ておらず、高精度のモデルの構築が必要とされていた。
や撚り線ケーブルといったように様々な種類があり、同
一種類でも様々な規格のものがある。しかるに、従来の
電磁界強度算出装置では、ケーブルに対して統一的な観
点からモデル化を行っておらず、これから、ケーブルの
種類が変わると、どのようにモデル化するのかをその都
度考えていかなくてはならないという問題点があった。
特に、撚り線ケーブルのモデル化は、実質的に構築され
ておらず、高精度のモデルの構築が必要とされていた。
【００１４】本発明はかかる事情に鑑みてなされたもの
であって、モーメント法を使って、電気回路装置の放射
する電磁界強度を算出する構成を採るときにあって、電
気回路装置の持つ２線ケーブルの放射する電磁界強度を
ユーザフレンドリーに、しかも高精度に算出できるよう
にする新たな電磁界強度算出装置及び方法の提供と、そ
の電磁界強度算出装置を実現するプログラムが記憶され
る新たなプログラム記憶媒体の提供とを目的とする。
であって、モーメント法を使って、電気回路装置の放射
する電磁界強度を算出する構成を採るときにあって、電
気回路装置の持つ２線ケーブルの放射する電磁界強度を
ユーザフレンドリーに、しかも高精度に算出できるよう
にする新たな電磁界強度算出装置及び方法の提供と、そ
の電磁界強度算出装置を実現するプログラムが記憶され
る新たなプログラム記憶媒体の提供とを目的とする。
【００１５】
【課題を解決するための手段】図１に本発明の原理構成
を図示する。図中、１は本発明を具備する電磁界強度算
出装置であって、モーメント法を使って、電気回路装置
の放射する電磁界強度を算出するものである。
を図示する。図中、１は本発明を具備する電磁界強度算
出装置であって、モーメント法を使って、電気回路装置
の放射する電磁界強度を算出するものである。
【００１６】本発明の電磁界強度算出装置１は、入力手
段１０と、ライブラリ手段１１と、設定手段１２と、モ
デル生成手段１３と、生成手段１４と、算出手段１５と
を備える。
段１０と、ライブラリ手段１１と、設定手段１２と、モ
デル生成手段１３と、生成手段１４と、算出手段１５と
を備える。
【００１７】入力手段１０は、解析対象となる電気回路
装置の構造情報を入力する。ライブラリ手段１１は、電
気回路装置に実装される２線ケーブルの属性情報を管理
する。設定手段１２は、入力手段１０の入力する構造情
報がケーブル接続箇所を持たないときには、対話処理に
従い、ケーブル接続箇所と、そのケーブル接続箇所を接
続する１つ又は複数のケーブル経由位置とを設定し、一
方、ケーブル接続箇所を持つときには、対話処理に従
い、そのケーブル接続箇所を接続する１つ又は複数のケ
ーブル経由位置を設定し、さらに、対話処理に従い、使
用する２線ケーブルの識別情報を設定する。
装置の構造情報を入力する。ライブラリ手段１１は、電
気回路装置に実装される２線ケーブルの属性情報を管理
する。設定手段１２は、入力手段１０の入力する構造情
報がケーブル接続箇所を持たないときには、対話処理に
従い、ケーブル接続箇所と、そのケーブル接続箇所を接
続する１つ又は複数のケーブル経由位置とを設定し、一
方、ケーブル接続箇所を持つときには、対話処理に従
い、そのケーブル接続箇所を接続する１つ又は複数のケ
ーブル経由位置を設定し、さらに、対話処理に従い、使
用する２線ケーブルの識別情報を設定する。
【００１８】モデル生成手段１３は、入力手段１０の入
力する電気回路装置の構造情報から、電気回路装置のモ
ーメント法モデルを生成する。生成手段１４は、モデル
生成手段１３に展開されて、ライブラリ手段１１の管理
する２線ケーブルの属性情報を使い、ケーブル接続箇所
を接続する２線ケーブルのモーメント法モデルを生成す
る。算出手段１４は、モデル生成手段１３の生成するモ
ーメント法モデルに基づき、モーメント法を用いて電気
回路装置の放射する電磁界強度を算出する。
力する電気回路装置の構造情報から、電気回路装置のモ
ーメント法モデルを生成する。生成手段１４は、モデル
生成手段１３に展開されて、ライブラリ手段１１の管理
する２線ケーブルの属性情報を使い、ケーブル接続箇所
を接続する２線ケーブルのモーメント法モデルを生成す
る。算出手段１４は、モデル生成手段１３の生成するモ
ーメント法モデルに基づき、モーメント法を用いて電気
回路装置の放射する電磁界強度を算出する。
【００１９】ここで、本発明の電磁界強度算出装置１の
持つ機能は具体的にはプログラムで実現されるものであ
り、このプログラムは、フロッピィディスクなどに記憶
されたり、サーバなどのディスクなどに記憶され、それ
らから電磁界強度算出装置１にインストールされてメモ
リ上で動作することで、本発明を実現することになる。
持つ機能は具体的にはプログラムで実現されるものであ
り、このプログラムは、フロッピィディスクなどに記憶
されたり、サーバなどのディスクなどに記憶され、それ
らから電磁界強度算出装置１にインストールされてメモ
リ上で動作することで、本発明を実現することになる。
【００２０】このように構成される本発明の電磁界強度
算出装置１では、入力手段１０が電気回路装置の構造情
報を入力すると、設定手段１２は、対話処理に従って、
ケーブル経由位置を設定しつつ、電気回路装置に実装さ
れる２線ケーブルの配置位置を設定するとともに、使用
する２線ケーブルの識別情報を設定する。
算出装置１では、入力手段１０が電気回路装置の構造情
報を入力すると、設定手段１２は、対話処理に従って、
ケーブル経由位置を設定しつつ、電気回路装置に実装さ
れる２線ケーブルの配置位置を設定するとともに、使用
する２線ケーブルの識別情報を設定する。
【００２１】この設定手段１２により設定された２線ケ
ーブルとその配置位置とを受けて、生成手段１４は、２
線ケーブルとして撚り線ケーブルが用いられるときに
は、２本の芯線の長さが同一となることを保証すべく、
ケーブル方向を変化させるケーブル経由位置に、撚り線
の腹位置が配置されることになるようにと撚り線ケーブ
ルの撚りピッチを変更するとともに、その変更した撚り
ピッチと、規定の回転角度との指定するケーブル芯線位
置に従って、ケーブルの方向が変化しないケーブル開始
点とケーブル終了点との間に線部材を生成するととも
に、その線部材の間にコンデンサを生成することで、ケ
ーブル接続箇所を接続する２線ケーブルのモーメント法
モデルを生成する。
ーブルとその配置位置とを受けて、生成手段１４は、２
線ケーブルとして撚り線ケーブルが用いられるときに
は、２本の芯線の長さが同一となることを保証すべく、
ケーブル方向を変化させるケーブル経由位置に、撚り線
の腹位置が配置されることになるようにと撚り線ケーブ
ルの撚りピッチを変更するとともに、その変更した撚り
ピッチと、規定の回転角度との指定するケーブル芯線位
置に従って、ケーブルの方向が変化しないケーブル開始
点とケーブル終了点との間に線部材を生成するととも
に、その線部材の間にコンデンサを生成することで、ケ
ーブル接続箇所を接続する２線ケーブルのモーメント法
モデルを生成する。
【００２２】ここで、生成手段１４は、ケーブルの方向
が変化しないケーブル開始点とケーブル終了点との間に
捩れが存在するときには、その捩れを吸収することによ
うにと回転角度を補正して、その補正した回転角度を使
って線部材を生成することになる。
が変化しないケーブル開始点とケーブル終了点との間に
捩れが存在するときには、その捩れを吸収することによ
うにと回転角度を補正して、その補正した回転角度を使
って線部材を生成することになる。
【００２３】このとき、生成手段１４が撚り線ケーブル
の撚りピッチを変更するのではなくて、設定手段１２
が、撚り線の腹位置がケーブル経由位置に配置されるこ
とになるようにと、設定されるケーブル経由位置をその
近傍で移動することもある。
の撚りピッチを変更するのではなくて、設定手段１２
が、撚り線の腹位置がケーブル経由位置に配置されるこ
とになるようにと、設定されるケーブル経由位置をその
近傍で移動することもある。
【００２４】また、設定手段１２により設定された２線
ケーブルとその配置位置とを受けて、生成手段１４は、
２線ケーブルとして平行線ケーブルが用いられるときに
は、ケーブル方向を変化させないケーブル開始点とケー
ブル終了点との間を等間隔に区切る規定長以下の長さの
指定するケーブル芯線位置に従って、そのケーブル開始
点とケーブル終了点との間の線部材を生成するととも
に、その線部材の間にコンデンサを生成することで、ケ
ーブル接続箇所を接続する２線ケーブルのモーメント法
モデルを生成する。
ケーブルとその配置位置とを受けて、生成手段１４は、
２線ケーブルとして平行線ケーブルが用いられるときに
は、ケーブル方向を変化させないケーブル開始点とケー
ブル終了点との間を等間隔に区切る規定長以下の長さの
指定するケーブル芯線位置に従って、そのケーブル開始
点とケーブル終了点との間の線部材を生成するととも
に、その線部材の間にコンデンサを生成することで、ケ
ーブル接続箇所を接続する２線ケーブルのモーメント法
モデルを生成する。
【００２５】ここで、生成手段１４は、ケーブルの方向
が変化しないケーブル開始点とケーブル終了点との間に
捩れが存在するときには、その捩れを考慮したケーブル
芯線位置に従って線部材を生成することになる。
が変化しないケーブル開始点とケーブル終了点との間に
捩れが存在するときには、その捩れを考慮したケーブル
芯線位置に従って線部材を生成することになる。
【００２６】この処理にあたって、生成手段１４は、規
定の公式に従って、２線ケーブルの芯線が空気中に存在
するときの単位長さ当たりのインダクタンス及び容量を
求め、それを使って、２線ケーブルの空気中での特性イ
ンピーダンスを算出し、その特性インピーダンスと、ラ
イブラリ手段１１の管理する２線ケーブルの属性情報で
指定される特性インピーダンスとから、２線ケーブルの
単位長さ当たりの容量を求めて、その容量と線部材の数
とから、コンデンサの容量を決定することができる。
定の公式に従って、２線ケーブルの芯線が空気中に存在
するときの単位長さ当たりのインダクタンス及び容量を
求め、それを使って、２線ケーブルの空気中での特性イ
ンピーダンスを算出し、その特性インピーダンスと、ラ
イブラリ手段１１の管理する２線ケーブルの属性情報で
指定される特性インピーダンスとから、２線ケーブルの
単位長さ当たりの容量を求めて、その容量と線部材の数
とから、コンデンサの容量を決定することができる。
【００２７】このように、本発明の電磁界強度算出装置
１によれば、モーメント法を使って、電気回路装置の放
射する電磁界強度を算出する構成を採るときにあって、
電気回路装置に実装される２線ケーブルの配設位置を自
由に変更できるようにする構成を採るので、実際のもの
に合ったシミュレーションが実行できるようになるとと
もに、ケーブルの配置を色々と変えてみて、その効果を
見るというようなことが簡単に実行できるようになる。
１によれば、モーメント法を使って、電気回路装置の放
射する電磁界強度を算出する構成を採るときにあって、
電気回路装置に実装される２線ケーブルの配設位置を自
由に変更できるようにする構成を採るので、実際のもの
に合ったシミュレーションが実行できるようになるとと
もに、ケーブルの配置を色々と変えてみて、その効果を
見るというようなことが簡単に実行できるようになる。
【００２８】そして、平行線ケーブルや撚り線ケーブル
を、金属線分とコンデンサとからなるはしご型回路で高
精度にモデル化するので、平行線ケーブルや撚り線ケー
ブルを実装する電気回路装置から放射される電磁界強度
を高精度に算出できるようになる。
を、金属線分とコンデンサとからなるはしご型回路で高
精度にモデル化するので、平行線ケーブルや撚り線ケー
ブルを実装する電気回路装置から放射される電磁界強度
を高精度に算出できるようになる。
【００２９】
【発明の実施の形態】以下、実施の形態に従って本発明
を詳細に説明する。図２に、本発明の電磁界強度算出装
置１の実行する処理フローの一実施例を図示する。ここ
で、図中、１００は入力データファイルであって、解析
対象となる電気回路装置の構造情報を格納するもの、２
００はケーブルライブラリであって、解析対象となる電
気回路装置に実装される２線ケーブルの属性情報を管理
するもの、３００は出力データファイルであって、算出
結果の電磁界強度を格納するものである。
を詳細に説明する。図２に、本発明の電磁界強度算出装
置１の実行する処理フローの一実施例を図示する。ここ
で、図中、１００は入力データファイルであって、解析
対象となる電気回路装置の構造情報を格納するもの、２
００はケーブルライブラリであって、解析対象となる電
気回路装置に実装される２線ケーブルの属性情報を管理
するもの、３００は出力データファイルであって、算出
結果の電磁界強度を格納するものである。
【００３０】ここで、この処理フローを実行するプログ
ラムは、フロッピィディスクなどに記憶されたり、サー
バなどのディスクなどに記憶され、それらから電磁界強
度算出装置１にインストールされてメモリ上で動作する
ことになる。
ラムは、フロッピィディスクなどに記憶されたり、サー
バなどのディスクなどに記憶され、それらから電磁界強
度算出装置１にインストールされてメモリ上で動作する
ことになる。
【００３１】また、図２では省略しているが、本発明の
電磁界強度算出装置１は、ユーザとの対話手段となる端
末装置を備えている。本発明の電磁界強度算出装置１
は、起動されると、図２の処理フローに示すように、先
ず最初に、ステップ１で、入力データファイル１００か
ら、解析対象となる電気回路装置の構造情報を読み込
み、続くステップ２で、この読み込んだ電気回路装置の
構造情報に２線ケーブルを配置し、続いて、ステップ３
で、ステップ１で読み込んだ構造情報と、ステップ２で
配置した２線ケーブルとに基づいて、解析対象の電気回
路装置をメッシュに分割してモーメント法を適用するた
めのモデルを生成する。
電磁界強度算出装置１は、ユーザとの対話手段となる端
末装置を備えている。本発明の電磁界強度算出装置１
は、起動されると、図２の処理フローに示すように、先
ず最初に、ステップ１で、入力データファイル１００か
ら、解析対象となる電気回路装置の構造情報を読み込
み、続くステップ２で、この読み込んだ電気回路装置の
構造情報に２線ケーブルを配置し、続いて、ステップ３
で、ステップ１で読み込んだ構造情報と、ステップ２で
配置した２線ケーブルとに基づいて、解析対象の電気回
路装置をメッシュに分割してモーメント法を適用するた
めのモデルを生成する。
【００３２】続いて、ステップ４で、処理対象となる周
波数の中から未処理のものを１つ選択し、続くステップ
５で、全周波数の選択を完了したのか否かを判断して、
全周波数の選択を完了したことを判断するときには、処
理を終了する。
波数の中から未処理のものを１つ選択し、続くステップ
５で、全周波数の選択を完了したのか否かを判断して、
全周波数の選択を完了したことを判断するときには、処
理を終了する。
【００３３】一方、全周波数の選択を完了していないこ
とを判断するとき、すなわち、ステップ４で周波数を選
択できたことを判断するときには、ステップ６に進ん
で、グリーン関数を使って、分割されたメッシュ間の相
互インピーダンスＺij（ｉ＝１〜ｎ，ｊ＝１〜ｎ）を算
出する。なお、説明の便宜上、ここでは、解析対象が金
属物であることを想定して、相互インピーダンスのみを
算出することを想定している。
とを判断するとき、すなわち、ステップ４で周波数を選
択できたことを判断するときには、ステップ６に進ん
で、グリーン関数を使って、分割されたメッシュ間の相
互インピーダンスＺij（ｉ＝１〜ｎ，ｊ＝１〜ｎ）を算
出する。なお、説明の便宜上、ここでは、解析対象が金
属物であることを想定して、相互インピーダンスのみを
算出することを想定している。
【００３４】続いて、ステップ７で、ステップ６で算出
した相互インピーダンスＺijと、構造情報で指定される
波源の電圧値Ｖi とを使い、モーメント法の連立方程式 〔Ｖi 〕＝〔Ｚij〕〔Ｉi 〕 を解くことで、各メッシュに流れる電流Ｉi を算出す
る。
した相互インピーダンスＺijと、構造情報で指定される
波源の電圧値Ｖi とを使い、モーメント法の連立方程式 〔Ｖi 〕＝〔Ｚij〕〔Ｉi 〕 を解くことで、各メッシュに流れる電流Ｉi を算出す
る。
【００３５】続いて、ステップ８で、予め指定される観
測点の中から未処理のものを１つ選択し、続くステップ
９で、全観測点の選択を完了したのか否かを判断して、
全観測点の選択を完了したことを判断するときには、次
の周波数を処理すべくステップ４に戻っていく。
測点の中から未処理のものを１つ選択し、続くステップ
９で、全観測点の選択を完了したのか否かを判断して、
全観測点の選択を完了したことを判断するときには、次
の周波数を処理すべくステップ４に戻っていく。
【００３６】一方、全観測点の選択を完了していないこ
とを判断するとき、すなわち、ステップ８で観測点を選
択できたことを判断するときには、ステップ１０に進ん
で、ステップ７で算出した各メッシュに流れる電流Ｉi
が、ステップ８で選択した観測点にもたらす電磁界強度
を算出して、その算出結果を出力データファイル３００
に格納してから、次の観測点を処理すべくステップ７に
戻っていく。
とを判断するとき、すなわち、ステップ８で観測点を選
択できたことを判断するときには、ステップ１０に進ん
で、ステップ７で算出した各メッシュに流れる電流Ｉi
が、ステップ８で選択した観測点にもたらす電磁界強度
を算出して、その算出結果を出力データファイル３００
に格納してから、次の観測点を処理すべくステップ７に
戻っていく。
【００３７】このようにして、本発明の電磁界強度算出
装置１は、モーメント法を用いて、解析対象となる電気
回路装置から放射される電磁界強度を算出していくよう
に処理するのである。
装置１は、モーメント法を用いて、解析対象となる電気
回路装置から放射される電磁界強度を算出していくよう
に処理するのである。
【００３８】図３に、ケーブルライブラリ２００のデー
タ構造を図示する。ケーブルライブラリ２００は、解析
対象となる電気回路装置に実装されるケーブルの属性情
報を管理するものであって、この図に示すように、ケー
ブルの名称をキーにして、その名称の指すケーブルの持
つ特性インピーダンスや、そのケーブルの種別（平行線
ケーブル／撚り線ケーブル／同軸ケーブル）や、平行線
ケーブルの場合は導体直径ｄ／導体間距離δ、撚り線ケ
ーブルの場合には導体直径ｄ／導体間距離δ／撚りピッ
チψ（正は時計回り、負は反時計回り）といったような
各種の属性情報を管理するものである。
タ構造を図示する。ケーブルライブラリ２００は、解析
対象となる電気回路装置に実装されるケーブルの属性情
報を管理するものであって、この図に示すように、ケー
ブルの名称をキーにして、その名称の指すケーブルの持
つ特性インピーダンスや、そのケーブルの種別（平行線
ケーブル／撚り線ケーブル／同軸ケーブル）や、平行線
ケーブルの場合は導体直径ｄ／導体間距離δ、撚り線ケ
ーブルの場合には導体直径ｄ／導体間距離δ／撚りピッ
チψ（正は時計回り、負は反時計回り）といったような
各種の属性情報を管理するものである。
【００３９】図４ないし図６に、図２の処理フローのス
テップ２で実行する「ケーブル配置決定処理」の詳細な
処理フロー、図７ないし図９に、図２の処理フローのス
テップ３で実行する「ケーブルモデル化処理」の詳細な
処理フローを図示する。次に、これらの処理フローに従
って、本発明の電磁界強度算出装置１の実行する処理に
ついて詳細に説明する。
テップ２で実行する「ケーブル配置決定処理」の詳細な
処理フロー、図７ないし図９に、図２の処理フローのス
テップ３で実行する「ケーブルモデル化処理」の詳細な
処理フローを図示する。次に、これらの処理フローに従
って、本発明の電磁界強度算出装置１の実行する処理に
ついて詳細に説明する。
【００４０】ここで、「ケーブル配置決定処理」は、ユ
ーザとの対話処理に従って実行する「ケーブル配置処
理」と、「ケーブル配置処理」の処理結果を受けて実行
する「リンク情報決定処理」とから構成されている。
ーザとの対話処理に従って実行する「ケーブル配置処
理」と、「ケーブル配置処理」の処理結果を受けて実行
する「リンク情報決定処理」とから構成されている。
【００４１】本発明の電磁界強度算出装置１は、図２の
処理フローのステップ２に従って「ケーブル配置決定処
理」に入ると、先ず最初に、図４の処理フローに従っ
て、「ケーブル配置処理」を実行する。
処理フローのステップ２に従って「ケーブル配置決定処
理」に入ると、先ず最初に、図４の処理フローに従っ
て、「ケーブル配置処理」を実行する。
【００４２】すなわち、ユーザからの要求に応答して
「ケーブル配置処理」に入ると、本発明の電磁界強度算
出装置１は、図４の処理フローに示すように、先ず最初
に、ステップ１で、ユーザと対話することで、２線ケー
ブルの入口の座標を設定し、続くステップ２で、ユーザ
と対話することで、２線ケーブルの出口の座標を設定す
る。ここで、この２線ケーブルの入口及び出口の座標
は、２線ケーブルの芯線（導体）の中心位置で定義され
るものであることから、２つの点の座標で設定されるこ
とになる。
「ケーブル配置処理」に入ると、本発明の電磁界強度算
出装置１は、図４の処理フローに示すように、先ず最初
に、ステップ１で、ユーザと対話することで、２線ケー
ブルの入口の座標を設定し、続くステップ２で、ユーザ
と対話することで、２線ケーブルの出口の座標を設定す
る。ここで、この２線ケーブルの入口及び出口の座標
は、２線ケーブルの芯線（導体）の中心位置で定義され
るものであることから、２つの点の座標で設定されるこ
とになる。
【００４３】続いて、ステップ３で、ユーザと対話する
ことで、ステップ１で設定した２線ケーブルの入口か
ら、ステップ２で設定した２線ケーブルの出口に至るま
でのケーブル経由点の座標を設定する。ここで、このケ
ーブル経由点の座標は、２線ケーブルの中心位置（２つ
の芯線の中間位置）で定義されている。なお、ケーブル
経由点が指定されないときには、ケーブルの入口から出
口まで一直線にケーブルが配置されることになる。
ことで、ステップ１で設定した２線ケーブルの入口か
ら、ステップ２で設定した２線ケーブルの出口に至るま
でのケーブル経由点の座標を設定する。ここで、このケ
ーブル経由点の座標は、２線ケーブルの中心位置（２つ
の芯線の中間位置）で定義されている。なお、ケーブル
経由点が指定されないときには、ケーブルの入口から出
口まで一直線にケーブルが配置されることになる。
【００４４】続いて、ステップ４で、ユーザと対話する
ことで、ケーブルライブラリ２００に登録されるケーブ
ルの名称を設定することで、使用する２線ケーブルを設
定し、続くステップ５で、次のケーブルの設定要求があ
るのか否かを判断して、設定要求があることを判断する
ときには、ステップ１に戻り、設定要求がないことを判
断するときには、「ケーブル配置処理」を終了する。
ことで、ケーブルライブラリ２００に登録されるケーブ
ルの名称を設定することで、使用する２線ケーブルを設
定し、続くステップ５で、次のケーブルの設定要求があ
るのか否かを判断して、設定要求があることを判断する
ときには、ステップ１に戻り、設定要求がないことを判
断するときには、「ケーブル配置処理」を終了する。
【００４５】なお、２線ケーブルの入口及び出口の座標
については、入力データファイル１００に予め登録され
ていることがあり、このときには、ユーザは、ケーブル
経由点の座標のみを設定することになる。
については、入力データファイル１００に予め登録され
ていることがあり、このときには、ユーザは、ケーブル
経由点の座標のみを設定することになる。
【００４６】このようにして、本発明の電磁界強度算出
装置１は、「ケーブル配置処理」に入ると、ユーザと対
話することで、図１０（ａ）に示すように、シミュレー
ション対象となる電気回路装置に２線ケーブルを配置す
るのである。この２線ケーブルは、実際には、図１０
（ｂ）に示すように、太さを持つものである。ここで、
図中に示すＡ及びＢは、ユーザの設定した２線ケーブル
の入口の座標、Ｃ及びＤは、ユーザの設定した２線ケー
ブルの出口の座標を示している。
装置１は、「ケーブル配置処理」に入ると、ユーザと対
話することで、図１０（ａ）に示すように、シミュレー
ション対象となる電気回路装置に２線ケーブルを配置す
るのである。この２線ケーブルは、実際には、図１０
（ｂ）に示すように、太さを持つものである。ここで、
図中に示すＡ及びＢは、ユーザの設定した２線ケーブル
の入口の座標、Ｃ及びＤは、ユーザの設定した２線ケー
ブルの出口の座標を示している。
【００４７】以下、図１０（ａ）に示すように、２線ケ
ーブルの中心を通り、２線ケーブルの各経由点を結ぶ線
分をリンクと呼ぶことにし、２線ケーブルの入口側のリ
ンク端点をリンクの始点と呼び、２線ケーブルの出口側
のリンク端点をリンクの終点と呼ぶことにする。
ーブルの中心を通り、２線ケーブルの各経由点を結ぶ線
分をリンクと呼ぶことにし、２線ケーブルの入口側のリ
ンク端点をリンクの始点と呼び、２線ケーブルの出口側
のリンク端点をリンクの終点と呼ぶことにする。
【００４８】図１１及び図１２に、ケーブル経由点の座
標設定に用いるディスプレイ画面の一実施例を図示す
る。このディスプレイ画面に示すように、本発明の電磁
界強度算出装置１では、設定された２線ケーブルの入口
／出口／経由点の３次元座標位置をユーザが正確に把握
できるようにするために、４つの方向（ＸＹ平面／ＸＺ
平面／ＹＺ平面／３次元図）からの見たときの図を使っ
て、これらの各点の座標を表示する構成を採っている。
ここで、図１１は、２線ケーブルの入口／出口の座標が
設定されるときに表示するディスプレイ画面であり、図
１２は、その後、ケーブル経由点の座標が設定されると
きに表示するディスプレイ画面である。
標設定に用いるディスプレイ画面の一実施例を図示す
る。このディスプレイ画面に示すように、本発明の電磁
界強度算出装置１では、設定された２線ケーブルの入口
／出口／経由点の３次元座標位置をユーザが正確に把握
できるようにするために、４つの方向（ＸＹ平面／ＸＺ
平面／ＹＺ平面／３次元図）からの見たときの図を使っ
て、これらの各点の座標を表示する構成を採っている。
ここで、図１１は、２線ケーブルの入口／出口の座標が
設定されるときに表示するディスプレイ画面であり、図
１２は、その後、ケーブル経由点の座標が設定されると
きに表示するディスプレイ画面である。
【００４９】このようにして、「ケーブル配置処理」に
従い、ユーザと対話することで、２線ケーブルの入口／
出口／経由点の座標を設定するとともに、２線ケーブル
の名称を設定すると、本発明の電磁界強度算出装置１
は、続いて、「リンク情報決定処理」を実行する。
従い、ユーザと対話することで、２線ケーブルの入口／
出口／経由点の座標を設定するとともに、２線ケーブル
の名称を設定すると、本発明の電磁界強度算出装置１
は、続いて、「リンク情報決定処理」を実行する。
【００５０】このリンク情報とは、リンクの両端の座標
で定義されるものであって、図１３に示すように、「ケ
ーブル配置処理」で設定されたケーブル経由点（ケーブ
ルの入口／出口は、その座標の中点）の座標と、その位
置での２線ケーブルの芯線位置（ケーブルの入口／出口
は、その座標）とで定義される。
で定義されるものであって、図１３に示すように、「ケ
ーブル配置処理」で設定されたケーブル経由点（ケーブ
ルの入口／出口は、その座標の中点）の座標と、その位
置での２線ケーブルの芯線位置（ケーブルの入口／出口
は、その座標）とで定義される。
【００５１】すなわち、リンク情報は、図１４に示すよ
うに、各ケーブル経由点の座標と、各ケーブル経由点で
リンクに直交する正負の２つの法線方向（正負はケーブ
ル進行方向に対して右手系で定義する）に設けられて、
ケーブル経由点からケーブルの導体間距離の１／２だけ
離れた位置として定義される２線ケーブルの芯線位置の
座標から求められることになる。「リンク情報決定処
理」とは、このリンク情報を求める処理である。
うに、各ケーブル経由点の座標と、各ケーブル経由点で
リンクに直交する正負の２つの法線方向（正負はケーブ
ル進行方向に対して右手系で定義する）に設けられて、
ケーブル経由点からケーブルの導体間距離の１／２だけ
離れた位置として定義される２線ケーブルの芯線位置の
座標から求められることになる。「リンク情報決定処
理」とは、このリンク情報を求める処理である。
【００５２】「ケーブル配置処理」に従って２線ケーブ
ルの入口／出口／経由点が設定されることで、「リンク
情報決定処理」に入ると、本発明の電磁界強度算出装置
１は、図５及び図６の処理フローに示すように、先ず最
初に、電気回路装置に配置されるケーブルの中からリン
ク情報の決定対象となるケーブルを１つ選択し、続くス
テップ２で、全てのケーブルの選択を判断するときに
は、「リンク情報決定処理」を終了する。
ルの入口／出口／経由点が設定されることで、「リンク
情報決定処理」に入ると、本発明の電磁界強度算出装置
１は、図５及び図６の処理フローに示すように、先ず最
初に、電気回路装置に配置されるケーブルの中からリン
ク情報の決定対象となるケーブルを１つ選択し、続くス
テップ２で、全てのケーブルの選択を判断するときに
は、「リンク情報決定処理」を終了する。
【００５３】一方、ステップ２で、全てのケーブルの選
択が完了していないことを判断するとき、すなわち、ス
テップ１でリンク情報の決定対象となるケーブルを選択
できたことを判断するときには、ステップ３に進んで、
選択したケーブルの名称をキーにしてケーブルライブラ
リ２００を参照することで、選択したケーブルが平行線
ケーブルか撚り線ケーブルのいずれかであるのかを判断
して、そのいずれでもないことを判断するとき、すなわ
ち、２線ケーブルでないことを判断するときには、別に
用意されるモデル化アルゴリズムに従ってモデル化処理
を実行することになるので、以下の処理を行わずにステ
ップ１に戻る。
択が完了していないことを判断するとき、すなわち、ス
テップ１でリンク情報の決定対象となるケーブルを選択
できたことを判断するときには、ステップ３に進んで、
選択したケーブルの名称をキーにしてケーブルライブラ
リ２００を参照することで、選択したケーブルが平行線
ケーブルか撚り線ケーブルのいずれかであるのかを判断
して、そのいずれでもないことを判断するとき、すなわ
ち、２線ケーブルでないことを判断するときには、別に
用意されるモデル化アルゴリズムに従ってモデル化処理
を実行することになるので、以下の処理を行わずにステ
ップ１に戻る。
【００５４】一方、ステップ３で、選択したケーブルが
２線ケーブル（平行線ケーブル／撚り線ケーブル）であ
ることを判断するときには、ステップ４に進んで、選択
した２線ケーブルの名称をキーにしてケーブルライブラ
リ２００を参照することで、選択した２線ケーブルの持
つ導体間距離δを読み取り、「ケーブル配置処理」で設
定されたその２線ケーブルの入口及び出口の座標から算
出される導体間距離が、その読み出した導体間距離δと
整合するのか否かを判断して、整合しないことを判断す
るときには、ステップ５に進んで、エラーメッセージを
出力してから、以下の処理を行わずにステップ１に戻
る。
２線ケーブル（平行線ケーブル／撚り線ケーブル）であ
ることを判断するときには、ステップ４に進んで、選択
した２線ケーブルの名称をキーにしてケーブルライブラ
リ２００を参照することで、選択した２線ケーブルの持
つ導体間距離δを読み取り、「ケーブル配置処理」で設
定されたその２線ケーブルの入口及び出口の座標から算
出される導体間距離が、その読み出した導体間距離δと
整合するのか否かを判断して、整合しないことを判断す
るときには、ステップ５に進んで、エラーメッセージを
出力してから、以下の処理を行わずにステップ１に戻
る。
【００５５】一方、ステップ４で、「ケーブル配置処
理」で設定された２線ケーブルの入口及び出口の座標か
ら算出される導体間距離に矛盾のないことを判断すると
きには、ステップ６に進んで、その入口の座標（２つの
座標からなっている）の中点を算出するとともに、その
出口の座標（２つの座標からなっている）の中点を算出
する。
理」で設定された２線ケーブルの入口及び出口の座標か
ら算出される導体間距離に矛盾のないことを判断すると
きには、ステップ６に進んで、その入口の座標（２つの
座標からなっている）の中点を算出するとともに、その
出口の座標（２つの座標からなっている）の中点を算出
する。
【００５６】続いて、ステップ７で、選択した２線ケー
ブルに設定されたケーブル経由点の数に“１”を加算す
ることで、その２線ケーブルの入口から出口に到達する
までのリンク数Ｎを算出し、続くステップ８で、そのリ
ンク数が“１”であるのか否かを判断して、“１”であ
ることを判断するとき、すなわち、ケーブル経由点が設
定されなかったことを判断するときには、ステップ９に
進んで、「ケーブル配置処理」で設定された２線ケーブ
ルの入口及び出口の座標と、ステップ６で算出したそれ
らの中点の座標とを、選択した２線ケーブルのリンク情
報として設定してステップ１に戻る。
ブルに設定されたケーブル経由点の数に“１”を加算す
ることで、その２線ケーブルの入口から出口に到達する
までのリンク数Ｎを算出し、続くステップ８で、そのリ
ンク数が“１”であるのか否かを判断して、“１”であ
ることを判断するとき、すなわち、ケーブル経由点が設
定されなかったことを判断するときには、ステップ９に
進んで、「ケーブル配置処理」で設定された２線ケーブ
ルの入口及び出口の座標と、ステップ６で算出したそれ
らの中点の座標とを、選択した２線ケーブルのリンク情
報として設定してステップ１に戻る。
【００５７】一方、ステップ８で、リンク数が“２”以
上であることを判断するときには、ステップ１０に進ん
で、２線ケーブルの入口に接続するリンクである開始リ
ンクを現在リンクとして設定する。続いて、ステップ１
１（図６の処理フロー）で、現在リンクから、現在リン
クの方向ベクトルと異なる方向ベクトルを持つリンクま
で辿ることで、現在リンクに直交する法線ベクトルを求
める。
上であることを判断するときには、ステップ１０に進ん
で、２線ケーブルの入口に接続するリンクである開始リ
ンクを現在リンクとして設定する。続いて、ステップ１
１（図６の処理フロー）で、現在リンクから、現在リン
クの方向ベクトルと異なる方向ベクトルを持つリンクま
で辿ることで、現在リンクに直交する法線ベクトルを求
める。
【００５８】すなわち、図１３及び図１４で説明したよ
うに、リンク情報とは、各ケーブル経由点の座標と、各
ケーブル経由点でリンクに直交する正負の２つの法線方
向に設けられて、ケーブル経由点からケーブルの導体間
距離の１／２だけ離れた位置として定義される２線ケー
ブルの芯線位置の座標とで定義されている。この法線方
向は、図１５に示すように、リンクｉの持つ方向ベクト
ルｌi と、それに接続するリンクｊの持つ方向ベクトル
ｌj とが同一方向とならないときに、その外積として決
定できる。これから、ステップ１１では、この法線ベク
トルが求められるまで、リンクを辿っていく処理を行う
のである。
うに、リンク情報とは、各ケーブル経由点の座標と、各
ケーブル経由点でリンクに直交する正負の２つの法線方
向に設けられて、ケーブル経由点からケーブルの導体間
距離の１／２だけ離れた位置として定義される２線ケー
ブルの芯線位置の座標とで定義されている。この法線方
向は、図１５に示すように、リンクｉの持つ方向ベクト
ルｌi と、それに接続するリンクｊの持つ方向ベクトル
ｌj とが同一方向とならないときに、その外積として決
定できる。これから、ステップ１１では、この法線ベク
トルが求められるまで、リンクを辿っていく処理を行う
のである。
【００５９】続いて、ステップ１２で、法線ベクトルの
求められるリンクまで辿れたのか否かを判断して、辿れ
たことを判断するときには、ステップ１３に進んで、現
在リンクと辿り先リンクとの間にリンクが存在するのか
否かを判断して、リンクが存在しないことを判断すると
きには、ステップ１４に進んで、法線ベクトルが求めら
れることで、現在リンクのリンク情報が求められること
になったことに対応して、現在リンクのリンク情報を求
める。
求められるリンクまで辿れたのか否かを判断して、辿れ
たことを判断するときには、ステップ１３に進んで、現
在リンクと辿り先リンクとの間にリンクが存在するのか
否かを判断して、リンクが存在しないことを判断すると
きには、ステップ１４に進んで、法線ベクトルが求めら
れることで、現在リンクのリンク情報が求められること
になったことに対応して、現在リンクのリンク情報を求
める。
【００６０】具体的には、現在リンクの進行先端点での
法線ベクトルを求めて、その法線ベクトルの方向に導体
間距離の１／２だけ離れた位置の座標を求めることで、
現在リンクの進行先端点でのリンク情報を求める。ここ
で、現在リンクの進行元端点でのリンク情報は、現在リ
ンクが開始リンクであるときには、２線ケーブルの入口
の座標とステップ５で求めた中点座標とをコピーし、現
在リンクが開始リンクでないときには、現在リンクに接
続するリンクの進行先端点でのリンク情報と同一となる
ので、それをコピーすることで求める。また、リンク情
報を構成する中点の座標は、ケーブル経由点の座標をコ
ピーすることで求める。
法線ベクトルを求めて、その法線ベクトルの方向に導体
間距離の１／２だけ離れた位置の座標を求めることで、
現在リンクの進行先端点でのリンク情報を求める。ここ
で、現在リンクの進行元端点でのリンク情報は、現在リ
ンクが開始リンクであるときには、２線ケーブルの入口
の座標とステップ５で求めた中点座標とをコピーし、現
在リンクが開始リンクでないときには、現在リンクに接
続するリンクの進行先端点でのリンク情報と同一となる
ので、それをコピーすることで求める。また、リンク情
報を構成する中点の座標は、ケーブル経由点の座標をコ
ピーすることで求める。
【００６１】一方、ステップ１３で、現在リンクと辿り
先リンクとの間にリンクが存在することを判断すると
き、すなわち、現在リンクと辿り先リンクとの間に、現
在リンクの方向ベクトルと同一の方向ベクトルを持つリ
ンクが存在することを判断するときには、ステップ１５
に進んで、法線ベクトルが求められることで、現在リン
クのリンク情報と、現在リンクと辿り先リンクとの間に
存在するリンクのリンク情報とが求められることになっ
たことに対応して、それらのリンク情報を求める。
先リンクとの間にリンクが存在することを判断すると
き、すなわち、現在リンクと辿り先リンクとの間に、現
在リンクの方向ベクトルと同一の方向ベクトルを持つリ
ンクが存在することを判断するときには、ステップ１５
に進んで、法線ベクトルが求められることで、現在リン
クのリンク情報と、現在リンクと辿り先リンクとの間に
存在するリンクのリンク情報とが求められることになっ
たことに対応して、それらのリンク情報を求める。
【００６２】具体的には、辿り先リンクの始点での法線
ベクトルを求めて、その法線ベクトルを、現在リンクの
進行先端点での法線ベクトルとして用いるとともに、現
在リンクと辿り先リンクとの間に存在するリンクの持つ
２つの端点での法線ベクトルとして用いて、その法線ベ
クトルの方向に導体間距離の１／２だけ離れた位置の座
標を求めることで行う。ここで、現在リンクの進行元端
点でのリンク情報と、リンク情報を構成する中点の座標
は、上述したコピー処理により行う。
ベクトルを求めて、その法線ベクトルを、現在リンクの
進行先端点での法線ベクトルとして用いるとともに、現
在リンクと辿り先リンクとの間に存在するリンクの持つ
２つの端点での法線ベクトルとして用いて、その法線ベ
クトルの方向に導体間距離の１／２だけ離れた位置の座
標を求めることで行う。ここで、現在リンクの進行元端
点でのリンク情報と、リンク情報を構成する中点の座標
は、上述したコピー処理により行う。
【００６３】このようにして、ステップ１３／ステップ
１５でリンク情報を求めると、ステップ１６に進んで、
辿り先リンクを現在リンクとして設定し、続くステップ
１７で、現在リンクが２線ケーブルの出口に接続するリ
ンクである終了リンクであるのか否かを判断して、終了
リンクでないことを判断するときには、ステップ１１に
戻り、終了リンクであることを判断するときには、ステ
ップ１８に進んで、終了リンクに接続するリンクの進行
先端点でのリンク情報をコピーするとともに、２線ケー
ブルの出口の座標とステップ５で求めた中点座標とをコ
ピーすることで終了リンクのリンク情報を求めて、ステ
ップ１に戻る。
１５でリンク情報を求めると、ステップ１６に進んで、
辿り先リンクを現在リンクとして設定し、続くステップ
１７で、現在リンクが２線ケーブルの出口に接続するリ
ンクである終了リンクであるのか否かを判断して、終了
リンクでないことを判断するときには、ステップ１１に
戻り、終了リンクであることを判断するときには、ステ
ップ１８に進んで、終了リンクに接続するリンクの進行
先端点でのリンク情報をコピーするとともに、２線ケー
ブルの出口の座標とステップ５で求めた中点座標とをコ
ピーすることで終了リンクのリンク情報を求めて、ステ
ップ１に戻る。
【００６４】一方、ステップ１２で、法線ベクトルの求
められるリンクまで辿れないことを判断するとき、すな
わち、終了リンクまで方向ベクトルが変化しないことを
判断するときには、ステップ１９に進んで、２線ケーブ
ルの出口の座標の規定する法線ベクトルに従って、現在
リンクのリンク情報と、終了リンクのリンク情報と、現
在リンクと終了リンクとの間に存在するリンクのリンク
情報を求めて、ステップ１に戻る。
められるリンクまで辿れないことを判断するとき、すな
わち、終了リンクまで方向ベクトルが変化しないことを
判断するときには、ステップ１９に進んで、２線ケーブ
ルの出口の座標の規定する法線ベクトルに従って、現在
リンクのリンク情報と、終了リンクのリンク情報と、現
在リンクと終了リンクとの間に存在するリンクのリンク
情報を求めて、ステップ１に戻る。
【００６５】具体的には、２線ケーブルの出口の座標の
規定する法線ベクトルを求めて、その法線ベクトルを、
現在リンクの進行先端点での法線ベクトルとし、終了リ
ンクの進行元端点での法線ベクトルとして用いるととも
に、現在リンクと終了リンクとの間に存在するリンクの
持つ２つの端点での法線ベクトルとして用いて、その法
線ベクトルの方向に導体間距離の１／２だけ離れた位置
の座標を求めることで行う。ここで、終了リンクの進行
先端点でのリンク情報は、２線ケーブルの出口の座標と
ステップ５で求めた中点座標とをコピーすることで行
う。また、現在リンクの進行元端点でのリンク情報と、
リンク情報を構成する中点の座標は、上述したコピー処
理により行う。
規定する法線ベクトルを求めて、その法線ベクトルを、
現在リンクの進行先端点での法線ベクトルとし、終了リ
ンクの進行元端点での法線ベクトルとして用いるととも
に、現在リンクと終了リンクとの間に存在するリンクの
持つ２つの端点での法線ベクトルとして用いて、その法
線ベクトルの方向に導体間距離の１／２だけ離れた位置
の座標を求めることで行う。ここで、終了リンクの進行
先端点でのリンク情報は、２線ケーブルの出口の座標と
ステップ５で求めた中点座標とをコピーすることで行
う。また、現在リンクの進行元端点でのリンク情報と、
リンク情報を構成する中点の座標は、上述したコピー処
理により行う。
【００６６】このようにして、本発明の電磁界強度算出
装置１は、「リンク情報決定処理」に入ると、２線ケー
ブルの入口から出口の間を接続するリンクの持つリンク
情報を求めるのである。
装置１は、「リンク情報決定処理」に入ると、２線ケー
ブルの入口から出口の間を接続するリンクの持つリンク
情報を求めるのである。
【００６７】次に、図２の処理フローのステップ３で実
行する「ケーブルモデル化処理」のについて説明する。
本発明の電磁界強度算出装置１は、「ケーブルモデル化
処理」に入ると、図７の処理フローに示すように、先ず
最初に、ステップ１で、モーメント法の区分長Ｌを決め
る。この区分長Ｌは、シミュレーション対象となる最大
周波数ｆと、光速度ｃとから、 Ｌ＝（１／１０）×（ｆ／ｃ） の式に従って求めることになる。
行する「ケーブルモデル化処理」のについて説明する。
本発明の電磁界強度算出装置１は、「ケーブルモデル化
処理」に入ると、図７の処理フローに示すように、先ず
最初に、ステップ１で、モーメント法の区分長Ｌを決め
る。この区分長Ｌは、シミュレーション対象となる最大
周波数ｆと、光速度ｃとから、 Ｌ＝（１／１０）×（ｆ／ｃ） の式に従って求めることになる。
【００６８】続いて、ステップ２で、ユーザと対話する
ことで、撚り線ケーブルの１ピッチ当たりの区分数ｍを
決める。モーメント法を用いる場合、電気回路装置の各
部分を区分長Ｌの大きさでメッシュ化していくことにな
るが、撚り線ケーブルの撚りピッチψは、通常、この区
分長Ｌよりもかなり小さい。従って、区分長Ｌに従って
撚り線ケーブルをメッシュ化すると、撚り線の形態が消
されてしまうことになる。これから、撚り線ケーブルに
ついては、１ピッチ当たり２４分割というような区分数
ｍに従って、メッシュ化していくことになるので、この
ステップ２では、ユーザと対話することでこの区分数ｍ
を決定する。ここで、この区分数ｍは、偶数の形態で決
定する。
ことで、撚り線ケーブルの１ピッチ当たりの区分数ｍを
決める。モーメント法を用いる場合、電気回路装置の各
部分を区分長Ｌの大きさでメッシュ化していくことにな
るが、撚り線ケーブルの撚りピッチψは、通常、この区
分長Ｌよりもかなり小さい。従って、区分長Ｌに従って
撚り線ケーブルをメッシュ化すると、撚り線の形態が消
されてしまうことになる。これから、撚り線ケーブルに
ついては、１ピッチ当たり２４分割というような区分数
ｍに従って、メッシュ化していくことになるので、この
ステップ２では、ユーザと対話することでこの区分数ｍ
を決定する。ここで、この区分数ｍは、偶数の形態で決
定する。
【００６９】続いて、ステップ３で、ケーブルライブラ
リ２００から、２線ケーブルの属性情報を読み込む。続
いて、ステップ４で、電気回路装置に配置される２線ケ
ーブルの中からモデル化対象となるものを１つ選択し、
続くステップ２で、全ての２線ケーブルの選択を判断す
るときには、「ケーブルモデル化処理」を終了する。
リ２００から、２線ケーブルの属性情報を読み込む。続
いて、ステップ４で、電気回路装置に配置される２線ケ
ーブルの中からモデル化対象となるものを１つ選択し、
続くステップ２で、全ての２線ケーブルの選択を判断す
るときには、「ケーブルモデル化処理」を終了する。
【００７０】一方、ステップ４で、全ての２線ケーブル
の選択が完了していないことを判断するとき、すなわ
ち、ステップ１でモデル化対象となる２線ケーブルを選
択できたことを判断するときには、ステップ５に進ん
で、選択した２線ケーブルが平行線ケーブルであるの
か、撚り線ケーブルであるのかを判断する。
の選択が完了していないことを判断するとき、すなわ
ち、ステップ１でモデル化対象となる２線ケーブルを選
択できたことを判断するときには、ステップ５に進ん
で、選択した２線ケーブルが平行線ケーブルであるの
か、撚り線ケーブルであるのかを判断する。
【００７１】このステップ６の判断処理により、選択し
た２線ケーブルが平行ケーブルであることを判断すると
きには、ステップ７に進んで、その平行線ケーブルの持
つ非平行リンク（法線方向が変化しないリンク、法線方
向が変化しないリンクが連続するときには、それをまと
めたリンク）の中から、非平行リンクを１つ選択する。
た２線ケーブルが平行ケーブルであることを判断すると
きには、ステップ７に進んで、その平行線ケーブルの持
つ非平行リンク（法線方向が変化しないリンク、法線方
向が変化しないリンクが連続するときには、それをまと
めたリンク）の中から、非平行リンクを１つ選択する。
【００７２】そして、続いて、ステップ８で、全ての非
平行リンクの選択が完了したのか否かを判断して、完了
したことを判断するときには、次の２線ケーブルのモデ
ル化を行うべくステップ４に戻り、非平行リンクを選択
できたことを判断するときには、ステップ９に進んで、
平行線ケーブルのモデル化アルゴリズムを使って、選択
した非平行リンクのモーメント法モデルを生成してか
ら、ステップ４で選択した平行線ケーブルの持つ次の非
平行リンクのモデル化を行うべくステップ７に戻ってい
く。
平行リンクの選択が完了したのか否かを判断して、完了
したことを判断するときには、次の２線ケーブルのモデ
ル化を行うべくステップ４に戻り、非平行リンクを選択
できたことを判断するときには、ステップ９に進んで、
平行線ケーブルのモデル化アルゴリズムを使って、選択
した非平行リンクのモーメント法モデルを生成してか
ら、ステップ４で選択した平行線ケーブルの持つ次の非
平行リンクのモデル化を行うべくステップ７に戻ってい
く。
【００７３】一方、このステップ６の判断処理により、
選択した２線ケーブルが撚り線ーブルであることを判断
するときには、ステップ１０に進んで、その撚り線ケー
ブルの持つ非平行リンクの中から、非平行リンクを１つ
選択する。
選択した２線ケーブルが撚り線ーブルであることを判断
するときには、ステップ１０に進んで、その撚り線ケー
ブルの持つ非平行リンクの中から、非平行リンクを１つ
選択する。
【００７４】そして、続いて、ステップ１１で、全ての
非平行リンクの選択が完了したのか否かを判断して、完
了したことを判断するときには、次の２線ケーブルのモ
デル化を行うべくステップ４に戻り、非平行リンクを選
択できたことを判断するときには、ステップ１２に進ん
で、撚り線ケーブルのモデル化アルゴリズムを使って、
選択した非平行リンクのモーメント法モデルを生成して
から、ステップ４で選択した撚り線ケーブルの持つ次の
非平行リンクのモデル化を行うべくステップ１０に戻っ
ていく。
非平行リンクの選択が完了したのか否かを判断して、完
了したことを判断するときには、次の２線ケーブルのモ
デル化を行うべくステップ４に戻り、非平行リンクを選
択できたことを判断するときには、ステップ１２に進ん
で、撚り線ケーブルのモデル化アルゴリズムを使って、
選択した非平行リンクのモーメント法モデルを生成して
から、ステップ４で選択した撚り線ケーブルの持つ次の
非平行リンクのモデル化を行うべくステップ１０に戻っ
ていく。
【００７５】次に、図８及び図９に示す処理フローにつ
いて説明する。ここで、図８の処理フローは、平行線ケ
ーブルのモデル化アルゴリズムの処理フロー（図７の処
理フローのステップ９の詳細な処理フロー）であり、図
９の処理フローは、撚り線ケーブルのモデル化アルゴリ
ズムの処理フロー（図７の処理フローのステップ１２の
詳細な処理フロー）である。
いて説明する。ここで、図８の処理フローは、平行線ケ
ーブルのモデル化アルゴリズムの処理フロー（図７の処
理フローのステップ９の詳細な処理フロー）であり、図
９の処理フローは、撚り線ケーブルのモデル化アルゴリ
ズムの処理フロー（図７の処理フローのステップ１２の
詳細な処理フロー）である。
【００７６】最初に、図８の処理フローについて説明す
る。すなわち、図７の処理フローのステップ９に従っ
て、平行線ケーブルの持つ非平行リンクのモデル化処理
に入ると、図８の処理フローに示すように、先ず最初
に、ステップ１で、非平行リンクのリンク長をモーメン
ト法の区分長Ｌで割り算することで、偶数を示す区分数
Ｍを求める。このとき、割り算値が奇数となるときに
は、メッシュを小さくする方がモーメント法の精度が高
まるので、その割り算値に“１”を加算する形で偶数を
示す区分数Ｍを求める。
る。すなわち、図７の処理フローのステップ９に従っ
て、平行線ケーブルの持つ非平行リンクのモデル化処理
に入ると、図８の処理フローに示すように、先ず最初
に、ステップ１で、非平行リンクのリンク長をモーメン
ト法の区分長Ｌで割り算することで、偶数を示す区分数
Ｍを求める。このとき、割り算値が奇数となるときに
は、メッシュを小さくする方がモーメント法の精度が高
まるので、その割り算値に“１”を加算する形で偶数を
示す区分数Ｍを求める。
【００７７】続いて、ステップ２で、非平行リンクの両
端の法線ベクトルのなす捩れ角度θを求める。すなわ
ち、図１６に示すように、非平行リンクの始点（Ａ，
Ｂ）での法線ベクトルと、終点（Ｃ，Ｄ）での法線ベク
トルとのなす捩れ角度θを求めるのである。
端の法線ベクトルのなす捩れ角度θを求める。すなわ
ち、図１６に示すように、非平行リンクの始点（Ａ，
Ｂ）での法線ベクトルと、終点（Ｃ，Ｄ）での法線ベク
トルとのなす捩れ角度θを求めるのである。
【００７８】続いて、ステップ３で、「θ＝９０度」で
あるのか否かを判断して、「θ＝９０度」であることを
判断するときには、ステップ４に進んで、前回の「θ＝
９０度」のときに、右回りに曲げているときには、今回
は左回りになるようにと曲げ方向を決定し、前回の「θ
＝９０度」のときに、左回りに曲げているときには、今
回は右回りになるようにと曲げ方向を決定する。
あるのか否かを判断して、「θ＝９０度」であることを
判断するときには、ステップ４に進んで、前回の「θ＝
９０度」のときに、右回りに曲げているときには、今回
は左回りになるようにと曲げ方向を決定し、前回の「θ
＝９０度」のときに、左回りに曲げているときには、今
回は右回りになるようにと曲げ方向を決定する。
【００７９】すなわち、「θ＜９０度」であるときに
は、図１７に示すように、平行線ケーブルの曲げが小さ
くなるようにと、ＡはＤに接続され、ＢはＥに接続され
るというように、始点（Ａ，Ｂ）と終点（Ｄ，Ｅ）との
対応をとって曲げることになるが、「θ＝９０度」のと
きには、右方向に曲げることもできるし、左方向に曲げ
ることもできる。このとき、どちらか一方に曲げ方向を
規定すると、平行線ケーブルがどんどん捩れていくこと
になるので、ここでは、前回の曲げ方向と逆方向に曲げ
ることを決定するのである。
は、図１７に示すように、平行線ケーブルの曲げが小さ
くなるようにと、ＡはＤに接続され、ＢはＥに接続され
るというように、始点（Ａ，Ｂ）と終点（Ｄ，Ｅ）との
対応をとって曲げることになるが、「θ＝９０度」のと
きには、右方向に曲げることもできるし、左方向に曲げ
ることもできる。このとき、どちらか一方に曲げ方向を
規定すると、平行線ケーブルがどんどん捩れていくこと
になるので、ここでは、前回の曲げ方向と逆方向に曲げ
ることを決定するのである。
【００８０】続いて、ステップ５で、法線ベクトルの捩
れ角度θを区分数Ｍで割り算することで、捩れ角度θを
補正するための補正角度φを求める。すなわち、区分数
Ｍで区分するときに、どの位ずつ回転したいったらよい
のかということを示す補正角度φを算出するのである。
れ角度θを区分数Ｍで割り算することで、捩れ角度θを
補正するための補正角度φを求める。すなわち、区分数
Ｍで区分するときに、どの位ずつ回転したいったらよい
のかということを示す補正角度φを算出するのである。
【００８１】続いて、ステップ６で、区分数Ｍと補正角
度φとから、平行線ケーブルの芯線の区分位置を決定す
るとともに、コンデンサの挿入位置を決定する。この芯
線の区分位置の決定処理は、非平行リンクのリンク長を
Ｘで表すならば、始点から終点に向けて、最初はＸ／
（２Ｍ）、その後はＸ／Ｍずつ進めていきながら、それ
と同期させて、最初はφ／２、その後はφずつ回転させ
ていって、その回転方向の指す方向に向けて導体間距離
の１／２だけ離れた位置の座標を特定することで行う。
また、コンデンサの挿入位置の決定処理は、図１８に示
すように、最初はＸ／（２Ｍ）の所、その後はＸ／Ｍず
つ進めた所を挿入位置として特定することで行う。ここ
で、図１８の黒丸は芯線の区分位置を表している。
度φとから、平行線ケーブルの芯線の区分位置を決定す
るとともに、コンデンサの挿入位置を決定する。この芯
線の区分位置の決定処理は、非平行リンクのリンク長を
Ｘで表すならば、始点から終点に向けて、最初はＸ／
（２Ｍ）、その後はＸ／Ｍずつ進めていきながら、それ
と同期させて、最初はφ／２、その後はφずつ回転させ
ていって、その回転方向の指す方向に向けて導体間距離
の１／２だけ離れた位置の座標を特定することで行う。
また、コンデンサの挿入位置の決定処理は、図１８に示
すように、最初はＸ／（２Ｍ）の所、その後はＸ／Ｍず
つ進めた所を挿入位置として特定することで行う。ここ
で、図１８の黒丸は芯線の区分位置を表している。
【００８２】続いて、ステップ７で、平行線ケーブルの
単位長当たりの容量を求める。この容量は、ケーブルラ
イブラリ２００に登録されているときには、それを参照
することで実現されることになるが、ケーブルライブラ
リ２００に登録されていないときには、数式処理に従っ
て求めることになる。
単位長当たりの容量を求める。この容量は、ケーブルラ
イブラリ２００に登録されているときには、それを参照
することで実現されることになるが、ケーブルライブラ
リ２００に登録されていないときには、数式処理に従っ
て求めることになる。
【００８３】すなわち、電磁気学で知られているよう
に、直径２ａの導線が距離Ｄで平行に配設されていると
きに、空気中での単位長当たりのインダクタンスＬ
0 は、 Ｌ0 ＝（μ0 ／π）ｌｎ〔（Ｄ＋（Ｄ2 −４ａ）1/2 ）／２ａ〕 ≒（μ0 ／π）ｌｎ（Ｄ／ａ） であり、空気中での単位長当たりの容量Ｃ0 は、 Ｃ0 ＝（πε0 ）／ｌｎ〔（Ｄ＋（Ｄ2 −４ａ）1/2 ）／２ａ〕 ≒（πε0 ）／ｌｎ（Ｄ／ａ） である。これから、空気中での特性インピーダンスＺ0
は、 Ｚ0 ≒（Ｌ0 ／Ｃ0 ）1/2 ＝（１／π)(μ0 ／ε0 ）1/2 ｌｎ（Ｄ／ａ） となる。
に、直径２ａの導線が距離Ｄで平行に配設されていると
きに、空気中での単位長当たりのインダクタンスＬ
0 は、 Ｌ0 ＝（μ0 ／π）ｌｎ〔（Ｄ＋（Ｄ2 −４ａ）1/2 ）／２ａ〕 ≒（μ0 ／π）ｌｎ（Ｄ／ａ） であり、空気中での単位長当たりの容量Ｃ0 は、 Ｃ0 ＝（πε0 ）／ｌｎ〔（Ｄ＋（Ｄ2 −４ａ）1/2 ）／２ａ〕 ≒（πε0 ）／ｌｎ（Ｄ／ａ） である。これから、空気中での特性インピーダンスＺ0
は、 Ｚ0 ≒（Ｌ0 ／Ｃ0 ）1/2 ＝（１／π)(μ0 ／ε0 ）1/2 ｌｎ（Ｄ／ａ） となる。
【００８４】一方、ケーブルライブラリ２００に登録さ
れる特性インピーダンスＺは、芯線を被覆する誘電体を
考慮した実測値である。この誘電体の実効誘電率εreff
は、空気中での特性インピーダンスＺ0 及び実測値Ｚと
の間に、 Ｚ／Ｚ0 ＝（εreff）1/2 の関係を持つ。これから、誘電体の実効誘電率ε
reffは、計算で算出される空気中での特性インピーダン
スＺ0 と、ケーブルライブラリ２００に登録される特性
インピーダンスＺとから、求めることが可能である。こ
の実効誘電率εreffが求まると、誘電体で被覆されたケ
ーブルの単位長さ当たりの容量Ｃは、 Ｃ＝εreffＣ0 −Ｃ0 ＝Ｃ0(εreff−１） に従って求めることができる。
れる特性インピーダンスＺは、芯線を被覆する誘電体を
考慮した実測値である。この誘電体の実効誘電率εreff
は、空気中での特性インピーダンスＺ0 及び実測値Ｚと
の間に、 Ｚ／Ｚ0 ＝（εreff）1/2 の関係を持つ。これから、誘電体の実効誘電率ε
reffは、計算で算出される空気中での特性インピーダン
スＺ0 と、ケーブルライブラリ２００に登録される特性
インピーダンスＺとから、求めることが可能である。こ
の実効誘電率εreffが求まると、誘電体で被覆されたケ
ーブルの単位長さ当たりの容量Ｃは、 Ｃ＝εreffＣ0 −Ｃ0 ＝Ｃ0(εreff−１） に従って求めることができる。
【００８５】これから、このステップ７では、先ず最初
に、公式に従って、空気中での単位長当たりの容量Ｃ0
と、空気中での特性インピーダンスＺ0 を算出し、次
に、この特性インピーダンスＺ0 と、ケーブルライブラ
リ２００に登録される特性インピーダンスＺとから誘電
体の実効誘電率εreffを算出し、続いて、その実効誘電
率εreffと、空気中での単位長当たりの容量Ｃ0 とか
ら、平行線ケーブルの単位長当たりの容量を求めること
になる。
に、公式に従って、空気中での単位長当たりの容量Ｃ0
と、空気中での特性インピーダンスＺ0 を算出し、次
に、この特性インピーダンスＺ0 と、ケーブルライブラ
リ２００に登録される特性インピーダンスＺとから誘電
体の実効誘電率εreffを算出し、続いて、その実効誘電
率εreffと、空気中での単位長当たりの容量Ｃ0 とか
ら、平行線ケーブルの単位長当たりの容量を求めること
になる。
【００８６】なお、導線の形状が円形でなくて任意の形
状を持つようなときには、このような簡単な式で表せな
いことになるが、そのような場合にも、２次元解析の計
算手法を用いることで、空気中での単位長当たりのイン
ダクタンスＬ0 と、空気中での単位長当たりの容量Ｃ0
とを求めることができるので、同様の数式処理に従っ
て、誘電体で被覆されたケーブルの単位長さ当たりの容
量Ｃを求めることができるものである。
状を持つようなときには、このような簡単な式で表せな
いことになるが、そのような場合にも、２次元解析の計
算手法を用いることで、空気中での単位長当たりのイン
ダクタンスＬ0 と、空気中での単位長当たりの容量Ｃ0
とを求めることができるので、同様の数式処理に従っ
て、誘電体で被覆されたケーブルの単位長さ当たりの容
量Ｃを求めることができるものである。
【００８７】ステップ７で、平行線ケーブルの単位長当
たりの容量を求めると、続いて、ステップ８で、非平行
リンクに挿入するコンデンサの容量を決定する。このコ
ンデンサの容量の決定処理は、その求めた平行線ケーブ
ルの単位長当たりの容量と、非平行リンクのリンク長と
を乗算し、その乗算値を挿入するコンデンサの個数で割
り算することで求める。図１８の例で示すならば、その
乗算値を“４”で割り算することで求めるのである。
たりの容量を求めると、続いて、ステップ８で、非平行
リンクに挿入するコンデンサの容量を決定する。このコ
ンデンサの容量の決定処理は、その求めた平行線ケーブ
ルの単位長当たりの容量と、非平行リンクのリンク長と
を乗算し、その乗算値を挿入するコンデンサの個数で割
り算することで求める。図１８の例で示すならば、その
乗算値を“４”で割り算することで求めるのである。
【００８８】このようにして、図８の処理フローに従っ
て、図１８に示すような金属線分（太さは平行線ケーブ
ルの芯線の直径）とコンデンサとからなるはしご型回路
形態のモーメント法モデルが生成されることになる。
て、図１８に示すような金属線分（太さは平行線ケーブ
ルの芯線の直径）とコンデンサとからなるはしご型回路
形態のモーメント法モデルが生成されることになる。
【００８９】次に、図９の処理フロー（図７の処理フロ
ーのステップ１２で実行する撚り線ケーブルのモデル化
アルゴリズムの処理フロー）について説明する。すなわ
ち、図７の処理フローのステップ１２に従って、撚り線
ケーブルの持つ非平行リンクのモデル化処理に入ると、
図９の処理フローに示すように、先ず最初に、ステップ
１で、非平行リンクのリンク長と、撚り線ケーブルの撚
りピッチとから、非平行リンクを構成する撚り線がπ×
ｎ回転となるようにと撚りピッチを最も近いものに変更
する。例えば、撚り線ケーブルが１０ｍｍで１回転を示
すという撚りピッチを持つときに、非平行リンクのリン
ク長が４７ｍｍであるときには、本来は４.7回転するの
であるが、これをπ×ｎ回転の形態の最も近い４.5回転
になるようにと撚りピッチを変更するのである。
ーのステップ１２で実行する撚り線ケーブルのモデル化
アルゴリズムの処理フロー）について説明する。すなわ
ち、図７の処理フローのステップ１２に従って、撚り線
ケーブルの持つ非平行リンクのモデル化処理に入ると、
図９の処理フローに示すように、先ず最初に、ステップ
１で、非平行リンクのリンク長と、撚り線ケーブルの撚
りピッチとから、非平行リンクを構成する撚り線がπ×
ｎ回転となるようにと撚りピッチを最も近いものに変更
する。例えば、撚り線ケーブルが１０ｍｍで１回転を示
すという撚りピッチを持つときに、非平行リンクのリン
ク長が４７ｍｍであるときには、本来は４.7回転するの
であるが、これをπ×ｎ回転の形態の最も近い４.5回転
になるようにと撚りピッチを変更するのである。
【００９０】具体的には、非平行リンクのリンク長Ｘを
撚りピッチψで割り算することで、商ｒと剰余Ｒを求め
て、０≦Ｒ≦０.5であるならば、ｒ回転となるようにと
撚りピッチψを“Ｘ／ｒ”に変更し、０.5＜Ｒ＜１であ
るならば、（ｒ＋０.5）回転となるようにと撚りピッチ
ψを“Ｘ／（ｒ＋０.5）”に変更することで行う。
撚りピッチψで割り算することで、商ｒと剰余Ｒを求め
て、０≦Ｒ≦０.5であるならば、ｒ回転となるようにと
撚りピッチψを“Ｘ／ｒ”に変更し、０.5＜Ｒ＜１であ
るならば、（ｒ＋０.5）回転となるようにと撚りピッチ
ψを“Ｘ／（ｒ＋０.5）”に変更することで行う。
【００９１】但し、このような一定のアルゴリズムに従
って撚りピッチψを変更すると、ケーブル全体で、撚り
ピッチψが小さくなる方に偏るとか、大きくなる方に偏
るとかいうことが起こる可能性がある。そこで、例え
ば、前回の非平行リンクで撚りピッチψを小さくすると
きには、今回の非平行リンクでは撚りピッチψを大きく
するというようにすることで、ケーブル全体でバランス
させていくような撚りピッチψの変更方法を採ってもよ
い。
って撚りピッチψを変更すると、ケーブル全体で、撚り
ピッチψが小さくなる方に偏るとか、大きくなる方に偏
るとかいうことが起こる可能性がある。そこで、例え
ば、前回の非平行リンクで撚りピッチψを小さくすると
きには、今回の非平行リンクでは撚りピッチψを大きく
するというようにすることで、ケーブル全体でバランス
させていくような撚りピッチψの変更方法を採ってもよ
い。
【００９２】この変更処理は、図１９に示すように、非
平行リンクの端点に、撚り線の腹位置がくるようにと撚
りピッチを補正することを意味する。このようにする
と、撚り線ケーブルの行きの長さと帰りの長さとが、同
一となるように実現できるからである。
平行リンクの端点に、撚り線の腹位置がくるようにと撚
りピッチを補正することを意味する。このようにする
と、撚り線ケーブルの行きの長さと帰りの長さとが、同
一となるように実現できるからである。
【００９３】ステップ１で撚りピッチを変更すると、続
いて、ステップ２で、非平行リンクの両端の法線ベクト
ルのなす捩れ角度θを求める。すなわち、図１６に示す
ように、非平行リンクの始点（Ａ，Ｂ）での法線ベクト
ルと、終点（Ｃ，Ｄ）での法線ベクトルとのなす捩れ角
度θを求めるのである。
いて、ステップ２で、非平行リンクの両端の法線ベクト
ルのなす捩れ角度θを求める。すなわち、図１６に示す
ように、非平行リンクの始点（Ａ，Ｂ）での法線ベクト
ルと、終点（Ｃ，Ｄ）での法線ベクトルとのなす捩れ角
度θを求めるのである。
【００９４】続いて、ステップ３で、図７の処理フロー
のステップ２で設定した区分数ｍで規定される回転角度
の補正処理に入って、その区分数ｍと、ステップ１で変
更した非平行リンクでの撚り回転数と、ステップ２で求
めた捩れ角度θとから、区分数ｍで規定される回転角度
を補正する。
のステップ２で設定した区分数ｍで規定される回転角度
の補正処理に入って、その区分数ｍと、ステップ１で変
更した非平行リンクでの撚り回転数と、ステップ２で求
めた捩れ角度θとから、区分数ｍで規定される回転角度
を補正する。
【００９５】この補正は、捩れ角度θを吸収するために
行うものである。例えば、偶数で設定される区分数ｍが
“１２”であるときには、「３６０度÷１２個＝３０
度」で規定される３０度ずつ回転したサンプリング位置
で、撚り線ケーブルの芯線位置を求めることを意味す
る。しかし、このまま３０度ずつ回転すると、非平行リ
ンクの両端の法線ベクトルのなす捩れ角度θが吸収でき
ない。
行うものである。例えば、偶数で設定される区分数ｍが
“１２”であるときには、「３６０度÷１２個＝３０
度」で規定される３０度ずつ回転したサンプリング位置
で、撚り線ケーブルの芯線位置を求めることを意味す
る。しかし、このまま３０度ずつ回転すると、非平行リ
ンクの両端の法線ベクトルのなす捩れ角度θが吸収でき
ない。
【００９６】そこで、例えば、ステップ１で変更した撚
り回転数が４.5回転である場合に、非平行リンクでの総
サンプリング数が「４.5回転×１２個＝５４個」となる
ことを考慮して、例えば、捩れ角度が４０度であるとき
には、この５４個で４０度を吸収すべく、３０度に「４
０度÷５４個＝０.7４度」を加算したり、３０度から
「４０度÷５４個＝０.7４度」を減算するようと、区分
数ｍで規定される回転角度を補正するのである。
り回転数が４.5回転である場合に、非平行リンクでの総
サンプリング数が「４.5回転×１２個＝５４個」となる
ことを考慮して、例えば、捩れ角度が４０度であるとき
には、この５４個で４０度を吸収すべく、３０度に「４
０度÷５４個＝０.7４度」を加算したり、３０度から
「４０度÷５４個＝０.7４度」を減算するようと、区分
数ｍで規定される回転角度を補正するのである。
【００９７】このとき、加算するか減算するかは、撚り
線の回転方向と、捩れ角度θの回転方向とで決定される
ことになる。具体的に説明するならば、撚り線が進行方
向に対して時計回りに回転するときにあって、捩れ角度
θが時計方向の回転角度を示すときには加算処理を行
い、捩れ角度θが反時計方向の回転角度を示すときには
減算処理を行う。また、撚り線が進行方向に対して反時
計回りに回転するときにあって、捩れ角度θが時計方向
の回転角度を示すときには減算処理を行い、捩れ角度θ
が反時計方向の回転角度を示すときには加算処理を行う
ことになる。
線の回転方向と、捩れ角度θの回転方向とで決定される
ことになる。具体的に説明するならば、撚り線が進行方
向に対して時計回りに回転するときにあって、捩れ角度
θが時計方向の回転角度を示すときには加算処理を行
い、捩れ角度θが反時計方向の回転角度を示すときには
減算処理を行う。また、撚り線が進行方向に対して反時
計回りに回転するときにあって、捩れ角度θが時計方向
の回転角度を示すときには減算処理を行い、捩れ角度θ
が反時計方向の回転角度を示すときには加算処理を行う
ことになる。
【００９８】ステップ３で区分数ｍで規定される回転角
度を補正すると、続いて、ステップ４で、ステップ１で
変更した撚りピッチと、ステップ３で補正した回転角度
とから、撚り線ケーブルの芯線の区分位置を決定すると
ともに、コンデンサの挿入位置を決定する。
度を補正すると、続いて、ステップ４で、ステップ１で
変更した撚りピッチと、ステップ３で補正した回転角度
とから、撚り線ケーブルの芯線の区分位置を決定すると
ともに、コンデンサの挿入位置を決定する。
【００９９】この芯線の区分位置の決定処理は、ステッ
プ３で補正した回転角度をαで表すならば、始点から終
点に向けて、最初はα／２度、その後はα度ずつ回転さ
せていって、その回転角度と撚りピッチとで決まる位置
の座標を特定することで行う。また、また、コンデンサ
の挿入位置の決定処理は、図２０に示すように、最初は
α／２度、その後はα度ずつ回転する所を挿入位置とし
て特定することで行う。
プ３で補正した回転角度をαで表すならば、始点から終
点に向けて、最初はα／２度、その後はα度ずつ回転さ
せていって、その回転角度と撚りピッチとで決まる位置
の座標を特定することで行う。また、また、コンデンサ
の挿入位置の決定処理は、図２０に示すように、最初は
α／２度、その後はα度ずつ回転する所を挿入位置とし
て特定することで行う。
【０１００】続いて、ステップ５で、撚り線ケーブルの
単位長当たりの容量を求める。この容量は、上述した平
行線ケーブルと同じ手法に従い、先ず最初に、公式に従
って、空気中での単位長当たりの容量Ｃ0 と、空気中で
の特性インピーダンスＺ0 を算出し、次に、この特性イ
ンピーダンスＺ0 と、ケーブルライブラリ２００に登録
される特性インピーダンスＺとから誘電体の実効誘電率
εreffを算出し、続いて、その実効誘電率εreffと、空
気中での単位長当たりの容量Ｃ0 とから求めることで行
うことになる。
単位長当たりの容量を求める。この容量は、上述した平
行線ケーブルと同じ手法に従い、先ず最初に、公式に従
って、空気中での単位長当たりの容量Ｃ0 と、空気中で
の特性インピーダンスＺ0 を算出し、次に、この特性イ
ンピーダンスＺ0 と、ケーブルライブラリ２００に登録
される特性インピーダンスＺとから誘電体の実効誘電率
εreffを算出し、続いて、その実効誘電率εreffと、空
気中での単位長当たりの容量Ｃ0 とから求めることで行
うことになる。
【０１０１】続いて、ステップ６で、非平行リンクに挿
入するコンデンサの容量を決定する。このコンデンサの
容量の決定処理は、ステップ５で求めた撚り線ケーブル
の単位長当たりの容量と、非平行リンクのリンク長とを
乗算し、その乗算値を挿入するコンデンサの個数で割り
算することで求める。図２０の例で示すならば、その乗
算値を“１２”で割り算することで求めるのである。
入するコンデンサの容量を決定する。このコンデンサの
容量の決定処理は、ステップ５で求めた撚り線ケーブル
の単位長当たりの容量と、非平行リンクのリンク長とを
乗算し、その乗算値を挿入するコンデンサの個数で割り
算することで求める。図２０の例で示すならば、その乗
算値を“１２”で割り算することで求めるのである。
【０１０２】このようにして、図９の処理フローに従っ
て、図２０に示すような金属線分（太さは撚り線ケーブ
ルの芯線の直径）とコンデンサとからなるはしご型回路
形態のモーメント法モデルが生成されることになる。
て、図２０に示すような金属線分（太さは撚り線ケーブ
ルの芯線の直径）とコンデンサとからなるはしご型回路
形態のモーメント法モデルが生成されることになる。
【０１０３】すなわち、図８及び図９の処理フローに従
って、平行線ケーブルや撚り線ケーブルは、図２１に示
すようなモーメント法モデルに変換されることになる。
このモーメント法モデルに対してモーメント法を適用す
る場合には、図２２（ａ）に示すような電流の流れを想
定すべく、図２２（ｂ）に示すような区分電流分布を想
定して、それを使って、相互インピーダンスを求めてい
くことになる。このとき、コンデンサの容量は自己イン
ピーダンス成分に加わることになる。
って、平行線ケーブルや撚り線ケーブルは、図２１に示
すようなモーメント法モデルに変換されることになる。
このモーメント法モデルに対してモーメント法を適用す
る場合には、図２２（ａ）に示すような電流の流れを想
定すべく、図２２（ｂ）に示すような区分電流分布を想
定して、それを使って、相互インピーダンスを求めてい
くことになる。このとき、コンデンサの容量は自己イン
ピーダンス成分に加わることになる。
【０１０４】このモーメント法モデルの精度を更に高め
る必要があるときには、図２３に示すように、芯線部分
から切り出される金属の線分部分にインダクタンスを持
たせることになる。
る必要があるときには、図２３に示すように、芯線部分
から切り出される金属の線分部分にインダクタンスを持
たせることになる。
【０１０５】最後に、本発明者が提案した平行線ケーブ
ル及び撚り線ケーブルのモデルの有効性を検証するため
に行った検証実験の結果について説明する。この検証実
験は、図２４（ａ）に示すように、プリント板に平行線
ケーブル（撚り線ケーブル）を接続し、プリント板から
この平行線ケーブル（撚り線ケーブル）に高周波信号を
供給するときに、この平行線ケーブル（撚り線ケーブ
ル）が放射する電磁界強度を上述のモデルを使ってシミ
ュレーションして、そのシミュレーション値を実測値と
比較することで行った。
ル及び撚り線ケーブルのモデルの有効性を検証するため
に行った検証実験の結果について説明する。この検証実
験は、図２４（ａ）に示すように、プリント板に平行線
ケーブル（撚り線ケーブル）を接続し、プリント板から
この平行線ケーブル（撚り線ケーブル）に高周波信号を
供給するときに、この平行線ケーブル（撚り線ケーブ
ル）が放射する電磁界強度を上述のモデルを使ってシミ
ュレーションして、そのシミュレーション値を実測値と
比較することで行った。
【０１０６】平行線ケーブルとしては、図２４（ｂ）に
示すように、９４５ｍｍの長さを持つものを使用し、先
端に１８０Ωを接続した。また、撚り線ケーブルとして
は、図２４（ｃ）に示すように、３００ｍｍの長さを持
つものを使用し、先端に１２０Ωを接続した。導体外径
等のデータは図２５に示してある。
示すように、９４５ｍｍの長さを持つものを使用し、先
端に１８０Ωを接続した。また、撚り線ケーブルとして
は、図２４（ｃ）に示すように、３００ｍｍの長さを持
つものを使用し、先端に１２０Ωを接続した。導体外径
等のデータは図２５に示してある。
【０１０７】図２６及び図２７に、検証実験の結果を図
示する。図２６は、平行線ケーブルについての実験結果
であり、図２７は、撚り線ケーブルについての実験結果
である。
示する。図２６は、平行線ケーブルについての実験結果
であり、図２７は、撚り線ケーブルについての実験結果
である。
【０１０８】平行線ケーブルについては、観測点を１０
ｍ先に設定する１０ｍ法に従った。周波数は５０ＭＨｚ
とした。図２６の上段の図は垂直方向の実験結果、下段
の図は水平方向の実験結果であり、図中に示す実線はシ
ミュレーション結果、点は実測値である。
ｍ先に設定する１０ｍ法に従った。周波数は５０ＭＨｚ
とした。図２６の上段の図は垂直方向の実験結果、下段
の図は水平方向の実験結果であり、図中に示す実線はシ
ミュレーション結果、点は実測値である。
【０１０９】また、撚り線ケーブルについては、観測点
を１０ｍ先に設定する３ｍ法に従った。周波数は５０Ｍ
Ｈｚとした。図２６の上段の図は垂直方向の実験結果、
下段の図は水平方向の実験結果であり、図中に示す実線
はシミュレーション結果、点は実測値である。
を１０ｍ先に設定する３ｍ法に従った。周波数は５０Ｍ
Ｈｚとした。図２６の上段の図は垂直方向の実験結果、
下段の図は水平方向の実験結果であり、図中に示す実線
はシミュレーション結果、点は実測値である。
【０１１０】この図２６及び図２７から、シミュレーシ
ョン結果と実測値とがよい一致することが分かる。これ
から、本発明の有効性が検証された。この実施例では、
撚り線ケーブルをモデル化するときに、ユーザにより指
定されたケーブル経由点で、撚り線の腹位置が配置され
ることになるようにと、撚りピッチを変更する構成を採
ったが、撚り線ピッチを変更するのではなくて、ユーザ
により指定されたケーブル経由点を変更する（方向ベク
トルの方向に長くしたり、短くすることで変更する）こ
とで、ケーブル経由点で、撚り線の腹位置が配置される
ことになるようにする構成を採ることも可能である。こ
の構成を採るときに、自動的にケーブル経由点を変更す
るのではなくて、ユーザにその旨を通知して、了解が得
られることを条件にして変更するようにしてもよい。
ョン結果と実測値とがよい一致することが分かる。これ
から、本発明の有効性が検証された。この実施例では、
撚り線ケーブルをモデル化するときに、ユーザにより指
定されたケーブル経由点で、撚り線の腹位置が配置され
ることになるようにと、撚りピッチを変更する構成を採
ったが、撚り線ピッチを変更するのではなくて、ユーザ
により指定されたケーブル経由点を変更する（方向ベク
トルの方向に長くしたり、短くすることで変更する）こ
とで、ケーブル経由点で、撚り線の腹位置が配置される
ことになるようにする構成を採ることも可能である。こ
の構成を採るときに、自動的にケーブル経由点を変更す
るのではなくて、ユーザにその旨を通知して、了解が得
られることを条件にして変更するようにしてもよい。
【０１１１】
【発明の効果】以上説明したように、本発明の電磁界強
度算出装置によれば、モーメント法を使って、電気回路
装置の放射する電磁界強度を算出する構成を採るときに
あって、電気回路装置に実装される２線ケーブルの配設
位置を自由に変更できるようにする構成を採るので、実
際のものに合ったシミュレーションが実行できるように
なるとともに、ケーブルの配置を色々と変えてみて、そ
の効果を見るというようなことが簡単に実行できるよう
になる。
度算出装置によれば、モーメント法を使って、電気回路
装置の放射する電磁界強度を算出する構成を採るときに
あって、電気回路装置に実装される２線ケーブルの配設
位置を自由に変更できるようにする構成を採るので、実
際のものに合ったシミュレーションが実行できるように
なるとともに、ケーブルの配置を色々と変えてみて、そ
の効果を見るというようなことが簡単に実行できるよう
になる。
【０１１２】そして、平行線ケーブルや撚り線ケーブル
を、金属線分とコンデンサとからなるはしご型回路で高
精度にモデル化するので、平行線ケーブルや撚り線ケー
ブルを実装する電気回路装置から放射される電磁界強度
を高精度に算出できるようになる。
を、金属線分とコンデンサとからなるはしご型回路で高
精度にモデル化するので、平行線ケーブルや撚り線ケー
ブルを実装する電気回路装置から放射される電磁界強度
を高精度に算出できるようになる。
【図１】本発明の原理構成図である。
【図２】本発明の実行する処理フローの一実施例であ
る。
る。
【図３】ケーブルライブラリのデータ構造図である。
【図４】本発明の実行する処理フローの一実施例であ
る。
る。
【図５】本発明の実行する処理フローの一実施例であ
る。
る。
【図６】本発明の実行する処理フローの一実施例であ
る。
る。
【図７】本発明の実行する処理フローの一実施例であ
る。
る。
【図８】本発明の実行する処理フローの一実施例であ
る。
る。
【図９】本発明の実行する処理フローの一実施例であ
る。
る。
【図１０】ケーブル配置の説明図である。
【図１１】ディスプレイ画面の一実施例である。
【図１２】ディスプレイ画面の一実施例である。
【図１３】リンク情報の説明図である。
【図１４】リンク情報の説明図である。
【図１５】法線ベクトルの説明図である。
【図１６】法線ベクトルの捩れ角度の説明図である。
【図１７】始点と終点の対応関係の説明図である。
【図１８】コンデンサ挿入位置の説明図である。
【図１９】撚りピッチ変更処理の説明図である。
【図２０】コンデンサ挿入位置の説明図である。
【図２１】モーメント法モデルの説明図である。
【図２２】区分電流の説明図である。
【図２３】モーメント法モデルの説明図である。
【図２４】検証実験の説明図である。
【図２５】検証実験の説明図である。
【図２６】検証実験の説明図である。
【図２７】検証実験の説明図である。
１ 電磁界強度算出装置
１０ 入力手段
１１ ライブラリ手段
１２ 設定手段
１３ モデル生成手段
１４ 生成手段
１５ 算出手段
─────────────────────────────────────────────────────
フロントページの続き
(56)参考文献 特開 平７−302278（ＪＰ，Ａ)
(58)調査した分野(Int.Cl.7，ＤＢ名)
G01R 29/08
G06F 17/00 - 17/50
Claims (8)
- 【請求項１】 モーメント法を使って、電気回路装置の
放射する電磁界強度を算出する電磁界強度算出装置にお
いて、 解析対象となる電気回路装置の構造情報を入力する入力
手段と、 上記構造情報がケーブル接続箇所を持たないときには、
対話処理に従い、ケーブル接続箇所と、該ケーブル接続
箇所を接続する１つ又は複数のケーブル経由位置とを設
定し、一方、ケーブル接続箇所を持つときには、対話処
理に従い、該ケーブル接続箇所を接続する１つ又は複数
のケーブル経由位置を設定し、さらに、対話処理に従
い、使用する撚り線ケーブルの識別情報を設定する設定
手段と、ケーブル方向を変化させるケーブル経由位置に撚り線の
腹位置が配置されることになるようにと、上記撚り線ケ
ーブルの撚りピッチを変更するとともに、その変更した
撚りピッチと規定の回転角度との指定するケーブル芯線
位置に従って、ケーブルの方向が変化しないケーブル開
始点とケーブル終了点との間の線部材を生成すること
で、上 記ケーブル接続箇所を接続する上記撚り線ケーブ
ルのモーメント法モデルを生成する生成手段と、 上記モーメント法モデルに基づき、モーメント法を用い
て電気回路装置の放射する電磁界強度を算出する算出手
段とを備えることを、 特徴とする電磁界強度算出装置。 - 【請求項２】 請求項１記載の電磁界強度算出装置にお
いて、 上記生成手段は、撚り線ケーブルの撚りピッチの変更に
あたって、撚りピッチを短くする変更と撚りピッチを長
くする変更とがバランスするように、撚りピッチを変更
することを、 特徴とする電磁界強度算出装置。 - 【請求項３】 モーメント法を使って、電気回路装置の
放射する電磁界強度を算出する電磁界強度算出装置にお
いて、 解析対象となる電気回路装置の構造情報を入力する入力
手段と、 上記構造情報がケーブル接続箇所を持たないときには、
対話処理に従い、ケーブル接続箇所と、該ケーブル接続
箇所を接続する１つ又は複数のケーブル経由位置とを設
定し、一方、ケーブル接続箇所を持つときには、対話処
理に従い、該ケーブル接続箇所を接続する１つ又は複数
のケーブル経由位置を設定し、さらに、対話処理に従
い、使用する撚り線ケーブルの識別情報を設定し、そし
て、ケーブル方向を変化させるケーブル経由位置に撚り
線の腹位置が配置されることになるようにと、設定した
ケーブル経由位置をその近傍で移動させることでケーブ
ル経由位置を最終的に設定する設定手段と、上記撚り線ケーブルの撚りピッチと規定の回転角度との
指定するケーブル芯線位置に従って、ケーブルの方向が
変化しないケーブル開始点とケーブル終了点との間の線
部材を生成することで、上 記ケーブル接続箇所を接続す
る上記撚り線ケーブルのモーメント法モデルを生成する
生成手段と、 上記モーメント法モデルに基づき、モーメント法を用い
て電気回路装置の放射する電磁界強度を算出する算出手
段とを備えることを、 特徴とする電磁界強度算出装置。 - 【請求項４】 請求項１ないし３のいずれか１項に記載
の電磁界強度算出装置において、 上記生成手段は、ケーブル開始点とケーブル終了点との
間に捩れが存在するときには、その捩れを吸収すること
になるようにと回転角度を補正して、その補正した回転
角度を使って線部材を生成することを、 特徴とする電磁界強度算出装置。 - 【請求項５】 モーメント法を使って、電気回路装置の
放射する電磁界強度を算出する電磁界強度算出装置にお
いて、 解析対象となる電気回路装置の構造情報を入力する入力
手段と、 上記構造情報がケーブル接続箇所を持たないときには、
対話処理に従い、ケーブル接続箇所と、該ケーブル接続
箇所を接続する１つ又は複数のケーブル経由位置とを設
定し、一方、ケーブル接続箇所を持つときには、対話処
理に従い、該ケーブル接続箇所を接続する１つ又は複数
のケーブル経由位置を設定し、さらに、対話処理に従
い、使用する平行線ケーブルの識別情報を設定する設定
手段と、ケーブル方向を変化させないケーブル開始点とケーブル
終了点との間を等間隔に区切る規定長以下の長さの指定
するケーブル芯線位置に従って、該ケーブル開始点と該
ケーブル終了点との間の線部材を生成し、さらに、この
とき、該ケーブ ル開始点と該ケーブル終了点との間に捩
れが存在するときには、その捩れを考慮したケーブル芯
線位置に従って線部材を生成することで、上 記ケーブル
接続箇所を接続する上記平行線ケーブルのモーメント法
モデルを生成する生成手段と、 上記モーメント法モデルに基づき、モーメント法を用い
て電気回路装置の放射する電磁界強度を算出する算出手
段とを備えることを、 特徴とする電磁界強度算出装置。 - 【請求項６】 モーメント法を使って、電気回路装置の
放射する電磁界強度を算出する電磁界強度算出装置にお
いて、 解析対象となる電気回路装置の構造情報を入力する入力
手段と、 上記構造情報がケーブル接続箇所を持たないときには、
対話処理に従い、ケーブル接続箇所と、該ケーブル接続
箇所を接続する１つ又は複数のケーブル経由位置とを設
定し、一方、ケーブル接続箇所を持つときには、対話処
理に従い、該ケーブル接続箇所を接続する１つ又は複数
のケーブル経由位置を設定し、さらに、対話処理に従
い、使用する平行線ケーブルの識別情報を設定する設定
手段と、ケーブル方向を変化させないケーブル開始点とケーブル
終了点との間を等間隔に区切る規定長以下の長さの指定
するケーブル芯線位置に従って、該ケーブル開始点と該
ケーブル終了点との間の線部材を生成し、さらに、この
とき、該ケーブル開始点と該ケーブル終了点との間に捩
れが存在するときには、上記平行線ケーブルの捩れが集
積しないようにと、該ケーブル開始点と該ケーブル終了
点との間の位置関係を決定することで、上 記ケーブル接
続箇所を接続する上記平行線ケーブルのモーメント法モ
デルを生成する生成手段と、 上記モーメント法モデルに基づき、モーメント法を用い
て電気回路装置の放射する電磁界強度を算出する算出手
段とを備えることを、 特徴とする電磁界強度算出装置。 - 【請求項７】 モーメント法を使って、電気回路装置の
放射する電磁界強度を算出する電磁界強度算出方法にお
いて、 解析対象となる電気回路装置の構造情報を入力する過程
と、 上記構造情報がケーブル接続箇所を持たないときには、
対話処理に従い、ケーブル接続箇所と、該ケーブル接続
箇所を接続する１つ又は複数のケーブル経由位置とを設
定し、一方、ケーブル接続箇所を持つときには、対話処
理に従い、該ケーブル接続箇所を接続する１つ又は複数
のケーブル経由位置を設定し、さらに、対話処理に従
い、使用する撚り線ケーブルの識別情報を設定する過程
と、ケーブル方向を変化させるケーブル経由位置に撚り線の
腹位置が配置されることになるようにと、上記撚り線ケ
ーブルの撚りピッチを変更するとともに、その変更した
撚りピッチと規定の回転角度との指定するケーブル芯線
位置に従って、ケーブルの方向が変化しないケーブル開
始点とケーブル終了点との間の線部材を生成すること
で、上 記ケーブル接続箇所を接続する上記撚り線ケーブ
ルのモーメント法モデルを生成する過程と、 上記モーメント法モデルに基づき、モーメント法を用い
て電気回路装置の放射する電磁界強度を算出する過程と
を備えることを、 特徴とする電磁界強度算出方法。 - 【請求項８】 モーメント法を使って、電気回路装置の
放射する電磁界強度を算出する電磁界強度算出装置の実
現に用いられるプログラムが記憶されるプログラム記憶
媒体であって、 解析対象となる電気回路装置の構造情報を入力する処理
と、 上記構造情報がケーブル接続箇所を持たないときには、
対話処理に従い、ケーブル接続箇所と、該ケーブル接続
箇所を接続する１つ又は複数のケーブル経由位置とを設
定し、一方、ケーブル接続箇所を持つときには、対話処
理に従い、該ケーブル接続箇所を接続する１つ又は複数
のケーブル経由位置を設定し、さらに、対話処理に従
い、使用する撚り線ケーブルの識別情報を設定する処理
と、ケーブル方向を変化させるケーブル経由位置に撚り線の
腹位置が配置されることになるようにと、上記撚り線ケ
ーブルの撚りピッチを変更するとともに、その変更した
撚りピッチと規定の回転角度との指定するケーブル芯線
位置に従って、ケーブルの方向が変化しないケーブル開
始点とケーブル終了点との間の線部材を生成すること
で、上 記ケーブル接続箇所を接続する上記撚り線ケーブ
ルのモーメント法モデルを生成する処理と、 上記モーメント法モデルに基づき、モーメント法を用い
て電気回路装置の放射する電磁界強度を算出する処理と
をコンピュータに実行させるプログラムが記憶されるこ
とを、 特徴とするプログラム記憶媒体。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP21758397A JP3405905B2 (ja) | 1997-08-12 | 1997-08-12 | 電磁界強度算出装置及び方法並びにプログラム記憶媒体 |
| US09/045,125 US6285957B1 (en) | 1997-08-12 | 1998-03-20 | Apparatus and method for calculating the strength of an electromagnetic field radiated from an electric device |
| DE19816097A DE19816097B4 (de) | 1997-08-12 | 1998-04-09 | Verfahren, Vorrichtung und Programmspeichermedium zum Berechnen der Stärke eines elektromagnetischen Feldes |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP21758397A JP3405905B2 (ja) | 1997-08-12 | 1997-08-12 | 電磁界強度算出装置及び方法並びにプログラム記憶媒体 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH1164412A JPH1164412A (ja) | 1999-03-05 |
| JP3405905B2 true JP3405905B2 (ja) | 2003-05-12 |
Family
ID=16706565
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP21758397A Expired - Fee Related JP3405905B2 (ja) | 1997-08-12 | 1997-08-12 | 電磁界強度算出装置及び方法並びにプログラム記憶媒体 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US6285957B1 (ja) |
| JP (1) | JP3405905B2 (ja) |
| DE (1) | DE19816097B4 (ja) |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3488629B2 (ja) * | 1998-04-07 | 2004-01-19 | 富士通株式会社 | 放射電磁界耐性算出装置及び方法並びにプログラム記録媒体 |
| JP2000214191A (ja) * | 1999-01-27 | 2000-08-04 | Fujitsu Ltd | シミュレ―ション装置及び方法並びにプログラム記録媒体 |
| JP3980811B2 (ja) * | 1999-10-20 | 2007-09-26 | 富士通株式会社 | 電磁界強度算出装置、電磁界強度算出方法およびその方法をコンピュータに実行させるプログラムを記録したコンピュータ読み取り可能な記録媒体 |
| US6711511B2 (en) | 2000-07-06 | 2004-03-23 | Fujitsu Limited | Electromagnetic wave analyzing apparatus and computer readable storage medium |
| US6625554B2 (en) * | 2001-06-22 | 2003-09-23 | Hitachi, Ltd. | Method and apparatus for determining a magnetic field |
| 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 |
| US7260797B2 (en) * | 2004-09-07 | 2007-08-21 | Synopsys, Inc. | Method and apparatus for estimating parasitic capacitance |
| US8352887B2 (en) | 2010-12-03 | 2013-01-08 | Synopsys, Inc. | High performance design rule checking technique |
| US8453103B2 (en) | 2010-12-03 | 2013-05-28 | Synopsys, Inc. | Real time DRC assistance for manual layout editing |
| US8677297B2 (en) | 2010-12-03 | 2014-03-18 | Synopsys, Inc. | Low-overhead multi-patterning design rule check |
| US8448097B2 (en) | 2010-12-03 | 2013-05-21 | Synopsys, Inc. | High performance DRC checking algorithm for derived layer based rules |
| DE102014005242B3 (de) * | 2014-04-08 | 2015-07-09 | SLE quality engineering GmbH & Co. KG | Verfahren und Vorrichtung zum Bestimmen einer Winkellage von Einzelleitungen an einer vorbestimmten Querschnittsstelle in einer mehradrigen Mantelleitung |
| CN106019987A (zh) * | 2016-07-19 | 2016-10-12 | 四川九洲空管科技有限责任公司 | 一种三维交互式模拟航迹生成方法及系统 |
| JP6906718B1 (ja) * | 2020-04-03 | 2021-07-21 | 三菱電機株式会社 | 伝送路長検出装置及びネットワーク通信装置 |
| CN112100809B (zh) * | 2020-08-10 | 2022-06-14 | 华南理工大学 | 一种基于多物理场耦合仿真三相同轴超导电缆设计的方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3091815B2 (ja) * | 1994-02-25 | 2000-09-25 | 富士通株式会社 | 電磁界強度算出装置 |
| DE19616772A1 (de) * | 1995-04-28 | 1996-10-31 | Fujitsu Ltd | Berechnungsgerät für eine elektromagnetische Feldintensität |
| JP3428232B2 (ja) * | 1995-06-16 | 2003-07-22 | 富士通株式会社 | 電磁界強度算出装置 |
| EP0778533B1 (en) * | 1995-11-16 | 2003-10-22 | Fujitsu Limited | Electromagnetic field intensity calculation apparatus and method |
| US5903477A (en) * | 1996-04-10 | 1999-05-11 | Fujitsu Limited | Simulation apparatus and simulation method for electromagnetic field intensity using moment method |
| JPH1115814A (ja) * | 1997-06-26 | 1999-01-22 | Fujitsu Ltd | モーメント法を用いたシミュレーション装置及び方法並びにプログラム記憶媒体 |
-
1997
- 1997-08-12 JP JP21758397A patent/JP3405905B2/ja not_active Expired - Fee Related
-
1998
- 1998-03-20 US US09/045,125 patent/US6285957B1/en not_active Expired - Lifetime
- 1998-04-09 DE DE19816097A patent/DE19816097B4/de not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| DE19816097B4 (de) | 2005-06-02 |
| DE19816097A1 (de) | 1999-02-18 |
| JPH1164412A (ja) | 1999-03-05 |
| US6285957B1 (en) | 2001-09-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3405905B2 (ja) | 電磁界強度算出装置及び方法並びにプログラム記憶媒体 | |
| Theodoulidis et al. | Eddy current coil interaction with a right-angled conductive wedge | |
| US10338117B2 (en) | Induced field determination using diffuse field reciprocity | |
| US6182023B1 (en) | Electromagnetic field intensity computing device | |
| CN108960206A (zh) | 视频帧处理方法和装置 | |
| US6879942B1 (en) | Apparatus for calculating immunity from radiated electromagnetic field, method for achieving calculation, and storage medium storing programs therefor | |
| CN110377994A (zh) | 一种基于屏蔽效能的编织屏蔽层半径设计方法 | |
| US20160170047A1 (en) | Method for predicting electromagnetic radiation characteristics, computer-readable recording medium and simulator | |
| Wang et al. | An integrated approach for electromagnetic compatible commercial aircraft engine cable harnessing | |
| JP3980811B2 (ja) | 電磁界強度算出装置、電磁界強度算出方法およびその方法をコンピュータに実行させるプログラムを記録したコンピュータ読み取り可能な記録媒体 | |
| Hampton et al. | Closed-form expressions for the magnetic fields of rectangular and circular finite-length solenoids and current loops | |
| CN111737885A (zh) | 复杂线束电磁耦合效应分析方法和装置 | |
| EP1846833A1 (en) | Apparatus and method for determining intersections | |
| Bai et al. | Comparison of Calculation Methods of Braided Shield Cable Transfer Impedance Using FSV Method. | |
| JP3703812B2 (ja) | Ｆｄｔｄ法を用いた電磁界解析方法、電磁界解析における媒質表現方法、シミュレーション装置、及びプログラム | |
| CN117709050A (zh) | Hemp辐照环境中线缆脉冲耦合特征仿真方法及系统 | |
| US11893320B2 (en) | Method for backshell components in 3D routing harness and flattening route harness | |
| Arlou et al. | Computationally-effective worst-case model of wire radiation in the frequency range 1 Hz–40 GHz | |
| Mocker et al. | Comparison of electromagnetic solvers for antennas mounted on vehicles | |
| Gansen et al. | An effective 3D leapfrog scheme for electromagnetic modelling of arbitrary shaped dielectric objects using unstructured meshes | |
| US20150186590A1 (en) | Wiring board design system and wiring board design method | |
| JP4868412B2 (ja) | 解析装置、解析方法、解析プログラムおよび解析プログラムが格納された記録媒体 | |
| EP4184374A1 (en) | Method for maintaining 3d orientation of route segments and components in route harness flattening | |
| Helmers et al. | Analyzing electromagnetic pulse coupling by combining TLT, MoM, and GTD/UTD | |
| JP2005291948A (ja) | 電磁波解析装置及びその解析方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20021029 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20030212 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080307 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090307 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100307 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100307 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110307 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110307 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120307 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130307 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130307 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140307 Year of fee payment: 11 |
|
| LAPS | Cancellation because of no payment of annual fees |