一种宽带阵列幅相与互耦误差联合校正算法

宽带数字阵列雷达接收通道中存在随频率变化的幅相误差和互耦误差,会严重影响雷达性能。针对这一问题,提出了一种宽带数字阵列幅相与互耦误差联合校正算法。首先选取通带内多个离散频点,对于每个频点,将幅相误差和互耦误差作为一个整体,采用基于子空间原理的窄带校正算法估计其阵列失真参数,并计算校正矩阵;然后将其组合起来,构成频域离散校正矩阵;最后基于最小二乘准则,设计了宽带有限长脉冲响应(FIR)校正滤波器矩阵。利用该矩阵,可实现通带范围内任意带宽入射信号的校正。计算机仿真实验验证了该算法的有效性。对实测数据的处理结果表明该算法在宽带数字阵列雷达系统中具有实用价值。     

Chebyshev模拟滤波器离散化问题分析与解决

针对高阶次、 低截止频率的Chebyshev模拟滤波器离散化发散问题,从滤波器阶次、采样步长、 截止频率、 离散化方法、 滤波器种类等方面进行了细致分析,得出定性原因,并针对该问题给出了通用解决方法.通过仿真实例和试验,验证了所提方法的有效性.     

基于卡尔曼滤波器的主燃油计量装置故障诊断

为了建立主燃油计量电液执行机构(EHA)的故障诊断算法,分析了主燃油计量液压机械控制装置的工作原理,提出了1种建立主燃油计量线性变参数(LPV)模型的方法,设计并得到基于离散卡尔曼滤波器的主燃油计量电液执行机构位移输出的残差生成器。仿真结果表明:该残差生成器具有对故障敏感而对噪声和建模误差不敏感、对残差的动态响应快,幅值变化较大的特点,有效降低故障误报警率和漏报警率。     

机动飞行中旋翼尾迹畸变效应实验

为准确计算机动飞行时的旋翼他轴响应,气动模型应计入尾迹畸变效应。体现此效应的关键在于推出时变的尾迹弯曲参数KR。本文通过烟流实验,得出随旋翼前进比、拉力系数和机动角速率变化的KR表达式,并分析指出,国外文献采用固定KR值不足以体现机动效应。为了给理论预测提供参考,本文还测量了机动旋翼的空气动力时间历程,表明旋翼轴接受准阶跃输入后,同轴及他轴的空气动力响应有过冲,且俯仰速率愈大,过冲幅值愈大;本实验中他轴响应为主响应的15%。     

基于矩阵运算的直升机旋翼/机体耦合动力学模型

针对直升机旋翼/机体耦合动力学建模中,小角度线性简化带来的计算精度下降以及适用范围缩小等问题,在计算桨叶速度及加速度时,基于矩阵运算,保留了运动的非线性。对不同机体运动幅值和频率情况下的小角度线性简化对桨叶速度的影响进行了仿真分析,结果发现,随着机体运动幅值的增加,线性简化影响越来越大,滚转幅值 15°时,简化带来的误差超过 5%;俯仰幅值 15°时,简化带来的误差超过 6%,且随着机体运动频率的减小,误差越来越大。对不同桨距的情况进行仿真,结果表明,所建模型可以很好地模拟直升机的动态响应,采用该模型可进行直升机旋翼/机体耦合动态响应仿真分析。     

针栓变推力固体火箭发动机动态响应特性研究

为了解针栓式固体发动机的动态响应特性,基于可变等效喉部面积的调控内弹道模型,研究了针栓型面及其调节过程对发动机动态响应的影响。结果表明:针栓式变推力固体发动机响应时间以及推力响应的动态特性与针栓移动速度、针栓型面以及推进剂压力指数密切相关。发动机响应时间由针栓移动时间和针栓移动造成的压力延迟时间组成;随针栓移动速度增大,针栓移动过程与压力延迟过程相对比重发生变化,使得响应时间先快速缩短后趋向于定值,由32ms缩短至11.6ms,而推力过冲逐渐增大,但在针栓打开过程中响应时间更短,推力过冲更大,并且计算表明压力和等效喉部面积的相对变化率决定了推力过冲;随着正压力指数的增大,压力调节范围增大,推力过冲和响应时间均增加,而负压力指数的推进剂明显缩短响应时间,抑制推力过冲;凸型面针栓能够缩短响应时间同时降低推力过冲,可根据等效喉部面积相对变化率获得最佳型面曲率,使得响应时间和推力过冲最小。     

柔性悬臂梁的振动特性与等效线性化方法的局限性

等效线性化方法是大变形柔性梁振动特性分析的一种高效近似处理方法,其成立前提是必须满足小振幅假设.传统完全非线性分析方法通常难以定量表征几何非线性效应对结构参量的影响,不利于开展定性分析,也无法满足降阶建模的需要.以大变形柔性悬臂梁为研究对象,将静态大变形对梁振动特性的影响直接等效为对截面惯性矩和质量密度分布变化的影响,构造了一种物理指示意义更加明确的等效线性化方法,并通过对比验证了该方法的有效性.为进一步讨论等效线性化方法的局限性,设计制作了一根大变形柔性悬臂梁,采用定响应幅的定频稳态激励法测量了悬臂梁在3种不同振幅水平下的频率响应函数(FRF)曲线.计算与试验结果的对比分析表明:静态大变形条件下,当悬臂尖端的动态位移振幅与静态变形幅值的比值低于10％时,结构的FRF曲线变化非常小,系统可视做线性系统,等效线性化方法的计算误差≤6％,算法有效;当悬臂尖端的动态位移振幅与静态变形幅值的比值达到20％时,测试FRF的幅值在共振点明显降低,等效线性化理论的计算结果误差偏高50％以上,等效线性化理论将不再适用;试验悬臂梁的振动非线性效应主要受阻尼的非线性特性控制.     

基于嵌入式数字信号处理器的数字滤波器设计和算法实现

介绍了用ＩＮＴＥＬ公司新推出的８０Ｃ１９６ＵＮ单片机来设计数字滤波器。用此嵌入式数字信号处理器设计的滤波器的幅值、相位和响应较为精确,消除了多个器件所引起的噪声、电压等误差。同时,本文还对Ｃ１９６ＵＮ硬件结构、软件设计及执行信号处理指令集中的关键语句做了概述     

实时干涉测量中对流层延迟与钟差精修正建模

由于中国深空干涉测量系统无法采用短时差分标校,使得实时测量精度较差,误差来源主要包括对流层延迟模型误差和钟差模型误差。为了提高实时测量精度,基于最大值最小准则,对对流层天顶延迟进行精确修正;引入邻近估计与线性回归模型,实现了高精度的钟差非线性预报。经任务数据验证,对流层延迟精修正模型预报值与实测值的差异最大值优于0.33m(仰角≥10°);相对EGNOS和"EGNOS+GMF"模型,该误差减小约1个数量级;相对线性预报,非线性预报误差减小约半个数量级。可为我国后续深空探测任务提供高精度对流层延迟和钟差建模及预报。     

拉伸板中心孔干涉配合对载荷幅值的影响

为研究干涉配合中被连接件最小截面（SCS）所传递载荷（截面载荷）的幅值得以降低的原理,按一定规律对销钉与板进行分区,并将其用同等刚度的弹簧代替,建立了销钉干涉配合的弹簧分析模型,进而从弹性角度解释了干涉配合中板的最小截面所传递载荷的幅值得以降低的原因。基于该模型分析了静载下截面载荷随外载的变化规律,得到了循环载荷下截面载荷幅值与外载幅值间的关系,其与有限元分析（FEA）结果之间的最大误差不到20%。分析发现：连接件与被连接件间的接触力随外载变化的自动调节作用是截面载荷幅值降低的原因;对于特定的绝对干涉量,干涉配合只对一定范围内的外载具有调节作用。当外载在干涉配合的有效影响区时,截面载荷幅值的降低程度完全决定于连接件与被连接件的相对刚度;若外载过大使得钉孔发生分离,干涉配合不仅对载荷幅值不再具有调节作用,还容易引起钉孔接触面的微动磨损。     

基于切比雪夫多项式展开的转子动力特性区间分析

