MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2015/630135 630135 Research Article Nonlinear Dynamic Analysis of Railway Wagon End Wall under the Action of Granular Media Wang Jing 1 Zhang Qichang 1 Yu Yuebin 2 Han Jianxin 1 Li Dong 3 Clayton John D. 1 School of Mechanical Engineering Tianjin University Tianjin 300072 China tju.edu.cn 2 Qiqihar Railway Rolling Stock Co., Ltd. Qiqihar 161992 China qrrs.chinacnr.com 3 Tianjin Key Laboratory of High Speed Cutting and Precision Machining Tianjin University of Technology and Education Tianjin 300222 China tute.edu.cn 2015 21102015 2015 12 05 2015 27 07 2015 21102015 2015 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.

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 wall-granular 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 stress-strain relationship, development of strength and deformation, and response of nonlinear dynamic and complicated coupling property of granular media.

With the rapid development of high-speed and heavy-loading 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  and related area. In Guaita et al.’s works , an elastoplastic constitutive law was applied to cylindrical silos with eccentric hoppers and rigid walls, using the Drucker-Prager criterion. Meng et al.  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.  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.  used FEM models to represent the various problems in silo analysis and the influence of different parameters was discussed. Zhang et al.  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.  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 soil-wall systems were carried out. To represent the unbounded half-space soil in a soil-structure interaction analysis, Wolf and Somaini  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  considered the vertical cross-sectional 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  discussed the effect of fundamental frequency of reinforced-soil 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 amplitude-frequency 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 high-speed and heavy-loading 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 . The wall was regarded as a pile with a diameter of a. The added mass of the whole pile is(1)Ma0=ρaπ27r2ba+16r3-ba9A=ρaπ27r2ba+16r3-baπa236,where r=1/2(2/3Ltanφ/2+d) is the diffusion radius, φ is the friction angle, and ba 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 is(2)M=M0+Ma=hρ+2ba-yhba2ρaπ127r2ba+16r3-a2ba36,where ρ is the density of the wall, ρa is the density of granular media, and M0 is the mass of the wall per unit area. Ma 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. a is the width of the wall, b is the height of the wall, h is the thickness of the wall, and ba 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:(3)D4wx4+24wx2y2+4wy4+cwt+M2wt2+v=px,y,t,where v=-{Eh/2(1-μ2)2w/x2[(w/x)2+μ(w/y)2] + Eh/2(1-μ2)2w/y2[(w/y)2+μ(w/x)2]+Eh/(1+μ)w/xw/y2w/xy}, in which w is the deflection of the wall, E is the elasticity modulus of the wall, μ is Poisson’s ratio of the wall, ρ is the density of the wall, D=Eh3/12(1-μ2) is the bending rigidity of the wall, and c is the damping coefficient.

Model of the rectangular plate with three clamped edges and a free edge.

The pressure  of granular media to the wall can be described by(4)px,y,t=p0+p1sinωt,where p0 is the static lateral pressure of the wall, p1 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 is(5)p0=K-γbαK-2b-ybαK-1+γb-yαK-2,where K=cosθcosδsin(θ-φ)/sinθcos(θ-φ+δ), α=cos(θ-φ-δ)tanθ/sin(θ-φ)cosδ, K 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 first-order modal of out-of-plane vibration, the vibration model functions can be written as(6)wx,y,t=Wx,yφt=XxYyφt.

The plate is with one free end on the top and three other fixed edges, so the boundary condition is as follows:

free end:(7)Myy=0=Vyy=0=0;

that is,(8)2wy2+μ2wx2y=0=0,3wy3+2-μ3wx2yy=0=0;

fixed edges:(9)wx=0=wx=a=wy=b=0;

that is,(10)wxx=0=wxx=a=wyy=b=0.

Given the vibration model function  as W(x,y)=(cos2πx/a-1)(1-cosπy/2b), then (6) can be described as(11)wx,y,t=cos2πxa-11-cosπy2b·φt.Substituting (11) into (3), the residual can be obtained by Galerkin method(12)R=4Wx,yφt+1DVx,yφ3t+cDWx,yφt+MDWx,yφt-p0+p1sinΩtDbecause(13)0a0bRWijx,ydxdy=0.Substituting (12) into (13), it can be obtained that(14)φt+a2a1φt+a3a1φt+a4a1φ3t=a5a1+a6a1sinΩt,where(15)a1=MD0a0bW2x,ydxdy,a2=cD0a0bW2x,ydxdy,a3=0a0bWx,y4Wx,ydxdy,a4=1D0a0bVx,yWx,ydxdy,a5=1D0a0baWx,yp0dxdy,a6=p1D0a0bWx,ydxdy.

For convenience, introduce the following nondimensional variables:(16)a3a1=ω02,Ψ=φh,τ=ω0t.

Consider ψ=u+φ0; then the equation of vibration becomes(17)u+u=εc1sinΩτω0+c2u+c3u3+3c3φ0u2+3c3φ02u,where(18)c1=a6a11hω02p1,c2=-a2a11ω0,c3=-a4a1h21ω02.

2.3. Perturbation Solution

The method of multiple scale is used to investigate the dynamic behavior of (17) under the primary resonance condition; namely, Ω/ω01.

To describe the nearness of the primary resonance condition, a detuning parameter σ is introduced and defined by(19)Ωω0=1+εσ.Introduce time-scale Tn=εnt, n=0,1,2,; then the approximate solution of (17) can be written as(20)ut=u0T0,T1+εu1T0,T1.Substituting (19) and (20) into (17) and equating the coefficients of ε0 and ε1 on both sides, we can obtain(21)D02u0+u0=0,(22)D02u1+u1=-2D0D1u0+c1sinT0+εσT0+c2D0u0+c3u03+3c3φ0u02+3c3φ02u0,where Dn=/Tn, n=0,1.

The solution of (21) is given in the form(23)u0=AT1eiT0+A-T1e-iT0;then (22) becomes(24)D02u1+u1=eiT0-2iD1A+c12ieiσT1+c2iA+3c3A2A-+3c3φ02A+c3A3e3iT0+3c3φ0A2e2iT0+3c3φ0AA-+cc,where cc denotes the complex conjugate of the preceding terms.

To eliminate secular terms, we must put(25)-2iD1A+c12ieiσT1+c2iA+3c3A2A-+3c3φ02A=0;then(26)2u1T02+u1=c3A3e3iT0+3c3φ0A2e2iT0+3c3φ0AA-+cc.So we obtain(27)u1=-18c3A3e3iT0-c3φ0A2e2iT0+3c3φ0AA-+cc.

It is convenient to write A in the polar form(28)AT1=12a0T1eiθT1.

Substituting (28) into (25), we obtain(29)da0dT1=12c2a0-c12cosγ,a0dγdT1=-a0σ-3c3a03-32c3φ02a0+c12sinγ,where γ=θ-σT1.

The amplitude-frequency response equation can be obtained as follows by imposing the condition da0/dT1=dγ/dT1=0:(30)2a0σ+6c3a03+3c3φ02a02+c2a02=c12,and the first approximated solution of system can be expressed as(31)wx,y,t=-164c3a03T1e3iθT1+3iT0-14c3φ0a02T1e2iθT1+2iT0+34c3φ0a02T1+cc·hcos2πxa-11-cosπy2b.

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.

Relationship between resonance frequency and height of granular media.

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.

Parameter table of granular media.

Particle size/(cm) ρ a /(kg/m3) φ /(°)
0.5 1900 32
0.6 1840 33
0.7 1780 34
0.8 1720 35
0.9 1660 36
1 1600 37
1.1 1540 38
1.2 1480 39
1.3 1420 40
1.4 1360 41
1.5 1300 42

Relationship between resonance frequency and particle size of granular media.

3.2. Amplitude-Frequency Response

Recall the amplitude-frequency 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 Amplitude-Frequency Response Curve

