一种推力杆球铰的优化设计方法技术

本发明专利技术涉及一种推力杆球铰的优化设计方法，包括如下工艺步骤：1)，选取优化变量及基因编码：2)，确定约束条件：3)，确定适应度函数(目标函数)：橡胶衬套体积优化的优化子目标定义为F1，球铰刚度优化的优化子目标定义为F2，硫化面积的优化子目标为F3；4)，优化计算：将F1、F2和F3作为各子种群的适应度函数，并在MATLAB遗传算法工具箱中采用混合法来进行优化计算，即可得到较优解的分布状态，并根据较优解的分布状态得到最优方案。本发明专利技术的推力杆球铰的优化设计方法不但没有大幅度改变推力杆的原有结构设计方案，控制了改进成本，而且使球铰径向刚度和轴向刚度的配比落入合理的区间，兼顾了整车性能，适宜工程应用。

An optimal design method of the ball joint of the thrust bar

【技术实现步骤摘要】
一种推力杆球铰的优化设计方法
本专利技术涉及一种橡胶衬套的设计方法，具体涉及一种推力杆球铰的优化设计方法，其应用于汽车悬架及轨道客车转向架的推力杆球铰中，属于推力杆球铰

技术介绍
推力杆是商用车底盘系统或轨道客车转向架系统中广泛应用的关键零部件，用于传递车桥与车架之间的纵向、侧向载荷及部分扭矩，并衰减车桥传来的部分振动能量。在推力杆服役过程中，其球铰中的橡胶衬套经常发生疲劳破坏的现象，这将导致商用车底盘系统的导向功能部分失效，并严重威胁行车安全。要显著提高推力杆球铰的疲劳寿命，除了要提升橡胶材料的力学性能外，还需对其球铰结构进行有针对性的优化设计。然而，推力杆的球铰结构复杂，球铰中的橡胶衬套不但具有变截面特征，而且其变形也具有非线性。此外，橡胶衬套被金属外套包裹，因此无法观察橡胶衬套的实际失效过程，导致推力杆球铰的性能优化变得非常困难和棘手。目前，国内外学者已对推力杆球铰应力分析，推力杆纵向刚度优化方法及修正方法以及疲劳寿命预测方法进行了深入研究。然而，相关研究仍然无法为推力杆的球铰结构优化提供理论和方法指导，导致推力杆的球铰设计具有一定的盲目性。要科学地优化推力杆的球铰结构，需解决如下关键问题：(1)准确理解橡胶衬套的失效机理；(2)确定球铰结构参数与其性能参数之间的关系；(3)针对性地建立球铰结构优化模型，并通过优化设计显著提高球铰的疲劳寿命。由于橡胶部件的刚度不仅影响其疲劳寿命，而且影响整车性能。因此，应在优化球铰应力应变分布状态的同时，兼顾其刚度特性。此外，橡胶结构具有非线性和变截面特征。因此，有限元分析是实现橡胶结构优化的基本途径。随着研究的深入，综合应用有限元法及遗传算法来完成橡胶部件的优化设计取得了成功，这为推力杆球铰的优化设计指明了方向。推力杆球铰的主要疲劳破坏形式为球铰橡胶衬套的疲劳失效。推力杆球铰的内部结构如图1所示。在球铰内部，橡胶衬套1通过硫化的方式与端盖2、芯轴3、塑料层4和金属外套固连。塑料层与芯轴之间为间隙配合，用于给橡胶衬套提供用于应力释放的旋转自由度。由于压缩预应力能促进橡胶材料疲劳裂纹的闭合，对橡胶衬套疲劳寿命有利，因此橡胶衬套自身也具有一定的过盈量，从而在装配过程中产生了一定的预压缩量。因此，影响橡胶衬套疲劳寿命的因素主要包括外部载荷、橡胶材料力学性能、环境温度、球铰结构和橡胶衬套预压缩量。受安装环境及成本限制，只有优化球铰结构及调整预压缩量在工程实际中是可行的。根据有限元球铰失效机理的分析结果，提高橡胶衬套疲劳寿命的关键在于限制橡胶衬套端部的延展变形并降低橡胶衬套环面的变形量。由于橡胶材料具有非线性且橡胶衬套为不规则结构，因此无法通过常规的数值计算方法准确量化橡胶衬套的变形量。然而，在球铰内部空间和外部载荷一定的前提下，橡胶衬套端部延展变形量的降低宏观表现为推力杆球铰轴向刚度的增加，而橡胶衬套环面变形量的降低宏观表现为推力杆纵向刚度的增加。因此，可用推力杆球铰轴向刚度和推力杆纵向刚度来间接量化橡胶衬套的变形量。
技术实现思路
为解决上述问题，本专利技术的目的在于提供一种推力杆球铰的优化设计方法，其能提高推力杆球铰的刚度，延长了使用寿命。为实现上述目的，本专利技术采取的技术方案为：一种推力杆球铰的优化设计方法，其包括如下工艺步骤：1)，选取优化变量及基因编码：选取塑料层的厚度d，宽度l、塑料层圆角半径r，以及橡胶衬套的轴向预压缩量(Sa)和径向预压缩量(Sr)这5个优化变量；选取的基因型是长度为5的双精度浮点数；2)，确定约束条件：3)，确定适应度函数(目标函数)：其中，橡胶衬套体积优化的优化子目标定义为：其中，R为芯轴半径；球铰刚度优化的优化子目标定义为：F2＝(kr-Kr)2+λk(ka-Ka)2其中，Kr和Ka分别是球铰径向刚度和轴向刚度的设计平均值；硫化面积的优化子目标为：其中，L为橡胶衬套与芯轴硫化区域的切面弧长；4)，优化计算：将F1、F2和F3作为各子种群的适应度函数，并在MATLAB遗传算法工具箱中采用混合法来进行优化计算，即可得到较优解的分布状态，并根据较优解的分布状态得到最优方案。本专利技术的推力杆球铰的优化设计方法进一步为：所述步骤1)中，塑料层表面形状为圆柱形，且在边缘设置圆角。本专利技术的推力杆球铰的优化设计方法进一步为：所述步骤2)的约束条件具体为：其中，kr和ka分别为球铰的径向刚度和轴向刚度；a和b分别为球铰径向刚度的合理上下限值；c和e分别为球铰轴向刚度的合理上下限值；f和g分别为基于经验得到的球铰径向刚度和轴向刚度配比的合理上下限值；h和i，j和k以及M和N分别是具体球铰结构参数塑料层的厚度d、宽度l和塑料层圆角半径r的合理上下限值。本专利技术的推力杆球铰的优化设计方法进一步为：所述步骤3)中的球铰刚度优化的优化子目标中，kr＝-30.459+10.353Sr+10.518Sa-6.885d+0.396l；ka＝-43.586+1.946Sr+3.733Sa-0.336d+0.758l。且kr和ka是对大量有限元计算结果进行多项式拟合得到的关系式。本专利技术的推力杆球铰的优化设计方法还为：所述步骤4)中，在MATLAB软件中按照上述优化模型建立遗传算法的适应度函数，具体按如下代码生成适应度函数的m文件：functionfit＝Thrust_rod_optim(vars)Sr＝vars(1)；Sa＝vars(2)；l＝vars(3)；d＝vars(4)；r＝vars(5)；R＝A1；％芯轴半径；fit1＝1./((40+d)*l+(pi/2-2)*r.^2-l/2*sqrt(40^2-l^2/4)-40^2*asin(l/80))；％％Kr＝vars(6)；Ka＝vars(7)；lambda_k＝vars(8)；Kr＝-30.459+10.353*Sr+10.518*Sa-6.855*d+0.396*l；Ka＝-43.586+1.946*Sr+3.733*Sa-0.336*d+0.758*l；lambda_k＝Kr./Ka；Kr_goal＝A2；％设计刚度Ka_goal＝A3；fit2＝(Kr-Kr_goal)^2+abs(lambda_k)*(Ka-Ka_goal)^2；％％L＝A4；fit3＝1./(L+l+2*d+2*R-4*r+pi*r-2*R*asin(l/(2*R))-sqrt(4*R^2-l^2))；％％fit＝[fit1；fit2；fit3]；其中，A1为芯轴半径R的已知量，A2和A3分别是球铰径向刚度和轴向刚度的设计理想值Kr_goal和Ka_goal的已知量，A4为橡胶衬套与芯轴硫化区域的切面弧长L的已知量。在设计要求一定的前提下，A1、A2、A3和A4均可指本文档来自技高网...

