Advertisement
andreacostabile

trasmissibilita_vibraz

Dec 5th, 2022 (edited)
855
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 1.57 KB | Science | 0 0
  1. close all; clear
  2.  
  3. om=linspace(0,5,5000);
  4. T=@(x)sqrt((1+(2*x*om).^2)./((1-om.^2).^2+4*(x^2)*om.^2));
  5. Ta=@(x,y)sqrt((1+(2*x*y)^2)./((1-y^2)^2+4*(x^2)*y^2));
  6.  
  7. %zita=0.01
  8. T1=T(0.01);
  9. pt1=plot(om,T1,"g","LineWidth",1.5);
  10. hold on
  11. %zita=0.05
  12. T2=T(0.05);
  13. pt2=plot(om,T2,"r","LineWidth",1.5);
  14.  
  15. %zita=0.1
  16. T3=T(0.1);
  17. pt3=plot(om,T3,"b","LineWidth",1.5);
  18.  
  19. %non smorzato
  20. T0=T(0);
  21. pt0=plot(om,T0,"k","LineWidth",1.5);
  22.  
  23. plot(om,zeros(1,length(om)),"k","LineWidth", 0.5)
  24. plot(sqrt(2)*(ones(1,100)),linspace(-20,20),"r--")
  25.  
  26. om_acc=25.133/11.1;     % w/wn(acciaio)
  27. pt_acc=plot(om_acc,Ta(0,om_acc),"+k","MarkerSize",10,"LineWidth",1.5);
  28.  
  29. om_zav=25.133/9.27;     % w/wn(acciaio con zavorra)
  30. pt_zav=plot(om_zav,Ta(0,om_zav),"ok","MarkerSize",10,"LineWidth",1.5);
  31.  
  32. om_prog=25.133/7.578;   % w/wn(progetto)
  33. pt_prog=plot(om_prog,Ta(0,om_prog),"xk","MarkerSize",10,"LineWidth",1.5);
  34.  
  35. om_aria=25.133/6.76;    % w/wn(aria)
  36. pt_aria0=plot(om_aria,Ta(0,om_aria),"*k","MarkerSize",10,"LineWidth",1.5);
  37. pt_aria1=plot(om_aria,Ta(0.01,om_aria),"*g","MarkerSize",10,"LineWidth",1.5);
  38. pt_aria2=plot(om_aria,Ta(0.05,om_aria),"*r","MarkerSize",10,"LineWidth",1.5);
  39. pt_aria3=plot(om_aria,Ta(0.1,om_aria),"*b","MarkerSize",10,"LineWidth",1.5);
  40.  
  41.  
  42. legd=legend([pt1 pt2 pt3 pt0 pt_prog pt_acc pt_zav pt_aria0], "z=0.01","z=0.05","z=0.1","z=0","Progetto","5xB40 305", "5xB40 305 con zavorra" ,"WBZ 500","Location","NorthEast");
  43. legd.FontSize=18; legd.FontName="FixedWidth";
  44. title("Trasmissibilità"), xlabel("w/wn"), ylabel("T")
  45. axis([0 5 0 6])
  46. grid off
  47.  
  48. pause
  49.  
  50. axis([3.7 3.735 0.06 0.11])
  51. grid on
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement