凸轮机构的弹性动力学

-.z.凸轮机构的弹性动力学分析〔附MATLAB代码〕【问题】一凸轮系统，欲使其考虑弹性因素后从动件的真实运动规律按照余弦加速度运动规律运动，建立该凸轮系统的弹性动力学模型，分析其未考虑弹性因素时从动件的运动规律，并绘制出从动件的理论运动规律及考虑弹性因素后的真实运动规律。凸轮系统的运动及动力参数自定。程序代码需提供电子版，并说明运行环境。【解答】建立动力学模型取图1所示的凸轮机构为研究对象，图2为其所对应的动力学模型。图1：凸轮机构运动简图图2：凸轮机构的动力学模型为使得问题简化，力学模型中忽略了凸轮轴的扭转变形、弯曲变形以及回位弹簧的阻尼作用。图2中为系统等效弹簧的刚度，为凸轮机构从动组件的阻尼系数，为回位弹簧的刚度，为回位弹簧的预紧力，为凸轮机构在从动件侧的当量质量，为与凸轮廓线有关的等效凸轮升程〔图中所示的凸轮并非真正的凸轮，其廓线对应的升程与真实凸轮廓线对应的升程具备关系，其中为摇臂比。因为与仅相差一个比例系数，为了便于表达，后文将只注重分析与从动件输出的关系，而不再专门区别与的差异〕，为从动件的实际升程。建立动力学方程该机构的自由度为1，利用牛顿第二定律建立运动微分方程：〔式1〕设凸轮转动的角速度为，它与时间微分、凸轮转角微分具有关系：〔式2〕将〔式2〕代入〔式1〕并整理可得：〔式3〕微分方程〔式3〕有两层含义：=1\*GB3①假设从动件的真实运动规律，可求解出凸轮在高速运转条件下考虑弹性变形影响的理论轮廓；=2\*GB3②假设凸轮廓线，可求解考虑弹性变形的从动件的动力学响应。运动方程的求解〔一〕凸轮轮廓的设计条件如下：，凸轮的转速，，，，；为防止余弦加速度运动产生的冲击，取凸轮的推程运动角和回程运动角均为，远休止角和近休止角均为，从动件的最大升程。根据条件，可以确定从动件的位移方程将上式代入〔式3〕可得：〔式4〕由于〔式4〕对应的常微分方程难以求出解析解，这里利用MATLAB求解出其数值解并与位移方程比拟如下列图：假设需要考虑机构的弹性变形，凸轮的轮廓应按照上图中的红色曲线进展设计。由廓线求解从动件的真实运动规律由于系统的弹性变形，从动件的真实运动与等效凸轮升程不再相等，当然，从动件的真实运动速度、加速度与对应的理论值也不相等。由于弹簧的预紧力为常数，它只影响系统振动的初始平衡位置，故在分析从动件的运动规律时不再考虑，从而〔式1〕被简化为：〔式5〕根据振动理论，系统自由振动的固有频率其中阻尼比代入相关数值计算可得，如果从动件按照前述的余弦加速度规律运动，则将上式代入〔式5〕可得：〔式6〕利用MATLAB求解〔式6〕的数值解，结果如下〔图中的〕：从上图中的位移、速度、加速度分析看，当考虑到系统的弹性和阻尼后，工作端的运动规律发生改变，。只有当很大时，也就是说当系统的固有频率很大时〔刚度大〕，而且凸轮的角速度很小时，才接近值。程序附录：%filename:tulun.m%function:cooperatewithjisuanlunkuo.mtocalculatethecurveofthecamfunctiond*=tulun(a,*)%%凸轮机构的构造参数h=6.2*10^(-3);%升程m=0.08459;%凸轮机构的当量质量n=1200;w=n/60*2*pi;%凸轮转速及角速度c=55.7148;%阻尼k1=10400;%回位弹簧刚度k=3194800;%系统等效弹簧刚度F=400;%回位弹簧预紧力%%理论轮廓的微分方程d*=(h/2*(m*w^2-k-k1)*cos(a)+h/2*c*w*sin(a)+0.5*(k+k1)*h+F-k**)/(c*w);%filename:jisuanlunkuo.m%function:cooperatewithtulun.mtocalculatethecurveofthecam%%计算考虑弹性变形的凸轮的理论轮廓曲线[a,*]=ode45('tulun',[02*pi],0);plot(a,*,'r');holdon;%%作出未考虑弹性变形的凸轮的理论轮廓曲线a1=linspace(0,2*pi);h=6.2*10^(-3);y=h/2*(1-cos(a1));plot(a1,y,'-b');%%添加标注gridon;legend('考虑弹性变形后凸轮的理论轮廓','未考虑弹性变形的凸轮理论轮廓');*label('凸轮的转角/rad');ylabel('升程/m');%filename:yundongfen*i.m%function:cooperatewithfen*i.mtocalculatethedisplacement,velocity%andaccelerationofthecammechanismwhentheelasticdeformationof%thecamisnotconsideredindesignfunctiondy=yundongfen*i(a,y)%%凸轮机构的构造参数h=6.2*10^(-3);m=0.08459;c=55.7148;k1=10400;w=1537;%不同情况下将凸轮的角速度值分别设为1537,2049,3074k=3194800;%%运动分析微分方程dy=zeros(2,1);dy(1)=y(2);dy(2)=(c*w*h/2*sin(a)+k*h/2*(1-cos(a))-(k+k1)*y(1)-c*w*y(1))/(m*w*w);%filename:fen*i.m%function:cooperatewithyundongfen*i.mtocalculatethedisplacement,velocity%andaccelerationofthecammechanismwhentheelasticdeformationof%thecamisnotconsideredindesign%%绘制从动件的位移曲线[a,y]=ode45('yundongfen*i',[02*pi],[00]);n=length(a);figure(1);plot(a,y(:,1),'--b');%不同情况下将线型分别设为--b，--g，--kholdon;%%绘制从动件的速度曲线figure(2);a1=zeros(n-1,1);fori=1:(n-1)a1(i)=a(i);endy1=diff(y(:,1))./diff(a);plot(a1,y1,'--b');%不同情况下将线型分别设为--b，--g，--kholdon;%%绘制从动件的加速度曲线figure(3);a2=zeros(n-2,1);fori=2:(n-1)a2(i-1)=a1(i);endy2=diff(y1(:,1))./diff(a1);plot(a2,y2,'--b');%不同情况下将线型分别设为--b，--g，--kholdon;%fliename:zuihou.m%function:givetheidealdisplacement,velocityandaccelerationofthe%cammechanism%%绘制理想位移曲线a0=linspace(0,2*pi);n=length(a0);y0=3.1*10^(-3)*(1-cos(a0));figure(1);plot(a0,y0,'r');holdon;%%绘制理想速度曲线figure(2);a01=zeros(1,n-1);fori=1:(n-1)a01(i)=a0(i);endy01=diff(y0)./diff(a0);plot(a01,y01,'r');holdon;%%绘制理想加速度曲线figure(3);a02=zeros(1,n-2);fori=2:(n-1)a02(i-1)=a01(i);endy02=diff(y01)./diff(a01);plot(a02,y02,'r');holdon;%%添加标注figure(1);*label('凸轮的转角/rad');ylabel('位移/m');legend('n=4','n=3','n=2','标准输出');figure(2);*label('凸轮的转角/rad');ylabel('速度/〔m/s〕');legend('n=4','n=3','n=2','标准输出');figure(3);*label('凸轮的转角/rad');ylabel('加速度/〔m/s^2〕');legend('n=4','n=3','n=2','标准输出');程序运行说明：在MATLAB主窗口中运行jisuanlunkuo.m文件可以得到考虑弹性变形后凸轮的理论轮廓曲线和未考虑弹