dettimo.m

利用matalab编制的转子动力学计算软件。其中模型利用有限元建模

function t0=dettimo(lambda,L,A,Gs,E,Iy,rho)
% t0=dettimo(lambda,L,A,Gs,E,Iy,rho)
%
  
      s2 = lambda*lambda*exp(-sqrt(2.0)*L*(sqrt(-E*Iy*Gs*(lambda*rho*Iy*Gs+ ...
lambda*rho*Iy*E-sqrt(Iy*rho*(rho*Iy*lambda*lambda*Gs*Gs-2.0*rho*Iy*lambda*...
lambda*Gs*E+rho*Iy*lambda*lambda*E*E+4.0*E*Gs*Gs*A)))*lambda)+sqrt(-E*Iy*Gs*(...
lambda*rho*Iy*Gs+lambda*rho*Iy*E+sqrt(Iy*rho*(rho*Iy*lambda*lambda*Gs*Gs-2.0*...
rho*Iy*lambda*lambda*Gs*E+rho*Iy*lambda*lambda*E*E+4.0*E*Gs*Gs*A)))*lambda))/E/...
Iy/Gs);
      s3 = (exp(sqrt(2.0)*L*(2.0*sqrt(-E*Iy*Gs*(lambda*rho*Iy*Gs+lambda*rho*Iy*...
E+sqrt(Iy*rho*(rho*Iy*lambda*lambda*Gs*Gs-2.0*rho*Iy*lambda*lambda*Gs*E+rho*Iy*...
lambda*lambda*E*E+4.0*E*Gs*Gs*A)))*lambda)+sqrt(-E*Iy*Gs*(lambda*rho*Iy*Gs+...
lambda*rho*Iy*E-sqrt(Iy*rho*(rho*Iy*lambda*lambda*Gs*Gs-2.0*rho*Iy*lambda*...
lambda*Gs*E+rho*Iy*lambda*lambda*E*E+4.0*E*Gs*Gs*A)))*lambda))/E/Iy/Gs/2)-exp(...
1/E/Iy/Gs*sqrt(2.0)*sqrt(-E*Iy*Gs*(lambda*rho*Iy*Gs+lambda*rho*Iy*E-sqrt(Iy*rho...
*(rho*Iy*lambda*lambda*Gs*Gs-2.0*rho*Iy*lambda*lambda*Gs*E+rho*Iy*lambda*lambda...
*E*E+4.0*E*Gs*Gs*A)))*lambda)*L/2))*rho;
      s1 = s2*s3;
      s2 = (rho*Iy*lambda*lambda*Gs*Gs-2.0*rho*Iy*lambda*lambda*Gs*E+rho*Iy*...
lambda*lambda*E*E+4.0*E*Gs*Gs*A)*(-exp(sqrt(2.0)*L*(2.0*sqrt(-E*Iy*Gs*(lambda*...
rho*Iy*Gs+lambda*rho*Iy*E-sqrt(Iy*rho*(rho*Iy*lambda*lambda*Gs*Gs-2.0*rho*Iy*...
lambda*lambda*Gs*E+rho*Iy*lambda*lambda*E*E+4.0*E*Gs*Gs*A)))*lambda)+sqrt(-E*Iy...
*Gs*(lambda*rho*Iy*Gs+lambda*rho*Iy*E+sqrt(Iy*rho*(rho*Iy*lambda*lambda*Gs*Gs...
-2.0*rho*Iy*lambda*lambda*Gs*E+rho*Iy*lambda*lambda*E*E+4.0*E*Gs*Gs*A)))*lambda...
))/E/Iy/Gs/2)+exp(1/E/Iy/Gs*sqrt(2.0)*sqrt(-E*Iy*Gs*(lambda*rho*Iy*Gs+lambda*...
rho*Iy*E+sqrt(Iy*rho*(rho*Iy*lambda*lambda*Gs*Gs-2.0*rho*Iy*lambda*lambda*Gs*E+...
rho*Iy*lambda*lambda*E*E+4.0*E*Gs*Gs*A)))*lambda)*L/2))/Iy/(E*E)/(Gs*Gs);
      t0 = s1*s2;