平面四杆机构的运动仿真模型分析
- 格式:doc
- 大小:185.00 KB
- 文档页数:7
平行四杆垂直升降机构的运动仿真与结构优化摘要:应用solidworks软件对平行四杆垂直起升机构进行了建模和运动分析,得到了起升摇臂的作用力图解,并通过有限元技术对起升摇臂进行了仿真分析,在仿真结果的基础上进行了结构优化,优化后,起升摇臂侧板质量减轻了19.7kg,降低了42%,优化结果为生产提供了理论依据。
关键词:垂直升降;运动仿真;有限元;结构优化引言:料车是一种在各种工业用炉前对物料进行搬运的专用设备,它可以实现对需处理物料的转运、装炉和出炉等功能。
料车垂直起升机构主要用以完成加料小车的垂直升降功能,是料车各组成部分中最关键的部分。
图1平行四杆垂直起升机构简图平行四杆垂直起升机构简图如图1所示,起升摇臂的下部与机架在o1处铰接,其上部在b点与液压缸的活塞杆铰接,其中部与拉杆的左端在a点铰接;拉杆的右端与起升连杆的上部铰接,并在c 点通过滚轮与导向槽板形成滚轮滑块机构;导向槽板的右端与机架在d点形成滚轮滑块机构;液压缸、起升连杆的下部均与机架铰接。
平行四杆垂直起升机构工作时,由液压缸形成原动力,推动起升摇臂绕o1点逆时针旋转,并通过铰接的a点带动拉杆向左移动,拉杆带动起升连杆绕着o2点逆时针旋转,并通过c点的滚轮迫使其上部的导向槽板产生左上方的斜向移动趋势,右端的滚轮滑块机构用于消除导向槽板的水平移动。
此时导向槽板即可带动固定其上部的重物在垂直方向上平移,以实现垂直起升目的。
一、建模与动力学仿真利用solidworks建立平行四杆垂直起升机构的三维模型,为了便于仿真,在不影响仿真结果的前提下将中部的起升连杆部分进行了简化,且右端滑块机构也使用配合功能以于实现,图2中只给出部分模型。
这样既可减少仿真模型中零、部件数量,提高仿真效率,也可减少仿真时的冗余自由度[1]。
图2平行四杆垂直起升机构模型用solidworks motion中的配合和马达功能,使每个零件的位置、速度、加速度等参数在各时间点都完全确定,确保没有冗余自由度和欠约束的存在。
基于UG软件的四连杆运动仿真分析UG软件是一款常用的CAD(计算机辅助设计)软件,它能够帮助工程师进行各种模型的建立、装配和分析。
在机械领域,UG软件被广泛应用于各类机械零部件的设计和仿真。
本文将就UG软件的四连杆运动仿真分析进行探讨,并详细介绍其原理、步骤及应用场景。
一、四连杆的基本概念四连杆是一种机械传动机构,由四条杆件和四个旋转副构成。
其中两条较长的杆件在一端旋转固定,称为地杆,另外两条较短的杆件同样旋转固定,称为摇杆。
四连杆的动作主要靠摇杆的运动驱动,使机械系统完成各种工作。
四连杆的工作原理强调套路重复的动作,即摇杆先向一个方向运动,然后再向另一个方向运动,执行往复的动作。
二、四连杆的运动仿真分析原理在使用UG软件进行四连杆运动仿真分析之前,我们需要了解一些基本原理。
首先,我们需要清楚地知道四连杆的各个参数,包括地杆长度、摇杆长度、连杆长度和摇杆旋转轴的位置等。
其次,我们还需要明确四连杆运动的动力学方程,即四个杆件的位置和速度之间的关系。
最后,我们需要掌握运动分析的方法,以便根据四连杆的参数和动力学方程,计算出各个杆件的位置和速度。
三、四连杆运动仿真分析的步骤1. 创建机械结构模型我们首先需要在UG软件中创建四连杆的机械结构模型,包括四连杆的杆件和旋转副等。
在创建过程中,需要设置结构的初始参数,如地杆长度、摇杆长度、连杆长度、摇杆旋转轴的位置等。
此外,还需要定义四连杆的运动路径和工作条件。
2. 定义杆件约束与运动学关系在创建四连杆的模型后,需要对杆件进行约束和位移关系的定义。
我们需要选择恰当的杆件,对其进行约束设置,确定其运动的自由度,以达到正确的运动效果。
同时,还需要定义杆件之间的运动学关系,解决各个杆件之间的相互作用问题。
3. 进行四连杆运动仿真完成约束和位移关系的设置后,我们就可以开始进行四连杆运动仿真。
在进行仿真前,我们需要确定仿真方案和仿真参数,如仿真时间、仿真速度和仿真环境等。
第1讲四连杆机构运动仿真一、建立连接1.设置工作目录选择【文件】→【设置工作目录】打开工作目录选取面板,如图1所示,选择如图所示2的文件夹为工作目录。
图1设置工作目录图2 选择文件夹2.建立新的装配文件打开PROE软件,点击'文件',选择‘新建’,有如下对话框弹出(如图3所示),在类型项选择‘组件’,子类型项选择‘设计’,名称改为‘2009109120’,不使用缺省模板,点击‘确定’。
有下对话框弹出(如图4所示),在模板中选择‘mmns -asm -design’,直接点击‘确定’开始进入制图过程。
图3 新建组件图4 选择单位二、装配文件1.机架的放置(1)进入PROE的主界面,点击右下角图标‘’,有如下对话框弹出(如图5所示),选择运动仿真四连杆中1ground.prt,单击打开。
图5 载入文件在主界面出现一行任务栏,在‘自动’选项中选择,再在右边单击‘’,如图6所示。
图6 机架1(2)再点击右下角图标‘’,选择运动仿真四连杆中1-ground-prt,单击‘打开’,则在主界面中出现一行任务栏,如图7所示。
图7 机架2用鼠标左键选择两平面对齐,如图8所示。
图8 平面对齐在选择两侧面对齐,在任务栏中选择,如图9所示,再单击右边''。
图9 侧面对齐2.曲柄的装配在单击右下角‘’,在运动仿真四连杆中选择‘2-crank-prt',单击‘打开’。
在主界面出现一行任务栏:,在用户定义栏中有选择'’,在操作区中选择曲柄的轴线与机座的轴线重合,如图10所示。
图10 曲柄面匹配再选择曲柄与该机座的一端面配对,如图11所示。
在任务栏中点击‘’,和‘’,完成该次联结。
图11 轴对齐3.摇杆的装配单击右下角‘’,运动仿真四连杆中选择‘4rocker-prt’,单击‘打开’。
任务栏:,同理在用户定义中选择‘’把第4摇杆与另一机座的轴线重合连结,如图12所示。
图12 轴对齐再将该摇杆与机座的端面配对连结,如图13所示。
实验六平面四杆机构实验一、实验目的1.了解曲柄摇杆机构和曲柄滑块机构位移、速度、加速度的变化规律;2.了解位移、速度、加速度的测定方法;3.增加对两种机构的运动规律的感性认识;4.比较曲柄(导杆)摇杆机构与曲柄(导杆)滑块两种机构的性能差别。
二、实验设备及工具1.QY-I曲柄(导杆)摇杆机构实验台;2.QH-I曲柄(导杆)滑块机构实验台;3.活动扳手,固定扳手,内六角扳手,螺丝刀,钢直尺。
三、实验台机械结构1.曲柄(导杆)摇杆机构实验台:图1如图1所示,其机构由曲柄、导杆、连杆、摇杆机构组成,其尺寸均可调。
可拼装曲柄摇杆机构,其底板在水平方向与机架构成一弹性系统,通过对水平方向振动变化的测试,可了解机构惯性力对机架振动的影响。
各构件长度可调范围如下:曲柄:20——60 ㎜导杆:0——50 ㎜连杆:0——220 ㎜摇杆:0——150 ㎜2.曲柄(导杆)滑块机构实验台图2如图2所示:其机构由曲柄、导杆、连杆、滑块组成,长度尺寸均可调节,可拼曲柄(导杆)滑块机构。
其联接结构与曲柄(导杆)摇杆机构相同,底板与机架的支承方式也相同。
各构件长度调节范围如下:曲柄:0——60 ㎜导杆:0——150 ㎜连杆:0——220 ㎜滑块偏心距:0——10 ㎜四、实验内容与步骤(一)系统联接及启动1.连接通讯线本实验必须通过计算机来完成。
将计算机串行口,通过标准的通讯线,连接到实验仪背面的接口。
2.打开实验软件,熟悉软件界面及各项操作的功能。
(可参阅操作系统软件简介)(二)组合机构实验操作(1)观察曲柄(导杆)摇杆机构,记录机构参数。
(2)打开实验仪上的电源,此时带有LED 数码管显示的面板上将显示"0"。
(3)启动曲柄(导杆)摇杆机构,在机构电源接通前应将电机调速电位器逆时针旋转至最低速位置,然后接通电源,并顺时针转动调速电位器,使转速逐渐加至每分 40 转,注意显示面板上实时显示曲柄轴的转速。
(4)机构运转正常后,在软件界面右侧的采样参数选择区内选择相应的采样方式和采样常数。
栏杆机四杆机构运动学分析之马矢奏春创作1 四杆机构运动学分析1.1 机构运动分析的任务、目的和方法曲柄摇杆机构是平面连杆机构中最基本的由转动副组成的四杆机构, 它可以用来实现转动和摆动之间运动形式的转换或传递动力.对四杆机构进行运动分析的意义是:在机构尺寸参数已知的情况下, 假定主动件(曲柄)做匀速转动, 撇开力的作用, 仅从运动几何关系上分析从动件(连杆、摇杆)的角位移、角速度、角加速度等运动参数的变动情况.还可以根据机构闭环矢量方程计算从动件的位移偏差.上述这些内容, 无论是设计新的机械, 还是为了了解现有机械的运动性能, 都是十分需要的, 而且它还是研究机械运动性能和动力性能提供需要的依据.机构运动分析的方法很多, 主要有图解法和解析法.当需要简捷直观地了解机构的某个或某几个位置的运动特性时, 采纳图解法比力方便, 而且精度也能满足实际问题的要求.而当需要精确地知道或要了解机构在整个运动循环过程中的运动特性时, 采纳解析法并借助计算机, 不单可获得很高的计算精度及一系列位置的分析结果, 并能绘制机构相应的运动线图, 同时还可以把机构分析和机构综合问题联系起来, 以便于机构的优化设计.1.2 机构的工作原理在平面四杆机构中, 其具有曲柄的条件为:a.各杆的长度应满足杆长条件, 即:最短杆长度+最长杆长度≤其余两杆长度之和.b.组成该周转副的两杆中必有一杆为最短杆, 且其最短杆为连架杆或机架(当最短杆为连架杆时, 四杆机构为曲柄摇杆机构;当最短杆为机架时, 则为双曲柄机构).三台设备测绘数据分别如下:, , ,最短杆长度+最长杆长度(125.36+73.4)<其余两杆长度之和(103.4+103.52)最短杆为连架杆, 四杆机构为曲柄摇杆机构图1-1 II-1型栏杆机机构测绘及其运动位置图, , ,最短杆长度+最长杆长度(125.36+50.1)<其余两杆长度之和(109.8+72.85)最短杆为连架杆, 四杆机构为曲柄摇杆机构图1-2 II-2型栏杆机机构测绘及其运动位置图, , L3=150mm,最短杆长度+最长杆长度(163.2+6)<其余两杆长度之和(150+90.1)最短杆为连架杆, 四杆机构为曲柄摇杆机构图1-3 III型栏杆机机构测绘及其运动位置图在如下图1所示的曲柄摇杆机构中, 构件AB为曲柄, 则B点应能通过曲柄与连杆两次共线的位置.曲柄摇杆机构死点情况分析:在曲柄摇杆机构中, 一般两连架杆一为主动件,一为从动件, 我们知道, 当从动件连架杆与连杆处于共线( 拉直共线或重叠共线) 位置时, 机构的传动角为0, 即机构处于死点位置, 机构在死点位置上无法启动且具有运动不确定性,主动时曲柄摇杆机构有两个死点位置, 而对曲柄主动时, 有否死点位置的问题, 基本没有涉及. 有的资料上则直接说, 曲柄主动时无死点位置. 本文对此问题进行了分析研究, 发现:曲柄主动时, 最短杆长度+最长杆长度<其余两杆长度之和,此时无死点位置.图1-4曲柄摇杆机构表1 曲柄摇杆机构的死点个数及位置情况表1.3 机构的数学模型的建立图1-5 曲柄摇杆机构数学模型简图在用矢量法建立机构的位置方程时, 需将构件用矢量来暗示, 并作出机构的封闭矢量多边形.如图1所示, 先建立一直角坐标系.设各构件的长度分别为L1 、L2 、L3 、L4 , 其方位角为、、、 .以各杆矢量组成一个封闭矢量多边形, 即ABCDA.其个矢量之和必即是零.即:式1式1为图1所示四杆机构的封闭矢量位置方程式.对一个特定的四杆机构, 其各构件的长度和原动件2的运动规律, 即为已知, 而 =0, 故由此矢量方程可求得未知方位角、 .角位移方程的分量形式为:式2闭环矢量方程分量形式对时间求一阶导数(角速度方程)为:式3其矩阵形式为:式4联立式3两公式可求得:式5式6闭环矢量方程分量形式对时间求二阶导数(角加速度方程)矩阵形式为:式7由式7可求得加速度:式8式9注:式1~式9中, Li(i=1,2,3,4)分别暗示机架1、曲柄2、连杆3、摇杆4的长度;(i=1,2,3,4)是各杆与x轴的正向夹角, 逆时针为正, 顺时针为负, 单元为 rad; 是各杆的角速度, , 单元为 rad/s; 为各杆的角加速度, 单元为.(1)求导中应用了下列公式:式10(2)在角位移方程分量形式(式2)中, 由于假定机架为参考系, 矢量1与x轴重合, =0, 则有非线性超越方程组:式11可以借助牛顿-辛普森数值解法或Matlab自带的fsolve函数求出连杆3的角位移和摇杆4的角位移.(3)求解具有n个未知量(i=1,2,…,n)的线性方程组:式12式中, 系列矩阵是一个阶方阵:式13的逆矩阵为 ;常数项b是一个n维矢量:式14因此, 线性方程组解的矢量为:式15式11是求解连杆3和摇杆4角速度和角加速度的依据.基于MATLAB法式设计四连杆机构的解析法同样可以用MATLAB 的计算工具来求值, 并结合MATLAB 的可视化手段, 把各点的计算值拟合成曲线, 获得四连杆机构的运动仿真轨迹.1.4.1 法式流程图图1-6 Matlab运动分析法式流程.2 M文件编写首先创立函数FoutBarPosition, 函数fsolve通过他确定 .function t=fourbarposition(th)%求解其他两杆的θ_3, θ_4L1=163.2mm;L2=64.25mm;L3=150mm;L4=90.1mm;%给定已知量, 各杆长L1,L2,L3,L4global th21 %给定初始θ_2t(1)= L2*cos(th21)+L3*cos(th(1))-L4*cos(th(2))-L1;t(2)=L2*sin(th21)+L3*sin(th(1))-L4*sin(th(2));主法式如下:disp ' * * * * * * 平面四杆机构的运动分析 * * * * * *' L1=;L2=;L3=;L4=;%各杆长度globalth21 %θ_2h=30; %给出转角步长30度th2=[0:h:360]*pi/180;%曲柄输入角度从0至360度, 步长为pi/6th34=zeros(length(th2),2);%建立一个N行2列的零矩阵, 第一列寄存options=optimset('display','off');%θ_3,第二列寄存θ_4form=1:length(th2)%建立for循环, 求解θ_3, θ_4th21= th2(m);y3=fsolve('fourbarposition',[1 1]); %的非线性超越方程, 结果保管在th34中th34(m,:)=y3;endy=L2*sin(th2)+L3*sin(th34(:,1)');%连杆3的C端点Y坐标值x=L2*cos(th2)+L3*cos(th34(:,1)');%连杆3的C端点X坐标值xx=[L2*cos(th2)];%连杆3的B端点X坐标值yy=[L2*sin(th2)];%连杆3的B端点Y坐标值figure(1)plot([x;xx],[y;yy],'k',[0 L1],[0 0], 'k--^', x,y,'ko', xx,yy,'ks') %绘制连杆3的几个位置点title('连杆3的几个位置点')xlabel('水平方向(m)')ylabel('垂直方向(m)')axisequal%XY坐标均衡h=5;%重新细分曲柄输入角度θ_2, 步长为5度th2=[0:h:360] *pi/180;th34=zeros(length(th2),2);options=optimset('display','off');for m=1:length(th2)%建立for循环, 求解θ_3, θ_4th21= th2(m);y3=fsolve('fourbarposition',[1 1]);th34(m,:)=y3;endfigure(2)%plot(th2*180/pi,th34(:,1),th2*180/pi,th34(:,2)) %绘制连杆3的角位移关于曲柄2的角位移图plot(th2*180/pi,th34(:,1)*180/pi,th2*180/pi,th34(:,2)* 180/pi)%绘制摇杆4的角位移关于曲柄2的角位移图axis([0 360 0170])%确定XY鸿沟值grid%图形加网格xlabel('主动件转角\theta_2(度)')%横坐题目目ylabel('从动件角位移(度)')%纵坐题目目title('角位移线图')text(120,120,'摇杆4角位移')%设定显示范围text(150,40,'连杆3角位移')w2=;%设定曲柄角速度for i=1:length(th2)A=[-L3*sin(th34(i,1)),L4*sin(th34(i,2));L3*cos(th34(i,1)), -L4*cos(th34(i,2))];B=[w2*L2*sin(th2(i)); -w2*L2*cos(th2(i))];w=inv(A)*B;w3(i)=w(1);%求解杆3角速度w4(i)=w(2);%求解杆4角速度%w3(i)=(w2*L2*sin(th34(i,2)-th2(i)))/(L3*sin(th34(i,1)-th34(i,2)));%w4(i)=(w2*L2*sin(th2(i) - th34(i,1)))/(L4*sin(th34(i,2)-th34(i,1)));endfigure(3)plot(th2*180/pi,w3,th2*180/pi,w4);%绘制角速度线图axis([0 360 -55])text(50,,'摇杆4角速度(\omega_4)')text(220,3,'连杆3角速度(\omega_3)')gridxlabel('主动件转角\theta_2(度)')ylabel('从动件角速度(rad\cdot s^{-1})')title('角速度线图')for i=1:length(th2)C=[L3*sin(th34(i,1)),-L4*sin(th34(i,2));L3*cos(th34(i,1)),-L4*cos(th34(i,2))];D=[w4(i)^2*L4*cos(th34(i,2))-w3(i)^2*L3*cos(th34(i,1))-w2^2*L2*cos(th2(i));w2^2*L2*sin(th2(i))+w3(i)^2*L3*sin(th34(i,1))-w4(i)^2*L4*sin(th34(i,2))];a=inv(C)*D;%s43=th34(:,2)-th34(:,1);%s23= th2'-th34(:,1);%a3=(-L2*w2^2.*cos(s23)-w3.^2*L3.*cos(-s43)+ L4.*w4.^2)./(L4.*sin(s43));%a4=(L2*w2^2.*cos(s23)-w4.^2*L4.*cos(s43)+L3.*w3.^2)./(L4.*sin(s43)); a3(i)=a(1);%求解杆3角加速度a4(i)=a(2);%求解杆4角加速度endfigure(4)plot(th2*180/pi,a3,th2*180/pi,a4);%绘制角加速度线图axis([0 360 -3060])text(30,18,'摇杆4角加速度(\alpha_4)')text(180,7,'连杆3角加速度(\alpha_3)')gridxlabel('主动件转角\theta_2(度)')ylabel('从动件角加速度(rad\cdot s^{-2})')title('角加速度线图')disp '曲柄转角连杆转角-摇杆转角-连杆角速度-摇杆角速度-连杆加速度-摇杆加速度'ydcs=[th2'*180/pi,th34(:,1)*180/pi,th34(:,2)*180/pi,w3 ',w4',a3',a4'];disp(ydcs)% 重新细分曲柄输入角度θ_2, 步长为1度h=1;%重新细分曲柄输入角度θ_2, 步长为度th2=[20:h:210] *pi/180;th34=zeros(length(th2),2);options=optimset('display','off');for m=1:length(th2)%建立for循环, 求解θ_3, θ_4th21= th2(m);y3=fsolve('fourbarposition',[1 1]);th34(m,:)=y3;endfigure(2)%plot(th2*180/pi,th34(:,1),th2*180/pi,th34(:,2)) %绘制连杆3的角位移关于曲柄2的角位移图plot(th2*180/pi,th34(:,1)*180/pi,th2*180/pi,th34(:,2)*180/pi)%绘制摇杆4的角位移关于曲柄2的角位移图axis([20210 0 180])%确定XY鸿沟值grid%图形加网格xlabel('主动件转角\theta_2(度)')%横坐题目目ylabel('从动件角位移(度)')%纵坐题目目title('角位移线图')text(120,120,'摇杆4角位移')%设定显示范围text(150,40,'连杆3角位移')w2=;%设定曲柄角速度for i=1:length(th2)A=[-L3*sin(th34(i,1)),L4*sin(th34(i,2));L3*cos(th34(i,1)), -L4*cos(th34(i,2))];B=[w2*L2*sin(th2(i)); -w2*L2*cos(th2(i))];w=inv(A)*B;w3(i)=w(1);%求解杆3角速度w4(i)=w(2);%求解杆4角速度%w3(i)=(w2*L2*sin(th34(i,2)-th2(i)))/(L3*sin(th34(i,1)-th34(i,2)));%w4(i)=(w2*L2*sin(th2(i) - th34(i,1)))/(L4*sin(th34(i,2)-th34(i,1)));endfigure(3)plot(th2*180/pi,w3,th2*180/pi,w4);%绘制角速度线图axis([20 210 -35])text(50,,'摇杆4角速度(\omega_4)')text(100,-1,'连杆3角速度(\omega_3)')gridxlabel('主动件转角\theta_2(度)')ylabel('从动件角速度(rad\cdot s^{-1})')title('角速度线图')for i=1:length(th2)C=[L3*sin(th34(i,1)),-L4*sin(th34(i,2));L3*cos(th34(i,1)),-L4*cos(th34(i,2))];D=[w4(i)^2*L4*cos(th34(i,2))-w3(i)^2*L3*cos(th34(i,1))-w2^2*L2*cos(th2(i));w2^2*L2*sin(th2(i))+w3(i)^2*L3*sin(th34(i,1))-w4(i)^2*L4*sin(th34(i,2))];a=inv(C)*D;%s43=th34(:,2)-th34(:,1);%s23= th2'-th34(:,1);%a3=(-L2*w2^2.*cos(s23)-w3.^2*L3.*cos(-s43)+ L4.*w4.^2)./(L4.*sin(s43));%a4=(L2*w2^2.*cos(s23)-w4.^2*L4.*cos(s43)+L3.*w3.^2)./(L4.*sin(s43)); a3(i)=a(1);%求解杆3角加速度a4(i)=a(2);%求解杆4角加速度endfigure(4)plot(th2*180/pi,a3,th2*180/pi,a4);%绘制角加速度线图axis([20210 -25 40])text(45,20,'摇杆4角加速度(\alpha_4)')text(160,5,'连杆3角加速度(\alpha_3)')gridxlabel('从动件角加速度')ylabel('从动件角加速度(rad\cdot s^{-2})')title('角加速度线图')disp '曲柄转角-连杆转角-摇杆转角-连杆角速度-摇杆角速度-连杆加速度-摇杆加速度'ydcs=[th2'*180/pi,th34(:,1)*180/pi,th34(:,2)*180/pi,w3 ',w4',a3',a4'];disp(ydcs)%RRR杆组各点约束力(动力学分析)Rbcd=zeros(length(th2),6);for m=1:length(th2)%求bcd三点约束反力M(1)=th34(m,1);M(2)=th34(m,2);M(3)= w3(m);M(4)= w4(m);M(5)= a3(m);M(6)= a4(m);M(7)=-L2*w2*w2*cos(th2(m));M(8)= -L2*w2*w2*sin(th2(m));M(9)=0;M(10)=0;M(11)=-100;Y1=RRRdy(M);Rbcd(m,:)=Y1;end%主动杆组各点约束力和力矩Ram=zeros(length(th2),3);for m=1:length(th2)N(1)=th2(m);N(2)=10;N(3)=0;N(4)=Rbcd(m,1);N(5)= Rbcd(m,2);Y1=crankdy(N);Ram(m,:)=Y1;endplot(th2*180/pi,Ram(:,3)); %绘制曲柄力矩线图axis([0 360 -5050])gridxlabel('曲柄角度')ylabel('曲柄力矩(N*m)')title('曲柄力矩线图') plot(th2*180/pi,Ram(:,1)); %绘制A 点约束力水平分力axis([0 360 -250250])gridxlabel('曲柄角度')ylabel('A 点水平分力(N)')title('A 点约束力水平分力')2栏杆机各机型的分析结果2.1 2代1机构尺寸参数, , , r4=125.36mm ;质心为rc1= mm, rc2= mm .rc3= mm质量为m1= kg, m2= kg .m3= kg ;转动惯量为J1= kg •m2, J2= kg •m2, J3= kg •m2, 构件3的工作阻力矩M3= N •m .顺时针方向,其他构件所受外力和外力矩(弹簧拉力年夜小及位置) 构件1以等角速度5.326 rad /s 逆时针方向回转曲柄角速度()tπθω*180/12=θ12为曲柄两极限点的转角范围1.θ12-=352.218θ31为摇杆两极限点转角范围θ31下表为曲柄转动一周,各参数变动量,角度间隔5度曲柄转角- 连杆转角-摇杆转角-连杆角速度-摇杆角速度-连杆加速度-摇杆加速度下表为曲柄转开工作区间30—225度,各参数变动量,角度间隔1度曲柄转角- 连杆转角-摇杆转角-连杆角速度-摇杆角速度-连杆加速度-摇杆加速度图2-1 连杆3的空间位置点图2-2 连杆3和摇杆4的角位移曲线图2-3 工作区间内连杆3和摇杆4的角位移曲线曲柄两极限点的转角范围1.θ=12-352.218摇杆两极限点转角范围θ31图2-4 连杆3和摇杆4角速度曲线图2-5 工作区间内连杆3和摇杆4角速度曲线工作区间内摇杆角速度最年夜值:工作区间内摇杆角速度最小值:图2-6 连杆3和摇杆4角加速度曲线图2-7工作区间内连杆3和摇杆4角加速度曲线工作区间内摇杆角加速度最年夜值:工作区间内摇杆角加速度最小值:2.2 2代2机构尺寸参数 各构件的尺寸为r1=50.1mm,r2=109.8mm,r3=72.85mm,r4=125.36mm ; 质心为rc1= mm, rc2= mm .rc3= mm 质量为m1= kg, m2= kg .m3= kg ;转动惯量为J1= kg •m2, J2= kg •m2, J3= kg •m2, 构件3的工作阻力矩M3= N •m .顺时针方向, 其他构件所受外力和外力矩(弹簧拉力年夜小及位置) 构件1以等角速度3.38594 rad /s 逆时针方向回转曲柄角速度()tπθω*180/12=θ12为曲柄两极限点的转角范围4.2620112-=θ θ31为摇杆两极限点转角范围θ31下表为曲柄转动一周,各参数变动量,角度间隔5度曲柄转角连杆转角-摇杆转角-连杆角速度-摇杆角速度-连杆加速度-摇杆加速度下表为曲柄转开工作区间20—210度,各参数变动量,角度间隔1度曲柄转角连杆转角-摇杆转角-连杆角速度-摇杆角速度-连杆加速度-摇杆加速度图2-8 连杆3的空间位置点图2-9 连杆3和摇杆4角位移曲线图2-10 工作区间内连杆3和摇杆4角位移曲线曲柄两极限点的转角范围4.θ=20112-26摇杆两极限点转角范围θ31图2-11 连杆3和摇杆4角速度曲线图2-12工作区间内连杆3和摇杆4角速度曲线工作区间内摇杆角速度最年夜值:曲柄转角116度, 摇杆转角, 摇杆角速度2.4237工作区间内摇杆角速度最小值:曲柄转角26度, 摇杆转角, 摇杆角速度角速度变动量:2.4237 +图2-13 连杆3和摇杆4角加速度曲线图2-14 工作区间内连杆3和摇杆4角加速度曲线工作区间内摇杆角加速度最年夜值:曲柄转角26度, 摇杆转角, 摇杆角加速度工作区间内摇杆角加速度最小值:曲柄转角181度, 摇杆转角, 摇杆角速度角加速度变动量:+2.3 3代机构尺寸参数 各构件的尺寸为r1=64.25mm,r2=150mm,r3=90.1mm,r4=163.2mm ; 质心为rc1= mm, rc2= mm .rc3= mm 质量为m1= kg, m2= kg .m3= kg ;转动惯量为J1= kg •m2, J2= kg •m2, J3= kg •m2, 构件3的工作阻力矩M3= N •m .顺时针方向, 其他构件所受外力和外力矩(弹簧拉力年夜小及位置) 构件1曲柄以等角速度逆时针方向回转 3代型运动时间曲柄角速度()tπθω*180/12=θ12为曲柄两极限点的转角范围2320312-=θ=180 θ31为摇杆两极限点转角范围θ31下表为3代栏杆机曲柄转动一周,各参数变动量,角度间隔5度(运动时间0.6s)曲柄转角连杆转角-摇杆转角-连杆角速度-摇杆角速度-连杆加速度-摇杆加速度下表为3代栏杆机曲柄转动20-210度,各参数变动量,角度间隔1度(运动时间0.6s)曲柄转角-连杆转角-摇杆转角-连杆角速度-摇杆角速度-连杆加速度-摇杆加速度图2-15 连杆3的空间位置点图2-16 连杆与摇杆的角位移曲线(0.6s) 图2-17 工作区间内连杆与摇杆的角位移曲线(0.6s)曲柄两极限点的转角范围23θ=180=12-203摇杆两极限点转角范围θ31图2-18 连杆与摇杆的角速度曲线(0.6s) 图2-19 工作区间内连杆与摇杆的角速度曲线(0.6s)工作区间内摇杆角速度最年夜值:曲柄转角110度, 摇杆转角, 摇杆角速度工作区间内摇杆角速度最小值:曲柄转角203度, 摇杆转角, 摇杆角速度角速度变动量:+图2-20 连杆与摇杆的角加速度曲线(0.6s) 图2-21工作区间内连杆与摇杆的角加速度曲线(0.6s)工作区间内摇杆角加速度最年夜值:曲柄转角23度, 摇杆转角, 摇杆角加速度工作区间内摇杆角加速度最小值:曲柄转角181度, 摇杆转角, 摇杆角速度角加速度变动量:+图2-22 工作区间内连杆与摇杆的角速度曲线(0.9s)工作区间内摇杆角速度最年夜值:曲柄转角110度, 摇杆转角, 摇杆角速度工作区间内摇杆角速度最小值:曲柄转角203度, 摇杆转角, 摇杆角速度角速度变动量:+图2-23 工作区间内连杆与摇杆的角加速度曲线(0.9s)工作区间内摇杆角加速度最年夜值:曲柄转角23度, 摇杆转角, 摇杆角加速度工作区间内摇杆角加速度最小值:曲柄转角181度, 摇杆转角, 摇杆角速度角加速度变动量:+下表为3代栏杆机曲柄转动20-210度,各参数变动量,角度间隔1度(运动时间0.9s)曲柄转角连杆转角-摇杆转角-连杆角速度-摇杆角速度-连杆加速度-摇杆加速度3代型栏杆机(运动时间1.3S)运动分析图2-24 工作区间内连杆与摇杆的角速度曲线(1.3s)工作区间内摇杆角速度最年夜值:曲柄转角110度, 摇杆转角, 摇杆角速度工作区间内摇杆角速度最小值:曲柄转角203度, 摇杆转角, 摇杆角速度角速度变动量:+图2-25 工作区间内连杆与摇杆的角加速度曲线(1.3s)工作区间内摇杆角加速度最年夜值:曲柄转角23度, 摇杆转角, 摇杆角加速度工作区间内摇杆角加速度最小值:曲柄转角181度, 摇杆转角, 摇杆角速度角加速度变动量:+下表为3代栏杆机曲柄转动20-210度,各参数变动量,角度间隔1度(运动时间1.3s)曲柄转角-连杆转角-摇杆转角-连杆角速度-摇杆角速度-连杆加速度-摇杆加速度。
平面四杆机构的运动仿真模型分析
1前言
平面四杆机构是是平面连杆机构的基础,它虽然结构简单,但其承载能力大,而且同样能够实现多种运动轨迹曲线和运动规律,因而在工程实践中得到广泛应用。
平面四杆机构的运动分析, 就是对机构上某点的位移、轨迹、速度、加速度进行分析, 根据原动件的运动规律, 求解出从动件的运动规律。
平面四杆机构的运动设计方法有很多,传统的有图解法、解析法和实验法。
随着计算机技术的飞速发展,机构设计及运动分析已逐渐脱离传统方法,取而代之的是计算机仿真技术。
本文在UG NX5环境下对平面四杆机构进行草图建模,通过草图中的尺寸约束、几何约束及动画尺寸等功能确定各连杆的尺寸,之后建立相应的连杆、运动副及运动驱动,对建立的运动模型进行运动学分析,给出构件上某点的运动轨迹及其速度和加速度变化规律曲线,文章最后简要分析几个应用于工程的平面四杆机构实例。
2平面四杆机构的建模
问题的提出
平面四杆机构因其承载能力大,可以满足或近似满足很多的运动规律,所以其应用非常广泛,本文以基于曲柄摇杆机构的物料传送机构为例,讨论其建模及运动分析。
如图1所示,ABCD为曲柄摇杆机构,曲柄AB为主动件,机构在运动中要求连杆BC的延伸线上E 点保持近似直线运动,其中直线轨迹为工作行程,圆弧轨迹为回程或空程,从而实现物料传送的功能。
平面四杆机构的建模
由于物料传送机构为曲柄摇杆机构,所以它符合曲柄存在条件。
根据机械原理课程中的应用实例[1],选取AB=100,BC=CD=CE=250,AD=200,单位均为毫米。
在UG NX5的Sketch环境里,创建如图2所示的草图,并作相应的尺寸约束和几何约束,其中EE'为通过E点的水平轨迹参考线,用以检验E点的工作行程运动轨迹。
现通过草图里的尺寸动画功能,令AB与AD 的夹角从0°到360°变化,可看到E点的变化轨迹为直线和圆弧,如图3所示为尺寸动画的四个截图,其中图3(a)中的E点为水平轨迹的起点,图3(b)中的E点为水平轨迹的中点,图3(c)中的E点为水平轨迹的终点,而图3(d)中的E点为圆弧轨迹(图中未画出)即回程的中点。
如E点轨迹不符合设计要求,则可适当调整各杆件的尺寸,再通过尺寸动画功能检验。
3平面四杆机构运动模型的建立及其分析
运动模型的建立
在UG NX5的Motion环境里,分别建立3个Link(连杆)即AB、BCE和CD以及4个Revolute Joint (旋转运动副),其中J001(即A点)和J004(即D点)为AB杆和CD杆的旋转运动副,而且J001还有一个常数驱动,其初始速度为30°/Sec,如图4所示。
需要指出的是,为了使四杆机构中的E点在后处理中能够生成运动轨迹,需要在E点添加一个关联点(Associate Point),在建立连杆BCE时将该点添加到连杆中,然后选择该点作为运动轨迹(Trace)点,在后处理中将其生成轨迹;同时为了使四杆机构中的E 点在后处理中能够生成运动规律曲线,需要在E点添加一个记号点(Marker),然后在后处理时选择该记
号点作为运动对象(Motion Object)以生成相关的变化规律曲线。
在运动模型建立后,即可对其进行解算。
解算参数设置如下:
解算类型:选择常态运行(Normal Run);
分析类型:选择运动学;
时长(Time):由于AB杆的初始速度为30(度/秒),为了使其能够旋转一周,选择时长为12
(秒),即12秒×30度/秒 =360度;
步数(Step):在速度和时长确定后,步数决定该时段内运动的快慢,这里选择步数为200。
运动分析及后处理
经过解算,即可对平面四杆机构进行运动仿真显示及其相关的后处理。
首先通过动画(Animation)可以观察机构的运动过程,并可随时暂停、回退等,如图5(a)所示为机构运动过程中某一位置。
选择动画中的轨迹(Trace)选项,继续观察机构的运动过程,即可看到E点所形成的运动轨迹,如图5(b)所示。
平面四杆机构的运动仿真显示后,还可生成指定记号点的位移、速度、加速度等规律曲线,如图6、图7和图8所示为指定记号点E点的相关曲线(包括X轴和Y轴的分量)。
4工程实例中的平面四杆机构运动模型
以上所述的平面四杆机构的建模及其运动模型的建立以及运动分析和后处理方法,完全可应用于工程中的平面四杆机构[2],以下仅简要分析其中3例,因篇幅所限不讨论有关规律曲线。
如图9所示为卡车自动卸料机构,该机构为曲柄摇块机构,其中2为车架(机架),摆动油缸3是
绕固定铰链C摆动的摇块,导杆4的一端与活塞固连,A、B、C处为转动副,3与4之间为滑动副(Slider)。
油从油缸下部流入时,推动活塞,使与车斗固连的曲柄1绕固定铰链B转动,从而达到卸料的目的,双点划线位置为卸料的初始位置。
注意这里的两个极限位置是通过UG NX5/Motion后处理中的整体机构轨迹(Trace Entire Mechanism)功能得到的,并非由几何作图获得(下同)。
如图10所示为汽车车门开关机构,该机构为偏置摇杆滑块机构,当摇杆1由动力驱动后绕固定铰链A转动时,带动连杆2(与车门固连)使滑块3沿固定导路4作往复移动,从而实现汽车车门开关的目的,图中双点划线位置为汽车车门开启状态。
这里对滑块3定义了两个运动副,即3与2之间的转动副以及3与4之间的滑动副。
如图11所示为飞机起落架机构,该机构为双摇杆机构,A、D为固定铰链,摇杆1由动力驱动,摇杆3与飞机着陆轮相连,当摇杆1转动时通过连杆2使摇杆3放下或收起,从而实现飞机起落架收放的目的,图中双点划线位置为飞机起落架收起状态。
5结论
通过UG NX5,对平面四杆机构模型进行草图建模,并建立其运动仿真模型,在运动分析后可生成机构上指定点的轨迹曲线,也可输出机构上指定点的位移、速度、加速度等变化规律曲线。
通过机构的简
要分析,证明以上方法完全可应用于工程实例中的各种平面四杆机构。