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

shouce

%CalpaMEF.m
%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
! e  o- k, D7 I* l3 gfunction [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
  }6 X: @2 j' M3 q# S' k1 dR1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
0 P" z. X' k. Z" Q- ^R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
7 F0 r0 D$ \5 l0 H%t=linspace(0,t,200);
; k( R$ Y; l& }- e6 u3 m# o%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
3 \# h. Q' {0 {; j, m& Y%plot(x1,y1)
1 a' M. T+ d% o% G
/ t- S- K) V( v0 v. C/ b& O# o) M
%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
( n/ J( w' U+ \) Z* d. `7 {%t1=0;
. g- K+ r% y' a) I9 M7 }9 u: e2 M%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
& h+ Y. X& Y3 I) Z# a' K%t1=linspace(0,t1,100);
8 S0 m1 e' p0 |; J. y" Z6 x%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));%转角参数
3 V" i  ]7 }1 T+ o) B%q=linspace(q1,q2,100);
9 m: s/ ]1 ^# S. G2 o( Zk=i+1;
, E) o2 \0 g" m  A4 [- v%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
  R* v% G' Q4 j. a( t. x( E6 d! {%plot(x22,y22)
9 C  |! ^7 ^$ w8 z7 N
# _8 n; f0 W$ {. @4 O5 I7 z
%t21=acos((2*R2^2-R^2)/(2*R2^2));
' O2 P1 g3 u0 E( t" C%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
& ^5 z/ i) l- T3 Q! F/ R%cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
% g6 g: O+ G* ?6 `' g%t22=acos((2*R2^2-cp^2)/(2*R2^2));
%P001=(A^2+R2^2-2*A*R2)^(1/2);
. X3 X2 L, B7 Z7 ?; P%P002=b1;
7 ]- l! A1 a7 ~5 L2 i( m%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);
" q8 Q/ o! b/ F%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
) u3 z3 o, H4 T. u+ d1 {: c/ L4 [%plot(x11,y11)
. ]. P3 H3 q. U4 y  j: s0 Q

%第二曲线方程 EF  EF  EF
t21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
p004=R2;
%PP=linspace(p003,p004,100);
" z: f/ Y5 n* s2 R& O! H  Zqm03=1/i*(acos(k*p003/A)-t21);
qm04=1/i*(acos(k*p004/A)-t21);
5 C/ o9 H/ |# Y9 Pqm1=linspace(qm03,qm04,100);
x12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
z12=0*qm1;
' t4 Y: U8 C% K( [% t# K, T" y6 [plot(x12,y12)

+ Q8 k7 C2 Z0 ~; o
: k( ~6 ^. W; w& y# {) u$ I8 u* DEF=[x12',y12',z12']
%save('EF.txt')
8 K$ N; d& I$ _0 Tend
( H' Z" d# z5 A0 o# ~
! V; D# S: N; v' \' Q
%CalpaMFG.m
%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
2 b. c! z( {7 F8 |& Zfunction [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
: g. u) w) a% pt=linspace(0,t,200);
x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
! y! U7 L4 Q! K( a3 V2 k- B* Oz1=0*t;
% \: E- H! K) F% cplot(x1,y1)
%
7 v" F9 T2 z9 SFG=[x1',y1',z1']
%save('FG.txt')
end
. r3 a) G  ~  S6 d
$ z; u) K# Z- _+ w+ Q2 v! F
% [ x1, y1] = CalpMFG(110, 5, 6,30)

* k6 c9 l& t7 K6 h- J6 ~5 ^
%CalpaMGH.m
%原始不对称型线计算程序
8 o3 c! d0 w! ]8 l) Q2 ifunction [ x11, y11] = CalpMGH(A, Z1, Z2, R)
* V3 `7 o4 l, Q% i" H/ E5 _i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
9 H3 P; t2 F. E4 F" ?%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
%t=linspace(0,t,200);
# `% b% |. B6 F* b$ `; |4 D+ Q%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
7 R3 G: p# l! g& u$ A%plot(x1,y1)

  I  Q0 k' Z2 T- `+ _
%第二曲线方程 GH  GH GH
b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
  j" Z3 x9 q+ Q) R$ D* K$ E9 bt1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
! u8 v2 G* N* }" \, v9 f. gq1=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));%转角参数
! D4 ^- F) {: c! V6 r%q=linspace(q1,q2,100);
  Q9 L4 \* e7 U9 C' m, o$ 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)
* V: u7 ]- U/ I: V% ~
7 O) X- q( W6 n7 W8 l: y
) u( f" W; m( `2 T7 K% L' j
* g* \$ W4 R7 L4 k%第三段曲线
& W7 t+ l9 ?- a3 o%t21=acos((2*R2^2-R^2)/(2*R2^2));
# ]5 z5 M1 H7 D6 o5 rx0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
1 x, {& q8 U* c/ \" r0 b) Pcp=((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);
0 V7 M6 ~( d# _$ v4 [0 J$ }P002=b1;
qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
3 m( r" d7 `. ?9 V* ]qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
8 w  a/ c9 Y8 R* Fqm=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;
/ ]) w# f" d3 \plot(x11,y11)
%
3 J, l/ j* A  C$ ^" X" wGH=[x11',y11',z11']
1 Q+ f; l7 @! ^- i. d2 e; x+ h6 B2 |1 S%save('GH.txt')
. r" X7 ?; f7 y$ @9 eend
/ a4 R4 J, d  y7 k# l$ h4 d

! z1 b+ f2 m3 i7 s# L. w
/ S. n' X5 i& R
* ]/ q( {5 g/ v$ c, l. c9 m