|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
l0=177.8;R=175;l=482.6;%xitamax=40;%xc=R*l/xitamax
xc0=1000;yc0=200;xg0=900;yg0=500;M=1100;g=9.8;
xd=800;yd=150;ad=900;ac=850;ab=400;ye=-150;xe=900;
for xita=40*pi/180%0:5*pi/180:35*pi/180
xc=xc0*cos(xita)-yc0*sin(xita);
yc=xc0*sin(xita)+yc0*cos(xita);
xg=xg0*cos(xita)-yg0*sin(xita);
yg=xg0*sin(xita)+yg0*cos(xita);
syms xa ya;
[xa,ya]=solve('(xa-xd)^2+(ya-yd)^2-ad^2=0','(xa-xc)^2+(xa-yc)^2-ac^2=0');
%[eval('(xa-xd).^2+(ya-yd).^2-ad.^2'),eval('(xa-xc).^2+(xa-yc).^2-ac.^2')]
syms xb yb;
[xb yb]=solve('(xb-xa)^2+(yb-ya)^2-ab^2=0','(xa-xc)^2+(ya-yc)^2-ac^2=0');
syms xf yf;
[xf,yf]=solve('(yd-ya)*xf+(xa-xd)*yf+ya*(xd-xa)-(yd-ya)*xa=0',...
'(ye-yb)*xf+yf*(xb-xe)+yb*(xe-xb)-xb*(ye-yb)=0 ');
Dofc=abs(yf*(xc-xf)-xf*(yc-yf))/sqrt((yf-yc)^2+(xf-xc)^2);
Wxg=M*g*yg;
Ffc=Wxg/Dofc
Dabe=abs(xa.*(ye-yb)+ya.*(xb-xe)+yb.*(xe-xb)-xb.*(ye-yb))./sqrt((yb-ye).^2+(xb-xe).^2);
Dafc=abs(xa.*(yc-yf)+ya.*(xf-xc)+ye.*(xc-xf)-xf.*(yc-yf))./sqrt((yf-ya).^2+(xf-xc).^2);
Fbe=Ffc*Dafc/Dabe;
end
以上是我的程序,我想算出Fbe和Ffc这两个力来,但是就没得我想要的结果,这是怎么回事啊?
高手来帮帮忙啊!
[ 本帖最后由 eight 于 2008-4-14 20:45 编辑 ] |
|