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

shouce

%CalpaMEF.m
) ~  e; A$ w, m" O%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
2 h2 N# t4 F& _  H( [function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
  }% ~( ]& D2 l& ^6 I9 MR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
4 R5 l; K* d( ]7 p" O. b%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
%t=linspace(0,t,200);
/ f+ c2 E  M* G  Z/ t%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
( j8 g3 W" F$ w: l%plot(x1,y1)
6 t# ~3 B' H0 J2 w% ]( o7 w2 Y7 F

%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
8 `0 T, _* q+ J# ^%t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
4 C( Y8 A1 {/ l' n%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
. A) z6 z) `) j9 {3 H  D%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
& O9 z/ Q- P. C6 u%q=linspace(q1,q2,100);
% Y  c& `8 A0 _0 L) \% Uk=i+1;
%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
3 H0 b$ i4 T+ Y! [# ]1 z%plot(x22,y22)


" R  w2 m6 T# M/ o% {  S) \2 f, [  P( k0 i( M%t21=acos((2*R2^2-R^2)/(2*R2^2));
%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
9 p2 k0 u! F* C, b/ \% p+ T%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));
5 O" }6 U1 g" `$ P! j) H1 J9 Z4 j%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); %方程
%plot(x11,y11)

9 z! ^+ I* C1 x  P4 t2 L% Q* V' S
%第二曲线方程 EF  EF  EF
t21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
p004=R2;
%PP=linspace(p003,p004,100);
, ^: z2 H1 d8 H: X% l/ R, O  rqm03=1/i*(acos(k*p003/A)-t21);
0 ^1 b, z/ m$ a7 w1 [qm04=1/i*(acos(k*p004/A)-t21);
5 ]" M+ W3 f- vqm1=linspace(qm03,qm04,100);
! I+ r( P$ a+ _5 _x12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
( ]  H. R2 I, _# y! t$ Mz12=0*qm1;
plot(x12,y12)
# r% B3 c8 d+ |5 b1 B: v" g
4 l" Z3 j/ r8 J, w) [' d
: U! {% @7 W5 D: R8 D7 r% F7 K$ l# R9 bEF=[x12',y12',z12']
5 m2 C0 i1 S8 t' u) h' L7 @; T%save('EF.txt')
# J# q) Z3 \" \0 ~8 r9 dend
* H$ ]* ^4 W1 U- x3 h; B. L& O
7 M8 V' J! h: \  o. V9 Q
%CalpaMFG.m
4 N, i! s% l' d%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
2 L$ d- ]8 E/ Q8 Y7 o) Q5 L. gfunction [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
8 Q3 {5 o$ C5 [+ RR1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
3 D( a1 t; |! g6 M, u' I( nt=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
t=linspace(0,t,200);
# d5 @# y8 f2 e& J8 ax1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
z1=0*t;
plot(x1,y1)
) E6 @( j; \6 f6 J%
/ v. y# }$ l  d- X# T+ l, `FG=[x1',y1',z1']
! L. h+ \, I) r0 F1 t/ j%save('FG.txt')
end
/ K( j& [! N  H
. a$ m  E3 h/ a+ l  c  ]6 `  q
% [ x1, y1] = CalpMFG(110, 5, 6,30)

0 q' k+ m  ]! o1 i4 t4 G
%CalpaMGH.m
%原始不对称型线计算程序
# F# n* C2 F4 w. u  [3 r6 n( W* Efunction [ x11, y11] = CalpMGH(A, Z1, Z2, R)
4 Z  k$ m- f$ o/ H' Z' l1 h! ai=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
4 Y$ E# t* Z3 H+ I! q& ]! NR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
5 o/ K( y+ l. Q% M$ q2 @& J( N* n%t=linspace(0,t,200);
: k3 N4 x3 J, [) _2 U( |& l0 O%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)
  N4 y$ t4 v1 a" m

%第二曲线方程 GH  GH GH
b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
t1=0;
8 x% g' M. O8 H9 y6 s3 X1 B9 e, o%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
  Z( R* @: j2 U) k4 }8 Cq1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
8 T3 f% ?% M% n8 m! w6 _%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)
8 n( n: C- s' \2 d6 {5 D+ D7 o& e

: r' V9 a- d5 E5 g
- r% m; C1 h5 ]1 V% ~2 ]& X6 |0 D%第三段曲线
%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点纵坐标
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);
8 u# x# c# l7 N% n8 y+ ]P002=b1;
% n* j1 ~  M% l3 m6 ]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); %方程
" G- Q* n( @1 Yz11=0*qm;
0 K5 Z( n) Y6 d' [( T2 g# X2 V) Mplot(x11,y11)
' C1 r. O& J! x, t. O3 S9 r: i* N%
0 q, y# X7 E1 w7 ]( _1 bGH=[x11',y11',z11']
2 z4 k6 q1 O6 Z% y% B%save('GH.txt')
end




# c& _+ p6 V0 s! i+ e