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

shouce

%CalpaMEF.m
%原始不對稱型線計(jì)算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
4 `2 P+ j9 b+ ifunction [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齒數(shù)比
. @) ?6 B% B  R/ A1 m2 H  aR1=(Z1/(Z1+Z2))*A; %陽轉(zhuǎn)子節(jié)圓半徑
R2=(Z2/(Z1+Z2))*A; %陰轉(zhuǎn)子節(jié)圓半徑
2 ~+ |5 P. K' g. L2 g%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數(shù)范圍 在等腰三角形中求
# E2 V1 F$ q- S% u% _& a%t=linspace(0,t,200);
' ]7 a8 D4 D* ?4 I%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
%plot(x1,y1)


* Y4 }& r# b1 v5 v%第二曲線方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個地方第一次弄錯了
% f7 S+ k3 n9 v/ \* F$ _. m) H4 {%t1=0;
) q7 p" F/ n6 D) T' u7 F%x11=b1*cos(t1);y11=b1*sin(t1);%陽轉(zhuǎn)子方程
%t1=linspace(0,t1,100);
- n( I4 N: W! B! s, q$ R%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉(zhuǎn)角參數(shù)
+ I0 p  L: Y6 s! Z; x%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%轉(zhuǎn)角參數(shù)
%q=linspace(q1,q2,100);
k=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)
5 W2 g: }) U7 J% }% ^. y2 R7 k
# u( h5 z# |5 O9 K
" t( T/ f8 P' l* H" s" E4 H%t21=acos((2*R2^2-R^2)/(2*R2^2));
%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C點(diǎn)橫坐標(biāo)   
( V7 R8 f: I9 h4 w%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C點(diǎn)縱坐標(biāo)
7 X: [, g2 k' O& I8 @1 S+ v2 _%cp=((x0-R2)^2+y0^2)^(1/2);%計(jì)算線段長度
7 c) L4 y0 w, T, }% _2 G%t22=acos((2*R2^2-cp^2)/(2*R2^2));
; W9 p9 J8 U$ C# P3 h5 ^) @%P001=(A^2+R2^2-2*A*R2)^(1/2);
7 K: h9 F- k2 n. W2 ]%P002=b1;
%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號輸錯
# ]7 @) y6 ^- ~0 R" V0 i+ @%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
2 L( `: w+ N/ e* ?%qm=linspace(qm01,qm02,100);
+ n" x' ]+ M( g) l%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
" m* q* F$ c/ V5 `& E%plot(x11,y11)

( o, e4 Z) u/ \6 `% }+ q
0 D& \! Y) Q8 y% o* n0 j) G4 Y%第二曲線方程 EF  EF  EF
t21=acos((2*R2^2-R^2)/(2*R2^2));
" `( d6 J3 g% F+ Q/ c: X9 |p003=R2*cos(t21);    %有點(diǎn)問題%  為什么是這個樣的,？   
. C" E% r3 ^0 T2 L. s/ Np004=R2;
' {4 u/ c9 ~7 I) c%PP=linspace(p003,p004,100);
! R$ K/ J6 m( `, aqm03=1/i*(acos(k*p003/A)-t21);
# S3 O0 G7 G1 r9 `qm04=1/i*(acos(k*p004/A)-t21);
qm1=linspace(qm03,qm04,100);
3 |  |1 [6 D; y. d" Z: A7 Ix12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
z12=0*qm1;
plot(x12,y12)
* x9 b0 v, t' R: U
+ ?. Z9 {6 w; r) ?
1 Y, [1 v' ~9 t  i; nEF=[x12',y12',z12']
%save('EF.txt')
6 S# Q4 P6 x" \end
; u3 K& }4 F: h2 x& I& _/ U2 p% p8 I
, W& L- w& K$ R
%CalpaMFG.m
%原始不對稱型線計(jì)算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
function [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齒數(shù)比
R1=(Z1/(Z1+Z2))*A; %陽轉(zhuǎn)子節(jié)圓半徑
' P) o: t  T! }# kR2=(Z2/(Z1+Z2))*A; %陰轉(zhuǎn)子節(jié)圓半徑
3 g+ v6 @& {3 L) \t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數(shù)范圍 在等腰三角形中求
t=linspace(0,t,200);
- i6 t0 W. Y$ M: d0 @# ax1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
z1=0*t;
* w9 o4 W  c6 a' n5 E# v) R# mplot(x1,y1)
# F3 }3 V/ A( Z& t%
  x. S' U% T- k, B9 CFG=[x1',y1',z1']
%save('FG.txt')
# h; d, Z; \. Jend


7 q' U1 p3 ?% T$ w% y) @! O1 z$ M% [ x1, y1] = CalpMFG(110, 5, 6,30)
! `" r5 i2 c! Z& A6 S
9 Q3 o6 c) V0 U0 K* Y# |; c! f
%CalpaMGH.m
. E! m2 x2 C- M8 M" i: x* o1 o! r9 u%原始不對稱型線計(jì)算程序
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
& Y9 K. w0 u  ~! i! d" ]0 _i=Z1/Z2; %齒數(shù)比
" p% z. W9 |! s6 [4 g  mR1=(Z1/(Z1+Z2))*A; %陽轉(zhuǎn)子節(jié)圓半徑
R2=(Z2/(Z1+Z2))*A; %陰轉(zhuǎn)子節(jié)圓半徑
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數(shù)范圍 在等腰三角形中求
6 C. `6 @* m" H9 e+ W%t=linspace(0,t,200);
0 X4 I8 L" d0 L6 H9 x%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
: a; g& w1 {, X3 O! q%plot(x1,y1)


%第二曲線方程 GH  GH GH
8 I" Z' d5 O4 w. m7 D3 Y8 Qb1=(R^2+R1^2+2*R*R1)^(1/2);  %這個地方第一次弄錯了
3 O# b# Y4 y! ?8 C8 M$ Nt1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%陽轉(zhuǎn)子方程
%t1=linspace(0,t1,100);
q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉(zhuǎn)角參數(shù)
6 _- {  U8 h$ u%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%轉(zhuǎn)角參數(shù)
! x9 k, }9 H* ~0 \7 g%q=linspace(q1,q2,100);
- u8 q' V) O0 O, z5 [k=i+1;
%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
% V  z% m# p$ S0 X- w0 ~( v%plot(x22,y22)



%第三段曲線
%t21=acos((2*R2^2-R^2)/(2*R2^2));
x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C點(diǎn)橫坐標(biāo)   
y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C點(diǎn)縱坐標(biāo)
cp=((x0-R2)^2+y0^2)^(1/2);%計(jì)算線段長度
t22=acos((2*R2^2-cp^2)/(2*R2^2));
& s) |7 C) K( Q0 J/ Z/ \P001=(A^2+R2^2-2*A*R2)^(1/2);
- D' e7 Y+ c9 e* cP002=b1;
qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號輸錯
3 P, t. N9 p5 U4 u# h" xqm02=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); %方程
- e6 b( E. k% {, Y  W2 q& Iz11=0*qm;
, V9 o- i* R; _4 U' cplot(x11,y11)
%
  Z1 r; {: W- X/ Q' v/ f1 VGH=[x11',y11',z11']
1 @1 j9 A; O- T9 d0 D%save('GH.txt')
end
3 [2 ?+ C9 a6 ~2 Y/ V+ M; h
6 }. V( |3 @! w; Z" t: w% K( C

, J4 |# y5 G1 @6 e! S) l. _
! u3 Q  `, g. f9 ]( I