穿甲燃烧弹侵彻陶瓷复合装甲和玻璃复合装甲的FEM-SPH耦合计算模型

为了提高小口径穿甲燃烧弹侵彻陶瓷复合装甲和玻璃复合装甲（透明装甲）的仿真分析精度,本文将传统的FEM（finite element method）-SPH（smooth particle hydrodynamics）耦合计算模型中穿甲燃烧弹弹芯的有限元模型和JC（Johnson-Cook）材料模型分别替换为SPH模型和JH2（Johnson-Holmquist-ceramics）材料模型,提出了新型FEM-SPH耦合计算模型。研究表明,新型FEM-SPH耦合计算模型可以有效模拟弹芯碎裂现象,减少SPH粒子和有限元耦合计算量,进而显著提高仿真模型的计算精度和计算效率,并给出了新型FEM-SPH耦合计算模型的有限元/粒子尺度和建模尺寸的优选结果。     

三维自适应FE-SPH耦合算法在多层间隔金属靶侵彻问题中的应用

鉴于有限元算法不能有效地模拟侵彻过程所产生的金属碎片, 本文中基于三维自适应FE-SPH耦合算法的基本理论, 自主开发了模拟多层间隔金属靶侵彻问题的三维FE-SPH耦合计算程序。该程序采用四面体单元对多层间隔金属靶侵彻模型进行初始离散, 计算过程中, 当四面体单元等效塑性应变超过某一设定值时, 单元自动转化为SPH粒子, 并引入有限单元-粒子接触算法和耦合算法, 实现大变形和破碎区域采用SPH方法计算, 克服有限元法单元畸变存在的问题。多层间隔靶侵彻算例分析表明, 三维FE-SPH耦合计算程序采用等效塑性应变作为转化判据计算结果较稳定, 并且能够有效地再现侵彻过程中所产生的碎片, 能够模拟侵彻碎片对后层靶的毁伤效应。     

近场动力学与有限元方法耦合求解热传导问题

求解含裂纹等不连续问题一直是计算力学的重点研究课题之一,以偏微分方程为基础的连续介质力学方法处理不连续问题时面临很大的困难. 近场动力学方法是一种基于积分方程的非局部理论,在处理不连续问题时有很大的优越性. 本文提出了求解含裂纹热传导问题的一种新的近场动力学与有限元法的耦合方法. 结合近场动力学方法处理不连续问题的优势以及有限元方法计算效率高的优势,将求解区域划分为两个区域,近场动力学区域和有限元区域. 包含裂纹的区域采用近场动力学方法建模,其他区域采用有限元方法建模. 本文提出的耦合方案实施简单方便,近场动力学区域与有限元区域之间不需要设置重叠区域. 耦合方法通过近场动力学粒子与其域内所有粒子(包括近场动力学粒子和有限元节点)以非局部方式连接,有限元节点与其周围的所有粒子以有限元方式相互作用. 将有限元热传导矩阵和近场动力学粒子相互作用矩阵写入同一整体热传导矩阵中,并采用Guyan缩聚法进一步减小计算量. 分别采用连续介质力学方法和近场动力学方法对一维以及二维温度场算例进行模拟,结果表明,本文的耦合方法具有较高的计算精度和计算效率. 该耦合方案可以进一步拓展到热力耦合条件下含裂纹材料和结构的裂纹扩展问题.      

一种新型SPH-FEM耦合算法及其在冲击动力学问题中的应用

为了充分发挥光滑粒子流体动力学方法(smoothed particle hydrodynamics,SPH)在处理大变形和有限元(finite element method,FEM)问题时计算精度高的优势,提出了一种新型SPH-FEM耦合算法.该耦合算法在大变形区域使用SPH粒子离散,其余区域使用有限元离散.在耦合界面...     

粘接-映射混合算法及其对颗粒材料中板的振动特性分析