Height ratio of granular media is defined as ba/b. It is used to express the loading capacity. Figure 4 shows the amplitude-frequency response curve for ba=0.1,0.2,0.25; that is, ba/b=2/5,4/5,1, 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.

Amplitude-frequency response curve of system under different height ratio of granular media.

3.2.2. Effect on Amplitude-Frequency Response Curve by Height of the End Wall

Figure 5 shows the amplitude-frequency response curves when the height of the end wall b=0.23, 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 height-width ratio of the end wall has to be considered during the design of the end wall. Reasonable height-width 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 height-width ratio is known. So, the key point of designing the carriage is to determine the optimal height-width ratio of the carriage.

Amplitude-frequency response curve of system under different height of the wall.

3.2.3. Effect of Particle Size of Granular Media on Amplitude-Frequency Response

We further discuss the influence of gravel particle size loaded in the carriage on the amplitude-frequency 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.

Amplitude-frequency response curve of system under different particle size of granular media.

3.2.4. Effect of Thickness of the End Wall on Amplitude-Frequency Response Curve

The thickness of the end wall is a key parameter in design. Figure 7 shows amplitude-frequency 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 amplitude-frequency 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.

Amplitude-frequency response curve of system under different thickness of the wall.

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 δ=ba/b, by introducing the theory of finite element modeling of the retaining wall structure; we considered the special contact relationship (e.g., see ) between granular media and continuum and set up finite element models with different δ; the FEM model is shown in Figure 8.

FEM model of the wall system with different δ.

δ = 1

δ = 0.8

δ = 0.4

δ = 0

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. Drucker-Prager yield criterion is a correction form of the Mohr-Coulomb yield criterion. It fits for granular media. So this paper adopts Drucker-Prager 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 δ=0.8, for example, Figures 9 and 10 show the first-order mode shape solved by finite element method and theoretical method.

Material properties.

Type Density (kg/m3) Elastic modulus (Pa) Poisson ratio Cohesion (Pa) Frictional angle (°) Flow angle (°)
Wall 7850 2 × 1 0 11 0.3
Granular media 1636.74 2 × 1 0 7 0.3 0 30 15

FEM calculation results and theoretical calculation results.

Height ratio Theoretical calculation (Hz) FEM calculation result (Hz) Relative error
0 127.806 135.444 5.64%
0.4 127.782 134.786 5.20%
0.8 126.805 133.993 5.36%
1 125.767 133.499 5.79%

The first-order mode by FEM method.

The first-order mode by 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 heavy-loading 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:

In theoretical methods, this paper sets up nonlinear vibration equations which consider the affection of granular media and got the first-order approximate solution and amplitude-frequency response equations in the case of taking account of first primary resonance.

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.

Nonlinear vibration of the structure can be minimized by setting an appropriate height-width ratio. So the optimum height-width ratio of the end wall is the key factor during the design stage of the wagon.

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.

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.

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.

