航天器舱体热密封结构多物理场耦合数值分析方法研究

综合考虑航天器舱体外围高超声速流动、密封带及舱体结构传热以及密封结构内部的空腔流动传热,提出航天器高温热密封结构的瞬态多物理场耦合分析方法。利用改进型Van Driest变换方法进行高超声速流动环境预测,基于高斯－赛德尔分块迭代耦合方法完成一体化耦合计算方法。采用包含不同材料结构部件的复合结构,计算了0～200 s内有无密封塞两种情况下的各部分结构的瞬态热传导过程。结果表明,密封塞的使用可显著降低空腔内的最高温度,瞬态变化情况的考虑更加准确地反映了各部分结构部件及内部空腔的温度变化情况。该文的计算方法可广泛应用于航天器热密封结构的传热特性分析,可为火箭等航天器上的高温密封部件设计提供有效的数值分析工具。     

ICVI工艺致密化C/C复合材料的数值模拟

为了研究等温化学气相渗透( ICVI)工艺制备C/C复合材料过程中预制体的致密化过程及流场和温度场的分布,利用COMSOL软件建立预制体致密化过程中传质、传热和孔隙率变化的多场耦合模型。以甲烷为前驱体,将动量和能量守恒方程进行耦合计算,计算结果表明,在传质最初始阶段,前驱体温度迅速升高至设定的沉积温度,且整个反应器内部温度分布均匀。根据以上计算结果,设定温度为定值,耦合质量、动量守恒方程和孔隙率变化方程,通过计算得到在开始致密化阶段预制体最大密度分布在预制体中部,随着致密化进行,该区域向外侧移动。致密化150 h后,不同时间预制体整体平均密度的计算值与实验值吻合较好,验证了致密化模型的可靠性。     

出口宽高比对S形二元收敛喷管雷达散射截面的影响

在S形二元收敛喷管进出口面积、偏心距、面积变化规律及中心线变化规律不变的条件下,采用自适应超椭圆方法设计不同出口宽高比的S形二元收敛喷管.基于迭代物理光学法与等效边缘电磁流法自主开发计算腔体部件雷达散射截面（RCS）的程序;然后通过文献中的实验数据验证了计算程序的准确性和可靠性.通过计算程序分别分析了水平、垂直两种极化方式下不同S形二元收敛喷管的边缘绕射场与总散射场的电磁散射特性.结果表明:在水平、垂直两种极化方式下,喷管出口宽高比的变化对S形二元收敛喷管边缘绕射场的RCS影响较小,不同出口宽高比的S形二元收敛喷管边缘绕射场的RCS相差不超过4dB.喷管出口宽高比的变化对总散射场的RCS影响较大;正探测角时,宽高比为1.5时,在大部分探测角范围内总散射场有较低的RCS;负探测角时,宽高比为3.5时有较低的RCS.      

基于自适应迭代的机器人曲面恒力跟踪

针对利用机器人进行打磨、抛光、去毛刺等场合时末端执行器对曲面工件轮廓跟踪时难以得到恒定接触力的问题,对机器人末端执行器和工件轮廓接触时的接触力进行研究,建立了实际跟踪过程中机器人末端执行器的接触力和已知传感器坐标系的映射关系,提出了一种基于自适应迭代学习算法的机器人力/位混合曲面恒力跟踪控制方法。该方法由两部分组成:基于机器人和环境接触时的阻抗模型设计了迭代学习控制律,在PD反馈控制的基础上通过迭代项克服机器人的未知参数和不确定性,并构建Lyapunov能量函数证明所提控制律的收敛性;将迭代学习控制律和力/位混合曲面恒力跟踪控制方法结合起来设计了用于曲面工件轮廓跟踪的控制方法。实验结果显示,经过15次迭代,接触力的波动范围逐渐变小并稳定在±3 N之内,验证了所提方法的有效性。      

轴对称推力矢量控制伺服机构空间运动解耦研究

针对由发动机摆动喷管矢量控制的两个伺服机构间的牵连效应导致的伺服机构运动耦合问题,提出了一种新型空间解耦计算方法。该方法利用欧拉角描述喷管姿态,通过空间齐次坐标变换矩阵原理对两路喷管伺服机构进行解耦计算,得到了两路伺服机构的伸长量解析解,确保了控制指令的精确,并通过ADAMS虚拟样机仿真验证了该计算方法的正确性。仿真实验结果表明,此计算方法与虚拟样机的仿真结果完全吻合,计算精度高,可以对喷管任意姿态进行运动精确控制。     

空间充气薄膜衍射成像结构的热力耦合分析

随着空间光学技术的发展,衍射成像系统成为空间光学领域的研究热点。提出了一种充气薄膜衍射成像结构设计方案,针对其在空间热载荷作用下的热力耦合效应开展研究。建立了充气薄膜衍射成像结构的热力耦合顺序分析框架,基于I DEAS和ANSYS软件二次开发技术,实现了热力耦合数值分析,得到了结构在轨运行时的瞬态温度场、热应力、热变形等。分析结果表明,充气支撑结构产生了较大热变形,但最大热应力远小于聚酰亚胺薄膜材料的抗拉强度,而镜面中心处没有出现较大的热变形响应,为充气薄膜衍射成像结构的防热设计提供了技术支持。     

热力耦合问题数学均匀化方法的物理意义

针对复合材料周期结构热力耦合问题,通过构造各阶摄动项的全解耦格式,推导了高阶数学均匀化方法（MHM）的数学表达式,并使用加权残量方法将其转换为易于编程实现的矩阵列式。将弹性影响函数和热影响函数分别比拟为弹性虚拟位移和热虚拟位移,通过弹性虚拟载荷和热虚拟载荷的自平衡特性、量纲分析及几何直观等角度揭示了各阶影响函数和摄动位移的物理意义,并指出二阶摄动位移对于细观结构分析的必要性。数值计算结果验证了高阶MHM矩阵列式及物理意义分析的正确性。      

燃气源控制气流非稳态耦合传热问题研究(Ⅰ)--计算模型和方法

针对某燃气源控制气流与导热套环之间的非稳态耦合传热问题,建立了数学模型,并采用整体求解的方法进行求解,即把流体区域的流动和对流换热与固体区域的导热过程组合起来,作为一个统一换热过程来处理。为确保计算结果在物理上的真实性,采用“假密度”法求解以温度T为求解变量的能量方程。对于非稳态控制方程离散后的离散方程采用了PISO方法进行求解。     

南北位保后东西耦合修正量计算

地球同步三轴稳定卫星在南北保持控制过程中,存在较大的东西方向耦舍量。实践中,在南北控制结束后,直接进行东西向轨道修正。给出了卫星寿命初期和中期东西修正量的计算模型,并在实际控制中取得了良好的效果。     

推力协同的大柔性飞行器纵向姿态控制

为了研究大柔性飞行器飞行/结构耦合动力学特性,提出了改进的面向控制的大柔性飞行器多体模型,开展了大柔性飞行器纵向动力学耦合特性分析与推力协同下纵向姿态控制律设计。采用二面角动态近似描述大柔性飞行器结构动力学特征,并推导了纵向耦合动力学模型。根据改进模型在配平点处的线性化模型,分析了飞行/结构耦合系统的纵向稳定性与结构变形量之间的关系。针对大柔性飞行器姿态稳定与跟踪,设计了纵向姿态控制器。与常规飞行器相比,大柔性飞行器飞行过程中会发生大变形,当载荷较大时,配平构型近似“U”形,此时纵向动力学具有长周期不稳定的特征。分析结果表明:大柔性飞行器各模态之间的耦合程度随着变形的增大而增大。此外,纵向姿态控制需要升降舵与推力协同控制速度和俯仰角并且考虑结构动力学的影响,否则飞行/结构的耦合作用会导致姿态跟踪误差衰减缓慢甚至发散。