平面连杆机构设计方案

平面连杆机构设计方案

1平面连杆机构的运动分析

1.1 机构运动分析的任务、目的和方法

曲柄摇杆机构是平面连杆机构中最基本的由转动副组成的四杆机构，它可以用来实现转动和摆动之间运动形式的转换或传递动力。

对四杆机构进行运动分析的意义是：在机构尺寸参数已知的情况下，假定主动件（曲柄）做匀速转动，撇开力的作用，仅从运动几何关系上分析从动件（连杆、摇杆）的角位移、角速度、角加速度等运动参数的变化情况。还可以根据机构闭环矢量方程计算从动件的位移偏差。上述这些容，无论是设计新的机械，还是为了了解现有机械的运动性能，都是十分必要的，而且它还是研究机械运动性能和动力性能提供必要的依据。

机构运动分析的方法很多，主要有图解法和解析法。当需要简捷直观地了解机构的某个或某几个位置的运动特性时，采用图解法比较方便，而且精度也能满足实际问题的要求。而当需要精确地知道或要了解机构在整个运动循环过程中的运动特性时，采用解析法并借助计算机，不仅可获得很高的计算精度及一系列位置的分析结果，并能绘制机构相应的运动线图，同时还可以把机构分析和机构综合问题联系起来，以便于机构的优化设计。

1.2 机构的工作原理

在平面四杆机构中，其具有曲柄的条件为：

a.各杆的长度应满足杆长条件，即：

最短杆长度+最长杆长度≤其余两杆长度之和。

b.组成该周转副的两杆中必有一杆为最短杆，且其最短杆为连架杆或机架（当最短杆为连架杆时，四杆机构为曲柄摇杆机构；当最短杆为机架时，则为双曲柄机构）。

在如下图1所示的曲柄摇杆机构中，构件AB为曲柄，则B点应能通过曲柄与连杆两次共线的位置。

1.3 机构的数学模型的建立

1.3.1建立机构的闭环矢量位置方程

在用矢量法建立机构的位置方程时，需将构件用矢量来表示，并作出机构的封闭矢量多边形。如图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.3.2求解方法

（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 是Mathworks 公司推出的交互式计算分析软件，具有强大的运算分析功能，具有集科学计算、程序设计和可视化于一体的高度集成化软件环境，是目前国际上公认的最优秀的计算分析软件之一，被广泛应用于自动控制、信号处理、机械设计、流体力学和数理统计等工程领域。通过运算分析，MATLAB 可以从众多的设计方案中寻找最佳途径，获取最优结果，大大提高了设计水平和质量。四连杆机构的解析法同样可以用MATLAB 的计算工具来求值，并结合MATLAB 的可视化手段，把各点的计算值拟合成曲线，得到四连杆机构的运动仿真轨迹。

2.1 程序流程图

2.2 M文件编写

首先创建函数FoutBarPosition，函数fsolve通过他确定。

function t=fourbarposition(th,th2,L2,L3,L4,L1)

t=[L2*cos(th2)+L3*cos(th(1))-L4*cos(th(2))-L1;…

L2*sin(th2)+L3*sin(th(1))-L4*sin(th(2))];

主程序如下：

disp ' * * * * * * 平面四杆机构的运动分析 * * * * * *'

L1=304.8;L2=101.6;L3=254.0;L4=177.8; %给定已知量，各杆长L1,L2,L3,L4

th2=[0:1/6:2]*pi;

%曲柄输入角度从0至360度，步长为pi/6

th34=zeros(length(th2),2);

%建立一个N行2列的零矩阵，第一列存放options=optimset('display','off'); %θ_3,第二列存放θ_3

for

m=1:length(th2)

%建立for循环，求解θ_3，θ_4

th34(m,:)=fsolve('fourbarposition',[1 1],… %调用fsove函数求解关于θ_3，θ_4

options,th2(m),L2,L3,L4,L1); %的非线性超越方程，结果保存在th34中

end

y=L2*sin(th2)+L3*sin(th34(:,1)'); %连杆3的D端点Y坐标值

x=L2*cos(th2)+L3*cos(th34(:,1)'); %连杆3的D端点X坐标值

xx=[L2*cos(th2)];

%连杆3的C端点X坐标值

yy=[L2*sin(th2)];

%连杆3的C端点Y坐标值

figure(1)

plot([x;xx],[y;yy],'k',[0 L1],[0 0],… %绘制连杆3的几个位置点

'k--^',x,y,'ko',xx,yy,'ks')

title('连杆3的几个位置点')

xlabel('水平方向')

ylabel('垂直方向')

axis

equal

%XY坐标均衡

th2=[0:2/72:2]*pi;

%重新细分曲柄输入角度θ_2，步长为5度

th34=zeros(length(th2),2);

options=optimset('display','off');

for m=1:length(th2)

th34(m,:)=fsolve('fourbarposition',[1 1],…

options,th2(m),L2,L3,L4,L1);

end

figure(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 0 170]) %确定XY边界值

grid

%图形加网格

xlabel('主动件转角\theta_2(度)')