哈工程机械动力学大作业
机械动力学大作业

机械动力学大作业
一、问题及要求
建立单自由度杆机构(有无滑块均可)动力学模型,由静止启动,选择一个固定驱动力矩,绘制原动件在一周内的运动关系线图,具体
机构及参数自拟。
、建立模型
建立如图一所示机构:已知驱动力矩M=20Nm阻力矩M=10 Nm 杆1长120mm转动惯量为J i=0.16kgm2 ;杆3质量为1kg,转动惯量为
j3=0.16kgm2;杆2质量、转动惯量忽略,其他参数如图所示。
求:建立系统运动方程。
二、求解过程
M3图一
方法:利用等效力学模型法进行动力学分析
解:取杆1为等效件,有
f w3
M v= Mi - M3—(1)
jv=ji+02+j3)(—y ⑵
< XWj/
w3120 3
—二—=一=0.75 Wi 160 4
又有(2)可知J V为常数则可知:
Ki
M v=J v i
则错误!未找到引用源。
四、采用ADAMS软件或Matlab/Simulink 环境,建立机械系统的动力学模型,借助软件进行求解计算和结果分析。
(1)利用Adams软件,建模后如图:
图2
图3
(2)当杆1由图1所示位置开始运动一周,机构运动时间为0.03秒, 利用Adam漱件分析杆件1角加速度错误!未找到引用源。
随时间的变化关系图,如图4所示:
图 线系关动运的内周一在件动原
4
图
uoqB 」曰|总3。
\/」e[n6uv
0/SS3
y
s s C6唾
二 N 2
I 」善E。
哈工大机械原理大作业1-16

Harbin Institute of Technology机械原理大作业设计说明书(一)课程名称:机械原理设计题目:连杆运动分析(16)院系:能源科学与工程学院班级:1102201设计者:学号:指导教师:赵永强唐德威设计时间:2013年6月8 日哈尔滨工业大学1 连杆机构运动分析题目16:如图所示机构,已知机构各构件的尺寸为AC l =CE l =100mm ,BC l =CD l =200mm ,90BCD ∠=,构件1的角速度为10/rad s ,试求构件5的角位移、角速度和角加速度,并对计算结果进行分析。
2 分析过程2.1 建立坐标系建立以点E 为原点的固定平面直角坐标系x-E-y,如图所示:图2 机构坐标系2.2结构分析将构件BCD 分为杆3和杆4。
该机构由2个Ⅰ级杆组RR (杆1和杆5)和两个Ⅱ级杆组RRP (杆3、杆4和滑块B 、D )。
其中原动件为杆1。
现将杆组分为如下两部分:图1 机构运动简图RRPRR图3 各级杆组2.3 建立数学模型2.3.1构件1、2、3的分析原动件杆1的转角:1θ=0—360。
原动件杆1的角速度:1ω=.1θ=10/rad s原动件杆1的角加速度:..1αθ==0运动副A 的坐标:0200A A x y mm =⎫⎬=⎭运动副A 的速度及加速度都为零。
构件1为BC (RRP Ⅱ级杆组)上滑块B 的导路 滑块B 的位置为:132cos cos B A C x x s x l θθ=+=+ 132sin sin B A C y y s x l θθ=+=+消去s,得:212arcsinA l θθ=+式中:011()sin ()cos C A C A A x x y y θθ=---构件3的角速度i ω和滑块B 沿导路的移动速度D υ:.211213(Q sin Q cos )/Q ωϕθθ==-+ 1322323(Q cos Q sin )/Q D s l l υθθ⋅==-+式中:..11111211321212Q sin ;Q cos ;Q sin sin cos sin l l l θθθθθθθθ=-==+构件3的角加速度和滑块B 沿导路移动的加速度:..241513(Q sin Q cos )/Q αθθθ==-+..4325323(Q cos Q sin )/Q B s l l υθθ==-+式中:122......21142211111Q cos sin cos 2sin l l l s θθθθθθθθ=---- 122......21152211111Q sin cos sin 2cos l l l s θθθθθθθθ=+-+2.3.2 构件3,4,5的分析构件3,4,5,由1个Ⅰ级基本杆组和一个RRP Ⅱ级杆组组成,与构件1,2,3结构相同,只运动分析过程与其相反。
哈工大机械原理大作业凸轮机构第四题

Harbin Institute of Technology机械原理大作业二课程名称:机械原理设计题目:凸轮机构设计姓名:李清蔚学号:1140810304 班级:1408103指导教师:林琳一.设计题目设计直动从动件盘形凸轮机构,其原始参数见表1表一: 凸轮机构原始参数升程(mm )升程运动角(º) 升程运动规律升程许用压力角(º)回程运动角(º) 回程运动规律 回程许用压力角(º)远休止角 (º) 近休止角(º) 4090等加等减速30 50 4-5-6-7多项式60 100 120二.凸轮推杆运动规律(1)推程运动规律(等加速等减速运动)推程F0=90°①位移方程如下:②速度方程如下:③加速度方程如下:(2)回程运动规律(4-5-6-7多项式)回程00240190≤≤ϕ,F0=90°,F s=100°,F0’=50°其中回程过程的位移方程,速度方程,加速度方程如下:三.运动线图及凸轮s d ds -φ线图本题目采用Matlab 编程,写出凸轮每一段的运动方程,运用Matlab 模拟将凸轮的运动曲线以及凸轮形状表现出来。
代码见报告的结尾。
1、程序流程框图开始输入凸轮推程回程的运动方程输入凸轮基圆偏距等基本参数输出ds,dv,da 图像 输出压力角、曲率半径图像 输出凸轮的构件形状结束2、运动规律ds图像如下:速度规律dv图像如下:加速度da规律如下图:3.凸轮的基圆半径和偏距以ds/dfψ-s图为基础,可分别作出三条限制线(推程许用压力角的切界限D t d t,回程许用压力角的限制线D t'd t',起始点压力角许用线B0d''),以这三条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点,即可满足压力角的限制条件。
得图如下:得最小基圆对应的坐标位置O点坐标大约为(13,-50)经计算取偏距e=13mm,r0=51.67mm.2.绘制理论轮廓线上的压力角曲线和曲率半径曲线针对凸轮转向及推杆偏置,令N1=1凸轮逆时针转;N2=1偏距为正。
哈工大机械原理大作业

H a r b i n I n s t i t u t e o f T e c h n o l o g y机械原理大作业一课程名称:机械原理设计题目:连杆机构运动分析院系:机电学院班级: 1208105分析者:殷琪学号:指导教师:丁刚设计时间:哈尔滨工业大学设计说明书1 、题目如图所示机构,一只机构各构件的尺寸为AB=100mm,BC=4.28AB,CE=4.86AB,BE=8.4AB,CD=2.14AB,AD=4.55AB,AF=7AB,DF=3.32AB,∠BCE=139?。
构件1的角速度为ω1=10rad/s,试求构件2上点E的轨迹及构件5的角位移、角速度和角加速度,并对计算结果进行分析。
2、机构结构分析该机构由6个构件组成,4和5之间通过移动副连接,其他各构件之间通过转动副连接,主动件为杆1,杆2、3、4、5为从动件,2和3组成Ⅱ级RRR基本杆组,4和5组成Ⅱ级RPR 基本杆组。
如图建立坐标系3、各基本杆组的运动分析数学模型1) 位置分析2) 速度和加速度分析 将上式对时间t 求导,可得速度方程:将上式对时间t 求导,可得加速度方程:RRR Ⅱ级杆组的运动分析如下图所示 当已知RRR 杆组中两杆长L BC 、L CD 和两外副B 、D 的位置和运动时,求内副C的位置、两杆的角位置、角运动以及E 点的运动。
1) 位置方程由移项消去j ϕ后可求得i ϕ:式中,可求得j ϕ:E 点坐标方程:其中2) 速度方程两杆角速度方程为式中,点E 速度方程为3) 加速度方程两杆角加速度为式中,点E 加速度方程为RPR Ⅱ级杆组的运动分析(1) 位移方程(2)速度方程其中(3)加速度方程4、 计算编程利用MATLAB 软件进行编程,程序如下:% 点B 和AB 杆运动状态分析>>r=pi/180;w 1=10;e 1=0;l 1=100;Xa=0;Ya=0;Vax=0;Vay=0;aax=0;aay=0;f1=0:1: 360;% B 点位置Xb=Xa+l1*cos(r*f1);Yb=Ya+l1*sin(r*f1);% B点速度Vbx=Vax-w1*l1*sin(r*f1);Vby=Vay+w1*l1*cos(r*f1);% B点加速度abx=aax-l1*w1.^2.*cos(r*f1);aby=aay-l1*w1.^2.*sin(r*f1);% RRR2级杆组运动分析% 输入D点参数l2=428;l3=214;Xd=455;Yd=0;Vdx=0;Vdy=0;adx=0;ady=0;% 计算E点、2杆、3杆运动参数lbe=840;lce=486;a0=2*l2*(Xd-Xb);b0=2*l2*(Yd-Yb);c0=l2^2+(Xb-Xd).^2+(Yb-Yd).^2-l3^2;f2=2*atan((b0+sqrt(a0.^2+b0.^2-c0.^2))./(a0+c0)); % C点位置Xc=Xb+l2*cos(f2);Yc=Yb+l2*sin(f2);% 2杆、3杆运动参数计算dX=Xc-Xd;dY=Yc-Yd;for n=1:length(dX)if dX(n)>0&dY(n)>=0f3(n)=atan(dY(n)/dX(n));elseif dX(n)==0&dY(n)>0f3(n)=pi/2;elseif dX(n)<0&dY(n)>=0f3(n)=pi+atan(dY(n)/dX(n));elseif dX(n)<0&dY(n)<0f3(n)=pi+atan(dY(n)/dX(n));elseif dX(n)==0&dY(n)<0f3(n)=1.5*pi;elseif dX(n)>0&dY(n)<0f3(n)=2*pi+atan(dY(n)/dX(n));endendC2=l2*cos(f2);C3=l3*cos(f3);S2=l2*sin(f2);S3=l3*sin(f3);G1=C2.*S3-C3.*S2;w2=(C3.*(Vdx-Vbx)+S3.*(Vdy-Vby))./G1;w3=(C2.*(Vdx-Vbx)+S2.*(Vdy-Vby))./G1;G2=adx-abx+(w2.^2).*C2-(w3.^2).*C3;G3=ady-aby+(w2.^2).*S2-(w3.^2).*S3;e2=(G2.*C3+G3.*S3)./G1;% E点位置w=acos((l2^2+lbe^2-lce^2)/(2*l2*lbe));Xe=Xb+lbe*cos(f2-w);Ye=Yb+lbe*sin(f2-w);Vex=Vbx-lbe*w2.*sin(f2-w);Vey=Vby+lbe*w2.*cos(f2-w);aex=abx-lbe*(e2.*sin(f2-w)+w2.^2.*cos(f2-w));aey=aby+lbe*(e2.*cos(f2-w)-w2.^2.*sin(f2-w));% 计算杆5运动参数Xf=646.2912088;Yf=-268.9008617;l5=sqrt((Xe-Xf).^2+(Ye-Yf).^2);dX=Xe-Xf;dY=Ye-Yf;for n=1:length(dX)if dX(n)>0&dY(n)>=0f5(n)=atan(dY(n)/dX(n));elseif dX(n)==0&dY(n)>0f5(n)=pi/2;elseif dX(n)<0&dY(n)>=0f5(n)=pi+atan(dY(n)/dX(n));elseif dX(n)<0&dY(n)<0f5(n)=pi+atan(dY(n)/dX(n));elseif dX(n)==0&dY(n)<0f5(n)=1.5*pi;elseif dX(n)>0&dY(n)<0f5(n)=2*pi+atan(dY(n)/dX(n));endendw5=(-Vex.*sin(f5)+Vey.*cos(f5))./l5;a5=(-aex.*sin(f5)+aey.*cos(f5))./l5;% 画出各参数曲线figure(1);plot(Xe,Ye,'k');xlabel('Xe/\mm');ylabel('Ye/mm');grid on;title('E点位置');figure(2);plot(f1,f5,'k');xlabel('f/\circ');ylabel('f5/\circ');grid on;title('5杆角位移');figure(3);plot(f1,w5,'k');xlabel('f/\circ');ylabel('w5/rad/s');grid on;title('5杆角速度');figure(4);plot(f1,a5,'k');xlabel('f/\circ');ylabel('a5/rad/s2');gridon;title('5杆角加速度');Warning: Unable to interpret TeX string "Xe/\mm"5、计算结果图一:E点的运动轨迹图二:5杆角位移图三:5杆角速度图四:5杆角加速度6、计算结果分析由E点位置图像可看出,构件4做周期往复运动,由图二、三、四可看出,构件5的角位移、角速度、角加速度均成周期性变化。
哈工大机械原理大作业连杆

