机械社区

 找回密码
 注册会员

QQ登录

只需一步,快速开始

搜索
查看: 2891|回复: 0

求助:matlab里球面坐标系下的非圆齿轮渐开线命令

[复制链接]
发表于 2010-1-19 20:30:19 | 显示全部楼层 |阅读模式
我想根据一篇文献里的关系式重现齿轮圆锥曲线,前面公式及命令如下:
; u+ U8 u, T# R9 v6 v2 e" d7 yfi1=[0:0.1:2*pi];
7 a" x/ w) f* \$ G0 b& m: ec=0.16;
8 R& l( R/ ^: b& C) e( vR=43.6;
( m( W! ~5 d- V' b+ D( Tfi2=0.75*fi1+(c/3)*0.75*cos(3*fi1+3*pi/2);      %0.75=z1/z2=3/4,我直接换成0.75代入关系式了。, Q- X5 V) n3 ?) {4 U3 @, ~
delta1=atan(0.75*(1-c*sin(3*fi1+3*pi/2)));$ V8 s: D* ~( c5 \' r" w- v' s
delta2=acot(0.75*(1-c*sin(3*fi1+3*pi/2)));0 ^# e/ H- m) [( a
1 P% d0 R7 W2 V' f9 w
球面方程如下:
/ v# Q7 M% M) cx1=R*(sin(delta1).*cos(fi1));1 L. ]+ `" Y. u) e0 `
y1=R*sin(delta1).*sin(fi1);+ n2 [; N; K: R: ^
z1=R*cos(delta1);
* F, K6 y: r  a. V2 a
r1=x1*i+y1*j+z1*k     %这里的i,j,k前的系数是节曲线上的点在三个坐标上的坐标。写到程序里面的不是这样的,下面也是一样。这样写出来大家理解直观一些。
  |" f/ q3 {1 t1 ~3 L" [8 w  nx2=R*(sin(delta2).*cos(fi2));# s9 \1 K0 `2 h
y2=R*sin(delta2).*sin(fi2);
7 x" V& j: J0 F0 t5 g2 Nz2=R*cos(delta2);

$ l( E% j9 `- F: br2=x2*i+y2*j+z2*k
  ]* z( H/ c( U3 L
. N: P' ]+ L0 D( Y7 Hr1和r2的值验算出来的都是43.6,没有问题。但是我是不知道用什么样的命令得到下面的球面图形。由于本人对matlab不很了解,到这里用命令老是出错,各位英雄女侠帮帮忙。感激不尽!!!! k( Z6 L- ?! Y- c) W" @9 n! g
9 [* B# I, B" f, y: @& ?. u9 a

1 U" v: }  D) d5 T

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?注册会员

x
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册会员

本版积分规则

小黑屋|手机版|Archiver|机械社区 ( 京ICP备10217105号-1,京ICP证050210号,浙公网安备33038202004372号 )

GMT+8, 2024-11-28 02:34 , Processed in 0.055191 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4 Licensed

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表