切口零件冲击破坏的SPH模拟

简要介绍了光滑粒子动力学方法的基本思想和基本方程，利用该方法编制程序对带长方形切口式样的冲击破坏过程进行了数值模拟。给出了用这种方法模拟三维冲击破坏的一个算例，并对有关问题进行相关的讨论。     

SPH方法在溃坝数值模拟中的应用

光滑粒子流体动力学(SPH)方法是一种纯拉格朗日无网格方法,适用于模拟波浪破碎、高速水流、高速冲击碰撞等瞬时极大变形的问题。介绍了SPH方法的基本原理、核函数及离散格式的控制方程、边界处理方法等,建立了数值水槽模型,验证了模型的有效性。应用SPH方法模拟了二维溃坝问题,考虑了添加无孔障碍物和有孔障碍物的2种情形,将模拟得到的结果进行对比分析。结果表明,SPH方法能够很好地捕捉流体自由面变形、飞溅及融合现象,在一定程度上有孔障碍物有更好的消能效果。     

自由表面流动问题的SPH方法数值模拟

利用光滑粒子流体动力学（SPH）方法对两个自由表面流动问题－水坝坍塌及涌波进入静止水塘作了数值模拟,尝试了一种处理固壁边界条件的方法－－边界力方法。数值模拟结果表明,SPH方法在处理流体复杂大变形问题方面具有优越笥,而边界力方法则是SPH中处理固壁边界的一种行之有效的方法。     

压铸充型过程的SPH方法建模及数值模拟

基于欧拉描述的数值方法在求解压铸充型过程自由表面问题时存在一定困难.为提高自由表面数值模拟的准确性,建立了基于纯拉格朗日描述的压铸充型过程光滑粒子流体动力学(SPH)方法计算程序.引入Monaghan边界模型建立了型腔壁面边界条件,通过划分入流区域粒子和流体粒子满足了入流边界条件.采用所编制的程序计算了Schmid验证...     

基于SPH方法的鸟撞飞机风挡的数值模拟

风挡抗鸟撞是飞机安全飞行的重要保证.文中基于飞机圆弧风挡受鸟体撞击的实验观察,建立了国产某型军用飞机圆弧风挡及鸟体的有限元分析模型,利用光滑粒子流体动力学法(SPH)耦合有限元法对圆弧风挡受鸟撞击的过程进行了数值模拟,计算得到风挡结构的变形、位移、应变、撞击力、应力、临界撞速、发生破坏的可能位置及其破坏方式等几方面的数据,并考察了SPH粒子疏密对计算结果的影响.研究表明,数值模拟结果与实验结果基本吻合;鸟撞整个过程约4ms,撞击中点、前1/3处和后1/3处,风挡发生破坏(包括安全破坏)的临界撞速分别约为(540±5)、(600±5)和(470±5)km/h;鸟撞过程中,风挡的位移与其厚度是同一量级,风挡的最大应变已达到10~(-2)量级;风挡首先发生破坏的位置在后弧框附近,然后向与风挡中线成45°角的方向发展;SPH粒子数越多,鸟体变形模态越好.     

基于SPH方法高速射流冲刷的模拟

针对有限单元法在处理大变形中的网格重构、网格畸变等问题的困难,采用SPH方法模拟高速射流冲刷过程。SPH方法是一种纯拉格朗日型的无网格粒子法,能够有效地避免大变形中的各种网格问题。采用该种方法主要模拟不同工况下的冲坑发展情况,冲刷对象选LS-DYNA中的MAT147材料模型。计算结果表明采用SPH方法可以形象、直观地反映高速射流冲刷的过程,并且能够清楚地描述粒子速度场和砂床内部压力场的演变,所得结果可以为高速淹没射流的研究提供参考。     

基于SPH方法的不同材质射流毁伤性能研究

为研究不同材质射流的毁伤性能,使用AUTODYN有限元软件,采用光滑粒子流体动力学（SPH）方法对Cu、PTFE、PTFE-Cu三种材料药型罩形成射流的成型及侵彻靶板过程进行了数值仿真,并通过实验进行验证.研究结果表明:Cu材料药型罩在爆轰波的作用下形成凝聚的射流,而PTFE和PTFE-Cu材料药型罩则形成飞散的粒子流;三种材料射流侵彻靶板过程中,Cu射流头部速度最低,侵彻深度最深,开孔最小;PTFE粒子流头部速度最高,侵彻深度最浅,开孔大小居中;PTFE-Cu射流的头部速度和侵彻深度都居中,而开孔最大;PTFE-Cu射流克服了PTFE射流侵彻性能不足的缺点,其开孔能力较之铜射流有所提高.      

SPH方法后处理的新途径

重点研究并挖掘对化合物分子三维立体结构可视化显示软件RasMol的功能。在用实心球体粒子集合显示连续体,并提出了一套比较完整的基于Rasmol软件的光滑粒子流体动力学（SmoothedParticleHydrodynamics,简称SPH）后处理的新途径。为sPH和其它无网格方法的后处理问题提供可行的方法和新途径。     

基于SPH方法与GPU并行计算的阶梯流数值模拟

为了更好地模拟分析地下空间阶梯流动力学规律和提升数值计算效率,本文利用光滑粒子流体动力学方法(SPH)并借助GPU并行计算技术对地下空间阶梯流进行数值建模和计算。基于GPU并行化加速的SPH地下空间阶梯流模型的数值计算结果表明：该模型不仅计算获得较好的阶梯流流态,而且计算效率大幅度提高,比串行SPH模型和移动粒子半隐方法(MPS)模型快数百倍,这为研究大尺度城市街区洪水提供了一个比较有前景的研究途径。     

面向光滑粒子流体动力学(SPH)流体的湍流细节模拟

研究基于光滑粒子流体动力学（smoothed particle hydrodynamics,SPH）方法的恢复涡流和湍流细节的方法,缓解因采用人工黏度而产生的数值耗散.根据由黏性力引起的动能损失率对离散点的涡度进行修正.通过流函数将涡度变化恢复到速度场,避免求解Biot-Savart积分这一耗时过程.该方法基于旋度计算,不仅可以增强现有的涡流,还可在潜在位置产生新的湍流,并且其易于集成到现有SPH方法中.实验证明,该方法可以逼真地模拟湍流细节,并实现不同强度水平的增强效果.      

SPH法数值仿真三维切削破岩和切削力估算

以Mohr-Coulomb模型描述岩石,在LS-DYNA中用SPH(Smootheed Particle Hydrodynamics)法模拟三维刀齿切削破岩,探讨切削参数和岩石性质对切削力的影响.基于切削力随各切削参数和岩石性质的变化规律,建立切削力估算公式.用所得公式计算切削力,并与实验值进行比较.仿真结果表明:切削厚度、切削前角、切削宽度、岩石的内聚力和内摩擦角对切削力影响显著,切削速度在较低范围内变化时对切削力影响非常小.三维切削仿真能够较好地反映刀齿两侧岩石颗粒对切削的影响,比二维理论模型更接近实际.所得切削力估算公式能准确计算仿真中的平均切向力,且其计算值能与实验吻合较好.验证了基于SPH法的数值仿真用于研究刀齿切削破岩的可靠性.     

棱形液舱内液体晃荡问题的SPH数值模拟

应用改进的SPH方法模拟棱形液舱在不同充液比、横摇激励频率工况下的液体晃荡行为。该方法能模拟出液体大幅度非线性晃荡产生的翻卷、破碎等现象,模拟结果与实验结果吻合较好。通过改变舱体结构,对带有隔板的棱形液舱进行数值模拟,分析隔板对液体晃荡特性的影响。结果表明:充液比及激励频率对液体的晃荡特性有重要影响,且适当改变箱体结构可以有效地抑制晃荡现象。     

液固高速撞击时材料表面损伤的数值模拟

为了研究液固撞击的机理,采用光滑粒子流体动力学方法（SPH）与有限单元法（FEM）建立了考虑流固耦合效应的高速液固撞击数值模型,详细分析了直径为2mm、撞击速度为1000m／s的液滴和射流对有机玻璃（PMMA）的三维撞击和破坏状况．分析表明：射流与液滴在撞击初始时刻的前缘变形和内部压力分布几乎是完全相同的;液滴撞击固体的最大压力值出现在0．20μs时,但此时材料内部最大等效应力只有104MPa,材料还不足以发生破坏;产生于液固撞击瞬时后0．32μs、速度高达2925m／s的侧向射流是使固体表面产生破坏的主要原因,因此撞击最初的破坏位于以撞击中心为圆心的一个圆环区域处．所得材料表面损伤情况与Brunton的实验数据吻合良好,证明了数值模型的可行性和精确性．     

SPH方法在两相流动问题中的典型应用

应用光滑粒子法(SPH)对具有间断面的不同介质的流动问题进行模拟分析.在现有的SPH方法的理论基础上进行修正,通过对密度近似方程进行正则化处理,在理想气体状态方程中引入Van Der Waals修正项以及对物质点运动速度进行修正,模拟了密度相差较大的两相流动问题,并且使用图解分析了两相交界面在模拟过程中的清晰程度.通过气(液)泡上升以及两相溃坝问题对SPH方法在两相流动问题模拟中的可行性进行了分析.最后应用修正后的SPH方法有效再现了两相流动的物理过程,并且能清晰分辨不同介质之间的交界面,从而说明SPH方法能有效对两相流动进行模拟.     

基于SPH方法的渤海海冰动力学数值模拟

简要介绍了光滑质点流体动力学 (smoothed particle hydrodynamics,简称SPH) 方法的基本原理,并将此方法应用到渤海海冰动力学数值模拟中. 在对海冰动力学方程和粘弹塑性本构方程进行SPH处理的基础上,对渤海海冰的动力演化过程进行了48 h的数值模拟,得到了海冰厚度、密集度、速度和主应力的演化过程及分布规律,同时对模拟的海冰厚度等值线分布与卫星遥感图像进行了对比分析. 数值模拟结果表明,SPH能够精确地模拟海冰分布规律以及冰缘线的位置,具有很高的计算精度和稳定性,是一种有效的海冰动力学数值模拟方法.     

一种光滑粒子流体动力学-有限元法转换算法及其在冲击动力学中的应用

为了解决冲击动力学中的大变形问题,提出了一种光滑粒子流体动力学-有限元法(SPH-FEM)转换算法,以等效Mises应力作为转换判据,将冲击过程中局部大变形区域的有限元网格转换为SPH粒子.该算法在大变形区域使用具有优势的SPH,在小变形区域使用精度和效率更高的FEM,为冲击动力学问题的数值计算提供了一条有效途径.使用SPH-FEM转换算法对圆柱形钢弹正冲击钢板发生冲塞破坏的过程进行了三维数值计算,计算结果与实验吻合较好,显示了该算法在计算精度方面的优势.在实际工程中,需要根据具体材料的失效模式,选择更加合适的转换判据.     

应力波和层裂计算中的光滑粒子法

采用光滑粒子法,数值模拟了平板碰撞下一维弹塑性应变波的传播及其在自由面和固定端的反射.通过引入含损伤的应变率相关的本构方程,对铝-锂合金的层裂实验进行了数值模拟.计算结果表明,光滑粒子法能很好地模拟应力波的传播和相互作用,层裂的计算结果与实验结果也比较吻合.