叶片曲面整体曲率光顺方法
【技术领域】
[0001] 本发明涉及微分几何以及叶轮叶片曲面造型方法，具体地，涉及一种叶片曲面整 体曲率光顺方法。
【背景技术】
[0002] 叶轮类零件在航空航天、能源动力等行业中均有广泛应用。叶片的三维造型是叶 轮设计的难点，叶片曲面一般设计为自由曲面，型面结构非常复杂，而且叶片曲面需要满足 几何连续性的要求，因为叶片曲面的光顺性对叶轮本身的效率有很大的影响。
[0003] 目前叶片曲面的设计通常是先设计出叶片各流面上的流线，然后积叠流线得到叶 片曲面。流线往往可以设计得很光顺，但由于没有对叶片叶高方向的光顺性进行控制，如曲 率或曲率变化的控制，叶片沿叶高方向可能会出现皱褶，不光顺的情况。另外，叶片曲面的 光顺性也会影响到后续的叶片数控加工工艺，如曲面不光顺会造成加工刀轨计算失败，或 者生成的刀轨不光滑，造成加工过程中机床的振动。
【发明内容】
[0004] 针对现有技术中的缺陷，本发明的目的是提供一种叶片曲面整体曲率光顺方法。
[0005] 根据本发明提供的叶片曲面整体曲率光顺方法，包括如下步骤：
[0006] 步骤Sl:输入初始流线，插值初始流线得到叶片曲面；
[0007] 步骤S2:计算叶片曲面在曲面采样点处沿流线方向和叶展方向的法曲率导数；
[0008] 步骤S3:计算叶片曲面沿流线方向和叶展方向的法曲率导数平方总和E;
[0009] 步骤S4:针对叶片曲面的点集{pkGR3，I < k < m}中的每个点，计算每个点到原叶片 曲面的距离d;其中，pk为点集中第k个点，R3为3维实向量空间，m为点集中点总个数；
[0010] 步骤S5:若距离d大于给定的允许误差值dmax，则不调整该点，并在当前迭代次数未 超过最大迭代次数T max时返回步骤S4以进行迭代计算；若距离d小于给定的允许误差值 dmax，则利用黄金分割法搜索该点沿叶片圆周方向的最优调整角度Θ，使叶片的法曲率导数 平方总和E取最小值，并记SE new;
[0011] 若Ene3W<E，则计算调整后该点的坐标，并替换原来的点，令E = Ene3w，并在当前迭代 次数未超过最大迭代次数Tmax时返回步骤S4以进行迭代计算;否则不调整该点，并在当前迭 代次数未超过最大迭代次数T max时返回步骤S4以进行迭代计算。
[0012] 优选地，所述步骤2包括如下步骤：
[0013]步骤S201:根据欧拉公式，叶片曲面在点P处沿切向X的法曲率计算公式Kp(X)为
[0015] 其中K1表示曲面在点p处的最大主曲率，K2表示曲面在点p处最小主曲率，P表示切 向X与主曲率^所在的主曲率方向的有向夹角；
[0016] 步骤S202:由于叶片曲面是通过插值初始流线得到的，则认为叶片曲面的u参数线 是叶片曲面的流线，V参数线方向是叶片的叶高方向；
[0017]步骤5203:设?1，风两点是在叶片曲面的同一1!参数线上，且两点的弧长距离为 8，叶 片曲面在点?1处沿流线方向的法曲率导数的计算公式为，
[0019] 其中tl表示叶片曲面在点P1处的流线方向，t2表示叶片曲面在点? 2处的流线方向， ⑷表示叶片曲面在点Pi处沿tI方向的法曲率，\ )表示叶片曲面在点P2处沿t2方向的 法曲率。
[0020] 优选地，所述步骤3具体为：
[0021] 根据沿叶片曲面流线方向和叶高方向给定的采样点数目，首先计算叶片在每点处 沿流线方向和叶高方向的法曲率导数，进而计算叶片曲面沿流线方向和叶展方向的法曲率 导数平方总和E，
[0023]
表示叶片在点？1处沿流线方向的法曲率导数，
P1处沿叶高方向的法曲率导数;w i，W 2分别表示叶片沿流线方向、叶高方向的法曲率导数平 方总和权重;η为曲面采样点个数。
[0024]优选地，所述步骤S5中叶片曲面上点的调整方向被限定在圆周方向上，保证调整 后的点与原来的点处于相同的流面上。
[0025]优选地，所述步骤5包括如下内容：
[0026]步骤S501:给定所述点的最大调整距离dmax，根据数据点的回转半径计算该点的最 大调整角度Qmax;在区间[_9max，0max]中，利用黄金分割法搜索最优调整角度Θ，使得叶片曲面 的法曲率导数平方总和E最小。
[0027] 与现有技术相比，本发明具有如下的有益效果：
[0028] 1、本发明解决了同时光顺优化叶片曲面沿流线方向和叶高方向的法曲率导数的 问题，不仅能够使叶片曲面沿流线方向光顺，而且能够使叶片曲面沿叶高方向也光顺，使叶 片曲面沿这两个方向的法曲率变化均匀，曲面变光滑；
[0029] 2、本方法适用于轴流式和离心式叶轮叶片曲面的整体曲率光顺造型方法。
【附图说明】
[0030] 通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、 目的和优点将会变得更明显：
[0031] 图1为本发明的流程图；
[0032] 图2为本发明中离心式叶轮叶片曲面在子午面上的流线模型示意图；
[0033] 图3为本发明中计算叶片曲面在采样点处沿流线方向和叶高方向法曲率导数示意 图；
[0034] 图4为本发明中沿圆周方向调整叶片曲面数据点示意图；
[0035] 图5为本发明中优化后的叶片曲面示意图。
[0036] 图中：
[0037] 1为在子午面上的叶片曲面流线；
[0038] 2为叶片曲面上的采样点；
[0039] 3为叶片曲面在同一流线上与采样点的弧长距离为s的点；
[0040] 4为叶片曲面沿叶高方向与采样点的弧长距离为s的点。
【具体实施方式】
[0041] 下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术 人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技 术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进。这些都属于本发 明的保护范围。
[0042] 在本实施例中，如图1、图2、图3、图4以及图5所示，本发明利用整体优化法使叶片 曲面沿流线方向和叶高方向的法曲率导数平方总和减少，使曲面变光顺。
[0043] 输入初始流线，通过插值初始流线得到叶片曲面，根据欧拉公式计算叶片曲面在 采样点P处沿流线方向和叶高方向的法曲率，
[0045] 其中K1,K2表示叶片曲面在点p处的主曲率，P表示切向X与主曲率^所在的主曲率 方向的有向夹角。
[0046] 计算叶片曲面在采样点处沿流线方向和叶高方向的法曲率导数。计算叶片曲面沿 流线方向的法曲率导数公式如下，
[0048]其中Pl，p2两点是在叶片曲面的同一 u参数线上，且两点的弧长距离为s，s是一个很 小的数。tl表示叶片曲面在点Pl处的流线方向，t2表示叶片曲面在点P2处的流线方向。计算 叶片沿叶高方向的法曲率导数的方法与计算沿流线方向的类似。
[0049]根据给定的沿叶片曲面流线方向和叶高方向的采样点个数，计算叶片沿这两个方 向的法曲率导数平方总和，
[0051]
表示叶片在点？1处沿流线方向的法曲率导数，
表示