参数化计算 原始不对称型线

shouce

%CalpaMEF.m
) o8 m) U! G, W) V' s5 Q%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
6 f5 N' }, O& ]; ~. Pi=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
6 g  Y: P3 B) M# KR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
! f6 k5 \# i4 f3 n2 E: u+ @; Y%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
) P5 z- W) Q& w2 B) {9 J2 p* m%plot(x1,y1)

3 v" e- }: D9 ~. n1 s, I
5 a7 V* L3 I: K%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
%t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
' x' }6 G7 c, O/ v+ h5 `& |7 \1 L%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)
+ n) |$ j  |: w+ i* L. i

%t21=acos((2*R2^2-R^2)/(2*R2^2));
  X/ j0 x! S4 U! D, J6 l%x0=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);%计算线段长度
%t22=acos((2*R2^2-cp^2)/(2*R2^2));
4 ], H' D6 e, K  y! \  ?6 k2 F%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))); %第一次在这儿括号输错
- |* l# Y$ a2 y' N. u# Q%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); %方程
0 [  y% G( S% s2 f5 X%plot(x11,y11)
4 S: b& N- C1 q% D2 J& d

2 f+ C6 ~; _) ]. v9 Q3 D9 K%第二曲线方程 EF  EF  EF
t21=acos((2*R2^2-R^2)/(2*R2^2));
$ K' K; R- c: X% op003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
. J+ Z( d7 y, `, S* Gp004=R2;
8 r8 |8 q, h8 L5 n# A( y; L%PP=linspace(p003,p004,100);
9 R% O- G' J5 d; f2 bqm03=1/i*(acos(k*p003/A)-t21);
( g  N* t: K5 G9 w. lqm04=1/i*(acos(k*p004/A)-t21);
qm1=linspace(qm03,qm04,100);
5 U+ `  \8 N- s* X2 w, Jx12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
1 u! f3 J: U- \; |8 K! f1 mz12=0*qm1;
plot(x12,y12)

. [3 Z2 S0 k# ]* k8 ?2 e+ Z
+ C! K* e1 {+ u0 bEF=[x12',y12',z12']
5 q2 T1 Q" J6 f! I, G%save('EF.txt')
0 K7 M# ~* H8 {& Rend
' t- o+ [6 d* G4 Y- W4 j% `

$ r7 {( Z+ K6 K& h# B/ p/ `%CalpaMFG.m
%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
" e3 Z) r1 j" q2 dfunction [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
8 n% J( V2 e/ h- e9 ]9 N7 |R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
2 T" x' n4 m. j' Z4 j% UR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
* K# M/ P8 I5 {' x- _4 b0 I3 jt=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
/ P3 k. Q! F7 t( f. N0 yt=linspace(0,t,200);
) t0 Z. u; B! a6 I  i( lx1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
0 c1 _8 X* t# K3 I8 `1 E, k: Hz1=0*t;
: W0 a& I3 q+ i! h; vplot(x1,y1)
% R  t) L1 t% k( F9 e% `# X%
FG=[x1',y1',z1']
%save('FG.txt')
end
7 ^) c! y4 v3 f

- r# ^& O( i& Y2 k" ?9 i% [ x1, y1] = CalpMFG(110, 5, 6,30)


%CalpaMGH.m
' U( x  h+ B; J% [, ]) m; l%原始不对称型线计算程序
" R( [/ M: g- z: U1 a" e4 @5 A; j- yfunction [ x11, y11] = CalpMGH(A, Z1, Z2, R)
5 m2 z: L+ `7 t0 v" m; `- F0 ni=Z1/Z2; %齿数比
" s' b' K* U9 S4 m$ ?R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
$ x' N  G6 W1 X) @! }; y1 x%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
# ~/ E5 `, K$ {( s: N%plot(x1,y1)


/ H* T6 W8 c/ F( F& J: h0 s5 K! \%第二曲线方程 GH  GH GH
" x4 `; m. ]- t8 s# q) [' b) qb1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
) |8 v: O. o) T9 wt1=0;
3 n$ G0 b& l" Z2 H8 ]%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
9 B9 A# t% z" rq1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
9 z$ ^* M2 R* M9 t" D( |%q=linspace(q1,q2,100);
8 E6 A4 ^6 k5 m8 Y+ hk=i+1;
%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
+ Z2 H) g' b% s. S  ]2 D- l7 v%plot(x22,y22)


. `" a  x' e' F  f
/ H$ ^3 t' ]) ?0 @/ x%第三段曲线
%t21=acos((2*R2^2-R^2)/(2*R2^2));
/ }3 u+ _3 o* ^; h9 }# nx0=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);%计算线段长度
3 [' }0 B- _" P5 ^t22=acos((2*R2^2-cp^2)/(2*R2^2));
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))); %第一次在这儿括号输错
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); %方程
. j" w' A0 ?/ {7 Nz11=0*qm;
( X: \/ C( i5 U, c9 j0 Splot(x11,y11)
%
GH=[x11',y11',z11']
%save('GH.txt')
end

  B" P4 u/ d& y3 d( F0 g
  |* `3 @: G5 o. _( F8 P3 k9 V; B
+ n7 T5 p6 Y) u6 C, v$ R9 q5 I& E" @
+ S3 z( H0 z" K5 e