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

shouce

%CalpaMEF.m
%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
: v2 r& |1 M; E; T8 ^$ Y8 @R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
: r" U$ y& R8 n6 x; y) `! L6 JR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
%t=linspace(0,t,200);
. S0 \  ?: V7 X7 R8 l%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
$ f1 B* A, x. U8 P! t) M3 b, m4 c%plot(x1,y1)
* N2 h0 U& P0 U0 U

  I' S" ?* T6 C, T0 ]" p7 q%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
%t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
; {5 S5 G( M( _/ {, \% f%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
; G) u2 F9 n7 ^8 z# t8 c' Q%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
8 C) e7 a% I! G3 Y: l5 I( W( f5 O%q=linspace(q1,q2,100);
6 t+ X, D, m  p  Hk=i+1;
. l# ]9 X# d1 w3 e; C; j/ o: i) j%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
%plot(x22,y22)
; p9 k! i0 Z& v
8 u+ m3 {: O, e- G+ M
6 }3 q5 \0 o- O% f! H' M/ W) I%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点纵坐标
& h) {4 O- X7 u1 S: M4 h' T7 Z%cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
7 b1 G: R3 k1 j%t22=acos((2*R2^2-cp^2)/(2*R2^2));
' @- T9 z8 s- m%P001=(A^2+R2^2-2*A*R2)^(1/2);
%P002=b1;
5 z7 t8 A4 i3 q* d) d# F. u%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
- f* c) L+ Q9 Z3 Q2 u) v- j3 h%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
%qm=linspace(qm01,qm02,100);
4 H8 D6 {" C0 a1 l  F) b%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
! r7 F) A6 j& \  m" ~%plot(x11,y11)
8 L( ~9 U# C3 L. Q
6 W- |: v7 s+ `8 M' m; p& \
%第二曲线方程 EF  EF  EF
t21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
/ z0 B; ^+ K" }4 D' q5 Xp004=R2;
$ Q4 f/ g' J9 y7 m* k* X8 h%PP=linspace(p003,p004,100);
( ~/ ?; _4 y9 r6 Jqm03=1/i*(acos(k*p003/A)-t21);
2 m; @$ R& f5 F7 V0 r- e! dqm04=1/i*(acos(k*p004/A)-t21);
( x0 g' F0 q4 o$ Y2 U. \qm1=linspace(qm03,qm04,100);
x12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
: ]- L: e8 D! }0 gz12=0*qm1;
# j) @  I  w3 ~5 a. {* `- ^plot(x12,y12)
9 g& ~  B2 U) v( G2 g' n' w
/ T% f9 u5 }4 ?  o- D
4 X) R- c/ f8 A$ HEF=[x12',y12',z12']
* S! F3 o2 z" T1 h%save('EF.txt')
0 r+ p+ {# O, \- z  h" H' _end
' Y; Y' E+ v; n6 u+ g! [

  o1 Q6 l2 {" [5 `+ r%CalpaMFG.m
%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
% |( j+ K4 q0 r% i8 l! ?7 U. xfunction [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
- B, a7 E0 h. f3 r4 W. MR1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
% V) ~3 ?* G3 `9 _t=linspace(0,t,200);
0 F4 W1 s1 J) ~+ L/ }  q# Ux1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
z1=0*t;
plot(x1,y1)
%
FG=[x1',y1',z1']
%save('FG.txt')
) N9 x6 r/ a8 |, D( {: Mend

+ e) g- {/ K9 f4 v! {3 t  C
: e( i) m, M" F  ]% [ x1, y1] = CalpMFG(110, 5, 6,30)
4 s0 v+ a0 r! {% f0 h& e

%CalpaMGH.m
%原始不对称型线计算程序
- N5 J: s7 i" E& }! }6 pfunction [ x11, y11] = CalpMGH(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
' Z3 F# u" ]! I5 v* h0 }8 `2 RR1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
2 l" N- X9 ]1 [7 w0 P1 `1 I%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
%t=linspace(0,t,200);
. ]* W- ^/ G' {+ f' U, Y. H%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)
; y* Z5 O2 |* X
" [# }/ L: L5 X7 O' C* n) J+ w
%第二曲线方程 GH  GH GH
b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
t1=0;
% ]3 R9 {8 P/ P* C3 C%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
0 B3 O* i- @8 U( |0 X" K6 u1 m) l%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;
4 d) F) q9 c0 k6 c%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
%plot(x22,y22)
- D' ]" P* g7 k7 c3 i

2 s& y% @1 ^/ {" C
%第三段曲线
6 n% a5 [& P& f% F7 `/ ~%t21=acos((2*R2^2-R^2)/(2*R2^2));
& Z- s" y, _. {# h$ Z2 ~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);%计算线段长度
$ y9 x5 b" Y- w1 k9 \t22=acos((2*R2^2-cp^2)/(2*R2^2));
P001=(A^2+R2^2-2*A*R2)^(1/2);
! _' N+ q1 x0 N) e" y. Y& x2 jP002=b1;
qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
3 e" P* I. a/ I( O/ xqm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
8 R  W$ [) b' @, N  Iqm=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;
) m; e  U/ q9 Y1 tplot(x11,y11)
  h) l, g1 ^! \( {+ R" k/ c%
) u! t* |! O' N# y2 [; [GH=[x11',y11',z11']
+ N2 U. k. X8 v- g%save('GH.txt')
8 V; c6 b4 n) j" S: g+ Q0 Qend


: j8 B( `: D9 e/ t