跨音速流中的一种边界元技术

边界元的名称源于C.A.Brebbian,边界元技术的特点是节点未知量均位于物面边界上,从而大大简化数值求解过程和降低对计算机容量的要求。本文将边界元技术推广应用于三维机翼的跨音速绕流的非线性问题,建立了一种适用于跨音速范围的带权余数表示式,并将其直接应用于全速位方程。为了保证冲波的不可逆特性,并使解在超临界流区域容易收敛,我们在全速位方程中增加有人工粘性项,而且将流场区域划分成许多小子域,并对每个子域应用格林定理。如果对每个子域我们所选用的插值函数对速位φ具有C~1连续性,则相邻元素交界面的面积分将相互抵消,使剩余的元素面积分(因而节点未知量)均位于机翼边界面上,从而导得边界元的积分方程。     

低雷诺数粘性流动的数值解

本文采用格林函数方法,使问题降低一个维度,计算了低雷诺数粘性流动,给出边界积分中几何系数的解析表示式,并讨论了计算区域选取对计算的影响。     

用口诀,巧算重积分

在重积分的计算中,一般解题步骤是:先画出积分区域图,由区域的边界方程来确定积分变量的上、下限,将重积分化为累次积分。以二重积分例:在直角坐标系下,若区域D为a≤x≤b,V_1(x)≤y≤甲_2(x) 如图一所示,则采用公式一,     

航空板结构强度的边界元法

边界元即边界积分方程方法,可方便地用来解板弯曲问题。本文在局部极坐标下建立公式,采用具有较高协调性的函数插值方案,利用外点技术完全避免边界奇异积分。列举的一系列算例表明该方法输入数据少、机时省、精度高,是板结构强度分析的一个有效方法。     

考虑地面影响的半空间气动噪声数值计算方法

含非紧致边界的半空间气动噪声问题普遍存在于航空领域,在数值计算过程中不仅要考虑非紧致边界的散射效应,还需要考虑地面对声波的反射。基于Lighthill声模拟理论和广义格林函数方法,引入边界声压积分项和可压缩修正项可以考虑非紧致边界以及流体可压缩性对声波散射的影响。然而在半空间声场的计算过程中,包含沿无限大地面的积分项,为了克服沿地面积分所带来的巨大工作量,根据无穷大平板对声波的反射特性,推导得出频域下基于半空间格林函数的声学积分方程。同时开展了半空间中应力源噪声与圆柱绕流噪声的数值预测与分析,验证了该方法的准确性与高效性。     

不规则边界上的近似超奇性边界积分方程

提出角点、棱边处不规则边界上的近似超奇性边界积分方程，其依据是边界上任一Cauchy或Hadamard主值积分可以由2个相应的近奇性边界积分的平均值近似表达。近似方程与光滑边界上的超奇性边界积分方程具有统一的表达形式，并且不需要正则化过程，有利于简化程序设计。给出的数值算例说明了近似方程的可行性与有效性。     

二维可压湍流边界层改进的积分方法

采用含高阶项的边界层动量积分方程,同时对常规边界层卷吸方程及延迟方程进行相应的高阶影响修正,得到改进的边界层积分方程组。应用此方程组对多种高亚音速及跨音速翼型边界层流动作了计算,并与一阶Ｇｒｅｅｎ方法计算结果及实验结果作了比较。结果表明,在边界层积分方程中保留法向压力梯度项及雷诺法向应力项,明显改进了翼面边界层接近分离区域处参数的计算精度。     

背腔式缝隙天线阵列的分析与设计

采用等效磁流分析方法,将波导激励的单背腔缝隙单元阵场求解问题离散,得到缝隙表面的边界积分方程。应用矩量法求解该边界积分方程,得到单背腔缝隙单元阵的场结构。再应用阵列天线理论仿真计算了一种易于共形的缝隙天线阵的方向图,与实际测量结果进行了对比,取得了较好的效果。     

实现边界元并行分布计算的一种新方法

 给出了实现边界元并行分布计算的一种新方法。它将边界元离散节点分为与计算网络所拥有的处理机台数相对应的若干组, 以此实现任务分配, 独立地计算形成每组节点对应的离散边界元积分方程, 最后利用迭代方法行求解分布的边界元方程组。以弹性问题为例, 在Inmos T800 Transputer 系统上实现了边界元计算。验证了本文方法的可行性与有效性。     

Hoff型夹层板振动分析的边界积分方程法

 本文利用静力问题的基本解建立了Hoff型夹层板振动分析的边界积分方程。在数值实现过程中,边界上和区域内的未知函数均釆用线性插值以提高计算精度。数值结果表明本文方法具有未知量少和计算精度高等优点。对于低阶频率,本文结果与精确解的误差不大于1％。     

直升机“地面共振”力学模型的探讨

 <正> I_B 桨叶绕垂直铰的转动惯量I_(?),I_(?),I_(?) 机身绕过其重心O的x,y,z 轴的转动惯量ι 桨叶长度m_B 单片桨叶的质量     

用杂交应力元及等参元计算应力强度因子

 本文介绍一种用奇异性杂交应力元及等参元计算平面裂纹的应力强度因子和J积分的方法。实际计算表明,此种方法的优点是精度高,计算时间短,便于实际应用。文中给出三个算例。各例分别采用216～256个自由度。计算结果同解析解及其他奇异元的解进行了比较。各例算得的应力强度因子值,均同解析解相差仅1～2％。每例中由4～6条路径求得的J积分值,其分散度均在±O.7％以内。 计算在TQ-16机上进行,程序编译约6～7分钟,每例的计算时间约6～8分钟。     

二维机翼二自由度跨音速颤振分析

 <正> On the assumption that airfoils take small amplitude ha-rmonic oscillating and middle reduced frequency, the two-dimensional, inviscid, small-disturbance nonlinear unsteady transonic flow equation can be split into two parts. One is the nonlinear steady small-disturbance equation, the other is the time-linearized unsteady small-disturbance equation. The steady equation is solved by Carlson code. The unsteady equation is solved by integral method. The standard v-g method is used to solve the flutter eigenvalue equations. A transonic flutter analysis is performed for a NACA 64A006 airfoil with pitching and plunging degrees of freedom.The Mach numbers considered are 0.7, 0.8, and 0.85. The aerodynamic coefficients are obtained by the above method. For each Mach number, the flutter speed and the corresponding flutter reduced frequencies are achieved by varying the airfoil-airmass ratio, plunge-pitch frequency ratio. The results are compared with other computational method and good agreement has been observe     

一个n阶微分方程的两点边值问题

有一些n阶非线性微分方程的两点边值问题已被人们讨论过,大多是在微分方程的右端函数有界、满足李普希兹条件或在相应边值问题存在上下解的情况下讨论的。[3]、[4]在这篇文章里,讨论了一个n阶非线性微分方程的两点线性边值问题,是首先通过n-2次累次积分将原方程化成了一个二阶微积分方程边值问题,并用拓扑度理论讨论了解的存在性,同时给出了解的唯一性条件。推广了文献[1]的结果。     

时间推进积分方程法求解非定常跨声速流动

本文提出了一种求解非定常跨声速流动的新方法——时间推进积分方程法,此法克服了时间线化积分方程法的限制,能较好地模拟激波的运动。本文首先用一维波(?)问题——模型问题阐明此法的基本思想,然后将它应用于二维非定常跨声速流动中。本文还首次引入拟速度位的概念,使时间推进积分方程式得到简化,尾涡条件和Kutta条件更易处理。数值计算表明时间推进积分方程法是合理可靠的。     

Two freedom linear parameter varying μ synthesis control for flight environment testbed

To solve the problem of robust servo performance of Flight Environment Testbed(FET)of Altitude Ground Test Facilities(AGTF) over the whole operational envelope, a two-degree-offreedom μ synthesis method based on Linear Parameter Varying(LPV) schematic is proposed, and meanwhile a new structure frame of μ synthesis control on two degrees of freedom with double integral and weighting functions is presented, which constitutes a core support part of the paper. Aimed at the problem of reference command's rapid change, one freedom feed forward is adopted, while another freedom output feedback is used to meet good servo tracking as well as disturbance and noise rejection; furthermore, to overcome the overshoot problem and acquire dynamic tuning,the integral is introduced in inner loop, and another integral controller is used in outer loop in order to guarantee steady errors; additionally, two performance weighting functions are designed to achieve robust specialty and control energy limit considering the uncertainties in system. As the schedule parameters change over large flight envelope, the stability of closed-loop LPV system is proved using Lyapunov inequalities. The simulation results show that the relative tracking errors of temperature and pressure are less than 0.5% with LPV μ synthesis controller. Meanwhile, compared with non-LPV μ synthesis controller in large uncertainty range, the proposed approach in this research can ensure robust servo performance of FET over the whole operational envelope.     

梁/液体耦合振动的边界积分对称列式法

 利用格林积分导出的液体内部扰动压力函数作为测试函数,通过变分得到一组具有对称系数的代数方程,将该方程同梁的有限元振动方程联立后,即可求解梁在液体中的耦合振动问题。本方法具有未知数目少,数据存储和计算效率高的特点。     

硅基内衬喷管扩张段烧蚀和温度场的耦合计算方法

本文叙述固体火箭喷管扩张段硅基内衬烧蚀和温度场的耦合计算方法.烧蚀采用非定常的液体层烧蚀模型;温度场计算采用圆柱坐标系下的一维径向瞬态导热模型,考虑材料在高温下形成的炭化层和热解层,而由烧蚀产生的边界移动问题用坐标变换的方法处理.利用烧蚀和导热计算方程组中二者的结合部,通过充分耦合,获得了更准确的预示计算结果.     

可压缩紊流分离边界层的计算

本文介绍了一种求解二维可压缩紊流分离边界层的正反积分方法。该方法以边界层动量积分方程和平均流动能积分方程作为基本方程,采用适用于贴附和分离紊流的Swafford速度型解析式,对比计算了四种代数涡流粘性模型对沿流向各点的耗散积分及诸边界层特征参数的影响。计算结果和存在激波边界层较强干扰的跨音速翼型边界层的测量数据吻合。该方法可扩展用于计算翼型的尾迹特性。     

定常可压缩无粘流非线性方程解的完全边界积分表示式

 本文以ρV为因变量将可压缩等熵无粘流非线性方程的解表示成完全边界积分形式,从而将问题的计算维数降低一维。