配气机构的动力学分析

1、配气机构动力学分析课程设计目 录一、配气机构的机构简图1二、配气机构运动学计算分析11)配气机构中间参数法的代数分析12）运初始值的设定及简化计算3三、配气机构动力学计算分析81)受力分析及微分方程的建立82）配气机构质量的换算及方程参数的计算103)动力学微分方程的求解12四、配气机构动力学优化比较16参考文献：23附件：241配气机构的运动学和动力学分析1、 配气机构的机构简图其自由度为主动件为凸轮轴，输出件为气门。2、 配气机构的运动学计算分析 1、配气机构中间参数法的代数分析 由上面的机构简图可以得到，摇臂轴与凸轮轴的竖直位移为： 化简得到： （1）摇臂轴与凸轮轴的水平位移：化简得到：

2、 （2）上面（1）（2）两式对时间求导得到解得 其中，分别为摇臂和推杆的角速度，两式对时间求导得到摇臂和推杆的角加速度为：同理，得到推杆的角加速度为其中即为挺柱和推杆长度比根据机构简图上的几何关系， 对时间求导可以得到将摇臂的角速度，角加速度带入可以得到：气门传动机构的传动比对中间参数进行线性近似可以得到 2、运动初始值的设定及运动学计算的简化计算初始参数的设定：凸轮轴转速：1000r/min 故运动开始时推杆与竖直位置成5度角，摇臂水平且摇臂轴两端摇臂成一条直线（即机构简图中所示和在一条直线上），故 ，。由此可得气门传动机构的传动比计算可简化为 气门传动机构的传动比：1.385挺柱最大升程7

3、.42mm初始时刻挺柱升程气门最大升程可计算得为10.28mm根据摇臂的空间尺寸我们计算出摇臂两端的长度分别为：由于上文的计算分析采用了中间参数（即摇臂轴和数值方向的夹角）的分析方法，因而应首先近似计算该参数值。（1） 式左侧可以近似为故根据初始设定值并画图调试得到较精确的中间参数计算值的简化计算式为 （3）上式即为中间参数的计算式下面进行配气机构各部件运动参数的简化计算，由于推杆的偏转角度一般比较小，因而可以近似认为。因，所以为一个小量，又接近于直角，所以由推杆角速度的表达式知道很小（表达式中存在高阶小量）根据高次方凸轮型线的设计得到在凸轮上升段、下降段挺柱的运动参数分别为：A、 上升缓冲段

4、（） 上升基本工作段（）B、下降基本工作段（） 下降缓冲段（）下面进行各运动件的运动学简化计算，这里均取以简化计算A、推杆角速度和角加速度的简化计算代入上面的挺柱升程与凸轮轴转角关系，并用Matlab计算绘图得到可见摆杆的角速度很小，几乎为零，因而在复杂计算中可以忽略推杆的角速度从而简化计算。推杆的角加速度也不大的，但要考虑其往复摆动的惯性力并进行校核计算。B、摇臂角速度和角加速度的简化计算摇臂由于一直做高速往复摆动，角速度和角加速都比较高。其加速度峰值出现在工作段起始和终止部分，需要做强度校核。C、气门的加速度和角加速度的简化计算气门的加速度峰值同样出现在工作段附近，在设计气门弹簧时需要考虑

5、其最大负惯性力，即最大负加速度。三、配气机构的动力学计算分析 配气机构是一个多质量的振动系统，由于只需着重研究气门的运动情况，故考虑阻尼、燃气作用力、气门座刚度、气门间隙等因素建立单质量系统。根据作用在质量M上的力平衡条件，即可建立气门运动的微分方程。其动力学计算模型如下图所示：上图中整个配气机构换算到气门端的质量为M，整个气门驱动机构的刚度为，两端分别连接凸轮轴和气门。气门弹簧刚度为，气门座刚度为，是阻尼系数。1、受力分析及微分方程的建立作用在质量M上的力描述如下：1）、气门弹簧力：其中为气门弹簧的预紧变形量2）、配气机构弹性变形力3）、气门落座时气门做的反作用力其中为气门座刚度，为气门座在

