二维均匀电场直流电阻率法的边界单元法.docVIP

二维均匀电场的直流电阻率法的边界单元法 一、边界元法计算均匀场中水平地形条件下二维不均匀体的异常 设在均匀场中。有一电阻率为的二维不均匀体，由于地面水平，故可由镜像法原理求空间任一点的电场（如图(1.1)所示） 上述电场的边值问题满足： 地表 (1.1) 地表 （1.2) (1.3) (1.4) (1.5) （一）积分方程的建立 由格林公式 图.1.1区域与边界 (1.6) 设为二维拉普拉斯方程的基本解 (1.7) 式中r为中任意点至电场计算点P的距离 则 (1.8) 为以p为中心的函数，(1.1)、（1.3)式代入(1.6)式左边 (1.3) 式中是p点的电位，为p对区域的张角。格林公式的右边为 而 由于边界在无穷远处，我们可近似将p点视作的圆心，因而在上，下列关系式成立 (1.10) 式中为p点的x坐标 (r,y)为失径与坐标y的夹角 (n,x)为法向与坐标x坐标的夹角 (1.10)式代入 式中为点p点的初始电位 于是当p点在内时，则(1.6)式变为 (1.11) 当P在地面时 (1.12) 由此可见，只要知道边界上的u和，便可由积分方程(3.1.12)求出地表的u.。 （二）边界上的u和的求取 当p点位于边界上时，(1.12)式变为 (1.11’) 该积分方程中包含函数u和，求解时，方程数为n个，而未知数有2n个，故必须导出另一组积分方程 在区域，从出发，可写出如下格林公式 (1.13) 式中出现负号是因为这里定义的法向方向，与前面定义区域时同向所致。 仍定义：当p位于边界上时，p对区域的张角为，则p对区域的张角为。于是与，中同样方法，得上的积分方程为 将边界条件(1.3)、(1.4)代入，得 (1.14) 联立(1.12)、(1.14)式，即可求得边界上的u及。 （三）边界单元法 用边界单元法求解(1.11)、（1.14)式。 用n个节点对边界进行剖分，上的剖分与完全对称，将(1.11)中的边界积分分解为诸单元的积分和，当p点位于第i个节点上时，则有 式中为第i节点的正常场值： 为i节点的x坐标 1、单元分析 为求解边界积分，我们必须求解在及边界上各单元的积分。我们先看在边界上单元两端节点编号为j,k,其坐标为，的单元（见图1.2）。 图1.2 首先我们作如下坐标变换，用形函数，表示单元中任意点的坐标x,y, 其中，是的数， 且 在j点, k点, , 可见，，是x,y的线性函数。 由于单元一般取得很小，故可假定电位u在各单元上是线性变化的，即单元上的u可表示为 式中u,j为节点j上的u; 为节点k上的 于是(1.15)式中第一项在边界上的积分 (1.16) 式中 (1.17) (1.17)式积分可用m点高斯求积公式计算 式中高斯求积公式所用的点数 i点至单元q点的距离 加权系数 L 单元长度 Q点的坐标为 同样(1.15)式中第二项是在的镜像上进行的 图1.3 (1.18) 式中 (1.15)式中第三项可写为 (1.19) 式中 (1.15)中第四项为 (1.20) 式中 2、总体结合 (1.16)、（1.18)、(1.19)、(1.20)代入(1.15)中，得 (1.21) 式中是j节点两侧单元的之和。 对上每个节点都有如上一个方程，由全部n个节点得一线性方程组 (1.22) 式中 对(1.14)式作同样处理，得 (1.23) 联立(1.22)、(1.23)式得