CN112257174A - 一种飞机供电系统智能化集成设计系统及方法 - Google Patents
一种飞机供电系统智能化集成设计系统及方法 Download PDFInfo
- Publication number
- CN112257174A CN112257174A CN202011092419.6A CN202011092419A CN112257174A CN 112257174 A CN112257174 A CN 112257174A CN 202011092419 A CN202011092419 A CN 202011092419A CN 112257174 A CN112257174 A CN 112257174A
- Authority
- CN
- China
- Prior art keywords
- power supply
- network topology
- model
- multidisciplinary
- supply system
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000013461 design Methods 0.000 title claims abstract description 110
- 238000000034 method Methods 0.000 title claims abstract description 66
- 238000005457 optimization Methods 0.000 claims abstract description 54
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 36
- 238000010801 machine learning Methods 0.000 claims abstract description 9
- 239000000446 fuel Substances 0.000 claims description 20
- 230000008878 coupling Effects 0.000 claims description 19
- 238000010168 coupling process Methods 0.000 claims description 19
- 238000005859 coupling reaction Methods 0.000 claims description 19
- 230000008569 process Effects 0.000 claims description 19
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 7
- 238000012549 training Methods 0.000 claims description 7
- 238000007418 data mining Methods 0.000 claims description 6
- 239000003795 chemical substances by application Substances 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 4
- 230000009467 reduction Effects 0.000 claims description 3
- 238000010845 search algorithm Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 2
- 239000004576 sand Substances 0.000 claims description 2
- 238000003860 storage Methods 0.000 claims description 2
- 230000010354 integration Effects 0.000 abstract description 8
- 230000008859 change Effects 0.000 abstract description 6
- 230000005291 magnetic effect Effects 0.000 description 43
- 238000004804 winding Methods 0.000 description 28
- 238000004146 energy storage Methods 0.000 description 18
- 230000005284 excitation Effects 0.000 description 10
- 230000006870 function Effects 0.000 description 9
- 238000012216 screening Methods 0.000 description 9
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 8
- 229910000831 Steel Inorganic materials 0.000 description 8
- 239000010959 steel Substances 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 7
- 229910000976 Electrical steel Inorganic materials 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000002790 cross-validation Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000004907 flux Effects 0.000 description 5
- 230000006698 induction Effects 0.000 description 5
- 238000005452 bending Methods 0.000 description 4
- 239000004020 conductor Substances 0.000 description 4
- 230000005415 magnetization Effects 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 3
- 239000010949 copper Substances 0.000 description 3
- 238000007599 discharging Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 229910052742 iron Inorganic materials 0.000 description 3
- 238000007637 random forest analysis Methods 0.000 description 3
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 2
- 239000003990 capacitor Substances 0.000 description 2
- 229910052802 copper Inorganic materials 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 238000003066 decision tree Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005538 encapsulation Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 239000000295 fuel oil Substances 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- 229910052744 lithium Inorganic materials 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000010349 pulsation Effects 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 239000011149 active material Substances 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 230000002500 effect on skin Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000002330 electrospray ionisation mass spectrometry Methods 0.000 description 1
- 238000013401 experimental design Methods 0.000 description 1
- 230000005294 ferromagnetic effect Effects 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- -1 machinery Substances 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000011241 protective layer Substances 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
Abstract
本发明公开一种飞机供电系统智能化集成设计系统及方法,在多学科建模模块中建立不同部件模型;在代理建模模块中设定设计参数变化区间,可以获得多学科模型的数据样本,可以分析设计参数在其变化区间内对设计目标的影响程度,提取多学科模型的关键参数,使用机器学习算法建立多学科模型的近似模型,满足系统集成设计精度和速度的需要;在网络拓扑优化设计模块中输入系统已知的电源、汇流条和负载信息,添加和删除供电系统网络拓扑应满足的约束条件,并自动生成网络拓扑可行解;在容量优化配置模块中输入系统功率需求,根据对工作状态进行能量优化规划的结果,优化配置系统各部件容量,实现高效率、低成本的飞机供电系统智能化集成设计。
Description
技术领域
本发明属于飞机供电系统智能化集成设计领域,具体涉及一种飞机供电系统智能化集成 设计方法。
背景技术
飞机供电系统智能化集成设计是一个复杂而庞大的工程设计问题,具有多层次、高耦合 强度、高维设计空间,高维多目标、多约束等特点。在目前的飞机供电系统设计过程中,在 获得飞机分解的性能指标后,设计人员使用不同的方法和工具对各个系统开展独立设计,系 统设计所需的假设条件和设计点通常来自于前期产品的经验;在集成过程中,来自不同供应 商的设备和系统之间通过人工和文档的形式实现信息的交互。这样的设计方式具有以下缺陷: ①设计人员根据经验,通过解耦和分解获得不同设备和系统的设计要求,往往具有片面性, 难以保证设计的覆盖率,也不适合用于新技术的应用;②根据逆向工程依次从负载、配电到 发电系统进行串行设计,一方面可能会带来大量人工迭代,提高了设计的难度和复杂度;③ 另一方面,无法最大化利用系统性能,还可能导致系统设计的保守性;④基于文档的设计受 到文件语言描述的歧义性、文件传递的静态非结构性等限制;⑤此外,随着电气化技术的发 展,飞机中原本独立分散、互不影响的气压、液压系统替换为电能系统,大大提高了飞机能 源系统的集成化程度,这种参数之间的复杂耦合关系也给系统的集成设计带来了困难。因此, 现行的系统设计方法已无法满足未来的发展需求。
发明内容
本发明针对现有技术问题,提出一种飞机供电系统智能化集成设计方法,可以高效、低 成本地完成飞机供电系统集成设计,获得综合性能最佳的飞机供电系统方案,具有重要的理 论和工程应用价值。
本发明具体采用如下技术方案:
本发明首先提供一种飞机供电系统智能化集成设计方法,其特征在于包括以下步骤:
步骤一:针对飞机供电系统不同部件内部参数的多物理耦合关系建立飞机供电系统各部 件的多学科模型,并设定各部件多学科模型的参数;
步骤二:对各部件多学科模型进行试验设计,提取表征部件多学科模型设计空间特征的 输入、输出参数样本,使用数据挖掘方法提取部件多学科模型的关键特征参数,利用机器学 习算法建立满足系统集成设计需求的各部件近似模型;
步骤三:根据已知的设计条件建立飞机供电系统网络拓扑模型,并根据各部件的近似模 型建立系统网络拓扑设计目标和约束条件,使用多目标优化算法求解获得满足设计要求的最 优网络拓扑方案;
步骤四:根据步骤三获得的系统最优网络拓扑方案生成系统可能存在的工作状态集合, 在此基础上建立每个工作状态下系统能量调度模型,借助优化算法对系统进行能量优化调度 的计算求解,而后对不同工作状态下系统能量调度方案进行综合,获得对应网络拓扑满足不 同工作状态工作需求的容量配置方案,并建立以部件容量、节点电压为输入的容量配置模型, 借助多目标优化算法对网络拓扑、部件容量、节点电压进行优化,获得性能最佳的系统容量 配置方案。
本发明的另一目的是提供一种飞机供电系统智能化集成设计系统,其特征在于包括:
多学科建模模块,针对飞机供电系统不同部件内部参数的多物理耦合关系建立飞机供电 系统各部件的多学科模型,并设定各部件多学科模型的参数;
代理建模模块,对各部件的多学科模型进行试验设计,提取表征部件多学科模型设计空 间特征的输入、输出参数样本,使用数据挖掘方法提取部件多学科模型的关键特征参数,利 用机器学习算法建立满足系统集成设计需求的各部件近似模型;
网络拓扑优化设计模块,根据已知的设计条件建立飞机供电系统网络拓扑模型,并根据 各部件的近似模型建立系统网络拓扑设计目标和约束条件,使用多目标优化算法求解获得满 足设计要求的最优网络拓扑方案;
容量优化配置模块,根据系统最优网络拓扑方案生成系统可能存在的工作状态集合,在 此基础上建立每个工作状态下系统能量调度模型,借助多目标优化算法对系统进行能量优化 调度的计算求解,对不同工作状态下系统能量调度方案进行综合,获得对应网络拓扑满足大 部分工作状态工作需求的容量配置方案并建立以部件容量、节点电压为输入的容量配置模型, 借助多目标优化算法对网络拓扑、部件容量、节点电压进行优化,获得性能最佳的系统容量 配置方案。
其中,飞机供电系统部件包括发电机、燃料电池、锂电池、超级电容、电感、电容、固态开关、变压器和线缆等。
多物理耦合具体包括电、气、液、机械、材料、结构、热等物理领域之间参数的耦合,多学科模型具体是指以多物理耦合参数为输入,以重量、体积、损耗、成本、可靠性、维修 性等一系列性能参数为输出的模型。
本发明技术效果如下:
1、本发明功能全面、操作简单。在多学科建模模块中建立不同部件模型;在代理建模模 块中设定设计参数变化区间,可以获得多学科模型的数据样本,可以分析设计参数在其变化 区间内对设计目标的影响程度,提取多学科模型的关键参数,使用机器学习算法建立多学科 模型的近似模型,满足系统集成设计精度和速度的需要;为了满足不同飞机的设计需求,在 网络拓扑优化设计模块中输入系统已知的电源、汇流条和负载信息,添加和删除供电系统网 络拓扑应满足的约束条件,并自动生成网络拓扑可行解;在容量优化配置模块中输入系统功 率需求,自动完成工作状态的生成、缩减和综合,自动对每个工作状态进行能量优化规划, 并根据对工作状态进行能量优化规划的结果,优化配置系统各部件容量,自动生成供电系统 可靠性最高、重量最轻、损耗最小的方案。
2、本发明高效率、低成本地完成从部件层面到系统层面供电系统智能化设计,在实现程 序中可以清晰展现部件和系统模型、优化设计的输入参数、优化过程的进度、优化设计的结 果。
附图说明
图1为一种飞机供电系统智能化集成设计方法;
图2为部件多学科模型框架示意图;
图3为代理建模模块框架示意图;
图4为网络拓扑优化设计模块框架示意图;
图5为供电系统容量优化配置模块示意图;
图6为飞机供电系统智能化集成设计的信号流程和关系图;
图7为同步发电机简化剖面;
图8为最优拉丁超立方试验设计方法;
图9为电机多学科模型试验设计结果;
图10为多学科模型特征参数提取方法;
图11为封装法原理;
图12为k折交叉验证原理;
图13为电机重量的近似模型精度;
图14为电机损耗的近似模型精度;
图15为网络拓扑优化设计帕累托图;
图16为网络拓扑优化设计方案之一;
图17为容量优化配置案例。
具体实施方式
一种飞机供电系统智能化集成设计系统,整体框架如图1-5所示,包括多学科建模模块、 代理建模模块、网络拓扑优化设计模块、容量优化配置模块,模块具体的功能和关系如图6 所示。
下面以某型飞机供电系统设计为例对本发明设计方法进行详细描述,具体步骤如下:
A:针对飞机供电系统不同部件内部参数的多物理耦合关系建立供电系统发电机、燃料电 池、蓄电池、变换器、开关、电缆等各多学科模型。多学科建模模块利用各部件多物理耦合 原理建立各部件的多学科模型,多物理耦合指多种物理领域之间参数的耦合,多学科模型具 体是指以多物理耦合参数为输入,以性能参数为输出的模型。以同步发电机为例:
在工程经验对电磁负荷、尺寸比等参数的取值范围内,对电机各部分尺寸进行计算。
计算功率P由额定功率PN根据下式计算得到
式中ηN为额定负载时的效率;为额定负载时的功率因数;Kg为考虑率发电机电枢 压降和绕组电流而引入的系数;Km为电动机电枢压降和绕组电流而引入的系数;KE为满载电 势标幺值,即满载时感应电势和端电压的比值,可由下式确定
KE=1-εL=0.931+0.0108lnPN-0.013p 2)
轴的工作能力一般取决于强度和刚度,转速高时还取决于震动稳定性。假设轴为圆柱体, 则其关键设计参数为其轴向和径向长度。
电机轴传递扭矩的强度条件为
即
式中[τ]为许用扭转剪切应力,MPa;T为轴传递的扭矩,Nmm;W为轴的抗扭截面系数, 0.2d3(1-α4)mm3;P为轴传递的功率,kW;D为轴的外直径,mm;d为轴的内直径,mm;α 为轴的内外直径之比。
根据轴的刚度计算,
式中系数B值查表获得。
在轴的结构初步设计完成后,通常要对其进行弯扭合成强度校核。
根据第三强度理论,钢的强度条件为
可以得到
式中σe为当量应力,N/mm2;Me为当量扭矩,Nmm;M为危险截面上的合成弯矩,Nmm; W为轴危险截面弯曲截面系数,圆截面W≈0.1D3(1-α4);a为折合系数,不变扭矩的a≈0.3, 循环脉动扭矩的a≈0.59,频繁正反转的a=1,若扭矩变化规律不清则按循环脉动处理。
电机主要尺寸是指其电枢直径和有效长度,可以根据电磁负荷和尺寸比初步计算得到。
线负荷A是指电枢圆周单位长度电流导体数,即
式中Nt为电枢绕组并联数,I为电枢电流,D为电枢直径。
电机常数CA指,单位转矩所耗用的有效材料的体积,反映材料的消耗量,即
式中n为电机转速,α′为极弧系数,kW为绕组因数,Bδ为气隙磁密。
电机尺寸比λ定义为电机定子或电枢的有效长度与极距之比。
定子外径D5可根据经验估算为
D5=D+Kτ 11)
式中K为比值经验系数,2极电机取0.35~0.4,4极电机取0.52~0.57,6极电机取0.6~0.72。
以图7为例说明磁路计算过程。
气隙磁感应强度为
齿距离齿根1/3齿高处的齿宽b1/3为
式中bc为槽宽。
齿的磁密(距离齿根1/3齿高处的磁感应强度)为
查硅钢片磁化特性曲线可获得对应的磁场强度Hrt。
转子齿部磁势为
Frt=Hrtlrt 16)
转子轭磁感应强度为
查硅钢片磁化特性曲线可获得对应的磁场强度Hry。
转子轭磁势为
Fry=Hrylryξ 18)
式中ξ为轭部磁路长度因数,可查表确定。
电极磁通为
Φep=kσΦ0 19)
式中kσ为主极漏磁系数,直流电机kσ为1.15~1.2,交流电机kσ=1+10δ/τ<1.35。
电极磁感应强度为
式中bep为电极极身宽度。
查硅钢片磁化特性曲线可获得对应的磁场强度Hep。
主极磁势为
Fep=Heplep 21)
定子轭磁感应强度为
查硅钢片磁化特性曲线可获得对应的磁场强度Hsy。
定子轭磁势为
Fsy=Hsylsyξ 23)
最终得到
F=Fδ+Frt+Fry+Fep+Fsy 24)
气隙系数kδ是由槽口影响使气隙磁阻增加而引入的系数,其计算方法为
式中bZ为齿宽。
若定、转子均开槽,则
kδ=kδ1·kδ2 26)
极弧系数决定于励磁磁势分布曲线的形状、气隙的均匀程度和磁路饱和程度。若励磁磁 势正弦分布、气隙均匀、磁路不饱和,则α′=2/π=0.637;磁路越饱和,α′越大。
绕组的设计参数包括单匝导体长度以及绕组串并联数和导线串并绕数。其中,单匝导体 长度在电机尺寸确定的情况下基本确定,因此需要求解单相绕组串并联数N1、a和单匝绕组 串并绕数N2、a2。
电枢绕组设计和校核步骤如下:①根据输出电压,计算获得单相绕组串数和单匝绕组串 并绕数的乘积N1N2;②根据输出功率和电压,计算得到电枢电流;③根据GB/T6108.1-2003, 选择线规AWG 4以下的电磁线组合,而后考虑电压等级的绝缘需求,得到电磁线的备选组合; ④设定电枢绕组电流密度Ja,计算获得单相绕组并联数和单匝绕组并联数的乘积aa2;⑤在槽 满率小于80%和N2、a、a2为正整数的条件下,枚举校核N1、N2、a、a2的可能性;⑥若不存 在N1、N2、a、a2的组合则说明该电机设计参数无法设计出合理的电机;⑦若存在N1、N2、a、 a2的组合,则进一步计算绕组参数。
绕组等效直流电阻为
式中ρc为某线规线缆单位长度的等效电阻,其随温度的变化关系可表示为
ρc=ρ0(1+αT) 28)
由于集肤效应,绕组通以交流电时电阻值将增大,用系数KF表示电阻增加系数,则
Re=KFR 29)
每相绕组的主电抗可表示为
式中λm为主磁路比磁导。
漏抗则为
式中槽比漏磁导λs为
式中h为绕组的径向长度,h0为槽的剩余径向长度,bs0为槽口宽。
谐波比漏磁导λδ为
齿顶比漏磁导λt为
式中λtd为槽口面对机靴时的齿顶比漏磁导,λtd为槽口面对极间区域时的齿顶比漏磁导, bt为齿顶宽度。
电机端部比漏磁导λE为
励磁绕组的设计包括确定励磁绕组匝数、线规及励磁系统对额定励磁电流和电压的要求。
电枢漏抗和电阻的标幺值为
额定负载时,电动势的标幺值为
每极电枢反应磁动势为
直轴、交轴电枢反应磁动势为
式中Kad、Kaq为直轴、交轴电枢反应磁势折算系数,表示为
式中Ad1、Aq1为直轴、交轴电枢反应磁场幅度系数,a1为基波磁场幅度系数。
式中Ad、Aq、Bd、Bq、A1、B1可以通过查表获取。
额定负载时,各部分磁密变为原来的Ei倍,据此算得每极磁压降为
Fa=1.1(FδN+FrtN+FryN+FepN+FsyN) 44)
然后获得额定负载时励磁磁动势为
根据全电流定律,
NfIf=F 46)
式中Nf为一对极励磁线圈匝数。
励磁绕组总长度为
lf=2pNf(L+bep) 47)
设励磁绕组径向方向有N1匝,切向方向有N2匝,励磁导线半径为rf,则有如下约束条 件
电机的设计目标可以是重量、损耗及其他电气性能。耦合约束层可以描述电机的大多电 气性能,而在确定材料和尺寸的前提下电机重量可以由材料密度计算得到,在此不做赘述, 而主要对损耗进行建模。
电机的损耗包括机械损耗、铜损、铁损和杂散损耗四个方面:机械损耗包括电刷及轴承 摩擦损耗、通风损耗等;铜损为绕组的损耗;铁损为定子和转子铁芯的损耗。
根据能量转换平衡方法,轴承摩擦损耗Pf可以简化为
Pf=1.05FNs×10-4 49)
式中F为轴承受力矩,Nm;Ns为电机转速,r/min。
轴承承受力矩F包括转子重量G和单边磁拉力P0所引起的弯曲应力Fw以及转矩产生的 剪应力Ft,即
Ft即电机的额定转矩,为
Ft=9.55P2/Ns 51)
转子重量Gr和单边磁拉力P0所引起的弯曲应力Fw为
Fw=9.8(Gr+P0)L1L2/L 52)
式中L为两轴承间的距离(或转子铁心长度),m;L1、L2分别为左、右轴承到转子重心点的距离,m。
单边磁拉力P0产生的主要原因是转子相对于定子偏心,有
P0=πDL(Bδ/7000)2e/δ 53)
式中D为转子外径,cm;L,cm;δ为气隙长度,cm;e为定转子偏心值,cm;Bδ为气 隙磁密平均值,Gs。
空气摩擦损耗主要分为圆柱形的转子表面和轴向磁力轴承转子的盘状结构两侧空气摩擦 损耗两部分。
圆柱形转子的表面空气摩擦损耗可以表示为
PM=kMρairu3Am 54)
式中Am为空气与圆柱形的接触面积。
圆盘的表面空气摩擦损耗可以表示为
PS=kSρairD2u3/8 55)
总摩擦损耗为
Pair=ρairDu3(kMπL+kSD/8) 56)
式中摩擦系数kM和kS分别为圆柱和圆盘表面的摩擦系数,与电机几何形状和结构设计有 关;ρair为空气密度;D、u分别为转子的直径和圆周速度。
气压和温度对空气密度ρair的影响公式可表示为
式中p为大气压(mmHg),T为风温。
摩擦系数在很大范围内时雷诺数Re的函数,圆柱体在流体中旋转的Re为
式中r为圆柱体半径,m;ω为圆柱体旋转角速度,rad/s;ν为流体的运动粘度,m2/s。
根据Re的大小,可以求得摩擦系数k为
开路或者空载铁心损耗包活磁滞和涡流损耗,是仅在主励磁绕组激励的情况下,由电机 铁心中时变磁通密度所引起的损耗。
正常运行条件下,单位重量的钢片的涡流损耗pe可以近似表示为
式中ρ为钢片电阻率,ρFe为钢片密度,Bmax为磁密最大值,f为频率,Δ为钢片厚度。
在厚度一定的情况下,有
基于经验公式,单位重量的铁磁物质的磁滞损耗ph可表示为
式中Kh为比例系数,取决于铁芯特性和体积以及所采用的单位;n在1.5~2.5之间,电 机估算时经常取2。
单位重量的钢的铁损pFe为
对于不同含硅量的硅钢片,可以通过查表获得Ke、Kh的数值。
为了计算方便,如果一般情况下,不同的f和不同的Bmax,损耗系数为
式中p10/50表示当Bmax=1T,f=50Hz时钢单位损耗系数,可通过查表获得。
钢的基本铁耗应为
PFe=kapFeGFe 65)
式中ka为经验系数,GFe为钢的重量。
根据上述数学表达式,电机多学科模型的输入、输出接口参数如表1所示。
B:对不同部件多学科模型进行试验设计,提取表征多学科模型设计空间特征的输入、 输出参数样本,而后使用数据挖掘方法提取多学科模型的关键特征参数,利用机器学习算法 建立满足系统集成设计需求的近似模型。以电机多学科模型为例:
确定电机的设计空间如表1所示,采用基于最优拉丁超立方的试验设计技术生成电机多 学科模型的数据样本,步骤如图8所示,其结果如图9所示。
表1电机多学科模型输入输出参数及设定的设计空间
使用数据挖掘方法分析和提取电机多学科模型的关键特征参数,步骤如图10所示。
封装法将参数筛选视为搜索寻优问题,根据一定的算法和交叉验证获得的目标函数给出 最优的参数组合,基本框架如图11所示。以前向贪婪搜索(Greedy Search,GS)、粒子群搜 索(Particle Swarm Optimization,PSO)为搜索算法,以M5P树、随机森林(RandomForest, RF)和基于BP训练的多层感知机(Back-Propagation Multi-Layer Perceptron,BP-MLP)为预 测算法,以交叉验证获得的归一化均方根值作为寻优目标,对参数子集进行优化,获得筛选 结果如表2所示。
表2基于封装法的电机多学科模型关键特征参数提取结果
嵌入法利用某些机器学习算法本身即具有特征打分排序功能的特点,对数据样本进行训 练,获得各个特征的权值系数,据此筛选出对模型有重要意义的参数组合。使用二次回归 (Quadratic Regression,QR)、决策树(Decision Tree,DT)、RF、快速有监督(RapidSupervised method,RS)以及梯度(Gradient Boosting,GB)方法对电机多学科模型的参数的重要程度 开展分析,获得电机重量和损耗的关键参数如表3所示。
表3基于嵌入法的电机多学科模型关键特征参数提取结果
采用多数投票规则对表2和表3所示筛选结果进行综合,具体原理为
定义wv(i)为第i个参数的权重,wa(j)为第j个参数筛选方法的权重,则第i个参数的投票 得分Vti可以表示为
式中当第j个参数筛选方法排除了第i个参数时,Vij=0,否则Vij=1;N为需要的参数筛 选结果的组数。
部件参数的权重wv的选取取决于该参数对部件性能的影响程度;参数对部件性能越敏 感,说明该参数的重要程度越高。参数筛选方法的权重wa则取决于各方法对多物理模型的拟 合效果;拟合效果越好,说明该方法筛选的参数组合越精确。
因此,根据多数投票规则,当Vti满足下士时,第i个参数被确认为多物理模型的关键参 数之一。
设置各参数和各算法的权重均为1,则电机多学科模型的关键特征参数筛选过程如表4 所示,筛选结果如表5所示。
表4基于无偏好的多数投票对电机多学科模型关键特征参数的投票结果
表5电机多学科模型最终的关键特征参数筛选结果
利用机器学习算法建立满足系统集成设计需求的近似模型,步骤如图12所示。
采用k折交叉验证建立电机多学科模型的近似模型。k折交叉验证是指,将训练集分为k 份,轮流使用其中一份作为验证集、另外(k-1)份作为训练集,使用以上k种情况训练得到模 型,并使用整个训练集训练误差最小的模型。
采用表6所示预测算法,以表5中的关键参数分别建立电机重量和损耗的近似模型,经 过对这些模型决定系数、均方根误差、平均误差、最大误差以及消耗计算资源的综合对比, 选择径向基神经网络作为最终的预测算法。最后,使用一组新的数据样本对最终模型的性能 进行测试验证,电机重量的近似模型精度如图13所示,电机损耗的近似模型精度如图14所 示。
表6电机近似模型选用的预测算法
C:根据飞机功能及已知的设计条件(包括电源的类型和数目、汇流条的类型和数目、 负载的类型和数目等),基于平台的设计方法定义系统各部件功能,采用“自顶向下”和“自 底向上”的思想,使用定义的部件库自动生成系统网络拓扑,建立基于设计、功能及适航标 准的网络拓扑模型;而后根据设计经验和适航标准确定系统网络拓扑需要满足的约束条件(包 括系统网络拓扑设计经验产生的约束、需要满足的功能性约束、需要满足的可靠性性约束等), 建立约束条件关于网络拓扑的函数关系,并借助约束满足算法求解满足这些约束条件的网络 拓扑可行方案;接着,根据建立的网络拓扑模型计算网络的基本可靠度和各负载节点对应的 任务可靠度,形成基本可靠度和任务可靠度关于网络拓扑的函数关系,作为网络拓扑可靠性 的衡量依据,使用搜索算法求解获得满足可靠性要求的最优网络拓扑方案。
网络拓扑优化设计以2个230V 400Hz的发电机、1个400V辅助动力装置燃料电池、3个230V 400Hz交流汇流条、2个270V直流汇流条和3个28V直流汇流条以及8个分别挂在 各个汇流条上的储能组成的分布式系统为例进行说明:
定义飞机供电系统平台如表7所示。
表7飞机供电系统网络拓扑设计平台
部件 | 输入/输出接口数目 |
发电机 | 1 |
燃料电池 | 1 |
锂电池 | 1 |
功率变换器 | 2 |
开关 | 2 |
电缆 | 2 |
汇流条 | 1-6 |
根据图论,将系统部件分为两类:节点包括电源、汇流条和储能,有向边为变换器、开 关、电缆构成的支路。则可以根据各个节点的重要程度和供电体制使用有向边将其连接城一 个系统,该系统可以使用节点邻接矩阵表示节点之间的连接关系。
根据GJB181A、CCAR25-R4中国民用航空规章第25部等标准,建立系统约束条件的数 学模型。
航空规定,每个电源都必须只连接一个汇流条,实际电源连接参数数目等于电源数目Ns, 因此电源在邻接矩阵中所在的行/列均只能有Ns个非零元素,并且每行/列必须有一个非零元 素,即
为了保证系统的供电裕度,多个电源不能连接到同一个汇流条上,因此
汇流条连接参数主要指汇流条节点之间的连接关系和能量可以流动的方向,满足
根据分布式电源的基本思想,假设在满足其他约束的条件下,每个汇流条均挂载一个储 能系统;若容量配置过程中优化得到某处储能容量为零,则说明该处无储能。
根据飞机电网约定俗成的设计经验,电源节点和负载节点不能够直接相连,即
负载或汇流条节点必须存在输入,即
网络拓扑结构满足故障重构后的负载的供电可靠性要求,保证系统发生故障时负载的可 靠工作,即
式中No_of_failure为系统中发生故障的部件的数目;imax为系统供电裕度要求,根据适航 标准取imax=3。
系统网络拓扑的供电能力约束是一个动态约束,设计人员需要依次对所有可能的故障情 况下系统供电路径的有无进行考察。
在系统故障重构过程中,系统满足:
系统工作时,任意节点i和j有
cij+cji=0 74)
交流电源不能并联供电,因此根据电网络理论,源节点和交流汇流条节点组成的网络支 路数不能大于交流汇流条数目Nla,即
当飞机主电源均无故障时,若每个汇流条节点均有电,则辅助电源不启动。
飞机中还存在可以正常工作的主电源和辅助电源时,应急电源不启动。
为充分利用系统支路的容错能力,系统各支路在遍历的过程中必须至少出现一次,从而 保证每条支路在故障重构中都能够起到作用。
式中NB为电网中除储能支路外其他支路的条数;Bri为电网第i个支路;ScenPj为由故障 产生的第j个工作状态;No_Scen为该系统可能存在的工作状态数目。
根据适航标准,系统必须满足:
起飞前,供电系统的失效率必须不高于1.33×10-3,飞机才能够派签。因此,系统基本可 靠度对应的失效率λbasic满足
λbasic≤1.33×10-3 79)
飞机出现任何降低飞机能力的失效情况的几率对应的失效率不大于10-4。因此,重要汇 流条的供电可靠度对应的失效率λESS满足
λESS≤10-4 80)
飞机出现次要的、会在一定程度上降低功能能力和安全余度的故障对应的失效率不大于 10-3。因此,一般汇流条的供电可靠度对应的失效率λGEN满足
λGEN≤10-3 81)
使用智能优化算法对该问题进行求解,可以获得96组可行解。
建立系统基本可靠度和任务可靠度对应的失效率关于网络邻接矩阵的数学模型。使用遍 历法从96组可行解中获得可靠性的帕累托最优解如图15所示,图16给出了其中一个解。
D:根据前面获得的系统最佳网络拓扑方案,制定网络重构策略,生成各个方案可能存在 的工作状态集合,删除其中冗余和不具有代表性的状态;然后制定系统能量管理策略,建立 飞机供电系统动态功率模型,根据系统的功率需求,使用多目标优化算法对不同工作状态下 系统的燃油消耗、系统损耗和储能重量进行优化,获得能满足功率需求的最佳能量调度方案; 最后对同一网络拓扑方案下多个工作状态的优化结果进行综合,获得能够满足大部分工作状 态工作需求的容量配置方案。
采用多目标优化算法对飞机供电系统不同工作状态在整个飞行过程中的动态能量调度进 行优化规划,即
式中No_Scen为飞机供电系统工作状态的数目;PwrSci为系统第i个工作状态在整个飞 行过程中能量优化调度的设计结果;pi(t)为系统第i个工作状态系统在t时刻的支路功率矩阵, 包含负载节点的动态需求;wF i为系统第i个工作状态的系统燃油损耗;wES i为系统第i个工 作状态的系统在整个飞行过程中所需的储能的重量;pls i为系统第i个工作状态的系统损耗; Pupper和Plower为能量分配策略给各支路设置的功率上、下限;KVL/KCL为系统在不同工作状 态下需要遵守的基尔霍夫电压/电流定律。
下面以图17所示2个230V 400Hz的主发电机、1个400V辅助动力装置燃料电池、2个270V直流汇流条、2个28V直流汇流条、1个900V直流汇流条以及2个240V、2个24V、1 个540V的储能组成的分布式系统为例对供电系统容量优化配置进行说明:
首先以燃油消耗和系统损耗为输出,建立单个工作状态下系统的能量优化调度模型。
系统第i个工作状态的燃油损耗wF i包括发电机消耗的航空燃油的重量wAo i和燃料电池消 耗的氢燃料的重量wH i,可以表示为
式中wAo_j i和wH_j i分别为第i个工作状态下第j个发电机消耗的航空燃油的重量和第j个 燃料电池消耗的氢燃料的重量。
wAo_j i可以根据发电机轴功率和发动机效率ηEng关于时间的积分进行行计算,即
式中T为整个飞行时间;pG_j i(t)和plG_j i(t)为分别为第j个发电机在第i个工作状态下t时 刻的输出功率和损耗;be为发动机燃油率。
以整个飞行过程的平均损耗表示工作状态i下系统损耗pls i,即
式中T为飞行总时间;支路损耗pljk(t)可表示为
当j=k时pljk(t)表示为节点i设备的损耗;当j≠k时,pljk(t)表示为支路jk的损耗。若cjk(t)=0, 两节点无连接,pljk(t)=0。否则,若两节点电压Vj=Vk,则支路jk的损耗包括由开关损耗plSwi_jk(t)、 线缆损耗plCa_jk(t);若两节点电压Vj≠Vk,则支路jk的损耗包括变换器损耗plCon_jk(t)、开关损 耗plSwi_jk(t)、线缆损耗plCa_jk(t)组成。
式中plSwi_j,k、plCa_j,k分别为支路jk上节点j端开关、线缆损耗,plSwi_k,j、plCa_k,j分别为支 路jk上节点k端开关、线缆损耗。线缆在不同位置产生的重量和损耗的差异是由其电压等级 的不同导致的。
在系统动态运行过程中,根据各支路实时功率,可计算得到第i个储能系统的动态输入/ 出功率pc_i(t),结合充放电器的动态损耗plC_i(t)。单次飞行过程中,储能的容量的变化capi(t) 为
式中ΔSOC为储能最大充/放电深度;Capi(0)为飞行初始时刻储能的能量,由于储能在飞 行刚开始时主要作用是吸收多余功率,因此设其为零。
则整个飞行过程对储能容量的需求Cap则可以通过式90)计算得到。
然后考虑飞机供电系统不同的故障情况,生成系统工作状态集合。
设定系统故障重构过程必须满足
式中R(pathl)为路径pathl的供电可靠度;pl(pathl)为路径pathl上的损耗;No_Swi(pathl) 为路径pathl上的接触器开关次数。
任一供电路径的可靠度为该路径上所有支路可靠度之积,考虑各支路可靠度相差不大(一 般最多相差一个变换器或者开关的可靠度),假设所有支路可靠度对应的失效率近似相等于 λBr,则R(pathl)可表示为
式中pathl_i为pathl上第i个支路;Ri为pathl上第i个支路的可靠度;Num为pathl上支 路的数目。
为简化计算,在生成工作状态时忽略供电路径上开关和变换器的损耗,仅考虑l节点供 电路径的线路损耗,则
式中ρCu为铜导线的电阻率;i_s和i_e分别为路径pathl第i个支路pathl_i起点和终点的 编号;Ii_s,i_e为支路pathl_i上传输的功率;Si为支路pathl_i线缆的截面积,根据线缆载流量进 行设计,与载流量IC和Ii_s,i_e近似有Ii_s,i_e=ICSi的关系。
对于一个确定的系统而言,导线载流量IC仅取决于温度,因此可认为所有导线的载流量 相同,因此
因此,将损耗pl分解为供电路径长度和供电路径上节点的其他支路数目,即
故
接着对产生的大量工作状态进行缩减,采用基于支路负荷极限和信息量的方法,选择系 统每个支路的负荷极限情况的工作状态中包含信息量最多、最具代表性的工作状态作为该支 路的工作状态。则最能代表支路i的工作状态求解为
式中Ii_max和Ii_min分别为无储能时系统中支路i可能的最大、最小负荷对应的工作状态集 合;Nrdi_j为支路i的极限负荷工作状态集合中第j个工作状态在所有支路极限负荷工作状态 集合中出现的次数,则Nrdi_j越大,工作状态j涵盖的极限负荷支路越多,其包含的信息量越 大;Nbi_j为支路i在第j个工作状态下系统运行的支路数,则Nbi_j越小,支路i在该工作状 态中的地位越重要,因此该工作状态对于支路i越具有代表性。
式中Pwrij为无储能时系统中支路i在工作状态j下工作的功率,即该矩阵第i行、第j 列元素,包括根据功率包线最大值计算的Pwrmax_ij和根据功率包线最小值计算的Pwrmin_ij,应 选用Pwrmax_ij;Imax为Ii_max构成的矩阵,表示支路与负荷极限状态的对应关系;Nrdi_j即统计 Imax中出现次数最多的非零元素,使用count表示。
最后,以系统可靠性、重量和损耗为目标,建立以部件容量、节点电压为输入的系统容 量配置模型,采用多目标优化算法优化系统网络拓扑、容量和节点电压,保证系统可靠性最 高、重量最小、损耗最小,即
式中Rbasic为系统基本可靠度;Rtask为系统任务可靠度;C为系统网络拓扑邻接矩阵;WS和PlS分别为系统重量和损耗;P为系统各支路部件容量;V为系统各节点电压;No_Scen为 飞机供电系统工作状态的数目;PwrSci为系统第i个工作状态在整个飞行过程中能量优化调 度的设计结果。
可靠性包括基本可靠性和任务可靠性,以其对应的失效率进行计算。
飞机供电系统基本可靠度对应的失效率λbasic可表示为系统所有部件失效率之和。
其中,λij为支路ij的失效率,为连接节点i和节点j的支路上部件失效率之和,即
其中,λSwi为开关的失效率;λCon为变换器的失效率。当i=j时,λij为节点i所代表的部件 的失效率;当i≠j时,根据节点i和j的类型和电压等级Vi和Vj,确定支路ij上的部件(开关、 变换器)。不计电缆的失效率,并默认相同类型的装置设备的失效率相同。
飞机供电系统任务即保证每个负载节点的供电,因此任务可靠性根据系统中每个负载节 点的供电路径的可靠性进行计算。基于深度优先遍历搜索,求取负载节点l的供电路径集合 pathl,并根据全概率公式,计算获得负载节点l的任务可靠性Rtask_l。
其中,Npath为最小路集合pathl中路径的条数;pathl_i为pathl中第i条路径。
但全概率公式计算量巨大,为简化计算,采用Aggrawal提出的“删去留下”法对最小路 集进行不交化处理,具体过程不再赘述。
重量Wb包括发动机当量重量和电网各支路的重量,可以表示为
式中WEngin为发动机当量重量,可表示为
式中NG为系统中发电机的数目;ηG_i为第i个发电机的效率;ηEngin为发动机的效率,取 60%;SpPEngin为发动机的比功率,取11kW/kg。
根据节点i和j的不同类型,支路ij的重量Wij可表示为
式中当i=j时,Wij表示为节点i设备的重量,可能是发电机WG、燃料电池WFC、汇流条负载WEL或者储能的重量WES。当i≠j时,Wij表示为支路ij的重量。若cij=0,两节点无连接,Wij=0。否则,当节点i、j均非储能节点时,若两节点电压Vi=Vj,则支路ij的重量由开关重量WSwi_ij、线缆重量WCa_ij组成;若两节点电压Vi≠Vj,则支路ij的重量由变换器重量WCon_ij、开关重量WSwi_ij、线缆重量WCa_ij组成。当节点i、j有储能节点时,支路ij的重量由充放电器构成。
线缆和开关的布局需要分具体情况讨论,有
式中WSwi_i,j和WCa_i,j分别为支路ij上节点i端开关和线缆重量;WSwi_j,i和WCa_j,i分别为支 路ij上节点j端开关和线缆重量。当节点i和j供电体制相同时,WCa_i,j=WCa_j,i;否则,需要 对电缆重量与供电体制的关系进行分析,选择重量较小的电缆敷设方案。
燃料和储能系统的重量WFuel和WESS需要根据系统不同工作状态下能量优化调度的结果 确定,即
式中wESS i和wF i为系统第i个优化工作状态所需要的储能系统和燃料的重量。
储能系统的重量wESS i则根据储能的容量需求进行计算,即
由于不同工作状态下系统效率不同,以不同工作状态系统损耗的加权平均进行衡量,即
式中plS k为工作状态k下系统的总损耗;Posk为工作状态k发生的概率。
假设工作状态k是由于p条支路brf_1、…、brf_p发生故障时,系统重构后由q条支路brn_1、…、brn_q进行供电的系统,则
式中Rbrf_i为第i个故障支路的可靠度;Rbrn_i为第i个工作支路的可靠度。
在上述模型的基础上,采用图5所示优化框架,两个优化器均选择NSG-II算法,开展运 算,获得支路容量设计结果如表8所示,储能容量设计结果如表9所示。
表8飞机供电系统支路容量优化配置结果
支路编号 | 支路容量/kW | 支路编号 | 支路容量/kW |
1-4 | 27.58 | 2-5 | 27.62 |
3-8 | 18.7 | 4-5 | 27.5 |
4-6 | 3.1 | 4-7 | 2.7 |
4-8 | 32 | 5-6 | 3.4 |
5-7 | 3.4 | 5-8 | 59 |
8-6 | 9 | 8-7 | 8.8 |
表9飞机供电系统储能容量优化配置结果
Claims (8)
1.一种飞机供电系统智能化集成设计系统,其特征在于包括:
多学科建模模块,针对飞机供电系统不同部件内部参数的多物理耦合关系建立飞机供电系统发电机、燃料电池、蓄电池、变换器、开关、电缆等各部件的多学科模型,并设定各部件多学科模型的参数;
代理建模模块,采用试验设计方法生成多学科模型的数据样本,采用数据挖掘方法分析和提取多学科模型的关键特征参数,采用机器学习算法对数据样本进行训练,建立各部件的近似模型;
网络拓扑优化设计模块,根据已知的设计条件建立飞机供电系统网络拓扑模型,并根据各部件的近似模型建立系统网络拓扑设计目标和约束条件,使用多目标优化算法求解获得满足设计要求的最优网络拓扑方案;
容量优化配置模块,根据系统最优网络拓扑方案生成系统可能存在的工作状态集合,在此基础上建立每个工作状态下系统能量调度模型,借助多目标优化算法对系统进行能量优化调度的计算求解,对不同工作状态下系统能量调度方案进行综合,获得对应网络拓扑满足大部分工作状态工作需求的容量配置方案并建立以部件容量、节点电压为输入的容量配置模型,借助多目标优化算法对网络拓扑、部件容量、节点电压进行优化,获得性能最佳的系统容量配置方案。
2.根据权利要求1所述的飞机供电系统智能化集成设计系统,其特征在于:多学科建模模块利用各部件多物理耦合原理建立各部件的多学科模型,多物理耦合指多种物理领域之间参数的耦合,多学科模型具体是指以多物理耦合参数为输入,以性能参数为输出的模型。
3.一种飞机供电系统智能化集成设计方法,其特征在于包括以下步骤:
步骤一:针对飞机供电系统不同部件内部参数的多物理耦合关系建立飞机供电系统各部件的多学科模型,并设定各部件多学科模型的参数;
步骤二:采用试验设计方法生成多学科模型的数据样本,采用数据挖掘方法分析和提取多学科模型的关键特征参数,采用机器学习算法对数据样本进行训练,建立各部件的近似模型;
步骤三:根据已知的设计条件建立飞机供电系统网络拓扑模型,并根据各部件的近似模型建立系统网络拓扑设计目标和约束条件,使用多目标优化算法求解获得满足设计要求的最优网络拓扑方案;
步骤四:根据步骤三获得的系统最优网络拓扑方案生成系统可能存在的工作状态集合,在此基础上建立每个工作状态下系统能量调度模型,借助优化算法对系统进行能量优化调度的计算求解,而后对不同工作状态下系统能量调度方案进行综合,获得对应网络拓扑满足不同工作状态工作需求的容量配置方案,并建立以部件容量、节点电压为输入的容量配置模型,借助多目标优化算法对网络拓扑、部件容量、节点电压进行优化,获得性能最佳的系统容量配置方案。
4.根据权利要求3所述的飞机供电系统智能化集成设计方法,其特征在于步骤三中,以图论、电网络理论为基础,以基于平台的设计方法建立飞机供电系统网络拓扑模型,而后根据设计经验和适航标准确定网络拓扑需要满足的约束条件,建立约束条件关于网络拓扑的函数关系,并使用约束满足算法求解满足这些约束条件的网络拓扑可行方案;根据建立的网络拓扑模型计算该网络的基本可靠度和各负载节点对应的任务可靠度,形成基本可靠度和任务可靠度关于网络拓扑的函数关系,作为网络拓扑可靠性的衡量依据,使用搜索算法求解获得满足可靠性要求的最优网络拓扑方案。
5.根据权利要求3所述的飞机供电系统智能化集成设计方法,其特征在于步骤四中,基于飞机供电余度要求遍历飞机供电系统所有可能的故障,根据指定的网络重构策略生成这些故障发生时系统的供电路径,生成供电系统网络拓扑的工作状态集合。
7.根据权利要求3所述的飞机供电系统智能化集成设计方法,其特征在于,采用多目标优化算法对飞机供电系统不同工作状态在整个飞行过程中的动态能量调度进行优化规划,即
For i=1:No_Scen
End
式中No_Scen为飞机供电系统工作状态的数目;PwrSci为系统第i个工作状态在整个飞行过程中能量优化调度的设计结果;pi(t)为系统第i个工作状态系统在t时刻的支路功率矩阵,包含负载节点的动态需求;wF i为系统第i个工作状态的系统燃油损耗;wES i为系统第i个工作状态的系统在整个飞行过程中所需的储能的重量;pls i为系统第i个工作状态的系统损耗;Pupper和Plower为能量分配策略给各支路设置的功率上、下限;KVL/KCL为系统在不同工作状态下需要遵守的基尔霍夫电压/电流定律。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011092419.6A CN112257174A (zh) | 2020-10-13 | 2020-10-13 | 一种飞机供电系统智能化集成设计系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011092419.6A CN112257174A (zh) | 2020-10-13 | 2020-10-13 | 一种飞机供电系统智能化集成设计系统及方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112257174A true CN112257174A (zh) | 2021-01-22 |
Family
ID=74243205
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011092419.6A Pending CN112257174A (zh) | 2020-10-13 | 2020-10-13 | 一种飞机供电系统智能化集成设计系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112257174A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114615149A (zh) * | 2022-05-12 | 2022-06-10 | 南昌航空大学 | 一种飞行器多动力系统数据交互网络结构优化方法 |
CN114880959A (zh) * | 2022-07-11 | 2022-08-09 | 广东电网有限责任公司佛山供电局 | 一种用于建筑能耗混合模型的输入变量获取方法和系统 |
US11515584B1 (en) * | 2021-11-10 | 2022-11-29 | Beta Air, Llc | Systems and methods for aircraft power management and distribution |
CN115470602A (zh) * | 2022-09-23 | 2022-12-13 | 中国人民解放军国防科技大学 | 一种考虑节点差异化抗毁需求的供电网络拓扑优化方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107317874A (zh) * | 2017-07-26 | 2017-11-03 | 中国航空工业集团公司西安飞机设计研究所 | 一种机载多拓扑双交换结构的网络设计方法 |
-
2020
- 2020-10-13 CN CN202011092419.6A patent/CN112257174A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107317874A (zh) * | 2017-07-26 | 2017-11-03 | 中国航空工业集团公司西安飞机设计研究所 | 一种机载多拓扑双交换结构的网络设计方法 |
Non-Patent Citations (2)
Title |
---|
ZEHUA DAI等: "Multi-level modeling methodology for optimal design of electric machines based on multi-disciplinary design optimization", 《ENERGIES》, 1 November 2019 (2019-11-01), pages 1 - 26 * |
王莉等: "电气化飞机电力系统智能化设计研究综述", 《航空学报》, 25 February 2019 (2019-02-25), pages 5 - 19 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11515584B1 (en) * | 2021-11-10 | 2022-11-29 | Beta Air, Llc | Systems and methods for aircraft power management and distribution |
CN114615149A (zh) * | 2022-05-12 | 2022-06-10 | 南昌航空大学 | 一种飞行器多动力系统数据交互网络结构优化方法 |
CN114615149B (zh) * | 2022-05-12 | 2022-08-02 | 南昌航空大学 | 一种飞行器多动力系统数据交互网络结构优化方法 |
CN114880959A (zh) * | 2022-07-11 | 2022-08-09 | 广东电网有限责任公司佛山供电局 | 一种用于建筑能耗混合模型的输入变量获取方法和系统 |
CN114880959B (zh) * | 2022-07-11 | 2022-12-30 | 广东电网有限责任公司佛山供电局 | 一种用于建筑能耗混合模型的输入变量获取方法和系统 |
CN115470602A (zh) * | 2022-09-23 | 2022-12-13 | 中国人民解放军国防科技大学 | 一种考虑节点差异化抗毁需求的供电网络拓扑优化方法 |
CN115470602B (zh) * | 2022-09-23 | 2023-06-16 | 中国人民解放军国防科技大学 | 一种考虑节点差异化抗毁需求的供电网络拓扑优化方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112257174A (zh) | 一种飞机供电系统智能化集成设计系统及方法 | |
Wheeler et al. | Electric/hybrid-electric aircraft propulsion systems | |
Xu et al. | Multi-objective dynamic VAR planning against short-term voltage instability using a decomposition-based evolutionary algorithm | |
Sarigiannidis et al. | Fault tolerant design of fractional slot winding permanent magnet aerospace actuator | |
Zhang et al. | Steady‐state performance evaluations of three‐phase brushless asynchronous excitation system for aircraft starter/generator | |
US8036858B2 (en) | Statistical method for electric machine optimization used for electric drives | |
Fletcher et al. | Modeling and simulation enabled UAV electrical power system design | |
Liu et al. | Mechanism and characteristic analyses of subsynchronous oscillations caused by the interactions between direct‐drive wind turbines and weak AC power systems | |
Madonna et al. | Thermal analysis of fault‐tolerant electrical machines for more electric aircraft applications | |
Moayed et al. | Presenting a stochastic model of simultaneous planning problem of distribution and subtransmission network development considering the reliability and security indicators | |
CN106026083A (zh) | 一种基于matlab的电网暂态稳定分析方法 | |
Darabi et al. | Employing adaptive particle swarm optimization algorithm for parameter estimation of an exciter machine | |
Liu et al. | Start-up sequence of generators in power system restoration avoiding the backtracking algorithm | |
Dong et al. | Reduction and modelling method of large‐scale alternating current/direct current power systems for electromagnetic transient simulation | |
Scheidler et al. | Methodology for electromagnetic optimization of a partially superconducting 1.4 MW electric machine for electrified aircraft propulsion | |
Ismagilov et al. | Topology selection of the high-speed high-voltage PMSM for aerospace application | |
CN108539737A (zh) | 一种工程实用的电力系统动态等值优化方法 | |
Chomat et al. | Analysis of magnetic field distribution in induction machine with combined star‐delta stator winding | |
Dalali et al. | Modified Thevenin‐based voltage instability indicator and load shedding approach for MCF connected network | |
Zhou et al. | Dynamic equivalent method of interconnected power systems with consideration of motor loads | |
Li et al. | Economic dispatch methods for smart grid based on improved SPEA2 and improved NSGA2 | |
DALCALI et al. | Electromagnetic Energy Harvester Design for Power Transmission Line | |
Ismagilov et al. | Investigation of an energy converter with an amorphous alloy magnetic core for a 270 V generation channel of an aircraft, taking into account climatic, vibration and electromagnetic conditions | |
Yan et al. | Influence of different time series on load recovery following faults in the vessel regional distribution system | |
Afzalan et al. | DG ALLOCATION AND SIZING USING MSFLA TO REDUCE LOSSES AND IMPROVE VOLTAGE PROFILE IN DISTRIBUTION NETWORK. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |