二元矢量喷管流场的N－S方程计算

对二元矢量喷管在不同工况下的流场进行了数值模拟。求解的控制方程为三维非定常Ｎ－Ｓ方程。空间方向用有限体积离散，时间方向用Ｒｕｎｇｅ－Ｋｕｔｔａ法推进。粘性项的计算上采用了一种有限体积框架下新的离散化方法。湍流模型采用了广泛使用的Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数模型。计算表明，本文方法能准确地计算二元矢量喷管、非矢量喷管在设计与大设计落压比下的流场。     

有限元的一种并行算法及实现

当前,由于有限元素法大都使用串行算法,处理时间较长。本文根据有限元素法的特点和采用的并行计算机结构,提出了一种并行有限元算法,能大大提高其处理速度,加快处理周期。且算法简单,实现方便。     

非线性机翼极限环颤振的研究

在定常空气动力的作用下,对含立方非线性刚度的二元机翼颤振系统的极限环颤振进行了研究.应用中心流形理论将原四维系统降为二维,采用后继函数法对分岔点类别进行了定性的分析,从而确定平衡点的性质,并应用范式理论对分岔点处中心流形约化方程进行化简,进而研究了系统参数对极限环颤振的稳定性以及幅值的影响.     

钝体边界层失稳特性分析

推导了适合可压缩二维/轴对称边界层的线性抛物化稳定性方程。在法向采用4阶中心格式,流向采用1阶向后差分,对线性抛物化稳定性方程进行离散。通过来流马赫数为4.5的平板与曲面边界层稳定性问题验证了程序的正确性。讨论了来流马赫数6的球钝锥边界层内,扰动波沿端头向下游演化的失稳特性。发现了一种新的失稳机制。     

单列叶栅紊流流场有限分析数值研究

本文将有限分析方法用于曲线座标系上紊流 N-S方程数值计算,研究了高雷诺数时叶栅粘性紊流流场。有限分析方法在网格单元上对非线性偏微分程进行线化处理,在解析边界条件下求出线化偏微分方程的解析解,以此解构造离散方程。有限分析方法能够根据对流的方向和大小自动改变格式系数,并具有数值扩散小、精度高和稳定性好等优点。本文以 k-ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。      

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

考虑预应力和液体影响的圆柱壳颤振分析

 为了研究预应力和液体对圆柱壳颤振的影响,采用杂交元方法,建立了圆柱壳的气动弹性方程。由Sanders薄壳理论和经典有限元理论,从壳的精确解推导节点位移函数入手,并由一阶活塞理论得到结构的气动阻尼矩阵和刚度矩阵,最终推导出考虑预应力和内部液体影响的圆柱壳的混合有限元公式和气动弹性方程。通过特征值法验证了此种有限元方法的正确性,并重点研究了预应力和内部注有液体对圆柱壳稳定性的影响。计算结果表明:预应力和液体对圆柱壳的颤振特性有显著的影响。     

导电圆柱壳体磁弹性二阶谐波共振响应分析

 研究了磁场环境中受机械载荷作用圆柱壳体的二阶谐波共振响应问题。首先给出了电磁场环境中导电圆柱壳体的非线性磁弹性振动方程和电磁场方程,并应用伽辽金积分法导出了相应的无量纲化非线性振动微分方程。然后采用多尺度法对系统的二阶超谐波和亚谐波共振问题进行求解,得到了稳态运动下关于定常解的幅值响应方程及解的稳定性判别式。最后通过数值算例,给出了共振幅值随磁感应强度和外激励力等参量变化的特征曲线图,分析了解的稳定性、奇点性态及分岔特性。结果表明,通过电磁或机械参数的适当选取,能够达到激发或抑制系统共振现象的目的。     

Stress–strain analysis of a thin-shell part of fuselage using a triangular finite element with lagrange multipliers

We developed a triangular discretization element of the fuselage fragment in the form of a thin shell compatible with the adjacent elements based on Lagrange multipliers. Solution of several test problems showed a high efficiency of using the modification of a triangular finite element being proposed.     

二维涡格法的超收敛性

本文提出并研究了二维涡格法的超收敛性。在离散化误差分析中,与薄翼理论的二维平板精确解进行比较,根据Chebychev多项式理论,导出了涡格法中具有超收敛性的离散化格式。并证明了这种离散化格式对二维抛物弧板和三次曲线弧板等问题也具有超收敛性。     

一种计算非定常二维流动的无网格算法

主要目的是发展一套求解非定常流动的无网格算法。计算区域的离散方面,提出了一种按区域进行填充布点的点云自动生成方法;发展了一种点云的运动技术来实现离散点云对物面边界的随体运动;在点云离散的基础上,采用最小二乘法求解矛盾方程的方法来求取空间导数,进而获得数值通量;采用双时间方法进行时间离散推进,其中物理时间迭代采用二阶隐式格式,伪时间迭代采用四步龙格一库塔显式格式,为了加速收敛,在伪时间迭代中采用了当地时间步长和隐式残值光顺等加速收敛措施。最后,利用本文算法模拟了NACA0012翼型和NACA64A010翼型的跨音速非定常流动,并将计算结果与实验测量结果进行了对比分析,验证了上述方法的正确性和实用性。     

用N－S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二级时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型￣［１］。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

极坐标系下谱元方法求解不可压缩Navier-Stokes方程