结构与颗粒材料相互作用广泛存在于各工程领域,其研究过程中涉及的连续-离散耦合计算方法面对诸多挑战.本文提出了粘接-映射混合算法来研究连续体与离散介质耦合动力学问题.将连续体模型划分为内部区域及与颗粒接触的边界区域.边界区域采用粘接算法模拟连续体外部形状并使用高效的球形接触判断准则;提出一种包含Rayleigh阻尼映射的有限元映射质点弹簧算法来精确计算连续体内部区域内力和变形.二者相结合构成粘接-映射混合算法,并引入计算机集群和GPU(图形处理器)并行技术,对埋没于颗粒材料中受激振动固支方板的连续-离散耦合动力学问题进行了数值仿真研究.结果表明,粘接-映射混合算法有利于双层级并行算法的程序实现及优化,并在连续-离散耦合界面进行快速接触判断的同时实现对颗粒材料中方板位移、变形、振动形态等参数的研究.通过定幅扫频和定频变幅方式考察激振力频率和幅值对振动板非线性动力学行为的影响并观察到二倍周期现象,同时给出了该连续-离散耦合系统中颗粒体系的能量耗散特性.      

使用时间插值的特征线法在光滑粒子法拉氏边界条件中的应用

在计算流体动力学的实际应用中无反射边界条件是一个重要的研究课题. 文中应用时间插值的特征线方法构造了一种新式的拉氏边界条件,并应用在光滑粒子法中. 该方法与使用特征线方法的欧拉边界条件方法相比,不需要区分超声速和亚声速流的不同,并且在入流和出流中具有相同的形式,因而更加简便易行. 数值结果表明,采用时间插值特征线法的拉氏边界条件方法在稀疏波、激波以及爆轰波的模拟中都能够得到较好的无反射效果.      

SPH-FEM接触算法在冲击动力学数值计算中的应用

为了充分发挥光滑粒子流体动力学方法(Smoothed Particle Hydrodynamics,SPH)易于处理大变形以及有限元(Finite Element Method,FEM)计算精度和效率高的优势,论文基于无网格粒子接触算法,在有限元节点处设置背景粒子,通过接触力的方式计算SPH粒子和有限单元之间的接触问题...     

径向基点插值法在旋转柔性梁动力学中的应用

将无网格径向基点插值法用于旋转柔性梁的动力学分析. 利用无网格方法对柔性梁的变形场进行离散,考虑梁的纵向拉伸变形和横向弯曲变形,并计入横向弯曲变形引起的纵向缩短,即非线性耦合项,运用第二类拉格朗日方程推导得到系统刚柔耦合动力学方程. 将无网格径向基点插值法的仿真结果有限元法和假设模态法进行比较分析,说明假设模态法的局限性,并表明其作为一种柔性体离散方法在刚柔耦合多体系统动力学的研究中具有可推广性,并讨论了径向基形状参数的影响. 同时运用3 种求解系统动力学方程的方法:纽马克方法、4阶龙格库塔法、亚当姆斯预报校正法,并比较各方法的计算效率, 结果表明纽马克方法最快.      

DEM与FEM动态耦合算法研究

离散单元法作为一种有效的数值分析方法,能够模拟脆性材料的裂纹扩展及碎片飞散等破坏特性,但是无法从根本上克服计算效率低下的诟病;传统有限单元法具有计算高效稳定的优点,却难以描述脆性材料冲击破坏过程中连续体向非连续体的转化。本文首先提出一种基于罚函数法的改进型离散单元和有限单元耦合方法,以提高耦合分析精度。在此基础上提出了动态耦合算法:即在初始阶段,模型全部为有限单元,当局部即将发生破坏时,仅使即将发生破坏的有限单元及相邻单元自动转化为离散单元,在离散单元区域研究破坏问题。这种算法充分利用有限单元法计算高效的优点,同时最大限度克服了离散单元法计算效率的不足。最后,通过两个简单算例验证了改进型耦合算法和动态耦合算法的有效性。     

一种考虑界面不连续的改进的有限粒子法

有限粒子法（finite particle method，FPM）作为SPH（smoothed particle hydrodynamics）方法的重要改进，有效提高了边界区域粒子的近似精度，但是当FPM处理多物理场时，在不连续界面附近的计算精度会大大降低，并且FPM必须满足的矩阵非奇异性也提高了对界面处理的要求。本文中基于DSPH（discontinuous SPH）方法，提出了一种考虑界面不连续的改进FPM—DSFPM（discontinuous special FPM）法，旨在改善FPM在界面不连续处的计算精度，从而进一步提高其计算效率和稳定性。首先，分析了DSFPM的核近似精度。其次，根据不同的工程问题，给出DSFPM处理小变形和大变形问题的算法流程。利用DSFPM、DSPH和FPM等3种方法对弹性铝块小变形碰撞冲击算例进行了模拟，通过对比分析铝块的速度和应力以及计算时间验证了DSFPM算法在非连续界面处计算精度和计算效率的优势。最后，通过结合DSFPM和DFPM（discontinuous FPM）实现了对于大变形问题的模拟。     

