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

shouce

%CalpaMEF.m
# k9 _: s) i( z4 u; ~%原始不對(duì)稱型線計(jì)算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
6 s; M) O- u8 J: p) F1 H& M2 Y6 `function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齒數(shù)比
, h" Q! r4 i. B* N9 d* R9 XR1=(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;%銷齒圓弧的參數(shù)范圍 在等腰三角形中求
%t=linspace(0,t,200);
3 I) k+ q  \% K  p0 [%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
%plot(x1,y1)

  ^* s' D% h% l; g8 V
%第二曲線方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個(gè)地方第一次弄錯(cuò)了
* w' a$ |( m% f0 v%t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%陽(yáng)轉(zhuǎn)子方程
3 M. }, {( ]6 ]% F' b9 l%t1=linspace(0,t1,100);
%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉(zhuǎn)角參數(shù)
5 n9 e6 T2 d4 Z  x5 [%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%轉(zhuǎn)角參數(shù)
+ A8 A; b4 Z, A%q=linspace(q1,q2,100);
k=i+1;
4 S4 w1 c: e: {# U3 [2 _1 I/ ?* Q8 \%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
+ `) G5 T* j) H6 f7 m- X2 N. _%plot(x22,y22)

- u' [$ l, d/ I; |; ]
%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)   
* B+ {, \1 Q* D: r, k  Q' ]% O! x4 z%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C點(diǎn)縱坐標(biāo)
& M' u7 S% U# j/ O5 s) u%cp=((x0-R2)^2+y0^2)^(1/2);%計(jì)算線段長(zhǎng)度
# Z% [/ ]/ X* C( ~3 t4 d- s%t22=acos((2*R2^2-cp^2)/(2*R2^2));
%P001=(A^2+R2^2-2*A*R2)^(1/2);
%P002=b1;
1 ?! g3 x) j+ A4 h& R# T6 U%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號(hào)輸錯(cuò)
' k) V: A1 @- T2 Z$ q9 m6 W%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
%qm=linspace(qm01,qm02,100);
, m; E+ m) \0 u3 D%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
%plot(x11,y11)
: K3 ~8 j; {# H% z4 _4 `+ F

%第二曲線方程 EF  EF  EF
  n; \. t3 m% n# k9 bt21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有點(diǎn)問題%  為什么是這個(gè)樣的？   
p004=R2;
%PP=linspace(p003,p004,100);
qm03=1/i*(acos(k*p003/A)-t21);
qm04=1/i*(acos(k*p004/A)-t21);
/ }) E) H+ P' Xqm1=linspace(qm03,qm04,100);
# ?! L) o/ I8 z# a7 Ux12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
/ ^5 q. y! A) n6 uz12=0*qm1;
plot(x12,y12)


. O7 `/ w' W: M8 L9 i* FEF=[x12',y12',z12']
%save('EF.txt')
; l& v: H1 Z4 B; Xend


( z# I; }  Q+ r%CalpaMFG.m
%原始不對(duì)稱型線計(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; %陽(yáng)轉(zhuǎn)子節(jié)圓半徑
R2=(Z2/(Z1+Z2))*A; %陰轉(zhuǎn)子節(jié)圓半徑
- X2 V- z' ]' _  {t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數(shù)范圍 在等腰三角形中求
+ T. Q# `; F( p3 h# i8 dt=linspace(0,t,200);
x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
- m; {/ f1 h+ H) i; d% kz1=0*t;
, ^7 f" s( h, @* f# Oplot(x1,y1)
%
3 r, k9 I- \4 [  ^) ?FG=[x1',y1',z1']
%save('FG.txt')
end
: \; A, G; j, V" w
* k6 A, B: f( W( B) U2 [
% [ x1, y1] = CalpMFG(110, 5, 6,30)
6 A# v) R2 k/ E/ |1 M4 c' \

%CalpaMGH.m
%原始不對(duì)稱型線計(jì)算程序
* j' t+ d3 u% n( U: [9 @) I) a- F! yfunction [ x11, y11] = CalpMGH(A, Z1, Z2, R)
$ p* @  f+ x1 X9 qi=Z1/Z2; %齒數(shù)比
R1=(Z1/(Z1+Z2))*A; %陽(yáng)轉(zhuǎn)子節(jié)圓半徑
, X' J! w: @# A$ R: p2 tR2=(Z2/(Z1+Z2))*A; %陰轉(zhuǎn)子節(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曲線段
%plot(x1,y1)

. i* X' w2 h5 s, O2 x; K! _
%第二曲線方程 GH  GH GH
+ L  e' A7 Z9 @4 c8 s0 Q; ]/ zb1=(R^2+R1^2+2*R*R1)^(1/2);  %這個(gè)地方第一次弄錯(cuò)了
t1=0;
1 u1 S. r0 K7 r( \- x%x11=b1*cos(t1);y11=b1*sin(t1);%陽(yáng)轉(zhuǎn)子方程
& z& X( N: T9 R%t1=linspace(0,t1,100);
q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉(zhuǎn)角參數(shù)
$ M& m0 Y; s: |0 Y9 U1 _%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%轉(zhuǎn)角參數(shù)
* x7 q% `6 k7 `( O, v%q=linspace(q1,q2,100);
! e7 X0 r5 u8 v* c! gk=i+1;
3 f( \2 C7 i: O( S: N%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
7 v. H  s, H: a* F& i% a9 V%plot(x22,y22)


  \! K. g, N( s# V, `) n  x
%第三段曲線
%t21=acos((2*R2^2-R^2)/(2*R2^2));
: |* s9 }. p5 J. [, Ox0=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));
4 W$ t) ~/ C- ]! [. S0 M5 U8 CP001=(A^2+R2^2-2*A*R2)^(1/2);
P002=b1;
7 M; d7 n2 M: ^6 p; Q. ^5 X- zqm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號(hào)輸錯(cuò)
# r+ x8 B* w& T! x% ~2 L3 Xqm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
0 \2 n- ?4 q4 H3 b! [qm=linspace(qm01,qm02,100);
7 W5 \, z; [" G8 Px11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
z11=0*qm;
plot(x11,y11)
%
GH=[x11',y11',z11']
%save('GH.txt')
end
3 x9 }' F9 t1 V; u


7 K' I' @. o4 _% y( n
. U, W4 x- H& L( I