CN103630378B - Bullet train comfort of passenger detection method - Google Patents

Abstract

The present invention relates to a kind of bullet train comfort of passenger detection method, including: A. uses pseudo-excitation method to generate the virtual rail irregularity for flexible car body random dynamic analysis；B. flexible car body random vibration frequency domain is analyzed, draws the random response of flexible car body；C. according to the random response drawn, different ride comfort standard is utilized to carry out the comfort of passenger of flexible car body calculating assessment.The present invention uses the large-scale FEM model of flexible car body to be estimated the comfort level of passenger on flexible car body diverse location, evaluation process use pseudo-excitation method and the Symplectic Algebra theoretical, there is provided provide powerful support for for promoting high-speed train body structure designed capacity, and for further vehicle and track structure being optimized, the research such as control established theory and the technology of numerical simulation of high efficient and reliable.

Description

Bullet train comfort of passenger detection method

Technical field

The present invention relates to the detection method of a kind of bullet train flexible car body comfort of passenger, for promoting height Speed train body structure designed capacity provides and provides powerful support for, and is to tie vehicle and track further Reliable theory and technology of numerical simulation have been established in the researchs such as structure is optimized, control.

Background technology

Railway is a kind of ground transportation system economic, efficient.When train orbits, vehicle Track Coupling System can be encouraged by from track irregularity, this excitation have significantly with Machine, is the main driving source causing vibration of rolling stock to generally acknowledge, to passenger during train operation Comfort level produces strong influence.And bullet train comfort of passenger is estimated, elastic The calculating of the computational analysis of car body random vibration, especially its power spectrum, occupies core status.

The most fewer to the research of stock rail coupled system random vibration, computation model is the most only Being limited to fairly simple Rigid Body Dynamics Model, its main difficulty is the free degree of huge structure Poor efficiency with tradition random vibration computational methods.Use conventional finite element method to set up and have limit for length's rail During road model, comprise before and after vehicle tens and just can substantially eliminate rail boundary across the length of sleeper The impact of condition.Therefore, even if vehicle uses the Rigid Body Dynamics Model that the free degree is relatively low, coupling The free degree of assembly system is the highest.And for such system, use tradition random vibration reason During opinion direct solution random response, need to solve frequency response function matrix, and carry out the company of big matrix Multiplication；And use Monte-Carlo method to be simulated, then need to choose substantial amounts of irregularity sample Basis and less integration step, and longer time history is to guarantee the precision solved.In addition real The body construction on border is also considerably complicated, uses rigid model to carry out calculating and can cause bigger Error, and the situation that car body diverse location vibrates cannot be reflected exactly.If car to be considered The comfort level of interior diverse location passenger, then need to set up detailed car body FEM model, its meter Calculation amount is difficult to accept.

Summary of the invention

Present invention is primarily aimed at solution the problems referred to above and deficiency, it is provided that a kind of bullet train passenger is relaxed Appropriateness detection method, can be quickly obtained bullet train elasticity car by efficient computational methods Every comfort level index of body passenger.

For achieving the above object, the technical scheme is that

A kind of bullet train comfort of passenger detection method, described method includes:

A. pseudo-excitation method is used to generate the virtual rail irregularity for flexible car body random dynamic analysis ；

B. flexible car body random vibration frequency domain is analyzed, draws the random response of flexible car body；

C. according to the random response drawn, the different ride comfort standard occupant comfort to flexible car body is utilized Degree carries out calculating assessment；

Further, described employing pseudo-excitation method generates virtual rail irregularity, its mathematic(al) representation For:

${\tilde{r}}_j(\mathrm{\ω},t)={\left\{\begin{array}{cccc}{e}^{-\mathrm{i\ω}{t}_1}& e^{-\mathrm{i\ω}{t}_2}& e^{-\mathrm{i\ω}{t}_3}& e^{-\mathrm{i\ω}{t}_4}\end{array}\right\}}^T\sqrt{S_j\left(\mathrm{\ω}\right)}{e}^{\mathrm{i\ωt}};(j=v,a,c)$

