如何计算ComsolMultiphysics®中的电容矩阵

作者图像

经过Sven Friedel

2017年6月7日

ComsolMultiphysics®软件中的电容计算似乎很容易。如果您只有两个导体,则配方很简单:拿一个导体并将其设置为接地,将另一个导体设置为终端,然后计算溶液。然后,一个内置的变量可提供电容。但是,如果您有两个以上的导体,例如触摸屏,传输线和电容传感器,该怎么办?如果您丢失了标准教科书术语,请遵循计算电容矩阵的工作示例。

什么是自我掌握?

电容是系统存储电荷的能力。它可以通过人体需要的充电量来定义,以便在接地参考电位上升高1伏。在线性系统中,这是

Q = C \ CDOT V

如果Q是电荷,则V是相对于地面的电势差,而C是电容。

在进入多个导管系统之前,请记住,从定义上讲,即使是一个孤立的导体也具有电容,相对于无穷大的接地球形壳的定义。对于指挥领域的情况,这种自我电力是

c = 4 \ pi \ epsilon_0 r

我们可以使用此公式来计算地球的自par体 - 大约是710 microfarad。

人体也可以充电,如这里所示

因此,人体也表现出自par(也称为身体电容)。根据姿势和周围区域的不同,身体电容在100 picofarad的范围内,甚至可以在人类中产生刺痛感。例如,您可以在早上梳理头发时轻松地将身体电容充满几千伏。在开始一天之前,请确保扎根!

比较相互和麦克斯韦电容矩阵

在典型的电气系统中,几个导体之间的电容是最感兴趣的。相互电容,也称为寄生或流浪电容,是所需的或不希望的电容(电荷的堆积),发生在两个电荷持有物之间。如果将带电的物体带到另一个对象附近,则由于静电诱导过程(不要与电磁诱导混淆),因此第一个对象上的电荷分布将发生变化。特别是在传输系统中,线之间的电容耦合通常是意想不到的,而且很麻烦,因为它会产生噪音。

相互电容示例的两个示意图。
屏蔽的三导电缆(左)以及地面板上的微带传输线(右)之间相互电容的典型示例(右)。从连续场模型到具有离散电容器的集总模型的过渡意味着将导体缩小到点,同时将电荷在其表面上移动到它们之间显示的电容器板的板。

为了方便起见,可以安排n个导体系统的相互电容和矩阵形式的另一种基础:

\ begin {bmatrix}
c_ {m,11}&c_ {m,12}&\ dots&c_ {m,1n} \\
c_ {m,21}&c_ {m,22}&\ dots&c_ {m,2n} \\
\ vdots&\ vdots&\ ddots&\ vdots \\
c_ {m,n1}&c_ {m,n2}&\ dots&c_ {m,nn}
\ end {bmatrix}。

当您将物理系统降低到离散元素网络时,该矩阵的系数(也称为部分电容或总电容)被使用在电路模拟器中。

在现场理论中,另一种矩阵形式更常见:麦克斯韦电容矩阵。由于名称是如此相似,并且系数不完全相同,因此重要的是要了解互助和麦克斯韦电容矩阵之间的关系。麦克斯韦电容矩阵描述了i的电荷之间的关系Th导体到系统中所有导体的电压。

\ begin {pmatrix}
q_1 \\
q_2 \\
\ vdots \\
q_n
\ end {pmatrix} =
\ begin {bmatrix}
c_ {11}&c_ {12}&\ dots&c_ {1n} \\
c_ {21}&c_ {22}&\ dots&c_ {2n} \\
\ vdots&\ vdots&\ ddots&\ vdots \\
c_ {n1}&c_ {n2}&\ dots&c_ {nn}
\ end {bmatrix}
\ begin {pmatrix}
v_1 \\
v_2 \\
\ vdots \\
v_n
\ end {pmatrix}

麦克斯韦电容矩阵系数C_ {11}可以通过测量导体1上的电荷来确定,而只有电势V_1 = 1并且所有其他电极均接地。因此,矩阵通常也称为接地电容矩阵。它的倒数c^{ - 1},称为弹性矩阵。

我们还可以通过总结自相互电容的贡献来计算导体1上的总电荷,如下所示。

