Research Article  Open Access
Nonlinear Dynamic Analysis of Railway Wagon End Wall under the Action of Granular Media
Abstract
In wagon end wall system under the action of granular media, the model of the wall is the result of multiple factors. To figure out how the key factors act, the large amplitude vibration together with the contribution from the action of granular media to build the nonlinear dynamics equation of the wagon end wall has been taken into account. In this paper, the effects of varying damping, added mass, and force due to the action of granular media on the patterns of movement and response performance were discussed in detail. Furthermore, the finite element model of end wallgranular media was established. It proves that the theoretical result is correct compared with the FEM results. It indicates that the resonance frequency of the wall declines with the increase of the height of granular media, while the frequency increases with the increase of the particle size of granular media. The nonlinear dynamics of the wall is less affected by the height ratio and particle size of granular media, while it is greatly affected by the height and the thickness of wall, which should be seriously considered in structural design of the wagon.
1. Introduction
70% of the total railway freight transportation is with granular media loaded in, such as coal and rubble. Research on dynamics of the wagon is acknowledged difficult for the special stressstrain relationship, development of strength and deformation, and response of nonlinear dynamic and complicated coupling property of granular media.
With the rapid development of highspeed and heavyloading trains, the security of the railroad running becomes more and more important. Thus, it is essential to design some new types of heavy haul freight cars. However, as the car speed increases, the intensity of excitation and the frequency variation increase. Then, wagon will show a complex dynamic response under the action of granular media, so the deformation of bending and torsion will increase and the probability of local damage will also increase. The mechanical properties increase the difficulty of the design and limit the development of cars. So, it is significantly meaningful to investigate the nonlinear dynamics of the bulk carrier.
A group of international and domestic academics launched a series of research on mechanics of granular media [1] and related area. In Guaita et al.’s works [2], an elastoplastic constitutive law was applied to cylindrical silos with eccentric hoppers and rigid walls, using the DruckerPrager criterion. Meng et al. [3] dealt with the application of a secant constitutive relationship for bulk solids in a transient finite element analysis of stresses, displacements, and velocities during discharge of soybean from a silo. Goodey et al. [4] outlined the development and verification of a finite element model for the filling pressure distribution in a square platform silo with flexible walls. Ayuga et al. [5] used FEM models to represent the various problems in silo analysis and the influence of different parameters was discussed. Zhang et al. [6] studied the mechanical behavior of granular media in a wagon by discrete element method (DEM). The mechanical behavior of granular media was studied in different conditions by He et al. [7] and they drew the conclusion that lateral static pressure is under linear distribution along height division of wagon.
However, the above works lacked studies on dynamics of the bulk carrier, and their theories do not perform very well when applied to wagon end wall system, since they fail to consider the coupling relationship between granular media and bulk carrier. Several researches on dynamic behavior of soilwall systems were carried out. To represent the unbounded halfspace soil in a soilstructure interaction analysis, Wolf and Somaini [8] used a discrete model that consists of a mass which is attached to a rigid support with a spring and a damper. Ghanbari et al.’s formula [9] considered the vertical crosssectional width change, thus calculating the natural frequency of retaining walls. As the increase of width to length ratio, the frequency of wall decreases. Hatami and Bathurst [10] discussed the effect of fundamental frequency of reinforcedsoil retaining walls on the influence of some design parameters including wall height, backfill width, reinforcement stiffness, reinforcement length, backfill friction angle, and toe restraint condition. But there are few researches focusing on wagon end wall. The coupling relationship could be measured by resonant frequency and amplitudefrequency response of system, and these parameters will significantly influence the performance of the wagon end wall. To fill this research gap, this paper proposes using theoretical and FEM method based on geotechnical mechanics to study the relationship between the key structural designing parameter and the dynamic response of the wagon end wall after highspeed and heavyloading trains are loaded with granular media. The numerical study results clearly show the influences of the key structural designing parameters on the dynamic response of the wagon end wall and demonstrate our conclusions. This paper provides a foundation for the design of the wagon and the structural optimization and improves the effect of working performance.
2. Theoretical Analysis
2.1. The Derivation Process of Added Mass about Wagon End Wall
We first need to consider the added mass caused by granular media when computing the dynamic characteristics of wall contact with granular media. Because the nature of granular media lies between solid and fluid, it has complicated mechanics characteristics, though the particle is of simple structure.
A semiempirical formula was introduced according to the results of experiment and the empirical derivation [11]. The wall was regarded as a pile with a diameter of . The added mass of the whole pile iswhere is the diffusion radius, is the friction angle, and is the height of granular media.
Distribution of added mass was simplified into triangular distribution. The added mass on the top was 0, so the total added mass per unit area iswhere is the density of the wall, is the density of granular media, and is the mass of the wall per unit area. is the added mass per unit area.
2.2. The Establishment of Dynamics Equation
The dynamical model of the wagon end wall was simplified as a plate with one free edge on the top and three other fixed edges. This model has a good performance on the dynamic characteristic of the system and can describe it accurately, as shown in Figure 1. is the width of the wall, is the height of the wall, is the thickness of the wall, and is the height of granular media. The nonlinear dynamic equation of a wall is derived based on Von Karman theory, with the added mass and damping effect taken into account:where + , in which is the deflection of the wall, is the elasticity modulus of the wall, is Poisson’s ratio of the wall, is the density of the wall, is the bending rigidity of the wall, and is the damping coefficient.
The pressure [12] of granular media to the wall can be described bywhere is the static lateral pressure of the wall, is the pressure amplitude of vibration provided by granular media, and is the vibration frequency of granular media.
According to Coulomb’s earth pressure theory, the lateral pressure of the wall iswhere , , is the lateral pressure coefficient of granular media, is the angle between slip surface and horizontal surface, is the internal friction angle of granular media, is the angle of friction between the wall and granular media, and is the force density of granular media.
Based on firstorder modal of outofplane vibration, the vibration model functions can be written as
The plate is with one free end on the top and three other fixed edges, so the boundary condition is as follows: free end: that is, fixed edges: that is,
Given the vibration model function [13] as , then (6) can be described asSubstituting (11) into (3), the residual can be obtained by Galerkin methodbecauseSubstituting (12) into (13), it can be obtained thatwhere
For convenience, introduce the following nondimensional variables:
Consider ; then the equation of vibration becomeswhere
2.3. Perturbation Solution
The method of multiple scale is used to investigate the dynamic behavior of (17) under the primary resonance condition; namely, .
To describe the nearness of the primary resonance condition, a detuning parameter is introduced and defined byIntroduce timescale , ; then the approximate solution of (17) can be written asSubstituting (19) and (20) into (17) and equating the coefficients of and on both sides, we can obtainwhere , .
The solution of (21) is given in the formthen (22) becomeswhere denotes the complex conjugate of the preceding terms.
To eliminate secular terms, we must putthenSo we obtain
It is convenient to write in the polar form
Substituting (28) into (25), we obtainwhere .
The amplitudefrequency response equation can be obtained as follows by imposing the condition :and the first approximated solution of system can be expressed as
3. Theoretical Analysis
3.1. Resonance Frequency
Resonance frequency is an important design index. The effect of height and particle size of granular media on resonance frequency is discussed in this section.
3.1.1. Effect of Height of Granular Media on Resonance Frequency
The height of granular media directly influences the system vibration. Figure 2 shows that the resonance frequency decreases with the increase of the height of granular media. The added mass increases with the increment of height, which induces the reduction of frequency. At the same time, the lateral pressure on the end wall increases as the granular media increase. It leads to the increment of stiffness and then increment of system frequency. Both factors have an influence together and consequently lead to the reduction of frequency. Furthermore, it should be noted that the added mass plays a greater influence on resonance frequency than the lateral pressure.
3.1.2. Effect of Particle Size of Granular Media on Resonance Frequency
The granular media transported by heavy railway freight car often show different particle sizes, and some typical types of cargoes include coal, ore, wood, steels, and other bulk cargoes. Granular media in different particle size have a direct impact on the vibration of the end wall system. Table 1 shows the density and friction angle of different particle size of gravel. As the particle size increases, density reduces and friction angle increases. And Figure 3 also illustrates that the system frequency increases with the increment of particle size.

3.2. AmplitudeFrequency Response
Recall the amplitudefrequency response of the end wall in (30). The nonlinear controlling factors include height ratio of granular media, height of the end wall, particle size of granular media, and thickness of the end wall. This section discusses the influence of the above parameters on the vibration characteristics of the system.
3.2.1. Effect of Height Ratio of Granular Media on AmplitudeFrequency Response Curve
Height ratio of granular media is defined as . It is used to express the loading capacity. Figure 4 shows the amplitudefrequency response curve for ; that is, , with red, blue, and green curve, respectively. The resonance frequency of the end wall decreases as the height ratio of granular media increases. However, the changing amplitude of the curves is negligibly small, and the stable amplitudes together with the resonance floating excursion are nearly the same. It suggests that the stiffness of the end wall contributed by the interaction between the end wall and granular media on working condition is negligible, and the amount of the granular media has little effect on the nonlinear characteristics of the end wall. Therefore, it is not necessary to consider whether the carriage is in working condition at design time. It will reduce the number of the working conditions to be checked and then light the designing workload.
3.2.2. Effect on AmplitudeFrequency Response Curve by Height of the End Wall
Figure 5 shows the amplitudefrequency response curves when the height of the end wall , 0.25, 0.3 with wagon fully loaded. With the increment of the height of end wall, the nonlinear resonance frequency and the resonance floating excursion of the end wall increase and the nonlinearity becomes stronger. Furthermore, the nonlinearity of structural vibration has the same trend. This indicates that the higher the end wall is, the greater the impact of the granular media on the nonlinear characteristics of the end wall is. So the heightwidth ratio of the end wall has to be considered during the design of the end wall. Reasonable heightwidth ratio can reduce the nonlinear vibration of the structure. The height of the carriage should be within a fixed range because the track width is a fixed value in a country. The size of the end wall can be definite if the heightwidth ratio is known. So, the key point of designing the carriage is to determine the optimal heightwidth ratio of the carriage.
3.2.3. Effect of Particle Size of Granular Media on AmplitudeFrequency Response
We further discuss the influence of gravel particle size loaded in the carriage on the amplitudefrequency response. Three particle sizes, 0.5 cm, 1 cm, and 1.5 cm, are considered. As shown in Figure 6, the nonlinear resonance frequency of the end wall increases with the increment of the particle size of the granular media. However, the increasing sizes of the curves are quite small and the resonance floating excursion is nearly the same. Therefore we can conclude that the particle size of the granular media has little impact on nonlinear characteristics of structure. The end wall of the carriage has good versatility. The carriage can carry different kinds of cargoes in working condition. So, it does not need to design different carriages for different cargoes. This point can reduce the design workload greatly and has significant implications for the design work.
3.2.4. Effect of Thickness of the End Wall on AmplitudeFrequency Response Curve
The thickness of the end wall is a key parameter in design. Figure 7 shows amplitudefrequency response curve of the system under different thickness, 1.92 mm, 2 mm, and 3 mm, of the wall. It can be seen clearly that the resonance floating excursion decreases greatly as the thickness increases. The resonance floating excursion decreases to nearly 0 when the thickness goes larger than 3 mm. Then, in this case, the end wall shows linear vibration characteristics and steady amplitudefrequency response. It proves that the thickness of end wall can affect the vibration characteristics of the structure. But when the thickness is larger than a specific value, stiffness should be regarded as the key factor which can affect the nonlinear characteristics of the structure instead of thickness. It will reduce the design workload to some extent.
4. Numerical Analysis
This section sets up finite element model to calculate the system for studying the vibration characteristics of the end wall. It validates the rationality of the theory and method above and testifies the correctness of the calculation result.
The key issues to consider are the nonlinear constitutive relation of the granular media, the simulation of contact relationship between granular media and end wall, and the simulation of the boundary condition to set up the finite element model of the end wall. We define the height ratio of the granular media , by introducing the theory of finite element modeling of the retaining wall structure; we considered the special contact relationship (e.g., see [14–22]) between granular media and continuum and set up finite element models with different ; the FEM model is shown in Figure 8.
(a)
(b)
(c)
(d)
Deformation of granular media relates to current stress state and loading process because granular media are anisotropic material and do not follow Hook’s law. So setting up increment model can reflect the loading process better. DruckerPrager yield criterion is a correction form of the MohrCoulomb yield criterion. It fits for granular media. So this paper adopts DruckerPrager elastoplastic model in granular media simulation and finds the contact relationship between granular media and the end wall from simulated contact elements. The length of the granular media area in the model should be much more than the height of the end wall to diminish the influence of the boundary condition. The granular media area is with one free edge on the top and three other fixed edges. The boundary condition of the end wall is with one free edge on the top and three other fixed edges. Table 2 shows the material properties of the end wall and granular media. The FEM model of the end wall system was calculated by Block Lanczos law. Table 3 compares the FEM calculation result and theoretical calculation results by solving the FEM model of the end wall system according to Block Lanczos law. Taking , for example, Figures 9 and 10 show the firstorder mode shape solved by finite element method and theoretical method.


