matlab ode45求解齿轮动力学,使用ODE45求解齿轮系统动力学方程后结果发散

EDA365欢迎您登录！

您需要 登录 才可以下载或查看，没有帐号？注册

x

1.1参数定义及动力学方程降阶

: g) c6 b" p' |  p7 _8 Gfunction [dx,ff1,ff2]=myfun(t,x)

! r  a6 [+ c' St

) k/ o7 {- i2 fbeita=26;

% ^- K% \' g# L, @" p4 q* u9 Umn=0.004;

' a5 ]3 t; N" v$ m+ n0 t3 mz1=46;

' A/ Z  R- B; K$ J3 P2 l, \/ p+ cz2=43;                       " p$ b' N. S- c, k' `3 e( ^5 D& H/ q% @

z3=122;

. l( l9 R8 q2 e  @" AT_in=200;

3 v; }) f% Z( F1 E: i) _( j: oT_out=80;

" L& g0 S! k2 Aroug1=7.8E3;            ) y# u$ V. z- S0 j

roug2=7.8E3;

% o- y( ]0 s7 e% Broug3=7.8E3;            / e% ]1 c' [1 j1 ?0 J; M3 [

alphan=20;                                                   7 f- R1 L9 y: H  j

alphat=atand(tand(alphan)/cosd(beita));      & F2 j" c* [8 w7 j

d1=z1*mn/cosd(beita)/1000;

0 E- w# H( n  K  fdb1=d1*cosd(alphat)/1000;

: X7 q4 Y- _2 }8 G' Ld2=z2*mn/cosd(beita)/1000;

2 _+ M6 c' Y% {0 bdb2=d2*cosd(alphat)/1000;

- r0 A  X( W/ L: D2 m6 a, R5 vd3=z3*mn/cosd(beita)/1000;

4 D2 Y4 @. d9 X6 u6 \  idb3=d3*cosd(alphat)/1000;

5 T/ C4 a/ ^' h! w; e. n; ]4 Jbp1=116/1000;

; u  V. ?9 I" g0 ], n7 X5 o6 F9 ~bp2=116/1000;