光滑粒子法与有限元的耦合算法及其在冲击动力学中的应用

扼要讨论了光滑粒子法的离散思想,充分利用光滑粒子法和有限元方法各自的优点,提出了一种初始时刻用有限元建模,计算过程中大变形单元自动转换为光滑粒子的耦合算法。高速碰撞的系列算例说明,耦合算法不但适宜于计算大变形冲击动力学问题,而且由于集两种方法的优点于一身,可以更高效地模拟一些高速碰撞问题,提高计算效率。     

近场动力学最小二乘和有限元耦合方法研究

准确高效地对损伤和断裂问题进行建模是计算力学中的关键研究课题之一。将近场动力学最小二乘在处理含裂纹等非连续问题上的优势和有限元计算效率高及便于施加边界条件的优势结合，提出了近场动力学最小二乘和有限元耦合方法。将裂纹及其可能扩展区域划分为近场动力学区域，边界及其他区域划分为有限元区域，并将其中的结点类型分为近场动力学结点和有限元结点。有限元结点仅与同单元中的其他结点产生作用，近场动力学结点则与其族内的所有结点产生作用。将以上的单元刚度矩阵和质量矩阵进行组装得到整体刚度矩阵和整体质量矩阵。本文的耦合方法数值实现简单有效，相对于键基和常规态基近场动力学，该耦合方法包含了应力和应变的概念，同时不受零能模式的影响。一维和二维静态和动态问题的研究，验证了本文的耦合方法的有效性和准确性。     

固体介质中SPH方法的拉伸不稳定性问题研究进展

光滑粒子流体动力学法(smoothed particle hydrodynamics, SPH)是一种基于核估计的无网格Lagrange数值方法.它用粒子方程离散流体动力学的连续方程, 既可以处理有限元难于处理的大变形和严重扭曲问题, 又可以处理有限差分法不易处理的自由边界和材料界面的问题, 在固体力学中的冲击、爆炸和裂纹模拟中具有广阔的发展前景.但是, 该算法的拉伸不稳定性(tensile instability)问题是它在固体力学领域中应用的最大障碍.对SPH稳定性分析表明, 算法不稳定性的条件仅与应力状态和核函数的2阶导数有关.目前, 应力点法(stress points)、Lagrange核函数法、人工应力法(artificialstress)、修正光滑粒子法(corrective smoothed particle method, CSPM)和守恒光滑法(conservativesmoothing)以及其他一些方法成功地改善了SPH的拉伸不稳定性, 但是每一种方法都不能彻底解决SPH的拉伸不稳定性问题.本文介绍了SPH法的方程和Von Neumann稳定性分析的思想, 以及国内外在这几个方面的研究成果及其最新进展, 同时指出目前研究中存在的问题和研究的方向.      

Coupling between finite element method and material point method for problems with extreme deformation

As a Lagrangian meshless method, the material point method (MPM) is suitable for dynamic problems with extreme deformation, but its efficiency and accuracy are not as good as that of the finite element method (FEM) for small deformation problems. Therefore, an algorithm for the coupling of FEM and MPM is proposed to take advantages of both methods. Furthermore, a conversion scheme of elements to particles is developed. Hence, the material domain is firstly discretized by finite elements, and then the distorted elements are automatically converted into MPM particles to avoid element entanglement. The interaction between finite elements and MPM particles is implemented based on the background grid in MPM framework. Numerical results are in good agreement with experimental data and the efficiency of this method is higher than that of both FEM and MPM.     

关于弹体冲击和贯穿混凝土的三维数值模拟

