[发明专利]基于单纯形样条函数的航空发动机稳态模型的建模方法在审

专利信息
申请号: 201510251444.7 申请日: 2015-05-15
公开(公告)号: CN104834785A 公开(公告)日: 2015-08-12
发明(设计)人: 郑前钢;华伟;李永进;孙丰勇;叶志锋;张海波 申请(专利权)人: 南京航空航天大学
主分类号: G06F17/50 分类号: G06F17/50
代理公司: 南京经纬专利商标代理有限公司 32200 代理人: 许方
地址: 210016 江*** 国省代码: 江苏;32
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 基于 单纯 形样条 函数 航空发动机 稳态 模型 建模 方法
【权利要求书】:

1.一种基于单纯形样条函数的航空发动机稳态模型的建模方法,根据给定发动机的飞行参数和控制参数,对航空发动机稳定状态进行实时估计,建立航空发动机的稳态模型,其特征在于包括以下步骤:

步骤一、获取航空发动机稳态模型的训练和测试数据;

步骤二、进行发动机工作区域的三角划分,并计算工作点所在单纯形的重心坐标;

步骤三、计算单纯形样条的基函数;

步骤四、求解单纯形样条函数的系数;

步骤五、利用所述的测试数据对模型泛化能力进行验证,若精度差则返回上述步骤二,重新进行三角划分,并调节单纯形样条的基函数的阶数;若精度好,则建立航空发动机稳态模型。

2.根据权利要求1所述的基于单纯形样条函数的航空发动机稳态模型的建模方法,其特征在于:在所述步骤一中,所述的获取航空发动机稳态模型的训练和测试数据是指:根据所述发动机的飞行参数和控制参数,进行发动机试车试验或根据航空发动机部件模型得到数据。

3.根据权利要求1所述的基于单纯形样条函数的航空发动机稳态模型的建模方法,其特征在于:

在所述步骤二中,所述的发动机工作区域的三角划分采用Delaunay三角划分算法进行;

所述的计算工作点所在单纯形的重心坐标是指单纯形内的点x对应单纯形的重心坐标:b(x)=(b0,b1,…,bn),其中,bi为x对应单纯形顶点的重心坐标,因此,x可以表示为:

x=Σi=0nbivpi]]>

其中,为单纯形的顶点,pi为顶点指数的排序即pi<pi+1

4.根据权利要求1所述的基于单纯形样条函数的航空发动机稳态模型的建模方法,其特征在于:在所述步骤三中,所述的求解单纯形样条函数的系数,其计算公式如下:

Bkd(b)=d!k!bk]]>

其中,k称为多维系k=(k0,k1,…,kn)∈Nn+1,k!为多维系数的阶乘积k!=k0!k1!…kn!,b为单纯形样条的重心坐标,bk等于d表示基函数多项式的阶数。

5.根据权利要求1所述的基于单纯形样条函数的航空发动机稳态模型的建模方法,其特征在于:在所述步骤四中,所述的单纯形样条函数的系数的求解为线性回归问题:

y(i)=Σj=1JΣ|k|=dcktjBkd(b)+r(i)]]>

式中,b(i)为发动机工作点所对应的单纯形区域tj的重心坐标,为单纯形样条的基函数,单纯形样条函数的系数的求解可进一步转换为求解残差平方和最小的最优化问题即

min J=(Y-Xc)T(Y-Xc)

Hc=0

式中,c为单纯形样条函数的系数,X为单纯形样条的基函数,Y为航空发动机输出数据,J=(Y-Xc)T(Y-Xc)为广义最小二乘的目标函数,Hc=0为单纯形样条函数的光滑条件,上式称为带有等式约束的广义最小二乘问题;对于上述问题,运用格朗日乘子得到Karush-Kuhn-Tucher条件:

XTXHTH0cv=XTY0]]>

式中,v为拉格朗日乘子,最终把单纯形样条函数的系数求解问题转化为线性方程的求解。

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于南京航空航天大学,未经南京航空航天大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/201510251444.7/1.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top