參數化計算 原始不對稱型線

shouce

%CalpaMEF.m
$ {6 M3 h4 C2 b& }1 H%原始不對稱型線計算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齒數比
1 f; a. z7 {( U' `0 FR1=(Z1/(Z1+Z2))*A; %陽轉子節圓半徑
3 A: j+ q2 P; u5 FR2=(Z2/(Z1+Z2))*A; %陰轉子節圓半徑
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數范圍 在等腰三角形中求
%t=linspace(0,t,200);
, S5 W. q' M/ s5 N3 n%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數方程  GF曲線段
% F: U+ ?$ V& p9 p; P%plot(x1,y1)


9 p, c3 s% W0 [  F1 W3 K$ g/ p%第二曲線方程 GH  GH GH
6 n) ]) Q$ {$ g% q! w; t3 l%b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個地方第一次弄錯了
%t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%陽轉子方程
3 L# G! U4 r7 Y% @& Z& ?, x% l%t1=linspace(0,t1,100);
, g! o( V/ a5 A; A3 W/ E, h3 K%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉角參數
! f% X2 [, o, @; R* |%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;
1 @& O% E: y! Z%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
%plot(x22,y22)
. `9 e1 {* W5 V, {6 v% C# _/ j

' B- g1 x/ l* `! H%t21=acos((2*R2^2-R^2)/(2*R2^2));
%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C點橫坐標   
6 ?! ?3 W7 [1 {% _; U%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C點縱坐標
%cp=((x0-R2)^2+y0^2)^(1/2);%計算線段長度
%t22=acos((2*R2^2-cp^2)/(2*R2^2));
%P001=(A^2+R2^2-2*A*R2)^(1/2);
' r2 Y- b: |( P- M: \5 a1 N%P002=b1;
%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號輸錯
6 _' [6 s- [$ j* c! i%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
$ ~: t+ Q, O& E2 f& z1 b%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)
! ]% E  E' y9 n* w

; N* N! y6 r7 y( B%第二曲線方程 EF  EF  EF
) r" Y% \( R  r8 j0 p! ct21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有點問題%  為什么是這個樣的？   
p004=R2;
%PP=linspace(p003,p004,100);
qm03=1/i*(acos(k*p003/A)-t21);
qm04=1/i*(acos(k*p004/A)-t21);
qm1=linspace(qm03,qm04,100);
x12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
z12=0*qm1;
plot(x12,y12)
9 P/ @8 Q! n, l
3 }1 B, t* e9 o; v  d' I
EF=[x12',y12',z12']
%save('EF.txt')
end
+ G$ t2 M5 P2 ^- l6 J9 ]! e
8 m0 q+ {; d7 b+ `# I
%CalpaMFG.m
% b' L6 k( H& v: Z0 G; H%原始不對稱型線計算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
function [ x1, y1] = CalpMFG(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;%銷齒圓弧的參數范圍 在等腰三角形中求
t=linspace(0,t,200);
x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數方程  GF曲線段
7 N: _8 n& E8 C0 c" B' V  L4 [z1=0*t;
) ]0 f. O+ S7 p2 K4 Dplot(x1,y1)
%
  G& q' _; R. D. }) CFG=[x1',y1',z1']
%save('FG.txt')
end

( Q6 x+ P) _* o% U* b
( s+ r+ @# c$ _2 O) `1 l" m% [ x1, y1] = CalpMFG(110, 5, 6,30)
3 N' U5 N- H& {$ p6 G2 c+ F

%CalpaMGH.m
' t8 |7 z: @( D  g  m! h/ Y%原始不對稱型線計算程序
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
7 ]2 F. y7 I+ \8 T: Oi=Z1/Z2; %齒數比
R1=(Z1/(Z1+Z2))*A; %陽轉子節圓半徑
* W" x* [4 ?( \R2=(Z2/(Z1+Z2))*A; %陰轉子節圓半徑
/ Q% w% ^# W; E1 r, Y/ n: V5 U%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數范圍 在等腰三角形中求
%t=linspace(0,t,200);
2 n, i$ z, z% a+ `( F, Z%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數方程  GF曲線段
%plot(x1,y1)
! ?4 W1 d3 C' z7 w' @  i

/ W- ^, O' m! J5 Z; \%第二曲線方程 GH  GH GH
: z& v6 }6 ]4 t/ ^9 k4 s; Gb1=(R^2+R1^2+2*R*R1)^(1/2);  %這個地方第一次弄錯了
; K; i0 g' H. ?6 A$ {# _t1=0;
! J, V7 N' Y. q" z%x11=b1*cos(t1);y11=b1*sin(t1);%陽轉子方程
%t1=linspace(0,t1,100);
- W+ Y5 d! X! a+ ^q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉角參數
%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%轉角參數
%q=linspace(q1,q2,100);
3 E. t9 D( d+ I" uk=i+1;
) J( w- g5 y& X- w%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
%plot(x22,y22)

' w! g7 R& T  i# c

3 m$ S. G9 J6 K4 k6 U+ s%第三段曲線
. G: C3 E: U, M- \%t21=acos((2*R2^2-R^2)/(2*R2^2));
) e$ @4 a% q, J- n" ?; a( yx0=A*cos(i*q1)-b1*cos(t1-k*q1);%C點橫坐標   
y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C點縱坐標
. g/ Y! c% U& x- w1 @cp=((x0-R2)^2+y0^2)^(1/2);%計算線段長度
/ M+ w- t9 f, B  @- w, `4 R9 p! ct22=acos((2*R2^2-cp^2)/(2*R2^2));
: E$ g3 U6 |" `P001=(A^2+R2^2-2*A*R2)^(1/2);
P002=b1;
3 }% m4 i" L2 z+ O0 e: N& j- q: Cqm01=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)));
5 _# f, I* f/ t+ Qqm=linspace(qm01,qm02,100);
x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
z11=0*qm;
8 {# V4 A6 `/ g$ b$ aplot(x11,y11)
%
GH=[x11',y11',z11']
0 C& Z- h: i1 ^9 \* R9 k+ @; v/ H%save('GH.txt')
) K) x' U0 P; l. F7 lend
4 t1 Z7 \: n* f6 G5 a7 W0 Q' T, _

, \6 _! y) w0 b3 T- p9 M
  @- f- Q7 n  {& d8 x; d8 f