• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 20
    • 下载费用:30 金币  

    怎样玩好重庆时时彩: 一种涡轴发动机动态实时模型建模方法及系统.pdf

    关 键 词:
    一种 发动机 动态 实时 模型 建模 方法 系统
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201410050351.3

    申请日:

    2014.02.13

    公开号:

    CN103838920A

    公开日:

    2014.06.04

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G06F 17/50申请日:20140213|||公开
    IPC分类号: G06F17/50 主分类号: G06F17/50
    申请人: 南京航空航天大学
    发明人: 鲁峰; 黄金泉; 张冬冬; 姚文荣; 邢耀东
    地址: 210016 江苏省南京市秦淮区御道街29号
    优先权:
    专利代理机构: 江苏圣典律师事务所 32237 代理人: 贺翔
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201410050351.3

    授权公告号:

    ||||||

    法律状态公告日:

    2016.08.17|||2014.07.02|||2014.06.04

    法律状态类型:

    授权|||实质审查的生效|||公开

    摘要

    本发明公开了一种涡轴发动机动态实时模型建模方法及系统,该方法包括:建立涡轴发动机稳态实时模型;采用剩余燃油流量动态加减系数法建立涡轴发动机动态实时模型。本发明解决了现有的发动机动态实时模型建模过程中在只有试车数据而缺乏部件特性数据情况下的建模问题,适用于发动机数控系统半物理仿真,亦可为传感器解析余度技术提供模型保证,对于缩短研制周期、降低试验风险及成本有着积极促进作用。

    权利要求书

    权利要求书
    1.  一种涡轴发动机动态实时模型建模方法,其特征在于,包括以下步骤:
    步骤A),建立涡轴发动机稳态实时模型;
    步骤B),采用剩余燃油流量动态加减速系数法建立涡轴发动机动态实时模型。

    2.  根据权利要求1所述的涡轴发动机动态实时模型建模方法,其特征在于,步骤A)的具体步骤如下:
    步骤A1),根据涡轴发动机工作条件,将试车数据进行相似归一化的换算;
    步骤A2),对换算后的数据进行野点剔除及平滑滤波处理;
    步骤A3),根据试车数据中稳态过程的数据,对工作状态相近的稳态数据进行聚类合并,以此构造稳态工作点;
    步骤A4),根据稳态工作点,采用多项式拟合法获得稳态关系曲线,利用稳态关系曲线插值获得稳态基点,并以此构造稳态基点插值表;
    步骤A5),根据构造的稳态基点插值表,采用线性插值法建立包含起动机在内的涡轴发动机的稳态实时模型。

    3.  如权利要求2所述的涡轴发动机动态实时模型建模方法,其特征在于,步骤A4)中利用稳态关系曲线插值获得稳态基点的详细步骤如下:
    将燃气涡轮转速换算值与预先设定的转速换算阈值进行比较,如果燃气涡轮转速换算值大于等于预先设定的转速换算阈值,将稳态关系曲线以约2%的间隔进行插值获得稳态基点;果燃气涡轮转速换算值小于预先设定的转速换算阈值,将稳态关系曲线以3%-5%的间隔进行插值获得稳态基点。

    4.  如权利要求2所述的涡轴发动机动态实时模型建模方法,其特征在于,步骤B)的具体步骤如下:
    步骤B1),根据试车数据中动态过程的数据,利用剩余燃油量动态系数法的定义,分别计算出不同燃气涡轮转子转速换算值下的燃气涡轮转速动态系数、燃气涡轮出口温度动态系数、压气机出口压力动态系数和动力涡轮输出功率动态系数,以获得动态系数插值表;
    步骤B2),建立慢车以上状态燃气发生器动态实时模型,该模型的输入参数为燃油流量、进口总压、进口总温,输出参数为燃气涡轮转速、压气机出口压力、燃气涡轮出口温度,采用剩余燃油流量动态系数法计算出模型输出参数;
    步骤B3),建立慢车以上状态动力涡轮及负载动态实时模型,该模型输入参数为燃气涡轮转速、燃油流量、进口总压、进口总温、负载杆角度,模型输出参数为动力涡轮转速、动力涡轮输出功率,分别利用剩余燃油流量动态系数法和动力涡轮转子功率平衡求解模型输出参数;
    步骤B4),建立涡轴发动机起动阶段数学模型,该模型包括起动机单独带转阶段的模型、起动机与燃气涡轮共同带转阶段的模型、以及燃气涡轮单独带转阶段的模型;
    步骤B5),建立涡轴发动机动态实时模型,该模型包含步骤A5) 建立的稳态实时模型、步骤B2) 建立的慢车以上燃气发生器动态实时模型、步骤B3)建立的动力涡轮及负载动态实时模型、以及步骤B4)建立的起动阶段数学模型,根据各个模型的输入参数解算各个模型的输出参数。

    5.  如权利要求4所述涡轴发动机动态实时模型建模方法,其特征在于,所述动态系数包含加速系数和减速系数。

    6.  如权利要求4所述涡轴发动机动态实时模型建模方法,其特征在于,步骤B2)中采用剩余燃油流量动态系数法计算出模型输出参数时,燃气涡轮转速以前一时刻的动态值为当前时刻的初始值进行累加,压气机出口压力和燃气涡轮出口温度以当前时刻的稳态值作为初始值进行累加。

    7.  如权利要求4所述涡轴发动机动态实时模型建模方法,其特征在于,步骤B3)中采用剩余燃油流量动态系数法计算出模型输出参数时,动力涡轮输出功率以当前时刻的稳态值作为初始值进行累加。

    8.  如权利要求4所述涡轴发动机动态实时模型建模方法,其特征在于,步骤B2)中采用剩余燃油流量动态系数法计算出燃气涡轮出口温度后,作为参数输入至一阶惯性温度传感器模型,以获得能够反映延迟效应的出燃气涡轮出口温度输出参数。

    9.  基于权利要求4所述涡轴发动机动态实时模型建模方法的系统,其特征在于,包括试车数据预处理???、稳态工作区域判定逻辑???、动态工作区域判定逻辑???、稳态基点自动生成???、动态加减速系数自动生成??楹投凳蹦P图扑阍诵心?椋?BR>所述试车数据预处理??橛糜诙栽际猿凳萁刑蕹暗愫推交瞬ù?;
    所述稳态工作区域判定逻辑??橛糜诟菥鍪猿凳菰ご砟?榇砗蟮氖猿凳葜形忍痰氖?,获得稳态工作点数据;
    所述动态工作区域判定逻辑??橛糜谂卸ň鍪猿凳菰ご砟?榇砗蟮氖猿凳葜械亩ぷ鞣段?,获得不同燃气涡轮转子转速换算值下燃气涡轮转速动态系数、燃气涡轮出口温度动态系数、压气机出口压力动态系数以及动力涡轮输出功率的加速系数与减速系数;
    所述稳态基点自动生成??橛糜诟菸忍ぷ髑蚺卸呒?橹猩傻奈忍ぷ鞯闶?、采用聚类合并法和拟合插值法计算出稳态基点,生成稳态基点插值表;
    所述动态加减速系数自动生成??橛糜诟菟龆ぷ髑蚺卸呒?橹猩傻牟煌计新肿幼倩凰阒迪氯计新肿俣凳?、燃气涡轮出口温度动态系数、压气机出口压力动态系数以及动力涡轮输出功率的加速、减速系数,采用拟合插值法获得其在稳态基点对应下的加速、减速系数,生成动态加减速系数插值表;
    所述动态实时模型计算运行??橛糜诟菟鑫忍阕远赡?樯傻奈忍悴逯当砗退龆蛹跛傧凳远赡?橹猩傻亩蛹跛傧凳逯当?,采用动态加减速系数法、根据涡轮发动机动态实时模型的输入参数计算出发动机动态实时模型的输出参数。

    说明书

    说明书一种涡轴发动机动态实时模型建模方法及系统
    技术领域
    本发明属于航空发动机建模与仿真领域,具体涉及一种涡轴发动机动态实时模型建模方法及系统。
    背景技术
    控制用传感器发生故障将直接导致闭环控制系统反馈信息不可用,使得涡轴发动机不能稳定工作,因此,涡轴发动机数控系统的传感器能否稳定可靠工作,对于保证涡轴发动机FADEC系统正常工作至关重要。传感器解析余度技术是一种典型的基于模型的传感器故障诊断方法,涡轴发动机动态实时模型是该方法的重要组成部分,同时该模型的精度决定了传感器故障诊断性能,因此,建立某型涡轴发动机动态实时模型对传感器故障诊断和容错控制具有重要的理论意义和工程实用价值。
    国内外用于发动机控制系统半物理试验的动态实时模型大都采用部件级建模的方法,但是当缺乏发动机部件级建模用的部件特性数据时,只能根据试验数据建模,同时部件级模型由于计算耗时长、占用资源多的特点,使得其应用范围受到限制,不能满足机载需求。因此,需要研究根据试验数据直接建立机载动态实时模型的方法。目前发动机动态实时模型建模方法主要集中在依赖部件特性数据的部件级建模方法上,有限的动态系数方法研究没有考虑加、减速的系数差异,以及没有在更为复杂的涡轴发动机上应用,更没有具有动态实时模型自动生成的通用软件开发的相关报道。
    发明内容
    本发明所要解决的技术问题是针对背景技术的缺陷,提供一种涡轴发动机动态实时模型建模方法及系统,用于解决缺乏部件特性数据、只提供试验数据时传统的涡轴发动机动态实时模型机载实时模型建模问题,该方法和系统满足机载实时性要求,能为传感器解析余度技术提供模型保证。
    本发明为解决上述技术问题采用以下技术方案:
    一种涡轴发动机动态实时模型建模方法,包括以下步骤:
    步骤A),建立涡轴发动机稳态实时模型;
    步骤B),采用剩余燃油流量动态加减速系数法建立涡轴发动机动态实时模型。
    作为本发明一种涡轴发动机动态实时模型建模方法进一步的优化方案,步骤A)的具体步骤如下:
    步骤A1),根据涡轴发动机工作条件,将试车数据进行相似归一化的换算;
    步骤A2),对换算后的数据进行野点剔除及平滑滤波处理;
    步骤A3),根据试车数据中稳态过程的数据,对工作状态相近的稳态数据进行聚类合并,以此构造稳态工作点;
    步骤A4),根据稳态工作点,采用多项式拟合法获得稳态关系曲线,利用稳态关系曲线插值获得稳态基点,并以此构造稳态基点插值表;
    步骤A5),根据构造的稳态基点插值表,采用线性插值法建立包含起动机在内的涡轴发动机的稳态实时模型。
    作为本发明一种涡轴发动机动态实时模型建模方法进一步的优化方案,步骤A4)中利用稳态关系曲线插值获得稳态基点的详细步骤如下:
    将燃气涡轮转子转速换算值与预先设定的转子转速换算阈值进行比较,如果燃气涡轮转子转速换算值大于等于预先设定的转子转速换算阈值,将稳态关系曲线以2%的间隔进行插值获得稳态基点;果燃气涡轮转子转速换算值小于预先设定的转子转速换算阈值,将稳态关系曲线以3%-5%的间隔进行插值获得稳态基点。
    作为本发明一种涡轴发动机动态实时模型建模方法进一步的优化方案,步骤B)的具体步骤如下:
    步骤B1),根据试车数据中动态过程的数据,利用剩余燃油量动态系数法的定义,分别计算出不同燃气涡轮转子转速换算值下的燃气涡轮转速动态系数、燃气涡轮出口温度动态系数、压气机出口压力动态系数和动力涡轮输出功率动态系数,以获得动态系数插值表;
    步骤B2),建立慢车以上状态燃气发生器动态实时模型,该模型的输入参数为燃油流量、进口总压、进口总温,输出参数为燃气涡轮转速、压气机出口压力、燃气涡轮出口温度,采用剩余燃油流量动态系数法计算出模型输出参数;
    步骤B3),建立慢车以上状态动力涡轮及负载动态实时模型,该模型输入参数为燃气涡轮转速、燃油流量、进口总压、进口总温、负载杆角度,模型输出参数为动力涡轮转速、动力涡轮输出功率,分别利用剩余燃油流量动态系数法和动力涡轮转子功率平衡求解模型输出参数;
    步骤B4),建立涡轴发动机起动阶段模型,该模型包括起动机单独带转阶段的模型、起动机与燃气涡轮共同带转阶段的模型、以及燃气涡轮单独带转阶段的模型;
    步骤B5),建立涡轴发动机动态实时模型,该模型包含步骤A1)建立的稳态实时模型、步骤B2)建立的慢车以上燃气发生器动态实时模型、,步骤B3)建立的动力涡轮及负载动态实 时模型、以及步骤B4)建立的起动阶段数学模型,根据各个模型的输入参数解算各个模型的输出参数。
    作为本发明一种涡轴发动机动态实时模型建模方法进一步的优化方案,所述动态系数包含加速系数和减速系数。
    作为本发明一种涡轴发动机动态实时模型建模方法进一步的优化方案,步骤B2)中采用剩余燃油流量动态系数法计算出模型输出参数时,燃气涡轮转速以前一时刻的动态值为当前时刻的初始值进行累加,压气机出口压力和燃气涡轮出口温度以当前时刻的稳态值作为初始值进行累加。
    作为本发明一种涡轴发动机动态实时模型建模方法进一步的优化方案,步骤B3)中采用剩余燃油流量动态系数法计算出模型输出参数时,动力涡轮输出功率以当前时刻的稳态值作为初始值进行累加。
    作为本发明一种涡轴发动机动态实时模型建模方法进一步的优化方案,步骤B2)中采用剩余燃油流量动态系数法计算出燃气涡轮出口温度后,作为参数输入至一阶惯性温度传感器模型,以获得能够反映延迟效应的出燃气涡轮出口温度输出参数。
    一种基于所述涡轴发动机动态实时模型建模方法的系统,包括试车数据预处理???、稳态工作区域判定逻辑???、动态工作区域判定逻辑???、稳态基点自动生成???、动态加减速系数自动生成??楹投凳蹦P图扑阍诵心?椋?
    所述试车数据预处理??橛糜诙栽际猿凳萁刑蕹暗愫推交瞬ù?;
    所述稳态工作区域判定逻辑??橛糜诟菥鍪猿凳菰ご砟?榇砗蟮氖猿凳葜形忍痰氖?,获得稳态工作点数据;
    所述动态工作区域判定逻辑??橛糜谂卸ň鍪猿凳菰ご砟?榇砗蟮氖猿凳葜械亩ぷ鞣段?,获得不同燃气涡轮转子转速换算值下燃气涡轮转速动态系数、燃气涡轮出口温度动态系数、压气机出口压力动态系数以及动力涡轮输出功率的加速系数与减速系数;
    所述稳态基点自动生成??橛糜诟菸忍ぷ髑蚺卸呒?橹猩傻奈忍ぷ鞯闶?、采用聚类合并法和拟合插值法计算出稳态基点,生成稳态基点插值表;
    所述动态加减速系数自动生成??橛糜诟菟龆ぷ髑蚺卸呒?橹猩傻牟煌计新肿幼倩凰阒迪氯计新肿俣凳?、燃气涡轮出口温度动态系数、压气机出口压力动态系数以及动力涡轮输出功率的加速、减速系数,采用拟合插值法获得其在稳态基点对应下的加速、减速系数,生成动态加减速系数插值表;
    所述动态实时模型计算运行??橛糜诟菟鑫忍阕远赡?樯傻奈忍悴逯当砗退龆蛹跛傧凳远赡?橹猩傻亩蛹跛傧凳逯当?,采用动态加减速 系数法、根据涡轮发动机动态实时模型的输入参数计算出发动机动态实时模型的输出参数。
    本发明采用以上技术方案与现有技术相比,具有以下技术效果:
    (1)本发明提出的涡轴发动机动态实时模型构建方法,能适用于在缺乏部件特性数据只利用试车数据情况下的涡轴发动机建模,是一种有效的数据驱动的涡轴发动机模型建模方法。
    (2)本发明提出的涡轴发动机动态实时模型构建方法,为改善动态模型精度,建模时引入了加速过程和减速过程两组动态系数,该模型能模拟包含起动、慢车以上状态的涡轴发动机全状态工况,且模型计算耗时小,能满足涡扇发动机机载实时需求;
    (3)本发明设计的涡轴发动机动态实时模型系统,能自动生成稳态基点插值表和动态系数插值表,将其存储在该软件相应路径下即可获得基于该试验数据的发动机动态实时模型,降低劳动强度,节约设计开发成本,且并不局限于某一型号涡轴发动机,具有一定的通用性。
    附图说明
    图1是涡轴发动机动态实时模型建模方法的流程图;
    图2是某型涡轴发动机慢车以上状态稳态区域聚类合并效果对比图;
    图3是某型涡轴发动机动态实时模型(含起动)稳态基点插值表;
    图4是某型涡轴发动机动态实时模型(含起动)动态加减速系数插值表;
    图5是某型涡轴发动机动态实时模型起动阶段燃气涡轮转速仿真效果图;
    图6是某型涡轴发动机动态实时模型起动阶段压气机出口压力仿真效果图;
    图7是某型涡轴发动机动态实时模型起动阶段燃气涡轮出口温度仿真效果图;
    图8是某型涡轴发动机动态实时模型的燃气发生器的燃气涡轮转速仿真效果图;
    图9是某型涡轴发动机动态实时模型的燃气发生器的压气机出口压力仿真效果图;
    图10是某型涡轴发动机动态实时模型的燃气发生器的燃气涡轮出口温度仿真效果图;
    图11是某型涡轴发动机动态实时模型的动力涡轮及负载仿真效果图。
    具体实施方式
    下面结合附图对本发明的技术方案做进一步的详细说明:
    本发明的思路是针对先进涡轴发动机在缺少部件特性数据情况下的发动机实时模型需求,在现有航空发动机通用模型建模方法基础上进行扩展,利用剩余燃油流量动态加减速系数法,建立起包括起动模型、慢车以上状态燃气发生器模型和动力涡轮及负载模型在内的涡轴发动机动态实时模型,相比于非线性部件级模型该模型计算耗时明显减小,能满足机载实时性要求,适用于缺少部件特性只依据试车数据的建模,同时设计的发动机动态实时模型自动生成软件能自动获取稳态基点插值表和动态加减速系数表,以固定格式与路径保存于编制 的动态模型系统中即可获得涡轴发动机动态模型,该方法是涡轴发动机实时模型的数据驱动建模的有益尝试,能显著降低劳动强度,减少设计开发成本。
    本发明的具体实施方式以某型涡轴发动机动态实时模型构建与模型自动生成为例,图1是涡轴发动机动态实时模型建模计算流程图,该动态实时模型及系统的建立方法包括以下步骤:
    步骤A),建立涡轴发动机稳态实时模型;
    步骤B),采用剩余燃油流量动态加减速系数法建立涡轴发动机动态实时模型。
    其中步骤A)的详细步骤如下:
    步骤A1)、涡轴发动机稳态实时模型参数相似归一化换算,由于试车数据是在不同的大气条件下进行试车试验获得的,因此需对原始试车数据进行相似变换,将其换算到标准大气状况下(P1=101325Pa,T1=288.15K),以便将其用于适用于不同大气条件的发动机简化模型建模中。在环境温度为T1,环境压力为P1的大气条件下,燃气涡轮转速Ng、动力涡轮输出功率Ne、燃油流量Wf、燃气涡轮出口温度T45和压气机出口压力P3的换算公式如下
    Ngc=Ng288.15T1+273.15---(1)]]>
    Nec=Ne101325P1288.15T1+273.15]]>
    Wfc=Wf101325P1288.15T1+273.15]]>
    T45c=T45288.15T1+273.15-273.15]]>
    P3c=P3101325P1]]>
    将试车数据按以下次序:P1、T1、Ng、Mp、Np、Wf、LDL存入指定路径下的data.txt文本文件中,参考步骤公式(1)计算出各物理参数的换算参数。
    步骤A2)、建立涡轴发动机稳态模型前应对稳态数据进行预处理。试车数据具有较好的一致连续性,但是试车数据中存在个别野点,根据设计值的±6%为其容许变化范围,当某一参数时间序列中某一采样点参数大于其前后连续3个点区域范围的容许变化范围,则判定野点,舍弃该值,取其前一采样点参数值为该时刻点参数值。
    动态试车数据预处理,发动机原始试车数据具有较严重的毛刺和噪声,为了建模需要,首先应对原始动态数据进行平滑处理。利用移动平均平滑方法对原始数据进行平滑,该方法可表述为:
    y^0(k)=y(k)+y(k-1)+...+y(k-N-1)N---(2)]]>
    式中N为取均值的个数。N值越大,得到的曲线越光滑;但是在数据真实值变化较大的地方,N的取值太大则会造成平滑以后的信号失真。通过试取发现N=20对原始数据进行平滑效果较好。
    对发动机原始试车数据进行平滑滤波处理,采用VC++6.0和matlab混合编程方法,即首先在matlab中编写的滤波程序的M文件,然后将其封装成动态链接库dll,再在VC++6.0中根据一定的规范调用此dll来对原始试车数据进行滤波,平滑滤波的窗口参数在VC++6.0模型软件中设置。
    步骤A3),根据试车数据自动确立稳态区域,对工作状态相近的稳态数据进行聚类合并,以此构造稳态工作点。
    预处理后的试车数据中存在多段相对稳定的区域,需从这些相对稳定的区域中获取稳态值,经过数据分析后发现如果燃气涡轮转速Ng达到相对稳定后,其它参数如燃油流量Wf、T45、P3也达到相对稳定,故只需考虑Ng达到相对稳定后,发动机即达到相对稳定状态,其他物理参数基本保持状态一致。
    依据试车数据中稳定状态保持时间可设置稳态区域判定窗口长度,在此设定该窗口长度为80,设置稳定状态区域的变化阈值为±0.03%,待判定的数据间隔为3个采样点,则稳定状态判定表达式如下
    |Ngc(i+3)-Ngc(i)|<0.03%  i=1,…,80   (3)
    若试车数据中连续80步采样数据都满足上式关系,则可判定系统处于稳态区域,式中的窗口长度,变化阈值和间隔都可在软件中设置,将稳定区域的索引值输出,如图2(a)所示,图中第一列为起始地址,第二列为终止地址,即第一行代表一个稳定区域,从试车数据中第267各点到765个点,以此类推,在这些稳定区域内对Wf、T45、P3、Mkp、Np分别求和取平均值作为稳态值,生成文本文件value.txt,如图2(b)所示。
    在初步确定稳态区域后,需对相近稳态区域值合并,通过均值处理获得相近稳态区域的一个稳态值,图2(c)给出了稳态区域相近的索引值,图中数字0-11分别代表图2(b)中各稳态区域,从上到下以此为0、1、2…11,一共12个稳态区域。通过聚类算法分析有以下稳定区域相近,分别是第1段和第12段稳定区域值,第4段、第6段和第9段稳定区域值,第5段、第7段、第8段、第10段和第11段稳定区域值,对这些相近的稳态区域分别对各参数求和取平均值作为新稳态值,经过升序排序后,可得如图2(d)所示的合并后的稳态插值表,共包含5段稳态区域,其Ngc值分别为59.36、80.55、83.02、88.93和97.10。
    步骤A4,根据步骤A3)中获得的稳态工作点,采用多项式拟合法获得稳态关系曲线,利用稳态关系曲线均匀插值获得稳态基点,并以此构造稳态基点插值表。
    稳态插值基点太少会引起插值的误差加大,根据已得到的5个稳态工作点利用多项式拟合法,将慢车以上状态稳态插值表中的插值基点扩展为19个,其中在较高工作状态下稳态基点选取间隔需比高工作状态的稳态基点间隔小,所以这里需要人为地添加。
    由于起动阶段没有稳定工作点,故需人为定义Ngc的稳态插值基点值,针对该型涡轴发动机选取稳态基点值对应的Ngc为:0、11.7、14、17、22、27、32、37、42、47、52。根据定义的Ng的稳态插值基点值,直接从预处理后的试车数据中对应的Wf、T45、P3数值作为稳态插值基点。图3给出了利用动态实时模型通用自动生成软件获得的涡轴发动机(含起动)稳态基点插值表。
    由于涡轴发动机主要工作在空中慢车以上状态(以燃气涡轮转子转速换算值为发动机工作状态的表征),为保证模型精度,此例中取换算的空中慢车的燃气涡轮转速为81%,即Ngc大于81%时以较小的间隔选点作为稳态基点插值数据,而在涡轴发动机较低工作状态下,以较大的间隔选点作为稳态插值数据,以减小插值时的计算量。
    步骤A5、稳态建模时根据稳态基点插值数据表,采用一维插值法来建立发动机的稳态实时模型。定义两个数组,分别用于存放自变量x和因变量y的插值样本。则插值算法如下
    y=x-x1x2-x1y2+x2-xx2-x1y1(x<x1)x-xixi+1-xiyi+1+xi+1-xxi+1-xiyi(xi<x<xi+1,1in)x-xn-1xn-xn-1yn+xn-xxn-xn-1yn-1(x>xn)---(4)]]>
    式中y为一维插值输出变量,x为插值变量,n为插值样本点的个数。涡轴发动机稳态实时模型实现如下:首先对模型的输入参数Wf进行相似换算,然后根据公式(4)对应关系进行稳态插值计算,最后通过相似换算的逆运算计算出各物理量的实际值。
    Ngc=f1(Wfc)   (5)
    T45c=f2(Ngc)
    P3c=f3(Ngc)
    Nec=f4(Ngc)
    式中f1,f2,f3,f4为不同物理量换算参数之间的稳态关系,如步骤A4所述方法构建的稳态基点插值数据表。
    步骤B)的详细步骤如下:
    步骤B1)、基于剩余燃油流量的动态系数法,发动机稳态工作过程中,每个Ngc都对应着一个需求燃油量。当实际燃油流量Wfc与需求燃油量不匹配时,就不能维持发动机稳态平衡工作,发动机进入动态工作过程。
    假设某稳定状态下燃气涡轮转速Ngc对应的需求油量Wfcs=f-1(Ngc),定义此时实际燃油流量Wfc与需求燃油量Wfcs的差为剩余供油量
    ΔWfc=Wfc-Wfcs   (6)
    燃气发生器动态过程中各状态参数的变化与ΔWfc存在对应关系,这种对应关系可以通过动态系数表示:当ΔWfc>0时,即实际供油量大于需求油量时,燃气涡轮转子存在正的剩余功率,Ngc增大;反之当ΔWfc<0时,Ngc减小。并且ΔWfc越大,Ngc变化越快。定义第i时刻Ngc、T45c、Nec的动态系数分别如下
    Kng(i)=Ngc(i+1)-Ngc(i)ΔWfc(i)&CenterDot;T]]>
    KT45(i)=T45c(i)-T45cs(i)ΔWfc(i)]]>(7)]]>
    KP3(i)=P3c(i)-P3cs(i)ΔWfc(i)]]>
    KNe(i)=Nec(i)-Necs(i)ΔWfc(i)]]>
    式中T为仿真步长,取试车数据的采样步长为0.02s,T45cs(i)、P3cs(i)、Necs(i)分别为Ngc(i)对应的第i步稳态燃气涡轮出口温度、燃气涡轮出口压力、动力涡轮输出功率,可利用公式(5)给出的稳态插值关系计算出。
    在动态系数Kng、KT45、KP3和Kne关于Ngc的动态系数插值表前需对试车数据进行分析,判断其是否反映加速或减速过程。针对慢车以上过程,确定稳态区域后即可得到其区间索引值如图2(a)所示,如果区间的末端索引值i的Ngc值小于末端索引值后十个点的Ngc值,即 式(8-a)成立,则判断此稳态区域到下一稳态区域间为加速过程;如果区间的末端索引值i的Ngc值大于末端索引值后十个点的Ngc值,即式(8-b)成立,则判断此稳态区域到下一稳态区域间为减速过程。
    Ngc(i)<Ngc(i+j) j=1,…,10  (a)   (8)
    Ngc(i)>Ngc(i+j) j=1,…,10  (b)
    采用公式(7)的动态系数计算公式可求得此稳态点下的加速系数和减速系数,此时获得了m个稳态点的加速系数和剩余n个稳态点的减速系数??悸堑蕉凳逯当碇忻扛鑫忍愣夹瓒杂σ桓黾铀傧凳鸵桓黾跛傧凳?,这里采用根据减速段的n个稳态点的减速系数进行一维插值得到m个稳态点的减速系数,n个稳态点的加速系数同样可以根据m个稳态点的加速系数插值得到。与扩展稳态基点方法相同,分别根据加速系数和减速系数采用多项式拟合法,构造出每个稳态基点下加速系数和减速系数,共计60个。
    在起动阶段,将T45c、P3c、Nec的动态加速系数和减速系数都默认为0,只考虑Ngc的动态加、减速系数。同样由于没有稳态区域,只需在步骤A3中定义的稳态插值基点值,从试车数据中根据其后一点的Ngc和Wfc数值,利用公式(7)计算得到Ngc的动态加、减速系数。图5给出了利用动态实时模型通用自动生成软件获得的涡轴发动机(含起动)动态加速、减速系数插值表。
    步骤B2)、建立慢车以上状态的燃气发生器动态实时模型,确定该模型的输入参数为Wf、P1、T1,输出参数为Ngc、燃气涡轮出口压力P3c、燃气涡轮出口温度T45c。
    (1)燃气涡轮转速动态计算
    给定燃气发生器初始换算转速Ngc(0)和供油规律,依据式(6)求出ΔWfc,则由式(7)可得燃气发生器换算转速变化率及转速
    N&CenterDot;gc(i)=Kng(i)ΔWfc(i)---(9)]]>
    Ngc(t)=Ngc(0)+&Integral;0tN&CenterDot;gc(i)dt]]>
    式中Kng按照动态系数表由Ngc插值获得。
    (2)燃气涡轮出口温度动态计算
    换算燃气涡轮出口温度T45c按下式计算
    T45c(t)=T45cs(t)+KT45(t)·ΔWfc(t)   (10) 式中T45cs与KT45分别根据动态系数表由Ngc插值求得。由于热电偶的滞后影响,使T45测量值与发动机实际值差异明显,故将温度传感器近似等效为一个惯性环节,则燃气涡轮出口温度传感器量测输出的T45cr表达式为
    T45cr(t)=T45cr(t-1)+[T45c(t)-T45cr(t-1)]·δ
    δ=1TΔT+1---(11)]]>
    式中T45cr为传感器真实量测值,δ为温升惯性常数,T为时间常数(根据该型传感器取值为3),采样步长ΔT为0.15s。
    (3)压气机出口压力动态计算
    燃气涡轮出口压力P3c
    P3c(t)=P3cs(t)+Kp3(t)·ΔWfc(t)   (12)
    式中KP3根据动态系数表由Ngc插值求得。
    步骤B3)、建立慢车以上状态的动力涡轮及负载动态实时模型,确定该模型的输入参数为Ng、Wf、P1、T1、负载杆角度LDL,模型输出参数为Np、Ne。
    (1)动力涡轮输出功率动态计算
    与T45c计算类似,采用动态系数法求解Nec,如下式所示
    Nec(t)=Necs(t)+KNe(t)·ΔWfc(t)   (13)
    (2)动力涡轮转速动态计算
    利用负载角LDL和Npc=100%时的动力涡轮转速得到Npc=100%时的换算需求功率,此时需求功率Nv0与输出功率平衡,据此对Nv0时的负载特性进行缩放修正,得到LDL和Nv0的插值表。当Npc=100%时,可直接利用LDL从插值表中插值出需求功率Nv0,而当Npc≠100%时,依据以下公式计算Nvc
    Nvc=Nv0&CenterDot;(Npc100)3---(14)]]>
    动力涡轮转子的转子动力学方程式
    Nec(i)=Jpt·Npc(i)·(Npc(i)-Npc(i-1))/T+Nvc(i)   (15)
    式中Jpt为动力涡轮转子的转动惯量,根据式(15)求解动力涡轮转子加速度,再由下式计算动力涡轮转速
    Npc=Npc(0)+&Integral;0tN&CenterDot;pcdt---(16)]]>
    步骤B4)、涡轴发动机的起动过程一般包括起动机单独带转、起动机和燃气涡轮共同带转、起动机脱开燃气涡轮带转这三个阶段。第一阶段下发动机由起动机单独带转,起动机扭矩Mg使发动机加速,直至达到点火转速n1;第二阶段下起动机继续带转,同时燃烧室点火涡轮开始产生扭矩,两者共同使燃气发生器转速Ngc加速。随着Ngc的增加,涡轮功率迅速增加,当涡轮功率大于压气机功率,这时涡轮能够独立带动压气机加速。然而,为了起动可靠,通常都转速上升到一定转速n2,涡轮功率已经大大超过压气机功率,起动机才停止工作;第三阶段下起动机脱开,涡轮带动压气机运转,依靠剩余功率继续加速直至达到地面慢车状态。
    (1)起动机单独带转阶段(n1=11.7%,即Ngc<11.7%)
    利用起动阶段的试车数据获得起动机特性,通过多项式拟合法燃气涡轮转子剩余扭矩ΔMg与Ngc的函数关系式:
    ΔMg=-0.01559·Ngc3+0.58955·Ngc2-6.97774·Ngc+27.28422   (17)
    燃气发生器换算转速Ngc:
    N&CenterDot;gc=ΔMg/Jgt]]>
    Ngc=&Integral;0tN&CenterDot;gcdt---(18)]]>
    这一阶段P3c、T45由Ngc进行一维插值获得。燃气涡轮不做功,Nec和Npc均为零。
    (2)起动机涡轮共同带转阶段(n2=42%,即11.7%<Ngc<42%)
    燃气发生器转子加速度由起动机扭矩和剩余供油量共同决定:
    N&CenterDot;gc=ΔMg/Jgt+KngΔWfc=1+KngΔWfc---(19)]]>
    式中第一项为起动机扭矩产生的燃气涡轮转子加速度分量,依据经验将其设为1。起动过程根据起动的T45c,P3c直接由Ngc插值得到,Nec和Npc均为零。
    起动机脱开阶段由燃气涡轮单独带转,故建模方法参照步骤B2慢车以上加减速动态实时模型建模方法获得。
    步骤B5)中动态实时模型系统中包括已建立的稳态实时模型、慢车以上燃气发生器动态实时模型、动力涡轮及负载动态实时模型,以及起动阶段数学模型。
    在Visual C++开发平台上以上述功能??槔炊ㄒ逵敕庾岸杂?,模型系统中稳态基点插值表和动态加减速系数插值表自动生成的步骤主要有:试车数据预处理??橛糜诙栽际猿凳萁刑蕹暗愫推交瞬ù?,对预处理后的试车数据进行相似换算;在稳态工作区域判定逻辑??橹薪卸ㄔご砗蟮氖猿凳莸奈忍ぷ鞣段?,获得稳态范围索引值和稳态工作点数据;在稳态基点自动生成??橹?,利用聚类合并法计算稳态基点,生成稳态基点插值表;在动态工作区域判定逻辑??橹信卸ㄔご砗蟮氖猿凳葜械亩ぷ鞣段?,获得动态范围索引值,及其加速系数与减速系数;动态加减速系数自动生成??橛糜诟菀延械募铀?、减速系数,采用拟合插值法获得稳态基点对应下的加速、减速系数,生成动态加减速系数插值表;而动态实时模型计算运行??橛糜诟萆傻奈忍悴逯当砗投蛹跛傧凳逯当?,采用动态加减速系数法解释发动机动态实时模型。
    将生成的涡轴发动机稳态基点插值表、动态系数插值表按指定路径存入涡轴发动机动态实时模型系统中,在模型系统中通过模型逻辑功能??橛肽P褪淙肓拷馑隳P褪涑霾问?。当输入飞行条件参数和控制参数后,该模型系统将控制参数进行相似归一化变换,换算至标况后采用剩余燃油流量动态加减系数法,解算模型输出参数,将其进行相似归一化反变换获得涡轴发动机物理量的实际值。
    图5-图11分别给出了涡轴发动机动态实时模型中起动阶段、燃气发生器和动力涡轮及负载各物理参数的输出参数,其中飞行条件为T1=31.84℃,P1=100.94Kpa。由上述仿真结果可以看出,根据试车数据自动生成的涡轴发动机动态实时模型具有较高的稳态、动态精度,其中稳态误差不大于2%,动态误差一般在5%以内。
    需要指出的是,以上所述,仅为本发明的具体实施方式,但本发明的?;し段Р⒉痪窒抻诖?,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化和替换,都应涵盖在本发明的?;し段?。因此,本发明的?;し段вσ运鋈ɡ蟮谋;し段?。

    关于本文
    本文标题:一种涡轴发动机动态实时模型建模方法及系统.pdf
    链接地址://www.4mum.com.cn/p-6147548.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 任九投注技巧 11选5计划软件免费版 分分彩组选包胆 彩名堂免费计划软件2.0.6 北京pk10直播开奖 快三大小单双计划预测软件 凤凰③肖默认论坛 稳赚六肖是什么网址 大乐透玩法中奖规则表 安徽时时哪里买 体彩大乐透手机扫码 手机体彩app 三分pk10稳赚技巧 北京pk10计划看不懂 足彩只买平局能赚6 幸运飞艇双面玩法怎么玩