Wherein, t₁, t₂, t₃, t₄Be four wheel to time lag；

J=v, a, c represent high and low direction and level three class track irregularity respectively；

ω is time circular frequency.

Further, it is considered to 4 wheels by the arbitrary excitation of same track, four wheels to time lag mathematics Expression formula is as follows:

$t_1=0;t_2=2l_1/v;t_3=2l_2/v;t_4=2(l_1+l_2)/v$

Wherein, l₁, l₂It is respectively semifixed wheelbase and half bogie pivot center of vehicle.

Further, in stepb, virtual rail irregularity is acted on stock rail coupled system fortune Dynamic equation, solves the random response of coupled system.

Further, in described stock rail coupled system, flexible car body uses large-scale FEM model , track structure uses the Symplectic Algebra Method Modeling.

Further, according to vehicle and the actual geometry of track structure and physical parameter, coupled system is set up Steady motion equation, its mathematic(al) representation is: $K_d{\tilde{u}}_j=r_j\sqrt{S_j}{e}^{\mathrm{i\ωt}}$

Wherein, K_dDynamic Stiffness Matrix for stock rail coupled system；

For the virtual displacement response vector of coupled system under dummy excitation effect；

r_jFor corresponding loading coefficient vector.

Further, according to dummy excitation law theory, coupled system is solved by following mathematic(al) representation Displacement and acceleration power spectrum:

$s_{\mathrm{uu}}\left(\mathrm{\ω}\right)=\underset{j=v,a,c}{\mathrm{\Σ}}{{\tilde{u}}_j}^{*}(\mathrm{\ω},t){{\tilde{u}}_j}^T(\mathrm{\ω},t);s_{\stackrel{..}{u}\stackrel{..}{u}}\left(\mathrm{\ω}\right)={\mathrm{\ω}}^4{s}_{\mathrm{uu}}\left(\mathrm{\ω}\right)$

Therefrom choose the corresponding free degree of flexible car body, obtain the random response of flexible car body.

Further, the seat chosen inside chassis of vehicle body on diverse location on flexible car body FEM model The corresponding node in chair installation place is as research object.

Further, different ride comfort standard is utilized to refer specifically to be utilized respectively China GB5 described in step C 599-85 marks Standard, railway systems TB/T2630-93 standard, the ISO2631-1 standard of International Organization for standardization and The UIC513 standard of International Union of Railways, the comfort level index to flexible car body diverse location passenger Calculate.

Further, the position calculating comfort of passenger index is chosen position near the window, centre position and leans on away Position, corridor.

To sum up content, bullet train comfort of passenger detection method of the present invention, with prior art Compare, have the advantage that

Use flexible car body large-scale FEM model that the comfort level of passenger on flexible car body diverse location is entered Row assessment, uses pseudo-excitation method to be converted into accordingly by three class track irregularities in evaluation process Simple harmonic quantity dummy excitation；The Symplectic Algebra theory is used to set up the low degree-of-freedom of stock rail coupled system The equation of motion；GB5599-85, TB/T2630-93, ISO2631-1 and UIC is used under frequency domain Comfort of passenger on flexible car body diverse location is referred in territory by domestic and international ride comfort standard carrier frequency such as 513 grade Mark calculates.

The present invention provides powerful support for for promoting the offer of high-speed train body structure designed capacity, and is further Be optimized vehicle and track structure, control etc. studies theory and the number having established high efficient and reliable Value emulation technology.

Accompanying drawing explanation

Fig. 1 is that stock rail coupled system analyzes model；

Fig. 2 is track structure model；

Fig. 3 is infinite long period chain structure；

Fig. 4 is domestic CR H type bullet train head car car body；

Fig. 5 is to calculate the vertical of three groups of focuses and horizontal comfort index with " GB5599-85 standard "；

