PTT黏弹性流体的光滑粒子动力学方法模拟

运用光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)方法对基于PHan-Thien-Tanner (PTT)模型的黏弹性流动进行了数值模拟. 首先, 利用SPH方法模拟了基于PTT模型的平板 Poiseuille流, 通过与文献结果的比较, 验证了SPH方法模拟黏弹性流动的准确性和有效性; 随后, 基于PTT模型对黏弹性自由表面流-液滴碰撞问题进行了SPH模拟, 研 究了PTT模型中拉伸参数对碰撞过程的影响. 为了解决张力不稳定问题, 采用简化的 人工应力公式. 数值结果表明, SPH方法可有效而灵活地模拟黏弹自由表面流问题.     

溃坝问题的粒子方法模拟

溃坝问题包括水体自由表面变形、翻卷、破碎、复杂的紊流和漩涡,涉及非线性复杂水动力学问题,含有丰富而复杂的流动机理。本文改进了传统的物质点方法,将物质点法从模拟固体问题拓展到了流体问题。采用改进的光滑粒子流体动力学和物质点方法,对两种长宽比的水坝坍塌过程进行了数值模拟,得到了波浪前沿及自由表面顶面随时间的演化过程。计算结果表明,两种粒子方法模拟的水体发展过程与实验结果吻合较好,较VOF更接近实验值。对涌波与竖直壁面的作用进行了计算,结果表明,SPH方法与其它方法相比能更好地捕捉水腔形态。最后分析了两种粒子方法在处理流体大变形问题时的特点以及计算的准确性和精度。     

非等温黏弹性复杂流动的改进SPH方法模拟

非等温黏弹性流体广泛存在于自然界和工业生产中,准确预测黏弹性流体的非等温流动机理和复杂流变特性有着重要的应用价值.文章提出一种改进的光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法对非等温黏弹性复杂流动进行了数值模拟,其中流体的黏弹特性通过eXtended Pom-Pom本构模型来表征.为了提高模拟结果的精度,采用了一种核函数梯度的修正算法;为了灵活地施加边界条件,发展了边界粒子和虚拟粒子相联合的边界处理方法;为了消除流动过程中的拉伸不稳定性,施加了粒子迁移技术.运用改进SPH方法数值模拟了液滴撞击固壁和F型腔注塑成型问题,通过与Basilisk软件得到的结果进行比较验证了改进SPH方法求解非等温黏弹性流体的有效性.通过利用不同粒子初始间距进行计算,评价了改进SPH方法的数值收敛性.研究了非等温流动相较于等温流动的不同流动特征,深入分析了不同热流变参数对流动过程的影响.数值结果表明,文章提出的改进SPH方法可稳定、准确地描述非等温黏弹性复杂流动的传热机理、复杂流变特性和自由面变化特性.     

Oldroyd-B黏弹性液滴弹跳行为的改进SPH模拟

基于光滑粒子流体动力学SPH(Smoothed Particle Hydrodynamics)方法对Oldroyd-B黏弹性液滴撞击固壁面产生的弹跳行为进行了模拟与分析。首先,为了解决SPH模拟黏弹性自由表面流出现的张力不稳定性问题,联合粒子迁移技术提出了一种改进SPH方法。然后,对Oldroyd-B黏弹性液滴撞击固壁面产生的铺展行为进行了改进SPH模拟,与文献结果的比较验证了方法的有效性。最后,通过降低Reynolds数捕捉到了液滴的弹跳行为;并在此基础上,分析了液滴黏度比、Weissenberg数和Reynolds数对液滴弹跳行为的影响。结果表明,改进SPH方法可有效地模拟黏弹性自由表面流问题;液滴黏度比、Weissenberg数和Reynolds数对液滴最大回弹高度均有显著的影响。     

DSPH方法的有效性验证及应用

讨论了不连续的光滑粒子流体动力学(discontinuoussmoothedparticlehydrodynamics,DSPH)方 法对粒子均匀性的敏感度,同时通过对DSPH 计算公式的详细推导,分析了该方法在界面处提高计算精度的 本质原因。通过简单的对碰模型,将DSPH 方法应用于实际问题中,研究了使用DSPH 方法计算密度、应变 率和能量时对模拟结果的影响,将多维的DSPH 方法推广到实际的物理计算中。     

基于SPH-FEM耦合方法的柔性导爆索分离装置爆炸分离过程数值模拟

为深入研究柔性导爆索在爆炸分离装置中的作用过程和机理,提出一种改进的光滑粒子流体动力学方法（smoothed particle hydrodynamics, SPH）与有限单元法（ finite element method, FEM）耦合算法。新方法中不仅包含导爆索模拟的SPH方法与分离装置模拟的FEM方法之间的接触算法,同时将完全损伤失效后的单元采用转化算法动态转化成SPH粒子继续参与计算,转化后的粒子与未转化的有限单元之间采用接触算法计算。采用该方法对环型和平板型两种爆炸分离结构的分离过程进行了数值模拟,验证了新方法的准确性与问题适用性;分析了分离板的变形断裂及损伤碎片的飞溅过程,得到了分离装置表面不同时刻的应力分布、损伤因子的变化趋势、von Mises应力的变化趋势;探讨了炸药在不同比内能情况下单元的屈服损伤速度、碎片的飞溅位移速度。     

流固耦合破坏分析的多分辨率PD-SPH方法

流固耦合破坏是一类涉及结构变形与破坏以及复杂自由表面现象的强非线性力学问题. 结合近场动力学(peridynamics, PD)与光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)各自的优势并考虑其计算效率问题, 提出一种适用于分析流?固耦合破坏问题的多分辨率PD-SPH混合方法. 分别采用SPH和PD方法以不同的空间和时间分辨率对流体和结构进行离散与求解, 利用具有与流体粒子相同光滑长度的虚粒子处理流?固界面, 以高精度满足界面边界条件. 通过两个经典算例: 液柱静压力下弹性板的变形和溃坝流体冲击弹性闸门的变形问题, 表明提出的多分辨率PD-SPH方法兼具较高的计算精度和计算效率; 对含裂缝的Koyna重力坝水力劈裂问题进行模拟计算, 所得裂缝扩展路径与文献结果吻合, 说明该方法适用于涉及结构破坏的流固耦合问题仿真. 最后尝试采用该方法进行流体冲击作用下含裂纹混凝土板崩塌过程数值仿真, 准确描述混凝土板的断裂破坏和全过程中的流体运动. 多分辨率PD-SPH混合方法或可为流?固耦合作用下的结构损伤破坏仿真提供一种新选择.      

基于修正光滑粒子流体动力学算法的低能量耗散数值波浪水池开发

目前, 无网格光滑粒子流体动力学SPH粒子法在波浪与结构物相互作用研究方面得到广泛应用, 但该方法模拟波浪远距离传播时, 常常面临严重的能量耗散问题, 导致波高非物理性降低, 给大范围海域、长时间作用下的波-物耦合作用研究带来一定困难. 对此, 本文采用一种核函数修正算法, 在确保粒子间相互作用对称性的同时, 改进压力梯度离散项的计算精度, 设法解决SPH方法中能量非物理性耗散的难题. 相较于前人减缓能量非物理性衰减的方法, 本文的修正SPH算法避免了自由液面搜索等复杂处理过程, 并能保证动量守恒特性. 数值结果中, 采用振荡液滴、规则波、不规则波等算例, 验证本修正SPH算法的准确性和有效性. 结果表明, 该修正SPH算法能准确模拟振荡液滴形态变化, 且动能保持较好守恒性. 通过数值水池与物理水池两者规则波与不规则波结果的对比分析表明, 基于本文修正SPH算法建立的数值波浪水池具有较好的抗能量衰减效果, 能实现长时间、远距离波浪传播的准确模拟. 此外, 本算法能在低光滑长度系数条件下, 实现精确模拟, 将极大缩减三维SPH模拟的时间, 从而节约计算成本.      

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

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

