液Ar冲击压缩特性的理论计算

用exp-6有效两体势模型和液体变分微扰理论计算了液Ar冲击压缩曲线,在35GPa以下的压力范围内计算的冲击压缩曲线与Thiel及Nellis等人的实验数据及其它理论的计算结果符合较好。计算结果表明文中所选的势较为准确地反映了液体分子间的相互作用。也对较高冲击压力下理论计算的冲击曲线和实验结果之间的偏差作了分析,结合不透明度实验的结果,我们认为当压力超过35GPa,温度在12000K以上时,液Ar体系电子激发对系统热力学状态有较大影响。     

多功能含能结构材料冲击压缩特性的理论计算

基于混合物冷能叠加原理,由各组分Hugoniot数据计算了密实材料的冲击压缩特性。再从等压路径出发,结合Wu-Jing模型由热力学关系得到了具有一定孔隙率多功能含能结构材料的冲击压缩特性计算方法。以W/Cu、Al/Ni、Ni/Ti和Al/Fe2O3/epoxy等典型颗粒金属材料及含能金属材料为例,计算了其冲击压缩过程中相关Hugoniot参数。计算结果与已有实验结果吻合较好,多功能含能结构材料冲击压缩特性受材料孔隙率、材料配比等影响明显。     

甲烷气体的冲击化学反应流研究

用二级轻气炮加载技术将W89Mo9Ni1Fe1飞片加速到约5.0km/s,撞击封装有接近常态的甲烷气体的LY12铝靶。利用六通道高温计和示波器记录到冲击压缩的甲烷气体的光辐射历史曲线,并获得甲烷气体的冲击波速度与温度。利用一维Euler方程与化学反应方程的全流场耦合,采用甲烷体系的13组元40反应步的化学反应模型,并用3阶WENO数值格式对空间进行离散。为避免刚性过强,化学反应源项用点隐方法处理。由此得到了甲烷的冲击反应流场和波后热力学参数。通过比较,数值结果与实验结果符合较好。最后,解释了甲烷用作电探针保护气体的物理机制。     

超高压范围内金属铝的Hugoniot线和等熵压缩线计算

由状态方程数据库出发,结合热力学关系,计算得到超高压范围内金属铝的冲击绝热线和等熵压缩线,与国外发表的文献数据基本一致,进一步说明沿着冲击绝热线对材料进行压缩会遇到一个压缩极限,对金属铝的极限压缩比约为4．97,而沿着等熵压缩线则不会遇到压缩极限问题,通过等熵压缩加载可以实现比冲击加栽更大的压缩比。     

不同应变率下冰破坏特性的试验研究

为研究寒区流冰对水工/船舶结构的影响,需要获得不同流速冰载荷的特征,准静态下不同应变率条件下冰材料变形至破坏的特征就成为一项重要的基础性问题。为此,开展了低温环境应变率10-2s-1至10-4s-1下淡水冰的单轴压缩试验。试验发现:随着应变率减小,极限应力由18.51MPa降低至8.44MPa,达到极限应力后,试件承载能力由瞬间消失变为逐渐降低至稳态,试件破坏形貌由劈裂转变为周向膨胀,即由脆性转变为韧性破坏,转变应变率约为10-3s-1;在双对数坐标系中,单轴压缩强度随应变率的增加近似呈线性增大。通过对应力-应变曲线进行积分,给出了不同破坏形式下冰变形至破坏的临界应变能密度,发现脆/韧转变状态下加载至破坏所需能量最大,该能量特征主要由冰中微裂纹的萌生、断面摩擦、再结晶引起。     

气液两相涡街的数值计算

本文利用有限差分法,建立了模拟气液两相涡街的物理模型和数值计算方法,成功地利用数值方法模拟出了气液两相涡街及涡街中含气率分布等主要特征,该方法简便,理论计算出的结果与实验值相符较好     

吸能包装模型结构的冲击响应

