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

原创 2020-07-21 09:57  阅读

  轴流风机机翼型叶片参数化建模方法_专业资料。本文通过创建翼型模板,结合Matlab与UG软件,探讨了风机翼型叶片参数化建模的方法,给出了翼型中线为圆弧时的翼型坐标算法、各截面安装角和站位的处理方法以及Matlab实现程序。并提出了叶片在UG建模时应注意的问题。文中提出的方法,减少了风机建模的工作量,缩短了风机CFD前处

  轴流风机机翼型叶片参数化建模方法 马静 王振亚 同济大学汽车学院 上海 (201804) Email: 摘 要:本文通过创建翼型模板,结合 Matlab 与 UG 软件,探讨了风机翼型叶片参数化建模 的方法,给出了翼型中线为圆弧时的翼型坐标算法、各截面安装角和站位的处理方法以及 Matlab 实现程序。并提出了叶片在 UG 建模时应注意的问题。文中提出的方法,减少了风机 建模的工作量,缩短了风机 CFD 前处理周期,提高了风机流场 CFD 分析计算的效率和质量。 关键词:叶片;参数化设计;UG;Matlab 1. 前言 随着CFD技术的迅速发展,对风机流场计算分析的要求越来越多。风机仿真计算的前期 工作量相当大,主要表现在机翼型叶片的建模,其中包括风机叶轮的机翼型叶片,机翼型前 导流叶片和叶轮后的止旋片建模。通常在UG软件中输入大量的翼型坐标点是相当麻烦的, 而使用*.dat文件导入这些数据的方法要方便的多, 但是对不同的叶片计算截面采用*.dat文件 手工导入翼型坐标点的工作量仍然非常大,并且修改起来也不方便。通过分析可知,叶片不 同计算截面的翼型曲线是相似的, 同种翼型只因弧长以及中线形状不同而不同, 因此完全可 以考虑采用参数化建模的设计方法。采用这种方法可以缩短建模时间,节省大量的工作量, 且所建的模型也易于修改。 因为在对风机流场进行CFD分析计算时改变风机叶片翼型是对风 机模型的重大修改需要花费大量的时间, 有了这种方法可以较轻松的完成修改。 本文就是基 于这种思想,介绍了用Matlab与UG两个软件结合进行风机叶片参数化建模的方法,本方法 利用Matlab强大的数据处理能力处理翼型离散点[1], 用UG强大的三维曲面建模能力构建叶片 复杂曲面。 2. 翼型离散点的参数化处理 2.1 翼型模板的建立 翼型模板的建立是实现参数化设计的第一步, 建立翼型模板库是一个积累的过程, 需要 将每次用到的翼型和收集到的有价值的翼型参数通过手工输入, 建立起翼型模板库, 在进行 风机叶片建模时就可以非常方便的从翼型模板库里直接调出所需要的翼型。 在 Matlab 中可以通过一个两列矩阵建立起翼型模板, 第一列输入原始翼型的 x / l 值, 第 1 二列输入原始翼型的 y / l 值,再存成以翼型名称命名的.mat 文件即可。有了翼型模板,在进 行叶片建模时就可以根据各计算截面的翼型弦长,直接得到原始翼型各离散点的坐标值。 2.2 实际翼型的处理 实际翼型有很多种,大部分都是在原始翼型的基础上通过改变翼型中线的形状得到的, 其中最常见的是中线为圆弧或多段圆弧组合的翼型, 因此这里以圆弧翼型为例, 给出翼型坐 标点的变换算法及 Matlab 程序,多段圆弧可作相似处理。 此算法有三个目的: 1. 使原始翼型变为中线. 使翼型的几何中心与坐标原点重合,以便于后期叶片拉伸建模方便。 3. 使翼型的安装角与设计相同,以便于后期叶片建模。 算法思想:先在圆柱坐标系下进行变换,以中线圆弧圆心为原点,用半径 r 和角度 θ 坐标确 定翼型各离散点位置,再转换到笛卡尔坐标系中,以便于 UG 导入。 算法实现: 1. 将原始翼型的 X 坐标转换为圆柱坐标系下的 θ 坐标: θ =α × 其中: α --翼型中线圆弧所对应的圆心角 X l l --翼型中线. 将原始翼型的 Y 坐标转换为圆柱坐标系下的 r 坐标: r = R +Y 其中: R --翼型中线. 再将翼型安装角转为零: θ1 = θ ? 4. 将翼型转至需要的安装角 β : α。 2 θ 2 = θ1 + β 5. 将圆柱坐标系下坐标 ( r , θ ) 转化到相同原点的笛卡尔坐标系下的坐标 ( x, y ) : x = r × cos(θ 2 ) ; y = r × sin(θ 2 ) 6. 最后将翼型几何中心移到坐标原点: 2 其中: xc , yc --翼型几何中心坐标 xr , yr --转换完的翼型离散点坐标 2.3 叶片各截面翼型参数化设计的Matlab实现 基于以上翼型模版和实际翼形处理算法的思想,在 Matlab 中编写 m 文件,对叶片所需 各截面翼形离散点进行批处理。下面以 NACA 叶型为例给出 Matlab 程序的框架。此段程序 是可以同时处理叶片的所有各截面翼形的,各截面的处理后的最终坐标存放在 Results 单元 中。本程序中使用结果保存的 Results 采用 cell 格式,可以很好的管理、存放各不同截面翼 形的数据,其中每一个小单元都是一个三列矩阵,存放一个截面翼型的 ( x, y , z ) 坐标,这对 于之后的输出*.dat 文件是很有帮助的,得到的截面翼型的 ( x, y , z ) 坐标直接复制到记事本 中,以截面站位坐标为文件名,保存为.dat 文件以供 UG 导入。 NACA=[ ; ]; %翼型模板建立 %初始化参数 C=[ ]; %弧长 c=0.4; %型心位置与弧长比 theta=[ ]; %圆心角 r=[ ]; %截面站位 beta=[ ]; %安装角 deta=0.5.*theta; alpha=(c-0.5).*theta+beta; %开始处理 RR=R(C,theta); %RR 中线圆弧半径 CC=RR.*theta.*pi./180; n=size(r,2); Results=cell(,); %结果保存单元创建 for i=1:n N=CC(i).*NACA(:,2); T=theta(i).*NACA(:,1)-deta(i)+beta(i); RTop=RR(i)+N; xTop=RTop.*sin(T/180*pi); yTop=RTop.*cos(T/180*pi); RBottom=RR(i)-N; xBottom=RBottom.*sin(T./180*pi); yBottom=RBottom.*cos(T./180*pi); one=ones(,1); RRX(i)=-sign(alpha(i)).*RR(i).*sin(alpha(i)/180 *pi); RRY(i)=RR(i).*cos(alpha(i)/180*pi); Bottom(:,1,i)=xBottom-RRX(i); Bottom(:,2,i)=yBottom-RRY(i); Bottom(:,3,i)=r(i).*one; Top(:,1,i)=xTop-RRX(i); Top(:,2,i)=yTop-RRY(i); Top(:,3,i)=r(i).*one; Results(i,1)={Bottom(:,:,i)}; Results(i,2)={Top(:,:,i)}; end 3.风机叶轮三维几何造型 3.1 叶片曲面生成 把包含各个截面站位的翼形离散点数据的*.dat 文件导入 UG, 就可以用三次样条曲线 合这些离散点,生成各截面的翼形曲线,然后再通过 UG 中网格曲面生成的功能,使曲面通 过各截面的翼形曲线,这样就可以得到三维的叶片曲面了。具体操作如下: 首先,打开 UG 新建一个部件,然后进入建模模式,选择工具菜单中的插入,然后在下 拉菜单下选择曲线中的“样条”,在弹出的对话框中选择样条生成方式为“通过点”,在曲线类 型中选择“多段”,曲线”次,之后选择“文件中的点”选项,这时找到你已经按上一 章方式创建好的.dat 文件打开,然后在选择约束方式的对话框中,你可以在对翼型斜率或曲 率有特殊的地方,加入斜率或曲率的限制条件,如果没有则直接按“确定”,这时就生成了一 条通过翼形离散点的三次样条曲线。 按照此方法, 导入一个叶片的各个控制站位截面的翼形 离散点的.dat 文件,生成各个控制站位截面的翼形曲线 然后在插入菜单中选择网格曲面中的“通过曲线”选项,然后依次选则各截面曲线,这里选择 曲线是有方向性的,应尽量选择每个曲线的相同部分以保证各个曲线) ,这 点要特别注意, 不然就生不成高质量的曲面了。 然后在弹出的曲面参数中按照需要定义参数, 这里按默认参数进行生成(如图3) 。 3.2 叶片曲面的检查与修整 通常情况下这样构造出来的曲面由于在过渡部分曲率变化比较剧烈,会出现扭曲现象, 可以使用UG的外观造型模式中使用“曲率梳”和“光影分析”工具进行 曲面的质量不是很好[2]。凯旋门娱乐国际, 分析检查,在不合要求的曲面区域可以用“剪裁片体”的命令将其切除,然后再将这个面补齐 并和周边的面进行匹配,以达到要求。在这里如果叶片建模的目的仅仅是用于CFD分析,则 不需要十分精确的建模, 因为叶片网格数量的限制, 网格近似模型的精度远远要粗于我们现 在建模的精度, 因此在这里除非曲面有很严重的缺陷, 否则不须用过多精力在修改几何模型 上。 4 图2 图3 3.3 叶轮几何模型的生成 有了一个修整好的叶片模型,则可以通过 UG 中的“特征引用”功能中的“环形矩阵”工具 或使用“编辑”菜单中的“变换”命令手动旋转复制生成其它叶片,然后再对轮毂(Hub)进行 几何建模。最后将叶片与轮毂的几何实体进行“合”的布尔运算,使其生成一个完整实体,这 样就处理好了叶片与轮毂之间的交接线问题, 这点对后期的叶轮面网格生成是很重要的, 否 则叶片面网格与轮毂面网格的衔接会带来很大的麻烦。 最后使用“剪裁体”工具中的圆柱面剪 裁,将叶顶多余部分剪除,留出叶顶间隙,最后生成风机外罩(shroud) ,完成叶轮几何建 模(如图4) 。由于许多划分网格的工具,几何控制能力都不是很强,因此为了后期网格划 分的方便,模型建立时的注意事项:首先是叶片与轮毂之间的交线处理,划分网格的叶片曲 面区域要明确,这之前已经提到;其次是,如果只算一个或几个叶栅的话,则最好留好分割 线;最后,注意叶片与轮毂之间的位置关系要精确,不要出现裂口或突出一块的现象,这些 都将给后期网格划分带来很烦。 4.结束语 本文介绍了一种使用 Matlab 和 UG 结合起来进行叶片参数化建模的方法,并提出了构 建翼型坐标函数库思想, 提出了当翼型中线为圆弧时的翼型坐标算法和各截面安装角和站位 的处理方法,以及它们的 Matlab 实现(当然也可使用其它计算软件实现,如 Excel 等) 。通 过一个实例, 演示了一个叶轮三维几何模型生成的全过程, 并提出为了后期网格划分的方便, 建模时应注意的事项。希望通过本文的介绍,可以缩短风机叶片建模的时间,提高 CFD 前 处理的质量和效率。 5 图4 参考文献 [1] 薛定宇,陈阳泉 基于MATLAB/Simulink的系统仿真技术与应用[M].清华大学出版社.2004.2 [2] 罗强. 三维CAD系统叶片造型方法研究[J].机械设计.2004,21(12):57 The parametric modeling method for wing shape vane of axial fan Ma Jing Wang Zhen Ya College of automotive engineering, Tongji University, Shanghai, 201804, China Abstract: In the preprocessing of fan’s CFD analysis, a parametric modeling method for fan’s wing shape vane is presented, through creating template for wing shape, using Matlab and UG. A processing means about discrete coordinate and setting angel of wing shape is put forward when the wing shape’s centerline is an arc, and offers the program in Matlab to carry out the function. Some problems also should be paid attention to, when using UG to create model. The method which is put forward by this paper, reduces working time of the preprocessing of fan’s CFD analysis, and improves the efficiency and quality of preprocessing. Key words: Vane; Parametric design; UG; Matlab 6

版权声明:本文为原创文章,版权归 凯旋门娱乐国际 所有,欢迎分享本文,转载请保留出处!
上一篇:机翼型管道轴流风机 YFIAD
下一篇:凯旋门娱乐国际上海市质监局抽查28批次风机产品