准双曲面小轮的齿面方程出问题了，我就是按照李特文的书编写的，帮我看看谢谢

飞龙在天88866

syms a t pi  rt f g j q e01 xb1 m1 x1  p  io1   s1  n    r1  s  b c d
pi=3.1415926
a=10*pi/180
rt= 96.52 ;       %sh(407)
f=(15+21.4723)*pi/180   ;     %sh(396)刀转角
g=350.111*pi/180  ;      %sh(417)刀倾角
q=96.444*pi/180    ;       %sh(431)小轮精切刀位
e01=40.2217   ;    %sh(430)
xb1=9.0661   ;       %sh(429)
m1= 357*pi/180 ;   %sh(402)齿陪安装角
x1=-50.8138    ;     %sh(425)%ffffffff安装距
    io1=  3.7051;        %sh(436)
    s1=  102.902 ;          %sh(4;16)
r1=[];
for p=0:0.1:0.8;
  for t= 0:pi/30:2*pi;
syms s real
rr1=[(rt+s*sin(a))*cos(t);((rt+s*sin(a)))*sin(t);-s*cos(a);1];
nt=[-cos(a)*cos(t);-cos(a)*sin(t);-sin(a);0];
mbt=[cos(f) 0 sin(f) 0;0 1 0 0;-sin(f) 0 cos(f) 0;0 0 0 1];
mcb=[-sin(g) -cos(g) 0  s1;cos(g) -sin(g) 0 0;0 0 1 0;0 0 0 1];
moc=[cos(q) sin(q) 0 0;-sin(q) cos(q) 0 0;0 0 1 0;0 0 0 1];
ro=moc*mcb*mbt*rr1;
no=moc*mcb*mbt*nt;
wc=-io1.*[0;0;1;1];%%%%%%%%%%%%%%%%%
wp=-[cos(m1);0;sin(m1);1];
ooon=[0;-e01;xb1;1];%%%%%%%%%%%%%%%%%%
wc1=wc(1:3,:);
wp1=wp(1:3,:);
no1=no(1:3,:);
ooon1=ooon(1:3,:);
ro1=ro(1:3,:);
no11=no1.';
s2=no11*(cross((wc1-wp1),ro1)+cross(ooon1,wp1));
s=solve(s2,s);   
s=vpa(s,5);
mnm=[1 0 0 0;0 1 0 e01;0 0 1 -xb1;0 0 0 1];
mqn=[cos(m1) 0 sin(m1) -x1;0 1 0 0;-sin(m1) 0 cos(m1) 0;0 0 0 1];
mpq=[1 0 0 0;0 cos(p) -sin(p) 0;0 sin(p) cos(p) 0;0 0 0 1];
rrr1=[((rt+s*sin(a)))*cos(t) ;((rt+s*sin(a)))*sin(t); -s*cos(a); 1];
r_1=mpq*mqn*mnm*moc*mcb*mbt*rrr1;
r_1=vpa(r_1,5);
r_1=double(r_1);
r1=[r1 r_1]
   end
end
  r1=r1';
  r1(:,4)=[];
xlswrite('1234.xls',r1,'Sheet1','A1');

q364228550

你不能照抄书呀，，得根据自己的实际情况做调整才行。。

mcray23

楼主用的是matlab写的吗

变量不用申明

meng198599

得根据自己的实际情况做调整才行

飞龙在天88866

q364228550 发表于 2012-3-12 11:11

                               
登录/注册后可看大图

你不能照抄书呀，，得根据自己的实际情况做调整才行。。

我都实验N次了     哎郁闷呀

飞龙在天88866

meng198599 发表于 2012-3-12 19:00

                               
登录/注册后可看大图

得根据自己的实际情况做调整才行

我都做N次了   能帮帮我不？谢谢哈

nieshaowu

我写了一个齿面求解程序文件（弧齿的），你看一下，可以参考此格式，具体的程序代码要根据自己情况编写，由于涉密，部分程序段没有给出，望见谅，可参考有关资料自己思考一下。另外，齿面求解需要解非线性方程组，初值选取很关键。 ToothFlank_Qiujie.rar (2.48 KB, 下载次数: 106)

2012-3-13 13:01 上传

点击文件名下载附件

飞龙在天88866

nieshaowu 发表于 2012-3-13 13:01

                               
登录/注册后可看大图

我写了一个齿面求解程序文件（弧齿的），你看一下，可以参考此格式，具体的程序代码要根据自己情况编写，由 ...

您好!可以告诉您我的QQ吗？我就是非线性方程组不会解，现在也挺迷茫的，小轮做了一个月了也没有什么进展。谢谢您

guyueguangda

楼主您好！能加您好有友么？我想和请教您下齿面方程的问题？还有您看的那本书全名叫什么啊？我的QQ375449107

天gao云淡

我QQ305303790，也想学习学习，那本书叫什么名字？