平面连杆机构优化设计Word格式.docx

上传人:wj 文档编号:8699415 上传时间:2023-05-13 格式:DOCX 页数:6 大小:96.23KB
下载 相关 举报
平面连杆机构优化设计Word格式.docx_第1页
第1页 / 共6页
平面连杆机构优化设计Word格式.docx_第2页
第2页 / 共6页
平面连杆机构优化设计Word格式.docx_第3页
第3页 / 共6页
平面连杆机构优化设计Word格式.docx_第4页
第4页 / 共6页
平面连杆机构优化设计Word格式.docx_第5页
第5页 / 共6页
平面连杆机构优化设计Word格式.docx_第6页
第6页 / 共6页
亲,该文档总共6页,全部预览完了,如果喜欢就下载吧!
下载资源
资源描述

平面连杆机构优化设计Word格式.docx

《平面连杆机构优化设计Word格式.docx》由会员分享,可在线阅读,更多相关《平面连杆机构优化设计Word格式.docx(6页珍藏版)》请在冰点文库上搜索。

平面连杆机构优化设计Word格式.docx

~135°

范围内变化。

ψ0

l3

φ0

90°

l4

l1

l2

图1机构运动简图

二、基本思路

四杆机构的设计要求可归纳为三类,即满足预定的连杆位置要求、满足预定的运动规律要求、满足预定的轨迹要求。

本案例中,要求曲柄作等速转动时,摇杆的转角满足预定运动规律。

优化设计时,通常无精确解,一般采用数值方法得到近似解。

本案例将机构预定的运动规律与实际运动规律观测量之间的偏差最小设为目标,由此建立优化设计数学模型,并运用MATLAB优化工具箱的相关函数进行求解。

三、要点分析

优化设计数学模型的三要素包括设计变量、目标函数和约束条件。

依次确定三要素后,编写程序进行计算。

1.设计变量的确定

通常将机构中的各杆长度,以及摇杆按预定运动规律运动时,曲柄所处的初始位置角φ0列为设计变量,即

(1)

考虑到机构各杆长按比例变化时,不会改变其运动规律,因此在计算可取l1为单位长度,而其他杆长则按比例取为l1的倍数。

若曲柄的初始位置对应摇杆的右极限位置,则φ0及ψ0均为杆长的函数,即

(2)

(3)

因此,设计变量缩减为3个独立变量,即

(4)

2.目标函数的建立

以机构预定的运动规律观测量ψEi与实际运动规律观测量ψi之间的偏差平方和最小为指标来建立目标函数,即

(5)

式中,m为输入角的等分数;

ψEi为预期输出角,ψEi=ψE(φi);

ψi为实际输出角。

由图2可知:

(6)

(7)

(8)

(9)

αi

βi

ψi

φi

(a)0≤φi<

π(b)π≤φi<

图2曲柄摇杆机构的运动学关系

3.约束条件的确定

(1)曲柄摇杆机构应满足曲柄存在条件,可得

(10)

(11)

(12)

(13)

(14)

(15)

(2)连杆与摇杆的夹角应在γmin和γmax之间,即

(16)

(17)

四、具体步骤

1.选择设计变量

已知l1=100mm,l4=500mm,且φ0和ψ0不是独立参数,它们可由下式

(2)、式(3)求出,即

所以该问题只有两个独立参数l2和l3,故设计向量为

2.建立目标函数

将输入角分成30等分,并依次取30个观测点ψ1,ψ2,...,ψ30,得目标函数

式中:

3.确定约束条件

约束函数按曲柄存在条件及对传动角的限制来建立,得

4.MATLAB程序及优化结果

这是一个具有2个设计变量、7个不等式约束条件的优化设计问题。

应用MATLAB软件的优化工具箱的fmincon函数对上述优化问题求解。

(1)编写m文件Objfun.m定义目标函数。

functionf=objfun(x)

l1=100;

l4=500;

th0=acos(((100+x

(1))^2-x

(2)^2+250000)/(1000*(100+x

(1))));

ps0=acos(((100+x

(1))^2-x

(2)^2-250000)/(1000*x

(2)));

f=0;

forth=th0:

pi/2/30:

th0+pi/2

r=(10000+250000-2*100*500*cos(th))^0.5;

a=acos((r^2+x

(2)^2-x

(1)^2)/(2*r*x

(2)));

b=acos((r^2+240000)/(1000*r));

ps=pi-a-b;

pse=ps0+2/(3*pi)*(th-th0)^2;

f=f+(ps-pse)^2;

end

(2)编写m文件confun.m定义约束。

function[c,ceq]=confun(x)

c

(1)=100-x

(1);

c

(2)=100-x

(2);

c(3)=600-x

(1)-x

(2);

c(4)=x

(1)-x

(2)-400;

c(5)=x

(2)-x

(1)-400;

c(6)=x

(1)^2+x

(2)^2-1.414*x

(1)*x

(2)-160000;

c(7)=360000-x

(1)^2-x

(2)^2-1.414*x

(1)*x

(2);

ceq=[];

(3)编写m文件run.m求解计算。

x0=[400400];

options=optimset('

LargeScale'

,'

off'

);

[x,fval]=fmincon(@objfun,x0,[],[],[],[],[],[],@confun)

(4)运行m文件run.m,得最优解X*=(412.8926mm,232.2417mm),f(X*)=0.0076mm2。

五、问题拓展

满足预定运动轨迹的优化设计,要求机构在运行过程中,连杆上的某点(分析点)尽可能沿着给定的曲线运动。

设计时,连杆分析点坐标可由机构杆长和夹角表示。

以分析点的预定轨迹观测点坐标值与实际轨迹观测点坐标值之间的偏差平和最小为指标来建立目标函数,并列出传动角要求、曲柄存在条件以及杆长尺寸限制等约束条件。

展开阅读全文
相关资源
猜你喜欢
相关搜索
资源标签

当前位置:首页 > 高等教育 > 军事

copyright@ 2008-2023 冰点文库 网站版权所有

经营许可证编号:鄂ICP备19020893号-2