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

shouce

%CalpaMEF.m
8 d2 c) d* a, u- r7 L7 X% P! h%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
: n( o3 h& v" b- D+ M( r8 Xfunction [ x12, y12] = CalpMEF(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;%销齿圆弧的参数范围 在等腰三角形中求
1 I* n; C4 g7 Q: \) g% j! D%t=linspace(0,t,200);
( Z# v3 N  C3 {* V5 r) a7 N3 e%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)


: r1 n: D6 X/ \%第二曲线方程 GH  GH GH
5 b* F- [) f' }%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));%转角参数
%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
%q=linspace(q1,q2,100);
  _# F  `0 f7 b) Dk=i+1;
9 @9 h7 n" S/ V1 u. N: ?* v) \% L%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
) B$ b6 l% Z* K2 ?$ e4 ~%plot(x22,y22)


%t21=acos((2*R2^2-R^2)/(2*R2^2));
%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
- e% ~# t. `9 j4 l# X%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));
, C! v8 a. L, `7 b. I%P001=(A^2+R2^2-2*A*R2)^(1/2);
  n: j3 f+ I9 z; f2 }3 `%P002=b1;
%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
& f5 [% c/ D- r' a. v4 H! Z+ m%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
# o/ y. a/ v( u9 M; V* ]$ @%qm=linspace(qm01,qm02,100);
  F" z- A* \2 T. S%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
%plot(x11,y11)

1 a4 H, T9 F, k) [, `# X- x
5 i- E' C! W# e4 x) K$ B%第二曲线方程 EF  EF  EF
! F3 \- m; \9 D0 p' z* \9 F* a& tt21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
p004=R2;
" o6 T, `, w2 u7 X( i( J1 P% w1 `%PP=linspace(p003,p004,100);
! u9 p/ J7 H* y& U' eqm03=1/i*(acos(k*p003/A)-t21);
: e; y! r" j4 S4 Cqm04=1/i*(acos(k*p004/A)-t21);
qm1=linspace(qm03,qm04,100);
/ w7 y& R1 v0 X5 `0 G: sx12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
z12=0*qm1;
plot(x12,y12)

, y; j9 O1 f8 E% X
EF=[x12',y12',z12']
%save('EF.txt')
( L5 y# m$ P) R) I) B5 z' zend

, e8 o  B1 c- V6 @( I) E! O! c! d
%CalpaMFG.m
. L6 l' Z6 }9 d' X# n4 P- m%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
/ L/ U& E4 G: c3 i. U) yfunction [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
6 Z8 F  G# |4 K2 K$ J& S' YR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
0 x; w' y& l. Z7 I; i/ Gt=linspace(0,t,200);
x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
+ R7 J5 q: Z+ q# k5 l8 `! hz1=0*t;
plot(x1,y1)
. U! l9 w+ n2 \) Y" a0 C%
FG=[x1',y1',z1']
%save('FG.txt')
end

. i' x, y1 ^* k
% [ x1, y1] = CalpMFG(110, 5, 6,30)

9 z. M  m0 M. L2 \: ]. [( G
%CalpaMGH.m
& b3 L, S4 |2 {3 X7 w: Y, o2 b%原始不对称型线计算程序
function [ x11, y11] = CalpMGH(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;%销齿圆弧的参数范围 在等腰三角形中求
9 x- O- S/ ^) y/ q%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)
3 I/ [" T0 ?8 L: Y% m
! L% ^. n' }1 o' g  S! {/ @
%第二曲线方程 GH  GH GH
b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
2 Y' w% q1 _) q7 ^2 g& L%t1=linspace(0,t1,100);
q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
( R/ L1 ?9 A: T0 P4 v% w; j- ~  z%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
%q=linspace(q1,q2,100);
% ~! e* w, U$ T8 [$ B' V$ Ok=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)


; y0 F8 M5 S& u: {. }
%第三段曲线
%t21=acos((2*R2^2-R^2)/(2*R2^2));
x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
. G0 j5 j4 O4 ^0 A  ~y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
. X  d' Z+ a) T8 a. N6 ht22=acos((2*R2^2-cp^2)/(2*R2^2));
& p6 o: t2 w( [/ e" TP001=(A^2+R2^2-2*A*R2)^(1/2);
1 |5 z- p" i# j2 E1 xP002=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)));
# B) o+ P3 G5 hqm=linspace(qm01,qm02,100);
x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
3 P/ g) J% G& r) a5 Xz11=0*qm;
) ]5 R! m* ^  j3 N& Uplot(x11,y11)
( S/ G8 G1 m8 L  J7 R' d' X% }9 Q%
GH=[x11',y11',z11']
%save('GH.txt')
9 q, i0 h$ y8 P7 o7 y- `0 Lend




% K' w! y3 c; x' A' }/ W& i1 Q