轮盘振动特性计算

本文提出了一种计算轮盘振动特性的数值解法。采用薄板和厚板两种理论,将问题化成一阶微分方程组进行数值计算。编制了计算机程序,得出了令人满意的结果。由于没有涉及到复杂的数学问题,计算机程序大为简化,计算时间少。      

压气机轮盘孔边应力的边界元分析

 本文根据压气机轮盘的轴对称问题,取出轮盘包含孔的一个扇形部分,以有限元法计算的结果作为扇形部分的边界条件,采用边界元法计算轮盘孔边的应力。这与直接用有限元法计算结果对比表明,边界元法具有数据准备工作量小、方便,计算点选择灵活,计算精度高等优点,用于计算孔边的应力集中问题是有效的。文中还推导丁常单元的二阶差分逼近式,用以求解边界上的应力。     

螺旋桨滑流对带后缘襟翼机翼气动特性影响的数值分析

采用动量理论结合有限体积法求解N-S方程的方法,进行了螺旋桨滑流对带后缘襟翼机翼气动特性干扰的数值模拟分析.螺旋桨模型简化为一个激励盘,即一个无厚度可穿透的圆盘,以给定的圆盘拉力载荷,即压强差来模拟螺旋桨滑流的效果.采用嵌套网格技术处理桨盘与机翼间复杂的几何关系.算例计算表明,该计算方法能够较好的模拟滑流对带后缘襟翼机翼气动特性的影响.     

任意坐标系中的弹塑性矩阵

 本文推导了三维问题及平面应力、平面应变、轴对称问题在任意曲线坐标系中的弹塑性矩阵。由此最一般的表达式可以方便地得到任何特定坐标系中的弹塑性矩阵,供各类不同特点问题的有限元计算时使用。 文中推导得到的任意曲线坐标系中的本构方程,也同时为曲线差分法求解弹塑性问题提供了刚度系数显式。     

处于非均匀温度场中含裂纹变厚度旋转盘的应力强度因子广义J积分解法

本文首先从直法线假设出发导出了处于非均匀温度场中变厚度各向同性钣元素本构方程并证明该方程与等厚度正交各向异性钣元素的本构方程相同。然后本文从上述本构方程出发,引用自由能与熵函数以及固体力学所有基本方程导出了与积分路径无关的广义J积分J_8,该广义J积分可用于处在非均匀温度场中的含裂纹变厚度旋转盘。最后本文应用这个广义J积分与有限元素法计算了涡轮盘应力强度.因子K_1并得到了令人满意的结果。     

错频转子系统动力特性研究

根据某级压气机叶／盘转子的实际结构，建立了叶／盘带短轴的有限元模型。采用NASTRAN结构分析软件，得到了模化协调转子所需的数据，构建了多自由度叶／盘转子系统动力学模型，并在此基础上通过改变多自由度模型中叶片刚度实现错频。采用基尔法求解错频转子系统的运动方程。获得了此多自由度系统的动态响应，并对不同错频方案进行了对比分析。     

用虚拟黏性法构造Navier-Stokes方程黏性项的隐式格式

提出的虚拟黏性法是一种关于Navier-Stokes方程黏性项(黏性应力项和传热项)的隐式新方法。通过引入虚拟时间和虚拟黏性项,将隐式格式的构造大大简化,从而避免了大型复杂隐式差分方程组的常规求解。在虚拟时间推进过程中所需求解的方程组的系数矩阵是一个三对角矩阵,它具有计算简单且计算量小的优点。用模型方程和Navier-Stokes方程进行了数值仿真,研究各种参数对计算的影响,并在精度和效率上与显式方法进行比较,证实了算法的正确性和优势。最后分析了本文格式的适用范围。     

电弧喷射推力器化学非平衡数值模拟

为了准确把握电弧喷射推力器工作过程物理机理与特征,采用化学非平衡模型对不同压缩室直径下低功率氮电弧喷射推力器工作情况进行了数值模拟。模型考虑了工作过程中的分解反应、电离反应和复合反应,化学动力学模型为4组分,4反应的有限速率化学反应模型。采用二阶精度NND格式数值求解耦合电磁源项和化学反应源项的N S方程组,采用有限控制容积积分方法求解由麦克斯韦方程组推导出的电磁场方程。数值模拟的结果揭示了推力器内部电弧能量输入作用和高温电离气体的离解电离状况,给出了不同压缩室直径下推力器的推力、比冲和推进效率。结果分析表明,压缩室直径对推力器性能具有较大影响。     

计算机翼跨声速压力分布的格林函数方法

本文给出一种用格林定理将跨声速小扰动速势方程化成积分方程式,然后进行数值求解的方法,文中对场元强度进行松弛迭代,一般只需20次左右即达收敛。算例表明,该方法结果与差分法结果一致。     

螺旋桨滑流与机翼气动干扰数值模拟研究

分别采用非定常方法和激励盘理论研究螺旋桨滑流对机翼的干扰作用,并对计算结果进行对比分析。非定常方法基于动态面搭接结构化网格,运用有限体积法求解非定常RANS方程。对于激励盘理论,首先采用点对点结构化网格,随单独螺旋桨旋转的方式进行非定常模拟得到螺旋桨桨叶压力分布,然后将桨叶压强分布转化为激励盘载荷分布,最后将载荷分布作为边界条件对干扰流场进行定常模拟。数值结果表明:除了在较大来流攻角时有所差异,两种方式得到的机翼气动力吻合良好;非定常方法可以得到螺旋桨与机翼的瞬时干扰流场,而激励盘理论则可以大大缩短数值模拟的时间并降低计算要求。     

轻负荷磁头的动态特性

本文采用有限差分法求解非定常修正Ｒｅｙｎｏｌｄｓ方程，采用四阶Ｒｕｎｇｅ－Ｋｕｔｔａ法求发磁头滑块的动力方程组，分别计算出磁头滑块受到冲击或飞越各种障碍物时的动态响应，较全面地反映了磁头滑块的动态时域特性。     

跟踪问题最优控制律精细积分

构成有限时间最优跟踪系统的控制律需要求解 Riccati微分方程及外部控制输入向量满足的微分方程 ,前者是非线性矩阵微分方程 ,后者是变系数线性微分方程。在结构力学与最优控制的模拟理论基础上所发展的精细积分方法借鉴了计算结构力学中的算法 ,可以精确有效地求解这些微分方程。这种方法的特点之一在于步长幅度变化较大时 ,Riccati微分方程的数值解仍可以保持很高的精度 ,并且变系数线性微分方程的求解亦可纳入其体系而不必用通常的差分方法。本文介绍了用精细积分方法求解这些方程的过程 ,并给出了数值算例。     

Small Signal Compensation of Magnetic Fields Resulting from Aircraft Maneuvers

Tolles and Lawson identified three permanent, five induced, and eight eddy-current fields as sources of magnetic interference associated with airframe maneuvers. Small signal approximations are used here to separate the eddy-current terms and thus decouple the sixteen equations into two sets of eight equations. It was found that a singularity exists in the small signal equations for the permanent and induced terms. This causes an ambiguity amoung three of the coefficients which can be resolved mathematically by resorting to large maneuvers. Flight test data exhibit a large amount of magnetic hysteresis and the magnetic anomality detector (MAD) equipped aircraft will not remain compensated from one flight to the next. This complicates the problem of resolving the dip angle ambiguity. Furthermore, significant differences exist in the compensation terms as determined from manual pilot maneuvers and the terms when the maneuvers are performed by the automatic pilot. This has been attributed to differences in altitude stability and the frequency content of the maneuvers. It was found that both optimal frequency filtering and altitude compensation can be used to improve the figure of merit (FOM) resulting from pilot maneuvers.     

具有随机参数的周期性结构振动分析的谱方法

为了探讨整圈连接的涡轮叶片组内各个叶片的制造偏差对叶片组振动的影响,本文利用一个具有周期性随机参数的结构模型来近似地模拟叶片组结构,并提出了分析这种结构的振动的一种谱方法。假定结构参数的标准差为微小量,因而可应用摄动法。将周期性随机结构参数展开成付氏级数,从而求解结构的自由振动和受迫振动,求得结构的频率和振型,以及共振振幅及其方差估计。论证了主振型的正交性,分析了此结构共振的特殊条件。算例表明,分析结果与实验结果具有相同的量级。     

一种基于叶盘结构几何失谐的降阶分析方法

叶片几何失谐严重影响叶盘结构的振动特性。为了准确且高效地对存在叶片几何失谐的叶盘结构进行振动分析,基于扇区建模方法,由动力学方程推导出降阶模型的表达式,利用该降阶方法对整个几何失谐叶盘结构分别进行模态分析和受迫响应分析,并与有限元方法进行固有频率和受迫响应结果的对比。结果表明:该降阶方法得到的结果与有限元方法得到的结果基本一致,验证了该方法的有效性和实用性。     

