一种自适应轴对称FEM-SPH 耦合算法及其在高速冲击模拟中的应用

为了准确、高效地模拟高速冲击问题,提出了一种自适应轴对称有限元(FEM)-光滑粒子流体动力学(SPH)耦合算法。该算法在初始时刻全部采用FEM 计算,在动态变形过程中自动将畸变单元转化为粒子,采用SPH 计算。该算法采用一种新的耦合算法实现单元与粒子间的高精度耦合,并应用最小内角转化准则和单元分组转化方式实现单元向粒子的自动转化。计算了几个典型的高速冲击问题:首先,通过计算应力波传播测试了新的单元-粒子耦合算法的精度;然后,通过计算泰勒杆问题验证了自适应耦合算法及相应程序的正确性;最后,计算了弹体侵彻铝板和混凝土板。结果表明:自适应耦合算法计算精度好且效率高,适合模拟高速冲击问题。     

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

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

高速碰撞数值计算中的SPH自适应粒子分布法

针对传统光滑粒子法在计算高速碰撞问题时会出现近邻粒子逸出核函数影响域而产生数值破坏这一缺陷,提出了一种根据粒子间距变化自动添加、合并粒子的SPH自适应粒子分布算法。采用该方法对Taylor碰撞和超高速碰撞问题进行了数值模拟,结果表明,该方法可以有效地消除计算中出现的数值破坏,提高计算精度。     

光滑粒子模拟方法在超高速碰撞现象中的应用

简要介绍了基于黎曼解的光滑粒子法,并将改进的SPH方法应用于超高速碰撞,对二维轴对称条件下的弹丸超高速碰撞薄板问题进行了数值模拟,研究了靶板厚度、弹丸速度、弹丸形状等因素对形成碎片云的影响。通过与实验数据比较,该算法模拟的碎片云的形状及特征与实验相吻合,验证了光滑粒子法对冲击动力学问题数值模拟的有效性。     

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

为了充分发挥光滑粒子流体动力学方法(Smoothed Particle Hydrodynamics,SPH)易于处理大变形以及有限元(Finite Element Method,FEM)计算精度和效率高的优势,论文基于无网格粒子接触算法,在有限元节点处设置背景粒子,通过接触力的方式计算SPH粒子和有限单元之间的接触问题.使用SPH-FEM接触算法分别对球头钢弹斜冲击钢板和平头钢弹正冲击钢板的情况进行了三维数值计算,SPH采用完全变光滑长度算法,FEM采用矩阵向量积的EBE(Element-By-Element)算法.将SPH-FEM接触算法的计算结果与LS-DYNA的计算结果以及实验结果进行了对比验证.     

光滑粒子法中的一种新的核函数

分析了传统的核函数产生压缩失稳现象的原因,提出了消除这种压缩失稳现象的一种新的核函数。采用改进的光滑粒子法,对几种常用的核函数进行了一维应变波的对比计算。结果表明:所提出新的核函数在应力波计算中不但保证了计算精度,还能有效地消除压缩失稳。     

弹性力学的重构核粒子边界无单元法与有限元的耦合法

将重构核粒子边界无单元法(RKP-BEFM)与有限元法(FEM)耦合,形成求解具有区域特征的弹性力学问题的重构核粒子边界无单元与有限元的耦合方法RKP-BEF/FE.推导了重构核粒子边界无单元与有限元耦合方法的离散化公式,建立了节点未知量的耦合方程.重构核粒子边界无单元法和有限单元法的较高精度保证了这一直接耦合方法的成功实现与求解精度.最后给出了平面问题的数值算例,验证了提出的耦合方法RKP-BEF/FE的有效性.     

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

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

弹塑性波计算中的光滑粒子法

采用改进的光滑粒子法 ,对脉冲应力载荷下板中弹塑性一维应变波的传播进行了数值计算 ,比较了人工粘性法和通量修正法处理冲击波间断面的效果。结果表明 ,改进的光滑粒子法在应力波数值计算中有良好的精度。     

基于光滑粒子动力学的煤岩刨刀刨削过程数值模拟