Harbin Institute of Technology机械原理大作业一课程名称: 机械原理 设计题目: 连杆机构运动分析 院 系: 机电工程学院 班 级: 设 计 者: 学 号: 指导教师: 设计时间:1.运动分析题目(11)在图所示的六杆机构中,已知:AB l =150mm, AC l =550mm, BD l =80mm, DE l =500mm,曲柄以等角速度1w =10rad/s 沿逆时针方向回转,求构件3的角速度、角加速度和构件5的位移、速度、加速度。
2.机构的结构分析2.1建立以点A 为原点的固定平面直角坐标系A-x, y,如下图: 2.2机构结构分析该机构由Ⅰ级杆组RR (原动件1)、Ⅱ级杆组RPR (杆2及滑块3)和Ⅱ级杆组RRP (杆4及滑块5)组成。
3.建立组成机构的各基本杆组的运动分析数学模型3.1原动件1(Ⅰ级杆组RR )由图所示,原动件杆1的转角a=0-360°,角速度1w =10rad/s ,角加速度1a =0,运动副A 的位置坐标A x =A y =0,速度(A ,A),加速度(A ,A ),原动件1的长度AB l =150mm 。
求出运动副B 的位置坐标(B x , B y )、速度(B ,B )和加速度(B ,B )。
3.2杆2、滑块3杆组(RPR Ⅱ级杆组)已出运动副B 的位置(B x , B y )、速度(B ,B )和加速度(B ,B ),已知运动副C 的位置坐标C x =0, C y =550mm,速度,加速度,杆长AC l =550mm 。
求出构件2的转角b,角速度2w 和角加速度2a . 3.3构件二上点D 的运动已知运动副B 的位置(B x , B y )、速度(B ,B )、加速度(B ,B ),已经求出构件2的转角b ,角速度2w 和角加速度2a ,杆BD 的长度BD l =80mm 。
根据Ⅰ级杆组RR 的运动分析数学模型求出点D 的位置坐标(D x ,D y )、速度(D ,D )和加速度(D ,D )。
哈工大机械原理大作业3