\ begin {matrix}
q_1&=&c_ {m,11} v_1+c_ {m,12}(v_1-v_2)+c_ {m,13}(v_1-v_3)\\
&=&(c_ {m,11}+c_ {m,12}+c_ {m,13})v_1-c_ {m,12} v_2-c_ {m,13} v_3)
\ end {matrix}

对于具有N导体的系统,相互关系和麦克斯韦电容矩阵之间的关系是

\ begin {pmatrix}
q_1 \\
q_2 \\
\ vdots \\
q_n
\ end {pmatrix} =
\ begin {bmatrix}
\ sum _ {} {^n_ {i = 1}} c_ {m,1i}&-c_ {m,12}&\ dots&-c_ {m,1n}
-c_ {m,21}&\ sum _ {} {^n_ {i = 1}} c_ {m,2i}&\ dots&-c_ {m,2n}
\ vdots&\ vdots&\ ddots&\ vdots \\
-c_ {m,n1}&-c_ {m,n2}&\ dots&\ sum _ {} {^n_ {i = 1}} c_ {m,ni}
\ end {bmatrix}
\ begin {pmatrix}
v_1 \\
v_2 \\
\ vdots \\
v_n
\ end {pmatrix}

您可以通过其负面的非对角元素来轻松地告诉麦克斯韦电容矩阵。

工作示例:两个球的相互电容

