一种时域颤振判定新方法研究.pdfVIP

2014年 l0月 西 北 工 业 大 学 学 报 0ct． 2014 第32卷第5期 JournalofNorthwestern PolytechnicalUniversity Vo1．32 No．5 一 种时域颤振判定新方法研究 秦之轩，史爱明 (西北工业大学 航空学院，陕西 西安 710072) 摘 要：采用离散傅里叶变换方法对颤振时间响应历程计算结果进行分析，提出了一种从多个模态分 支计算结果中归纳耦合主颤振分支的方法。分别从频率移动理论和能量变化观点出发，提出了针对 颤振边界的临界颤振动压判定频率移动判据和能量因子判据。通过对 国际标准跨声速颤振算例 AGARD445．6机翼的时域颤振计算结果进行分析，验证了方法在耦合模态分支确定、临界颤振动压判 定、主颤振分支判定3个方面的高有效性。此外，将该方法应用于高速中等展弦比气动翼面的跨声速 颤振特性研究，成功地根据较复杂的响应 曲线判断出了颤振边界，表明该方法具有 良好的工程应用 前景。 关 键 词：时域判颤，幅频分析，AGARD445．6机翼，频率移动判据，能量因子判据 中图分类号：V215．3+4 文献标志码：A 文章编号：1000．2758(2014)05—0700—07 从非线性力学角度来看，颤振属于一种复杂的 最后将本文的方法应用于实际的算例进行检验。 非线性 自激振荡 问题_】J。颤振 问题的研究手段包 括地面实验 。、颤振计算和颤振试飞。其中颤振 1 时域颤振方程 计算是设计初期成本较低的方法。时域方法 与 传统的频域方法相比其具有气动力求解精度较高， 时域颤振计算的振动控制方程为： 能够直观给出振动响应历程等优点。但是，时域方 脚 + +Kx=F() (1) 法也存在一些缺陷：无法直接提供振动的频率特性， 式中：M、C、K依次是振动系统的质量、阻尼和刚度 不易在大量的模态分支中快速筛选出耦合模态分 矩阵，、、依次是三维笛卡尔物理坐标下振动系 支，对于临界颤振点的判断缺乏有效的、定量的判 统的加速度、速度和位移列向量，F(t)是非定常气 据，主要依赖人工经验，主观因素过强。 动力。对于J7、r个 自由度的振动系统，这是一个包含 本文将针对时域计算方法的不足开展研究工 Ⅳ个相互耦合的二阶常微分方程的方程组，其求解 作，具体研究思路如下： 比较困难，因此需要将其从物理坐标转化到模态坐 提出一种临界颤振点的判断方法，该方法包含 标下，对其进行解耦。 2个方面：①使用傅里叶变换 对时域颤振计算得 设模态坐标下，广义位移用 表示，则广义位移 到的广义位移响应进行变换，得到各阶模态的主频 和实际位移 的关系是： 并将其作为振动的频率，绘制出各分支模态的频率 = (2) 随来流动压变化的曲线，频率有互相接近趋势的两