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

shouce

%CalpaMEF.m
%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
8 z& b4 k: \  H1 {6 R/ mi=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
# z7 X5 W. q* [8 rR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
%t=linspace(0,t,200);
, j- T3 R. L) w7 |( X. t5 }%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)
; X( n6 O4 J. H8 Z! `
. V6 ^. N+ o$ v. a- p" y; E" \
" y* ]) a$ A- ?( B' E9 n# m%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
%t1=0;
2 N* w, g: M3 E$ k1 f" p%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));%转角参数
" q( S2 g) j8 K2 L2 e%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
: N9 c: @0 G- i( Q- C9 |%q=linspace(q1,q2,100);
. c3 s, Z. h6 R, u3 e! ~k=i+1;
%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
%plot(x22,y22)
# _( [9 Q" i. c4 A! Q4 [
4 t- z9 o8 H* y  y% y
+ v2 w8 G6 `4 V9 i; {9 f* |%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点纵坐标
%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;
%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
$ m5 g2 }" C/ V. y%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
1 N; F3 c! n. R3 m, S%qm=linspace(qm01,qm02,100);
! _$ b3 z( M7 H%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
% w: s8 R- _$ l+ x: M%plot(x11,y11)

# i; y! c9 W! j" h* P
%第二曲线方程 EF  EF  EF
t21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
p004=R2;
, [" j: M/ C: E6 i" P%PP=linspace(p003,p004,100);
qm03=1/i*(acos(k*p003/A)-t21);
1 L/ V: H% F. oqm04=1/i*(acos(k*p004/A)-t21);
qm1=linspace(qm03,qm04,100);
# O2 c; }2 d; r+ R8 z$ {# N" dx12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
z12=0*qm1;
0 A( p/ l2 R: h! v$ E1 Cplot(x12,y12)
$ p+ M( h( w/ q. q! k
; L  f( ]/ r( [9 ~: d, I# Y+ N3 m
+ t; r( _& u2 yEF=[x12',y12',z12']
( F  q+ ?8 [& @. o3 c! i  G%save('EF.txt')
end
9 g% _5 Z- h/ p1 E" }5 ~; M
" O; y. O3 @5 F3 r. \; \( v0 r6 m
& v$ \5 G9 @$ h9 X%CalpaMFG.m
7 ?! O4 [! j9 E6 j7 [+ @%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
function [ x1, y1] = CalpMFG(A, Z1, Z2, R)
* ]2 B9 ~0 D& K7 A+ di=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
' B+ [/ T8 W: H& \; B! eR2=(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曲线段
& h. a, r; e4 @' R7 {  Y, [z1=0*t;
7 m/ ?* y( z+ U) a, H' e; Dplot(x1,y1)
%
. m0 e: D) g. J, R5 F1 [FG=[x1',y1',z1']
%save('FG.txt')
+ f8 f# Z/ f5 I; q- {! o; Tend


% [ x1, y1] = CalpMFG(110, 5, 6,30)
9 O1 I$ f; n& P# l; y' {
; m) }. E# `# E; D) I. u* i
%CalpaMGH.m
%原始不对称型线计算程序
5 z, T- h) k! p6 r% zfunction [ x11, y11] = CalpMGH(A, Z1, Z2, R)
- R8 O. Y# K5 D( Mi=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
! ^/ U1 K1 D7 X( i  Z%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曲线段
% K' x) U  J1 E7 p%plot(x1,y1)
, @- A/ [4 x5 R: S, q5 o

%第二曲线方程 GH  GH GH
b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
t1=0;
, W( I/ P' }& L+ z0 v1 u' n4 k, m%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
5 m( O8 ~5 G3 y) J3 I%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));%转角参数
. G  E6 Z& a' X3 r%q=linspace(q1,q2,100);
k=i+1;
3 @9 W; _. j' v) q) X%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
%plot(x22,y22)
) \. ^1 I' S! ?6 F& m, D  r


%第三段曲线
9 j2 L/ G2 q" T' s& V" ~, s5 w%t21=acos((2*R2^2-R^2)/(2*R2^2));
, B" c' q( C4 d/ ~+ vx0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
- ~3 R) B5 C$ i6 T8 iy0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
+ L0 V& O" Y$ i% A, N+ Z5 f$ ?t22=acos((2*R2^2-cp^2)/(2*R2^2));
* l4 Z* j6 O, p# ]6 R7 ^, qP001=(A^2+R2^2-2*A*R2)^(1/2);
; T0 l. h( Q# U" b4 kP002=b1;
qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
9 O5 n+ l! q$ Fqm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
: K. N- E* T6 dqm=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)
! _$ D- v  D+ i) |! B6 g, f0 {/ f%
GH=[x11',y11',z11']
%save('GH.txt')
9 J& S# S& J0 V" a. h: Zend
" A7 _4 Y1 q* R2 V
; _1 H1 U5 }9 @

  ^4 ~7 N- }+ `% x' l/ I2 [