黏弹性接触界面端附近的奇异应力场

研究蠕变加载条件下线黏弹性材料接触界面端附近的奇异应力场问题.考虑接触界面的摩擦,假设界面端的滑移方向不改变,相对滑移量微小,且其与位移同量级,由此线性化局部边界条件,根据对应原理得到Laplace变换域中的界面端应力场,导出时域中奇异应力场的卷积积分表达式.对卷积积分核函数进行数值反演,考虑接触材料的两类组合,一是持久模量具有量级上的差异,另一是持久模量接近相同.算例结果证实核函数可以用准弹性法求得的解析式较准确地近似.在此基础上,利用积分中值定理,并引入各应力分量的修正系数,得到黏弹性奇异应力场的简化式.结合核函数的数值反演结果分析修正系数表达式的取值范围,得到如下结论,若两相接触材料的持久模量相差很大,可以采用准弹性解的解析式较准确地描述界面端的奇异应力场;一般情况下,应力场不存在统一的奇异值和应力强度系数,当采用类似于准弹性解的表达式近似给出黏弹性应力场时,可以估计此近似描述的误差限.文中最后采用有限元分析黏弹性板端部嵌入部位的应力场,算例包括了黏弹性板与弹性金属支承、黏弹性板与黏弹性垫层所形成的滑移接触界面端,利用黏弹性有限元的数值结果验证理论分析所得结论的有效性.      

HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics,CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度. 其次,通过NASATurbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶). 最后,通过NASA CommonResearchModel标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.      

圆柱形汇聚激波诱导Richtmyer-Meshkov不稳定的SPH模拟

汇聚激波诱导不同物质界面的Richtmyer-Meshkov(RM)不稳定现象在惯性约束核聚变领域有重要的学术意义和工程背景.基于网格离散的宏观流体力学方法由于数值扩散问题往往需要高阶精度算法才能准确追踪界面演化,且对大变形和破碎合并等复杂界面追踪也极为困难.光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法采用纯拉格朗日算法,可以有效克服上述难点.但经典SPH算法需采用人工黏性处理强间断,在激波间断处往往会出现严重的非物理振荡,对于涉及强冲击不稳定性问题,很难达到理想的模拟效果.本文采用基于HLL黎曼求解器的SPH算法,实现了对强激波和大密度比物质界面的有效分辨和追踪.一维数值校核证明了代码的可靠性、健壮性,并进一步模拟了二维圆柱形汇聚冲击波冲击四边形轻/重气界面诱导的RM不稳定性问题,与已有实验结果进行了对比,发现模拟结果与实验结果吻合.通过分析界面演化过程中的密度及压力变化,发现本文所采用的方法可准确地追踪激波与界面作用的复杂界面和波系演化规律.研究结果为进一步理解和解释汇聚冲击条件下的RM不稳定性机理奠定了基础.     

具有恒定冲击载荷的梯度泡沫金属材料设计

多胞材料可通过大变形大量地吸收冲击能量，引入密度梯度可进一步提高其耐撞性。梯度多胞材料的宏观力学响应对材料密度分布极为敏感，不同类型的细观构型的影响也极为不同。已有的研究工作主要局限在对给定的密度梯度分析其动态响应，较少对耐撞性设计方法进行研究。本文针对梯度闭孔泡沫金属材料，基于非线性塑性冲击波模型发展了耐撞性反向设计方法，以维持冲击物受载恒定为目标，运用级数法获得了简化模型和渐近解。利用变胞元尺寸法构建了连续梯度变化的三维Voronoi细观有限元模型，并利用ABAQUS/Explicit有限元软件对理论设计进行数值验证。结果表明，反向设计理论简化模型的渐近解对于梯度闭孔泡沫金属材料的耐撞性设计是有效的，所提出的耐撞性设计方法在控制冲击吸能过程和冲击物受载方面具有指导意义。     

旋转射弹高速倾斜入水多相流场与弹道数值模拟