现在,我们已经对术语有了明确的定义,让我们探索计算Comsol多物理学中指导器的任意导体系统的电容矩阵的容易。为了使我们的脚下保持坚实的地面,我们将从具有已知分析解决方案的系统开始。(我是否提到我喜欢分析解决方案?

在我们的情况下,我们可以使用由两个Radii的导电球组成的系统Ab,距离分开C和无穷大的地面参考。

两个导电球的例证。
自麦克斯韦时代以来,就已经知道了这种系统的封闭式表达式。我指的是de Queiroz(2003)和Lekner(2011)的两份出版物。三个麦克斯韦电容矩阵的表达式是

c_ {11} = fab \ sinh u \ sum_ {n = 0}^\ infty [a \ sinh nu+ b \ sinh(n+ 1)u]^{ - 1}
c_ {22} = fab \ sinh u \ sum_ {n = 0}^\ infty [b \ sinh nu+ a \ sinh(n+ 1)u]^{ - 1}
c_ {12} = - \ frac {fab} {c} \ sinh u \ sum_ {n = 1}^\ infty [\ sinh nu]^{ - 1}

在哪里

\ cosh u = \ frac {c^2-a^2-b^2} {2ab} \ quad和\ quad f = 4 \ pi \ epsilon_0 \ cdot 1 [m]

您可以使用总和运算符轻松将这些表达式声明为comsol多物理中的变量:

屏幕截图,显示表达式为comsol多物理学中的变量。

在n上进行参数扫描时,我们发现当球体彼此之间不太接近时,该系列会迅速收敛。对于给定的一组参数,我们可以将n安全设置为10:a = 0.1,b = 0.3,c = 0.5。

为了计算电容矩阵静电界面,我们将终端条件设置为一个势为1 V的球体。

带注释的屏幕截图显示了Comsol多物理学中的终端条件设置。

接下来,我们复制该功能,将其应用于第二个球体,然后将终端名称设置为2。为了计算电容矩阵,我们需要在端子上应用不同的电压或电荷模式。出于教学原因,我们将讨论传统的手动终端扫描,然后再引入comsol®软件5.3版发布的新的,更快的技术。尽管新技术在一系列非常常见的问题中更快,但手动方法更为笼统。

手动终端直接在静电界面。

屏幕截图显示了Comsol多物理学中的静电接口。

宣布扫描参数名称在里面全局参数部分 (portname默认情况下),然后您可以运行参数扫描portname

带有参数扫描节点的屏幕截图突出显示了模型树中的屏幕截图。

在模型中,comsol多物理软件将一个终端设置为1 V,并且在扫描期间将所有其他终端设置为接地,从而导致以下两个解决方案:

Comsol多物理学中电容矩阵模型的解决方案。

您可以使用结果>全球矩阵评估在不同符号中提取电容矩阵,包括麦克斯韦电容矩阵和相互电容矩阵。

带注释的屏幕快照,显示Comsol多物理中的电容矩阵。

在这个简单的例子中,他们的关系是

\ begin {bmatrix}
c_ {m,11}&c_ {m,12} \\
c_ {m,21}&c_ {m,22}
\ end {bmatrix}
=
\ begin {bmatrix}
C_ {11}+C_ {12}&-C_ {12} \\
-c_ {21}&c_ {22}+c {21}
\ end {bmatrix}

如果要设置电荷端子而不是电压端子,则主要解决方案是逆电容矩阵。有一组转换可以帮助您将电荷转换为上述矩阵。

使用固定源扫描和边界元素方法增长速度

在Comsol Multiphysics版本5.3中,我们引入了许多强大的新建模方法。与电容矩阵计算特别相关的一项功能是新的固定来源清扫学习步骤。

与手动终端扫描相反,该终端使用参数扫描portname,这项新技术解释了以下事实:在静电系统上应用不同的电荷或电压模式不会改变基础FEM方程的系统矩阵,而只会改变其负载。这意味着矩阵只需要一次倒置一次,并且可以在所有其他负载情况下重复使用。这种方法可以大大减少计算时间,尤其是当需要大量的终端数量或其他参数扫描(例如,对于几何形状)时。

即使对于中等数量的终端,触摸屏仿真器仿真应用程序,速度增长令人惊讶:对于此型号,我的机器上可以达到7.8倍的DOF!

新的电容位置传感器教程模型使用固定来源清扫学习步骤。

固定的源扫描也更容易设置。无需激活手动终端并定义一个portname变量和参数扫描。您需要做的就是选择一个学习步骤。默认情况下,该研究将贯穿所有终端。另外,您可以定义要涵盖的指定资源。

comsol多物理学中固定源扫描设置的屏幕截图。

如果固定源扫描如此强大,那么为什么我们要保持传统方法呢?

在某些情况下,我们感谢对系统矩阵的重新计算;例如,在非线性或多物理问题中,或者应针对每个终端配置对网格进行调整。在这种情况下,手动终端扫描是一种更好的方法。

comsol®软件的5.3版带有的另一个功能功能是静电中的边界元素方法(BEM)。与FEM相比,在所有域中需要网格(包括周围空气域)的网格,BEM避免在无限的空隙中划分,从而减少了DOF的数量。您可以在静电和电容矩阵计算中了解更多有关将BEM组合在使用BEM教程对电容位置传感器进行建模

等不及要了解有关新BEM实施的更多信息?阅读有关该方法如何有助于简化腐蚀建模的信息以前的博客文章

关于电容矩阵计算的总结想法

在这篇博客文章中,我们研究了电容矩阵的计算,讨论了不同的术语,并为具有封闭形式解决方案的众所周知问题提供了数值解决方案。虽然此处介绍的简单模型可以为您自己的模型提供参考,但版本5.3中发布的新功能可以帮助您更有效地创建具有许多终端和其他参数扫描的大型模型。

其他资源

参考

  1. De Queiroz,A.C.M。,2003年,“电容计算“。
  2. Lekner,J。,2011年,“两个球体的电容系数”,静电学杂志69(1):11-14。

评论(11)

发表评论
乐动体育app
加载...
李·温宁(Li Wenzheng)
李·温宁(Li Wenzheng)
2018年3月5日

嗨,弗里德尔,
感谢这个有用的博客!我有一个问题:
有什么办法可以按时间相关模式计算电容矩阵?因为我的模型是时间依赖的其他东西(神经活动)。我已经尝试过,但是当研究时间依赖时,没有电容矩阵的表达。
期待您的帮助!
问候,
Wenzheng

Sven Friedel
Sven Friedel
2018年4月3日

亲爱的温宁,
电容计算仅用于固定和频率依赖性研究。
想象一下RC情况。R和C可以从阻抗或入场矩阵的真实和虚构部分中唯一确定。在时间相关的情况下,此信息隐藏在u(t)和i(t)曲线中。看:https://www.comsol.ch/model/transient-modeling-of-a-capacitor-in-a-a-circuit-12695
在系统缓慢变化的情况下,您可以通过对参数p的依赖性替换对时间t的依赖性,然后可能提取C(p)。
斯文

大卫·里斯(David Reens)
大卫·里斯(David Reens)
2018年5月1日

亲爱的Sven,

多么出色的模型,对我的研究非常有用,谢谢。关于您从集团/香料/相互电容矩阵到Maxwellian的转换,我对我有些奇怪。我认为,由于地面的处理方式,前者的维度应该比后者少。我相信Maxwellian矩阵应该具有总和为零的行和列,不是吗?还是在此处留下麦克斯韦矩阵的最后一行和列的约定,以便给出vec {q} = matrix(c)*vec {v},其中电压和电荷向量排除了地面?

谢谢,
戴夫

Sven Friedel
Sven Friedel
2018年5月2日

亲爱的戴夫,
实际上,我们遵循公约vec {q} = matrix(c)*vec {v},其中排除了地面。如果您查看图。1上面,您会看到一个由3个导体组成的系统(+一个无穷大的地面或周围的盾牌)。这使得总共n = 4杆。电容矩阵具有N*(N-1)/2系数,在这种情况下6,可以通过对称(N-1)X(N-1)矩阵表示。
此致
斯文

岩岩
岩岩
2018年6月25日

我有四个核心电缆,我想在频率上获得电容矩阵。
我将公共屏幕设置为接地,并将四个内核设置为不同的端子,将它们列出1至4。并将它们设置为电压端子类型,电压为1V。
我想使用PortName的参数扫描。但是,我对如何在全局参数定义中定义portname感到困惑。我可以给portname的名字,但是那是什么表达式?现在,我将其设置为portname,并得到黄色警告。仿真似乎运行了,但我一无所获。
另一个问题是如何在频率范围内获得结果?
谢谢

Gleb Kuman
Gleb Kuman
2019年5月17日

亲爱的,斯文
非常感谢您提供如此不错的博客文章!

我对您的模型有一个问题:
- 根据等式Q1 = CM11V1 + CM12(V1 - V2) + ETC,如果我们从1st旁边的模型中删除所有导体,则获得单个导体Q1 = CM11*V1容量的简单方程。
- 因此,第一球CM11的容量等于5.4pf(形式相互矩阵)应等于单个球体的容量,根据c = 4*pi*pi*eps0*r等于11.3pf。
- 您能说出这种差异的原因吗?为什么CM11不等于单个球的容量?
谢谢

Sven Friedel
Sven Friedel
2019年5月17日

相互电容矩阵的对角元素确实与各自身体的自相射率并不相同。让我们参考原始定义:“电容描述了人体需要在1伏上向接地参考潜力提高电势所需的电荷。”
如果除了无穷大的地面外,还有另一个机构(也接地),那么为了将其潜力提高到1 V而需要存储的充电量较小。电容矩阵的对角元素是在另一个机构的“存在”中定义的,这也很明显,从分析描述C11 C11是A和B的函数。对于b –> 0,您可以获得半径为a的球体。

lim_b-> 0(f*a*b*sinh(u)*sum(1/(a*sinh(n*u)+b*sinh((n+1)*u)),n,0,n))= 4*pi*epsilon0_const*a

一探究竟!
斯文

Gleb Kuman
Gleb Kuman
2019年5月17日

亲爱的,斯文
非常感谢您的全面和快速的回​​应。
现在一切都清楚了!

问候,格莱布

奥斯卡·迪亚兹(Oscar Diaz)
奥斯卡·迪亚兹(Oscar Diaz)
2021年2月4日

嗨,弗里德尔,
关于相互电容计算的非常有趣的博客。我运行两种方法都解释了:1)手动终端扫描和2)固定源扫描。使用第一个,所得的求解模型更重,并且可以在没有问题的情况下获得相互电容矩阵。对于第二个,我无法提取矩阵,因为每次尝试进行全局矩阵评估时,输出都只会返回“ NAN”单元。
你能建议我想念什么吗?
谢谢!

Sven Friedel
Sven Friedel
2021年2月4日 Comsol员工

尊敬的奥斯卡 - 我会补充您创建一个常规的支持案例。在那里我们可以为您提供最好的帮助。

Prasanna Kumar Routray
Prasanna Kumar Routray
2021年7月29日

嗨,Sven,
您有本教程的.mph文件吗?

探索comsol乐动体育赛事播报博客