JP4377437B2 - Substrate model creation device, and substrate noise analysis device and method - Google Patents

Substrate model creation device, and substrate noise analysis device and method Download PDF

Info

Publication number
JP4377437B2
JP4377437B2 JP2008059869A JP2008059869A JP4377437B2 JP 4377437 B2 JP4377437 B2 JP 4377437B2 JP 2008059869 A JP2008059869 A JP 2008059869A JP 2008059869 A JP2008059869 A JP 2008059869A JP 4377437 B2 JP4377437 B2 JP 4377437B2
Authority
JP
Japan
Prior art keywords
substrate
cell
impedance
cell region
model
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
Application number
JP2008059869A
Other languages
Japanese (ja)
Other versions
JP2008165830A (en
Inventor
淳平 野中
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NEC Electronics Corp
Original Assignee
NEC Electronics Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by NEC Electronics Corp filed Critical NEC Electronics Corp
Priority to JP2008059869A priority Critical patent/JP4377437B2/en
Publication of JP2008165830A publication Critical patent/JP2008165830A/en
Application granted granted Critical
Publication of JP4377437B2 publication Critical patent/JP4377437B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Design And Manufacture Of Integrated Circuits (AREA)

Description

本発明は、基板モデル作成装置及び方法、並びに基板ノイズ解析装置及び方法に関し、更に詳しくは、半導体集積回路の基板領域の寄生素子の影響を含めた回路シミュレーションを可能とする基板モデル作成装置及び方法、並びに基板ノイズ解析装置及び方法に関する。   The present invention relates to a substrate model creation apparatus and method, and a substrate noise analysis apparatus and method, and more particularly, a substrate model creation apparatus and method capable of circuit simulation including the influence of parasitic elements in a substrate region of a semiconductor integrated circuit. And a substrate noise analyzing apparatus and method.

半導体集積回路は、1つの半導体基板上に、複数の回路素子や配線が作りこまれ、所定の機能を有するブロックの集合として構成される。半導体集積回路では、半導体基板の回路パターンが形成されていない部分が寄生素子として働き、あるブロックで発生したノイズが、他のブロックに伝播し、このブロックの動作に影響を与えることが知られている。特に、アナログ回路とディジタル回路とが混載される半導体集積回路では、ディジタル回路ブロックで発生したノイズが、アナログ回路ブロックに伝わって影響を与え、動作の不具合を引き起こすことがある。このため、アナログ・ディジタル混載の半導体集積回路を設計する際には、半導体基板を伝達するノイズの影響を計算し、ノイズの影響で不具合が生じないように設計する必要がある。   A semiconductor integrated circuit is configured as a set of blocks having a predetermined function in which a plurality of circuit elements and wirings are formed on one semiconductor substrate. In a semiconductor integrated circuit, it is known that a portion of a semiconductor substrate where a circuit pattern is not formed functions as a parasitic element, and noise generated in one block propagates to another block and affects the operation of this block. Yes. In particular, in a semiconductor integrated circuit in which an analog circuit and a digital circuit are mixedly mounted, noise generated in the digital circuit block may be transmitted to the analog circuit block and affect it, thereby causing a malfunction. For this reason, when designing an analog / digital mixed semiconductor integrated circuit, it is necessary to calculate the influence of noise transmitted through the semiconductor substrate and to design such that no trouble is caused by the influence of noise.

半導体基板を伝播する基板ノイズの影響を計算するものとして、文献(Balshaz R. Stanisic, Nishath K. Verghese, RobA. Rutenbar, L.Richard Carleyand, David J. Allistot, "Addressing Substrate Coupling Mixed-Mode IC's: Simulation and Power Distribution Synthesis, "IEEE Journal of Solid-State Circuit, Vol.9, No.3, pp.226-238, March 1994)には、半導体集積回路の半導体基板を複数のメッシュに分割し、メッシュ分割された各セルを抵抗成分と静電容量成分とで模擬して得られた基板モデルを、“SPICE”に代表される回路シミュレータで分析する技術が記載されている。   The literature (Balshaz R. Stanisic, Nishath K. Verghese, RobA. Rutenbar, L. Richard Carleyand, David J. Allistot, "Addressing Substrate Coupling Mixed-Mode IC's: Simulation and Power Distribution Synthesis, "IEEE Journal of Solid-State Circuit, Vol.9, No.3, pp.226-238, March 1994), divides a semiconductor substrate of a semiconductor integrated circuit into a plurality of meshes. A technique is described in which a substrate model obtained by simulating each divided cell with a resistance component and a capacitance component is analyzed by a circuit simulator represented by “SPICE”.

図21は、上記した従来のシミュレーションに用いる半導体基板の基板部分をメッシュ状に分割した基板モデルを模式的に示している。また、図22は、図21のように分割された基板モデルの単位セルモデル201の詳細を示している。図21に示す例では、基板モデルは、x方向に5分割、y方向に4分割、z方向に3分割と、等分に分割され、合計60個の単位セルモデル201の集合として表現される。それぞれの単位セルモデル201は、図22に示すように、7つのノードN0〜N6と、6つの抵抗要素203とからなる。   FIG. 21 schematically shows a substrate model obtained by dividing the substrate portion of the semiconductor substrate used in the above-described conventional simulation into a mesh shape. FIG. 22 shows details of the unit cell model 201 of the substrate model divided as shown in FIG. In the example shown in FIG. 21, the board model is equally divided into five parts in the x direction, four parts in the y direction, and three parts in the z direction, and is expressed as a set of 60 unit cell models 201 in total. . Each unit cell model 201 includes seven nodes N0 to N6 and six resistance elements 203 as shown in FIG.

7つのノードは、単位セルモデルの中心のノードN0と、立方体形状の6つの表面のそれぞれの中心点のノードN1〜N6から構成される。6つの抵抗要素203は、ノードN0と、ノードN1〜N6のそれぞれとの間の等価的な抵抗成分及び容量成分が並列に接続された構造からなる。ここで、抵抗要素203は、基板の抵抗率と誘電率との関係から定まる誘電緩和時間が、回路の扱う信号の速度に比べて十分に短い場合には、その容量成分を省略して抵抗成分のみで近似することが可能である。   The seven nodes are composed of a node N0 at the center of the unit cell model and nodes N1 to N6 at the center points of the six cube-shaped surfaces. The six resistance elements 203 have a structure in which equivalent resistance components and capacitance components between the node N0 and each of the nodes N1 to N6 are connected in parallel. Here, when the dielectric relaxation time determined from the relationship between the resistivity and the dielectric constant of the substrate is sufficiently shorter than the signal speed handled by the circuit, the resistance element 203 omits the capacitance component and has a resistance component. It is possible to approximate with only.

抵抗要素203は、その抵抗成分のみを考えると、半導体基板の抵抗率をρとし、単位セルモデル201の各辺の長さをdx、dy、dzとすると、次式より算出される。
Rx=ρ・dx/(2・dy・dz)
Ry=ρ・dy/(2・dx・dz) 式(1)
Rz=ρ・dz/(2・dx・dy)
ただし、式(1)のRx,Ry,Rzは、それぞれ図21のx、y、z軸方向の抵抗成分を表し、RxはノードN0とノードN2、N4との間のそれぞれの抵抗要素203の抵抗値を、RyはノードN0とノードN3、N5との間の抵抗要素203の抵抗値を、RzはノードN0とノードN1、N6との間の抵抗要素203の抵抗値を表している。半導体基板は、ウェル、拡散層、又は、エピタキシャル層などによって構成されるため、場所によってその抵抗値が異なる。従って、単位セルモデル201は、その場所の抵抗値に応じて、抵抗要素203の抵抗値が計算される。
The resistance element 203 is calculated by the following equation when the resistivity of the semiconductor substrate is ρ and the length of each side of the unit cell model 201 is dx, dy, dz, considering only the resistance component.
Rx = ρ · dx / (2 · dy · dz)
Ry = ρ · dy / (2 · dx · dz) Equation (1)
Rz = ρ · dz / (2 · dx · dy)
However, Rx, Ry, and Rz in Expression (1) represent resistance components in the x-, y-, and z-axis directions in FIG. 21, respectively, and Rx is the resistance element 203 between the node N0 and the nodes N2 and N4. Ry represents the resistance value of the resistance element 203 between the node N0 and the nodes N3 and N5, and Rz represents the resistance value of the resistance element 203 between the node N0 and the nodes N1 and N6. Since the semiconductor substrate is composed of a well, a diffusion layer, an epitaxial layer, or the like, its resistance value varies depending on the location. Therefore, in the unit cell model 201, the resistance value of the resistance element 203 is calculated according to the resistance value at that location.

図23は、図21の半導体基板の等価回路を示している。半導体基板を抵抗成分からなる単位セルモデル201の集合として表現することで、図23に示すような抵抗成分がメッシュ接続された等価回路が得られる。この等価回路が、ディジタル回路ブロックと、アナログ回路ブロックとから構成される半導体装置の等価回路である場合には、ディジタル回路ブロック内のあるノードで発生したノイズが、アナログ回路ブロック内のノードへどのように伝播するかを、回路シミュレータを用いて計算するで、ディジタル回路のノイズが、アナログ回路に及ぼす影響を計算によって予測することが可能となる。   FIG. 23 shows an equivalent circuit of the semiconductor substrate of FIG. By expressing the semiconductor substrate as a set of unit cell models 201 composed of resistance components, an equivalent circuit in which the resistance components are mesh-connected as shown in FIG. 23 is obtained. When this equivalent circuit is an equivalent circuit of a semiconductor device composed of a digital circuit block and an analog circuit block, noise generated at a certain node in the digital circuit block is transferred to a node in the analog circuit block. Thus, it is possible to predict the influence of noise of the digital circuit on the analog circuit by calculation using the circuit simulator.

上記した基板ノイズの解析モデルは、精度の高い計算ができるが、ノード数が多いほど、シミュレーションにかかる時間が長くなるという欠点がある。この欠点を克服するものとして、特願2001−010359(インピーダンス網モデルの作成方法)では、基板モデル作成時の単位セルの大きさを局所的に変化させることにより、ノード数を削減する方法が提案されている。   The above-described analysis model for substrate noise can be calculated with high accuracy, but there is a drawback that the time required for simulation increases as the number of nodes increases. In order to overcome this drawback, Japanese Patent Application No. 2001-010359 (method for creating an impedance network model) proposes a method for reducing the number of nodes by locally changing the size of a unit cell when creating a substrate model. Has been.

図24は、シミュレーションの対象の半導体集積回路の断面を示している。半導体集積回路は、基板113上に、ウェル111、拡散層112が形成されている。例えば、上記先願の方法を使用して、同図に示す半導体集積回路の基板をモデル化する場合について考える。図25は、図24の半導体集積回路をモデル化した例であり、(a)はセル分割の様子を、(b)は等価回路を示している。同図(a)に示すように、上記先願では、ウェル111及び拡散層112を含む部分は、基板領域のセルよりも小さなセルでモデル化し、基板領域113の領域では、ウェル及び拡散層の領域よりも大きなセルでモデル化する。   FIG. 24 shows a cross section of a semiconductor integrated circuit to be simulated. In the semiconductor integrated circuit, a well 111 and a diffusion layer 112 are formed on a substrate 113. For example, consider the case where the substrate of the semiconductor integrated circuit shown in FIG. 25 shows an example in which the semiconductor integrated circuit of FIG. 24 is modeled. FIG. 25A shows a state of cell division, and FIG. 25B shows an equivalent circuit. As shown in FIG. 5A, in the above-mentioned prior application, the portion including the well 111 and the diffusion layer 112 is modeled by a cell smaller than the cell in the substrate region, and in the region of the substrate region 113, the well and diffusion layer are formed. Model with a cell larger than the region.

図25(b)は、同図(a)のようにメッシュ分割したセルのノード間を抵抗素子で表した等価回路を示している。この等価回路を、回路シミュレータ等で解析することで、あるノードが他のノードに対して与える影響を、計算によって求めることができる。この例では、細かな解析を要する箇所は小さなセルでモデル化し、細かな解析を要しない箇所は大きいセルでモデル化する。このため、解析するノード数を減少して、回路シミュレータの計算を簡易にしながらも、所望の箇所では細かな解析が可能になる。   FIG. 25B shows an equivalent circuit in which the nodes of the cells divided into meshes as shown in FIG. By analyzing this equivalent circuit with a circuit simulator or the like, the influence of a certain node on other nodes can be obtained by calculation. In this example, a part that requires detailed analysis is modeled by a small cell, and a part that does not require detailed analysis is modeled by a large cell. Therefore, the number of nodes to be analyzed is reduced and the calculation of the circuit simulator is simplified, but detailed analysis can be performed at a desired location.

ところで、上記先願では、半導体集積回路の箇所によって、モデル化するセルのサイズが異なるため、サイズの異なるセルが相互に接する部分においては、例えば図25(b)に示すノードAからノードBに至るルートを考えると、ノードAから拡散層内のノードを通過してBに至るルートと、ノードAからノードPを経由してBに至るルートができる。これは、サイズの大きなセルの1つに対して、そのセルに接するサイズの小さなセルが複数あるために、それらのセル間のノードが、同図(b)のノードPからノードA及びBなどに向かって斜め方向に接続されるためである。   By the way, in the above-mentioned prior application, since the size of the cell to be modeled differs depending on the location of the semiconductor integrated circuit, in the portion where the cells having different sizes contact each other, for example, from node A to node B shown in FIG. Considering the route to reach, there can be a route from node A through a node in the diffusion layer to B, and a route from node A to node B via node P. This is because, for one large cell, there are a plurality of small cells in contact with the cell, so the nodes between those cells are changed from node P to nodes A and B in FIG. This is because they are connected in an oblique direction toward.

メッシュ分割によっては、ノードAからノードPを経由してノードBに至るルートが、低インピーダンスな電流経路(以降このような電流経路をショートカットパスと呼ぶ)となり、ショートカットパスの影響が無視できなくなる。ショートカットパスは、主にウェル構造のシート抵抗に影響を与え、回路解析の精度を低下させる。従来のモデル化方法においては、このショートカットパスの影響を予測する手段は知られていなかった。   Depending on mesh division, the route from node A to node B via node P becomes a low-impedance current path (hereinafter, such a current path is referred to as a shortcut path), and the influence of the shortcut path cannot be ignored. The shortcut path mainly affects the sheet resistance of the well structure and reduces the accuracy of circuit analysis. In the conventional modeling method, means for predicting the influence of this shortcut path has not been known.

本発明は、半導体集積回路を複数のサイズを有するセルにメッシュ分割して得られた基板モデルを解析する際に、ショートカットパスが解析精度に与える影響を予測して、解析精度の低下が発生するか否かを判定し、解析の精度を高める基板モデル作成装置及び方法、並びに基板ノイズ解析装置及び方法を提供することを目的とする。   According to the present invention, when analyzing a substrate model obtained by dividing a semiconductor integrated circuit into cells having a plurality of sizes, the influence of the shortcut path on the analysis accuracy is predicted, and the analysis accuracy is reduced. It is an object of the present invention to provide a substrate model creation device and method, and a substrate noise analysis device and method for determining whether or not to improve analysis accuracy.

本発明の基板モデル作成装置は、半導体集積回路における基板部分を、複数種類のサイズを持つセルにメッシュ分割し、各セルを1個以上の抵抗素子から成る、又は、抵抗素子及び容量素子から成るインピーダンス要素に模擬することで、前記基板部分を等価回路でモデル化する基板モデル作成装置であって、第1のサイズを有する複数の第1セル領域と、該第1セル領域に隣接し前記第1のサイズよりも大きな第2のサイズを有する1つの第2セル領域との間の境界がpn接合を形成する際には、前記第1のセル領域の夫々と前記第2のセル領域との間の接続を、前記第1のセルのインピーダンス要素で定まるインピーダンスと、pn接合で定まる容量素子と、前記第2のセル領域のインピーダンス要素で定まるインピーダンスとの直列接続で表現すると共に、前記容量素子が第2セル領域に直接的に接続されるように素子の配置ないし移動を行うことを特徴とする。   In the substrate model creation apparatus of the present invention, a substrate portion in a semiconductor integrated circuit is mesh-divided into cells having a plurality of sizes, and each cell is composed of one or more resistance elements, or is composed of a resistance element and a capacitance element. A board model creation device for modeling the board portion with an equivalent circuit by simulating an impedance element, wherein a plurality of first cell areas having a first size and the first cell area adjacent to the first cell area When a boundary between one second cell region having a second size larger than one size forms a pn junction, each of the first cell regions and the second cell region A series connection of an impedance determined by an impedance element of the first cell, a capacitive element determined by a pn junction, and an impedance determined by an impedance element of the second cell region. Together they represent in, and performing arrangement or movement of the device so that the capacitive element is directly connected to a second cell region.

本発明の基板モデル作成方法は、データ処理装置を用い、半導体集積回路における基板部分を、複数種類のサイズを持つセルにメッシュ分割し、各セルを1個以上の抵抗素子から成る、又は、抵抗素子及び容量素子から成るインピーダンス要素に模擬することで、前記基板部分を等価回路でモデル化する基板モデル作成方法であって、前記データ処理装置が、第1のサイズを有する複数の第1セル領域と、該第1セル領域に隣接し前記第1のサイズよりも大きな第2のサイズを有する1つの第2セル領域との間の境界がpn接合を形成する際には、前記第1のセル領域の夫々と前記第2のセル領域との間の接続を、前記第1のセルのインピーダンス要素で定まるインピーダンスと、pn接合で定まる容量素子と、前記第2のセル領域のインピーダンス要素で定まるインピーダンスとの直列接続で表現すると共に、前記容量素子が第2セル領域に直接的に接続されるように素子の配置ないし移動を行うことを特徴とする。   The substrate model creation method of the present invention uses a data processing apparatus to divide a substrate portion in a semiconductor integrated circuit into cells having a plurality of types of sizes, and each cell is composed of one or more resistive elements, or a resistor. A substrate model generation method for modeling the substrate portion with an equivalent circuit by simulating an impedance element including an element and a capacitive element, wherein the data processing device includes a plurality of first cell regions having a first size. And a boundary between the second cell region adjacent to the first cell region and having a second size larger than the first size forms a pn junction, the first cell Connection between each of the regions and the second cell region is determined by impedance determined by the impedance element of the first cell, capacitive element determined by a pn junction, and impedance of the second cell region. Together they represent a series connection of an impedance determined by the dance element, and performing arrangement or movement of the device so that the capacitive element is directly connected to a second cell region.

本発明の基板モデル作成装置及び基板モデル作成方法では、異なるセルサイズにメッシュ分割したときに、複数の小さなサイズのセルと、その複数の小さなサイズのセルに接する、1つの大きなサイズのセルとがpn接合を形成する際には、小さなサイズのセルのインピーダンス要素と、大きなサイズのインピーダンス要素を小さなサイズのセルの数で並列に分割したインピーダンス要素とを、静電容量素子を介して、小さなサイズのセルと大きなサイズのセルとをそれぞれ接続したモデルを作成する。この場合、複数の小さなサイズのセルのインピーダンス要素を並列接続したものと、大きなサイズのセルの1つのインピーダンス要素とを、1つの静電容量素子を介して接続するのに比べて、ショートカットパスの影響が小さくなる。   In the board model creation device and the board model creation method of the present invention, when the mesh is divided into different cell sizes, a plurality of small sized cells and one large sized cell in contact with the plurality of small sized cells are provided. When forming a pn junction, an impedance element of a small-sized cell and an impedance element obtained by dividing a large-sized impedance element in parallel by the number of small-sized cells are connected via a capacitive element to a small size. A model is created by connecting the cells and large cells. In this case, the shortcut path has a short path compared to the case where the impedance elements of a plurality of small sized cells are connected in parallel and one impedance element of a large sized cell are connected via one capacitance element. The impact is reduced.

本発明の基板ノイズ解析装置は、基板モデル作成装置が作成した基板モデルに対し、基板上に作りこまれた機能回路のネットリストを結合して回路シミュレーションを行うことを特徴とする。   The board noise analyzing apparatus of the present invention is characterized in that a circuit simulation is performed by combining a net list of functional circuits built on a board with a board model created by the board model creating apparatus.

本発明の基板ノイズ解析方法は、データ処理装置を用いて回路シミュレーションを行い、基板ノイズ解析を行う方法であって、前記データ処理装置が、本発明の基板モデル作成方法で作成した基板モデルに対し、基板上に作りこまれた機能回路のネットリストを結合して回路シミュレーションを行うことを特徴とする。   The substrate noise analysis method of the present invention is a method of performing circuit noise analysis by performing circuit simulation using a data processing device, wherein the data processing device applies the substrate model created by the substrate model creation method of the present invention. The circuit simulation is performed by combining the net lists of the functional circuits built on the substrate.

本発明の基板ノイズ解析装置及び基板ノイズ解析方法では、本発明のいずれかの基板モデル作成装置又は基板モデル解析方法により作成された基板モデルを使用して、回路シミュレーションを行う。このため、精度の良い回路シミュレーションを簡易に実現する。   In the substrate noise analysis device and the substrate noise analysis method of the present invention, circuit simulation is performed using the substrate model created by any one of the substrate model creation device or the substrate model analysis method of the present invention. For this reason, an accurate circuit simulation is easily realized.

本発明の基板モデル作成装置では、前記判定手段が、所望の解析精度サイズを満足しないと判定すると、メッシュサイズを変更して再度メッシュ分割を行うことが好ましい。この場合、メッシュ分割の修正を再度行うことで、ショートカットパスの影響を小さくすることができる。   In the board model creation device of the present invention, it is preferable that if the determination unit determines that the desired analysis accuracy size is not satisfied, the mesh size is changed and mesh division is performed again. In this case, the effect of the shortcut path can be reduced by correcting the mesh division again.

また、本発明の基板モデル作成装置は、前記インピーダンスの比率と前記所定値との違いに依存してメッシュ分割サイズを修正することが好ましい。この場合、抵抗比はセルのサイズに依存して定まるため、メッシュ分割サイズを、所望の基板抵抗減少率となる抵抗比にあわせて設定する。このため、ショートカットパスの影響が小さい基板モデルを簡易に作成することができる。   Moreover, it is preferable that the board | substrate model creation apparatus of this invention corrects a mesh division | segmentation size depending on the difference between the ratio of the said impedance and the said predetermined value. In this case, since the resistance ratio is determined depending on the cell size, the mesh division size is set in accordance with the resistance ratio that provides a desired substrate resistance reduction rate. For this reason, it is possible to easily create a board model that is less affected by the shortcut path.

本発明の基板モデル作成装置では、前記メッシュ分割に際して分割の妥当性を判定する判定手段を備え、該判定手段は、前記第1セル領域と前記第2セル領域との間の境界と平行なインピーダンスを、前記第1のセル領域を表現する第1のインピーダンスとして算出し、前記第1のセル領域と前記第2のセル領域との間の接続を表現する第2のインピーダンス要素を算出し、前記第1のセル領域と前記第2のセル領域の間のpn接合容量を表現するインピーダンスを前記第2のインピーダンス要素から減算して第3のインピーダンス要素を算出し、前記第1のインピーダンスと前記第3のインピーダンスとの比率と、モデル精度に関連づけられた所定値とを比較することにより、基板モデルが所望の解析精度を満足するか否かの判定を行うことが好ましい。この場合、基板抵抗減少率が所望の値以下となる抵抗比の値をあらかじめ求めておき、メッシュ分割によって得られた抵抗比と、予め求めておいた抵抗比とを比べることで、ショートカットパスの影響を判断する。このため、得られた基板モデルの解析精度が所望の精度を満たすか否かを、簡易に判定することができる。   The substrate model creation apparatus of the present invention further includes a determination unit that determines the appropriateness of the division when the mesh is divided, and the determination unit includes an impedance parallel to a boundary between the first cell region and the second cell region. Is calculated as a first impedance representing the first cell region, a second impedance element representing a connection between the first cell region and the second cell region is calculated, A third impedance element is calculated by subtracting an impedance representing a pn junction capacitance between the first cell area and the second cell area from the second impedance element, and the first impedance and the first impedance are calculated. 3 is compared with a predetermined value associated with the model accuracy to determine whether or not the board model satisfies the desired analysis accuracy. It is preferred. In this case, a resistance ratio value at which the substrate resistance reduction rate is equal to or less than a desired value is obtained in advance, and the resistance ratio obtained by mesh division is compared with the resistance ratio obtained in advance, so that the shortcut path Determine the impact. For this reason, it is possible to easily determine whether or not the analysis accuracy of the obtained substrate model satisfies a desired accuracy.

本発明の基板モデル作成装置は、前記抵抗比と前記所定値との違いに依存して、サイズが異なるセル相互の接続部分を模擬するインピーダンス要素のインピーダンスを重み付けする補正を行うことが好ましい。
ショートカットパスの影響が大きいと判断されたときには、ショートカットパスを構成するインピーダンス要素に重み付けをすることもできる。この場合、重み付け後の基板モデルでは、同じメッシュ分割であっても、重み付け前と比較して、基板抵抗減少率は低い値となる。このため、メッシュ分割を再度行わなくても、ショートカットパスの影響を低くすることができる。
The substrate model creation apparatus of the present invention preferably performs correction for weighting the impedance of an impedance element that simulates a connection portion between cells having different sizes, depending on the difference between the resistance ratio and the predetermined value.
When it is determined that the influence of the shortcut path is large, the impedance elements constituting the shortcut path can be weighted. In this case, in the weighted board model, even if the same mesh division is used, the board resistance reduction rate is lower than that before weighting. For this reason, it is possible to reduce the influence of the shortcut path without performing mesh division again.

以上説明したように、本発明の基板モデル作成装置及び方法、並びに基板ノイズ解析装置及び方法は、異なるセルサイズに半導体集積回路の基板部分をメッシュ分割して基板モデルを作成する場合でも、抵抗比を計算し、基板抵抗減少率を求めることで、簡易にショートカットパスの影響が大きいか否かを判断できる。ショートカットパスの影響が大きいと判断されたときには、メッシュサイズの変更や、重み付け計算を行うことで、ショートカットパスの影響が所望の値以下となる基板モデルを作成する。このため、基板ノイズの解析を精度よく行うことができる。   As described above, the substrate model creation apparatus and method, and the substrate noise analysis apparatus and method according to the present invention provide a resistance ratio even when creating a substrate model by dividing the substrate portion of the semiconductor integrated circuit into different cell sizes. And calculating the substrate resistance reduction rate, it is possible to easily determine whether or not the influence of the shortcut path is large. When it is determined that the influence of the shortcut path is large, a board model in which the influence of the shortcut path is equal to or less than a desired value is created by changing the mesh size or performing weighting calculation. For this reason, it is possible to accurately analyze the substrate noise.

以下、図面を参照し、本発明の実施形態例に基づいて、本発明を更に詳細に説明する。   Hereinafter, with reference to the drawings, the present invention will be described in more detail based on exemplary embodiments of the present invention.

(第1の実施形態例)
本発明の第1の実施形態例の基板モデル作成方法及び基板ノイズ解析方法の説明に先立って、本実施形態例の基板モデル作成方法及び基板ノイズ解析方法がどのような手法により、ショートカットパスの影響を判断するのかについて説明する。図1は、半導体集積回路の基板モデルであり、(a)は基板モデルのセル分割の様子を、(b)、(c)、(d)はそれぞれセル分割における等価回路を示している。便宜上、l方向を横方向、w方向を奥行き方向、h方向を高さ方向と呼ぶ。また、直方体の端面のh方向の面を上面及び下面と呼ぶ。
(First embodiment)
Prior to the description of the substrate model creation method and substrate noise analysis method according to the first embodiment of the present invention, the method of the substrate model creation method and substrate noise analysis method according to the present embodiment determines the effect of the shortcut path. Will be described. FIG. 1 shows a substrate model of a semiconductor integrated circuit, where (a) shows a state of cell division of the substrate model, and (b), (c), and (d) show equivalent circuits in cell division, respectively. For convenience, the l direction is referred to as the horizontal direction, the w direction is referred to as the depth direction, and the h direction is referred to as the height direction. Also, the h-direction surfaces of the end faces of the rectangular parallelepiped are referred to as an upper surface and a lower surface.

半導体基板のモデルは、例えば、抵抗率がρ1のウェル領域と、抵抗率がρ2の基板領域を模しているものとする。基板モデルのセル分割においては、同図(a)に示すように、比較的微細な構造を含むウェル領域は、それぞれが(w×h1×l/5)のサイズを有する、5つの小さなセル101に分割され、基板領域は、(w×h2×l)のサイズを有する大きなセル102に分割される。セル101は、紙面横方向に直列に接続され、セル102の上面は、5つのセル101の下面とそれぞれ接している。 It is assumed that the model of the semiconductor substrate simulates, for example, a well region having a resistivity ρ 1 and a substrate region having a resistivity ρ 2 . In the cell division of the substrate model, as shown in FIG. 5A, the well region including a relatively fine structure has five small cells each having a size of (w × h 1 × 1/5). The substrate area is divided into large cells 102 having a size of (w × h 2 × l). The cells 101 are connected in series in the horizontal direction of the paper, and the upper surface of the cell 102 is in contact with the lower surfaces of the five cells 101.

図1(a)のようにセル分割された半導体集積回路の基板モデルを、等価回路に変換すると、同図(b)のような等価回路が得られる。セル101及び102は、それぞれが図21に示す単位セルモデル201と同様の構成を持つ。各セル101の中心のノード間を接続する横方向の抵抗値R1は、次式によって計算される。

Figure 0004377437
When the substrate model of the semiconductor integrated circuit divided into cells as shown in FIG. 1A is converted into an equivalent circuit, an equivalent circuit as shown in FIG. 1B is obtained. Each of the cells 101 and 102 has the same configuration as the unit cell model 201 shown in FIG. The lateral resistance value R 1 connecting the nodes at the center of each cell 101 is calculated by the following equation.
Figure 0004377437

縦方向の、セル102の中心のノードPとそれぞれのセル101の中心のノードを接続する斜め方向の抵抗R2は、セル102の領域の成分と、それぞれのセル101の領域の成分とに分けて考えるとよい。このとき、図1(c)に示すように、セル101の抵抗成分R2aを並列に接続し、これとセル102の抵抗成分R2bとを直列に接続するモデルを考えると、端子Aから端子Bに至る電流経路に、ノードPを経由しない経路ができる。この経路は、2つの抵抗R2aを介して端子Aと端子Bを接続するため、端子AB間を低インピーダンスな電流経路で接続するショートカットパスとなる。このような等価回路を用いると、ショートカットパスの影響が大きくなり、解析精度が低下する。そこで、同図(d)に示すように、セル102の抵抗成分R2bを複数にわけ、それらを、それぞれの抵抗成分R2aに直列に接続するモデルを考える。このモデルでは、端子Aから端子Bに至る経路に、ノードPを経由しない経路が存在しないため、同図(c)示すモデルよりも、ショートカットパスの影響が小さくなる。 The resistance R 2 in the diagonal direction connecting the node P at the center of the cell 102 and the node at the center of each cell 101 in the vertical direction is divided into a component of the region of the cell 102 and a component of the region of the cell 101. It is good to think about it. At this time, as shown in FIG. 1C, considering a model in which the resistance component R 2a of the cell 101 is connected in parallel and the resistance component R 2b of the cell 102 is connected in series, the terminal A to the terminal A path that does not pass through the node P can be formed in the current path reaching B. Since this path connects the terminal A and the terminal B via the two resistors R 2a , it becomes a shortcut path that connects the terminals AB with a low-impedance current path. When such an equivalent circuit is used, the influence of the shortcut path is increased, and the analysis accuracy is lowered. Therefore, as shown in FIG. 6D, a model is considered in which the resistance component R 2b of the cell 102 is divided into a plurality and connected in series to each resistance component R 2a . In this model, since there is no route that does not pass through the node P in the route from the terminal A to the terminal B, the influence of the shortcut path is smaller than that of the model shown in FIG.

図1(b)の抵抗R2は、セル102の領域の成分R2bと、セル101の抵抗成分R2aの和となり、次式により計算される。

Figure 0004377437
ただし、nは、セル101の個数であり、図1の例では、n=5である。このとき、ノードPに対して、図1(b)に示す等価回路以外に外部からの電流の出入りが無ければ、ショートカットパスが存在する場合の端子A−B間の合成抵抗RSは次式によって計算される。
Figure 0004377437
ここで、γは、
Figure 0004377437
である。なお、式(1.3)の算出については、後で説明する。 The resistance R 2 in FIG. 1B is the sum of the component R 2b in the region of the cell 102 and the resistance component R 2a of the cell 101, and is calculated by the following equation.
Figure 0004377437
However, n is the number of the cells 101, and n = 5 in the example of FIG. At this time, if there is no external current in / out of the node P other than the equivalent circuit shown in FIG. 1B, the combined resistance R S between the terminals A and B when the shortcut path exists is Is calculated by
Figure 0004377437
Where γ is
Figure 0004377437
It is. The calculation of equation (1.3) will be described later.

一方、ショートカットパスが無い場合には、端子Aから端子Bに至る経路にはノードPを経由する経路がないため、端子AB間の合成抵抗RABは、

Figure 0004377437
となる。ここで、ショートカットパスがある場合とない場合との、端子AB間の合成抵抗の比をとり、ショートカットパスによる基板抵抗の減少率をEとすると、基板抵抗減少率Eは次式により計算される。
Figure 0004377437
式(1.6)からショートカットパスがある場合の抵抗値RSと、ない場合の抵抗値RABとの比が1に近ければ、基板抵抗減少率Eが十分に小さいとなり、ショートカットパスの影響が小さいと言える。 On the other hand, when there is no shortcut path, there is no path via the node P in the path from the terminal A to the terminal B, so the combined resistance R AB between the terminals AB is
Figure 0004377437
It becomes. Here, when the ratio of the combined resistance between the terminals AB with and without the shortcut path is taken and the reduction rate of the substrate resistance due to the shortcut path is E, the substrate resistance reduction rate E is calculated by the following equation. .
Figure 0004377437
If the ratio of the resistance value R S when there is a shortcut path and the resistance value R AB when there is no shortcut path is close to 1 from Equation (1.6), the substrate resistance reduction rate E becomes sufficiently small, and the influence of the shortcut path Can be said to be small.

図2は、図1(b)の等価回路のうち、ノードPに接続された抵抗R2のみを示している。ショートカットパスの影響が大きくなるのは、端子AB間の抵抗値に比べて、ショートカットパスを構成する領域である、図2に点線で示される面Cと面Dとの間の領域の抵抗値が低い場合であると予想される。面CD間の抵抗値は、

Figure 0004377437
で計算される。また、ショートカットパスが無い場合の端子A−B間の抵抗RABは前述のように式(1.5)で求められる。式(1.5)をR1について書き直し、また、式(1.2)と式(1.7)の関係から次式が得られる。
1=RAB/n (1.8)
2=n・RCD (1.9)
つまり、RCDが小さくなるとR2も小さくなるために、ショートカットパスの影響が大きくなる。 FIG. 2 shows only the resistor R 2 connected to the node P in the equivalent circuit of FIG. The influence of the shortcut path is greater than the resistance value between the terminals AB. The resistance value of the area between the plane C and the plane D shown in FIG. Expected to be low. The resistance value between the surface CDs is
Figure 0004377437
Calculated by Further, the resistance R AB between the terminals A and B when there is no shortcut path is obtained by the equation (1.5) as described above. Equation (1.5) is rewritten for R 1 , and the following equation is obtained from the relationship between Equation (1.2) and Equation (1.7).
R 1 = R AB / n (1.8)
R 2 = n · R CD (1.9)
That is, as R CD becomes smaller, R 2 also becomes smaller, so the influence of the shortcut path becomes larger.

ここで、式(1.3)と(1.8)とを式(1.6)に代入して、ショートカットパスによる基板抵抗減少率Eを求めると、

Figure 0004377437
となる。また、このとき抵抗RCDとRABとの比を抵抗比xとおくと、式(1.5)と式(1.7)から
Figure 0004377437
となり、抵抗比xはセルのサイズ及び抵抗率ρ1、ρ2に依存することがわかる。γは、抵抗比xを用いて、
Figure 0004377437
と表せる。式(1.10)及び式(1.12)から、基板抵抗減少率Eは、抵抗比xとセル数nのみに依存することがわかる。このため、ショートカットパスの影響は、抵抗比xとセル数nのみで算出できると考えてよい。 Here, substituting Equations (1.3) and (1.8) into Equation (1.6) to determine the substrate resistance reduction rate E by the shortcut path,
Figure 0004377437
It becomes. At this time, when the ratio of the resistances R CD and R AB is set to the resistance ratio x, from the formulas (1.5) and (1.7),
Figure 0004377437
It can be seen that the resistance ratio x depends on the cell size and the resistivity ρ 1 , ρ 2 . γ is the resistance ratio x,
Figure 0004377437
It can be expressed. From the expressions (1.10) and (1.12), it can be seen that the substrate resistance reduction rate E depends only on the resistance ratio x and the number of cells n. For this reason, it may be considered that the influence of the shortcut path can be calculated only by the resistance ratio x and the number of cells n.

図3は、基板抵抗減少率Eと抵抗比x、セル数nとの関係を求めた結果をグラフとして示している。同図に示すように、基板抵抗減少率Eは、セル数nが変化してもほとんど影響を受けない。このため、基板抵抗減少率Eは、ほぼ抵抗比xのみによって算出できると言える。セル分割に当たって、抵抗比xは、前述のように、セルのサイズ及び抵抗率ρ1、ρ2に依存するが、抵抗率ρ1、ρ2は解析対象の領域に固有なものであるため、セルのサイズのみに依存すると言える。図4は、図3のグラフから、基板抵抗の減少率が、それぞれ5%,10%,20%,30%以下になる抵抗比xの最小値xminを求め、これを表として示している。例えば、ショートカットパスによる基板抵抗減少率Eを10%以下にしたい場合には、図4を参照するとこのときのxminは0.74であるため、
x>xmin (1.13)
となる抵抗比xを有するセルに分割すればよい。
FIG. 3 is a graph showing the results of determining the relationship between the substrate resistance reduction rate E, the resistance ratio x, and the number of cells n. As shown in the figure, the substrate resistance reduction rate E is hardly affected even when the number of cells n changes. For this reason, it can be said that the substrate resistance reduction rate E can be calculated almost only by the resistance ratio x. In the cell division, the resistance ratio x depends on the cell size and the resistivity ρ 1 , ρ 2 as described above, but the resistivity ρ 1 , ρ 2 is specific to the region to be analyzed. It can be said that it depends only on the cell size. FIG. 4 shows, as a table, the minimum value x min of the resistance ratio x at which the reduction rate of the substrate resistance is 5%, 10%, 20%, and 30% or less, respectively, from the graph of FIG. . For example, when it is desired to reduce the substrate resistance reduction rate E by the shortcut path to 10% or less, referring to FIG. 4, x min at this time is 0.74.
x> x min (1.13)
What is necessary is just to divide | segment into the cell which has the resistance ratio x used as follows.

なお、式(1.5)及び式(1.7)では、RAB及びRCDは、容量成分を省略し、その抵抗成分のみで近似しているが、図21に示す単位セルモデル201で説明したように、本来抵抗要素は、抵抗素子と容量素子とが並列に接続された構造である。回路が扱う信号が、半導体基板の誘電緩和時間よりも高速な場合には、上記容量素子の影響が顕著になる。このため、RABとRCDのそれぞれに上記容量素子を並列接続させた構造のインピーダンスZABとZCDとを求め、ZABとZCDとのインピーダンス比xからショートカットパスの影響を求めることで、より高精度に基板抵抗減少率、つまり、解析誤差を求めることが可能になる。 In equations (1.5) and (1.7), R AB and R CD are approximated only by their resistance components, omitting the capacitance component, but in unit cell model 201 shown in FIG. As described above, the resistive element originally has a structure in which a resistive element and a capacitive element are connected in parallel. When the signal handled by the circuit is faster than the dielectric relaxation time of the semiconductor substrate, the influence of the capacitor element becomes significant. Therefore, the impedance Z AB and Z CD having a structure in which the capacitive element is connected in parallel to each of R AB and R CD are obtained, and the influence of the shortcut path is obtained from the impedance ratio x between Z AB and Z CD. Thus, the substrate resistance reduction rate, that is, the analysis error can be obtained with higher accuracy.

図5は、本発明の第1の実施形態例の基板ノイズ解析装置を、その動作手順を示すフローチャートを兼用して示している。以下、同図を参照して、本発明の第1の実施形態例の基板モデル作成方法、基板モデル作成装置、基板ノイズ解析方法、及び、基板ノイズ解析装置について説明する。本実施形態例の基板モデル作成装置及び基板ノイズ解析装置は、少なくとも、キーボード等の入力装置11と、プログラム制御により動作するデータ処理装置12と、情報を記憶する記憶装置13と、ディスプレイ装置や印刷装置等の出力装置14とを備えるコンピュータシステムにより構成される。   FIG. 5 shows a substrate noise analyzing apparatus according to the first embodiment of the present invention by using a flowchart showing an operation procedure thereof. Hereinafter, a substrate model creation method, a substrate model creation device, a substrate noise analysis method, and a substrate noise analysis device according to a first embodiment of the present invention will be described with reference to FIG. The board model creation apparatus and board noise analysis apparatus according to this embodiment include at least an input device 11 such as a keyboard, a data processing device 12 that operates under program control, a storage device 13 that stores information, a display device, and a printing device. And a computer system including an output device 14 such as a device.

記憶装置13は、レイアウトデータ記憶部15と、デバイスパラメータ記憶部16と、基板モデルネットリスト記憶部17とを備え、データ処理装置12にデータを入出力する。データ処理装置12は、メッシュ分割手段18と、メッシュ分割の妥当性判定手段19と、メッシュ分割修正手段20と、ネットリスト変換手段21と、回路シミュレーション手段22とを備え、図示しないソフトウェアプログラムに基づいて動作する。入力装置11は、データ処理装置12に解析者等が与えるべきコマンド等を入力する。出力装置は、シミュレーションによって得られた結果などを出力する。   The storage device 13 includes a layout data storage unit 15, a device parameter storage unit 16, and a board model netlist storage unit 17, and inputs / outputs data to / from the data processing device 12. The data processing device 12 includes a mesh division unit 18, a mesh division validity determination unit 19, a mesh division correction unit 20, a netlist conversion unit 21, and a circuit simulation unit 22, and is based on a software program (not shown). Works. The input device 11 inputs a command or the like to be given by an analyst or the like to the data processing device 12. The output device outputs a result obtained by the simulation.

メッシュ分割手段18は、レイアウトデータ記憶部15に記憶されている、解析対象の半導体集積回路のレイアウトデータと、デバイスパラメータ記憶部16に記憶されている、半導体基板各部の抵抗率や容量パラメータなどのデバイスパラメータとを入力し、それらをもとに基板をメッシュ分割する。メッシュ分割の妥当性判定手段19は、上述の手法に基づいて、ショートカットパスが解析精度低下の原因となるか否かを判定し、メッシュ分割パターンが妥当なものか否かの判定を行う。   The mesh dividing means 18 includes the layout data of the semiconductor integrated circuit to be analyzed stored in the layout data storage unit 15 and the resistivity and capacitance parameters of each part of the semiconductor substrate stored in the device parameter storage unit 16. The device parameters are input, and the substrate is meshed based on them. Based on the above-described method, the mesh division validity determination unit 19 determines whether or not the shortcut path causes a decrease in analysis accuracy, and determines whether or not the mesh division pattern is appropriate.

メッシュ分割修正手段20は、メッシュ分割パターンが妥当でないと判定された場合には、メッシュ分割パターンを修正する。ネットリスト変換手段21は、メッシュ分割された基板を、基板上に構成される回路のネットリストへと変換し、メッシュ分割パターンとネットリストとを、基板モデルネットリスト記憶部17に出力する。シミュレーション手段22は、生成された基板モデルを用いて回路シミュレーションを行い、基板ノイズを解析する。   When it is determined that the mesh division pattern is not valid, the mesh division correction unit 20 corrects the mesh division pattern. The netlist conversion means 21 converts the mesh-divided board into a netlist of a circuit configured on the board, and outputs the mesh division pattern and netlist to the board model netlist storage unit 17. The simulation means 22 performs circuit simulation using the generated board model and analyzes board noise.

次に、本実施形態例の基板モデル作成方法及び基板ノイズ解析方法について説明する。メッシュ分割手段18は、レイアウトデータとデバイスパラメータとを入力し、基板のメッシュ分割パターンを得る(ステップS1)。ステップS1は、コンピュータシステムが自動的に行ってもよいし、解析者が手動で行ってもよい。メッシュ分割の妥当性判定手段19は、ステップS1で得られたメッシュ分割パターンの抵抗比xと、要求される解析精度によって、例えば図4の関係で、定まるxminとを比較して、式(1.13)の条件を満たすか否かを判断する(ステップS2)。 Next, a substrate model creation method and a substrate noise analysis method according to this embodiment will be described. The mesh dividing means 18 inputs layout data and device parameters, and obtains a board mesh dividing pattern (step S1). Step S1 may be performed automatically by a computer system or manually by an analyst. The mesh division validity determination means 19 compares the resistance ratio x of the mesh division pattern obtained in step S1 with x min determined in the relationship shown in FIG. It is determined whether or not the condition of 1.13) is satisfied (step S2).