基于VOF多相流模型和有限体积法求解水、汽、气多相流动的RANS方程，结合重叠网格技术和six DOF算法对某一型号舰载射弹倾斜入水过程进行数值模拟研究。首先基于该方法研究了射弹旋转效应对射弹运动特性及流体动力特性的影响，然后对不同入水角下倾斜入水过程进行分析，得到不同倾角下旋转射弹入水空泡形态发展规律、弹体运动特征及流体动力特性变化规律。研究结果表明:射弹的旋转有利于弹体在初始对称面内的弹道稳定性，但会降低弹体侧向稳定性，使射弹受到的阻力系数、俯仰力矩系数变小；入水角越小，形成的空泡越不对称，由射弹运动状态的改变引起的空泡形态变化越明显，在超空泡航行阶段，弹体运动较稳定，不同角度下流体动力系数差别很小，当弹体下表面刺破空泡壁沾湿时，弹体运动状态发生较大变化，流体动力系数迅速增大，此时入水角度过小，弹体容易失稳；弹体的沾湿对空泡形态、弹体运动稳定性和流体动力特性有着重要的影响。     

等热流密度加热垂直下降液膜的研究

利用等热流密度加热条件下降膜流动的三维模型方程进行线性稳定性分析和数值模拟。线性稳定性分析表明，模型方程在小到中等Reynolds数下都适用，并且流向不稳定性增长率随着Reynolds数和Marangoni数增加而增加，展向不稳定性增长率则随着Marangoni数增加而增加，随着Reynolds数增加而减小，流向和展向对扰动波数都存在一个不稳定区间。三维数值模拟表明，在等热流密度加热条件下，液膜在随机扰动的情况下最终会形成带孤立波的三维溪流状结构，液膜与气体的换热也因溪流状结构的出现而加强；在随机扰动的基础上引入占优势地位的展向最不稳定扰动会使得换热增强，液膜会提前破裂；在随机扰动的基础上引入占优势地位的流向最不稳定扰动时，液膜的换热会增强，但不会提前破裂；在随机扰动的基础上同时引入占优势地位的流向和展向最不稳定扰动时，换热会加强且液膜会提前破裂。     

一种考虑界面相物性影响的界面单元法

界面是由复杂的界面相简化而成的,界面破坏实际是界面相材料的破坏。数值计算为了方便,如经典模型和内聚力模型等,都把很薄的界面相作无厚度化处理。导致只能考虑界面的面力,而无法考虑界面相内的应力(平行于界面方向的应力)。使界面失效准则先天性地排除了界面相内部应力的影响,从界面相材料失效机理的角度来看这是不够严谨的。本文将界面相材料等效为一种弹性连续体,由界面本构关系推导得到了一种新的界面单元。该单元具有界面参数易确定、对界面相物性可以进行等效描述等优点。通过商用有限元软件ABAQUS和用户子程序UEL实现了数值分析,并与直接物理模型的数值模拟结果进行对比,证明了本方法的简便及准确性。通过对不同界面相厚度结构的进一步分析,探讨了本文方法的可行范围。     

含参广义集值强向量平衡问题的稳定性

本文借助集合极限的性质和弱f-性假设证明了含参广义集值强向量平衡问题解集映射的下半连续性,其方法不同于最近文献(Zhao,2016和Meng,2018).此外,建立了含参广义集值强向量平衡问题解集连通性的充分条件,并举例验证了所得结果的正确性.本文得结果推广和改进了已有文献(Gong,2008,Xu,2009,Chen,2010,Xu,2013和Zhao,2013)中相应结果.     

Euler-Bernoulli梁的高阶二次摄动解及收敛性讨论

首次用解析的方式给出了Euler-Bernoulli梁后屈曲与非线性弯曲问题的高阶二次摄动解答.假定梁的中线不可伸长，用精确曲率公式与能量变分原理导出了非线性Euler-Bernoulli梁的模型.通过与精确解或高阶摄动解的比较，讨论了二次摄动解答的收敛性及适用域.得到主要结论如下:低阶摄动解适用于描述梁的初始后屈曲阶段及初始非线性弯曲阶段；更高阶次的摄动解适用于描述梁的深度后屈曲以及深度非线性弯曲.从这个意义上去说，该文不仅仅指出某些文献上的部分结果不精确是由于摄动解答超出了其特定的适用域，并且还进一步发展与完善了二次摄动法.     

求解矩阵方程AXB+CXD=F参数迭代法的最优参数分析

本文针对求矩阵方程AXB+CXD=F唯一解的参数迭代法，分析当矩阵A，B，C，D均是Hermite正（负）定矩阵时，迭代矩阵的特征值表达式，给出了最优参数的确定方法，并提出了相应的加速算法.