天然差分方法在地下水渗流问题中的应用

在地下水渗流问题中,对不规则网格剖分后的结点控制元进行水量均衡分析,得到与伽辽金有限元法结果形式相同,但性质更优且自动满足质量守恒条件的稳定的计算格式.并可以简便的解决非线性的潜水问题.最后对泰斯问题进行了计算比较.     

一种改进的无单元伽辽金方法

使用单位分解积分,对传统的无单元伽辽金方法进行改进.有限覆盖和单位分解是单位分解积分的数学基础,对单位分解积分进行了严格证明,并指出使用Shepard函数作为单位分解函数是一个很好的选择.数值实例表明,使用单位分解积分进行数值求积的无单元伽辽金方法是一种真正的无网格方法,与经典的背景网格积分相比具有更高的精度.     

无网格法剖析及基扩充EFG法在电磁计算中的应用

主要论述基扩充的无网格法(MLM)用于2D电磁问题计算时的具体算法及编程问题。以独特的分步骤操作方法，介绍了无网格方法；从数值拟合的角度，对无网格伽辽金法(EFG)的核心技术——移动最小二乘法进行了深入剖析；严格按照加权余量法原理，利用偏微分方程的余量加权在节点支持域上的积分，导出了基扩充的EFG离散格式；应用FEM和基扩充的EFG两种方法对一些实例进行了计算验证。     

非线性热传导方程MLPG/RBF-FD无网格数值模拟

结合无网格局部彼得洛夫-伽辽金(MLPG)方法和径向基函数有限差分(RBF-FD)无网格方法求解非线性热传导问题。MLPG方法属于弱式无网格方法,具有处理边界条件方便的优点,然而因其要做大量的插值、积分运算而计算效率偏低;RBF-FD无网格方法属于强式无网格方法,直接对微分算子进行数值离散,计算效率高,然而其边界条件的处理较复杂。将二者相结合,在求解域边界附近采用MLPG方法,其它区域采用RBF-FD无网格方法,则能扬长避短。介绍了MLPG方法和RBF-FD无网格方法的基本原理,将该混合方法用来求解非线性热传导方程,数值算例显示了方法的正确性和高效性。     

数值模拟输运问题的一种破开算子方法

本文通过破开算子方法,把二维输运问题的控制方程破开为对流问题和扩散问题。在任意四边形网格的离散下,用特征线法解对流问题,并采用伽辽金加权余量法,从而有效地减少插值所引起的数值阻尼,提高计算精度。用有限单元法和迭代计算格式解扩散问题。由于采用了辛普生积分公式,在每个时间步长都不需要求逆矩律,节省了计算时间。算例表明,本文数值模拟结果与精确的理论解吻合较好。     

轴对称构件受力分析的插值粒子法

面对土木工程与机械工程中广泛存在的轴对称力学问题, 采用具有离散点插值特性的无网格方法形函数, 结合弹性力学空间轴对称问题的最小势能原理, 建立了轴对称构件力学分析的插值粒子法. 本文无网格法方法构造形函数不依赖网格, 也具有像有限元法一样可直接施加边界条件的优点. 本方法能直接获得全域连续应力场, 避免了有限元法应力后处理二次拟合带来的计算误差. 最后通过实例分析, 验证了所建立的无网格方法的有效性.     

基于伽辽金级数的空间电荷场计算模型

利用局部伽辽金矩量法,建立了电子注空间电荷场的2.5维计算模型. 将空间电荷场在径向展开为各个模式的叠加,利用伽辽金级数表示,通过求解其系数方程组, 最终求得空间电荷场.利用本文的模型,计算分析了电子注空间电荷场的特性, 并研究了模型中仿真参量对空间电荷场计算结果的影响.本文的空间电荷场计算模型利用粒子模拟基本方法, 可以应用到速调管注波互作用计算模型中.     

运动介质中电磁场计算的一种迎风有限元方法

本文提出了运动介质中正弦稳态电磁场问题的一种迎风有限元解法。用伽僚金法求解这类问题,当离散网格的Peclet数大于1时,计算结果会出现伪振荡。为了抑制这种振荡,引入了采用在迎风面与背风面具有不同迎风参数的权函数的迎风有限元法。该方法对一维问题,在均匀网格下能在节点上给出问题的精确解,在一维结果的基础上,提出了相应的二维解法,并用一个二维模型进行了验证。     

弹性力学的复变量无网格方法

