CN1403966A - 一种用于虚拟多介质电容提取中最优切割数的生成方法 - Google Patents
一种用于虚拟多介质电容提取中最优切割数的生成方法 Download PDFInfo
- Publication number
- CN1403966A CN1403966A CN 02130850 CN02130850A CN1403966A CN 1403966 A CN1403966 A CN 1403966A CN 02130850 CN02130850 CN 02130850 CN 02130850 A CN02130850 A CN 02130850A CN 1403966 A CN1403966 A CN 1403966A
- Authority
- CN
- China
- Prior art keywords
- medium
- unit
- boundary element
- virtual
- length
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 127
- 238000005520 cutting process Methods 0.000 title claims abstract description 54
- 239000003990 capacitor Substances 0.000 title claims abstract description 35
- 238000000605 extraction Methods 0.000 title claims abstract description 26
- 230000008569 process Effects 0.000 title description 12
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 239000004020 conductor Substances 0.000 claims description 85
- JEIPFZHSYJVQDO-UHFFFAOYSA-N ferric oxide Chemical compound O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 claims description 33
- 239000000758 substrate Substances 0.000 claims description 31
- 239000002184 metal Substances 0.000 claims description 26
- 239000000284 extract Substances 0.000 claims description 17
- 238000010586 diagram Methods 0.000 claims description 13
- 238000002360 preparation method Methods 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 10
- 230000003287 optical effect Effects 0.000 abstract 1
- 239000002609 medium Substances 0.000 description 92
- 230000005684 electric field Effects 0.000 description 11
- 230000014509 gene expression Effects 0.000 description 6
- 230000008878 coupling Effects 0.000 description 5
- 238000010168 coupling process Methods 0.000 description 5
- 238000005859 coupling reaction Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 239000013598 vector Substances 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000011960 computer-aided design Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 230000003071 parasitic effect Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000001965 increasing effect Effects 0.000 description 2
- 230000001939 inductive effect Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000004323 axial length Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 230000005686 electrostatic field Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
Images
Landscapes
- Fixed Capacitors And Capacitor Manufacturing Machines (AREA)
Abstract
一种用于虚拟多介质电容提取中最优切割数的生成方法,属于集成电路计算机辅助设计中的三维互连电容提取技术领域。采用该方法可自动生成最优切割数,使电容提取速度达到最优或者近似最优。该方法包括由计算机依次执行的以下步骤:根据模拟结构的几何参数,确定候选切割数(m,n),这些(m,n)的值构成集合S;依次取集合S中的每个候选切割数(m,n),计算采用此切割数进行电容提取时Z参数的值(Z参数为系数矩阵A中非零元数目),得到集合Z(S);根据候选切割数集合S及与其对应的Z值集合Z(S),选定最优切割数。在虚拟多介质电容提取中,采用本发明,可在保证计算精度的情况下,提高自动化程度,缩短电容提取计算时间。
Description
技术领域
一种用于虚拟多介质电容提取中最优切割数的生成方法属于集成电路计算机辅助设计(IC-CAD)中的三维互连电容提取技术领域。
背景技术
集成电路(Integrated Circuit,IC)是当前电子工业、乃至信息产业的基石。随着半导体集成电路制造技术的发展,电路中金属互连线越来越窄(达到0.13微米以下),线与线的间距也越来越小。这使得互连线之间的电磁场寄生效应已成为影响诸如延迟、功耗和可靠性等电路性能的主要因素。因此,在集成电路设计中必须考虑互连寄生效应。
当前集成电路的设计流程如图1所示。要设计一块芯片,首先要提出功能描述,然后经过逻辑设计、版图设计得到描述半导体工艺结构的版图。这时需要进行一个“版图验证”的步骤,来验证电路设计是否能达到当初设定的性能要求,如果满足要求就可以进行生产制造,否则还要回到逻辑设计进行修正,重复一个迭代的过程。在版图验证中,一个重要的环节称为“寄生参数提取”,它包括计算金属互连线之间的寄生电阻、电容、电感等参数。只有得到这些互连寄生参数后,才能进行电路模拟以判断电路是否达到性能要求。
互连线间的寄生电容(简称“互连电容”)对电路性能的影响比较广泛,而且互连电阻、数模混合电路的衬底耦合和微电子机械***(Micro-Electro-Mechanical Systems,MEMS)等领域的研究也与之十分类似。互连电容的提取(计算)最受学术界和工业界的关注。在当前工艺特征下,为进行精确的电容提取,需要对集成电路版图的三维模型进行数值模拟,其主要方法包括区域型解法(有限差分方法和有限元方法)、直接边界元方法(Boundary ElementMethod,BEM)、间接边界元方法等。与区域型解法相比,边界元法的优点是精度高,较少的离散变量和较强处理复杂边界形状的能力。
来自实际版图的互连结构通常定义在有限区域内,其中静电场的电势可由带混合边界条件的拉普拉斯(Laplace)方程描述。由于直接边界元方法所离散的边界积分方程包含电势和法向电场强度两个变量,它比间接边界元法更适合这种有限域、多介质寄生电容结构的模拟。虽然直接边界元方法有上述优点,但由于工业实际中计算需求越来越多、其中考虑的三维结构也越来越大,如何提高三维电容提取方法的计算速度就成了当务之急。
本人在文献“喻文健,王泽毅,侯劲松,‘一种快速的三维VLSI互连电容提取方法:虚拟多介质方法’,电子学报,第29卷,第11期,pp.1526-1529,2001年”和“Wenjian Yu andZeyi Wang,‘An efficient quasi-multiple medium algorithm for the capacitance extraction of actual3-D VLSI interconnects’,in Proc.IEEE ASP-DAC 2001,Yokohama,Japan,Jan.2001,pp.366-371.”中提出虚拟多介质(Quasi-Multiple Medium,QMM)方法并用于三维多介质电容提取的计算中,它已成为一种重要的快速三维互连电容提取方法。
QMM电容提取的输入(通常为遵循一定格式的文本文件)是描述三维多介质互连结构的数据,其中包括几何形体数据和电学参数(导体上电压,介质的介电常数等)。互连结构整体是一个长方体区域,如图2所示,其中介质3分层排列,一些金属导体嵌入在某些介质层中。在计算中,假设其中相互连接的一块或多块导体的电压为1伏特,称为主导体1,而其他导体电压均为0伏特,称为环境导体2。金属衬底4与第1层介质下底面相贴,偏压为0伏特。
通过计算可得到各导体上的电量,设主导体上电量为Qm,其他第j块导体上电量为Qj,则导体j与主导体的耦合电容为:Cjm=-Qj/Vmj
由于电压差Vmj为1伏特,所以导体j与主导体间耦合电容在数值上即为-Qj。另外,主导体上的电量Qm在数值上等于其与环境导体间的总电容。总电容和耦合电容就是电容提取计算的输出结果。
QMM方法在直接边界元电容提取的过程中,增加了一个“虚拟切割”的步骤,即将分层介质结构垂直于底面均匀地切割为m×n块,如图3所示。这里的m,n为两个正整数,我们称二元整数对(m,n)为虚拟切割数。
为研究方便,从需要进行电容提取计算的集成电路版图中切出的一个长方体形状的多层介质互连结构作为一个模拟结构,其中包含多层介质以及嵌入其中的若干导体。在虚拟多介质电容提取方法中,首先要对模拟结构的各层介质进行切割,经过这种“虚拟切割”得到的这些介质区域称为虚拟介质区域。然后要对各介质区域(即虚拟介质区域)的边界表面进行离散化(划分为边界元)。这里的边界元划分采用的是非均匀边界元划分方法,即区分各种边界表面的类型采用不同大小的边界元划分间隔。对边界进行离散化后,就可以在每个介质区域中列出一系列离散的边界积分方程,其中包含了各边界元上的电势和法向电场强度变量,带入边界条件(导体边界上电势已知,模拟结构外边界上法向电场强度为零),再通过介质区域交界面上电势和电位移连续性方程,可将上述所有离散边界积分方程耦合成一个整体线性方程组。求解该方程组得到导体表面各边界元的法向电场强度,进而可计算出各导体的带电量,也即得到所需的互连电容值。
如图4所示,虚拟多介质电容提取具体包括如下六个步骤:
1.处理输入的模拟结构,为后面虚拟切割和非均匀边界元划分进行数据准备。首先将要计算的模拟结构(长方体形状)置于一个三维直角坐标系中,并使长方体的条边都对应地平行于坐标系的x、y、z三个轴。输入的信息包括各个介质和导体块的位置坐标和几何长度、各介质的相对介电常数、以及“主导体”的编号。处理上述几何信息,可以获得介质对导体的空间包含关系(用于判断各介质区域的边界),邻近主导体的导体(简称“邻近导体”,表示距离主导体较近,并与主导体之间无其他导体阻隔的导体),以及用于非均匀边界元划分方法中的参数,它们包括:用于含“内孔”表面(也称有孔表面)的边界元划分间隔MGx,MGy,MGz及EGx,EGy,EGz,和各形体(介质、导体块和金属衬底)的参考边界元划分份数ESx,ESy,ESz(由于对各形体计算出的数值可能不同,所以对每个形体需单独存储这三个参数)。
“内孔”指导体块嵌入介质表面,给原来完整的介质表面造成的“孔”。如图2所示,第3层介质中包含一些导体块,这些导体块的下底面与第3层介质的下底面相贴,相当于嵌在第三层介质的下底面。图5显示了第3层介质的下底面,对该面进行边界元划分时,应去掉这些导体块形成的内孔5,将其余部分,即有效的介质表面6划分为边界元。
在文献(古江春,王泽毅,洪先龙,“多孔平面的快速边界元划分”,计算机辅助设计与图形学学报,Vol.12,No.3,pp.211-215,2000年)中公开了对有孔表面划分边界元的切片法。切片法首先将有孔表面划分为母元,其主要步骤是:
首先,经过该面每个顶点和每个内孔多边形的顶点沿x轴方向作扫描线,若多个顶点在同一直线上,可以只做一条扫描线,如图6所示。每对相邻扫描线间的区域在内孔之外的部分被分为若干个梯形,其中,每个梯形称为切片单元7;然后,考虑紧挨着多个切片单元,判断它们是否可以合并形成一个梯形,如果可以就进行合并;最后,如图7所示,表面上的切片单元经过判断、合并后,得到的不可能进一步合并的新切片单元区域称为“母元”。
母元又可分为主(Master)母元和一般母元两种。主母元是在有孔介质外表面上的或有孔介质交界面上的紧挨主导体的母元,有孔表面的其他母元的就是一般母元。图7为对图5的有孔表面划分母元后的结果,紧挨主导体形成的内孔8分布着主母元9~14。
MGx,MGy,MGz分别为用于划分主母元的沿x、y、z轴方向的边界元划分间隔,
EGx,EGy,EGz分别为用于划分一般母元的沿x、y、z轴方向的边界元划分间隔,
它们的计算公式如下:
其中,上述公式中xlen为模拟结构x方向长度,ylen为模拟结构y方向长度,单位均为微米(下同)。由于主导体可能由互相衔接的多个导体块组成,故设MstLengthx,MstLengthy,MstLengthz分别表示组成主导体的各导体块沿x、y、z轴方向长度的最大值;
MaxMedZ为介质层高度(沿z轴方向)的最大值;
Min{}是取最小值函数,Max{}是取最大值函数,ln()为取自然对数函数。
ESx,ESy,ESz分别表示形体沿x、y、z轴方向的参考边界元划分份数,所述的形体根据分析对象的不同可为导体块,金属衬底或介质三种形体中的一种,下面区分不同形体,介绍ESx,ESy,ESz的计算:
其中,Lengthx,Lengthy,Lengthz分别为该形体沿x、y、z轴方向的长度,所述的形体根据分析对象的不同可为导体块,金属衬底或介质三种形体中的一种(下同)。
组成主导体的各导体块:
邻近导体块:
设整个主导体沿x轴方向长度为MasterLengthx,沿y轴方向长度为MasterLengthy。
2.确定虚拟切割数(m,n)。需要人为指定虚拟切割数,或者根据某种经验公式得到。经验公式为:
3.虚拟切割处理。垂直于xoy坐标平面,将整个模拟结构均匀切割为m×n块虚拟介质区域,嵌在介质中的导体块也可能被切成若干小块,生成表示这些新产生的形体的数据,同时,得到各虚拟介质块包含导体块的信息以及描述各介质区域有效边界表面的信息。
4.边界元划分。将每个介质区域(虚拟介质块)的有效边界表面划分为四边形边界元,边界元的信息包括四边形各顶点坐标、边界元类型(介质外表面、导体表面或介质交界面)、所属介质区域的编号。该步骤采用一种非均匀边界元划分方法,区别不同的情况用不同的边界元密度划分表面,从而减少边界元的数量。对一个有效的边界表面划分边界元具体包括如下几步:判断该面是否含内孔,若是,需先将它划分为“母元”。若该表面为介质无孔外表面、或介质无孔交界面、或导体表面(包括金属衬底),利
用参考边界元划分份数ES(包括ESx,ESy,ESz)确定该表面沿两个正交方向(如果是平行于xoy坐标平面的表面,则为沿x轴、y轴方向,依次类推)上的边界元划分间隔G1和G2;若该表面为介质有孔外表面、或介质有孔交界面,利用边界元划分间隔MG、EG确定该表面每个母元的两个正交方向上的边界元划分间隔G1和G2。按照G1和G2将该表面或该表面上所有的母元生成边界元的数据。
5.形成方程。根据边界元的信息,将介质i边界上的某个边界元Гk中心点作为配置点,可列离散边界积分方程:
其中:Ni为介质i边界上边界元的总数,Гj(j=1,…,Ni)代表第j号边界元,uj(j=1,…,Ni)为第j号边界元上点的电势(假设每个边界元上电势、电场强度均为常量),qj(j=1,…,Ni)为第j号边界元上点的电场强度的外法向(由该介质区域内指向外)分量。另外,积分号中的两个函数表达式为:
式中,r(k)为边界元k中心点到边界元Γj上积分点dГ的距离,n为边界元Γj上单位外法向矢量,r(k)为从边界元k中心点到边界元Γj上积分点dΓ的矢量,“(,)”为两矢量求内积符号。
方程(2)是关于离散变量uj,qj(j=1,…,Ni)的线性方程。然后,取遍所有配置点,即介质i边界上边界元Гk的中心点(i=1,…,M,k=1,…,Ni),其中M为介质区域的总个数,可得到一个完整的线性方程组。该线性方程组的离散变量是各边界元上的电势和法向电场强度,带入已知的边界条件(导体边界上电势为给定偏压,介质外边界上法向电场为零),整理后可得如下形式:
Ax=f (4)
其中x为各边界元上未知的电势和法向电场强度组成的向量。
6.求解方程。用预条件GMRES算法求解线性方程组(4),所有边界元上电势和法向电场强度均可得到。最后使用下述公式计算各导体所带电量:
其中Гci是导体i的表面,Nci为导体i表面的边界元数目,qj(j=1,…,Nci)为各边界元上的法向电场强度,ε为包含导体的介质的介电常数。
在上述虚拟多介质电容提取方法的第2个步骤中需要人为指定切割参数m和n(或者采用某种经验公式),这使得QMM电容提取的整体计算性能很大程度上依赖于用户的经验,并且电容提取的计算速度较慢。
发明内容
本发明的目的是在虚拟多介质电容提取方法中,采用一种新方法自动生成QMM切割数(m,n),使虚拟多介质电容提取的计算速度达到最优或者近似最优,从而可以高效提取三维电容。
本发明公开了一种用于虚拟多介质电容提取中最优切割数的生成方法,这种方法又称为“最小Z值的最优解”(OSMZ:optimal selection with minimal Z-value)方法,Z参数是系数矩阵A中非零元的数目;其特征在于此方法包括由计算机依次执行的以下步骤:
1)将从集成电路版图中切出的一个长方体形状的多层介质互连结构作为模拟结构,根据模拟结构的几何参数,取m为闭区间[mmin,mmax]内的任意整数,n为闭区间[nmin,nmax]内的任意整数,这样得到的所有整数对(m,n)构成集合S的初始值;设定判断模拟结构xoy截面长宽比是否失调的界值RL,利用RL来进一步筛选集合S中的元素;
其中,m:从模拟结构的xoy截面看,其沿y方向虚拟介质块的切割份数;
n:从模拟结构的xoy截面看,其沿x方向虚拟介质块的切割份数;
2)选取筛选后的集合S中的每个整数对作为候选切割数(m,n),分别计算采用此切割数进行虚拟多介质电容提取时Z参数的值,其中Z参数为系数矩阵A中非零元数目,最终得到与候选切割数集合S相对应的非零元数目的集合Z(S),其步骤依次如下:
2.1)对取出的一个候选切割数(m,n),计算其虚拟切割出的介质区域i的总数M,M=m×n×L,i=1,2…M;
其中,L为模拟结构中的介质层的数目;
2.2)将介质区域上的边界表面划分为两类:第一类为导体表面、金属衬底、介质无孔外表面和介质有孔外表面,第二类为介质无孔交界面和介质有孔交界面;对所有的介质区域,采用非均匀边界元划分方法分别计算其所有边界表面上的边界元数目,用数组单元a[i]存储介质区域i上所有的第一类边界表面上的边界元数目,用数组单元b[i]存储介质区域i上所有的第二类边界表面上的边界元数目,i=1,2,…,M;
2.3)根据公式
,计算与集合S中取出的候选切割数(m,n)相对应的Z值;
2.4)反复执行步骤2.1)至2.3),直到对集合S中的每个候选切割数进行了相应的处理;
3)根据集合S中各切割数对应的Z值的最小值确定阈值Zg;若集合S中的某元素(m,n)对应的Z值小于阈值Zg,并且与集合S中的其它元素相比,该元素的m和n乘积最小,则该元素即为最优切割数。
其中,xlen为模拟结构x方向长度,ylen为模拟结构y方向长度,单位均为微米;
在该方法步骤1)中的利用RL来进一步筛选集合S中元素的方法如下:
a.如果xlen/ylen大于RL,则保留集合S中满足“n>m”的元素,删除其他元素;否则进行步骤b;
其中,判断模拟结构xoy截面长宽比是否失调的界值RL大于1;
b.如果ylen/xlen大于RL,则保留集合S中满足“m>n”的元素,删除其他元素;否则保留集合S中满足“|n-m|<D”的元素,删除其他元素,其中,D为对于长宽比不失调的模拟结构,m和n两数取值之差的上限。
该方法中所述的RL的推荐值为1.4。
该方法中所述的m和n两数取值之差的上限D的推荐值为2。
该方法步骤3)中的选取最优切割数的方法如下:
首先,求集合Z(S)中的最小值,记为Z0;
其次,将Z0乘以大于1的系数Rz得到一个阈值Zg,即Zg=Z0×Rz;
再次,初始化Q为无限大整数;
最后,对集合S中的给定元素(m,n),如果其对应值Z(m,n)<Zg且m×n<Q,则mO=m,nO=n,Q=m×n;对集合S中所有的元素(m,n)都进行上述操作,最终得到集合S中的最优切割数(mO,nO)。
该方法中所述的系数Rz的推荐值为1.25。
该方法中所述的采用非均匀边界元划分方法计算导体表面上边界元数目ea1的方法如下:
首先,将已知的未进行虚拟切割前的导体块的参考边界元划分份数按长度比例分配到虚拟切割得到的相应子导体上,得到子导体块沿x、y、z轴方向的实际边界元划分份数Ex,Ey,Ez: Ez=ESz;
其中,“[]”为中取整符号,例如,[a]为最接近a的整数,若小于1,则取值1;
ESx,ESy,ESz,:当前表面所属的未进行虚拟切割前的形体的沿x,y,z轴方向的参考边界元划分份数,所述的形体根据分析对象的不同可为导体块,金属衬底或介质三种形体中的一种,它们的具体数值已由前序的模拟结构参数准备步骤计算出来;
由于虚拟多介质切割可能使原来的一块导体被切成多块子导体,故设:
Lengthx′,Lengthy′分别表示该面所属的进行虚拟切割后的子导体块沿x、y方向的长度;
Lengthx,Lengthy分别表示该面所属的导体块沿x、y方向的长度;
然后,区分该面的位置情况,分别计算边界元数目ea1:
若该面平行于xoy平面,ea1=Ex×Ey;
若该面平行于yoz平面,ea1=Ey×Ez;
若该面平行于zox平面,ea1=Ez×Ex。
该方法中所述的采用非均匀边界元划分方法计算金属衬底表面上边界元数目的方法如下:
金属衬底表面是一种平行于xoy平面的特殊导体表面,虚拟切割前的金属衬底在虚拟切割后变成m×n块子衬底,每个子衬底沿x,y轴方向的边界元划分份数由原金属衬底的参考边界元划分份数和切割数决定,其上边界元数目为
该方法中所述的采用非均匀边界元划分方法计算介质无孔外表面上边界元数目的方法如下:
首先,将已知的介质沿x,y轴方向的参考边界元划分份数均匀分配到虚拟介质块上,介质沿z轴方向的划分份数保持不变,得到虚拟介质块沿x、y、z轴方向的实际边界元划分份数Ex,Ey,Ez: Ez=ESz;
然后,区分虚拟介质块上该无孔介质表面的位置情况,分别计算边界元数目ea3:
若该面平行于xoy平面,ea3=Ex×Ey;
若该面平行于yoz平面:ea3=Ey×Ez;
若该面平行于zox平面:ea3=Ez×Ex。
该方法中采用非均匀边界元划分方法计算介质有孔外表面上边界元数目时,要使用切片法对此有孔表面划分母元,分别计算各母元的边界元数目,它们的和即为此有孔外表面上的边界元数目ea4,其计算母元边界元数目的方法如下:
首先,设当前母元沿x、y、z方向的划分间隔分别为Gx,Gy,Gz;
当前母元为主母元,则Gx=MGx,Gy=MGy,Gz=MGz;
否则当前母元为一般母元,Gx=EGx,Gy=EGy,Gz=EGz;
MGx,MGy,MGz:分别为用于划分主母元的沿x、y、z轴方向的边界元划分间隔,已由前序的模拟结构参数准备步骤计算出来;
EGx,EGy,EGz:分别为用于划分一般母元的沿x、y、z轴方向的边界元划分间隔,已由前序的模拟结构参数准备步骤计算出来;
然后,设Lengthmx,Lengthmy,Lengthmz为该母元沿x、y、z方向的长度,用它们分别除以当前母元沿x、y、z轴方向的划分间隔则得到该母元上实际的边界元划分份数Ex,Ey,Ez:
最后,区分该母元的位置情况,计算母元的边界元数目em:
若该母元平行于xoy平面,em=Ex×Ey;
若该母元平行于yoz平面,em=Ey×Ez;
若该母元平行于zox平面,em=Ez×Ex。
该方法中所述的采用非均匀边界元划分方法计算介质无孔交界面上边界元数目的方法如下:
首先,将已知的介质沿x,y轴方向的参考边界元划分份数均匀分配到虚拟介质块上,介质沿z轴方向的划分份数保持不变,得到虚拟介质块沿x、y、z轴方向的实际边界元划分份数Ex,Ey,Ez: Ez=ESz;
然后区分该无孔交界面的位置情况,分别计算边界元数目eb1值:
若该面平行于xoy平面,eb1=Ex×Ey;
对于平行yoz平面或zox平面的虚拟介质交界面的边界元划分,还需根据其离主导体的远近进行区分,离主导体越远,其上边界元就划分得越少,具体的公式如下:
若该面平行于yoz平面:
若该面所属介质层为主导体所在层,eb1=Ey×Ez;
若该面所属介质层为主导体所在层的上、下相邻层,
若该面所属介质层不是主导体所在层及其相邻层,
若该面平行于zox平面:
若该面所属介质层为主导体所在层,eb1=Ez×Ex;
若该面所属介质层为主导体所在层的上、下相邻层,
若该面所属介质层不是主导体所在层及其相邻层,
该方法中采用非均匀边界元划分方法计算介质有孔交界面上边界元数目时,要使用切片法对此有孔表面划分母元,分别计算各母元的边界元数目,它们的和即为此有孔交界面上的边界元数目eb2,其计算母元边界元数目的方法如下:
首先,设当前母元沿x、y、z方向的划分间隔分别为Gx,Gy,Gz;
当前母元为主母元,则Gx=MGx,Gy=MGy,Gz=MGz;
否则当前母元为一般母元,Gx=EGx,Gy=EGy,Gz=EGz;
MGx,MGy,MGz:分别为用于划分主母元的沿x、y、z轴方向的边界元划分间隔,已由前序的模拟结构参数准备步骤计算出来;
EGx,EGy,EGz:分别为用于划分一般母元的沿x、y、z轴方向的边界元划分间隔,已由前序的模拟结构参数准备步骤计算出来;
然后,设Lengthmx,Lengthmy,Lengthmz为该母元沿x、y、z方向的长度,用它们分别除以当前母元沿x、y、z轴方向的划分间隔则得到该母元上实际的边界元划分份数Ex,Ey,Ez:
最后,区分该母元的位置情况,计算母元的边界元数目em:
若该母元平行于xoy平面,em=Ex×Ey;
若该母元平行于yoz平面,em=Ey×Ey;
若该母元平行于zox平面,em=Ez×Ex。
虚拟多介质电容提取方法中,采用OSMZ方法可以在保持计算精度和较小内存用量的情况下,使虚拟多介质计算的时间近似达到最短,并较少地依赖于用户经验,使得自动化程度高。
附图说明
图1为集成电路设计流程。
图2为实施例的三维互连模拟结构示意图。
图3为采用虚拟多介质方法将图2中的每个介质层切割为3×2结构的示意图。
图4为QMM电容提取流程图。
图5为图2所示模拟结构的第3层介质下底面。
图6为利用切片法划分图5所示有孔表面的示意图。
图7为图5所示有孔表面经切片法划分后,合并生成的母元的示意图。
图8为OSMZ方法的流程框图。
图9为确定切割数取值范围S步骤的流程图。
图10为计算系数矩阵非零元数目Z步骤的流程图。
图11为求最优切割数步骤的流程图。
图12为本发明实施例的主视图。
图13为本发明实施例第三层介质的俯视图。
图14为本发明实施例第二层介质的俯视图。
图15为第一层介质各虚拟介质区域的编号示意图(俯视图)。
图16a为第一层介质中的虚拟介质区域1的主视图。
图16b为第一层介质中的虚拟介质区域1的俯视图。
具体实施方式
下面结合一个具体实例说明含OSMZ方法的虚拟多介质电容提取的执行过程。
在Sun Ultra Enterprise 450工作站上运行嵌入了含有OSMZ方法的QMM电容提取软件,来提取本实施例中主导体的总电容以及它与其他导体间的耦合电容。
QMM电容提取方法步骤如下:
OSMZ方法之前的信息处理和参数准备:
首先将互连模拟结构置于三维直角坐标系中,并使模拟结构区域的每条边都对应地平行于坐标系的x、y、z三个轴。
图12显示了本实施例的正视图,这是个含5层介质(D1,D2,D3,D4,D5)的模拟结构,模拟结构的大小是11×9×5.8(单位为微米,下同),从下到上各介质层高度为0.3,0.9,1,4,1.2,2.0。各层介质的相对介电常数分别为3.2,3.9,3.2,3.9,3.2。从下向上数第二、三层介质中包含导体块,导体块底面与介质层底面相贴,这两层导体块的高度分别为0.5,0.8,它们的位置和大小见图13和图14。
图12中黑方块表示主导体m7(1V偏压),灰色块表示环境导体(0V偏压)在整个结构的最底下是一层金属衬底平面(相当于接地导体,0V偏压),用粗线表示。
图13为从下往上数第三层介质的俯视图,其中显示了5块导体在xoy平面上的大小和位置。导体m1、m2、m3的尺寸均为1.5×0.7,导体m4的尺寸为0.5×1.5,导体m5的尺寸为11×1.5,其他几何尺寸已在图中标出。
图14为从下往上数第二层介质的俯视图,其中显示了4块导体在xoy平面上的大小和位置。四块导体m6,m7,m8和m9在x方向上宽度均为0.3,在y方向上长度分别2,4,1和1,其他几何尺寸已在图中标出。
计算MG、EG和ES等参数:
在本例中,主导体m7沿三个方向长度为:MstLengthx=0.3,MstLengthy=4,MstLengthz=0.5;x、y方向模拟结构的长度分别为xlen=11、ylen=9;介质层高度最大值MaxMedZ=2。根据公式(a),得到:
MGx=0.3,MGy=0.638,MGz=0.25;
EGx=1.1,EGy=1.275,EGz=0.5;
本例中邻近导体是m2、m5、m6、m8。根据各形体的尺寸和虚拟多介质电容提取步骤1中描述的方法,得到各形体(包括介质、导体和金属衬底)相应的ES参数,如表1所示。
表1 实施例中介质、导体和金属衬底相应的ES参数
D1 | D2 | D3 | D4 | D5 | m1 | m2 | m3 | m4 | m3 | m6 | m7 | m8 | m9 | 衬底 | |
ESx | 10 | 10 | 10 | 10 | 10 | 2 | 6 | 2 | 1 | 20 | 1 | 4 | 2 | 1 | 17 |
ESy | 9 | 9 | 9 | 9 | 9 | 1 | 1 | 1 | 2 | 2 | 2 | 8 | 1 | 1 | 14 |
ESz | 3 | 3 | 3 | 2 | 2 | 1 | 2 | 1 | 1 | 2 | 3 | 3 | 3 | 1 | / |
利用OSMZ方法确定最优虚拟切割数(mO,nO),图8是OSMZ方法的流程框图,分下面三个步骤执行:
步骤1)确定最优虚拟切割数的取值范围S,流程如图9所示;
1.2)RL取推荐值1.4;因为xlen/ylen<1.4,进行步骤1.3);
1.3)又由于ylen/xlen<1.4,所以模拟结构的长宽比不失调。D取推荐值2,故取S={(m,n)|m,n为整数,3≤m≤7,3≤n≤8,且|n-m|<2},即S={(3,3)(3,4)(4,3)(4,4)(4,5)(5,4)(5,5)(5,6)(6,5)(6,6)(6,7)(7,6)(7,7)(7,8)},共14个(m,n)取值。
步骤2)对集合S中的每个候选切割数(m,n),计算Z值,流程如图10所示;
以切割数取(3,3)时的情况为例,
2.1)介质区域数M=m×n×L=3×3×5=45。
2.2)将介质区域按1到45进行编号,最底层的介质被切割为9个介质区域,它们的编号如图15所示。然后从下往上顺序对各层的介质区域进行编号。初始化数组a[],b[]均有45个单元,每个单元内容清零。
2.3)、2.4)以第一层介质中的虚拟介质区域1为例,考虑它的每个表面,进行计算。
如图16a和图16b所示,这个介质区域共包括7个有效边界表面,分别为:
◆介质无孔外表面(x=0)
如表1所示,介质层D1的ES值为:ESx=10,ESy=9,ESz=3;对应的此虚拟介质块的参考划分数为
Ez=ESz=3.(对后面几个面也有用)。
此面平行于yoz平面,边界元数目ea3=Ey×Ez=9,为第一类边界元,a[1]=a[1]+ea3=9.
◆介质无孔外表面(y=0)
此面平行于zox平面,边界元数目ea3=Ez×Ex=9,为第一类边界元,a[1]=a[1]+ea3=18.
◆介质无孔交界面(x=3.67)
该面平行于yoz平面且该介质区域为主导体所在层的相邻层,则:
此面的边界元数目
为第二类边界元,b[1]=b[1]+eb1=4.
◆介质无孔交界面(y=3)
此面平行于zox平面,边界元数目
为第二类边界元,b[1]=b[1]+eb1=8.
◆金属衬底表面(z=0)
如表1所示,金属衬底的ES值为:ESx=17,ESy=14;
此面的边界元数目
为第一类边界元,a[1]=a[1]+ea2=48.
◆导体表面(m6的下底面,z=0.3)
如表1所示,导体块m6的ES值为:ESx=1,ESy=2,ESz=3;再计算:
此面平行于xoy平面,边界元数目ea1=Ex×Ey=2,为第一类边界元,a[1]=a[1]+ea1=50.
◆介质有孔交界面(z=0.3)
用切片法将该面划分为三个母元me1、me2、me3,如图16b所示,它们均为一般母元,划分间隔为:Gx=EGx=1.1,Gy=EGy=1.275,Gz=EGz=0.5;这些母元均平行于xoy平面,计算其边界元数目的过程如下:
母元me1:
母元me2:
母元me3:
该有孔交界面共有边界元eb2=12,为第二类边界元,b[1]=b[1]+eb2=20。
最后,得到介质区域1的两类边界元数目为a[1]=50,b[1]=20。
2.5)反复执行2.3)和2.4),得到所有45个介质区域的两类边界元数目,用数组a[],b[]存储,如表2所示。
表2数组a[],b[]的内容
i | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 |
a[i] | 50 | 63 | 51 | 39 | 38 | 39 | 48 | 39 | 48 | 51 | 120 | 56 | 9 | 32 | 9 |
b[i] | 20 | 58 | 21 | 21 | 74 | 21 | 17 | 21 | 17 | 37 | 108 | 39 | 45 | 132 | 45 |
i | 16 | 17 | 18 | 19 | 2a | 21 | 22 | 23 | 24 | 25 | 26 | 27 | 28 | 29 | 30 |
a[i] | 18 | 9 | 18 | 57 | 69 | 63 | 23 | 14 | 23 | 18 | 9 | 18 | 12 | 6 | 12 |
b[i] | 36 | 45 | 36 | 21 | 23 | 22 | 28 | 32 | 28 | 26 | 30 | 26 | 20 | 21 | 20 |
i | 31 | 32 | 33 | 34 | 35 | 36 | 37 | 38 | 39 | 40 | 41 | 42 | 43 | 44 | 45 |
a[i] | 6 | 0 | 6 | 12 | 6 | 12 | 21 | 15 | 21 | 15 | 9 | 15 | 21 | 15 | 21 |
b[i] | 21 | 22 | 21 | 20 | 21 | 20 | 11 | 12 | 11 | 12 | 13 | 12 | 11 | 12 | 11 |
2.6)计算相应的Z值,
结果为326555,即:Z(3,3)=326555。
2.7)反复执行2.1)-2.6)的步骤,得到与候选切割数集合S相对应的非零元数目的集合Z(S),结果如下:
Z(3,3)=326555; Z(3,4)=259006; Z(4,3)=286184;
Z(4,4)=230010; Z(4,5)=270737; Z(5,4)=184303;
Z(5,5)=215074; Z(5,6)=208758; Z(6,5)=216683;
Z(6,6)=215788; Z(6,7)=180822; Z(7,6)=149421;
Z(7,7)=134335; Z(8,7)=126451;
步骤3)根据集合S中各切割数对应的Z值大小,选定最优切割数,流程如图11所示;
3.1)求出最小值Z0=126451;
3.2)阈值Zg=126451×1.25=158063.75;
3.3)初始化Q为无穷大,从集合S中依次取值,判断是否满足条件:Z(m,n)<Zg且m×n<Q;第一个满足条件的是(7,6),更新Q=42;mO=7,nO=6;S中其他元素均不满足条件,最终得到最优切割数为(7,6)。
OSMZ方法之后的处理:
垂直于xoy坐标平面,将整个模拟结构均匀切割为7×6块虚拟介质区域,嵌在介质中的导体也可能被切成若干小块,生成表示这些新产生的形体的数据,同时,得到各虚拟介质块包含导体的信息以及描述各介质区域有效边界表面的信息。
将每个介质区域(虚拟介质块)的有效边界表面划分为四边形边界元,最后得到2355个边界元,得出的线性方程组“Ax=f”含3500个变量,矩阵A中共有149421个非零的元素。使用带预条件的GMRES方法求解线性方程组“Ax=f”,迭代17次后收敛。最终各导体上电量(单位为10-18库仑)为:
衬底: -415.720237
导体m6: -6.069719
导体m8: -0.699613
导体m9: -0.035718
导体m1 -1.335993
导体m2 -68.465130
导体m3 -12.096716
导体m4 -0.119530
导体m5 -202.377256
导体m7(主导体):718.467608
下面就计算本实施例的情况,通过与用现有的经验公式生成切割数的QMM电容提取方法比较,说明使用OSMZ方法进行QMM电容提取的优点。下表显示了两种QMM电容提取方法计算本实施例的有关数据:
表3两种提取方法的比较
经验公式生成切割数 | OSMZ方法 | |
切割数生成时间(秒) | 0 | 0.20 |
得到的虚拟切割数 | (3,3) | (7,6) |
系数矩阵非零元数目(Z参数) | 326555 | 149421 |
虚拟多介质电容提取总时间(秒) | 3.28 | 2.68 |
总电容值(10-18法拉) | 706.0 | 718.5 |
总内存用量(兆字节) | 2.38 | 4.56 |
从表3中看出,经验公式生成切割数的方法的运行时间几乎为零,而OSMZ方法由于要计算各种情况对应的Z值,需要一定的CPU时间,如本例为0.20秒。但两者得到的虚拟切割数是不同的,也导致整体的虚拟多介质电容提取的效率有很大差别。采用了OSMZ方法,整体电容提取时间为2.68秒,比采用经验公式方法缩短了0.6秒(速度提高约22%)。而从计算结果(总电容)看,两者的差别在2%以内,说明采用OSMZ方法保证了计算精度(通常差别5%为精度的合理范围)。而从内存的使用来看,使用OSMZ方法的电容提取用了较多内存,但4.56兆的大小仍然比很多其他电容提取软件小得多,在实际应用中影响不大。
与采用经验公式方法相比,虽然采用OSMZ方法本身要增加一些计算时间,但由它得到的最优切割数使整个虚拟多介质电容提取有较大计算速度的提高。另外根据OSMZ方法的原理,它比经验公式方法自动化程度更高,能对更广泛的各种结构都产生很好的计算效果。
Claims (13)
1.一种用于虚拟多介质电容提取中最优切割数的生成方法,其特征在于此方法包括由计算机依次执行的以下步骤:
1)将从集成电路版图中切出的一个长方体形状的多层介质互连结构作为模拟结构,根据模拟结构的几何参数,取m为闭区间[mmin,mmax]内的任意整数,n为闭区间[nmin,nmax]内的任意整数,这样得到的所有整数对(m,n)构成集合S的初始值;设定判断模拟结构xoy截面长宽比是否失调的界值RL,利用RL来进一步筛选集合S中的元素;
其中,m:从模拟结构的xoy截面看,其沿y方向虚拟介质块的切割份数;
n:从模拟结构的xoy截面看,其沿x方向虚拟介质块的切割份数;
2)选取筛选后的集合S中的每个整数对作为候选切割数(m,n),分别计算采用此切割数进行虚拟多介质电容提取时Z参数的值,其中Z参数为系数矩阵A中非零元数目,最终得到与候选切割数集合S相对应的非零元数目的集合Z(S),其步骤依次如下:
2.1)对取出的一个候选切割数(m,n),计算其虚拟切割出的介质区域i的总数M,M=m×n×L,i=1,2…M;
其中,上为模拟结构中的介质层的数目;
2.2)将介质区域上的边界表面划分为两类:第一类为导体表面、金属衬底、介质无孔外表面和介质有孔外表面,第二类为介质无孔交界面和介质有孔交界面;对所有的介质区域,采用非均匀边界元划分方法分别计算其所有边界表面上的边界元数目,用数组单元a[i]存储介质区域i上所有的第一类边界表面上的边界元数目,用数组单元b[i]存储介质区域i上所有的第二类边界表面上的边界元数目,i=1,2,…,M;
2.3)根据公式
,计算与集合S中取出的候选切割数(m,n)相对应的Z值;
2.4)反复执行步骤2.1)至2.3),直到对集合S中的每个候选切割数进行了相应的处理;
3)根据集合S中各切割数对应的Z值的最小值确定阈值Zg;若集合S中的某元素(m,n)对应的Z值小于阈值Zg,并且与集合S中的其它元素相比,该元素的m和n乘积最小,则该元素即为最优切割数。
3.根据权利要求1所述的用于虚拟多介质电容提取中最优切割数的生成方法,其特征在于所述步骤1)中的利用RL来进一步筛选集合S中元素的方法如下:
a.如果xlen/ylen大于RL,则保留集合S中满足“n>m”的元素,删除其他元素;否则进行步骤b;
其中,判断模拟结构xoy截面长宽比是否失调的界值RL大于1;
b.如果ylen/xlen大于RL,则保留集合S中满足“m>n”的元素,删除其他元素;否则保留集合S中满足“|n-m|<D”的元素,删除其他元素,其中,D为对于长宽比不失调的模拟结构,m和n两数取值之差的上限。
4.根据权利要求1或3所述的用于虚拟多介质电容提取中最优切割数的生成方法,其特征在于:所述的RL的推荐值为1.4。
5.根据权利要求3所述的用于虚拟多介质电容提取中最优切割数的生成方法,其特征在于:所述的m和n两数取值之差的上限D的推荐值为2。
6.根据权利要求1所述的用于虚拟多介质电容提取中最优切割数的生成方法,其特征在于所述步骤3)中的选取最优切割数的方法如下:
首先,求集合Z(S)中的最小值,记为Z0;
其次,将Z0乘以大于1的系数Rz得到一个阈值Zg,即Zg=Z0×Rz;
再次,初始化Q为无限大整数;
最后,对集合S中的给定元素(m,n),如果其对应值Z(m,n)<Zg且m×n<Q,则mO=m,nO=n,Q=m×n;对集合S中所有的元素(m,n)都进行上述操作,最终得到集合S中的最优切割数(mO,nO)。
7.根据权利要求6所述的用于虚拟多介质电容提取中最优切割数的生成方法,其特征在于:所述的系数Rz的推荐值为1.25。
8.根据权利要求1所述的用于虚拟多介质电容提取中最优切割数的生成方法,其特征在于所述的采用非均匀边界元划分方法计算导体表面上边界元数目ea1的方法如下:
首先,将已知的未进行虚拟切割前的导体块的参考边界元划分份数按长度比例分配到虚拟切割得到的相应子导体上,得到子导体块沿x、y、z轴方向的实际边界元划分份数Ex,Ey,Ez: Ez=ESz;
其中,“[]”为中取整符号,例如,[a]为最接近a的整数,若小于1,则取值1;
ESx,ESy,ESz,:当前表面所属的未进行虚拟切割前的形体的沿x,y,z轴方向的参考边界元划分份数,所述的形体根据分析对象的不同可为导体块,金属衬底或介质三种形体中的一种,它们的具体数值已由前序的模拟结构参数准备步骤计算出来;
由于虚拟多介质切割可能使原来的一块导体被切成多块子导体,故设:
Lengthx′,Lengthy′分别表示该面所属的进行虚拟切割后的子导体块沿x、y方向的长度;
Lengthx,Lengthy分别表示该面所属的导体块沿x、y方向的长度;
然后,区分该面的位置情况,分别计算边界元数目ea1:
若该面平行于xoy平面,ea1=Ex×Ey;
若该面平行于yoz平面,ea1=Ey×Ez;
若该面平行于zox平面,ea1=Ez×Ex。
9.根据权利要求1所述的用于虚拟多介质电容提取中最优切割数的生成方法,其特征在于所述的采用非均匀边界元划分方法计算金属衬底表面上边界元数目的方法如下:
金属衬底表面是一种平行于xoy平面的特殊导体表面,虚拟切割前的金属衬底在虚拟切割后变成m×n块子衬底,每个子衬底沿x,y轴方向的边界元划分份数由原金属衬底的参考边界元划分份数和切割数决定,其上边界元数目为
10.根据权利要求1所述的用于虚拟多介质电容提取中最优切割数的生成方法,其特征在于所述的采用非均匀边界元划分方法计算介质无孔外表面上边界元数目的方法如下:
首先,将已知的介质沿x,y轴方向的参考边界元划分份数均匀分配到虚拟介质块上,介质沿z轴方向的划分份数保持不变,得到虚拟介质块沿x、y、z轴方向的实际边界元划分份数Ex,Ey,Ez: Ez=ESz;
然后,区分虚拟介质块上该无孔介质表面的位置情况,分别计算边界元数目ea3:
若该面平行于xoy平面,ea3=Ex×Ey;
若该面平行于yoz平面:ea3=Ey×Ez;
若该面平行于zox平面:ea3=Ez×Ex。
11.根据权利要求1所述的用于虚拟多介质电容提取中最优切割数的生成方法,其特征在于采用非均匀边界元划分方法计算介质有孔外表面上边界元数目时,要使用切片法对此有孔表面划分母元,分别计算各母元的边界元数目,它们的和即为此有孔外表面上的边界元数目ea4,其计算母元边界元数目的方法如下:
首先,设当前母元沿x、y、z方向的划分间隔分别为Gx,Gy,Gz;
当前母元为主母元,则Gx=MGx,Gy=MGy,Gz=MGz;
否则当前母元为一般母元,Gx=EGx,Gy=EGy,Gz=EGz;
MGx,MGy,MGz:分别为用于划分主母元的沿x、y、z轴方向的边界元划分间隔,已由前序的模拟结构参数准备步骤计算出来;
EGx,EGy,EGz:分别为用于划分一般母元的沿x、y、z轴方向的边界元划分间隔,已由前序的模拟结构参数准备步骤计算出来;
最后,区分该母元的位置情况,计算母元的边界元数目em:
若该母元平行于xoy平面,em=Ex×Ey;
若该母元平行于yoz平面,em=Ey×Ez;
若该母元平行于zox平面,em=Ez×Ex。
12.根据权利要求1所述的用于虚拟多介质电容提取中最优切割数的生成方法,其特征在于所述的采用非均匀边界元划分方法计算介质无孔交界面上边界元数目的方法如下:
首先,将已知的介质沿x,y轴方向的参考边界元划分份数均匀分配到虚拟介质块上,介质沿z轴方向的划分份数保持不变,得到虚拟介质块沿x、y、z轴方向的实际边界元划分份数Ex,Ey,Ez: Ez=ESz;
然后区分该无孔交界面的位置情况,分别计算边界元数目eb1值:
若该面平行于xoy平面,eb1=Ex×Ey;
对于平行yoz平面或zox平面的虚拟介质交界面的边界元划分,还需根据其离主导体的远近进行区分,离主导体越远,其上边界元就划分得越少,具体的公式如下:
若该面平行于yoz平面:
若该面所属介质层为主导体所在层,eb1=Ey×Ez;
若该面所属介质层为主导体所在层的上、下相邻层,
若该面所属介质层不是主导体所在层及其相邻层,
若该面平行于zox平面:
若该面所属介质层为主导体所在层,eb1=Ez×Ex;
若该面所属介质层为主导体所在层的上、下相邻层,
若该面所属介质层不是主导体所在层及其相邻层,
13.根据权利要求1所述的用于虚拟多介质电容提取中最优切割数的生成方法,其特征在于采用非均匀边界元划分方法计算介质有孔交界面上边界元数目时,要使用切片法对此有孔表面划分母元,分别计算各母元的边界元数目,它们的和即为此有孔交界面上的边界元数目eb2,其计算母元边界元数目的方法如下:
首先,设当前母元沿x、y、z方向的划分间隔分别为Gx,Gy,Gz;
当前母元为主母元,则Gx=MGx,Gy=MGy,Gz=MGz;
否则当前母元为一般母元,Gx=EGx,Gy=EGy,Gz=EGz;
MGx,MGy,MGz:分别为用于划分主母元的沿x、y、z轴方向的边界元划分间隔,已由前序的模拟结构参数准备步骤计算出来;
EGx,EGy,EGz:分别为用于划分一般母元的沿x、y、z轴方向的边界元划分间隔,已由前序的模拟结构参数准备步骤计算出来;
最后,区分该母元的位置情况,计算母元的边界元数目em:
若该母元平行于xoy平面,em=Ex×Ey;
若该母元平行于yoz平面,em=Ey×Ez;
若该母元平行于zox平面,em=Ez×Ex。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB021308500A CN1153164C (zh) | 2002-10-11 | 2002-10-11 | 一种用于虚拟多介质电容提取中最优切割数的生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB021308500A CN1153164C (zh) | 2002-10-11 | 2002-10-11 | 一种用于虚拟多介质电容提取中最优切割数的生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1403966A true CN1403966A (zh) | 2003-03-19 |
CN1153164C CN1153164C (zh) | 2004-06-09 |
Family
ID=4746486
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB021308500A Expired - Fee Related CN1153164C (zh) | 2002-10-11 | 2002-10-11 | 一种用于虚拟多介质电容提取中最优切割数的生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1153164C (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100452063C (zh) * | 2006-06-07 | 2009-01-14 | 清华大学 | 硅集成电路衬底多频率点下综合耦合参数的快速提取方法 |
CN102024083A (zh) * | 2010-12-15 | 2011-04-20 | 中国科学院微电子研究所 | 一种提取含有冗余金属的互连结构的电容的方法 |
CN102063528A (zh) * | 2010-12-20 | 2011-05-18 | 西安电子科技大学 | 基于查找表算法的菱形冗余填充寄生电容提取方法 |
CN102651047A (zh) * | 2012-04-11 | 2012-08-29 | 清华大学 | 集成电路设计中基于随机行走的电容参数提取计算方法 |
CN103198177A (zh) * | 2013-03-11 | 2013-07-10 | 清华大学 | 基于gpu的集成电路电容参数提取***及方法 |
CN103473402A (zh) * | 2013-08-30 | 2013-12-25 | 清华大学 | 面向集成电路互连电容参数提取的空间管理数据生成方法 |
CN104723392A (zh) * | 2013-12-18 | 2015-06-24 | 株式会社御牧工程 | 边界确定方法和介质切割方法 |
CN110428506A (zh) * | 2019-08-09 | 2019-11-08 | 成都景中教育软件有限公司 | 一种基于参数的动态几何三维图形切割实现方法 |
-
2002
- 2002-10-11 CN CNB021308500A patent/CN1153164C/zh not_active Expired - Fee Related
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100452063C (zh) * | 2006-06-07 | 2009-01-14 | 清华大学 | 硅集成电路衬底多频率点下综合耦合参数的快速提取方法 |
CN102024083A (zh) * | 2010-12-15 | 2011-04-20 | 中国科学院微电子研究所 | 一种提取含有冗余金属的互连结构的电容的方法 |
CN102063528A (zh) * | 2010-12-20 | 2011-05-18 | 西安电子科技大学 | 基于查找表算法的菱形冗余填充寄生电容提取方法 |
CN102063528B (zh) * | 2010-12-20 | 2012-10-24 | 西安电子科技大学 | 基于查找表算法的菱形冗余填充寄生电容提取方法 |
CN102651047A (zh) * | 2012-04-11 | 2012-08-29 | 清华大学 | 集成电路设计中基于随机行走的电容参数提取计算方法 |
CN103198177A (zh) * | 2013-03-11 | 2013-07-10 | 清华大学 | 基于gpu的集成电路电容参数提取***及方法 |
CN103473402A (zh) * | 2013-08-30 | 2013-12-25 | 清华大学 | 面向集成电路互连电容参数提取的空间管理数据生成方法 |
CN103473402B (zh) * | 2013-08-30 | 2016-08-10 | 清华大学 | 面向集成电路互连电容参数提取的空间管理数据生成方法 |
CN104723392A (zh) * | 2013-12-18 | 2015-06-24 | 株式会社御牧工程 | 边界确定方法和介质切割方法 |
CN110428506A (zh) * | 2019-08-09 | 2019-11-08 | 成都景中教育软件有限公司 | 一种基于参数的动态几何三维图形切割实现方法 |
CN110428506B (zh) * | 2019-08-09 | 2023-04-25 | 成都景中教育软件有限公司 | 一种基于参数的动态几何三维图形切割实现方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1153164C (zh) | 2004-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1306563C (zh) | 模拟装置 | |
CN1881562A (zh) | 一种自动适应工艺特征尺寸的互连寄生电容提取方法 | |
CN1154950C (zh) | 制造具有优化响应特性的物体的方法和设备 | |
CN1206722C (zh) | 基于等效电路的集成电路电源网络瞬态分析求解的方法 | |
CN1212642C (zh) | 掩模图形制成方法及半导体装置的制造方法 | |
TWI803690B (zh) | 推估方法、推估裝置、模型之產生方法及學習裝置 | |
CN1860489A (zh) | 使用最小阻力路径来特征化储层模型中的连通性 | |
CN1779686A (zh) | 确定缓冲器***的方法和*** | |
CN1369114A (zh) | 修改集成电路的方法 | |
CN1804827A (zh) | 海底热液活动探测数据处理和信息管理方法 | |
CN1205546C (zh) | 山区公路平面线形自动仿真设计方法 | |
CN1403966A (zh) | 一种用于虚拟多介质电容提取中最优切割数的生成方法 | |
CN1862546A (zh) | 一个快速的集成电路可布性分析方法 | |
CN1279480C (zh) | 考虑耦合效应进行时延优化的标准单元总体布线方法 | |
CN1681085A (zh) | 电子束描绘***、方法、程序及直接描绘制造半导体器件方法 | |
CN1796942A (zh) | 城际导航路径的计算方法 | |
CN1702466A (zh) | 记录了蓄电元件的等效电路模型的记录介质、导出程序 | |
CN100347713C (zh) | 基于最小自由度优先原则的提高集装箱装箱效率的方法 | |
CN1300731C (zh) | 半导体集成电路及其设计方法 | |
CN1503341A (zh) | 半导体器件的布局检验方法 | |
CN1485769A (zh) | 掩膜数据加工装置 | |
CN1643526A (zh) | 边界数据的内外判定方法及其程序 | |
CN1547252A (zh) | 集成电路布图规划与缓冲器规划集成的布局方法 | |
CN1240015C (zh) | 标准单元总体布线时障碍下时延驱动直角斯坦纳树方法 | |
CN1867920A (zh) | 物体的奇异向量物理量的数值建模过程、相关***和程序产品 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |