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

shouce

%CalpaMEF.m
: c2 r' h8 q8 ~3 K& a%原始不對稱型線計算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齒數比
R1=(Z1/(Z1+Z2))*A; %陽轉子節圓半徑
2 Q( M; j6 d1 Z/ G3 f. l8 C3 mR2=(Z2/(Z1+Z2))*A; %陰轉子節圓半徑
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數范圍 在等腰三角形中求
, ]: A+ F- O) C1 {6 C4 f6 G%t=linspace(0,t,200);
( E8 T/ ^& E) e% V% q: I%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數方程  GF曲線段
%plot(x1,y1)
( N  `* d9 H: g6 y4 v8 p" ]" ~3 P

%第二曲線方程 GH  GH GH
9 N" z- m( ?5 V- U2 J3 \%b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個地方第一次弄錯了
%t1=0;
- X" x; q$ C1 v5 E%x11=b1*cos(t1);y11=b1*sin(t1);%陽轉子方程
%t1=linspace(0,t1,100);
%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉角參數
; m* F- C) L7 P$ P%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%轉角參數
%q=linspace(q1,q2,100);
1 ~8 X0 [- u; u; d# \4 i3 U6 Nk=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)


%t21=acos((2*R2^2-R^2)/(2*R2^2));
5 H. o4 [5 I3 b- ?: R* }" B6 P3 ^%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C點橫坐標   
% c, ]8 v. d9 R3 Y2 P4 y% z0 s%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C點縱坐標
%cp=((x0-R2)^2+y0^2)^(1/2);%計算線段長度
; _5 S' M9 d9 A%t22=acos((2*R2^2-cp^2)/(2*R2^2));
%P001=(A^2+R2^2-2*A*R2)^(1/2);
2 F8 w% A( h0 {% N. Z%P002=b1;
%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號輸錯
" H% s" N7 d; t1 d+ g%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
2 e! M! V  j  R) 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)
6 o4 ]$ A# s8 B+ ]
$ e) p1 k' I3 K$ J
: {$ ~/ {5 ]" Q%第二曲線方程 EF  EF  EF
0 K) k& ?1 ?( A4 a- jt21=acos((2*R2^2-R^2)/(2*R2^2));
+ o6 G% D3 n! A; x7 t# A5 Qp003=R2*cos(t21);    %有點問題%  為什么是這個樣的？   
/ ~2 n( T' d$ p5 S# y% sp004=R2;
$ Z/ a) i, w- K, q: z4 Y% F%PP=linspace(p003,p004,100);
  O2 S$ `0 j0 P8 k8 L! Vqm03=1/i*(acos(k*p003/A)-t21);
, d8 p3 {; E1 G( o! @2 G+ rqm04=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;
( U: F4 q  F; h# y6 }' D, Q# }plot(x12,y12)
; v: _: Z! c$ L; k
# ]" l0 a6 W0 H$ v
. y6 y  q, }9 G, Q+ sEF=[x12',y12',z12']
%save('EF.txt')
end
$ U! U. S5 q. ]% A" |* \
' {9 `, h* X+ b  U
%CalpaMFG.m
) {* p( A: [9 m8 H%原始不對稱型線計算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
function [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齒數比
' i& N3 G2 \5 i6 b- Q  ~% [% xR1=(Z1/(Z1+Z2))*A; %陽轉子節圓半徑
R2=(Z2/(Z1+Z2))*A; %陰轉子節圓半徑
, T9 b' {. E, P- w' |t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數范圍 在等腰三角形中求
" h$ P% q3 m) i( v8 v/ x- Et=linspace(0,t,200);
, a8 ~6 {+ V8 q+ i; M* I! |. w, S3 ]x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數方程  GF曲線段
' s7 F3 [8 o+ h1 V0 N1 q  fz1=0*t;
9 I0 |4 S& c0 h3 z' i( `plot(x1,y1)
9 c4 s! v5 g( u( H# W%
FG=[x1',y1',z1']
  Q" E8 C7 j, R1 g1 [%save('FG.txt')
end
+ z2 z9 J1 O" L  F

% [ x1, y1] = CalpMFG(110, 5, 6,30)

- S( ]7 f4 x7 L3 W! V
%CalpaMGH.m
" J" i8 l1 _# Q%原始不對稱型線計算程序
0 Y6 u5 f, y  Afunction [ x11, y11] = CalpMGH(A, Z1, Z2, R)
i=Z1/Z2; %齒數比
  X" F. o3 x- \R1=(Z1/(Z1+Z2))*A; %陽轉子節圓半徑
R2=(Z2/(Z1+Z2))*A; %陰轉子節圓半徑
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數范圍 在等腰三角形中求
! H0 ]' q9 r+ p( e%t=linspace(0,t,200);
& F2 e- U. ~' Z+ k( h* c0 M%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數方程  GF曲線段
" k- a" k; M! w( W, J%plot(x1,y1)

6 r& r- X" U( G
%第二曲線方程 GH  GH GH
b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個地方第一次弄錯了
t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%陽轉子方程
%t1=linspace(0,t1,100);
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);
& _) r/ `, v/ j5 o  |' }k=i+1;
6 D9 `4 y8 {8 @" ^, P2 M1 o" e%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
%plot(x22,y22)
* S' f* ]% _$ S- t, C! M: _

& R% F7 _* r$ j1 e: F! @% n5 N
7 V* p2 I$ a- G8 c' L/ I6 }%第三段曲線
%t21=acos((2*R2^2-R^2)/(2*R2^2));
8 w! t6 ^% g' Q  y) O. Tx0=A*cos(i*q1)-b1*cos(t1-k*q1);%C點橫坐標   
; S- ~* Q0 g) A8 L/ V" u! d1 jy0=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);
P002=b1;
4 R  g9 c) S8 D* M1 p; {) eqm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號輸錯
  r' O0 K! b. h9 x; Qqm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
qm=linspace(qm01,qm02,100);
' f) U4 U  D& `; b) o# k% P9 Mx11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
% U( s- |4 N7 w* Z' F6 \1 Lz11=0*qm;
0 v9 O4 h2 H6 `1 {& d, J1 D, bplot(x11,y11)
; s# i& U6 _; l6 [0 L- V%
GH=[x11',y11',z11']
%save('GH.txt')
& S+ q9 i% [9 V- w! h: B+ h0 e* r; ~6 Nend
: v1 p9 E' e! `/ t- R. U
1 h/ A* S; L( ?: u4 H  V
: S! @1 q/ Q1 W! v2 }

( O4 r9 Z; O3 p) S  N