Harbin Institute of Technology机械原理大作业三课程名称:设计题目:院系:班级:设计者:学号:指导教师:设计时间:哈尔滨工业大学大作业3 齿轮传动设计 1、设计题目1.1机构运动简图1.2机械传动系统原始参数2、传动比的分配计算由已知条件,电动机转速n=1450r/min ,输出转速n 1=27 r/min ,n 2=31 r/min ,n 3=37 r/min ,带传动最大传动比max p i =2.5,滑移齿轮传动最大传动比=4,定轴齿轮传动最大 传动比=4。
可求得:传动系统的总传动比为:11n ni == 1450/27=53.70322n ni == 1450/31=46.774 33n ni == 1450/37=39.189 传动系统的总传动比分别由带传动、滑移齿轮传动和定轴齿轮传动三部分实现。
设带传的传送比为其最大传送比5.2max =p i ,滑移齿轮的传动比为321,,v v v i i i ,定轴齿轮传动的传动比为f i ,则总传动比由于1i > 2i > 3i ,故取1max 4v v i i ==则定轴齿轮传动部分的传动比为1max max5.37f p v i i i i ==滑移齿轮传动的传动比22max3.49v f p i i i i ==33m a x2.92v f p i i i i ==定轴齿轮传动由3对齿轮传动组成, 每对齿轮的传动比为:1.754d i ==≤3、齿轮齿数的确定滑移齿轮齿数3=v i 65622.9521z z == 2=v i 8766 3.4719z z ==1=v i 10967 3.9417z z ==齿轮7,齿轮8:719z = 866z =781()852a m z z =+=齿轮9,齿轮10:917z =1067z =此时已知条件为'a =85mm ,910211()842()ni i a m z z X X ==+=-∑mm ''arccos(cos )21.78a aαα==总变位系数:'910()0.552tan z z x inv inv ααα∑+=-=根据x ∑值和1093.94 3.0z uz ==>,按选择变位系数线图左部斜线⑤分配变位系数,得90.45x =齿轮5,齿轮6:5=21z662z =此时已知'a =85mm ,561()832a m z z =+= ''arccos(cos )23.42a a αα==100.10x =总变位系数:'65() 1.082tan z z x inv inv ααα∑+=-=根据x ∑值和652.953z u z ==>,按选择变位系数线图左部斜线④分配变位系数,得 50.5x =60.58x =定轴圆柱齿轮齿数=d i 1214111326 1.5317z z z z ===齿轮11,齿轮12:角度变位正传动。
哈工大机械原理大作业——凸轮——5号