条件を満たさないと判断した場合、つまり抵抗比xがxmin未満のため、ショートカットパスによる基板抵抗の減少が無視できなくなり、解析精度が低下すると判断した場合には、メッシュ分割修正手段20は、メッシュ分割パターンの修正を行い(ステップS3)、ステップS2に戻って、修正したメッシュ分割パターンで判断をやり直す。ステップS2で、条件を満たすと判断されると、ネットリスト変換手段21は、メッシュ分割パターンを、ネットリストに変換し、メッシュ分割パターンとネットリストとを、基板モデルとして基板モデルネットリスト記憶部17に出力する(ステップS4)。 If it is determined that the condition is not satisfied, that is, the resistance ratio x is less than x min , the decrease in the substrate resistance due to the shortcut path cannot be ignored, and if it is determined that the analysis accuracy is reduced, the mesh division correcting unit 20 The mesh division pattern is corrected (step S3), the process returns to step S2, and the determination is performed again using the corrected mesh division pattern. If it is determined in step S2 that the condition is satisfied, the net list conversion unit 21 converts the mesh division pattern into a net list, and uses the mesh division pattern and the net list as a board model as a board model net list storage unit 17. (Step S4).

ステップS1からステップS4までの処理によって、基板モデルが作成され、作成された基板モデルを用いて基板ノイズの解析が可能になる。シミュレーション手段22は、任意のときに、基板モデルネットリスト記憶部17からメッシュパターン及びネットリストを読み込んで、メッシュパターンに、基板上に構成される回路のネットリストを結合して回路シミュレーションを行い、基板ノイズが回路の動作に与える影響を求める(ステップS5)。   A substrate model is created by the processing from step S1 to step S4, and the substrate noise can be analyzed using the created substrate model. The simulation means 22 reads the mesh pattern and the net list from the board model net list storage unit 17 at any time, and performs a circuit simulation by combining the mesh pattern with the net list of the circuit configured on the board. The influence of substrate noise on circuit operation is determined (step S5).

ステップS3でのメッシュ分割の修正は、図1のメッシュ分割で考えると、前述の式(1.11)で示したように抵抗比xはセルのサイズに依存するため、l又はh2を調節することで、抵抗比xをxmin未満の値となるようにする。例えばh2を調節する場合で考えると、hは、式(1.11)のxをxminに置き換えてh2について解くと、

Figure 0004377437
となる。また、lを調節する場合を考えると、同様にして、lについて解くと、
Figure 0004377437
となる。l又はh2をこのように調節することで、基板抵抗減少率を所望の値以下にするメッシュ分割が得られる。 In the modification of the mesh division in step S3, when considering the mesh division in FIG. 1, the resistance ratio x depends on the cell size as shown in the above formula (1.11), so l or h 2 is adjusted. As a result, the resistance ratio x is set to a value less than x min . For example, consider the case of adjusting the h 2, h 2 is solved for h 2 by replacing the x in equation (1.11) to x min,
Figure 0004377437
It becomes. Also, considering the case of adjusting l, similarly, if l is solved,
Figure 0004377437
It becomes. By adjusting l or h 2 in this way, mesh division can be obtained in which the substrate resistance reduction rate is less than or equal to a desired value.

ステップS4では、多端子F行列(村坂 住隆,永田 真,森江 隆,岩田 穆,“F行列を用いたチップレベルの基板雑音解析法” 信学技報ICD99-147,1999年9月)や、AWE法(Balshaz R. Stanisic, Nishath K. Verghese,Rob A. Rutenbar, L.Richard Carleyand, David J. Allistot, "AddressingSubstrate Coupling Mixed-Mode IC's: Simulation and Power DistributionSynthesis," IEEE Journal of Solid-State Circuit, Vol.9, No.3, pp.226-238,March 1994)等のアルゴリズムを用いて、ネットリストの縮退を行うことができる。ステップS3では、メッシュ分割の修正を行って、そのメッシュサイズが変化するが、メッシュ構造自体が変わるものではないので、ステップS4でネットリストを縮退するアルゴリズムを適用することができる。   In step S4, a multi-terminal F matrix (Sumitaka Murasaka, Makoto Nagata, Takashi Morie, Satoshi Iwata, “A Chip Level Substrate Noise Analysis Method Using F Matrix” IEICE Tech. Bulletin ICD99-147, September 1999) And AWE (Balshaz R. Stanisic, Nishath K. Verghese, Rob A. Rutenbar, L. Richard Carleyand, David J. Allistot, "Addressing Substrate Coupling Mixed-Mode IC's: Simulation and Power Distribution Synthesis," IEEE Journal of Solid-State. Circuit, Vol. 9, No. 3, pp. 226-238, March 1994) can be used to reduce the netlist. In step S3, the mesh division is corrected to change the mesh size. However, since the mesh structure itself does not change, an algorithm for degenerating the netlist in step S4 can be applied.

本実施形態例では、半導体集積回路の基板部分を、メッシュ状の、サイズが異なるセルを組み合わせて分割することでモデル化する基板モデルにおいて、大きい方のセル側と小さい方のセル側の抵抗比を用いて、ショートカットパスが解析精度に与える影響を予測する。ショートカットパスの影響が大きいと予想されるときには、基板モデルのメッシュ分割が適当でないと判断して、抵抗比が所定の値未満になるように、セルの大きさを修正する。このため、サイズが異なるセルに分割することにより問題となるショートカットパスの影響を低く抑えた基板モデルを、簡易に作成することができる。このような基板モデルを、回路シミュレータなどで解析すると、計算規模を縮小しても、解析の精度に大きな問題が生じない。   In the present embodiment example, in a substrate model in which a substrate portion of a semiconductor integrated circuit is modeled by combining mesh-shaped cells of different sizes, the resistance ratio between the larger cell side and the smaller cell side Is used to predict the effect of the shortcut path on the analysis accuracy. When the influence of the shortcut path is expected to be large, it is determined that the mesh division of the substrate model is not appropriate, and the cell size is corrected so that the resistance ratio is less than a predetermined value. For this reason, it is possible to easily create a substrate model in which the influence of a shortcut path which is a problem is reduced by dividing the cell into cells having different sizes. When such a substrate model is analyzed by a circuit simulator or the like, even if the calculation scale is reduced, no serious problem occurs in the accuracy of the analysis.

(第2の実施形態例)
図6は、半導体集積回路の基板モデルを示しており、(a)は基板モデルのセル分割を、(b)セル分割の等価回路を示している。図6に示す基板モデルは、セル101を高さ方向に3分割している点で、図1の基板モデルと相違する。本発明の第2の実施形態例の基板モデル作成方法及び基板ノイズ解析方法は、基板モデルのメッシュ分割の仕方を除いて、第1の実施形態例の基板モデル作成方法及び基板ノイズ解析方法と同様である。
(Second Embodiment)
FIG. 6 shows a substrate model of a semiconductor integrated circuit. (A) shows cell division of the substrate model, and (b) shows an equivalent circuit of cell division. The substrate model shown in FIG. 6 is different from the substrate model of FIG. 1 in that the cell 101 is divided into three in the height direction. The substrate model creation method and the substrate noise analysis method according to the second embodiment of the present invention are the same as the substrate model creation method and the substrate noise analysis method according to the first embodiment except for the method of dividing the mesh of the substrate model. It is.

セル101の横方向のセル数をn、高さ方向のセル数をnvとすると、本実施形態例では、半導体集積回路は、lの方向にn分割、h1の方向にnv分割された小さなセル101の群と、それらと垂直方向に接続された1つの大きなセル102とでモデル化される。図6(a)に示す基板モデルは、n=5、nv=3の例である。この基板モデルから、同図(b)に示す等価回路が得られる。ここで、端子AB間の抵抗値と、図2に示したCD面の抵抗値の比である抵抗比xについて考えると、第1の実施形態例と同様に、式(1.11)が得られる。ここで、式(1.11)中のh1、h2、l,wは、セルの分割数にかかわらず一定である。従って、抵抗比xもセル数にかかわらず一定であり、ここではx=1.0であったとする。 In this embodiment, the semiconductor integrated circuit is divided into n in the l direction and n v in the h1 direction, where n is the number of cells in the horizontal direction and n v is the number of cells in the height direction. It is modeled by a group of small cells 101 and one large cell 102 connected vertically to them. The substrate model shown in FIG. 6A is an example where n = 5 and n v = 3. From this substrate model, an equivalent circuit shown in FIG. Here, considering the resistance ratio x, which is the ratio of the resistance value between the terminals AB and the resistance value of the CD surface shown in FIG. 2, the equation (1.11) is obtained as in the first embodiment. It is done. Here, h 1 , h 2 , l, and w in the formula (1.11) are constant regardless of the number of divided cells. Therefore, the resistance ratio x is also constant regardless of the number of cells, and it is assumed here that x = 1.0.

図7は、図6(b)に示す等価回路を用いて回路のネットリストを作成し、SPICEシミュレーションを行って、基板抵抗減少率Eを求めた結果をグラフとして示している。図7に示すように、高さ方向のセル数nvを、nv=1から9までの間で変化させても、グラフのプロット点はほぼ重なり、基板抵抗減少率Eにはほとんど影響を与えない。つまり、図6(a)に示すように、ウェルの高さ方向を複数のセルでモデル化した場合であっても、ショートカットパスの影響は、ウェルの高さ方向を1つのセルでモデル化した第1の実施形態例と同様に判断できる。 FIG. 7 is a graph showing the results of creating a circuit netlist using the equivalent circuit shown in FIG. 6B and performing SPICE simulation to obtain the substrate resistance reduction rate E. As shown in FIG. 7, even if the number of cells n v in the height direction is changed between n v = 1 and 9, the plot points of the graph almost overlap, and the substrate resistance reduction rate E is hardly affected. Don't give. That is, as shown in FIG. 6A, even if the well height direction is modeled by a plurality of cells, the effect of the shortcut path is that the well height direction is modeled by one cell. The determination can be made in the same manner as in the first embodiment.

本実施形態例では、図5に示すステップS1において、メッシュ分割手段18は、第1の実施形態例で小さなセルとして構成される領域を、その高さ方向を更に複数のセルに分割する。この場合、第1の実施形態例と同様の効果に加えて、小さなセルとして構成される領域を更に細かく解析できる。   In the present embodiment example, in step S1 shown in FIG. 5, the mesh dividing means 18 further divides the area configured as small cells in the first embodiment example into a plurality of cells. In this case, in addition to the same effects as those of the first embodiment, a region configured as a small cell can be analyzed in more detail.

(第3の実施形態例)
図8は、半導体集積回路の基板モデルを示しており、(a)は基板モデルのセル分割を、(b)(c)(d)はそれぞれ基板モデルの面を示している。図8に示す基板モデルは、ウェルを模した領域が、抵抗率が異なる2つの領域として構成される点で、図1の基板モデルと相違する。本発明の第3の実施形態例の基板モデル作成方法及び基板ノイズ解析方法は、基板モデルのメッシュ分割、及び、2つのショートカットパスによる影響を判断する点を除いて、図5に示す基板モデル作成方法及び基板ノイズ解析方法と同様である。
(Third embodiment)
FIG. 8 shows a substrate model of the semiconductor integrated circuit. (A) shows cell division of the substrate model, and (b), (c), and (d) show the surface of the substrate model. The substrate model shown in FIG. 8 is different from the substrate model shown in FIG. 1 in that a region simulating a well is configured as two regions having different resistivities. The board model creation method and board noise analysis method of the third exemplary embodiment of the present invention are the board model creation shown in FIG. 5 except that the mesh division of the board model and the influence of two shortcut paths are judged. This is the same as the method and the substrate noise analysis method.

半導体基板のモデルは、例えば、抵抗率がρ1及びρ3のウェル領域と、抵抗率がρ2の基板領域を模しているものとする。半導体集積回路は、基板モデルのセル分割において、図8(a)に示すように、抵抗率がρ1のウェル領域はそれぞれが(w1×h1×l/5)のサイズを有する、5つのセル101に、抵抗率がρ1のウェル領域は(w3×h1×l)のサイズを有するセル103に、基板領域は、(w1×h2×l)及び(w3×h2×l)のサイズを有する、2つのセル102にそれぞれメッシュ分割される。セル101は、それぞれ、その下面で前面側のセル102と、その背面側の側面でセル103と接している。 The model of the semiconductor substrate is assumed to mimic a well region having resistivity ρ 1 and ρ 3 and a substrate region having resistivity ρ 2 , for example. In the semiconductor integrated circuit, in the cell division of the substrate model, as shown in FIG. 8A, each of the well regions having the resistivity ρ 1 has a size of (w1 × h 1 × 1/5). In the cell 101, the well region having a resistivity of ρ 1 has a size of (w3 × h 1 × l), and the substrate regions have (w1 × h 2 × l) and (w3 × h 2 × l). Are divided into two cells 102 each having a size of Each of the cells 101 is in contact with the cell 102 on the front surface on the lower surface and the cell 103 on the side surface on the rear surface.

図8(a)に示す基板モデルでは、複数のセル101に対して、1つのセル102及び103が接している。ショートカットパスの影響の計算に際しては、図中Zで示すセル101の下面側のショートカットパスのみでなく、図中Yで示すショートカットパスの影響も考慮する必要がある。このため、本実施形態例では、基板モデルの作成に際して、図5に示すステップS2で、ショートカットパスZに加えて、ショートカットパスYが及ぼす影響を、基板抵抗減少率として算出し、メッシュ分割が妥当であるかを判断する。   In the substrate model shown in FIG. 8A, one cell 102 and 103 are in contact with a plurality of cells 101. In calculating the influence of the shortcut path, it is necessary to consider not only the shortcut path on the lower surface side of the cell 101 indicated by Z in the figure but also the influence of the shortcut path indicated by Y in the figure. For this reason, in this embodiment, when creating a board model, in step S2 shown in FIG. 5, the influence of the shortcut path Y in addition to the shortcut path Z is calculated as the board resistance reduction rate, and mesh division is appropriate. It is judged whether it is.

図8(a)の例では、基板抵抗減少率は、ショートカットパスY及びZによる基板抵抗減少率の和として求められる。具体的には、ショートカットパスYによる基板抵抗減少率E2と、ショートカットパスZによる基板抵抗減少率E1とを個別に求め、それらを合計し、ウェル領域全体の基板抵抗減少率Eとする。 In the example of FIG. 8A, the substrate resistance reduction rate is obtained as the sum of the substrate resistance reduction rates by the shortcut paths Y and Z. Specifically, the substrate resistance reduction rate E 2 due to the shortcut path Y and the substrate resistance reduction rate E 1 due to the shortcut path Z are individually obtained and summed to obtain the substrate resistance reduction rate E of the entire well region.

ショートカットパスZによる基板抵抗減少率をE1は、第1の実施形態例と同様に計算され、

Figure 0004377437
Figure 0004377437
Figure 0004377437
となる。ここで、RABは図8(b)に示す面AB間の抵抗、RCDは同図(c)に示す面CD間の抵抗であり、それぞれ、
Figure 0004377437
Figure 0004377437
として求められる。 The substrate resistance reduction rate E 1 by the shortcut path Z is calculated in the same manner as in the first embodiment,
Figure 0004377437
Figure 0004377437
Figure 0004377437
It becomes. Here, R AB is the resistance between the planes AB shown in FIG. 8B, and R CD is the resistance between the planes CD shown in FIG.
Figure 0004377437
Figure 0004377437
As required.

一方、ショートカットパスYによる基板抵抗減少率E2については、面の方向は異なるが、上記と同様に、

Figure 0004377437
Figure 0004377437
Figure 0004377437
と計算される。ここで、REFは、図8(d)に示す面EF間の抵抗であり、
Figure 0004377437
として求められる。ウェル領域全体の基板抵抗減少率Eは、式(1.18)と式(1.23)との和をとり、
E=E+E (1.25)
となる。 On the other hand, for the substrate resistance reduction rate E 2 by the shortcut path Y, the direction of the surface is different.
Figure 0004377437
Figure 0004377437
Figure 0004377437
Is calculated. Here, R EF is the resistance between the planes EF shown in FIG.
Figure 0004377437
As required. The substrate resistance reduction rate E of the entire well region is the sum of Equation (1.18) and Equation (1.23),
E = E 1 + E 2 (1.25)
It becomes.

本実施形態例では、基板モデルの作成及び基板ノイズの解析に際して、2つのショートカットパスの合計の基板抵抗減少率を求め、その基板抵抗減少率が所望の値以下になるようにセルのサイズを調整して、メッシュ分割をやり直す。このため、ショートカットパスが2つある場合でも、簡易にショートカットパスの影響を回避した基板モデルを得ることができる。   In this embodiment, when creating a substrate model and analyzing substrate noise, the total substrate resistance reduction rate of the two shortcut paths is obtained, and the cell size is adjusted so that the substrate resistance reduction rate is less than the desired value. Then redo the mesh division. Therefore, even when there are two shortcut paths, it is possible to easily obtain a board model that avoids the influence of the shortcut paths.

(第4の実施形態例)
図9は、半導体集積回路の基板モデルを示しており、(a)は基板モデルのセル分割を、(b)は基板モデルのセルの面を示している。図9に示す基板モデルは、セル101を格子状に分割している点で、図1の基板モデルと相違する。本発明の第4の実施形態例の基板モデル作成方法及び基板ノイズ解析方法は、基板モデルのメッシュ分割、及び、2つのショートカットパスによる影響を判断する点を除いて、図5に示す基板モデル作成方法及び基板ノイズ解析方法と同様である。
(Fourth embodiment)
FIG. 9 shows a substrate model of a semiconductor integrated circuit, where (a) shows cell division of the substrate model and (b) shows a cell surface of the substrate model. The substrate model shown in FIG. 9 is different from the substrate model shown in FIG. 1 in that the cells 101 are divided in a lattice shape. The board model creation method and the board noise analysis method of the fourth embodiment of the present invention are the board model creation shown in FIG. 5 except that the influence of the mesh division of the board model and two shortcut paths is determined. This is the same as the method and the substrate noise analysis method.

本実施形態例では、本導体集積回路は、セル101を、2次元的に格子状に並べたセルとしてモデル化する。図9(a)の例では、ウェル領域を5×3にメッシュ分割している。この場合には、図中に示すショートカットパスZの他に、ショートカットパスXができる。横方向のショートカットパスZは、端子AB間のショートカットパスとなり、奥行き方向のショートカットパスXは、同図(b)に示す端子GH間のショートカットパスになる。   In the present embodiment example, the conductor integrated circuit models the cell 101 as a two-dimensionally arranged cell. In the example of FIG. 9A, the well region is divided into 5 × 3 meshes. In this case, in addition to the shortcut path Z shown in FIG. The horizontal shortcut path Z is a shortcut path between terminals AB, and the shortcut path X in the depth direction is a shortcut path between terminals GH shown in FIG.

ここで、直交する電流経路は互いに干渉しあわないので、ショートカットパスXは、端子AB間の抵抗値に影響を与えない。また、ショートカットパスZは、端子GH間の抵抗値に影響を与えない。つまり、ショートカットパスXとショートカットパスZとによる基板抵抗減少率は、互いに独立となり、図5のステップS2では、それぞれの端子間の基板抵抗減少率が、所望の精度を満足しているかを確認する。   Here, since the orthogonal current paths do not interfere with each other, the shortcut path X does not affect the resistance value between the terminals AB. Further, the shortcut path Z does not affect the resistance value between the terminals GH. That is, the substrate resistance reduction rates by the shortcut path X and the shortcut path Z are independent from each other. In step S2 of FIG. 5, it is confirmed whether the substrate resistance reduction rates between the terminals satisfy the desired accuracy. .

ショートカットパスZによる、端子AB間の基板抵抗減少率E1は、第1の実施形態例と同様に、式(1.10)で計算される。また、ショートカットパスXによる、端子GH間の基板抵抗減少率E3は、

Figure 0004377437
Figure 0004377437
Figure 0004377437
となる。ここで、RGHは、図9(b)に示す面G−H間の抵抗値である。何れの方向でも、基板抵抗減少率は、セルのサイズに依存する。 The substrate resistance reduction rate E 1 between the terminals AB due to the shortcut path Z is calculated by the equation (1.10) similarly to the first embodiment. The substrate resistance reduction rate E 3 between the terminals GH due to the shortcut path X is
Figure 0004377437
Figure 0004377437
Figure 0004377437
It becomes. Here, R GH is a resistance value between the planes GH shown in FIG. In either direction, the substrate resistance reduction rate depends on the cell size.

本実施形態例では、2つの方向にショートカットパスができる場合に、それぞれの方向の基板抵抗減少率を求める。求められた基板抵抗減少率が、それぞれの端子間において所望の値以下であるか否かを判断し、セル分割をやり直してセルのサイズを調節することで、ショートカットパスによる基板ノイズ解析の精度の低下を防ぐことができる。   In this embodiment, when the shortcut path can be formed in two directions, the substrate resistance reduction rate in each direction is obtained. It is judged whether the obtained substrate resistance reduction rate is less than the desired value between the terminals, and the cell size is adjusted by re-dividing the cell, thereby improving the accuracy of the substrate noise analysis by the shortcut path. Decline can be prevented.

(第5の実施形態例)
図10は、半導体集積回路の基板モデルを示し、(a)〜(d)は、それぞれ基板モデルの等価回路を示している。図10に示す基板モデルは、サイズの異なるセルの接合部が、pn接合である点で、図1の基板モデルと相違する。本発明の第5の実施形態例の基板モデル作成方法及び基板ノイズ解析方法は、等価回路に静電容量の成分が含まれる点を除いて、図5に示す基板モデル作成方法及び基板ノイズ解析方法と同様である。
(Fifth embodiment)
FIG. 10 shows a substrate model of a semiconductor integrated circuit, and (a) to (d) show equivalent circuits of the substrate model, respectively. The substrate model shown in FIG. 10 is different from the substrate model of FIG. 1 in that the junctions of cells having different sizes are pn junctions. The board model creation method and board noise analysis method of the fifth embodiment of the present invention are the same as the board model creation method and board noise analysis method shown in FIG. 5 except that the equivalent circuit includes a capacitance component. It is the same.

例えば図24に示す半導体集積回路の基板で、ウェル領域111と基板領域113とのpn極性が異なる場合には、ウェル領域−基板領域間の接続は、接合容量を介したものになる。このような半導体集積回路をメッシュ分割して基板モデルを作成すると、図10(a)に示すようなセル分割と、等価回路とが得られる。等価回路では、セル101の高さ方向の抵抗R2a、及び、セル102の中心のノードとセル101とをそれぞれ接続する抵抗R2bが、静電容量Cpを介して接続されている。 For example, in the case of the substrate of the semiconductor integrated circuit shown in FIG. 24, when the pn polarity of the well region 111 and the substrate region 113 is different, the connection between the well region and the substrate region is via a junction capacitance. When such a semiconductor integrated circuit is divided into meshes to create a substrate model, cell division as shown in FIG. 10A and an equivalent circuit are obtained. In the equivalent circuit, a resistance R 2a in the height direction of the cell 101 and a resistance R 2b that connects the node at the center of the cell 102 and the cell 101 are connected via a capacitance C p .

多端子F行列を用いた基板モデル作成方法では、抵抗素子のみで構成されたネットリストであれば、行列演算により電流の出入りが無いノードを削除してネットリストを縮退することが可能である。しかし、容量素子は、F行列演算で縮退が不可能であり、基板モデルに多数の容量素子を含むと、基板モデルのネットリスト規模が増大し、回路解析に時間がかかる。このため、可能であれば、F行列によるネットリスト縮退の前に、容量素子の数を減らしておくのが望ましい。   In the substrate model creation method using the multi-terminal F matrix, if the net list is composed only of resistance elements, it is possible to degenerate the net list by deleting a node where current does not enter and exit by matrix calculation. However, the capacitive element cannot be degenerated by F matrix calculation. If the board model includes a large number of capacitive elements, the netlist size of the board model increases, and circuit analysis takes time. For this reason, if possible, it is desirable to reduce the number of capacitive elements before the net list degeneration by the F matrix.

セル102の抵抗素子R2bと容量素子Cpとを、図10(b)に示すように、単純に1つの素子で置き換えて、容量素子数を減らそうとすると、図1(c)と同様に、セル101同士が低インピーダンスで接続されるために、ショートカットパスの影響が大きい。図10(b)に示す等価回路では、ショートカットパスの影響で解析精度が目標に到達しない場合には、同図(a)示す等価回路のR2bを、容量素子Cpよりもセル101側に移動して、同図(c)に示す等価回路に変換する。同図(c)に示す等価回路の容量素子Cpを1つにまとめて、同図(d)に示す等価回路が得られる。同図(d)に示す等価回路は、同図(b)に示す等価回路よりも、ショートカットパスの影響が小さくなる。 As shown in FIG. 10B, when the resistance element R 2b and the capacitive element C p of the cell 102 are simply replaced with one element as shown in FIG. In addition, since the cells 101 are connected with low impedance, the influence of the shortcut path is large. In the equivalent circuit shown in FIG. 10 (b), when the analysis accuracy due to the influence of short-cut path does not reach the goal, the R 2b of the equivalent circuit shown FIG. (A), the cell 101 side of the capacitive element C p Move to convert to the equivalent circuit shown in FIG. By combining the capacitive elements C p of the equivalent circuit shown in FIG. 4C into one, the equivalent circuit shown in FIG. The equivalent circuit shown in FIG. 4D is less affected by the shortcut path than the equivalent circuit shown in FIG.