为研究刨削参数对刨刀载荷特性的影响规律,对刨刀破煤过程进行仿真,为得到更可靠的仿真数据,采用光滑粒子动力学(SPH)与有限元进行耦合(FEM)对刨刀破煤过程进行模拟,并将仿真结果与传统有限元算法进行比较证明此算法的优越性。通过模拟不同刨削参数下的刨刀的破煤过程,得到刨刀载荷数据。研究结果表明：随着刨削速度的增加,刨刀在刨削过程中的整体受载变化不大,但刨刀所受载荷峰值和波动范围随着随刨削速度的增加而增大;随着刨削深度的增加,刨刀在刨削过程中的整体受载有着明显的变化,具体表现为：刨刀所受载荷峰值和均值随着随刨削深度的增加而增大,但刨削深度对刨刀载荷波动范围影响不明显。     

基于拟流体模型的SPH新方法及其在弹丸超高速碰撞薄板中的应用

引入颗粒动力学理论(拟流体模型)建立了适用于超高速碰撞的SPH新方法。将超高速碰撞中处于损伤状态的碎片等效为拟流体,在描述其运动过程中引入了碎片间相互作用和气体相对碎片的作用。采用该方法对球形弹丸超高速碰撞薄板形成碎片云的过程进行了数值模拟,得到了弹坑直径、外泡碎片云和内核碎片云的形状、分布,并与使用传统SPH方法、自适应光滑粒子流体动力学(ASPH)方法的模拟结果进行对比,结果显示：新方法在内核碎片云形状和分布上计算结果更加准确。同时对Whipple屏超高速碰撞问题进行了研究,分析了不同撞击速度下防护屏弹坑尺寸及舱壁损伤特性等特性,计算结果与实验吻合较好且符合Whipple防护结构的典型撞击极限曲线。

     

基于SPH方法的表面张力修正算法及其应用

针对传统SPEI方法中基于CSF模型的表面张力算法,在计算边界、尖角等粒子缺失部位的曲率时存在偏差较大,且粒子秩序较差,对大变形问题表面张力计算精度较低的问题,在Morris提出的表面张力SPH方法基础上,通过引入CSPM方法对边界法向的计算和曲率的计算进行修正,得到了表面张力修正方程组.应用本文方法模拟了水溶液中初始...     

近场动力学与有限元的混合建模方法

综述了近场动力学与有限元混合建模方法的研究进展,阐明了各种混合建模方法的基本原理与特点,并重点介绍本课题组在近场动力学与有限元方法混合建模方面的研究工作。现有近场动力学与有限元混合建模方法包括位移协调约束、力耦合、混合函数方法以及子模型方法等,除子模型方法外,都可归结为并行式多尺度分析方法,其基本思想是将计算结构划分为近场动力学子域、有限元子域以及两者的交界区域(或重叠区域、或界面单元、或过渡区域)。子模型方法可归结为显-显分析方法,先采用显式有限元进行整体分析,后采用近场动力学方法对重点区域进行分析。混合建模方法需要着重提高交界区域的计算精度,并且消除虚假力和虚假应力波问题。提出了通过力耦合的近场动力学与有限元混合建模的隐式分析方法,该方法不再设置重叠区,通过杆单元连接近场动力学子域与有限元子域,其中界面上的有限元结点不仅与其所在单元的其他结点发生作用,还通过杆单元与以其为圆心、一定半径的圆域内的其他物质点相互作用。研究表明,本文提出的混合模型和求解方法既能有效解决裂纹扩展等不连续问题,又可提高计算效率,为工程结构破坏问题的计算分析提供一种有效方法。     

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

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

随机有限元方法在断裂分析中的应用

在幂律非线性随机有限元基础上,以单边裂纹板为例给出计算含量钢继裂参数,J（J积分）,δ（裂纹张开位移）,Δ（由裂纹引起的裂纹板上下底面相对位移）,θ（由裂纹引起的裂纹板上下底在相对转角）及其对基本随机变量变化率的方法和分析算例。     

超高速撞击数值仿真结果分析

为了分析超高速撞击过程的宏观现象和内在机理,对9.53 mm铝球以6.64 km/s的速度撞击2.2 mm厚的铝靶的SPH仿真结果进行了量化分析。结果表明：采用SPH方法以及Steinberg弹塑性模型和Mie-Grneisen状态方程,可获得与试验相符的仿真结果;球形破片开坑或穿孔直径遵循初始快速增加、然后缓慢增加,直至稳定的变化规律;破片/靶板界面的最大撞击压力比材料强度大两个量级以上;靶板阻抗力最大值发生在破片最大直径侵入靶板时刻;碎片云的运动过程具有自相似演化特征,其运动范围不会超出碎片云的包络圆锥范围。