基于ONERA失速气动力模型,建立二元大攻角失速翼型非线性颤振运动的状态方程,通过龙格库塔法进行数值积分,获得二元翼型的非线性颤振时域响应曲线,并通过二分法寻找到失速翼型的颤振临界速度。仿真预测了半弦长对失速翼型气动弹性的影响规律,计算发现,随着半弦长的值变大,颤振速度会发生突跃,当弦长为0.5 m时,颤振速度最大,发生颤振的难度最大。