压气机叶片设计程序及使用说明.docVIP

程序说明 姓名：唐天全 学号：2015201171 使用说明： （BladeDesign.exe放置于BladeDesign文件夹的distrib文件夹中。 若需要在MATLAB中运行程序，需要修改MATLAB工作地址，并在命令窗口中输入“blade_design”既可运行） 默认已知进气角：50°，出气角：13°； 程序输入： 攻角，单位 “°” 落后角，单位“°” 中弧线由几段组成（安照老师要求，只需在此处输入“2”既可） 输入相邻两段中弧线的中间载荷βk，单位“°”（当输入的中弧线组成段数大于1时） 中弧线各段的线型（此处提供3种线性，分别为抛物型，指数型及圆弧型） 程序输出： 给出一个满足输入的2D叶型图。 功能说明： 该程序主要功能是实现了在给定攻角和落后角时计算出一个符合一定要求的2D叶型图。其中，中弧线可以由抛物型、指数型及圆弧型任意搭配组合而成（由用户指定）；此外中弧线的组成段数理论上可以达到任意大（由用户指定），但是当段数大于1时需要而外给出相邻两段中弧线中间载荷βk的值。 叶片厚度可以在子程序add_thickdistribution.m中通过修改“x”和“R”变量（均为矩阵）进行修改。 在绘制叶盆及叶背过程中，为了使得表面光滑，使用了三次样条插值进行光滑，在计算三次样条时为了节省时间，使用SOR迭代进行求解。 另外，给出绘制前后小圆时算法： 以叶型前缘为例，在叶型的近前缘或尾缘处的压力面上的适当位置,选取一点a（程序中取a_x=0.07）,过点作压力面叶型线的垂线,它与吸力面叶型线交于点b,将线段ab分为N个小等分di=ab入,计算机程序计算时确定o点可从a点加di开始。c点在吸力面上从前缘点开始。如果直线co的斜率与吸力面叶型线在c点的斜率相等,而且线段co的值与线段ao的值相等,则认为o点即为所找的内切圆的圆心,程序就可以o点为圆心,线段co的长度为半径做出圆及相应的圆弧。如果线段co的值与线段ao的值不相等,则o点在原有的基础上再增加di,吸力面上c点从前缘点开始重复以上过程。尾缘圆弧的求解与前缘类似。 举例： 打开BladeDesign.exe需要大概1分钟的时间缓冲。 对应于：攻角为2°，落后角为7°；中弧线由一段抛物线组成。 对应于：攻角为2°，落后角为7°；中弧线分别有一段抛物线和一段指数型曲线组成；两段中弧线间给定35°的载荷。 对应于：攻角为2°，落后角为7°；中弧线分别由抛物线、指数型曲线和圆弧曲线三段组成；每两段中弧线间的载荷分别为40°和20°。 代码 1.主程序：blade_design.m tic %% while 1 clear beita1=50; %turning_angle=60; beita2=13; i=input(incidence (unity: ??;in general,between -3?? and 5??)= ); def=input(deflect angle (unity: ??;in general,between 1?? and 15??)= ); beita1_k=beita1+i; beita2_k=beita2-def; turning_angle=abs(beita1_k-beita2_k); while 1 %n=input(how many kinds of lines you want to consist the mid-line(in general its 2): ); n=input(??ê?è??D????óé?àéù??×é3é￡¨í????é????￡?ê?è? 2￡?: ); n=fix(n); fprintf(\n); ii=0; while 1 ii=ii+1; %fprintf(the shape of line %d :\n,ii); fprintf(μú %d ??μ???Dí:\n,ii); line(ii)=input(input: 1--parabolas; 2--index curve; 3--arc curve: ); line(ii)=ceil(line(ii)); if line(ii)3 || line(ii)0 ii=ii-1; continue; end if ii==n break; end end if n1 fprintf(\n????2?·