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

shouce

%CalpaMEF.m
%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
1 w6 t0 M1 [0 Lfunction [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
& m6 k1 s1 |. B, Y( G4 y& D/ fR1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
% |# V. e- f! H9 oR2=(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曲线段
%plot(x1,y1)

$ ~' d/ s! k1 h9 p6 b
%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
$ D  i5 G2 G) X( W$ q4 m" T9 }%t1=0;
. T$ t' E7 Z" L  |: I) L. t& n) O%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
+ N2 d" K' l3 g" }) J%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
8 L. n0 P7 y5 ?4 ]9 j5 n%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
%q=linspace(q1,q2,100);
$ q2 Q  g6 T  L  w9 Mk=i+1;
- g5 p% [0 g5 l0 _& K0 S! [- h%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
%plot(x22,y22)
" C1 I7 |! ~4 z+ I1 W( h
" [: W! V& w/ N8 G3 N3 K: Y4 N
# i5 T: I" C7 ]9 q* t8 p0 l5 `%t21=acos((2*R2^2-R^2)/(2*R2^2));
1 E2 s- F8 o- L# U! u: i9 j4 b) Z9 i%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
; z) X% f" W+ L%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
%cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
2 B! ?3 f% b+ y: v( f/ Z+ b% N%t22=acos((2*R2^2-cp^2)/(2*R2^2));
6 I' @# R/ G5 y# b%P001=(A^2+R2^2-2*A*R2)^(1/2);
%P002=b1;
; w- F) S+ f7 X% n2 E%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);
  S  N  q$ A: T5 P2 @9 _* q  f2 g( V%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
%plot(x11,y11)


5 U/ ]* R* T" L* }9 H! f%第二曲线方程 EF  EF  EF
/ a' l0 u6 {, S& ^$ wt21=acos((2*R2^2-R^2)/(2*R2^2));
* z& Z% [3 C) @0 y# V- Z5 r; ^p003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
p004=R2;
- E& {) [7 B- P6 z0 o' f%PP=linspace(p003,p004,100);
" V. e" I5 x7 @3 W! U* Hqm03=1/i*(acos(k*p003/A)-t21);
2 [& G& ?1 ]# d3 P3 hqm04=1/i*(acos(k*p004/A)-t21);
; q3 i8 a: t, J4 gqm1=linspace(qm03,qm04,100);
x12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
3 o7 |! D8 Q9 n- P9 V5 U9 {z12=0*qm1;
plot(x12,y12)


' p2 F+ @& r' V5 d/ FEF=[x12',y12',z12']
0 M$ D  U% f2 F! G- F0 q' U%save('EF.txt')
* `6 r1 k8 j& K- P' Y7 Xend

) v  Q( \% v# @5 Z2 G0 S
, M* v6 s0 k* u/ F- c%CalpaMFG.m
%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
; [4 E8 h+ B$ ?5 e" E2 T. E0 Kfunction [ x1, y1] = CalpMFG(A, Z1, Z2, R)
' g3 Z& L, I" p& z" pi=Z1/Z2; %齿数比
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);
% s/ s! m' k: [! g8 S' Bx1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
z1=0*t;
1 O7 ^& D2 [5 R( }5 q: W5 q4 fplot(x1,y1)
%
* A1 v# u) N2 b7 [$ r4 YFG=[x1',y1',z1']
: [, b9 U5 x  z4 Q3 V" ~%save('FG.txt')
5 G; L% v/ W1 A* i$ X+ Z+ B) Pend
6 d0 G# D+ Q; b

% [ x1, y1] = CalpMFG(110, 5, 6,30)


%CalpaMGH.m
%原始不对称型线计算程序
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
* a0 b0 @+ Y4 I" a1 T  ui=Z1/Z2; %齿数比
+ b9 g( U- d5 w; G+ z/ @R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
/ k' |  R  h0 v8 Y$ jR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
; w, U& A& M5 o) s1 ]%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
%t=linspace(0,t,200);
% ^& J0 J+ _  j& N0 J( t2 ?$ l%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)
: B1 Z7 _) J+ }
8 n. N! U4 Z. T
%第二曲线方程 GH  GH GH
b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
t1=0;
: ^: v& w" F" M3 B%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
" `, E! y! F+ O; L# R3 W% wq1=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);%曲线方程
3 T/ h3 n2 |& Y3 l7 c( t%plot(x22,y22)



8 u, ?" I) k& Q5 K%第三段曲线
/ l: u% M( }- W- z% J" s%t21=acos((2*R2^2-R^2)/(2*R2^2));
0 P" v% L' b& s/ h' _  G4 ax0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
- {& b8 o. I8 E" F" V+ by0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
/ d' y; N9 d: C  Dcp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
t22=acos((2*R2^2-cp^2)/(2*R2^2));
( R$ }2 x3 G) U; jP001=(A^2+R2^2-2*A*R2)^(1/2);
7 g* H- ~: k: p1 u$ p6 JP002=b1;
qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
1 K& `) X( Z& w9 o/ oqm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
qm=linspace(qm01,qm02,100);
& e9 {  I' m, }; Z- G6 l! W: C, [6 @x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
z11=0*qm;
$ J4 X( G" r5 J2 l3 [, {5 rplot(x11,y11)
# K5 m  i& ]) [% J1 C+ B%
GH=[x11',y11',z11']
4 W% y% C) O- d%save('GH.txt')
end

/ A2 W9 w/ P0 u+ e8 n" _


+ K- z* A* f" D, I( b