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

shouce

%CalpaMEF.m
& S: H% z! w1 c1 W# M4 |$ N%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
6 F* V' z% ^7 A7 a/ v%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
, f  t; |$ h, [) h, H1 j%t=linspace(0,t,200);
, S& c& P. I6 P2 D%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)
' n1 b0 t$ X0 V8 a
) C: w: b/ W4 e! ]6 s! }' R# x0 H! j5 V
%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
5 ]$ q- T* w. @) M7 c%t1=0;
& I; J: S0 {% i* S& _7 a" @%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
$ _8 e5 O: K) h% O%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);
( l- _* m; S( i8 z, F% Tk=i+1;
6 b9 n1 Q/ T  T%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
%plot(x22,y22)


%t21=acos((2*R2^2-R^2)/(2*R2^2));
%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
; T8 O: t6 ?' \* C+ q/ N8 f%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
# @1 W% u, ~7 ^& h+ j%cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
$ \, e  v% d8 @  }' M%t22=acos((2*R2^2-cp^2)/(2*R2^2));
/ V  Q8 u( q- y9 e. a: w  W- T& }%P001=(A^2+R2^2-2*A*R2)^(1/2);
%P002=b1;
2 L" L# ?- ~. t: A( T%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);
%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
%plot(x11,y11)

4 g, [: Y1 p! ^; z. A
1 _# N7 i6 ?$ E9 R( F( @. h%第二曲线方程 EF  EF  EF
t21=acos((2*R2^2-R^2)/(2*R2^2));
, }1 l- O/ {0 u" kp003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
* r: }8 S9 {. i+ }% e$ bp004=R2;
$ U, R) ^, t  H3 ~%PP=linspace(p003,p004,100);
, b) m2 M! \" S. Tqm03=1/i*(acos(k*p003/A)-t21);
, ], o& w4 O6 W' aqm04=1/i*(acos(k*p004/A)-t21);
qm1=linspace(qm03,qm04,100);
" k2 ?; b: H  {$ r0 yx12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
6 l1 V7 G' |5 n, x) d+ Fz12=0*qm1;
plot(x12,y12)
1 ]" r& a: b7 g. M7 u
! ?( {2 z0 ~) y
EF=[x12',y12',z12']
3 M( e- S5 S: e2 I/ o%save('EF.txt')
6 @; ~7 k8 {5 p) p) c- C  Vend
- B) l. G6 R- u( ]
' ~; B' l# {/ [. r7 ]; a+ O
%CalpaMFG.m
%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
, s  R3 H$ }! ?1 ~% Q; x1 E( ifunction [ x1, y1] = CalpMFG(A, Z1, Z2, R)
5 l9 k# ]& ^' Vi=Z1/Z2; %齿数比
" Q8 L3 f8 E" L) }R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
t=linspace(0,t,200);
1 I/ k, W7 }; ^9 z. J6 bx1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
z1=0*t;
* F! ~5 ]$ d: J' Uplot(x1,y1)
%
FG=[x1',y1',z1']
$ c4 Z* d* I4 j( n( D%save('FG.txt')
# @1 Z9 O% J+ r, y' Dend
( P* v0 ^+ S9 l1 R* [4 }# k' u
6 x5 I$ ]1 Q7 j( B$ ]$ n! X
% [ x1, y1] = CalpMFG(110, 5, 6,30)

6 `; B5 r3 K: C! s: U1 D
1 i$ @0 N) f# t8 s  X%CalpaMGH.m
%原始不对称型线计算程序
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
! i9 `, A# d/ a2 }) nR2=(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曲线段
4 T; n5 _# Z1 O, Z, W%plot(x1,y1)


%第二曲线方程 GH  GH GH
' S- S7 {; X/ o3 jb1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
% p2 G% W- V- D3 |, O( h2 `* xt1=0;
, R! A& W: p5 m; w: c" j+ U%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
4 p9 `: Q. {; k% D; F" c) _# Gq1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
/ a" _( s: {: R8 u0 }5 Z%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
%q=linspace(q1,q2,100);
5 M% {4 o" Z+ r0 q+ Hk=i+1;
5 `2 K6 Y: C0 L& ]6 G* m%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
%plot(x22,y22)

5 w$ X/ |4 V5 ~/ c8 ?, T: a
# z* s- q: Y) W+ }9 x, b: ~3 j# d! P: c
9 G" R9 C9 `  ]6 Z. G%第三段曲线
8 }$ B% R% y4 x8 N4 _%t21=acos((2*R2^2-R^2)/(2*R2^2));
$ q, J4 F: C, p5 h- C6 o! Ox0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
4 f8 `( V7 o  H1 \) |9 A- q% S0 [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);
, h- f+ G/ K* M0 e# J, x: o7 N: ~P002=b1;
1 D9 H2 j$ q+ l# s1 nqm01=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);
. E/ w2 i) w, F" t3 {+ m- ^+ Bx11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
z11=0*qm;
plot(x11,y11)
" p8 K; T& U2 ?4 j+ f3 @4 T%
GH=[x11',y11',z11']
%save('GH.txt')
end


5 {" T; o5 S5 q6 n7 I