CN109815439B - 一种减速器设计协同优化方法 - Google Patents

一种减速器设计协同优化方法 Download PDF

Info

Publication number
CN109815439B
CN109815439B CN201811622395.3A CN201811622395A CN109815439B CN 109815439 B CN109815439 B CN 109815439B CN 201811622395 A CN201811622395 A CN 201811622395A CN 109815439 B CN109815439 B CN 109815439B
Authority
CN
China
Prior art keywords
level
optimization
objective function
design
subject
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
Application number
CN201811622395.3A
Other languages
English (en)
Other versions
CN109815439A (zh
Inventor
郑松
徐柯
葛铭
郑小青
魏江
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Dianzi University
Original Assignee
Hangzhou Dianzi University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN201811622395.3A priority Critical patent/CN109815439B/zh
Publication of CN109815439A publication Critical patent/CN109815439A/zh
Application granted granted Critical
Publication of CN109815439B publication Critical patent/CN109815439B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明提供了一种减速器设计协同优化方法,本发明将减速器的设计问题分解为一个系统级和两个学科级问题,系统级为减速器体积最小的优化问题,齿轮和轴的约束分别形成两个子学科优化问题。按照如下步骤进行优化:(1)系统级传递设定的初始参数给子学科,学科级目标函数中加入了系统级目标函数来减小初始点的影响,实现学科级的优化。(2)将各个子学科的优化结果传递给系统级的约束条件,在系统级的约束条件中加入动态松弛因子进行优化。(3)经过系统级和各个学科级之间的多次迭代,满足收敛条件,得到优化的结果。本发明有效的解决传统协同优化算法可行域不存在、优化收敛速度慢、初始点选取不当的缺陷。

Description

一种减速器设计协同优化方法
技术领域
本发明属于优化算法应用技术领域,具体是一种减速器设计协同优化方法。
背景技术
减速器设计优化其目标是在满足转轴和齿轮大量约束的同时,使减速器的体积最小(质量最轻)。协同优化(Collaborative Optimization,CO)方法是斯坦福大学KROO等提出的一种解决复杂优化问题的有效算法,其基本原理是将复杂优化设计问题分解为一个系统级优化问题和几个并行的学科级优化问题。CO方法的各学科寻优计算都有很好的独立性,而不同学科所得优化结果的不一致性则由系统级优化来协调,经过多次迭代获得最优的结果。但CO仍存在不完善之处,其系统级采用的一致性等式约束会导致优化结果陷人局部最优解,出现收敛速度慢等问题。另外,CO的初始点选取不当也会造成优化结果陷入局部极值点的问题。
发明内容
本发明的目的在于提供一种减速器设计协同优化方法,有效解决协同优化算法易陷入局部最优、收敛速度慢、初始点选取不当的缺陷对协同优化算法进行改进,提高减速器设计的优化精度和速度,以解决上述背景技术中提出的问题。
本发明为了实现上述目的,采用以下方案:
本发明提供一种减速器设计协同优化方法,将减速器的设计问题分解为一个系统级和两个学科级问题,系统级为减速器体积最小的优化问题,齿轮和轴的约束分别形成两个子学科优化问题,其特征在于,包括以下步骤:
步骤1:初始化减速器的设计参数,该参数包括设计变量的初始点Z0、学科级目标函数的控制系数D、系统级目标函数的权重系数d和控制系数K,将初始点Z0作为系统级的优化结果;
步骤2:系统级的优化结果传递给学科级,在学科级目标函数中加入系统级目标函数,满足自身约束条件下进行学科级优化求解;
步骤3:将学科级优化结果传递给系统级进行优化,在系统级的约束条件中加入松弛因子ε,增大寻优范围;
步骤4:更改系统级目标函数权重系数d和控制系数K,对减速器协同优化算法的优化结果判断是否收敛,若是,则输出全局最优解的位置,作为优化问题的解;若否,则返回步骤2继续进行迭代;
其中系统级目标函数权重系数d和控制系数K的更改如下:
Figure BDA0001927198860000021
式中,
Figure BDA0001927198860000022
表示第n次迭代系统级的优化结果,
Figure BDA0001927198860000023
表示第n+1次迭代系统级的优化结果,
Figure BDA0001927198860000024
表示系统级前后两次优化设计点的距离;
Figure BDA0001927198860000025
式中,d表示系统级目标函数权重系数。
作为本发明进一步的方案:所述步骤2中学科级优化中加入系统级目标函数如下:
Figure BDA0001927198860000026
式中,Ji(xi)表示第i个学科级的目标函数,xi表示第i个学科级的设计变量点,D表示学科级目标函数的控制系数,d和K分别表示系统级目标函数的权重系数和控制系数,s表示该学科级设计变量的个数,xij表示第i个学科级的第j个设计变量,
Figure BDA00019271988600000212
表示系统级分配给学科级的第j个设计变量的期望值,F(xi)表示系统级的目标函数。
作为优选,所述步骤3中在系统级的约束条件中加入松弛因子ε具体如下:
Figure BDA0001927198860000027
式中,
Figure BDA0001927198860000028
表示第i个系统级的约束条件,si表示系统级设计变量的个数,zj表示第j个系统级变量,
Figure BDA0001927198860000029
表示第i个学科级优化结果传递给系统级的第j个设计变量,ε表示松弛因子由学科间不一致信息来确定;
ε=(λ×Δ)2
式中,λ表示常数取值范围0.5≤λ≤1,Δ代表了学科级之间的不一致性;
Figure BDA00019271988600000210
式中,||*||表示2范数,
Figure BDA00019271988600000211
表示任意两个学科级之间的2范数,max表示取最大值。
与现有技术相比,本发明的有益效果是:
本发明建立了一种减速器设计协同优化方法,与现有协同优化算法相比,本发明方法针对初始点选取、可行域不存在、收敛速度慢等问题,用协同优化系统级的约束中加入松弛因子和学科级目标函数中加入系统级目标函数的改进方法,将复杂优化设计问题分解为一个系统级优化问题和几个并行的学科级优化问题,减小初始点的影响和提高计算精度,同时又能减少迭代次数提高减速器设计的优化精度和速度。
附图说明
图1为本发明的流程示意图。
图2为本发明实施例涉及的改进协同优化迭代过程示意图。
图3为本发明实施例涉及的标准协同优化迭代过程示意图。
具体实施方式
以下结合附图对本发明涉及的一种减速器设计协同优化方法具体实施方案进行详细地说明。
将减速器的设计问题分解为一个系统级和两个学科级问题,系统级为减速器体积最小的优化问题,齿轮和轴的约束分别形成两个子学科优化问题。如图1所示,本实施例所提供的一种减速器设计协同优化方法包括以下步骤:
步骤1:初始化减速器的设计参数,该参数包括设计变量的初始点Z0、学科级目标函数的控制系数D、系统级目标函数的权重系数d和控制系数K,将初始点Z0作为系统级的优化结果;
步骤2:系统级的优化结果传递给学科级,在学科级目标函数中加入系统级目标函数,满足自身约束条件下进行学科级优化;
步骤3:将学科级优化结果传递给系统级进行优化,在系统级的约束条件中加入松弛因子ε,增大寻优范围;
步骤4:更改系统级目标函数权重系数d和控制系数K,对减速器协同优化算法的优化结果判断是否收敛,若是,则输出全局最优解的位置,作为优化问题的解;若否,则返回步骤2继续进行迭代。
所述步骤1中初始化减速器的参数如下:
Z0=[z01,z02,z03,z04,z05,z06,z07]
式中,Z0表示系统级的初始点,z01,z02,…,z07分别表示系统级设计变量的初始值;
在第一次迭代过程中为减小初始点的影响,学科级目标函数的控制系数D设为0,系统级目标函数的权重系数d和控制系数K都设为1。
所述步骤2中学科级优化中加入系统级目标函数具体如下:
Figure BDA0001927198860000031
式中,Ji(xi)表示第i个学科级的目标函数,xi表示第i个学科级的设计变量点,D表示学科级目标函数的控制系数,d和K分别表示系统级目标函数的权重系数和控制系数,s表示该学科级设计变量的个数,xij表示第i个学科级的第j个设计变量,
Figure BDA0001927198860000032
表示系统级分配给学科级的第j个设计变量的期望值,F(xi)表示系统级的目标函数;
在学科级优化过程中学科级目标函数和系统级目标函数共同作用下向全局最优设计点逼近,增加寻优范围,从而减小初始点选择的依赖,利于寻找全局最优设计点。
所述步骤3中在系统级的约束条件中加入松弛因子ε具体如下:
Figure BDA0001927198860000033
式中,
Figure BDA0001927198860000035
表示第i个系统级的约束条件,si表示系统级设计变量的个数,zj表示第j个系统级变量,
Figure BDA0001927198860000034
表示第i个学科级优化结果传递给系统级的第j个设计变量,ε表示松弛因子由学科间不一致信息来确定;
ε=(λ×Δ)2
式中,λ表示常数取值范围0.5≤λ≤1,Δ代表了学科级之间的不一致性;
Figure BDA0001927198860000041
式中,||*||表示2范数,
Figure BDA0001927198860000042
表示任意两个学科级之间的2范数,max表示取最大值;
在系统级的约束条件中加入基于学科间不一致信息的动态松弛因子,在优化的初期,Δ学科不一致性较大,相应的ε松弛因子较大,决定了系统级可行域范围较大,从而加快收敛,确保了系统级优化可行域的存在。
所述步骤4中更改系统级目标函数权重系数d和控制系数K具体如下:
Figure BDA0001927198860000043
式中,
Figure BDA0001927198860000044
表示第n次迭代系统级的优化结果,
Figure BDA0001927198860000045
表示第n+1次迭代系统级的优化结果,
Figure BDA0001927198860000046
表示系统级前后两次优化设计点的距离,随着迭代的进行d会逐渐变小从而改变系统级目标函数在学科级优化设计时的权重;
Figure BDA0001927198860000047
式中,d表示系统级目标函数的权重系数,当d<0.05时,K为0,去除了学科级中系统级目标函数的作用,避免学科级在寻找设计点时的震荡,使之尽快趋于收敛。
根据如上内容,本实施所提供的一种减速器设计协同优化方法,有效的解决了初始点选取、可行域不存在、收敛速度慢等问题,提高了计算精度,加快了计算速度。
为验证上述协同优化算法的正确性和有效性,本发明对减速器设计在不同初始点条件下进行了仿真验证,与使用动态松弛因子的标准协同优化在优化结果精度、迭代次数方面进行了对比。本实验的开发环境为MATLABR2017a。不同初始点的数值如表1所示。标准协同优化和改进协同优化的优化结果对比如表2所示。
表1不同初始点的数值
Figure BDA0001927198860000048
表2标准协同优化和改进协同优化的优化结果对比
Figure BDA0001927198860000049
Figure BDA0001927198860000051
由表1和表2可知,该改进协同优化算法的方法与标准动态松弛协同优化相比。标准协同优化算法易受初始点的影响,不同的初始点迭代计算所得优化结果均不同。而改进的协同优化算法能够避免初始点选取的问题,不同的初始点迭代计算所得优化结果都极为接近,并且优化结果均优于标准协同优化,加快了收敛计算的速度。
在不同初始点条件下,两种优化算法的迭代过程如图2和图3所示。标准协同优化在不同初始点情况下,每次迭代的结果都不同,原因在于学科级的目标函数以系统级优化结果为参考,在满足自身约束的情况下尽可能的向其靠近。如果初始点选取不当,就会使学科级的最优解向初始点靠近,经过多次迭代之后得到的最优解受到影响从而陷入局部最优解。本发明的协同优化改变传统协同优化算法单方面考虑系统级优化结果的方法。在迭代初期,将初始点拉向全局最优点附近,减少初始点的影响,提高减速器设计的优化精度和速度。

Claims (3)

1.一种减速器设计协同优化方法,将减速器的设计问题分解为一个系统级和两个学科级问题,系统级为减速器体积最小的优化问题,齿轮和轴的约束分别形成两个子学科优化问题,其特征在于,包括如下步骤:
步骤1:初始化减速器的设计参数,该参数包括设计变量的初始点Z0、学科级目标函数的控制系数D、系统级目标函数的权重系数d和控制系数K,将初始点Z0作为系统级的优化结果;
步骤2:系统级的优化结果传递给学科级,在学科级目标函数中加入系统级目标函数,满足自身约束条件下进行学科级优化求解;
步骤3:将学科级优化结果传递给系统级进行优化,在系统级的约束条件中加入松弛因子ε,增大寻优范围;
步骤4:更改系统级目标函数权重系数d和控制系数K,对减速器协同优化算法的优化结果判断是否收敛,若是,则输出全局最优解的位置,作为优化问题的解;若否,则返回步骤2,继续进行迭代;
其中系统级目标函数权重系数d和控制系数K的更改如下:
Figure FDA0001927198850000011
式中,
Figure FDA00019271988500000110
表示第n次迭代系统级的优化结果,
Figure FDA0001927198850000012
表示第n+1次迭代系统级的优化结果,
Figure FDA0001927198850000013
表示系统级前后两次优化设计点的距离;
Figure FDA0001927198850000014
式中,d表示系统级目标函数权重系数。
2.根据权利要求1所述的一种减速器设计协同优化方法,其特征在于:所述步骤2中学科级优化中加入系统级目标函数如下:
Figure FDA0001927198850000015
式中,Ji(xi)表示第i个学科级的目标函数,xi表示第i个学科级的设计变量点,D表示学科级目标函数的控制系数,d和K分别表示系统级目标函数的权重系数和控制系数,s表示该学科级设计变量的个数,xij表示第i个学科级的第j个设计变量,
Figure FDA0001927198850000016
表示系统级分配给学科级的第j个设计变量的期望值,F(xi)表示系统级的目标函数。
3.根据权利要求1所述的一种减速器设计协同优化方法,其特征在于:所述步骤3中在系统级的约束条件中加入松弛因子ε具体如下:
Figure FDA0001927198850000017
式中,
Figure FDA0001927198850000018
表示第i个系统级的约束条件,si表示系统级设计变量的个数,zj表示第j个系统级变量,
Figure FDA0001927198850000019
表示第i个学科级优化结果传递给系统级的第j个设计变量,ε表示松弛因子由学科间不一致信息来确定;
ε=(λ×Δ)2
式中,λ表示常数取值范围0.5≤λ≤1,Δ代表了学科级之间的不一致性;
Figure FDA0001927198850000021
式中,||*||表示2范数,
Figure FDA0001927198850000022
表示任意两个学科级之间的2范数,max表示取最大值。
CN201811622395.3A 2018-12-28 2018-12-28 一种减速器设计协同优化方法 Active CN109815439B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811622395.3A CN109815439B (zh) 2018-12-28 2018-12-28 一种减速器设计协同优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811622395.3A CN109815439B (zh) 2018-12-28 2018-12-28 一种减速器设计协同优化方法

Publications (2)

Publication Number Publication Date
CN109815439A CN109815439A (zh) 2019-05-28
CN109815439B true CN109815439B (zh) 2022-11-25

Family

ID=66601803

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811622395.3A Active CN109815439B (zh) 2018-12-28 2018-12-28 一种减速器设计协同优化方法

Country Status (1)

Country Link
CN (1) CN109815439B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111461387B (zh) * 2019-12-30 2022-08-23 杭州电子科技大学 一种环保岛系统设计的径向基协同优化方法
CN116306383B (zh) * 2023-05-22 2024-01-30 华东交通大学 展向波纹杆件协同优化方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015180196A1 (zh) * 2014-05-27 2015-12-03 中国空间技术研究院 一种基于多学科优化的卫星总体方案确定系统及实现方法
CN106294972A (zh) * 2016-08-05 2017-01-04 四川理工学院 一种车桥多学科可靠性设计优化方法
CN106650999A (zh) * 2016-10-25 2017-05-10 杭州电子科技大学 一种啤酒生产调度优化方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015180196A1 (zh) * 2014-05-27 2015-12-03 中国空间技术研究院 一种基于多学科优化的卫星总体方案确定系统及实现方法
CN106294972A (zh) * 2016-08-05 2017-01-04 四川理工学院 一种车桥多学科可靠性设计优化方法
CN106650999A (zh) * 2016-10-25 2017-05-10 杭州电子科技大学 一种啤酒生产调度优化方法

Also Published As

Publication number Publication date
CN109815439A (zh) 2019-05-28

Similar Documents

Publication Publication Date Title
CN109815439B (zh) 一种减速器设计协同优化方法
CN110826791A (zh) 一种基于长短时记忆神经网络的混合风功率预测方法
CN108319147B (zh) 一类具有短时延和数据丢包的网络化线性参数变化系统的h∞容错控制方法
CN107943919A (zh) 一种面向会话式实体搜索的查询扩展方法
CN107121929A (zh) 基于线性协方差模型预测控制的鲁棒再入制导方法
CN104217041A (zh) 一种多约束在线高斯伪谱末制导方法
CN106294418A (zh) 检索方法和检索系统
CN109389222A (zh) 一种快速的自适应神经网络优化方法
CN104616072A (zh) 一种基于区间优化的提高谷氨酸发酵产物浓度的方法
Liu et al. Parameter calibration in wake effect simulation model with stochastic gradient descent and stratified sampling
CN107039972A (zh) 一种基于afsa‑bfo算法的电力系统pss参数整定方法
CN109325288A (zh) 一种基于不确定性优化的固体运载器总体参数确定方法及系统
CN107220233A (zh) 一种基于高斯混合模型的用户知识需求模型构建方法
CN110516231A (zh) 基于注意力机制的膨胀卷积实体名识别方法
CN103809446B (zh) 飞行器多回路模型簇颤振抑制复合频率鲁棒控制器设计方法
CN110826182B (zh) 基于顶点法和序贯优化策略的飞行器结构气动弹性设计法
CN105976057B (zh) 一种基于改进灰色理论的中长期负荷预测方法
CN104991447B (zh) 一种基于小波神经网络的火电机组变负荷速率预测方法
CN115358167B (zh) 一种考虑发动机参数的飞发一体气动伴随优化设计方法
CN105095960A (zh) 一种梯度粒子群混合优化方法
CN116975192A (zh) 一种基于树形结构的大语言模型运算链装置及方法
CN116070768A (zh) 一种基于数据重构和TCN-BiLSTM的短期风电功率预测方法
Dong et al. Application of Adam-BP neural network in leveling fitting
CN112560326A (zh) 压力场的确定方法及装置
CN103809453B (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