CN111666630B - 一种煤矿井下无人机结构的多目标智能优化设计方法 - Google Patents
一种煤矿井下无人机结构的多目标智能优化设计方法 Download PDFInfo
- Publication number
- CN111666630B CN111666630B CN202010458424.8A CN202010458424A CN111666630B CN 111666630 B CN111666630 B CN 111666630B CN 202010458424 A CN202010458424 A CN 202010458424A CN 111666630 B CN111666630 B CN 111666630B
- Authority
- CN
- China
- Prior art keywords
- aerial vehicle
- unmanned aerial
- coal mine
- underground
- model
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 78
- 238000005457 optimization Methods 0.000 title claims abstract description 74
- 238000013461 design Methods 0.000 title claims abstract description 73
- 239000003245 coal Substances 0.000 title claims abstract description 69
- 238000004458 analytical method Methods 0.000 claims abstract description 19
- 238000004088 simulation Methods 0.000 claims abstract description 15
- 230000006870 function Effects 0.000 claims abstract description 12
- 238000011156 evaluation Methods 0.000 claims abstract description 11
- 238000012795 verification Methods 0.000 claims abstract description 8
- 230000007246 mechanism Effects 0.000 claims abstract description 4
- 238000004519 manufacturing process Methods 0.000 claims description 13
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 8
- 239000000463 material Substances 0.000 claims description 8
- 230000001052 transient effect Effects 0.000 claims description 8
- 238000012549 training Methods 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 4
- 230000009471 action Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 230000035772 mutation Effects 0.000 claims description 3
- 239000002994 raw material Substances 0.000 claims description 3
- 238000010801 machine learning Methods 0.000 abstract 1
- 238000012545 processing Methods 0.000 abstract 1
- 239000003795 chemical substances by application Substances 0.000 description 6
- 239000002131 composite material Substances 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical group C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 229920000049 Carbon (fiber) Polymers 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000004917 carbon fiber Substances 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000011152 fibreglass Substances 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000007704 transition Effects 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种煤矿井下无人机结构的多目标智能优化设计方法,属于无人机结构设计方法。无人机设计:根据井下无人机飞行环境特点和巡航要求,选取无人机结构设计的决策变量,确定优化目标和约束条件,构建井下无人机结构设计的多目标优化模型;考虑到无人机飞行稳定性这一优化目标很难建立数值模型,基于决策变量取值,获得代表个体及其飞行稳定性的数值模拟评价,进而采用机器学习方法,构建其代理模型,形成飞行稳定性优化目标评价函数;采用NSGA‑II多目标优化方法,结合约束处理机制,实现无人机结构的寻优,并对获得的最优结构进行仿真验证分析。本方法用于煤矿井下无人机结构设计,有效保证了无人机结构的合理性与科学性,保障了井下作业的安全性和经济性。
Description
技术领域
本发明涉及一种无人机结构设计方法,特别是一种煤矿井下无人机结构的多目标智能优化设计方法。
背景技术
目前,我国煤炭资源的分布存在浅层埋藏量不足特性,煤炭开采逐步向深部开发利用发展,开采难度越来越高,出现瓦斯突出、冲击地压等灾害的可能性也随之升高。当出现作业事故时,必须对灾后的井下环境和人员进行及时准确探测。然而,煤矿井下复杂地生产条件,使探测救援的难度大大增加;传统的人工施救策略,产生二次意外事故的可能性较高,对矿井安全隐患的危险倍增。
无人机作为一款高度智能机器人,是一种通过无线电技术、设备和自控程序实现远端智能调控的非载人飞行器。无人机具有远程操控或自主能力,不存在执行人员的伤亡隐患,成为有效的救援探测设备。文献(张军红,李国民.轻型无人机复合材料机翼设计分析与验证.西安航空技术高等专科学校学报,2013,31(3):33-36)基于有限元模型,根据机翼的静强度和应力分布,选取轻型无人机的机翼构件布局。该方法对中小型复合材料无人机结构设计和静强度试验具有一定的参考价值。文献(Lee M K,Cho C M,Jang S Y.HALEUAV.Composite wing structure design.Advanced Materials Research,2010,123-125:105-108) 选取具有最佳无人机质量的复合材料,并对所设计的无人机机翼进行有限元分析。文献(王一凡,孙刚,张淼.考虑转捩点约束的自然层流翼型变弯度设计.航空计算技术,2018,48(2):33-36)针对大型民用航空飞机的机翼,采用代理模型和遗传算法,提出针对自然层流翼型的变弯度优化设计方法。实验分析表明在巡航范围内的不同升力系数下,随着翼型弯度的变化,翼型表面仍然能够保持一定程度的层流区域,使整个巡航段中翼型的气动效率得到提升。
针对煤矿井下四旋翼无人机的智能优化方法研究较少,大多数研究工作主要是采用数值模拟的方法,对四旋翼无人机结构进行有限元分析与模态分析。文献(刘峰,高鸿渐,喻辉,等.基于有限元的四旋翼无人机碳纤维结构优化设计与固有模态分析.玻璃钢/复合材料, 2017,4:17-23)针对消费级的无人机,采用有限元分析方法,基于最大应力强度准则,设计质量轻、结构简单的碳纤维四旋翼无人机,分析了所设计结构的稳定性,计算了初始结构失稳临界载荷和失稳模态,基于最小重量要求,对无人机初始结构的碳纤维铺层进行了优化设计。
发明内容
本发明所要解决的技术问题是:提出一种煤矿井下无人机结构的多目标智能优化设计方法,旨在获得一组兼顾无人机质量、无人机升力和无人机成本的最优无人机制造方案。
本发明的技术解决方案:无人机多目标智能优化设计方法:建立煤矿井下无人机结构的多目标优化模型,并采用多目标优化方法求解该结构优化设计模型,为煤矿井下巷道狭闭空间内实现巡航的无人机设计方案,其步骤如下:
Step1:建立煤矿井下无人机结构方案的升力代理模型。
井下无人机由机架和动力系统组成,整机具有对称性;在机架和动力系统的共同作用下,无人机的升力由单一桨叶上提供的拉力差决定;无人机机架结构与动力系统包含多个参数,其与无人机升力之间存在高度非线性关系,且无法建立机理模型,采用辨识方法,构建无人机升力的代理模型,用于后续智能优化过程中的个体适应度评价。
Step2:构建兼顾无人机质量、成本和升力的多目标优化设计模型。
为了满足煤矿井下探测和巡航工作需要,同时满足无人机的飞行安全性和制造经济性,无人机结构方案的评价指标包括两个部分;其一,采用无人机质量f1(x)和升力f2(x),评估无人机的飞行平稳性,其中,无人机升力采用代理模型,无人机质量通过机架结构参数和机架材料加以描述;其二,采用无人机成本f3(x)评价所设计的无人机结构经济性,选取无人机机架和动力系统所需材料的总价作为无人机成本;确定决策变量及其约束范围,兼顾无人机质量、升力和成本三个优化目标,建立煤矿井下无人机结构的三目标优化设计模型。
Step3:煤矿井下无人机结构的多目标优化设计方法。
采用多目标智能优化算法,求解煤矿井下无人机结构设计三目标优化模型,得到Pareto 最优解集,实现井下无人机结构参数的优化设计。
Step4:数值模拟验证获得的最优结构方案性能。
为了验证获得的最优结构方案是否能够满足实际工程需要,从Pareto最优解集中选取满足成本预算需求的若干结构方案,采用Fluent数值模拟软件,对无人机升力进行分析验证,并对无人机机架进行模态分析,验证机架稳定性和飞行安全性。
所述的步骤1中,建立煤矿井下无人机结构方案的升力代理模型,采用以下步骤获得:
Step1:选取无人机机架和动力系统关键参数,确定煤矿井下无人机结构优化的决策变量,记为x={x1,x2,x3,x4,x5,x6,x7,x8};其中,x1表示电机座直径,x2表示机架直径,x3表示框架宽度,x4表示起落架高度,x5表示厚度,x6表示高度,x7表示支架个数,x8表示桨叶尺寸。
Step2:根据无人机原材料制造限制,x5,x7,x8采用离散整数值;其余变量采用连续值;确定各个决策变量的取值范围,记为 x5∈{1,1.5,2,2.5},x7∈{2,3,4,5,6,7,8},x8∈{127.0,152.4}。
Step3:采用拉丁超立方抽样法,产生无人机结构样本数据,记为D=[xi];其中,xi,i=1,2,…,N为第i个样本,N为样本规模。
Step4:采用Fluent数值模拟软件,对相应结构方案建立三维网格和有限元分析模型,获得每种结构方案的无人机升力,记为U(xi),xi∈D。
Step6:采用由Kriging法和径向基函数构成的集成学习模型,训练样本数据,辨识无人机结构参数和升力之间的非线性映射关系,构建高精度的无人机升力代理模型f2(x);其中,各个子模型的输出通过加权方法加以集成,权值根据其训练误差动态修正;令为第l个子模型的训练误差,为无人机升力估计值,则第k个子模型权重wk记为:
所述的步骤2中,构建兼顾无人机质量、成本和升力的多目标优化设计模型,采用以下步骤获得:
Step1:记ρt为无人机机架密度,ρs为无人机支撑架密度,d为无人机支架宽度,M0为桨叶的单位质量,ct为无人机机架造价,cs为无人机支架造价,kf为桨叶造价,建立无人机质量f1(x)和无人机成本f3(x)的目标函数分别为:
Step2:构建煤矿井下无人机结构的三目标优化设计模型,记为
F(x)={minf1(x),maxf2(x),minf3(x)}。
所述的步骤3中,煤矿井下无人机结构的多目标优化设计方法,采用以下步骤获得:
Step1:设定种群规模为N,最大进化代数为Gen;设置优化算法的初始进化代数t=0,设置初始种群P(0)。
Step2:对个体进行适应度评价,并根据个体的适应度值进行非支配排序。
Step3:依据个体的序值和拥挤距离,选择非占优个体。
Step4:通过交叉和变异操作,生成下一代种群P(t+1)。
Step5:判断是否达到最大迭代次数,若满足,则输出Pareto最优解集;否则,t=t+1,跳转到Step2。
所述的步骤4中,数值模拟验证获得的最优结构方案性能,采用以下步骤获得:
Step1:从Pareto最优解集中,分别选取具有最小成本、最小质量、最大升力和具有折中目标值的4种无人机结构方案。
Step2:结合SolidWorks和Fluent仿真软件,仿真分析无人机结构方案中单一桨叶的拉力,并将其作为机架的反向作用力,求解无人机的瞬态模态响应。
Step3:选取可以满足无人机质量、升力和成本实际需求的最优无人机结构方案。
有益效果,由于采用了上述方案,本发明面向煤矿井下这一特有狭闭空间,采用多目标优化设计方法,寻优获得可以满足井下平稳飞行的小型四旋翼无人机最佳设计方案,实现矿井环境的探测和信息的获取,提高煤矿井下灾后救援的有效性。四旋翼无人机通常具有结构对称性,为井下四旋翼无人机方案设计提供了有利条件。所设计的四旋翼无人机不仅可以避免救援人员在灾难现场靠近未知区域采集数据时发生二次灾害的危险,降低井下盲区救援的危险灾害系数,为救援决策提供科学有利的依据。
本发明有机融合数据驱动的代理模型技术和智能优化方法,综合考虑煤矿井下巷道用于巡航和探测的无人机飞行安全性、平稳性和制造经济性,兼顾无人机质量、升力和成本三个目标,采用先进的多目标进化优化方法,实现无人机结构方案的寻优。
构建的煤矿井下无人机升力代理模型可以替代耗时的数值模拟评价过程,对优化过程中的个体性能进行评价,有效降低计算成本;多目标优化模型兼顾无人机质量、成本和升力,基于此模型获得的无人机结构方案符合煤矿井下实际运行环境需求,有利于保障其在煤矿井下巡航和探测任务的顺利实施,还可以为煤矿井下灾后救援工作的安全推进和井下环境状态探测提供有力的实现基础。
无人机结构设计目的是在保证无人机结构稳定性的前提下,尽可能的降低无人机质量和成本,提高无人机的拉力,从而实现煤矿井下无人机的长续航、低成本与高灵活性,提升探测和巡航效率。
基于此,本发明采用无人机作为一种自主移动设备,实现及时有效的灾后环境和危险源探测,为构建高效可靠的灾难应急救援体系提供必要的实现基础。
优点:
(1)煤矿井下无人机结构参数直接影响无人机的飞行平稳性、安全性和制造经济性;并且用于度量上述性能的无人机质量、升力和成本之间存在彼此制约关系。若要获得兼顾以上目标的最佳结构方案,是一个NP-hard问题。本发明以无人机质量、升力和成本作为目标函数,将煤矿井下无人机结构设计建模为一个多目标优化问题,进而采用多目标优化方法进行求解,从而获得一组综合考虑上述三个目标的无人机结构方案集,通过对最优结构方案进行瞬态模态响应仿真验证,获得满足实际工程需求的无人机结构最优设计方案。
(2)通过建立无人机升力代理模型,避免对每个结构方案进行耗时的数值模拟,从而降低结构方案评价的计算成本。本发明采用拉丁超立方抽样方法和数值模拟软件,获得无人机结构方案的初始样本集;采用Kriging方法和径向基函数构成集成学习模型,显著提升了代理模型的整体准确性,提高了模型评价可靠性。
附图说明
图1为本发明的煤矿井下无人机结构的多目标智能优化设计方法流程图。
图2为本发明的煤矿井下无人机结构设计的决策变量物理含义示意图。
图2-1为图2的正视结构图。
图3为本发明的煤矿井下无人机结构Pareto最优解集。
图4为本发明的单一桨叶无人机升力仿真图。
图5为本发明的方案一的无人机机架瞬态模态分析。
图6为本发明的方案二的无人机机架瞬态模态分析。
图7为本发明的方案三的无人机机架瞬态模态分析。
图中,1、电机座直径;2、机架直径;3、框架宽度;4、起落架高度;5、厚度;6、高度;7、支架个数;8、桨叶尺寸。
具体实施方式
无人机多目标智能优化设计方法:建立煤矿井下无人机结构的多目标优化模型,并采用多目标优化方法求解该结构优化设计模型,为煤矿井下巷道狭闭空间内实现巡航的无人机设计方案,其步骤如下:
Step1:建立煤矿井下无人机结构方案的升力代理模型。
井下无人机由机架和动力系统组成,整机具有对称性;在机架和动力系统的共同作用下,无人机的升力由单一桨叶上提供的拉力差决定;无人机机架结构与动力系统包含多个参数,其与无人机升力之间存在高度非线性关系,且无法建立机理模型,采用辨识方法,构建无人机升力的代理模型,用于后续智能优化过程中的个体适应度评价。
Step2:构建兼顾无人机质量、成本和升力的多目标优化设计模型。
为了满足煤矿井下探测和巡航工作需要,同时满足无人机的飞行安全性和制造经济性,无人机结构方案的评价指标包括两个部分;其一,采用无人机质量f1(x)和升力f2(x),评估无人机的飞行平稳性,其中,无人机升力采用代理模型,无人机质量通过机架结构参数和机架材料加以描述;其二,采用无人机成本f3(x)评价所设计的无人机结构经济性,选取无人机机架和动力系统所需材料的总价作为无人机成本;确定决策变量及其约束范围,兼顾无人机质量、升力和成本三个优化目标,建立煤矿井下无人机结构的三目标优化设计模型。
Step3:煤矿井下无人机结构的多目标优化设计方法。
采用多目标智能优化算法,求解煤矿井下无人机结构设计三目标优化模型,得到Pareto 最优解集,实现井下无人机结构参数的优化设计。
Step4:数值模拟验证获得的最优结构方案性能。
为了验证获得的最优结构方案是否能够满足实际工程需要,从Pareto最优解集中选取满足成本预算需求的若干结构方案,采用Fluent数值模拟软件,对无人机升力进行分析验证,并对无人机机架进行模态分析,验证机架稳定性和飞行安全性。
所述的步骤1中,建立煤矿井下无人机结构方案的升力代理模型,采用以下步骤获得:
Step1:选取无人机机架和动力系统关键参数,确定煤矿井下无人机结构优化的决策变量,记为x={x1,x2,x3,x4,x5,x6,x7,x8};其中,x1表示电机座直径,x2表示机架直径,x3表示框架宽度,x4表示起落架高度,x5表示厚度,x6表示高度,x7表示支架个数,x8表示桨叶尺寸。
Step2:根据无人机原材料制造限制,x5,x7,x8采用离散整数值;其余变量采用连续值;确定各个决策变量的取值范围,记为 x5∈{1,1.5,2,2.5},x7∈{2,3,4,5,6,7,8},x8∈{127.0,152.4}。其中,xi ,i=1,2,3,4,6 表示相应的决策变量取值下限,表示相应的决策变量取值上限。
Step3:采用拉丁超立方抽样法,产生无人机结构样本数据,记为D=[xi];其中,xi,i=1,2,…,N为第i个样本,N为样本规模。
Step4:采用Fluent数值模拟软件,对相应结构方案建立三维网格和有限元分析模型,获得每种结构方案的无人机升力,记为U(xi),xi∈D。
Step6:采用由Kriging法和径向基函数构成的集成学习模型,训练样本数据,辨识无人机结构参数和升力之间的非线性映射关系,构建高精度的无人机升力代理模型f2(x);其中,各个子模型的输出通过加权方法加以集成,权值根据其训练误差动态修正;令为第l个子模型的训练误差,为无人机升力估计值,则第k个子模型权重wk记为:
所述的步骤2中,构建兼顾无人机质量、成本和升力的多目标优化设计模型,采用以下步骤获得:
Step1:记ρt为无人机机架密度,ρs为无人机支撑架密度,d为无人机支架宽度,M0为桨叶的单位质量,ct为无人机机架造价,cs为无人机支架造价,kf为桨叶造价,建立无人机质量f1(x)和无人机成本f3(x)的目标函数分别为:
Step2:构建煤矿井下无人机结构的三目标优化设计模型,记为
F(x)={minf1(x),maxf2(x),minf3(x)}。
所述的步骤3中,煤矿井下无人机结构的多目标优化设计方法,采用以下步骤获得:
Step1:设定种群规模为N,最大进化代数为Gen;设置优化算法的初始进化代数t=0,设置初始种群P(0)。
Step2:对个体进行适应度评价,并根据个体的适应度值进行非支配排序。
Step3:依据个体的序值和拥挤距离,选择非占优个体。
Step4:通过交叉和变异操作,生成下一代种群P(t+1)。
Step5:判断是否达到最大迭代次数,若满足,则输出Pareto最优解集;否则,t=t+1,跳转到Step2。
所述的步骤4中,数值模拟验证获得的最优结构方案性能,采用以下步骤获得:
Step1:从Pareto最优解集中,分别选取具有最小成本、最小质量、最大升力和具有折中目标值的4种无人机结构方案。
Step2:结合SolidWorks和Fluent仿真软件,仿真分析无人机结构方案中单一桨叶的拉力,并将其作为机架的反向作用力,求解无人机的瞬态模态响应。
Step3:选取可以满足无人机质量、升力和成本实际需求的最优无人机结构方案。
下面结合具体附图和实例对本发明所提方法的实施方式进行详细说明。
实施例1:煤矿井下无人机结构的多目标智能优化设计方法,其特征是,包括如下步骤:
步骤1)选取无人机机架和动力系统关键参数,确定煤矿井下无人机结构优化的决策变量,确定其变量类型和取值范围。
步骤2)采用拉丁超立方抽样法,结合Fluent数值模拟软件,产生无人机结构及其升力的样本数据集。
步骤3)采用由Kriging法和径向基函数构成的集成学习模型,建立煤矿井下无人机结构方案的升力代理模型。
步骤4)兼顾无人机质量、升力和成本三个优化目标,建立煤矿井下无人机结构的三目标优化设计模型。
步骤5)采用多目标智能优化算法,求解煤矿井下无人机结构设计三目标优化模型,得到Pareto最优解集,实现井下无人机结构参数的优化设计。
步骤6)分别选取具有最小成本、最小质量、最大升力和具有折中目标值的4种无人机结构方案;结合SolidWorks和Fluent仿真软件,仿真分析无人机结构方案中单一桨叶的拉力,并将其作为机架的反向作用力,求解无人机的瞬态模态响应;选取可以满足无人机质量、升力和成本实际需求的最优无人机结构方案。
所述步骤1)中,确定无人机结构参数的具体步骤为:
Step1:决定无人机质量、无人机升力和无人机成本的关键因素是无人机的机架与动力系统,其具体参数包括:电机座直径x1,机架直径x2,框架宽度x3,起落架高度x4,厚度x5,高度x6,支架个数x7,桨叶尺寸x8;
Step2:根据煤矿井下无人机的设计要求,确定上述结构参数的取值范围为:x1∈[24,60], x2∈[380,450],x3∈[0,40],x4∈[10,100],x5∈{1,1.5,2,2.5},x6∈[260,290],x7∈{2,3,4,5,6,7,8}, x8∈{127.0,152.4}。
所述步骤2)中,产生无人机结构及其升力的样本数据集的具体步骤为:
Step1:根据结构参数的取值范围,产生20组无人机结构样本数据;
Step2:通过Fluent对无人机的整体结构进行非结构网格划分,获得无人机结构参数的单一桨叶升力,构成样本数据库。
表1无人机结构参数样本数据
所述步骤3)中,建立煤矿井下无人机结构方案的升力代理模型具体步骤为:
Step1:高斯径向基函数中的参数c选取为1、0.1、0.05、0.01,Multiquadric径向基函数中的参数c选取为[0.1,1,2,4,4.5,5],采用CV验证分别获取其最佳值;构建Kriging模型,其关联函数中的系数θ取值如下表2:
表2 Kriging模型中θ取值
Step2:采用判定系数R2和均方根误差RMSE评价模型,分析各个子模型精度如表3。
表3模型精度评价
所述步骤4)中,建立煤矿井下无人机结构的三目标优化设计模型的具体步骤为:
根据无人机其内部参数表,建立无人机质量f1(x)和无人机成本f3(x)的目标函数分别为:
表4无人机机架和动力系统的内部参数
由此,建立煤矿井下无人机结构的三目标优化设计模型为:
F(x)=(minf1(x),maxf2(x),minf3(x))T。
所述步骤6)中,4种无人机结构方案的模拟数值分析步骤如下:
Step1:从图3所示的Pareto最优结构集中选取4组设计方案,其结构参数和目标值如表5所示:
表5无人机最佳结构设计方案
Step2:导入创建的实体模型或直接利用节点和单元生成有限元模型,定义模型的单元类型与材料属性;
Step3:设置物理环境、边界条件与接触条件,对结构设计方案进行求解验证,验证无人机的稳定性。无人机有限元分析仿真验证如图4-图7所示。
Claims (5)
1.一种煤矿井下无人机结构的多目标智能优化设计方法,其特征在于:无人机多目标智能优化设计方法:建立煤矿井下无人机结构的多目标优化模型,并采用多目标优化方法求解该结构优化设计模型,为煤矿井下巷道狭闭空间内实现巡航的无人机设计方案,其步骤如下:
Step1:建立煤矿井下无人机结构方案的升力代理模型;
井下无人机由机架和动力系统组成,整机具有对称性;在机架和动力系统的共同作用下,无人机的升力由单一桨叶上提供的拉力差决定;无人机机架结构与动力系统包含多个参数,其与无人机升力之间存在高度非线性关系,且无法建立机理模型,采用辨识方法,构建无人机升力的代理模型,用于后续智能优化过程中的个体适应度评价;
Step2:构建兼顾无人机质量、成本和升力的多目标优化设计模型;
为了满足煤矿井下探测和巡航工作需要,同时满足无人机的飞行安全性和制造经济性,无人机结构方案的评价指标包括两个部分;其一,采用无人机质量f1(x)和升力f2(x),评估无人机的飞行平稳性,其中,无人机升力采用代理模型,无人机质量通过机架结构参数和机架材料加以描述;其二,采用无人机成本f3(x)评价所设计的无人机结构经济性,选取无人机机架和动力系统所需材料的总价作为无人机成本;确定决策变量及其约束范围,兼顾无人机质量、升力和成本三个优化目标,建立煤矿井下无人机结构的三目标优化设计模型;
Step3:煤矿井下无人机结构的多目标优化设计方法;
采用多目标智能优化算法,求解煤矿井下无人机结构设计三目标优化模型,得到Pareto最优解集,实现井下无人机结构参数的优化设计;
Step4:数值模拟验证获得的最优结构方案性能;
为了验证获得的最优结构方案是否能够满足实际工程需要,从Pareto最优解集中选取满足成本预算需求的若干结构方案,采用Fluent数值模拟软件,对无人机升力进行分析验证,并对无人机机架进行模态分析,验证机架稳定性和飞行安全性。
2.根据权利要求1所述的煤矿井下无人机结构的多目标智能优化设计方法,其特征在于:步骤1中,建立煤矿井下无人机结构方案的升力代理模型,采用以下步骤获得:
Step1:选取无人机机架和动力系统关键参数,确定煤矿井下无人机结构优化的决策变量,记为x={x1,x2,x3,x4,x5,x6,x7,x8};其中,x1表示电机座直径,x2表示机架直径,x3表示框架宽度,x4表示起落架高度,x5表示厚度,x6表示高度,x7表示支架个数,x8表示桨叶尺寸;
Step2:根据无人机原材料制造限制,x5,x7,x8采用离散整数值;其余变量采用连续值;确定各个决策变量的取值范围,记为x5∈{1,1.5,2,2.5},x7∈{2,3,4,5,6,7,8},x8∈{127.0,152.4};
Step3:采用拉丁超立方抽样法,产生无人机结构样本数据,记为D=[xi];其中,xi,i=1,2,…,N为第i个样本,N为样本规模;
Step4:采用Fluent数值模拟软件,对相应结构方案建立三维网格和有限元分析模型,获得每种结构方案的无人机升力,记为U(xi),xi∈D;
Step6:采用由Kriging法和径向基函数构成的集成学习模型,训练样本数据,辨识无人机结构参数和升力之间的非线性映射关系,构建高精度的无人机升力代理模型f2(x);其中,各个子模型的输出通过加权方法加以集成,权值根据其训练误差动态修正;令为第l个子模型的训练误差,为无人机升力估计值,则第k个子模型权重wk记为:
4.根据权利要求1所述的煤矿井下无人机结构的多目标智能优化设计方法,其特征在于:步骤3中,煤矿井下无人机结构的多目标优化设计方法,采用以下步骤获得:
Step1:设定种群规模为N,最大进化代数为Gen;设置优化算法的初始进化代数t=0,设置初始种群P(0);
Step2:对个体进行适应度评价,并根据个体的适应度值进行非支配排序;
Step3:依据个体的序值和拥挤距离,选择非占优个体;
Step4:通过交叉和变异操作,生成下一代种群P(t+1);
Step5:判断是否达到最大迭代次数,若满足,则输出Pareto最优解集;否则,t=t+1,跳转到Step2。
5.根据权利要求1所述的煤矿井下无人机结构的多目标智能优化设计方法,其特征在于:步骤4中,数值模拟验证获得的最优结构方案性能,采用以下步骤获得:
Step1:从Pareto最优解集中,分别选取具有最小成本、最小质量、最大升力和具有折中目标值的4种无人机结构方案;
Step2:结合SolidWorks和Fluent仿真软件,仿真分析无人机结构方案中单一桨叶的拉力,并将其作为机架的反向作用力,求解无人机的瞬态模态响应;
Step3:选取可以满足实际需求的无人机结构方案。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010458424.8A CN111666630B (zh) | 2020-05-26 | 2020-05-26 | 一种煤矿井下无人机结构的多目标智能优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010458424.8A CN111666630B (zh) | 2020-05-26 | 2020-05-26 | 一种煤矿井下无人机结构的多目标智能优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111666630A CN111666630A (zh) | 2020-09-15 |
CN111666630B true CN111666630B (zh) | 2023-05-09 |
Family
ID=72384762
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010458424.8A Active CN111666630B (zh) | 2020-05-26 | 2020-05-26 | 一种煤矿井下无人机结构的多目标智能优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111666630B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113516368A (zh) * | 2021-06-11 | 2021-10-19 | 北京大学 | 城市和社区不确定性风险的预测方法、装置、设备及介质 |
CN117313237B (zh) * | 2023-09-23 | 2024-04-16 | 哈尔滨工业大学 | 一种基于机器学习的特种无人机组构方案优选方法 |
CN118070678B (zh) * | 2024-04-15 | 2024-06-21 | 海澜智云科技有限公司 | 一种基于流程模拟和机器学习的煤气化装置智能优化方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106021808A (zh) * | 2016-06-07 | 2016-10-12 | 西北工业大学 | 考虑三维效应的小展弦比机翼翼型设计方法 |
CN106529093A (zh) * | 2016-12-15 | 2017-03-22 | 北京航空航天大学 | 一种针对大展弦比机翼的气动/结构/静气弹耦合优化方法 |
-
2020
- 2020-05-26 CN CN202010458424.8A patent/CN111666630B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106021808A (zh) * | 2016-06-07 | 2016-10-12 | 西北工业大学 | 考虑三维效应的小展弦比机翼翼型设计方法 |
CN106529093A (zh) * | 2016-12-15 | 2017-03-22 | 北京航空航天大学 | 一种针对大展弦比机翼的气动/结构/静气弹耦合优化方法 |
Non-Patent Citations (5)
Title |
---|
Low-Cost Approaches to UAV Design Using Advanced Manufacturing Techniques;Eric P. Flynn;3rd IEEE Integrated STEM Education Conference;全文 * |
基于Kriging模型机翼平面外形气动优化设计;于方圆;高永;王允良;韩维;;海军航空工程学院学报(第05期);全文 * |
基于Modelica的四旋翼飞行器多领域统一建模及多目标优化设计;马小辉;硕士电子期刊全文数据库;全文 * |
基于代理模型的翼型多目标集成优化设计;龙腾;刘莉;孟令涛;朱华光;;系统仿真学报(第07期);全文 * |
基于计算试验设计与代理模型的飞行器 近似优化策略探讨;龙腾 等;机械工程学报;第52卷(第14期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111666630A (zh) | 2020-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111666630B (zh) | 一种煤矿井下无人机结构的多目标智能优化设计方法 | |
CN110276479B (zh) | 飞机质量变化的巡航阶段燃油消耗预测方法 | |
CN109190283A (zh) | 一种考虑高湍流自由来流效应的风力机翼型气动优化方法 | |
Benini et al. | High-lift multi-element airfoil shape and setting optimization using multi-objective evolutionary algorithms | |
CN106126860B (zh) | 一种考虑加工误差的高超声速机翼鲁棒优化设计方法 | |
CN103488869A (zh) | 一种最小二乘支持向量机的风力发电短期负荷预测方法 | |
CN107480335A (zh) | 一种高超声速飞行器模型迭代设计方法 | |
CN109307845A (zh) | 一种用于对梯次利用电池储能系统进行健康状态评估的方法及系统 | |
CN113344408B (zh) | 民航交通管制运行多尺度态势感知流程的处理方法 | |
CN111045080A (zh) | 一种基于pso-bp模型和地震属性参数的煤层含气量预测方法 | |
CN105005878A (zh) | 一种坚强智能电网的综合评估方法 | |
CN114357606A (zh) | 一种基于模型驱动的飞机机电系统效能权衡决策方法 | |
CN102968553A (zh) | 一种飞行器着陆风险评测方法 | |
CN101898644A (zh) | 用来设计翼型的系统与方法 | |
CN105138766B (zh) | 基于模糊聚类的高超声速气动热降阶模型的加点方法 | |
Fitas | Optimal design of composite structures using the particle swarm method and hybridizations | |
CN117313235A (zh) | 一种基于机器学习的特种无人机优化翼型桨叶设计方法 | |
CN104462787A (zh) | 基于rst/iahp的航空维修人因可靠性评估方法 | |
CN117057190A (zh) | 一种轨道车辆制动盘温升预测方法 | |
Yadlin et al. | Icing simulation research supporting the ice-accretion testing of large-scale swept-wing models | |
Kirz | Surrogate-Based Low-Boom Low-Drag Nose Design for the JAXA S4 Supersonic Airliner | |
Liang et al. | Plithogenic multi‐criteria decision making approach on airspace planning scheme evaluation based on ATC‐flight real‐time simulation | |
Trevisi | Configuration optimisation of kite-based wind turbines | |
OLADEJO et al. | Aerodynamic Performance Analysis of an Optimized Aerofoil for Unmanned Aerial Vehicles | |
CN117313237B (zh) | 一种基于机器学习的特种无人机组构方案优选方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |