JP2004246503A - Device and method for calculating impedance of circuit board - Google Patents
Device and method for calculating impedance of circuit board Download PDFInfo
- Publication number
- JP2004246503A JP2004246503A JP2003034243A JP2003034243A JP2004246503A JP 2004246503 A JP2004246503 A JP 2004246503A JP 2003034243 A JP2003034243 A JP 2003034243A JP 2003034243 A JP2003034243 A JP 2003034243A JP 2004246503 A JP2004246503 A JP 2004246503A
- Authority
- JP
- Japan
- Prior art keywords
- impedance
- circuit board
- noise source
- passive element
- mounting distance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
Description
【0001】
【発明の属する技術分野】
本発明は、能動素子のようなノイズ源と受動素子とが実装された回路基板のインピーダンスを計算する装置および方法に関する。
【0002】
【従来の技術】
プリント回路基板上に実装されている素子(部品)を駆動したとき、基板の電源層−GND(グランド)層間に電流が流れる。この電流が他の素子に流れ込むとき、基板のインピーダンスに比例した電圧ノイズが発生し、他の素子の誤動作の原因となる。基板の設計段階でこのノイズを正確に計算するには、コンデンサ等の受動素子が実装された基板のインピーダンスを正確に把握する必要がある。
【0003】
従来より基板のモデル化に使用されている方法としては、メッシュ分割法や有限要素法がある。また、能動素子から見たデカップリングコンデンサを含む簡易インピーダンスモデルを元に、インピーダンス特性を計算する手法もある(例えば、特許文献1参照)。
【0004】
【特許文献1】
特開2002−41594号公報
【0005】
【発明が解決しようとする課題】
しかしながら、上述した従来のインピーダンス計算手法には、次のような問題がある。
【0006】
メッシュ分割法や有限要素法は、正確な計算を目指すほどモデル形状が複雑になり、計算に膨大な時間を要するという欠点がある。また、現在使用されている素子の中には数100MHzで動作するものが増えてきており、このような周波数領域では能動素子からコンデンサ等の受動素子までの実装距離により、基板全体の共振点(共振周波数)が大きく変化する。しかし、特許文献1の計算法では、実装距離に対する等価回路インピーダンスの変化を等式化するという手法は用いていない。
【0007】
したがって、数100MHz以上の周波数で動作する素子を含む回路基板において、受動素子の実装距離に対するインピーダンスの変化を高速かつ正確に計算する手法は知られていない。
【0008】
本発明の課題は、能動素子のようなノイズ源と受動素子とが実装された回路基板において、ノイズ源と受動素子の間の実装距離に対するインピーダンスを高速かつ正確に計算する装置および方法を提供することである。
【0009】
【課題を解決するための手段】
図1は、本発明のインピーダンス計算装置の原理図である。図1のインピーダンス計算装置は、入力手段11、解析手段12、および出力手段13を備え、ノイズ源と受動素子が実装された回路基板のインピーダンスを計算する。
【0010】
入力手段11は、回路基板上におけるノイズ源と受動素子の間の実装距離を入力する。解析手段12は、入力された実装距離をパラメータとして、ノイズ源と受動素子が実装された回路基板のモデルのインピーダンスを計算し、計算されたインピーダンスを用いて回路解析を行う。そして、出力手段13は、回路解析の結果を出力する。
【0011】
ノイズ源と受動素子が実装された回路基板のモデルとしては、例えば、回路基板および受動素子をそれぞれRLCの直列回路とみなし、回路基板の回路と受動素子の回路を直並列回路とみなした簡易モデルが用いられる。このような簡単なモデルを用いて実装距離に対するインピーダンスの変化を等式化することで、ノイズ源と受動素子が実装された回路基板のインピーダンスを高速かつ正確に計算することができる。
【0012】
また、計算されたインピーダンスによりモデルを特定して回路解析を行えば、受動素子の実装距離に対する電源ノイズの変化等を効率よく解析することができる。
【0013】
図1の入力手段11は、例えば、後述する図6の入力装置41に対応し、図1の解析手段12は、例えば、図6のデータ処理装置42に対応し、図1の出力手段13は、例えば、図6の出力装置44に対応する。
【0014】
【発明の実施の形態】
以下、図面を参照しながら、本発明の実施の形態を詳細に説明する。
本実施形態では、受動素子実装時の回路基板の周波数特性に注目し、これを簡単な直並列回路で表現することにより基板をモデル化する。このようなモデルを使用することで、実装基板のインピーダンスや基板に発生する電源ノイズを高速かつ正確に計算することが可能となる。
【0015】
まず、図2から図5までを参照しながら、基板のモデル化の方法について説明する。
図2は、プリント回路基板にパスコン(バイパスコンデンサ)を含む部品(素子)を実装した様子を示している。一般に、パスコンは、回路上のノイズを除去するために用いられる。
【0016】
DD(direct current−direct current )コンバータ21は、電圧値および電源の供給源であり、低周波パスコン22は、低周波領域(1MHz以下)で効果のあるバイパスコンデンサである。また、LSI(large−scale integration)23は、負荷となる能動素子であり、駆動されることによりノイズを発生する。高周波パスコン24は、高周波領域(1MHz以上)で効果のあるバイパスコンデンサである。
【0017】
図3は、プリント回路基板上に実装されている能動素子とパスコンの距離(実装距離)を変化させたときの、素子から見た基板インピーダンスの変化の様子を示している。図3の横軸および縦軸は、それぞれ周波数および基板インピーダンスを表し、曲線31、32、および33は、それぞれ異なる実装距離に対応するインピーダンス特性を表している。このように、パスコンの実装距離、実装数等の実装条件の変化により、基板インピーダンスの特性が変化する。
【0018】
図3を見ると、曲線31、32、および33のいずれにおいても3つの共振点(f1、f2、f3)が確認できる。このような特性は、プリント回路基板およびパスコンをRLCの直列回路とみなしたときのプリント回路基板の回路とパスコン回路の直並列共振特性とみなすことができる。したがって、プリント回路基板上に、電源、低周波パスコン、高周波パスコン、および素子(負荷LSI)が実装されたビヘイビアモデル(簡易モデル)は、図4のように表現することができる。図4のモデルにおける各構成要素の意味は以下の通りである。
▲1▼電源
電圧値および電源の供給源であるDDコンバータの等価インピーダンスモデル。
▲2▼パスコン(低周波)
低周波パスコンの等価インピーダンスモデル。パスコンが複数のときは、それらが並列接続されているものとしてパラメータ(容量値、抵抗値、インダクタンス値)を算出する。
▲3▼パスコン(高周波)
高周波パスコンの等価インピーダンスモデル。パスコンが等距離に複数実装されているときは、容量値および抵抗値を並列接続として算出する。インダクタンス値は素子からの実装距離によって変化する。
▲4▼PCB
プリント回路基板の等価インピーダンスモデル。抵抗値以外は素子の電源端子(電源/GNDピン)と高周波パスコンの距離によって変化する。
▲5▼負荷LSI
LSI素子におけるVdd端子とGND端子の間の等価インピーダンスモデル。
【0019】
ここで、▲3▼の高周波パスコンと▲4▼のPCBの等価インピーダンスモデルのパラメータは、パスコンの実装数および高周波パスコンから▲5▼の素子の電源/GNDピンまでの距離によって変化する。この変化量を定式化することで、パスコン実装条件に応じて図4のモデルパラメータを決定することができる。したがって、解析時間が速く、かつ精度の高い図4のモデルで、パスコン実装条件の変化に対するパラメータの変化を検証することができ、より速く正確にパスコンの最適解を得ることが可能となる。
【0020】
図3の3つの共振点の周波数を、低い方から順にf1、f2、f3とする。電源および低周波パスコンのインピーダンスによる共振点は、高周波パスコンおよびプリント回路基板による共振点に比べて十分に小さい。したがって、これらはプリント回路基板のどの位置に実装されていても、プリント回路基板の共振点に対する影響は小さい。しかしながら、高周波パスコンの共振点とプリント回路基板の共振点は近いため、高周波パスコンの実装条件により、図3の共振点の位置が変化する。
【0021】
また、抵抗値は共振周波数には影響を与えないので、容量値およびインダクタンス値のみを用いて、プリント回路基板と高周波パスコンの並列回路を図5のように表現し、これらのパラメータの算出方法を説明する。図5において、高周波パスコンの容量値およびインダクタンス値を、それぞれC1およびL1とし、プリント回路基板の容量値およびインダクタンス値を、それぞれC2およびL2とする。
【0022】
図5の4つのパラメータのうち、高周波パスコンの容量値C1のみが、パスコンの実装距離を変えても変化しない。また、図3の第1周波数f1は、C1およびL1により決定されるから、C1およびf1よりL1を求めることができる。パスコンの実装数、実装距離、およびパスコン1つ当りの寄生誘導値等を変化させることにより、素子の電源/GNDピンから高周波パスコンまでの実装距離dとL1との間の関係をグラフ化し、これを近似式で表現すると、次のようになる。
L1=aL1log10d+bL1 (1)
また、aL1およびbL1は、グラフより次のような近似式で表現できる。
aL1=a1nb1
bL1=b2La+b3
ここで、nはパスコン実装数、Laはパスコン1つ当りの寄生誘導値であり、a1およびb1は、基板の電源層−GND層間の誘電体の誘電率により一意的に決まる値である。(1)式を用いてL1を計算する場合、あらかじめ誘電率毎のa1およびb1を計算し、データベースとして保存しておくと、処理を高速化することができる。また、b2およびb3は、グラフより次のような近似式で表現できる。
b2=m1nm2
b3=m3nm4
m1、m2、m3、およびm4は、基板の電源層−GND層間の、誘電体の誘電率により一意的に決まる値であり、あらかじめ誘電率毎のm1、m2、m3、およびm4を計算し、データベースとして保存しておくと、処理を高速化することができる。
【0023】
以上の定式化により、実装距離d、パスコン実装数n、パスコン1つ当りの寄生誘導値La、および誘電率が決まれば、(1)式からL1を求めることができる。
【0024】
次に、図5のC2およびL2は、f2およびf3より一意的に決められる。C2およびL2とf2およびf3の関係は、図5より次のように表される。
【0025】
【数1】
【0026】
したがって、f2およびf3が決まれば、L1と(2)および(3)式よりL2およびC2を求めることができる。実装距離dに対するf2の変化をグラフ化すると、次のような曲線で近似することができる。
f2=af2d+bf2 (4)
したがって、f2は(4)式より求められる。ここで、af2およびbf2はグラフからの近似により、次式のようになる。
af2=a5La+a6
bf2=a7La+a8
ここで、a5、a6、a7、およびa8のパスコン実装数nに対する変化をグラフ化すると、次のような曲線で近似することができる。
ai=kin+mi (i=5,6,7,8) (5)
したがって、a5、a6、a7、およびa8は(5)式より求められる。ここで、kiおよびmiは基板の電源層−GND層間の誘電体の誘電率により一意的に決まる値であり、あらかじめ誘電率毎のkiおよびmiを計算し、データベースとして保存しておくと、処理を高速化することができる。
【0027】
以上の定式化により、実装距離d、パスコン実装数n、パスコン1つ当りの寄生誘導値La、および誘電率が決まれば、(4)式からf2を求めることができる。
【0028】
また、実装距離dに対するf3の変化をグラフ化すると、次のような曲線で近似することができる。
f3=af3d+bf3 (6)
したがって、f2は(4)式より求められる。ここで、af3およびbf3は次式のようになる。
af3=a9La+a10
bf3=b4La b5
ここで、a9およびa10のLaに対する変化をグラフ化すると、次のような曲線で近似することができる。
ai=kin+mi (i=9,10) (7)
したがって、a9およびa10は(7)式より求められる。ここで、ki、mi、b4、およびb5は、基板の電源層−GND層間の誘電体の誘電率により一意的に決まる値であり、あらかじめ誘電率毎のki、mi、b4、およびb5を計算し、データベースとして保存しておくと、処理を高速化することができる。
【0029】
以上の定式化により、実装距離d、パスコン実装数n、パスコン1つ当りの寄生誘導値La、および誘電率が決まれば、(6)式からf3を求めることができる。こうして得られたf2およびf3を用いれば、(2)および(3)式よりL2およびC2を求めることができ、得られたパラメータを用いて図5の回路を解くことにより、高い精度を保ったまま、高速にパスコンの効果(例えば、図3のようなインピーダンス特性)を計算することができる。
【0030】
ところで、上述した数式は、基板の電源−GND間の誘電体の厚みおよび基板の面積を一定と仮定した場合の数式である。L2は基板の誘電体の厚みに比例するので、例えば、図4の各パラメータを求めるときの基板厚を100μmとし、実際は200μmの基板に対するパラメータを取得したい場合は、得られたL2を2倍すればよい。また、基板の面積については、実装距離dを一定としたときの基板面積Sに対するf2およびf3をグラフ化し、それを近似式により近似すると、次のようになる。
【0031】
【数2】
【0032】
(8)式において、asiは、実装距離dによって変化するパラメータであり、bsおよびcは定数である。実装距離dに対するasiの変化をグラフ化した曲線を近似すると、f2とf3に応じてそれぞれ次のような近似式となる。
as2=af2d+bf2
as3=af3lnd+bf3
ここで、lnは自然対数であり、af2、bf2、af3、およびbf3は定数である。(8)式を用いてf2およびf3の変化量を求め、(2)および(3)式を計算するときに、この変化量を考慮したf2およびf3の値を使用すれば、基板面積が変化しても図2の各パラメータを求めることが可能となる。
【0033】
また、以上の定式化では受動素子の一例としてコンデンサを用いて説明したが、インダクタおよび抵抗のような受動素子の場合でも、コンデンサの場合と同様に上記の操作を行うことで、対応する近似式を導くことができる。
【0034】
次に、図6から図8までを参照しながら、上述した近似式により図4の各パラメータを求める装置および処理について説明する。
図6は、本実施形態のインピーダンス計算装置の構成図である。図6のインピーダンス計算装置は、例えば、コンピュータを用いて構成され、入力装置41、データ処理装置42、記憶装置43、および出力装置44を備える。データ処理装置42は、データ処理部45および演算処理部46を含む。
【0035】
記憶装置43は、あらかじめ計算しておいたパラメータ(a1、b1、b4、b5、m1〜m10、k5〜k10等)をデータベース(パラメータテーブル47)として記憶している。入力装置41は、プリント回路基板の物性情報および各素子の実装条件に関する情報を、解析条件として入力する。プリント回路基板の物性情報としては、基板の形状(電源層−GND層間の距離、基板面積等)や基板の誘電体の誘電率が入力される。
【0036】
データ処理装置42のデータ処理部45は、入力装置41から入力された解析条件から、パラメータテーブル47から取得すべきパラメータに関する条件(誘電率等)を抽出し、記憶装置43へ渡す。演算処理部46は、記憶装置43からデータ処理装置42へ渡されたテーブルデータを元に、上述した数式の演算を行い、演算結果を元にシミュレーションを実行する。そして、出力装置44は、データ処理部45および演算処理部46から、それぞれ、入力された解析条件およびシミュレーション結果を受け取って出力する。
【0037】
図7は、図6のインピーダンス計算装置によるインピーダンス計算処理のフローチャートであり、図8は、出力装置44のインタフェース画面の例を示している。
【0038】
まず、ユーザは、図8のインタフェース画面から基板の物性情報(PCBの条件55)を入力し(ステップS1)、受動素子の実装条件(電源の出力特性54およびパスコンの条件56)を入力する(ステップS2)。ステップS1では、基板の幅および長さ、電源層−GND層間の距離等が入力され、ステップS2では、基板上に実装されているDDコンバータおよびパスコンの等価回路パラメータ、パスコンの実装数、LSIの電源/GNDピンから高周波パスコンまでの距離等が入力される。
【0039】
次に、データ処理装置42は、簡易等価回路を求めるのに必要なパラメータをパラメータテーブル47から読み込み(ステップS3)、読み込んだパラメータを元に上述した数式の演算を行って、簡易等価回路を求める(ステップS4)。そして、トランジェント解析および/または周波数解析を行う。
【0040】
トランジェント解析は、パスコン実装による効果を時間領域で確認するシミュレーションである。このトランジェント解析では、まず、ユーザが、図8のインタフェース画面から簡易素子モデルの電源電圧、立上がり時間tr/立下り時間tf、ピーク値等のパラメータ(装置の条件53)を入力する(ステップS5)。この簡易素子モデルとしては、例えば、先願の特願2001−335416号に記載されたモデルを用いることができる。
【0041】
次に、演算処理部46が、SPICE(simulation program with integrated circuit emphasis )等の電気回路シミュレータにより、素子モデルを波源とした簡易等価回路の解析を行い(ステップS6)、解析結果である電源ノイズの時間波形のグラフ52をインタフェース画面に表示する(ステップS7)。
【0042】
周波数解析では、演算処理部46が、簡易等価回路全体の周波数に対するインピーダンス特性を、回路理論を元に計算し(ステップS8)、計算したインピーダンス特性のグラフ51をインタフェース画面に表示する(ステップS9)。
【0043】
図6のインピーダンス計算装置は、例えば、図9に示すような情報処理装置(コンピュータ)を用いて構成される。図9の情報処理装置は、CPU(中央処理装置)61、メモリ62、入力装置63、出力装置64、外部記憶装置65、媒体駆動装置66、およびネットワーク接続装置67を備え、それらはバス68により互いに接続されている。図6の入力装置41は、記憶装置43、および出力装置44を備える。
【0044】
メモリ62は、例えば、ROM(read only memory)、RAM(random access memory)等を含み、処理に用いられるプログラムとデータを格納する。CPU61は、メモリ62を利用してプログラムを実行することにより、必要な処理を行う。
【0045】
図6のデータ処理装置42は、CPU61およびメモリ62に対応し、図6のデータ処理部45および演算処理部46は、メモリ62に格納されたプログラムに対応する。
【0046】
入力装置63は、例えば、キーボード、ポインティングデバイス、タッチパネル等であり、図6の入力装置41に対応する。入力装置63は、ユーザからの指示や情報の入力に用いられる。
【0047】
出力装置64は、例えば、ディスプレイ装置およびスピーカを含み、図6の出力装置44に対応する。出力装置64は、ユーザへの問い合わせや処理結果の出力に用いられる。
【0048】
外部記憶装置65は、例えば、磁気ディスク装置、光ディスク装置、光磁気ディスク装置、テープ装置等である。情報処理装置は、この外部記憶装置65に、上述のプログラムとデータを保存しておき、必要に応じて、それらをメモリ62にロードして使用する。外部記憶装置65は、図6の記憶装置43としても用いられる。
【0049】
媒体駆動装置66は、可搬記録媒体69を駆動し、その記録内容にアクセスする。可搬記録媒体69としては、メモリカード、フレキシブルディスク、CD−ROM(compact disk read only memory )、光ディスク、光磁気ディスク等、任意のコンピュータ読み取り可能な記録媒体が用いられる。ユーザは、この可搬記録媒体69に上述のプログラムとデータを格納しておき、必要に応じて、それらをメモリ62にロードして使用する。
【0050】
ネットワーク接続装置67は、インターネット等の任意の通信ネットワークに接続され、通信に伴うデータ変換を行う。情報処理装置は、上述のプログラムとデータをネットワーク接続装置67を介して他の装置から受け取り、必要に応じて、それらをメモリ62にロードして使用する。
【0051】
図10は、図9の情報処理装置にプログラムとデータを供給することのできるコンピュータ読み取り可能な記録媒体を示している。可搬記録媒体69やサーバ70のデータベース71に保存されたプログラムとデータは、メモリ62にロードされる。このとき、サーバ70は、プログラムとデータを搬送する搬送信号を生成し、ネットワーク上の任意の伝送媒体を介して情報処理装置に送信する。そして、CPU61は、そのデータを用いてそのプログラムを実行し、必要な処理を行う。
【0052】
(付記1) ノイズ源と受動素子が実装された回路基板のインピーダンスを計算するインピーダンス計算装置であって、
前記回路基板上における前記ノイズ源と受動素子の間の実装距離を入力する入力手段と、
入力された実装距離をパラメータとして、前記ノイズ源と受動素子が実装された回路基板のモデルのインピーダンスを計算し、計算されたインピーダンスを用いて回路解析を行う解析手段と、
前記回路解析の結果を出力する出力手段と
を備えることを特徴とするインピーダンス計算装置。
【0053】
(付記2) ノイズ源と受動素子が実装された回路基板のインピーダンスを計算するインピーダンス計算方法であって、
前記ノイズ源と受動素子が実装された回路基板のモデルを生成し、
前記回路基板上における前記ノイズ源と受動素子の間の実装距離をパラメータとして、前記モデルのインピーダンスを計算する
ことを特徴とするインピーダンス計算方法。
【0054】
(付記3) ノイズ源と受動素子が実装された回路基板のインピーダンスを計算するコンピュータのためのプログラムであって、
前記回路基板上における前記ノイズ源と受動素子の間の実装距離を入力し、
入力された実装距離をパラメータとして、前記ノイズ源と受動素子が実装された回路基板のモデルのインピーダンスを計算する
処理を前記コンピュータに実行させることを特徴とするプログラム。
【0055】
(付記4) 前記受動素子は、コンデンサ、インダクタ、および抵抗のいずれかであることを特徴とする付記3記載のプログラム。
(付記5) 前記コンピュータは、前記ノイズ源と受動素子が実装された回路基板を該回路基板と受動素子の直並列回路とみなしたモデルに基づき、前記インピーダンスを計算することを特徴とする付記3記載のプログラム。
【0056】
(付記6) 前記プログラムは、前記回路基板の電源層とグランド層の間の誘電体の誘電率を入力する処理を前記コンピュータにさらに実行させ、前記コンピュータは、誘電率毎のパラメータ値を格納したパラメータテーブルから、入力された誘電率に応じたパラメータ値を読み出し、読み出したパラメータ値を用いて前記インピーダンスを計算することを特徴とする付記3記載のプログラム。
【0057】
(付記7) ノイズ源と受動素子が実装された回路基板のインピーダンスを計算するコンピュータのためのプログラムを記録した記録媒体であって、該プログラムは、
前記回路基板上における前記ノイズ源と受動素子の間の実装距離を入力し、
入力された実装距離をパラメータとして、前記ノイズ源と受動素子が実装された回路基板のモデルのインピーダンスを計算する
処理を前記コンピュータに実行させることを特徴とするコンピュータ読み取り可能なプログラム記録媒体。
【0058】
(付記8) ノイズ源と受動素子が実装された回路基板のインピーダンスを計算するコンピュータにプログラムを搬送する搬送信号であって、該プログラムは、
前記回路基板上における前記ノイズ源と受動素子の間の実装距離を入力し、
入力された実装距離をパラメータとして、前記ノイズ源と受動素子が実装された回路基板のモデルのインピーダンスを計算する
処理を前記コンピュータに実行させることを特徴とする搬送信号。
【0059】
【発明の効果】
本発明によれば、能動素子のようなノイズ源と受動素子とが実装された回路基板において、受動素子の実装距離に対するインピーダンスを高速かつ正確に計算することができる。したがって、このような回路基板上で発生する電源ノイズの解析に寄与するところが大きい。
【図面の簡単な説明】
【図1】本発明のインピーダンス計算装置の原理図である。
【図2】プリント回路基板に部品を実装した様子を示す図である。
【図3】パスコンの実装条件による基板インピーダンスの変化を示す図である。
【図4】回路部品が実装された基板の簡易モデルを示す図である。
【図5】プリント回路基板と高周波パスコンとの並列等価回路を示す図である。
【図6】インピーダンス計算装置の構成図である。
【図7】インピーダンス計算処理のフローチャートである。
【図8】インタフェース画面を示す図である。
【図9】情報処理装置の構成図である。
【図10】記録媒体を示す図である。
【符号の説明】
11 入力手段
12 解析手段
13 出力手段
21 DDコンバータ
22 低周波パスコン
23 LSI
24 高周波パスコン
31、32、33 曲線
41、63 入力装置
44、64 出力装置
42 データ処理装置
43 記憶装置
45 データ処理部
46 演算処理部
47 パラメータテーブル
51、52 グラフ
53 装置の条件
54 電源の出力特性
55 PCBの条件
56 パスコンの条件
61 CPU
62 メモリ
65 外部記憶装置
66 媒体駆動装置
67 ネットワーク接続装置
68 バス
69 可搬記録媒体
70 サーバ
71 データベース[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to an apparatus and a method for calculating the impedance of a circuit board on which a noise source such as an active element and a passive element are mounted.
[0002]
[Prior art]
When an element (component) mounted on a printed circuit board is driven, a current flows between a power supply layer and a GND (ground) layer of the board. When this current flows into another element, a voltage noise proportional to the impedance of the substrate is generated, which causes a malfunction of the other element. In order to accurately calculate this noise at the board design stage, it is necessary to accurately grasp the impedance of the board on which passive elements such as capacitors are mounted.
[0003]
Conventionally, methods used for modeling a substrate include a mesh division method and a finite element method. There is also a method of calculating impedance characteristics based on a simple impedance model including a decoupling capacitor viewed from an active element (for example, see Patent Document 1).
[0004]
[Patent Document 1]
JP, 2002-41594, A
[Problems to be solved by the invention]
However, the above-described conventional impedance calculation method has the following problems.
[0006]
The mesh division method and the finite element method have the drawback that the model shape becomes more complicated as accurate calculation is aimed, and the calculation takes an enormous amount of time. In addition, the number of currently used devices that operate at several hundred MHz is increasing, and in such a frequency range, the resonance point of the entire substrate (depending on the mounting distance from an active device to a passive device such as a capacitor) is considered. Resonance frequency) greatly changes. However, the calculation method of
[0007]
Therefore, in a circuit board including an element operating at a frequency of several hundred MHz or more, a method of calculating a change in impedance with respect to a mounting distance of a passive element at high speed and accurately is not known.
[0008]
An object of the present invention is to provide an apparatus and a method for quickly and accurately calculating an impedance with respect to a mounting distance between a noise source and a passive element on a circuit board on which a noise source such as an active element and a passive element are mounted. That is.
[0009]
[Means for Solving the Problems]
FIG. 1 is a diagram illustrating the principle of the impedance calculating apparatus according to the present invention. 1 includes an
[0010]
The input means 11 inputs the mounting distance between the noise source and the passive element on the circuit board. The
[0011]
As a model of a circuit board on which a noise source and a passive element are mounted, for example, a simple model in which the circuit board and the passive element are each regarded as a series circuit of RLC, and the circuit of the circuit board and the circuit of the passive element are regarded as a series-parallel circuit Is used. By equalizing the change in impedance with respect to the mounting distance using such a simple model, the impedance of the circuit board on which the noise source and the passive element are mounted can be calculated quickly and accurately.
[0012]
Further, if a circuit is analyzed by specifying a model based on the calculated impedance, it is possible to efficiently analyze a change in power supply noise with respect to a mounting distance of a passive element.
[0013]
The
[0014]
BEST MODE FOR CARRYING OUT THE INVENTION
Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.
In the present embodiment, the frequency characteristics of the circuit board at the time of mounting the passive element are noted, and the board is modeled by expressing this with a simple series-parallel circuit. By using such a model, the impedance of the mounting board and the power supply noise generated in the board can be calculated quickly and accurately.
[0015]
First, a method of modeling a substrate will be described with reference to FIGS.
FIG. 2 shows a state in which components (elements) including decaps (bypass capacitors) are mounted on a printed circuit board. Generally, a bypass capacitor is used to remove noise on a circuit.
[0016]
The DD (direct current-direct current)
[0017]
FIG. 3 shows how the substrate impedance changes when the distance (mounting distance) between the active element mounted on the printed circuit board and the decap is changed. The horizontal axis and the vertical axis in FIG. 3 represent frequency and substrate impedance, respectively, and
[0018]
Referring to FIG. 3, three resonance points (f 1 , f 2 , f 3 ) can be confirmed in any of the
(1) Power supply voltage value and equivalent impedance model of the DD converter that is the power supply source.
(2) Decap (low frequency)
Equivalent impedance model of low frequency decap. When there are a plurality of decaps, the parameters (capacitance, resistance, inductance) are calculated assuming that they are connected in parallel.
(3) Decap (high frequency)
Equivalent impedance model of high frequency decap. When a plurality of decaps are mounted equidistantly, the capacitance value and the resistance value are calculated as a parallel connection. The inductance value changes depending on the mounting distance from the element.
(4) PCB
Equivalent impedance model of the printed circuit board. Other than the resistance value, it changes depending on the distance between the power supply terminal (power supply / GND pin) of the element and the high-frequency decap.
(5) Load LSI
5 is an equivalent impedance model between a Vdd terminal and a GND terminal in an LSI element.
[0019]
Here, the parameters of the equivalent impedance model of the high-frequency decaps of (3) and the PCB of (4) vary depending on the number of mounted decaps and the distance from the high-frequency decap to the power supply / GND pin of the element (5). By formulating the amount of change, the model parameters in FIG. 4 can be determined according to the mounting conditions of the bypass capacitor. Therefore, it is possible to verify the change of the parameter with respect to the change of the decapacitor mounting condition by using the model of FIG. 4 having a short analysis time and high accuracy, and it is possible to more quickly and accurately obtain the optimal decaptor solution.
[0020]
The frequencies of the three resonance points in FIG. 3, and f 1,
[0021]
In addition, since the resistance value does not affect the resonance frequency, a parallel circuit of the printed circuit board and the high-frequency decap is expressed as shown in FIG. 5 using only the capacitance value and the inductance value, and a method for calculating these parameters is described. explain. 5, the capacitance value and the inductance value of the high frequency bypass capacitor, respectively and C 1 and L 1, a capacitance value and the inductance value of the printed circuit board, respectively and C 2 and L 2.
[0022]
Of the four parameters in FIG. 5, only the capacitance value C 1 of the high-frequency bypass capacitor does not change even by changing the mounting distance bypass capacitor. The first frequency f 1 in FIG. 3, since being determined by the C 1 and L 1, can be obtained L 1 from C 1 and f 1. The number of mounted bypass capacitor implementation distance, and by changing the bypass capacitor one per parasitic induction values, etc., to graph the relationship between the mounting distance d and L 1 from the power / GND pin of the device to a high frequency bypass capacitor, This can be expressed as an approximate expression as follows.
L 1 = a L1 log 10 d + b L1 (1)
Further, a L1 and b L1 can be represented by the following approximate expressions from the graph.
a L1 = a 1 n b1
b L1 = b 2 L a + b 3
Here, n is the number of bypass capacitor implementation, L a is a parasitic inductive value per one
b 2 = m 1 n m2
b 3 = m 3 nm 4
m 1 , m 2 , m 3 , and m 4 are values uniquely determined by the dielectric constant of the dielectric between the power supply layer and the GND layer of the substrate, and are previously determined as m 1 , m 2 , m 3 for each dielectric constant. , and m 4 calculates, when stored as a database, it is possible to speed up the process.
[0023]
The formulation described above, mounting the distance d, the bypass capacitor implementation the number n, the bypass capacitor one per parasitic inductive value L a, and if the dielectric constant is Kimare, can be obtained L 1 from (1).
[0024]
Next, C 2 and L 2 in FIG. 5 are uniquely determined from f 2 and f 3 . Relationship C 2 and L 2 and f 2 and f 3 are expressed as from FIG. 5 below.
[0025]
(Equation 1)
[0026]
Therefore, once the f 2 and f 3, L 1 and (2) and (3) can be obtained L 2 and C 2 from the equation. When graphing the change in f 2 with respect to the mounting distance d, it can be approximated by the following curve.
f 2 = a f2 d + b f2 (4)
Therefore, f 2 is obtained from the equation (4). Here, a f2 and b f2 are as follows by approximation from a graph.
a f2 = a 5 L a + a 6
b f2 = a 7 L a + a 8
Here, if the changes of a 5 , a 6 , a 7 , and a 8 with respect to the number n of decaps mounted are graphed, they can be approximated by the following curves.
a i = k i n + m i (i = 5,6,7,8) (5)
Therefore, a 5 , a 6 , a 7 , and a 8 can be obtained from equation (5). Here, k i and m i is a value uniquely determined by the dielectric constant of the dielectric of the power supply layer -GND layers of the substrate, to calculate the k i and m i for each advance dielectric constant, and stored as a database By doing so, the processing can be sped up.
[0027]
The formulation described above, mounting the distance d, the bypass capacitor implementation the number n, the bypass capacitor one per parasitic inductive value L a, and if the dielectric constant is Kimare, can be obtained f 2 from equation (4).
[0028]
Further, when graphed changes in f 3 with respect to the mounting distance d, can be approximated by the following curve.
f 3 = a f3 d + b f3 (6)
Therefore, f 2 is obtained from the equation (4). Here, a f3 and b f3 are as follows.
a f3 = a 9 L a + a 10
b f3 =
Here, when graphed changes to L a of a 9 and a 10, it can be approximated by the following curve.
a i = k inn + m i (i = 9,10) (7)
Therefore, a 9 and a 10 can be obtained from equation (7). Here, k i, m i, b 4, and b 5 is a value uniquely determined by the dielectric constant of the dielectric of the power supply layer -GND layers of the substrate, k i for each advance permittivity, m i, If b 4 and b 5 are calculated and stored as a database, the processing can be speeded up.
[0029]
The formulation described above, mounting the distance d, the bypass capacitor implementation the number n, the bypass capacitor one per parasitic inductive value L a, and if the dielectric constant is Kimare can be obtained f 3 from (6). The use of f 2 and f 3 thus obtained, by solving the circuit (2) and (3) than can be obtained L 2 and C 2 type, using the obtained
[0030]
By the way, the above formula is a formula on the assumption that the thickness of the dielectric between the power supply and GND of the substrate and the area of the substrate are constant. Since L 2 is proportional to the thickness of the dielectric material of the substrate, for example, the substrate thickness when obtaining each parameter in FIG. 4 is 100 μm, and when it is desired to actually obtain a parameter for a substrate of 200 μm, the obtained L 2 is 2 You only have to double it. Also, the area of the substrate, when graphed f 2 and f 3 to the substrate area S when the mounting distance d constant, is approximated by the approximate equation it becomes as follows.
[0031]
(Equation 2)
[0032]
(8) In the formula, a si is a parameter that varies with mounting distance d, b s and c are constants. When a change in a si with respect to the mounting distance d approximates the graphed curve, respectively the following approximate expression according to f 2 and f 3.
as2 = af2d + bf2
a s3 = a f3 lnd + b f3
Here, ln is a natural logarithm, and a f2 , b f2 , a f3 , and b f3 are constants. The change amounts of f 2 and f 3 are obtained by using the expression (8), and when calculating the expressions (2) and (3), the values of f 2 and f 3 in consideration of the change amounts are used. Each parameter in FIG. 2 can be obtained even if the substrate area changes.
[0033]
Also, in the above formulation, a capacitor was described as an example of a passive element. However, in the case of a passive element such as an inductor and a resistor, by performing the above operation in the same manner as in the case of a capacitor, a corresponding approximate expression is obtained. Can be led.
[0034]
Next, an apparatus and a process for obtaining each parameter in FIG. 4 by the above-described approximate expression will be described with reference to FIGS.
FIG. 6 is a configuration diagram of the impedance calculation device of the present embodiment. The impedance calculation device in FIG. 6 is configured using, for example, a computer, and includes an
[0035]
The
[0036]
The data processing unit 45 of the
[0037]
FIG. 7 is a flowchart of the impedance calculation process by the impedance calculation device of FIG. 6, and FIG. 8 shows an example of an interface screen of the
[0038]
First, the user inputs physical property information of the board (PCB condition 55) from the interface screen of FIG. 8 (step S1), and inputs mounting conditions of the passive element (output characteristic 54 of the power supply and
[0039]
Next, the
[0040]
The transient analysis is a simulation for confirming the effect of mounting a decap in a time domain. In this transient analysis, first, the user inputs parameters (device conditions 53) such as the power supply voltage, rise time tr / fall time tf, and peak value of the simple element model from the interface screen of FIG. 8 (step S5). . For example, a model described in Japanese Patent Application No. 2001-335416 can be used as the simplified element model.
[0041]
Next, the
[0042]
In the frequency analysis, the
[0043]
The impedance calculation device in FIG. 6 is configured using, for example, an information processing device (computer) as illustrated in FIG. 9 includes a CPU (Central Processing Unit) 61, a
[0044]
The
[0045]
The
[0046]
The
[0047]
The
[0048]
The
[0049]
The medium driving device 66 drives the
[0050]
The
[0051]
FIG. 10 illustrates a computer-readable recording medium that can supply a program and data to the information processing apparatus in FIG. 9. The programs and data stored in the
[0052]
(Supplementary Note 1) An impedance calculation device that calculates the impedance of a circuit board on which a noise source and a passive element are mounted,
Input means for inputting a mounting distance between the noise source and the passive element on the circuit board,
With the input mounting distance as a parameter, an analysis unit that calculates the impedance of the model of the circuit board on which the noise source and the passive element are mounted, and performs a circuit analysis using the calculated impedance,
Output means for outputting a result of the circuit analysis.
[0053]
(Supplementary Note 2) An impedance calculation method for calculating an impedance of a circuit board on which a noise source and a passive element are mounted,
Generate a model of the circuit board on which the noise source and passive elements are mounted,
An impedance calculation method, wherein the impedance of the model is calculated using a mounting distance between the noise source and the passive element on the circuit board as a parameter.
[0054]
(Supplementary Note 3) A computer program for calculating the impedance of a circuit board on which a noise source and a passive element are mounted,
Input the mounting distance between the noise source and the passive element on the circuit board,
A program causing the computer to execute a process of calculating an impedance of a model of a circuit board on which the noise source and the passive element are mounted, using the input mounting distance as a parameter.
[0055]
(Supplementary Note 4) The program according to
(Supplementary Note 5) The computer calculates the impedance based on a model in which the circuit board on which the noise source and the passive element are mounted is regarded as a series-parallel circuit of the circuit board and the passive element. The program described.
[0056]
(Supplementary Note 6) The program causes the computer to further execute a process of inputting a dielectric constant of a dielectric between a power supply layer and a ground layer of the circuit board, and the computer stores a parameter value for each dielectric constant. The program according to
[0057]
(Supplementary Note 7) A recording medium recording a program for a computer for calculating the impedance of a circuit board on which a noise source and a passive element are mounted, the program comprising:
Input the mounting distance between the noise source and the passive element on the circuit board,
A computer-readable program recording medium, which causes the computer to execute a process of calculating an impedance of a model of a circuit board on which the noise source and the passive element are mounted, using the input mounting distance as a parameter.
[0058]
(Supplementary Note 8) A carrier signal for carrying a program to a computer for calculating the impedance of a circuit board on which a noise source and a passive element are mounted, the program comprising:
Input the mounting distance between the noise source and the passive element on the circuit board,
A carrier signal, which causes the computer to execute a process of calculating an impedance of a model of a circuit board on which the noise source and the passive element are mounted, using the input mounting distance as a parameter.
[0059]
【The invention's effect】
According to the present invention, on a circuit board on which a noise source such as an active element and a passive element are mounted, the impedance with respect to the mounting distance of the passive element can be calculated quickly and accurately. Therefore, it greatly contributes to the analysis of power supply noise generated on such a circuit board.
[Brief description of the drawings]
FIG. 1 is a diagram illustrating the principle of an impedance calculation device according to the present invention.
FIG. 2 is a diagram showing a state where components are mounted on a printed circuit board.
FIG. 3 is a diagram showing a change in substrate impedance depending on mounting conditions of a decap.
FIG. 4 is a diagram showing a simplified model of a board on which circuit components are mounted.
FIG. 5 is a diagram showing a parallel equivalent circuit of a printed circuit board and a high-frequency decap.
FIG. 6 is a configuration diagram of an impedance calculation device.
FIG. 7 is a flowchart of an impedance calculation process.
FIG. 8 is a diagram showing an interface screen.
FIG. 9 is a configuration diagram of an information processing apparatus.
FIG. 10 is a diagram showing a recording medium.
[Explanation of symbols]
DESCRIPTION OF
24 High-
62
Claims (4)
前記回路基板上における前記ノイズ源と受動素子の間の実装距離を入力する入力手段と、
入力された実装距離をパラメータとして、前記ノイズ源と受動素子が実装された回路基板のモデルのインピーダンスを計算し、計算されたインピーダンスを用いて回路解析を行う解析手段と、
前記回路解析の結果を出力する出力手段と
を備えることを特徴とするインピーダンス計算装置。An impedance calculator for calculating the impedance of a circuit board on which a noise source and a passive element are mounted,
Input means for inputting a mounting distance between the noise source and the passive element on the circuit board,
With the input mounting distance as a parameter, an analysis unit that calculates the impedance of the model of the circuit board on which the noise source and the passive element are mounted, and performs a circuit analysis using the calculated impedance,
Output means for outputting a result of the circuit analysis.
前記ノイズ源と受動素子が実装された回路基板のモデルを生成し、
前記回路基板上における前記ノイズ源と受動素子の間の実装距離をパラメータとして、前記モデルのインピーダンスを計算する
ことを特徴とするインピーダンス計算方法。An impedance calculation method for calculating an impedance of a circuit board on which a noise source and a passive element are mounted,
Generate a model of the circuit board on which the noise source and passive elements are mounted,
An impedance calculation method, wherein the impedance of the model is calculated using a mounting distance between the noise source and the passive element on the circuit board as a parameter.
前記回路基板上における前記ノイズ源と受動素子の間の実装距離を入力し、
入力された実装距離をパラメータとして、前記ノイズ源と受動素子が実装された回路基板のモデルのインピーダンスを計算する
処理を前記コンピュータに実行させることを特徴とするプログラム。A program for a computer that calculates an impedance of a circuit board on which a noise source and a passive element are mounted,
Input the mounting distance between the noise source and the passive element on the circuit board,
A program causing the computer to execute a process of calculating an impedance of a model of a circuit board on which the noise source and the passive element are mounted, using the input mounting distance as a parameter.
前記回路基板上における前記ノイズ源と受動素子の間の実装距離を入力し、
入力された実装距離をパラメータとして、前記ノイズ源と受動素子が実装された回路基板のモデルのインピーダンスを計算する
処理を前記コンピュータに実行させることを特徴とするコンピュータ読み取り可能なプログラム記録媒体。A recording medium for recording a computer program for calculating the impedance of a circuit board on which a noise source and a passive element are mounted, the program comprising:
Input the mounting distance between the noise source and the passive element on the circuit board,
A computer-readable program recording medium, which causes the computer to execute a process of calculating an impedance of a model of a circuit board on which the noise source and the passive element are mounted, using the input mounting distance as a parameter.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003034243A JP4216088B2 (en) | 2003-02-12 | 2003-02-12 | Apparatus and method for calculating the impedance of a circuit board |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003034243A JP4216088B2 (en) | 2003-02-12 | 2003-02-12 | Apparatus and method for calculating the impedance of a circuit board |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004246503A true JP2004246503A (en) | 2004-09-02 |
JP4216088B2 JP4216088B2 (en) | 2009-01-28 |
Family
ID=33019992
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003034243A Expired - Fee Related JP4216088B2 (en) | 2003-02-12 | 2003-02-12 | Apparatus and method for calculating the impedance of a circuit board |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4216088B2 (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007293725A (en) * | 2006-04-26 | 2007-11-08 | Nec Corp | Resonance frequency calculation device and resonance frequency calculation method |
JP2009163449A (en) * | 2007-12-28 | 2009-07-23 | Fujitsu Ltd | Technique for creating analysis model and technique for creating circuit board model |
WO2013038511A1 (en) * | 2011-09-13 | 2013-03-21 | 富士通株式会社 | Design method for semiconductor integrated circuit and design program for semiconductor integrated circuit |
WO2022239275A1 (en) * | 2021-05-10 | 2022-11-17 | 日本たばこ産業株式会社 | Power supply unit for aerosol generation device |
-
2003
- 2003-02-12 JP JP2003034243A patent/JP4216088B2/en not_active Expired - Fee Related
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007293725A (en) * | 2006-04-26 | 2007-11-08 | Nec Corp | Resonance frequency calculation device and resonance frequency calculation method |
JP4697445B2 (en) * | 2006-04-26 | 2011-06-08 | 日本電気株式会社 | Resonance frequency calculation device and resonance frequency calculation method |
JP2009163449A (en) * | 2007-12-28 | 2009-07-23 | Fujitsu Ltd | Technique for creating analysis model and technique for creating circuit board model |
WO2013038511A1 (en) * | 2011-09-13 | 2013-03-21 | 富士通株式会社 | Design method for semiconductor integrated circuit and design program for semiconductor integrated circuit |
WO2022239275A1 (en) * | 2021-05-10 | 2022-11-17 | 日本たばこ産業株式会社 | Power supply unit for aerosol generation device |
Also Published As
Publication number | Publication date |
---|---|
JP4216088B2 (en) | 2009-01-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wu et al. | Overview of signal integrity and EMC design technologies on PCB: Fundamentals and latest progress | |
US8850375B2 (en) | Integrated circuit design and simulation | |
US10521533B2 (en) | Inductor simulation method and inductor nonlinear equivalent circuit model | |
JP5773101B2 (en) | Capacitor simulation method, capacitor simulation apparatus and method of using the same | |
US20070283299A1 (en) | A Method And Apparatus To Target Pre-Determined Spatially Varying Voltage Variation Across The Area Of The VLSI Power Distribution System Using Frequency Domain Analysis | |
JP2010204869A (en) | Circuit constant analysis method for equivalent circuit model of laminate chip inductor and circuit simulation method | |
US9507906B2 (en) | Metal interconnect modeling | |
de Jesús Leal-Romo et al. | Surrogate-based analysis and design optimization of power delivery networks | |
Kovacevic-Badstuebner et al. | Broadband circuit-oriented electromagnetic modeling for power electronics: 3-D PEEC solver vs. RLCG-solver | |
US20110301922A1 (en) | Equivalent circuit simulation system and method | |
JP4216088B2 (en) | Apparatus and method for calculating the impedance of a circuit board | |
JP4169755B2 (en) | Simulated measurement device for generated noise on electronic substrate and simulated measurement method for generated noise | |
JP5304460B2 (en) | Printed wiring board power circuit design apparatus, printed wiring board power circuit design method and program | |
Sanaullah et al. | A new real pole delay model for RLC interconnect using second order approximation | |
Schierholz et al. | Engineering-informed design space reduction for PCB based power delivery networks | |
CN112180179B (en) | Electromagnetic compatibility analysis method and device | |
JP4614094B2 (en) | Resonance frequency calculation device and resonance frequency calculation method | |
Knockaert et al. | Laguerre-based bandlimited reduced-order modeling | |
Liang et al. | Small-signal fractional-order model of PN junction long-base diode | |
Yin et al. | Switching power supply EMC optimisation design method | |
Li et al. | Development and validation of a microcontroller model for EMC | |
Archambeault et al. | Power and ground-reference plane impedance determination as decoupling capacitor distance increases | |
WO2023032786A1 (en) | Method for creating equivalent circuit model of multi-terminal capacitor, equivalent circuit model creation program, storage medium storing equivalent circuit model creation program, simulation method, and simulation device | |
US10621387B2 (en) | On-die decoupling capacitor area optimization | |
Erdin et al. | Placement of decoupling capacitors on power transmission lines |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060202 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20080701 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080708 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080903 |
|
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: 20081104 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20081105 |
|
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: 20111114 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111114 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121114 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121114 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131114 Year of fee payment: 5 |
|
LAPS | Cancellation because of no payment of annual fees |