没有合适的资源?快使用搜索试试~ 我知道了~
温馨提示
采用离散傅里叶变换方法对颤振时间响应历程计算结果进行分析,提出了一种从多个模态分支计算结果中归纳耦合主颤振分支的方法。分别从频率移动理论和能量变化观点出发,提出了针对颤振边界的临界颤振动压判定频率移动判据和能量因子判据。通过对国际标准跨声速颤振算例AGARD445. 6机翼的时域颤振计算结果进行分析,验证了方法在耦合模态分支确定、临界颤振动压判定、主颤振分支判定3个方面的高有效性。此外,将该方法应用于高速中等展弦比气动翼面的跨声速颤振特性研究,成功地根据较复杂的响应曲线判断出了颤振边界,表明该方法具有良好
资源推荐
资源详情
资源评论
!"#$
年
#"
月
第
%!
卷第
&
期
西北工业大学学报
’()*+,- (. /(*0123403*+ 5(-60371+87,- 9+8:3*4806
;70<
=(-< %!
!"#$
/(< &
收稿日期
:!"#$>"%>#?
基金项目
:
国家自然科学基金
(#"E"!"$E)
与航天科技创新基金重点项目
(IDPI"###)
资助
作者简介
:
秦文轩
(#A??—)
,
西北工业大学硕士研究生
,
主要从事颤振计算研究
。
一种时域颤振判定新方法研究
秦
之轩
,
史爱明
(
西
北工业大学 航空学院
,
陕西 西安
@#""@!)
摘 要
:
采用离散傅里叶变换方法对颤振时间响应历程计算结果进行分析
,
提出了一种从多个模态分
支计算结果中归纳耦合主颤振分支的方法
。
分别从频率移动理论和能量变化观点出发
,
提出了针对
颤振边界的临界颤振动压判定频率移动判据和能量因子判据
。
通过对国际标准跨声速颤振算例
DfDR_$$&B E
机翼的时域颤振计算结果进行分析
,
验证了方法在耦合模态分支确定
、
临界颤振动压判
定
、
主颤振分支判定
%
个方面的高有效性
。
此外
,
将该方法应用于高速中等展弦比气动翼面的跨声速
颤振特性研究
,
成功地根据较复杂的响应曲线判断出了颤振边界
,
表明该方法具有良好的工程应用
前景
。
关 键 词
:
时域判颤
,
幅频分析
,DfDR_$$&B E
机翼
,
频率移动判据
,
能量因子判据
中图分类号
:=!#&B % X $
文献标志码
:D
文章编号
:#""">!@&?(!"#$)"&>"@"">"@
从非线性力学角度来看
,
颤振属于一种复杂的
非线性自激振荡问题
[#]
。
颤振问题的研究手段包
括
地面实验
[!>%]
、
颤振计算和颤振试飞
。
其中颤振
计算是设计初期成本较低的方法
。
时域方法
[$>&]
与
传统的频域方法相比其具有气动力求解精度较高
,
能够直观给出振动响应历程等优点
。
但是
,
时域方
法也存在一些缺陷
:
无法直接提供振动的频率特性
,
不易在大量的模态分支中快速筛选出耦合模态分
支
,
对于临界颤振点的判断缺乏有效的
、
定量的判
据
,
主要依赖人工经验
,
主观因素过强
。
本文将针对时域计算方法的不足开展研究工
作
,
具体研究思路如下
:
提出一种临界颤振点的判断方法
,
该方法包含
!
个方面
:
①使用傅里叶变换
[E]
对时域颤振计算得
到的广义位移响应进行变换
,
得到各阶模态的主频
并将其作为振动的频率
,
绘制出各分支模态的频率
随来流动压变化的曲线
,
频率有互相接近趋势的两
阶模态即为频率耦合模态
,
将频率最靠近点作为临
界颤振点
;
②计算基于广义位移和广义速度的弹性
势能和振动动能
,
绘制出总能量随时间变化的曲线
,
以能量因子的正负和大小表征振动收敛和发散的快
慢
,
然后计算出能量因子的零点
,
就是临界颤振点
。
最后将本文的方法应用于实际的算例进行检验
。
1
时域颤振方程
时域颤振计算的振动控制方程为
:
M¨x + C x + Kx = F(t)
(#)
式中
: M、C、K
依
次是振动系统的质量
、
阻尼和刚度
矩阵
,¨x、x、x
依次是三维笛卡尔物理坐标下振动系
统的加速度
、
速度和位移列向量
,F(t)
是非定常气
动力
。
对于
N
个自由度的振动系统
,
这是一个包含
N
个相互耦合的二阶常微分方程的方程组
,
其求解比
较困难
,
因此需要将其从物理坐标转化到模态坐标
下
,
对其进行解耦
。
设模态坐标下
,
广义位移用 ξ 表示
,
则广义位移
ξ 和实际位移
x
的关系是
:
x =
Φξ
(!)
式中
:
Φ 是振动系统的模态矩阵
,
将方程
(!)
代入方
程
(#)
,
然后在方程
(#)
两边同时左乘 Φ
T
得
:
Φ
T
M
Φ
¨
ξ
+
Φ
T
C
Φ
ξ
+
Φ
T
K
Φξ
=
Φ
T
F(t) (%)
令
:
.
M =
Φ
T
M
Φ
第
&
期
秦之轩
,
等
:
一种时域颤振判定新方法研究
/
C =
Φ
T
C
Φ
/
K =
Φ
T
K
Φ
/
F(t) =
Φ
T
F(t)
.
M、
/
C、
/
K
分别称为系统的模态质量矩阵
、
模
态阻尼矩
阵和模态刚度矩阵
,
/
F(t)
称为广义力列向量
。
由
此
,
振动方程转化为模态坐标下的形式
:
.
M
¨
ξ
+
/
C
ξ
+
/
K
ξ
=
/
F(t) ($)
其中的
.
M、
/
C、
/
K
都是对角矩阵
,
因此在模态坐标下
振动方程变成了
#
组包含
N
个互相独立的二阶常微
分方程的方程组
,
这就是时域颤振计算使用的控制
方程
。
2
颤振耦合分支确定
本节的内容是对时域颤振计算得到的响应曲线
进行傅里叶变换
,
确定颤振耦合分支
。
研究对象为
颤振计算标模
DfDR_$$&B E
机翼
[@]
,
计
算
状
态为
:
来流马赫数
"B A,
图
#
给出了不同来流动压下机翼
第一阶模态
(
一阶弯曲
)
广义位移的响应曲线
,
图中
收敛
、
等幅
、
发散的曲线分别对应来流动压为
# A#!
UL ^(H·4
!
)、$ @?@ UL ^(H·4
!
)、& @$# UL ^ (H·4
!
)
的流动状态
。
这里之所以只给出第一阶模态的响应
曲线
,
是为了在能够反映振动状态随来流动压增加
的变化趋势前提下
,
尽可能地使图片清晰易读
,
其它
模态广义位移的变化趋势与第一阶模态是一致的
。
本文确定耦合分支的方法是基于频率移动原
理
,
相互耦合的分支频率在颤振点前后都会有靠近
的现象
,
因此只要所研究的状态在颤振点附近即可
。
研究中分别对广义位移
、
广义力
、
实位移
、
物面
压力的时间响应进行了傅里叶变换
。
图
!
给出了上述
%
个来流动压下广义位移在进
行离散傅里叶变换后的正弦分量的幅频特性曲线
。
为了能更加清楚地观察到峰值的位置
,
图中只显示
了频率范围在
" C $" bO
之间的曲线
。
图中
,Q S #
A#! UL ^ (H·4
!
)
时
(
振动收敛状态
),
第一阶和第二
阶模态的主频分别为
##B A bO
和
%!B A bO;
当来流动
压增大到
$ @?@ UL ^( H · 4
!
)
时
(
等幅振动状态
),
一
、
二
阶模态的主频分别变为
#@B $ bO
和
#?B % bO,
出现了明显的频率靠近现象
,
可以断定第一和第二
阶模态为颤振耦合模态分支
。
随着来流动压继续增
大
,
这两阶模态的频率一直保持着接近状态
,
同时振
幅不断增大
,
这符合颤振耦合模态分支的特性
。
余
弦分量的分析结果与正弦分量完全一致
,
如图
%
所示
。
图
# Ma S "B A
时不同来流动压下
DfDR_$$&B E
图
!
广义位移正弦分量 图
%
广义位移余弦分量
机翼广义位移响应曲线
广义力响应的频率特性与广义位移完全一致
,
傅里叶变换的结果也印证了这一点
。
广义力的分析
结果依然显示第一阶和第二阶模态为耦合模态分
支
,
且给出的主频与广义位移的结果完全相同
。
3
颤振边界定量判据
确定耦合模态分支后
,
本节的主要内容是分别
·#"@·
剩余6页未读,继续阅读
资源评论
weixin_38637272
- 粉丝: 4
- 资源: 935
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- 基于Django和HTML的新疆地区水稻产量影响因素可视化分析系统(含数据集)
- windows conan2应用构建模板
- 3_base.apk.1
- 基于STM32F103C8T6的4g模块(air724ug)
- 基于Java技术的ASC学业支持中心并行项目开发设计源码
- 基于Java和微信支付的wxmall开源卖票商城设计源码
- 基于Java和前端技术的东软环保公众监督系统设计源码
- 基于Python、HTML、CSS的crawlerdemo软件工程实训爬虫设计源码
- 基于多智能体深度强化学习的边缘协同任务卸载方法设计源码
- 基于BS架构的Java、Vue、JavaScript、CSS、HTML整合的毕业设计源码
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功