二维有限差分析是求解两个变量的拉普拉斯方程的一种近似方法,这种方法的要点如下:
在平面场中,将平面划分成若干正方形格子,每个格子的边长都等于h ,图13-10表示其中的一部分,设0点的电位为V 0,0点周围方格顶点的电位分别为V 1、V 2、V 3和V 4。
现在来推导一个用V 1、V 2、V 3和V 4表示V 0的公式:
图13-10
已知平面场的电位满足两个变量的拉普拉斯方程:
0222
2=∂∂+∂∂y
V
x V 其中
h
x
V x
V x V x x V c
a
∂∂-
∂∂≈⎪⎭⎫ ⎝⎛∂∂∂∂=
∂∂0
22
但是
h V V x V
h V V x
V c a
30
01 ,-≈∂∂-≈
∂∂ 所以
2
30013
0010
2
2h V V V V h h V V h V V x V +--≈--
-≈∂∂
同理
2
4
0020
2
2h
V V V V y V
+--≈
∂∂ 将上面两个方程相加一起得:
042
43212222=-+++≈∂∂+∂∂h
V V V V V y V x V 由上面方程推出:)(4
1
43210V V V V V +++≈
(13.47)
该式说明0点的电位近似等于相互垂直的方向上和0点等距离的四个点上的电位平均值,距离h 愈小则结果愈精确,方程(13.47)是用近似法求解两个变量拉普拉斯方程的依据。
然而,V 0和V 1、V 2、V 3、V 4都是未知值,这种情况下需要按照方程(13.47)写出每一点的电位方程,然后求这些方程的联立解。
求解时较简便的方法是选代法,这种方法可求出平面场中各点电位的近似值。
图13-11表示一个截面为正方形的导体槽,槽的顶面与侧面相互绝缘,顶面的电位为
V 0,侧面与底面的电位都等于零。
为了求出槽中各点的电位,将槽分成十六个相同的方格,这些方格在槽中共有九个顶点。
用V 1、V 2,…,V 9表示各顶点的电位。
求解步骤如下:
图13-11
第一步,假设某点的电位为某值,称为某点的原始电位,原始电位等于多少并不影响最后的结果。
如果原始电位选择得当,则计算步骤会得到简化。
第二步,根据原始电位,利用式(13.47)求出每点周围四个点电位的平均值,电位平均值一般不等于电位的原始值,将平均值代替原始值就得到每点电位的第一次选代值。
然后根据第一次选代值求出每点周围四个点电位的平均值,如果平均值不等于第一次选代值,就将平均值代替第一次选代值,得到每点电位的第二次选代值。
第三步,利用式(13.47)对每点电位进行选代,一直到每点的电位与它的周围四个点的电位平均值相差在允许范围内为止。
【例13.1】在图13-12中,设V=100,试用选代法求方格顶点上的电位。
图13-12
解:设九个顶点的电位分别用V 1、V 2、…、V 9来表示。
第一步:设每点的原始电位都等于零。
第二步:根据原始电位利用公式,)(4
1
43210V V V V V +++≈
,求出各点的周围电位的平均值为:25)000100(4
1
321=+++===
V V V 。
其余各点周围电位的平均值都等于零。
然后将所得的平均值代替原始值,得到第一次选代值。
第三步,根据第一次选代值,求出各点周围电位的平均值为:
3.31)0025100(41
1=+++=V 5.37)25025100(41
2
=+++=V 3.313
=V
3.6)00025(4
1
654
=+++===V V V
0987
===V V V
然后将所得的平均值代替第一次选代值,得到第二次选代值。
第四步:根据第二次选代值,求出各点周围电位的平均值为: 36)03.65.37100(41
1=+++=V 2.42)3.63.313.31100(41
2
=+++=V 363
=V 4.9)003.63.31(41
4=+++=V 5.12)3.603.65.37(41
5
=+++=V 4.96
=V
6.1)0003.6(4
1
987
=+++===V V V
然后将所得的平均值代替第二次选代值,得到第三次选代值。
按照同样方法对每一点进行选代,结果如表13-1,可以看出,步骤18以后,各点的电位收敛于某固定值。
利用有限差分法求解电位方程时,需要进行大量的计算,本题解仅求九个点的电位,计
算工作量已可观、如果求电位的点数目很大,则必须用电子计算机进行计算。
【例13.2】
如图13-13所示表示四个不同形状的电极围成一个不规则槽,各电极的电位如图所示。
槽的截面共分成14个相同的方格,试用选代法求出每个方格顶点的电位。
图13-13
解:第一步,设每点的原始电位都等于零。
第二步,根据原始电位,求出各点的周围电位的平均值。
10)002020(41
1=+++=V 5)00200(41
2=+++=V 5)00200(41
3=+++=V 5)00020(4
1
4
=+++=V 065
==V V
然后将所得的平均值代替原始值,得到第一次选代值。
第三步,根据第一次选代值,求出各点周围电位的平均值为: 5.125041
)552020(411==+++=V 75.83541
)052010(412==+++=V 25.52541
)00205(413==+++=V 25.73041
)001020(414==+++=V 5.21041
)0055(415==+++=V
25.1)0050(4
1
6
=+++=V 然后将所得的平均值代替第一次选代值,得到第二次选代值。
第四步,根据第二次选代值,求出各点周围电位的平均值为:
145641
)25.775.82020(411==+++=V 1025.4041
)5.225.55.1220(412==+++=V 5.73041
)25.102075.8(413==+++=V 75.83541
)05.25.1220(414==+++=V 3.425.1741
)025.175.825.7(415==+++=V
94.175.74
1
)0025.55.2(416
==+++=V 然后将所得的平均值代替第二次选代值得到第三次选代值,这样继续重复下去,一直计算到
满意的程度为止。