1. 引言
连杆在发动机中承担连接和传动的功能,所以承受较大的周期性交变应力作用,此外为了限制连杆运动惯性力,提高连杆运动性能,需要将连杆的质量尽量降低,同时还必须满足刚度和强度的要求,防止连杆在工作中断裂失效。因此,一个重量轻且具有足够强度的连杆在现代汽车设计中备受重视 [1] 。有学者利用Matlab优化工具箱对连杆杆身的工字形结构进行可靠性优化设计并给出定量描述 [2] 。有研究人员以ANSYS软件为平台,建立了连杆的有限元模型,在满足屈服强度条件下对连杆进行结构优化设计,得出的连杆质量比传统设计的质量要小 [3] 。此外也有学者在灵敏度分析的基础上,选择合理的设计参数以连杆质量最轻为目标对连杆进行了优化设计 [4] 。在国外有学者对连杆进行拓扑优化,操纵并比较拓扑优化前后连杆的最大主应力、等效应力、最大主应变和总变形值 [5] 。也有对锻钢连杆模型进行优化,并进行一些减重活动。在优化的帮助下实现了轻量化的结果 [6] 。
参考国内外的研究现状,本文将运用优化理论以及有限元分析的理论,以连杆的质量作为目标进行最优化求解,使连杆质量更加小,同时对比优化前后连杆的最大等效应力和最大变形,保证其有足够强度,不影响发动机正常运行。
2. 优化理论
机械优化设计即在满足约束前提下,寻找一组设计参数,使机械产品单项设计指标达到最优的过程。其方法主要包括解析法和数值计算方法,解析法主要是利用微分学和变分学的理论,用于解决小型和简单的问题;而数值计算方法是利用已知的信息,通过迭代计算来逼近最优化问题的解。机械优化的数学模型如下:
(1)
其中
为设计变量,是在优化中决定的参数,在优化设计过程中优选并且动态变化。
为目标函数,
为约束条件,分别为不等式约束和等式约束。求解优化设计问题即为在由所有约束条件构成的可行域内,求使目标函数取得最小值的一个设计点。本设计的设计变量是结构的尺寸,目标函数是连杆质量的最小值,约束条件为应力约束、变形约束和结构边界约束。
3. 连杆模型建立及优化对象参数化
3.1. 模型建立
本次优化的模型是在solid works中所建立的,其连杆大头以及连杆小头的内径参照曲轴半径以及活塞销半径进行设计,其余参数均根据汽油发动机设计开发手册 [7] 中的连杆设计模块进行设计。由于完全按照实物模型进行建模,到有限元分析步骤会非常困难,所以本连杆模型在不影响主要力学性能的前提下,适当的对连杆结构进行了简化,将连杆整合成了整体式结构,删掉一些螺栓结构以及倒角倒圆。整体建模如图1所示。
Figure 1. Three dimensional model of connecting rod
图1. 连杆的三维模型图
3.2. 参数化设计
优化的三要素,设计变量、约束条件、目标函数都需要将实体化的模型转换为参数化的模型,所以这一步也是优化中较为关键的步骤。首先是设计变量,设计变量即为自变量,需要通过改变设计变量的大小来获得优化结果。本次优化对象的设计变量都是选择的连杆的尺寸大小,根据以往的研究文献,共选择六个设计变量,如图2所示,分别为连杆小头外圆直径P1 (X1),连杆大头外圆直径P2 (X2),连杆小头处工字形弧半径P6 (X3),连杆大头处工字形弧半径P5 (X4),连杆大头结构长P11 (X5),连杆大头结构宽P12 (X6)。
Figure 2. Parameterization of design variables
图2. 设计变量示意图
同样的,目标函数也需要进行参数化,本次优化的目标函数是连杆的质量,目标函数的表达形式如下:
(2)
其中m1、m2、……、mn为连杆各部分质量。本次优化的发动机连杆材料使用40 Mn2S,其参数性能如表1所示,将连杆材料定义完成后,在ANSYS中即可获得连杆的质量M,计算结果显示,连杆的整体质量为5.5716 kg,而后将此数据参数化。
Table 1. Performance table of 40 Mn2S
表1. 40 Mn2S性能表
约束条件参数化同上面一样,连杆的约束主要是应力约束条件,变形约束条件以及边界约束条件。
1) 应力约束条件:
(3)
式中:σs为材料屈服极限,n为安全系数,本设计n取1.5。
2) 变形约束条件:
(4)
式中:ΔRmax表示大头内孔半径的允许最大变形量,D1为连杆的轴承与曲柄销进行初始配合的间隙,Δrmax表示小头内孔半径的最大变形量,D2表示活塞销与小头衬套进行初始配合时的间隙。
3) 边界约束条件:
为保证结构的合理性,本研究控制各设计变量在一定范围内变动。
(5)
4. 优化过程及结果
4.1. 优化过程
将连杆的优化要素参数化之后,就可以在直接优化模块进行优化。优化过程中首先要对目标参数进行设置,根据图3所示,优化名称为几何结构质量最小化,优化对象即目标参数为连杆的质量M,限制类型为连杆质量小于等于连杆质量上边界,而上边界的值为5.57 Kg。通过以上设置对连杆质量进行轻量化设计。
后续进行设计变量参数范围的设定,根据边界条件的设置,其范围如图4所示。
最后设置优化点数,本次优化选择在变量范围内设置30组设计点,软件对每组设计点数据进行计算,并且通过求解得到30组连杆的质量大小,由于需要求解的变量及生成的点数较多,所以需要长时间的运行。图5节选了30组数据中的24组,分别为设计变量X1~X6在各自约束条件内的部分取值点,并且对这些取值点进行排列,每组数据都将获得各自的质量值。
Figure 5. Connecting rod design point data diagram
图5. 连杆设计点数据图
4.2. 优化结果
通过30组数据的计算,最终会得到30个连杆的质量数值,我们可以通过对比得到连杆质量最小的点,此外直接优化模块给出了3组最优参考点供我们选择,如图6所示。
由于本优化的目的是轻量化,所以选择连杆质量最小的一组解,即参考点1,并且依据数据进行重新建模。连杆杆身各设计变量优化前后如表2所示。
Table 2. Design variable optimization results
表2. 设计变量优化结果
从表中我们可以看出,优化前后数据有增有减,增加的是连杆大头结构长,其余尺寸均有减小,并且设计变量的尺寸越大,其变化幅度也越大,连杆在优化前的质量为5.5716 kg,优化后质量减小为5.1147 kg,降低了8.2%。
5. 优化前后静力学性能
将模型导入静力学分析模块后,依次进行材料定义,网格划分,施加约束以及载荷。连杆载荷的施加主要考虑两种情况,最大拉升工况和最大压缩工况。最大拉伸工况是指出现在活塞接近排气行程上止点时,连杆所处的工作状况;最大压缩工况取气缸最大爆发压力下,连杆所处的工作状况 [8] ,因为压缩工况的力比较大,所以只对连杆受压进行分析,气体压力计算公式为:
(6)
式中Pmax为燃气最大爆发压力,D为气缸直径,通过计算及资料查阅,将压力施加到连杆小头孔的180度范围内,所施加的压力为10 MPa,设置如图7所示。并在连杆大头孔设置固定支撑。
通过计算得出连杆所受的最大等效应力为131.5 MPa,发生在连杆小头与连杆杆身的交界处,最大位移为0.07 mm,发生在连杆小头上部。图8为其应力应变云图。
同样,利用优化后的结构数据重新生成模型,导入到ANSYS中进行静力学分析,在载荷和约束条件不变的情况下得出了优化后模型的最大等效应力为139 MPa,同样发生在连杆小头与杆身的交界处,最大变形为0.005 mm。图9为优化后的有限元云图。
同原始模型相比,优化后模型的最大等效应力略微增加,数值上增加了7.5 MPa,通过应力约束条件可知,最大应力应小于413 MPa,优化后的模型应力值远小于许用应力,所以符合条件。而优化后模型的变形相较于优化前降低了0.065 mm,所以受力性能较好。
6. 结论
1) 通过workbench对发动机连杆进行优化,使其质量降低了8.2%,符合轻量化要求。
2) 优化后连杆所受最大应力增加,满足许用应力条件,在轻量化的前提下符合设计要求。最大变形减小,表明优化结构良好,静力学分析结果表明优化后连杆的静力学性能较好。
3) 使用solid work及workbench优化工具联合进行优化,在发动机连杆优化方面有一定的实用价值。