zjg2000 发表于 2006-5-23 18:08

新手求助装载机举升机构优化问题,大虾前辈们帮我!!

那位前辈能帮我编一下这个程序:关于装载机举升机构优化问题
目标函数:举升机构的传力比要求最大:UF2=R9*R10*sin((UG3+UG-UG5)/180*pi)/(HI4*(H4+XA));
其中R9、R10、(杆长)UG3、UG5(角度)是需要优化的变量。
其中初始给定值:R9=165;R10=70;UG3=0;UG5=12;<BR>UG=UG4:0.5:UG6;
给定常数及关系式:H1=115;
H2=300;
H3=133;
Y1=210;
Y3=27;
Xj=80;
R11=80.75;
C1=sqrt(Y3^2+H3^2);
D1=atan(Y3/H3);
L=sqrt((Xj+R11+H1-C1*cos(D1+pi/4))^2+(H2-Y1+C1*sin(D1+pi/4))^2);
H4=H3-10;<BR>UG4=(acos((Y1-Y3)/L))*180/pi;
A1=H2-Y1+C1*sin(D1+pi/4);
A2=Xj+R11+H1-C1*cos(D1+pi/4);
A3=atan(A1/A2);
UG6=(A3+pi/2)*180/pi;
HI4=sqrt(R9^2+R10^2-2*R9*R10*cos(UG3+UG4-UG5));
XA=L*sin (UG/180*pi);
约束条件:
50-R9*R10*sin(UG3+UG6-UG5)/sqrt(R9^2+R10^2-2*R9*R10*cos(UG3+UG6-UG5))〈0
900-sqrt(R9^2+R10^2-2*R9*R10*cos(UG3+UG4-UG5))〈0
400-(sqrt(R9^2+R10^2-2*R9*R10*cos(UG3+UG6-UG5))-sqrt(R9^2+R10^2-2*R9*R10*cos(UG3+UG4-UG5)))〈0
可是用fmincon等函数编写,我写的程序老是有错误,实在没办法,着急啊,请各位大虾前辈帮帮我。

[ 本帖最后由 xinyuxf 于 2007-3-14 11:33 编辑 ]

frogfish 发表于 2006-5-24 21:04

回复:(zjg2000)新手求助,大虾钱 辈们帮我!!谢谢...

<P>什么错误贴一下,我的现在的机器里没装matlab </P>

zjg2000 发表于 2006-5-27 09:44

谢谢你的关心

我慢慢的调试好了,是关于我毕业设计的,不过老师要求我们用基于遗传算法的工具箱进行优化,可是我不知道这个工具箱 ,你老人家能不能指教一下

zjg2000 发表于 2006-5-28 09:11

博士给我看看是什么错误?

<P>我的目标函数是:<BR>function f=Mubiaohanshu(x,UG,H4,XA)<BR>H1=115;<BR>H2=300;<BR>H3=133;<BR>Y1=210;<BR>Y3=27;<BR>Xj=80;<BR>R11=80.75;<BR>C1=sqrt(Y3^2+H3^2);<BR>D1=atan(Y3/H3);<BR>L=sqrt((Xj+R11+H1-C1*cos(D1+pi/4))^2+(H2-Y1+C1*sin(D1+pi/4))^2);<BR>UG4=(acos((Y1-Y3)/L))*180/pi; <BR>A1=H2-Y1+C1*sin(D1+pi/4);<BR>A2=Xj+R11+H1-C1*cos(D1+pi/4);<BR>A3=atan(A1/A2);<BR>UG6=(A3+pi/2)*180/pi;<BR>UG=(UG4:0.5:UG6);<BR>H4=H3-10;<BR>XA=285.*sin (UG./180.*pi);<BR>f=x(1)*x(2)*sin((x(3)+UG-x(4))/180*pi)/((sqrt(x(1)^2+x(2)^2-2*x(1)*x(2)*cos(x(3)+UG4-x(4)))*(H4+XA)))<BR>end<BR>约束条件是:function =Yueshutiaojian(x)<BR>H1=115;<BR>H2=300;<BR>H3=133;<BR>Y1=210;<BR>Y3=27;<BR>Xj=80;<BR>R11=80.75;<BR>C1=sqrt(Y3^2+H3^2);<BR>D1=atan(Y3/H3);<BR>L=sqrt((Xj+R11+H1-C1*cos(D1+pi/4))^2+(H2-Y1+C1*sin(D1+pi/4))^2);<BR>UG4=(acos((Y1-Y3)/L))*180/pi; <BR>A1=H2-Y1+C1*sin(D1+pi/4);<BR>A2=Xj+R11+H1-C1*cos(D1+pi/4);<BR>A3=atan(A1/A2);<BR>UG6=(A3+pi/2)*180/pi;<BR>C(1)=50-x(1)*x(2)*sin(x(3)+UG6-x(4))/sqrt(x(1)^2+x(2)^2-2*x(1)*x(2)*cos(x(3)+UG6-x(4)));<BR>C(2)=900-sqrt(x(1)^2+x(2)^2-2*x(1)*x(2)*cos(x(3)+UG4-x(4)));<BR>C(3)=1400-(sqrt(x(1)^2+x(2)^2-2*x(1)*x(2)*cos(x(3)+UG6-x(4)))-sqrt(x(1)^2+x(2)^2-2*x(1)*x(2)*cos(x(3)+UG4-x(4))));<BR>Ceq=[];<BR>主函数是:x0=;<BR>LBnd=;<BR>UBnd=;<BR>=fmincon('Mubiaohanshu',x0,[],[],[],[],LBnd,UBnd,'Yueshutiaojian')<BR>R9=x(1),R10=x(2),UG3=x(3),UG5=x(4)<BR>但是结果R9 =200.2440、R10 =100.2440、UG3 = -9.7700UG5 =1.7560,超出了我的限制范围为什么?</P>

yanpi88 发表于 2007-3-14 11:12

求助

你好 我也是研究类似课题的 请问:约束条件你是怎么建立的啊?我参考了王国彪的《装载机工作装置优化设计》里的,建立了一些约束条件,采用了VB编程,但是在满足约束条件的情况下,得出的优化结构数据是不对的(要么越界,要么为0)。因此,我觉得可能是我的约束条件列错了。但是有许多是无法知道的。例如:在卸料性中的CC(最小安装距),DD(最大行程)等。如果知道,麻烦赐教。我的邮箱是:yanpi88@yahoo.com.cn,万分感谢。

[ 本帖最后由 mjhzhjg 于 2007-3-18 18:58 编辑 ]

jli 发表于 2007-3-14 12:37

可以在isight中间方便的集成,你说的详细一点,我可以做个例子给你。

mjhzhjg 发表于 2007-3-18 19:05

楼主还不尽快提供啊,你遇到了热心人

jli 发表于 2007-3-19 15:12

UG6=(A3+pi/2)*180/pi;
HI4=sqrt(R9^2+R10^2-2*R9*R10*cos(UG3+UG4-UG5));
XA=L*sin (UG/180*pi);

UG6和XA的值小于零,这个物理意义对吗?
页: [1]
查看完整版本: 新手求助装载机举升机构优化问题,大虾前辈们帮我!!