CN112131774B - 用于集成电路三棱柱网格剖分的混合阶有限元方法及装置 - Google Patents
用于集成电路三棱柱网格剖分的混合阶有限元方法及装置 Download PDFInfo
- Publication number
- CN112131774B CN112131774B CN202011334620.0A CN202011334620A CN112131774B CN 112131774 B CN112131774 B CN 112131774B CN 202011334620 A CN202011334620 A CN 202011334620A CN 112131774 B CN112131774 B CN 112131774B
- Authority
- CN
- China
- Prior art keywords
- layer
- mesh
- polygons
- triangular
- integrated circuit
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Design And Manufacture Of Integrated Circuits (AREA)
- Semiconductor Integrated Circuits (AREA)
Abstract
本发明提供用于集成电路三棱柱网格剖分的混合阶有限元方法及装置：收集并设置集成电路版图多边形的层信息，将多边形投影到同一层并进行初始网格剖分；简化与对齐集成电路版图多边形，从xy方向对集成电路版图多边形进行三角形网格剖分和自适应细分；将三角形网格在厚度方向扩展为三棱柱网格，从z方向对三棱柱网格进行自适应细分；采用一致剖分规则实现三棱柱网格到四面体网格的剖分，基于四面体网格求解三维超大规模集成电路电磁场的初始值，根据电磁场的初始值变化速度超过和不超过预设值的四面体网格以及集成电路版图小尺度区域的四面体网格建立混合阶单元，实现混合阶有限元技术。实现快速高质量的三棱柱网格剖分；提高电磁响应的计算精度。
Description
技术领域
本发明涉及集成电路技术领域，特别涉及一种用于集成电路三棱柱网格剖分的混合阶有限元方法及装置。
背景技术
随着集成电路产业的飞速发展，传统的二维平面的封装技术已经无法满足日益增长的芯片性能和封装密度的要求，三维封装技术已经成为高密度芯片封装的一个不可逆转的趋势。三维集成电路就是将多块平面芯片从垂直方向上进行堆叠，使用金线键合或者其他方式进行互连，从而形成一种立体的封装结构。由于三维技术最大可能地利用了芯片的立体空间，以堆叠的形式使得在单位面积内能封装更多的晶体管，有效地降低了芯片的延时和功耗，在提高芯片性能的同时降低了芯片的制造成本。随着芯片特征尺寸的不断缩小，芯片上单位面积内可集成越来越多的晶体管。而为了获得更高的芯片性能，芯片的工作频率不断提升，电源电压不断降低，噪声容限也不断缩小。因而所产生的电源完整性和信号完整性问题，已经成为芯片设计中的关键性问题，因此针对三维集成电路进行电磁仿真，分析其电源完整性问题、信号完整性问题已成为迫在眉睫需要解决的问题。
早期，针对结构简单的集成电路对其电源完整性和信号完整性问题主要采用简化的传输线法或有限差分法，以上方法不需要针对复杂结构的集成电路进行复杂的网格剖分，计算速度快，但由于对集成电路版图进行了很大的近似，在处理早期结构简单、规则的版图计算结果是准确的，但对于近年来结构越来越复杂版图则产生不可比对的误差。近年来，针对带有结构复杂的集成电路版图的多层集成电路采用基于场的方法计算，对各层集成电路版图采用二维非结构网格剖分，在网格剖分的过程中考虑了版图的复杂结构，因此其计算结果更为准确，但这种方法采用的是二维网格对计算场域进行离散，其假设的前提是集成电路版图尺寸远远大于集成电路金属层的厚度以及金属层之间的介质层的厚度，对于目前正在发展的三维集成电路来说，其电源、信号不光在各层的平面结构传递，还在垂直方向产生大量的互连，且集成电路电源部分的供电金属线的宽度越来越窄，已经可以和金属层或介质层厚度相比拟，此时，这个假设已经显现出越来越大的误差，因此有必要直接采用三维电磁场数值计算方法如三维有限元法直接对整个三维集成电路进行电磁仿真。然而，由于这种发展的三维集成电路结构上具有非常明显的多尺度特征，其尺度范围从互连线宽的纳米级到电源层版图区域的厘米级，层间距、过孔尺寸也从微米级到纳米级，如果直接对这样的多尺度的复杂的三维集成电路进行三维非结构的四面体网格剖分，将花费大量的CPU时间，并且，可能会在小尺度区域产生过于密集的网格，导致最终产生数量庞大的网格，对三维非结构的四面体网格剖分技术以及由此产生的超大规模稀疏矩阵的求解技术将是一个极大的挑战。
为了提高三维超大规模集成电路电磁仿真的计算精度，采用有限元法计算时通常采用两种方法：一是对网格剖分进行加密，生成更加密集的网格；二是增加有限元网格单元的阶数，采用更高阶的插值函数逼近待计算的场在单元内的分布。但对于三维多尺度结构的超大规模集成电路来说，对网格加密通常导致网格密度在大量小尺度区域异常增加，一方面使得网格密度分布不合理，增加的网格量并没有起到增加计算精度的作用，导致计算效率低下；同时，由于对网格加密导致大量区域网格密度异常增加，使得总体网格数量有上千万的增加，对后面稀疏矩阵的求解提出挑战。因此，对三维超大规模集成电路电磁仿真，常采用增加有限元的阶数提高计算精度，如将一阶单元提高到二阶单元，但如果将所有单元都提高到二阶单元，同样也会导致计算的未知量急剧增大，计算效率仍然很低。为避免将所有单元都提高到二阶单元导致计算的未知量急剧增大，提出采用混合阶有限元，基于后验误差估计对误差大的区域采用二阶单元，其他地方采用一阶单元，一阶单元与二阶单元相邻的区域采用过渡单元，但目前针对三棱柱单元，仍未找到过渡单元的实现方法。
本发明提出一种用于集成电路三棱柱网格剖分的混合阶有限元方法及装置。实现了对三维多尺度结构超大规模集成电路进行快速高质量的三棱柱网格剖分，又实现了基于后验误差估计选择性的提高某些单元的阶数，使得三维多尺度结构超大规模集成电路的电磁响应计算精度大幅提高。
发明内容
（一）发明目的
为克服上述现有技术存在的至少一种缺陷，本发明提供了一种用于集成电路三棱柱网格剖分的混合阶有限元方法及装置，能够实现对三维多尺度结构超大规模集成电路进行快速高质量的三棱柱网格剖分，能够实现基于后验误差估计选择性的提高某些单元的阶数，使得三维多尺度结构超大规模集成电路的电磁响应计算精度大幅提高。
（二）技术方案
作为本发明的第一方面，本发明公开了一种用于集成电路三棱柱网格剖分的混合阶有限元方法，包括以下步骤：
步骤1，收集并设置所有层的集成电路版图多边形的层信息，将所有层的所述多边形投影到同一层并进行初始网格剖分；简化与对齐所述集成电路版图多边形，从xy平面方向对所述集成电路版图多边形进行三角形网格剖分和自适应细分；
步骤2，将三角形网格在厚度方向扩展为三棱柱网格，从z方向对三棱柱网格进行自适应细分；
步骤3，采用一致剖分规则实现所述三棱柱网格到四面体网格的网格剖分，基于所述四面体网格求解三维超大规模集成电路电磁场的初始值，根据所述电磁场的初始值变化速度超过和不超过预设值的所述四面体网格以及所述集成电路版图小尺度区域的所述四面体网格，建立混合阶单元，实现混合阶有限元技术，使得集成电路版图的电磁响应计算精度得到提高。
一种可能的实施方式中，所述步骤1包括：采用边交换法恢复所有所述多边形的边，在不同层的所述多边形的边的交点插入新的网格节点；
所述收集并设置所有层的集成电路版图多边形的层信息，将所有层的所述多边形投影到同一层并进行初始网格剖分；以及，所述采用边交换法恢复所有所述多边形的边，在不同层的所述多边形的边的交点插入新的网格节点；包括：
步骤1.1，获取多层所述集成电路版图包含多个顶点的多个所述多边形，增加两个能覆盖所有层的所述集成电路版图多边形的矩形多边形，分别作为顶层空气层和底层空气层的所述多边形；
步骤1.2，设置各层所述多边形的层信息，并用二进制数字代表各自的层；
步骤1.3，将各层的多个所述多边形垂直投影到同一层，根据Delaunay三角剖分算法形成以多边形顶点为网格节点的Delaunay三角形网格，其中，所述多边形的各个边包含预先设定的所在多边形的所述多边形信息和所在层的层信息；
步骤1.4，合并投影后重合的多个所述多边形的边的所述多边形信息和所述多边形的层信息；
步骤1.5，根据所述边交换法将所述Delaunay三角形网格对齐到多个所述多边形的各个边，同时计算多个不同层的所述多边形的边的交点并将所述交点新增为所述多边形的顶点和所述Delaunay三角形网格的节点，形成第一三角形网格。
一种可能的实施方式中，在所述步骤1中，所述简化与对齐所述集成电路版图多边形，包括：
对投影到同一层的所述多边形的边进行简化与对齐，并对简化与对齐后的所述多边形及整个集成电路版图区域进行三角形网格剖分；
所述对投影到同一层的所述多边形的边进行简化与对齐，并对简化与对齐后的所述多边形及整个集成电路版图区域进行三角形网格剖分；以及，所述从xy平面方向对所述集成电路版图多边形进行三角形网格剖分和自适应细分；包括：
步骤1.6，基于所述第一三角形网格，在每个所述多边形的内外分别形成夹住所述多边形的内外辅助多边形，并通过设定的距离阈值控制所述内外辅助多边形与所述多边形的距离；
步骤1.7，对落在所述内外辅助多边形之间的各层所述多边形的边进行对齐和简化处理，并根据各个所述多边形的边所包含的多边形编号信息将投影到同一层的多层所述多边形还原到各层中，同时更新所述三角形网格及其层信息，形成第二三角形网格；
步骤1.8，基于网格尺寸和单元质量控制技术对所述第二三角形网格进行自适应网格细分，形成更新的所述第二三角形网格。
一种可能的实施方式中，所述步骤2包括：
步骤2.1：根据包括金属层和介质层的每层纵向上的排列顺序，记录所述集成电路版图每层唯一的区域编号；
步骤2.2：设定纵向上的坐标原点，根据所述集成电路版图每层的厚度及区域编号顺序，记录每层的位置信息；
步骤2.3：基于所述步骤1获得的所述第二三角形网格，将所有所述三角形网格在厚度方向扩展为一个所述三棱柱，所述三棱柱的上下底面三角形为所述第二三角形网格对应的所述三角形，上下底面位置为层界面最大值和最小值，从而形成第一三棱柱网格；
步骤2.4：基于所述步骤2.3获得的所述第一三棱柱网格，将所有所述三棱柱沿z方向按zinterface进行分割，每个所述三棱柱被分割为N-1个所述三棱柱，形成第二三棱柱网格，其中N为所述集成电路版图包括上下空气层的层界面个数，zinterface为层界面的位置。
步骤2.5：基于所述步骤1.7获得的所述第二三角形网格，根据所述三角形网格的层信息，还原所述三角形网格被哪些层界面的所述多边形共用；
步骤2.6：根据所述步骤2.5确定的所述三角形网格被哪些层界面的所述多边形共用，以此推算出其依次包含的区域，进而确定所述第二三棱柱网格中每个所述三棱柱所在的区域；
步骤2.7，实现所述三棱柱网格在z方向统一的由细到粗的渐进变化的细分过程。
一种可能的实施方式中，在所述步骤3中，所述采用一致剖分规则实现所述三棱柱网格到四面体网格的网格剖分，包括：
将所有所述三棱柱网格的顶点统一编号；
对每个所述三棱柱网格，由上层所述编号较小的点开始，向上层所述编号较大的点对应的下层点连接，将所述三棱柱网格剖分为所述四面体网格；或者，由上层所述编号较大的点开始，向上层所述编号较小的点对应的下层点连接，将所述三棱柱网格剖分为所述四面体网格。
作为本发明的第二方面，本发明公开了一种用于集成电路三棱柱网格剖分的混合阶有限元装置，包括：
第一模块，用于收集并设置所有层的集成电路版图多边形的层信息，将所有层的所述多边形投影到同一层并进行初始网格剖分；简化与对齐所述集成电路版图多边形，从xy平面方向对所述集成电路版图多边形进行三角形网格剖分和自适应细分；
第二模块，用于将三角形网格在厚度方向扩展为三棱柱网格，从z方向对三棱柱网格进行自适应细分；
第三模块，用于采用一致剖分规则实现所述三棱柱网格到四面体网格的网格剖分，基于所述四面体网格求解三维超大规模集成电路电磁场的初始值，根据所述电磁场的初始值变化速度超过和不超过预设值的所述四面体网格以及所述集成电路版图小尺度区域的所述四面体网格，建立混合阶单元，实现混合阶有限元技术，使得集成电路版图的电磁响应计算精度得到提高。
一种可能的实施方式中，所述第一模块包括：采用边交换法恢复所有所述多边形的边，在不同层的所述多边形的边的交点插入新的网格节点；
所述收集并设置所有层的集成电路版图多边形的层信息，将所有层的所述多边形投影到同一层并进行初始网格剖分；以及，所述采用边交换法恢复所有所述多边形的边，在不同层的所述多边形的边的交点插入新的网格节点；包括：
第一单元，用于获取多层所述集成电路版图包含多个顶点的多个所述多边形，增加两个能覆盖所有层所述集成电路版图多边形的矩形多边形，分别作为顶层空气层和底层空气层的所述多边形；
第二单元，用于设置各层所述多边形的层信息，并用二进制数字代表各自的层；
第三单元，用于将各层的多个所述多边形垂直投影到同一层，根据Delaunay三角剖分算法形成以多边形顶点为网格节点的Delaunay三角形网格，其中，所述多边形的各个边包含预先设定的所在多边形的所述多边形信息和所在层的层信息；
第四单元，用于合并投影后重合的多个所述多边形的边的所述多边形信息和所述多边形的层信息；
第五单元，用于根据所述边交换法将所述Delaunay三角形网格对齐到多个所述多边形的各个边，同时计算多个不同层的所述多边形的边的交点并将所述交点新增为所述多边形的顶点和所述Delaunay三角形网格的节点，形成第一三角形网格。
一种可能的实施方式中，在所述第一模块中，所述简化与对齐所述集成电路版图多边形，包括：
对投影到同一层的所述多边形的边进行简化与对齐，并对简化与对齐后的所述多边形及整个集成电路版图区域进行三角形网格剖分；
所述对投影到同一层的所述多边形的边进行简化与对齐，并对简化与对齐后的所述多边形及整个集成电路版图区域进行三角形网格剖分；以及，所述从xy平面方向对所述集成电路版图多边形进行三角形网格剖分和自适应细分；包括：
第六单元，用于基于所述第一三角形网格，在每个所述多边形的内外分别形成夹住所述多边形的内外辅助多边形，并通过设定的距离阈值控制所述内外辅助多边形与所述多边形的距离；
第七单元，用于对落在所述内外辅助多边形之间的各层所述多边形的边进行对齐和简化处理，并根据各个所述多边形的边所包含的多边形编号信息将投影到同一层的多层所述多边形还原到各层中，同时更新所述三角形网格及其层信息，形成第二三角形网格；
第八单元，用于基于网格尺寸和单元质量控制技术对所述第二三角形网格进行自适应网格细分，形成更新的所述第二三角形网格。
一种可能的实施方式中，所述第二模块，
用于根据包括金属层和介质层的每层纵向上的排列顺序，记录所述集成电路版图每层唯一的区域编号；
用于设定纵向上的坐标原点，根据所述集成电路版图每层的厚度及区域编号顺序，记录每层的位置信息；
用于基于所述第一模块获得的所述第二三角形网格，将所有所述三角形网格在厚度方向扩展为一个所述三棱柱，所述三棱柱的上下底面三角形为所述第二三角形网格对应的所述三角形，上下底面位置为层界面最大值和最小值，从而形成第一三棱柱网格；
用于所述第三单元获得的所述第一三棱柱网格，将所有所述三棱柱沿z方向按zinterface进行分割，每个所述三棱柱被分割为N-1个所述三棱柱，形成第二三棱柱网格，其中N为所述集成电路版图包括上下空气层的层界面个数，zinterface为层界面的位置；
用于基于所述第七单元获得的所述第二三角形网格，根据所述三角形网格的层信息，还原所述三角形网格被哪些层界面的所述多边形共用；
用于根据确定的所述三角形网格被哪些层界面的所述多边形共用，以此推算出其依次包含的区域，进而确定所述第二三棱柱网格中每个所述三棱柱所在的区域；
用于实现所述三棱柱网格在z方向统一的由细到粗的渐进变化的细分过程。
一种可能的实施方式中，在所述第三模块中，所述采用一致剖分规则实现所述三棱柱网格到四面体网格的网格剖分，包括：
将所有所述三棱柱网格的顶点统一编号；
对每个所述三棱柱网格，由上层所述编号较小的点开始，向上层所述编号较大的点对应的下层点连接，将所述三棱柱网格剖分为所述四面体网格；或者，由上层所述编号较大的点开始，向上层所述编号较小的点对应的下层点连接，将所述三棱柱网格剖分为所述四面体网格。
（三）有益效果
本发明提供的一种用于集成电路三棱柱网格剖分的混合阶有限元方法及装置，通过步骤1收集并设置所有层的集成电路版图多边形的层信息，将所有层的多边形投影到同一层并进行初始网格剖分；简化与对齐所述集成电路版图多边形，从xy平面方向对所述集成电路版图多边形进行所述三角形网格剖分和自适应细分；通过步骤2将三角形网格在厚度方向扩展为三棱柱网格，从z方向对所述三棱柱网格进行自适应细分；步骤3，采用一致剖分规则实现所述三棱柱网格到四面体网格的网格剖分，基于所述四面体网格求解三维超大规模集成电路电磁场的初始值，根据所述电磁场的初始值变化速度超过和不超过预设值的所述四面体网格以及所述集成电路版图小尺度区域的所述四面体网格，建立混合阶单元，实现混合阶有限元技术，使得集成电路版图的电磁响应计算精度得到提高。实现了对三维多尺度结构超大规模集成电路进行快速高质量的三棱柱网格剖分，又实现了基于后验误差估计选择性的提高某些单元的阶数，使得三维多尺度结构超大规模集成电路的电磁响应计算精度大幅提高。
附图说明
以下参考附图描述的实施例是示例性的，旨在用于解释和说明本发明，而不能理解为对本发明的保护范围的限制。
图1是本发明提供的一种用于集成电路三棱柱网格剖分的混合阶有限元方法的流程图。
图2是本发明提供的一种用于集成电路三棱柱网格剖分的混合阶有限元方法及装置的简单的多层集成电路版图金属层与介质层示意图。
图3是本发明提供的用于集成电路三棱柱网格剖分的混合阶有限元方法及装置的多层集成电路版图的区域编号与层界面位置关系示意图。
图4(a)~4(f)是本发明提供的用于集成电路三棱柱网格剖分的混合阶有限元方法及装置的三棱柱网格到四面体网格的剖分示意图。
图5是本发明提供的用于集成电路三棱柱网格剖分的混合阶有限元方法及装置的相邻三棱柱网格到四面体网格错误的剖分示意图。
图6是本发明提供的用于集成电路三棱柱网格剖分的混合阶有限元方法及装置的相邻三棱柱网格示意图。
图7（a）是本发明提供的用于集成电路三棱柱网格剖分的混合阶有限元方法及装置的采用规则1的剖分结果示意图。
图7(b)是本发明提供的用于集成电路三棱柱网格剖分的混合阶有限元方法及装置的采用规则2的剖分结果示意图。
图8是本发明提供的用于集成电路三棱柱网格剖分的混合阶有限元装置的结构示意图。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行更加详细的描述。
需要说明的是：在附图中，自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本发明一部分实施例，而不是全部的实施例，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，均仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明保护范围的限制。
下面参考图1-7详细描述本发明提供的一种用于集成电路三棱柱网格剖分的混合阶有限元方法的第一实施例。如图1-7所示，本实施例提供的混合阶有限元方法主要包括有：步骤1、步骤2和步骤3。
本发明直接针对尺度范围为厘米级到纳米级的复杂的三维集成电路进行三维非结构的网格剖分，但与传统三维非结构的四面体网格剖分不同，本专利将基于平行平板场域识别技术，对这样的多尺度的复杂的三维集成电路进行平面非结构的三维三棱柱网格剖分。
本专利在形成的第二三角形网格和第二三棱柱网格的基础上，分别在xy平面方向和z方向进行自适应网格细分，以此形成三维超大规模集成电路的自适应细分的三棱柱网格。
本发明一方面采用快速网格剖分技术，对三维多尺度结构超大规模集成电路（即集成电路版图）进行快速高质量的三棱柱网格剖分，避免直接采用三维非结构四面体网格剖分导致网格剖分进行不下去的问题；另一方面，为提高其电磁响应的计算精度，采用混合阶有限元技术，基于后验误差估计对误差大的区域采用二阶单元，其他地方采用一阶单元，一阶单元与二阶单元相邻的区域采用过渡单元。针对目前三棱柱单元（即三棱柱网格）仍未找到过渡单元实现方法的问题，本发明提出将剖分的三棱柱单元采用一致的规则剖分成四面体网格，再针对四面体网格采用混合阶有限元的技术。该技术既实现了对三维多尺度结构超大规模集成电路进行快速高质量的三棱柱网格剖分，又实现了基于后验误差估计选择性的提高某些单元的阶数，使得三维多尺度结构超大规模集成电路的电磁响应计算精度大幅提高。
首先通过步骤1收集并设置所有层的集成电路版图多边形的层信息，将所有层的所述多边形投影到同一层并进行初始网格剖分；简化与对齐所述集成电路版图多边形，从xy平面方向对所述集成电路版图多边形进行所述三角形网格剖分和自适应细分；进而通过步骤2将三角形网格在厚度方向扩展为三棱柱网格，从z方向对所述三棱柱网格进行自适应细分。
从上述可以看出，本专利实现了通过二维网格剖分的时间复杂度对多尺度的复杂的三维集成电路进行了平面方向非结构的三维三棱柱网格剖分，并且采取分别在xy平面方向和z方向的自适应细分的方式，完全避免在小尺度区域产生过于密集的三棱柱网格。
通过步骤3采用一致剖分规则实现所述三棱柱网格到四面体网格的网格剖分，基于所述四面体网格求解三维超大规模集成电路电磁场的初始值，根据所述电磁场的初始值变化速度超过和不超过预设值的所述四面体网格以及所述集成电路版图小尺度区域的所述四面体网格，建立混合阶单元，实现混合阶有限元技术，使得集成电路版图的电磁响应计算精度得到提高。
所述步骤3包括：
步骤3.1，基于四面体网格，建立一阶单元及其有限元刚度矩阵，并根据场路耦合求解三维超大规模集成电路电磁场的初始值；
步骤3.2，基于四面体网格，将所述电磁场的初始值变化速度不超过预设值的四面体网格标记为一阶单元，将所述电磁场的初始值变化速度超过预设值的四面体网格以及集成电路版图小尺度区域的四面体网格标记为二阶单元，根据所述一阶单元和二阶单元建立混合阶单元；其中，集成电路版图小尺度区域包括集成电路版图中各层之间的过孔处、各版图之间的互连线处以及用户设定的端口处；
步骤3.3，根据混合阶单元建立有限元刚度矩阵，计算所述三维超大规模集成电路电磁场。
所述步骤1可以为三维超大规模集成电路版图多边形收集与带层信息的三角形网格剖分。
其中，所述步骤1包括：采用边交换法恢复所有所述多边形的边，在不同层的所述多边形的边的交点插入新的网格节点；
所述收集并设置所有层的集成电路版图多边形的层信息，将所有层的所述多边形投影到同一层并进行初始网格剖分；以及，所述采用边交换法恢复所有所述多边形的边，在不同层的所述多边形的边的交点插入新的网格节点(可以为先试图恢复多边形的边，如果两个相交的边都为多边形的边（即无法恢复时），则在交点插入新的网格节点)；包括：
步骤1.1，获取多层所述集成电路版图包含多个顶点的多个所述多边形，除此之外在顶层之上和底层之下再增加两个能覆盖所有层的所述集成电路版图多边形的矩形多边形，分别作为顶层空气层和底层空气层的所述多边形；
步骤1.2，设置各层所述多边形的层信息，这种多边形的层信息应该既能唯一代表多边形所在的层，又能进行任意层之间的合并等布尔操作；如对于图2所示的多层集成电路版图，其为有4个金属层和3个介质层的集成电路模型，则共有2个空气层，4个金属层及3个介质层，即共有9层，为10个层界面；并用二进制数字代表各自的层，如果多边形位于哪个层界面，则该层界面所在的二进制的位为1。由于需要考虑金属层的厚度，因此位于金属层上的多边形的层信息同时包含金属层上下层界面。对图2所示的多层集成电路版图，可设置多边形的层信息如表1所示。
表1 图2所示多层集成电路各层多边形的层信息
步骤1.3，将各层的多个所述多边形垂直投影到同一层，根据Delaunay三角剖分算法形成以多边形顶点为网格节点的Delaunay三角形网格，其中，所述多边形的各个边包含预先设定的所在多边形的所述多边形信息和所在层的层信息；
步骤1.4，合并投影后重合的多个所述多边形的边的所述多边形信息和所述多边形的层信息。
例如图2所示的多层集成电路版图，某个金属层1和金属层2上的多边形重合，进行合并后，这两个多边形合并为一个多边形，其层信息为：6+24=30。
步骤1.5，根据所述边交换法将所述Delaunay三角形网格对齐到多个所述多边形的各个边，同时计算多个不同层的所述多边形的边的交点并将所述交点新增为所述多边形的顶点和所述Delaunay三角形网格的节点，形成第一三角形网格；
其中，在所述步骤1中，所述简化与对齐所述集成电路版图多边形，包括：
对投影到同一层的所述多边形的边进行简化与对齐，并对简化与对齐后的所述多边形及整个集成电路版图区域进行三角形网格剖分；
所述对投影到同一层的所述多边形的边进行简化与对齐，并对简化与对齐后的所述多边形及整个集成电路版图区域进行三角形网格剖分；以及，所述从xy平面方向对所述集成电路版图多边形进行三角形网格剖分和自适应细分；包括：
步骤1.6，基于所述第一三角形网格，在每个所述多边形P的内外分别形成夹住所述多边形P的内外辅助多边形P 0与P 9，并通过设定的距离阈值控制所述内外辅助多边形与所述多边形P的距离；
步骤1.7，对落在所述内外辅助多边形之间的各层所述多边形的边进行对齐和简化处理，并根据各个所述多边形的边所包含的多边形编号信息将投影到同一层的多层所述多边形还原到各层中，同时更新所述三角形网格及其层信息，形成第二三角形网格；
步骤1.8，直接基于网格尺寸和单元质量控制技术对所述包含对齐简化后的多边形的第二三角形网格进行自适应网格细分，形成更新的所述第二三角形网格。
所述步骤2可以为三维超大规模集成电路进行快速的三棱柱网格剖分。
上述步骤1，已经将所有层（包括顶层空气层和底层空气层）的集成电路版图多边形投影到同一层，并且经过对齐和简化处理，再在此基础上形成经过自适应细分的第二三角形网格，步骤2将在以上步骤1的基础上，直接快速形成三维多层超大规模集成电路的三棱柱网格剖分。
基于以上更新的第二三角形网格，再完成下述步骤S2.1~S2.6，即完成了根据层界面对三棱柱在z方向进行初步分割的操作。
其中，所述步骤2包括：
步骤2.1：根据包括金属层和介质层的每层纵向上的排列顺序（从上到下或从下到上），记录所述集成电路版图每层唯一的区域编号；例如，对于图2所示的多层集成电路版图，其共有2个空气层，4个金属层及3个介质层，可设置其区域编号如表2所示；
表2 图2所示多层集成电路各层区域编号信息
步骤2.2：设定纵向上的坐标原点，根据所述集成电路版图每层的厚度及区域编号顺序，记录每层的位置信息，包括顶层空气层和底层空气层、金属层以及介质层，例如，对于图2所示的集成电路版图，其为有4个金属层和3个介质层的集成电路模型，则共有2个空气层，4个金属层及3个介质层，即共有9层，为10个层界面，记录这些层界面的位置，记为zinterface={z1, z2, …, z10}，其中z1~z10如图3所示，其排列顺序与步骤2.1中的区域编号顺序一致；
步骤2.3：基于所述步骤1获得的第二三角形网格，将所有所述三角形网格在厚度方向扩展为一个所述三棱柱，所述三棱柱的上下底面三角形为所述第二三角形网格对应的所述三角形，上下底面位置为层界面最大值和最小值，从而形成第一三棱柱网格；由于步骤1.1在所有层集成电路版图多边形的基础上增加了覆盖所有层集成电路版图多边形的矩形多边形作为顶层空气层和底层空气层的多边形，所以第二三角形网格的所有三角形的层均包含顶层层界面（即最上层层界面）和底层层界面（即最下层层界面），因此，所有三角形扩展的三棱柱上下底面位置相同，均为顶层层界面（即最上层层界面）和底层层界面（即最下层层界面）所在的位置。
步骤2.4：基于所述步骤2.3获得的所述第一三棱柱网格，将所有所述三棱柱沿z方向（即在厚度方向）按zinterface进行分割，每个所述三棱柱被分割为N-1个所述三棱柱，形成第二三棱柱网格，其中N为所述集成电路版图包括上下空气层（即顶层空气层和底层空气层）的层界面个数，zinterface为层界面的位置。
步骤2.5：基于所述步骤1.7获得的所述第二三角形网格，根据所述三角形网格的层信息，还原所述三角形网格被哪些层界面的多边形共用；如对于图2所示的集成电路示意图，假设第二三角形网格中的某个三角形的层编码为543，则543=29+23+24+21+22+20，即该三角形网格被第0层（即底层空气层，用编码式中的20表示）、金属层1（用编码式中的21+22表示）、金属层2（用编码式中的23+24表示）及第9层（即顶层空气层，用编码式中的29表示）共用；
步骤2.6：根据所述步骤2.5确定的所述三角形网格被哪些层界面的多边形共用，以此推算出其依次包含的区域，进而确定所述第二三棱柱网格中每个所述三棱柱所在的区域。如对于图2所示的集成电路示意图，假设第二三角形网格中的某个三角形的层编码为543，其被第0层（即底层空气层）、金属层1、金属层2及第9层（即顶层空气层）共用，则依据图3所示的多层集成电路版图的区域编号与层界面位置关系示意图，可确定该三角形网格对应的第二三棱柱网格中每个三棱柱所在的区域从下到上依次为：区域1、区域2、区域3、区域4、区域5、区域7、区域7、区域9和区域9。
可通过步骤2.7实现所述三棱柱网格在z方向统一的由细到粗的渐进变化的细分过程。
通常情况下，空气层厚度远大于介质层厚度，介质层厚度远大于金属层厚度，显然在z方向不进行网格细分是不合理的，为此本专利提出一种z方向渐变的细分方法。其思路在于，设置一个小尺寸到大尺寸的渐变过渡的近似比例r 0，使得z方向的尺寸从薄层到厚层是渐进变化的。针对本专利讨论的三维超大规模集成电路模型，薄层与厚层的位置关系可归结为以下两种情况：
1) 一边为薄层，一边为厚层：介质层/金属层-空气层
设薄层上三棱柱网格z方向上的长度为t 1, 厚层厚度为t 2，则应找到这样一个渐变比例r，使得
式中n为厚层上细分的段数，n应满足以下公式
显然，如果
则r=r 0，否则，可通过二分法计算r:
首先，设置n=1，逐渐增大n，使得n满足式（2），
然后，基于计算的n，设r left=r 0, 
反复进行以下计算：
Step 1: 设
r=(r left+ r right )/2
Step 2: 计算
如果
设
r
right
=r
转入step 1，否则，如果
设
r
left
=r
转入step 1，否则，r即为计算的正确结果，结束。式中
为设置的值为足够小的误差精度，如取
。
一种实施例为：设最上层金属层上等分3份，即将金属层在z方向分为3份，每个三棱柱网格z方向长度为5 nm，空气层厚度为1 cm，设置薄层到厚层的渐变过渡比例为r 0=1.5。经计算，可得厚层（空气层）上细分的段数n=33，r=1.50143756，三棱柱网格z方向长度由小到大依次如表3所示。
表3 金属层过渡到空气层时空气层的细分网格的长度渐进变化
| 层编号 | 长度（m） |
| 0 | 5E-9 |
| 1 | 7.5071880E-9 |
| 2 | 1.1271574E-8 |
| 3 | 1.6923564E-8 |
| 4 | 2.5409676E-8 |
| 5 | 3.8151041E-8 |
| 6 | 5.7281407E-8 |
| 7 | 8.6004456E-8 |
| 8 | 1.2913032E-7 |
| 9 | 1.9388112E-7 |
| 10 | 2.9110039E-7 |
| 11 | 4.3706907E-7 |
| 12 | 6.5623192E-7 |
| 13 | 9.8529125E-7 |
| 14 | 1.4793533E-6 |
| 15 | 2.2211566E-6 |
| 16 | 3.3349280E-6 |
| 17 | 5.0071862E-6 |
| 18 | 7.5179775E-6 |
| 19 | 1.1287774E-5 |
| 20 | 1.6947888E-5 |
| 21 | 2.5446195E-5 |
| 22 | 3.8205874E-5 |
| 23 | 5.7363734E-5 |
| 24 | 8.6128066E-5 |
| 25 | 1.2931591E-4 |
| 26 | 1.9415977E-4 |
| 27 | 2.9151877E-4 |
| 28 | 4.3769724E-4 |
| 29 | 6.5717508E-4 |
| 30 | 9.8670735E-4 |
| 31 | 1.4814795E-3 |
| 32 | 2.2243490E-3 |
| 33 | 3.3397211E-3 |
其中，E-9指的是10的-9次方，E-8、E-7等以此类推。
2) 上下为薄层，中间为厚层：金属层-介质层-金属层
设上下薄层（金属层）上三棱柱网格z方向上的长度为t 1,中间厚层（介质层）厚度为t 2，则应找到这样一个渐变比例r，使得
式中2n-1为厚层上细分的段数，n应满足以下公式
显然，如果
则r=r 0，否则，可通过二分法计算r:
首先，设置n=1，逐渐增大n，使得n满足式（5），
然后，基于计算的n，设r left=r 0, 
反复进行以下计算：
Step 1: 设
r=(r left+ r right )/2
Step 2: 计算
如果
设
r
right
=r
转入step 1，否则，如果
设
r
left
=r
转入step 1，否则，r即为计算的正确结果，结束。式中
为设置的值为足够小的误差精度，如取
。
另一种实施例为：设上下层金属层上等分3份，即将金属层在z方向分为3份，每个三棱柱网格z方向长度为5 nm，中间介质层厚度为1 μm，设置薄层到厚层的渐变过渡比例为r 0=1.5。经计算，可得厚层（空气层）上细分的段数2n-1=17，r=1.51561265，三棱柱网格z方向长度由小到大依次如表4所示。
表4 金属层-介质层-金属层过渡时介质层的细分网格的长度渐进变化
以上三棱柱网格在z方向统一的由细到粗的渐进变化的细分，一方面保证了三棱柱网格的渐进变化满足场的渐进变化规律，使得电磁仿真的计算结果更加准确，另一方面，这种自适应细分保证了通过二维非结构网格自适应细分的速度来实现三维网格自适应细分的速度，且不会导致局部过于密集的三棱柱网格。
如图4所示，将三棱柱网格剖分为四面体网格有以下6种剖分方法：
图4(a)，三棱柱网格123456被分为四面体网格1234, 2564和2634；
图4(b)，三棱柱网格123456被分为四面体网格1234, 2534和5634；
图4(c)，三棱柱网格123456被分为四面体网格1235, 1645和1365；
图4(d)，三棱柱网格123456被分为四面体网格1235, 1345和3645；
图4(e)，三棱柱网格123456被分为四面体网格1236, 1526和1456；
图4(f)，三棱柱网格123456被分为四面体网格1236, 1426和2456。
将三棱柱网格剖分为四面体网格时，如果直接采用图4中的一种剖分方法进行剖分，可能因为相邻三棱柱网格的剖分不一致，导致公共面的剖分出现交叉的问题，如图5所示，这将使得公共面的定义出现二义性。
图5中，三棱柱网格123456被分为四面体网格1234, 2564和2634，三棱柱网格327658被分为四面体网格3278, 2358和3658，则四面体网格2564, 2634, 2358和3658在公共面2563出现了交叉，剖分错误。
为避免这种错误发生，本发明定义一种一致剖分规则，规则如下：
将所有所述三棱柱网格的顶点统一编号；
规则1：对每个所述三棱柱网格，由上层所述编号较小的点开始，向上层所述编号较大的点对应的下层点连接，将所述三棱柱网格剖分为3个所述四面体网格；或者，规则2：对每个所述三棱柱网格，由上层所述编号较大的点开始，向上层所述编号较小的点对应的下层点连接，将所述三棱柱网格剖分为3个所述四面体网格。
对采用快速三棱柱网格生成技术生成的三维超大规模集成电路的高质量网格，一致的采样规则1或规则2将所有三棱柱网格剖分成四面体网格，即可实现三棱柱网格到四面体网格的网格剖分。
如图6所示的相邻三棱柱网格123456和327658，采用规则1的剖分结果如图7(a)所示，剖分结果为：4132, 4523, 4536, 5237, 5673和5687。采用规则2的剖分结果如图7(b)所示，剖分结果为：6321, 6251, 6541, 8723, 8632和8256。
将三棱柱网格完全剖分为四面体网格后，即可针对四面体网格采用混合阶有限元的技术，实现基于后验误差估计选择性的提高某些单元的阶数，使得三维多尺度结构超大规模集成电路的电磁响应计算精度大幅提高。
本发明通过步骤1收集并设置所有层的集成电路版图多边形的层信息，将所有层的多边形投影到同一层并进行初始网格剖分；简化与对齐所述集成电路版图多边形，从xy平面方向对所述集成电路版图多边形进行所述三角形网格剖分和自适应细分；通过步骤2将三角形网格在厚度方向扩展为三棱柱网格，从z方向对所述三棱柱网格进行自适应细分；步骤3采用一致剖分规则实现所述三棱柱网格到四面体网格的网格剖分，基于所述四面体网格求解三维超大规模集成电路电磁场的初始值，根据所述电磁场的初始值变化速度超过和不超过预设值的所述四面体网格以及所述集成电路版图小尺度区域的所述四面体网格，建立混合阶单元，实现混合阶有限元技术，使得集成电路版图的电磁响应计算精度得到提高。实现了对三维多尺度结构超大规模集成电路进行快速高质量的三棱柱网格剖分，又实现了基于后验误差估计选择性的提高某些单元的阶数，使得三维多尺度结构超大规模集成电路的电磁响应计算精度大幅提高。
下面参考图2-8详细描述本发明提供的一种用于集成电路三棱柱网格剖分的混合阶有限元装置的第一实施例。如图2-8所示，本实施例提供的混合阶有限元装置主要包括有：第一模块、第二模块和第三模块。
第一模块，用于收集并设置所有层的集成电路版图多边形的层信息，将所有层的所述多边形投影到同一层并进行初始网格剖分；简化与对齐所述集成电路版图多边形，从xy平面方向对所述集成电路版图多边形进行所述三角形网格剖分和自适应细分；
第二模块，用于将三角形网格在厚度方向扩展为三棱柱网格，从z方向对所述三棱柱网格进行自适应细分。
从上述可以看出，本专利实现了通过二维网格剖分的时间复杂度对多尺度的复杂的三维集成电路进行平面方向非结构的三维三棱柱网格剖分，并且采取分别在xy平面方向和z方向的自适应细分的方式，完全避免了在小尺度区域产生过于密集的三棱柱网格。
第三模块，用于采用一致剖分规则实现所述三棱柱网格到四面体网格的网格剖分，基于所述四面体网格求解三维超大规模集成电路电磁场的初始值，根据所述电磁场的初始值变化速度超过和不超过预设值的所述四面体网格以及所述集成电路版图小尺度区域的所述四面体网格，建立混合阶单元，实现混合阶有限元技术，使得集成电路版图的电磁响应计算精度得到提高。
所述第三模块：
用于基于四面体网格，建立一阶单元及其有限元刚度矩阵，并根据场路耦合求解三维超大规模集成电路电磁场的初始值；
用于基于四面体网格，将所述电磁场的初始值变化速度不超过预设值的四面体网格标记为一阶单元，将所述电磁场的初始值变化速度超过预设值的四面体网格以及集成电路版图小尺度区域的四面体网格标记为二阶单元，根据所述一阶单元和二阶单元建立混合阶单元；其中，集成电路版图小尺度区域包括集成电路版图中各层之间的过孔处、各版图之间的互连线处以及用户设定的端口处；
用于根据混合阶单元建立有限元刚度矩阵，计算所述三维超大规模集成电路电磁场。
其中，所述第一模块包括：采用边交换法恢复所有所述多边形的边，在不同层的所述多边形的边的交点插入新的网格节点；
所述收集并设置所有层的集成电路版图多边形的层信息，将所有层的所述多边形投影到同一层并进行初始网格剖分；以及，所述采用边交换法恢复所有所述多边形的边，在不同层的所述多边形的边的交点插入新的网格节点(可以为先试图恢复多边形的边，如果两个相交的边都为多边形的边（即无法恢复时），则在交点插入新的网格节点)，包括：
第一单元，用于获取多层所述集成电路版图包含多个顶点的多个所述多边形，除此之外在顶层之上和底层之下再增加两个能覆盖所有层所述集成电路版图多边形的矩形多边形，分别作为顶层空气层和底层空气层的所述多边形；
第二单元，用于设置各层所述多边形的层信息，这种多边形的层信息应该既能唯一代表多边形所在的层，又能进行任意层之间的合并等布尔操作；如对于图2所示的多层集成电路版图，其为有4个金属层和3个介质层的集成电路模型，则共有2个空气层，4个金属层及3个介质层，即共有9层，为10个层界面；并用二进制数字代表各自的层，如果多边形位于哪个层界面，则该层界面所在的二进制的位为1。由于需要考虑金属层的厚度，因此位于金属层上的多边形的层信息同时包含金属层上下层界面。对图2所示的多层集成电路版图，可设置多边形的层信息如表1所示。
表1 图2所示多层集成电路各层多边形的层信息
第三单元，用于将各层的多个所述多边形垂直投影到同一层，根据Delaunay三角剖分算法形成以多边形顶点为网格节点的Delaunay三角形网格，其中，所述多边形的各个边包含预先设定的所在多边形的所述多边形信息和所在层的层信息；
第四单元，用于合并投影后重合的多个所述多边形的边的所述多边形信息和所述多边形的层信息。
例如图2所示的多层集成电路版图，某个金属层1和金属层2上的多边形重合，进行合并后，这两个多边形合并为一个多边形，其层信息为：6+24=30。
第五单元，用于根据所述边交换法将所述Delaunay三角形网格对齐到多个所述多边形的各个边，同时计算多个不同层的所述多边形的边的交点并将所述交点新增为所述多边形的顶点和所述Delaunay三角形网格的节点，形成第一三角形网格；
其中，在所述第一模块中，所述简化与对齐所述集成电路版图多边形，包括：
对投影到同一层的所述多边形的边进行简化与对齐，并对简化与对齐后的所述多边形及整个集成电路版图区域进行三角形网格剖分；
所述对投影到同一层的所述多边形的边进行简化与对齐，并对简化与对齐后的所述多边形及整个集成电路版图区域进行三角形网格剖分；以及，所述从xy平面方向对所述集成电路版图多边形进行三角形网格剖分和自适应细分；包括：
第六单元，用于基于所述第一三角形网格，在每个所述多边形P的内外分别形成夹住所述多边形P的内外辅助多边形P 0与P 9，并通过设定的距离阈值控制所述内外辅助多边形与所述多边形P的距离；
第七单元，用于对落在所述内外辅助多边形之间的各层所述多边形的边进行对齐和简化处理，并根据各个所述多边形的边所包含的多边形编号信息将投影到同一层的多层所述多边形还原到各层中，同时更新所述三角形网格及其层信息，形成第二三角形网格；
第八单元，用于直接基于网格尺寸和单元质量控制技术对所述包含对齐简化后的多边形的第二三角形网格进行自适应网格细分，形成更新的所述第二三角形网格。
上述第一模块，已经将所有层（包括顶层空气层和底层空气层）的集成电路版图多边形投影到同一层，并且经过对齐和简化处理，再在此基础上形成经过自适应细分的第二三角形网格，第二模块将在以上第一模块的基础上，直接快速形成三维多层超大规模集成电路的三棱柱网格剖分。
其中，所述第二模块：
用于根据包括金属层和介质层的每层纵向上的排列顺序（从上到下或从下到上），记录所述集成电路版图每层唯一的区域编号；例如，对于图2所示的多层集成电路版图，其共有2个空气层，4个金属层及3个介质层，可设置其区域编号如表2所示；
表2 图2所示多层集成电路各层区域编号信息
用于设定纵向上的坐标原点，根据所述集成电路版图每层的厚度及区域编号顺序，记录每层的位置信息，包括顶层空气层和底层空气层、金属层以及介质层，例如，对于图2所示的集成电路版图，其为有4个金属层和3个介质层的集成电路模型，则共有2个空气层，4个金属层及3个介质层，即共有9层，为10个层界面，记录这些层界面的位置，记为zinterface={z1, z2, …, z10}，其中z1~z10如图3所示，其排列顺序与上述的区域编号顺序一致。
用于基于所述第一模块获得的所述第二三角形网格，将所有所述三角形网格在厚度方向扩展为一个所述三棱柱，所述三棱柱的上下底面三角形为所述第二三角形网格对应的所述三角形，上下底面位置为层界面最大值和最小值，从而形成第一三棱柱网格；由于第一单元在所有层集成电路版图多边形的基础上增加了覆盖所有层集成电路版图多边形的矩形多边形作为顶层空气层和底层空气层的多边形，所以第二三角形网格的所有三角形的层均包含顶层层界面（即最上层层界面）和底层层界面（即最下层层界面），因此，所有三角形扩展的三棱柱上下底面位置相同，均为顶层层界面（即最上层层界面）和底层层界面（即最下层层界面）所在的位置。
用于所述第三单元获得的所述第一三棱柱网格，将所有所述三棱柱沿z方向（即在厚度方向）按zinterface进行分割，每个所述三棱柱被分割为N-1个所述三棱柱，形成第二三棱柱网格，其中N为所述集成电路版图包括上下空气层（顶层空气层和底层空气层）的层界面个数，zinterface为层界面的位置。
用于基于所述第七单元获得的所述第二三角形网格，根据所述三角形网格的层信息，还原所述三角形网格被哪些层界面的所述多边形共用；如对于图2所示的集成电路示意图，假设第二三角形网格中的某个三角形的层编码为543，则543=29+23+24+21+22+20，即该三角形网格被第0层（即底层空气层，用编码式中的20表示）、金属层1（用编码式中的21+22表示）、金属层2（用编码式中的23+24表示）及第9层（即顶层空气层，用编码式中的29表示）共用。
用于根据确定的所述三角形网格被哪些层界面的所述多边形共用，以此推算出其依次包含的区域，进而确定所述第二三棱柱网格中每个所述三棱柱所在的区域。如对于图2所示的集成电路示意图，假设第二三角形网格中的某个三角形的层编码为543，其被第0层（即底层空气层）、金属层1、金属层2及第9层（即顶层空气层）共用，则依据图3所示的多层集成电路版图的区域编号与层界面位置关系示意图，可确定该三角形网格对应的第二三棱柱网格中每个三棱柱所在的区域从下到上依次为：区域1、区域2、区域3、区域4、区域5、区域7、区域7、区域9和区域9。
可通过第二模块实现所述三棱柱网格在z方向统一的由细到粗的渐进变化的细分过程。
通常情况下，空气层厚度远大于介质层厚度，介质层厚度远大于金属层厚度，显然在z方向不进行网格细分是不合理的，为此本专利提出一种z方向渐变的细分方法。其思路在于，设置一个小尺寸到大尺寸的渐变过渡的近似比例r 0，使得z方向的尺寸从薄层到厚层是渐进变化的。针对本专利讨论的三维超大规模集成电路模型，薄层与厚层的位置关系可归结为以下两种情况：
1) 一边为薄层，一边为厚层：介质层/金属层-空气层
设薄层上三棱柱网格z方向上的长度为t 1, 厚层厚度为t 2，则应找到这样一个渐变比例r，使得
式中n为厚层上细分的段数，n应满足以下公式
显然，如果
则r=r 0，否则，可通过二分法计算r:
首先，设置n=1，逐渐增大n，使得n满足式（2），
然后，基于计算的n，设r left=r 0, 
反复进行以下计算：
Step 1: 设
r=(r left+ r right )/2
Step 2: 计算
如果
设
r
right
=r
转入step 1，否则，如果
设
r
left
=r
转入step 1，否则，r即为计算的正确结果，结束。式中
为设置的值为足够小的误差精度，如取
。
一种实施例为：设最上层金属层上等分3份，即将金属层在z方向分为3份，每个三棱柱网格z方向长度为5 nm，空气层厚度为1 cm，设置薄层到厚层的渐变过渡比例为r 0=1.5。经计算，可得厚层（空气层）上细分的段数n=33，r=1.50143756，三棱柱网格z方向长度由小到大依次如表3所示。
表3 金属层过渡到空气层时空气层的细分网格的长度渐进变化
| 层编号 | 长度（m） |
| 0 | 5E-9 |
| 1 | 7.5071880E-9 |
| 2 | 1.1271574E-8 |
| 3 | 1.6923564E-8 |
| 4 | 2.5409676E-8 |
| 5 | 3.8151041E-8 |
| 6 | 5.7281407E-8 |
| 7 | 8.6004456E-8 |
| 8 | 1.2913032E-7 |
| 9 | 1.9388112E-7 |
| 10 | 2.9110039E-7 |
| 11 | 4.3706907E-7 |
| 12 | 6.5623192E-7 |
| 13 | 9.8529125E-7 |
| 14 | 1.4793533E-6 |
| 15 | 2.2211566E-6 |
| 16 | 3.3349280E-6 |
| 17 | 5.0071862E-6 |
| 18 | 7.5179775E-6 |
| 19 | 1.1287774E-5 |
| 20 | 1.6947888E-5 |
| 21 | 2.5446195E-5 |
| 22 | 3.8205874E-5 |
| 23 | 5.7363734E-5 |
| 24 | 8.6128066E-5 |
| 25 | 1.2931591E-4 |
| 26 | 1.9415977E-4 |
| 27 | 2.9151877E-4 |
| 28 | 4.3769724E-4 |
| 29 | 6.5717508E-4 |
| 30 | 9.8670735E-4 |
| 31 | 1.4814795E-3 |
| 32 | 2.2243490E-3 |
| 33 | 3.3397211E-3 |
2) 上下为薄层，中间为厚层：金属层-介质层-金属层
设上下薄层（金属层）上三棱柱网格z方向上的长度为t 1,中间厚层（介质层）厚度为t 2，则应找到这样一个渐变比例r，使得
式中2n-1为厚层上细分的段数，n应满足以下公式
显然，如果
则r=r 0，否则，可通过二分法计算r:
首先，设置n=1，逐渐增大n，使得n满足式（5），
然后，基于计算的n，设r left=r 0, 
反复进行以下计算：
Step 1: 设
r=(r left+ r right )/2
Step 2: 计算
如果
设
r
right
=r
转入step 1，否则，如果
设
r
left
=r
转入step 1，否则，r即为计算的正确结果，结束。式中
为设置的值为足够小的误差精度，如取
。
另一种实施例为：设上下层金属层上等分3份，即将金属层在z方向分为3份，每个三棱柱网格z方向长度为5 nm，中间介质层厚度为1 μm，设置薄层到厚层的渐变过渡比例为r 0=1.5。经计算，可得厚层（空气层）上细分的段数2n-1=17，r=1.51561265，三棱柱网格z方向长度由小到大依次如表4所示。
表4 金属层-介质层-金属层过渡时介质层的细分网格的长度渐进变化
以上三棱柱网格在z方向统一的由细到粗的渐进变化的细分，一方面保证了三棱柱网格的渐进变化满足场的渐进变化规律，使得电磁仿真的计算结果更加准确，另一方面，这种自适应细分保证了通过二维非结构网格自适应细分的速度来实现三维网格自适应细分的速度，且不会导致局部过于密集的三棱柱网格。
如图4所示，将三棱柱网格剖分为四面体网格有以下6种剖分方法：
图4(a)，三棱柱网格123456被分为四面体网格1234, 2564和 2634；
图4(b)，三棱柱网格123456被分为四面体网格1234, 2534和5634；
图4(c)，三棱柱网格123456被分为四面体网格1235, 1645和1365；
图4(d)，三棱柱网格123456被分为四面体网格1235, 1345和3645；
图4(e)，三棱柱网格123456被分为四面体网格1236, 1526和1456；
图4(f)，三棱柱网格123456被分为四面体网格1236, 1426和2456。
将三棱柱网格剖分为四面体网格时，如果直接采用图4中的一种剖分方法进行剖分，可能因为相邻三棱柱网格的剖分不一致，导致公共面的剖分出现交叉的问题，如图5所示，这将使得公共面的定义出现二义性。
图5中，三棱柱网格123456被分为四面体网格1234, 2564和2634，三棱柱网格327658被分为四面体网格3278, 2358和3658，则四面体网格2564, 2634, 2358和3658。在公共面2563出现了交叉，剖分错误。
为避免这种错误发生，在所述第三模块中，所述实现所述三棱柱网格到四面体网格的网格剖分，采用一致剖分规则：
将所有所述三棱柱网格的顶点统一编号；
规则1：对每个所述三棱柱网格，由上层所述编号较小的点开始，向上层所述编号较大的点对应的下层点连接，将所述三棱柱网格剖分为3个所述四面体网格；或者，规则2：对每个所述三棱柱网格，由上层所述编号较大的点开始，向上层所述编号较小的点对应的下层点连接，将所述三棱柱网格剖分为3个所述四面体网格。
对采用快速三棱柱网格生成技术生成的三维超大规模集成电路的高质量网格，一致的采样规则1或规则2将所有三棱柱网格剖分成四面体网格，即可实现三棱柱网格到四面体网格的网格剖分。
如图6所示的相邻三棱柱网格123456和327658，采用规则1的剖分结果如图7(a)所示，剖分结果为：4132, 4523, 4536, 5237, 5673和5687。采用规则2的剖分结果如图7(b)所示，剖分结果为：6321, 6251, 6541, 8723, 8632和8256。
将三棱柱网格完全剖分为四面体网格后，即可针对四面体网格采用混合阶有限元的技术，实现基于后验误差估计选择性的提高某些单元的阶数，使得三维多尺度结构超大规模集成电路的电磁响应计算精度大幅提高。
本发明通过第一模块收集并设置所有层的集成电路版图多边形的层信息，将所有层的多边形投影到同一层并进行初始网格剖分；简化与对齐所述集成电路版图多边形，从xy平面方向对所述集成电路版图多边形进行所述三角形网格剖分和自适应细分；通过第二模块将三角形网格在厚度方向扩展为三棱柱网格，从z方向对所述三棱柱网格进行自适应细分；第三模块采用一致剖分规则实现所述三棱柱网格到四面体网格的网格剖分，基于所述四面体网格求解三维超大规模集成电路电磁场的初始值，根据所述电磁场的初始值变化速度超过和不超过预设值的所述四面体网格以及所述集成电路版图小尺度区域的所述四面体网格，建立混合阶单元，实现混合阶有限元技术，使得集成电路版图的电磁响应计算精度得到提高。实现了对三维多尺度结构超大规模集成电路进行快速高质量的三棱柱网格剖分，又实现了基于后验误差估计选择性的提高某些单元的阶数，使得三维多尺度结构超大规模集成电路的电磁响应计算精度大幅提高。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。
Claims (6)
1.一种用于集成电路三棱柱网格剖分的混合阶有限元方法，其特征在于，包括以下步骤：
步骤1，收集并设置所有层的集成电路版图多边形的层信息，将所有层的所述多边形投影到同一层并进行初始网格剖分；简化与对齐所述集成电路版图多边形，从xy平面方向对所述集成电路版图多边形进行三角形网格剖分和自适应细分；
步骤2，将三角形网格在厚度方向扩展为三棱柱网格，从z方向对三棱柱网格进行自适应细分；
步骤3，采用一致剖分规则实现所述三棱柱网格到四面体网格的网格剖分，基于所述四面体网格求解三维超大规模集成电路电磁场的初始值，根据所述电磁场的初始值变化速度超过和不超过预设值的所述四面体网格以及所述集成电路版图小尺度区域的所述四面体网格，建立混合阶单元，实现混合阶有限元技术，使得所述集成电路版图的电磁响应计算精度得到提高；
所述步骤1包括：采用边交换法恢复所有所述多边形的边，在不同层的所述多边形的边的交点插入新的网格节点；
所述收集并设置所有层的集成电路版图多边形的层信息，将所有层的所述多边形投影到同一层并进行初始网格剖分；以及，所述采用边交换法恢复所有所述多边形的边，在不同层的所述多边形的边的交点插入新的网格节点；包括：
步骤1.1，获取多层所述集成电路版图包含多个顶点的多个所述多边形，增加两个能覆盖所有层的所述集成电路版图多边形的矩形多边形，分别作为顶层空气层和底层空气层的所述多边形；
步骤1.2，设置各层所述多边形的层信息，并用二进制数字代表各自的层；
步骤1.3，将各层的多个所述多边形垂直投影到同一层，根据Delaunay三角剖分算法形成以多边形顶点为网格节点的Delaunay三角形网格，其中，所述多边形的各个边包含预先设定的所在多边形的所述多边形信息和所在层的层信息；
步骤1.4，合并投影后重合的多个所述多边形的边的所述多边形信息和所述多边形的层信息；
步骤1.5，根据所述边交换法将所述Delaunay三角形网格对齐到多个所述多边形的各个边，同时计算多个不同层的所述多边形的边的交点并将所述交点新增为所述多边形的顶点和所述Delaunay三角形网格的节点，形成第一三角形网格；
在所述步骤1中，所述简化与对齐所述集成电路版图多边形，包括：
对投影到同一层的所述多边形的边进行简化与对齐，并对简化与对齐后的所述多边形及整个集成电路版图区域进行三角形网格剖分；
所述对投影到同一层的所述多边形的边进行简化与对齐，并对简化与对齐后的所述多边形及整个集成电路版图区域进行三角形网格剖分；以及，所述从xy平面方向对所述集成电路版图多边形进行三角形网格剖分和自适应细分；包括：
步骤1.6，基于所述第一三角形网格，在每个所述多边形的内外分别形成夹住所述多边形的内外辅助多边形，并通过设定的距离阈值控制所述内外辅助多边形与所述多边形的距离；
步骤1.7，对落在所述内外辅助多边形之间的各层所述多边形的边进行对齐和简化处理，并根据各个所述多边形的边所包含的多边形编号信息将投影到同一层的多层所述多边形还原到各层中，同时更新所述三角形网格及其层信息，形成第二三角形网格；
步骤1.8，基于网格尺寸和网格单元质量控制对所述第二三角形网格进行自适应网格细分，形成更新的所述第二三角形网格。
2.根据权利要求1所述的混合阶有限元方法，其特征在于，所述步骤2包括：
步骤2.1：根据包括金属层和介质层的每层纵向上的排列顺序，记录所述集成电路版图每层唯一的区域编号；
步骤2.2：设定纵向上的坐标原点，根据所述集成电路版图每层的厚度及区域编号顺序，记录每层的位置信息；
步骤2.3：基于所述步骤1获得的所述第二三角形网格，将所有所述三角形网格在厚度方向扩展为一个所述三棱柱，所述三棱柱的上下底面三角形为所述第二三角形网格对应的所述三角形，上下底面位置为层界面最大值和最小值，从而形成第一三棱柱网格；
步骤2.4：基于所述步骤2.3获得的所述第一三棱柱网格，将所有所述三棱柱沿z方向按zinterface进行分割，每个所述三棱柱被分割为N-1个所述三棱柱，形成第二三棱柱网格，其中N为所述集成电路版图包括上下空气层的层界面个数，zinterface为层界面的位置；
步骤2.5：基于所述步骤1.7获得的所述第二三角形网格，根据所述三角形网格的层信息，还原所述三角形网格被哪些层界面的所述多边形共用；
步骤2.6：根据所述步骤2.5确定的所述三角形网格被哪些层界面的所述多边形共用，以此推算出其依次包含的区域，进而确定所述第二三棱柱网格中每个所述三棱柱所在的区域；
步骤2.7，实现所述三棱柱网格在z方向统一的由细到粗的渐进变化的细分过程。
3.根据权利要求1所述的混合阶有限元方法，其特征在于，在所述步骤3中，所述采用一致剖分规则实现所述三棱柱网格到四面体网格的网格剖分，包括：
将所有所述三棱柱网格的顶点统一编号；
对每个所述三棱柱网格，由上层所述编号较小的点开始，向上层所述编号较大的点对应的下层点连接，将所述三棱柱网格剖分为所述四面体网格；或者，由上层所述编号较大的点开始，向上层所述编号较小的点对应的下层点连接，将所述三棱柱网格剖分为所述四面体网格。
4.一种用于集成电路三棱柱网格剖分的混合阶有限元装置，其特征在于，包括：
第一模块，用于收集并设置所有层的集成电路版图多边形的层信息，将所有层的所述多边形投影到同一层并进行初始网格剖分；简化与对齐所述集成电路版图多边形，从xy平面方向对所述集成电路版图多边形进行三角形网格剖分和自适应细分；
第二模块，用于将三角形网格在厚度方向扩展为三棱柱网格，从z方向对三棱柱网格进行自适应细分；
第三模块，用于采用一致剖分规则实现所述三棱柱网格到四面体网格的网格剖分，基于所述四面体网格求解三维超大规模集成电路电磁场的初始值，根据所述电磁场的初始值变化速度超过和不超过预设值的所述四面体网格以及所述集成电路版图小尺度区域的所述四面体网格，建立混合阶单元，实现混合阶有限元技术，使得所述集成电路版图的电磁响应计算精度得到提高；
所述第一模块包括：采用边交换法恢复所有所述多边形的边，在不同层的所述多边形的边的交点插入新的网格节点；
所述收集并设置所有层的集成电路版图多边形的层信息，将所有层的所述多边形投影到同一层并进行初始网格剖分；以及，所述采用边交换法恢复所有所述多边形的边，在不同层的所述多边形的边的交点插入新的网格节点；包括：
第一单元，用于获取多层所述集成电路版图包含多个顶点的多个所述多边形，增加两个能覆盖所有层所述集成电路版图多边形的矩形多边形，分别作为顶层空气层和底层空气层的所述多边形；
第二单元，用于设置各层所述多边形的层信息，并用二进制数字代表各自的层；
第三单元，用于将各层的多个所述多边形垂直投影到同一层，根据Delaunay三角剖分算法形成以多边形顶点为网格节点的Delaunay三角形网格，其中，所述多边形的各个边包含预先设定的所在多边形的所述多边形信息和所在层的层信息；
第四单元，用于合并投影后重合的多个所述多边形的边的所述多边形信息和所述多边形的层信息；
第五单元，用于根据所述边交换法将所述Delaunay三角形网格对齐到多个所述多边形的各个边，同时计算多个不同层的所述多边形的边的交点并将所述交点新增为所述多边形的顶点和所述Delaunay三角形网格的节点，形成第一三角形网格；
在所述第一模块中，所述简化与对齐所述集成电路版图多边形，包括：
对投影到同一层的所述多边形的边进行简化与对齐，并对简化与对齐后的所述多边形及整个集成电路版图区域进行三角形网格剖分；
所述对投影到同一层的所述多边形的边进行简化与对齐，并对简化与对齐后的所述多边形及整个集成电路版图区域进行三角形网格剖分；以及，所述从xy平面方向对所述集成电路版图多边形进行三角形网格剖分和自适应细分；包括：
第六单元，用于基于所述第一三角形网格，在每个所述多边形的内外分别形成夹住所述多边形的内外辅助多边形，并通过设定的距离阈值控制所述内外辅助多边形与所述多边形的距离；
第七单元，用于对落在所述内外辅助多边形之间的各层所述多边形的边进行对齐和简化处理，并根据各个所述多边形的边所包含的多边形编号信息将投影到同一层的多层所述多边形还原到各层中，同时更新所述三角形网格及其层信息，形成第二三角形网格；
第八单元，用于基于网格尺寸和网格单元质量控制对所述第二三角形网格进行自适应网格细分，形成更新的所述第二三角形网格。
5.根据权利要求4所述的混合阶有限元装置，其特征在于，所述第二模块，
用于根据包括金属层和介质层的每层纵向上的排列顺序，记录所述集成电路版图每层唯一的区域编号；
用于设定纵向上的坐标原点，根据所述集成电路版图每层的厚度及区域编号顺序，记录每层的位置信息；
用于基于所述第一模块获得的所述第二三角形网格，将所有所述三角形网格在厚度方向扩展为一个所述三棱柱，所述三棱柱的上下底面三角形为所述第二三角形网格对应的所述三角形，上下底面位置为层界面最大值和最小值，从而形成第一三棱柱网格；
用于所述第三单元获得的所述第一三棱柱网格，将所有所述三棱柱沿z方向按zinterface进行分割，每个所述三棱柱被分割为N-1个所述三棱柱，形成第二三棱柱网格，其中N为所述集成电路版图包括上下空气层的层界面个数，zinterface为层界面的位置；
用于基于所述第七单元获得的所述第二三角形网格，根据所述三角形网格的层信息，还原所述三角形网格被哪些层界面的所述多边形共用；
用于根据确定的所述三角形网格被哪些层界面的所述多边形共用，以此推算出其依次包含的区域，进而确定所述第二三棱柱网格中每个所述三棱柱所在的区域；
用于实现所述三棱柱网格在z方向统一的由细到粗的渐进变化的细分过程。
6.根据权利要求4所述的混合阶有限元装置，其特征在于，在所述第三模块中，所述采用一致剖分规则实现所述三棱柱网格到四面体网格的网格剖分，包括：
将所有所述三棱柱网格的顶点统一编号；
对每个所述三棱柱网格，由上层所述编号较小的点开始，向上层所述编号较大的点对应的下层点连接，将所述三棱柱网格剖分为所述四面体网格；或者，由上层所述编号较大的点开始，向上层所述编号较小的点对应的下层点连接，将所述三棱柱网格剖分为所述四面体网格。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011334620.0A CN112131774B (zh) | 2020-11-24 | 2020-11-24 | 用于集成电路三棱柱网格剖分的混合阶有限元方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011334620.0A CN112131774B (zh) | 2020-11-24 | 2020-11-24 | 用于集成电路三棱柱网格剖分的混合阶有限元方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112131774A CN112131774A (zh) | 2020-12-25 |
| CN112131774B true CN112131774B (zh) | 2021-02-26 |
Family
ID=73852387
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011334620.0A Active CN112131774B (zh) | 2020-11-24 | 2020-11-24 | 用于集成电路三棱柱网格剖分的混合阶有限元方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112131774B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113094960B (zh) * | 2021-04-20 | 2021-09-24 | 北京航空航天大学 | 一种基于矩量法的多芯航空连接器电磁模型的快速生成方法 |
| CN114357941B (zh) * | 2022-02-22 | 2022-09-23 | 北京智芯仿真科技有限公司 | 一种集成电路版图电流密度超标区域优化方法及系统 |
| CN114970443B (zh) * | 2022-06-14 | 2024-08-23 | 上海格宇软件有限公司 | 一种基于区域分解和网格压印的芯片网格生成方法及系统 |
| CN117390639B (zh) * | 2023-10-12 | 2024-03-29 | 巨霖科技(上海)有限公司 | 一种非结构化网格加密方法和系统 |
| CN118334276B (zh) * | 2024-06-12 | 2024-08-09 | 芯瑞微(上海)电子科技有限公司 | 基于pcb模型的混合元共形网格生成方法及生成装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102930071A (zh) * | 2012-08-29 | 2013-02-13 | 电子科技大学 | 基于非匹配网格的周期结构的三维电磁场仿真模拟方法 |
| CN104318618A (zh) * | 2014-10-20 | 2015-01-28 | 南京师范大学 | 一种广义三棱柱空间数据模型的三维剖切方法 |
-
2020
- 2020-11-24 CN CN202011334620.0A patent/CN112131774B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102930071A (zh) * | 2012-08-29 | 2013-02-13 | 电子科技大学 | 基于非匹配网格的周期结构的三维电磁场仿真模拟方法 |
| CN104318618A (zh) * | 2014-10-20 | 2015-01-28 | 南京师范大学 | 一种广义三棱柱空间数据模型的三维剖切方法 |
Non-Patent Citations (5)
| Title |
|---|
| 《A Straightforward Updating Criterion for 2-D/3-D Hybrid Discontinuous Galerkin Time-Domain Method Controlling Comparative Error》;Wending Mai等;《IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES》;20180430;第66卷(第4期);正文第1713-1722页 * |
| 《An Efficient and Stable 2-D/3-D Hybrid Discontinuous Galerkin Time-Domain Analysis With Adaptive Criterion for Arbitrarily Shaped Antipads in Dispersive Parallel-Plate Pair》;Wending Mai 等;《IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES》;20171031;第65卷(第10期);正文第3671-3681页 * |
| 《基于三棱柱网格的三维矢量有限元方法研究》;邓攀;《中国硕士学位论文全文数据库（电子期刊）》;20090131;正文第I138-1200页 * |
| 《改进的二维三维混合时域不连续伽辽金方法》;买文鼎等;《电波科学学报》;20180228;第33卷(第1期);正文第33-40页 * |
| 《时域不连续伽略金方法基函数研究与应用》;买文鼎;《中国博士学位论文全文数据库（电子期刊）》;20200131;正文第A002-245页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112131774A (zh) | 2020-12-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112131774B (zh) | 用于集成电路三棱柱网格剖分的混合阶有限元方法及装置 | |
| CN112132973B (zh) | 三维集成电路电磁仿真全三维网格快速生成方法及装置 | |
| CN112149336B (zh) | 三维集成电路电磁仿真高质量网格快速生成方法及装置 | |
| CN111931458B (zh) | 基于混合阶有限元的三维集成电路电磁场计算方法及装置 | |
| TWI390707B (zh) | 半導體晶片、巨集擺置的方法、多向密集堆積巨集擺置器與多尺寸混合擺置設計方法 | |
| CN110807289B (zh) | 基于后验误差估计的集成电路自适应有限元网格细分方法 | |
| US7844939B2 (en) | Mask pattern correction program and system | |
| CN115859899A (zh) | 一种多驱动能力的集成电路标准单元版图迁移的方法 | |
| KR20230109646A (ko) | 혼합 높이 셀 라이브러리들에 대한 커스텀-타일 배치 패브릭들에 대한 적응적 행 패턴들 | |
| CN109551768A (zh) | 一种基于stl的3d打印文件的数据处理方法 | |
| CN111767688B (zh) | 集成电路版图多边形非结构网格自适应细分方法和系统 | |
| CN114330215A (zh) | 一种快速处理集成电路版图走线的电磁仿真方法及装置 | |
| CN111243094A (zh) | 一种基于点灯法的三维模型精确体素化方法 | |
| US8751975B2 (en) | RC corner solutions for double patterning technology | |
| CN111814422B (zh) | 集成电路电磁仿真多重网格法的插值矩阵构造方法及装置 | |
| CN112290955B (zh) | 基于集成电路阻抗网络提取的网格节点编码方法及系统 | |
| Mattison | A high quality, low cost router for MOS/LSI | |
| CN114330214B (zh) | 一种包含走线的集成电路快速高精度计算的方法及装置 | |
| CN112052641B (zh) | 大规模集成电路版图非结构网格偏心中点生成方法和系统 | |
| US20120304134A1 (en) | Exposure data generation method | |
| CN114330213B (zh) | 一种基于版图多边形走线的自动识别转换方法及装置 | |
| JP3554479B2 (ja) | 自動配置配線方法及び自動配置配線装置 | |
| JP3288190B2 (ja) | Ｌｓｉのレイアウト設計方法およびその装置 | |
| KR100336163B1 (ko) | 수치 해석을 위한 3차원 구조물 생성 시스템 및 방법 | |
| Xiang et al. | An algorithm for simultaneous pin assignment and routing |
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 |