第8章_边界元法学案.pptVIP

（2）线性单元 1）主对角线元素 采用线性单元离散曲线边界，相当于用N条直线段去逼近曲线。设定节点为线段的端点，显然节点处曲线不再光滑，如图8-8所示。根据式（8-51）和式（8-52），对角线元素为 和 式中，ci=1-θ/2π，θ是节点i所张的平面角。由于单元中r0与en相互垂直，经过与式（8-54）和式（8-56）相类似的推导，可以得出 和 2）非对角线元素 此时 根据式（8-47）和式（8-48），并注意到l=Lj（1+ξ）/2，所以有 同理可得 和 式中，r是ξ的函数， 如上所述，应用高斯求积公式，即可算得以上各系数矩阵的元素值。 8.5 典型算例 仍以长直接地金属槽中电位和电场强度分布的计算为例，具体阐述二维边界元法的实施及其计算精度的评价。该金属槽如图8-9所示，侧壁和底面电位为零，顶盖电位的相对值为100。基于场分布的对称性，此二维电场的求解区域可进而压缩而归结为D/2。这样，对应于本问题的数学模型为 采用8.4节中线性单元的离散化模式，沿计算场域边界按逆时针方向编写单元的节点号，如图8-10所示。然后，基于线性插值函数，建立系数矩阵［H］和［G］，由所得边界元方程（8-38）解出边界上的未知量后，即可计算场域内任一点的电位与场强。以上采用线性单元边界元法的数值求解过程，可应用列于附录8.1的线性单元边界元法的通用计