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

shouce

%CalpaMEF.m
9 {! H% j# c) A+ s%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
* @( S, c1 e+ u) X: ]8 ^R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
' n1 Z0 B- n! ~  e%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
%t=linspace(0,t,200);
$ P5 z6 K8 Z6 K7 C+ b%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
1 s& x& |0 ?1 D%plot(x1,y1)
: C' m% S/ m% j8 X2 x

8 s& ^7 n! E6 Q7 c8 \3 j%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
: S1 Z  p+ `& ~6 b4 B%t1=0;
& T% }4 e5 Y3 w+ G+ x$ \$ r%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));%转角参数
%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);%曲线方程
' S1 G* B3 e8 o* g%plot(x22,y22)
% ~+ w" y! ?% N; v, \" y
2 G$ p1 ]' B, D9 g
%t21=acos((2*R2^2-R^2)/(2*R2^2));
%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
0 z2 e. }. E: L/ H%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
" e1 u0 q+ b' Z7 g. K%cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
%t22=acos((2*R2^2-cp^2)/(2*R2^2));
6 j- k2 W, B7 l%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))); %第一次在这儿括号输错
2 z2 {+ j! j2 ^! n0 D) B%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); %方程
. K: U. I' @0 H3 F: [%plot(x11,y11)

, m4 w+ f! R) i' E6 P+ x9 Z
$ V- p! k+ f  ?$ Z& t: m! `! @%第二曲线方程 EF  EF  EF
+ R  c% w3 l2 a' Nt21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
5 \7 h4 G4 n" zp004=R2;
* H' T/ {* }" Z: E* @. y0 X' L4 f%PP=linspace(p003,p004,100);
+ w$ K8 Z' ~! S2 U% iqm03=1/i*(acos(k*p003/A)-t21);
qm04=1/i*(acos(k*p004/A)-t21);
qm1=linspace(qm03,qm04,100);
. e3 A5 Y2 a9 L% A8 O. ?" jx12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
z12=0*qm1;
plot(x12,y12)
; Z0 W' h/ u8 ^& M4 A- x

3 f# W3 k1 w5 M% \4 Y' e' mEF=[x12',y12',z12']
% p! a3 J6 o6 }5 N: `: ^; x! ]% a/ p* y8 o%save('EF.txt')
end
3 k' s  a) _( p! Z. D0 N/ j
+ }3 ?1 A; l* M8 f( D- ^
%CalpaMFG.m
%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
2 L, g/ x  }  q1 e! Afunction [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
* z- K3 T4 B4 Q2 FR1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
- d* f- u- @5 N7 e' ^. k6 r1 B: zR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
. E- V% p( E9 vt=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
5 _8 j1 J+ m0 Q) e( Dt=linspace(0,t,200);
/ a0 t7 D( m& }2 e. Fx1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
7 ?5 V- k/ P. R  [, e# k+ ~6 Sz1=0*t;
! a/ @/ h* J1 Y$ Kplot(x1,y1)
%
4 a) R& ]8 |- T0 ?FG=[x1',y1',z1']
* ~) _$ H7 w# ?4 v5 J* V%save('FG.txt')
5 i3 ?! B' Z* s1 e0 G! iend
( s$ A/ y# y: O! @2 D5 N8 A
% r: D( K* h- Y% `
5 `+ S! `* T- R$ F3 J. ?+ R% [ x1, y1] = CalpMFG(110, 5, 6,30)

5 S# U3 _# }: Y, H4 `% |
: ^1 o# }" [, z%CalpaMGH.m
%原始不对称型线计算程序
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
( I. V- L+ |. ]6 ~i=Z1/Z2; %齿数比
0 |- h$ ~5 i5 N. Y/ a6 qR1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
0 y/ m7 _' _, ^, t, v. ?% WR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
. ]. d4 O6 Q6 C%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
%t=linspace(0,t,200);
3 y2 d: F: @) q7 C  _" M+ q) l; c%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
  i7 V# A$ E* ^" Z! P5 V%plot(x1,y1)


%第二曲线方程 GH  GH GH
$ Z8 E% P! T' N5 r6 u7 Vb1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
. @/ g( s3 ?: U' U4 k: x6 R% Z/ st1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
. i- q# t. U+ X1 F- V/ C. b%t1=linspace(0,t1,100);
) ^2 z# o2 h5 V# Aq1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
2 Z* p) m- X# K4 G: Z* {" Z%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;
- h9 _$ p9 F, q3 z5 l%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
" \3 B6 g* l1 Z. ]- i& {( ~%plot(x22,y22)

3 i( \( W1 D+ j: @( i/ _
* O, V+ _4 a  Q4 D
9 j% J. |8 n! F- o" K7 @%第三段曲线
%t21=acos((2*R2^2-R^2)/(2*R2^2));
  k3 P' ^  G& D* j$ ]& X4 Fx0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
5 \1 F( e, A6 Q4 ?' e7 `; jy0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
2 h: w; O5 ]) h! Ocp=((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);
! h' k6 z3 ]& Z# eP002=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)));
  E% v! P- u$ o8 Y6 N" j" bqm=linspace(qm01,qm02,100);
% H3 ?6 a/ A& E6 z$ Nx11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
z11=0*qm;
plot(x11,y11)
6 k2 x% k  z* D%
! `" q% i3 ^/ }- o5 h5 }GH=[x11',y11',z11']
%save('GH.txt')
; G# {% B6 B, `* {* ]. oend

& _1 i# [7 \5 b; ?: T& e3 Y$ p