|

楼主 |
发表于 2011-4-4 11:32:07
|
显示全部楼层
V1。2版本发布:
本程序采用更好的算法,直接采用电磁场论公式计算,速度快。
电感量计算.rar
(4.97 KB, 下载次数: 819)
快速算法说明
椭圆积分的数值计算:
·第一类椭圆积分 K(k)=∫0π/2sqrt(1-k2sin2θ)-1dθ
·第二类椭圆积分 E(k)=∫0π/2sqrt(1-k2sin2θ)dθ
·K(k)和E(k)可以利用“几何算术平均值”法计算。迭代过程如下:
a0 = 1, b0 = sqrt(1-k2), c0 = k
an = (an-1 + bn-1)/2
bn = sqrt(an-1·bn-1)
cn = (an-1 - bn-1)/2
C = 20c02 + 21c12 + 22c22 + … + 2NcN2,式中n=1,2,3,…N
最后得:K(k) = π/2/an, E(k) = K(k)*(1-C/2)
以上每迭代计算一次,精度的有效数字的位数提高一倍,收敛速度很快。一般只需几次迭代计算就可以得到足够精度。
·利用纽曼公式计算线圈C1与C1之间的互感
M = ∫C1∫C2dl2·dl1/|r2-r1|
对轴对称两圆圈积分后得到:
M = uRD [ (1-k2/2)K(k) - E(k) ]
式中:R=大圈半径,b=小圈半径/R,h=两圈距离/R,u = 4π1e-7
D = sqrt( (1+b)2 + h2 ),k = 2*sqrt(b)/D
关键要计算出出椭圆积分K(k)和E(k)才能计算出电感量M
如果有N圈,那么需计算任意两圈之前的互感量的总和,总共有N*(N-1)个互感。
然后,计算N线圈的各圈单圈自感,共N个单圈自感。
再计算出导线的内自感,其值为0.05*导线长度,单位是uH
那么N圈的总电感为 L = 各圈互感 + 各圈自感 + 导线内自感 |
|