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

shouce

%CalpaMEF.m
$ ?% o0 n5 R7 Y1 ?4 C- f%原始不對(duì)稱(chēng)型線計(jì)算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
  @7 @; v: n& Y' yfunction [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齒數(shù)比
9 c2 j- J- q+ aR1=(Z1/(Z1+Z2))*A; %陽(yáng)轉(zhuǎn)子節(jié)圓半徑
! _3 u0 O# n, J2 h" }7 VR2=(Z2/(Z1+Z2))*A; %陰轉(zhuǎn)子節(jié)圓半徑
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷(xiāo)齒圓弧的參數(shù)范圍 在等腰三角形中求
8 H  v3 j) a7 r# n* ]) B( F%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%銷(xiāo)齒圓弧的參數(shù)方程  GF曲線段
  f0 ^/ e8 i3 X( x; V%plot(x1,y1)
% {& z: t" ~* ~3 q) z

: }. E* j$ c- S. ?%第二曲線方程 GH  GH GH
3 Q  T- w( b/ k. L/ {3 k%b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個(gè)地方第一次弄錯(cuò)了
%t1=0;
7 h2 s7 i. d$ H6 h%x11=b1*cos(t1);y11=b1*sin(t1);%陽(yáng)轉(zhuǎn)子方程
%t1=linspace(0,t1,100);
: p0 m4 `2 D! Z8 B%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉(zhuǎn)角參數(shù)
%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%轉(zhuǎn)角參數(shù)
9 T2 m. N0 _" A+ A3 {' ]* l& i%q=linspace(q1,q2,100);
% Z9 w4 B' u9 M' s5 ]k=i+1;
%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
% r3 d) ?# N' R9 C% k; @+ \%plot(x22,y22)
$ D) _) d- W0 k" J9 x. x( p

%t21=acos((2*R2^2-R^2)/(2*R2^2));
$ k& U2 \) [6 R2 E%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ì)算線段長(zhǎng)度
%t22=acos((2*R2^2-cp^2)/(2*R2^2));
- j. j) c# s: ]; q  _%P001=(A^2+R2^2-2*A*R2)^(1/2);
: S+ z% \/ h* y* q%P002=b1;
& ]" C+ Q$ |3 [0 v7 h4 S%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號(hào)輸錯(cuò)
6 R+ ~9 Q& c8 m4 ^%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
: C8 ~: w: T+ A& Q$ ~. X! U3 _%qm=linspace(qm01,qm02,100);
+ K# G7 x6 I$ r: V* c) K% f) h%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
! R% @. _( r+ `8 d%plot(x11,y11)

# j: h7 ?+ A. o5 Y3 [+ U( I" p3 W
' q9 w; N2 e0 K6 I7 K( S4 h%第二曲線方程 EF  EF  EF
t21=acos((2*R2^2-R^2)/(2*R2^2));
- a% c: d2 Q( U9 z' x6 S3 }2 n- tp003=R2*cos(t21);    %有點(diǎn)問(wèn)題%  為什么是這個(gè)樣的,？   
p004=R2;
7 e$ j7 L3 t. C1 a0 E%PP=linspace(p003,p004,100);
qm03=1/i*(acos(k*p003/A)-t21);
( V- r6 f  c6 t  @2 C: s6 O. Zqm04=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;
7 r3 s+ e& _5 z. Dplot(x12,y12)
& q8 c& j- G" q/ h+ ]+ O

EF=[x12',y12',z12']
%save('EF.txt')
end


%CalpaMFG.m
%原始不對(duì)稱(chēng)型線計(jì)算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
, n3 h) t5 \6 h+ Afunction [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齒數(shù)比
R1=(Z1/(Z1+Z2))*A; %陽(yáng)轉(zhuǎn)子節(jié)圓半徑
R2=(Z2/(Z1+Z2))*A; %陰轉(zhuǎn)子節(jié)圓半徑
t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷(xiāo)齒圓弧的參數(shù)范圍 在等腰三角形中求
" W1 e% j/ n( T- m; It=linspace(0,t,200);
x1=R2-R*cos(t);y1=-R*sin(t);%銷(xiāo)齒圓弧的參數(shù)方程  GF曲線段
5 ?/ \. @5 e- P& ~( Nz1=0*t;
plot(x1,y1)
; y6 m; H+ w7 N5 g7 k%
& |. l; m- C7 V1 }8 X4 aFG=[x1',y1',z1']
6 @* ^: N( U  l%save('FG.txt')
end


% [ x1, y1] = CalpMFG(110, 5, 6,30)
. p% U2 |! r* }' B, a

%CalpaMGH.m
5 K4 ^" x9 D' `- J2 {8 e%原始不對(duì)稱(chēng)型線計(jì)算程序
: g# F2 B! r7 B7 z7 N+ k! D3 L8 Y0 O/ Pfunction [ x11, y11] = CalpMGH(A, Z1, Z2, R)
( o/ k5 R2 i8 m4 M0 t/ F% p3 \- _i=Z1/Z2; %齒數(shù)比
. A8 Q0 O# ~5 Z0 o- bR1=(Z1/(Z1+Z2))*A; %陽(yáng)轉(zhuǎn)子節(jié)圓半徑
R2=(Z2/(Z1+Z2))*A; %陰轉(zhuǎn)子節(jié)圓半徑
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷(xiāo)齒圓弧的參數(shù)范圍 在等腰三角形中求
4 p6 J* V+ U  @7 X& k, p5 e% t%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%銷(xiāo)齒圓弧的參數(shù)方程  GF曲線段
" T! I  b5 N. a* K8 x% Q%plot(x1,y1)


%第二曲線方程 GH  GH GH
b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個(gè)地方第一次弄錯(cuò)了
0 Q( ~0 h, ^! E9 e/ G2 ]t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%陽(yáng)轉(zhuǎn)子方程
# Q# o- b& g/ N' ^%t1=linspace(0,t1,100);
q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉(zhuǎn)角參數(shù)
%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%轉(zhuǎn)角參數(shù)
6 H% b% G" m: e7 G$ i2 r- F%q=linspace(q1,q2,100);
k=i+1;
5 Z: h1 `/ p$ y1 l1 t/ r%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
%plot(x22,y22)
+ j0 p8 f# R1 K3 Y: k
, y$ d0 X% I, i

+ \% B0 n& @9 a- }& ^$ s# o8 v%第三段曲線
%t21=acos((2*R2^2-R^2)/(2*R2^2));
* q7 ]) x2 z3 B; Gx0=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ì)算線段長(zhǎng)度
t22=acos((2*R2^2-cp^2)/(2*R2^2));
P001=(A^2+R2^2-2*A*R2)^(1/2);
8 a, R. [$ n& E" T) _P002=b1;
. I: S  f4 [* a$ Uqm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號(hào)輸錯(cuò)
3 r& |- s: G# c1 f( [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); %方程
! q; \4 B( `0 n1 K! H$ g0 \z11=0*qm;
plot(x11,y11)
2 J* Q" Z' @4 D% v/ a%
GH=[x11',y11',z11']
9 B* m; ^  r' w# E6 @' p%save('GH.txt')
/ n) Q, T1 d( B# G# I+ w8 G# Lend


4 U  {+ p8 q8 f