matlab ode45求解齿轮动力学,Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断 - Powered by Discuz!...

function doty=GearTorqueVibratory(t,y)

%转动惯量kg.m2

doty=zeros(10,2);

Jin=4.217823;

J1=0.008952085;

J2=0.190956155;

J3=0.020671561;

J4=0.766038806;

J5=0.090060;

J6=4.24048334;

J7=0.6044254;

J8=27.9071658;

Jout=58.764562;

%齿轮基圆半径m

R1=(63.171e-3)/2;

R2=(234.749e-3)/2;

R3=(97.732e-3)/2;

R4=(321.941e-3)/2;

R5=(137.755e-3)/2;

R6=(451.530e-3)/2;

R7=(213.959e-3)/2;

R8=(601.758e-3)/2;

%轴扭转刚度Nm.rad

k1=252851;

k2=1936800;

k3=6186712;

k4=33459824;

k5=10190124;

Tin=819.4;%输入扭矩N.m

Tout=91675;%输出轴扭矩N.m

% 传动轴阻尼计算

lanmda=0.01;

c1=2*lanmda*sqrt(k1/(1/Jin+1/J1));

<