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

shouce

%CalpaMEF.m
$ v. h3 n7 s* P" D" N%原始不对称型线计算程序  [ x12, y12] = CalpMEF(100, 4, 6, 25)
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
) k8 }0 @) Z+ }i=Z1/Z2; %齿数比
0 a' }( H' O$ m3 R2 D) [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);
; I- W2 Y4 z& O# T- `, a6 Z%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
4 B  ~3 W$ f' l0 c* h%plot(x1,y1)
- z3 ]2 }" ^8 `
3 _  w0 m. C2 O7 v) a3 N
%第二曲线方程 GH  GH GH
%b1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
. ]3 m! x- O3 O4 K%t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
%t1=linspace(0,t1,100);
6 J) n! o) J" N+ q1 N% Y9 L: {4 [%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
* a  U$ o5 u1 j+ Z/ U3 e9 {) T%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)
3 L& f# o+ V+ r( u
$ h6 D6 n3 F. ~! n1 p9 X" x/ ^
%t21=acos((2*R2^2-R^2)/(2*R2^2));
%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
! ~  q. W0 J4 v+ ?$ i3 P%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C点纵坐标
# y' z% J3 z5 l%cp=((x0-R2)^2+y0^2)^(1/2);%计算线段长度
%t22=acos((2*R2^2-cp^2)/(2*R2^2));
8 R! l, z: Q- a7 ^& m2 W* q0 r%P001=(A^2+R2^2-2*A*R2)^(1/2);
% B! Q2 L" [0 t$ y+ y. e4 C) Y2 T9 G+ N%P002=b1;
1 F- n) l8 T' S* ^- O6 n%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在这儿括号输错
$ Q  U" w( M- a%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
%qm=linspace(qm01,qm02,100);
! F4 u# ]+ c3 C" L2 E4 V6 x%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
0 K. y# c1 [4 ?) w%plot(x11,y11)
' p6 b7 L9 Y  l+ c) G) l
- F7 l# K( n3 h5 j& q. B% N9 z9 K; J
" U4 n; J& ~0 B+ ]2 [* ~%第二曲线方程 EF  EF  EF
, k/ T9 I3 O# Q% [. m8 ]t21=acos((2*R2^2-R^2)/(2*R2^2));
p003=R2*cos(t21);    %有点问题%  为什么是这个样的？   
2 Z1 D  Y% d' I6 Q, ap004=R2;
, k" w# r2 l5 Q2 ?1 q) E%PP=linspace(p003,p004,100);
8 w, }# u! k8 z  O3 _qm03=1/i*(acos(k*p003/A)-t21);
/ d! k, X1 f# }  B6 @+ m; hqm04=1/i*(acos(k*p004/A)-t21);
# e2 X; O  F- f6 Vqm1=linspace(qm03,qm04,100);
x12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
& g! S. X, S, Iz12=0*qm1;
5 k9 x, \2 ^' [# uplot(x12,y12)

5 Z  r: z/ o9 h4 k5 C8 g) `
3 o. `1 ~+ p' P/ T4 X& h: T/ PEF=[x12',y12',z12']
%save('EF.txt')
! O, p% m: x, Qend
7 v' L! l5 G& H
2 \- a' R( A  L1 R
2 c/ u4 ]7 v) I6 ]6 h%CalpaMFG.m
%原始不对称型线计算程序  [ x1, y1] = CalpMFG(100, 4, 6, 25)
$ t9 B& w; y8 p6 l. ^+ b, Dfunction [ x1, y1] = CalpMFG(A, Z1, Z2, R)
i=Z1/Z2; %齿数比
9 k" X* [' ~7 ER1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
2 m3 x7 i7 R) tR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
. |& I4 P) `) e/ I0 B: t: z2 ut=linspace(0,t,200);
. b9 i' E2 S1 ]0 J( A0 B; Rx1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
% ^( @; _! A8 D! I# qz1=0*t;
4 C; x! r+ [1 `6 Lplot(x1,y1)
6 o4 P/ {* c; a%
FG=[x1',y1',z1']
%save('FG.txt')
1 ^+ B* d1 D/ @9 Zend

7 g+ d9 o& C7 g0 g
% [ x1, y1] = CalpMFG(110, 5, 6,30)
8 j5 I3 M* t) {2 @3 o: ^" I* A

%CalpaMGH.m
%原始不对称型线计算程序
" l4 \3 C! a0 ofunction [ x11, y11] = CalpMGH(A, Z1, Z2, R)
- j7 m% v9 R  m  a$ `i=Z1/Z2; %齿数比
R1=(Z1/(Z1+Z2))*A; %阳转子节圆半径
% e4 t0 [, \2 ?" MR2=(Z2/(Z1+Z2))*A; %阴转子节圆半径
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%销齿圆弧的参数范围 在等腰三角形中求
%t=linspace(0,t,200);
7 k/ q6 u1 w2 |' X  F5 u%x1=R2-R*cos(t);y1=-R*sin(t);%销齿圆弧的参数方程  GF曲线段
%plot(x1,y1)

& i! j* e  Y; F  d% l" L
0 C  C+ E8 ?; s%第二曲线方程 GH  GH GH
9 F) X8 v0 i. ?/ ^6 C: Zb1=(R^2+R1^2+2*R*R1)^(1/2);  %这个地方第一次弄错了
9 T: M" W" {2 d4 {t1=0;
%x11=b1*cos(t1);y11=b1*sin(t1);%阳转子方程
( A0 w/ d# N  F1 K%t1=linspace(0,t1,100);
q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%转角参数
# b/ x+ o9 h1 Y0 v1 K%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%转角参数
%q=linspace(q1,q2,100);
$ @9 @2 p& f1 e+ T% Mk=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)
7 G* M% j& I) A* @


2 s3 p6 T! ?6 q%第三段曲线
%t21=acos((2*R2^2-R^2)/(2*R2^2));
x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C点横坐标   
, c7 Z% g" I8 s4 R' by0=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));
P001=(A^2+R2^2-2*A*R2)^(1/2);
! o/ B: O% l2 m* b5 |7 RP002=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)));
5 k8 ]' G  @& {) \& L8 d; rqm=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 l/ C# |7 i8 i4 kplot(x11,y11)
%
GH=[x11',y11',z11']
% j( B/ J% @2 T5 c7 M%save('GH.txt')
end
# T6 ]/ d9 u5 P0 ]- k
" y. u1 p, v7 L* L# f