CN112290955B - Grid node coding method and system based on integrated circuit impedance network extraction - Google Patents
Grid node coding method and system based on integrated circuit impedance network extraction Download PDFInfo
- Publication number
- CN112290955B CN112290955B CN202011513386.8A CN202011513386A CN112290955B CN 112290955 B CN112290955 B CN 112290955B CN 202011513386 A CN202011513386 A CN 202011513386A CN 112290955 B CN112290955 B CN 112290955B
- Authority
- CN
- China
- Prior art keywords
- mesh
- grid
- triangular
- port
- nodes
- 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.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M13/00—Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
- H03M13/25—Error detection or forward error correction by signal space coding, i.e. adding redundancy in the signal constellation, e.g. Trellis Coded Modulation [TCM]
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M13/00—Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
- H03M13/03—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words
- H03M13/23—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words using convolutional codes, e.g. unit memory codes
Landscapes
- Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Design And Manufacture Of Integrated Circuits (AREA)
Abstract
The invention relates to a grid node coding method and a grid node coding system based on integrated circuit impedance network extraction. The method comprises the steps of carrying out mesh subdivision on a plurality of acquired multilayer integrated circuit layouts to obtain triangular meshes; sequentially and continuously coding ports in the triangular meshes; acquiring the distance from a mesh node to a port in a triangular mesh; coding the corresponding mesh nodes according to the distances from the mesh nodes to the ports in the triangular mesh; according to the information of the grid nodes on the coded triangular grid, column writing a finite element equation set of a potential field of the computing integrated circuit to obtain a finite element sparse matrix; according to the coding of the grid nodes of the finite element grid, repeatedly utilizing triangle-star transformation to eliminate the non-port grid nodes of the finite element sparse matrix to obtain a multi-port network admittance matrix; and extracting an impedance network model according to the multiport network admittance matrix. The sequential elimination method provided by the invention improves the accuracy of the extraction of the impedance network model of the multilayer very large scale integrated circuit.
Description
Technical Field
The invention relates to the field of integrated circuits, in particular to a grid node coding method and a grid node coding system based on integrated circuit impedance network extraction.
Background
The extraction of the multilayer impedance network model of the very large scale integrated circuit is an important work of the back end verification of the integrated circuit, and the conventional extraction method is that an S parameter matrix of the very large scale integrated circuit is calculated by a simplified transmission line method or a field-based calculation method, and then the S parameter matrix is converted into the impedance network model based on the S parameter matrix. The transmission line method has high calculation speed, but because the integrated circuit layout is greatly approximated, the layout calculation result which has a simple and regular structure at the early stage of processing is accurate, but the layout which has more and more complex structure and a scale range from centimeter to nanometer in recent years generates an unmatchable error. The method comprises the steps of calculating based on a field method, considering a complex structure of a layout in the mesh generation process, and therefore the calculation result is accurate, but the method needs to completely solve the whole sparse matrix to obtain the field distribution, then calculating the size (such as current, voltage, power and the like) of the port quantity aiming at a predefined port according to the field distribution, thereby calculating an S parameter matrix, and finally converting the S parameter matrix into an impedance network model. The method can provide a very accurate calculation result aiming at the complex layout structure of the integrated circuit, but the time for solving the whole sparse matrix is very long, and after the field distribution is obtained, an additional process is needed for calculating the S parameter matrix and then converting the S parameter matrix into the impedance network model, so that the method has long calculation time and complex process, and inevitable errors can be generated in the processes of calculating the S parameter matrix, converting the S parameter matrix into the impedance network model and the like.
Based on the above problems, in order to reduce the calculation time and reduce the calculation complexity, the field distribution is not calculated according to the conventional method, the S parameter matrix of the multiple ports is calculated by calculating the port quantity, and finally the S parameter matrix is converted into an impedance matrix, but the admittance matrix of the ports is obtained by directly eliminating the internal nodes of the finite element mesh of the non-ports by the admittance network represented by the finite element sparse matrix through a elimination method. However, in the elimination process, the elimination sequence of the internal nodes affects the accuracy of the finally obtained multiport network impedance matrix, and an improper elimination sequence of the internal nodes may cause inaccurate calculation results, so that the impedance network model cannot be effectively extracted.
Disclosure of Invention
The invention aims to provide a grid node coding method and a grid node coding system based on integrated circuit impedance network extraction, which improve the accuracy of multilayer very large scale integrated circuit impedance network model extraction.
In order to achieve the purpose, the invention provides the following scheme:
a mesh node coding method based on integrated circuit impedance network extraction comprises the following steps:
acquiring a multilayer integrated circuit layout; the multilayer integrated circuit layout comprises a plurality of polygons with a plurality of vertexes;
mesh subdivision is carried out on the multilayer integrated circuit layout to obtain a non-structural triangular mesh for dividing the multilayer integrated circuit layout;
sequentially and continuously coding ports in the triangular meshes; the ports are a plurality of ports in a multi-port network impedance matrix to be calculated, which is defined by a user in advance;
acquiring the distance from the mesh node in the triangular mesh to the port;
coding the corresponding mesh nodes according to the distances from the mesh nodes in the triangular mesh to the ports; the distance from the mesh node to the port in the triangular mesh is positively correlated with the code of the mesh node, namely the larger the distance is, the larger the code of the mesh node is;
according to the information of the grid nodes on the coded triangular grid, column writing a finite element equation set of a potential field of the computing integrated circuit to obtain a finite element sparse matrix; the nodes of the finite element mesh are associated with the finite element sparse matrix; the codes of the grid nodes of the finite element grid are the same as the codes of the grid nodes of the triangular grid after the codes are obtained; the encoding of the finite element mesh is the same as that of the triangular mesh after encoding; the information of the grid nodes comprises voltage, current, power, conductivity, potential, dielectric constant, permeability, electric field intensity and magnetic field intensity;
according to the codes of the grid nodes of the finite element grid, repeatedly utilizing triangle-star transformation to eliminate the non-port grid nodes of the finite element sparse matrix to obtain a multi-port network admittance matrix;
and extracting an impedance network model according to the multiport network admittance matrix.
Optionally, the mesh subdivision is performed on the multilayer integrated circuit layout to obtain a non-structural triangular mesh for dividing the multilayer integrated circuit layout, and then the method further includes:
and encoding the triangular meshes.
Optionally, the obtaining the distance from the mesh node in the triangular mesh to the port specifically includes:
initializing mesh nodes in the triangular mesh; the codes of the initialized grid nodes are all 0;
determining a triangular mesh associated with the encoded port from the encoded port; the triangle mesh associated with the encoded port is a triangle mesh whose triangle vertices include the port;
obtaining unprocessed mesh nodes in the triangular mesh associated with the encoded port;
sequentially calculating the distance between the unprocessed grid nodes in the triangular grid associated with the coded port and the port;
sorting the distances between the unprocessed grid nodes in the triangular grid associated with the coded port and the port from small to large to obtain a distance sequence;
acquiring a neighbor triangular mesh having a common edge with the related triangular mesh;
sequentially calculating the distance between the unprocessed grid nodes in the neighbor triangular grid and the port;
inserting the distance between the unprocessed grid nodes in the neighbor triangular grid and the port into the original distance sequence to form an updated distance sequence from small to large;
acquiring a new neighbor triangular mesh having a common edge with the neighbor triangle, and returning to the step of sequentially calculating the distances between the unprocessed mesh nodes in the neighbor triangular mesh and the ports until the mesh nodes of all the triangular meshes are processed, so as to obtain an updated distance sequence arranged from small to large;
continuously coding the triangular mesh nodes outside all the ports according to the distance sequence arranged from small to large; the codes with small distance are small, and the codes with large distance are large.
Optionally, the eliminating the non-port mesh nodes of the sparse finite element matrix by repeatedly using triangle-star transformation according to the encoding of the mesh nodes of the finite element mesh to obtain a multi-port network admittance matrix, specifically including:
and sequentially eliminating the grid nodes except the ports of the finite element grid according to the sequence of grid node codes from large to small by utilizing the triangular-star transformation to obtain the admittance matrix of the multiport network.
A mesh node encoding system based on integrated circuit impedance network extraction, comprising:
the multilayer integrated circuit layout acquisition module is used for acquiring a multilayer integrated circuit layout; the multilayer integrated circuit layout comprises a plurality of polygons with a plurality of vertexes;
the mesh generation module is used for carrying out mesh generation on the multilayer integrated circuit layout to obtain a non-structural triangular mesh for dividing the multilayer integrated circuit layout;
the port coding module is used for sequentially and continuously coding the ports in the triangular meshes; the ports are a plurality of ports in a multi-port network impedance matrix to be calculated, which is defined by a user in advance;
a mesh node-to-port distance obtaining module, configured to obtain a distance from a mesh node in the triangular mesh to the port;
the mesh node coding module is used for coding the corresponding mesh nodes according to the distances from the mesh nodes in the triangular mesh to the ports; the distance from the mesh node to the port in the triangular mesh is positively correlated with the code of the mesh node, namely the larger the distance is, the larger the code of the mesh node is;
the finite element sparse matrix determining module is used for writing a finite element equation set of a potential field of the computing integrated circuit in a row mode according to the information of the grid nodes on the coded triangular grid to obtain a finite element sparse matrix; the nodes of the finite element mesh are associated with the finite element sparse matrix; the codes of the grid nodes of the finite element grid are the same as the codes of the grid nodes of the triangular grid after the codes are obtained; the encoding of the finite element mesh is the same as that of the triangular mesh after encoding; the information of the grid nodes comprises voltage, current, power, conductivity, potential, dielectric constant, permeability, electric field intensity and magnetic field intensity;
the admittance matrix determination module of the multiport network is used for eliminating the non-port grid nodes of the finite element sparse matrix by repeatedly utilizing triangle-star transformation according to the codes of the grid nodes of the finite element grid to obtain a multiport network admittance matrix;
and the impedance network model extraction module is used for extracting an impedance network model according to the multiport network admittance matrix.
Optionally, the method further includes:
and the triangular mesh coding module is used for coding the triangular mesh.
Optionally, the grid node-to-port distance obtaining module specifically includes:
the initialization unit is used for initializing the mesh nodes in the triangular mesh; the codes of the initialized grid nodes are all 0;
an associated triangle mesh determination unit for determining a triangle mesh associated with the encoded port from the encoded port; the triangle mesh associated with the encoded port is a triangle mesh whose triangle vertices include the port;
an unprocessed mesh node acquisition unit configured to acquire an unprocessed mesh node in the triangular mesh associated with the encoded port;
a first distance determining unit, configured to sequentially calculate distances between unprocessed mesh nodes in the triangular mesh associated with the encoded port and the port;
a distance sequence determining unit, configured to sort distances between unprocessed mesh nodes in the triangular mesh associated with the encoded port and the port from small to large to obtain a distance sequence;
a neighbor triangular mesh acquisition unit for acquiring a neighbor triangular mesh having a common edge with the related triangular mesh;
the second distance determining unit is used for sequentially calculating the distance between the unprocessed grid nodes in the neighbor triangular grid and the port;
a distance sequence updating unit, configured to insert the distance between an unprocessed mesh node in the neighboring triangular mesh and the port into an original distance sequence to form an updated small-to-large distance sequence;
the updated distance sequence determining unit is used for acquiring a new neighbor triangular mesh which has a common edge with the neighbor triangle, and returning to the step of sequentially calculating the distances between the unprocessed mesh nodes in the neighbor triangular mesh and the ports until the mesh nodes of all the triangular meshes are processed, so as to obtain an updated distance sequence which is arranged from small to large;
the continuous coding unit is used for continuously coding the triangular mesh nodes outside all the ports according to the distance sequence arranged from small to large; the codes with small distance are small, and the codes with large distance are large.
Optionally, the admittance matrix determining module of the multiport network specifically includes:
and the admittance matrix determining unit of the multiport network is used for sequentially eliminating the grid nodes except the ports of the finite element grid according to the sequence of grid node coding from large to small by utilizing the triangular-star transformation to obtain the admittance matrix of the multiport network.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the grid node coding method and system based on the integrated circuit impedance network extraction provided by the invention are characterized in that the grid nodes in the triangular grid are coded according to the distance between the grid nodes in the triangular grid and the port, and the grid nodes of the finite element sparse matrix non-port are eliminated according to the coding and by repeatedly utilizing triangle-star transformation, so that the multiport network admittance matrix is obtained. The phenomenon that the accuracy of the impedance network model is low due to the fact that more important grid nodes closer to the port are eliminated in the grid node elimination process is avoided, unimportant grid nodes far away from the port are eliminated at first, then important grid nodes close to the port are eliminated step by step, the situation that larger errors are amplified continuously in the elimination process due to the fact that grid nodes close to the port are eliminated at first is avoided, and the extraction accuracy of the multilayer very large scale integrated circuit impedance network model is improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
Fig. 1 is a schematic flow chart of a method for encoding a mesh node based on an integrated circuit impedance network extraction according to the present invention;
FIG. 2 is a diagram illustrating the encoding of mesh nodes and the encoding of a triangular mesh in the prior art;
FIG. 3 is a schematic diagram of the encoding of mesh nodes of a triangular mesh and the encoding of the triangular mesh provided by the present invention;
fig. 4 is a schematic structural diagram of a mesh node coding system based on integrated circuit impedance network extraction according to the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a grid node coding method and a grid node coding system based on integrated circuit impedance network extraction, which improve the accuracy of multilayer very large scale integrated circuit impedance network model extraction.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
Fig. 1 is a schematic flow chart of a method for encoding a mesh node based on an integrated circuit impedance network extraction, as shown in fig. 1, the method for encoding a mesh node based on an integrated circuit impedance network extraction provided by the present invention includes:
s101, obtaining a multilayer integrated circuit layout; the multi-level integrated circuit layout includes a plurality of polygons with a plurality of vertices.
Before S101, the method further includes:
determining a three-dimensional model of the multi-layer integrated circuit; the three-dimensional model refers to the fact that the distribution of electromagnetic parameters such as conductivity, permeability and dielectric constant in the model is a function of a three-dimensional space coordinate, and the distribution of field quantities such as potential, magnetic potential, electric field intensity and magnetic field intensity in the model is a function of the three-dimensional space coordinate;
and simplifying the three-dimensional model of the multilayer integrated circuit into a two-dimensional model.
And S102, performing mesh subdivision on the multilayer integrated circuit layout to obtain a non-structural triangular mesh for dividing the multilayer integrated circuit layout. Wherein, S102 is mesh generation performed in a two-dimensional model.
S102 specifically comprises the following steps:
for a direct-current field analysis model of a multilayer integrated circuit, forming a Delaunay triangular mesh with polygon vertexes as mesh nodes for a plurality of polygons of each layer of an integrated circuit layout according to a Delaunay triangulation algorithm, and simplifying layout polygons without losing precision based on the Delaunay triangular mesh; and performing self-adaptive subdivision on the integrated circuit calculation region defined by the polygon based on the simplified layout polygon to form an unstructured triangular mesh.
For an alternating electromagnetic field analysis model of a multilayer integrated circuit, a plurality of polygons of each layer are vertically projected to the same layer, a Delaunay triangular mesh with polygon vertexes as mesh nodes is formed according to a Delaunay triangulation algorithm, wherein each side of each polygon comprises preset polygon information of the polygon and layer information of the layer.
And merging the polygon information and the layer information of the projected and superposed polygon edges.
And aligning the Delaunay triangular mesh to each side of the plurality of polygons according to a side exchange method, simultaneously calculating intersection points of the sides of the polygons, and newly adding the intersection points as vertexes of the polygons and nodes of the Delaunay triangular mesh to form a first triangular mesh.
And on the basis of the first triangular mesh, overlapping the layer information of each polygon edge to all triangles in each polygon on the basis of Boolean operation.
And identifying and collecting the triangular and polygonal sides contained in each parallel flat plate field by a parallel flat plate field identification method according to the layer information of all the triangles and each polygonal side.
According to the calculation precision requirement and the public areas of different parallel flat plate fields, carrying out self-adaptive mesh subdivision processing on the triangles in each parallel flat plate field to finally form the non-structural triangular mesh.
After S102, further comprising:
and encoding the triangular meshes.
S103, sequentially and continuously coding ports in the triangular meshes; the ports are a plurality of ports in a multi-port network impedance matrix to be calculated, which is defined by a user in advance.
S104, obtaining the distance from the mesh node in the triangular mesh to the port.
S104 specifically comprises the following steps:
initializing mesh nodes in the triangular mesh; and the codes of the initialized grid nodes are all 0.
Determining a triangular mesh associated with the encoded port from the encoded port; the triangle mesh associated with the encoded port is the triangle mesh whose vertices include the port.
Unprocessed mesh nodes in the triangular mesh associated with the encoded port are obtained.
And sequentially calculating the distance between the unprocessed grid nodes in the triangular grid associated with the coded port and the port.
And sequencing the distances between the unprocessed grid nodes in the triangular grid associated with the coded ports and the ports from small to large to obtain a distance sequence.
And acquiring a neighbor triangular mesh with a common edge with the related triangular mesh.
And sequentially calculating the distance between the unprocessed grid nodes in the neighbor triangular grid and the port.
And inserting the distance between the unprocessed grid nodes in the neighbor triangular grid and the port into the original distance sequence to form an updated distance sequence from small to large.
And acquiring a new neighbor triangular mesh having a common edge with the neighbor triangle, and returning to the step of sequentially calculating the distances between the unprocessed mesh nodes in the neighbor triangular mesh and the ports until the mesh nodes of all the triangular meshes are processed, so as to obtain an updated distance sequence arranged from small to large.
Continuously coding the triangular mesh nodes outside all the ports according to the distance sequence arranged from small to large; the codes with small distance are small, and the codes with large distance are large.
S105, coding the corresponding grid node according to the distance from the grid node in the triangular grid to the port; the distance from the mesh node to the port in the triangular mesh is positively correlated with the code of the mesh node, i.e., the larger the distance, the larger the code of the mesh node.
According to the distribution and attenuation law of the field, the farther away from the position of the user-set port, the smaller the influence of the grid node and medium change on the port, and conversely, the closer to the position of the user-set port, the larger the influence of the grid node and medium change on the port.
Fig. 2 is a schematic diagram of coding of mesh nodes of a triangular mesh and coding of the triangular mesh in the prior art, fig. 3 is a schematic diagram of coding of mesh nodes of the triangular mesh and coding of the triangular mesh provided by the present invention, and as shown in fig. 2 and fig. 3, a specific coding process of mesh nodes in the triangular mesh is as follows:
(1) the mesh nodes corresponding to the 3 user-defined ports are numbered 1-3, as shown in fig. 3.
And numbering the rest grid nodes on the basis of port node coding. The numbering method is that the coded ports are taken as centers and are gradually transmitted to the unprocessed grid nodes through the triangle neighbors.
(2) Setting the coding state of all grid nodes to be 0, wherein the coding state indicates that all nodes are not processed, namely, defining an array encoded = 0; setting the initial distance sequence as an empty set: dlist = { }.
(3) The transfer status of all mesh cells (triangles) is set to 0, indicating that all triangle cells are not transferred, i.e. the array fed =0 is defined.
(4) And coding the ports defined by the user in sequence and setting the coding state of the ports. Suppose that the user-defined 3 ports are originally encoded as i in turn1，i2，i3The 3 ports are sequentially encoded as 1, 2, 3, and encoded (i) is set1)=1，encoded(i2)=1，encoded(i3)=1。
(5) And finding the triangles associated with the port nodes which are already coded to form a set T. The triangle associated with a node is defined as the triangle containing the node, and the triangle associated with the node is T as shown in port 1 in FIG. 11，T2，T3The triangle associated with port 2 is T6，T11，T12The triangle associated with port 3 is T16，T17，T18Then the set T = { T = { (T)1，T2，T3，T6，T11，T12，T16，T17，T18}; and marks the corresponding value of the fed array as 1, e.g., set fed (1) =1, fed (2) =1, fed (3) =1, fed (6) =1 ….
(6) And finding out unprocessed vertexes of all the triangles in the set T, sequentially calculating the weight distance from each vertex to a user set port node, and inserting the weight distances into the existing distance sequence from small to large. If the number of the found unprocessed vertices is 0, the process proceeds to step 9. Wherein the weight distance is defined as follows:
in the formulad a For mesh nodes in the triangular meshaThe weighted distance to the port or ports,mdefining the number of ports for the user, N the number of mesh nodes,w b is as followsbThe weight of each port is set to be,d ab is as followsaFrom mesh node to mesh nodebDistance of each port. The sorted distance sequence is as follows:
dlist={d12，d15，d2，d6，d4，d13，d5，d16，d7，d9}。
(7) and adding the neighbor which is not transmitted by the element of the set T into the set T to form a new set T. For the current set T, T2Is T7，T3Is T4…, the final set T is updated as:
T={T1，T2，T3，T6，T11，T12，T16，T17，T18，T7，T4，T5，T10，T15}
here the cell neighbors are defined as: the neighboring triangle is a triangle having a common side with the triangle.
(8) The set T is updated. If the element of the set T is passed, or all the neighbors of the element of the set T are passed, or the neighbors are also in the set T, the element is deleted from the set T. T is updated as:
T={T7，T4，T5，T10，T15and f, turning to the step (6).
(9) And recoding the grid nodes according to the sorted distance sequence order. Assume that the distance sequence of the final update is:
dlist={d12，d15，d2，d6，d4，d13，d5，d16，d7，d9，d1，d14，d3the encoded mesh nodes are shown in table 1, where table 1 is as follows:
| original node encoding | |
| 8 | 1 |
| 11 | 2 |
| 10 | 3 |
| 12 | 4 |
| 15 | 5 |
| 2 | 6 |
| 6 | 7 |
| 4 | 8 |
| 13 | 9 |
| 5 | 10 |
| 16 | 11 |
| 7 | 12 |
| 9 | 13 |
| 1 | 14 |
| 14 | 15 |
| 3 | 16 |
S106, writing a finite element equation set of a potential field of the computing integrated circuit in a row according to the information of the grid nodes on the coded triangular grid to obtain a finite element sparse matrix; the nodes of the finite element mesh are associated with the finite element sparse matrix; the codes of the grid nodes of the finite element grid are the same as the codes of the grid nodes of the triangular grid after the codes are obtained; the encoding of the finite element mesh is the same as that of the triangular mesh after encoding; the information of the grid nodes includes voltage, current, power, conductivity, potential, permittivity, permeability, electric field strength, and magnetic field strength.
S106 specifically comprises:
based on the simplified two-dimensional field problem, writing a finite element equation set of a potential field of the computing integrated circuit based on the split grid columns to obtain a finite element sparse matrix.
For the direct current field model, the three-dimensional model of the multilayer integrated circuit refers to the conductivity in the direct current field model
Potential of the electrodeuAll the distributions of (A) and (B) are three-dimensional space coordinatesx，y，z) I.e.:
，
which satisfies the following equation (1):
boundary condition (2):
in the formula (I), the compound is shown in the specification,
is a boundary of the first type and is,nis normal to the boundary of the second type,
is an electric potentialuAt the first kind boundary
Value of above, using
It is shown that,
is the bulk current density of the external circuit.
The dimension of an actual PCB or a chip packaged board in the multilayer ultra-large scale integrated circuit is far larger than the thickness of the metal layer, and the three-dimensional direct current field problem of the multilayer integrated circuit is simplified into a two-dimensional direct current field problem.
The field solving equation set established by the finite element method for the direct current electric field two-dimensional model of each layer of integrated circuit is an equation set (3):
in the formula (I), theI(u)In order to be a functional function,his the thickness of the metal layer or layers,
is the conductivity of grid cell e;
as a grid celleThe potential vector of the node of (a);
is the face of the grid cell e and,
as the density of the surface current, the current density,
representing grid cellseThe edge of (2). And (4) taking an extreme value of the equation set (3), namely forming a finite element sparse matrix for solving the potential field.
For the alternating electromagnetic field model, the three-dimensional model of the multilayer integrated circuit refers to a multilayer super-large gaugeDielectric constant in three-dimensional model of electromagnetic response characteristics in frequency domain simulation of modular integrated circuit
Magnetic permeability of
Electric field intensityEMagnetic field intensityHAll the distributions of (A) and (B) are three-dimensional space coordinatesx，y，z) I.e.:
，
，
，
. The function of the three-dimensional model satisfies the following equation:
in the formulaJFor the purpose of the applied current density distribution,
for the angular frequency simulated for the integrated circuit,
indicating the strength of the magnetic fieldHThe degree of rotation of the screw is reduced,
indicates the electric field intensityEThe degree of rotation of the screw is reduced,jis the unit of an imaginary number,j 2 =-1。
when the actual PCB or chip package size in the multilayer VLSI is far larger than the metal layer, the three-dimensional model of the electromagnetic response characteristic of a frequency point in the frequency domain simulation of the multilayer VLSI can be simplified into a two-dimensional model, and the dielectric constant in the model is at the moment
Magnetic permeability of
Electric field intensityEMagnetic field intensityHAll the distributions are two-dimensional plane coordinates (x，y) I.e.:
，
，
，
distribution thereof andzis irrelevant. And the potential in the fielduAnd surface current densityJ s Satisfies the following conditions:
in the formula
Respectively representx，y，zUnit vector of direction, EzOf electric field strengthzThe direction component of the light beam is,H x andH y respectively of magnetic field strengthxAndythe direction component of the light beam is,his the metal layer spacing.
Through the equivalence from the three-dimensional model to the two-dimensional model, the two-dimensional finite element functional extreme value formula corresponding to the two-dimensional model is obtained as follows:
in the formula (I), the compound is shown in the specification,
the functional is a functional, the minimum value of the functional is corresponding to the variation of the functional as 0,
the variation of the functional is represented by the functional,
as a grid celliThe surface admittance of the first and second electrodes,kis referred to askThe number of the boundaries is such that,
is a boundary
Opening boundary condition of (u)kIs a boundary
The distribution of the electric potential on the upper side,
indicating a position to the right of the boundary and infinitely close to the boundary,
indicating a position to the left of the boundary and infinitely close to the boundary,
presentation unitiThe area of (a) is,
as a grid celliThe current density of (a) is,
as a grid celliThe surface resistance of the glass substrate is higher than the surface resistance of the glass substrate,
as a grid celliThe potential of (2).
By using the extreme condition of the formula (3), a finite element sparse matrix for solving the potential field can be formed.
And S107, according to the codes of the grid nodes of the finite element grid, repeatedly utilizing triangle-star transformation to eliminate the non-port grid nodes of the finite element sparse matrix to obtain the multiport network admittance matrix.
S107 specifically comprises the following steps:
and sequentially eliminating the grid nodes except the ports of the finite element grid according to the sequence of grid node codes from large to small by utilizing the triangular-star transformation to obtain the admittance matrix of the multiport network.
And S108, extracting an impedance network model according to the multiport network admittance matrix.
Fig. 4 is a schematic structural diagram of a mesh node coding system extracted based on an integrated circuit impedance network, as shown in fig. 4, the mesh node coding system extracted based on the integrated circuit impedance network provided by the present invention includes: the system comprises a multilayer integrated circuit layout acquisition module 401, a grid subdivision module 402, a port coding module 403, a grid node-to-port distance acquisition module 404, a grid node coding module 405, a finite element sparse matrix determination module 406, an admittance matrix determination module 407 of a multiport network and an impedance network model extraction module 408.
The multilayer integrated circuit layout obtaining module 401 is configured to obtain a multilayer integrated circuit layout; the multi-level integrated circuit layout includes a plurality of polygons with a plurality of vertices.
The mesh generation module 402 is configured to perform mesh generation on the multilayer integrated circuit layout to obtain a non-structural triangular mesh for dividing the multilayer integrated circuit layout.
The port coding module 403 is configured to sequentially and continuously code the ports in the triangular mesh; the ports are a plurality of ports in a multi-port network impedance matrix to be calculated, which is defined by a user in advance.
The mesh node-to-port distance obtaining module 404 is configured to obtain distances from mesh nodes in the triangular mesh to the ports.
The mesh node encoding module 405 is configured to encode a corresponding mesh node according to a distance from the mesh node in the triangular mesh to the port; the distance from the mesh node to the port in the triangular mesh is positively correlated with the code of the mesh node, i.e., the larger the distance, the larger the code of the mesh node.
The finite element sparse matrix determining module 406 is configured to column write a finite element equation set for calculating a potential field of the integrated circuit according to the information of the mesh nodes on the encoded triangular mesh to obtain a finite element sparse matrix; the nodes of the finite element mesh are associated with the finite element sparse matrix; the codes of the grid nodes of the finite element grid are the same as the codes of the grid nodes of the triangular grid after the codes are obtained; the encoding of the finite element mesh is the same as that of the triangular mesh after encoding; the information of the grid nodes includes voltage, current, power, conductivity, potential, permittivity, permeability, electric field strength, and magnetic field strength.
The admittance matrix determination module 407 of the multiport network is configured to repeatedly utilize triangle-star transformation to eliminate non-port mesh nodes of the sparse matrix of the finite element according to the coding of the mesh nodes of the finite element mesh, so as to obtain an admittance matrix of the multiport network.
The impedance network model extraction module 408 is configured to extract an impedance network model from the multiport network admittance matrix.
The invention provides a grid node coding system based on integrated circuit impedance network extraction, which further comprises: and a triangular mesh coding module.
The triangle mesh coding module is used for coding the triangle mesh.
The mesh node-to-port distance obtaining module 404 specifically includes: the device comprises an initialization unit, an associated triangular mesh determination unit, an unprocessed mesh node acquisition unit, a first distance determination unit, a distance sequence determination unit, a neighbor triangular mesh acquisition unit, a second distance determination unit, a distance sequence updating unit, an updated distance sequence determination unit and a continuous coding unit.
The initialization unit is used for initializing the mesh nodes in the triangular mesh; and the codes of the initialized grid nodes are all 0.
The associated triangular mesh determining unit is used for determining a triangular mesh associated with the coded port according to the coded port; the triangle mesh associated with the encoded port is the triangle mesh whose vertices include the port.
The unprocessed mesh node acquisition unit is used for acquiring unprocessed mesh nodes in the triangular mesh associated with the coded port.
The first distance determining unit is used for sequentially calculating the distance between the unprocessed grid nodes in the triangular grid associated with the coded port and the port.
The distance sequence determining unit is used for sequencing the distances between the unprocessed grid nodes in the triangular grid associated with the coded ports and the ports from small to large to obtain a distance sequence.
The neighbor triangular mesh acquisition unit is used for acquiring a neighbor triangular mesh having a common edge with the related triangular mesh.
The second distance determining unit is used for sequentially calculating the distance between the unprocessed grid nodes in the neighbor triangular grid and the ports.
And the distance sequence updating unit is used for inserting the distance between the unprocessed grid nodes in the neighbor triangular grid and the ports into the original distance sequence to form an updated distance sequence from small to large.
The updated distance sequence determining unit is used for acquiring a new neighbor triangular mesh having a common edge with the neighbor triangle, and returning to the step of sequentially calculating the distances between the unprocessed mesh nodes in the neighbor triangular mesh and the ports until the mesh nodes of all the triangular meshes are processed, so as to obtain an updated distance sequence arranged from small to large.
The continuous coding unit is used for continuously coding the triangular mesh nodes outside all the ports according to the distance sequence arranged from small to large; the codes with small distance are small, and the codes with large distance are large.
The admittance matrix determination module 407 of the multiport network specifically includes: an encoding ordering unit and an admittance matrix determination unit of the multiport network.
And the admittance matrix determination unit of the multiport network is used for sequentially eliminating the grid nodes except the ports of the finite element grid according to the sequence of grid node coding from large to small by utilizing the triangular-star transformation to obtain the admittance matrix of the multiport network.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (8)
1. A mesh node coding method based on integrated circuit impedance network extraction is characterized by comprising the following steps:
acquiring a multilayer integrated circuit layout; the multilayer integrated circuit layout comprises a plurality of polygons with a plurality of vertexes;
mesh subdivision is carried out on the multilayer integrated circuit layout to obtain a non-structural triangular mesh for dividing the multilayer integrated circuit layout;
sequentially and continuously coding ports in the triangular meshes; the ports are a plurality of ports in a multi-port network impedance matrix to be calculated, which is defined by a user in advance;
acquiring the distance from the mesh node in the triangular mesh to the port;
coding the corresponding mesh nodes according to the distances from the mesh nodes in the triangular mesh to the ports; the distance from the mesh node to the port in the triangular mesh is positively correlated with the code of the mesh node, namely the larger the distance is, the larger the code of the mesh node is;
according to the information of the grid nodes on the coded triangular grid, column writing a finite element equation set of a potential field of the computing integrated circuit to obtain a finite element sparse matrix; the grid nodes of the finite element grid are associated with the finite element sparse matrix; the codes of the grid nodes of the finite element grid are the same as the codes of the grid nodes of the triangular grid after the codes are obtained; the encoding of the finite element mesh is the same as that of the triangular mesh after encoding; the information of the grid nodes comprises voltage, current, power, conductivity, potential, dielectric constant, permeability, electric field intensity and magnetic field intensity;
according to the codes of the grid nodes of the finite element grid, repeatedly utilizing triangle-star transformation to eliminate the non-port grid nodes of the finite element sparse matrix to obtain a multi-port network admittance matrix;
and extracting an impedance network model according to the multiport network admittance matrix.
2. The method of claim 1, wherein the mesh subdivision is performed on the multilayer integrated circuit layout to obtain unstructured triangular meshes for dividing the multilayer integrated circuit layout, and then the method further comprises:
and encoding the triangular meshes.
3. The method according to claim 1, wherein the obtaining of the distances from the mesh nodes in the triangular mesh to the ports specifically comprises:
initializing mesh nodes in the triangular mesh; the codes of the initialized grid nodes are all 0;
determining a triangular mesh associated with the encoded port from the encoded port; the triangle mesh associated with the encoded port is a triangle mesh whose triangle vertices include the port;
obtaining unprocessed mesh nodes in the triangular mesh associated with the encoded port;
sequentially calculating the distance between the unprocessed grid nodes in the triangular grid associated with the coded port and the port;
sorting the distances between the unprocessed grid nodes in the triangular grid associated with the coded port and the port from small to large to obtain a distance sequence;
acquiring a neighbor triangular mesh having a common edge with the associated triangular mesh;
sequentially calculating the distance between the unprocessed grid nodes in the neighbor triangular grid and the port;
inserting the distance between the unprocessed grid nodes in the neighbor triangular grid and the port into the original distance sequence to form an updated distance sequence from small to large;
acquiring a new neighbor triangular mesh having a common edge with the neighbor triangular mesh, and returning to the step of sequentially calculating the distances between the unprocessed mesh nodes in the neighbor triangular mesh and the ports until the mesh nodes of all triangular meshes are processed, so as to obtain an updated distance sequence arranged from small to large;
continuously coding the triangular mesh nodes outside all the ports according to the distance sequence arranged from small to large; the codes with small distance are small, and the codes with large distance are large.
4. The method according to claim 1, wherein the repeatedly eliminating non-ported mesh nodes of the sparse finite element matrix by using triangle-star transformation according to the encoding of the mesh nodes of the finite element mesh to obtain the multiport network admittance matrix comprises:
and sequentially eliminating the grid nodes except the ports of the finite element grid according to the sequence of grid node codes from large to small by utilizing the triangular-star transformation to obtain the admittance matrix of the multiport network.
5. A trellis node coding system based on integrated circuit impedance network extraction, comprising:
the multilayer integrated circuit layout acquisition module is used for acquiring a multilayer integrated circuit layout; the multilayer integrated circuit layout comprises a plurality of polygons with a plurality of vertexes;
the mesh generation module is used for carrying out mesh generation on the multilayer integrated circuit layout to obtain a non-structural triangular mesh for dividing the multilayer integrated circuit layout;
the port coding module is used for sequentially and continuously coding the ports in the triangular meshes; the ports are a plurality of ports in a multi-port network impedance matrix to be calculated, which is defined by a user in advance;
a mesh node-to-port distance obtaining module, configured to obtain a distance from a mesh node in the triangular mesh to the port;
the mesh node coding module is used for coding the corresponding mesh nodes according to the distances from the mesh nodes in the triangular mesh to the ports; the distance from the mesh node to the port in the triangular mesh is positively correlated with the code of the mesh node, namely the larger the distance is, the larger the code of the mesh node is;
the finite element sparse matrix determining module is used for writing a finite element equation set of a potential field of the computing integrated circuit in a row mode according to the information of the grid nodes on the coded triangular grid to obtain a finite element sparse matrix; the grid nodes of the finite element grid are associated with the finite element sparse matrix; the codes of the grid nodes of the finite element grid are the same as the codes of the grid nodes of the triangular grid after the codes are obtained; the encoding of the finite element mesh is the same as that of the triangular mesh after encoding; the information of the grid nodes comprises voltage, current, power, conductivity, potential, dielectric constant, permeability, electric field intensity and magnetic field intensity;
the admittance matrix determination module of the multiport network is used for eliminating the non-port grid nodes of the finite element sparse matrix by repeatedly utilizing triangle-star transformation according to the codes of the grid nodes of the finite element grid to obtain a multiport network admittance matrix;
and the impedance network model extraction module is used for extracting an impedance network model according to the multiport network admittance matrix.
6. The system of claim 5, further comprising:
and the triangular mesh coding module is used for coding the triangular mesh.
7. The system of claim 5, wherein the grid node to port distance obtaining module specifically comprises:
the initialization unit is used for initializing the mesh nodes in the triangular mesh; the codes of the initialized grid nodes are all 0;
an associated triangle mesh determination unit for determining a triangle mesh associated with the encoded port from the encoded port; the triangle mesh associated with the encoded port is a triangle mesh whose triangle vertices include the port;
an unprocessed mesh node acquisition unit configured to acquire an unprocessed mesh node in the triangular mesh associated with the encoded port;
a first distance determining unit, configured to sequentially calculate distances between unprocessed mesh nodes in the triangular mesh associated with the encoded port and the port;
a distance sequence determining unit, configured to sort distances between unprocessed mesh nodes in the triangular mesh associated with the encoded port and the port from small to large to obtain a distance sequence;
a neighbor triangular mesh acquisition unit for acquiring a neighbor triangular mesh having a common edge with the associated triangular mesh;
the second distance determining unit is used for sequentially calculating the distance between the unprocessed grid nodes in the neighbor triangular grid and the port;
a distance sequence updating unit, configured to insert the distance between an unprocessed mesh node in the neighboring triangular mesh and the port into an original distance sequence to form an updated small-to-large distance sequence;
the updated distance sequence determining unit is used for acquiring a new neighbor triangular mesh which has a common edge with the neighbor triangular mesh, and continuously and sequentially calculating the distance between the unprocessed mesh nodes in the neighbor triangular mesh and the port by using the second distance determining unit until the mesh nodes of all the triangular meshes are processed, so as to obtain an updated distance sequence which is arranged from small to large;
the continuous coding unit is used for continuously coding the triangular mesh nodes outside all the ports according to the distance sequence arranged from small to large; the codes with small distance are small, and the codes with large distance are large.
8. The system of claim 5, wherein the admittance matrix determination module of the multiport network comprises:
and the admittance matrix determining unit of the multiport network is used for sequentially eliminating the grid nodes except the ports of the finite element grid according to the sequence of grid node coding from large to small by utilizing the triangular-star transformation to obtain the admittance matrix of the multiport network.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011513386.8A CN112290955B (en) | 2020-12-21 | 2020-12-21 | Grid node coding method and system based on integrated circuit impedance network extraction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011513386.8A CN112290955B (en) | 2020-12-21 | 2020-12-21 | Grid node coding method and system based on integrated circuit impedance network extraction |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112290955A CN112290955A (en) | 2021-01-29 |
| CN112290955B true CN112290955B (en) | 2021-03-30 |
Family
ID=74425935
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011513386.8A Active CN112290955B (en) | 2020-12-21 | 2020-12-21 | Grid node coding method and system based on integrated circuit impedance network extraction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112290955B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113591260B (en) * | 2021-08-19 | 2023-09-05 | 上海电机系统节能工程技术研究中心有限公司 | Method and device for generating electromagnetic field grid model of motor and electronic equipment |
| CN115600550B (en) * | 2022-11-30 | 2023-03-10 | 北京智芯仿真科技有限公司 | Accurate detection method for determining layout design defects of integrated circuit based on port resistance |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9465896B1 (en) * | 2015-03-04 | 2016-10-11 | Cadence Design Systems, Inc. | Systems and methods for testing integrated circuit designs |
| CN110096823A (en) * | 2019-05-08 | 2019-08-06 | 深圳职业技术学院 | Based on binary-coded digit integrated circuit wiring method and terminal device |
| CN111737947A (en) * | 2020-08-06 | 2020-10-02 | 北京智芯仿真科技有限公司 | Integrated circuit full-wave IBIS model extraction method and device based on field-circuit coupling |
| CN111737943A (en) * | 2020-08-06 | 2020-10-02 | 北京智芯仿真科技有限公司 | Integrated circuit IBIS model extraction method and system based on equivalent circuit model |
-
2020
- 2020-12-21 CN CN202011513386.8A patent/CN112290955B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9465896B1 (en) * | 2015-03-04 | 2016-10-11 | Cadence Design Systems, Inc. | Systems and methods for testing integrated circuit designs |
| CN110096823A (en) * | 2019-05-08 | 2019-08-06 | 深圳职业技术学院 | Based on binary-coded digit integrated circuit wiring method and terminal device |
| CN111737947A (en) * | 2020-08-06 | 2020-10-02 | 北京智芯仿真科技有限公司 | Integrated circuit full-wave IBIS model extraction method and device based on field-circuit coupling |
| CN111737943A (en) * | 2020-08-06 | 2020-10-02 | 北京智芯仿真科技有限公司 | Integrated circuit IBIS model extraction method and system based on equivalent circuit model |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112290955A (en) | 2021-01-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112257372B (en) | Method and system for extracting impedance network model of integrated circuit | |
| CN112290955B (en) | Grid node coding method and system based on integrated circuit impedance network extraction | |
| CN111737947B (en) | Integrated circuit full-wave IBIS model extraction method and device based on field-circuit coupling | |
| CN111931458B (en) | Three-dimensional integrated circuit electromagnetic field calculation method and device based on mixed order finite element | |
| TW432304B (en) | Integrated circuit design method, database apparatus for designing integrated circuit and integrated circuit design support apparatus | |
| CN111737943B (en) | Integrated circuit IBIS model extraction method and system based on equivalent circuit model | |
| CN107315848B (en) | Method and system for integrated circuit design | |
| CN115600549B (en) | Accurate detection method for determining integrated circuit layout design defects based on grid subdivision | |
| Yu et al. | Fast capacitance extraction of actual 3-D VLSI interconnects using quasi-multiple medium accelerated BEM | |
| CN114065673A (en) | Bessel integral self-adaptive segmentation method and system in integrated circuit rapid calculation | |
| CN112149336B (en) | Method and device for quickly generating three-dimensional integrated circuit electromagnetic simulation high-quality grid | |
| JP2658917B2 (en) | 3D wiring inductance calculation method | |
| CN112131774A (en) | Mixed-order finite element method and device for integrated circuit triangular prism mesh subdivision | |
| CN112287629B (en) | Impedance network model extraction method and system for very large scale integrated circuit | |
| CN114330215A (en) | Electromagnetic simulation method and device for rapidly processing layout routing of integrated circuit | |
| CN111523160B (en) | Method and apparatus for emulating an integrated circuit and computer readable medium | |
| Stine et al. | A simulation methodology for assessing the impact of spatial/pattern dependent interconnect parameter variation on circuit performance | |
| Lee et al. | Accuracy investigation of a neuromorphic machine learning system due to electromagnetic noises using PEEC model | |
| CN112989735B (en) | Iterative solution method and device for interlayer coupling of multilayer very large scale integrated circuit | |
| Kumashiro et al. | Asymptotic waveform evaluation for transient analysis of 3-D interconnect structures | |
| CN113868863A (en) | Statistical electromagnetic DGTD (differential gradient time division) calculation method for uncertain structure target | |
| CN110703038B (en) | Harmonic impedance estimation method suitable for fan access power distribution network | |
| CN114330214B (en) | Method and device for fast and high-precision calculation of integrated circuit containing routing | |
| WO2023158579A1 (en) | Integrated circuit interconnect shape optimizer | |
| CN112818633B (en) | Iterative method and device for integrated circuit current distribution dynamically applied by interlayer coupling |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |