参数化计算 原始不对称型线

shouce

%CalpaMEF.m
%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
, O7 E+ a$ |- `$ C) T* ]0 s%t=linspace(0,t,200);
9 O0 o+ @" s% [5 Z% K%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)
9 s0 B$ S' v0 L: b: c

%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
%t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
, G  S, O: f  p%t1=linspace(0,t1,100);
%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
: t% Y5 n. O7 g5 _: o* @%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
3 ]2 G7 t7 y4 ?0 R- b%q=linspace(q1,q2,100);
# D' {3 G! G% s* r- T# Z9 d; p5 `; kk=i+1;
%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
%plot(x22,y22)

6 q$ F- u6 `2 x5 V! e* Y& i
1 b( D) P! o! e* C# s% I%t21=acos((2*R2^2-R^2)/(2*R2^2));
  p5 X2 f: l1 g%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
%cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
# \' i1 u7 n& y1 `%t22=acos((2*R2^2-cp^2)/(2*R2^2));
, P) b2 C) l8 J! }%P001=(A^2+R2^2-2*A*R2)^(1/2);
) D/ _7 r9 A6 T$ n0 i$ L  _7 Q%P002=b1;
, U# y/ R5 x. |%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
, X( N* `8 O  U: G' m; b- L/ f# W%qm=linspace(qm01,qm02,100);
%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
%plot(x11,y11)

3 Q  G% @% {4 t# y7 F, ]
3 N% A) z/ T+ \0 i3 ~# g%第二曲线方程 EF  EF  EF
( z  q( k3 C5 I* z7 [2 v! at21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
p004=R2;
%PP=linspace(p003,p004,100);
2 M* w1 ]& G- _4 mqm03=1/i*(acos(k*p003/A)-t21);
qm04=1/i*(acos(k*p004/A)-t21);
qm1=linspace(qm03,qm04,100);
' M8 W; L* v" r1 K- Bx12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
z12=0*qm1;
plot(x12,y12)


EF=[x12',y12',z12']
+ |9 o, x% J7 v% p  e%save('EF.txt')
end
; _; y: ~" _% E6 N

%CalpaMFG.m
%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
function [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
4 ~+ ]) y2 \" J: ~9 \7 SR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
5 C( L! I$ H: f: A" }" w# f# n2 {t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
& d& s/ m  ~* S/ t1 r: d9 Bt=linspace(0,t,200);
% r7 _0 Z8 e# o0 ]5 {x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
z1=0*t;
plot(x1,y1)
%
- }! V+ ^( t$ I% v7 P& I" ZFG=[x1',y1',z1']
%save('FG.txt')
end
4 }0 a+ ?# t2 h6 Z8 ]- b
+ U6 m; V) I* N6 [7 r5 Q* ~
% [ x1, y1] = CalpMFG(110, 5, 6,30)


%CalpaMGH.m
# e  s$ J& j" t( Y%原始不对称型线计算程序
1 M% |- M( {) K- I2 S; ~function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
# B6 ?- j2 G# N, W& N%t=linspace(0,t,200);
" I; V2 v7 E* w! X%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)

; B+ U6 Z' M6 n$ @- x- U! C) @" N
; V' f7 I. Q' K5 x) a%第二曲线方程 GH  GH GH
b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
4 t' n" u4 p% `# i! K0 _%t1=linspace(0,t1,100);
q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
8 _; Y8 d- _' p3 y5 m* A0 M4 X%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
%q=linspace(q1,q2,100);
k=i+1;
8 C! N1 d; u$ w3 E, p. M%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
+ H: a' m: k' U; z! [* I%plot(x22,y22)
' H2 |9 f: H8 _2 l$ W4 _/ ]/ i
9 U. N* X9 _. c( d" ?- ~
: {- F* d: t  u3 n
%第三段曲线
7 p% D1 F" V  E%t21=acos((2*R2^2-R^2)/(2*R2^2));
x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
( o" q$ ^" v) |3 Et22=acos((2*R2^2-cp^2)/(2*R2^2));
P001=(A^2+R2^2-2*A*R2)^(1/2);
P002=b1;
6 ~% J9 r9 H# y' @* [3 N5 a1 h' O- gqm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
6 @# M: X8 R' \# r+ vqm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
qm=linspace(qm01,qm02,100);
x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
1 A% m& O3 I& O. @6 Oz11=0*qm;
3 {8 L' n3 U. Q3 X+ s! E: R1 \plot(x11,y11)
%
: _% {1 m8 i0 J* m. B6 `! [GH=[x11',y11',z11']
$ X1 b2 l7 D+ V0 O9 z& t0 X5 I%save('GH.txt')
9 w* H4 E8 u' e# Jend

8 i+ j1 f2 e% V4 Q2 E
* y& K- y) [+ t$ m

1 n. ]& \- j0 f