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

shouce

%CalpaMEF.m
4 Q% f6 \0 C" a5 m# @  N& a  t%原始不對(duì)稱型線計(jì)算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
6 F6 x% ?8 y  E0 [+ efunction [ x12, y12] = CalpMEF(A, Z1, Z2, R)
7 G8 N) q3 X1 M. Q3 `' @3 xi=Z1/Z2; %齒數(shù)比
1 K3 `  B! n/ F2 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;%銷齒圓弧的參數(shù)范圍 在等腰三角形中求
%t=linspace(0,t,200);
# B+ W1 W. S/ ]. i( C% Q8 Y%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
%plot(x1,y1)

! s: x$ o! ~- B7 J3 b
3 V8 _$ J& x/ u+ |1 l$ L%第二曲線方程 GH  GH GH
9 `) n. ^" ~1 `( Z+ v3 {%b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個(gè)地方第一次弄錯(cuò)了
( P" M$ e1 H# D; F) ~4 z%t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%陽(yáng)轉(zhuǎn)子方程
%t1=linspace(0,t1,100);
2 _, v7 q. ^8 K* y%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉(zhuǎn)角參數(shù)
' ]/ g+ ?' a. j6 m- K: B" D. [%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)

# s6 ~: ?, z: V, a% V+ K
%t21=acos((2*R2^2-R^2)/(2*R2^2));
- ~& }! T. Z% O) D* ?%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C點(diǎn)橫坐標(biāo)   
8 M# G5 t8 L9 N9 F%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)度
1 p6 l) D$ {1 H$ [. m/ }: e%t22=acos((2*R2^2-cp^2)/(2*R2^2));
6 O# `* B9 C" x%P001=(A^2+R2^2-2*A*R2)^(1/2);
%P002=b1;
* m+ V! f- p8 `$ x) }/ E6 I5 Q%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號(hào)輸錯(cuò)
- I" ~) M( R3 \; R9 w! e7 I' ~' R%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
5 M- ^/ }& q+ K. v7 W%qm=linspace(qm01,qm02,100);
: X: u& K# a! G: i# s+ M; ^9 [+ c%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
%plot(x11,y11)


2 O0 ?8 K: w; ^7 N* ^5 m%第二曲線方程 EF  EF  EF
! ]1 l! {% V$ \  x8 W4 m! St21=acos((2*R2^2-R^2)/(2*R2^2));
, w, b2 b: M* k8 x' U- @p003=R2*cos(t21);    %有點(diǎn)問(wèn)題%  為什么是這個(gè)樣的？   
p004=R2;
%PP=linspace(p003,p004,100);
qm03=1/i*(acos(k*p003/A)-t21);
4 T. l  P4 s  A1 _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); %方程
1 I3 ?. b5 X" x3 Zz12=0*qm1;
plot(x12,y12)

( K) Y0 b1 Q! I6 o0 h1 M
EF=[x12',y12',z12']
%save('EF.txt')
8 S8 \4 X/ n& _7 a+ hend

) L4 @  x4 K( Y- j
) {) I# l8 f8 K%CalpaMFG.m
( g0 L. Q; p3 g2 T, q0 T5 M%原始不對(duì)稱型線計(jì)算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
& t& g3 f( l& I+ m; V! @+ [function [ x1, y1] = CalpMFG(A, Z1, Z2, R)
: k' k# m& E8 `" _8 |: i, U$ Y: ^i=Z1/Z2; %齒數(shù)比
3 O; K' ~4 u; HR1=(Z1/(Z1+Z2))*A; %陽(yáng)轉(zhuǎn)子節(jié)圓半徑
+ v" a' Y( M8 j4 WR2=(Z2/(Z1+Z2))*A; %陰轉(zhuǎn)子節(jié)圓半徑
t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數(shù)范圍 在等腰三角形中求
! @2 o# w4 c6 _7 T" c! z; S* rt=linspace(0,t,200);
5 h" X) G0 ^$ K/ u) Tx1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
z1=0*t;
plot(x1,y1)
0 q& }% L  Q6 u* V$ }  P%
" t1 Z) j3 X* Q2 AFG=[x1',y1',z1']
  e! Y- H3 T$ U* t%save('FG.txt')
2 l& {3 l- `, A( D4 ?2 Q+ eend


% [ x1, y1] = CalpMFG(110, 5, 6,30)

5 |- n! Y! z; \0 ~' T
& U0 w5 R8 D: S5 ~: q; P' |%CalpaMGH.m
%原始不對(duì)稱型線計(jì)算程序
function [ x11, y11] = CalpMGH(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;%銷齒圓弧的參數(shù)范圍 在等腰三角形中求
6 z) G# m8 `, E/ Z6 `%t=linspace(0,t,200);
- l  }  X5 B) L8 J( k8 E1 K2 ]%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
" t1 u6 t) I3 e- T" l9 {%plot(x1,y1)


+ H5 C) ^+ t" T- s$ Z/ `%第二曲線方程 GH  GH GH
8 h; x, K4 b; d  K1 L2 ~b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個(gè)地方第一次弄錯(cuò)了
# X. D+ ^2 H( s/ q. }- G+ I; N* ut1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%陽(yáng)轉(zhuǎn)子方程
  A) i) p0 s9 m* T%t1=linspace(0,t1,100);
) V! t2 a" b0 c( yq1=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;
# H, P! B5 F" |; K%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
$ ]3 `3 i6 x( N9 T%plot(x22,y22)

& w' {8 p  r3 P

; U4 P( u0 G: [5 c%第三段曲線
%t21=acos((2*R2^2-R^2)/(2*R2^2));
- c; o$ A6 X# B7 K& Z8 \. Fx0=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)
& e" r) \) H; {% s5 g3 {2 scp=((x0-R2)^2+y0^2)^(1/2);%計(jì)算線段長(zhǎng)度
t22=acos((2*R2^2-cp^2)/(2*R2^2));
: a9 ~- @/ R: z- }. [P001=(A^2+R2^2-2*A*R2)^(1/2);
9 c9 W% B5 B9 z5 `) IP002=b1;
! _) L( d. N* L  V: p0 }qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號(hào)輸錯(cuò)
+ Q3 Y2 J/ r# o( zqm02=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); %方程
" A8 W1 X1 H, y* {z11=0*qm;
$ }. e; ?" F9 |3 Lplot(x11,y11)
/ ^& J4 w% k0 U* J%
GH=[x11',y11',z11']
: J7 K5 B1 |* Z# `2 h1 B%save('GH.txt')
end

" {' ^+ |7 t: O  T8 [. ]
& Y/ m% O0 P& m, \) u: g
8 `! L* k, T* K9 f
; r0 T3 U1 G9 Q! H) x