JP3853602B2 - 電磁界強度算出装置への入力データ作成方法および作成装置 - Google Patents
電磁界強度算出装置への入力データ作成方法および作成装置 Download PDFInfo
- Publication number
- JP3853602B2 JP3853602B2 JP2001089735A JP2001089735A JP3853602B2 JP 3853602 B2 JP3853602 B2 JP 3853602B2 JP 2001089735 A JP2001089735 A JP 2001089735A JP 2001089735 A JP2001089735 A JP 2001089735A JP 3853602 B2 JP3853602 B2 JP 3853602B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- mesh
- electromagnetic field
- field intensity
- points
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/0864—Measuring electromagnetic field characteristics characterised by constructional or functional features
- G01R29/0892—Details related to signal analysis or treatment; presenting results, e.g. displays; measuring specific signal features other than field strength, e.g. polarisation, field modes, phase, envelope, maximum value
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
【発明の属する技術分野】
本発明は電気回路装置から放射される電磁界強度の算出方式に係り、更に詳しくはモーメント法を用いて電気回路装置から放射される電磁界の強度を算出する電磁界強度算出装置への入力データ作成方法、および作成装置に関する。
【0002】
【従来の技術】
任意の形状の物体、例えば電気回路装置から放射される電磁波の電磁界強度は、物体の各部分に流れる電流が分かれば容易に計算することができる。この電流値は、理論的にはマクスウェルの電磁界方程式を与えられた境界条件の下で解くことによって与えられるが、任意の形状の物体を対象とした複雑な境界条件の下での理論的な解法は存在しない。
【0003】
従って、現在の電磁界強度算出装置で用いられている電流の計算法は全て近似的なものである。この近似的な解法として、微小ループアンテナ近似法と、分布定数線路近似法、およびモーメント法の3つが知られている。
【0004】
微小ループアンテナ近似法は、波源回路と負荷回路との間を接続する配線をループアンテナとして扱い、ループ上の電流は平坦なものと仮定して、その電流を集中定数回路の計算手法で求める方法である。この方法の計算は最も簡単であるが、ループの寸法が電磁波の波長に比べて無視できない条件では精度が悪く、現実にはほとんど用いることができない。
【0005】
次に分布定数線路近似法は、一次元の構造物として近似できる物体に対して分布定数線路の方程式を適用することによって電流を求める方法である。計算は比較的簡単であり、電流を高速、かつ高精度に求めることができるが、解析対象の物体が一次元の構造物として近似できない場合、例えばプリント板や金属筐体などの任意の形状の物体については、電流を求めることができない。
【0006】
これに対してモーメント法は、マクスウェルの電磁界方程式から導かれる積分方程式の解法の1つであり、三次元の任意形状の物体を扱うことができる。モーメント法では物体は小さな要素(パッチ、またはメッシュ)に分割され、電流の計算が行われる。
【0007】
モーメント法を用いて電気回路装置から放射される電磁界の強度を算出する場合には、電気回路装置内のプリント板、ケーブル、リード、および金属筐体など全ての構成要素を組合わせて、かつ正確にモデル化することによって、電磁界強度を高精度で計算することができる。
【0008】
プリント板、ケーブル、およびリードなどのデータは、CADから三次元データや配線のデータを取り込んで、比較的容易に解析モデルを作成することができる。しかしながら電気回路装置の金属筐体などの三次元構造データに対しては、CADからデータを取り込んだ後で、モーメント法に適したメッシュ(パッチ)への分割を行ってメッシュモデルを作成する必要があり、そのモデルを正確に、かつ効率的に作成する必要がある。
【0009】
【発明が解決しようとする課題】
このように電気回路装置の金属筐体などをモーメント法に適するメッシュに分割したい場合にも、従来は自動的な分割法としては3点で構成される三角形の形状のメッシュ(パッチ)への分割が行われていた。後述するようにモーメント法に最も適したメッシュ分割法は四角形へのメッシュ分割であり、従来は金属筐体などを4点で構成される四角形のメッシュへの分割を自動化することが困難であった。
【0010】
すなわち従来は、四角に切れないようなところや、穴のある場所については、CADのデータを手動で修正して、四角形メッシュへの分割に適した形式にCADモデルを修正する必要があり、そのための作業に手間がかかるという問題点があった。
【0011】
図35、および図36はモーメント法を用いた電流計算におけるメッシュ形状の違いの影響の説明図である。図35はメッシュ(パッチ) 内での電流計算の相違の説明図であり、三角形メッシュの場合には頂点から対辺方向の電流 (方向を矢印で示す)が計算されるのに対して、四角形メッシュの場合には対辺方向の電流が計算されることになる。
【0012】
図36はこのようなメッシュを組合わせた場合の全体としての電流方向の説明図である。同図(a)のように三角形のメッシュに分割した場合には、電流の流れがジグザグ、すなわち不均一となり、電流の流れる経路の長さが長くなるために、電流の伝搬遅延が発生し、モーメント法による解析精度が低下することになる。これに対して(b)のような四角形メッシュへの分割の場合には、電流がスムーズに流れ、モーメント法による解析の精度が高くなる。
【0013】
本発明の課題は、上述の問題点に鑑み、電気回路装置の金属筐体などをモーメント法の電流計算に適する四角形形状のメッシュに自動分割し、分割されたデータを電磁界強度算出装置に与えることによって、全体として電磁界強度算出の解析精度を高めると共に、その処理を効率化することである。
【0014】
【課題を解決するための手段】
図1は本発明の原理的な機能ブロック図である。同図は金属筐体を有する電気回路装置から放射される電磁界の強度を算出する電磁界強度算出装置への入力データの作成方法の機能ブロック図である。
【0015】
図1において、本発明の入力データ作成方法では、1で電気回路装置の三次元データから金属筐体の面データを抽出し、2でその面データに対応する面を四角形の形状のメッシュに分割し、3でそのメッシュ分割データを電磁界強度算出装置に出力する。
【0016】
発明の実施の形態においては、金属筐体の面データが金属筐体を構成する金属板の表面、および裏面を含む複数の面データから構成される時には、その複数の面データから金属板の表面、または裏面に存在する複数の面から構成される同一面のデータを抽出し、その同一面を構成する各面を四角形の形状のメッシュに分割することもできる。
【0017】
またその同一面を構成する各面の面データとしての制御点のデータを用いて、その各面を四角形状に近似させる時に、その四角形状の頂点に対応する制御点を抽出することもでき、更にその各面を近似する四角形状の頂点に対応する制御点を4つの頂点とし、制御点のデータによって決定される平面、または曲面を対辺同志で等分割し、四角形形状メッシュに分割することもできる。
【0018】
また発明の実施の形態においては、金属筐体の面データが複数の面に対応するデータである時に、その複数の面のそれぞれを四角形の形状のメッシュに分割し、その複数の面のうち隣接する面で共通辺となる可能性がある2つの辺上の分割点の座標を比較し、その分割点の座標があらかじめ定められた許容誤差内で一致すると判定された時に、その分割点のデータを隣接する面における共有データとし、その共有データを含むメッシュ分割データを電磁界強度算出装置に出力することもできる。
【0019】
更に実施の形態においては、金属筐体の面データに対応する面の四角形形状へのメッシュ分割後に、その金属筐体と異なる材質であって、対応する面に重畳し、対応する面より小さい面積の面が存在する時、その小さい面積の面の形状に目合わせして四角形形状のメッシュを再分割し、再分割後のメッシュ分割データを電磁界強度算出装置に出力することもできる。
【0020】
次に本発明における記憶媒体として、金属筐体を有する電気回路装置から放射される電磁界の強度を算出する電磁界強度算出装置への入力データを作成する計算機によって使用される記憶媒体であって、金属筐体を構成する面のうちで1つ以上の、面の指定をユーザから受取るステップと、電気回路装置の三次元データから指定された面の面データを抽出するステップと、その面データに対応する面を四角形の形状のメッシュに分割するステップと、そのメッシュ分割データを電磁界強度算出装置に出力するステップとを計算機に実行させるためのプログラムを格納した計算機読出し可能可搬型記憶媒体が用いられる。
【0021】
また本発明の電磁界強度算出装置への入力データを作成する装置において、電気回路装置の三次元データから金属筐体の面データを抽出する面データ抽出手段と、その面データに対応する面を四角形の形状のメッシュに分割するメッシュ分割手段と、そのメッシュ分割データを電磁界強度算出装置に出力する作成データ出力手段とを備える。
【0022】
以上のように本発明によれば、金属筐体を有する電気回路装置から放射される電磁界の強度を算出する電磁界強度算出装置への入力データを作成するために、金属筐体の面データに対応する面がモーメント法の計算に適した四角形の形状のメッシュに分割される。
【0023】
【発明の実施の形態】
図2は、本発明における電磁界強度算出装置と、電磁界強度算出装置への入力データを作成する入力データ作成装置との関係の説明図である。入力データ作成装置11の出力は電磁界強度算出装置12に与えられ、例えば電気回路装置から放射される電磁界の強度が算出される。
【0024】
本発明の実施形態としては、入力データ作成装置11は例えば電気回路装置の金属筐体をモーメント法の計算に適する四角形形状のメッシュに分割し、そのメッシュ分割データを電磁界強度算出装置12に与えるものとして実施形態を説明するが、入力データ作成装置11は金属筐体だけでなく、電気回路装置内のプリント板、ケーブル、リードなどの構成要素に対応する入力データを作成することも当然可能である。
【0025】
電磁界強度算出装置は、例えば金属筐体に対するメッシュ分割データの入力に対応して、分割された各メッシュ (パッチ)の自己インピーダンスや、メッシュの間における相互インピーダンスなどを計算するインピーダンス計算部13、インピーダンス計算部13の計算結果と、例えばユーザから与えられる波源の電圧を用いて電気回路装置の各部における電流を計算する電流分布計算部14、求められた電流分布を用いて、例えば電気回路装置からある距離だけ離れた点における電界、磁界などを計算する電磁界計算部15によって構成されている。
【0026】
図3は本発明の実施形態における電磁界強度算出装置への入力データ作成装置の構成ブロック図である。
図3において、CADデータファイル20から三次元構造データ読み込み部21によって金属筐体の三次元構造データが読み込まれ、メモリ部22に面データとして全ての面のデータが順次格納され、かつディスプレイ表示部23によって各面の形状が表示される。例えばユーザによって、キーボード/マウス入力部24を介してメッシュ分割の対象となる面がクリック選択され、同一面の抽出処理部25に与えられる。
【0027】
同一面の抽出処理部25は、CADデータファイル20に格納されている各面のデータがブロックデータである時に、同一面の抽出を行うものである。このブロックデータでは、金属筐体を構成する各金属板について表面、裏面、および側面のデータが格納されている。この場合、金属板の表面、または裏面のいずれかのみをメッシュ分割すればよく、板厚方向の反対側の面、および側面についてはメッシュ分割が不必要であることから、表面、または裏面を構成する複数の面を同一面として抽出するものである。CADデータファイル20に格納されているデータが金属板の重みを無視したサーフェスデータである場合には、この抽出処理は不必要となる。同一面の抽出処理の詳細については、更に後述する。
【0028】
同一面の抽出処理部25の処理結果は、 4点ブロック面の抽出とメッシュ分割処理部26に与えられる。本実施形態では、1つの面をできるだけ四角形で近似し、その四角形を対辺等分割で四角形形状のメッシュに分割するが、4点ブロック面抽出はこのメッシュ分割前の面の形状をできるだけ四角形に近い形になるように4つの頂点を抽出するものである。そして抽出された4点ブロック面については、例えばユーザによってキーボード/マウス入力部24から入力された分割長を用いて、対辺等分割が行われる。
【0029】
次に隣接面ポイント共有処理部27によって、隣接面の共有辺に対応する分割点の座標の比較が行われ、例えばあらかじめ定められた許容誤差内で一致すると判定された分割点が隣接面の間で共有データとされ、隣接面がつながったメッシュが作成される。
【0030】
次にメッシュ形状の目合わせ処理部28によって、材料の異なる面の重畳部分に対応して、メッシュの再分割が行われる。例えば金属板に穴があるような場合にも、その穴の形状に合わせてメッシュの再分割が行われる。
【0031】
作成されたメッシュデータはメッシュデータ書き込み部29によってメッシュデータファイル30に格納され、このファイルの内容が電磁界強度算出装置への入力データとして出力され、筐体の透き間などからの電波のもれや、シールド解析などに用いられる。三次元構造データ読み込み部21からメッシュデータ書き込み部29までの処理は、必要に応じてメモリ部22の作業領域を用いて実行される。
【0032】
図3の入力データ作成装置の動作について、図4〜図11を用いて更に説明する。図4および図5はCADデータファイル20の格納データの例である。図4は金属筐体を構成する面、一般的には曲面の例であり、曲面はBスプライン曲面として、曲面の制御点によって表現される。各制御点はT方向のパラメータuと、S方向のパラメータv、および重みによって表現される。すなわちx,y,z座標を持つ基準点が存在し、各制御点は基準点からの面単位の(変位)座標u,v を持つ。
【0033】
図5は、図4に対応するCADデータのデータ構造である。このデータはBスプライン曲面の面番号から、各制御点の空間座標値X、Y、Zなどを含んでいる。
【0034】
図6は、4点ブロック面の抽出とメッシュ分割処理部26によるメッシュ分割の説明図である。CADデータファイル20に格納されている三次元構造データから金属筐体の面データが抽出され、その面データに対応する面が筐体面に沿って四角形形状にメッシュ分割される。この四角形としては正方形又は長方形が最適であるが、台形やひし形などであってもよい。
【0035】
図7は同一面の抽出処理部25による処理の説明図である。図3のキーボード/マウス入力部24を介して、ユーザから金属筐体を構成する複数の面のうちいずれかの面が指定(選択)されるが、指定された面が金属筐体の表面、または裏面のいずれにあるかによって表面、または裏面に存在する同一面のデータが、同一面の抽出処理部25によって抽出される。
【0036】
図8、および図9は4点ブロック面の抽出とメッシュ分割処理部26による処理の説明図である。図8においては、例えば図7で抽出された同一面を構成する各面のそれぞれを表現する制御点から、その面を四角形に近似するための4つの頂点となる制御点が抽出される。
【0037】
図9においては、図8のようにして抽出された4つの制御点が四角形の4つの頂点とされ、この4つの頂点によって構成される平面、または曲面を対辺同志で等分割する処理が行われる。
【0038】
図10は隣接面ポイント共有処理部27による処理の説明図である。隣接面が接続される場合の共有辺となると考えられる、それぞれの辺上にあるメッシュ分割点の座標値が許容誤差以内で同じと判断される場合、それらの面はその辺で接続されたものとみなし、分割点の共有が行われる。
【0039】
図11は材質の異なる2つの面が重畳する場合に、メッシュ形状の目合わせ処理部28によって行われる処理の説明図である。同図において、異なる材料の面が重畳している場合に、小さい面積の面の形状に目合わせしてメッシュの再分割が行われる。例えば図11の右側で、9つの大きな正方形のメッシュのうちで最も左上の正方形のメッシュは、3つの長方形のメッシュに再分割されることになる。
【0040】
図12は、図3の同一面の抽出処理部25によって実行される同一面抽出処理の詳細フローチャートである。同図において処理が開始されると、ステップS1で例えばユーザによって金属筐体を構成する1つの面が基準面としてマウスクリックされ、ステップS2で基準面の選択が行われ、ステップS3で図3のメモリ部22から面データが出力されて、ステップS4で基準面、すなわち最初の処理対象面、または処理対象面に隣接する面が抽出される。
【0041】
ここではまず基準面が抽出され、ステップS5でその面を重心と金属面の外側方向への法線の算出が行われ、ステップS6でこの法線と逆方向のベクトルが算出され、ステップS7で面データを用いてその逆方向に交差する面が探索される。
【0042】
このような処理について、図13〜図22を用いて説明する。図13は金属筐体の曲面の例であり、図18は平面の例である。図13の例で、上側の左側の面がユーザによって基準面としてマウスクリックされると、図14に示すように選択された基準面Aに対する法線ベクトルと、その逆方向ベクトルとが算出される。
【0043】
図12のステップS7でこの逆方向ベクトルに交差する面が探索されるが、図14では基準面Aの下側の面が探索されることになる。図12のステップS8で、その交差する面との距離があらかじめ定められた許容値以内か否かが判定され、ここでは距離は金属板の厚さであり、許容値以内と判定されるとすると、ステップS9で交差する面がすでに同一面と判断された面であるか否かが判定される。
【0044】
図13において同一面と判断されるべき面は基準面と基準面の右側の面すなわち上側の右側の面であるが、ここではまだ処理が開始されたばかりであり、交差する面は同一面と判断された面でないと判定され、ステップS10ですでに除去された面であるか否かが判定される。ここでも交差する面はすでに除去された面ではないと判断され、ステップS11でその面に対して除去フラグがセットされ、基準面Aの下側の面が除去される。
【0045】
そしてステップS12で処理対象面、ここでは基準面Aに対して同一面フラグがセットされ、ステップS13でその面、ここでは基準面Aがディスプレイ上でハイライト表示されて、ステップS14で基準面に対して次の隣接面があるか否かが判定される。
【0046】
図15に示すように基準面Aに対して隣接面Bが存在し、ステップS4以降の処理がその隣接面を処理対象面として実行される。処理対象面Bの法線ベクトルと逆方向のベクトルと交差する面が同様に探索され、ステップS8〜ステップS10の処理を経て、ステップS11でこの交差する面に対しても除去フラグがセットされ、処理対象面、すなわち面Bに対してステップS12で同一面フラグがセットされ、ステップS13でその処理対象面がハイライト表示される。
【0047】
続いてステップS14で次の隣接面があるか否かが判定される。ここで隣接面とは、直前の処理対象面に必ずしも隣接している必要はなく、まだ処理が行われていない面と考えることもでき、図16に示すように例えば側面Cが選択され、ステップS4以降の処理が行われる。
【0048】
前述と同様にしてステップS7で逆方向ベクトルに交差する面が探索されるが、ここではその交差する面はユーザによって選択された基準面としての面Aとなる。この面Aに対してはすでに同一面フラグがセットされており、ステップS9ですでに同一面と判断された面と判定され、ステップS15で処理対象面、すなわち側面C自身に対して除去フラグがセットされ、側面Cは除去される。
【0049】
続いてステップS14で、処理が残っている面として図17に示すように側面Dが選択され、ステップS4以降の処理が実行される。ステップS7で逆方向ベクトルに交差する面として探索される面は図15で説明したようにすでに除去された面であり、ステップS10でそのことが判定されて、ステップS15で処理対象面、すなわち側面Dも除去されて、全ての面に対する処理を終了する。前述のように基準面と同一面と判断される面は、図13の右側に示すように上側の2つの面となる。
【0050】
図18の平面に対する処理の説明が図19〜図22であるが、図14〜図17と比較して異なる点は図21、および図22の処理である。図21では左側の側面Cを処理対象面として処理が行われるが、図12のステップS8で逆方向ベクトルと交差する面、ここでは面Dとの距離が許容値より大きいと判断され、ステップS15で処理対象面、すなわち面C自体が除去されることになる。
【0051】
図12でステップS8の処理が行われないか、または交差する面との距離が金属板の厚さよりも大きく、かつ幅 (側面の間の距離)より小さく設定されていない場合には、ステップS9〜S11の処理を経てステップS12で処理対象面、ここでは側面Cに同一面フラグがセットされ、側面Cが基準面Aと同一面であると判定されてしまうことになる。
【0052】
図22においては、右側の側面Dを処理対象面とする処理が実行されるが、側面Cに対すると同様に、ステップS8で交差する面との距離が許容値より大きいと判定されて、ステップS15で処理対象面、すなわち側面Dに対して除去フラグがセットされることになる。
【0053】
なお金属筐体板が図13のように曲面を有する場合には、本実施形態で説明した同一面抽出処理が正しい結果を常に生ずるとは必ずしも言えない。例えば、図14で基準面Aに対して同一面フラグがセットされた後に,図15で面Bではなく、側面Cが処理対象面とされ、かつ面Cに対する逆方向ベクトルが基準面Aでなく、面Bに交差する場合には側面Cも同一面とされる可能性がある。このように曲面の形状や、同一面を構成する各面の配置、図12のステップS8における交差する面との距離に対する許容値の設定などによっては正しい結果が得られないことに注意を要する。
【0054】
図23は図3の4点ブロック面の抽出とメッシュ分割処理部26の処理の詳細フローチャートである。同図において処理が開始されると、ステップS20におけるユーザからのメッシュの分割長入力と、ステップS21における同一面データ出力とを用いて、ステップS22で同一面から同一面を構成する複数の面の1つが取り出され、ステップS23で同一面を構成する複数の面の全てに対して処理が終了したか否かが判定され、終了していない場合には、ステップS24で取り出された1つの面に対してその面の重心からその面を表現する複数の制御点までの距離が算出され、ステップS25で距離が1番大きい制御点が基準点とされる。
【0055】
次のステップS26〜S30が4点ブロック面の抽出処理に相当するが、この処理について図24〜図27の例を用いて説明する。
ステップS26で重心を中心として基準点、図24では点Aから180度の範囲内にある全ての制御点が1つずつ順次抽出される。そしてステップS27で抽出された制御点、重心、および基準点の3点から構成される三角形の面積が算出される。ステップS28で180度の範囲内にある全ての制御点が抽出されたか否かが判定され、全てが抽出されていない場合にはステップS26からの処理が繰返される。
【0056】
180度の範囲内にある全ての制御点が抽出されたと判定されると、ステップS29で面積が1番大きい三角形に対応する制御点が次の基準点とされ、ステップS30で4点ブロック面の4点、すなわち四角形の4つの頂点が基準点として抽出されたか否かが判定され、まだ抽出されていない場合にはステップS26からの処理が繰返される。
【0057】
図24ではステップS25で基準点とされた点、すなわち重心からの距離が最も遠いと判定された点Aから180度の範囲内にある全制御点として点B,C、およびDがステップS26からS28のループで順次抽出され、ステップS27でそれぞれの三角形の面積が算出され、ステップS29で面積が1番大きい三角形に対応する制御点、ここではBが次の基準点とされる。なお図24で重心を点Oとする時、三角形AOBの面積の方が、三角形AOCの面積より大きいものとする。
【0058】
図25は点Bを基準点としたステップS26〜S30の処理の説明図である。重心Oを中心として基準点Bから180度の範囲内にある制御点としてC,D、およびEが抽出され、面積の最も大きい三角形に対応する制御点Eが次の基準点とされる。
【0059】
図26においては基準点Eから180度の範囲内にある制御点の抽出が行われる。このように重心を中心として基準点から180度の範囲内にある制御点を抽出する場合、すでに基準点とされた制御点は抽出しないことにする。そこでS29で点Fが次の基準点とされるが、ステップS30ですでに4点が抽出されたと判定され、ステップS31で4つの辺の分割数が算出される。
【0060】
この分割数はステップS20でユーザから入力されたメッシュの分割長、すなわちメッシュの辺の長さに対応して算出されるが、4点ブロックとしての四角形の各辺の長さを指定された分割長で割った時に剰余が生ずる場合には、分割数に1が加算される。
【0061】
なおユーザから入力されるメッシュの分割長(メッシュの辺の長さ)はモーメント法の計算に適するように、例えば最大の周波数における波長の1/20程度とされる。また対辺の分割数が異なる場合には長い方の辺の分割数に合わせて短い辺の分割も行われる。
【0062】
そしてステップS32で、図9で説明したような対辺等分割が行われる。図9では一方の対辺は5分割され、また他方の対辺はそれぞれ3分割されいる。そしてステップS33でメッシュデータの保存が行われる。メッシュデータが保存された後にステップS22で同一面のデータから残っている1つの面が取り出され、ステップS23以降の処理が実行されるが、ステップS22で残りの面が取り出せない時、すなわちステップS23で全ての面に対する処理が終了したと判定された時点で処理を終了する。
【0063】
図28、および図29は図3のメッシュデータファイル30に格納されるメッシュデータのデータ構造の説明図である。図28は金属筐体の曲面に対応する四角形メッシュ(パッチ)の例を示し、それぞれのメッシュはメッシュ分割点に相当するポイントによって表現されている。
【0064】
図29はメッシュデータ構造の例である。メッシュデータはポリゴン、すなわち多角形の頂点の座標指定データと、ポリゴン構成点の指定データとからなっている。座標指定データは、図28のポイントのそれぞれに対するポイント番号、x,y,zの座標値から構成され、またポリゴン構成点の指定データは、1つのポリゴン (四角形メッシュ,パッチ)に対するパッチ番号と、そのパッチの4つの頂点となるポイントのポイント番号から構成されている。
【0065】
図30は隣接面の間におけるメッシュ分割点共有処理、すなわち図3の隣接面ポイント共有処理部27によると処理の詳細フローチャートである。同図において処理が開始されると、ステップS41でメッシュデータからメッシュの分割点の座標値が取り出され、ステップS42で全ての分割点の座標値に対する処理が終了したか否かが判定され、まだ終了していない場合にはステップS43で隣接面を構成する分割点の座標値との間で間隔が算出され、ステップS44でその間隔が許容値以内であるか否かが判定され、許容値より大きい場合には隣接面の他の分割点の座標値との間隔算出以降の処理が繰返される。この座標値の間の間隔の算出処理においては、図29で説明したように各分割点に対応するポイントのx,y,z座標値がすでに求められており、この座標値を用いて間隔が算出される。
【0066】
ステップS44で座標値の間隔が許容値以内と判定されると、ステップS45でその2つの分割点が共有点としてセットされ、メッシュデータとして格納される。続いてステップS41に戻り、他のメッシュ分割点の座標値が取り出され、ステップS42以降の処理が繰返される。ステップS41で取り出されるべきメッシュ分割点の座標値がない場合には、ステップS42で全ての分割点に対する処理が終了したと判定されて、ポイント共有処理を終了する。
【0067】
図31はポイント共有処理の結果の例の説明図である。同図のように、2つの面のそれぞれの辺の間で、2つの辺の分割点の間の間隔が許容値以内であると判定されると、その2つの辺は2つの面の接続辺と判定され、その2つの辺の分割点は共有点としてセットされることになる。
【0068】
図23において4点ブロック面の4つの頂点を基準点として次々に抽出する処理を説明したが、隣接する4点ブロック面との間に共有点が存在する場合には、次の基準点としてその共有点が優先的に抽出される。
【0069】
図32は図3のメッシュ形状の目合わせ処理部28による処理の詳細フローチャートである。この処理では、異なる材料の面が重畳している場合に、それぞれの面のメッシュが目合わせされるようにメッシュの再分割が行われる。
【0070】
図33はこの目合わせの例の説明図である。同図においては、対辺がそれぞれ3分割された細い実線のメッシュの面と、その面とは材質が異なり、面積が小さく、太い実線で示される面とが重畳している。
【0071】
図33で目合わせを行うために、太い実線の面のエッジを大きな面のメッシュの端まで延長し、延長線が通るメッシュを延長線に対応して再分割することにより、目合わせを実行する。これによって再分割する前のメッシュは削除される。例えば図33で再分割する前の細い実線のメッシュのうち、最も右上のメッシュは6個のメッシュに再分割されることになる。このような目合わせの方法は、例えば金属面に穴や、すき間などがあいている場合にも適用することができ、穴やすき間が重畳するメッシュは、目合わせした後でその穴やすき間に対応するメッシュを電気伝導度0として扱うことが可能となる。
【0072】
図32において処理が開始されると、まずステップS51でメッシュデータと全ての面データから、メッシュ領域内で重畳する異なる材料の面が取り出され、ステップS52でそのような全ての面に対する処理が終了したか否かが判定され、まだ終了していない場合にはステップS53で重畳する面のエッジが延長され、ステップS54でエッジの延長線上にあるメッシュが再分割され、ステップS55で再分割されたメッシュがメッシュデータに追加され、ステップS56で全てのエッジに対する処理が終了したか否かが判定され、まだ終了していない場合には次のエッジに対してステップS53以降の処理が繰返される。
【0073】
ステップS56で取り出された面の全てのエッジに対する処理が終了したと判定されると、ステップS51に戻って重畳している次の面が取り出され、ステップS52以降の処理が繰返される。ステップS51で重畳する次の面が取り出されなかった場合には、ステップS52で全ての面に対する処理が終了したと判定されて、目合わせ処理を終了する。
【0074】
以上説明したような電磁界強度算出装置への入力データの作成は、当然一般的なコンピュータによって実行することができる。図34はそのようなコンピュータシステムの構成ブロック図である。同図においてコンピュータ51は、本体52とメモリ53とによって構成されている。メモリ53としては、ランダムアクセスメモリ(RAM)、ハードディスク、磁気ディスクなどの各種の記憶装置を使用することができ、このようなメモリに本発明の特許請求の範囲の請求項7、および8のプログラムや、図12、図23、図30、および図32のフローチャートに示したプログラムなどが格納され、そのプログラムが本体52によって実行されることにより本実施形態における電磁界強度算出装置への入力データを作成すること、すなわち電気回路装置の金属筐体の金属面ををモーメント法の計算に適する四角形形状メッシュに分割し、分割後のメッシュデータを電磁界強度算出装置に出力することが可能となる。
【0075】
このようなプログラムは、プログラム提供者側からネットワーク54を介して送られ、そのプログラムがコンピュータ51によって実行されることも可能である。また更にそのようなプログラムは市販され、流通している可搬型記憶媒体55に格納され、そのような可搬型記憶媒体55がコンピュータ51にロードされることによっても実行可能である。可搬型記憶媒体55としてはCD−ROM、フロッピーディスク、光ディスク,光磁気ディスクなど様々な形式の記憶媒体を使用することができ、そのような記憶媒体に前記のようなプログラムが格納され、そのプログラムがコンピュータ51により実行されることによって、本実施形態で説明した電磁界強度算出装置への入力データの作成が可能となる。
【0076】
【発明の効果】
以上詳細に説明したように、本発明によればCADの三次元構造データからモーメント法の計算に適した金属筐体モデル、すなわち四角形形状のメッシュに分割されたモデルを、自動的、かつ高速に作成することが可能となり、そのようなモデルを電磁界強度算出装置に出力することによって、電気回路装置の金属筐体を含めた装置全体から放射される電磁界の強度を正確に計算することが可能となり、電磁界強度算出装置の実用性の向上に寄与するところが大きい。
【図面の簡単な説明】
【図1】本発明の原理的な機能ブロック図である。
【図2】本発明の入力データ作成装置と電磁界強度算出装置の関係を示す図である。
【図3】入力データ作成装置の構成を示すブロック図である。
【図4】CADデータファイルにおける曲面の表現方法を説明する図である。
【図5】CADデータファイルのデータ構造を示す図である。
【図6】金属筐体の筐体面の四角形形状へのメッシュ分割を説明する図である。
【図7】金属筐体の表面、または裏面に存在する同一面を説明する図である。
【図8】同一面を構成する各面を表現する制御点から各面を四角形形状に近似するための4つの頂点を抽出する方法の説明図である。
【図9】図8で求められた4つの点を頂点とする平面、または曲面を対辺同志で等分割する方法の説明図である。
【図10】隣接する面の接続辺において分割点を共有する方法を説明する図である。
【図11】異なる材料の面が重畳する場合にメッシュを再分割する方法を説明する図である。
【図12】同一面抽出処理の詳細フローチャートである。
【図13】曲面に対する同一面を説明する図である。
【図14】曲面を持つ金属板の面データから同一面を抽出する処理の経過(その1)を説明する図である。
【図15】曲面を持つ金属板の面データから同一面を抽出する処理の経過(その2)を説明する図である。
【図16】曲面を持つ金属板の面データから同一面を抽出する処理の経過(その3)を説明する図である。
【図17】曲面を持つ金属板の面データから同一面を抽出する処理の経過(その4)を説明する図である。
【図18】平面状の金属板における同一面を説明する図である。
【図19】平面状の金属板における同一面抽出処理の経過(その1)を説明する図である。
【図20】平面上の金属板における同一面抽出処理の経過(その2)を説明する図である。
【図21】平面上の金属板における同一面抽出処理の経過(その3)を説明する図である。
【図22】平面上の金属板における同一面抽出処理の経過(その4)を説明する図である。
【図23】4点ブロック面の抽出とメッシュ分割処理の詳細フローチャートである。
【図24】4点ブロック面の抽出例を説明する図(その1)である。
【図25】4点ブロック面の抽出例を説明する図(その2)である。
【図26】4点ブロック面の抽出例を説明する図(その3)である。
【図27】4点ブロック面の抽出例を説明する図(その4)である。
【図28】メッシュ分割データにおける曲面の表現方法を説明する図である。
【図29】メッシュ分割データのデータ構造を示す図である。
【図30】隣接面ポイント共有処理の詳細フローチャートである。
【図31】隣接面ポイント共有処理の具体例の説明図である。
【図32】メッシュ形状の目合わせ処理の詳細フローチャートである。
【図33】メッシュ形状の目合わせの例を説明する図である。
【図34】本実施形態におけるプログラムのコンピュータへのローディングを説明する図である。
【図35】メッシュ形状による電流計算の相違を説明する図である。
【図36】三角形メッシュと四角形メッシュとにおける電流方向の相違を説明する図である。
【符号の説明】
11 入力データ作成装置
12 電磁界強度算出装置
13 インピーダンス計算部
14 電流分布計算部
15 電磁界計算部
20 CADデータファイル
21 三次元構造データ読み込み部
22 メモリ部
23 ディスプレイ表示部
24 キーボード/マウス入力部
25 同一面の抽出処理部
26 4点ブロック面の抽出とメッシュ分割処理部
27 隣接面ポイント共有処理部
28 メッシュ形状の目合わせ処理部
29 メッシュデータ書き込み部
30 メッシュデータファイル
Claims (4)
- 金属筐体を有する電気回路装置から放射される電磁界強度を算出する電磁界強度算出装置への入力データを作成する装置において、
前記電気回路装置の三次元データから、前記金属筐体の面データを抽出する面データ抽出手段と、
該面データに対応する面を、該面データが、該金属筐体を構成する金属板の表面、および裏面を含む複数の面データから構成される時、該復数の面データから、該金属板の表面、または裏面に存在する複数の面から構成される同一面のデータを抽出し、該同一面を構成する複数の面の各面の面データとしての制御点のデータを用いて、該各面を四角形状に近似させる時、該四角形状の頂点に対応する制御点を抽出し、該同一面を構成する該金属板の表面、または裏面に存在する複数の面の各面を四角形の形状のメッシュに分割するメッシュ分割手段と、
該メッシュ分割データを前記電磁界強度算出装置に出力する作成データ出力手段と
を備えることを特徴とする電磁界強度算出装置への入力データ作成装置。 - 前記各面を近似する四角形状の頂点に対応する制御点を4つの頂点とし、前記制御点のデータによって決定される平面、または曲面を、2つの制御点が規定する辺と、この辺に対向する他の2つの制御点で規定される対辺について、対辺同士で同じ個数に等分割し、前記四角形状メッシュに分割することを特徴とする請求項1記載の電磁界強度算出装置への入力データ作成装置。
- 前記金属筐体の面データが複数の面に対応するデータである時、該複数の面のそれぞれを四角形の形状のメッシュに分割し、該複数の面のうち隣接する面で共通辺となる可能性がある2つの辺上の分割点の座標を比較し、該分割点の座標があらかじめ定められた許容誤差内で一致すると判定された時、該分割点のデータを隣接する面における共有データとし、該共有データを含むメッシュ分割データを前記電磁界強度算出装置に出力することを特徴とする請求項1記載の電磁界強度算出装置への入力データ作成装置。
- 前記金属筐体の面データに対応する面の四角形形状へのメッシュ分割後に、該全体筐体と異なる材質であり、該対応する面に重畳し、該対応する面より小さい面積の面が存在する時、該小さい面積の面の形状に目合わせして前記四角形形状のメッシュを再分割し、該再分割後のメッシュ分割データを前記電磁界強度算出装置に出力することを特徴とする請求項1記載の電磁界強度算出装置への入力データ作成装置。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001089735A JP3853602B2 (ja) | 2001-03-27 | 2001-03-27 | 電磁界強度算出装置への入力データ作成方法および作成装置 |
US10/057,887 US6766261B2 (en) | 2001-03-27 | 2002-01-29 | Method and device generating data input to an electromagnetic field intensity calculating device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001089735A JP3853602B2 (ja) | 2001-03-27 | 2001-03-27 | 電磁界強度算出装置への入力データ作成方法および作成装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2002288241A JP2002288241A (ja) | 2002-10-04 |
JP3853602B2 true JP3853602B2 (ja) | 2006-12-06 |
Family
ID=18944623
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001089735A Expired - Fee Related JP3853602B2 (ja) | 2001-03-27 | 2001-03-27 | 電磁界強度算出装置への入力データ作成方法および作成装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US6766261B2 (ja) |
JP (1) | JP3853602B2 (ja) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003296395A (ja) * | 2002-04-01 | 2003-10-17 | Fujitsu Ltd | 相互インミッタンス計算プログラム、相互インミッタンス計算装置、相互インミッタンス計算方法および電磁界強度計算プログラム |
JP2004094674A (ja) * | 2002-08-30 | 2004-03-25 | Fujitsu Ltd | メッシュデータ生成方法をコンピュータに実行させるためのプログラム及びメッシュデータ生成装置 |
US7222033B1 (en) | 2003-08-18 | 2007-05-22 | Steven Lynn Newson | Electromagnetic emissions and susceptibility calculating method and apparatus |
US7260810B2 (en) * | 2003-10-16 | 2007-08-21 | International Business Machines Corporation | Method of extracting properties of back end of line (BEOL) chip architecture |
JP4528024B2 (ja) | 2004-04-28 | 2010-08-18 | 富士通株式会社 | 回路解析方法を実行させるためのプログラム |
JP4654020B2 (ja) * | 2004-12-21 | 2011-03-16 | 富士通株式会社 | メッシュデータ作成装置,メッシュデータ作成プログラム,同プログラムを記録したコンピュータ読取可能な記録媒体及びメッシュデータ作成方法並びに構造解析装置 |
JP4450751B2 (ja) | 2005-03-17 | 2010-04-14 | 富士通株式会社 | メッシュモデル作成方法、シミュレーション装置及びプログラム |
JP4879163B2 (ja) | 2005-03-22 | 2012-02-22 | 富士通株式会社 | モデリング方法及び装置、プログラム及び記憶媒体 |
JP4499781B2 (ja) | 2005-03-28 | 2010-07-07 | 富士通株式会社 | 電磁界強度算出方法、電磁界強度算出装置、制御プログラム |
WO2007043160A1 (ja) * | 2005-10-07 | 2007-04-19 | Mitsubishi Denki Kabushiki Kaisha | 電磁界分布計算方法、電磁界分布計算装置、および磁界分布計算プログラム |
JP4957080B2 (ja) * | 2006-05-31 | 2012-06-20 | パナソニック株式会社 | 数値解析用メッシュ作成装置 |
JP4920470B2 (ja) * | 2007-03-28 | 2012-04-18 | 三菱電機株式会社 | 交番磁界解析方法、および交番磁界解析プログラム |
WO2009004675A1 (ja) * | 2007-06-29 | 2009-01-08 | Fujitsu Limited | 設計支援装置、設計支援プログラム、設計支援方法 |
JP5089478B2 (ja) * | 2008-04-28 | 2012-12-05 | キヤノン株式会社 | 解析モデル作成装置及び解析モデル作成方法 |
US8949102B2 (en) | 2011-02-24 | 2015-02-03 | Cadence Design Systems, Inc. | Method and system for power delivery network analysis |
CN103136404B (zh) * | 2011-12-02 | 2016-01-27 | 深圳光启高等理工研究院 | 一种获得超材料折射率分布的方法和装置 |
CN103164555B (zh) * | 2011-12-14 | 2016-03-09 | 深圳光启高等理工研究院 | 一种人工电磁材料设计数据的优化方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04222061A (ja) | 1990-12-25 | 1992-08-12 | Babcock Hitachi Kk | メッシュ生成方法 |
JP3091815B2 (ja) * | 1994-02-25 | 2000-09-25 | 富士通株式会社 | 電磁界強度算出装置 |
JP2768900B2 (ja) | 1994-05-10 | 1998-06-25 | 富士通株式会社 | 電磁界強度算出装置 |
JP3825133B2 (ja) * | 1997-04-30 | 2006-09-20 | 富士通株式会社 | 電磁界強度算出装置及び方法並びにプログラム記録媒体 |
US5966524A (en) * | 1997-07-24 | 1999-10-12 | Lucent Technologies Inc. | 3-D electromagnetic infinite element |
-
2001
- 2001-03-27 JP JP2001089735A patent/JP3853602B2/ja not_active Expired - Fee Related
-
2002
- 2002-01-29 US US10/057,887 patent/US6766261B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
US6766261B2 (en) | 2004-07-20 |
US20020140435A1 (en) | 2002-10-03 |
JP2002288241A (ja) | 2002-10-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3853602B2 (ja) | 電磁界強度算出装置への入力データ作成方法および作成装置 | |
JP3428232B2 (ja) | 電磁界強度算出装置 | |
JP2768900B2 (ja) | 電磁界強度算出装置 | |
US20030204343A1 (en) | Electromagnetic field analysis method based on FDTD method, medium representation method in electromagnetic field analysis, simulation device, and storage medium | |
JP3091815B2 (ja) | 電磁界強度算出装置 | |
JP2006260285A (ja) | メッシュモデル作成方法、シミュレーション装置及びプログラム | |
Tirkas et al. | Finite-difference time-domain method for electromagnetic radiation, interference, and interaction with complex structures | |
US8065101B2 (en) | Electromagnetic field intensity calculating method and apparatus | |
Monje-Real et al. | Electric field integral equation fast frequency sweep for scattering of nonpenetrable objects via the reduced-basis method | |
US9568624B2 (en) | Method for predicting electromagnetic radiation characteristics, computer-readable recording medium and simulator | |
Strodthoff et al. | Automatic decomposition of 3D solids into contractible pieces using Reeb graphs | |
JP2004094674A (ja) | メッシュデータ生成方法をコンピュータに実行させるためのプログラム及びメッシュデータ生成装置 | |
JP3825133B2 (ja) | 電磁界強度算出装置及び方法並びにプログラム記録媒体 | |
JP2001188075A (ja) | 電磁界強度算出装置、電磁界強度算出方法およびその方法をコンピュータに実行させるプログラムを記録したコンピュータ読み取り可能な記録媒体 | |
JPH11161690A (ja) | 電磁界強度算出装置、電磁界強度算出方法、及び電磁界強度算出手段を有するプログラムを格納した記録媒体 | |
JP3443162B2 (ja) | 電磁界強度算出装置 | |
US7952578B2 (en) | Support apparatus, design support program, and design support method | |
CN111832183B (zh) | 基于八叉树结构mlfma区域分解的电磁仿真方法 | |
CN113987652A (zh) | 墙模型的生成方法和装置 | |
JP2004004054A (ja) | Fdtd法を用いた電磁界解析方法、電磁界解析における媒質表現方法、シミュレーション装置、及びプログラム | |
US12061858B1 (en) | Method and system for determining initial mesh for conductor in electronic circuit design | |
Giniyatova et al. | Parallel Algorithm for Solving Problem of Electromagnetic Wave Diffraction by a Tooth-Shaped Plate | |
CN118334276B (zh) | 基于pcb模型的混合元共形网格生成方法及生成装置 | |
JP2001013184A (ja) | 電磁界強度算出装置および記録媒体 | |
CN113987651A (zh) | 板模型的生成方法、装置和计算机设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060613 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060810 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20060905 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060906 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090915 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100915 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100915 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110915 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120915 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120915 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130915 Year of fee payment: 7 |
|
LAPS | Cancellation because of no payment of annual fees |