本実施形態例では、サイズの異なるセル同士の接合面が、pn接合の接合面である場合には、半導体集積回路の基板部分を、図10(c)又は図(d)のようにモデル化する。同図(d)に示すような基板モデルの等価回路を用いる場合には、容量素子Cpはショートカットパスの経路にはならないので、Cpの存在は無視して、抵抗素子に起因するショートカットパスのみを考慮することで、モデル全体のショートカットパスの影響を判断できる。このため、第1の実施形態例と同様に、ショートカットパスを構成する抵抗比xを算出することで、ショートカットパスの影響を判断することができる。 In the present embodiment, when the joint surface between cells of different sizes is a joint surface of a pn junction, the substrate portion of the semiconductor integrated circuit is modeled as shown in FIG. 10C or FIG. To do. When the equivalent circuit of the substrate model as shown in FIG. 6D is used, since the capacitive element C p does not serve as a shortcut path, the presence of C p is ignored and the shortcut path caused by the resistive element is ignored. By considering only this, the effect of the shortcut path of the entire model can be determined. Therefore, as in the first embodiment, the influence of the shortcut path can be determined by calculating the resistance ratio x constituting the shortcut path.

(第6の実施形態例)
上記第1〜第5の実施形態例では、例えば図1(b)に示すように、ノードPとそれぞれのセル101とを接続する抵抗R2は、セル101の位置にかかわらず、同じ抵抗値を持っている。しかし、実際には、ノードPから各セル101の中心のノードまでの距離は一様ではなく、従って、抵抗値は、ノードPとセル101の距離に応じて大きくなる。
(Sixth embodiment)
In the first to fifth embodiments, for example, as shown in FIG. 1B, the resistance R 2 that connects the node P and each cell 101 has the same resistance value regardless of the position of the cell 101. have. However, in practice, the distance from the node P to the center node of each cell 101 is not uniform, and therefore the resistance value increases according to the distance between the node P and the cell 101.

図11は、半導体集積回路の基板モデルを示し、(a)及び(b)はそれぞれ基板モデルの等価回路を示している。本発明の第6の実施形態例の基板モデル作成方法及び基板ノイズ解析方法は、ノードPからセル101までの距離に応じて抵抗R2に重み付けを行う点で、第1〜第5の実施形態例と相違する。 FIG. 11 shows a substrate model of a semiconductor integrated circuit, and (a) and (b) show equivalent circuits of the substrate model, respectively. The substrate model creation method and substrate noise analysis method of the sixth embodiment of the present invention are the first to fifth embodiments in that the resistance R 2 is weighted according to the distance from the node P to the cell 101. Different from the example.

図11(a)に示す等価回路は、抵抗R2がR2i(i=1からnまでの整数)に置き換わる点で、図1(b)に示す等価回路と相違する。図11(a)のR2iを、セル101の成分とセル102の成分とに分けると、同図(b)に示す等価回路が得られる。抵抗R2iと抵抗R2a、R2biとの関係は、
2i=R2a+R2bi (1.29)
である。抵抗R2aは、セル101の垂直方向の抵抗成分のため、セルの位置によって変化せず、その抵抗値は、

Figure 0004377437
となる。一方、セル102の抵抗R2biは、接続するセル102によってノード間の距離が異なるので、次式に示すような重みαi/βを持つものとする。
Figure 0004377437
The equivalent circuit shown in FIG. 11A is different from the equivalent circuit shown in FIG. 1B in that the resistor R 2 is replaced with R 2i (i = 1 to an integer from 1 to n). When R 2i in FIG. 11A is divided into a component of the cell 101 and a component of the cell 102, an equivalent circuit shown in FIG. 11B is obtained. The relationship between the resistor R 2i and the resistors R 2a and R 2bi is
R 2i = R 2a + R 2bi (1.29)
It is. Since the resistance R 2a is a resistance component in the vertical direction of the cell 101, it does not change depending on the position of the cell, and its resistance value is
Figure 0004377437
It becomes. On the other hand, the resistance R 2bi of the cell 102 has a weight α i / β as shown in the following equation because the distance between nodes differs depending on the cell 102 to be connected.
Figure 0004377437

αiの重み付けは、例えば、i=1を端子A側からみて1つ目のセル101としたとき、

Figure 0004377437
のように定義する。このように定義することで、重みαiは、ノードPと対象セルの下面のノードとの距離に応じたものとなり、対象セルがn個のセル101の中央のセルのときに、つまり、i=(n+1)/2のときに、重みαiが最小となり、端子A又はBに近いほど、重みαiが大きくなる。 The weighting of α i is, for example, when i = 1 is the first cell 101 when viewed from the terminal A side,
Figure 0004377437
Define as follows. By defining in this way, the weight α i corresponds to the distance between the node P and the node on the lower surface of the target cell, and when the target cell is the center cell of the n cells 101, that is, i When = (n + 1) / 2, the weight α i is minimum, and the closer to the terminal A or B, the larger the weight α i .

式(1.31)中のβは、R2b1〜R2bnの合成抵抗値が、式1.2中のR2bの項で示される値になるように、設定された定数であり、以下のようにして求める。R2biのアドミタンスYiは、

Figure 0004377437
となる。図11(b)に示す等価回路ではR2b1〜R2bnは並列接続されるため、その合成アドミタンスは、
Figure 0004377437
である。一方、式1.2中のR2bから、アドミタンスYは、
Figure 0004377437
で表される。その合成アドミタンスは、
Figure 0004377437
である。式(1.34)=式(1.36)とおいて、βを算出すると、
Figure 0004377437
となる。式(1.37)を、式(1.31)に代入すると、
Figure 0004377437
が得られる。ただし、αiは式(1.33)で与えられる。 Β in the formula (1.31) is a constant set so that the combined resistance value of R 2b1 to R 2bn becomes the value shown by the term of R 2b in the formula 1.2. So ask. R 2bi admittance Y i is
Figure 0004377437
It becomes. In the equivalent circuit shown in FIG. 11B, since R 2b1 to R 2bn are connected in parallel, the combined admittance is
Figure 0004377437
It is. On the other hand, from R 2b in Formula 1.2, the admittance Y is
Figure 0004377437
It is represented by Its synthetic admittance is
Figure 0004377437
It is. When β is calculated by formula (1.34) = formula (1.36),
Figure 0004377437
It becomes. Substituting equation (1.37) into equation (1.31),
Figure 0004377437
Is obtained. However, α i is given by equation (1.33).

図12は、図11の等価回路の具体的数値例を示している。図11(a)に示す各値を、l=5、w=5、h1=1、h2=5、ρ1=1、及び、ρ2=10にし、抵抗R2iを式(1.38)に設定したときには、図12に示すように、R2a=0.02、R2b1=R2b5=5.66、R2b2=R2b4=4.76、R2b3=4.42になった。5つのセル101の中心である3番目のセルの抵抗R2b3が一番小さく、端子A又はBを有するセルの抵抗が大きいことがわかる。 FIG. 12 shows a specific numerical example of the equivalent circuit of FIG. Each value shown in FIG. 11A is set to l = 5, w = 5, h 1 = 1, h 2 = 5, ρ 1 = 1, and ρ 2 = 10, and the resistance R 2i is expressed by the equation (1. when set to 38), as shown in FIG. 12, R 2a = 0.02, R 2b1 = R 2b5 = 5.66, R 2b2 = R 2b4 = 4.76, becomes R 2b3 = 4.42 . It can be seen that the resistance R 2b3 of the third cell which is the center of the five cells 101 is the smallest, and the resistance of the cell having the terminal A or B is large.

図13は、抵抗R2bに重み付けをした場合と、重み付けをしない場合との、抵抗比Xに対する基板抵抗減少率の変化をグラフとして示している。抵抗比xが同じときで比較すると、重みづけ無しの場合よりも、重みづけを行った場合のほうが、基板抵抗減少率Eは小さな値となる。また、ショートカットパスによる基板抵抗の減少率を10%以下にしたい場合には、重みづけ無しの場合はx>0.74とする必要があったが、重みづけを行った場合にはx>0.6となる。このため、メッシュ分割の柔軟性を高くすることができる。 FIG. 13 is a graph showing changes in the substrate resistance reduction rate with respect to the resistance ratio X when the resistance R 2b is weighted and when the resistance R 2b is not weighted. Comparing when the resistance ratio x is the same, the substrate resistance reduction rate E is smaller when weighting is performed than when weighting is not performed. Further, when it is desired to reduce the substrate resistance reduction rate by the shortcut path to 10% or less, it is necessary to set x> 0.74 when there is no weighting, but when weighting is performed, x> 0. .6. For this reason, the flexibility of mesh division can be increased.

図14は、半導体集積回路の基板モデルを示している。サイズの異なるセル間を接続する抵抗に重み付けを行った場合には、接合面に垂直な方向の抵抗成分への影響が懸念される。そこで、同図に示す基板モデルで、重み付け有りの場合と無しの場合とについて、面C-D間の合成抵抗値を求めた。図15は、面C−D間の合成抵抗と抵抗比xとの関係をグラフとして示している。同図に示すように、抵抗比xが0.1以上であれば、重み付けの有無による抵抗成分の差は発生しないといえる。一方、ショートカットパスの影響を小さくするためには、図13を参照すると、基板抵抗減少率を10%以下にする場合には、抵抗比xは0.6〜0.7よりも大きな値とする必要がある。つまり、ショートカットパスの影響を小さくするようなメッシュ分割をおこなっていれば、サイズの異なるセル間を接続する抵抗に重みづけを行った場合でも、接合面に垂直な方向の抵抗成分へは大きな影響は与えない。   FIG. 14 shows a substrate model of a semiconductor integrated circuit. When weights are applied to resistors that connect cells of different sizes, there is concern about the effect on the resistance component in the direction perpendicular to the joint surface. Therefore, in the substrate model shown in the figure, the combined resistance value between the planes CD was obtained for cases with and without weighting. FIG. 15 is a graph showing the relationship between the combined resistance between the planes CD and the resistance ratio x. As shown in the figure, if the resistance ratio x is 0.1 or more, it can be said that there is no difference in resistance component due to weighting. On the other hand, in order to reduce the influence of the shortcut path, referring to FIG. 13, when the substrate resistance reduction rate is 10% or less, the resistance ratio x is set to a value larger than 0.6 to 0.7. There is a need. In other words, if mesh division is performed to reduce the influence of the shortcut path, even if weighting is applied to the resistance connecting cells of different sizes, the resistance component in the direction perpendicular to the joint surface is greatly affected. Will not give.

図16は、本発明の第6の実施形態例の基板ノイズ解析装置を、その動作手順を示すフローチャートを兼用して示している。図16に示す基板ノイズ解析装置は、接続モデル適用手段51を備える点で、図5に示す第1の実施形態例の基板ノイズ解析装置とは相違する。以下、図16を参照して、本実施形態例の基板モデル作成方法、基板モデル作成装置、基板ノイズ解析方法、及び、基板ノイズ解析装置について説明する。   FIG. 16 shows a substrate noise analyzing apparatus according to the sixth embodiment of the present invention, which is also used as a flowchart showing its operation procedure. The board noise analyzing apparatus shown in FIG. 16 is different from the board noise analyzing apparatus of the first embodiment shown in FIG. Hereinafter, a substrate model creation method, a substrate model creation device, a substrate noise analysis method, and a substrate noise analysis device according to this embodiment will be described with reference to FIG.

接続モデル適用手段51は、サイズが異なるセル間を接続する抵抗に対して、式1.38に示すような重み付けを行う。本実施形態例では、ステップS2で、ショートカットパスが解析精度を低下させると判定されたときに、ステップS3に進んでメッシュ分割をやり直す代わりに、ステップS6に進んでサイズの異なるセル間を接続する抵抗に重み付けを行う。このため、メッシュ分割をやり直すことなく、ショートカットパスの影響を小さくすることができる。   The connection model application unit 51 performs weighting as shown in Expression 1.38 on the resistance that connects cells having different sizes. In this embodiment, when it is determined in step S2 that the shortcut path lowers the analysis accuracy, instead of proceeding to step S3 and performing mesh division again, the process proceeds to step S6 to connect cells of different sizes. Weight resistance. For this reason, the influence of the shortcut path can be reduced without redoing the mesh division.

本実施形態例では、サイズが異なるセル間を接続する抵抗に重みづけを行って、半導体集積回路の基板部分をモデル化する。これにより、メッシュ分割構造及びメッシュサイズを変化させること無く、ショートカットパスの影響を小さくすることができる。また、サイズが異なるセル間を接続する抵抗に重み付けをしただけでは、ショートカットパスの影響を十分に小さくすることが出来ない場合には、メッシュ分割を修正することもできる。   In the present embodiment example, the substrate portion of the semiconductor integrated circuit is modeled by weighting the resistors connecting the cells having different sizes. Thereby, the influence of the shortcut path can be reduced without changing the mesh division structure and the mesh size. Further, if the influence of the shortcut path cannot be reduced sufficiently by simply weighting the resistances connecting the cells having different sizes, the mesh division can be corrected.

ここで、式(1.3)の導出について説明する。図17は、図1と同様の基板モデルであり、(a)、(b)はそれぞれ等価回路を示している。図17を参照して、セル101の数nが奇数の場合について抵抗RABの導入について説明する。セル101の数が奇数の場合には、セル102の中心ノードであるノードPに対して外部からの電流の供給が無いと仮定すると、図中に示す中央のセル101のノードQと、ノードPとは、等電位になる。従って、ノードPとQとを接続する抵抗R2は除去しても、端子AB間のインピーダンスRABには影響を与えない。 Here, the derivation of Expression (1.3) will be described. FIG. 17 shows a substrate model similar to that in FIG. 1, and (a) and (b) each show an equivalent circuit. Referring to FIG. 17, the introduction of resistor R AB will be described for the case where the number n of cells 101 is an odd number. When the number of cells 101 is an odd number, assuming that no external current is supplied to the node P, which is the central node of the cell 102, the node Q and the node P of the center cell 101 shown in the figure. Becomes equipotential. Therefore, even if the resistor R 2 connecting the nodes P and Q is removed, the impedance R AB between the terminals AB is not affected.

ノードQとノードPとを接地して、ノードP、Qよりも端子B側の抵抗を全て削除すると、図17(b)に示す等価回路が得られる。回路の対称性を考慮すると、この等価回路の端子A−グランド間のインピーダンスを求めて、これを2倍することで、同図(a)に示す回路の端子A−B間のインピーダンスRABを求めることができる。なお、説明の簡略化のために、以後、R1=1、R2=Rとして説明する。 When the nodes Q and P are grounded and all the resistors on the terminal B side from the nodes P and Q are deleted, an equivalent circuit shown in FIG. 17B is obtained. Considering the symmetry of the circuit, the impedance between the terminal A and the ground of the equivalent circuit is obtained and doubled to obtain the impedance R AB between the terminals A and B of the circuit shown in FIG. Can be sought. For simplification of description, R 1 = 1 and R 2 = R will be described below.