Harbin Institute of Technology机械原理大作业二课程名称:机械原理设计题目:凸轮机构设计院系:机电工程学院班级:设计者:学号:指导教师:设计时间:凸轮机构设计一. 设计题目表1 凸轮机构原始参数二. 凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图1. 推杆升程运动方程(ϕ为转角)(203ϕπ≤≤) 令1/rad s ω= 位移:001212040[sin()]sin(3)2s h ϕπϕϕϕφπφππ=-=- 速度:1002120120[1cos()]cos(3)h v ωπϕϕφφππ=-=- 加速度:2120022360[sin()]sin(3)h a πωπϕϕφφπ==2. 推杆回程运动方程(ϕ为转角)7563πϕπ≤≤ 位移:''00128001604014[1sin()]sin(4)]233T s h T πϕϕπφπφππ=-+=-+- 速度:1''00216016014[1cos()]cos(4)]3h v T ωπϕϕπφφππ=--=-+- 加速度:21'2'002264014sin()sin(4)]3h a T πωπϕπφφπ=-=-- 式中:027()()326s T πππϕϕϕϕϕ=-+=-+=-由MATLAB编程得到线位移图像:线速度图像:线加速度图像:三.绘制机构的ds sdϕ-线图由dsds vdtdd wtϕϕ==可知线图即vsw-线图,由MATLAB编程后得图四.确定基圆半径和偏距1.画切线图如上页图,得其切线图2.求基圆半径和偏距从图中读取e=50mmr0=111.80mm五.画理论及实际轮廓附:MATLAB编程1.推杆位移线图代码%升程阶段t=0:0.001:2*pi/3;s=120.*t./pi-40./pi.*sin(3.*t);hold onplot(t,s);%远休止阶段t=2*pi/3:0.001:7*pi/6;s=80;hold onplot(t,s);%回程阶段t=7*pi/6:0.001:5*pi/3;s=800/3-160.*t./pi+40./pi.*sin(4.*t-14*pi/3) hold onplot(t,s);%近休止阶段t=5*pi/3:0.001:2*pi;s=0;hold onplot(t,s);grid onhold off2.推杆速度线图代码%升程阶段t=0:0.001:2*pi/3;v=120/pi-120/pi*cos(3.*t);hold onplot(t,s);%远休止阶段t=2*pi/3:0.001:7*pi/6;v=0;hold onplot(t,s);%回程阶段t=7*pi/6:0.001:5*pi/3;v=-160/pi+160/pi*cos(4.*t-14/3*pi) hold onplot(t,s);%近休止阶段t=5*pi/3:0.001:2*pi;v=0;hold onplot(t,s);grid onhold off3.推杆加速度线图代码%升程阶段t=0:0.001:2*pi/3;a=360/pi.*sin(3.*t)hold onplot(t,a);%远休止阶段t=2*pi/3:0.001:7*pi/6;a=0;hold onplot(t,a);%回程阶段t=7*pi/6:0.001:5*pi/3;a=-640/pi.*sin(4.*t-14/3*pi) hold onplot(t,a);%近休止阶段t=5*pi/3:0.001:2*pi;a=0;hold onplot(t,a);grid onhold off4.dssdϕ-线图代码%升程阶段t=0:0.001:2*pi/3;s=120.*t./pi-40./pi.*sin(3.*t);v=120/pi-120/pi*cos(3.*t);hold onplot(v,s);%远休止阶段t=2*pi/3:0.001:7*pi/6;s=80;v=0;hold onplot(v,s);%回程阶段t=7*pi/6:0.001:5*pi/3;s=800/3-160.*t./pi+40./pi.*sin(4.*t-14*pi/3);v=-160/pi+160/pi*cos(4.*t-14/3*pi)hold onplot(v,s);5.最终轮廓线图代码h=80;w=1;e=50;rr=20;s0=100;m=120*pi/180;ms=(120+90)*pi/180;m1=(120+90+90)*pi/180; for i=1:1:120mm(i)=i*pi/180.0;s1=h.*( mm(i)./m -(1/(2*pi)).*sin(2.*pi.*mm(i)./m));v1=(w*h./m)*(1-cos(2*pi*mm(i)/m));x(i)=(s0+s1)*sin(mm(i))+e*cos(mm(i));y(i)=(s0+s1)*cos(mm(i))-e*sin(mm(i));a(i)=(s0+s1)*cos(mm(i))-e*sin(mm(i))+v1/w*sin(mm(i)); b(i)=-(s0+s1)*sin(mm(i))-e*cos(mm(i))+v1/w*cos(mm(i));xx(i)=x(i)+rr*b(i)/sqrt(a(i)*a(i)+b(i)*b(i));yy(i)=y(i)-rr*a(i)/sqrt(a(i)*a(i)+b(i)*b(i));endfor i=121:1:210mm(i)=i*pi/180;s2=h;v2=0;x(i)=(s0+s2)*sin(mm(i))+e*cos(mm(i));y(i)=(s0+s2)*cos(mm(i))-e*sin(mm(i));a(i)=(s0+s2)*cos(mm(i))-e*sin(mm(i))+v2/w*sin(mm(i));b(i)=-(s0+s2)*sin(mm(i))-e*cos(mm(i))+v2/w*cos(mm(i));xx(i)=x(i)+rr*b(i)/sqrt(a(i)*a(i)+b(i)*b(i));yy(i)=y(i)-rr*a(i)/sqrt(a(i)*a(i)+b(i)*b(i));endfor i=211:1:300mm(i)=i*pi/180;s3=800/3-160/pi.*mm(i)+40/pi.*sin(4.*mm(i)-14/3*pi);v3=-160/pi+160/pi.*cos(4*mm(i)-14/3*pi)x(i)=(s0+s3)*sin(mm(i))+e*cos(mm(i));y(i)=(s0+s3)*cos(mm(i))-e*sin(mm(i));a(i)=(s0+s3)*cos(mm(i))-e*sin(mm(i))+v3/w*sin(mm(i));b(i)=-(s0+s3)*sin(mm(i))-e*cos(mm(i))+v3/w*cos(mm(i));xx(i)=x(i)+rr*b(i)/sqrt(a(i)*a(i)+b(i)*b(i));yy(i)=y(i)-rr*a(i)/sqrt(a(i)*a(i)+b(i)*b(i));endfor i=301:1:360mm(i)=i*pi/180;x(i)=(s0+0)*sin(mm(i))+e*cos(mm(i));y(i)=(s0+0)*cos(mm(i))-e*sin(mm(i));a(i)=(s0+0)*cos(mm(i))-e*sin(mm(i))+v3/w*sin(mm(i));b(i)=-(s0+0)*sin(mm(i))-e*cos(mm(i))+v3/w*cos(mm(i));xx(i)=x(i)+rr*b(i)/sqrt(a(i)*a(i)+b(i)*b(i));yy(i)=y(i)-rr*a(i)/sqrt(a(i)*a(i)+b(i)*b(i));endplot(x,y,'r',xx,yy,'g')text(0,20,'实际轮廓线')text(120,100,'理论轮廓线')hold on。
哈工大机械原理大作业连杆机构9

