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

1、动力弹塑性分析的无网格自然单元法*_无网格法论文导读:：基于无网格自然单元法，提出了结构动力弹塑性响应分析的一条新途径自然单元法是一种新兴的无网格数值计算方法，其实质是基于自然邻近插值的伽辽金法自然单元法在本质边界条件的施加上较采用移动最小二乘法的无网格法具有明显的优势在空间域上采用自然单元法离散，并运用加权余量法推导了动力弹塑性分析的离散控制方程然后，采用预校正形式的Newmark法在时间域上进行求解最后给出了数值算例，并验证了所提方法的有效性和正确性论文关键词：无网格法，自然单元法，动力响应，弹塑性0引言在结构的动力响应过程中，通常总是既有弹性变形，又有塑性变形，而且这两种变形以及它们之间

2、的分界面都随时间变化1.正是由于弹塑性动力响应问题的复杂性，往往得不到结构弹塑性动力响应问题的解析解，需要借助数值方法进行求解. 在结构弹塑性动力响应的模拟方面，不论是有限元法2还是边界元法3, 4，都是基于网格的数值方法不合适的网格形状严重影响计算的精度，而且复杂结构的网格生成也是极具挑战性的问题. 相反，近年来发展迅速的无网格法具有近似函数不依赖于网格，且计算精度高，前后处理简单的优点5.目前发展的无网格方法主要有无单元Galerkin法6, 7、无网格局部Petrov-Galerkin法 8、边界无单元法9, 10、杂交边界点法11, 12、无网格流形方法13以及自然单元法14, 15等

3、. 与大多数无网格法不同，自然单元法的近似函数具有插值性且在边界结点间是线性变化的无网格法，从而可以方便地施加本质边界条件. 此外，自然单元法的形函数计算不仅具有不涉及矩阵求逆运算以及计算量较小的优点，而且也没有任何人为参数的选择问题. 因此，自然单元法是一种非常有发展前景的数值方法，并且已引起了国内外许多学者1622的极大关注.近年来，许多学者都致力于弹性动力学的无网格法研究. 程玉民等10将弹性动力学的边界积分方程方法与改进的移动最小二乘法结合，提出了弹性动力学的边界无单元法. 苗雨等12将双互易法同杂交边界点法相结合，提出了求解弹性动力学问题的双互易杂交边界点法. 李树忱等13运用无网格

4、流形方法求解了弹性动力学问题. 刘应华等23发展了弹性动力学分析的基于Voronoi结构的无网格局部Petrov-Galerkin法. 虽然朱合华等19，张英新等20和江涛等21在无网格自然单元法应用于弹塑性力学问题方面做了一些有益的工作，但目前尚未见到动力弹塑性分析中无网格自然单元法的研究成果.本文尝试将自然单元法应用于动力弹塑性分析问题的求解计算. 首先采用加权残值法详细推导了动力弹塑性分析的自然单元法理论公式，然后给出了其详细的数值实现过程. 最后，通过典型算例的计算和对比分析验证了方法的有效性和合理性1 控制方程考虑二维动力弹塑性分析问题，其计算域为，边界为，则平衡方程为：(1)式中，

5、为质量密度；为结点加速度，且有，为结点位移；为应力；为应变；为体力矢量；和分别为位移和面力已知的边界.应力和位移满足如下的边界条件(2a)(2b)和初始条件(3a)(3b)式中，为上点处的外法线的方向余弦；和分别为已知的位移和面力分量；和分别为初始位移和初始速度.弹塑性材料在进入塑性状态后，继续加载时的应力应变关系可以表示为(4)式中，为总应变增量. 采用关联的流动法则以及各向同性的加工强化模型，可得到弹塑性刚度张量为 2(5)其中(6a)(6b)(6c)这里，和分别为等效应力和等效塑性应变，为单向屈服应力，而和分别为剪切模量和泊松比.2 动力弹塑性分析的无网格自然单元法2.1 自然邻近插值1423自然邻近插值是一种多变量的插值方案，目前应用于无网格法的有Sibson插值和Laplace插值.本文采用Sibson插值.对区域内任一结点，其Voronoi结构定义为:(7)式中,为点与结点的距离.此外无网格法，二次Voronoi结构可定义为:(8)图1所示为平面7个结点的Vor