參數(shù)化計算 原始不對稱型線

shouce

%CalpaMEF.m
%原始不對稱型線計算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齒數(shù)比
R1=(Z1/(Z1+Z2))*A; %陽轉子節(jié)圓半徑
R2=(Z2/(Z1+Z2))*A; %陰轉子節(jié)圓半徑
! d2 S: e7 K3 D0 Y6 ^' C%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數(shù)范圍 在等腰三角形中求
%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
%plot(x1,y1)

6 h' W: j: _- C: }* g- u
%第二曲線方程 GH  GH GH
4 w- U0 z6 |5 X% Z4 T$ K$ t%b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個地方第一次弄錯了
1 P9 v$ b- B& q4 d8 X: B%t1=0;
9 u- \$ J2 i0 Y: K% j+ j1 W%x11=b1*cos(t1);y11=b1*sin(t1);%陽轉子方程
% k( O1 p1 _! Q% H%t1=linspace(0,t1,100);
%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉角參數(shù)
%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%轉角參數(shù)
  F6 E% y6 _) J& k%q=linspace(q1,q2,100);
k=i+1;
0 ?2 M' p* ~6 F' D8 S%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
) p3 L4 T1 K3 s% z& X7 o0 K& W, m%plot(x22,y22)

, m' A7 H: f& A  P2 M' }
, K) q% K4 S9 h( V7 w%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點縱坐標
4 ~% y  W1 Q3 ~% X7 V4 B%cp=((x0-R2)^2+y0^2)^(1/2);%計算線段長度
  u) l9 `: {2 ?# G6 D( `%t22=acos((2*R2^2-cp^2)/(2*R2^2));
3 t8 o# x1 m4 n" ?# x! R% G8 |%P001=(A^2+R2^2-2*A*R2)^(1/2);
4 @, f5 o+ Z, c' H. ~%P002=b1;
* k& h5 j1 s1 R. k/ q% s. ]%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)));
%qm=linspace(qm01,qm02,100);
%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
9 P) Y* z3 V0 d% t# `. f- b, b" @%plot(x11,y11)
% W6 v( W2 o, U" ^) o2 K- ?
- ]8 [2 R/ v  m- x
% {/ Z. C& k" s* a8 \%第二曲線方程 EF  EF  EF
* M3 Q* f; ~5 x1 Et21=acos((2*R2^2-R^2)/(2*R2^2));
# M7 N4 p6 f- `p003=R2*cos(t21);    %有點問題%  為什么是這個樣的,？   
: U* Q2 k2 R; [. b& `/ g- \p004=R2;
+ d: {1 W& w1 a" `( Q2 k, V%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;
( h# I0 Q  w8 Xplot(x12,y12)
1 b+ }+ C0 c' [# n; D+ D
4 ?( x5 }% b8 n; H. s9 l. t
5 w+ B; V7 H5 @: s7 c' V8 VEF=[x12',y12',z12']
$ n) ]2 f" y! E- C5 Z%save('EF.txt')
end

# C" |" ^, t: ~$ p% v# g
, d2 b& Q  F, Z: t0 f7 {%CalpaMFG.m
%原始不對稱型線計算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
, R0 b9 Q0 j. J- P5 @- Xfunction [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齒數(shù)比
. |2 M+ ?1 Z2 Z9 `R1=(Z1/(Z1+Z2))*A; %陽轉子節(jié)圓半徑
; @/ u0 T$ A3 B- r3 \4 d0 TR2=(Z2/(Z1+Z2))*A; %陰轉子節(jié)圓半徑
t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數(shù)范圍 在等腰三角形中求
t=linspace(0,t,200);
3 o% X1 ^4 ^5 r% v% Hx1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
z1=0*t;
: Z- X& }0 D& n0 m% P# r3 V; uplot(x1,y1)
%
) [2 N; w1 H# m: K, \FG=[x1',y1',z1']
2 p2 k. G# |# A9 ?- V%save('FG.txt')
end
8 c; e/ P( L/ a+ |( d2 G0 e2 f

7 m" w- H7 j7 q+ B1 T4 c% [ x1, y1] = CalpMFG(110, 5, 6,30)
' f" \5 `# e( h' g+ f
+ o( m2 b& E& E
%CalpaMGH.m
%原始不對稱型線計算程序
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
i=Z1/Z2; %齒數(shù)比
R1=(Z1/(Z1+Z2))*A; %陽轉子節(jié)圓半徑
8 Q4 _+ E4 {, o1 dR2=(Z2/(Z1+Z2))*A; %陰轉子節(jié)圓半徑
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數(shù)范圍 在等腰三角形中求
%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
/ {. a7 v  J+ x& f( n+ x%plot(x1,y1)

5 F! C$ f; }! r5 q
%第二曲線方程 GH  GH GH
b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個地方第一次弄錯了
t1=0;
2 n+ V7 T0 W+ ~. |; E" ]%x11=b1*cos(t1);y11=b1*sin(t1);%陽轉子方程
%t1=linspace(0,t1,100);
2 B% q7 n* r/ a7 Nq1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉角參數(shù)
, k3 n6 v1 H/ {" j" s' H2 l%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%轉角參數(shù)
%q=linspace(q1,q2,100);
/ D. p. A) n/ Z( L7 m$ ^$ s, ~% qk=i+1;
; l" c( |1 p1 g% S. q%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
: s* I+ X% m; R8 _  ?%plot(x22,y22)


- T" b! ]$ a- k( l2 M8 X
* M5 |' O$ u- ?4 [3 g; M%第三段曲線
%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點縱坐標
% U( q$ d, I. c, X$ m! L- kcp=((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;
) d" X" X! D$ L" sqm01=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)));
qm=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;
plot(x11,y11)
1 \. S9 m7 d' X) W/ [9 {%
  y' n* s& U8 ^: ?GH=[x11',y11',z11']
  n5 n+ a) {$ H: K%save('GH.txt')
end
/ R; d: ]: w9 t" d" a8 {
, I& a' h5 \  B8 [  x' m# t. F0 R