在线工博会

轴流式三元叶片的计算机辅助成型研究
山东工业大学 刘日良 苑国强
为节省流量,手机版未显示文章中的图片,请点击此处浏览网页版
摘要 对轴流式三元叶片成型的关键技术进行了研究,给出了一种从有限个锥形截面出发,利用计算机技术描述叶片空间形状、展开叶片以及生成简便压模型线的方法。
关键词:叶片 三元流动 计算机辅助设计
一、引言
三元流动叶轮接近气体在旋转叶轮中流动的真实情况,因而性能优越,获得越来越广泛的应用。目前制造三元叶轮最简便、经济的工艺方法是钢板焊接,主要过程为:放样展开—板材下料—压制成型—焊接—修型加工。因此,对叶片的空间描述、展开以及压模的设计制造是直接关系着本工艺方法效率和叶轮性能的关键。由于三元叶片本身的复杂性和特殊性,气动设计中一般只给定有限个基元叶栅和叶型的参数,对叶片缺乏统一有效的数学描述,叶片的展开和模具的制造困难。因此,该类风机设计制造的周期长,性能难以保证。基于这种现实,本文给出了一种利用计算机技术解决上述关键问题的一般性方法。
二、叶片型面的坐标计算
1.叶片型面的平面坐标
根据风机的气动设计理论可知,叶片由一系列积迭在一起的叶型面组成。为了便于计算叶栅参数和引用平面叶栅试验数据,叶型是在平面(圆锥面展开)上进行设计的,型面中弧线一般为圆弧形,而轮廓线则由沿中弧线的叶型厚度分布确定。
以型面中弧线的弦长方向作为X1轴方向,Y1轴过中弧线的圆心建立坐标系O1、X1、Y1,如图1。设型面中弧线的半径为r,弦长为B,型面轮廓线上任一点P处的叶片厚度为δ,所在径向线与Y1轴形成的 圆心角为t(顺时针为正),则P点在O1X1Y1中的坐标可表示为

(图片)

图1 描述型面的平面坐标系

(图片)

其中:
(1)对于叶片外表面(吸力面)取正号;内表面(压力面)取负号;
(2)对于薄板型叶片δ取常数(暂不考虑边缘修形);
(3)-T≤t≤T,T=crcsin(B/2R)。
设该型面的形心为O2(x0,y0),安装角为θ,以O2为原点建立与型面安装角无关的坐标系O2X2Y2(X2轴与圆锥面的素线重合),则P点在坐标系O2X2Y2中的坐标为

(图片)

其中,对于图中圆弧板型叶片有

(图片)

2.叶片型面的 空间坐标
上述的叶型是用圆锥面截切叶片,将截面展开后得到的平面图形。为了得到真实形状,需将截面卷回到圆锥面上。设圆锥面的端面半径为R1、R2(R1<R2),长度为L,建立如图2所示的右手坐标系O3X3Y3Z3,其中,X3轴与圆锥轴线重合,Z3轴与过积迭点的径向线重合。则有如下的映射关系。

(图片)

图2 空间坐标系的建立

(图片)

其中,z0为积迭点的高度,

(图片)

三、叶片的成型
1.叶片曲面的形成
根据公式(1)、(2)、(3),对于气动设计中给定的任一叶型面,确定了积迭线(一般采用径向线)后便可以准确地求得空间位置。型面轮廓线上任意点的空间位置可表示为
Pi(x,y,z)=P(r1±δ/2,t1,z01) (i=1,2,…,n) (4)
将上述n个截面积迭起来,便得到整个叶片曲面的骨架。为了得到光滑的叶片,还需要进行曲面拟合。下面是采用三次B样条曲线拟合叶片的具体方法。
(1)求型值点Pi
对于截面i(i=1,2,…,n),令ti=wTi(-1≤w≤1),固定w,根据公式(4),得到叶片上处于不同截面上的n个对应点P1,P2,…,Pn。
(2)构造三次B样条曲线,反求控制点Qi
过上述点构造一条三次B样条曲线,表示为

(图片)

其中 i=1,2,…,n,令u=0,则有

(图片)

令Pi(0)=Pi,Q0=Q1,Qn+1=Qn,得

(图片)

至此,用追赶法可解得Q1,Q2,…Qn+1。
(3)生成曲面
令w连续变化,即控制P1,P2,…,Pn沿截面连续移动则得到光滑曲面,以[B]表示三次B样条系数矩阵,可记为
Pi(u,w)=[u3 u2 u 1][B]
[Qi-1(w) Qi(w) Qi+1(w) Qi+2(w)]T (6)
式中,0≤u≤1,-1≤w≤1,i=1,2,…,n-1。
2.叶片的展开
三元叶片属于不可展曲面,只能进行近似的展开。令δ=1,根据公式(4)、(5)、(6)可得到叶片中面,将之划分为若干微小 单元后可求得任意曲面单元的四个顶点p(u,w),p(u,w+Δw),P(u+Δu,w),p(u+Δu,w+Δw)。至此利用三角线法可得到较为准确的叶片展开图。值得一提的是,除叶片的构型精度和划分方式及曲面单元的数量以外,展开基点的选择对展开精度的影响也很大。习惯上沿一端顺次展开的方法会产生较大的误差。建议以根截面(i=1)的中点(t=0)为基点。
3.平面型线的生成
(1)任意平面型线的求解
对于式(6)确定的叶片曲面,与任意平面ax+by+cz+d=0联立求交,可得交线的方程为
[u3 u2 u 1][B](α[X(w)]T
+b[Y(w)]T+c[Z(w)]T)+d=0 (7)
其中,[X(w)]、[Y(w)]、[Z(w)]为控制矩阵[Q(w)]的坐标分量。
上式中,给定w=w*,得到关于u的三次方程,解方程后得u*,则(u*,w*)对应着一个交点(B样条曲线与平面的交点)。按合适的步长,连续求得一系列交点,则得到需要的平面型线。
(2)等x值型线的简便求解方法
将坐标系O3X3Y3Z3绕Z3轴旋转β角(取β约等于叶弦角的平均值)得新坐标系O4X4Y4Z4,如图3a,将式(4)表示的各叶型截面转换到该坐标系中,得

(图片)

如图3 b所示,式(4)表示的曲线以及边缘线(样条曲线),构成了叶片曲面的基本骨架,。取特殊位置平面x4=x*与之求交,可得叶片上的一系列点A1,A2,…,An;B1,B2,…,Bn。分别以A1,A2,…,An和B1,B2,…,Bn作为型值点可生成与上、下压模对应的两条型线。根据制造精度和工艺的要求,取一系列合适的x*,则得到一组简便的压模型线。

(图片)

图3 坐标转换及型线求解

四、结束语
根据气动设计提供的若干扭变截面,采用合适的数学模型拟合叶片曲面的方法简便易行,叶片造型较为准确。在此基础上生成的展开图和平面型线图可作为板材下料和设计压模时的重要依据。经过软件编制和样机试制证明,本方法简便可靠,叶片成型准确,焊接及修型加工方便,是提高该类风机生产效率和性能的有效措施。
6/5/2004


电脑版 客户端 关于我们
佳工机电网 - 机电行业首选网站