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

shouce

%CalpaMEF.m
, Z3 W$ C7 M. P2 d) t& a- T%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
+ K0 D) y8 t7 Y: g& L! a; k0 K+ `i=Z1/Z2; %齿数比
6 |  F) v- M, G  }R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
5 {4 }  A2 b: g! ^" IR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
8 t1 j+ L; F8 a  t5 B3 j%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
%t=linspace(0,t,200);
* H1 l2 U% o2 {: I7 n2 z9 _6 u2 }%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
4 U  @; B, n& y/ C& s! |%plot(x1,y1)


%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
%t1=0;
" O7 d# \+ p/ f- X! E8 U! O%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
7 }+ Q. i: A; N" w6 V. N%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));%转角参数
1 h$ P  f/ i' v$ Q! u7 d%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)


- P  p2 n/ y, ^4 f2 `* |%t21=acos((2*R2^2-R^2)/(2*R2^2));
%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
3 F+ n8 n" ?7 `/ ~%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;
) Y# M- M4 K) k# e+ g% j. i%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
  k, s6 O. R2 L& [3 p! g( p%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
%qm=linspace(qm01,qm02,100);
$ U. y1 `( j; Q! Y, Y%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
%plot(x11,y11)

- _- W4 T5 p; ^  {5 G* o% ^, G
%第二曲线方程 EF  EF  EF
t21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
p004=R2;
& Y' n+ d+ L1 k2 n%PP=linspace(p003,p004,100);
3 m9 H, Y( x0 s# p: }8 T5 xqm03=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); %方程
. j! E4 T+ {3 N, zz12=0*qm1;
plot(x12,y12)
. j$ m2 D; H  e) m

EF=[x12',y12',z12']
+ y: o7 S* c; Y, s. x. S% e2 w' I%save('EF.txt')
* Z" j! h" ^" v% Gend

& {, r; p+ j- }
) L$ T9 M0 i( j; O%CalpaMFG.m
%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
function [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
* y) i, _8 O0 d/ W; L! nR1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
( U+ ]( H/ ^1 A; R0 a! Wt=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
& a3 ^- |2 t7 u9 ?& t: St=linspace(0,t,200);
x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
z1=0*t;
0 ^& n+ ~3 V6 Q6 [# w: {; Pplot(x1,y1)
%
FG=[x1',y1',z1']
/ g9 a. W# `  u7 q  g%save('FG.txt')
end
9 s$ Q" y3 }) Q. ~: B$ t( X# N5 @

; b6 \) }; U9 d8 q, p7 l& c' P9 Q& Y% [ x1, y1] = CalpMFG(110, 5, 6,30)
2 R3 G) @) r/ y1 i  n8 K& S3 p
- o8 \% b: D# o7 J5 I; R
0 O9 ~9 i9 e- v$ @" J# Z9 h. d%CalpaMGH.m
%原始不对称型线计算程序
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
! P5 X6 A: V( A' q  V0 y  mi=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
, A; `; y: s9 v! f" A: c& j%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曲线段
- o- @$ F5 x6 a$ W! f3 {& I. _%plot(x1,y1)
  f/ m  v3 O1 |) I1 M, t

" G% T; N" @* k1 u/ ]  A%第二曲线方程 GH  GH GH
6 v3 n6 k! a* Y. D; J* Nb1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
t1=0;
! n3 w% m& Q9 _%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
: [8 W# r" X5 k& Z6 x9 J2 yq1=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));%转角参数
! |& F2 E. E2 G* r6 i%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);%曲线方程
: }+ {, K* p2 d; w2 q0 [7 I# R%plot(x22,y22)



9 F# Q8 \) T: p; `# A8 ?3 x%第三段曲线
1 w% Q. t# ]9 W3 }$ O%t21=acos((2*R2^2-R^2)/(2*R2^2));
x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
* U0 J( v8 V2 M8 u: R  C3 {0 Jy0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
5 z  f( n1 ]1 Bt22=acos((2*R2^2-cp^2)/(2*R2^2));
* g5 J6 A- ~) {" K" k( VP001=(A^2+R2^2-2*A*R2)^(1/2);
3 A. J( O2 E: Z4 AP002=b1;
1 P  G+ b, T0 ]- Uqm01=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);
: @& ?2 h$ M% U! W5 C$ w, Ux11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
z11=0*qm;
plot(x11,y11)
%
GH=[x11',y11',z11']
) e: d3 @5 h( `%save('GH.txt')
end
, L) \, ?3 V$ W8 W  Z6 v4 J