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文件：function fit＝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。

7.如权利要求6所述的推力杆球铰的优化设计方法，其特征在于：所述步骤4-2)中，所述代码如下：clear

clc

options＝gaoptimset('PopulationSize',400,'Generations',Inf,'display','diagnose')；％,'OutputFcn',@outputfcn LB＝[0,0,j,h,M,a,c,f]；％Sr,Sa,l,d,r,Kr,Ka,lambda_kUB＝[10,10,k,i,N,b,e,g]；

[F,Fval]＝gamultiobj(@Thrust_rod_optim,5,[],[],[],[],LB,UB,[],options)；

％for i＝1:length(Fval)％plot3(Fval(i,1),Fval(i,2),Fval(i,3),'o')％hold on

％end

plot3(Fval(:,1),Fval(:,2),Fval(:,3),'o')xlabel('塑料层截面面积')

ylabel('刚度优化目标')

zlabel('硫化面积')。

8.如权利要求7所述的推力杆球铰的优化设计方法，其特征在于：所述步骤4)中，在计算完成后，根据较优解的分布状态及趋势，找到三个子目标函数的数值均最小的那个较优解，此较优解即为问题的最优解。