CN101261658B - 基于电压延迟行为的铁电电容行为仿真用等效电路 - Google Patents
基于电压延迟行为的铁电电容行为仿真用等效电路 Download PDFInfo
- Publication number
- CN101261658B CN101261658B CN2008101034476A CN200810103447A CN101261658B CN 101261658 B CN101261658 B CN 101261658B CN 2008101034476 A CN2008101034476 A CN 2008101034476A CN 200810103447 A CN200810103447 A CN 200810103447A CN 101261658 B CN101261658 B CN 101261658B
- Authority
- CN
- China
- Prior art keywords
- voltage
- circuit
- ferroelectric capacitor
- model
- behavior
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 239000003990 capacitor Substances 0.000 title claims description 55
- 230000008859 change Effects 0.000 claims abstract description 20
- 239000004615 ingredient Substances 0.000 claims description 4
- 229910052793 cadmium Inorganic materials 0.000 claims description 3
- 229910052802 copper Inorganic materials 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000000034 method Methods 0.000 abstract description 8
- 238000004088 simulation Methods 0.000 abstract description 8
- 238000013461 design Methods 0.000 abstract description 7
- 235000013599 spices Nutrition 0.000 abstract description 5
- 230000001052 transient effect Effects 0.000 description 7
- 230000010287 polarization Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 239000000463 material Substances 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000005684 electric field Effects 0.000 description 4
- 230000005294 ferromagnetic effect Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000005291 magnetic effect Effects 0.000 description 4
- 230000003068 static effect Effects 0.000 description 2
- 101000878457 Macrocallista nimbosa FMRFamide Proteins 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000002902 bimodal effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000013341 scale-up Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000007306 turnover Effects 0.000 description 1
Images
Landscapes
- Design And Manufacture Of Integrated Circuits (AREA)
Abstract
本发明涉及一种基于电压延迟行为的铁电电容行为模型,属于集成电路设计技术领域。所述模型具有两个端口节点,正节点和负节点;在正负节点间并联连接一个线性电容,一个非线性电容和一个压控电压源;压控电压源的正输出端连接一个电压变化率跟踪电路,负输出端接地;所述模型通过正负端口节点与外电路相连。该模型利用SPICE程序中的表达式功能建立非线性电容原件,并配合电压延迟行为元件作为控制电路,来模拟铁电电容的物理行为,具有结构简单仿真速度快,仿真精度高等优点,适合各种基于SPICE的电路仿真程序。
Description
技术领域
本发明属于集成电路设计技术领域,特别涉及一种基于电压延迟行为的铁电电容行为仿真用等效电路。
背景技术
铁电电容是一种利用铁电材料制作的电容元件,在当今的集成电路设计尤其是存储器电路设计中,具有越来越重要的地位。对于规范而精确的电路设计,基于EDA(电子设计辅助)软件的电路仿真是必不可少的。然而铁电电容作为一种新型的元件,还没有在各种EDA软件中得到独立的模拟模型,给集成电路设计带来了很大困难。因此目前适合EDA软件的铁电电容模型设计仍然是一个很有挑战性的工作。铁电电容模型大体上可分为两种:物理模型和行为模型。物理模型能直观反应铁电材料的物理特性,但构造物理模型需要较深的物理基础,且缺少直观的电路描述,很难在EDA软件中实现。行为模型则抛开铁电理论的细节,只是单纯地从电路的角度观察铁电材料的宏观电学特征,并进行描述,不需要深厚的铁电学物理知识背景,因而比较容易实现。以下简要介绍一下已有的几种行为模型。
(1)双电容模型:这是一个早期的建立在电滞回线基础上的模型。这个模型将电滞回线在纵轴右边的两支用两条直线来近似,每条直线对应一个线性电容,电容值即直线的斜率,模型参数可以通过从实验曲线中提取得到,模型结构参照图3。但是由于电滞回线不但跟铁电电容大小有关,还跟所加的信号波形有关,因此当改变信号后,参数需要重新提取。在实际存储器电路中,加在铁电电容上的信号幅度和波形是无法预知的,这样不同情形下的模型参数只能设定得很粗略,在实际模拟中存在很大误差。因此这个模型是一个很粗糙的电容模型。
(2)数值模型:用一个双曲正切函数来模拟回线。数值模型为铁电电容静态特性分析提供了一种很好的方法,但不适合作瞬态分析,因为它不包含任何和瞬态行为有关的参数。
(3)分布式反转场模型:在偶极子上加上足够的电场时,偶极子会从一种状态翻转为另一种状态,使偶极子发生极化方向翻转的临界电场即矫顽电场。对于单畴的铁电体,这个矫顽电场跟这种材料的电滞回线的矫顽场一致;对于多畴铁电体,由于电畴的方向各不相同,不同电畴在不同的外加电压下发生反转。当电畴数目非常大时,不同电压下发生反转的电畴数目对电压形成高斯分布,这样,假定每个电畴的极化对总极化强度的贡献相同,则每增加单位电压所导致的极化强度的增加对电压也符合高斯分布。这个模型和数值模型相比,本质上是类似的。
(4)基于铁磁的模型:铁电体和铁磁体有许多类似之处,因此可以利用已有的铁磁核模型来仿真铁电模型。基于铁磁基础的模型将加在铁电电容上的电压按比例转换为电流源,加在磁芯上,产生基于磁滞回线的电压,再将电压按比例转换为通过铁电电容的电流,模型结构参照图4,其中K1、K2为两个比例系数。如果磁芯模型精度提高,铁电电容类磁模型的精度也可以进一步提高,但这个模型不适于作瞬态分析,因为这两者在瞬态特性上没有共性。
(5)基于电流的激化反转模型:这个模型既可以作静态分析也可以作瞬态分析。模型将铁电电容看作电容和一个电流源并联。电容由铁电材料的小信号介电系数来确定,电流源代表电容的极化反转电流,通过对电容进行脉冲测试得到。对于每种不同的脉冲幅度,电流源参数需要重新确定,这增加了模型的复杂性,不利于模型的实际应用,因为加在铁电电容上的电压事先是不知道的。另外模型参数对电路的时间常数也具有依赖型。
(6)零反转时间瞬态模型(ZSTT):模型假定铁电电容的极化反转时间为零,当所研究电路的RC常数远远大于铁电电容的极化反转时间时,假定所带来的误差是可以忽略的。反转时间为零意味着铁电电容极化电荷的增加是瞬时的,电荷的增量只是外加电压合铁电电容初始状态的函数,它与时间无关。ZSTT铁电电容模型的电路主要由两部分组成,一是控制铁电电容状态转换的开关,二是C0(V0)和C1(V1)两个非线性电容,模型结构参见图5,其中S01、S02、S11和S12为电容状态转换开关。电容开关的逻辑控制需要复杂的电路,不利于模拟。同时,两个非线性电容的值通常采用分段线性的数据列表,对不同的电容,要更换大量数据才能仿真。
以上的各种模型或是不够精确,或是使用不便,或是结构复杂难于控制,因此应用在EDA工具中时都有所不便。SPICE是目前最常用的EDA电路仿真工具,针对这一工具的铁电电容模型变得非常重要。
发明内容
本发明的目的是提供一种基于电压延迟行为的铁电电容行为仿真用等效电路,便于应用在SPICE EDA工具中,可以实现快速、高精度的仿真,同时结构简单、控制方便。为实现所述目的,本发明所采用的技术方案为:
一种基于电压延迟行为的铁电电容行为仿真用等效电路,所述仿真用等效电路具有正节点和负节点;在正、负节点间并联连接一个线性电容,一个非线性电容和一个压控电压源;压控电压源的正输出端连接一个电压变化率跟踪电路,负输出端接地;所述仿真用等效电路通过正、负节点接电压V;
所述电压变化率跟踪电路,其组成部分包括一个电压延迟行为元件和一个差分方程电路;所述差分方程电路具有n+,n-和out三个端口节点;所述电压延迟行为元件的两端分别连接所述压控电压源的正输出端和所述差分方程电路的端口节点n+;所述差分方程电路的端口节点n-与所述压控电压源的正输出端相连;所述差分方程电路的端口节点out端的信号V(out)是电压V的变化率,是仿真用等效电路的控制信号;所述差分方程电路利用表达式(V(N+)-V(n-))/Δt计算电压V的变化率;所述电压延迟行为元件将电压V延迟时间Δt;
所述非线性电容,其电容值的计算表达式为:
其 中:
C1(V)表示描述以电压V为参数的铁电电容行为曲线的右半支、C2(V)表示描述以电压V为参数的铁电电容行为曲线的左半支;α、Cu、Cd表示描述铁电电容行为曲线左右两支形状特征的参数;V表示施加在铁电电容行为仿真用等效电路两端的电压,Vc表示铁电电容的矫顽电压,也是描述铁电电容特征的参数;
V(out)是电压V的变化率;sgn()是符号函数,用来判断变量的符号,如果是正号,则值是1;如果是负号,则值是-1;V(n+1)的值是电压V延迟了时间Δt后的电压值,V(n-)的值是电压V当前的电压值。
本发明中铁电电容的电学行为完全由非线性电容Cx表现,可以用表达式进行精确的赋值,精确度提高很多。此模型的工作控制完全由微分电路完成,控制过程简单了许多,有利于提高仿真速度。由于控制信号是输入电压对时间的微分,即电压变化率,因此更适合瞬态仿真。
附图说明
图1为一种基于电压延迟行为的铁电电容行为模型;
图2为电容模型中的线性电容的电容曲线与电滞回线。
图3为双电容模型结构;
图4为基于铁磁的模型结构;
图5为ZSTT模型结构。
具体实施方式
本发明提供了一种基于电压延迟行为的铁电电容行为模型,下面结合附图对本发明作进一步详细的说明。
图1所示为一种基于电压延迟行为的铁电电容行为模型,其组成部分包括:一个线性电容Cy,一个非线性电容Cy,一个压控电压源VCVS和一个电压变化率跟踪电路。所述电容模型包括1、2两个端点,对电容模型施加的信号电压V加在这两个端点上,在这两个端点上连接着描述电容行为的主要元件:线性电容Cy和非线性电容Cx;压控电压源VCVS的增益为1,其输入端V+,V-端口间电压为施加到1、2端点间的电压V,输出端O+,O-端口间的电压等于输入端端口的电压,即等于V。从输出端O+,O-端口出来的电压施加到电压变化率跟踪电路上,跟踪电路输出端out的电压值V(out)即是电压V的变化率。
线性电容Cy是一个普通电容元件,其大小是一个模型参数;非线性电容Cx的大小由表达式的值决定。图2为电容模型中线性电容的电容曲线与电滞回线,横轴表示施加电压V,Vm和-Vm为V的正向和负向最大值;纵轴对于电容曲线是电容值,对于电滞回线是电量值。电容曲线是V轴上方的双峰曲线,电滞回线是中心对称的方形回线。如图所示电滞回线和电容曲线都包含两个分支,电滞回线的ACB分支和电容曲线的右侧峰对应施加电压V从-Vm向Vm变化的过程;电滞回线的BDA分支和电容曲线的左侧峰对应施加电压V从Vm向-Vm变化的过程;两个过程如图中箭头方向所示。非线性电容Cx的表达式分为描述右侧峰的C1(V)和描述左侧峰的C2(V),他们的大小由施加电压V决定,具体表达式为:
C1(V)表示描述以电压V为参数的铁电电容行为曲线的右半支、C2(V)表示描述以电压V为参数的铁电电容行为曲线的左半支;α、Cu、Cd表示描述铁电电容行为曲线左右两支形状特征的参数,可以从电容测试数据中提取得到;V表示施加在电容模型两端的电压,Vc表示铁电电容的矫顽电压,也是描述铁电电容特征的参数。
由于右侧峰对应一个电压上升的过程,左侧峰对应一个电压下降的过程,因此C1(V)对应的施加电压变化率是正值,C2(V)对应的施加电压变化率是负值,通过判断V(out)的正负号就可以确定该选择哪一个表达式来计算Cx的值,因此:
式中sgn()是符号函数,用来判断变量的符号,如果是正号,则值是1;如果是负号,则值是-1。通过这种简单的计算,当V(out)是正号时Cx=C1(V),V(out)是负号时Cx=C2(V)。右侧峰C1(V)的其形状描述参数为Cu,下标u表示电压上升;左侧峰C2(V)的形状描述参数为Cd,下标d表示电压下降;由于对于实际电容来说左右两支的形状通常不是完全对称的,因此这两个参数一般不相等.
电压变化率跟踪电路的组成部分包括:一个电压延迟行为元件Δt以及一个差分方程电路。如图中所示:V(n+)的值是施加电压V延迟了时间Δt后的电压值,V(n-)的值是施加电压V当前的电压值,V(out)=(V(n+)-V(n-))/Δt,因此V(out)表示的是施加电压V的变化率,二者符号相同。电压延迟行为元件和差分方程电路都可以用SPICE软件中的行为描述电压源来描述。
对比本模型和普通电容,二者都是具有两个端点,因此本模型在使用时也可以像普通电容一样直接连接到所需连接的端口上。利用本模型可以完成各种利用铁电电容的电路模拟。
Claims (1)
1.一种基于电压延迟行为的铁电电容行为仿真用等效电路,其特征在于:所述仿真用等效电路具有正节点和负节点;在正、负节点间并联连接一个线性电容,一个非线性电容和一个压控电压源;压控电压源的正输出端连接一个电压变化率跟踪电路,负输出端接地;所述仿真用等效电路通过正、负节点接电压V;
所述电压变化率跟踪电路,其组成部分包括一个电压延迟行为元件和一个差分方程电路;所述差分方程电路具有n+,n-和out三个端口节点;所述电压延迟行为元件的两端分别连接所述压控电压源的正输出端和所述差分方程电路的端口节点n+;所述差分方程电路的端口节点n-与所述压控电压源的正输出端相连;所述差分方程电路的端口节点out端的信号V(out)是电压V的变化率,是仿真用等效电路的控制信号;所述差分方程电路利用表达式(V(n+)-V(n-))/Δt计算电压V的变化率;所述电压延迟行为元件将电压V延迟时间Δt;
所述非线性电容,其电容值的计算表达式为:
其中:
C1(V)表示描述以电压V为参数的铁电电容行为曲线的右半支、C2(V)表示描述以电压V为参数的铁电电容行为曲线的左半支;α、Cu、Cd表示描述铁电电容行为曲线左右两支形状特征的参数;V表示施加在铁电电容行为仿真用等效电路两端的电压,Vc表示铁电电容的矫顽电压,也是描述铁电电容特征的参数;
V(out)是电压V的变化率;sgn()是符号函数,用来判断变量的符号,如果是正号,则值是1;如果是负号,则值是-1;V(n+)的值是电压V延迟了时间Δt后的电压值,V(n-)的值是电压V当前的电压值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101034476A CN101261658B (zh) | 2008-04-03 | 2008-04-03 | 基于电压延迟行为的铁电电容行为仿真用等效电路 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101034476A CN101261658B (zh) | 2008-04-03 | 2008-04-03 | 基于电压延迟行为的铁电电容行为仿真用等效电路 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101261658A CN101261658A (zh) | 2008-09-10 |
CN101261658B true CN101261658B (zh) | 2010-10-13 |
Family
ID=39962112
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008101034476A Expired - Fee Related CN101261658B (zh) | 2008-04-03 | 2008-04-03 | 基于电压延迟行为的铁电电容行为仿真用等效电路 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101261658B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101694675B (zh) * | 2009-10-22 | 2014-12-31 | 中兴通讯股份有限公司 | 电路仿真方法和系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1303126A (zh) * | 2000-12-08 | 2001-07-11 | 清华大学 | 一种用于模拟铁电电容的等效电路 |
CN1888923A (zh) * | 2006-07-17 | 2007-01-03 | 大连理工大学 | 铁电薄膜电滞回线的测量补偿方法 |
-
2008
- 2008-04-03 CN CN2008101034476A patent/CN101261658B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1303126A (zh) * | 2000-12-08 | 2001-07-11 | 清华大学 | 一种用于模拟铁电电容的等效电路 |
CN1888923A (zh) * | 2006-07-17 | 2007-01-03 | 大连理工大学 | 铁电薄膜电滞回线的测量补偿方法 |
Non-Patent Citations (4)
Title |
---|
CN 1888923 A,全文. |
JP特開2003-218345A 2003.07.31 |
任天令,张武全,李春晓,陈宏毅,朱钧,刘理天.一种改进的用于铁电存储器设计的铁电电容宏模型.电子学报29 8.2001,29(8),1135-1137. |
任天令,张武全,李春晓,陈宏毅,朱钧,刘理天.一种改进的用于铁电存储器设计的铁电电容宏模型.电子学报29 8.2001,29(8),1135-1137. * |
Also Published As
Publication number | Publication date |
---|---|
CN101261658A (zh) | 2008-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102623062B (zh) | 一种模拟忆阻器的电路 | |
CN102236736B (zh) | 一种十字形cmos集成霍尔磁传感器的电路模型的仿真方法 | |
CN105553459A (zh) | 浮地压控忆阻器仿真器电路 | |
CN206991310U (zh) | 一种对数型忆容器等效模拟电路 | |
CN103336866B (zh) | 一种电磁暂态仿真中含负电阻支路的处理方法 | |
CN106896327A (zh) | 分数阶KiBaM‑等效电路综合特征电池模型及其参数辨识方法 | |
Bebikhov et al. | The application of mathematical simulation for solution of linear algebraic and ordinary differential equations in electrical engineering | |
Zhu et al. | Accurate and scalable IO buffer macromodel based on surrogate modeling | |
Fontana et al. | A new simulation program for analog circuits using symbolic analysis techniques | |
CN109543339A (zh) | 三电平换流器的定导纳建模与实时仿真方法 | |
CN101800512A (zh) | 动态幅度线性可调的混沌信号源 | |
CN110399677A (zh) | 基于偏磁状态下改进j-a公式的变压器直流偏磁仿真模拟方法 | |
CN101295454B (zh) | 一种无电感蔡氏电路 | |
CN105224755A (zh) | 一种通用的金属互联层电容spice预测模型的构建方法 | |
CN100570618C (zh) | 应用于spice电路仿真程序中的铁电电容行为模型 | |
CN101261658B (zh) | 基于电压延迟行为的铁电电容行为仿真用等效电路 | |
CN102930762A (zh) | 一种三维混沌电路 | |
CN205453724U (zh) | 一种多功能通用混沌电路系统 | |
CN112865109A (zh) | 一种数据驱动型电力系统潮流计算方法 | |
CN110417022A (zh) | 矩阵乘法运算提取雅可比元素的配电网三相潮流计算方法 | |
CN102521466A (zh) | 一种扩散电阻的仿真方法及系统 | |
CN211506501U (zh) | 一种浮地型磁控忆阻模拟器 | |
CN106326537A (zh) | 一种计算时域同时开关噪声的方法以及装置 | |
CN103390251A (zh) | 一种电力系统状态估计的量测量权重设置方法 | |
CN107658880B (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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20101013 Termination date: 20130403 |