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

shouce

%CalpaMEF.m
%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
% P1 M% J2 H* b' Q+ o3 tfunction [ x12, y12] = CalpMEF(A, Z1, Z2, R)
, \. W+ U- h$ u% s3 B6 V+ X  ii=Z1/Z2; %齿数比
& `8 W  A3 Z% B; L4 T. TR1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
* d' V6 |4 ]; u* f0 e  g7 |R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
2 f+ @0 \( i# K* O/ t%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)

7 Q  n  b8 G. R- {
%第二曲线方程 GH  GH GH
) V' }4 ^8 w& M+ O%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);
%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
% E% b! b( C" Y8 t9 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);%曲线方程
$ F* E" l2 O3 a+ |/ R* T0 ^%plot(x22,y22)
; ]- f- S& S1 Z6 Z! H5 _

: ~; o, D" E' \%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点纵坐标
. f  l) m+ [& C* }1 I' R8 K%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);
%P002=b1;
0 U$ E- t( V: }* S% n%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);
  s3 |( i4 h, ^" ^3 g/ V4 o  U2 o%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
- o8 ^" i9 g/ b* |' k%plot(x11,y11)


9 K, f1 M. n7 y; M; x3 A%第二曲线方程 EF  EF  EF
5 j0 D: g! M1 v  ^2 u* \$ V$ C! e! q  st21=acos((2*R2^2-R^2)/(2*R2^2));
  Y0 v: b" ]3 f, m$ Y) pp003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
5 ?  T6 D6 b7 I6 f( kp004=R2;
%PP=linspace(p003,p004,100);
qm03=1/i*(acos(k*p003/A)-t21);
qm04=1/i*(acos(k*p004/A)-t21);
qm1=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;
plot(x12,y12)


/ d( G$ k5 e1 z. q, A; AEF=[x12',y12',z12']
# T$ Z! u: E& _& [%save('EF.txt')
; D' `9 q  w$ m) A/ W" d3 Wend

* Y8 D1 U' N8 Z3 F0 }5 @( k8 c2 l
%CalpaMFG.m
%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
) h7 _1 \0 w( `0 H3 B! @function [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
' k7 u1 P. W: {+ ~$ SR1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
8 y! s6 w6 I# ?: h. L2 Bt=linspace(0,t,200);
) a: E" K* Q% g: y! w. vx1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
z1=0*t;
plot(x1,y1)
%
/ j" W" [% Z, G& h3 JFG=[x1',y1',z1']
%save('FG.txt')
end
# Z. ]8 \" L/ |( ^$ ?" S/ u

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

' j, e( U& H6 f% i9 }) o
/ f4 R% l* G# _' K/ Q& i%CalpaMGH.m
& p: u' S3 K. l* ]%原始不对称型线计算程序
9 j# n7 p* r& s  o- vfunction [ x11, y11] = CalpMGH(A, Z1, Z2, R)
! N+ X" S# J/ n. [i=Z1/Z2; %齿数比
! n( x8 s5 f) [3 oR1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
+ w: ~. U* n# a& m8 q. N" vR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
4 l) |  v7 s$ J% R& u2 j8 E0 {%t=linspace(0,t,200);
' ?2 Y9 G, B' R1 {2 ~%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)
3 z& e3 \* Z; X! k

( r, {4 r0 @* C/ m%第二曲线方程 GH  GH GH
# u' }9 |. Y9 l6 i& P3 x3 Z1 Ab1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
' o/ ?2 @- f5 \: U%t1=linspace(0,t1,100);
q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
7 Y5 Y9 r3 P; q8 J1 p3 R+ t# a%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
%q=linspace(q1,q2,100);
; _8 n( {& u' W# J  A6 R5 Nk=i+1;
%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
, G; R) u. e7 I5 X%plot(x22,y22)
1 l( m  {  Q& Y+ @' b, \1 X
2 Z6 X% X6 l$ c, I8 g/ [
% r  C& N2 c8 V" a$ ^; M
%第三段曲线
0 C5 ~( o. @' G' f- D9 G  _; z%t21=acos((2*R2^2-R^2)/(2*R2^2));
) u$ ]- v4 D: Q% p" e( q2 Wx0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
6 o5 a+ k2 i/ i& Ey0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
3 m9 \6 `% e$ k8 dcp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
4 V" C. j  A# v2 l7 {4 H/ H! rt22=acos((2*R2^2-cp^2)/(2*R2^2));
+ I$ E5 ?! ]. K1 B: RP001=(A^2+R2^2-2*A*R2)^(1/2);
P002=b1;
  r/ N7 L# ^) U) zqm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
/ P# B' [8 d* U" }1 g7 |6 Fqm02=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); %方程
z11=0*qm;
plot(x11,y11)
%
* u0 Q( u# m  aGH=[x11',y11',z11']
# U1 U. V6 w* ?8 |+ l& j* S%save('GH.txt')
' ?$ m3 J7 `8 t7 S4 iend
9 u6 ]9 m  z9 q


! ^7 s- Q' w9 v/ G; m
9 H) h  V0 W6 ]7 F- r, U7 d