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

shouce

%CalpaMEF.m
4 Q" S. E: ~. l& }%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
& T. y8 H, e% a' {9 c2 o0 @( {R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
( r" C+ O: F# x0 P%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
: L+ O/ s( g# w/ p% e% X%plot(x1,y1)


%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
%t1=0;
% X& X& @6 T0 M2 z7 A%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
( I; k5 i; D  `; [9 c9 G/ w, y; B%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));%转角参数
+ w! w8 W# ]7 R7 K%q=linspace(q1,q2,100);
; c( m9 Q# X" N5 Tk=i+1;
%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
1 r  {5 |  |1 N- S! W( N%plot(x22,y22)
0 N. }+ H7 K* p) y" G
3 i% ]- U" L6 J' G3 x6 s
%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);%计算线段长度
%t22=acos((2*R2^2-cp^2)/(2*R2^2));
* ?( c, K; z7 k+ g5 [0 h" Z%P001=(A^2+R2^2-2*A*R2)^(1/2);
5 s1 A* f  o* Y9 @: ]( P7 k%P002=b1;
* |8 Q" z6 [; w( _$ g! [: e%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)));
8 U3 V3 u) j5 V3 C6 d& P%qm=linspace(qm01,qm02,100);
%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
! w) ~# Q: g* g! G%plot(x11,y11)


  q- t( L; o1 S( b# Z6 e# t+ g%第二曲线方程 EF  EF  EF
6 [# G* {! C8 s2 R( n$ ]t21=acos((2*R2^2-R^2)/(2*R2^2));
! @1 F, `( Z: r7 zp003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
p004=R2;
%PP=linspace(p003,p004,100);
qm03=1/i*(acos(k*p003/A)-t21);
8 N( j6 ?6 q3 Q  ^# q8 y' W0 kqm04=1/i*(acos(k*p004/A)-t21);
" H4 G+ V+ p0 y, R5 r8 Wqm1=linspace(qm03,qm04,100);
6 y# ^# a3 v- j3 n) h' Y' T) _x12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
: P9 F2 p) D+ Q4 k8 X. S/ h. t3 pz12=0*qm1;
0 ]$ a- u/ u5 v# y9 Z# _4 uplot(x12,y12)
) M& t- t5 ?8 ~  q( _' K
; p6 [  X- ~$ ]0 U1 l  n% t+ j: k
( \+ {  X& V+ ?2 U* G7 Y/ c* mEF=[x12',y12',z12']
%save('EF.txt')
end


8 J! ^0 H" d2 ~! ]  W%CalpaMFG.m
0 G& l; ^5 [% {; l%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
* _  t' b& z9 t9 H* Gfunction [ x1, y1] = CalpMFG(A, Z1, Z2, R)
6 R5 H4 L4 m+ X  _i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
R2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
5 l2 o( y5 Z# y4 E1 V6 N( Bt=(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曲线段
) B7 c& p0 F. k* o  @z1=0*t;
" j* Q+ e( t1 `% f) dplot(x1,y1)
9 X% N0 _! C5 a5 E%
FG=[x1',y1',z1']
# R$ g9 {; ]( f$ [" d%save('FG.txt')
end


% [ x1, y1] = CalpMFG(110, 5, 6,30)
5 L8 x; G4 k3 M' Z. u1 s7 a
( t& M9 |+ K  r
%CalpaMGH.m
%原始不对称型线计算程序
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
3 p  ^$ t8 _5 N2 |& r4 e9 n% di=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
" ]4 b1 X8 [7 b; m' VR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
" C. c" a$ t+ C0 \6 N) ?%t=linspace(0,t,200);
%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
- i. g. u. p" R) M1 c' W8 v& s%plot(x1,y1)
+ N+ Y) V7 S: \9 j4 Y
) D. s4 r' B' ]  f# n& O+ F0 A3 m
9 ^6 D8 m% H- J. H%第二曲线方程 GH  GH GH
b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
% Z+ g$ W1 ?) I( T: M8 G* ft1=0;
0 y0 P/ z  [3 I* ~%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));%转角参数
: g6 b1 |) [& l! v) o6 l%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);%曲线方程
%plot(x22,y22)


* i/ p- j! t0 c$ l
! `' Z( C7 p" k7 O2 L, ~# T6 c%第三段曲线
; S* m3 e* E$ S% D. S% ?%t21=acos((2*R2^2-R^2)/(2*R2^2));
3 Q* t% j( ]4 T5 Nx0=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);%计算线段长度
t22=acos((2*R2^2-cp^2)/(2*R2^2));
; i8 j% K3 ?) e+ G. dP001=(A^2+R2^2-2*A*R2)^(1/2);
7 Q0 T. s  S# w7 w  KP002=b1;
! I- ~: \# ^/ D3 Fqm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
+ m/ r/ U6 _: lqm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
" T" X+ z) E; H* x1 Y8 z" [, Q6 Cqm=linspace(qm01,qm02,100);
! t! _( d' c# P2 C- \1 wx11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
' _2 J7 [" g' vz11=0*qm;
% j1 W5 c$ B4 b, |7 kplot(x11,y11)
%
GH=[x11',y11',z11']
: P" e8 V5 W/ V' x! Y3 i%save('GH.txt')
* T1 [2 B/ {; a$ dend
3 x/ \/ ^/ e5 ^. O. f

% G2 |; w' ]: E/ T, S6 D  F