非定常空气动力学大作业
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
非定常空气动力大作业
一、问题要求
1、采用非线性代数模型建模。
2、样本数据为某飞机模型单自由度滚转运动风洞试验中测得的滚转力矩系数（对应数据文件中“Cl ”列。数据文件名为cb0.dat-cb7.dat ，运动规律为：40cos(2)ft φπ=-︒，
分别对应运动频率0.0Hz-0.7Hz 。“φ”对应数据文件中“phi ”列。试验风速v=25m/s ，模型展长（参考长度）0.75m 。
3、要求编写建模程序（语言不限），给出源程序。
4、根据建模精度，调整系数个数，给出系数矩阵。
5、根据建模结果，计算运动规律为40cos(2)ft φπ=-︒，f ＝0.35Hz ，滚转力矩迟滞环；计算运动规律分别为2010cos(2)ft φπ=-︒-︒，2010cos(2)ft φπ=︒-︒，和10cos(2)ft φπ=-︒，f ＝0.4Hz ，滚转力矩迟滞环。
6、给出计算曲线。
实验数据
-0.06
-0.04-0.0200.02
0.040.060.08-40
-30
-20
-10
010
20
30
40
phi
C l
cb0
cb1cb2cb3cb4cb5cb6cb7
图1 原始实验数据曲线
二、模型建立
考虑一般的非线性运动规律
()1cos m a eff eff k t αααφ=-+ （1）
式中
2eff b
k f v
π=⋅⋅ （2）
其中，f 为非定常运动的频率（单位Hz ），b 为模型展长（单位m ），v 试验风速（单位m/s ）。
对于一般的非线性运动，可以建立横向非定常气动力的非线性代数模型如下：
23012345678Ca C C C C C C C C C αααααααααααα=++++++++&&&&& （3） 其中，α即为方程（1）中的1α，α
&由方程（1）求导可得 ()sin eff a eff k t ααφ=-+g
（4） 系数i c 是减缩频率eff k 的函数, 其定义如下： 与α有关的系数为
231234 0,1,2,3,4i i i i i C a a k a k a k i =+++=
（5） 与α
&有关的系数为 231234log() 5,6,7,8i i i i i C a k a k a k a k i =+++= （6）
因此，对于该模型共有36个待定系数。
三、模型求解
尽管求解模型方程（3）中的系数可以转化为一个线性最小二乘问题，但为了适用于任意形式的模型方程，本文采用非线性最小二乘逼近的方法来计算模型
方程（3）中的系数。
程序基于的软件平台为MATLAB 。源程序由包含：main.m ，myfun.m ，myfun_output.m 。
1. 计算模型中的系数并将拟合的结果与样本数据进行比较
系数 a1
a2
a3
a4
C0 0.011076794 -0.019088882 0.129860051 0.227670455 C1 -0.002145374 -0.000482186 -0.002171056 -0.009108143 C2 -4.04489E-06 -7.5293E-05 0.001403595 -0.007146387 C3 -3.3296E-06 7.25361E-05 -1.24556E-05 -0.00172726 C4 5.30213E-07 -2.12203E-06 1.00441E-05 1.5538E-05 C5 0.000749592 0.033686433 -0.190526344 -0.405792502 C6 -5.26448E-06 -0.000363143 0.004059039 -0.011643338 C7 1.34121E-05 0.000941117 -0.019001616 0.099957474 C8
0.001431945
0.015076089
-0.109972314
0.322587438
表1、模型方程（3）的系数矩阵
f=0.0
-0.06
-0.04-0.0200.02
0.040.060.08-40
-30
-20
-1001020
30
40
phi（deg）
C l
实验结果拟合结果
-0.06
-0.04-0.0200.02
0.040.060.08-40
-30
-20
-10
010203040
phi（deg）C l
实验结果拟合结果
f=0.2
-0.06
-0.04-0.0200.02
0.040.060.08-40
-30
-20
-10
01020
30
40
phi（deg）
C l
实验结果拟合结果
-0.06
-0.04-0.0200.02
0.040.060.08-40
-30
-20
-10
010203040
phi（deg）C l
实验结果拟合结果
f=0.4
-0.06
-0.04-0.0200.02
0.040.060.08-40
-30
-20
-10
01020
30
40
phi（deg）
C l
实验结果拟合结果