Argyroudis S. Kaynia A. M. Pitilakis K. Development of fragility functions for geotechnical constructions: application to cantilever retaining walls Soil Dynamics and Earthquake Engineering 2013 50 106 116 10.1016/j.soildyn.2013.02.014 2-s2.0-84876305049 Guaita M. Couto A. Ayuga F. Numerical simulation of wall pressure during discharge of granular material from cylindrical silos with eccentric hoppers Biosystems Engineering 2003 85 1 101 109 10.1016/S1537-5110(03)00037-0 2-s2.0-0037691526 Meng Q. Jofriet J. C. Negi S. C. Finite element analysis of bulk solids flow: part 2, application to a parametric study Journal of Agricultural Engineering Research 1997 67 2 151 159 10.1006/jaer.1997.0162 2-s2.0-0005068318 Goodey R. J. Brown C. J. Rotter J. M. Verification of a 3-dimensional model for filling pressures in square thin-walled silos Engineering Structures 2003 25 14 1773 1783 10.1016/j.engstruct.2003.07.003 2-s2.0-0242658759 Ayuga F. Guaita M. Aguado P. Static and dynamic silo loads using finite element models Journal of Agricultural Engineering Research 2001 78 3 299 308 10.1006/jaer.2000.0640 2-s2.0-0012627816 Zhang C. S. Yang Q. W. Xu Y. DEM simulation to the behaviors of bulk granule within a truck body during a sudden breaking The Chinese Journal of Process Engineering 2002 2 supplement 79 83 He H. Tian H. Q. Jiang F. Granule movement and influents to gondola car safe Proceedings of the International Conference on Mechanical Engineering and Mechanics (ICMEM '07) November 2007 Wuxi, China 1917 1921 Wolf J. P. Somaini D. R. Approximate dynamic model of embedded foundation in time domain Earthquake Engineering & Structural Dynamics 1986 14 5 683 703 10.1002/eqe.4290140502 Ghanbari A. Hoomaan E. Mojallal M. An analytycal method for calculating the natural frequency of retaining walls International Journal of Civil Engineering 2013 11 1 1 9 2-s2.0-84878336422 Hatami K. Bathurst R. J. Effect of structural design on fundamental frequency of reinforced-soil retaining walls Soil Dynamics and Earthquake Engineering 2000 19 3 137 157 10.1016/s0267-7261(00)00010-5 2-s2.0-0034024043 Zhang Y. A. Liu S. A. Chen J. G. Influence and its application of soil on the mode of cantilever retaining wall Journal of Civil, Architectural & Environmental Engineering 2009 31 3 72 77 2-s2.0-70350482963 Liu L. B. Zhang Y. X. Chen J. G. Dynamic-detection model of soil-wall system and the identification method of virtual parameter Journal of Northwest A & F University (Natural Science Edition) 2012 40 10 223 227 Cao G. X. Vibration of Thin Rectangular Elastic Plate 1983 Beijing, China China Building Industry Press Liu F. Q. Lateral earth pressures acting on circular retaining walls International Journal of Geomechanics 2014 14 3 10.1061/(asce)gm.1943-5622.0000291 Tanchaisawat T. Bergado D. T. Voottipruex P. Numerical simulation and sensitivity analyses of full-scale test embankment with reinforced lightweight geomaterials on soft Bangkok clay Geotextiles and Geomembranes 2008 26 6 498 511 10.1016/j.geotexmem.2008.05.005 2-s2.0-55749103896 Kitsabunnarat A. Alsaleh M. Helwany S. Capturing strain localization in reinforced soils Acta Geotechnica 2008 3 3 175 190 10.1007/s11440-008-0078-8 2-s2.0-51749103727 Hsia S.-Y. Chou Y.-T. Traffic noise propagating from vibration of railway wagon Mathematical Problems in Engineering 2013 2013 7 651797 10.1155/2013/651797 2-s2.0-84893834980 Bao X.-H. Ye G.-L. Ye B. Sago Y. Zhang F. Seismic performance of SSPQ retaining wall-centrifuge model tests and numerical evaluation Soil Dynamics and Earthquake Engineering 2014 61-62 63 82 10.1016/j.soildyn.2014.01.019 2-s2.0-84894453075 Zhang F. Ye B. Ye G. L. Unified description of sand behavior Frontiers of Architecture and Civil Engineering in China 2011 5 2 121 150 10.1007/s11709-011-0104-z 2-s2.0-79958036420 Athanasopoulos -Zekkos A. Vlachakis V. S. Athanasopoulos G. A. Phasing issues in the seismic response of yielding, gravity-type earth retaining walls—overview and results from a FEM study Soil Dynamics and Earthquake Engineering 2013 55 59 70 10.1016/j.soildyn.2013.08.004 2-s2.0-84885042427 Mao J. Q. A finite element approach to solve contact problems in geotechnical engineering International Journal for Numerical and Analytical Methods in Geomechanics 2005 29 5 525 550 10.1002/nag.424 2-s2.0-17744376859 Liang X. W. Liu Q. Finite element analysis of the cantilever retain wall Journal of Huazhong University of Science and Technology (Natural Science Edition) 2003 20 2 99 102