水波传播的无网格数值模拟

在水波传播的数值模拟中采用了一种基于配点和径向基函数的无网格方法.采用Laplace方程的基本解作为径向基函数,将源点布置在模拟波浪场之外,沿边界布置配点而不是划分网格,从而自动满足控制方程,且不存在奇点,不需要求解积分方程.数值造波采用给定入射波面和速度势的方法,数值消波综合采用阻尼层消波和Sommerfeld辐射条件,非线性自由面的演化追踪采用二阶Taylor级数展开式.对线性规则波和非线性三阶Stokes波的模拟显示,数值结果与理论解吻合良好.表明无网格方法不但形式简单、计算速度快,而且稳定性和准确性令人满意,有望成为水波模拟问题的一种有效的数值方法.     

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

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

四阶常微分方程的Birkhoff配点法

提出求解四阶常微分方程的Birkhoff配点法.通过构造满足边界条件的Birkhoff插值基函数,得到具有稳定条件数的代数方程组.在数值算例中,通过与一类Legendre 配点法的数值结果进行比较.结果表明:Birkhoff配点法的有效性和高精度.     

样条虚边界元法的数值稳定性与误差估计

样条虚边界元法是针对传统间接奇异边界元法存在的问题而提出的一种半解析半数值方法。它既保留了边界元法的优点，也避开了求解奇异积分方程的问题，在试函数和权函数的选取方面也作出了改进，具有精度好、效率高等优点。本文主要针对弹性力学平面问题样条虚边界元法在数值稳定性与误差估计方面的问题展开讨论，获得了虚边界的布设规律及方法误差的直观度量，为该法的实际应用打下了更好的基础。     

笛卡尔切割单元法在湖泊二维流场模拟中的应用

针对湖泊几何边界复杂的特点,为了高效并准确地计算流场,采用了基于有限体积法的笛卡尔切割单元法.该方法在流场主要区域采用笛卡尔网格,而对与固体边界相交处的单元采用特殊的处理.这样既保留了笛卡尔网格求解速度快的优点,同时在边界处又保证了精度,避免了传统差分法处理该类问题需采用阶梯形网格逼近而导致的失真问题.通过与一个标准算例的贴体网格计算结果对比,显示计算结果一致.最后运用笛卡尔切割单元法,求解了杭州西湖的二维流场,结果表明该方法适用于具有复杂几何边界的湖泊流场求解,且精度较高.     

解三维温度场问题的边界单元法

本文用加权余量法,建立了各向同性体三维稳定温度场的边界积分方程和离散型方程,导出了边界单元法在稳定温度场中出现的两个奇异积分的具体解析表达式。此外,对全域积分和系数矩阵的对称化也作了一些工作。     

无网格局部Petrov-Galerkin法分析板弯曲的剪切自锁问题

用径向基函数构造无网格局部Petrov-Galerkin方法的形函数,插值函数具有Kronecker delta函数性质,因此可以很方便地施加本质边界条件.分析了板弯曲时剪切自锁现象产生的原因,利用无网格局部Petrov-Galerkin方法对两对边固支另对边简支中厚板的弯曲进行了分析和计算.发现无网格方法相对于有限元...     

利用点插值函数作为试函数的最小二乘配点法

采用径向基函数与多项式基函数作为耦合的基函数,并利用点插值法构造加权残值法中的近似试函数,试函数中的形函数具有狄拉克-δ函数性质,因此可以直接施加本质边界条件．利用这种试函数和采用最小二乘配点法求解了一维二阶微分方程和薄板的弯曲问题,并与理论结果进行对比;同时还检验了配点数以及节点支持域半径对计算精度的影响．数值结果表明：这是一种与单元划分无关的无网格方法,具有模拟简单,计算精度高,收敛快的优点．     

计算含多处裂纹结构基频的一种近似方法

将计算单裂纹结构基频的方法推广到求含多处裂纹梁弯曲振动的近似基频。通过在未损伤梁位移的基础上叠加一个多项式来构造裂纹梁的位移函数,它满足梁的边界条件和运动条件;应用Rayleigh方法,求出含多处裂纹简支梁的弯曲振动基频,得到了梁的近似基频的闭式表达式。与Ansys计算结果比较表明,该方法具有较高的精度。     

基于2.5维有限元算法的边界处理及网格划分

由于2.5维周期化技术能有效地处理列车振动波传播模拟过程中计算尺度与计算时间的矛盾,因而国内外学者在此理论基础上将FE(有限元)和粘滞边界以及IFE(无限元)等边界处理方法进行耦合.为讨论上述算法中有限元网格划分区域及单元大小对计算精度影响,通过运用软件Matlab编制程序读取Ansys中的模型信息,并将2.5维有限元和无限元边界耦合等算法在其中得以实现.在此基础上发现在相同的网格划分条件下,采用2.5维无限元边界较粘滞阻尼边界具有更高的精度,同时通过选取不同尺寸的有限元模拟区域及最大网格大小,将计算结果与经典解析解进行比较,并对上述网格划分参数的合理选择提出了建议.     

