Wheelrail_profiles.zip_S1002_railway vehicle_vehicle and rail_wh

% Rail Profile % reference from % The Dynamics of a Railway Vehicle on a Disturbed Track % page 80. Appendix A clc; %close all; clear; % function=z=s1002(y); wp=zeros(1001,3); for i=1:1001 y=-70+66/500*(i-1); if (y>60) z=-(y-60)-2.6357; elseif ((y<=60) & (y>32.15796)) z=1.364323640-0.066666667*y; elseif (y>-26.0) z=-3.358537058e-2*y + 1.565681624e-3*y^2 - 2.810427944e-5*y^3+... 5.844240864e-8*y^4 - 1.562379023e-8*y^5 + 5.309217349e-15*y^6-... 5.957839843e-12*y^7 + 2.646656573e-13*y^8; elseif (y>-35.0) z=-4.320221063e3 - 1.038384026e3*y - 1.065501873e2*y^2 -... 6.051367875*y^3 - 2.054332446e-1*y^4 - 4.169739389e-3*y^5 -... 4.687195829e-5*y^6 - 2.252755540e-7*y^7; elseif (y>-38.426669071) z=16.446 - sqrt(13*13-(y+26.210665)^2); elseif (y>-39.3764473993) z=-93.576667419 - 2.747477419*y; elseif (y>-49.662510381) z=8.834924130 + sqrt(400-(y+58.558326413)^2); elseif (y>-62.764705882) z=16 + sqrt(144-(y+55)^2); elseif (y>=-70) z=9.519259302 + sqrt(20.5^2-(y+49.5)^2); else z=0; end wp(i,1)=i; wp(i,2)=y; wp(i,3)=z; end Y=wp(:,2); Z=wp(:,3); Z2=-Z; Y2=-Y; figure(20) plot(Y,Z,'r.'); axis([-80 80 -10 30]); Title('Wheel profile using mathematical definition [S1002 wheel]'); legend('S1002 wheel profile') xlabel('Lateral distance[mm]'); ylabel('Profile height [mm]'); grid; figure(30) plot(Y,Z2,'r.'); axis([-80 80 -30 10]); Title('Wheel profile using mathematical definition [S1002 wheel]'); legend('S1002 wheel profile') xlabel('Lateral distance[mm]'); ylabel('Profile height [mm]'); grid; figure(40) plot(Y2,Z2,'r.'); axis([-80 80 -30 10]); Title('Wheel profile using mathematical definition [S1002 wheel]'); legend('S1002 wheel profile') xlabel('Lateral distance[mm]'); ylabel('Profile height [mm]'); grid;