顶置凸轮配气机构气门升程的精确计算
- 格式:doc
- 大小:90.00 KB
- 文档页数:10
内燃机凸轮机构最大气门升程影响的实验分析Suat Saridemir,Hamit SaruhanDepartment of Manufacturing Engineering, Faculty of Technology, University of Düzce, Düzce, TurkeyDepartment of Mechanical Engineering, Faculty of Engineering, University of Düzce, Düzce, Turkey摘要:凸轮机构气门行程的设计对改善内燃及的动态性能是非常重要的,高速发动机不正常的振动导致随动件跳动,碰撞阀门表面和中心,进而导致凸轮系统的碰撞。
碰撞的程度取决于阀升程值和阀门关闭速度,凸轮和随动件碰撞时会引入大的力和压力。
由于不正常的振动,这可能会导致系统的早期失效。
本实验使用时域和频域分析方法研究,主要目的是研究和分析不同气门升程值和运转速度下凸轮系统的动态特性,两种最大的气门升程值为8和10毫米分别进行了运转速度为450,930,1440,1950和2430rpm条件下的实验。
从实验结果看,振幅随运转速度的增加而增加,统计分析获得的数据显示,10mm的气门升程比8mm气门升程产生更多的能量。
关键词:凸轮机构;气门行程的动态特性;振动分析;内燃机1.引言当今,越来越多的研究者热衷于提高内燃机(ICE))的性能,减少噪音及振动。
由于内燃机由许多承受大的负载的组件构成,内燃机的动态特性十分复杂[1]。
凸轮系统的运转是内燃机震动的一个震动源,凸轮从动件带动发动机气门运动,凸轮轴的旋转驱动随动件运动,定位弹簧使阀门始终压在凸轮表面。
由于弹簧力的作用,运转过程中系统会承受大的负载力。
阀弹簧的性能在凸轮从动装置的动态特性中起着重要作用,尤其是在高速运转的情况下。
因此,凸轮从动系统的动态特性主要影响因素是阀弹簧的性能。
LJ750发动机配气凸轮型线设计计算书配气机构的设计思想:配气机构是发动机的重要组成部分,一台发动机的经济性能是否优越,工作是否可靠,噪声与振动能否控制在较低的限度,常常与其配其机构的设计是否合理密切关系。
特别对于高速大功率发动机,因对其较高性能指标的设计要求,配气机构的设计及其零部件的设计和制造就更加重要和严格。
配气机构的设计涉及到如下诸多问题:1.如果根据具体的机型的要求选取合理的函数凸轮型线并进行有关特性参数的计算和分析;2.配气机构方案选择及其气门的结构布置;3.对于高速配气机构来说,必须进行配气机构的动力学计算;4.在设计过程中,需要考虑很多重要的因素,如配气相位,平稳性(包括气门速度,加速度,脉冲数值大小及有无飞脱,落座反跳等),充气性能,润滑性能,凸轮与挺拄之间的接触应力等问题;因此,进行配气机构的设计往往需要很大的工作量,特别是为了获得一个最优的设计方案。
为了简化工作量,在目前国内外的设计中,常常采用参考某些成功机型的方法,来获得设计一台新发动机所需要的配气机构。
为此,在本设计了,参考了CBR900发动机的配气机构,通过对该机型配气机构的特性进行分析,并对提出的几种方案的凸轮型线进行的对比,最终获得了LJ750发动机配气机构的有关数据。
对比的结果最终还是采用了和CBR900发动机相同的配气结构,因此涉及到配气结构的许多动力学计算在这里都略去了。
本章主要介绍了多项动力凸轮设计的基本原理及其准则,并对CBR900发动机的凸轮型线进行了细致的分析和研究,在此基础上提出了LJ750发动机的配气机构。
一、凸轮型线类型的选择配气机构是发动机的一个重要系统,其设计好坏对发动机的性能、可靠性和寿命有极大的影响。
其中凸轮型线设计是配气机构设计中最为关键的部分,在确定了系统参数后,重要的问题是根据发动机的性能和用途,正确选择凸轮型线类型及凸轮参数。
凸轮型线有多种,如复合正弦,复合摆线,低次方,高次方,多项动力,谐波凸轮等。
19.应⽤GT-POWER设计发动机⽓门升程应⽤GT-POWER设计发动机⽓门升程张⼩燕蒲运平(长安汽车⼯程研究院重庆 401120)摘要:本⽂利⽤GT-POWER计算,介绍了车⽤发动机的⽓门升程设计问题。
设计中,主要从发动机性能⽅⾯考虑⽓门升程的设计要求,利⽤DOE的⽅法,考察了⽓门开启持续⾓及⽓门升程对发动机充量系数的影响。
关键词:GT-POWER ⽓门升程⽓门开启持续⾓充量系数Design of Valve Lift of Automotive Engine by GT-POWERAbstract:A GT-POWER model is used to show the details of valve lift designed of automotive engine in this paper.The effection of valve opening and valve lift on engine’s volumetric coefficiency is studied by DOE(design of experiment),considering primarily the required profile for valve lift on performance of engine.Key words: GT-POWER Valve Lift Valve Opening Duration Volumetric Coefficiency1、介绍对于传统凸轮驱动的配⽓机构,由于受到结构的限制,⽓门不可能瞬间开启到最⼤升程的位置,其升程特性只能是连续变化的,这在换⽓过程中造成很⼤的流动损失,有损于发动机的动⼒性。
⽓门升程对⽓道流量系数有⾮常重要的影响,⽓道流量系数随⽓门升程的增⼤⽽增⼤,如图1所⽰,并且可以看出,⽓门升程开启达⼀定⾼度后,流量系数⼏乎不再增加。
因此合理设计⽓门的运动规律,也就是⽓门升程(Valve Lift或者Valve Profile),对提⾼充量系数Φc,改善发动机换⽓质量是相当重要的。
摘要 (2)一、可变气门正时技术 (3)(一)、可变气门正时系统的原理 (3)1、可变配气相位调整原理 (4)2、可变配气相位技术条件 (5)(二)、可变气门正时技术的现状 (5)(三)、可变气门正时技术的发展趋势 (6)二、国内外可变气门配气机构的现状和发展趋势 (7)(一)、可变配气机构分类 (7)(二)、可变气门技术的发展现状 (7)三、可变配气相位技术研究意义 (8)三、连续可变配气凸轮轴设计浅析 (9)(一)、连续可变凸轮轴作用 (9)(二)、连续可变配气凸轮轴的工作原理 (9)(三)、连续可变配气凸轮轴与传统可变配气技术凸轮轴优缺点比较 (9)(四)、可变配气相位技术条件 (11)四、可变气门正时技术的发展趋势 (11)参考文献 (13)本文介绍了从进气晚关角及进排气的动态效应几方面着手,不断改进发动机的配气相位以及进排气系统,使发动机的实际性能曲线逐步接近计算机仿真曲线。
配气相位、进气门间隙、排气门间隙、转速、负荷五个调整参数之间是相互影响的。
通过在配气机构多刚体模型中引入柔性体,描述了配气机构的动力学性能;建立了柔性体气门弹簧,分析了气门弹簧动刚度的非线性行为,并且依据模态技术计算得到其动态应力;该方法为优化设计配气机构等机械产品及对其进行疲劳性能研究提供了依据。
该仪器可检测各种汽、柴油发动机的启动性能、高压点火性能、燃油喷射性能、充电性能、动力性能、配气相位、发动机异响震动分析等30余种技术参数,并分析故障产生的原因、检测过程中,可随时显示各种波形及技术参数和结果并可随机打印,该仪器内存有一百多种国内外发动机技术参数,内容十分丰富,随时可以与检测结果对比。
Passat B5轿车有4缸和6缸两种发动机,4缸机有4G54与4G64两种型号,6缸机型号为6G72,其配气机构均采用顶置凸轮轴式配气机构。
介绍了气门间隙自动调整器的结构、工作原理,以及其维护与保养。
目前,汽车工业的发展正在面临着两个主要问题——能源的枯竭与环境的污染。
・文摘・顶置凸轮配气机构气门升程的精确计算苏军申屠淼(江苏理工大学力学系 , 镇江 212013为使发动机布置得更紧凑 , 提高配气机构的刚性与减轻运动件的质量 , 以适应高转速的要求 , 许多现代的四冲程发动机都采用顶置凸轮配气机构 , 其中一些汽车、摩托车发动机的顶置凸轮配气机构的结构形式如图 1所示。
图 1顶置凸轮摇臂机构现对以上典型结构的气门升程的计算方法与公式进行了推演 , 以供应用参考。
假定摇臂为绝对刚体。
(1 将凸轮升程表换算为极坐标升程表 h ( 。
(2 建立 2个坐标系 , 一个为定坐标系 x Oy ; 另一个为随凸轮转动的动坐标系 x Oy (详见图 2 。
两坐标系原点 O 均在凸轮轴心 ; 横坐标 Ox垂直于凸轮轴图 2顶置凸轮配气机构简图线 , 平行于气缸盖底平面 ; 横坐标 Ox 垂直于凸轮轴线 , 并令其通过凸轮升程廓线过渡段的起点 ; 认定摇臂初始位置 , 即 =0°。
此时 Ox 与 Ox 的夹角为 : 0=180-arctg1-x 1-arctg 1R +S 1(1 式中 :x 1和 y 1为摇臂中心 O 1在定坐标系中的坐标 ; L 1为 O 1点至 M 点的距离 ; M 点为摇臂零件图中水平中心线与 OP 0的交点 ; S 1为 OP 0线上 R 3圆弧面至 M 点的弦高 ; R 为凸轮基圆半径。
(3 计算出 R 3圆弧的圆心 P 0在定坐标系 x Oy 中的坐标位置。
x P 0=(R +R 3 ・ cos 0y P 0=(R +R 3 ・ sin 0(2式中 :R 3为摇臂凸轮侧臂圆弧的半径。
当凸轮转过角度时 , 凸轮廓线与摇臂 R 3圆弧面的接触点在 C (x C , y C 点。
在凸轮转动中 , R 3圆弧面上的接触点 C 在 C L 和 C R 之间变动 ; 而 C 点在凸轮廓线上则是在 x Oy 坐标系中逆旋转方向而变动。
令 OC 与 Ox 的夹角为。
凸轮推杆的运动规律可以通过一些例题来理解和掌握。
以下是一个关于凸轮推杆升程的例题:
例题:已知一偏置直动尖顶推杆盘形凸轮机构,其基圆半径为10mm,凸轮转速为10 r/min,要求推杆做一次往返运动,求推杆的升程和降程。
解:根据题目,凸轮转速为10 r/min,所以每分钟转过的角度为360°×10 = 3600°。
推杆做一次往返运动,即完成一个周期,需要转过的角度为360°。
因此,推杆完成一个周期所需的时间为3600°/ 360°= 10 min。
推杆的升程和降程可以根据基圆半径和推杆长度计算得出。
假设推杆长度为L,则升程为L -基圆半径,降程为L + 基圆半径。
根据题目,基圆半径为10mm,假设推杆长度L = 30mm,则推杆的升程和降程分别为20mm和40mm。
% ******** 偏置移动从动件盘形凸轮设计绘图和运动分析******** disp ' ######## 已知条件########'disp ' 凸轮作逆时针方向转动,从动件偏置在凸轮轴心的右边'disp ' 从动件在推程作等加速/等减速运动,在回程作余弦加速度运动' % 基圆半径;滚子半径;从动件偏距;从动件升程rb=40;rt=10;e=15;h=50;% 推程运动角;远休止角;回程运动角;推程许用压力角;凸轮转速ft=100;fs=60;fh=90;alpha_p=35;n=200;% 角度和弧度转换系数;机构尺度hd=pi/180;du=180/pi;se=sqrt(rb^2-e^2);w=n*pi/30; omega=w*du; % 凸轮角速度(°/s)fprintf(' 基圆半径rb = %3.4f mm \n',rb)fprintf(' 滚子半径rt = %3.4f mm \n',rt)fprintf(' 推杆偏距 e = %3.4f mm \n',e)fprintf(' 推程升程h = %3.4f mm \n',h)fprintf(' 推程运动角ft = %3.4f 度\n',ft)fprintf(' 远休止角fs = %3.4f 度\n',fs)fprintf(' 回程运动角fh = %3.4f 度\n',fh)fprintf(' 推程许用压力角alpha_p = %3.4f 度\n',alpha_p) fprintf(' 凸轮转速n = %3.4f r/min \n',n) fprintf(' 凸轮角速度(弧度) w = %3.4f rad/s \n',w)fprintf(' 凸轮角速度(度) omega = %3.4f 度/s \n',omega) disp ' 'disp ' @@@@@@ 计算过程和输出结果@@@@@@' disp ' '% (1)---校核凸轮机构的压力角和轮廓曲率半径'disp ' *** 计算凸轮理论轮廓的压力角和曲率半径***'disp ' 1 推程(等加速/等减速运动)'for f=1:ftif f<=ft/2s(f)=2*h*f^2/ft^2;s=s(f); % 等加速-位移方程ds(f)=4*h*f*hd/(ft*hd)^2;ds=ds(f);d2s(f)=4*h/(ft*hd)^2;d2s=d2s(f);vt(f)=4*h*omega*f/ft^2; % 等加速-速度方程elses(f)=h-2*h*(ft-f)^2/ft^2;s=s(f); % 等减速-位移方程ds(f)=4*h*(ft-f)*hd/(ft*hd)^2;ds=ds(f);d2s(f)=-4*h/(ft*hd)^2;d2s=d2s(f);vt(f)=4*h*omega*(ft-f)/ft^2; % 等减速-速度方程endalpha_t(f)=atan(abs(ds-e)/(se+s)); % 推程压力角(弧度)alpha_td(f)=alpha_t(f)*du; % 推程压力角(度)pt1=((se+s)^2+(ds-e)^2)^1.5;pt2=abs((se+s)*(d2s-se-s)-(ds-e)*(2*ds-e));rho_t(f)=pt1/pt2; % 推程曲率半径st(f)=s;endalpha_tm=max(alpha_td);fprintf(' 推程最大压力角alpha_tm = %3.4f 度\n',alpha_tm)for f=1:ftif alpha_td(f)==alpha_tm;ftm=f;break;endendfprintf (' 对应的位置角ftm = %3.4f 度\n',ftm)if alpha_tm>alpha_pfprintf(' * 凸轮推程压力角超过许用值,需要增大基圆!\n')endrho_tn = min(rho_t);fprintf (' 最小曲率半径rho_tn = %3.4f mm\n',rho_tn)for f=1:ftif rho_t(f)==rho_tn;ftn=f;break;endendfprintf(' 对应的位置角ftn = %3.4f 度\n',ftn)if rho_tn<rt+5fprintf(' * 凸轮推程轮廓曲率半径小于许用值,需要增大基圆或减小滚子!\n') enddisp ' 2 回程(余弦加速度运动-简谐运动)'d1=ft+fs;d2=ft+fs+fh; % 回程运动角范围for f=d1:d2k=f-d1;s(f)=0.5*h*(1+cos(pi*k/fh));s=s(f); % 简谐运动-位移方程ds(f)=-0.5*pi*h*sin(pi*k/fh)/(fh*hd);ds=ds(f);d2s(f)=-0.5*pi^2*h*cos(pi*k/fh)/(fh*hd)^2;d2s=d2s(f);alpha_h(f)=atan(abs(ds+e)/(se+s)); % 回程压力角(弧度)alpha_hd(f)=alpha_h(f)*du; % 回程压力角(度)ph1=((se+s)^2+(ds-e)^2)^1.5;ph2=abs((se+s)*(d2s-se-s)-(ds-e)*(2*ds-e));rho_h(f)=ph1/ph2; % 回程曲率半径sh(f)=s;vh(f)=-0.5*pi*h*omega*sin(pi*f/fh)/fh; % 简谐运动-速度方程ah(f)=-0.5*pi^2*h*omega^2*cos(pi*f/fh)/fh^2; % 简谐运动-加速度方程endalpha_hm = max(alpha_hd(d1:d2));fprintf(' 回程最大压力角alpha_hm = %3.4f 度\n',alpha_hm)for f=d1:d2if alpha_hd(f)==alpha_hm;fhm=f;break;endendfprintf(' 对应的位置角fhm = %3.4f 度\n',fhm)rho_hn=min(rho_h(d1:d2));fprintf(' 最小曲率半径rho_hn = %3.4f mm\n',rho_hn)for f=d1:d2if rho_h(f)==rho_hn;fhn=f;break;endendfprintf(' 对应的位置角fhn = %3.4f 度\n',fhn)if rho_hn<rt+5fprintf(' * 凸轮回程轮廓曲率半径小于许用值,需要增大基圆或减小滚子!\n') enddisp ' '% (2)---计算凸轮机构的从动件运动参数'disp ' *** 计算凸轮机构从动件的运动参数***'disp ' 1 推程(等加速/等减速运动)'disp ' 凸轮转角位移s(mm) 速度v(mm/s)'for f=10:10:ftydcs_t=[f st(f) vt(f)];disp(ydcs_t)endat_1=4*h*omega^2/ft^2;at_2=-4*h*omega^2/ft^2;fprintf(' 等加速上升的加速度at_1 = %3.4f (mm/s^2) \n',at_1)fprintf(' 等减速上升的加速度at_2 = %3.4f (mm/s^2) \n',at_2)disp ' 2 回程(余弦加速度运动-简谐运动)'disp ' 凸轮转角位移s(mm) 速度v(mm/s) 加速度a(mm/s^2)'for f=d1:10:d2ydcs_h=[f sh(f) vh(f) ah(f)];disp(ydcs_h)end% (3)---绘制凸轮机构的从动件运动线图figure(1);subplot(3,2,1) % 推程位移线图f=1:ft;plot(f,st);xlabel ('凸轮转角\it \phi / \rm( °)')ylabel ('\it s / \rm(mm)')title('从动件推程位移线图');subplot(3,2,2) % 回程位移线图f=d1:d2;plot(f,sh(d1:d2));xlabel ('凸轮转角\it \phi / \rm( °)')ylabel ('\it s / \rm(mm)')title('从动件回程位移线图');subplot(3,2,3) % 推程速度线图f=1:ft;plot(f,vt);xlabel ('凸轮转角\it \phi / \rm( °)')ylabel ('\it v / \rm(mm/s)')title('从动件推程速度线图');subplot(3,2,4) % 回程速度线图f=d1:d2;plot(f,-vh(d1:d2));xlabel ('凸轮转角\it \phi / \rm( °)')ylabel ('\it v / \rm(mm/s)')title('从动件回程速度线图');subplot(3,2,5) % 推程加速度线图line([0,ft/2],[at_1,at_1]);line([ft/2,ft/2],[at_1,at_2]); % 等加速等减速之间的突变垂线line([ft/2,ft],[at_2,at_2]);xlabel ('凸轮转角\it \phi / \rm( °)')ylabel ('\it a / \rm(mm/s^2)')title('从动件推程加速度线图');subplot(3,2,6) % 回程加速度线图f=d1:d2;plot(f,-ah(d1:d2));xlabel ('凸轮转角\it \phi / \rm( °)')ylabel ('\it a / \rm(mm/s^2)')title('从动件回程加速度线图');disp ' '% (4)---计算凸轮理论廓线与实际廓线的直角坐标和向径'disp ' ****** 凸轮理论轮廓与实际轮廓的直角坐标******'nd=360;for f=1:ndif f<=ft/2 % 等加速运动s(f)=2*h*f^2/ft^2;s=s(f);ds(f)=4*h*f*hd/(ft*hd)^2;ds=ds(f);elseif f>ft/2 & f<=ft % 等减速运动s(f)=h-2*h*(ft-f)^2/ft^2;s=s(f);ds(f)=4*h*(ft-f)*hd/(ft*hd)^2;ds=ds(f);elseif f>ft & f<=d1 % 远休止角s=h;ds=0;elseif f>d1 & f<=d2 % 简谐运动k=f-d1;s(f)=0.5*h*(1+cos(pi*k/fh));s=s(f);ds(f)=-0.5*pi*h*sin(pi*k/fh)/(fh*hd);ds=ds(f);elseif f>d2 & f<=nds=0;ds=0;endxx(f)=(se+s)*sin(f*hd)+e*cos(f*hd);x=xx(f); % 理论轮廓横坐标yy(f)=(se+s)*cos(f*hd)-e*sin(f*hd);y=yy(f); % 理论轮廓纵坐标dx(f)=(ds-e)*sin(f*hd)+(se+s)*cos(f*hd);dx=dx(f);dy(f)=(ds-e)*cos(f*hd)-(se+s)*sin(f*hd);dy=dy(f);xp(f)=x+rt*dy/sqrt(dx^2+dy^2);xxp=xp(f); % 实际轮廓横坐标yp(f)=y-rt*dx/sqrt(dx^2+dy^2);yyp=yp(f); % 实际轮廓纵坐标r(f)=sqrt(x^2+y^2); % 理论轮廓向径rp(f)=sqrt(xxp^2+yyp^2); % 实际轮廓向径enddisp ' 1 推程(等加速/等减速运动)'disp ' 凸轮转角理论x 理论y 实际x 实际y'for f=10:10:ftnu=[f xx(f) yy(f) xp(f) yp(f)];disp(nu)enddisp ' 2 回程(余弦加速度运动)'disp ' 凸轮转角理论x 理论y 实际x 实际y'for f=d1:10:d2nu=[f xx(f) yy(f) xp(f) yp(f)];disp(nu)enddisp '*** 凸轮理论轮廓与实际轮廓的向径***'disp ' 1 推程(等加速/等减速运动)'disp ' 凸轮转角理论r 实际r'for f=10:10:ftnu=[f r(f) rp(f)];disp(nu)enddisp ' 'disp ' 2 回程(余弦加速度运动)'for f=d1:10:d2nu=[f r(f) rp(f)];disp(nu)end% (5)---绘制凸轮的理论轮廓和实际轮廓figure(2);plot(xx,yy,'r-.') % 理论轮廓(红色,点划线)axis ([-(rb+h-10) (rb+h+10) -(rb+h+10) (rb+rt+10)]) % 横轴和纵轴的下限和上限axis equal % 横轴和纵轴的尺度比例相同text(rb+h+3,0,'X') % 标注横轴text(0,rb+rt+3,'Y') % 标注纵轴text(-5,5,'O') % 标注直角坐标系原点title('偏置移动从动件盘形凸轮轮廓') % 标注图形标题hold on; % 保持图形plot([-(rb+h) (rb+h)],[0 0],'k') % 横轴(黑色)plot([0 0],[-(rb+h) (rb+rt)],'k') % 纵轴(黑色)plot([e e],[0 (rb+rt)],'k--') % 初始偏置位置(黑色,虚线)ct=linspace(0,2*pi); % 画圆的极角变化范围plot(rb*cos(ct),rb*sin(ct),'g') % 基圆(绿色)plot(e*cos(ct),e*sin(ct),'c--') % 偏距圆(青色,虚线)plot(e + rt*cos(ct),se + rt*sin(ct),'m') % 滚子圆(品红色)plot(xp,yp,'b') % 实际轮廓(蓝色)******** 偏置移动从动件盘形凸轮设计绘图和运动分析********######## 已知条件########凸轮作逆时针方向转动,从动件偏置在凸轮轴心的右边从动件在推程作等加速/等减速运动,在回程作余弦加速度运动基圆半径rb = 40.0000 mm滚子半径rt = 10.0000 mm推杆偏距 e = 15.0000 mm推程升程h = 50.0000 mm推程运动角ft = 100.0000 度远休止角fs = 60.0000 度回程运动角fh = 90.0000 度推程许用压力角alpha_p = 35.0000 度凸轮转速n = 200.0000 r/min凸轮角速度(弧度) w = 20.9440 rad/s凸轮角速度(度) omega = 1200.0000 度/s@@@@@@ 计算过程和输出结果@@@@@@*** 计算凸轮理论轮廓的压力角和曲率半径***1 推程(等加速/等减速运动)推程最大压力角alpha_tm = 34.2666 度对应的位置角ftm = 50.0000 度最小曲率半径rho_tn = 35.2303 mm对应的位置角ftn = 51.0000 度2 回程(余弦加速度运动-简谐运动)回程最大压力角alpha_hm = 30.9248 度对应的位置角fhm = 213.0000 度最小曲率半径rho_hn = 30.3591 mm对应的位置角fhn = 250.0000 度*** 计算凸轮机构从动件的运动参数***1 推程(等加速/等减速运动)凸轮转角位移s(mm) 速度v(mm/s)10 1 24020 4 48030 9 72040 16 96050 25 120060 34 96070 41 72080 46 48090 49 240100 50 0等加速上升的加速度at_1 = 28800.0000 (mm/s^2)等减速上升的加速度at_2 = -28800.0000 (mm/s^2)2 回程(余弦加速度运动-简谐运动)凸轮转角位移s(mm) 速度v(mm/s) 加速度a(mm/s^2) 160 50 673 -33602170 48 358 -41220180 44 0 -43865190 37 -358 -41220200 29 -673 -33602210 21 -907 -21932220 12.5 -1031.3 -7617.1230 5.8 -1031.3 7617.1240 2 -907 21932250 0 -673 33602****** 凸轮理论轮廓与实际轮廓的直角坐标******1 推程(等加速/等减速运动)凸轮转角理论x 理论y 实际x 实际y 10.0000 21.3848 34.8977 18.7440 25.2527 20.0000 28.1459 33.4732 26.5660 23.5988 30.0000 36.0309 32.4073 34.7788 22.4860 40.0000 45.6105 31.0206 43.9004 21.1679 50.0000 57.1986 28.4142 54.4870 18.7889 60.0000 69.0579 22.5501 63.1030 14.5165 70.0000 78.5024 12.6099 70.2060 7.0270 80.0000 84.4235 -0.3453 74.7846 -3.008390.0000 86.0810 -15.0000 76.0894 -14.5890 100.0000 83.1533 -29.8936 73.7429 -26.51052 回程(余弦加速度运动)凸轮转角理论x 理论y 实际x 实际y 160.0000 15.6881 -86.9597 13.9127 -77.1185 170.0000 0.0875 -86.8780 1.9206 -77.0474 180.0000 -15.0000 -81.2321 -9.9808 -72.5829 190.0000 -27.7230 -70.8432 -20.2897 -64.1539 200.0000 -36.8131 -57.2861 -27.8219 -52.9092 210.0000 -41.8603 -42.5041 -32.0770 -40.4336 220.0000 -43.3607 -28.3394 -33.3609 -28.2733 230.0000 -42.5280 -16.1041 -32.6176 -17.4398 240.0000 -40.9188 -6.3040 -31.0634 -7.9985 250.0000 -39.9750 1.4129 -29.9813 1.0597*** 凸轮理论轮廓与实际轮廓的向径***1 推程(等加速/等减速运动)凸轮转角理论r 实际r10.0000 40.9287 31.449020.0000 43.7338 35.533930.0000 48.4609 41.414840.0000 55.1597 48.737350.0000 63.8674 57.635560.0000 72.6465 64.751270.0000 79.5088 70.556880.0000 84.4242 74.845190.0000 87.3781 77.4754100.0000 88.3634 78.36342 回程(余弦加速度运动)160.0000 88.3634 78.3634170.0000 86.8780 77.0714180.0000 82.6054 73.2660190.0000 76.0745 67.2859200.0000 68.0948 59.7783210.0000 59.6564 51.6121220.0000 51.8003 43.7302230.0000 45.4750 36.9872240.0000 41.4015 32.0766250.0000 40.0000 30.0000。
・文摘・
顶置凸轮配气机构气门升程的精确计算
苏军申屠淼
(江苏理工大学力学系 , 镇江 212013
为使发动机布置得更紧凑 , 提高配气机构的刚性与减轻运动件的质量 , 以适应高转速的要求 , 许多现代的四冲程发动机都采用顶置凸轮配气机构 , 其中一些汽车、摩托车发动机的顶置凸轮配气机构的结构形式如图 1
所示。
图 1顶置凸轮摇臂机构
现对以上典型结构的气门升程的计算方法与公式进行了推演 , 以供应用参考。
假定摇臂为绝对刚体。
(1 将凸轮升程表换算为极坐标升程表 h ( 。
(2 建立 2个坐标系 , 一个为定坐标系 x Oy ; 另一个为随凸轮转动的动坐标系 x Oy (详见图 2 。
两坐标系原点 O 均在凸轮轴心 ; 横坐标 Ox
垂直于凸轮轴
图 2顶置凸轮配气机构简图
线 , 平行于气缸盖底平面 ; 横坐标 Ox 垂直于凸轮轴
线 , 并令其通过凸轮升程廓线过渡段的起点 ; 认定摇臂初始位置 , 即 =0°。
此时 Ox 与 Ox 的夹角为 : 0=180-arctg
1-x 1-arctg 1
R +S 1
(1 式中 :x 1和 y 1为摇臂中心 O 1在定坐标系中的坐标 ; L 1为 O 1点至 M 点的距离 ; M 点为摇臂零件图中水平中心线与 OP 0的交点 ; S 1为 OP 0线上 R 3圆弧面至 M 点的弦高 ; R 为凸轮基圆半径。
(3 计算出 R 3圆弧的圆心 P 0在定坐标系 x Oy 中的坐标位置。
x P 0=(R +R 3 ・ cos 0
y P 0=(R +R 3 ・ sin 0
(2
式中 :R 3为摇臂凸轮侧臂圆弧的半径。
当凸轮转过角度时 , 凸轮廓线与摇臂 R 3圆弧面
的接触点在 C (x C , y C 点。
在凸轮转动中 , R 3圆弧面上的接触点 C 在 C L 和 C R 之间变动 ; 而 C 点在凸轮廓线上则是在 x Oy 坐标系中逆旋转方向而变动。
令 OC 与 Ox 的夹角为。
(4 气门升程计算。
计算的第 1步为由设定的来计算。
当在某一接触点 C 位置时 , 对应角度是。
此时 C 点在动坐标系中的坐标为 :x C =[R +h ( ]・ cos
y
C =[R +h ( ]・ sin
(3
式中 :h ( 为凸轮在角时的极坐标升程。
计算 C 点前后 -1和 +1点的坐标 , 得 (x
K -1, y
K -1 和 (x
K +1, y
K +1 。
据此用差商法即可求得 C 处凸轮升程廓线的切向斜率 k ( 和法向斜率k 1( :
k ( =
K +1
K -1
x K +1-x K -1k 1( =1/k (
(4
则 R 3圆弧的圆心 P 在动坐标系 x Oy 中的坐标为 :x P =x
C +R 3・ co s[arctg k 1( ]
y P =y
C +R 3・ sin[ar ctg k 1( ]
(5 在动坐标中某一位置时 , 摇臂中心 O 1点可视为以 O 为圆心 , 以 OO 1(即 R 1 为半径的圆周与以 P 为圆心 , 以 O 1P (即 R 2 为半径的圆周相交的交点 , 这 2
个圆第 17卷 (1999 第 2期内燃机学报Transactions of CSICE
Vo l . 17(1999 N o . 2
x 2+y 2=R 21
(x -x P 2+(y -y P 2=R 22
(6 将式 (6 整理化简 , 得 :
(1+A 2 y 2+2A By +(B 2-R 21 =0
(7
式中 :A = P x P
; B =2221 2P 2P
2x P。
求解式 (7 , 其根为 :y 1
=
22222
1(1+A 2
x 1=-(B +Ay
1
(8
上式中的“ ±” 号可通过判别取一个正确值。
于是可得 O 1点在 x Oy 坐标系中的坐标 , 进而求出 O 1O 与 Ox 的夹角 :
∠ O 1Ox =arctg y 1
x 1
而在固定坐标系 x Oy 中应有 :∠ O 1Ox =arctg 1
x 1,
由此可得 , 在某一角度时 : =∠ O 1Ox -∠ O 1Ox 。
于是第 1步计算完成 :
= 0- =180-ar ctg 1R +S 1-arctg 1x 1+2arctg 1
x 1
第 2步计算摇臂的摆动角 ( , 这可由 R 3的圆心 P 与初始位置 P 0求得。
P 点在定坐标系 x Oy 中的 (5 求出的 (x
P , y
P 经坐标旋转得到 :
x P =x P ・ cos -y P sin
y P =x P ・ sin -y
P cos
(10 于是有 :
( =arctg(
P 1x P -x 1 -arctg(P 01
x P 0-x 1
(11 第 3步计算气门升程 , 令摇臂气门侧臂 R 4圆弧的圆心为 Q , 其定坐标系中的坐标为 :
x Q =x 1-R 5・ cos
y Q =y 1-R 5・ sin
(12
式中 :R 5为 O 1点至 Q 的距离 , O 1Q 是定值 ; 角是基准水平线与 O 1Q 的夹角 , 当 =0时 , = 0, 当摇臂摆动时 , = 0+ ( 。
令 R 4圆弧与气门头的接触点为 D , 其坐标为 :
x D =x Q +R 4・ sin
y D =y Q -R 4・ cos
(13
=0°时 , D 点的初始位置 D 0为 :
x D 0=x Q
0+R 4・ sin
y D 0=y Q 0+R 4・ cos
(14
摇臂摆动中 D 点的变化距离为 :
L =DD 0=(y D -y D 0 2
+(x D -x D 0
2
气门的实际升程 h q ( 是 L 在气门中心线上的投影 ;
=上面 3个式中的角是气门中心线与垂直线的夹角。
按照上述计算过程编制了计算机程序 , 程序中可以直接应用凸轮图纸上的滚子或平面凸轮升程表 , 自动转换为凸轮极坐标转角 , 最终计算出凸轮升程与气门升程表之间的对应关系 , 且可以计算出任意位置的实际摇臂比。
上述计算是一种完全精确的计算方法。
(5 计算实例。
现对该类型气门进行了计算 , 表 1为计算所得的进气门某些特定点的数值对比。
由表 1可知 :
表 1计算实例中的进气门特定点对比表
凸轮轴转角 /°凸轮和摇臂实际接触点 /°滚子凸轮升程 /mm 气门实际
升程 /mm 备注 000. 00. 01. 1起始点 22230. 08630. 01. 115气门开启点 39. 5640. 331. 91. 52最左接触点
81885. 68736. 841. 34顶点 1471310. 6151. 1870. 98最右接触点 1821810. 080. 01. 1气门关闭点
198
198
0. 0
0. 0
1. 1
终点
(a 凸轮的转角与实际接触点的角度 , 除始末两点外 , 相差较大 , 因此形成了实际升程与摇臂比的变化都很大。
如凸轮转角为 39. 5°时 , 按滚子凸轮计算的升程是 0. 3m m , 而凸轮与摇臂的实际接触点已在 64°, 此时的实际升程已是 1. 25m m ; 摇臂比为 1. 52, 气门的实际升程为 1. 9m m 。
故“ 放大” 了气门的时间断面。
图 3为凸轮升程与气门升程对照图。
图 3计算结果
(b 上述凸轮升程与摇臂比变化最大的位置是在凸轮转角 30°~60°之间 , 故能在很短的几度内使气门迅速开至较大值 , 形成进气门速开。
对于反方向装的
排气摇臂 , 正好形成速关。
(c 经多次计算得图 2中摇臂在凸轮侧的长度 h 2, 对于已确定的 R 3等
其它结构参数有一个最佳值。
如偏离最佳值 0. 2mm 以上 , 则凸轮和摇臂的接触面之间将产生较大的滑移 , 同时升程曲线产生鼓出或下凹现象 , 从而使气门产生跳动 , 大大加剧了凸轮与摇臂接触面的磨损或损伤。
・
203・ 1999年 4月苏军等 :顶置凸轮配气机构气门升程的精确计算。