Fig. 6 is to calculate the vertical of three groups of focuses and horizontal comfort index with " TB/T2630-93 standard " ；

Fig. 7 is the comfort index calculating three groups of focuses with " ISO-2631-1 standard "；

Fig. 8 is the comfort index calculating three groups of focuses with " UIC513 standard ".

Detailed description of the invention

With detailed description of the invention, the present invention is described in further detail below in conjunction with the accompanying drawings:

As shown in Figure 1 to Figure 3, bullet train comfort of passenger detection method of the present invention, including Following steps:

Step A, employing pseudo-excitation method generate the void for bullet train flexible car body random dynamic analysis Intend track irregularity；

Consider 4 wheels arbitrary excitation by same track, four wheels to time lag mathematic(al) representation such as Under:

$t_1=0;t_2=2l_1/v;t_3=2l_2/v;t_4=2(l_1+l_2)/v---\left(1\right)$

Wherein l₁, l₂It is respectively semifixed wheelbase and half bogie pivot center of vehicle.

Assume that three class track irregularities are zero-mean Stationary Gauss Random process, and there is no correlation each other , its power spectrum is respectively S_v(ω)、S_a(ω) and S_c(ω).Wherein, ω is time circular frequency 。

According to dummy excitation law theory, generating corresponding virtual rail irregularity, its mathematic(al) representation is :

${\tilde{r}}_j(\mathrm{\ω},t)={\left\{\begin{array}{cccc}{e}^{-\mathrm{i\ω}{t}_1}& e^{-\mathrm{i\ω}{t}_2}& e^{-\mathrm{i\ω}{t}_3}& e^{-\mathrm{i\ω}{t}_4}\end{array}\right\}}^T\sqrt{S_j\left(\mathrm{\ω}\right)}{e}^{\mathrm{i\ωt}};(j=v,a,c)---\left(2\right)$

Wherein j=v, a, c represent high and low direction and level three class track irregularity respectively.

Virtual rail irregularity is acted on the stock rail coupled system equation of motion, coupling can be solved The virtual responsive of assembly system.

Step B, car body random vibration frequency domain is analyzed, draws the random response of car body；

Use stock rail coupling system model as shown in Figure 1, wherein flexible car body structure finite element Model uses commercial finite element software to generate, and calculates its Free Modal, synthesizes flexible car body mould The state equation of motion, its mathematic(al) representation is as follows:

$M_c{\stackrel{..}{q}}_c+C_c{\stackrel{.}{q}}_c+K_c{q}_c=C_a^{\mathrm{ct}}{\stackrel{.}{u}}_t+K_a^{\mathrm{ct}}{u}_t---\left(3\right)$

Wherein: q_cModal displacement vector for car body；

M_c、C_cAnd K_cIt is respectively the modal mass of car body, damping and stiffness matrix；

With For corresponding loading coefficient matrix.

Bogie uses Rigid Body Dynamics Model, any dynamic method can be used to model, its motion side Journey can be expressed as following form:

$M_t{\stackrel{..}{u}}_t+C_t{\stackrel{.}{u}}_t+K_t{u}_t=C_a^{\mathrm{ct}}{q}_t+Q_{\mathrm{tw}}^C{\stackrel{.}{u}}_w+Q_{\mathrm{tw}}^K{u}_w---\left(4\right)$

Wherein, For wheel to motion vector；

M_t, C_tThe respectively quality of bogie, damping and stiffness matrix；

、 、 With For corresponding loading coefficient matrix.

Car body is coupled with the bogie equation of motion, sets up the equation of motion of vehicle:

$\left[\begin{array}{cc}{\stackrel{\‾}{M}}_c& \\ & M_t\end{array}\right]\left\{\begin{array}{c}{\stackrel{..}{q}}_t\\ {\stackrel{..}{u}}_t\end{array}\right\}+\left[\begin{array}{cc}{\stackrel{\‾}{C}}_c& -{\stackrel{\‾}{C}}_a^{\mathrm{ct}}\\ -{\stackrel{\‾}{C}}_a^{\mathrm{tc}}& C_t\end{array}\right]\left\{\begin{array}{c}{\stackrel{.}{q}}_c\\ {\stackrel{.}{u}}_t\end{array}\right\}+\left[\begin{array}{cc}{\stackrel{\‾}{K}}_c& -{\stackrel{\‾}{K}}_a^{\mathrm{ct}}\\ -{\stackrel{\‾}{K}}_a^{\mathrm{tc}}& K_t\end{array}\right]\left\{\begin{array}{c}{q}_c\\ u_t\end{array}\right\}=\left[\begin{array}{c}0\\ Q_{\mathrm{tw}}^C\end{array}\right]{\stackrel{.}{u}}_w+\left[\begin{array}{c}0\\ Q_{\mathrm{tw}}^K\end{array}\right]u_w---\left(5\right)$

It is abbreviated as:

$M_v{\stackrel{..}{u}}_v+C_v{\stackrel{.}{u}}_v+K_v{u}_v=Q_{\mathrm{vt}}^C{\stackrel{.}{u}}_w+Q_{\mathrm{vt}}^K{u}_w---\left(6\right)$

And it is rewritten as form stable:

$(K_v+\mathrm{i\ω}{C}_v-{\mathrm{\ω}}^2{M}_v)u_v=K_v^d{u}_v=(\mathrm{i\ω}{Q}_{\mathrm{vt}}^C+Q_{\mathrm{vt}}^K)u_w---\left(7\right)$

Wherein With Represent respectively vehicle and wheel under jth class track irregularity effect to virtual displacement.

Track structure is considered as three leafing scatterplot supporting constructions as shown in Figure 2, and regards it as shown in Figure 3 Infinite Long cycle chain structure, with stress track minor structure as research object, sets up the equation of motion:

$M_e{\stackrel{..}{u}}_e+C_e{\stackrel{.}{u}}_e+K_e{u}_e=f_e\left(t\right)+f_b---\left(8\right)$

Wherein, u_eFor minor structure motion vector；

M_e、C_eAnd K_eIt is respectively track mega-structure mass, damping and stiffness matrix；

f_bThe load vectors of this minor structure it is applied to for closing on minor structure.

It is rewritten as form stable:

$(-{\mathrm{\ω}}^2{M}_e+\mathrm{i\ω}{C}_e+K_e)u_e=K_e^d{u}_e=f_e\left(t\right)+f_b---\left(9\right)$

Its block form mathematical expression is stated as follows:

$\left[\begin{array}{ccc}{K}_{\mathrm{aa}}^0& K_{\mathrm{ab}}^0& K_{\mathrm{ai}}^0\\ K_{\mathrm{ba}}^0& K_{\mathrm{bb}}^0& K_{\mathrm{bi}}^0\\ K_{\mathrm{ia}}^0& K_{\mathrm{ib}}^0& K_{\mathrm{ii}}^0\end{array}\right]\left\{\begin{array}{c}{u}_a\\ u_b\\ u_i\end{array}\right\}=\left\{\begin{array}{c}{f}_{\mathrm{ae}}\\ f_{\mathrm{be}}\\ f_{\mathrm{ie}}\end{array}\right\}+\left\{\begin{array}{c}{p}_a\\ -p_b\\ 0\end{array}\right\}---\left(10\right)$

Wherein, u_aAnd u_bRepresent the motion vector of the two side outlet frees degree respectively；

u_iRepresent the motion vector of internal degree of freedom；

f_ae,f_beAnd f_ieIt is respectively outer force vector suffered by both sides Egress node and internal node；

p_aAnd p_bIt is to close on minor structure to be applied to the nodal force vector of this minor structure.

According to following mathematic(al) representation synthesis the Symplectic Algebra matrix S (ω):