起伏地形下地下间断面被动源叠前逆时偏移成像

为了在起伏自由边界条件下应用逆时偏移成像, 提出一种边界元与有限差分混合的方法来计算地震波场的逆时传播。首先, 使用边界元法反传地表台阵数据, 可以得到地下一水平界面上的波场位移数据; 然后, 将此水平面作为有限差分方法的边界, 并反传此水平面处的波形数据, 获得此水平面下空间的地震波场。为了应用成像条件, 利用地震波的性质, 将地震波场分离为矢量P波和S波成分, 对分离出的P波和S波分量, 应用互相关成像条件可获得地下界面。通过对理论数据的测试, 验证了方法的有效性。     

闭曲面有限元网格生成的边界预调整方法

为使闭曲面相对其参数域具有边界完整性,确保波前法(AFT)生成网格的可靠性,提出了以参/虚边界补充为特征的闭曲面边界预调整方法,将闭曲面分割为一系列开曲面片,进而利用AFT法使网格在各开曲面片上分区域同时生成.该方法可使基于AFT法的曲面网格生成算法直接用于各种形式的闭曲面.     

汽车运动空气流场的有限容积法

介绍了汽车运动空气流场数值计算的有限容积法的主要内容，涉及基本方程，有限容积法的应用，交错网格系统边界条件 处理等，并给出了应用算例，建立了汽车运动流场数值计算的理论基础。     

The Boundary Contour Method for Half-Plane Piezoelectric Media

This paper presents a further development of the Boundary Contour Method (BCM) for half-plane piezoelectric media. Firstly, the divergence free property of the integrand of the half-plane piezoelectric boundary element is proven. Secondly, the boundary contour method formulation is derived and potential functions are obtained by introducing linear shape functions and Green's functions[1] for half-plane piezoelectric media. Finally, numerical solutions for illustrative example are compared with exact ones and that of conventional boundary element method (BEM) ones. The numerical results of BCM coincide very well with exact solution, and the feasibility and efficiency of the method are verified.     

瞬变电磁勘探中的人工神经网络反演法

为使瞬变电磁勘探的反演法简单和通用,提出了一种瞬变电磁勘探吵的人工神经网络反演法,基于人工神经网络的BP算法,设计和开发了一个自适应的人工神经网络反演系统,该系统避免了复杂的电磁场计算,只需经过学习训练就能够解决复杂的实际问题,而且具有记忆功能,从而使瞬变电磁法的反演工作具有延续性和可继承性,通过对实际的地质数据进行反演,结果表明,该方法是切实可行的。     

应用组合法计算三维非线性涡流场中的损耗

在样条积分方程法，有限元法和表面阻抗法的基础上，提出了一种求解三维非线性磁场中涡流损耗的新方法－－组合法，其特点是针对工程电磁场数值计算中几种方法各自的适用范围，将求解区域划分为积分方程求解区和有限元求解区，通过边界条件的处理将两种算法耦合，考虑到波在导体中传播的特性，应用表面阻抗的概念，把导体内部的涡流损耗作用反映到导体表面，并与有限元方程耦合求解出导体表面的涡流损耗，针对ＴＥＡＭ ＷＯＲＫＳＨ     

基于松耦合法的求解气动弹性问题的数值模拟方法

通过修改流体-结构界面处流体域的边界条件并引入松耦合方法，改进了一种求解气动弹性问题的数值计算方法．该数值计算方法的网格生成、流场求解、结构求解和数据交换四个过程相互独立，因而有可能利用现有的发展比较成熟的商用网格生成器、CFD（computational fluid dyllamics）求解器和CSD（computational structural dymrmcs）求解器分别完成这些功能．通过对折减阻尼比为0．0208的圆柱涡激振动问题在雷诺数为100条件下的计算结果进行分析，表明该方法是可行的．     

用随井追踪法确定油气藏边界

本文提出的随井追踪法是三极梯度法中的一种野外工作方法。大量的水槽模拟实验结果表明,一般的三极梯度法只适宜于确定油气藏短轴边界,油气藏长轴边界的确定还有赖于随井追踪法．等轴状油气藏边界的确定亦须采用随井追踪法。随井追踪法还可用于油气藏的垂向预测和开发监测。     

三维定常对流扩散方程的双方程方法

三维定常对流扩散方程的经典边界积分方程,其类型关于未知对流扩散势导数是第一类积分方程,关于未知对流扩散势是第二类积分方程。本文从格林公式出发,通过建立位势的单、双场守恒积分公式,推导出三维定常对流扩散方程新的边界积分方程,其类型与经典方程相反。对不同的边界采用不同的方程,由此把双方程边界元方法推广到三维空间。     

有限差分法计算多边形平面静电场

在直角坐标系中用有限差分法计算边界为五边形的平面静电场.提出了具有该边界条件的定解问题和相应算法,用fortran语言编写程序计算差分方程求出静电场的数值解,给出数值结果,将数值用origin画图软件得到五边形静电场分布情况.