基于偶应力理论的局部彼得洛夫—伽辽金无网格法：原理、构建与工程应用.docxVIP

基于偶应力理论的局部彼得洛夫—伽辽金无网格法：原理、构建与工程应用

一、引言

1.1研究背景与意义

在现代材料科学与工程领域，随着微纳制造技术的飞速发展，材料和结构的尺寸不断减小，进入微米甚至纳米尺度范围。在这些微观尺度下，大量实验研究发现材料的力学性能呈现出明显的尺寸效应，即材料的力学响应不再仅仅取决于宏观的应力应变关系，还与结构的特征尺寸密切相关。例如，在细铜丝的扭转实验中，当铜丝直径减小到微米量级时，无量纲的扭转硬化显著增加；在薄梁弯曲实验中，随着梁厚度的减小，无量纲的弯曲硬化也明显增强。这些现象无法用传统的经典连续介质力学理论进行解释，因为经典理论基于连续性假设，忽略了材料微结构的影响以及力偶的作用，认为物体是由无数没有体积的空间点组成，只考虑了平移特征承受力矢，而未考虑转动与改变转动特征的力偶。

为了解释材料在微米尺度下的尺寸效应现象，应变梯度理论应运而生。应变梯度理论通过将高阶应变梯度纳入支配材料行为的本构或演化方程，来引入尺度对结构或系统的弹、塑性变形等力学行为的影响。偶应力理论作为应变梯度理论的一种，在解释材料尺寸效应方面具有重要作用。它考虑了材料微结构中力偶的作用，引入了独立的微观转动矢量，使得系统能量不仅包括应力对应变做的功，还包含偶应力对旋转形变做的功。通过引入与材料微结构相关的特征长度参数，偶应力理论能够有效地描述材料在微观尺度下的力学行为，弥补了经典连续介质力学的不足。

在数值计算方法领域，传统的有限元法作为一种广泛应用的数值求解方法，在处理复杂工程问题时取得了巨大的成功。然而，有限元法依赖于网格的划分，在处理大变形、移动边界、裂纹扩展等问题时，会面临网格畸变、重新划分网格等困难，这不仅增加了计算的复杂性，还可能导致计算精度的下降。为了克服有限元法的这些局限性，无网格法应运而生。无网格法摒弃了传统的网格划分，采用一系列离散的节点来近似求解问题的域，从而解除了节点之间的网格约束，在处理复杂几何形状、大变形和移动边界等问题时具有天然的优势。

局部彼得洛夫—伽辽金无网格（MLPG）法作为一种“真正”的无网格法，在求解过程中不需要背景网格。它将控制方程的子域积分“弱”形式引入方法中，降低了对试探函数连续性的要求，使得试探函数与测试函数可以取自不同的函数空间，子域构造形式多样，具有灵活多样的实施方案。同时，MLPG法刚度矩阵的形成不需要装配过程，计算效率较高。将偶应力理论与MLPG法有机结合，构建基于偶应力理论的MLPG法，能够充分发挥两者的优势，为解决具有尺寸效应的复杂工程问题提供一种强有力的数值计算方法。这种方法不仅能够准确地描述材料在微观尺度下的力学行为，还能有效地克服传统数值方法在处理复杂问题时的局限性，对于推动微纳尺度下材料与结构的力学分析和工程应用具有重要的理论意义和实际应用价值。

1.2国内外研究现状

偶应力理论的发展可以追溯到一个多世纪前，Voigt最早提出了体力偶和面力偶的概念，并建议构建考虑作用在材料微粒界面上力偶的连续模型。随后，Cosserat兄弟根据相关假设建立了Cossera理论，对应的运动方程中出现了偶应力。但在当时，该理论并未得到广泛关注和深入研究。直到20世纪60年代左右，Toupin和Mindlin等学者开始尝试对Cossera理论进行改进和扩展工作，他们对Cossera连续体物质点的旋转施加一定约束，逐渐发展了一种更为普遍的偶应力理论。在传统偶应力弹性理论中，材料颗粒不仅在力的作用下作平动，还在力偶的作用下作转动，系统能量包括应力功和偶应力对旋转形变做的功，对于各向同性线弹性材料，本构方程中除了两个经典的拉梅系数外，还包含两个与材料微结构有关的附加常数。然而，这种理论在构建中仅用到传统的力和力矩平衡关系，对偶应力张量的球形部分和体力偶在应力张量本构关系中的表现形式无法确定，存在一定的局限性。

为了克服传统偶应力理论的不足，Yang等学者从引入高阶平衡关系角度出发，提出了一种修正偶应力理论。在该理论中，添加力偶矩平衡关系后，偶应力张量被约束成对称量，它对与之共轭张量的曲率张量的对称部分做功，并与应力对应变做的功一起转变为系统能量。这种理论下的本构方程仅包含一个附加常数，大大降低了非经典常数的确定难度，使得理论在实际应用中更加简便和可行。Park和Gao使用这种新理论计算弯曲微梁，发现当微梁厚度与材料内察长度相当时，呈现出明显的尺度效应，所求的无量纲弯曲刚度与弯曲实验测量值吻合得较好，进一步验证了修正偶应力理论在描述材料尺寸效应方面的有效性。

无网格法最早可以追溯到Lucy等解决天体物理问题中用到的光滑质点流体动力学（SPH）方法。自那以后，无网格法得到了迅速发展，根据