1. 引言
近年来,随着我国高铁、地铁、城市轻轨、城际铁路等重型载运交通工具的快速发展,在便捷交通、促进经济和社会快速发展的同时,产生的强振动对周围环境也造成了很大扰动,特别是地铁以及市域轨道交通等。此外,一些重型动力装备,是现代工业生产和科学研究的关键工艺装备和关键技术手段,如锻锤、冲压装备、振动台等,这些装备会对周围环境、人员舒适造成影响,甚至造成附属建筑结构破坏。
采取振动控制措施是解决重型载运工具以及动力装备对周围环境振动影响的重要手段,但许多工程实践表明,当对地铁等重型载运工具采取浮置板道床振动控制措施时,其对周围环境的扰动会显著减小,但车厢内的人员舒适可能会下降,是一个矛盾问题。另外,当对重型动力装备采取振动控制措施时,如基础间接隔振、装备直接隔振,动力装备对周边环境的扰动显著减小,但装备自身的振动可能会变大,不利于装备的安全长期使用,这也是一个矛盾问题。
如何采取有效的振动控制措施,对于重型载运工具,在振动控制效果以及车厢人员舒适两者间获得较为均衡的效果,对于重型动力装备,既能保证对环境影响的振动控制效果,也能兼顾装备自身的振动水平,是当前实际工程中面临的一个迫切优化问题。
Eberhart和Kennedy于1995年首次提出了一种新型的群智能优化算法——粒子群优化(Particle Swarm Optimization, PSO)算法[1],其主要思想是基于粒子间的合作和竞争以找出最优解(粒子),该算法具有简单、易实现、收敛快且可调参数少等优点,Coello C A C等[2]首次提出了多目标粒子群优化算法(Multi-Objective PSO, MOPSO),目前已广泛应用于工程领域的优化计算[3]-[5]。
本文将基于多目标粒子群优化算法,对重型载运工具及重型动力装备开展多目标振动控制研究,以期解决前述的矛盾问题。
2. 多目标粒子群优化算法简介
多目标粒子群优化算法(MOPSO)的主要计算方法[6],可表述如下:
第1步:初始化粒子群,计算各对应粒子的目标函数向量,并将其中的非劣解加入到外部档案之中;
第2步:初始化粒子群的局部最优解pbest和全局最优解gbest;
第3步:在搜索空间内,通过(1)和(2)式不断调整粒子的飞行速度和位置,形成新的pbest;
(1)
(2)
第4步:根据新的非劣解维护外部档案,并为每个粒子选取gbest (档案的内容将决定全局最优解的选取);
第5步:是否达到最大迭代次数,否的话继续计算,是则停止计算,并输出Pareto最优解集。
MOPSO与单目标粒子群优化算法的主要区别就是全局最优解的选取方式以及外部档案的设定和更新;其中,外部档案的更新将直接决定全局最优解,如图1示意,图中,NDs表示非支配解集,s1~s5表示一组解集。在多目标优化中,最直接的方法就是利用Pareto支配的概念,考虑档案中的所有非劣解,并从中确定一个“主导者”,从而形成一组Pareto前沿,如图2所示。
Figure 1. Update of the external file
图1. 外部档案的更新
Figure 2. Pareto frontier
图2. Pareto前沿
3. 重型动力装备多目标振动问题描述
Figure 3. Vibration control of heavy power equipment
图3. 重型动力装备振动控制
采取隔振措施的重型动力装备的振动控制示意如图3所示,其中,m为浮置板道床及上部轿厢、乘客等的总质量,k、c为隔振体系的刚度和阻尼。对于环境扰动的评判,可以动力装备传递至地面的力来评价;对于装备自身的振动水平,可以其加速度来评价。假设装备的振动位移、速度和加速度分别为
,则传递至地面的力可表示为,则多目标优化问题即可表示为对:目标1——以及目标2——
的多目标优化。
4. 重型载运工具多目标振动问题描述
采取浮置板道床的重型载运轨道交通工具的振动控制示意如图4所示,其中,m为浮置板道床及上部轿厢、乘客等的总质量,k、c为浮置板隔振体系的刚度和阻尼。对于环境扰动的评判,可以轨道交通传递至地面的力来评价;对于轿厢中的人体舒适评判,有多种评价指标,如振动频率、加速度以及振动强度等[7]-[11],本文以加速度来评价。假设浮置板道床及上部轿厢体系的振动位移、速度和加速度分别为
,则轨道交通传递至地面的力可表示为,其多目标优化问题同重型动力装备。
Figure 4. Vibration control of heavy-duty vehicles
图4. 重型载运工具振动控制
5. 算例
由于重型动力装备及重型载运工具的简化振动控制模型相同(计算框图如图5所示),本文考虑两种荷载工况,一种是幅值为5000 N、振动频率为25 Hz的简谐振动荷载,另一种是幅值为5000 N的随机振动荷载。
Figure 5. Vibration control calculation block diagram of heavy power equipment and heavy-duty vehicles
图5. 重型动力装备及重型载运工具振动控制计算框图
设总质量m为10000 kg,隔振体系的刚度k和阻尼c未知,根据上述描述的多目标振动问题,对隔振参数进行寻优计算。根据常规隔振设计经验,浮置板道床的隔振体系固有频率可按3~15 Hz范围寻优,阻尼比的寻优范围按固有频率范围作相应计算。
MOPSO的算法设置:最大迭代次数100,最大、最小惯性权重系数ωmax = 0.9、ωmin = 0.4,惯性权重系数按式
变化,i为第i次迭代,m为最大迭代次数,加速因子c1和c2分别取1.0、0.4。
1) 简谐振动荷载的计算结果:
Figure 6. Pareto frontier based on multi-objective calculation (harmonic vibration load)
图6. 多目标计算得到的Pareto前沿(简谐振动荷载)
Figure 7. Relation between designed natural frequency and damping ratio based on multi-objective calculation solution set (harmonic vibration load)
图7. 基于多目标计算解集得到的设计固有频率与阻尼比关系(简谐振动荷载)
多目标优化计算的Pareto前沿如图6所示,根据解集计算的固有频率与阻尼比的对应关系如图7所示。由图6中的Pareto前沿,A区域中的解集相对较为理想,在实际工程中,可以参考选取,其中,A区域对应图7的C区域。以A区域中的某点解为例,其传递力为652.7 N,装备加速度为0.5 m/s2,则对应的设计固有频率为3 Hz,设计阻尼比为0.236。
2) 随机振动荷载的计算结果:
Figure 8. Pareto frontier based on multi-objective calculation (random vibration load)
图8. 多目标计算得到的Pareto前沿(随机振动荷载)
Figure 9. Relation between designed natural frequency and damping ratio based on multi-objective calculation solution set (random vibration load)
图9. 基于多目标计算解集得到的设计固有频率与阻尼比关系(随机振动荷载)
多目标优化计算的Pareto前沿如图8所示,根据解集计算的固有频率与阻尼比的对应关系如图9所示。以Pareto前沿中的某点解为例,其传递力为7856.5 N,装备加速度为0.84 m/s2,则对应的设计固有频率为9.95 Hz,设计阻尼比为0.35。
6. 小结
本研究基于多目标粒子群优化算法对重型载运工具及重型动力装备中常见的一类“矛盾”振动问题进行了优化,对于重型载运交通工具,可以通过优化设计,兼顾对轿厢人员舒适以及环境扰动的影响,对于重型动力装备,可以实现对装备自身振动水平以及环境扰动影响的平衡考虑,是对振动控制优化设计的重要补充,对于实际工程具有较强的指导意义。
文中考虑了两种振动荷载,结果表示不同的振动荷载其多目标优化设计的结果和规律有明显不同,如果参数设计合理,不仅可以有效兼顾轿厢人员舒适(装备自身振动影响)以及传递至环境的振动水平,还可以将二者均控制在较为理想的水平,如果参数选择和设计不利,可能导致某个目标控制效果明显不足。故在实际工程中,对于不同振动荷载形式的装备或载运工具,应开展针对性设计,避免按照经验设计导致的振动控制效果不佳,甚至造成不利影响。
未来,笔者还将基于实际工程中面临的类似矛盾振动控制问题开展多目标优化设计,如精密工程中的微振动控制问题等。
项目基金
国家重点研发计划“地下城轨交通综合体振动及结构安全控制方法研究与应用”(课题编号:2024YFF0508104)。中国机械工业集团青年基金重点项目“大科学工程群微纳级环境振动控制关键技术研究与应用”。