弹塑性流体力学拉氏自适应网格和网格重分守恒重映数值模拟方法

拉氏自适应重分弹塑性流体力学有限元程序实现了网格完全自适应，具有良好、灵活的非结构自适应网格数据结构，实现了滑移界面两边(接触间断)网格动态调整，网格的细分和合并处理灵活，网格重分和网格自适应模块兼容、守恒重映，网格重分中采用多种方法控制新网格的质量，爆轰计算可采用Lee-Tarver的化学反应率模式。初步数值计算结果表明，弹塑性流体力学拉氏自适应重分数值模拟方法合理，计算结果正确，基本反映了流场的物理结构。     

钽的层裂实验数值模拟

 采用连续介质力学基唯象模型模拟分析了钽的平板撞击层裂行为。该模型包括了材料的非线性弹性（状态方程）、率相关塑性和孔洞的形核及生长等多种效应，并且采用一种对角隐式Runge-Kutta方法来求解本构率方程组，提高了热粘塑性本构关系计算的稳定性及精度。将数值模拟结果和相关实验数据进行了对比分析，结果表明，对于样品中的拉应力峰值明显高于材料层裂强度的实验（中、高速平板撞击实验），理论模型具有较好的预估能力，但对于临界层裂问题（低速平板撞击实验），该模型对材料损伤与失效过程的描述可能不够准确，需要进一步改进。     

一种岩石损伤本构模型在地下强爆炸中的应用

为了描述地下爆炸波的传播,建立了一种考虑塑性硬化、剪胀和损伤软化效应的岩石本构模型,并应用于地下强爆炸自由场的数值计算,获得的速度和位移波形以及峰值速度衰减曲线等与国外地下试验测量数据和同类计算相比十分接近,从而验证了本构模型的有效性。     

欧拉流体力学数值方法中的界面计算

本文所提算法适用于二维和三维多介质流体力学两步欧拉数值方法中输运计算的混合网格(包括自由面网格)界面处理。在一个混合网格中，界面被近似地看作直线(二维)或平面(三维)。整个方法分为三步：(1)第一步，用混合网格周围的八个网格的介质面积份额(二维)或二十六个网格的介质体积份额(三维)确定界面的法线方向；第二步，用混合网格的本身的介质面积份额(二维)或体积份额(三维)确定界面的方程(位置)；第三步，用此直线方程求出通过网格边界的流以及下一时刻网格的面积份额(二维)或体积份额(三维)。最后给出了用此方法所做的一些数值计算及与SLIC算法的比较。     

用浮阻力模型研究Richtmyer-Meshkov不稳定性诱导混合

采用浮阻力模型对激波管低压缩和激光加载高压缩情况下的Richtmyer-Meshkov不稳定性诱导混合现象进行了研究。通过与实验和理论分析结果进行比较发现:为了达到好的吻合, Richtmyer-Meshkov不稳定性情况下阻力系数的取值范围(2.0~5.36)比Rayleigh-Taylor不稳定性情况下的值(3.3~4.0)宽得多; 而在Richtmyer-Meshkov不稳定性情况下, 高压缩时阻力系数的不确定度(约为3.36)明显高于低压缩时的值(约为1.46), 模型的进一步完善还有待于更精确实验的验证。研究显示:指数律经验公式中指数随工况的不同而显著变化, 目前工程设计中采用指数律经验公式是粗糙的。     

三维光滑接触界面上的分配参数滑移面算法

针对三维滑移界面上的分配参数滑移算法,提出一种三维滑移面光滑目的 ,仅根据单元面顶点的坐标和法向等局部信息,构造出曲面片边界G¹连续(切平面连续),曲面片内C¹连续的三维光滑曲面,此目的 为实际的滑移面提供一个准确的几何表达.在此光滑曲面上构建了具有C¹连续的法向压力场和质量面密度场.然后在光滑滑移界面上通过Newton-Raphson迭代计算接触点,依据分配参数滑移算法施加接触约束和接触计算.数值算例显示光滑的接触滑移面可以减少滑移节点的振动,提高计算的收敛性.     

卵形杆弹对铝靶的斜侵彻

在Lagrange有限元基础上,扼要介绍了一种适用于三维斜侵彻数值模拟计算的滑移面处理技术。该方法放弃了传统滑移计算中从单元的设置,而代之以预设从节点,从而避免三维斜侵彻数值计算中主从单元的相交、滑移面的识别、修正与再定义的困难。在确保计算精度的同时,有效地提高了计算效率。卵形杆弹对铝靶侵彻的系列数值模拟计算表明,无论对于正撞击或斜撞击,本文中所获结果与实验结果都有良好的一致性,这说明本文中所述方法和所建程序的合理性和有效性,为侵彻贯穿过程的数值分析提供了一种实用和有效的手段。     

圆筒试验JWL状态方程参数的贝叶斯标定

在研究数值模拟的输入参数引入的不确定性时,通常需要人为给定每个输入参数的概率分布,且输入参数概率分布的选择可能会对分析结果产生直接影响。利用贝叶斯方法标定了圆筒试验JWL状态方程参数,得到了标定参数的估计值和后验分布,并研究了不同统计模型假设对标定参数的估计值和后验分布的影响。贝叶斯后验分布融合了基准试验的试验数据的信息,因此将其作为不确定度量化分析时输入参数的初始概率分布,可以尽量减少分布选择引入的认知不确定性。     

基于多项式混沌的全局敏感度分析

回顾基于多项式混沌和方差分解的全局敏感度分析方法,针对高维数随机空间和高阶多项式混沌展开面临的“维数灾难”问题,采用回归法、稀疏网格积分及基于l₁优化的稀疏重构技术（即压缩感知技术）来减少非嵌入式多项式混沌方法所需的样本配置点数目.针对几个典型响应面模型（包括Ishigami函数、Sobol函数、Corner peak函数和Morris函数）进行Sobol全局敏感度指标计算,展示多项式混沌方法在基于方差分解的全局敏感度分析中的有效性.     

拉氏数值模拟中的人工粘性方法

较好的人工粘性需要满足较小的计算开销、不能去除真实具有的涡运动等条件.提出一种应用于拉氏数值模拟中基于Lew人工粘性,同时增加了限制器的人工粘性方法.可以有效减少数值模拟结果对网格的依赖;采用特征值限制器控制施加的人工粘性大小,通过限制器能够区分激波压缩和等熵压缩;方便应用在二维、三维,结构网格或者非结构网格上.