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

shouce

%CalpaMEF.m
%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
& ?7 @& w5 Z- O3 i1 Dfunction [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
) M, f: i7 @. [) U% sR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
# t) o) A2 ^7 c) I8 s: i" x%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)

  X: ]0 L' e0 n9 ^7 {( \: K3 }
( X  f2 e& _/ h- h: l$ h# o. V%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
%t1=0;
+ O( B' o, A; ?& R9 o- _%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
: k2 l+ F+ w5 v5 a%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));%转角参数
%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)
- B" t8 j: o- j, I8 J7 f+ u
# g' ]+ y/ x# |& x& z* e
%t21=acos((2*R2^2-R^2)/(2*R2^2));
%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
1 O7 m9 I2 H7 ], h% a' P%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
7 o: U, L7 Y& g5 t( [; J- W/ z%cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
%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)));
5 B- s: Z2 V, x. p%qm=linspace(qm01,qm02,100);
& @- o/ k% y9 f1 x8 V" v%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
! y9 q" ?; y) B+ g4 j: n%plot(x11,y11)
0 S2 d% n1 z9 j' a
# ~9 s: h3 R0 }( M: H6 H
%第二曲线方程 EF  EF  EF
t21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
p004=R2;
%PP=linspace(p003,p004,100);
qm03=1/i*(acos(k*p003/A)-t21);
qm04=1/i*(acos(k*p004/A)-t21);
( T; U: R. Z6 e: E# R1 uqm1=linspace(qm03,qm04,100);
x12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
5 n+ R* w* ?- U- E& ]z12=0*qm1;
plot(x12,y12)

  l, _% ?' k3 x" r
EF=[x12',y12',z12']
%save('EF.txt')
8 n& G3 S, P+ @( s+ \) \+ ]end

! `8 F- Y" F3 w4 d
%CalpaMFG.m
! o# O3 r3 C' b: n1 n" D1 M%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
function [ x1, y1] = CalpMFG(A, Z1, Z2, R)
- v3 [$ S# d2 t+ V; di=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
& E. N  `3 g. n& Xt=linspace(0,t,200);
. ]* C1 v# I+ V; q* |/ @x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
4 a  n5 g3 K. ~$ _# Z9 a( O  Pz1=0*t;
% g- w  X+ ?/ Q6 x( x& zplot(x1,y1)
7 X2 E) h; {+ y! H5 w, z' ]%
FG=[x1',y1',z1']
%save('FG.txt')
* g7 k1 `, J1 V6 r" y( Kend


% [ x1, y1] = CalpMFG(110, 5, 6,30)
/ t8 N! o( q  v' f

%CalpaMGH.m
3 g4 b, I" T/ `$ B%原始不对称型线计算程序
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
! \) X$ k' ^9 ~7 M/ a. Q# R0 ^/ Z' y% [( ei=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
- ], k/ f2 f6 h7 b) P1 J  ~) G& P! uR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
; I3 q3 _0 b4 h& o2 A* S%plot(x1,y1)

& \! r0 S" A1 e* ^
%第二曲线方程 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));%转角参数
) @" }+ T/ f8 }( {6 r& f/ i- m- R%q=linspace(q1,q2,100);
k=i+1;
* L: H1 Z# X6 F6 F$ H8 a%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
3 T/ l8 Z# A) ^6 l/ p# M+ v%plot(x22,y22)
6 H3 _. L: i3 K7 h1 B, `
+ U  J, R9 Z. {: E( B/ T2 T$ {$ B

%第三段曲线
" H6 N0 B3 Q  C4 m%t21=acos((2*R2^2-R^2)/(2*R2^2));
x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
' H6 z& U) S. A+ gy0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
- Q" v) E% d  _/ E9 wt22=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);
  ^4 R! ^" z; j2 qx11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
4 ^  n9 O- a; V0 o/ G( ?. D6 E+ I4 az11=0*qm;
plot(x11,y11)
%
GH=[x11',y11',z11']
2 q; M/ A3 z/ W+ p%save('GH.txt')
end
& Z- y% h# ^4 S; t