$\begin{array}{c}{K}_{\mathrm{aa}}=K_{\mathrm{aa}}^0-K_{\mathrm{ai}}^0{\left(K_{\mathrm{ii}}^0\right)}^{-1}{K}_{\mathrm{ia}}^0;K_{\mathrm{ab}}=K_{\mathrm{ab}}^0-K_{\mathrm{ai}}^0{\left(K_{\mathrm{ii}}^0\right)}^{-1}{K}_{\mathrm{ib}}^0\\ K_{\mathrm{ba}}={\left(K_{\mathrm{ab}}\right)}^T;K_{\mathrm{bb}}=K_{\mathrm{bb}}^0-K_{\mathrm{bi}}^0{\left(K_{\mathrm{ii}}^0\right)}^{-1}{K}_{\mathrm{ib}}^0\\ S_{\mathrm{aa}}=-{\left(K_{\mathrm{ab}}\right)}^{-1}{K}_{\mathrm{aa}};S_{\mathrm{ab}}={\left(K_{\mathrm{ab}}\right)}^{-1}\\ S_{\mathrm{ba}}=-K_{\mathrm{ba}}+K_{\mathrm{bb}}{\left(K_{\mathrm{ab}}\right)}^{-1}{K}_{\mathrm{aa}};S_{\mathrm{bb}}=K_{\mathrm{bb}}^0-K_{\mathrm{bi}}^0{\left(K_{\mathrm{ii}}^0\right)}^{-1}{K}_{\mathrm{ib}}^0\\ S\left(\mathrm{\ω}\right)=\left[\begin{array}{cc}{S}_{\mathrm{aa}}& S_{\mathrm{ab}}\\ S_{\mathrm{ba}}& S_{\mathrm{bb}}\end{array}\right]\end{array}---\left(11\right)$

Solve the characteristic vector of the Symplectic Algebra matrix, and it write as block form:

Wherein, , i=1,2 ..., 2n is the characteristic vector of symplectic matrix S (ω), and n is the freedom of minor structure Degree.

Synthesize by following mathematic(al) representation:

$P_{\mathrm{\α}}=N_a{X}_a^{-1};P_{\mathrm{\β}}=-N_b{X}_b^{-1}---\left(13\right)$

Then stress minor structure steady motion equation is rewritten as following form:

$\left[\begin{array}{ccc}{K}_{\mathrm{aa}}^0+P_{\mathrm{\β}}& K_{\mathrm{ai}}^0& K_{\mathrm{ab}}^0\\ K_{\mathrm{ia}}^0& K_{\mathrm{ii}}^0& K_{\mathrm{ib}}^0\\ K_{\mathrm{ba}}^0& K_{\mathrm{bi}}^0& K_{\mathrm{bb}}^0+P_{\mathrm{\α}}\end{array}\right]\left\{\begin{array}{c}{u}_a\\ u_i\\ u_b\end{array}\right\}=\left\{\begin{array}{c}{f}_{\mathrm{ae}}\\ f_{\mathrm{be}}\\ f_{\mathrm{ie}}\end{array}\right\}---\left(14\right)$

Then the response of kth adjacent tracks minor structure uses following mathematic(al) representation to calculate:

$u_{\mathrm{ek}}=\left\{\begin{array}{cc}\left[\begin{array}{c}{X}_b{\mathrm{\&mu;}}^k{X}_b^{-1}\\ X_b{\mathrm{\&mu;}}^{k-1}{X}_b^{-1}\\ -{\left(K_{\mathrm{ii}}^0\right)}^{-1}{K}_{\mathrm{ia}}^0{X}_b{\mathrm{\&mu;}}^k{X}_b^{-1}-{\left(K_{\mathrm{ii}}^0\right)}^{-1}{K}_{\mathrm{ib}}^0{X}_b{\mathrm{\&mu;}}^{k-1}{X}_b^{-1}\end{array}\right]u_a& k<0\\ \left[\begin{array}{c}{X}_a{\mathrm{\&mu;}}^{k-1}{X}_a^{-1}\\ X_a{\mathrm{\&mu;}}^k{X}_a^{-1}\\ -{\left(K_{\mathrm{ii}}^0\right)}^{-1}{K}_{\mathrm{ia}}^0{X}_a{\mathrm{\&mu;}}^{k-1}{X}_a^{-1}-{\left(K_{\mathrm{ii}}^0\right)}^{-1}{K}_{\mathrm{ib}}^0{X}_a{\mathrm{\&mu;}}^k{X}_a^{-1}\end{array}\right]u_b& k>0\end{array}\right.---\left(15\right)$

The equation of motion of vehicle and four stress track minor structures is coupled together, and considers in above formula each Interaction between minor structure, obtains the steady motion equation of coupled system:

$K_d{\tilde{u}}_j=r_j\sqrt{S_j}{e}^{\mathrm{i\&omega;t}}---\left(16\right)$

Wherein, flexible car body uses large-scale FEM model；Track structure uses the Symplectic Algebra Method Modeling 。K_dDynamic Stiffness Matrix for stock rail coupled system； For coupled systemes under dummy excitation effect The virtual displacement response vector of system, r_jFor corresponding loading coefficient vector.

According to dummy excitation law theory, solve the displacement of coupled system by following mathematic(al) representation and add Speed-power is composed:

$s_{\mathrm{uu}}\left(\mathrm{\&omega;}\right)=\underset{j=v,a,c}{\mathrm{\&Sigma;}}{{\tilde{u}}_j}^{*}(\mathrm{\&omega;},t){{\tilde{u}}_j}^T(\mathrm{\&omega;},t);s_{\stackrel{..}{u}\stackrel{..}{u}}\left(\mathrm{\&omega;}\right)={\mathrm{\&omega;}}^4{s}_{\mathrm{uu}}\left(\mathrm{\&omega;}\right)---\left(17\right)$

The most therefrom choose the corresponding free degree of flexible car body, i.e. obtain the random response of flexible car body.

The random response that step C, basis draw, utilizes the different ride comfort standard occupant comfort to car body Degree carries out calculating assessment.

Flexible car body FEM model is chosen the seat installation place on diverse location inside chassis of vehicle body Corresponding node is as research object, the random response power spectrum of flexible car body step B obtained Substitute into corresponding ride comfort standard both at home and abroad, calculate every flexibility index of car body, including:

China's GB5599-85 national standard:

$W=7.08{\left(\underset{0.5}{\overset{80}{\&Integral;}}{\left(\frac{F\left(f\right)}{f}\right)}^{\frac{3}{2}}G\left(f\right)\mathrm{df}\right)}^{\frac{3}{20}}---\left(18\right)$

Wherein, G (f) is the average power spectral density of acceleration samples；F is vibration frequency, and unit is H z；F (f) is frequency correction factor.Its evaluation criterion:

Stationarity grade	Evaluation	Riding index W
			1 grade	Excellent	<2.5
2 grades	Well	2.5-2.75
			3 grades	Qualified	2.75-3.0

TB/T2630-93:

$A_w=\sqrt{\underset{1}{\overset{80}{\&Integral;}}{B}^2\left(f\right)G\left(f\right)\mathrm{df}}---\left(19\right)$

Wherein, B (f) is frequency weighting function.Its evaluation criterion:

ISO2631-1:

$a_w=\sqrt{\underset{0.5}{\overset{80}{\&Integral;}}{W}^2\left(f\right)G\left(f\right)\mathrm{df}}---\left(20\right)$

Wherein, W (f) is frequency weighting function.Its evaluation criterion:

Root mean square of weighed acceleration aw	The subjective sensation of people
		＜ 0.315	The most uncomfortable
0.315～0.63	There are some uncomfortable
		0.55～1.0	The most uncomfortable
0.85～1.6	Uncomfortable
		1.255～2.5	Very uncomfortable
＞ 2.0	Like a hen on a hot girdle

UIC513:

$N_{\mathrm{MV}}=6\sqrt{{\left(a_{\mathrm{XP}95}^{w_d}\right)}^2+{\left(a_{\mathrm{YP}95}^{w_d}\right)}^2+{\left(a_{\mathrm{ZP}95}^{w_d}\right)}^2}---\left(21\right)$

Or

$N_{\mathrm{VA}}=4\left(a_{\mathrm{ZP}95}^{w_b}\right)+2\sqrt{{\left(a_{\mathrm{YA}95}^{w_b}\right)}^2+{\left(a_{\mathrm{ZA}95}^{w_b}\right)}^2}+4\left(a_{\mathrm{XD}95}^{w_e}\right)---\left(22\right)$

Wherein, , , , , With The weighted acceleration being respectively on the relevant position that this standard specifies is equal Root value.Its evaluation criterion:

Comfort level index	Human body sensory
		N ＜ 1	Best comfort
1<N<2	Good comfort
		2<N<4	Middle equivalent comfort
4<N<5	Bad comfortableness
		N>5	Extreme difference comfortableness

As a example by domestic CR H type bullet train, the application present invention provides the comfort of passenger specifying position Testing result.Car body limited element calculation model is as shown in Figure 4.

In step A, track irregularity takes China's circuit actual measurement spectrum, and speed is 200km/h, provides height , low direction and the dummy excitation of level.

In step B, ANSYS carries out the model analysis of FEM model, provide free vibration mode , and intercept its front 35 rank and set up coupled system equation, according to the calculating of formula (17), closed The acceleration power spectrum response of note position.Consider the symmetry of car body, choose on the left of car body floor The three corresponding positions of row seat are as research object, and each column selection takes 18 focuses.

In step C, the acceleration power spectrum obtained according to step B, respectively according to GB5599-85 standard, It is the most near the window that car body calculated by TB/T2630-93 standard, ISO-2631-1 standard and UIC513 standard meter Position, centre position, by the comfort level index of the focus of position, corridor.

As seen from Figure 5, the vertical and horizontal comfort level index using GB5599-85 criterion calculation is equal Reach " excellent " grade.As seen from Figure 6, hanging down of TB/T2630-93 criterion calculation is used " excellent " grade has all been reached to horizontal comfort level index.As seen from Figure 7, ISO is used The comfort level index that-2631-1 criterion calculation draws meets optimum, i.e. " does not has uncomfortable ".By Fig. 8 is it can be seen that the comfort level index using UIC513 criterion calculation to draw meets good grade, all Meet the requirement that new car dispatches from the factory.

As it has been described above, combine the plan content given by accompanying drawing, similar technical scheme can be derived .In every case be the content without departing from technical solution of the present invention, according to the present invention technical spirit to Any simple modification, equivalent variations and the modification that upper embodiment is made, all still falls within skill of the present invention In the range of art scheme.

Claims (3)

1. a bullet train comfort of passenger detection method, it is characterised in that described method includes:

A. pseudo-excitation method is used to generate the virtual rail irregularity for flexible car body random dynamic analysis；

Consider 4 wheels arbitrary excitation by same track, four wheels to time lag mathematic(al) representation as follows:

t₁=0；t₂=2l₁/v；t₃=2l₂/v；t₄=2 (l₁+l₂)/v

Wherein l₁, l₂It is respectively semifixed wheelbase and half bogie pivot center of vehicle；

Assume that three class track irregularities are zero-mean Stationary Gauss Random process, and there is no correlation each other, its power spectrum It is respectively S_v(ω)、S_a(ω) and S_c(ω), wherein, ω is time circular frequency；

According to dummy excitation law theory, generating corresponding virtual rail irregularity, its mathematic(al) representation is:

${\tilde{r}}_j(\mathrm{\&omega;},t)={\left\{\begin{array}{cccc}{e}^{-{\mathrm{i\&omega;t}}_1}& e^{-{\mathrm{i\&omega;t}}_2}& e^{-{\mathrm{i\&omega;t}}_2}& e^{-{\mathrm{i\&omega;t}}_4}\end{array}\right\}}^T\sqrt{S_j\left(\mathrm{\&omega;}\right)}{e}^{i\mathrm{\&omega;}t};(j=v,a,c)$

Wherein j=v, a, c represent high and low direction and level three class track irregularity respectively；

B. flexible car body random vibration frequency domain is analyzed, draws the random response of flexible car body；

Specifically include: use stock rail coupling system model, synthesize flexible car body mode motion equation, turn to Frame uses Rigid Body Dynamics Model to form the equation of motion, then is coupled with the bogie equation of motion by car body, sets up vehicle The equation of motion:

$(K_v+{\mathrm{i\&omega;C}}_v-{\mathrm{\&omega;}}^2{M}_v)u_v=K_v^d{u}_v=({\mathrm{i\&omega;Q}}_{vt}^C+Q_{vt}^K)u_w$

WhereinWithRepresent respectively vehicle and wheel under jth class track irregularity effect to virtual displacement；

Track structure is considered as three leafing scatterplot supporting constructions, and regards it as the cycle chain structure of infinite length, with stress Track minor structure is that research object sets up the equation of motion, synthesizes the Symplectic Algebra matrix, solve the feature of the Symplectic Algebra matrix to Amount, synthesizes stress minor structure steady motion equation；

The equation of motion of vehicle and four stress track minor structures is coupled together the most at last, and considers each sub-knot in above formula Interaction between structure, obtains the steady motion equation of coupled system:

$K_d{\tilde{u}}_j=r_j\sqrt{S_j}{e}^{i\mathrm{\&omega;}t}$

Wherein, flexible car body uses large-scale FEM model；Track structure uses the Symplectic Algebra Method Modeling；K_dFor car The Dynamic Stiffness Matrix of Track Coupling System；For under dummy excitation effect coupled system virtual displacement response to Amount, r_jFor corresponding loading coefficient vector；

According to dummy excitation law theory, solved displacement and the acceleration power of coupled system by following mathematic(al) representation Spectrum:

$S_{uu}\left(\mathrm{\&omega;}\right)=\underset{j=v,a,c}{\&Sigma;}{{\tilde{u}}_j}^{*}(\mathrm{\&omega;},t){{\tilde{u}}_j}^T(\mathrm{\&omega;},t);S_{\stackrel{\&CenterDot;\&CenterDot;}{u}\stackrel{\&CenterDot;\&CenterDot;}{u}}\left(\mathrm{\&omega;}\right)={\mathrm{\&omega;}}^4{S}_{uu}\left(\mathrm{\&omega;}\right)$

The most therefrom choose the corresponding free degree of flexible car body, i.e. obtain the random response of flexible car body；

C. the seat installation place chosen on flexible car body FEM model inside chassis of vehicle body on diverse location is corresponding Node is as research object, and the random response power spectrum of flexible car body step B obtained substitutes into corresponding domestic and international Ride comfort standard, calculates every flexibility index of car body.

Bullet train comfort of passenger detection method the most according to claim 1, it is characterised in that: in step C In both at home and abroad ride comfort standard refer specifically to be utilized respectively China's GB5599-85 standard, railway systems TB/T2630-93 mark Standard, the ISO2631-1 standard of International Organization for standardization and the UIC513 standard of International Union of Railways, to flexible car body not The comfort level index of co-located passenger calculates.

Bullet train comfort of passenger detection method the most according to claim 2, it is characterised in that: calculate and take advantage of The position of visitor's comfort level index is chosen position near the window, centre position and leans on position, corridor.