一种跨声速定常流场求解加速方法

跨声速定常流场的隐式求解相当于使用牛顿迭代法求解一个非线性方程组。为满足牛顿迭代收敛性的要求,通常需要对所求解问题进行全局化处理。在同伦延拓的框架内,提出了一种基于拉普拉斯算子的方程延拓方法,提高了定常流场隐式求解收敛速度。针对定常流场通常初始化为均匀来流的特点,一方面利用拉普拉斯算子的椭圆性加快边界条件信息向流场内部的传播,另一方面利用拉普拉斯算子的线性和正定性改善延拓问题的正则性,综合两者增加拟牛顿算法的稳定性,提高可用CFL数,最终达到提高流场求解效率的目的。由于流场问题的复杂性和非线性,难以通过理论分析得出先验的最优非线性求解策略。因此,通过无黏NACA0012翼型、湍流RAE2822翼型和三维ONERA M6机翼等算例的数值实验,研究了拉普拉斯项参数对收敛效率的影响,给出了效率较优的参数组合,验证了本文方法在跨声速情况下相对于经典伪时间推进法可以节约20%以上的CPU计算时间。      

多刚体系统动力学的正则方程与辛算法

针对多刚体系统动力学数值计算精度的关键问题，首先采用Ｋａｎｅ方程导出了二阶形式的动力学方程，并研究了其中系数矩阵的恒等关系，进而建立了多刚体系统动力学的哈密顿体系并获得正则方程。在此基础上，分析了多刚体系统动力学数值积分的能量耗散现象，揭示其根源在于，数值算法用前一时刻值代替当前值。理论分析和计算结果表明，采用多刚体系统动力学的正则方程与辛算法的结合，可使多体系统动力学的数值计算性态得以改善，能有效消除“能量耗散”现象，从而使动力学计算的误差累积与违约得以控制，保证了计算精度。     

串联型MPPT空间电源系统稳定性技术

串联型最大功率点跟踪(MPPT)空间电源系统存在母线电压振荡或跌落等不稳定现象,文章以MPPT不调节母线电源拓扑为研究对象,建立了系统等效电路模型,通过小信号等效分析、状态空间方程特征根求解等手段,分别对MPPT控制域及电压控制域下电源系统的稳定性进行了分析.结果表明,在太阳电池阵及阵列功率调节器特性曲线一定的条件下,通过对二者间电感、电容参数合理配置,可以保证系统稳定.在MATLAB/SIMULINK中搭建了串联型MPPT不调节母线电源系统仿真模型并进行实例仿真,仿真结果验证了理论分析的正确性.     

摄动条件下椭圆轨道卫星相对运动研究

基于严格定义的相对轨道要素,考虑非球形引力、大气阻力和三体引力等摄动力的影响,用拟平均根数法添加相关摄动项,推导了椭圆轨道摄动条件下无奇点的卫星近距离相对运动模型。分析了考虑摄动力的影响时相对轨道要素的特性,并据此分析了各摄动力对卫星相对运动轨道预报精度的影响,比较了有无摄动的相对运动模型与精确解间的误差。算例验证了方法与结果的正确性和有效性。研究结果完善了相对轨道要素的基本理论,并将其适用范围拓展至摄动条件下的椭圆基准轨道卫星相对运动研究领域。     

基于量纲分析的弹塑性应力强度因子探讨

以Ⅲ型裂纹(反平面剪切)问题为例,对裂纹尖端弹塑性应力场进行了分析研究。首先,基于断裂力学基本理论,引入弹塑性本构关系,推导弹塑性裂纹尖端场;其次,依据线弹性场中应力强度因子的定义和基于J-积分的弹塑性裂纹尖端应力场的分布,并结合物理意义,定义了与线弹性场应力强度因子一致的弹塑性应力强度因子,来表征弹塑性裂纹尖端应力场的强度;最后,对本文所定义的弹塑性应力强度因子进行了合理性分析,说明了本文定义的适用性。本文所定义的弹塑性应力强度因子,相比于表征弹塑性场的J-积分,具有更明确的物理意义和更简明的表达形式;为弹塑性裂纹尖端场的分析和含裂纹材料弹塑性破坏判据的建立提供了理论方法。另外,本文方式也适合于定义Ⅰ型和Ⅱ型问题。     

复杂外形网格生成技术及应用

以某型双立尾战斗机为例,给出了实际复杂型号飞行器结构网格生成的一种分块方案及实现过程,用一种矢性三次多项式插值的方法生成相邻块公共交界面的网格,用Hilgenstock方法修正源项,求解椭圆型方程生成块内空间网格,使网格线与边界面夹角、与边界面间距控制为期望的值.在所生成的双立尾战斗机分块结构网格中求解Euler方程,获得了合理的全机气动力系数、翼面压力分布和绕流图谱结果.     

非结构三角形网格的一个ENO型有限体积方法

基于双曲型守恒律方程,对非结构三角形网格给出了一种ENO(Essentially Nonoscillatory Scheme)型有限体积格式,方法的主要思想是先对每一个三角形单元构造一个加权的二次插值多项式,而在计算交界面的流通量时采用了两点高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法.最后给出了该格式收敛的数值阶,并对前台阶问题进行了计算.      

强迫振动法求解偏航阻尼导数

求解了非定常薄层近似Ｎ－Ｓ方程模拟强迫振动条件下的非定常扰动流场,并在Ｅｔｋｉｎ定义下给出偏航阻尼导数Ｃηγ－Ｃηβｃｏｓα的计算方法和数值积分公式,还给出Ｍａ＝１０,α＝９．８４７°,１９．６８°下的计算结果。     

HTPB推进剂宽泛应变率下粘弹性本构模型研究

为研究HTPB推进剂的率相关性力学性能,采用材料万能试验机、液压试验机和分离式霍普金森压杆(SHPB),分别开展了低(1.67×10~(-4)~1.67×10~(-1)s~(-1))、中(1~100 s~(-1))、高(700~2 500 s~(-1))应变率的单轴压缩实验。实验结果表明,HTPB推进剂的压缩力学性能是率相关性,随应变率的升高,给定应变下的应力逐渐增大。采用广义非线性ZWT本构模型描述HTPB推进剂宽泛应变率下的压缩力学行为,模型预测与实验数据对比表明,模型中至少需要4个麦克斯韦元件。     

CRM翼身组合体模型高阶精度数值模拟

基于五阶空间离散精度的WCNS格式,开展了CRM翼身组合体模型的高阶精度数值模拟,以评估WCNS格式对复杂外形的模拟能力以及典型运输机巡航构型阻力预测的精度。首先依照DPW组委会提出的网格生成指导原则,利用ICEM软件生成了粗、中、细、极细四套网格,网格规模从"粗网格"的2 578 687个网格点逐渐扩展到"极细网格"的65 464 511个网格点。研究了设计升力系数下,网格规模对气动特性、压力分布和翼根后缘局部分离区的影响,采用"中等网格"开展了抖振特性的数值模拟研究。通过与二阶精度的计算结果、DPW V统计结果和部分试验结果的对比分析,高阶精度数值模拟结果表明,阻力系数计算结果与DPW V统计平均结果吻合较好;网格密度对机翼上表面的激波位置和翼身结合部后缘局部分离区略有影响;迎角为4°时,升力系数下降的主要原因是机翼上表面激波诱导分离区和翼身结合部后缘局部分离区的增加。