牛头刨床matlab方针代码.docVIP

下载本文档

18
0
约2.72千字
约 4页
2016-07-05 发布于河南
举报
版权申诉

牛头刨床matlab方针代码.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

方案一 clear all clc n = 360 fr = 4500 l1 = 0.35 * sin(15 / 180 * pi) l3 = 0.7727407 l4 = 0.22 l6 = 0.35 slm = 0.4196 m5 = 700/9.8 g5 = 700 js3 = 1.1 g3 = 200 m3 = 200/9.8 h = 0.76 w1 = 72*2*pi/60 fi = 75 / 180 * pi fi1 = fi:0.0175:(2*pi+fi) fi1t=1:1:360 lx = sqrt(l1^2+l6^2-2*l6*l1*cos(fi1)) for i=1:1:n lxd1(i) = (l6*l1*w1*sin(fi1(i)))/(sqrt(l1^2+l6^2-2*l6*l1*cos(fi1(i)))) lxd2(i) = (l1*l6*w1^2*cos(fi1(i)))/(l1^2 - 2*cos(fi1(i))*l1*l6 + l6^2)^(1/2) - (l1^2*l6^2*w1^2*sin(fi1(i))^2)/(l1^2 - 2*cos(fi1(i))*l1*l6 + l6^2)^(3/2) end for i=1:1:n fi3(i)=atan((l1*sin(fi1(i)))/(l6-l1*cos(fi1(i)))) end for i=1:1:n fi4(i)=asin((h-l3*cos(fi3(i)))/l4) end for i=1:1:n s(i) = slm-l4*cos(fi4(i))-l3*sin(fi3(i)) end for i = 1:1:n w3(i) = (l1*w1*cos(fi1(i))-lxd1(i)*sin(fi3(i)))/(lx(i)*cos(fi3(i))) end for i = 1:1:n w4(i) = (l3*w3(i)*sin(fi3(i)))/(l4*cos(fi4(i))) end for i=1:1:n v(i)=l4*w4(i)*sin(fi4(i))-l3*w3(i)*cos(fi3(i)) end for i=1:1:n e3(i) = (-l1*w1^2*sin(fi1(i))+lx(i)*w3(i)^2*sin(fi3(i))-lxd1(i)*w3(i)*cos(fi3(i))-lxd1(i)*w3(i)*cos(fi3(i))-lxd2(i)*sin(fi3(i)))/(lx(i)*cos(fi3(i))) end for i=1:1:n e4(i) = (l4*w4(i)^2*sin(fi4(i))+l3*e3(i)*sin(fi3(i))+l3*w3(i)^2*cos(fi3(i)))/(l4*cos(fi4(i))) end for i=1:1:n a(i)=l4*e4(i)*sin(fi4(i))+l4*w4(i)^2*cos(fi4(i))-l3*e3(i)*cos(fi3(i))+l3*w3(i)^2*sin(fi3(i)) end for i=43:1:170,b=1:1:170-43 ad(b) = (a(i)-a(i-1))/(2*pi/360) end for i=1:1:n if v(i)0(s(i)0.36s(i)0.04) re45x(i) = m5*a(i)+fr else re45x(i) = m5*a(i) end end for i=1:1:n re45(i) = re45x(i)/cos(fi4(i)) re45y(i) = re45(i)*sin(abs(fi4(i))) r65(i) = g5+re45y(i) end for i=1:1:n r12(i) = (re45x(i)*l3*cos(fi3(i))+re45y(i)*l3*sin(abs(fi3(i)))+js3*e3(i)+m3*(l3/2*e3(i)*cos(fi3(i))+w3(i)^2*l3/2*sin(abs(fi3(i)))+l3/2*e3(i)*sin(abs(fi3(i)))+w3(i)^2*l3/2*cos(fi3(i)))*l3/2+g3*l3/2*sin(fi3(i)))/lx(i) end for i=1:1:n mi(i) = r12(i)*sin(fi1(i)+fi3(i)-pi/2)*l1 end for i=1:1:n rc3y(i) = re45y(i)+m3*(l3/2*e3(i)*sin(abs(fi3(i)))+w3(i)^2*l3/2*cos(