図18は、セルの数nを変化させたときの、図17(b)に示す等価回路を示し、(a)はn=1、(b)はn=3、(c)はn=5、(d)はn=7のときの回路をそれぞれ示している。図17(b)に示す等価回路の、端子A−グランド間のインピーダンスを求める方法について説明する。まず、図18(b)に示す回路の端子A−グランド間のインピーダンスを求めると、

Figure 0004377437
となる。次に、図18(c)に示す回路の端子A−グランド間のインピーダンスを求めると、式(2.1)のZ2を利用することで、
Figure 0004377437
のように書くことができる。そして、図18に示す回路のような繰り返し構造の回路のインピーダンスは、漸化式で表現することができ、図18(a)に示す初期状態をZ1とすると、次式が得られる。
Figure 0004377437
FIG. 18 shows the equivalent circuit shown in FIG. 17B when the number n of cells is changed. FIG. 18A shows n = 1, FIG. 18B shows n = 3, and FIG. 18C shows n = 5. , (D) respectively show circuits when n = 7. A method for obtaining the impedance between the terminal A and the ground in the equivalent circuit shown in FIG. First, when the impedance between the terminal A and the ground of the circuit shown in FIG.
Figure 0004377437
It becomes. Next, when the impedance between the terminal A and the ground of the circuit shown in FIG. 18C is obtained, Z 2 in Expression (2.1) is used,
Figure 0004377437
Can be written as The impedance of the circuit of the repeating structure, such as the circuit shown in FIG. 18, the recurrence formula can be expressed, if the initial state shown in FIG. 18 (a) and Z 1, the following equation is obtained.
Figure 0004377437

式(2.3)を整理すると、

Figure 0004377437
となり、更に式(2.5)から特性方程式を求めてxについて解くと、
Figure 0004377437
Figure 0004377437
が得られる。このとき、α,βを次式のように定める。
Figure 0004377437
ここでbi
Figure 0004377437
のように定義し、これをZiについて解くと次式が得られる。
Figure 0004377437
Rearranging equation (2.3)
Figure 0004377437
Further, when a characteristic equation is obtained from Equation (2.5) and solved for x,
Figure 0004377437
Figure 0004377437
Is obtained. At this time, α and β are determined as follows.
Figure 0004377437
Where b i
Figure 0004377437
When this is solved for Z i , the following equation is obtained.
Figure 0004377437

式(2.11)を式(2.3)に代入すると、

Figure 0004377437
が得られるが、これをbi+1について解くと次式となる。
Figure 0004377437
一方、式(2.9)をRについて解くと、
R=α2+α、R=β2+β (2.14)
が得られる。 Substituting equation (2.11) into equation (2.3),
Figure 0004377437
If this is solved for bi + 1 , the following equation is obtained.
Figure 0004377437
On the other hand, when equation (2.9) is solved for R,
R = α 2 + α, R = β 2 + β (2.14)
Is obtained.

式(2.14)を式(2.13)に代入すると

Figure 0004377437
となる。一方、式(2.4)よりZ1=0だから、
Figure 0004377437
であり、biは次式のように一般化される。
Figure 0004377437
式(2.17)を式(2.11)に代入すると、次式のようにZiが定式化される。
Figure 0004377437
Substituting equation (2.14) into equation (2.13)
Figure 0004377437
It becomes. On the other hand, since Z 1 = 0 from equation (2.4),
Figure 0004377437
And b i is generalized as:
Figure 0004377437
By substituting equation (2.17) into equation (2.11), Z i is formulated as in the following equation.
Figure 0004377437

回路の対称性を考慮して、ZiからRABを求めると、

Figure 0004377437
が得られる。このとき、
Figure 0004377437
とすれば、
Figure 0004377437
となるが、ここではnは奇数なので、(−1)n=−1であり、式(2.20)は、
Figure 0004377437
となる。 In consideration of the symmetry of the circuit, R AB is obtained from Z i .
Figure 0004377437
Is obtained. At this time,
Figure 0004377437
given that,
Figure 0004377437
Here, since n is an odd number, (−1) n = −1, and equation (2.20) is
Figure 0004377437
It becomes.

図19は、図1と同様の基板モデルであり、(a)、(b)はそれぞれ等価回路を示している。図19を参照して、セル101の数nが偶数の場合の抵抗RABの導入について説明する。セル101の数nが、偶数の場合には、等価回路は図19(a)に示すような回路となる。回路の対称性を考慮すると、抵抗RABは、図19(b)に示す回路の端子A-グランド間のインピーダンスを算出し、これを2倍することで、求めることができる。 FIG. 19 shows a substrate model similar to that in FIG. 1, and (a) and (b) each show an equivalent circuit. Referring to FIG. 19, the introduction of resistor R AB when the number n of cells 101 is an even number will be described. When the number n of cells 101 is an even number, the equivalent circuit is a circuit as shown in FIG. Considering the symmetry of the circuit, the resistor RAB can be obtained by calculating the impedance between the terminal A and the ground of the circuit shown in FIG. 19B and doubling this.

1を、

Figure 0004377437
のように定義すれば、上記したnが奇数の場合と同じ方法でRABを定式化できる。式(2.22)を、式(2.10)に代入すると、
Figure 0004377437
が得られるが、更に式(2.14)を代入すると
Figure 0004377437
が得られる。式(2.24)と式(2.15)より、biは次式のように一般化される。
Figure 0004377437
式(2.25)を、式(2.11)に代入すると、次式のようにZiが定式化される。
Figure 0004377437
Z 1
Figure 0004377437
R AB can be formulated in the same way as in the case where n is an odd number. Substituting equation (2.22) into equation (2.10),
Figure 0004377437
Is obtained, but further substituting equation (2.14)
Figure 0004377437
Is obtained. From the equations (2.24) and (2.15), b i is generalized as the following equation.
Figure 0004377437
Substituting equation (2.25) into equation (2.11), Z i is formulated as in the following equation.
Figure 0004377437

回路の対称性を考慮してZiからRABを求めると、

Figure 0004377437
が得られる。このとき、
Figure 0004377437
とすれば、式(2.27)は、
Figure 0004377437
となるが、ここではnは偶数なので、(−1)n=1であり、式(2.28)は、
Figure 0004377437
となる。式(2.29)は、式(2.21)と同じ式であり、図17に示す回路の端子A−B間のインピーダンスRABは、セルの数nが偶数か奇数かにかかわらず、式(2.29)で計算できる。なお、式(2.29)は、説明の簡略化のためにR1=1、R2=Rとした場合の式であり、RABの抵抗値は、
Figure 0004377437
と表すことができる。 When R AB is obtained from Z i in consideration of the symmetry of the circuit,
Figure 0004377437
Is obtained. At this time,
Figure 0004377437
Then, the equation (2.27) becomes
Figure 0004377437
Here, since n is an even number, (−1) n = 1, and equation (2.28) is
Figure 0004377437
It becomes. Expression (2.29) is the same expression as Expression (2.21), and the impedance R AB between terminals AB of the circuit shown in FIG. 17 is equal to whether the number n of cells is even or odd. It can be calculated by equation (2.29). The expression (2.29) is an expression when R 1 = 1 and R 2 = R for simplification of explanation, and the resistance value of R AB is
Figure 0004377437
It can be expressed as.

なお、図8では、図面の簡略化のために、図中のウェル領域の奥行き方向の奥側のショートカットパスの影響のみを考慮したが、ウェル領域の奥行き方向の手前側にもセル103と同様なショートカットパスがある場合には、これを考慮すると良い。   In FIG. 8, for the sake of simplification of the drawing, only the influence of the shortcut path on the far side in the depth direction of the well region in the figure is taken into account. If there is a shortcut path, consider this.

また、第2から第6の実施形態例で示した基板モデルのセル分割は、組み合わせても良い。例えば、図20に示すように、第4の実施形態例で説明した格子状に配列されたセルを、第2の実施形態例で説明したセルの高さ方向を複数のセルに分割してモデル化することもできる。その際には、図5に示すメッシュ分割手段18及びメッシュ分割修正手段20は、図20に示すメッシュ分割に対応する手段として構成される。第2の実施形態例で説明したように、セル101の高さ方向を複数のセルでモデル化しても、ショートカットパスの影響は変化しない。従って、図20に示す基板モデルと、図9に示す基板モデルとは、同じアルゴリズムでショートカットパスの影響の算出が可能である。   Further, the cell division of the substrate model shown in the second to sixth embodiments may be combined. For example, as shown in FIG. 20, the cells arranged in the lattice shape described in the fourth embodiment are divided into a plurality of cells in the cell height direction described in the second embodiment. It can also be converted. In that case, the mesh division means 18 and the mesh division correction means 20 shown in FIG. 5 are configured as means corresponding to the mesh division shown in FIG. As described in the second embodiment, even if the height direction of the cell 101 is modeled by a plurality of cells, the influence of the shortcut path does not change. Therefore, the board model shown in FIG. 20 and the board model shown in FIG. 9 can calculate the influence of the shortcut path by the same algorithm.

メッシュ分割の修正については、計算機が自動で修正する方法と、解析者が前述のアルゴリズムに従って手動で修正する方法とが考えられる。第6の実施形態例では、ステップS2で、ショートカットパスが解析精度を低下させると判定されたときは、ステップS3又はステップS6の何れに進むかは、解析者が選択して決定することもできる。また、αiの重み付けは、前述のものに限られず、様々な重みの付け方の中から選ぶことができる。 Regarding the mesh division correction, there are a method in which the computer automatically corrects, and a method in which the analyst manually corrects according to the algorithm described above. In the sixth embodiment, when it is determined in step S2 that the shortcut path reduces the analysis accuracy, the analyzer can select and decide whether to proceed to step S3 or step S6. . Further, the weighting of α i is not limited to the above, and can be selected from various weighting methods.

以上、本発明をその好適な実施形態例に基づいて説明したが、本発明の基板モデル作成装置及び方法、並びに基板ノイズ解析装置及び方法は、上記実施形態例にのみ限定されるものでなく、上記実施形態例の構成から種々の修正及び変更を施した基板モデル作成装置及び方法、並びに基板ノイズ解析装置及び方法も、本発明の範囲に含まれる。   As mentioned above, although this invention was demonstrated based on the suitable embodiment example, the board | substrate model creation apparatus and method of this invention, and a board | substrate noise analysis apparatus and method are not limited only to the said embodiment example, Substrate model creation devices and methods, and substrate noise analysis devices and methods in which various modifications and changes have been made from the configuration of the above embodiment, are also included in the scope of the present invention.

半導体集積回路の基板モデルを示す模式立体図であり、(a)はメッシュ分割を、(b)、(c)、(d)は基板モデルの等価回路を共に示している。It is a schematic three-dimensional view showing a substrate model of a semiconductor integrated circuit, where (a) shows mesh division and (b), (c), and (d) show both equivalent circuits of the substrate model. 半導体集積回路の基板モデルを示す模式立体図。FIG. 3 is a schematic three-dimensional view showing a substrate model of a semiconductor integrated circuit. 抵抗比xと基板抵抗減少率Eとの関係を示すグラフ。The graph which shows the relationship between resistance ratio x and board | substrate resistance decreasing rate E. FIG. 基板抵抗減少率Eを所望の値にする抵抗比xminを示す表。The table | surface which shows resistance ratio xmin which makes the board | substrate resistance decreasing rate E a desired value. 本発明の第1の実施形態例の基板ノイズ解析装置を、その動作手順を示すフローチャートと共に示すブロック図。The block diagram which shows the board | substrate noise analysis apparatus of the 1st Example of this invention with the flowchart which shows the operation | movement procedure. 半導体集積回路の基板モデルを示す模式立体図であり、(a)はメッシュ分割を、(b)等価回路を共に示している。It is a schematic three-dimensional view showing a substrate model of a semiconductor integrated circuit, (a) shows mesh division and (b) shows an equivalent circuit together. 基板抵抗減少率と、メッシュ分割数との関係を示すグラフ。The graph which shows the relationship between a board | substrate resistance decreasing rate and the number of mesh divisions. 半導体集積回路の基板モデルを示す模式立体図であり、(a)、(b)、(c)、(d)はそれぞれメッシュ分割を示している。It is a schematic three-dimensional view showing a substrate model of a semiconductor integrated circuit, and (a), (b), (c), and (d) each show mesh division. 半導体集積回路の基板モデルを示す模式立体図であり、(a)、(b)はそれぞれメッシュ分割を示している。It is a schematic three-dimensional view showing a substrate model of a semiconductor integrated circuit, and (a) and (b) show mesh division, respectively. 半導体集積回路の基板モデルを示す模式立体図であり、(a)、(b)、(c)、(d)はそれぞれメッシュ分割を等価回路と共に示している。It is a schematic three-dimensional view showing a substrate model of a semiconductor integrated circuit, and (a), (b), (c), and (d) show mesh division together with an equivalent circuit. 半導体集積回路の基板モデルを示す模式立体図であり、(a)、(b)はそれぞれメッシュ分割を等価回路と共に示している。It is a schematic three-dimensional view showing a substrate model of a semiconductor integrated circuit, and (a) and (b) show mesh division together with an equivalent circuit. 半導体集積回路の基板モデルを等価回路と共に示す模式立体図。The model three-dimensional view which shows the board | substrate model of a semiconductor integrated circuit with an equivalent circuit. 基板抵抗減少率Eと抵抗比xとの関係を、重み付けがある場合及びない場合に付いて示すグラフ。The graph which shows the relationship between the board | substrate resistance decreasing rate E and resistance ratio x with and without weighting. 半導体集積回路の基板モデルを等価回路と共に示す模式立体図。The model three-dimensional view which shows the board | substrate model of a semiconductor integrated circuit with an equivalent circuit. 接続面に垂直な方向の基板抵抗と抵抗比xとの関係を示すグラフ。The graph which shows the relationship between the board | substrate resistance of a direction perpendicular | vertical to a connection surface, and resistance ratio x. 本発明の第6の実施形態例の基板ノイズ解析装置を、その動作手順を示すフローチャートと共に示すブロック図。The block diagram which shows the board | substrate noise analysis apparatus of the 6th Example of this invention with the flowchart which shows the operation | movement procedure. 半導体集積回路の基板モデルを示し、(a)は基板モデルのセル分割を等価回路と共に示す模式ブロック図、(b)は簡略化した等価回路を示す回路図。1A is a schematic block diagram illustrating a cell model of a semiconductor integrated circuit together with an equivalent circuit, and FIG. 2B is a circuit diagram illustrating a simplified equivalent circuit. 図17のセル分割数を変化させたときの図17(b)に示す回路を示す回路図であり、(a)はn=1、(b)はn=3、(c)はn=5、(d)はn=7のときの回路を示す。FIG. 18B is a circuit diagram showing the circuit shown in FIG. 17B when the number of cell divisions in FIG. 17 is changed, where FIG. 17A is n = 1, FIG. 17B is n = 3, and FIG. 17C is n = 5. , (D) shows a circuit when n = 7. 半導体集積回路の基板モデルを示し、(a)は基板モデルのセル分割を等価回路と共に示す模式ブロック図、(b)は簡略化した等価回路を示す回路図。1A is a schematic block diagram illustrating a cell model of a semiconductor integrated circuit together with an equivalent circuit, and FIG. 2B is a circuit diagram illustrating a simplified equivalent circuit. 半導体集積回路の基板モデルを示す模式立体図。FIG. 3 is a schematic three-dimensional view showing a substrate model of a semiconductor integrated circuit. 図22の単位セルモデル201を示す模式立体図。FIG. 23 is a schematic three-dimensional view showing the unit cell model 201 of FIG. 22. 半導体集積回路の基板モデルのメッシュ分割を示す模式立体図。FIG. 3 is a schematic three-dimensional view showing mesh division of a substrate model of a semiconductor integrated circuit. 図22の基板モデルの等価回路を示す回路図。The circuit diagram which shows the equivalent circuit of the board | substrate model of FIG. 半導体集積回路の一例を示す断面図。FIG. 6 is a cross-sectional view illustrating an example of a semiconductor integrated circuit. 異なるセルサイズを有するセルに分割する基板モデルを示す模式断面図であり、(a)はセル分割を、(b)は等価回路と共に示している。It is a schematic cross section which shows the board | substrate model divided | segmented into the cell which has a different cell size, (a) has shown cell division, (b) has shown with the equivalent circuit.

