哈工程机械动力学大作业
机械动力学大作业
机械动力学大作业
一、问题及要求
建立单自由度杆机构(有无滑块均可)动力学模型,由静止启动,选择一个固定驱动力矩,绘制原动件在一周内的运动关系线图,具体
机构及参数自拟。
、建立模型
建立如图一所示机构:已知驱动力矩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)同样地,速度、加速度通过求导即可得出算式,可以编出程序。
哈工程机械动力学大作业
机械动力学大作业含弹性摆杆的铰链四杆机构动力学仿真学号:院系名称:机电工程学院专业:机械工程学生姓名:本次进行设计和分析的对象为平面铰链四杆机构,在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、载入动画,进行图像处理,绘制曲线。
哈工大机械原理大作业——凸轮——2号
哈工大机械原理大作业——凸轮——2号————————————————————————————————作者:————————————————————————————————日期:Harbin Institute of Technology机械原理大作业课程名称:机械原理设计题目:凸轮机构设计一、设计题目(1)凸轮机构运动简图:(2)凸轮机构的原始参数序号升程升程运动角升程运动规律升程许用压力角回程运动角回程运动规律回程许用压力角远休止角近休止角14 90°120°余弦加速度35°90°3-4-5多项式65°80°70°(1) 推杆升程、回程运动方程如下:A.推杆升程方程:设为1rad sω=升程位移为:()()1cos451cos1.52hsπψψψ⎡⎤⎛⎫=-=-⎢⎥⎪Φ⎝⎭⎣⎦23ψπ≤≤升程速度为:()()1100sin67.5sin1.52hvπωπψψωψ⎛⎫==⎪ΦΦ⎝⎭23ψπ≤≤升程加速度为:()()2221100cos101.25cos1.52haπωπψψωψ⎛⎫==⎪ΦΦ⎝⎭23ψπ≤≤B.推杆回程方程:回程位移为:()()345111110156s h T T T ψ⎡⎤=--+⎣⎦1029918ψπ≤≤ 回程速度为:()()22111103012h v T T T ωψ=--+'Φ 1029918ψπ≤≤ 回程加速度为:()()221111260132h a T T T ωψ=--+'Φ 1029918ψπ≤≤其中:()010s T ψ-Φ+Φ='Φ1029918ψπ≤≤ (2) 利用Matlab 绘制推杆位移、速度、加速度线图 A. 推杆位移线图clcclearx1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300); x3=linspace(10*pi/9,29*pi/18,300); x4=linspace(29*pi/18,2*pi,300); T1=(x3-10*pi/9)/(pi/2); s1=45*(1-cos(1.5*x1)) s2=90;s3=90*(1-(10*T1.^3-15*T1.^4+6*T1.^5)); s4=0;plot(x1,s1,'r',x2,s2,'r',x3,s3,'r',x4,s4,'r') xlabel('角度ψ/rad'); ylabel('位移s/mm') title('推杆位移线图') gridaxis([0,7,-10,100]) 得到推杆位移线图:B.推杆速度线图clcclearx1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300);x3=linspace(10*pi/9,29*pi/18,300);x4=linspace(29*pi/18,2*pi,300);T1=(x3-10*pi/9)/(pi/2);v1=67.5*1*sin(1.5*x1);v2=0;v3=-30*90*1*T1.^2/(pi/2).*(1-2*T1+T1.^2);v4=0;plot(x1,v1,'r',x2,v2,'r',x3,v3,'r',x4,v4,'r') xlabel('角度ψ/rad');ylabel('速度v/(mm/s)')title('推杆速度线图')Grid得到推杆速度线图:C.推杆加速度线图clcclearx1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300);x3=linspace(10*pi/9,29*pi/18,300);x4=linspace(29*pi/18,2*pi,300);T1=(x3-10*pi/9)/(pi/2);a1=101.25*1.^2.*cos(1.5*x1);a2=0;a3=-60.*90.*T1./((pi/2).^2).*(1-3*T1+2*T1.^2); a4=0;plot(x1,a1,'r',x2,a2,'r',x3,a3,'r',x4,a4,'r') xlabel('角度ψ/rad');ylabel('加速度a/')title('推杆加速度线图')Grid得到推杆加速度线图:三、凸轮机构的ds/dψ-s线图,并依次确定凸轮的基圆半径和偏距.1、凸轮机构的ds/dψ--s线图:x1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300);x3=linspace(10*pi/9,29*pi/18,300);x4=linspace(29*pi/18,2*pi,300);T1=(x3-10*pi/9)/(pi/2);s1=45*(1-cos(1.5*x1))s2=90;s3=90*(1-(10*T1.^3-15*T1.^4+6*T1.^5));s4=0;v1=67.5*1*sin(1.5*x1);v2=0;v3=-30*90*1*T1.^2/(pi/2).*(1-2*T1+T1.^2);v4=0;plot(v1,s1,'r',v2,s2,'r',v3,s3,'r',v4,s4,'r')xlabel('ds/dψ');ylabel('(位移s/mm)')title('ds/dψ—s曲线')gridaxis([-120,80,-10,100])得到ds/dψ—s曲线:2、确定凸轮的基圆半径和偏距:在dssdϕ-线图中,右侧曲线为升程阶段的类速度-位移图,作直线D t d t与其相切,且与位移轴正方向呈夹角[α1]=350, 故该直线斜率:32sin2=tan5533cos2okϕϕ⨯=⨯通过编程求其角度。
哈工大机械原理大作业连杆机构
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连杆机构运动分析课程名称:机械原理设计题目:院系:班级:设计者:学号:指导教师:设计时间:哈尔滨工业大学一、连杆机构运动分析已知机构各构件的尺寸为AB=100mm,BC=2.73AB,CD=1.36AB,CG=2.36AB,BG=4.9AB,AF=2.36AB,AD=2.87AB,DF=2AB,GE=1.45AB,EF=2.82AB,GM=1.36AB,MK=1.91AB,KD=0.54AB,KF=2.81AB,HF=3.1AB,DH=3.63AB, =135 ,构件1的角速度为w1=10rad/s。
试求构件2上点G的轨迹以及构件4、构件6和构件8的角位移、角速度与角加速度。
并对计算结果进行分析。
二、机构的结构分析,组成机构的基本杆组划分该机构由一个I级杆组以及若干个II级杆组组成。
I 级杆组为:原动件1。
II级杆组为:杆2与杆CD(RRR杆组)、滑块G、杆4(RPR杆组)、杆5、杆6(RRR 杆组)、杆7与杆8(RRR杆组)。
I级杆组 1(原动件与机架)II 级杆组 1(RRR 杆2与杆CD )II 级杆组 2(RPR 杆4与滑块3)II 级杆组 3(RRR 杆5与杆6)II 级杆组 4(RRR 杆7与杆8)三、 各基本杆组的运动分析数学模型1) 根据机构运动简图,运动副A 为固定点,其位置、速度、加速度,,,,a a ax ay ax ayx y v v a a ,和杆AB 的角位移、角速度、角加速度,,ab ab ab ϕωα均已知,并且杆AB 长度为AB l 已知,使用同一构件上点的运动分析模型,可以计算出B 点的位置、速度以及加速度,,,,b b bx by bx byx y v v a a ,。
2) 已知运动副D 的位置、速度、加速度,,,,d d dx dy dx dyx y v v a a ,,以及杆CD 的长度CD l 、杆BC 的长度BC l ,B 点的参数上面已求出,使用RRRII 级杆组运动分析模型,可以求解出运动副C 的位置、速度、加速度,,,,c c cx cy cx cyx y v v a a ,,以及杆BC 的角位移、角速度、角加速度,,bc bc bc ϕωα。
哈工大机械原理大作业2号连杆
连杆分析说明书一、题目分析干草压缩机,当杆件1以w=5π/3(n=50r/min)转动时,带动RRR杆组中2和3杆件的摆动。
通过E点的连接使杆件4带动杆件5沿着导路水平运动。
二、杆组划分该装置由一个一级杆组RR、两个二级杆组(RRR和RRP)以及机架组成。
一级杆组:二级杆组包括一个RRR杆组以及一个RRP杆组。
RRR 杆组:RRP 杆组:三、各基本杆组的运动分析数学模型建立如图所示的坐标系。
一级杆组:B点的位置w1=5*3.14159/3;a1=0;xA=0;yA=0;vxA=0;vyA=0;axA=0;ayA=0;lAB=150;xB=xA+cos(f1)*lAB;yB=yA+sin(f1)*lAB; xyvxB=vxA-w1*lAB*sin(f1);vyB=vyA+w1*lAB*cos(f1);axB=0-w1^2*lAB*cos(f1)-0*lAB*sin(f1);ayB=ayA-w1^2*lAB*sin(f1);RRR杆组:xD=400;yD=500;vxD=0;vyD=0;axD=0;ayD=0;lBC=600lCD=500;lBE=480A0=2*lBC*(xD-xB)B0=2*lBC*(yD-yB)E=(xD-xB);F=(yD-yB);C0fang=E.^2.-F.^2;C0=C0fang.^(1/2);f2=2*atan((B0+(A0.^2+B0.^2-C0.^2).^(1/2))./(A0+C0)) xC=xB+lBC*cos(f2);yC=yB+lBC*sin(f2);f3=atan((yC-yD)./(xC-xD))+pi;C2=lBC*cos(f2);S2=lBC*sin(f2);C3=lCD*cos(f3);S3=lCD*sin(f3);G1=C2.*S3-C3.*S2;w2=(C3.*(vxD-vxB)+S3.*(vyD-vyB))./G1;w3=(C2.*(vxD-vxB)+S2.*(vyD-vyB))./G1;vxC=vxB-lBC*w2.*sin(f2);vyC=vyB+lBC*w2.*cos(f2);G2=axD-axB+w2.^2.*C2-w3.^2.*C3;G3=ayD-ayB+w2.^2.*S2-w3.^2.*S3;e2=(G2.*C3+G3.*S3)./G1;e3=(G2.*C2+G3.*S2)./G1;axC=axB-lBC*e2.*sin(f2)-lBC*w2.^2.*cos(f2);ayC=ayB+lBC*e2.*cos(f2)-lBC*w2.^2.*sin(f2);xE=xB+lBE*cos(f2);yE=yB+lBE*sin(f2);vxE=vxB-lBE*w2.*sin(f2);vyE=vyB+lBE*w2.*cos(f2);axE=axB-lBE*w2.^2.*cos(f2)-lBE*e2.*sin(f2);ayE=ayB-lBE*w2.^2.*sin(f2)+lBE*e2.*cos(f2);RRP杆组:lEF=600;yF=600;A02=600-yE;f4=asin(A02/lEF);Q2=-vyE;Q3=-lEF*cos(f4);w4=Q2/Q3;Q4=ayD-ayE+lEF*w4.^2.*sin(f4);xF=xE+lEF*cos(f4);e4=Q4/Q3;vFx=vxE-lEF*w4.*sin(f4);aFx=axE-lEF*w4.^2.*cos(f4)-lEF*e4.*sin(f4);s=xF-xA;四、计算编程f1=[0:pi/180:2*pi];w1=5*3.14159/3;a1=0;xA=0;yA=0;vxA=0;vyA=0;axA=0;ayA=0;lAB=150;xB=xA+cos(f1)*lAB;yB=yA+sin(f1)*lAB;vxB=vxA-w1*lAB*sin(f1);vyB=vyA+w1*lAB*cos(f1);axB=0-w1^2*lAB*cos(f1)-0*lAB*sin(f1);ayB=ayA-w1^2*lAB*sin(f1);xD=400;yD=500;vxD=0;vyD=0;axD=0;ayD=0;lBC=600lCD=500;lBE=480A0=2*lBC*(xD-xB)B0=2*lBC*(yD-yB)E=(xD-xB);F=(yD-yB);C0fang=E.^2.-F.^2;C0=C0fang.^(1/2);f2=2*atan((B0+(A0.^2+B0.^2-C0.^2).^(1/2))./(A0+C0)) xC=xB+lBC*cos(f2);yC=yB+lBC*sin(f2);f3=atan((yC-yD)./(xC-xD))+pi;C2=lBC*cos(f2);S2=lBC*sin(f2);C3=lCD*cos(f3);S3=lCD*sin(f3);G1=C2.*S3-C3.*S2;w2=(C3.*(vxD-vxB)+S3.*(vyD-vyB))./G1;w3=(C2.*(vxD-vxB)+S2.*(vyD-vyB))./G1; vxC=vxB-lBC*w2.*sin(f2);vyC=vyB+lBC*w2.*cos(f2);G2=axD-axB+w2.^2.*C2-w3.^2.*C3;G3=ayD-ayB+w2.^2.*S2-w3.^2.*S3;e2=(G2.*C3+G3.*S3)./G1;e3=(G2.*C2+G3.*S2)./G1;axC=axB-lBC*e2.*sin(f2)-lBC*w2.^2.*cos(f2); ayC=ayB+lBC*e2.*cos(f2)-lBC*w2.^2.*sin(f2); xE=xB+lBE*cos(f2);yE=yB+lBE*sin(f2);vxE=vxB-lBE*w2.*sin(f2);vyE=vyB+lBE*w2.*cos(f2);axE=axB-lBE*w2.^2.*cos(f2)-lBE*e2.*sin(f2); ayE=ayB-lBE*w2.^2.*sin(f2)+lBE*e2.*cos(f2); lEF=600;yF=600;A02=600-yE;f4=asin(A02/lEF);Q2=-vyE;Q3=-lEF*cos(f4);w4=Q2/Q3;Q4=ayD-ayE+lEF*w4.^2.*sin(f4);xF=xE+lEF*cos(f4);e4=Q4/Q3;vFx=vxE-lEF*w4.*sin(f4);aFx=axE-lEF*w4.^2.*cos(f4)-lEF*e4.*sin(f4); s=xF-xA;f=f1./pi*180;plot(f,s);plot(f,vFx);plot(f,aFx);五、计算结果分析0123456758060062064066068070001234567-400-300-200-100100200300400速度变化曲线01234567-1000-500500100015002000250030003500。
机械动力学大作业(哈尔滨工业大学)汇编
大作业:如下图所示,两端简支梁(自己也可设定其他支承形式,如悬臂梁等),梁的截面尺寸及材料属性如下。
x, UXy, UY Rot Z9.144 m横截面积A = 1.006 x 10-2 m2截面惯性矩I = 1.774 x 10-4 m 4梁的高度 h = 0.3063 m 弹性模量E = 207x109 Nm -2密度 ρ= 7860 kg/m 3 泊松比可设定为0.25,响应计算时阻尼自己设定(小于0.05即可),要求完成如下分析计算: 1. 梁的固有频率理论值(1000Hz 以内)2. 利用有限元方法,分别用梁单元及板单元计算其固有频率及振型,同时分析不同单元数对计算结果的影响(如单元数为2、5、10、20、50、100等);3. 在有限元软件中利用模态叠加法进行响应计算(取前10阶模态):激振力为正弦形式,作用点与幅值自己设定,频率5-250Hz 范围内,步长5Hz 。
4. 计算至少3个频响函数。
撰写报告时要注意:报告要有理论说明,步骤要尽量详细。
计算结果用图表形式表示。
1 梁的固有频率理论值(1000Hz 以内)解:由参考文献可知,假设等截面简支梁长度为l ,抗弯刚度为EI ,密度为ρ,横截面积为A ,r 表示第r 阶振型,则可以得到简支梁的自由振动的固有频率为:222r r EIw l A πρ=由Matlab 计算梁的前十阶固有频率如下图所示:模态阶数 1 2 3 4 5 6 7 8 9 10 固有频率12.8 51.2 115.2 204.8 320.1 460.9 627.3 819.4 1037 1280 计算程序为:function NatFre = Natural_frequency_10 ( A , h , Rou , L , E , I ) clc; %其中,A=横截面积,h=梁高度,Rou=材料密度,L为梁长度clear; % E=材料的弹性模量,I=截面惯性矩NatFre = randn(1,10,'double');for r = 1 : 10NatFre[r] = r*r*pi*pi*sqrt( E*I /(A*h*Rou*L^3) );NatFre = NatFre / (2*pi) ;end2 利用有限元方法,分别用梁单元及板单元计算其固有频率及振型,同时分析不同单元数对计算结果的影响(如单元数为2、5、10、20、50、100等);(1)梁单元模型a 单元数为2 时的固有频率:SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE1 12.840 1 1 12 56.601 1 2 23 141.47 1 3 34 143.93 1 4 45 255.98 1 5 5第1阶模态振型(单元数为2):第2阶模态振型(单元数为2):第3阶模态振型(单元数为2):第4阶模态振型(单元数为2):第5阶模态振型(单元数为2):b 单元数为5 时的固有频率:SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE1 12.791 1 1 12 51.083 1 2 23 115.06 1 3 34 140.88 1 4 45 206.14 1 5 5第1阶模态振型(单元数为5):第2阶模态振型(单元数为5):第3阶模态振型(单元数为5):第5阶模态振型(单元数为5):c 单元数为10时的固有频率:SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE1 12.789 1 1 12 51.004 1 2 23 114.22 1 3 34 140.45 1 4 45 201.85 1 5 5第2阶模态振型(单元数为10):第3阶模态振型(单元数为10):第5阶模态振型(单元数为10):d 单元数为20时的固有频率:SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE1 12.789 1 1 12 50.999 1 2 23 114.16 1 3 34 140.34 1 4 45 201.54 1 5 5第2阶模态振型(单元数为20):第3阶模态振型(单元数为20):第5阶模态振型(单元数为20):e 单元数为50时的固有频率:SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE1 12.789 1 1 12 50.999 1 2 23 114.16 1 3 34 140.31 1 4 45 201.51 1 5 5第2阶模态振型(单元数为50):第3阶模态振型(单元数为50):第5阶模态振型(单元数为50):f 单元数为100时的固有频率:SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE1 12.789 1 1 12 50.999 1 2 23 114.16 1 3 34 140.31 1 4 45 201.51 1 5 5第1阶模态振型(单元数为100):第2阶模态振型(单元数为100):第3阶模态振型(单元数为100):第4阶模态振型(单元数为100):第5阶模态振型(单元数为100):由以上可以看出,当选定为梁单元时,选择的划分网格的单元数越多,所得到的固有频率与理论值越精确。
哈工大机械原理大作业2凸轮机构设计
机械原理大作业(二)作业名称:凸轮机构设计设计题目:23题院系:班级:设计者:学号:指导教师:设计时刻:哈尔滨工业大学机械设计1.运动分析题目:设计直动从动件盘形凸轮机构,其原始参数见下表2.确信凸轮机构推杆升程、回程运动方程(设定角速度为ω=10 rad/s)升程:0°< Φ < 120°由公式可得:s=60-60*cos(3*Φ/2);v=90*ω*sin(3*Φ/2);a=135*ω2 *cos(3*Φ/2);远停止:120°< Φ < 200°由公式可得:s=120;v=0;a=0;回程:200°< Φ < 290°由公式可得:s=h[1-(10T23-15T24+6T25)]v=(-30hω1/Φ0')T22(1–2T2+T22)a=(-60hω12/Φ0'2)T2(1–3T2+2T22)式中:T2=(Φ-Φ0-Φs)/ Φ0'近停止: 290°< Φ < 360°由公式可得:s=0;v=0;a=0;3.绘制推杆位移、速度、加速度线图(设ω=10rad/s)1) 推拉位移曲线代码:%推杆位移曲线;x=0:(pi/1000):(2*pi/3);s1=60-60*cos(1.5*x);y=(2*pi/3):(pi/1000):(10*pi/9);s2=120;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;s3=120*(1-(10*T2.^3-15*T2.^4+6*T2.^5));m=(29*pi/18):(pi/1000):(2*pi);s4=0;plot(x,s1,'b',y,s2,'b',z,s3,'b',m,s4,'b'); xlabel('角度(rad)');ylabel('行程(mm)');title('推杆位移曲线');grid;2)推杆速度曲线代码:%推杆速度曲线;w=10;x=0:(pi/1000):(2*pi/3);v1=90*w*sin(1.5*x);y=(2*pi/3):(pi/1000):(10*pi/9);v2=0;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;v3=(-30*120*w/(pi/2))*T2.^2.*(1-2*T2.^2+T2.^2); % v3=-120*w*sin(2*z-20*pi/9);m=(29*pi/18):(pi/1000):(2*pi);v4=0;plot(x,v1,'r',y,v2,'r',z,v3,'r',m,v4,'r'); xlabel('角度(rad)');ylabel('速度(mm/s)');title('推杆速度曲线(w=10rad/s)');grid;3)凸轮推杆加速度曲线代码:%凸轮推杆加速度曲线;w=10;x=0:(pi/1000):(2*pi/3);a1=135*w^2*cos(3*x/2);y=(2*pi/3):(pi/1000):(10*pi/9);a2=0;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;a3=(-60*120*w^2/(pi/2)^2)*T2.*(1-3*T2.^2+2*T2.^2); m=(29*pi/18):(pi/1000):(2*pi);a4=0;plot(x,a1,'m',y,a2,'m',z,a3,'m',m,a4,'m');xlabel('角度(rad)');ylabel('加速度(mm/s^2)');title('凸轮推杆加速度曲线(w=10rad/s)');grid;4)绘制凸轮机构的dd/dd−d线图,并依次确信凸轮的基圆半径和偏距代码:%dd/dd−d线图,确信e,s0;x=0:(pi/1000):(2*pi/3);s1=60-60*cos(1.5*x);ns1=90*sin(1.5*x);y=(2*pi/3):(pi/1000):(10*pi/9);s2=120;ns2=0;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;s3n=120*(1-(10*T2.^3-15*T2.^4+6*T2.^5));ns3=-120*10*3*T2.^2+120*15*4*T2.^3-120*6*5*T2.^4 ;m=(29*pi/18):(pi/1000):(2*pi);s4=0;ns4=0;x1=0:pi/36000:pi/2;s1n=60-60*cos(1.5*x1);v1=90*sin(1.5*x1);m1=diff(s1n);%求切线1n1=diff(v1);z=m1./n1;for i=1:length(z);if abs(z(i)+tan(-55*pi/180))<0.001;breakendendb11=s1n(i)-z(i)*v1(i);x1=-300:200;y01=z(i)*x1+b11;%切线1k1=z(i);plot(x1,y01)x3=10*pi/9:pi/36000:14*pi/9;%求切线2s3n=120*(1-(10*T2.^3-15*T2.^4+6*T2.^5));v3=-120*10*3*T2.^2+120*15*4*T2.^3-120*6*5*T2.^4 ;m3=diff(s3n);n3=diff(v3);p=m3./n3;for o=1:length(p);if abs(p(o)-tan(-25*pi/180))<0.01;breakendendo;b33=s3n(o)-p(o)*v3(o);x3=-300:700;y03=p(o)*x3+b33;%切线2plot(x3,y03);sym uv[u,v]=solve('u= 1.4281*v-81.7665','u=-0.4663*v-59.6715');%v=11.66332347972972972972972972973 x%u=-65.110107738597972972972972972973 yplot(ns1,s1,'m',ns2,s2,'b',ns3,s3n,'b',ns4,s4,'b',x1,y01,'g',x3,y03,'g',v,u,'*'); xlabel('ds/d¦µ');ylabel('S');axis([-300,200,-300,300]);title('s0,e 的确信');grid;确信凸轮基圆半径与偏距:偏距e=90mm,d020mm;基圆半径为d0=150mm。
哈工大机械原理大作业_连杆35
Harbin Institute of Technology机械原理大作业一课程名称:机械原理设计题目:连杆机构设计院系:机电学院班级:1308102分析者:学号:指导教师:设计时间:2015年06月哈尔滨工业大学一、连杆机构运动分析题目以及坐标系的建立题目中的连杆机构可以简化为图示机构二、机构的结构分析,组成机构的基本杆组划分该机构由机架、一个原动件AB和两个II级杆组组成。
有题目分析易知B点的位置、速度和加速度是我们知道的由B点我们可以推导出c点的位移速度加速度从而得到D点的轨迹最终我们可以求得E点的各个参数。
五、计算编程利用vb软件进行编程,程序如下:Public q As SingleOption ExplicitPrivate Sub Command1_Click()Dim yd, xd, yc, xb, yb, a, b, t, q, i As Singleq = 0a1.ScaleWidth = 300a1.ScaleHeight = -300a1.ScaleLeft = -150a1.ScaleTop = 250a1.Line (-300, 0)-(300, 0)a1.Line (0, -300)-(0, 300)For i = -300 To 300If i Mod 50 = 0 Thena1.Line (-5, i)-(5, i)a1.Line (i, -5)-(i, 5)a1.CurrentX = 0a1.CurrentY = ia1.Print ia1.CurrentX = ia1.CurrentY = 0a1.Print iEnd IfNextFor t = 0 To 200q = q + t / 1000xb = 100 * Sin(q)yb = -100 * Cos(q)xd = 3 / 2 * xbyd = yb + Sqr(40000 - xb * xb) / 2 a1.Line (a, b)-(xd, yd)a = xdb = ydNextEnd SubPrivate Sub Command2_Click()Dim yd, xd, yc, xb, yb, a, b, t, q, i, ve, vc, ye, ye0, yc0, ve0, vc0, ae, ac, ae0, ac0, x1 As Singlea1.ScaleWidth = 10a1.ScaleHeight = -100a1.ScaleLeft = 150a1.ScaleTop = 50ve0 = 0vc0 = 0ye = 200x1 = 151For i = 1 To 15x1 = 3.14159 / 6 + x1a1.Line (x1, 1)-(x1, -1)a1.CurrentX = x1 - 3.14159 / 12a1.CurrentY = -2a1.Print i * 30NextFor t = 0 To 3000q = q + t / 10000xb = 100 * Sin(q)yb = -100 * Cos(q)xd = (3 / 2) * xbyd = yb + Sqr(40000 - xb * xb) / 2ye = yd + Sqr(40000 - xd * xd)yc = yb - Sqr(40000 - xb * xb)ve = (ye - ye0)vc = (yc - yc0)ae = ve - ve0ac = ae - ae0a1.Line (151, -100)-(151, 100)a1.Line (t - 50, 0)-(t, 0)If t Mod 5 = 0 Thena1.Line (150.9, t)-(151.1, t)a1.CurrentX = 151a1.CurrentY = ta1.Print t * 55a1.CurrentX = 151a1.CurrentY = -ta1.Print -t * 55End Ifa1.CurrentX = 151a1.CurrentY = 0a1.Line (q - 0.2, vc0)-(q, vc)a1.Line (q - 0.2, ve0)-(q, ve)yc0 = ycye0 = yeve0 = vevc0 = vcae0 = aeac0 = acNextEnd SubPrivate Sub Command4_Click()Dim yd, xd, yc, xb, yb, a, b, t, q, i, ve, vc, ye, ye0, yc0, ve0, vc0, ae, ac, ae0, ac0, x1 As Singlea1.ScaleWidth = 10a1.ScaleHeight = -30a1.ScaleLeft = 150a1.ScaleTop = 15ve0 = 0vc0 = 0ye = 200x1 = 151For i = 1 To 15x1 = 3.14159 / 6 + x1a1.Line (x1, 0.5)-(x1, -0.5)a1.CurrentX = x1 - 3.14159 / 12a1.CurrentY = -2a1.Print i * 30Nexta1.Line (151, -100)-(151, 100)a1.Line (0, 0)-(300, 0)For t = 0 To 3000q = q + t / 10000xb = 100 * Sin(q)yb = -100 * Cos(q)xd = (3 / 2) * xbyd = yb + Sqr(40000 - xb * xb) / 2 ye = yd + Sqr(40000 - xd * xd)yc = yb - Sqr(40000 - xb * xb)ve = (ye - ye0)vc = (yc - yc0)ae = ve - ve0ac = ae - ae0a1.Line (q - 0.2, ae0)-(q, ae)a1.Line (q - 0.2, ac0)-(q, ac)If t Mod 2 = 0 Thena1.Line (150.9, t)-(151.1, t)a1.CurrentX = 151a1.CurrentY = ta1.Print t * 300a1.CurrentX = 151a1.CurrentY = -ta1.Print -t * 300End Ifyc0 = ycye0 = yeve0 = vevc0 = vcae0 = aeac0 = acNextEnd SubPrivate Sub Command5_Click()Dim yd, xd, yc, xb, yb, a, b, t, q, i, ve, vc, ye, ye0, yc0, ve0, vc0, ae, ac, ae0, ac0, x1 As Singlea1.ScaleWidth = 8.5a1.ScaleHeight = -500a1.ScaleLeft = 25a1.ScaleTop = 300ye = 200x1 = 26For t = 0 To 3000q = q + t / 10000xb = 100 * Sin(q)yb = -100 * Cos(q)xd = (3 / 2) * xbyd = yb + Sqr(40000 - xb * xb) / 2ye = yd + Sqr(40000 - xd * xd)yc = yb - Sqr(40000 - xb * xb)a1.Line (q - 1 / 10, yc0 + 300)-(q, yc + 300) a1.Line (q - 1 / 10, ye0 - 150)-(q, ye - 150) yc0 = ycye0 = yeIf t Mod 27 = 0 Thena1.Line (25.8, t)-(26.2, t)a1.CurrentX = 26a1.CurrentY = ta1.Print 16500 * (t / 270 - 1)End IfNexta1.Line (0, 270)-(50, 270)a1.Line (26, 2270)-(26, -270)For i = 1 To 15x1 = 3.14159 / 6 + x1a1.Line (x1, 275)-(x1, 265)a1.CurrentX = x1 - 3.14159 / 12a1.CurrentY = 270a1.Print i * 30NextEnd SubPrivate Sub Command3_Click()a1.ClsTimer1.Enabled = FalseEnd SubPrivate Sub Command6_Click()Timer1.Enabled = TrueEnd SubPrivate Sub Timer1_Timer()Dim yd, xd, yc, xb, yb, a, b, i, ve, vc, ye, ye0, yc0, ve0, vc0, ae, ac, ae0, ac0, x1 As Singlea1.Clsa1.ScaleWidth = 500a1.ScaleHeight = -700a1.ScaleLeft = -250a1.ScaleTop = 350a1.Line (-350, 0)-(350, 0)a1.Line (0, -350)-(0, 350)For i = -7 To 7a1.Line (-2, i * 50)-(2, i * 50)a1.CurrentX = 0a1.CurrentY = i * 50a1.Print i * 50a1.Line (i * 50, -2)-(i * 50, 2)a1.CurrentY = 0a1.CurrentX = i * 50a1.Print i * 50Nextye = 200x1 = 26q = q + 1 / 100xb = 100 * Sin(q)yb = -100 * Cos(q)xd = (3 / 2) * xbyd = yb + Sqr(40000 - xb * xb) / 2ye = yd + Sqr(40000 - xd * xd)yc = yb - Sqr(40000 - xb * xb)yc0 = ycye0 = yea1.Line (0, 0)-(xb, yb)a1.Line (0, yc)-(xd, yd)a1.Line (10, yc + 10)-(-10, yc + 10) a1.Line (10, yc - 10)-(-10, yc - 10) a1.Line (-10, yc + 10)-(-10, yc - 10) a1.Line (10, yc + 10)-(10, yc - 10)a1.Line (0, ye)-(xd, yd)a1.Line (10, ye + 10)-(-10, ye + 10) a1.Line (10, ye - 10)-(-10, ye - 10) a1.Line (-10, ye + 10)-(-10, ye - 10) a1.Line (10, ye + 10)-(10, ye - 10) End Sub得到如下结果动画D点的行程曲线速度曲线位移曲线加速度曲线六、结果分析我们从可以看出计算出的曲线与标准答案一致我们可以看出我们的结果是正确的。
哈工大机械原理大作业2号凸轮
一、如图2-1所示直动从动件盘形凸轮机构,其原始参数见表2-1。
从表2-1中选择一组凸轮机构的原始参数,据此设计该凸轮机构。
图2-1表2-1 凸轮机构原始参数升程:(0<φ<2π/3)假设w=10;s1=25-25*cos(1.5*φ);v1=37.5*w*sin(1.5*φ);a1=56.25*w^2* cos(1.5φ);远休止:(2π/3<φ<π)s2=50;v2=0;a2=0;回程·1(π<φ<1.25π)s31=50-(100/(π^2/4))*( φ-π).^2;v31=-200*w/(π^2/4)*( φ-π);a31=-200*w^2/(π^2/4);回程·2(1.25π<φ<1.5π)s32=100/(π^2/4)*(3*π/2-φ).^2; v32=-200*w/(π^2/4)*( 3*π/2-φ); a32=200*w^2/(π^2/4);近休止(1.5π<φ<2π)s4=0;v4=0;a4=0;三、求ds/dφ与s的关系:x=0:(pi/1000):(2*pi/3);s1=25-25*cos(1.5*x);ns1=37.5*sin(1.5*x);y=(2*pi/3):(pi/1000):(pi);s2=50;ns2=0;z1=(pi):(pi/1000):(1.25*pi);s31=50-(100/(pi^2/4))*(z1-(pi)).^2;ns31=-(100/(pi^2/4))*(2*z1-2*pi)z2=(1.25*pi):(pi/1000):(1.5*pi);s32=100/(pi^2/4)*(3*pi/2-z2).^2;ns32=100/(pi^2/4)*(-3*pi+2*z2)m=(1.5*pi):(pi/1000):(2*pi);s4=0;ns4=0;plot(ns1,s1,'b',ns2,s2,'b',ns31,s31,'b',ns32,s32,'b',ns4,s4,'b')凸轮机构的s s -ϕd d 线图凸轮中心应该取在公共区下方。
- 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、本次建模,设置了多个对照组,可以进行横向和纵向对比,使得结果比较可靠。