北京国际力学中心无网格、质点类和扩展有限元法国际学术研讨会将在南京召开

近年来，无网格、质点类和扩展有限元法的发展突飞猛进，国内外学者在广义流体动力学、非连续的光滑粒子流体动力学方法、爆炸模拟、水下爆炸冲击模拟、与分子动力学耦合的多尺度模拟等应用领域作出了大量贡献．为了进一步推动相关科研和学术活动的开展，北京国际力学中心拟于2009年10月12～17日在南京举办“北京国际力学中心无网格、质点类和扩展有限元法国际学术研讨会”．     

一种改进的无网格局部Petrov-Galerkin方法

本文提出了一种改进的无网格局部Petrov-Galerkin方法来分析平面弹性力学问题．这种无网格方法采用移动最小二乘近似函数（MLS）来近似试函数,采用Heaviside函数作为加权残值法中的权函数,采用直接插值法来施加本质边界条件．最后通过数值实例表明平面弹性力学问题中改进的无网格Petrov-Galerkin方法具有收敛快、稳定性好、精度高和简单有效的特点．     

科技人才招聘

中国科学院计算技术研究所1．大规模数值模拟算法与软件博士后2～3人研究方向：大规模数值模拟算法与软件。要求：计算力学、计算数学等相关专业背景。对有限元，无网格等通用数值方法有深入了解和自己独到见解；有并行程序开发经验，懂MPI，C／C等；扎实的数学力学基础。     

有限元网格自动生成的组合处理方法

给出有限元网格自动生成的相互配合协调处理的新方案，它处理的是平面上任意形状的多连通区域．边界描述采用边界顶点序列法，生成有限元分析的单元以四边形元（八点元或四点元）为主；根据边界形状的复杂性辅以三角形元（六点元或三点元）．所建议的组合处理方法具有容易控制网格密度、确保有限元分析中的单元质量、对于指定点（如集中荷载作用点或支撑点）可以自动成为网格的结点、可实现无人工控制的有限元网格自动生成等特点，对于有限元重分析系统尤其表现出高效率．本方法可推广到相当广泛一类三维问题中，在并行处理方面也有很大的潜力．     

自适应多重网格有限元求解应力集中问题

基于Ｚｉｅｎｋｉｅｗｉｃｚ—Ｚｈｕ误差估计方法、自适应策略以及分子表结构，本文实现了自适应多重网格有限元程序．计算结果表明，自适应多重网格有限元求解应力集中问题具有精度高、速度快的特性．     

空间有限元网格自动划分的方法

本文介绍了利用计算机编程，对空间有限元计算所需网格进行划分的一种方法，并给出了有关计算公式，从而为空间有限元计算提供了网格、节点和单元的信息．作者使用该方法对内燃机某些零件进行了有限元网格划分，结果表明该方法简单可行．     

三维弹性力学问题中有限元方程的预处理方法

针对三维弹性问题中有限元方程的数值求解，建立了一类简单且实用的代数多重网格预处理共轭梯度法(AMG-CG法)，详细描述了相应代数多重网格方法的粗化技术及网格转移算子的构造．由于该预处理方法能有效地降低刚度矩阵的条件数，使刚度矩阵的谱分布更集中，从而大大提高了计算效率．数值结果表明，AMGCG法对求解三维弹性问题有限元方程是十分有效和健壮的。     

斜辊管材矫直机机架的强度有限元分析

利用有限元方法建立了三雏的有限元力学模型,施加在矫直机机架的前面计算的工作载荷,合理配置了边界条件,采用映射方法划分了实体有限元网格,创新性的完成了矫直机机架的力学模型的实体有限元化.     

线弹性力学问题局部间断有限元方法的后验误差估计

基于线弹性力学问题的局部间断有限元方法,对于无悬点型的三角剖分网格, 给出了应力张量场L2-范数下的一个残差型全局后验误差估计,并给出该误差估计式的可靠性和有效性分析.     

无网格法与有限元法的耦合在线弹性断裂力学中的应用

文章采用无网格Galerkin方法与有限元(EFG-FE)耦合的方法来计算裂纹问题,这种耦合的方法不仅解决了无网格Galerkin法力学边界条件施加的难点,而且还克服了无网格Galerkin法耗时较多的缺点;运用线弹性断裂力学理论,分别采用了线性基函数、二次基函数和部分扩展基,对有限板单边裂纹的应力强度因子和受拉中心斜裂纹方形板进行了分析,数值计算结果表明了方法的有效性.     

Meshless Least-Squares Method for Solving the Steady-State Heat Conduction Equation