【技术保护点】
1.一种推力杆球铰的优化设计方法，其特征在于：包括如下工艺步骤：/n1)，选取优化变量及基因编码：选取塑料层的厚度d，宽度l、塑料层圆角半径r，以及橡胶衬套的轴向预压缩量(Sa)和径向预压缩量(Sr)这5个优化变量；选取的基因型是长度为5的双精度浮点数；/n2)，确定约束条件：/n3)，确定适应度函数(目标函数)：其中，橡胶衬套体积优化的优化子目标定义为：/n

【技术特征摘要】
1.一种推力杆球铰的优化设计方法，其特征在于：包括如下工艺步骤：
1)，选取优化变量及基因编码：选取塑料层的厚度d，宽度l、塑料层圆角半径r，以及橡胶衬套的轴向预压缩量(Sa)和径向预压缩量(Sr)这5个优化变量；选取的基因型是长度为5的双精度浮点数；
2)，确定约束条件：
3)，确定适应度函数(目标函数)：其中，橡胶衬套体积优化的优化子目标定义为：



其中，R为芯轴半径；
球铰刚度优化的优化子目标定义为：
F2＝(kr-Kr)2+λk(ka-Ka)2
其中，Kr和Ka分别是球铰径向刚度和轴向刚度的设计平均值；
硫化面积的优化子目标为：



