动力弹塑性分析的无网格自然单元法

基于无网格自然单元法,提出了结构动力弹塑性响应分析的一条新途径.自然单元法是一种新兴的无网格数值计算方法,其实质是基于自然邻近插值的伽辽金法.自然单元法在本质边界条件的施加上较采用移动最小二乘法的无网格法具有明显的优势.在空间域上采用自然单元法离散,并运用加权余量法推导了动力弹塑性分析的离散控制方程.然后,采用预校正形式的Newmark法在时间域上进行求解.最后给出了数值算例,并验证了所提方法的有效性和正确性.     

三维弹塑性自然单元法算法实现

自然单元法是一种新兴的无网格数值计算方法,其实质是基于自然相邻插值(C∞)的伽辽金法。该方法计算精度与四边形或六面体单元有限元法相当,自然相邻插值函数比其他无网格法插值函数的计算速度快。由于自然相邻插值在凸域的边界上的相邻点之间是严格线性的,所以自然单元法在边界面的处理也相当简单。本文研究了在自然单元法中采用Von.Mises,Mohr-Coulomb和Drucker-Prager屈服准则解决三维弹塑性问题,并编制了相应计算程序,最后通过算例验证算法的正确性。     

弹塑性分析的线性互补问题

为了实现弹塑性分析的非迭代类算法，本文从一条简捷而严密的途径导出了弹塑性边值问题的线性互补方程，与其它非迭代类算法相比，这里所建立的线性互补问题的规模要小得多。文中讨论了解的适定性的问题，给出了求解工程问题的实用算法和若干算例。     

自然单元法的自适应研究

自然单元法(NEM)是一种新兴的无网格数值计算方法,具有前处理简单和易于准确施加本质边界条件等优点.本文基于Z-Z后验误差估计方法,给出了一种自然单元法的误差估计因子和自适应分析细化方案.采用结点处的光滑应变计算相应的恢复应力,并用于构造全域上的恢复应力场.通过结点Voronoi单胞内的能量范数相对误差指示需要进行结点...     

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

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

弹性力学的杂交自然单元法

为了解决自然单元法无法直接求解节点应力以及应力解精度不高的问题, 将应力杂交的思想引入自然单元法中, 与弹性问题的Hellinger-Reissner变分原理结合, 提出了弹性问题的杂交自然单元法. 算例表明: 杂交自然单元法的计算结果与解析解吻合, 证明该方法是可行的; 在求解应力方面, 杂交自然单元法比自然单元法有着更高的计算精度, 而且可以直接求解出节点的应力.     

直接增强自然单元法计算应力强度因子

自然单元法是一种新兴的无网格数值计算方法,但应用于裂纹问题计算时,其近似函数并不能准确反映裂纹尖端渐进应力场的奇异性,为获得足够的计算精度,需要在缝尖附近增大结点的布置密度。针对裂纹问题提出一种增强的自然单元法,将缝尖渐近位移场函数嵌入到自然单元法近似函数中,给出了增强试函数的构造方法,推导了总体刚度矩阵和荷载列阵的相关列式。应力强度因子可以作为附加未知量直接算得,也可用J积分或相互作用能量积分方法进行计算,对增强区域的选择和影响进行了分析。算例结果表明,基于增强自然单元法采用围线积分方法计算应力强度因子具有很高的精度,但直接以附加结点自由度形式计算则精度有所降低。     

弹塑性力学问题的虚单元法

具有有限差分法特征的虚单元法,可视为是有限元法向任意多边形单元的扩展.对于颗粒增强复合材料等非均质材料,在材料细观力学性能表征、非均质材料力学分析等非线性问题方面,传统的弹塑性有限元具有网格数目多、效率低下等不足之处,而虚单元法使网格划分更加灵活,更能细观反应材料真实结构,为材料的弹塑性力学分析等非线性问题提供了新的思路.基于增量法弹塑性力学原理和双线性投影算子,建立了弹塑性力学问题的虚单元法求解技术,提出了弹塑性力学问题虚单元法的应力更新方案,研究了弹性力学问题虚单元法的精度和收敛性,讨论了虚单元法求解弹塑性力学问题的网格依赖性.同时,开展了任意多边形和凹多边形单元的数值试验研究,结果表明,虚单元法无须分割多边形,仅需节点自由度便可求得单元刚度矩阵和应力等效荷载,程序实现简单,计算精度高,改善了传统有限元的网格依赖性和塑性区的网格奇异性.     

自然单元法研究进展

自然单元法是一种基于Voronoi图和Delaunay三角化几何结构,以自然邻点插值为试函数的一种新型数值方法.其既具有无网格方法和经典有限元方法的优点,又克服了两者的一些缺陷,是一种发展前景广阔的求解微分方程的数值方法.自然单元法的形函数满足插值性质,可以像有限元法一样直接施加本质边界条件,不存在基于移动最小二乘拟合的无网格方法不能直接施加本质边界条件的难题.由于自然单元法是无网格方法,可以方便处理有限元方法较难处理的一些问题,例如移动边界和大变形等问题.自然单元法与其他数值方法的最根本区别于其插值格式的不同.将自然邻点插值用于Galerkin过程,就得到基于Voronoi结构的自然单元Galerkin法.自然邻点插值有自然邻点Sibson插值和Laplace插值(非Sibson插值)两种.Laplace插值比Sibson插值在计算上要简单的多,并且不论对凸的或非凸的区域都能精确施加本质边界条件.以Laplace插值为试函数的自然单元法在数值实施上比以Sibson插值为试函数的自然单元法简单.本文对基于Voronoi结构的自然邻点插值和自然单元法的基本思想作了介绍,综述了国内外关于自然单元法的研究成果,总结了自然单元法的优点和尚需解决的问题.      