在移动最小二乘法的基础上，提出了复变量移动最小二乘法.复变量移动最小二乘法的优点是采用一维基函数建立二维问题的逼近函数，所形成的无网格方法计算量小.然后，将复变量移动最小二乘法应用于弹性力学的无网格方法，提出了复变量无网格方法，推导了复变量无网格方法的公式.与传统的无网格方法相比，复变量无网格方法具有计算量小、精度高的优点.最后给出了数值算例. 关键词： 移动最小二乘法 复变量移动最小二乘法 无网格方法 弹性力学 复变量无网格方法     

三维电磁扩散场数值模拟及磁化效应的影响

采用矢量有限元法实现了三维电磁扩散场数值模拟,并成功将其应用在大地电磁的正演研究中.为灵活精确地拟合起伏地形和地下不规则构造,采用由不规则四面体单元组成的非结构化网格,可根据模型设计的需要调整网格的大小.引入了基于二次场理论,将解析的一次场从总场中扣除,直接计算二次场,使得误差仅局限于相对较小的二次场,以提高总场计算精度.常规的节点有限元法不满足电性分界面上法向电场不连续和无源区单元内电流密度无散,违反麦克斯韦方程组.为克服节点有限元法的弊端,使用矢量有限元法求解基于二次电场的偏微分方程.另外,在算法设计中,考虑了磁导率参数的变化,可以模拟磁导率不均匀的模型.通过与COMMEMI模型已发表的结果对比,证明了本文算法的正确性和精确性.为突显非结构网格优势,计算了椭球异常体模型和任意地形模型的MT响应,并详细讨论了地形和磁化效应对三维数值模拟结果的影响.     

相变温度场中无网格伽辽金法的应用

简述了无单元法的基础理论,推导出相变温度场的无单元法计算公式,采用罚函数法引入了第一类边界条件,编制了相应的计算程序.通过经典相变的应用例子,和有限元计算结果及解析解的比较,说明了无单元法应用于相变温度场具有连续性好,精度高,前后处理简单等优点.     

Study on Combined Method Based on 3-D ESPI

1　Ｉｎｔｒｏｄｕｃｔｉｏｎ　　Ｔｈｅｃｏｍｂｉｎａｔｉｏｎｏｆｅｘｐｅｒｉｍｅｎｔｗｉｔｈｃａｌｃｕｌａｔｉｏｎｆｏｒｍｓａｎｅｗｍｅｔｈｏｄ ,ｔｈａｔｉｓｃｏｍｂｉｎｅｄｍｅｔｈｏｄ (ＣＭ ) ,ｗｈｉｃｈｈａｓｂｅｃｏｍｅａｎｅｗｅｆｆｅｃｔｉｖｅｔｏｏｌｆｏｒｓｏｌｖｉｎｇ ｐｒｏｂｌｅｍｓｉｎｍｅｃｈａｎｉｃｓｉｎｒｅｃｅｎｔｙｅａｒｓ .ＴｈｅｃｏｍｂｉｎａｔｉｏｎｏｆＦＥＭｗｉｔｈｔｒａｄｉｔｉｏｎａｌｍｏｉｒ啨ｉｎｔｅｒｆｅｒｏｍｅｔｒｙｈａｓｂｅｅｎｒｅｐｏｒｔｅｄ[1] .Ｉｎｔｈｅｃｏ…     

计算流管声场的三维有限元数值模型

发展了一种三维有限元数值模型和计算方法来对矩形流管声场进行整体的计算.与以往的二维方法相比,此种数值方法不仅全面反映了矩形流管内声波的传播情况,而且提高了网格精度,从而大大扩展了对铺设有声衬的流管的计算领域.结果表明,该数值模型是有效和准确的,与其它方法和文献的计算结果吻合得非常好.同时,在大大增加计算量的同时,也对程序代码进行了优化工作,提高了计算效率.     

Study on Combined Method Based on 3-D ESPI

The finite element method (FEM),whether the calculation is accurate or not,depends closely on object boundary condition.If the three dimensional displacement of the object obtained in experiment is regarded as its boundary condition,a new method combining the results of experiment and calculation,called combined method (CM),is formed.The combined method possess advantages of experiment and calculation.It can correct calculation and improve the accuracy of FEM.Accordingly it has more practicability.In this paper,the three dimensional displacement fields of a typical beam loaded at three points are tested by using 3-D electric speckle pattern interferometry (ESPI).Using the experimental results as boundary condition the whole three-dimensional displacement fields can be calculated by FEM.The beam′s three-dimensional displacement fields obtained by FEM agree very well with those obtained by experiment.This proves that the combined method is effective and practicable.     

高分辨电子显微象实空间象模拟方法的改进

在深入研究高分辨实空间象模拟方法的基础上,本文导出了一个通用的收敛判据。这个判据确定了实空间象模拟计算中,在x-y平面上的取样间距δ和在电子束入射方向上片层厚度δ的依赖关系,从而有效地克服了实空间法中存在的发散问题。此外,还给出了精确计算传播因子的新公式,并在动力学衍射计算中获得了和传统的FFT多片层法符合很好的结果。 关键词：     

A finite-element method using dispersion reduced spline elements for room acoustics simulation

This paper presents a finite element method (FEM) using hexahedral 27-node spline acoustic elements (Spl27) with low numerical dispersion for room acoustics simulation in both the frequency and time domains, especially at higher frequencies. Dispersion error analysis in one dimension is performed to increase the accuracy of FEM using Spl27 by modifying the numerical integration points of element stiffness and mass matrices. The basic accuracy and efficiency of the FEM using the improved Spl27, which uses modified integration points, are presented through numerical experiments using benchmark problems in both the frequency and time domains, revealing that FEM using the improved Spl27 in both domains provides more accurate results than the conventional method does, and with fewer degrees of freedom. Moreover, the effectiveness of FEM using the improved Spl27 over that using hexahedral 27-node Lagrange elements is shown for time domain analysis of the sound field in a practical sized room.     

卷积完全匹配层在两维声波有限元计算中的应用

将基于复坐标变换和复频移扩展坐标变量的卷积完全匹配层Convolution Perfectly Machted Layer(CPML)引入到两维声波方程的有限元(FEM)计算中,该匹配层作为一种吸收边界条件Absorbing Boundary Condition(ABC)应用在有限元计算的边界截断上。文中分别给出了频域和时域的CPML方程的表达形式,并在有限元计算软件COMSOL中完成数值计算。相对于经典的PML,CPML最大的优势在于它不需要把场分裂开,这使其具有更好的稳定性和更高的吸收性能,且更易于实现。数值计算结果表明,CPML边界层有着比PML更好的吸收效果,其更有效的吸收了进入其中的声场能量。      

Application of convolution perfectly matched layer in finite element method calculation for 2D acoustic wave

A method was presented to extend the Convolution Perfectly Matched Layer(CPML), which bases on the complex coordinates transformation and complex frequency shifted stretched-coordinate metrics,to the 2D acoustic equation calculated with the method of Finite Element Method(FEM).This non-physical layer is used at the computational edge of a FEM as an Absorbing Boundary Condition(ABC) to truncate unbounded media.In this paper,the CPML equations have been presented in frequency domain and in time domain,respectively,and the calculations have been realized in the FEM software of COMSOL.The main advantage of CPML over the classical PML layer is that it is based on the unsplit components of the wave field leading to a more stable,highly effective absorption and a more facility to realize.The results of numerical simulation demonstrate that CPML has better absorbability than PML and it absorbs the outgoing energy more effectively.     

Comparison of 1D transfer matrix method and finite element method with tests for acoustic performance of multi-chamber perforated resonator

Multi-chamber perforated resonator (MCPR) is a kind of typical silencer element which can both attenuate broadband noise and satisfy specific installation requirements. The one-dimensional transfer matrix method (TMM) and finite element method (FEM) are widely used to predict the transmission loss of the resonators. This paper mainly focuses on the comparison between 1D TMM and FEM in which detailed perforation modeling is applied for the acoustic modeling of MCPRs. Five resonators with different acoustic attenuation frequency ranges are built for simulation and test. In order to verify the results of the above methods, a transmission loss test facility is designed based on two-load method. Through adjusting the distance between microphones, the facility’s effective measurement frequency can be changed. The results show that despite of the complex modeling and calculation, FEM with detailed perforation modeling shows good consistency with test results in both frequency and amplitude within entire frequency range. In comparison, TMM is limited by the cut-off frequency when calculating transmission losses. Besides, accuracy of TMM in low frequency range is also affected by perforation conditions. However, TMM is time-saving in calculation and structure optimization. In MCPRs’ development process, TMM can be used to quickly design and optimize structure parameters while FEM can be used to verify the acoustic performance before prototyping.