机械社区

 找回密码
 注册会员

QQ登录

只需一步,快速开始

搜索
查看: 2875|回复: 0

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

[复制链接]
发表于 2010-1-19 20:30:19 | 显示全部楼层 |阅读模式
我想根据一篇文献里的关系式重现齿轮圆锥曲线,前面公式及命令如下:
3 G/ f0 z: s5 [* y, tfi1=[0:0.1:2*pi];
$ @6 |* D4 c2 ac=0.16;
" N& N) T; I6 R% h3 X3 ~R=43.6;
1 d4 J4 k) q$ q* \/ vfi2=0.75*fi1+(c/3)*0.75*cos(3*fi1+3*pi/2);      %0.75=z1/z2=3/4,我直接换成0.75代入关系式了。' I7 z9 J4 E& {1 B7 \& j/ _3 ]# C# a
delta1=atan(0.75*(1-c*sin(3*fi1+3*pi/2)));
; b) v9 o3 y# s. F$ {+ g/ L1 xdelta2=acot(0.75*(1-c*sin(3*fi1+3*pi/2)));" U4 k3 m; S9 X

- J4 G% W# e( w/ l- ^3 p球面方程如下:0 v( {' ?' R' l8 T% l$ l
x1=R*(sin(delta1).*cos(fi1));. J! ~% P# k8 `8 n% w1 G
y1=R*sin(delta1).*sin(fi1);
* u8 h# r5 A1 U5 Y/ N% ez1=R*cos(delta1);

1 f, `0 q3 n4 g% Q$ I4 Er1=x1*i+y1*j+z1*k     %这里的i,j,k前的系数是节曲线上的点在三个坐标上的坐标。写到程序里面的不是这样的,下面也是一样。这样写出来大家理解直观一些。
2 ]2 X* P  O7 k" F% jx2=R*(sin(delta2).*cos(fi2));
4 O2 \: O& c+ x) f2 K0 \y2=R*sin(delta2).*sin(fi2);$ ^6 F/ e$ E2 q$ \9 \* z
z2=R*cos(delta2);
) G3 c6 P  R5 v3 u
r2=x2*i+y2*j+z2*k
9 B, z7 Z4 E# d1 ^: y; i/ f5 u
. J2 U% a: [$ f+ C8 |5 b3 g& }r1和r2的值验算出来的都是43.6,没有问题。但是我是不知道用什么样的命令得到下面的球面图形。由于本人对matlab不很了解,到这里用命令老是出错,各位英雄女侠帮帮忙。感激不尽!!!
7 L* F/ c( h* P: v; S/ a # ~( f: R( E* d5 T% K3 [
, S; r2 }: k, i- ^, g/ {

本帖子中包含更多资源

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

x
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2024-11-1 12:33 , Processed in 0.053566 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4 Licensed

© 2001-2017 Comsenz Inc.

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