凸轮机构弹性动力学.docxVIP

高等机构学大作业 PAGE \* MERGEFORMAT PAGE \* MERGEFORMAT 7 凸轮机构的弹性动力学分析（附MATLAB代码） 【问题】已知一凸轮系统，欲使其考虑弹性因素后从动件的真实运动规律按照余弦加速度运动规律运动，建立该凸轮系统的弹性动力学模型，分析其未考虑弹性因素时从动件的运动规律，并绘制出从动件的理论运动规律及考虑弹性因素后的真实运动规律。凸轮系统的运动及动力参数自定。程序代码需提供电子版，并说明运行环境。 【解答】 建立动力学模型 取图1所示的凸轮机构为研究对象，图2为其所对应的动力学模型。 图1：凸轮机构运动简图 图2：凸轮机构的动力学模型 为使得问题简化，力学模型中忽略了凸轮轴的扭转变形、弯曲变形以及回位弹簧的阻尼作用。图2中为系统等效弹簧的刚度，为凸轮机构从动组件的阻尼系数，为回位弹簧的刚度，为回位弹簧的预紧力，为凸轮机构在从动件侧的当量质量，为与凸轮廓线有关的等效凸轮升程（图中所示的凸轮并非真正的凸轮，其廓线对应的升程与真实凸轮廓线对应的升程具备关系，其中为摇臂比。因为与仅相差一个比例系数，为了便于叙述，后文将只注重分析与从动件输出的关系，而不再专门区别与的差异），为从动件的实际升程。 建立动力学方程 该机构的自由度为1，利用牛顿第二定律建立运动微分方程： （式1） 设凸轮转动的角速度为，它与时间微分、凸轮转角微分具有关系： （式2） 将（式2）代入（式1）并整理可得： （式3） 微分方程（式3）有两层含义： = 1 \* GB3 \* MERGEFORMAT ①若已知从动件的真实运动规律，可求解出凸轮在高速运转条件下考虑弹性变形影响的理论轮廓； = 2 \* GB3 \* MERGEFORMAT ②若已知凸轮廓线，可求解考虑弹性变形的从动件的动力学响应。 运动方程的求解 （一）凸轮轮廓的设计 已知条件如下：，凸轮的转速，，，，；为避免余弦加速度运动产生的冲击，取凸轮的推程运动角和回程运动角均为，远休止角和近休止角均为，从动件的最大升程。 根据已知条件，可以确定从动件的位移方程 将上式代入（式3）可得： （式4） 由于（式4）对应的常微分方程难以求出解析解，这里利用MATLAB求解出其数值解并与位移方程比较如下图： 若需要考虑机构的弹性变形，凸轮的轮廓应按照上图中的红色曲线进行设计。 由已知廓线求解从动件的真实运动规律 由于系统的弹性变形，从动件的真实运动与等效凸轮升程不再相等，当然，从动件的真实运动速度、加速度与对应的理论值也不相等。 由于弹簧的预紧力为常数，它只影响系统振动的初始平衡位置，故在分析从动件的运动规律时不再考虑，从而（式1）被简化为： （式5） 根据振动理论，系统自由振动的固有频率 其中阻尼比 代入相关数值计算可得， 如果从动件按照前述的余弦加速度规律运动，则 将上式代入（式5）可得： （式6） 利用MATLAB求解（式6）的数值解，结果如下（图中的）： 从上图中的位移、速度、加速度分析看，当考虑到系统的弹性和阻尼后，工作端的运动规律发生改变，。只有当很大时，也就是说当系统的固有频率很大时（刚度大），而且凸轮的角速度很小时，才接近值。 程序附录： % filename: tulun.m % function: cooperate with jisuanlunkuo.m to calculate the curve of the cam function dx=tulun(a,x) %% 凸轮机构的结构参数 h=6.2*10^(-3); % 升程 m=0.08459; % 凸轮机构的当量质量 n=1200;w=n/60*2*pi; % 凸轮转速及角速度 c=55.7148; % 阻尼 k1=10400; % 回位弹簧刚度 k=3194800;