The meshless weighted least-squares (MWLS) method is a pure meshless method that combines the moving least-squares approximation scheme and least-square discretization. Previous studies of the MWLS method for elastostatics and wave propagation problems have shown that the MWLS method possesses several advantages, such as high accuracy, high convergence rate, good stability, and high computational efficiency. In this paper, the MWLS method is extended to heat conduction problems. The MWLS computational parameters are chosen based on a thorough numerical study of 1-dimensional problems.Several 2-dimensional examples show that the MWLS method is much faster than the element free Galerkin method (EFGM), while the accuracy of the MWLS method is close to, or even better than the EFGM These numerical results demonstrate that the MWLS method has good potential for numerical analyses of heat transfer problems.     

无网格伽辽金法(EFGM)求解接触问题

本文考虑到无网格伽辽金法只需节点信息,不需将节点连成单元,并且有精度高,后处理方便等优点,从而用它求解接触问题。这里采用的方法是将Ｋａｔｏｎａ界面单元引入ＥＦＧＭ,迭代求得两物体间的接触状态。算例表明,本文方法基本可行。     

采用无网格伽辽金法求解电容层析成像正问题

由于有限元法求解电容层析成像正问题的计算准备及后处理非常费时，对正问题的三维求解造成了瓶颈，为此，提出采用无网格伽辽金法求解电容层析成像正问题，获得正问题的弱变分形式，并用拉格朗日乘子法施加边界条件，从而得到数值解．在同样的仿真条件下，2种方法的计算时间分别为14．046S和5．078S．对5种典型流型进行仿真，结果表明，2种方法计算结果的最大相对误差为2．25％．因此，无网格伽辽金法与有限元法具有相当的精度，且计算速度有较大提高．     

无网格伽辽金法在二维结构问题中的应用研究

无网格伽辽金法(EFGM)是一种新型的求解偏微分方程的数值计算方法,不需对结构进行有限元网格的离散化,只需节点信息而不需将节点连成单元.本文论述和研究了EF-GM的基本原理与实现过程,主要包括用移动最小二乘法(MLS)构造形函数、用变分原理推导控制方程、用拉格朗日乘子法增强本征边界条件和域的高斯积分4个主要过程.基于MAT-LAB平台,实现了二维弹性结构的EFGM算法,并将典型算例的EFGM求解结果与有限元近似解、解析解结果进行了比较,结果表明了EFGM算法的正确性和有效性.     

一种自适应影响域半径无网格Galerkin法

针对某些力学问题的数值求解需要结点的局部加密,在采用背景积分网格积分方式的基础上,提出一种影响域半径随结点疏密程度而变化的自适应影响域半径无网格Galerkin法。在方法中,无网格结点与背景积分网格的结点重合,结点的影响域半径即可根据该结点周围的网格的最大边长来选取。算例显示,该文方法是可行而有效的。     

基于Taylor展开的无单元插值形函数及应用

在无单元伽辽金法的基础上,构造了基于Taylor展开的具有过点插值的无单元形函数,它可以和有限元法一样处理边界条件,克服了传统的无单元伽辽金法遇到的瓶颈问题;对非凸边界的处理,提出了新的准则--弧弦准则(arc-string criterion).这样,可大大减少了无单元法的计算工作量,提高了边界处理的精度,并且继承了无单元法及有限元法的优点.     

关于无单元法中的插值基函数选取的探讨

无单元法不需要单元信息，它采用了一种基于移动最小二乘（MLS）的插值函数。插值基函数对插值函数以及无单元法的计算精度影响很大。本文就不同的基函数对插值函数及无单元法的计算精度的影响作了分析比较，得出了一些有益的结论，并用算例说明了这些结论的正确性。     

模拟平板轧制的新方法——无网格伽辽金法

将一种新的数值方法无网格伽辽金法(EFGM)用于刚塑性可压缩材料稳态轧制过程的模拟,由于形函数不满足插值条件,采用罚函数法满足本质边界条件;为提高精度,选用矩形影响域的张量积核函数;利用有限元背景网格作为积分单元,对求解域内和边界上采用不同的高斯积分方案·数值计算结果与刚塑性有限元的计算结果和文献中的实验数据吻合较好,说明无网格伽辽金法用于刚塑性可压缩材料轧制过程的可行性和正确性·     

无网格伽辽金法在裂纹扩展计算中的应用

采用了一种基于t-分布的新型权函数，提高了无网格伽辽金法的计算精度；采用完全变换法处理本质边界条件，实现了本质边界条件在节点处的精确施加；针对裂纹扩展中的实际情况，对动态影响半径法作了进一步的补充和改进．算例验证了方法的正确性和有效性．     

分区粘弹性介质样条边界元反分析数值解法

通过改造非均质样条边界元列式方法,建立了分区粘弹性介质边界元反分析的数值求解格式,该数值解法与均匀弹性体边界元法相类似,不需作分区组合计算,简化了计算程序,减少了计算工作量、提高了计算精度和解题效率,为解决非均质围岩反分析中遇到的困难奠定了基础。