弹丸侵彻混凝土的SPH算法

给出了弹丸侵彻混凝土的数值计算,其中弹体作为刚体处理并划分成Lagrangian标准有限元网格,而混凝土划分成光滑粒子并经历大应变、高应变率和高压作用。为了描述混凝土的非线性变形及断裂特性,在计算中引入了Holmquist Johnson Cook本构模型及损伤模型。计算结果与实验结果的对比表明,利用光滑粒子流体动力学方法对混凝土材料进行大应变、高应变率的变形计算是有效的,并可避免网格重分或网格消蚀。此外,将光滑粒子流体动力学方法和有限元方法结合可以保持计算过程中材料界面的清晰。     

流体与结构相互作用问题的FE-SPH耦合模拟

应用有限元(FE)-光滑粒子流体动力学(SPH)耦合法模拟了具有自由表面的不可压流体与结构的相互作用问题.流体和结构分别采用SPH法和有限元法同时求解,两者在交界面处的相互作用通过接触算法进行处理.为了避免隐式计算压力,通过引入人工压缩率,将不可压流体近似为人工可压缩流体.采用FE-SPH耦合法对弹性板在随时间变化的水压作用下的变形以及倒塌水柱冲击弹性结构两个问题进行了模拟.模拟结果与实验结果以及其他已有数值结果符合良好,说明FE-SPH耦合法用于流体与结构相互作用问题的模拟是可行和有效的.     

NUMERICAL SIMULATION OF WIND-BLOWN SAND MOVEMENT BASED ON SPH

运用光滑粒子流体动力学(smoothed particles hydrodynamics, SPH)方法对沙粒和气流的相互耦合运动特性进行了分析,研究提出了风沙流的SPH数值方法并进行了数值模拟. 首先提出了风沙流的SPH建模方法和基本理论,建立了风沙流动的SPH数值模拟平台. 其次通过建立风沙流的SPH 模型并施加边界条件,对自然风作用下沙粒的运动情况进行了数值模拟,详细分析了沙粒运动轨迹及特性,最后通过与相关研究成果对比分析,验证了完善后的SPH方法有效性. 通过考虑气流场的可变性,在风沙流SPH计算模型中引入了加载(起风)和卸载(停风)方式,观察并对比分析了沙粒的运动轨迹和特性. 为进一步研究风沙流的实时动态非线性行为提供了SPH理论基础和数值分析方法.     

碎片云SPH方法数值模拟中的材料失效模型

光滑粒子流体动力学方法（smoothed particle hydrodynamics,SPH）被广泛应用于薄板超高速撞击碎片云的数值模拟。利用AUTODYN软件中的SPH模块,考察了无失效模型、Grady失效模型和最大拉应力失效模型3种方案下碎片云模拟结果,发现无失效模型时计算结果及材料表现与实验明显不符;相比于Grady失效模型,最大拉应力失效模型下材料更难失效,将小幅度减弱碎片云扩散程度,碎片总数减少,粒子聚集产生更大碎片,碎片云侵彻性能提高,增大模型失效阈值亦有上述表现。相比而言,Grady失效模型计算结果更符合实验,但2种失效模型间差异与撞击工况相关,材料破碎越充分差异越小。     

应用基于GPU的SPH方法模拟二维楔形体入水砰击问题

光滑粒子流体动力学(SPH)法是一种无网格的拉格朗日效值方法,广泛应用于计算流体领域模拟复杂自由表面流问题.SPH方法的主要缺点就是计算量过大,而基于GPU的并行计算方法可使SPH计算得到有效加速.本文应用基于GPU的SPH并行计算方法研究了二维楔形体的入水砰击问题.数值计算结果与文献中对应的解析解比较一致,验证了基于GPU的SPH方法的精度和可靠性.仿真结果同时显示基于GPU的并行计算方法可使SPH计算速度得到显著提高.     

一种SPH应力修正算法及自由表面流中的应用

提出了一种SPH应力修正算法,即模型中的拉应力和压应力分别采用不同的插值核函数和状态方程来处理,改善应力稳定性问题。介绍了一种改进的Quintic核函数,用于改善模型中压应力的稳定性。通过增加钟型核函数的光滑长度,改善模型中拉应力的稳定性。采用该应力修正算法模拟了无重力条件下方形液滴的震荡变形过程,对比分析了不同算法的模拟结果。此外,为进一步验证算法的适用性,模拟了溃坝算例。研究表明,改进的Quintic型核函数明显改善了粒子聚集现象,该SPH应力修正方法可以使液滴具有更均匀的粒子分布以及更光滑的自由表面,有效改善了SPH方法中的压应力不稳定作用以及自由表面流的模拟精度。     

基于光滑粒子流体动力学的滑油泵数值仿真

基于光滑粒子流体动力学方法,构建齿轮泵壳体及内部流体的粒子模型,对内啮合齿轮泵在不同工况下的流量特性进行数值模拟。首先设置了均匀分布的油泵出入口压强,计算得到的流量结果与试验结果吻合,且在中低转速下流量与转速呈线性关系;针对高转速下油泵流量降低的问题,通过适当减小周期模型内部的流体粒子数反映流量降低,获得了与试验值相吻合的连续转速流量的模拟结果;针对影响齿轮泵性能的间隙和空化现象,结合模型特点给出了相应的近似处理方法。通过以上研究,将SPH方法成功地应用于滑油泵问题的分析计算。     

RESEARCH DEVELOPMENTS OF SMOOTHED PARTICLE HYDRODYNAMICS METHOD AND ITS COUPLING WITH FINITE ELEMENT METHOD

拉格朗日型的有限元法和光滑粒子法在模拟材料大变形问题时各存优缺点, 而有限元与光滑粒子耦合算法实现了在小变形区域采用有限元法计算, 在局部的大变形区域采用光滑粒子法计算, 有效地综合了有限元法计算效率高和光滑粒子法能够自然地模拟材料大变形问题的特点.重点论述了有限元法、光滑粒子法以及有限元与光滑粒子耦合算法的研究现状及应用进展, 并讨论了各方法中需要进一步解决的问题.最后通过算例对3种方法的计算精度和计算效率进行了分析, 供研究人员参考.     

激光激励微型硅悬臂梁的振动特性研究

本文用光学探测方法研究了半导体硅悬臂梁的振动问题;运用基于外差干涉原理的实验装置得到了悬臂梁在激光激励下的振动响应（振动振幅和相位随调制激光频率的变化）;采用等离子波和热弹性波的数学模型,对悬臂梁的振动进行了理论分析;可看到实验与理论模拟结果吻合很好,同时通过分析可得振动相位与调制激光频率的平方根之间有线性关系。关键词词: 光学探测, 硅悬臂梁,振动.      

避免核函数导数的二阶导数核近似方法的修正

采用核近似光滑粒子流体力学SPH(Smooth Particle Hydrodynamics)方法计算一元函数二阶导数和多元函数二阶偏导数,对避免核函数导数算法进行了改进,提出了对核函数光滑长度处处具有o(h)精度的算法,并分别推导出一元和多元函数的修正公式。用不同的粒子间距和不同的光滑长度进行了二阶导数的计算和热传导问题的模拟,进行了修正方法与原方法的误差分析。结果表明,本文提出的修正措施在提高精度、减少误差及加快收敛速度等方面起到很大的作用。     

固体介质中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稳定性分析的思想, 以及国内外在这几个方面的研究成果及其最新进展, 同时指出目前研究中存在的问题和研究的方向.