求解二维弹性问题的径向基差分法

常规的配点型无网格法在求解弹性力学问题中,存在求解精度差和纽曼边界条件处理等局限．为解决这一问题,通过利用流体力学中基于径向基构造的差分格式(RBF-FD),来求解弹性力学平面问题．同时,为了进一步提高求解精度,对纽曼边界条件采用Hermite插值进行处理．数值算例表明,该方法具备良好的收敛性,并有着较高的精度,可有效解决传统配点型无网格法精度差的问题．同时,也表明该方法可以应用于弹性力学问题的求解．     

基于局部Petrov-Galerkin离散方案的无网格法

基于局部Petrov-Galerkin离散方案,选用自然邻近插值构造试函数,用Shepard函数作为权函数,提出了一种无网格方法(MNNPG),这种方法充分发挥了局部Petrov-Galerkin法的优势,并且结合了自然邻近插值的特点,方便引入边界条件,由于以Shepard函数的圆形支集作为积分子域,用分片中点插值来完成区域积分,无需额外背景网格,是一种真正的无网格法。本文将该无网格方法用于求解二维弹性力学边值问题,算例结果很好地吻合了精确解,表明该方法具有良好的数值精度和稳定性。     

配点型广义节点无网格法的基本原理及其应用

借鉴流形方法思想,引入广义节点的概念,对传统的无网格法进行了改进,建立了可具有任意高阶多项式插值函数的广义节点无网格方法。同时采用径向插值函数构造具有插值特性的逼近函数;采用配点法建立系统的离散方程。在阐述了这种方法基本原理的同时,针对线弹性力学问题给出了这种方法的数值计算列式。与传统无网格方法相比,这种方法更具有一般性;同时由于采用了配点法而不需要背景积分网格,所以可以认为这种方法是某种真正意义上的无网格法。当选取0阶广义节点位移插值函数时便可得到传统的无网格法;在不增加支持域内节点数目的条件下,通过选取高阶广义节点位移插值函数可以提高计算精度。最后通过算例分析,对0阶、1阶及2阶广义节点无网格法与现有的有关解答进行了对比,论证了其合理性。     

轴对称弹性体扭转问题的无网格自然邻接点Petrov-Galerkin法

本文将无网格自然邻接点Petrov-Galerkin 法应用于轴对称弹性体扭转问题的求解．无网格自然邻接点Petrov-Galerkin 法采用自然邻接点插值构造试函数，并且采用三角形线性单元的形函数作为加权残值法的加权函数．自然邻接点插值构造的试函数满足Kronecker delta 函数性质，因此本质边界条件的施加十分方便．由于几何形状和边界条件的轴对称特点，原来的空间问题简化为二维问题求解，因此计算时只需要横截面上离散节点的信息．数值算例结果表明，所提出的方法对求解轴对称弹性体扭转问题是行之有效的．     

基于局部移动Kriging无网格方法的层合板自由振动分析

利用基于局部移动Kriging插值无网格法对层合板自由振动进行了数值分析,基于一阶剪切层合理论导出了层合板振动的控制方程和边界条件,进一步得到了自由振动的离散化特征方程。由于Kriging插值函数具有Kronecker delta函数性质,可以直接施加本质边界条件。通过本文给出的方法,对不同边界条件、不同跨厚比、不同材料参数和铺设角度的层合板的振动频率进行了计算,均得到满意结果。最后用该方法对层合板的铺设角度进行优化设计,得到了与已有文献完全一致的优化结果。数值结果充分表明了无网格Kriging方法分析层合板自由振动问题的有效性和高精确度。     

弹性静力问题的无网格弱-强形式结合法

基于改进的移动最小二乘(MLS)二阶导数近似,建立了一种求解弹性静力问题的无网格弱-强形式结合法(MLS-MWS)。该方法采用节点离散求解域,通过MLS构造形函数,将求解域划分为边界域和内部域,并分别使用控制方程的局部弱形式和强形式来建立离散系统方程。对强形式中涉及的近似函数二阶导数计算,提出了一种将其转化为求两次一阶导数的方法,与传统方法相比,该方法计算简单、精度高。MLS-MWS法结合了弱、强形式无网格法的优点,Neumann边界条件容易满足,并且只需在边界区域进行积分。文中应用该方法分析了两个弹性力学平面问题,分析结果表明本文方法具有良好的精度和收敛性。     

轴对称结构极限下限分析的自然单元法

作为一种基于自然邻近插值的新型无网格法,自然单元法克服了大多数无网格法难以施加本质边界条件的困难．将自然单元法与减缩基技术相结合,建立了一种轴对称结构极限下限分析的数值格式和求解算法．通过不断修正自平衡应力场,轴对称结构极限下限分析可转化为一系列的非线性数学规划子问题,并由复合形法求解．在每个非线性规划子问题中,自平衡应力场表示为一组带有待定系数的自平衡应力场基矢量的线性组合,并且这些自平衡应力场基矢量可由弹塑性增量分析的平衡迭代结果得到．算例结果表明,本文所提的轴对称结构极限下限分析方法行之有效．     

伽辽金型无网格法的数值积分方法

无网格法直接通过节点信息构造形函数,不依赖于节点之间的有序单元连接,能够建立任意高阶连续的整体协调形函数。与传统的有限元法相比,无网格法对大变形问题、移动边界问题和高阶问题的求解方面有比较明显的优势。伽辽金型无网格法是目前应用最为广泛的一类无网格法。虽然无网格形函数本身不依赖于单元,但伽辽金型无网格法需要采取合适的方法进行弱形式的数值积分。由于无网格形函数一般不是多项式,具有非插值性且影响域与背景积分网格通常不重合,伽辽金型无网格法通常需要采用高阶的高斯积分进行数值积分,导致了计算效率低下,难于求解大型实际问题。因此,如何通过建立高效积分方法提高无网格法的计算效率成为无网格法研究领域的一个核心问题。本文对无网格法中强制边界条件的施加方法进行了简要归纳,总结了伽辽金型无网格法中若干常用的数值积分方法,并对伽辽金型无网格法的数值积分方法领域存在的一些问题进行了探讨。     

考虑尺寸效应的双材料轴对称界面端应力奇异性

提出了一种分析双材料轴对称界面端的应力奇异行为的特征值法.基于弹性力学空间轴对称问题的基本方程和一阶近似假设,利用分离变量形式的位移函数和无网格算法,导出了关于应力奇异性指数的离散形式的奇异性特征方程.由奇异性特征方程的特征值和特征向量,即可确定应力奇异性指数、位移角函数和应力角函数.数值求解了纤维/基体轴对称界面端模型的奇异性特征方程, 结果表明:尺寸效应参数δ(奇异点与轴对称轴的距离和应力奇异性支配区域大小的比值)影响着应力奇异性的强弱与阶次, 准一阶近似解析解只是δ>>1时的一个特例.      

无网格稳定配点法及其在弹性力学中的应用

伽辽金型无网格法具有精度高、稳定性好的优点,但是实现高阶准确积分过程复杂,计算效率低.配点型无网格法的计算效率高,但是其在求解复杂问题时往往会出现精度和稳定性较差的结果.本文介绍一种新的无网格法-无网格稳定配点法,采用重构核近似作为近似函数,在规则子域内非常容易实现高阶准确积分,既保留了配点型无网格法效率高的特点,又具备伽辽金型无网格法精度高和稳定性好的特点,而且还兼具有限体积法满足局域离散方程守恒的特点.通过弹性力学算例验证了该算法的优越性,未来可将其进一步应用于流体和流固耦合问题分析.     

基于自然单元法的轴对称结构极限上限分析

自然单元法是一种以自然邻近插值为试函数的新兴无网格数值方法,其形函数的计算不涉及矩阵求逆,也不需要任何人为参数。为了充分发挥自然单元法的优势,本文基于极限分析上限定理建立了轴对称结构极限上限分析的整套求解算法。轴对称结构的位移场由自然邻近插值构造,并且采用罚函数法处理材料的不可压条件。为了消除目标函数非光滑所引起的数值困难,采用逐步识别刚性区和塑性区,并对两者用不同方法进行处理。数值算例结果表明,本文提出的轴对称结构极限上限分析方法是行之有效的。     

用无网格LRPIM分析中厚板弯曲问题

用径向基函数构造无网格点插值法的形函数,插值函数具有Kronecker delta函数性质,因此可以很方便地施加本质边界条件.利用无网格局部径向点插值方法分别对一个对边固支另对边简支中厚板和一个悬臂中厚板的弯曲进行了分析计算.该方法不需要任何形式的网格划分,所有的积分都在规则形状的子域及其边界上进行,是一种真正的无网格方法.算例表明:将无网格局部径向点插值法应用于计算中厚板的弯曲问题,所求得的位移场和应力场都是光滑的;在径向基函数的基础上,附加多项式大大提高了插值精度;所得结果与弹性力学理论解以及有限元解都十分吻合.     

弹性动力学的双互易杂交边界点法

将双互易法同杂交边界点法相结合,提出了求解弹性动力问题的新型数值方法------双互易杂交边界点方法. 该算法在求解弹性动力问题时,将控制方程非齐次项的域内积分转化为边界积分. 该方法将问题的解分为通解和特解两部分,通解使用杂交边界点法求得,特解则使用局部径向基函数插值得到,从而实现了使用静力问题的基本解来求解动力问题. 计算时仅仅需要边界上离散点的信息,无论积分还是插值都不需要网格,域内节点仅用来插值非齐次项,因此该算法仍是一种边界类型的无网格方法. 数值算例表明,该方法后处理简单,计算精度高,适合于求解弹性动力问题.      

基于分区径向基函数配点法的大变形分析

无网格法因为不需要划分网格, 可以避免网格畸变问题,使得其广泛应用于大变形和一些复杂问题. 径向基函数配点法是一种典型的强形式无网格法,这种方法具有完全不需要任何网格、求解过程简单、精度高、收敛性好以及易于扩展到高维空间等优点,但是由于其采用全域的形函数, 在求解高梯度问题时 存在精度较低和无法很好地反应局部特性的缺点. 针对这个问题,本文引入分区径向基函数配点法来求解局部存在高梯度的大变形问题. 基于完全拉格朗日格式,采用牛顿迭代法建立了分区径向基函数配点法在大变形分析中的增量求解模式.这种方法将求解域根据其几何特点划分成若干个子域, 在子域内构建径向基函数插值, 在界面上施加所有的界面连续条件,构建分块稀疏矩阵统一求解. 该方法仍然保持超收敛性, 且将原来的满阵转化成了稀疏矩阵, 降低了存储空间,提高了计算效率. 相比较于传统的径向基函数配点法和有限元法, 这种方法能够更好地反应局部特性和求解高梯度问题.数值分析表明该方法能够有效求解局部存在高梯度的大变形问题.      

基于分区径向基函数配点法的大变形分析

无网格法因为不需要划分网格,可以避免网格畸变问题,使得其广泛应用于大变形和一些复杂问题.径向基函数配点法是一种典型的强形式无网格法,这种方法具有完全不需要任何网格、求解过程简单、精度高、收敛性好以及易于扩展到高维空间等优点,但是由于其采用全域的形函数,在求解高梯度问题时存在精度较低和无法很好地反应局部特性的缺点.针对这个问题,本文引入分区径向基函数配点法来求解局部存在高梯度的大变形问题.基于完全拉格朗日格式,采用牛顿迭代法建立了分区径向基函数配点法在大变形分析中的增量求解模式.这种方法将求解域根据其几何特点划分成若干个子域,在子域内构建径向基函数插值,在界面上施加所有的界面连续条件,构建分块稀疏矩阵统一求解.该方法仍然保持超收敛性,且将原来的满阵转化成了稀疏矩阵,降低了存储空间,提高了计算效率.相比较于传统的径向基函数配点法和有限元法,这种方法能够更好地反应局部特性和求解高梯度问题.数值分析表明该方法能够有效求解局部存在高梯度的大变形问题.     

求解弹性力学问题的有理单元法

传统的位移有限元法采用多项式形式的位移试函数,对于边数大于4的多边形单元,构造满足单元间协调性要求的多项式形式位移插值函数是一件困难的工作。本文利用逆距离权插值的思想并考虑到单元节点的分布,建立了边数大于4多边形单元上的有理函数形式的形函数。利用有理试函数,采用Galerkin法推导出求解平面弹性力学问题的有理单元法。采用有理单元法求解弹性力学问题,求解区域根据需要可以划分为任意多边形单元,极大地提高了网格划分的灵活性。有理单元法不依赖等参变换,不同单元的形函数表达形式统一,方便计算程序的编写。     

平面弹性问题的位移-应力混合重心插值配点法

提出数值分析平面弹性问题的位移-应力混合重心插值配点法。将弹性力学控制方程表达为位移和应力的耦合偏微分方程组,采用重心插值近似未知量,利用重心插值微分矩阵得到平面问题控制方程的矩阵形式离散表达式。使用重心插值离散位移和应力边界条件,采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,应用最小二乘法求解过约束方程组,得到平面弹性问题位移和应力数值解。数值算例结果表明,重心Lagrange插值方法的计算精度可达到10~(-10)量级。位移-应力混合重心插值配点法的计算公式简单、程序实施方便,是一种高精度的无网格数值分析方法。     

无网格法及其最新进展

无网格法具有许多独特的优点,因此有人认为无网格法将成为继有限元法之后新一代的数值方法. 到目前为止,已提出了不下十几种无网格法, 这些无网格法各有不同的优缺点.本文评述几种主要无网格法(着重应用于固体力学中的基于弱式的无网格法),将它们适当地分类, 论述典型的无网格形状函数的构造方法,介绍无网格法的发展现状, 评价已提出无网格法的优缺点,并比较典型无网格法的相同和不同点, 以及讨论无网格法发展所面临的问题等.最后就无网格法的发展趋势进行了展望.      

奇异边界法模拟二维弹性力学问题

奇异边界法是一种新的边界型无网格数值离散方法。该方法使用基本解作为插值基函数,在继承传统边界型方法优点的同时,不需要费时费力的网格划分和奇异积分,数学简单,编程容易,是一个真正的无网格方法。为避免配置点与插值源点重合时带来的基本解源点奇异性,该方法提出了源点强度因子的概念,从而将边界型强格式方法的核心归结为求解源点强度因子。本文首次将该方法应用于求解平面弹性力学问题。数值算例表明,本文算法稳定,效率高,并可达到很高的计算精度。     

考虑体积力时应力函数的求法

求解弹性力学平面问题时,常常引用应力函数解法.但在一般弹性力学教材中,对如何寻求应力函数论述不多.文献[1]、[2]采用了用边界上应力函数φ及其导数的力学意义来确定应力函数.这是一种行之有效的方法,它适用范围较大,力学意义明确,可为寻求应力函数指明方向.但[1]、[2]中有关公式只适于无体力的情况,这无疑地使其应用受到一定限制.本...