绕升力机身跨音速流动的有限元素法数值分析

 <正> 本文应用全速位方程最小压强积分的有限元素法解绕升力机身的跨音速流动。机身头部可以是尖的或带有进气道的。用头部伸出无穷长圆柱来模拟头部进气。尖头机身在元素形状处理上具有一定的复杂性,从数值实验角度考虑我们分别按尖头和头部伸出无穷长细圆柱来近似尖头头部两种方法对尖头机身进行计算。采用人工密度法捕获激波,但由于绕机身流线形状比较复杂,直接采用机翼问题中的人工粘性公式,解往往不收敛。我们推导出一种较精确的,适用于复杂流动情况的人工粘性计算公式,为了加速收敛,应用网格逐次加密技术,可使收敛速度提高2～3倍。本文还给出绕升力机身跨音速流动的远场速位解析式。     

基于CFD的电磁散射数值模拟

计算流体力学(CFD)技术中的时域有限体积法(FVTD)被推广到计算电磁学中,直接数值求解时变麦克斯韦方程组.FVTD使用电磁场矢量共置于网格单元中心、散射场积分形式的麦克斯韦方程组守恒形,空间离散使用基于特征值的近似黎曼解构建网格单元边界通量,时间推进采用四阶龙格-库塔法.TM和TE波极化下几种二维完全导电体的表面诱导电流密度和雷达散射截面(RCS)计算验证表明,时域有限体积法是一种高精度有效的时域方法.     

Motion equations of hemispherical resonator and analysis of frequency split caused by slight mass non-uniformity

The mass non-uniformity of hemispherical resonator is one of reasons for frequency split, and frequency split can cause gyroscope to drift. Therefore, it is of great significance to analyze the relationship between mass non-uniformity and frequency split, which can provide a theoretical basis for mass balance of imperfect resonator. The starting point of error mechanism analysis for gyroscope is the motion equations of resonator. Firstly, based on the Kirchhoff-Love hypothesis in the elastic thin shell theory, the geometric deformation equations of resonator are deduced. Secondly, the deformation energy equation of resonator is derived according to the vibration mode and relationship between the stress and strain of hemispherical thin shell. Thirdly, the kinetic energy equation of resonator is deduced by the Coriolis theorem. Finally, the motion equations of resonator are established by the Lagrange mechanics principle. The theoretical values of precession factor and natural frequency are calculated by the motion equations, which are substantially consistent with the ones by the finite element method and practical measurement, the errors are within a reasonable range. Simultaneously, the varying trend of natural frequency with respect to the geometrical and physical parameters of resonator by the motion equations is consistent with that by the finite element analysis. The above conclusions prove the correctness and rationality of motion equations. Similarly, the motion equations of resonator with mass non-uniformity are established by the same modeling method in case of ignoring the input angular rate and damping, and the state equations with respect to the velocity and displacement of vibration system are derived, then two natural frequencies are solved by the characteristic equation. It is concluded that one of reasons for frequency split is the 4th harmonic of mass non-uniformity, and thus much attention should be paid to minimizing the 4th harmonic of mass non-uniformity in the course of mass balancing for imperfect resonator.     

电弧喷射推力器电磁场模拟

为了揭示电弧喷射推力器内部电磁场特征及其与高温电离气体间的相互作用机理,建立了比较完善的电磁场模型,详细讨论了电磁场数值模拟方法及其计算效果.模型同时考虑了电场力、洛仑兹力和电子压强梯度对传导电流的贡献,由电子动力学方程和粒子微观碰撞理论推导出电流方程,由麦克斯韦方程组推导出电磁场控制方程,采用包括强隐过程(SIP)迭代法在内的十种迭代方法求解离散的电磁场方程.研究表明,SIP迭代法是适用于电弧喷射推力器电磁场模拟的快速有效的方法,模型能够准确地反映推力器内部电磁场与高温电离气体间的相互作用机理.     

旋成体轴对称跨声速全位势流的高效差分算法

本文用守恒型全位势方程,贴体坐标网格,对旋成体轴对称跨声速绕流的差分数值计算方法进行了研究;根据最佳收敛准则,提出了轴对称情形的AF2迭代算法,并将此算法应用于半球头柱体、弹体等各种外形的旋成体。与一般方法仅适用于亚声速自由流不同,本计算可从亚声速、跨声速自由流一直到低超声速自由流。计算结果表明,本文方法收敛快,与实验及其它方法的结果符合较好。