利用空气炮冲击实验对吸能包装结构的跌落过程进行模拟,进行了缩比模型的正撞和30°斜撞实验,针对模型实验进行了数值分析,获得了吸能包装结构模型在撞击过程中的应力分布和塑性变形,并将计算情况与实验结果进行了分析。结果表明：在撞击中吸能包装结构主要通过缓冲木材的塑性变形及外钢壳屈曲产生的塑性铰吸收能量,塑性变形主要集中于撞击端,而远离撞击端未见塑性变形;计算中木材本构参数采用顺纹方向压缩应力应变曲线具有一定的有效性。

     

轴向冲击载荷下开孔薄壁方管吸能特性的研究

利用有限元软件ABAQUS对开孔薄壁方管在轴向冲击载荷作用下的变形规律与吸能特性进行了数值仿真分析,并根据数值模拟结果进行了多目标优化。在薄壁方管的每个侧面上开六个长方形孔,通过改变长方形孔的长度和宽度来研究开孔薄壁方管的耐撞性。结果表明:在研究的范围内,除开孔尺寸为20mm×5mm的薄壁方管外,其余开孔薄壁方管的破坏模式为外延式破坏模式,其变形过程为动态渐进屈曲;当开孔长度增大时,开孔薄壁方管初始峰值力减小;当开孔长度和宽度同时增大时,开孔薄壁方管平均载荷均减小;开孔面积越大,管吸收的能量和比吸能越小;开孔长度不变时,随着宽度的增大压缩力效率降低,当宽度不变时,随着长度的增大压缩力效率增大。基于开孔方管的数值模拟结果构造了响应面模型并对其进行了多目标优化,给出了Pareto前沿图,可根据各目标的偏好确定方管的开孔尺寸。     

激光冲击强化过程中蒸气等离子体压力计算的耦合模型

首先基于系统能量守恒条件,提出了一种计算蒸气等离子体压力的一维耦合计算模型。模型中不仅考虑了蒸气等离子体界面压力与质点速度的非线性效应,同时也考虑了界面烧蚀所致的运动速度,将蒸气等离子体的膨胀与约束介质的变形耦合。在耦合模型的基础上,采用显式差分计算程序与显式有限元计算程序LS-DYNA互相迭代求解的方法,对不同激光功率密度分布下的蒸气等离子体压力进行了计算。结果表明,计算结果与实验测量结果具有很好的一致性,证明了计算模型的合理性。     

基于光纤陀螺的温度补偿模型改进、压缩和FPGA实现

为提高光纤陀螺（FOG）在变温环境下输出误差的补偿精度,在长短期记忆神经网络模型（LSTM）基础上,利用分段非线性粒子群算法（PN-PSO）对LSTM模型超参数寻优,建立PN-PSO-LSTM光纤陀螺温度补偿模型。为有效降低计算和存储开销,便于部署在资源受限的硬件环境中,提出一套适用于光纤陀螺应用场景的模型压缩方案,包括：知识蒸馏、剪枝、激活函数线性化、定点数量化等。最后基于Xilinx公司某芯片完成部署。对比实验结果表明,相较于传统反向传播（BP）模型和传统PSO-LSTM模型,采用所提模型补偿后,陀螺零偏输出均方误差分别降低74.4%和53.5%,模型压缩后在大小减小94.1%的同时,陀螺零偏输出均方误差仍然比传统全精度模型更低,在FPGA实现后对比PC端模型推理速度提升98.47%。     

考虑侧面摩擦的侧限压缩模型的数值解法

海上箱筒型基础结构下沉到位后,舱内土体的受力变形可以看作是侧限压缩模型问题,上部荷载引起的侧壁摩擦力与土体应力之间的关系较为复杂。本文通过简化模型,列出了侧限压缩模型的平衡微分方程,分析了方程解析计算时存在的难点,结合摩擦应力边界条件的特点,构造函数迭代法对该方程进行数值求解,展示了计算的收敛过程,得到接触侧面上水平应力和土体表面位移的非线性分布规律,并将其结果与ABAQUS数值模拟结果进行了对比。结果表明,函数迭代法收敛速度较快,计算过程稳定,对其他存在变量耦合边界条件的微分方程数值求解有一定的参考价值。     