机械原理大作业大作业一:连杆机构运动分析学生姓名:学号:指导教师:丁刚完成时间:机电工程学院机械设计系制二〇一八年四月连杆机构运动分析1题目(9)图1 设计题目在图1所示的机构中,已知l AB=60mm,l BC=180mm,l DE=200mm,l CD=120mm,l EF=300mm,h=80mm,h1=85mm,h2=225mm,构件1以等角速度ω1=100rad/s 转动。
求在一个运动循环中,滑块5的位移、速度和加速度曲线。
2分析结构1、杆1为RR主动件,绕A以ω1 转动,自由度1.2、4杆和滑块5为RRP II级杆组.,自由度0.3、2,3杆组成II级杆组RRR,自由度0.总共自由度为F=5*3-2*7=1 .由上述的杆组类型,确认出所需运动分析数学模型:同一构件上的点、RRP、RRR。
3.杆组法对平面连杆机构进行运动分析3.1对主动件杆1 RR I级构件的分析主动杆1转角:φ= [0°,360°) δ=0°,则φ’=ω1=100 rad/s角加速度φ’’=0.已知h2=225mm, h=80mm, l AB=60mm 所以A(225mm,80mm)A点速度(0,0),加速度(0,0)B点位置(x A+l AB*cos(φ), Y A+l AB*sin(φ))B点速度(-l AB*sin(φ), l AB*cos(φ)),加速度(-l AB*cos(φ), -l AB*sin(φ))3.2RRRII 级杆组分析(模型参考教材P37-38)图3 如图所示两个构件组成II 级杆组。
已知了B 的位置(x B ,y B )= (x A +l AB *cos(φ), Y A +l AB *sin(φ)),速度(x ’B ,y ’B ) 和加速度(x ’’B ,y ’’B ), 已知运动副D (0,0), 还可知,x ’D =y ’D =0, x ’’D =y ’’D =0. l BC =180 mm, l CD = 120mm所以,x c =x D +l CD *cos(φi)= x B +l BC *cos(φj) y c =x D +l CD *sin(φi)= x B +l BC *sin(φi) 对于φ的求解: A 0=2*l CD (x B -x D ) B 0=2*l BC (y B -y D ) C 0=l CD 2+ l BD 2- l BC 2为了保证机构的装配正常:l BD ≤l CD + l BC AND l BD ≥Abs (l CD - l BC )可求3杆的转角φi=2*arctan((B 0±sqrt (A 02 + B 02- C 02))/(A 02+ C 02)),角速度w3=φi ’和角加速度α3= φi ’’3.3 同一构件上的点(模型参考书P35-36)Φiφjφi已知D(0,0),速度(0,0),加速度(0,0),3杆转角φi 角速度φi’角加速度φi’’,Φi和它的导数在3.2都有体现LDE= 200mm可求出E的坐标,速度,加速度.x E =x c+lCE*cos(φi)y E =x C+lCE*sin(φi)同样地,速度、加速度通过求导即可得出算式,可以编出程序。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
机械动力学大作业含弹性摆杆的铰链四杆机构动力学仿真学号:院系名称:机电工程学院专业:机械工程学生:本次进行设计和分析的对象为平面铰链四杆机构,在Adams的环境下,通过对四杆机构进行建模以及运动仿真,绘制出摆杆的相关曲线图。
为了形成有效的对比,先建立含有刚性摆杆的四杆机构,进行运动仿真,绘制出摆杆的相关曲线。
再建立含有柔性摆杆的铰链四杆机构,所有参数设置均和刚性摆杆一样。
考虑到弹性摇杆可能发生较大的形变,不利于观测,绘制摇杆运动曲线时选择摇杆的质心作为参考点。
在Adams中主要有三种方法创建柔性构件,第一种是将刚性构件离散化后采用柔性梁连接;第二种是直接将刚体替换为柔性体;第三种是运用有限元分析的方法建立柔性构件。
本次建模,主要采用前两种方法建立柔性摆杆。
运用有限元建立柔性构件,等以后再进行深入研究。
同时两种方法建立的柔性杆可以形成对比。
通过本次设计,主要学习了Adams 软件建模以及运动仿真、图形处理、刚柔混合建模的操作方法,对自己也是一个很大锻炼和提升。
设计的为平面曲柄摇杆机构。
相关参数如:曲柄长L=200mm,宽W=60mm,高D=30mm;连杆长L=427mm,宽W=30mm,高D=20mm;摇杆长L=403mm,宽W=40mm,高D=20mm;机架长L=600mm,宽W=40mm,高D=20mm;曲柄角速度为40deg/sec。
经过验证,最短杆长度加上最长杆长度小于中间两根杆的长度之和,满足曲柄存在的条件,且最长杆为机架,故为曲柄摇杆机构。
一、建模过程1、建立四个标记点,这四个点依次连接就可以确定一个铰链四杆机构。
2、建立四根杆的模型3、在杆件之间添加转动副4、选择最长杆为机架并固定5、给曲柄添加驱动,使曲柄角速度为40deg/sec。
6、使模型的显示方式为实体显示7、进行运动仿真,主要通过时间和步长来控制仿真运动的快慢。
8、载入动画,进行图像处理,绘制曲线图。
9、将刚性摇杆换成柔性的摇杆并添加转动副和驱动10、进行运动仿真11、载入动画,进行图像处理,绘制曲线。
12、将连杆改为柔性体,添加转动副,并生成模型,进行运动仿真。
13、载入动画,生成曲线图。
14、将摇杆换成弹性体,生成模型,进行运动仿真。
15、载入动画,进行图像处理,生成数据。
二、运动学及动力学分析图1 刚性摆杆质心位移图2 柔性摆杆质心位移图3 柔性摆杆质心位移分析图1,图2,图3分析摆杆质心位置的变化曲线。
刚性和柔性摆杆质心的初始位置的X坐标均为500mm,Y坐标均为175mm,Z坐标均为0,说明两次运动过程起始位置一模一样,两个模型具有对比性。
刚性摆杆和柔性摆杆质心的位置均按类似余弦曲线变化,说明摇杆在一定角度围转动,与预期相符合。
由图中观察和数据分析,周期T大概为9s左右,从而验证了角速度 为40 。
其中刚性摆杆最低点Y坐标大约为50mm左右,说明摇杆没有运动到水平位置,最高点大约为200mm左右,说明当曲柄和连杆共线时,摇杆接近和机架垂直。
柔性摆杆质心最低位置Y坐标大约25mm,说明摇杆也没有运动到水平位置,最高位置Y坐标也约为200mm,说明曲柄和连杆共线时,摇杆也能和机架接近垂直。
柔性摆杆最低点Y坐标小于刚性摆杆,可能是由于杆发生了弯曲变形。
在任意时刻,柔性摆杆的X坐标均大于刚性摆杆的,说明杆发生弯曲了,使得X坐标增大。
Z坐标没有变化且为0,说明机构在XY平面运动,与预期相符合。
由图还看出曲柄摇杆机构具有急回特性。
无论在什么方向上,摆杆质心的坐标始终为正,说明摆杆在XY 平面的第一象限运动,并且摆角不超过180度。
图4 刚性摆杆质心速度曲线图5 柔性摆杆质心速度曲线图6 柔性摆杆质心速度曲线分析图4,图5,图6分析摆杆质心速度变化曲线。
在Z方向上,刚性杆和柔性杆质心速度始终为0,可见机构在XY平面运动。
在X方向上,刚性杆质心初始速度为-57mm/s,而柔性杆在X方向上速度为0,由此可以看出,对于刚性杆,任一点的速度都可以用质心的速度来代替,而柔性杆则不可以。
刚性杆运动时整个杆件同时运动,而柔性杆则是由一部分带动另一部分类似于波的传播过程。
从图中看出,无论刚性还是弹性的摆杆,质心的速度均有正有负,说明摇杆是往复摆动,符合实际。
在Y方向上,刚性摆杆质心初始速度为-33mm/s,而柔性杆初始速度为0,也说明了刚性构件和弹性构件在运动时的不同。
刚性摆杆质心最低速度大于25mm/sec,而柔性摆杆质心最低速度小于25mm/sec,说明柔性摆杆由于弯曲变形,使得运动速度变得缓慢。
刚性摆杆和柔性摆杆质心的最大速度几乎相同,说明此时虽然柔性摆杆发生了弹性形变,但是对速度的影响很低。
刚性摆杆速度在最大值附近瞬间变化,而柔性摆杆几乎是缓慢过渡。
在Adams中,对刚性摆杆质心的位移曲线对时间求导数可以得到质心的速度变化曲线,和直接生成的质心速度曲线一致。
从中可以看出,速度的变化周期为9s,和预期相符合。
而对柔性杆,由于是将构件离散化后采用柔性梁连接,属于不连续体,不能直接对质心的位移变化曲线求导来得到相应的速度曲线。
虽然X,Y方向速度有正有负,但是空间的速度始终为正。
图7 刚性摆杆质心加速度曲线图8 柔性摆杆质心X方向加速度图9 柔性摆杆质心加速度分析图7,图8,图9由图可以看出,刚性摆杆和柔性摆杆区别十分明显。
刚性摆杆呈现明显的周期性,而柔性摆杆只在X方向上有周期性。
对于刚性摆杆,加速度最小值为-50mm/s2,最大值为60mm/s2。
分析图9可知,对于柔性摆杆,在刚启动瞬间,加速度有一个突变,之后几乎一直为零,说明柔性体在Y方向上几乎不受力。
由图7和图8均可分析出,摇杆运动的周期为9s。
图10 刚性摆杆质心角速度图11 柔性摆杆质心角速度图12 柔性摆杆质心角速度分析角速度图10和图11说明,在绕X轴和绕Y轴,摆杆的角速度都为0,说明角速度为绕Z轴施加。
在绕Z轴,角速度有正有负,说明摆杆是往复摆动。
绕Z轴的角速度和摆杆在空间的角速度曲线在0—3秒完全重合,说明绕Z轴的角速度即为摆杆的角速度。
虽然绕Z轴角速度有正,有负,但是在空间,摆杆的角速度始终为正。
刚性摆杆最大角速度为22.5deg/sec,最小角速度为-22.5dec/sec。
说明刚性摆杆是左右对称的往复运动。
柔性摆杆最大角速度36dec/sec,最小角速度-33dec/sec,最大和最小角速度绝对值不等,说明摆杆在经过同一位置时往复运动速度不一样。
从图可以清楚的看出摆杆运动的周期为9s。
摆杆在空间的角速度始终为正,当绕Z轴的角速度为负时,将曲线沿着X轴翻转就得到了摆杆在空间的曲线。
图12说明,启动瞬间,摆杆的角速度有一个变化的过程,之后和图10,图11几乎一样,说明采用两种方法建立的柔性体十分类似。
图13 刚性摆杆质心角加速度图14柔性摆杆质心角加速度图15 柔性摆杆质心角加速度分析角加速度曲线对比分析图13和图14,在X,Y方向上,角加速度均为0,和角速度曲线符合。
在0到6s,刚性摆杆角加速度有两个最大值,分别为19.5,7.5,一个最小值点为4.5,在1.5秒和6秒的时候,角加速度为0.且在空间的角加速度为沿着轴的角速度沿着X轴翻转。
整个周期最大角加速度为24,最小角加速度为-19.5.角加速度的变化周期为9秒。
柔性摆杆有三个最值,两个最小值,两个0点。
变化周期为9秒。
而对于图15,空间角加速度在启动瞬间有一个突变,之后一直为0 。
图16 连杆和刚性摇杆之间的转动副受力曲线图17 连杆与柔性摇杆之间的转动副受力曲线图18 连杆与柔性摇杆之间的转动副受力曲线分析受力曲线对比分析可知,转动副受力明显不一样,刚性的时候转动副受力比较小,而柔性的时候转动副受力明显比较大。
在空间,转动副受到的力始终为正。
刚性摆杆在X方向,受力始终为负值,Y方向上有正值和负值,Z方向上不受力。
柔性摆杆在X方向上受力有正值和负值,在Y方向上有正值和负值,在Z方向上不受力。
图18说明,启动瞬间,转动副受力有突变。
图19 曲柄角速度曲线图19 说明,曲柄是匀速转动的,从而验证了其它曲线图的正确性。
从角速度和周期的角度来看,互相印证,说明了本次结果的可靠性。
三、结论1、绝对刚体在现实生活中是不存在的,平时我们在做有关杆件的计算时,如果不考虑应力和应变,可以将杆件看成刚体。
如果需要研究应力和应变,可先将物体当作刚体,用理论力学的方法求得加给它的各未知力,然后再用变形体力学,包括材料力学、弹性力学、塑性力学等的理论和方法进行研究。
2、作为传动构件,如果传递的载荷比较大,则必须要考虑到构件的变形,否则可能引起重大的事故发生。
3、在曲柄摇杆机构中,曲柄的刚性也很重要,如果曲柄很容易发生变形,很有可能无法带动整个机构实现运动。
4、对于无法直接获得的柔性体,可以考虑采用离散柔性连接件的方法来获得柔性体。
对于任何一个刚体,可以将其分成多个小块,理论上划分的小块越多,越接近柔性体。
这在数学方法上类似于以直代曲。
对于任意一条平面曲线,对曲线进行划分,如果划分的段数足够多,那么每一个小段近似都是一条直线。
对于刚体,也类似。
同样,对于弹性体来说,虽然整体发生了变形,但是把整体划分成足够多的小部分之后,每一部分发生的变形很小,可以看成刚体。
5、本次建模,设置了多个对照组,可以进行横向和纵向对比,使得结果比较可靠。