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

shouce

%CalpaMEF.m
* `: ~) C$ B# a/ e%原始不对称型线计算程序 [ x12, y12] = CalpMEF(100, 4, 6, 25)
\/ G( T% W8 S' C: L# b+ c* Efunction [ x12, y12] = CalpMEF(A, Z1, Z2, R)
- z$ p" K( J& ?4 {# ~" j9 di=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
$ r* Y( |. _6 O6 vR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
%t=linspace(0,t,200);
" a0 S$ b) |8 V%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程 GF曲线段
2 n. x9 J# E7 @# Q0 I2 H%plot(x1,y1)


5 b" D8 W) T; C& h8 k6 E%第二曲线方程 GH GH GH
7 T7 _# y* ~* y2 r; |%b1=(R^2+R1^2+2*R*R1)^(1/2); %这个地方第一次弄错了
%t1=0;
+ j- D: G* s3 G2 Q( a%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
) E& G' N5 |2 g% ]8 d%t1=linspace(0,t1,100);
* W: O k3 L1 I# A%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);
( C+ X( X6 x2 z) W# Fk=i+1;
%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲线方程
/ q4 Z5 `% @0 {%plot(x22,y22)
2 c1 g3 \% H4 T. s. h8 |2 M

E" [3 c% \- C4 K* L/ x. l%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);%计算线段长度
7 h& A" Y# q& ]( Q! l%t22=acos((2*R2^2-cp^2)/(2*R2^2));
; q' n. }" s: Q% O: L; Q9 [%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))); %第一次在这儿括号输错
4 L5 `2 S( O; j, a$ i. g%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
%qm=linspace(qm01,qm02,100);
2 d3 E [6 ?9 @% a8 h%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
%plot(x11,y11)
1 c! E x( M! \( x
4 m0 m& K: A7 ~/ }( E# z
; t% w5 g2 p' r: [4 j- Q4 p, s%第二曲线方程 EF EF EF
t21=acos((2*R2^2-R^2)/(2*R2^2));
' a' w0 ]$ R! X; A7 b1 ?p003=R2*cos(t21); %有点问题% 为什么是这个样的？
+ h7 W+ c( N( qp004=R2;
%PP=linspace(p003,p004,100);
/ @8 X. [$ `/ e& j4 Dqm03=1/i*(acos(k*p003/A)-t21);
qm04=1/i*(acos(k*p004/A)-t21);
qm1=linspace(qm03,qm04,100);
+ G. b3 @! J! ? r$ I2 @' sx12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
. ~ x- n' X0 l# W. L" oz12=0*qm1;
plot(x12,y12)

6 b7 s$ G( M p6 c, @' W
+ Y: q. u/ D+ {EF=[x12',y12',z12']
) e9 o) p% B5 a& G; v- K%save('EF.txt')
end
9 Q6 H/ O$ M7 d4 Y+ F/ E
! c' b% B# @# V' k6 J o: y# }$ P
. @6 y% Q/ a4 b& T, } n%CalpaMFG.m
: x/ Q" p A" l3 F%原始不对称型线计算程序 [ x1, y1] = CalpMFG(100, 4, 6, 25)
function [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
6 u: \7 \3 W+ f$ f, KR1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
" n. c0 `0 \- t$ H9 JR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
$ k4 V7 _ y' bt=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
. w' O( }2 |! j" h( ^7 ?7 w4 a Ut=linspace(0,t,200);
+ v9 s$ R. i5 J3 vx1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程 GF曲线段
3 e, o) I" G7 F+ {( N9 {) @3 pz1=0*t;
# i: z7 _ U/ z) tplot(x1,y1)
9 x; }$ H7 V/ ]: A, j* C2 i0 R, M3 w" f%
0 [4 n, A( c& R8 `6 }, q# t0 L' cFG=[x1',y1',z1']
& d! g: }' J2 C! n7 `) C, W%save('FG.txt')
end

% l; H1 A4 Y- c: M0 A/ L/ @% O
) M" ~0 S: T0 w% [ x1, y1] = CalpMFG(110, 5, 6,30)


%CalpaMGH.m
& V7 T9 E9 o! m7 i3 `%原始不对称型线计算程序
- n5 B* d. z/ R7 O) \1 \5 P5 O, mfunction [ x11, y11] = CalpMGH(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
3 G" I- P F* l# {. e$ }/ f" cR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
%t=linspace(0,t,200);
5 t( t9 B& y* B%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程 GF曲线段
%plot(x1,y1)


%第二曲线方程 GH GH GH
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);
% N+ J [* O0 u- C( wq1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
( I/ C4 o$ ?5 P2 f/ ]3 O%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
j' ?3 C! ~% P7 ?%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);%曲线方程
& w* ^% `7 R7 G; @" X%plot(x22,y22)
3 c1 z) x* a; [9 L& ~4 j3 A


# W' Z: b: _5 n3 K! Z$ r: E- T/ w%第三段曲线
2 p( W& ^; Z ]%t21=acos((2*R2^2-R^2)/(2*R2^2));
x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标
- B- T/ {: e8 ry0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
% E5 f% Q r' w8 z6 ]! xcp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
$ R8 @: m1 w, Y5 G- {5 a; P* g( @t22=acos((2*R2^2-cp^2)/(2*R2^2));
. L/ K9 m3 x4 B! a. hP001=(A^2+R2^2-2*A*R2)^(1/2);
% z! x# n8 O( R) w0 N2 Y/ }- FP002=b1;
qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
, v! y6 o- n8 f( wqm02=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); %方程
# K8 t+ f6 [( p; E+ m) h! H. S+ hz11=0*qm;
% M7 Y. n6 V, }( }plot(x11,y11)
! Y& y. N$ d3 `( x+ e* K( l( }5 K%
GH=[x11',y11',z11']
: J/ M* p6 e: O$ X4 R9 L$ a" ~%save('GH.txt')
end
) c4 u% z( g$ W8 B; ^



1 b) N5 t# \: P% G6 W