6、气门弹簧预紧力及燃气作用力作用下的初变形量。 4）、燃气作用力，即由于气缸中气体压力与排气管中气体压力差而作用在气门头部面积上的力。进气门可取为0，排气门可根据排气门开启前后气缸中气压与排气管气压的变化曲线求得。 5）、气门座摩擦力（忽略） 6）、配气机构外阻尼力，表现为配气机构中运动件与固定件之间的摩擦力 7）、配气机构内阻尼力，表现为配气机构运动件之间的摩擦根据以上各力的平衡条件，得到以下微分方程：以凸轮轴转角为自变量 代入上式得到对上面的微分方程补充初始条件：气门刚刚打开的瞬间凸轮轴转角，引入函数带入上面的微分方程得到一阶微分方程组 （*） 式中，是关于凸轮轴转角的已知函数。（*）式可运

7、用经典四阶Runge-Kutta 方法进行数值求解。初始条件为： 2、 配气机构质量的换算及方程参数的计算 a、质量换算 由于凸轮轴刚度远大于配气机构的整体刚度，故上式可简化为： 挺柱的刚度较大，可以忽略挺柱对系统刚度的影响。假设推杆在靠近挺柱一端是不动的，则推杆端质量，为推杆质量。气门端质量，为气门质量，忽略弹簧上座及锁夹质量，为弹簧质量，为摇臂的转动惯量。 b、系统刚度的计算 1）、凸轮轴刚度 其中凸轮轴材料弹性模数凸轮轴横断面直径惯性矩凸轮轴支撑轴颈内端面之间距离253mm挺柱中心线到支撑轴颈间的距离：进气凸轮， 排气凸轮 故 2）、推杆刚度 其中推杆材料的弹性模数推杆横断面积推杆长度故

8、 3）、摇臂、摇臂轴及摇臂座刚度 预估值定为 4）、单质量系统刚度进气系统换算的刚度为： 排气系统换算刚度为： 5）、气门座刚度 小组成员计算的配气机构各零件惯量参数如下： 挺柱质量135g 推杆质量 摇臂质量125.8g 摇臂转动惯量换算到气门端的当量质量 进气门质量 排气门质量 气门弹簧质量所以进气系统的当量质量：排气系统的当量质量：c、阻尼系数的确定 内阻尼一般很小，可取 外阻尼可用阻尼比来进行估算： 故进气系统外阻尼： 排气系统外阻尼：d、其他参数的确定 气门弹簧预紧力定为 气门弹簧预压缩量定为 排气门开启前后缸内和排气管压差约为4bar 所以燃气作用力 气门座初变形量： 弹簧刚度：

9、忽略气门座的摩擦力 气门间隙进气门，排气门3、 动力学微分方程的求解 初始值的设定：气门刚刚开启时刻，气门端的位移应刚刚抵消气门间隙和气门座在弹簧预紧力及燃气作用力作用下的初变形量，即排气门对于进气系统的初始值计算如下：此时的凸轮轴转角通过Matlab求解得到气门升程、气门速度和气门加速度随曲轴转角的变化关系如下：凸轮轴转角时摇臂末端并未接触气门顶部，此时仅是消除气门间隙和气门座的预变形。上图可以看出，进气门在气门开启和气门落座过程中均有速度波动，在凸轮型线中缓冲段和工作段的交界处也有速度波动，但波动量几乎不可见。下图为进气系统各类曲线的总览排气门系统的初始值计算如下：凸轮轴转角通过Matla

10、b求解得到气门升程、气门速度和气门加速度随曲轴转角的变化关系如下：同理，凸轮轴转角时排气系统的摇臂未接触排气门顶部，此时仅是消除气门间隙和气门座的预变形。上图可以看出，同样在排气门起座和落座的时，气门出现速度波动，由于燃气作用力的影响等，排气门的速度波动明显大于进气门。在凸轮型线缓冲段和工作段的交界处，也存在速度波动，但相对而言波动很小。下图为排气系统各类曲线：四、配气机构动力学优化比较进气门速度的理想运动学曲线和动力学分析的速度曲线比较进气门最大升程10.3496mm进气门最大速度1.693m/s进气门最大加速度830.6730m/s2动力学升程的响应慢于运动学，且由于气门间隙的影响，运动学

