无网格配置法与两类问题有限元：理论、应用与比较研究.docxVIP

无网格配置法与两类问题有限元：理论、应用与比较研究

一、引言

1.1研究背景与意义

1.1.1数值计算方法发展概述

数值计算方法的发展源远流长，其历史可以追溯到古代的数值逼近方法，如差分法和插值法。随着计算机科学与技术的迅猛发展，数值计算方法迎来了前所未有的变革与发展，逐渐涵盖了插值和逼近、微积分、常微分方程、偏微分方程、线性代数等众多领域，成为科学计算和工程技术中不可或缺的关键工具。

在众多数值计算方法中，有限元法（FiniteElementMethod，FEM）占据着举足轻重的地位。其起源可追溯到20世纪40年代，AlexanderHrennikoff在1941年提出将连续介质离散为格子结构以求解弹性问题的思想，这成为有限元法的雏形。随后，RichardCourant在1943年发表的论文中，利用Rayleigh–Ritz方法和分片多项式函数在三角形区域内逼近问题解，为有限元法奠定了重要的数学理论基础。到了1960年，美国加州大学伯克利分校的RayW.Clough教授首次正式提出“有限元方法”这一术语，并展示了其在飞机结构分析中的应用，标志着有限元法作为一种通用数值分析工具的正式诞生。此后，有限元法在理论和应用方面都取得了长足的发展。IvoBabu?ka和FrancoBrezzi提出的Babu?ka–Brezzi条件（又称LBB条件）为混合有限元方法提供了稳定性和收敛性的充分条件，Sobolev空间理论的引入则进一步完善了有限元法的误差估计和收敛性分析，使其在数学上建立了更加严格的理论基础。

然而，随着科学技术的不断进步，有限元法在处理一些复杂问题时逐渐暴露出局限性，如在分析大变形、不连续性、奇异性或裂纹动态扩展等问题时，有限元法面临着网格畸变、网格重构等挑战。例如，在处理材料加工成型过程中的大变形问题时，传统有限元法划分的网格会发生严重畸变，导致计算精度下降甚至计算失败。在求解裂纹动态扩展问题时，需要不断重新划分网格以适应裂纹的扩展，这不仅增加了计算的复杂性，还降低了计算效率。

为了克服有限元法的这些缺陷，无网格配置法（MeshlessMethod）应运而生。无网格配置法最早可追溯到20世纪70年代提出的光滑粒子流体动力学（SmoothedParticleHydrodynamics，SPH）方法，该方法最初用于天体物理领域的研究，随后在其他领域得到了应用和发展。无网格配置法的核心特点是不需要进行网格的生成和调整，而是直接在计算域内的离散点集上进行计算。这种方法能够有效避免有限元法中网格相关的问题，在处理复杂几何形状、动态变化边界以及大变形等问题时具有显著优势。经过多年的发展，无网格配置法已衍生出多种具体方法，如有限点法（FinitePointMethod，FPM）、再生核粒子法（ReproducingKernelParticleMethod，RKPM）等，在工程和科学计算领域展现出了巨大的潜力。

1.1.2无网格配置法和有限元法的重要性

有限元法凭借其对复杂几何构形的良好适应性、在各种物理问题中的广泛可应用性、建立在严格理论基础上的可靠性以及适合计算机实现的高效性，在工程学、物理学、生物学等众多领域得到了极为广泛的应用。在土木工程领域，有限元法被广泛用于建筑结构、桥梁、地下工程等的力学分析和设计。通过建立结构的有限元模型，可以准确预测结构在各种荷载作用下的应力、应变和变形情况，为结构的优化设计提供重要依据。在航空航天工程中，有限元法用于飞机、航天器等的结构强度分析、气动弹性分析等，对于保障飞行器的安全性能和飞行性能起着关键作用。在生物医学工程领域，有限元法可用于模拟人体骨骼、软组织等的力学行为，为医学诊断、手术规划和医疗器械研发提供支持。

无网格配置法由于其独特的无需网格化特性，在处理有限元法难以解决的问题时发挥着重要作用。在处理具有复杂几何形状的流场模拟问题时，无网格配置法能够更加灵活地适应复杂的边界条件，使得流体流动的模拟更加精细和真实。在材料科学领域，无网格配置法在研究材料的大变形、裂纹扩展等问题时具有显著优势，能够为材料的性能优化和失效分析提供更准确的结果。在生物组织模拟方面，无网格配置法可以更好地模拟生物组织的复杂形态和力学特性，为生物力学研究提供了新的手段。

1.1.3研究意义

从理论完善的角度来看，深入研究无网格配置法和两类问题有限元，有助于进一步拓展和完善数值计算方法的理论体系。无网格配置法中总体刚度矩阵的构建方法尚处于不断探索和发展阶段，研究其构建方法可以丰富无网格方法的理论内涵。对两类问题有限元在刚体动力学和弹性动力学中的深入研究，能够深化对有限元法在不同力学领域应用的理解，