轴流风机机翼型叶片参数化建模方法

1、LU也網仅I匕轴流风机机翼型叶片参数化建模方法马静王振亚同济大学汽车学院上海（201804）Email:basei摘 要：本文通过创建翼型模板，结合Matlab与UG软件,探讨了风机翼型叶片参数化建模的方法，给出了翼型中线为圆弧时的翼型坐标算法、各截面安装角和站位的处理方法以及 Matlab实现程序。并提出了叶片在UG建模时应注意的问题。文中提出的方法，减少了风机建模的工作量，缩短了风机CFD前处理周期，提高了风机流场CFD分析计算的效率和质量。关键词：叶片；参数化设计；UG Matlab1. 前言随着CFD技术的迅速发展，对风机流场计算分析的要求越来越多。风机仿真计算的前期工作量相当大，主要

2、表现在机翼型叶片的建模，其中包括风机叶轮的机翼型叶片，机翼型前导流叶片和叶轮后的止旋片建模。通常在UG软件中输入大量的翼型坐标点是相当麻烦的，而使用*.dat文件导入这些数据的方法要方便的多，但是对不同的叶片计算截面采用*.dat文件手工导入翼型坐标点的工作量仍然非常大，并且修改起来也不方便。通过分析可知，叶片不同计算截面的翼型曲线是相似的，同种翼型只因弧长以及中线形状不同而不同，因此完全可以考虑采用参数化建模的设计方法。采用这种方法可以缩短建模时间，节省大量的工作量， 且所建的模型也易于修改。因为在对风机流场进行CFD分析计算时改变风机叶片翼型是对风 机模型的重大修改需要花费大量的时间，有了

3、这种方法可以较轻松的完成修改。本文就是基于这种思想，介绍了用 Matlab与UG两个软件结合进行风机叶片参数化建模的方法，本方法 利用Matlab强大的数据处理能力处理翼型离散点1，用UG强大的三维曲面建模能力构建叶片复杂曲面。2. 翼型离散点的参数化处理2.1翼型模板的建立翼型模板的建立是实现参数化设计的第一步，建立翼型模板库是一个积累的过程， 需要将每次用到的翼型和收集到的有价值的翼型参数通过手工输入，建立起翼型模板库，在进行风机叶片建模时就可以非常方便的从翼型模板库里直接调出所需要的翼型。在Matlab中可以通过一个两列矩阵建立起翼型模板，第一列输入原始翼型的 x/l值，第 二列输入原始

4、翼型的y/l值，再存成以翼型名称命名的.mat文件即可。有了翼型模板，在进 行叶片建模时就可以根据各计算截面的翼型弦长，直接得到原始翼型各离散点的坐标值。2.2实际翼型的处理实际翼型有很多种，大部分都是在原始翼型的基础上通过改变翼型中线的形状得到的，其中最常见的是中线为圆弧或多段圆弧组合的翼型，因此这里以圆弧翼型为例，给出翼型坐标点的变换算法及 Matlab程序，多段圆弧可作相似处理。此算法有三个目的：1. 使原始翼型变为中线为圆弧的翼型2. 使翼型的几何中心与坐标原点重合，以便于后期叶片拉伸建模方便。3. 使翼型的安装角与设计相同，以便于后期叶片建模。算法思想：先在圆柱坐标系下进行变换，以中

5、线圆弧圆心为原点，用半径r和角度e坐标确定翼型各离散点位置，再转换到笛卡尔坐标系中，以便于UG导入。算法实现：1. 将原始翼型的x坐标转换为圆柱坐标系下的e坐标：Xe= aXAl其中：a -翼型中线圆弧所对应的圆心角l -翼型中线圆弧长2. 将原始翼型的Y坐标转换为圆柱坐标系下的r坐标：r = R+Y其中：R-翼型中线圆弧半径3. 再将翼型安装角转为零:24. 将翼型转至需要的安装角B :e2 = e + B5. 将圆柱坐标系下坐标 （r, e）转化到相同原点的笛卡尔坐标系下的坐标（x, y）:x = r Xcos（e）； y = r Xsin（e2）最后将翼型几何中心移到坐标原点：3宙禺龜技

6、逗兗在线xr = x- Xc yr = y - yc其中：xc, yc-翼型几何中心坐标xr, yr -转换完的翼型离散点坐标2.3叶片各截面翼型参数化设计的Matlab实现基于以上翼型模版和实际翼形处理算法的思想，在Matlab中编写m文件，对叶片所需各截面翼形离散点进行批处理。下面以NACA叶型为例给出 Matlab程序的框架。此段程序是可以同时处理叶片的所有各截面翼形的，各截面的处理后的最终坐标存放在Results单元中。本程序中使用结果保存的Results采用cell格式，可以很好的管理、存放各不同截面翼形的数据，其中每一个小单元都是一个三列矩阵，存放一个截面翼型的（x, y, z）坐

7、标，这对于之后的输出*.dat文件是很有帮助的，得到的截面翼型的（x, y,z）坐标直接复制到记事本中，以截面站位坐标为文件名，保存为.dat文件以供UG导入。NACA= ; ;%翼型模板建立xTop=RTop.*si n( T/180*pi);%初始化参数yTop=RTop.*cos(T/180*pi);C=；%弧长RBottom=RR(i)-N;c=0.4;%型心位置与弧长比xBottom=RBottom.*si n(T./180*pi);theta=;%圆心角yBottom=RBottom.*cos(T./180*pi);r=；%截面站位one=on es(,1);beta=;%安装角R

8、RX(i)=-sig n(alpha(i).*RR(i).*si n( alpha(i)/180deta=0.5.*theta;*pi);alpha=(c-0.5).*theta+beta;RRY(i)=RR(i).*cos(alpha(i)/180*pi);%开始处理Bottom(:,1,i)=xBottom-RRX(i);RR=R(C,theta);%RR中线圆弧半径Bottom(:,2,i)=yBottom-RRY(i);CC=RR.*theta.*pi./180;Bottom(:,3,i)=r(i).*o ne;n=size(r,2);Top(:,1,i)=xTop-RRX(i);Re

9、sults=cell(,);%结果保存单元创建Top(:,2,i)=yTop-RRY(i);for i=1: nTop(:,3,i)=r(i).*o ne;N=CC(i).*NACA(:,2);Results(i,1)=Bottom(:,:,i);T=theta(i).*NACA(:,1)-deta(i)+beta(i);Results(i,2)=Top(:,:,i);RTop=RR(i)+N;end3. 风机叶轮三维几何造型3.1叶片曲面生成把包含各个截面站位的翼形离散点数据的*.dat文件导入UG,就可以用三次样条曲线拟合这些离散点，生成各截面的翼形曲线，然后再通过UG中网格曲面生成的功能

10、，使曲面通过各截面的翼形曲线，这样就可以得到三维的叶片曲面了。具体操作如下：首先，打开UG新建一个部件，然后进入建模模式，选择工具菜单中的插入，然后在下拉菜单下选择曲线中的 样条”在弹出的对话框中选择样条生成方式为通过点”在曲线类型中选择多段”曲线次数为 “3次，之后选择 文件中的点”选项，这时找到你已经按上一 章方式创建好的.dat文件打开，然后在选择约束方式的对话框中，你可以在对翼型斜率或曲 率有特殊的地方，加入斜率或曲率的限制条件，如果没有则直接按确定”这时就生成了一条通过翼形离散点的三次样条曲线。按照此方法，导入一个叶片的各个控制站位截面的翼形离散点的.dat文件，生成各个控制站位截面

11、的翼形曲线如图1所示。然后在插入菜单中选择网格曲面中的通过曲线”选项，然后依次选则各截面曲线，这里选择曲线是有方向性的， 应尽量选择每个曲线的相同部分以保证各个曲线方向一致（如图2），这点要特别注意，不然就生不成高质量的曲面了。 然后在弹出的曲面参数中按照需要定义参数， 这里按默认参数进行生成（如图3）。3.2叶片曲面的检查与修整通常情况下这样构造出来的曲 面由于在过渡部分曲率变化比较剧烈，会出现扭曲现象， 曲面的质量不是很好2。可以使用UG的外观造型模式中使用曲率梳”和 光影分析”工具进行分析检查，在不合要求的曲面区域可以用剪裁片体”的命令将其切除，然后再将这个面补齐并和周边的面进行匹配，以

12、达到要求。在这里如果叶片建模的目的仅仅是用于CFD分析，则不需要十分精确的建模，因为叶片网格数量的限制，网格近似模型的精度远远要粗于我们现 在建模的精度，因此在这里除非曲面有很严重的缺陷，否则不须用过多精力在修改几何模型上。6in 田跑vfr2±LUSElHX fri±.' " rv-v »»r i#in 田跑vfr2±#in 田跑vfr2±图2图3#in 田跑vfr2±3.3叶轮几何模型的生成有了一个修整好的叶片模型 ，则可以通过UG中的特征引用”功能中的 环形矩阵”工具 或使用 编辑”菜单中的 变换”命

13、令手动旋转复制生成其它叶片，然后再对轮毂（Hub）进行几何建模。最后将叶片与轮毂的几何实体进行合”的布尔运算，使其生成一个完整实体，这样就处理好了叶片与轮毂之间的交接线问题，这点对后期的叶轮面网格生成是很重要的，否则叶片面网格与轮毂面网格的衔接会带来很大的麻烦。最后使用 剪裁体”工具中的圆柱面剪裁，将叶顶多余部分剪除，留出叶顶间隙，最后生成风机外罩（shroud），完成叶轮几何建模（如图4）。由于许多划分网格的工具，几何控制能力都不是很强，因此为了后期网格划 分的方便，模型建立时的注意事项： 首先是叶片与轮毂之间的交线处理，划分网格的叶片曲面区域要明确，这之前已经提到；其次是，如果只算一个或几

14、个叶栅的话，则最好留好分割 线；最后，注意叶片与轮毂之间的位置关系要精确，不要出现裂口或突出一块的现象，这些 都将给后期网格划分带来很大麻烦。4 结束语本文介绍了一种使用 Matlab和UG结合起来进行叶片参数化建模的方法，并提出了构 建翼型坐标函数库思想，提出了当翼型中线为圆弧时的翼型坐标算法和各截面安装角和站位 的处理方法，以及它们的 Matlab实现（当然也可使用其它计算软件实现，如Excel等）。通过一个实例，演示了一个叶轮三维几何模型生成的全过程，并提出为了后期网格划分的方便，建模时应注意的事项。希望通过本文的介绍，可以缩短风机叶片建模的时间，提高CFD前处理的质量和效率。图4参考文

15、献1 薛定宇，陈阳泉基于MATLAB/Simulink的系统仿真技术与应用M.清华大学出版社.2004.22 罗强.三维CA系统叶片造型方法研究J.机械设计.2004,21(12):57The parametric modeling method for wing shape vane ofaxial fanMa Jing Wang Zhe n YaCollege of automotive e ngineering, Tongji University, Shang hai, 201804, ChinaAbstract:In the preprocesshg of fan' s CF

16、D analysis, a parametric modeling method for fan' s wing shape vane isprese nted, through creati ng template for wing shape, using Matlab and UG. A process ing means aboutdiscrete coord in ate and sett ing an gel of wing shape is put forward whe n the wing shape' s cen terlian arc, and offers the program in Matlab to carry out the function. Some problems also should be paid attentio