11、峰值高出约0.2mm排气门最大升程10.1789mm排气门最大速度1.7m/s排气门最大加速度835.9052m/s2同样，动力学升程的响应慢于运动学，且由于气门间隙的影响，运动学峰值高出约0.3mm下面为运用2100发动机原有凸轮型线得到的进排气系统升程曲线、速度曲线和加速度曲线。进气系统：排气系统：通过对比原2100配气系统动力学循环曲线和运用高次方凸轮型线后的2100配气系统动力学循环曲线可以得到以下结论：1、 在进排气升程变化不大的情况下，采用多项式高次方凸轮型线可以显著的减少进排气系统的速度和加速度幅值，速度由2.5m/s左右降至1.7m/s左右，加速度幅值由1900m/s2左右降至

12、800m/s2左右。2、 采用高次方凸轮型线使得进排气系统在开启和落座时刻的速度、加速度波动幅度明显减少，速度几乎没有波动。这样在使用材料和加工工艺相同的情况下保证了配气系统具备较高的使用寿命，较低是噪声，同时可以降低了成本。参考文献1李佳, 刘震涛, 刘忠民, 等. 内燃机配气机构系统刚度的虚拟设计J. 农业工程学报, 2012, 28(6): 44-49.2赵冬青, 苏铁熊, 赵振锋, 等. 顶置凸轮轴配气机构运动学和动力学计算J. 车用发动机, 2004 (6): 1-3.3史绍熙. 柴油机设计手册J. 北京: 中国农业机械出版社, 1984.4高文志, 李明海, 袁文华, 等.内燃机课

13、程设计J. 北京: 中国水利水电出版社, 2010.85内燃机设计附件运动学计算的Matlab代码 x=linspace(0,174,1740); y=0.25*(1-cosd(3.67*x).*(x>=0&x<24.5)+(7.67-13.780*(87-x)/62.5).2+8.257*(87-x)/62.5).6-2.343*(87-x)/62.5).14+0.447*(87-x)/62.5).22).*(x>=24.5&x<87)+(7.67-13.780*(x-87)/62.5).2+8.257*(x-87)/62.5).6-2.343*(x-

14、87)/62.5).14+0.447*(x-87)/62.5).22).*(x>=87&x<149.5)+0.25*(1-cosd(3.67*(174-x).*(x>=149.5&x<174); y1=0.9175*sind(3.67*x).*(x>=0&x<24.5)+(25.265*(87-x)/62.5)-45.417*(87-x)/62.5).5+30.071*(87-x)/62.5).13-9.015*(87-x)/62.5).21).*(x>=24.5&x<87)+(-25.265*(x-87)/62.

15、5)+45.417*(x-87)/62.5).5-30.071*(x-87)/62.5).13+9.015*(x-87)/62.5).21).*(x>=87&x<149.5)+(-0.9175*sind(3.67*(174-x).*(x>=149.5&x<174); y2=3.367*cosd(3.67*x).*(x>=0&x<24.5)+(-23.161+208.18*(87-x)/62.5).4-358.372*(87-x)/62.5).12+173.551*(87-x)/62.5).20).*(x>=24.5&x&

16、lt;87)+(-23.161+208.18*(x-87)/62.5).4-358.372*(x-87)/62.5).12+173.551*(x-87)/62.5).20).*(x>=87&x<149.5)+(3.367*cosd(3.67*(174-x).*(x>=149.5&x<174); wt=104.72*y1./(248*tan(acos(cosd(85)-y/34.55); plot(x,wt,'b') et=-y2*104.722./(248*tan(acos(cosd(85)-y/34.55)+0.14*(104.72*y1/34.55).2./(sin(acos(cosd(85)-y/34.55).3; plot(x,et,'r') wy=y1*104.72./(34.55*sin(acos(cosd(85)-y/34.55); plot(x,wy,'b')ey=y2*104.722./(34.55*sin(acos(cosd(