CN111339611A - 一种非线性动力吸振器的参数优化方法 - Google Patents
一种非线性动力吸振器的参数优化方法 Download PDFInfo
- Publication number
- CN111339611A CN111339611A CN202010091796.1A CN202010091796A CN111339611A CN 111339611 A CN111339611 A CN 111339611A CN 202010091796 A CN202010091796 A CN 202010091796A CN 111339611 A CN111339611 A CN 111339611A
- Authority
- CN
- China
- Prior art keywords
- particle
- vibration absorber
- dynamic vibration
- iteration
- nonlinear dynamic
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种非线性动力吸振器的参数优化方法,采用动态变惯性权重的粒子群优化算法进行非线性动力吸振器的参数优化,惯性权重与迭代次数呈单调递减的二次函数关系。开始迭代时,惯性权重较大,具有较强的全局搜索能力;迭代将要结束时,惯性权重小,具有较强的局部搜索能力。惯性权重随迭代次数的变化关系决定了改进的粒子群优化算法搜索能力快,求解精度高的特点。此外,粒子群优化算法本身算法规则简单,目的性强,易于实现。
Description
技术领域
本发明属于机械振动领域,具体涉及一种非线性动力吸振器的参数优化方法。
背景技术
采用动力吸振器吸振是一种通过设计吸振器结构使其具有适当的阻尼,在此阻尼下动力吸振器频率略微小于结构的共振频率,来减小原结构共振频率下振动响应的控制方法。当动力吸振器的响应与激励不具有一次非线性关系时,该动力吸振器就称为“非线性动力吸振器”。为了获取适当的阻尼,抑制原结构共振频率下的振动响应,实现最好的吸振效果,需要对非线性动力吸振器的参数进行优化。
参数优化是为达到设计目标采取的一种方法,通过将设计目标参数化,采用优化方法,不断的调整设计变量,使得设计结果不断接近目标值。常用数值优化方法和智能优化算法进行参数优化,但现有的优化算法不能兼顾全局寻优和收敛速度,收敛快的易于陷入局部最优解,能够找到全局最优解的往往收敛速度慢,迭代次数较多。
发明内容
发明目的:针对上述现有技术,提出一种非线性动力吸振器的参数优化方法,提高参数优化过程的全局搜索能力,使参数优化过程能够以最快的速度找到全局最优解。
技术方案:一种非线性动力吸振器的参数优化方法,包括如下步骤:
1)基于非线性动力吸振器的结构,建立非线性动力吸振器的动力学模型,量化吸振效果;
2)将吸振效果的量化函数作为目标函数,采用变惯性权重粒子群优化算法进行全局寻优。
进一步的,步骤2)中采用变惯性权重粒子群优化算法进行全局寻优包括如下步骤:
21)初始化粒子群算法因子、粒子群位置和速度;
22)基于惯性权重的计算表达式和迭代次数,计算惯性权重;
23)基于粒子的位置,计算每个粒子的适应度函数值;
24)对于每个粒子,将粒子经历过的最优位置作为粒子个体最优点;
25)对于整个群体,将所有粒子个体所经历过位置中的最优位置作为种群最优点;
26)基于速度更新公式和位置更新公式更新每个粒子的速度和位置;
27)判断迭代次数是否达到迭代次数阈值,如未达到,重复步骤22)到步骤26),否则算法结束。
进一步的,所述步骤21)中,粒子群算法因子包括种群大小n、加速度因子c1、加速度因子c2以及迭代次数阈值H;粒子群位置和速度都是一个3维向量,位置为k=(k1,k2,k3),速度为v=(v1,v2,v3),一个位置代表一个潜在解。
进一步的,所述步骤22)中,惯性权重的计算表达式是关于迭代次数i的一个二次函数,具体为:
wi为惯性权重。
进一步的,所述步骤24)中,最优位置为适应度函数最小时对应的粒子位置。
进一步的,所述步骤26)中,速度更新公式和位置更新公式分别为:
vi+1=wivi+c1·r1·(pbesti-ki)+c2·r2·(gbesti-ki)
ki+1=ki+vi
其中,vi表示第i次迭代粒子的速度,ki表示第i次迭代粒子的位置,r1、r2是一个[0,1]之间的随机数,pbesti是第i次迭代粒子个体最优点,gbesti是第i次迭代种群最优点。
有益效果:采用动态变惯性权重的粒子群优化算法进行非线性动力吸振器的参数优化,惯性权重与迭代次数呈单调递减的二次函数关系。开始迭代时,惯性权重较大,具有较强的全局搜索能力;迭代将要结束时,惯性权重小,具有较强的局部搜索能力。惯性权重随迭代次数的变化关系决定了改进的粒子群优化算法搜索能力快,求解精度高的特点。此外,粒子群优化算法本身算法规则简单,目的性强,易于实现。
附图说明
图1为非线性动力吸振器剖面图,1是导杆、2是橡胶环、3是帘布层、4是簧片、5是支撑弹簧、6是导向筒。
图2为非线性动力吸振器等效力学模型图,ma是导杆1等效的上质量块质量,ms是导向筒6等效的下质量块质量,k1是18个周布簧片4等效成的弹簧刚度,k2是支撑弹簧5的刚度,k3是橡胶环2等效的弹簧元件刚度,c是橡胶环2等效的阻尼元件阻尼系数,F是施加在导向筒6上的激励,xa是导杆1的位移,xs是导向筒6的位移。
图3是变惯性权重的粒子群优化算法流程图。
具体实施方式
下面结合附图对本发明做更进一步的解释。
一种非线性动力吸振器的参数优化方法,包括如下步骤:
1)基于非线性动力吸振器的结构,建立非线性动力吸振器的动力学模型,量化吸振效果;
2)将吸振效果的量化函数作为目标函数,采用变惯性权重粒子群优化算法进行全局寻优。
其中,步骤1)包括如下具体步骤:
图1给出了所述的一种非线性动力吸振器的基本结构,导向筒6位于吸振器最下方,支撑弹簧5下端连接在导向筒6上,上端连接导杆1,导杆1和导向筒6之间还有18个均布的簧片4,簧片4被橡胶环2包裹,最外层的帘布层3裹住橡胶环2。
按照如图1所示的非线性动力吸振器的基本结构,将其等效为集中参数的动力学模型,导杆1等效成上质量块,导向筒6等效成下质量块,18个周布的簧片4等效成一个弹簧,橡胶环2等效成一个弹簧和阻尼元件,如图2所述。
根据动力学模型建立刚体平面运动微分方程:
利用拉普拉斯变换,在0初值条件下求解上面的微分方程,得到以激励为输入、以导杆1位移为输出的传递函数:
在已知激励F的情况下,可以根据拉普拉斯反变换求出导杆1的位移xa:
xa=L-1[G(s)·F(s)]
其中,L-1表示拉普拉斯反变换。
进一步地,对导杆1的位移xa求二阶导数可以得到导杆1的加速度;
基于加速度值计算导杆1一个周期T内的加速度均方根值:
这是吸振器的吸振效果评价函数,导杆1的加速度均方根值越小,说明吸振效果越好。
k3、c是橡胶环的等效元件参数,存在一定的函数关系,最终得到参数优化的目标函数为:
步骤2)中,将上述函数作为粒子群优化算法的目标函数,优化目标是找到目标函数的最小值以及函数取到最小值时对应的k1、k2、k3值。
粒子群优化算法为:首先在解空间内初始化一群粒子,每个粒子代表一个可行解;用粒子的速度和位置表示其特征,引入根据粒子位置计算的适应度函数,以此评价粒子的优劣;通过跟踪粒子个体最优点和种群最优点来更新个体位置;每更新一次粒子位置,就重新计算适应度值,以此更新粒子个体最优点和种群最优点。值得注意的是,采用粒子群优化算法进行参数优化,粒子群优化算法的惯性权重是随着优化迭代次数的增加而减小的。
如图3所示,步骤2)包括如下步骤:
21)初始化粒子群算法因子、粒子群位置和速度。其中,粒子群算法因子包括种群大小n、加速度因子c1、加速度因子c2以及迭代次数阈值H。粒子群位置和速度都是一个3维向量,位置为k=(k1,k2,k3),速度为v=(v1,v2,v3),一个位置代表一个潜在解,速度向量的每个分量表示每一次迭代中对应刚度的变化值。
22)基于惯性权重的计算表达式和迭代次数,计算惯性权重wi。惯性权重的计算表达式是关于迭代次数i的一个二次函数,具体为:
24)对于每个粒子,将粒子经历过的最优位置作为粒子个体最优点,最优位置为适应度函数最小时对应的粒子位置。
25)对于整个群体,将所有粒子个体所经历过位置中的最优位置作为种群最优点。
26)基于速度更新公式和位置更新公式更新每个粒子的速度和位置。速度更新公式和位置更新公式分别为:
vi+1=wivi+c1·r1·(pbesti-ki)+c2·r2·(gbesti-ki)
ki+1=ki+vi
其中,vi表示第i次迭代粒子的速度,ki表示第i次迭代粒子的位置,r1、r2是一个[0,1]之间的随机数,pbesti是第i次迭代粒子个体最优点,gbesti是第i次迭代种群最优点。
27)判断迭代次数是否达到迭代次数阈值,如未达到,重复步骤22)到步骤26),否则算法结束。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (7)
1.一种非线性动力吸振器的参数优化方法,其特征在于,包括如下步骤:
1)基于非线性动力吸振器的结构,建立非线性动力吸振器的动力学模型,量化吸振效果;
2)将吸振效果的量化函数作为目标函数,采用变惯性权重粒子群优化算法进行全局寻优。
2.根据权利要求1所述的非线性动力吸振器的参数优化方法,其特征在于,步骤2)中采用变惯性权重粒子群优化算法进行全局寻优包括如下步骤:
21)初始化粒子群算法因子、粒子群位置和速度;
22)基于惯性权重的计算表达式和迭代次数,计算惯性权重;
23)基于粒子的位置,计算每个粒子的适应度函数值;
24)对于每个粒子,将粒子经历过的最优位置作为粒子个体最优点;
25)对于整个群体,将所有粒子个体所经历过位置中的最优位置作为种群最优点;
26)基于速度更新公式和位置更新公式更新每个粒子的速度和位置;
27)判断迭代次数是否达到迭代次数阈值,如未达到,重复步骤22)到步骤26),否则算法结束。
3.根据权利要求2所述的非线性动力吸振器的参数优化方法,其特征在于,所述步骤21)中,粒子群算法因子包括种群大小n、加速度因子c1、加速度因子c2以及迭代次数阈值H;粒子群位置和速度都是一个3维向量,位置为k=(k1,k2,k3),速度为v=(v1,v2,v3),一个位置代表一个潜在解。
6.根据权利要求5所述的非线性动力吸振器的参数优化方法,其特征在于,所述步骤24)中,最优位置为适应度函数最小时对应的粒子位置。
7.根据权利要求6所述的非线性动力吸振器的参数优化方法,其特征在于,所述步骤26)中,速度更新公式和位置更新公式分别为:
vi+1=wivi+c1·r1·(pbesti-ki)+c2·r2·(gbesti-ki)
ki+1=ki+vi
其中,vi表示第i次迭代粒子的速度,ki表示第i次迭代粒子的位置,r1、r2是一个[0,1]之间的随机数,pbesti是第i次迭代粒子个体最优点,gbesti是第i次迭代种群最优点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010091796.1A CN111339611A (zh) | 2020-02-13 | 2020-02-13 | 一种非线性动力吸振器的参数优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010091796.1A CN111339611A (zh) | 2020-02-13 | 2020-02-13 | 一种非线性动力吸振器的参数优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111339611A true CN111339611A (zh) | 2020-06-26 |
Family
ID=71185195
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010091796.1A Pending CN111339611A (zh) | 2020-02-13 | 2020-02-13 | 一种非线性动力吸振器的参数优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111339611A (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103455728A (zh) * | 2013-09-13 | 2013-12-18 | 天津大学 | 一种基于机械加工过程的动力吸振器参数调谐优化方法 |
US20150375304A1 (en) * | 2014-06-30 | 2015-12-31 | Kennametal Inc. | Vibration Absorber with High Viscous Dampening |
CN106326561A (zh) * | 2016-08-25 | 2017-01-11 | 同济大学 | 一种新型外旋轮线型离心摆吸振器设计方法 |
CN108615068A (zh) * | 2018-03-24 | 2018-10-02 | 西安电子科技大学 | 一种混沌扰动与自适应惯性权重的粒子群优化方法 |
CN109255426A (zh) * | 2018-08-17 | 2019-01-22 | 长安大学 | 一种权重动态调整的地质断层参数粒子群优化搜索算法 |
CN109977520A (zh) * | 2019-03-20 | 2019-07-05 | 吉林大学 | 一种二自由度动力吸振器联合优化方法 |
-
2020
- 2020-02-13 CN CN202010091796.1A patent/CN111339611A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103455728A (zh) * | 2013-09-13 | 2013-12-18 | 天津大学 | 一种基于机械加工过程的动力吸振器参数调谐优化方法 |
US20150375304A1 (en) * | 2014-06-30 | 2015-12-31 | Kennametal Inc. | Vibration Absorber with High Viscous Dampening |
CN106326561A (zh) * | 2016-08-25 | 2017-01-11 | 同济大学 | 一种新型外旋轮线型离心摆吸振器设计方法 |
CN108615068A (zh) * | 2018-03-24 | 2018-10-02 | 西安电子科技大学 | 一种混沌扰动与自适应惯性权重的粒子群优化方法 |
CN109255426A (zh) * | 2018-08-17 | 2019-01-22 | 长安大学 | 一种权重动态调整的地质断层参数粒子群优化搜索算法 |
CN109977520A (zh) * | 2019-03-20 | 2019-07-05 | 吉林大学 | 一种二自由度动力吸振器联合优化方法 |
Non-Patent Citations (4)
Title |
---|
TIGLI O F 等: "Optimum vibration absorber (tuned mass damper) design for linear damped systems subjected to random loads", 《JOURNAL OF SOUND AND VIBRATION》 * |
张多: "动力吸振器参数优化及其主动控制技术研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 * |
贾安兵 等: "随机激励下动力吸振器参数优化方法", 《山东理工大学学报(自然科学版)》 * |
邱晓燕 等: "计及智能电网多不确定性的多目标储能优化", 《计及智能电网多不确定性的多目标储能优化》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sun et al. | Finite frequency $ H_ {\infty} $ control for vehicle active suspension systems | |
CN109725537B (zh) | 一种改进的线性最优半主动控制方法 | |
CN107977497B (zh) | 一种电动轮轮内减振系统参数优化方法 | |
CN110110380B (zh) | 一种压电执行器迟滞非线性建模方法及应用 | |
CN113591360B (zh) | 基于整车动力学模型的磁流变阻尼器结构参数优化方法 | |
CN115866621B (zh) | 一种基于鲸鱼算法的无线传感器网络覆盖方法 | |
CN108416178A (zh) | 一种机械弹性车轮结构参数优化设计方法 | |
CN110001336B (zh) | 一种基于add正实网络优化的车辆isd悬架一阶理想模型 | |
CN114936518A (zh) | 一种拉伸/压缩弹簧设计参数求解方法 | |
CN115342812A (zh) | 一种基于改进蝴蝶优化算法的无人机三维航迹规划方法 | |
CN111339611A (zh) | 一种非线性动力吸振器的参数优化方法 | |
CN111439086B (zh) | 一种基于粒子群优化的车辆isd悬架pdd控制理想模型 | |
Panda et al. | Reduction of linear time-invariant systems using Routh-approximation and PSO | |
CN112561024A (zh) | 一种跌倒检测方法及其系统、通信设备和存储介质 | |
Panda et al. | Model reduction of linear systems by conventional and evolutionary techniques | |
CN111452579A (zh) | 基于悬架不变点理论的车辆半主动自适应控制系统及方法 | |
Hada et al. | Optimisation of an active suspension force controller using genetic algorithm for random input | |
CN115826402A (zh) | 一种基于深度强化学习算法的主动悬架控制策略生成方法 | |
Ma et al. | Optimization design of semi-active controller for in-wheel motors suspension | |
CN107992671B (zh) | 一种基于生物遗传算法的智能机器人调频方法 | |
Zhang et al. | A hybrid max–min ant system by levy flight and opposition-based learning | |
CN117421999B (zh) | 一种建筑消能黏滞阻尼器参数优化设计方法 | |
CN117841591B (zh) | 基于改进模糊神经网络pid的isd悬架控制方法 | |
Liu et al. | Optimal control for automotive seat suspension system based on acceleration based particle swarm optimization | |
CN111923679B (zh) | 基于稳定性切换理论的车辆isd半主动悬架pid时滞切换补偿控制方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200626 |