基于CFD-DEM算法的固体火箭发动机气-固两相流模拟

针对固体火箭发动机中含金属推进体系中颗粒冲刷烧蚀及推力两相流损失的难题,基于连续相-离散元(CFD-DEM)耦合模型,考虑颗粒间的碰撞力与非碰撞力以及颗粒-气流作用力等,对固体发动机中气-固两相流进行了数值模拟研究。研究结果表明,相比于传统的双流体模型及轨道法,CFD-DEM能提供更丰富的粒子尺度信息包括粒子的运动轨迹、颗粒间的碰撞、颗粒受力情况等等;以及粒子相与气相相互作用过程,粒子相与壁面相互作用等。粒子的碰撞在喷管收敛段壁面、喉部区域以及喷管扩张段中心区域发生,与壁面的碰撞在收敛段壁面位置,且颗粒在该区域角速度较大,气相-颗粒曳力最大值出现在气相加速区域。单个粒子尺度的流场数据为固体发动机气固两相流流场的认识提供了更加丰富的信息,并为粒子聚集聚合及对烧蚀层冲刷等提供研究基础。     

复合材料宽弦风扇叶片模态仿真分析

针对复合材料宽弦风扇叶片,以铺层结构设计信息为基准建立了1种有限元建模方法,对某型复合材料风扇叶片模型进行了固有频率与模态振型的数值仿真计算,并采用激光测振仪对该风扇叶片的固有频率进行了试验测量。将试验测量与数值仿真结果进行对比分析,前3阶固有频率的误差在5%以内,第4阶固有频率的误差为6.6%。结果表明:该有限元建模方法建立的复合材料风扇叶片模型与实际叶片较吻合,基本满足工程分析要求,可以为复合材料风扇叶片的优化设计提供参考。     

光纤光栅传感器对非均匀应变场的响应特性

提出利用测量位置处光栅传感器的光谱特征反推该处对应的应变分布状态的方法。采用改进的传递矩阵方法,模拟了不同应变分布函数对应的光谱响应,分析了其基本特征参数同应变分布之间的变化规律。通过有限元仿真计算试验件在弯曲载荷下光栅粘贴处的应变分布,重构了传感器的光谱,并将其和应变分布函数的理论计算光谱图进行比较,验证了应变分布形式与反射光谱特征的关联性。构建带有不同半径圆孔的铝合金板孔口弯曲加载试验系统,测试试验件弯曲挠度和孔径大小与反射谱特性的关系。理论及试验结果表明,实测与仿真计算的谱图变化具有相同规律,光谱特征参数可作为预测铝合金孔径的有效指标,并为其他类型结构的应变场变化规律的监测提供参考依据。     

复合材料壳体固体火箭发动机温度载荷下力学响应分析

基于非线性有限元数值仿真分析方法,使用有限元计算软件ABAQUS,分别建立金属壳体与复合材料壳体两组发动机仿真模型,分析复合材料壳体固体发动机的结构完整性,并研究复合材料壳体的各项参数对于发动机装药和壳体的应力场的影响。分析结果表明:在温度载荷下,选用复合材料壳体的端燃装药固体发动机比选用金属材料壳体的发动机具有更好的结构完整性。随着复合材料壳体厚度的增大,装药的应力、应变值均增大,壳体的应力逐渐减小;随着复合材料壳体弹性模量增大或泊松比减小,装药的应力、应变逐渐减小,而壳体的应力、应变逐渐增大。     

AUTOMATIC FINITE ELEMENT MODELING OF A WING

研究了建立机翼有限元模型的自动化。用行列法标识结构的部件（肋、梁、蒙皮、立柱）,并用简单明了的表格界面来输入所有部件的几何、网格密度、材料、载荷、边界条件等参数。设计并用ＰＣＬ（ＰＡＴＲＡＮｃｏｍｍａｎｄｌａｎｇｕａｇｅｕｎｄｅｒＰＡＴＲＡＮ６．０）编制了相应的自动化建立整个机翼有限元模型的用户化程序模块。作为算例建立了ＶＦＷ６１４机翼的有限元模型并用ＮＡＳＴＲＡＮ６８进行了计算。结果表明,本方法有效和高效。     

不连续散粒体的离散单元法

建立了球体元三维离散单元法的单元接触模型。该模型把两个球体空间接触时的相互作用简化为四个弹簧阻尼器和三个滑动摩擦器。根据赫兹接触理论和非完全弹性碰撞理论,合理地确定了各接触参数。文中详细讨论了接触球体间相互作用力的分析计算方法和球体间产生接触的判断及检索方法,给出了球体元三维离散单元法的原理、算法和数据结构,并编制了相应的程序。利用该程序对球体自重下落和颗粒阻尼减振器分别进行了仿真计算,理论结果和实验结果有较好的一致性,验证了本文方法和程序的正确性。     

具有随机刚度薄板小挠度问题的数值分析

从薄板小挠度问题的控制方程出发,将结构的随机参数统一到刚度随机场中,采用局部平均法对刚度随机场进行离散,并将其融于二阶差分方法的列式中,给出了求解刚度随机时薄板上各点位移的数字特征的随机有限元法。该方法不但可以解决复杂边界下的薄板小挠度问题,而且进一步解决了结构的可靠度问题。文中最后对算例进行了计算,用解析解验证了本文方法的正确性。     

非确定结构系统区间分析的直接优化法

工程中的非确定性问题可以用区间分析、随机理论或模糊集理论进行求解。本文采用区间分析法来处理结构静力分析和设计中的不确定性问题。将结构系统中的不确定性参数用区间数来表示,用有限元法建立系统的控制方程。该控制方程是线性区间方程组。本文讨论了一些求解区间线性方程组的方法并提出了一种直接优化法。将方程组中的所有区间数都作为设计变量,区间量的变化区间作为相应的设计变量的边界约束,运用约束优化法求出方程组解的各元素的最大值和最小值。文中给出了两个算例,列出了本文算法与其他算法的结果比较     

深冲钢W08Al热轧模拟试验研究

在热模拟试验机上对 Ｗ ０８ Ａｌ深冲钢进行了热轧模拟试验,测量了材料在不同情况下的相变温度,考察了不同热轧工艺参数对晶粒尺寸和织构的影响。结果表明,铸坯经预先热变形后重新加热测量相变时,相变温度较低,而直接热变形后测量相变温度升高;在 Ａｒ３温度以上采用大变形量可以使晶粒细化,但在两相区终轧变形会造成混晶现象; Ａｒ３温度以上热变形所得的变形织构是弱的,进入两相区终轧变形则不利织构（１００）成分很快发展起来,并且织构的取向密度沿板厚方向明显不均匀。最后对 Ｗ ０８ Ａｌ热轧合理热轧工艺提出了建议。     

半穿甲数值计算初步探讨：正碰撞均质钢甲计算

半穿甲子弹,又名穿甲爆破弹,其在侵彻一定厚度装甲后爆炸以起到杀伤破坏作用,因此子弹的侵彻机理和弹体的动强度问题就显得尤为重要。利用ANSYS／LS-DYNA3D计算软件对半穿甲子弹正碰撞均质钢甲的侵彻过程进行了数值模拟,并与试验结果进行了对比。