气液两相管流水力瞬变的数值计算及实验研究

由气液两相管流的基本方程出发,通过引入矢通量分裂,对传统的特征线差分做了较大的改进,形成了基于矢通量分裂的特征线差分解法。该法首先将控制方程组的特征值分解成正、负两部分,进而将控制方程中的矢通量雅可比矩阵分裂成两个亚矢量矩阵,对其按各自的迎风格式差分,从而建立了稳定的差分求解格式。该计算法适合于计算声速变化较大且计及液流速度的气液管流的瞬变。计算求解得到的各种不同初始空隙比的压力变化曲线、声速曲线、波速变化曲线、空隙比变化曲线及气体释放影响曲线,通过与不同初始空隙比时气液管流水力瞬变的实验结果对比分析,结果表明两者吻合较好,说明本文方法对于低空隙比的气液两相管流具有较普遍的适用性。     

考虑环境温、湿度的球形装药爆炸冲击波参数计算模型

为获得实际温、湿度环境中冲击波参数计算模型,计算了考虑温、湿度的理想气体状态方程参数,利用SPEED软件,针对典型状态空气中球形装药爆炸过程进行数值模拟,得到了温、湿度对爆炸冲击波参数的影响规律。结果表明,温度和相对湿度对冲击波超压的影响较小,而正压作用时间和冲量随温度和相对湿度的升高均呈线性递减关系,在高温高湿和寒冷干燥条件下,冲击波正压作用时间和冲量相差分别达21.8%和18.4%。以经典工程计算模型为基础,通过引入含有温度、湿度和对比距离的修正因子,建立了考虑环境温、湿度的球形装药爆炸冲击波参数的计算模型。采用该模型计算得到的不同药量球形TNT爆炸冲击波参数与数值仿真结果吻合较好,可对装药在实际环境中威力评估提供参考。     

多元混合PBX炸药冲击起爆的多元Duan-Zhang-Kim反应速率模型研究

提出了多元混合PBX炸药孔隙塌缩热点模型新的处理方法,构建了新的细观反应速率模型,系列数值模拟结果与实验结果均一致,表明该细观反应速率模型可较好地描述和预测炸药组分配比及颗粒度对多元混合PBX炸药冲击起爆过程的影响。PBX炸药冲击起爆过程主要受热点点火过程和燃烧反应过程共同作用：HMX占主导成分的PBXC03炸药,起爆压力低,冲击起爆过程受热点点火影响较明显,热点点火后的燃烧反应速度较快,表现为加速反应特性;TATB占主导成分的钝感PBXC10炸药,起爆压力高,冲击起爆过程主要受点火后的燃烧反应过程控制,且点火后燃烧反应速度较慢,表现为稳定反应特性。

     

复合材料层合板低速冲击损伤分析的连续介质损伤力学模型

针对复合材料层合板低速冲击损伤问题,提出了一种各向异性材料连续介质损伤力学模型,模型涵盖损伤表征、损伤起始判定和损伤演化法则3 个方面. 通过材料断裂面坐标下的损伤状态变量矩阵完成损伤表征,并考虑断裂面角度的影响,建立了主轴坐标系下的材料损伤本构关系. 损伤起始由卜克(Puck) 失效准则预测,损伤演化由断裂面上的等效应变控制,服从基于材料应变能释放的线性软化行为. 模型区分了纤维损伤和基体损伤,并根据冲击载荷下层内产生多条基体裂纹继而扩展至界面形成层间裂纹(分层) 的试验观察,引入基体裂纹饱和密度参数表征层间分层. 以[0₃/45/-45]_S 和[45/0/-45/90]_4S 两种铺层的复合材料层合板为例,预测了不同冲击能量下复合材料层合板的低速冲击损伤响应参数,试验结果证明了连续介质损伤力学模型的有效性.模型在不同网格密度下的计算结果表明单元特征长度的引入可以在一定程度上降低损伤演化阶段对网格密度的依赖性.      

