跨声速颤振计算方法研究

采用改进的弹簧近似方法,以ALE(Arbitrary Lagrangian-Eulerian)描述的Euler方程为控制方程,耦合结构运动方程,计算了Isogai机翼模型的颤振边界。气动力求解采用双时间推进,外时间为物理时间,与结构运动方程同步;结构运动方程求解采用Rayleigh-Ritz方法,时间推进也采用双时间法,节省了大量计算时间。研究表明,计算结果与参考文献提供的结果吻合良好,证明了所提出的跨声速颤振计算方法的准确性。     

低速冲击与准静态压痕力下复合材料层合板的损伤等效性

低速冲击与集中准静态压痕(QSI)力对复合材料层合板所造成的损伤具有等效性。通过两类损伤方式对比试验,获得了冲击能量(或准静态压痕力)与层合板损伤面积、损伤宽度和凹坑深度3组对应关系。3组关系的对比分析表明,凹坑深度是表征损伤的最佳参数,因此将凹坑深度作为损伤参数来建立落锤冲击损伤与准静态压痕力损伤间的等效关系。发现拐点是层合板纤维和基体整体抵抗冲击(或准静态压痕力)能力的最大值,两类试验的拐点相差很小,且两类试验变化趋势相同,说明了用集中准静态压痕试验取代落锤低速冲击试验是可行的。对相同凹坑深度下的冲击能量与准静态压痕力进行了数值分析,并提出了冲击能量与准静态压痕力间的等效性公式。研究表明,文中建立的等效性公式能够很好地反映落锤冲击能量与准静态压痕力之间的关系,可为今后由集中准静态压痕试验代替落锤冲击试验提供理论支持。     

复杂前机身—进气道组合体的流场计算

本文提出了一种能够对复杂前机身-进气道组合体进行流场计算与气动分析的方法。通过求解Euler方程来模拟前机身-进气道组合体在亚、跨、超声速,不同攻角和发动机流量条件下的各类流场。     

全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

横向喷流与超音速主流干扰流场的数值研究

 采用欧拉方程与可压缩涡对模型相结合的数值方法,研究了超音速主流与横向喷流干扰流场的气动特性。喷流横截面用一对反向旋转的可压缩涡来模拟。采用有限元Osher格式及非结构网格数值求解二维欧拉方程以得到非等熵外流场的气动特性及弓波位置。将欧拉方程数值解结果与可压缩涡对流场解结果相结合以研究喷流的轨迹及喷流横截面的变化。最后给出了数值模拟结果并与工程估算结果进行了比较,表明本文所述方法是可行的。     

基于Riemann问题的解在自适应无结构网格上求解Euler方程

给出了一种根据 Riemann问题的解 ,计算网格单元边界处的守恒量通量的方法。该方法不受网格单元形状的限制 ,具有较好的通用性。实际计算表明 ,采用这一方法在自适应无结构网格上求解 Euler方程 ,对于复杂的流动细节具有较高的分辨率 ,并且使用多重网格法能够加速收敛     

用运动嵌套网格方法数值模拟旋翼前飞非定常流场

通过求解Ｅｕｌｅｒ方程数值模拟了直升机旋翼前飞非定常流场。为了模拟包括旋转、周期性变距和周期性挥舞的非定常运动,本文采用了一种能够快速完成重叠网格间流场信息交换的运动嵌套网格方法。空间采用中心平均的有限体积法进行离散,时间方向采用双时间法隐式推进求解。为了验证方法的正确性,数值模拟了悬停流场,数值结果和实验值吻合很好。尽管前飞状态的计算没有实验验证,但计算结果与其它文献计算结果吻合很好。     

高超声速化学非平衡欧拉方程数值模拟

本文给出了求解高超声速化学非平衡无粘流场的数值模拟方法。计算中采用守恒变量型的ＮＮＤ差分格式和全流场捕捉技术;建立了（数学）奇性轴上的控制方程,克服了通常在其上插值方法所引起的扰动,改善了收敛性和计算精度;化学反应模型由五组元五反应组成;温度和压力的计算是非迭代的,有别于国内外惯用的迭代算法;此外,还发展了背风过低压区的处理方法,从而提高了欧拉方程的应用范围,使无粘非平衡绕流的计算攻角达到３０°     

二维轴流压气机转—静非定常Euler方程数值计算

将Euler方程写在包括转子和静子的绝对坐标系下,通过在转—静之间设计一个可以随时间扭动的平行四边形网格计算了二维轴流压气机转—静非定常流场。由于整个方程统一写在动网格下,转—静之间没有使用插值。      

三维机翼气动结构多学科优化方法

采用Euler方程、连续伴随优化方法、自由形面变形技术（free-form deformation,FFD）以及一种四边形线性壳单元模型对三维机翼气动结构多学科优化方法开展了研究。目标函数为阻力系数和机翼翼盒结构质量的加权和,设计变量为FFD控制点以及10段蒙皮的厚度。气动方面升阻力系数对设计变量的梯度由伴随方法求得,结构方面翼盒质量对于蒙皮厚度的梯度可直接计算,应力的梯度采用有限差分方法获得。对跨声速大展弦比机翼进行气动结构多学科优化,并将该计算结果与机翼气动单学科优化设计后的静气弹计算结果进行对比,结果表明:该气动结构多学科优化方法能够在满足约束条件下,实现阻力系数和翼盒质量同时降低,保证优化结果的合理性。