wheelset_suspension.m

Matlab中实现轮对分析的程序以及相关的程序包,可以实现火车轮轨接触关系的分析

% M-file name: wheelset_suspension.m
% M-file type: Function file

% This function file computes the suspension forces acting on the wheelset due to the primary
% suspension components. The primary suspension consists of longitudinal and lateral springs
% and dampers. The inputs to the function are the wheelset lateral and yaw displacements and
% velocities. The outputs are the lateral suspension force and the vertical suspension moment
% on the wheelset.

% The outputs given by this function are called by the function file 'equations.m'.

function [Fsusp1,Msusp1]=wheelset_suspension(x1,x3,x2,x4)

% Parameters used for simulation

% kpx: Primary longitudinal stiffness coefficient (N/m)
% cpx: Primary longitudinal damping coefficient (N-sec/m)
% kpy: Primary lateral stiffness coefficient (N/m)
% cpy: Primary lateral damping coefficient (N-sec/m)
% dp: Half of lateral distance between primary longitudinal springs (m)

% Indicating the global nature of the variables This means that the value of the variables need
% not be specified in this function file. This value is automatically obtained from the main file
% 'single_wheelset.m'.

global kpx cpx kpy cpy dp;

% Computing the suspension forces and moments acting on the wheelset. These variables are
% passed to the function 'equations'.

Fsusp1=-2*kpy*x1-2*cpy*x3;
Msusp1=-2*dp^2*kpx*x2-2*dp^2*cpx*x4;