基于连续-非连续单元法的三维脆性颗粒冲击破碎特性分析

通过在连续-非连续单元法(CDEM)中引入考虑应变率效应的断裂能本构以及能量统计算法,实现了球体冲击破碎过程中损伤破裂程度及能量演化的定量分析。计算结果表明,冲击破碎过程分为接触蓄能阶段、损伤破碎阶段和碎块飞散阶段。首先,颗粒的部分动能转化为单元弹性变形能,随后这部分变形能和动能迅速转化为摩擦消耗、阻尼消耗及弹簧断裂能,破碎基本完全后碎块继续飞散。不同冲击速度下,颗粒分别出现了反弹、开裂、破碎和粉碎的现象。随冲击速度的增加,D50的变化速率逐渐放缓,破碎块度逐渐趋于稳定;破裂度、损伤度以及平均损伤因子的变化速率先增加后放缓,颗粒破坏以拉伸破坏为主。以上结论可为脆性材料冲击破碎工艺的优化设计提供依据。     

基于总循环应变能的疲劳寿命计算模型

从总循环应变能的角度出发, 提出了一个新的疲劳寿命计算模型, 通过对30CrMnSiA 钢和AISI 4340钢两种材料的的疲劳寿命计算, 并与相关实验数据对比, 证明了该模型计算疲劳寿命较为准确.     

考虑当地涡旋运动修正的湍流模型在非定常空化湍流流场计算中的应用

为了更精确地模拟当地涡旋运动对复杂非定常空化湍流流动的影响,考虑当地涡旋运动特性对标准k-ε湍流模型进行了修正;并采用该修正模型对二维Clark-Y水翼周围非定常空化湍流场进行了计算;通过与实验结果进行对比,分析并验证了该修正模型计算结果的有效性。研究表明:该修正湍流模型根据当地涡旋效应对湍动能输运方程中的生成项进行修正,可以有效解决标准k-ε模型对湍动能的过度预测问题,得到的空化核心区域时均湍动能和湍流黏性的预测结果降低了约30%,与实际值更为接近;标准k-ε模型计算所得的时均空泡长度不大于0.6倍弦长,而修正模型计算所得的时均空穴可以覆盖水翼吸力面;修正模型可以准确预测水翼吸力面表面空泡的非定常脱落过程,包括附着空穴的拉伸和断裂,以及脱落空泡的逐渐消失。     

基于Level Set方法的气-液-固三相流动模型与模拟

采用基于Level Set方法与离散颗粒模型相结合的方法,建立了一个用于描述气液固三相流动的新模型.模型耦合了颗粒与气泡、颗粒与液相以及气泡与液相之间的相互作用.应用该模型对液固悬浮液中的典型现象--气泡的单孔及多孔形成过程以及颗粒夹带现象进行了三维模拟,检验了其可行性.并进一步研究了颗粒的存在对气泡的形成与上升过程的影响以及气泡诱导的液相流动对颗粒行为的影响.研究结果表明,所提出的模型能够真实地预测三相流中气泡与颗粒分散相的特征,为研究多尺度的三相流动提供了一种新途径.     

用三项式点火增长反应速率模型计算JB-9014的反应区

在一维流体动力学编码SSS程序中,利用三项式点火增长反应速率模型对钝感炸药JB-9014的反应区进行了计算。计算中,未反应炸药采用固态HOM状态方程,产物采用气态JWL状态方程,计算得到了钝感炸药JB-9014化学反应区的峰值压力、CJ压力、反应区宽度和反应时间等参数,计算结果与实验结果符合较好。分析表明,三项式点火增长反应速率模型可用于研究钝感炸药的反应区结构。