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

shouce

%CalpaMEF.m
# W# P% X* h: E0 I6 @%原始不對(duì)稱型線計(jì)算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
7 }5 x7 j7 U! X' T9 |+ z7 ^function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齒數(shù)比
R1=(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ù)范圍 在等腰三角形中求
4 A. p& C. Q! T: k& H$ g$ e%t=linspace(0,t,200);
* O4 c9 z& {8 f& l; I7 q- [%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
%plot(x1,y1)


%第二曲線方程 GH  GH GH
8 \) O6 C5 Q7 K. i& G! K%b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個(gè)地方第一次弄錯(cuò)了
) A' i9 l9 @& H4 j+ x%t1=0;
2 H9 C2 N, m$ P) e; M%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ù)
%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)
" j& ~- c; d8 v1 }; Z1 G, O* e0 }- B

%t21=acos((2*R2^2-R^2)/(2*R2^2));
, c. Y% i$ y5 e4 c& U. r0 V9 a4 s, V%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)
- V3 J& q& H2 A2 n%cp=((x0-R2)^2+y0^2)^(1/2);%計(jì)算線段長度
%t22=acos((2*R2^2-cp^2)/(2*R2^2));
" T) R  R6 a; ^- E' ?+ w& l%P001=(A^2+R2^2-2*A*R2)^(1/2);
: a0 z( }% _* }" {%P002=b1;
%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號(hào)輸錯(cuò)
%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
# h/ T$ T5 V; s) f%qm=linspace(qm01,qm02,100);
%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
' x4 i7 }! @) ^( [%plot(x11,y11)
, h  s- k, ^, p3 ?

%第二曲線方程 EF  EF  EF
6 t! Q; L4 K, s0 L3 it21=acos((2*R2^2-R^2)/(2*R2^2));
8 N& q' Y# x0 q1 rp003=R2*cos(t21);    %有點(diǎn)問題%  為什么是這個(gè)樣的？   
- N# r& Q  w! c) o' a" D/ C; Tp004=R2;
%PP=linspace(p003,p004,100);
2 S3 T- J6 ^& B8 h/ E, ^qm03=1/i*(acos(k*p003/A)-t21);
qm04=1/i*(acos(k*p004/A)-t21);
1 o* ^( M/ g# z  d3 l( ~% l0 aqm1=linspace(qm03,qm04,100);
x12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
& E5 h& _: O& v! E! k8 R/ `z12=0*qm1;
plot(x12,y12)
# {3 l$ y) ?1 ^
! J: d0 V9 F. u1 w6 U0 @; L8 M( M- q
) H3 h2 l9 G+ |EF=[x12',y12',z12']
$ y" F! t1 e! c' n, T! o8 X/ C%save('EF.txt')
end
' U# j1 z0 |0 F! B: b& ~! l1 }/ X

%CalpaMFG.m
%原始不對(duì)稱型線計(jì)算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
0 N1 s. h0 L4 D) o( e# {function [ x1, y1] = CalpMFG(A, Z1, Z2, R)
5 d0 I! s( x) L3 vi=Z1/Z2; %齒數(shù)比
- `2 s1 S2 w+ ^9 r& U2 D; Y3 J* r1 cR1=(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ù)范圍 在等腰三角形中求
9 u- P7 H1 I' R0 e: ot=linspace(0,t,200);
& z3 X! z8 c" X* O* m6 gx1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
1 g% Q9 r- h3 o" M+ h4 U% [0 pz1=0*t;
plot(x1,y1)
%
6 z6 k9 F3 T/ OFG=[x1',y1',z1']
%save('FG.txt')
end
/ e9 k, G$ {* L

0 _# ^1 L* i  O% v% [ x1, y1] = CalpMFG(110, 5, 6,30)
0 X) [4 p% M: t3 J# J! q, J9 B
' D1 R4 M- ~5 M/ ]
9 s. r$ r8 f: }2 M%CalpaMGH.m
%原始不對(duì)稱型線計(jì)算程序
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
& Z5 o9 E+ F% x5 f/ Wi=Z1/Z2; %齒數(shù)比
( g9 j( B% [! o9 x- K- a+ jR1=(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ù)范圍 在等腰三角形中求
, ~4 p4 u7 q5 E7 U%t=linspace(0,t,200);
! I6 x8 E: e5 z9 i' ]& v%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
7 j2 f. `8 w  D& U$ |' S+ x) H; Q%plot(x1,y1)


%第二曲線方程 GH  GH GH
b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個(gè)地方第一次弄錯(cuò)了
* K! c# A7 N/ l5 S0 dt1=0;
5 h! ?# [) O7 M- W+ d%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ù)
%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);
. M( a1 y$ u' s' H: Gk=i+1;
& m$ A7 |1 s" i8 ~! a%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
%plot(x22,y22)


& S( \  y  ?$ `+ d0 z% k
%第三段曲線
8 D2 S2 K0 E4 H0 H# C%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));
' y- e( C# p+ K5 FP001=(A^2+R2^2-2*A*R2)^(1/2);
) @* t$ m! x! d$ U8 ^. O" N: J9 [P002=b1;
  z/ u4 Z" H9 ]) D" iqm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號(hào)輸錯(cuò)
$ u. {& S6 V7 }* oqm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
; k7 U" n( G+ ?1 N$ \9 \* V; Gqm=linspace(qm01,qm02,100);
$ l; d8 w0 ~3 {6 ]0 m+ Ex11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
& D) i! Q8 ?4 r0 @7 g, Z, I# Vz11=0*qm;
) k/ J% k" l* T: j3 zplot(x11,y11)
$ n- c& F  k+ @& U/ f' ~3 u%
* w$ n0 p" n+ N: i- e4 _" k+ jGH=[x11',y11',z11']
%save('GH.txt')
: }2 T) w: ~1 U( Mend

$ P2 l3 R9 ^2 D" ~