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

shouce

%CalpaMEF.m
# ]& ]; b! T9 ~. j) B# A( i%原始不對稱型線計算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
# H6 `4 ?& y( F* K5 E8 w  |: B# E* ufunction [ 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é)圓半徑
1 }' [. ^8 u; I& V% m& W( \2 v%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數(shù)范圍 在等腰三角形中求
%t=linspace(0,t,200);
9 N' |5 Y$ u" h+ [%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
%plot(x1,y1)
) E- D! M/ A0 s/ t! O6 n
3 \) a: b) }9 ^2 C% e, X; @" m
%第二曲線方程 GH  GH GH
- A9 X8 i9 A% F8 \7 d' r4 h7 v6 q' `%b1=(R^2+R1^2+2*R*R1)^(1/2);  %這個地方第一次弄錯了
" R) J1 W5 t/ I/ S, V9 d$ H% _4 J7 u4 C%t1=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ù)
2 v) F5 G" I/ P$ W- x/ _# c: N$ t%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%轉(zhuǎn)角參數(shù)
! A. l/ t  ?8 `, `% E%q=linspace(q1,q2,100);
k=i+1;
; v! Y+ I4 Y1 z%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
& d5 F7 v9 |9 y%plot(x22,y22)

- }  B; O" p- B' a$ }. l1 o0 a
%t21=acos((2*R2^2-R^2)/(2*R2^2));
%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C點橫坐標   
2 P. ^" o& Q# l/ F%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C點縱坐標
6 T! L2 W- m& E$ L4 s%cp=((x0-R2)^2+y0^2)^(1/2);%計算線段長度
%t22=acos((2*R2^2-cp^2)/(2*R2^2));
/ e- H3 l+ _* W& \3 I, H4 V%P001=(A^2+R2^2-2*A*R2)^(1/2);
%P002=b1;
%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號輸錯
% @9 _' H: h, `$ D: Q" ?8 f+ G%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
%qm=linspace(qm01,qm02,100);
" g, ]6 j# X& ]; q%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
8 V0 U& v" H0 R%plot(x11,y11)
, U  v& v8 D3 @. l) W

+ `( B4 G3 c' O0 P2 K) j%第二曲線方程 EF  EF  EF
! z1 J9 }% Q4 at21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有點問題%  為什么是這個樣的？   
p004=R2;
! f) I# {- c( Y. V%PP=linspace(p003,p004,100);
+ P/ e2 L! @5 uqm03=1/i*(acos(k*p003/A)-t21);
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); %方程
; Z* y9 k) A" r& s/ x2 O6 L' ^' Ez12=0*qm1;
plot(x12,y12)


EF=[x12',y12',z12']
%save('EF.txt')
8 @' ?' x2 o" x2 c/ D, f7 Qend
0 l8 u( }5 u8 u3 E9 i
1 s0 A" H" |* y8 {8 G) E1 n; j
2 q1 c( V/ E& \%CalpaMFG.m
* l; L- m: [5 N# N+ {: S%原始不對稱型線計算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
( M& z7 }4 L: m  ?' B! Efunction [ x1, y1] = CalpMFG(A, Z1, Z2, R)
- F$ Z1 V$ w+ b) D* j8 Xi=Z1/Z2; %齒數(shù)比
R1=(Z1/(Z1+Z2))*A; %陽轉(zhuǎn)子節(jié)圓半徑
R2=(Z2/(Z1+Z2))*A; %陰轉(zhuǎn)子節(jié)圓半徑
, T- t6 m' ?8 W/ {* {t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數(shù)范圍 在等腰三角形中求
; U' K- a' D3 ^t=linspace(0,t,200);
x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
. p3 ?7 q) k% w! P/ xz1=0*t;
- d& L: {, Y% ~' H; _/ Iplot(x1,y1)
%
FG=[x1',y1',z1']
%save('FG.txt')
, ]1 D6 ?, N) l  F& f* ?8 {  Lend
( o2 i* V3 V# p
1 y3 i* L& f, q
% [ x1, y1] = CalpMFG(110, 5, 6,30)
" m* g8 c$ s7 I4 K0 S

) g$ m- z  \2 p# l%CalpaMGH.m
! S1 a) X( k/ T. r%原始不對稱型線計算程序
$ }7 F3 S0 Y* }3 d) S% X2 B5 t( S9 Ufunction [ x11, y11] = CalpMGH(A, Z1, Z2, R)
i=Z1/Z2; %齒數(shù)比
R1=(Z1/(Z1+Z2))*A; %陽轉(zhuǎn)子節(jié)圓半徑
R2=(Z2/(Z1+Z2))*A; %陰轉(zhuǎn)子節(jié)圓半徑
4 S" G. q8 U' O7 a" v%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數(shù)范圍 在等腰三角形中求
4 h: Y# I4 s" V%t=linspace(0,t,200);
# q, @0 W. k2 k& s/ o7 G! A& D%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數(shù)方程  GF曲線段
+ h) L1 G8 c( T7 v# [9 i%plot(x1,y1)
0 T% Z9 h" L+ R$ ?1 t
' ]8 r, V/ A. H+ m
%第二曲線方程 GH  GH GH
( P' U# k& s- U/ l, q2 p; U  jb1=(R^2+R1^2+2*R*R1)^(1/2);  %這個地方第一次弄錯了
4 M( I/ B# X6 U3 k" S* B$ }- s7 S2 zt1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%陽轉(zhuǎn)子方程
7 V9 X, `! J1 c%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ù)
4 Q: z) j1 W* M& O: F- R%q=linspace(q1,q2,100);
3 u$ j. |; ]6 lk=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)
2 p9 E5 C' b$ w4 W  K


0 b5 w' `5 A& b3 ~1 z  i1 [%第三段曲線
%t21=acos((2*R2^2-R^2)/(2*R2^2));
4 [- t9 S  M% X% P6 yx0=A*cos(i*q1)-b1*cos(t1-k*q1);%C點橫坐標   
y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C點縱坐標
cp=((x0-R2)^2+y0^2)^(1/2);%計算線段長度
$ B8 ?( T- ~7 S! k+ T$ Vt22=acos((2*R2^2-cp^2)/(2*R2^2));
2 F3 V: A4 @8 L( }' H* zP001=(A^2+R2^2-2*A*R2)^(1/2);
. i6 c3 T/ W! @9 ]5 NP002=b1;
qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號輸錯
qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
qm=linspace(qm01,qm02,100);
- x0 D* g- \5 q) Fx11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
z11=0*qm;
plot(x11,y11)
%
, a# _6 o+ q* }$ \3 @GH=[x11',y11',z11']
+ j& Y1 ]* }: k7 S& x/ u%save('GH.txt')
end



5 h5 e6 e5 k9 ~6 y4 |" i- z  o; }