发动机叶片数学建模.doc

1.问题重述 航空发动机叶片及其截面线的形状如图1所示，它可以分为叶盆、叶背、前缘和后缘四个部分，其中前缘的形状对整个叶片的气动性能具有非常重要的作用。传统的方法是将前缘设计成一段圆弧，则这段圆弧的圆心坐标和半径（也称前缘半径）r就是重要的前缘参数。但是近年来越来越多的实验和数值研究显示，采用非圆弧形前缘，例如椭圆弧形前缘可以明显改善叶片的气动性能。如果是椭圆弧形的前缘，则椭圆弧的中心坐标、长短半轴a、b，倾角（椭圆长轴方向与y轴的夹角）就是重要的前缘参数。 图1 航空发动机叶片及其截面线 图2. 测量得到的截面线前缘附近的离散数据 在某型航空发动机的仿制过程中，如图2所示测量了截面线前缘附近的一组离散数据。它们的坐标值记录在附件qianyuan_data.txt中, 文件中的每一行代表一个数据点的坐标（第一个值为横坐标x，第二个值为纵坐标y）。由于在测量时前缘数据与相邻的叶盆叶背数据无法直接分离，故其中混有相邻的叶盆叶背数据，请建立数学模型，解决以下问题： 根据所给数据，通过计算分离出前缘数据，并判断出前缘的形状是一段圆弧还是一段椭圆弧，计算出相应的前缘参数。将整个由离散数据出发，经过分离前缘数据、判断前缘类型、计算前缘参数的过程用一个算法描述。 评价你所给出模型的精确度和稳定性（指如果给数据添加均值较小的随机扰动，是否还能得到较高的计算结果），给出计算实例。 能否对1）中计算前缘参数的方法进行改进提高计算精度，给出具体的方法并比较改进前后的计算结果。 2.模型的假设与符号说明 2.1.模型的假设： （1）.实际测量所得到的前缘附近数据时离散的，可以假设各点的连线是一段 光滑曲线。 （2） 2.2.符号说明： （x01,y01）为圆弧的圆心坐标，r为圆弧的半径； （x02,y02）为叶片前缘椭圆弧的中心坐标，a,b为叶片前缘椭圆弧的长短半轴； 倾角β叶片前缘椭圆长轴方向与y轴方向的夹角。 3.模型的分析 从混有相邻的叶盆叶背数据中，分离出前缘数据，必须找到只有前缘数据才满足的条件，如方程，将数据逐一代人，条件成立的数据，即为前缘数据，实现前缘数据的筛选。 因为叶片前缘是一段圆弧或椭圆弧，叶片前缘的数据都应满足圆或椭圆的方程，而叶盆叶背的数据则不满足此方程。就能做出判断数据是否属于叶片前缘。 求出叶片前缘数据所满足的方程是问题的关键。 假设叶片前缘是一段椭圆弧，求出椭圆一般方程需知道圆弧上的4个点的坐标。 4.对问题一模型的建立与计算 4.1建立： 由于圆是椭圆的特殊情况，所以可以直接设叶片前缘类型为椭圆弧，根据椭圆的一般方程建立模型：Ax2+By2+Cx+Dy+E=0(其中A,B,C,D,E为椭圆方程的系数). 如果A=B则为圆； 如果A≠B则为 椭圆。 4.2.计算： 第一步：分析截面线前缘附近数据点，可得到相对准确的四个点在椭圆弧上，分别是（-18.9572,-0.241143）,(-18.9982,-0.0403025),(-19.1117,0.699986),(-19.0947,0.886326),(-18.7235,2.00039)四个点。如果代入四个点坐标，可得一个五元一次线性方程组： 由于线性方程组中常数项全为0，故此方程组为齐次线性方程组，而系数行列式为 ： 由此行列式可由MATLAB软件编程可求得，故次方程组只有零解。 第二步：为了能求出合理的系数，可设方程组的解E=-K（其中K∈R）.则可得非齐次线性方程：Ax2+By2+Cx+Dy=K。 则在上述分析所得到的五组数据中任意找4点，即（-18.9752；-0.241143） ，代入方程：Ax2+By2+Cx+Dy=K，得到非齐次线性方程组： 常数项行列式： 则Gx=H. 该方程的解可由附录1中的程序在MALAB软件上运行求得，为了使数据更精确，此处表示为分数形式：A=…*k,B=…*k,把所求得系数代入方程则可得（0.0027x2+0.00014y2+0.1044x+0.00018y） 展开得：b2x2+a2y2-2b2x0x-2a2y0y+b2x02+a2y02-a2b2=0 由第二步可得到椭圆一般方程Ax2+By2+Cx+Dy=k 两边同除以k并移项得（Ax2+By2+Cx+Dy）-1=0