欢迎来到天天文库
浏览记录
ID:16641314
大小:92.50 KB
页数:4页
时间:2018-08-23
《机翼翼型的matlab编程》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、1.编程界面类似下图,再加一个图形显示框和一个生成三维坐标数据的按钮。该模块的功能是根据用户指定的中线形式以及相关的几何参数,选择相应计算公式来计算叶片截面上、下轮廓的坐标数据点,并将计算结果分别存入数组中。b——叶片弦长q叶片转折角叶片安装角和最大厚度不需要。2.原始翼型数据见txt文件,两种即可,文件名NACA4309和NACA4412,文件中存放的是翼型上缘和下缘的二维坐标数据。Matlab将从此文件中读取二维坐标数据。注意:其中x坐标为弦长b的百分数。3.翼型中线形式,圆弧形和抛物线型两种。翼型中线通常采用圆弧形或
2、抛物线形。圆弧形中线可由两段或一段圆弧构成。(1)对于图1所示的双圆弧中线,(03、翼型提供型值点数据代入计算;yz值根据不同形式的中线分别由式(1)、(2)、(3)确定;中线上任意点的切线方向与x轴之间的夹角φ可由下式计算:(5)4.我想要的程序是能够读入原始坐标参数,然后通过计算得到最终的坐标参数,然后存储为txt或者dat格式的文件。如果你感觉数据缺失也没有关系,你只需要帮我把整体程序做一下,然后我自己再修改。clcclear%x,ye的值可由所选原始翼型提供型值点数据代入计算,这里假设二者均已知%先计算yz的值,以式(1)为例,其他类似%假设θ、a、b(若未知,再讨论)sita=(θ具体的值);a4、=(a的具体值);b=(b的具体值);%计算yz的值symsxyzdyzfaiyz=((a/sin(0.6*sita)^2-(x-a)^2))^0.5-a/(tan(0.6*sita));yzvalue=subs(yz,x);%注意这里的x就指的是一个具体的x的值,由于不知道x的具体形式,这里只讨论了x为一个具体数值的情况,若x为数组编程稍加修改即可%下面计算φdyz=diff(yz,x);fai=subs(atan(dyz),x);%这里的x也跟上述一样是具体的数值,不是参数%下面叶片截面上轮廓xu,yuxu=x-ye*5、sin(fai);yu=yzvalue+ye*cos(fai);%叶片截面下轮廓与此类似数据的导出比较简单,若能调通,再编写建议取一个x及ye的值代入计算,看运行结果是否满足要求另我对题目的理解可能有问题,再联系
3、翼型提供型值点数据代入计算;yz值根据不同形式的中线分别由式(1)、(2)、(3)确定;中线上任意点的切线方向与x轴之间的夹角φ可由下式计算:(5)4.我想要的程序是能够读入原始坐标参数,然后通过计算得到最终的坐标参数,然后存储为txt或者dat格式的文件。如果你感觉数据缺失也没有关系,你只需要帮我把整体程序做一下,然后我自己再修改。clcclear%x,ye的值可由所选原始翼型提供型值点数据代入计算,这里假设二者均已知%先计算yz的值,以式(1)为例,其他类似%假设θ、a、b(若未知,再讨论)sita=(θ具体的值);a
4、=(a的具体值);b=(b的具体值);%计算yz的值symsxyzdyzfaiyz=((a/sin(0.6*sita)^2-(x-a)^2))^0.5-a/(tan(0.6*sita));yzvalue=subs(yz,x);%注意这里的x就指的是一个具体的x的值,由于不知道x的具体形式,这里只讨论了x为一个具体数值的情况,若x为数组编程稍加修改即可%下面计算φdyz=diff(yz,x);fai=subs(atan(dyz),x);%这里的x也跟上述一样是具体的数值,不是参数%下面叶片截面上轮廓xu,yuxu=x-ye*
5、sin(fai);yu=yzvalue+ye*cos(fai);%叶片截面下轮廓与此类似数据的导出比较简单,若能调通,再编写建议取一个x及ye的值代入计算,看运行结果是否满足要求另我对题目的理解可能有问题,再联系
此文档下载收益归作者所有