CN102945303A - 一种基于径向基网络算法获取集成电路成品率的方法 - Google Patents
一种基于径向基网络算法获取集成电路成品率的方法 Download PDFInfo
- Publication number
- CN102945303A CN102945303A CN2012104519757A CN201210451975A CN102945303A CN 102945303 A CN102945303 A CN 102945303A CN 2012104519757 A CN2012104519757 A CN 2012104519757A CN 201210451975 A CN201210451975 A CN 201210451975A CN 102945303 A CN102945303 A CN 102945303A
- Authority
- CN
- China
- Prior art keywords
- integrated circuits
- integrated circuit
- sampled point
- yield
- 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.)
- Pending
Links
Images
Abstract
本发明涉及一种基于径向基网络算法获取集成电路成品率的方法,属于集成电路技术领域,该方法包括:根据集成电路工艺厂商提供的工艺参数,采用径向基网络算法,建立一个替代电路仿真的替代模型,将工艺参数作为替代模型的自变量,电路性能指标作为替代模型的函数值;根据最小范数方法,获取最易使集成电路失效的工艺浮动值;得到的替代模型和最易失效的工艺浮动值,进行统计采样,获取采样点和电路性能指标;根据所述采样点及其电路性能指标,通过统计学方法得到该集成电路的成品率。该方法可降低成品率获取过程中的电路仿真次数,减少分析集成电路成品率所用的时间,缩短集成电路设计周期,加快集成电路生产,降低集成电路的成本,提高经济价值。
Description
技术领域
本发明属于集成电路技术领域,特涉及获取集成电路成品率的方法。
背景技术
集成电路的生产分为两个阶段。第一阶段,集成电路的设计者进行集成电路设计:设计者通过集成电路仿真工具对所设计的集成电路进行仿真分析,获得所设计的集成电路性能。当所设计的性能满足要求后,设计者将设计好的集成电路设计方案交给集成电路生产厂商。第二阶段,集成电路生产厂商进行集成电路制造:生产厂商根据设计者提供的集成电路设计方案,采用集成电路工艺技术,制造出符合设计者要求的集成电路成品。
工艺浮动是指集成电路制造过程中,由于工艺不一致所产生的偏差,包括阈值电压,氧化层厚度的浮动等等。随着集成电路工艺尺寸逐渐接近物理极限,工艺浮动对于集成电路制造的影响愈发严重。
集成电路成品率是指实际制造的全部集成电路中,符合设计者要求的集成电路占全部集成电路的比例。由于工艺浮动的存在,集成电路成品率不能达到100%。成品率的高低决定了集成电路的成本。成品率越高,单个集成电路的成本就越低,带来的经济效益就越高。所以,成品率成为了集成电路领域一个重要的指标,需要集成电路设计者和集成电路生产厂商共同努力,提高集成电路的成品率。
集成电路生产厂商需要改进集成电路制造工艺,减小集成电路制造上所产生的工艺偏差,以提高成品率。
集成电路设计者需要在设计集成电路时,不仅仅考虑所设计的集成电路性能,还需要获取所设计的集成电路成品率大小,如果成品率不符合要求,则重新设计集成电路,直到集成电路满足成品率要求为止。
综上所述:集成电路设计者在设计集成电路时需要预先获取成品率的要求,当设计集成电路时获得的成品率较低时,设计者需要通过改进集成电路设计,提高成品率,获得更大的经济效益。
现阶段,集成电路设计者广泛应用的成品率获取方法是蒙特卡洛方法(该方法由X.Li,J.Le,和L.T.Pileggi提出,题目为统计性能建模和优化Statistical performance modeling andoptimization,发表在Foundations and Trends in Electronic Design Automation,vol.1,no.4,pp.331-480,Apr.2006)。该方法的具体流程如图1所示。该方法获取成品率的过程具体包括:
第一步,根据集成电路制造厂商提供的工艺浮动参数(具体表示为X~p(X),其中,X为工艺浮动向量,p(X)为工艺浮动向量X所服从的统计分布函数)进行统计采样,得到一个采样点X1;
第二步,通过集成电路仿真工具获取该采样点X1的电路性能指标f1;
重复上述两个步骤,当采样点数量达到收敛要求时(收敛要求为统计理论中优值系数Figure of Merit等于0.1时,所需要的采样点数量),停止采样;这时,获得的采样点和对应的性能指标可以表示为{X1,f1},{X2,f2},…,{XN,fN},N为收敛时的采样数量;
第三步,根据所述N个采样点和对应的性能指标,通过统计学方法得到该集成电路的成品率。
该方法的一个主要问题是获取集成电路成品率的时间太长,有时需要进行上百万次采样,才能获得集成电路的成品率。上百万次的采样是非常耗时的,以每次采样需要1秒为例,为了获取一个集成电路的成品率,则需要大约10天的时间。因此,发明一个更加快速的成品率获取方法是十分必要的。
发明内容
本发明的目的在于克服已有技术的不足之处,提出一种新的基于径向基网络算法获取集成电路成品率的方法,该方法可降低成品率获取过程中的电路仿真次数,从而减少分析集成电路成品率所用的时间,缩短集成电路设计周期,加快集成电路生产,降低集成电路的成本,提高经济价值。
本发明区别于传统的蒙特卡洛方法的特点及有益效果:
本发明应用径向基网络算法建立一个替代模型,把获取集成电路成品率所需的电路仿真转移到替代模型上。替代模型的计算时间会比电路仿真工具的仿真时间快三到四个数量级。
同时,本发明不是直接根据工艺浮动参数X~p(X)进行采样,而是首先获取最易使集成电路失效的工艺浮动值Xopt,根据X~p(X-Xopt)进行采样(p(X-Xopt)表示对原有工艺参数的统计分布p(X)进行偏移,偏移向量为最易失效的工艺浮动值),使得采样点的数量减少一到两个数量级。
附图说明
图1为蒙特卡洛方法的流程图。
图2为本发明的方法流程图。
图3为实验所用的电路图。
具体实施方式
下面结合附图及实施例进一步详细说明本发明的具体内容。
本发明的基本流程如图2所示,具体包括以下步骤:
1)根据集成电路工艺厂商提供的工艺参数X,采用径向基网络算法,建立一个替代电路仿真的替代模型(替代模型的作用是形成从工艺参数X到电路性能指标f的解析映射),将工艺参数X作为替代模型的自变量,电路性能指标f作为替代模型的函数值,即f=f(X),从而替代电路仿真工具(该算法由S.-F.Su,C.-C.Chuang,C.Tao,J.-T.Jeng和C.-C.Hsiao提出,题目为用于区间型符号数据的线性区间回归权重的径向基网络算法Radial basisfunction networks with linear interval regression weights for symbolic interval data,发表在IEEE Transactions on Systems,Man,and Cybernetics,Part B:Cybernetics,vol.42,no.1,pp.69–80,Feb.2012);
2)根据最小范数方法,获取最易使集成电路失效的工艺浮动值Xopt(该最小范数方法由L.Dolecek,M.Qazi,D.Shah和A.Chandrakasan提出,题目为克服仿真局限:通过最小范数方法评估静态存储器Breaking the simulation barrier:Sram evaluation through normminimization,发表在IEEE/ACM International Conference on Computer-Aided Design,Nov.2008,pp.322–329.);
3)根据1)得到的替代模型和2)得到的最易失效的工艺浮动值Xopt,进行统计采样,获取采样点和电路性能指标,所述采样过程具体包括以下步骤:
31)根据X~p(X-Xopt)进行采样(p(X-Xopt)表示对原有工艺参数的统计分布p(X)进行偏移,偏移向量为最易失效的工艺浮动值),得到一个采样点X1′;
32)通过所述的替代模型获取该采样点X1′的电路性能f1′;
重复上述31)、32)两步,当采样点数量达到收敛要求时,停止采样,获得的采样点和对应的性能指标可以表示为{X1′,f1′},{X2′,f2′},…,{XM′,fM′},M为收敛时的采样数量;
4)根据所述M个采样点及其电路性能指标,通过统计学方法得到该集成电路的成品率。
本发明的实施例如下:
1)本实施例所需要获取成品率的集成电路为静态存储器,如图3所示,该集成电路包括六个场效应管(M1-M6),WL为存储器字线,BL,BLB为存储器位线,VDD为电源电压,GND为地线;
2)获取集成电路商用工艺厂商提供的六个场效应管的工艺浮动参数X,该工艺浮动参数X为服从独立正态统计分布的向量,即X~N(6,0,1),其中N()表示正态分布,6表示工艺浮动向量的维度,0表示均值,1表示方差。
3)根据所述工艺浮动参数向量采用径向基网络算法,建立替代模型:建立的替代模型如下:
其中f为电路性能指标,X为工艺浮动参数,ω1,ω2,....,ω100为径向基网络权值ω1=0.014,ω2=0.028,ω3=0.012,...,ω100=0.003,M1,M2,..,M100为径向基网络基向量;
M1=[0,0,0,0,0,0]T,M2=[-0.33,-0.93,-0.27,-2.13,1.05,2.01]T,
M3=[-0.39,0.81,0.87,-0.93,-0.45,0.57]T,…,M100=[2.91,-2.61,1.11,2.91,-2.49,1.05]T;
(本实施例建立替代模型的时间为3.21分钟)
4)根据最小范数方法,获取最易使集成电路失效的工艺浮动值:获得的最易失效的工艺浮动值为Xopt=[-0.0202,1.4024,0.9744,-0.3569,-0.5086,-0.1235]T。得到最易失效的工艺浮动值的时间为9.91分钟;
5)根据得到的替代模型和得到的最易失效的工艺浮动值Xopt,进行统计采样,具体包括以下步骤:
51)根据X~p(X-[-0.0202,1.4024,0.9744,-0.3569,-0.5086,-0.1235]T)进行采样,得到一个采样点;
52)通过所述的替代模型获取该采样点的电路性能;
重复上述51)、52)两步,当采样点数量达到1619次时,满足收敛要求,停止采样,采样时间为7.13分钟;
6)根据这1619个采样点及其电路性能指标,获取成品率:通过统计学方法得到该集成电路的成品率,获取成品率时间为1.36秒。
本实施例最终获取的成品率为99.97%,本实施例所消耗的全部时间为20.25分钟。
为了比较本发明所带来的益处,通过蒙特卡洛方法获取同样电路的成品率,最终获取的成品率同样为99.97%,所需要的采样点数量为33441次,所消耗的时间为9.29小时。
与蒙特卡洛方法相比,本发明在可以准确获取成品率的同时,效率更高,可以更快的获取集成电路的成品率。本发明较之蒙特卡洛方法可以将获取时间缩短27.5倍。
Claims (1)
1.一种基于径向基网络算法获取集成电路成品率的方法,其特征在于,该方法包括以下步骤:
1)根据集成电路工艺厂商提供的工艺参数X,采用径向基网络算法,建立一个替代电路仿真的替代模型,将工艺参数X作为替代模型的自变量,电路性能指标f作为替代模型的函数值,即f=f(X),从而替代电路仿真工具;
2)根据最小范数方法,获取最易使集成电路失效的工艺浮动值Xopt;
3)根据1)得到的替代模型和2)得到的最易失效的工艺浮动值Xopt,进行统计采样,获取采样点和电路性能指标,所述采样过程具体包括以下步骤:
31)根据X~p(X-Xopt)进行采样,p(X-Xopt)表示对原有工艺参数的统计分布p(X)进行偏移,偏移向量为最易失效的工艺浮动值,得到一个采样点X1′;
32)通过所述的替代模型获取该采样点X1′的电路性能f1′;
重复上述31)、32)两步,当采样点数量达到收敛要求时,停止采样,获得的采样点和对应的性能指标可以表示为{X1′,f1′},{X2′,f2′},…,{XM′,fM′},M为收敛时的采样数量;
4)根据所述M个采样点及其电路性能指标,通过统计学方法得到该集成电路的成品率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012104519757A CN102945303A (zh) | 2012-11-12 | 2012-11-12 | 一种基于径向基网络算法获取集成电路成品率的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012104519757A CN102945303A (zh) | 2012-11-12 | 2012-11-12 | 一种基于径向基网络算法获取集成电路成品率的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102945303A true CN102945303A (zh) | 2013-02-27 |
Family
ID=47728247
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012104519757A Pending CN102945303A (zh) | 2012-11-12 | 2012-11-12 | 一种基于径向基网络算法获取集成电路成品率的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102945303A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108021738A (zh) * | 2017-11-15 | 2018-05-11 | 南京邮电大学 | 一种具有稀疏特性的芯片单性能成品率预测方法 |
CN108614904A (zh) * | 2016-12-11 | 2018-10-02 | 复旦大学 | 基于仿真数据相关性的集成电路参数成品率快速分析方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1563966A (zh) * | 2004-03-25 | 2005-01-12 | 杭州电子工业学院 | 一种基于核估计的集成电路成品率测定方法 |
US7493574B2 (en) * | 2006-02-23 | 2009-02-17 | Cadence Designs Systems, Inc. | Method and system for improving yield of an integrated circuit |
US20120253775A1 (en) * | 2011-04-01 | 2012-10-04 | Taiwan Semiconductor Manufacturing Company, Ltd. | Multidimensional Monte-Carlo Simulation for Yield Prediction |
-
2012
- 2012-11-12 CN CN2012104519757A patent/CN102945303A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1563966A (zh) * | 2004-03-25 | 2005-01-12 | 杭州电子工业学院 | 一种基于核估计的集成电路成品率测定方法 |
US7493574B2 (en) * | 2006-02-23 | 2009-02-17 | Cadence Designs Systems, Inc. | Method and system for improving yield of an integrated circuit |
US20120253775A1 (en) * | 2011-04-01 | 2012-10-04 | Taiwan Semiconductor Manufacturing Company, Ltd. | Multidimensional Monte-Carlo Simulation for Yield Prediction |
Non-Patent Citations (6)
Title |
---|
LARA DOLECEK等: "Breaking the Simulation Barrier:SRAM Evaluation Through Norm Minimization", 《INTERNATIONAL CONFERENCE ON COMPUTER-AIDED DESIGN》, 31 October 2008 (2008-10-31), pages 322 - 329 * |
YU BEN: "Statistical Verification and Optimization of Integrated Circuits", 《PROQUEST DISSERTATIONS & THESES DATABASE》, 31 December 2011 (2011-12-31) * |
伊燕平等: "基于径向基函数神经网络的地下水数值模拟模型的替代模型研究", 《水土保持研究》, vol. 19, no. 4, 15 August 2012 (2012-08-15), pages 265 - 269 * |
杨茂峰: "基于代理模型的电路成品率计算与优化研究", 《中国优秀硕士学位论文全文数据库信息科技辑》, 15 January 2010 (2010-01-15) * |
王向东等: "基于神经网络的集成电路生产过程建模与优化", 《自动化学报》, vol. 27, no. 3, 22 May 2001 (2001-05-22), pages 289 - 295 * |
荆明娥: "集成电路参数成品率的预测与优化技术研究", 《中国博士学位论文全文数据库信息科技辑》, 15 June 2005 (2005-06-15) * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108614904A (zh) * | 2016-12-11 | 2018-10-02 | 复旦大学 | 基于仿真数据相关性的集成电路参数成品率快速分析方法 |
CN108614904B (zh) * | 2016-12-11 | 2021-10-12 | 复旦大学 | 基于仿真数据相关性的集成电路参数成品率快速分析方法 |
CN108021738A (zh) * | 2017-11-15 | 2018-05-11 | 南京邮电大学 | 一种具有稀疏特性的芯片单性能成品率预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107533576A (zh) | 针对使用电路模板的电路设计的提取的布局依赖效应的重用 | |
CN103577643A (zh) | 一种sram型fpga单粒子翻转效应仿真方法 | |
CN112257361B (zh) | 一种基于二次拟合模型的标准单元库构建方法 | |
CN103838905A (zh) | 可描述高低温失配特性的晶体管模型及仿真方法 | |
CN101706831B (zh) | 半导体设计仿真领域中电路容差测量方法 | |
US8813006B1 (en) | Accelerated characterization of circuits for within-die process variations | |
CN103645646B (zh) | 内含微小电感或电容的复杂电路系统建模与控制方法 | |
CN101196936B (zh) | Mos晶体管电学统计模型的快速建模方法 | |
CN102945303A (zh) | 一种基于径向基网络算法获取集成电路成品率的方法 | |
CN106649920A (zh) | 一种基于ibis的集成电路总剂量效应建模方法 | |
CN107292026A (zh) | 一种工艺参数波动引起mosfet性能变化的估计方法 | |
CN102542098A (zh) | 一种fpga互连线延时获取方法及其系统 | |
TWI824233B (zh) | 具有漏電流中的統計變化之動態隨機存取記憶體通路電晶體的設計 | |
CN108052727B (zh) | 一种金属栅功函数变化导致栅电容统计分布的估计方法 | |
CN108491606B (zh) | 一种材料强度分布获取方法 | |
US10755015B2 (en) | Agnostic model of semiconductor devices and related methods | |
Li et al. | A yield-enhanced global optimization methodology for analog circuit based on extreme value theory | |
Han et al. | Diagnosis method for analog circuit hard fault and soft fault | |
US11669667B2 (en) | Automatic test pattern generation (ATPG) for parametric faults | |
US20220171910A1 (en) | Modelling timing behavior using augmented sensitivity data for physical parameters | |
Ye et al. | Fast variational interconnect delay and slew computation using quadratic models | |
CN117252130B (zh) | 飞行器近场压力分布预测方法、系统及相关设备 | |
US11568117B2 (en) | Generating simulation-friendly compact physical models for passive structures | |
CN103455689B (zh) | 基于泰勒级数展开的区间模型修正方法 | |
CN103235865B (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 | ||
C05 | Deemed withdrawal (patent law before 1993) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20130227 |