基于切比雪夫多项式展开,建立了转子系统临界转速和稳态响应区间分析方法,能在仅具有参数上下界的情况下获得结果分布区间。针对临界转速区间分析问题,将支承刚度、密度和弹性模量视为区间变量,建立转子系统固有频率的切比雪夫多项式表达式,最终获得临界转速区间,临界转速误差不超过0.07%。在固有特性分析基础上,应用切比雪夫方法求解转子系统的稳态响应,并由模态叠加原理对响应区间进行修正,解决了固有特性不确定性较大时稳态响应求解精度不足的问题。与蒙特卡洛插值结果相比,响应上界误差不超过8%。基于切比雪夫展开的非确定分析方法,不需修改原结构的振动微分方程。结果表明该方法结合商业有限元软件后,应用简单,并具有理想的精度。     

Direct trajectory optimization based on a mapped Chebyshev pseudospectral method

 In view of generating optimal trajectories of Bolza problems, standard Chebyshev pseudospectral (PS) method makes the points' accumulation near the extremities and rarefaction of nodes close to the center of interval, which causes an ill-condition of differentiation matrix and an oscillation of the optimal solution. For improvement upon the difficulties, a mapped Chebyshev pseudospectral method is proposed. A conformal map is applied to Chebyshev points to move the points closer to equidistant nodes. Condition number and spectral radius of differentiation matrices from both methods are presented to show the improvement. Furthermore, the modification keeps the Chebyshev pseudospectral method's advantage, the spectral convergence rate. Based on three numerical examples, a comparison of the execution time, convergence and accuracy is presented among the standard Chebyshev pseudospectral method, other collocation methods and the proposed one. In one example, the error of results from mapped Chebyshev pseudospectral method is reduced to 5% of that from standard Chebyshev pseudospectral method.     

Maximizing the Usable Bandwidth of MTI Signal Processors

A technique is presented for maximizing the percentage of usable Doppler bandwidth throughout which a radar return can be detected while maintaining an acceptable clutter suppression. The technique employs the weighted Chebyshev approximation to the design of a transversal high-pass digital filter which has an optimal passband ripple for a given number of filter weights and associated integration gain consistent with the required increase in signal-to-noise ratio needed for acceptable probabilities of detection and false alarm. Conventional approaches to the design of a movingtarget arget indictor (MTI) filter which maximizes the improvement factor by clutter suppression typically improve the signal-to-background noise ratio over less than 50 percent of the range between dc and the pulse-repetition frequency fT. This technique can increase the usable bandwidth to 80 percent or more of fT. Two examples are included which utilize parameter values from the Army Missile Command's experimental radar and demonstrate the interactive influence of such filter parameters as the number of weights, passband ripple and bandedge, and stopband attenuation and cutoff.     

Interval analysis of rotor dynamic response based on Chebyshev polynomials

Uncertainty is extensively involved in the rotor systems of rotating machinery, which may cause an unstable vibrational response. To take the uncertainty into consideration for the uncertain rotor-bearing system, an improved unified interval analysis method based on the Chebyshev expansion is established in this paper. Firstly, the Chebyshev Interval Method (CIM) to calculate not only the critical speeds but also the dynamic response of rotor with uncertain parameters is introduced. Then, the numerical investigation is carried out based on the developed double disk rotor model and computation procedure, and the results demonstrate the validity. But when the uncertainty is sufficiently large to influence critical speeds, the upper and lower bounds are far from the actual bounds. In order to overcome the defects, a Bound Correction Interval analysis Method (BCIM) is proposed based on the Chebyshev expansion and the modal superposition. In use of the improved method, the bounds of the interval responses, especially the upper bound, are corrected, and the comparison with other methods demonstrates that the higher accuracy and a wider application range.     

旋转盘腔盘缘热流密度的敏感性分析

为保证涡轮盘满足适航规章的安全性要求,采用单向流固耦合数值方法,研究了转静系旋转盘腔以基比切夫数表示的盘缘加热热流密度的变化对冷却效果的影响,并依据旋转盘腔冷却问题的工程评价体系对旋转盘腔的冷却效果进行评价。研究结果表明:基比切夫数的变化对于旋转盘腔的流动结构和流动阻力基本没有影响,对盘面的换热效果影响也较微弱,仅引起转盘迎风面热流密度和温度的改变。同时,温度分布的改变导致了与温度梯度紧密相关的热应力水平发生变化。随着热流密度的增加,转盘整体应力水平上升,并且盘缘附近区域的等效应力提高的幅度大于中心区域。当基比切夫数高于临界值后,最大等效应力值从转盘中心转移到盘缘。基比切夫数的变化能够从部件承受能力和实际使用载荷两方面对涡轮盘的失效概率产生较大影响,因此,在涡轮盘腔的设计阶段,需要考虑基比切夫数对涡轮盘安全性的影响。     

