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

shouce

%CalpaMEF.m
. x* k; ?$ x$ N) N%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
5 F; B! i# J9 O5 YR1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
" b- G' h& B/ ?6 b%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曲线段
%plot(x1,y1)
5 s! Y  l" T7 Z8 M2 w2 u

%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
# }! L/ z5 B. B%t1=0;
7 ^1 G3 U" [3 n% q3 m$ ?+ b2 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));%转角参数
1 {$ N) m- b6 E3 q6 C%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);%曲线方程
! n! f) r! l7 x5 O4 `  \9 B%plot(x22,y22)
9 i5 W, N" |. Q9 H
$ u2 @. y* a- L2 |/ Y0 k& K8 F
8 l) L" u# Y: e/ O* s6 v%t21=acos((2*R2^2-R^2)/(2*R2^2));
  ^( e  K8 t* U6 n3 B: K0 W) E%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
4 U+ y/ t, ^& [) O$ U5 x' Y4 n%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
# p5 {  q6 J) B8 a6 k7 B%cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
%t22=acos((2*R2^2-cp^2)/(2*R2^2));
+ b. F) r2 d7 v$ V5 w" \( L%P001=(A^2+R2^2-2*A*R2)^(1/2);
  @, X' J% N& J! C, m. n8 E%P002=b1;
& @6 c! Z( R! c% x7 \! e%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
; ~4 B! |9 Z1 N; N' o%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
; v6 f8 K  F  X4 v0 L# U%qm=linspace(qm01,qm02,100);
%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
# q7 S1 M, R% J0 h%plot(x11,y11)


* i6 O$ X$ m& E0 Y- t' V%第二曲线方程 EF  EF  EF
: h  C8 Y. ]2 w# Gt21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
/ B  H2 m9 {; E! k/ p( b. Ip004=R2;
%PP=linspace(p003,p004,100);
qm03=1/i*(acos(k*p003/A)-t21);
3 l% d  |, M/ _6 o. n$ |/ {qm04=1/i*(acos(k*p004/A)-t21);
  I- D3 p3 l6 e2 `; S. v& Tqm1=linspace(qm03,qm04,100);
! I7 u9 X, p! ^' Cx12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
z12=0*qm1;
/ z* n$ g: N2 Rplot(x12,y12)
( o, C0 g. ?8 M0 M
+ _0 k& x! v6 Q) b- C
) Y9 C  v5 U4 j6 E0 zEF=[x12',y12',z12']
%save('EF.txt')
end
0 y4 \1 U6 U# M. u" b" _
) S. O; I+ \# }( W6 N* p
%CalpaMFG.m
%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
: ]5 A' C4 U. t/ V" U( p; Wfunction [ 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;%销齿圆弧的参数范围 在等腰三角形中求
: i* |- n" |% r# r9 _; w# c4 Ut=linspace(0,t,200);
x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
z1=0*t;
3 I( ^0 Q4 M+ C9 d$ ~plot(x1,y1)
' N$ n$ T. a3 u0 q7 U%
FG=[x1',y1',z1']
$ l7 g# |8 f* _4 d$ ?8 t! T5 {%save('FG.txt')
* I: P- U" ~/ Q2 Y8 u6 R: v  Hend
& a4 }( M, v3 F1 n4 {+ C

6 k0 M# C  n. J3 `; _) ?% [ x1, y1] = CalpMFG(110, 5, 6,30)

3 C+ `' S# t4 @! n6 H9 p" ]" e
/ X! d8 i) B# K( D# M0 ]%CalpaMGH.m
4 q* E$ Y3 o5 @/ S& W; S%原始不对称型线计算程序
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
% O, F/ O* v/ p9 a) z. ei=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
) W; Y! ]( g, qR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
2 [+ a4 r* Y/ V6 V%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
# u8 i" e8 Q  x# z" [. c# q%t=linspace(0,t,200);
  N' u' d( a7 X9 N: `( R%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)

( }7 [: _# C' Y3 l0 [/ o$ e
%第二曲线方程 GH  GH GH
( I- X$ a. c7 Cb1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
& t2 F2 ^! B/ C1 q7 a* ^%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);
. _1 g: Z! O/ p  ik=i+1;
%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
" m, ?- G  H/ U9 T9 n3 O%plot(x22,y22)
; S* Z( L% E' L# T; m6 t


%第三段曲线
8 V$ V7 e: e5 m5 O! e%t21=acos((2*R2^2-R^2)/(2*R2^2));
2 c( b$ L9 C6 x% m- s+ ux0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
% X) j5 U7 B5 c0 K0 Jy0=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));
$ S! a6 d3 t* W0 a  l1 I' |P001=(A^2+R2^2-2*A*R2)^(1/2);
P002=b1;
( L  H: H4 O9 _! Qqm01=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);
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)
3 |* B2 w0 _& f# x%
" S3 P' D. P$ L- T. sGH=[x11',y11',z11']
%save('GH.txt')
end


9 k9 w) F0 z& _  Q8 N

% t- k+ `& d* Z! M4 e5 ~7 a7 G