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

shouce

%CalpaMEF.m
) D. r" a2 H5 X* Y0 A- q%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
9 j7 h, f# s3 ^5 c" y+ Q* ufunction [ x12, y12] = CalpMEF(A, Z1, Z2, R)
3 g/ F( y: \% ~2 F- O9 B( ci=Z1/Z2; %齿数比
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);
%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)
* i) ~3 P8 M) [: t3 b3 Z
4 V4 b& ^' G3 k1 H% B+ y: x9 r
% }$ ~: d8 X' ]" E5 \" M. k" ?( J%第二曲线方程 GH  GH GH
7 B+ o/ ]; y$ `8 ~4 w: |%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
%t1=0;
! U# R1 h. D4 Q' U$ q% n3 t9 T%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
+ ?3 ]! l: Y% T+ q0 y2 X%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
* Z+ p1 S8 ~. c' J2 O% V: _, A%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);%曲线方程
( C9 r6 W9 j9 z8 \8 N, [0 L%plot(x22,y22)
  V' m' ^' u% w6 D0 g

5 L" l7 q! @9 E1 G' H7 X%t21=acos((2*R2^2-R^2)/(2*R2^2));
  \; s! k* K+ C% e( Z$ j. v%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
: B* R) o! y* O0 x5 {9 Z& y%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
8 k+ Z2 k# j& s%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);
' U* a9 n; z* G* s: a! g7 E- V! [, h' X%P002=b1;
) Z* p! x; j5 t  o* R%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)));
9 A! |/ J8 P& N& _& p$ v%qm=linspace(qm01,qm02,100);
%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
6 c5 \3 d+ m& X%plot(x11,y11)


& U) f6 m5 y' V( e/ F+ R6 h3 U%第二曲线方程 EF  EF  EF
+ ~- _5 I. z- F. ]+ q- u' P+ l$ |t21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
8 V. x+ t. Z) }- D0 cp004=R2;
, f5 Q- P' e+ u9 R2 y1 I; p' g/ G%PP=linspace(p003,p004,100);
0 |$ ]; L" s* w, n  o4 ^qm03=1/i*(acos(k*p003/A)-t21);
qm04=1/i*(acos(k*p004/A)-t21);
qm1=linspace(qm03,qm04,100);
x12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
z12=0*qm1;
2 p7 H& F" z+ `* wplot(x12,y12)

' i  I# m  F2 U% s4 E5 }
EF=[x12',y12',z12']
8 b9 U( D# Z. T' {%save('EF.txt')
end


6 d' o& f) t: n6 f  e4 [%CalpaMFG.m
%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
function [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
4 @/ [7 C* x% k' _5 u# _# GR2=(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曲线段
z1=0*t;
plot(x1,y1)
%
& O" \/ ?$ l3 \$ i( @+ p2 nFG=[x1',y1',z1']
+ X) p. G3 _6 j: ^  Y8 v# @* F%save('FG.txt')
end
+ Y% V$ \% H. X: a
. w; b; I2 N; {& Y3 ]
# [6 m1 A" ~- n4 ~% v! T% [ x1, y1] = CalpMFG(110, 5, 6,30)
4 Y& Z, L6 [2 X2 r: M8 o, I9 E7 i4 o

%CalpaMGH.m
- r& y( k0 O/ w4 \7 O1 M8 j' W: X: H%原始不对称型线计算程序
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
7 Y! ~" ]$ J% L# |) AR1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
0 e9 p) G$ i4 `- @/ H6 j6 r, MR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
! T/ D0 [6 Q2 O- @) a) H%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
%t=linspace(0,t,200);
* x. w+ [8 y: @%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
& ~1 [/ P/ }& H%plot(x1,y1)

% K) r: j6 x+ h* Y
%第二曲线方程 GH  GH GH
b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
3 a/ z; L* z$ u! Z7 Ct1=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);
k=i+1;
7 z8 \5 o$ E: z) T* L9 K" n%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
%plot(x22,y22)

: B3 C7 v/ t( }" N$ {* T. a

4 ^1 p8 s. m* B! ]%第三段曲线
! e- b- C4 T6 n. W  y%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);%计算线段长度
6 n% d: ?* ?& v. d3 o5 R0 bt22=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))); %第一次在这儿括号输错
qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
" ?/ |$ u' E; d. n; Kqm=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;
# \5 M, W0 j$ S7 S4 u: x) Pplot(x11,y11)
%
GH=[x11',y11',z11']
%save('GH.txt')
end




& }0 N& f+ U* _# E