光滑粒子法及其在冲击动力学中的应用-工程力学专业论文.docxVIP

摘 摘 要 本文系统地分析了光滑粒子法的数学基础及离散过程，详细地讨论了与该方法 有关的一些基本问题，如离散格式、核函数、光滑长度、激波振荡处理、计算精度 和稳定性等问题，并在一些方面作了改进．最后将光滑粒子法应用于应力波传播和 高速碰撞计算，通过一些对比计算说明该方法的特点和长处． 首先利用插值理论，给出函数及函数各阶导数的核估计，分析了传统光滑粒子 法离散连续介质力学守恒方程组的基本方法及相应的离散形式．在保证系统总动量 和总能量守恒的基础上，提出了离散方程的理想离散形式．此外还讨论了光滑粒子 法中的一些基本问题，比如核函数的性质、光滑长度的选取、人工粘性、本构关系 以及该方法计算的主要步骤和时间步长选取等． 接着讨论了传统光滑粒子法中存在的两个缺点：边界点计算精度低和拉应力引 起的计算不稳定现象．提出通过设立虚拟点的方法，可以在一定程度上改善边界点 的核估计，但是却无法改善边界点导数的核估计．通过Fourier稳定性分析，给出 了出现拉引力不稳定的条件．建立在Taylor展开思想基础上的改进的光滑粒子法， 可以有效地克服传统光滑粒子法中存在的这些问题，因此是一种更为理想的离散格 式。通过一些简单分析指出，在某些限制条件下，光滑粒子法可以退化为有限差分 法． ， 一一一【粒子搜索算法在光滑粒子法中起着非常重要的作用，直接决定了计算效率，以 此为目标，利用核函数影响域的局域性，提出了一种搜索简便、计算效率获得很大 提高的分区搜索算法．边界的确定一直是光滑粒子法的一个难点．本文通过引入粒 子边界圆的概念，提出了依据边界圆的覆盖方式区分边界点和内点，有效地解决了 边界难以确定的问题． 将光滑粒子法用于应力波计算是冲击力学数值方法研究中的重要内容之一。一 维应变波的算例指出，采用通量修正输运法处理激波阵面有很好的计算效果．分析 指出，当光滑长度过小时，传统的光滑粒子法对光滑长度有很强的依赖性．但是光 滑长度取得过大，无论是传统的还是改进的光滑粒子法都可能发生压缩失稳．本文 在分析出现压缩失稳现象原因的基础上，提出了一种新的核函数，消除了压缩失稳 现象．一维和二维应力波算例表明，在数值模拟弹塑性波传播方面，光滑粒子法完 针对高速碰撞变形大、存在滑移面等基本特点，提出并建立了一套防止粒子与粒子 和粒子与界面发生相互“侵入”的算法，为高速碰撞数值模拟计算打下基础．Taylor 碰撞考核计算表明，光滑粒子法可以获得比有限元方法更好的结果和更逼真的变形 图象．高速侵彻和超高速碰撞数值箅例则更充分地展示了光滑粒子法计算大变形和 强滑移面问题的能力和优势．因此，在冲击力学数值计算中，光滑粒子法是一种有潜力的新算法． 强滑移面问题的能力和优势．因此，在冲击力学数值计算中，光滑粒子法是一种有 潜力的新算法． 关键词：光滑粒子法 核函数 搜索算法 应力波 冲击动力学／谴b{筻 ABSTRACTIn ABSTRACT In this thesis the mathematical background and discretization method of Smoothed Particle Hydrodynamics(SPH)are described and some related problems of SPH，such as discretization forms，kernel functions，smoothing length，shock oscillation control，com— putation accuracy，numerical stability and SO on，are also discussed．Meanwhile】some aspects of the method about SPH are improved and developed．Finally,the method is applied to the study of stress wave propagation and hypervelocity impact phenomenon and the results show that the method has good features in the numerical study of impact dynamics． Kernel estimate of a function and its derivatives D2e the key point in SPH method and it can be obtained by kernel interpolation theory．Based on the analysis of SPH method， the discretiza