The result shows that the nonlinear resonance frequency decreases as the height of the granular media increases. But the decreasing magnitude is not significant, and the resonant frequencies and vibration mode shapes solved by theoretical method and finite element method are almost the same. However, the theoretical solution is a little lower than that of FEM result. The reason is that the frequency is influenced by stiffness and quality of the system. To be more specific, the length of the granular media area is unlimited in the theoretical method, but it is limited in the finite element method. Therefore, the increment of quality makes the result lower. But the error is within acceptable range. We will discuss the influence of the finite soil volume in the theoretical model in future work. At the same time, this model provides a basis for simulating the retaining wall structure by FEM.
5. Conclusion
This paper researched the dynamics of the end wall of heavyloading trains by theoretical methods and verified this theory by numerical simulations method. Finally, we can reach the conclusion that some key parameters are related to the design of end wall. The conclusions are as follows:(1)In theoretical methods, this paper sets up nonlinear vibration equations which consider the affection of granular media and got the firstorder approximate solution and amplitudefrequency response equations in the case of taking account of first primary resonance.(2)It certified that the amount of materials has little effect on the nonlinear characteristics of the end wall. So whether the carriage is in working condition can be neglected.(3)Nonlinear vibration of the structure can be minimized by setting an appropriate heightwidth ratio. So the optimum heightwidth ratio of the end wall is the key factor during the design stage of the wagon.(4)The particle size of the granular media has little effect on the nonlinear characteristics of the structure. It is not necessary to design different carriages for different kinds of cargo.(5)The thickness of the end wall affects the dynamics of the structure in a large part. However, it can just consider the influence of the strength instead of thickness when the thickness is greater than the particular value when studying the nonlinear characteristics.(6)The finite element model of the end wall was established which considers the influence of the granular media. By calculating the resonant frequency under different height ratio of granular media, the result validates the correctness of the theory in the paper.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work was financially supported by the Natural Science Foundation of Tianjin, China, Project no. 12JCZDJC28000.
References
 S. Argyroudis, A. M. Kaynia, and K. Pitilakis, “Development of fragility functions for geotechnical constructions: application to cantilever retaining walls,” Soil Dynamics and Earthquake Engineering, vol. 50, pp. 106–116, 2013. View at: Publisher Site  Google Scholar
 M. Guaita, A. Couto, and F. Ayuga, “Numerical simulation of wall pressure during discharge of granular material from cylindrical silos with eccentric hoppers,” Biosystems Engineering, vol. 85, no. 1, pp. 101–109, 2003. View at: Publisher Site  Google Scholar
 Q. Meng, J. C. Jofriet, and S. C. Negi, “Finite element analysis of bulk solids flow: part 2, application to a parametric study,” Journal of Agricultural Engineering Research, vol. 67, no. 2, pp. 151–159, 1997. View at: Publisher Site  Google Scholar
 R. J. Goodey, C. J. Brown, and J. M. Rotter, “Verification of a 3dimensional model for filling pressures in square thinwalled silos,” Engineering Structures, vol. 25, no. 14, pp. 1773–1783, 2003. View at: Publisher Site  Google Scholar
 F. Ayuga, M. Guaita, and P. Aguado, “Static and dynamic silo loads using finite element models,” Journal of Agricultural Engineering Research, vol. 78, no. 3, pp. 299–308, 2001. View at: Publisher Site  Google Scholar
 C. S. Zhang, Q. W. Yang, and Y. Xu, “DEM simulation to the behaviors of bulk granule within a truck body during a sudden breaking,” The Chinese Journal of Process Engineering, vol. 2, supplement, pp. 79–83, 2002. View at: Google Scholar
 H. He, H. Q. Tian, and F. Jiang, “Granule movement and influents to gondola car safe,” in Proceedings of the International Conference on Mechanical Engineering and Mechanics (ICMEM '07), pp. 1917–1921, Wuxi, China, November 2007. View at: Google Scholar
 J. P. Wolf and D. R. Somaini, “Approximate dynamic model of embedded foundation in time domain,” Earthquake Engineering & Structural Dynamics, vol. 14, no. 5, pp. 683–703, 1986. View at: Publisher Site  Google Scholar
 A. Ghanbari, E. Hoomaan, and M. Mojallal, “An analytycal method for calculating the natural frequency of retaining walls,” International Journal of Civil Engineering, vol. 11, no. 1, pp. 1–9, 2013. View at: Google Scholar
 K. Hatami and R. J. Bathurst, “Effect of structural design on fundamental frequency of reinforcedsoil retaining walls,” Soil Dynamics and Earthquake Engineering, vol. 19, no. 3, pp. 137–157, 2000. View at: Publisher Site  Google Scholar
 Y. A. Zhang, S. A. Liu, and J. G. Chen, “Influence and its application of soil on the mode of cantilever retaining wall,” Journal of Civil, Architectural & Environmental Engineering, vol. 31, no. 3, pp. 72–77, 2009. View at: Google Scholar
 L. B. Liu, Y. X. Zhang, and J. G. Chen, “Dynamicdetection model of soilwall system and the identification method of virtual parameter,” Journal of Northwest A & F University (Natural Science Edition), vol. 40, no. 10, pp. 223–227, 2012. View at: Google Scholar
 G. X. Cao, Vibration of Thin Rectangular Elastic Plate, China Building Industry Press, Beijing, China, 1983.
 F. Q. Liu, “Lateral earth pressures acting on circular retaining walls,” International Journal of Geomechanics, vol. 14, no. 3, 2014. View at: Publisher Site  Google Scholar
 T. Tanchaisawat, D. T. Bergado, and P. Voottipruex, “Numerical simulation and sensitivity analyses of fullscale test embankment with reinforced lightweight geomaterials on soft Bangkok clay,” Geotextiles and Geomembranes, vol. 26, no. 6, pp. 498–511, 2008. View at: Publisher Site  Google Scholar
 A. Kitsabunnarat, M. Alsaleh, and S. Helwany, “Capturing strain localization in reinforced soils,” Acta Geotechnica, vol. 3, no. 3, pp. 175–190, 2008. View at: Publisher Site  Google Scholar
 S.Y. Hsia and Y.T. Chou, “Traffic noise propagating from vibration of railway wagon,” Mathematical Problems in Engineering, vol. 2013, Article ID 651797, 7 pages, 2013. View at: Publisher Site  Google Scholar
 X.H. Bao, G.L. Ye, B. Ye, Y. Sago, and F. Zhang, “Seismic performance of SSPQ retaining wallcentrifuge model tests and numerical evaluation,” Soil Dynamics and Earthquake Engineering, vol. 6162, pp. 63–82, 2014. View at: Publisher Site  Google Scholar
 F. Zhang, B. Ye, and G. L. Ye, “Unified description of sand behavior,” Frontiers of Architecture and Civil Engineering in China, vol. 5, no. 2, pp. 121–150, 2011. View at: Publisher Site  Google Scholar
 A. Athanasopoulos Zekkos, V. S. Vlachakis, and G. A. Athanasopoulos, “Phasing issues in the seismic response of yielding, gravitytype earth retaining walls—overview and results from a FEM study,” Soil Dynamics and Earthquake Engineering, vol. 55, pp. 59–70, 2013. View at: Publisher Site  Google Scholar
 J. Q. Mao, “A finite element approach to solve contact problems in geotechnical engineering,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 29, no. 5, pp. 525–550, 2005. View at: Publisher Site  Google Scholar
 X. W. Liang and Q. Liu, “Finite element analysis of the cantilever retain wall,” Journal of Huazhong University of Science and Technology (Natural Science Edition), vol. 20, no. 2, pp. 99–102, 2003. View at: Google Scholar
Copyright
Copyright © 2015 Jing Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.