基于自然单元法的极限上限分析

自然单元法是一种基于离散点集的Voronoi图和Delaunay三角化几何信息,以自然邻近插值为试函数的新型数值方法.相对于一般无网格法中常采用的移动最小二乘近似而言,自然邻近插值不涉及到复杂的矩阵求逆运算,更不需要任何人为的参数,可以提高计算效率.采用该方法构造的形函数满足Delta函数的性质,可以像有限元一样准确地施加边界条件,可以方便处理场函数及其导数的不连续性的问题.论文将自然单元法应用到极限上限分析中,编制了相应的计算程序,通过极限分析的几个经典算例进行了验证,同时采用类似于分片应力磨平的方式,编制相应的磨平程序,由计算点上的塑性耗散功外推得到了节点上的塑性耗散功的值,从而画出了极限状态下结构的塑性耗散功的分布云图.计算结果表明采用自然单元法求解极限上限分析具有稳定性好,精度高,收敛快等优点.     

自然单元法计算裂纹与材料边界问题

提出一种新的非凸边界上自然单元法形函数计算方法,通过边界结点限制点对间的邻点关系,对包括裂纹和材料边界在内的各种类型的非凸边界具有统一的处理原则,所得到的近似函数在边界结点间具有线性插值性.     

用数字散斑相关法研究聚氨酯泡沫塑料的压缩力学性能

提出了一种改进的数字散班相关计算方法，使其能直接进行应变迭代，灵敏度可达１０００με，并对它进行了验证试验。对于两种不同密度的聚氨酯泡沫塑料，采用改进的数字散斑相关法对其压缩力学性能进行了测试，得到了其应力——应变关系及弹性模量等力学参数。     

On the Strain Hardening Stability of Initially Isotropic Metallic Materials under Natural Ageing

The strain hardening stability of initially isotropic structurally stable metallic materials under natural ageing is studied experimentally. Martensitic steel and titanium and aluminum–magnesium alloys, which harden in accordance with the isotropic–kinematic hypothesis, are used as an example to show that the shape of the loading surface does not change during ageing and its size and position may vary with time     

应变硬化薄板夹持型滚弯成形过程分析

现存文献对薄板的连续滚弯成形过程的数学模拟多限于讨论三辊轮解锥型弯板过程,且基于较多的假定和板的小变形情况。而开发有更多的功能和更高的生产效率的四辊轮连续弯板机需研究夹持型滚弯成形过程。本文将「１」的工作,理想塑性材料平面应变薄板的夹持型连续滚弯成形过程的弹塑性大变形弯曲的力学分析,推广至一般应变硬化材料的情形。得出了控制微分方程的解析解。采用固支边界条件模拟刚度很大的平持系统对板的夹持。通过对线     

High-order Asymptotic Analysis for the V-notch in a Power-law Hardening Material

The aim of this study is to determine the complete elastic-plastic stress asymp-totic solutions at the plane V-notch tip with the two edges being clamped-clampe...     

加锚体自然单元法分析与程序设计

建立了加锚体的自然单元法分析模型。结合自然相邻插值和一维轴力杆单元的劲度矩阵,推导了锚杆支护数值模拟的单元列式,给出了自然单元法分析加锚体的实施步骤和主要程序设计。结构的离散和锚杆的离散完全独立,因此可以自由地加设或拆除锚杆。算例分析表明自然单元法分析加锚体的精度与有限元的相当,但比有限元的实施要方便得多。     

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

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

线性蠕变问题的摄动边界元-有限元方法

本文将摄动、边界元、有限元方法结合起来,提出一种求解线性蠕变问题的新方法。该方法不采用一般增量法中在一个时段内各物理量保持不变或作线性变化的假设,加大了计算步长提高了精度。文中构造了边界元摄动格式,构造了包含钢筋在内的边界元有限元耦合摄动格式,并给出了满意的数值结果。     

线性船舶兴波阻力边界元新方法

采用常数边界元对船舶与流体界面进行离散，求解船舶兴波势及船舶兴波阻力。这种方法可避免在船舶与流体自由面交线上安置节点，因而避免了这些节点建立补充方程。因为满足自由面条件的Ｈａｖｅｌｏｃｋ源函数的源点和场点不能同时在自由面上，使得自由面上的节点无法用Ｈａｖｅｌｏｃｋ源函数的建立方程。如对自由面交线上的节点建立补充方程，则要对线性自由面条件中包含的未知势函数的二阶导数用差分形式表达，引入较大误差。     

Finite Element Analysis of Elasto-plastic Plate Bending Problems using Transition Rectangular Plate Elements

In this work, the finite element analysis of the elasto-plastic plate bending problems is carried out using transition rectangular plate elements. The shape functions of the transition plate elements are derived based on a practical rule. The transition plate elements are all quadrilateral and can be used to obtain efficient finite element models using minimum number of elements. The mesh convergence rates of the models including the transition elements are compared with the regular element models. To verify the developed elements, simple tests are demonstrated and various elasto-plastic problems are solved. Their results are compared with ANSYS results.The English text was polished by Keren Wang.