本文利用有限元的思想并结合谱方法的精度提出求解偏微分方程的谱元方法,并将谱元方法应用到极坐标系中;详细推导了在极坐标系下的谱元方法的具体计算公式,求解了极坐标系下的简单椭圆型二阶偏微分方程;并结合时间分裂方法应用于同心旋转圆筒间流体的流动问题,具体求解了原始变量速度和压力的不可压缩NavierStokes方程,均取得了满意的结果.     

TIME ACCURATE 2D NAVIER-STOKES CALCULATIONS WITH A DUAL-TIME STEPPING METHOD

ＴＩＭＥＡＣＣＵＲＡＴＥ２ＤＮＡＶＩＥＲＳＴＯＫＥＳＣＡＬＣＵＬＡＴＩＯＮＳＷＩＴＨＡＤＵＡＬＴＩＭＥＳＴＥＰＰＩＮＧＭＥＴＨＯＤＹａｎｇＹｏｎｇ，ＱｉａｏＺｈｉｄｅ（Ｆａｃｕｌｔｙ５０３，ＮｏｒｔｈｗｅｓｔｅｒｎＰｏｌｙｔｅｃｈｎｉｃａｌＵｎｉ...     

基于QSS的自适应多步校正算法及其在柔性航天器动力学中的应用

量化状态系统（QSS）算法是一种基于状态变量离散化的数值积分方法,该方法与基于时间离散的传统方法显著不同,QSS通过计算状态变量每次跃迁所需的时间来推进下一步的积分。在求解非刚性常微分方程时,QSS算法比传统算法更具优势,但它不适合求解刚性问题。为此提出一种基于QSS的自适应多步校正算法（AMCQSS）,该算法以QSS为基础、结合隐式多步法思想,在计算过程中可以自适应选择二步法或三步法,以有效提高求解刚性问题的精度和效率。通过对柔性航天器动力学的仿真求解,验证了算法的可行性。将该算法与ODE23tb、ODE15s、ODE45以及QSS等算法进行对比,结果表明AMCQSS算法既能保证求解的效率及精度,又具有较好的收敛性和稳定性。     

扫描有限元法求解二元喷管的三维可压粘性流场

发展了一种二元喷管流场计算的三维Ｎ－Ｓ方程扫描有限元方法,给出了Ｎ－Ｓ方程中二阶导数项新的离散方法,以及更完善、简洁的扫描有限元方程的表达式。开发了相应的计算程序,用所发展的方法对二维平面收扩喷管、两种大宽高比的二元喷管及其圆变方转接段的三维粘性流场进行了计算,计算结果表明本方法具有较宽广的适用范围,计算结果良好     

三维复杂区域对流换热问题的控制方程及计算

研究并发展了一套计算复杂区域对流换热问题的通用方法并开发了相应的可视化计算程序。针对所建立的一组以直角坐标速度分量表示的非正交曲线坐标系下的三维流动及传热控制方程,采用有限体积方法离散,SIMPLER算法求解,并利用网格界面速度插值方法解决采用非交错网格时出现的压力波动问题。本文利用所开发的程序计算了低散热发动机冷却通道自然对流换热问题。      

Analysis of nonlinear deformation and stability of oval cylindrical shells under combined loading by bending and torsion moments

A problem of geometrically nonlinear deformation and stability of cylindrical shells with noncircular cross-section contour is solved by the finite element method. The nonlinear system of algebraic equations to determine the nodal unknowns of finite elements is obtained with the use of the variational Lagrange principle. The system is solved by the method of successive loadings with the use of the Newton-Kantorovich method of linearization, the method of solving linear Craut’s equations and Sylvestor’s stability criterion. The nonlinear deformation and stability of shells with the oval and elliptical cross-sections under combined loading by bending and torsion moments are analyzed. The critical loads and buckling modes of shells are determined. The influence of deformation nonlinearity, shell ovality and ellipticity on the critical loads is clarified.     

直升机旋翼/机身耦合系统气动/机械稳定性分析

考虑了铰接式旋翼桨叶绕挥舞、摆振和变距铰的整体刚性运动与桨叶中等弹性变形之间的动力学耦合作用,将直升机机身和传动轴作为弹性体。通过构造一种特殊的24自由度刚柔混合单元得到旋翼/机身耦合系统的周期时变动力学方程,根据Floquet理论对稳态周期解的稳定性进行研究。计算结果表明,旋翼轴的截面刚度对旋翼/机身耦合系统的气动/机械稳定性有一定的影响。      

控制面间隙对非线性二元机翼气动弹性响应的影响

采用Theodorsen非定常气动力建立同时含有俯仰立方非线性和控制面间隙非线性二元机翼的运动方程,并以状态空间形式描述。基于状态依赖Riccati方程控制方法,设计了非线性颤振控制律。综合运用Runge-Kutta数值方法与Henon方法,研究了控制面间隙对系统开环/闭环响应的影响。其中Henon方法用以准确快速地确定间隙非线性的转折点,从而可以避免间隙非线性引起的数值不稳定现象。仿真结果显示,俯仰立方非线性可以使间隙非线性系统产生振幅稳定的极限环振荡;控制面间隙对开环响应的影响随着来流速度的升高而减弱;在速度较高的情况下,控制面间隙导致闭环系统响应产生过阻尼现象。