求解类拉普拉斯方程的径向基函数配点法

建立了径向基函数配点法求解类拉普拉斯方程定解问题的方法.将求解域依据系数张量为分片常量来分解为若干子域.在每个子域上分别利用所布置的中心点建立用径向基函数表达的近似待解函数.在每个子域内及子域边界与外边界重合部分的配置节点上分别利用类拉普拉斯方程和定解条件建立近似解函数的待定系数满足的配点方程组,在相邻子域的分界线的配...     

应用DDM/FEM方法研究开口腔体电磁散射特性

将非重叠区域分裂法（DDM）与矢量有限元法（FEM）相结合对无限大接地的三维开口腔体的电磁（EM）散射问题进行分析，将原尺寸较大的腔体分成若干较小的子域，在各子域内推导出此边值问题的等价泛函，在此基础上应用矢量有限元法进行分析。通过传输条件实现各子域间的耜合，在腔体口面应用边界积分方程（BIE）进行描述，这样，减少了程序运行时对计算机内存的需求，易于对电大尺寸腔体的电磁散射问题分析计算，所得数值结果与有关文献的结果一致。从而验证了本算法的正确性。     

应用DDM/FEM/BEM混合法计算各向异性介质柱电磁散射

应用重叠型区域分解法（DDM）结合有限元（FEM）和边界元法（BEM）计算二维各向异性介质柱电磁散射.对介质柱外的无限大区域采用边界元法分析,将介质柱所在区域分解为若干个重叠的子域,每个子域用有限元法分析,各子域间通过传输条件进行耦合.为了提高计算速度,引入了多波前法求解有限元方程,并用内观法结合多波前法解有限元和边界...     

带环形弹性顶盖圆柱形储液罐的流-固耦合特性

研究带有环形弹性顶盖的刚性圆柱罐在充满液体时的流-固耦合振动特性.将液体域分割成一个圆柱形子域和一个圆环柱形子域,从而保证在2个子域内的液体速度势函数为C1类连续函数并具有连续的边界条件.利用分离变量法和叠加原理求得各子域内液体速度势的解析解,利用环形顶盖的干模态来展开环形顶盖的湿模态.对液体子域间的界面连接条件和自由表面波条件作级数展开:沿液体深度方向作Fourier展开,沿径向作Bessel展开.截断级数得到环形顶盖与液体耦合振动的特征方程.比较了数值结果与有限元法结果,显示出很好的一致性.     

非线性机械系统动力分析的辛时间子域-迭代法

建立了用于非线性机械动力响应分析的辛时间子域-迭代法.首先,在任一时间子域内,将原非线性动力方程分解为线性和非线性两部分.对于线性方程,采用基于相空间非传统Hamilton型变分原理的辛时间子域法求解.将线性响应代入非线性项并视其为载荷,重复使用辛时间子域法迭代求解至一定精度,各解叠加即得原方程在当前时间子域的数值解.然后,将这个时间子域的末端值作为下个时间子域的初始值,进行下一个时间子域的计算,如此递推下去直到最后一个时间子域.结果表明,该方法具有较高的计算精度和计算效率.     

二阶时域波动方程的无网格方法求解

将径向基函数配点型无网格方法引入二阶时域波动方程的求解中,方程的空间导数采用径向基函数逼近,时间导数采用Crank-Nicolson方法离散,对应的边界条件直接施加在离散的边界数据点上.采用该方法对二维非规则求解域内的波传播问题进行了数值计算,并与有限元计算结果进行了对比分析.结果表明:基于径向基函数配点的无网格方法不但形式简单、易于实施,而且能够有效解决复杂求解域高维的波动问题.     

一种基于SVD-CBFM和RACA的单站RCS快速求解方法

在奇异值分解特征基函数法(SVD-CBFM)的基础上提出一种快速求解目标单站RCS的有效数值方法.为了降低入射激励的数目,该方法考虑了子域间的相互耦合作用,计算出各子域的次要特征基函数(SCBF).采用再压缩自适应交叉近似(RACA)算法对各子域的特征基函数(CBFs)进行压缩,加速CBFs的生成.同时运用RACA算法填充远场区阻抗矩阵,从而进一步提高SCBF和缩减矩阵构造过程中的矩阵矢量相乘速度.数值算例验证了该方法的精准性和有效性.     

平面黏弹性力学问题的插值型重构核粒子法

插值型重构核粒子法的形函数结合了包含Kronecker delta特性的简单函数和由基函数向量采用重构条件构造的增强函数,并且具有点插值特性和不低于核函数的高阶光滑性.该方法可直接施加本质边界条件,同时也保证了较高的计算精度.基于插值型重构核粒子法,文章提出了一种求解平面黏弹性力学问题的新方法.采用弹性-黏弹性对应原理和Laplace变换,将黏弹性问题转化为Laplace域内的准弹性问题,并采用插值型重构核粒子法进行求解,然后借助Laplace数值逆变换求得黏弹性问题的解.数值算例验证了本文所提方法的有效性.     

向量三阶拟线性边值问题的奇摄动

本文考虑一类向量三阶拟线性边值问题。在适当条件下,通过构造边界层函数,求得高阶渐近展开,然后利用对角化方法,证明了其解和高阶渐近解的误差估计。     

线性广义特征值问题在模态子空间中的摄动解

介绍一个求解形如（［Ａ０］＋［△Ａ］）［Ｚ］＝（［Ｂ０］＋［△Ｂ］）［Ｚ］［Ａ］的广义线性特征值问题的摄动方法，这一摄动解是在由特征值方程［Ａ０］［Ｘ］＝［Ｂ０］［Ｘ］［Ａ］的部分模态向量张成的子空间内求得．应用这一方法，可将特征值方程的求解转化为线性代数方程组的求解．文中还讨论了有关重特征值的摄动求解方法．     

Study on the EM Scattering of Rivets on the Plate

The electromagnetic (EM) scattering by rivets on the conducting plate is studied for the first time by using electric field integral equation (EFIE) in conjunction with the moment method. The surfaces of the rivets and the plate are partitioned into triangular cells, the current distribution on the patches is represented by sub-domain type basis function, the EFIE is translated into matrix equation by the Galerkin method, then the current coefficient is obtained. The results show the properties of radar cross section (RCS) varying with the incident angle when there are rivets on the plate.     

短峰波与圆弧型贯底式多孔介质防波堤相互作用研究

将比例边界有限元法(SBFEM)推广到短峰波与圆弧型贯底式多孔介质防波堤相互作用研究中.将圆弧延伸构建一个虚拟圆柱面,对于圆弧段与圆弧延伸段,分别采用不同的孔隙影响系数.该圆将整个流场划分成一个无限子域和一个有限子域,并利用变分原理在各个子域上推导了SBFEM的基本方程.然后,SBFEM针对无限域和有限域分别采用汉开尔函数和贝塞尔函数作为基函数来求解对应域的解,并在公共圆边界处加以匹配.数值计算验证了该方法是一种利用很少单元离散虚拟圆便能得到精确结果的高效方法.进一步给出了不同波向入射、相对波数、防波堤不同孔隙影响系数、不同张角等参数对防波堤波浪荷载、波浪爬升的影响,这为该类防波堤水动力分析和结构设计提供了有价值的参考.     

区域分裂方法应用于屏蔽门的电磁泄漏分析

电大尺寸和复杂物体的电磁泄漏分析是计算电磁学的一个重要的研究课题，区域分裂方法（DDM）作为微分方程数值求解的新技术十分适宜求解电大尺寸的电磁场问题。提出了一种基于区域分裂方法和有限元方法（FEM）的混合算法来分析电磁防护中的电磁屏蔽门问题，在屏蔽门的边界和区域分裂的虚拟边界上分别利用吸收边界条件和传输边界条件，具有良好的收敛性质，数值结果表明了这一混合算法的有效性，同时区域分裂方法十分适合于计算机的并行计算，所以这里给出的方法适合于计算电大尺寸物体的电磁计算问题。     

不连续位移法解弹性力学平面问题

本文讨论了间接边界元法中的不连续位移法的一些问题;给出了不连续位移基本解和开尔文基本解的关系;分析了传统的不连续位移法解弹性力学平面问题精度低的原因,提出了将不连续位移作用在域外的断续附设边界上的方法,理论分析和算例表明该方法能大大提高精度,同时还能应用于含有夹层的由多种材料组成的地基计算.     

TE波时域棱边有限元方法：线磁流辐射与散射分析

讨论有耗介质TE波时域棱边有限元方法,导出电场矢量波动方程边值问题的弱解形式,应用棱边基函数给出单元矩阵方程,通过组合获得时域全域矩阵方程,详细讨论棱边有限元组合中符号函数的作用和累加填充步骤.给出了激励矢量中线磁流的加入以及棱边有限元的定量验证,分析了线磁流照射下有耗介质物体散射.     

点配法应用于矩量法分析的若干问题及改进

在矩量法分析中,点配法由于简单、计算量小和物理意义明确而获得了广泛的应用,其缺点是解不稳定,收敛性差。本文从双位方程出发,讨论了线天线匹配点位置的选择对解的影响,介绍了一种改善计算精度和收敛特性的方法,给出一种节省计算机CPU时间的技术。     

二维开口腔体填充多层各向异性介质TM波散射特性分析

应用区域分解法(DDM)结合有限元(FEM)和边界积分法(BIM)分析了填充多层各向异性介质二维开口腔体TM波散射特性.对腔体外区域采用边界积分法进行分析,将腔体内的每层介质作为一个子域用有限元法分析,各子域间通过传输条件进行耦合.作为算例,分别计算了腔体中填充各向同性和各向异性介质时的散射截面,数值结果表明了该方法的有效性.采用这种技术,大大地减少了对计算机内存的需求.     

A-?方法在高频电磁场计算中的应用

基于节点的有限元方法具有网格剖分、构造高阶基函数容易的优点。由于在节点上定义场量,节点有限元更适用于多物理问题。但节点有限元方法直接求解电磁场会出现伪解,场量在不均匀介质中不连续等问题。基于节点的A-φ方法可以有效避免传统节点有限元方法存在的问题。本文研究A-φ方法的工程应用,研究开域和闭域问题中如何设置关于矢势和标势函数的边界条件,特别是波导问题和理想导体球的散射问题,讨论了端口边界条件,辐射边界条件的使用方法。对于理想导体边界条件采用了阻抗边界条件,与端口条件配合,克服方程的奇异性。数值卖验比较分析了A-φ法节点有限元和棱边法有限元的计算结果,验证了A-φ法节点有限元的正确性和有效性。