超声压气机叶型设计方法

在设计超声叶型时,为使得叶栅进口马赫数和气流角等于给定值,提出一种新的叶型参数化方法。该方法以经典唯一进气角计算方法为基础,将超声叶栅的唯一进气角和叶型几何形状关联起来,由进口马赫数和气流角确定吸力面进口段叶型;根据喉部面积、前后缘小圆半径、最大挠度和最大厚度等特征参数确定其他部分叶型。用此参数化方法设计了6个超声叶型,并用Fluent对设计结果进行了验证。结果表明,来流马赫数及进气角的设计值与Fluent求解结果基本一致,进气角最大误差仅为0.7°,进口马赫数最大误差仅为0.01;并且实现了多激波增压和减小激波损失等效果。     

适用于非周期流固耦合问题的时间谱方法

针对非定常流固耦合设计、分析问题所面临的计算效率与精度、鲁棒性之间的矛盾,采取将Chebyshev谱方法与基于非定常面元和几何非线性梁有限元模型的流固耦合分析方法相结合的方式,建立了针对大展弦比机翼非定常流固耦合优化设计问题的,可与伴随方法相结合的时间谱方法。Chebyshev谱方法直接对流固耦合的控制方程进行处理,利用Chebyshev算子替换系统状态变量,将非定常问题转化为Chebyshev控制点处耦合的定常问题。通过这种方式建立的流固耦合时间谱方法具有较高的计算精度、效率和足够的鲁棒性。验证算例及Goland机翼颤振速度计算实例表明,Chebyshev谱方法的计算精度随着Chebyshev控制点个数的增加而不断增大。只需选取较少的控制点,Chebyshev谱方法便可以达到满足精度要求的计算结果。与此同时,建立的流固耦合时间谱方法不仅适用于周期性非定常问题还适用于非周期性非定常问题。     

大弯角扩压叶栅吸力面全域压力分布的PSP实验

基于自行组建的PSP(pressure sensitive paint)测量系统,采用中国科学院化学研究所PSP,实验测量了来流马赫数分别为0.4和0.5时大弯角扩压叶栅叶片吸力面全域压力分布,并与传统测压技术所得结果进行了比较.研究结果表明:①采用PSP技术测量可以获得叶片表面全域的连续压力分布,具有空间分布率高、定量测量的特点;②叶片中部PSP测量压力值与压力扫描阀结果吻合得非常好,但受到叶片吸力面曲率的影响,叶片前缘与尾缘附近的误差较大,最大误差为4.48%;③随着来流速度的提高,PSP与压力扫描阀之间的误差逐渐减小.      

Multi-scale thermodynamic analysis method for 2D SiC/SiC composite turbine guide vanes

Ceramic Matrix Composite (CMC) turbine guide vanes possess multi-scale stress and strain with inhomogeneity at the microscopic scale. Given that the macroscopic distribution cannot reflect the microscopic stress fluctuation, the macroscopic method fails to meet the requirements of stress and strain analysis of CMC turbine guide vanes. Furthermore, the complete thermodynamic properties of 2D woven SiC/SiC-CMC cannot be obtained through experimentation. Accordingly, a method to calculate the thermodynamic properties of CMC and analyze multi-scale stress and strain of the turbine guide vanes should be established. In this study, the multi-scale thermodynamic analysis is investigated. The thermodynamic properties of Chemical Vapor Infiltration (CVI) processed SiC/SiC-CMC are predicted by a Representative Volume Element (RVE) model with porosity, leading to the result that the relative error between the calculated in-plane tensile modulus and the experimental value is 4.2%. The macroscopic response of a guide vane under given conditions is predicted. The relative error between the predicted strain on the trailing edge and the experimental value is 9.7%. The calculation of the stress distribution of micro-scale RVE shows that the maximum value of microscopic stress, which is located in the interlayer matrix, is more than 1.5 times that of macroscopic stress in the same direction and the microscopic stress distribution of the interlayer matrix is related to the pore distribution of the composite.