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

shouce

%CalpaMEF.m
%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
5 x# Z- {8 g& v0 i0 o+ ~* ~2 p- MR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
- @9 ~3 E' E5 _- f: l. O%t=linspace(0,t,200);
. I  V( r' g4 I2 E5 H3 U4 Z%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
% g3 N2 e& ]1 Y9 [' w+ @%plot(x1,y1)
7 \. F  Y' U% _, |+ S
7 T3 Y1 E# L+ f; X2 ?
%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
%t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
- N) [% Q. y# R( [  o%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));%转角参数
%q=linspace(q1,q2,100);
k=i+1;
0 K( p6 M- e7 }2 g6 B6 V" E/ K/ l%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
%plot(x22,y22)

/ P3 h! \1 b8 G9 `4 Q: T/ u) h# y( ^
%t21=acos((2*R2^2-R^2)/(2*R2^2));
+ _, C  c: Y" P' ^# v2 v%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
# ]) B, O7 ^4 I; I3 q%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
3 B5 m' `9 @  ^/ E0 b%cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
%t22=acos((2*R2^2-cp^2)/(2*R2^2));
2 X7 Q0 Q$ B" ~( P%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);
7 }. q' m/ O& Q7 g' n' l9 `0 ]%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
%plot(x11,y11)


%第二曲线方程 EF  EF  EF
8 D8 D6 u1 j5 G& V& ?4 _9 P7 Z" Mt21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
p004=R2;
7 o: ~* z- r3 }5 V%PP=linspace(p003,p004,100);
/ Z) V( j: E! W  ?: k% p/ @: bqm03=1/i*(acos(k*p003/A)-t21);
1 P4 e8 r) h/ t$ ^qm04=1/i*(acos(k*p004/A)-t21);
qm1=linspace(qm03,qm04,100);
/ {/ J5 b2 l* O# `x12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
z12=0*qm1;
plot(x12,y12)


. Q) h* T) c; P1 @: hEF=[x12',y12',z12']
  [. `6 w# [; v8 W6 S0 L/ |%save('EF.txt')
8 z1 s' I/ |( d+ _( ]end
/ I: e) B5 [1 W" E* }

) |3 i; A/ ^& r%CalpaMFG.m
2 j1 Y( \# n' {  e%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
function [ 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;%销齿圆弧的参数范围 在等腰三角形中求
t=linspace(0,t,200);
x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
z1=0*t;
plot(x1,y1)
8 r3 x# ~  S2 E( x1 Y3 q3 s) I%
FG=[x1',y1',z1']
%save('FG.txt')
( a2 `# O; r: t8 R3 T+ Z: tend
: h" p: ~' A; i! n2 r! f% O
9 o9 C- c  n0 f: P3 X
) o. h9 r# A1 l" K8 a$ h% [ x1, y1] = CalpMFG(110, 5, 6,30)
  f" ]9 F* B1 u7 l9 j. f: |
1 B4 C7 }1 r3 T0 X# {1 [
%CalpaMGH.m
% T1 I3 ]( t6 k%原始不对称型线计算程序
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
! ~+ A" V* T2 g% eR1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
4 q$ s% f9 T$ d' j%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)

0 n2 w% m. F+ }9 Q- j* Y
- L3 N$ G, M; S1 h. S& a& [! n  b2 `%第二曲线方程 GH  GH GH
) v7 X8 f3 f0 X: Y% Wb1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
- s0 F& a3 O- k" [t1=0;
+ A( ?1 Z0 f/ X* |0 ]8 Q9 M%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
6 j$ u3 H6 a$ C* c%t1=linspace(0,t1,100);
2 A9 d' G" n( ?: R: pq1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
( ?7 j6 F/ R/ s2 \/ m: ^1 ?. G! D%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);%曲线方程
  J" G+ ^( N) L+ A. }. [4 `%plot(x22,y22)

$ i2 g5 _, F! p! q

%第三段曲线
  K- \- ?3 k; G' j' W  h%t21=acos((2*R2^2-R^2)/(2*R2^2));
$ Q( i. v8 b% n7 e) D1 p% ?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);
9 c$ |1 K% C( v! ?$ |, mP002=b1;
- j; W3 v/ 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);
5 q4 H( t+ h  p$ Sx11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
/ `; l3 M9 o; k* }* I; A6 Ez11=0*qm;
5 Q  A: |) W/ Q4 Pplot(x11,y11)
%
: z# n3 L. A$ xGH=[x11',y11',z11']
7 A8 ^+ f6 }2 }%save('GH.txt')
% X! g( ~( p$ e0 \* k; L9 d: n; Mend
( e/ k5 {9 A0 @! p7 a! v
7 L& n% P8 [9 `7 V! [
) a8 T/ i( o& E6 F1 e5 n
1 G/ X/ |; t4 T5 Q5 A