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

shouce

%CalpaMEF.m
%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
5 u  K! i7 L. p' Ri=Z1/Z2; %齿数比
2 J7 h) r+ f; t3 @, _8 @R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
%t=linspace(0,t,200);
' t6 h, f7 _+ [) @5 M" b%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)

1 K. c+ v0 [0 Z! e
) y, U  `1 d. t  A%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
7 [! |( |" h5 u5 X%t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
$ c9 J) ?/ ~% C3 X& \6 v%t1=linspace(0,t1,100);
%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
* }9 b3 g3 U3 V; w+ d1 n% e% N5 p$ U%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
0 x6 U2 l0 J3 Q%q=linspace(q1,q2,100);
9 f" ]0 `' E; z5 S; zk=i+1;
8 U1 y' L5 S$ t" o$ i+ A%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
" u3 k" L  X2 |0 P# Y- m%plot(x22,y22)


7 t/ w( Y; h. s+ w& C%t21=acos((2*R2^2-R^2)/(2*R2^2));
4 t$ h8 ]+ k2 G, I# m( m2 W%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
2 e& U+ ?% t. w: C# J+ U& ^: u%cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
4 F1 M3 m+ q9 m) |/ N* u$ `! Z%t22=acos((2*R2^2-cp^2)/(2*R2^2));
%P001=(A^2+R2^2-2*A*R2)^(1/2);
3 m5 y- w5 |5 v4 r8 M%P002=b1;
9 W* X; d" ?4 u* z# |6 L%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
1 Q- }3 ?+ |' k%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
) f3 @4 i4 t" L# @3 E%qm=linspace(qm01,qm02,100);
6 Y  L! }" Y+ q- ~%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
%plot(x11,y11)
" O' ]- c( i3 F9 }3 c: y9 a
0 c- X6 h) f' Y- ~1 N" n5 q0 C3 ]/ x
%第二曲线方程 EF  EF  EF
' c+ z% \& z8 L: H9 W4 Kt21=acos((2*R2^2-R^2)/(2*R2^2));
6 O# R9 a3 D0 C5 h/ H3 L; Fp003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
) N$ x+ k$ P) N4 ?4 B. ?6 Cp004=R2;
%PP=linspace(p003,p004,100);
qm03=1/i*(acos(k*p003/A)-t21);
/ C% N: J- e! P+ Hqm04=1/i*(acos(k*p004/A)-t21);
# ~+ e/ b( i5 o  u  }qm1=linspace(qm03,qm04,100);
x12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
5 C1 Q: ~7 l- T- {, y0 Yz12=0*qm1;
1 V( P4 J: i( r" o# a3 gplot(x12,y12)

& h2 g" r8 u: N+ h/ D4 l
( ?3 p. g3 x7 O$ ^* }& X0 kEF=[x12',y12',z12']
%save('EF.txt')
end


  B- ~8 S' r9 v" e%CalpaMFG.m
- J1 s/ b( y! B( D6 ?" k! `- W7 w4 T%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
# J& J* k  m$ o; b7 ~1 T- d& Pfunction [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
2 Q  u% M6 e- j; R$ o8 J  Kt=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
7 \- x2 z& N; W! Ct=linspace(0,t,200);
+ U7 B$ n! j$ J9 ex1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
z1=0*t;
plot(x1,y1)
3 E6 l$ U: D9 x8 \0 H. R5 p3 h2 [%
- }. V5 M5 h  |( N! d) eFG=[x1',y1',z1']
%save('FG.txt')
* \( }- y* [6 \) _5 Eend

! c1 z; _) y4 w* p) v# W! c
% [ x1, y1] = CalpMFG(110, 5, 6,30)

  L5 M, k, S) E+ |: n3 h/ e. ~
%CalpaMGH.m
: a* y8 C+ |, u, P%原始不对称型线计算程序
- P4 V9 y3 J' S2 H# Ofunction [ x11, y11] = CalpMGH(A, Z1, Z2, R)
& K2 J9 }. l% V& ]i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
; a  O( h  s' t, F* uR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
$ y7 w* k/ v9 r%t=linspace(0,t,200);
$ g  E5 ?. E& l5 }%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
4 v$ X( u. e+ @9 G- g%plot(x1,y1)
9 c6 O" k- T" u/ O

%第二曲线方程 GH  GH GH
3 ?; L' N8 O2 S1 ?; M# |b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
t1=0;
# T$ w! O% }! I2 e8 b/ Z& l%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
( X+ |# d% A& @7 \%t1=linspace(0,t1,100);
q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
* V" Z! c; O( W%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
%q=linspace(q1,q2,100);
4 J+ |& u6 F( u5 pk=i+1;
- {* C# T5 Q6 V1 _5 F%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
%plot(x22,y22)


, T: @  t( R& @# t7 _6 T1 J  i- d5 Q
%第三段曲线
0 G' `$ P. e1 p' e1 \%t21=acos((2*R2^2-R^2)/(2*R2^2));
* p/ {1 `5 y- I) N, k  Wx0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
* O# e7 [  e7 f" X/ F- Q  ~y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
6 t' S3 |1 E5 U$ z1 Nt22=acos((2*R2^2-cp^2)/(2*R2^2));
) }" k/ g8 x  m- y+ FP001=(A^2+R2^2-2*A*R2)^(1/2);
2 V, {  @! M( UP002=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)));
/ T2 o: t. T8 r0 q2 d" Rqm=linspace(qm01,qm02,100);
x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
z11=0*qm;
plot(x11,y11)
+ n7 G- a3 W) h% g7 y3 I%
GH=[x11',y11',z11']
, ]: H4 f# e# I# A' B%save('GH.txt')
end

( o, r6 F5 s" O8 S. ^

4 Z" s  T0 T7 Q) J" R4 g7 z