符号の説明Explanation of symbols

11 入力装置
12 データ処理装置
13 記憶装置
14 出力装置
15 レイアウトデータ記憶部
16 デバイスパラメータ記憶部
17 基板モデルネットリスト記憶部
18 メッシュ分割手段
19 メッシュ分割の妥当性判定手段
20 メッシュ分割修正手段
21 ネットリスト変換手段
22 シミュレーション手段
62 抵抗素子
101、102、103:基板モデルのセル
111 ウェル領域
112 拡散層
113 基板領域
117 セル分割線
201 単位セルモデル
203 抵抗要素
11 Input device 12 Data processing device 13 Storage device 14 Output device 15 Layout data storage unit 16 Device parameter storage unit 17 Substrate model netlist storage unit 18 Mesh division means 19 Mesh division validity determination means 20 Mesh division correction means 21 Netlist Conversion means 22 Simulation means 62 Resistance elements 101, 102, 103: Cell 111 of substrate model Well region 112 Diffusion layer 113 Substrate region 117 Cell dividing line 201 Unit cell model 203 Resistance element

Claims (6)

半導体集積回路における基板部分を、複数種類のサイズを持つセルにメッシュ分割し、各セルを1個以上の抵抗素子から成る、又は、抵抗素子及び容量素子から成るインピーダンス要素に模擬することで、前記基板部分を等価回路でモデル化する基板モデル作成装置であって、
第1のサイズを有する複数の第1セル領域と、該第1セル領域に隣接し前記第1のサイズよりも大きな第2のサイズを有する1つの第2セル領域との間の境界がpn接合を形成する際に、前記第1セル領域の夫々と前記第2セル領域との間の接続を、前記第1セル領域のインピーダンス要素で定まるインピーダンスと、pn接合で定まる容量素子と、前記第2セル領域のインピーダンス要素で定まるインピーダンスとの直列接続で表現すると共に、前記容量素子が第2セル領域に直接的に接続されるように素子の配置ないし移動を行い、前記第1セル領域と前記第2セル領域との間に分岐点を設け、前記分岐点を前記第2セル領域のインピーダンス要素より第2セル領域側に設けた回路をネットリストに変換するネットリスト変換手段と、
前記メッシュ分割に際して分割の妥当性を判定する判定手段とを備え、
該判定手段は、
前記第1セル領域と前記第2セル領域との間の境界と平行なインピーダンスを、前記第
1セル領域を表現する第1のインピーダンスとして算出し、
前記第1セル領域と前記第2セル領域との間の接続を表現する第2のインピーダンス要
素を算出し、
前記第1セル領域と前記第2セル領域との間のpn接合容量を表現するインピーダンス
を前記第2のインピーダンス要素から減算して第3のインピーダンス要素を算出し、
前記第1のインピーダンスと前記第3のインピーダンスとの比率と、モデル精度に関連
づけられた所定値とを比較することにより、基板モデルが所望の解析精度を満足するか否
かの判定を行うことを特徴とする基板モデル作成装置。
The substrate portion in the semiconductor integrated circuit is mesh-divided into cells having a plurality of types of sizes, and each cell is composed of one or more resistance elements, or is simulated by an impedance element composed of a resistance element and a capacitance element. A board model creation device for modeling a board part with an equivalent circuit,
A boundary between a plurality of first cell regions having a first size and a second cell region adjacent to the first cell region and having a second size larger than the first size is a pn junction. Forming a connection between each of the first cell regions and the second cell region, an impedance determined by an impedance element of the first cell region, a capacitive element determined by a pn junction, and the second It is expressed by a series connection with an impedance determined by the impedance element of the cell region, and the element is arranged or moved so that the capacitive element is directly connected to the second cell region, and the first cell region and the first cell A net list conversion means for converting a circuit in which a branch point is provided between two cell regions, and a circuit in which the branch point is provided on the second cell region side from the impedance element of the second cell region into a net list;
Determination means for determining the validity of the division at the time of the mesh division,
The determination means includes
Calculating an impedance parallel to a boundary between the first cell region and the second cell region as a first impedance representing the first cell region;
Calculating a second impedance element representing a connection between the first cell region and the second cell region;
A third impedance element is calculated by subtracting an impedance representing a pn junction capacitance between the first cell area and the second cell area from the second impedance element;
By comparing the ratio between the first impedance and the third impedance and a predetermined value associated with the model accuracy, it is determined whether or not the substrate model satisfies a desired analysis accuracy. Characteristic board model creation device.
前記判定手段が、所望の解析精度を満足しないと判定すると、メッシュサイズを変更して再度メッシュ分割を行うメッシュ分割修正手段を更に備える、請求項1に記載の基板モデル作成装置。   The substrate model creation device according to claim 1, further comprising a mesh division correction unit that changes the mesh size and performs mesh division again when the determination unit determines that the desired analysis accuracy is not satisfied. 前記メッシュ分割修正手段は、前記インピーダンスの比率と前記所定値との違いに依存してメッシュ分割サイズを修正する、請求項2に記載の基板モデル作成装置。   The board model creation device according to claim 2, wherein the mesh division correction unit corrects the mesh division size depending on a difference between the impedance ratio and the predetermined value. 前記判定手段が、所望の解析精度を満足しないと判定すると、前記インピーダンスの比率と前記所定値との違いに依存して、サイズが異なるセル相互の接続部分を模擬するインピーダンス要素のインピーダンスを重み付けする補正を行う接続モデル適用手段を更に備える、請求項1〜3の何れかに記載の基板モデル作成装置。   If the determination means determines that the desired analysis accuracy is not satisfied, the impedance of the impedance element that simulates the connection portion of cells having different sizes is weighted depending on the difference between the impedance ratio and the predetermined value. The board model creation device according to claim 1, further comprising connection model application means for performing correction. 請求項1〜4の何れかに記載の基板モデル作成装置が作成した基板モデルに対し、基板上に作りこまれた機能回路のネットリストを結合して回路シミュレーションを行うシミュレーション手段を備えることを特徴とする基板ノイズ解析装置。   5. A simulation means for performing circuit simulation by combining a netlist of functional circuits built on a substrate with respect to the substrate model created by the substrate model creation device according to claim 1. Substrate noise analyzer. ネットリスト変換手段と、シミュレーション手段とを備えるデータ処理装置を用いて回路シミュレーションを行い、基板ノイズ解析を行う方法であって、
前記データ処理装置を用い、半導体集積回路における基板部分を、複数種類のサイズを持つセルにメッシュ分割し、各セルを1個以上の抵抗素子から成る、又は、抵抗素子及び容量素子から成るインピーダンス要素に模擬することで、前記基板部分を等価回路でモデル化することによって、基板モデルを作成するステップと、
前記シミュレーション手段が、前記基板モデルに対し、基板上に作りこまれた機能回路のネットリストを結合して回路シミュレーションを行うステップと、
を有し、
基板モデルを作成する前記ステップは、
前記データ処理装置の前記ネットリスト変換手段が、第1のサイズを有する複数の第1セル領域と、該第1セル領域に隣接し前記第1のサイズよりも大きな第2のサイズを有する1つの第2セル領域との間の境界がpn接合を形成する際に、前記第1セル領域の夫々と前記第2セル領域との間の接続を、前記第1セル領域のインピーダンス要素で定まるインピーダンスと、pn接合で定まる容量素子と、前記第2セル領域のインピーダンス要素で定まるインピーダンスとの直列接続で表現すると共に、前記容量素子が第2セル領域に直接的に接続されるように素子の配置ないし移動を行い、前記第1セル領域と前記第2セル領域との間に分岐点を設け、前記分岐点を前記第2セル領域のインピーダンス要素より第2セル領域側に設けた回路をネットリストに変換するステップを有することを特徴とする基板ノイズ解析方法。
A method of performing circuit simulation using a data processing device including a netlist conversion unit and a simulation unit, and performing substrate noise analysis,
Using the data processing apparatus, the substrate portion of the semiconductor integrated circuit is mesh-divided into cells having a plurality of sizes, and each cell is composed of one or more resistive elements, or impedance elements composed of resistive elements and capacitive elements. By creating a substrate model by modeling the substrate portion with an equivalent circuit,
It said simulation means, with respect to the substrate model, and performing a circuit simulation by combining a netlist of the functional circuit which is built on a substrate,
Have
The step of creating a substrate model includes:
The net list converting means of the data processing device includes a plurality of first cell areas having a first size and a second size adjacent to the first cell area and having a second size larger than the first size. When the boundary with the second cell region forms a pn junction, the connection between each of the first cell regions and the second cell region is an impedance determined by the impedance element of the first cell region. , A series connection of a capacitive element determined by a pn junction and an impedance determined by an impedance element of the second cell region, and the arrangement of the elements so that the capacitive element is directly connected to the second cell region. A circuit that moves, provides a branch point between the first cell region and the second cell region, and provides the branch point closer to the second cell region than the impedance element of the second cell region. Substrate noise analysis method characterized by comprising the step of converting the Ttorisuto.
JP2008059869A 2008-03-10 2008-03-10 Substrate model creation device, and substrate noise analysis device and method Expired - Fee Related JP4377437B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008059869A JP4377437B2 (en) 2008-03-10 2008-03-10 Substrate model creation device, and substrate noise analysis device and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008059869A JP4377437B2 (en) 2008-03-10 2008-03-10 Substrate model creation device, and substrate noise analysis device and method

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2002068220A Division JP4119140B2 (en) 2002-03-13 2002-03-13 Substrate model creation device and method, and substrate noise analysis device and method

Publications (2)

Publication Number Publication Date
JP2008165830A JP2008165830A (en) 2008-07-17
JP4377437B2 true JP4377437B2 (en) 2009-12-02

Family

ID=39695115

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008059869A Expired - Fee Related JP4377437B2 (en) 2008-03-10 2008-03-10 Substrate model creation device, and substrate noise analysis device and method

Country Status (1)

Country Link
JP (1) JP4377437B2 (en)

Also Published As

Publication number Publication date
JP2008165830A (en) 2008-07-17

Similar Documents

Publication Publication Date Title
Habal et al. Constraint-based layout-driven sizing of analog circuits
US9767240B2 (en) Temperature-aware integrated circuit design methods and systems
US20060036980A1 (en) Method and apparatus for jitter analysis and program therefor
JP2006209590A (en) Electromagnetic field analysis device, analysis method, and analysis program
JPH07501906A (en) Computer-aided design method and apparatus for multilevel interconnect technology
JP5165564B2 (en) Method and apparatus for supporting integrated circuit design
US6868374B1 (en) Method of power distribution analysis for I/O circuits in ASIC designs
US20220215152A1 (en) Generation of layout including power delivery network
Suda et al. A 65 nm programmable analog device array (PANDA) for analog circuit emulation
US8850374B2 (en) Method of reducing parasitic mismatch
CN112307702A (en) System, method and computer readable medium for developing electronic architectural design layouts
KR100704577B1 (en) Timing budget designing method
US20060190854A1 (en) Method for incorporating pattern dependent effects in circuit simulations
JP4325274B2 (en) Semiconductor device model creation method and apparatus
JP4119140B2 (en) Substrate model creation device and method, and substrate noise analysis device and method
JP5071081B2 (en) Simultaneous operation signal noise basic characteristic acquisition method and program in simultaneous operation signal noise estimation method for semiconductor device
Jiang et al. MIRID: Mixed-mode IR-drop induced delay simulator
JP4377437B2 (en) Substrate model creation device, and substrate noise analysis device and method
JP2008165831A (en) Substrate model preparing device and method, and substrate noise analysis device and method
US11042682B2 (en) Analog design tool having a cell set, and related methods
JP2009276822A (en) Semiconductor device design support device and semiconductor device design support method
JP5832252B2 (en) Noise analysis model and noise analysis method
JP4241582B2 (en) Layout creating apparatus and layout creating method
US20240111937A1 (en) Machine learning tool for layout design of printed circuit board
US12032896B2 (en) Generation of layout including power delivery network

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20080728

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20080926

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20081030

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20081226

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20090126

RD01 Notification of change of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7421

Effective date: 20090225

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20090226

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090326

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20090406

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090519

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090721

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: 20090811

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: 20090910

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: 20120918

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120918

Year of fee payment: 3

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120918

Year of fee payment: 3

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120918

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130918

Year of fee payment: 4

LAPS Cancellation because of no payment of annual fees