其中，L为橡胶衬套与芯轴硫化区域的切面弧长；
4)，优化计算：将F1、F2和F3作为各子种群的适应度函数，并在MATLAB遗传算法工具箱中采用混合法来进行优化计算，即可得到较优解的分布状态，并根据较优解的分布状态得到最优方案。


2.如权利要求1所述的推力杆球铰的优化设计方法，其特征在于：所述步骤1)中，塑料层表面形状为圆柱形，且在边缘设置圆角。


3.如权利要求1所述的推力杆球铰的优化设计方法，其特征在于：所述步骤2)的约束条件具体为：



其中，kr和ka分别为球铰的径向刚度和轴向刚度；a和b分别为球铰径向刚度的合理上下限值；c和e分别为球铰轴向刚度的合理上下限值；f和g分别为基于经验得到的球铰径向刚度和轴向刚度配比的合理上下限值；h和i，j和k以及M和N分别是具体球铰结构参数塑料层的厚度d、宽度l和塑料层圆角半径r的合理上下限值。


4.如权利要求1所述的推力杆球铰的优化设计方法，其特征在于：所述步骤3)中的球铰刚度优化的优化子目标中，
kr＝-30.459+10.353Sr+10.518Sa-6.885d+0.396l；
ka＝-43.586+1.946Sr+3.733Sa-0.336d+0.758l。


5.如权利要求1所述的推力杆球铰的优化设计方法，其特征在于：所述步骤4)的优化计算方法具体包括如下步骤：
4-1)，在MATLAB软件中按照上述优化模型建立遗传算法的适应度函数，
4-2)，通过gaoptimset指令调用上述适应度函数的m文件及遗传算法工具箱，并在代码中指定初始种群个体数量等相关初始参数，并制定输出图像的作图指令。


6.如权利要求5所述的推力杆球铰的优化设计方法，其特征在于：所述步骤4-1)中，按如下代码生成适应度函数的m文件：
functionfit＝Thrust_rod_optim(vars)
Sr＝vars(1)；
Sa＝vars(2)；

【专利技术属性】
技术研发人员：柯俊，高晋，徐敬恩，
申请(专利权)人：浙江理工大学，
类型：发明
国别省市：浙江;33