久久久国产一区二区_国产精品av电影_日韩精品中文字幕一区二区三区_精品一区二区三区免费毛片爱
機械社區
標題:
參數化計算 原始不對稱型線
[打印本頁]
作者:
shouce
時間:
2015-12-22 16:32
標題:
參數化計算 原始不對稱型線
%CalpaMEF.m
+ X4 D b0 C! b# `; P
%原始不對稱型線計算程序 [ x12, y12] = CalpMEF(100, 4, 6, 25)
% |! d8 P+ V. ]- ~. m6 X9 F
function [ x12, y12] = CalpMEF(A, Z1, Z2, R)
5 ^2 t+ A( [# P4 Y6 h. g- v7 i
i=Z1/Z2; %齒數比
* |4 G+ \7 \9 I: e- N
R1=(Z1/(Z1+Z2))*A; %陽轉子節圓半徑
. d, _- Q/ o6 ?! k$ _- V" {9 v
R2=(Z2/(Z1+Z2))*A; %陰轉子節圓半徑
7 `) a& t- x9 ~8 g* C* J
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數范圍 在等腰三角形中求
: A# k! V+ o) {
%t=linspace(0,t,200);
5 f f7 |7 z% ]" l1 D n
%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數方程 GF曲線段
1 [9 E5 S0 f- Y
%plot(x1,y1)
& z* U% l7 }/ S
4 [7 G! F4 O C7 Z
/ U, h* n0 D+ O6 g' [
%第二曲線方程 GH GH GH
) p: K7 h# g( z& z! }
%b1=(R^2+R1^2+2*R*R1)^(1/2); %這個地方第一次弄錯了
6 v2 @; ^! n, ?; A7 Y2 K
%t1=0;
9 \6 \3 V( W0 n" R
%x11=b1*cos(t1);y11=b1*sin(t1);%陽轉子方程
/ z! g+ ?( m( |/ c, l7 a. @
%t1=linspace(0,t1,100);
2 C7 l4 W5 M8 V% \" s0 r) u7 n0 ^
%q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉角參數
$ [+ U; m; s7 \
%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%轉角參數
0 p5 v N# I" Z) I, F. s! L
%q=linspace(q1,q2,100);
6 F/ [' p r9 ^& ]+ l
k=i+1;
3 W. X \% ^& G" n7 L
%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
; h3 z8 J/ I2 P& K% E1 b
%plot(x22,y22)
]* E1 j7 @3 n* s
3 _0 B1 R _( J
. o& W: a$ L; D$ c4 o$ h
%t21=acos((2*R2^2-R^2)/(2*R2^2));
. Y% b* k* C: K6 o/ z5 ?0 x
%x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C點橫坐標
# A- r! i C& k5 E5 X
%y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C點縱坐標
$ @9 p; e( r& Y0 r% V, Q0 Y
%cp=((x0-R2)^2+y0^2)^(1/2);%計算線段長度
# a/ h0 e/ S% [# t* }: Q1 }
%t22=acos((2*R2^2-cp^2)/(2*R2^2));
- S* Y/ y" X- ^0 ?' ^
%P001=(A^2+R2^2-2*A*R2)^(1/2);
2 o& a' [6 J" R$ K5 O/ R
%P002=b1;
6 p+ u6 q, |1 t8 Q
%qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號輸錯
8 J- c4 ?8 S w5 ?
%qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
+ B! C' }; b6 a7 X
%qm=linspace(qm01,qm02,100);
! w& o# L6 G3 N& N3 O+ M: K
%x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
* l# X$ f( N+ N/ N5 {* A
%plot(x11,y11)
5 n+ X. {) B1 x( t4 v. g3 y
6 y1 w" f, v! v9 Y& [* r: \
* B7 | X; ?; X# z6 t# M
%第二曲線方程 EF EF EF
: c8 }, L$ P z/ E. {* T# K
t21=acos((2*R2^2-R^2)/(2*R2^2));
4 K/ W3 m% d9 j$ H$ K3 G
p003=R2*cos(t21); %有點問題% 為什么是這個樣的？
9 @& }. J8 Y3 d* g
p004=R2;
! A$ R x3 ~& ]1 W
%PP=linspace(p003,p004,100);
. m2 S) l7 {8 T: T3 c: M$ m
qm03=1/i*(acos(k*p003/A)-t21);
7 V! [% H4 Q( c8 u5 g3 o% g
qm04=1/i*(acos(k*p004/A)-t21);
! ]$ V& C, C7 |$ n. ?1 u
qm1=linspace(qm03,qm04,100);
0 p. ]1 c7 \* I5 d. P, w4 R' H
x12=A-(A*cos(qm1)-R2*cos(t21+k*qm1));y12=A*sin(qm1)-R2*sin(t21+k*qm1); %方程
" n) T- O; g, |, ] C" I+ V
z12=0*qm1;
% K) N" e2 v, O! h
plot(x12,y12)
- \. s% F! ]1 T: F* I2 q
3 X: ~2 T: N' v
. t1 k0 G& {7 J: j% b
EF=[x12',y12',z12']
( ~! e2 \3 y6 x }
%save('EF.txt')
8 p: l$ I) ~4 Y
end
\6 |0 j/ ]: q+ X
7 g0 }- p* h- w: Z+ Y, \8 c8 r
' T% L1 l5 n& @& W
%CalpaMFG.m
9 u2 L2 V5 ~( u
%原始不對稱型線計算程序 [ x1, y1] = CalpMFG(100, 4, 6, 25)
" b7 L8 e% V) p+ t
function [ x1, y1] = CalpMFG(A, Z1, Z2, R)
: u& l) o {/ c& P
i=Z1/Z2; %齒數比
! ?1 d( ?3 F1 o: o, `# I5 P, y* v3 N
R1=(Z1/(Z1+Z2))*A; %陽轉子節圓半徑
3 z4 C4 `0 d) K5 f
R2=(Z2/(Z1+Z2))*A; %陰轉子節圓半徑
. d7 W Q( q3 C1 K+ w
t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數范圍 在等腰三角形中求
( A4 I) S) `/ H) a: K9 I3 g
t=linspace(0,t,200);
' H7 x: x! F! Q; d1 g/ p/ K( W
x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數方程 GF曲線段
0 s: `' w/ W" _; f- l
z1=0*t;
3 ` k% S% |8 K) z5 T
plot(x1,y1)
- S3 M, k; V1 y
%
3 k, m2 J) @, {7 k% U2 {
FG=[x1',y1',z1']
6 G" Y# R* _! V. b' M
%save('FG.txt')
5 d" j' y$ e, S) K6 q
end
8 J$ [& f. \" s' H
+ d7 U) F( E: r$ d1 s, _
2 \: w l" v/ z- H% q
% [ x1, y1] = CalpMFG(110, 5, 6,30)
; U; f; W: d" j! [% q
2 Z% ]' f1 b% L8 d) ?
# k& E% ?+ b; T" j8 u6 ]
%CalpaMGH.m
# {4 g! r& T, _: x
%原始不對稱型線計算程序
% ]) K( a* I |9 }# G8 q: J) w
function [ x11, y11] = CalpMGH(A, Z1, Z2, R)
+ ]3 ]& m- g& x% |
i=Z1/Z2; %齒數比
+ d+ i8 a' z0 ~, T
R1=(Z1/(Z1+Z2))*A; %陽轉子節圓半徑
1 \5 E: Q9 m) w
R2=(Z2/(Z1+Z2))*A; %陰轉子節圓半徑
. Y% I9 H7 l1 k+ C9 g4 \6 i
%t=(pi-acos((2*R2^2-R^2)/(2*R2^2)))/2;%銷齒圓弧的參數范圍 在等腰三角形中求
& f# n% z* g, k+ a- C- L& R
%t=linspace(0,t,200);
' C" ]7 P7 u( n; z
%x1=R2-R*cos(t);y1=-R*sin(t);%銷齒圓弧的參數方程 GF曲線段
' X* |) w' \% e$ D& n% Y7 a
%plot(x1,y1)
# I6 z6 A" c. I W+ l
/ O! _" G( _# K" i3 U' n/ `
0 }. n' n- s" V9 L1 i h
%第二曲線方程 GH GH GH
# [, O. P4 {( F" S# c) L. l
b1=(R^2+R1^2+2*R*R1)^(1/2); %這個地方第一次弄錯了
9 b2 v a5 a* h: z/ X
t1=0;
! u' g4 N, m9 `: f9 ?4 }9 P
%x11=b1*cos(t1);y11=b1*sin(t1);%陽轉子方程
7 w" a/ `! e9 F, r
%t1=linspace(0,t1,100);
# m2 {7 P& S7 t
q1=0-acos((A^2+b1^2-R2^2)/(2*A*b1));%轉角參數
$ K( y. w( q9 p, ?4 u5 P, \
%q2=0-acos((A^2+b1^2-(R^2+R2^2-2*R*R2))/(2*A*b1));%轉角參數
; {, ^) m& q2 ]: x8 D7 l
%q=linspace(q1,q2,100);
& h9 L, V, ?) f. l: ]
k=i+1;
F% }. N' A0 ~
%x22=A*cos(i*q)-b1*cos(t1-k*q);y22=A*sin(i*q)+b1*sin(t1-k*q);%曲線方程
' C z9 x9 ?, }0 k0 f. ]! t
%plot(x22,y22)
, i9 M0 w5 R% N) ~" y8 F2 R3 W
' S/ J/ g$ y3 Z" [+ e( H m3 O
+ b: y1 x: ~2 k( C
9 U" d0 {% a, n; \9 i
%第三段曲線
3 ?* U1 H0 r3 I9 A" A( U9 @
%t21=acos((2*R2^2-R^2)/(2*R2^2));
2 `1 ?! n% m: ~+ H& x. S+ b$ p
x0=A*cos(i*q1)-b1*cos(t1-k*q1);%C點橫坐標
( U; j. r2 S: l: h. E% z9 o5 \
y0=A*sin(i*q1)+b1*sin(t1-k*q1);%C點縱坐標
# F6 J/ s2 P4 E M
cp=((x0-R2)^2+y0^2)^(1/2);%計算線段長度
+ B- S6 g K% c2 d3 |
t22=acos((2*R2^2-cp^2)/(2*R2^2));
+ T# k' J" h* z; k0 h2 g1 T6 O, H
P001=(A^2+R2^2-2*A*R2)^(1/2);
( i* N- Z- W/ ?6 _# E
P002=b1;
# W9 S6 H7 A0 ^' ^ `. u8 k
qm01=1/i*(t22-acos((A^2+R2^2-P001^2)/(2*A*R2))); %第一次在這兒括號輸錯
/ h3 N9 t9 ]; w- |0 c
qm02=1/i*(t22-acos((A^2+R2^2-P002^2)/(2*A*R2)));
W. y/ {; T; i2 V
qm=linspace(qm01,qm02,100);
% u, }9 _! k& K! [) W, Q
x11=A-(A*cos(qm)-R2*cos(t22-k*qm));y11=A*sin(qm)+R2*sin(t22-k*qm); %方程
0 p: V3 B' z# z9 P2 x
z11=0*qm;
, C2 E( S! E. w% D5 h6 l# M
plot(x11,y11)
r/ Z5 q, [% G2 S
%
/ T9 |# H/ w) l& J0 A. g# |* D
GH=[x11',y11',z11']
8 w& n; ?. k- r& T9 O7 y
%save('GH.txt')
: {& O& F% y- L$ b J* @
end
5 S6 V, V1 b9 l: y$ Q
; ]' H$ j Z5 w. l7 m9 u
7 |9 P5 r- H5 W- _% c/ o
# \; ]4 R3 S/ h, y; H/ W! b
4 w( ^& @0 V Z1 ^( K* J1 Y
歡迎光臨 機械社區 (http://www.ytsybjq.com/)
Powered by Discuz! X3.5