利用FEM和SPH算法相结合给出弹体冲击和贯穿混凝土的三维数值模拟。将弹体划分为四面体有限元网格,混凝土划分为光滑粒子。为形成单元方便,将弹体划分成三个独立的模块并且这三个独立的模块具有拼凑功能。给出光滑粒子-有限元界面的滑移面算法。为了提高计算效率,在程序的前处理中预设界面的从属点和每个光滑粒子的关联点。计算结果与实验结果对比表明,FEM和SPH算法相结合及编制的程序能有效模拟弹体对混凝土的三维冲击和贯穿过程。同时给出在垂直冲击中二维计算与三维计算的差别。     

THE PARAMETRIC VARATIONAL PRINCIPLE AND NON-LINEAR FINITE ELEMENT METHOD FOR ANALYSIS OF ASTROMESH ANTENNA STRUCTURES

针对大型周边桁架式索网天线由拉索拉压模量不同引起的本构非线性和结构大变形引起的几何非线性问题,给出了基于参变量变分原理的几何非线性有限元方法. 首先针对含预应力索单元拉压模量不同分段描述的本构关系,通过引入参变量,导出了基于参变量及其互补方程的统一描述形式,避免了传统算法需要根据当前变形对索单元张紧/松弛状态的预测,提高了算法收敛性. 然后利用拉格朗日应变描述索网天线结构大变形问题,结合几何非线性有限元法,建立了基于参变量的非线性平衡方程和线性互补方程;并给出了牛顿-拉斐逊迭代法与莱姆算法相结合的求解算法. 数值算例验证了本文提出的算法比传统算法具有更稳定的收敛性和更高的求解精度,特别适合于大型索网天线结构的高精度变形分析和预测.     

SIZE OPTIMIZATION FOR CANTILEVER I-BEAM BASED ON VARIATIONAL METHOD

以变分方法为工具,研究了分布载荷下工字型悬臂梁的尺寸优化问题,给出了在强度和位移约束下梁缘条沿展向的宽度函数,并简要介绍了优化准则法和多岛遗传算法. 以平直机翼为工程实例,将之简化为一个工字型悬臂梁,分别采用该3 种方法进行优化设计. 结果表明变分法不仅效率高,而且给出了尺寸分布的函数,这对飞行器概念设计阶段的机翼结构设计具有较大的指导意义.     

基于径向基点插值法的旋转Mindlin板高次刚柔耦合动力学模型

将无网格径向基点插值法(radial point interpolation method,RPIM)用于中心刚体?旋转柔性板的动力学分析.基于浮动坐标系方法和一阶剪切变形理论即Mindlin板理论,考虑剪切变形的影响,并计入板面内变形的非线性耦合变形项,采用径向基点插值法描述板的变形场,保留动能中有关非线性耦合变形项...     

光滑粒子动力学方法的发展与应用

光滑粒子动力学(smoothed particle hydrodynamics,SPH)是一种拉格朗日型无网格粒子方法,已经成功地应用到了工程和科学的众多领域.SPH使用粒子离散及代表所模拟的介质,并且基于粒子体系估算和近似介质运动的控制方程.本文分析和综述了SPH模拟方法的发展历程、数值方法与应用进展.介绍了SPH方法的基本思想;从连续性、边界处理、稳定性和计算效率4个方面阐述了SPH方法的研究现状;介绍了SPH方法近年来在可压缩流动、不可压缩流动以及弹塑性材料高速变形与失效方面的一些典型应用;并对SPH方法的发展与应用进行了预测与展望.      

光滑节点插值法:计算固有频率下界值的新方法

将光滑节点插值法用于悬臂梁的静力学,并首次用于旋转柔性梁的频率分析. 采用梯度光滑技术,用线性插值形函数描述梁的位移场,求解4 阶微分方程. 在静力学分析中,将该方法所得梁中各点位移与假设模态法、有限元法及解析解的结果对比,可知该方法虽用简单的线性插值形函数描述梁的位移场,但精度却很高. 进一步研究表明,采用模态高于9 阶的假设模态法会使刚度阵条件数变差,导致结果发散. 在频率分析中,与有限元法、假设模态法和解析解对比,表明该方法一个重要特性:能提供固有频率的下界值,而有限元法和假设模态法只能提供固有频率的上界值,说明该方法结合有限元法在处理无解析解的问题时可以从上下界最大程度的逼近真实解,提高精度. 光滑节点插值法具有形函数结构简单、独立变量少且能提供固有频率下界值的特性,因此,具有较高的推广及应用价值.