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

shouce

%CalpaMEF.m
$ m5 ^/ p: [. Z: J%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
: ^" Q  m9 r. Ifunction [ 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;%销齿圆弧的参数范围 在等腰三角形中求
7 P3 L1 B1 m9 y+ P%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)


%第二曲线方程 GH  GH GH
! B$ S" z# M# R% D, C%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
%t1=0;
( F6 }- O/ S# P% Q. q, c) A%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
! ~! B$ W- U8 D9 b4 Q2 s%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
  C0 w' Q9 w1 ^6 u8 q. J%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);%曲线方程
7 ]# m: _) k+ @%plot(x22,y22)

' P3 c4 Q3 b& g( g. C
%t21=acos((2*R2^2-R^2)/(2*R2^2));
%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
; a8 w5 U( ~* G. e* O2 D%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
. v+ k9 `2 H! O* N4 e%cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
%t22=acos((2*R2^2-cp^2)/(2*R2^2));
% M; ]2 `$ `$ i%P001=(A^2+R2^2-2*A*R2)^(1/2);
; u6 Y/ M2 Z9 C0 C% I# t%P002=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)));
7 b9 U/ y/ U- O; p0 C' 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); %方程
/ y6 v0 M# T1 \5 {2 ^  E%plot(x11,y11)


9 O+ Z- {( @' Q%第二曲线方程 EF  EF  EF
' ]& h, [" P! {# f7 D5 u9 L$ it21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
p004=R2;
%PP=linspace(p003,p004,100);
qm03=1/i*(acos(k*p003/A)-t21);
, ^2 b+ z) K( W( U5 T/ ^qm04=1/i*(acos(k*p004/A)-t21);
qm1=linspace(qm03,qm04,100);
0 M) t( ~! T1 m) P3 b4 Zx12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
% u# f4 I* d, W3 `1 bz12=0*qm1;
plot(x12,y12)

; B" M5 ]$ l! j! L
EF=[x12',y12',z12']
- f- n# d* b9 K6 n5 I8 C; j% y%save('EF.txt')
end
4 X) P5 S  D7 g. Y- l" w
  g$ l$ m4 g1 J1 v
8 z. G, q! X% U' p%CalpaMFG.m
%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
0 \! c9 B4 D. k) T" ofunction [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
1 ]% p; y6 t' jR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
! E  `: p, [8 t1 s1 j/ y0 ct=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
, c! T; h8 K% q' G3 G# l# n2 \/ o% ?' xt=linspace(0,t,200);
4 H- q8 y! _6 q  Ix1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
z1=0*t;
( j1 {3 G9 M& k1 O  Zplot(x1,y1)
%
% m0 j; j3 m) K: mFG=[x1',y1',z1']
% ~. S5 Y" p# p5 {3 Y# v0 L% \%save('FG.txt')
end


$ ?: e. l( |# R8 c% [ x1, y1] = CalpMFG(110, 5, 6,30)
' h% o5 T# l* ~8 \& c% J
# n' j1 \& |, {' X/ B' B
& p  L. G! @% q3 P, Z%CalpaMGH.m
%原始不对称型线计算程序
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
- P& H6 f- |, O3 a/ o$ Fi=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
+ A" V2 @6 L- L%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
, z. Y- {: @, I3 P5 q2 ^+ ~3 P%t=linspace(0,t,200);
$ C! K9 j# j, ?" V%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)
$ N0 l1 e. J' L0 b; p: J; M
% ?$ A$ ?( W9 ]8 U& B+ ^
%第二曲线方程 GH  GH GH
+ ]+ D. S8 w9 _+ L- }b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
t1=0;
# r, N& r1 n! ~# m8 s2 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));%转角参数
6 a2 s( q, J* B9 t. v+ I0 ^) `%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
%q=linspace(q1,q2,100);
. `- k$ G$ e  W1 Lk=i+1;
3 G+ v9 y2 ?  o. D( [%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
%plot(x22,y22)

" q  R4 C' n; f" ~  n0 q

%第三段曲线
%t21=acos((2*R2^2-R^2)/(2*R2^2));
- p/ u& N; @8 O& lx0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
5 u/ \' @. ^; f- N) a$ ?7 Wy0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
, @. j8 C5 Z2 jcp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
t22=acos((2*R2^2-cp^2)/(2*R2^2));
& \8 `8 N. I, F  Z5 sP001=(A^2+R2^2-2*A*R2)^(1/2);
P002=b1;
& t6 f' v$ p8 p- e7 h$ Yqm01=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)));
5 O1 {8 n3 V  M0 O2 E8 sqm=linspace(qm01,qm02,100);
x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
/ j: D6 V( C- A% {z11=0*qm;
plot(x11,y11)
%
! l4 T. Q$ \$ Q( y( |! UGH=[x11',y11',z11']
4 ^; {  H6 S, L1 H( J; J) U%save('GH.txt')
end

% C" K% G8 O4 R8 s
' V4 F) c2 w1 [7 }, ]
. `5 c, w. y3 W5 n1 T
+ x2 O5 P2 N8 w& Q- {' S. @7 x- f