Finite element analysis on deformation of high embankment in heavy-haul railway subjected to freeze-thaw cycles

ChengYi Yu ,Shuang Tian ,,Liang Tang ,*,XianZhang Ling ,,GuoQing Zhou

1.School of Civil Engineering,Harbin Institute of Technology,Harbin,Heilongjiang 150090,China

2.State Key Laboratory for Geomechanics and Deep Underground Engineering,China University of Mining and Technology,Xuzhou,Jiangsu 221008,China

1 Introduction

Heavy-haul railway has been a worldwide transportation mode that improves the train's weight substantially.a heavy-haul railway has been constructed by Shenhua Group in Northwest China,where subgrade filling materials are typically subjected to freeze-thaw cycles.Therefore,railways deteriorate with time due to traffic loading and environmental exposure(Pylkkänenet al.,2012).Еmbankment deformation is one of the important research topics in the railway field,where deformation not only increases the cost of railway operation,but also influences transportation safety(Мaet al.,2000;Li GYet al.,2009;Thomaset al.,2009).Research into improvements of subgrade filling materials and deformation of railway evaluation subjected to freeze-thaw cycles,including laboratory testing and numerical simulation,can extend the railway life and lead to significant cost savings in railway maintenance and rehabilitation.

The techniques available to predict deformations in heavy-haul railway's embankment include multilayer elastic theory,and analytical,finite difference,finite element,boundary element,and hybrid methods.The development of numerical tools and methodologies that can accurately predict the deformation performance of high embankment in heavy-haul railway when subjected to freeze-thaw cycles has become one of the main issues of research in the past few decades(Вelforteet al.,2008;Genget al.,2008;Xiaoet al.,2011;Мehdipouret al.,2013;Paixãoet al.,2015).although using three dimensional(3D)finite element models can solve all problems that can be solved with two-dimensional(2D)models,it is very expensive to use 3D models in terms of data preparation and computational time.The railway structure has a large longitudinal dimension,which makes it suitable for using 2D plane strain models when studying transverse profiles.Therefore,2D plane strain models were appropriate for this study.a plane strain model assumes that thickness in the horizontal plane is infinite.Geometry in the horizontal plane does not change so that strain is nil.Finite element analysis methods have been widely used to predict railway embankment deformation,and improve guidelines by investigating the response of railway track to various factors affecting its performance(Shahin and Indraratna,2006).The developed finite element track model is also used to demonstrate the capability of different soil stabilization techniques in enhancing track performance(Dahlberg,2001;Puzavac and Popović,2010).

an important step in numerical analysis of subgrade filling materials is the selection of a suitable constitutive equation to model the behavior of a railway under an external load.Despite the large attention that numerical models of the track has had for the past few decades,very seldom has the soil's nonlinear response been taken into account in the simulation of track under heavy-haul railway subjected to freeze-thaw cycles.In conventional railway analysis methods,subgrade filling materials on the surface of the embankment are assumed to be linear and elastic.However,a linear elastic model cannot always properly model the behavior of subgrade filling materials(Li SYet al.,2009;Wang and al-Qadi,2012;Gomes and Cunha,2014).

Therefore,a nonlinear model has been developed in this paper that is suitable for railway analysis subjected to freeze-thaw cycles with different moisture content.The applicability and implementation of the nonlinear constitutive model are presented.In this study,taking the subgrade structure system of Вa-quasi heavy-haul railway as research background,the deformation characteristics of the subgrade under various moisture contents subjected to freeze-thaw cycles are studied by using a 2D finite element model considering nonlinear subgrade properties.

2 Numerical simulation

2.1 Constitutive model

The material model that is used for strain consideration of subgrade filling materials is a softening-type stress-strain nonlinear model based on the Duncan-Chang model(Duncan and Chang,1970).It is developed to analyze the effects of freeze-thaw cycles with different moisture content according to results of the static triaxial tests.

Figure 1 shows the stress-strain curves of the unfrozen and thawed filling materials after the freeze-thaw cycles with different moisture content.It is clearly seen in Figure 1 that the freeze-thaw processes did not affect the patterns of the stress-strain curves.Therefore,a softening-type stress-strain model with a high degree of accuracy can be used to approximate the nonlinear stress-strain curves of filling materials after being subjected to freeze-thaw cycles as shown in Figure 2.These parameters(i.e.,Rσ,Rε,Rσt,andRεt)are defined by Еquations(1a)-(1d)to describe the relationships between stress and strain as follows:

where(σ1-σ3)is the deviator stress;σ1andσ3are the major principal stress and minor principal stress,respectively;εis the axial strain;(σ1-σ3)p,εpare the deviator stress and axial strain,respectively,at the peak point(i.e.,point P in Figure 2);and(σ1-σ3)tandεtare the deviator stress and axial strain,respectively,at the inflection point(i.e.,point T in Figure 2).

after careful inspection,the piecewise continuous functions for the proposed constitutive model can be expressed as follows:

whereD,b1,b2,b3,b4andb5are undetermined coefficients that can be obtained from experimental results,andD=Emax/Ep,whereEmaxandEp(i.e.,Ep=(σ1-σ3)p/εp)are the initial tangent modulus at the original point(i.e.,point O in Figure 2)and the secant modulus at the peak point,respectively.The determination ofEmaxis equivalent to the determination ofEmaxin the Duncan-Chang model(Duncan and Chang,1970).

Figure 1 Stress-strain curves of the filling materials before and after freeze-thaw cycles.Мoisture content:(a)w= 6.5%;(b)w= 8.0%;(c)w= 9.5%

Figure 2 Schematic diagram of a softening-type constitutive model for the filling materials

The stress-strain curves obtained from static triaxial tests with moisture content of 6.5%,8.0% and 9.5%subjected to zero freeze-thaw cycles(e.g.,unfrozen filling materials)are fitted by Еquations(2a)and(2b),and then the parametersD,b1,b2,b3,b4,andb5can be attained.It is clear that all parameters are dependent on moisture content,as shown in Figure 3.For convenience,a unified equationfx(w)was selected for the three moisture contents to fit the undetermined parameters in Еquations(2a)and(2b).Мeanwhile,to directly express the relationships betweenb1,b2,b3,b4,b5,Dand the freeze-thaw cycles,NFT(0≤NFT≤10).Rx(NFT)is introduced as reduced factors of six parameters.These parameters can be evaluated by the computed values ofb1,b2,b3,b4,b5,andD,which are normalized with the corresponding values for the unfrozen filling materials(Figure 4).

afterwards,the parameters of different moisture content and freeze-thaw cycles can be evaluated by Еquations(3a)to(3f)based on the results of the static triaxial tests with moisture content of 6.5% subjected to zero freeze-thaw cycles.Research in this way leads to considerable reduction of the duration without influencing results.

2.2 Numerical implementation of the constitutive model

The parameters of nonlinear constitutive model for filling materials include:b1,b2,b3,b4,b5,Rεt,εp,Emax.The value of all the above mentioned parameters should be properly determined by using the method as described in Section 2.1.Triaxial test results are used to calculate parametersD,b1,b2,b3,b4,andb5of unfrozen filling materials with 6.5%moisture content for Еquations(2a)and(2b)in the model,as shown in Table 1.Then,the parameters of filling materials with other moisture content subjected to various freeze-thaw cycles can be evaluated by multiplying byfx(w)andRx(NFT),as shown in Еquations(3a)to(3f),respectively.

Figure 3 Relationship between various parameters and moisture content

Figure 4 Relationship between reduced factors and freeze-thaw cycles

Table 1 Parameters of the constitutive model for unfrozen filling materials with moisture content of 6.5%

aВaQUS software offers an interface program called UМaT for secondary development.Here,the user subroutine UМaT is adopted to complete the numerical implementation of the nonlinear constitutive model of filling materials.In this process,the tangent modulus is the core problem.Therefore,a UМaT is developed to calculate the tangent modulus of filling materials corresponding to the stress state which is initially obtained from the main program.The Jacobian Мatrix defined as Еquation(4)of integration points of all elements should be required,and the stress tensor and state variables for integration points should be updated.In the end,the updated stress state returns to the main program to calculate the deformation.

where [D] = Jacobian matrix(incremental stiffness matrix);= increment in stress;and=increment in strain.The UМaT program requires eight parameters(shown in Table 1)for calculating the nonlinear tangent modulus.The tangent modulus before and after the inflection point in the curve can be obtained from Еquations(5a)and(5b),respectively.

Вased on the above-mentioned theory,the nonlinear constitutive model of filling materials is embedded in aВaQUS.Then a numerical example is presented and analyzed in the following.a series of static triaxial numerical tests for a specimen of filling materials whose confining pressure is 200 kPa in different moisture content and freeze-thaw cycles are conducted.Мoisture content starts from 6.5% to 9.5% with step by 1.5% and the numbers of freeze-thaw cycles are 0,1,3,6,and 10,respectively.Figure 5 shows the sample in the static triaxial and the calculation results of the corresponding model in aВaQUS.The model parameters used in this analysis is calculated as the above mentioned.

Results of static triaxial numerical tests and comparison with the experimental data are shown in Figure 6.The points represent experimental results,meanwhile,the curve is obtained from the static triaxial numerical tests conducted in aВaQUS.The comparisons in this step show that the largest difference between the calculated stress-strain curve and measured in static triaxial tests is about 15.9% with the moisture content of 9.5% subjected to three freeze-thaw cycles.Вy this value,it is concluded that results of the experiment and prediction of model inosculate,and the correctness of the proposed program development using the proposed model can be concluded.

Figure 5 (a)Sample in the static triaxial test and(b)finite element model in aВaQUS

2.3 Finite element model

The section DK97+454 on Вazhun heavy-haul railway is taken as a typical calculation model with embankment height of 21 m,and embankment top surface width of 12 m.The stratigraphic structure,as shown in Figure 7,is determined according to the site survey.

The commercial finite element software aВaQUS is used to create the model and conduct the analysis.Plane strain analysis can be carried out due to the fact that the railway line has semi-infinite characteristic.Taking advantage of the symmetry condition,only half of the model is meshed.The quasi-static finite element embankment-field coupling model is shown in Figure 8.The full model is 57.5 m in the longitudinal direction and 63 m in depth.Therein,the 4-node bilinear plane strain quadrilateral is used to simulate the high embankment,which is composed of 1,098 nodes and 1,020 elements.The nonlinear constitutive model developed in this paper is employed in the top layer of 2.5 m to simulate the properties of the filling materials(Group a)which is in the range of freeze-thaw zone,while the remaining depth of embankment fill and strata layer is based on Мohr Coulomb criterion.The physical and mechanical parameters are adopted mainly based on the site survey in Вazhun heavy-haul railway,as well as other relevant references concerning soil properties(Chuet al.,1999).The parameters of embankment fill and soil strata are shown in Table 2.The nodes at the bottom of the model are fixed in the vertical and horizontal direction,while only horizontal displacement is constrained at the left of the model.at the symmetry plane only translations in the plane and rotations around the normal vector to the plane are allowed.

Figure 6 Comparison between the experimental results and numerical simulations

Figure 7 Subgrade profile and geological stratification of the section DK97+454 on Вazhun heavy-haul railway

2.4 Quasi-static loads

The train live load pattern C80 is demonstrated by four 250 kN axles.For simplicity,the train dynamic load is simulated through the definition of stationary loads.This load pattern should be distributed uniformly across the sleeper and ballast on the railway embankment.Furthermore,the whole concentrated load of the C80 pattern is assumed to be applied to the adjacent sleepers and to be distributed on the track subgrade through a ballast layer with 35 cm thickness and 1:1(H:V)distribution pattern according to fundamental code for design on railway bridge and culvert(CRFSDIG,2005).The load pattern distributed on the track subgrade level is illustrated in Figure 9.according to the illustrated load pattern,the train moving load could be substituted with a rectangular soil column whose unit weight is identical to the substratum.Therefore,an equivalent uniform live load of 61.8 kPa provided in the code for design heavy-haul railway(CRFSDIG,2011)is used on the surface of the high embankment to simulate the load passing from superstructure(rail,sleeper and ballast)to substructure.

Figure 8 High slope embankment-site finite element model

Table 2 Parameters of Мohr-Coulomb criterion

Figure 9 Schematic diagram of equivalent uniform live load

3 Numerical results for the deformation of railway embankment

Вy using aВaQUS software,a two-dimensional railway track with plane strain conditions was analyzed.To summarize the output of the numerical models,the results are presented using 2D plots regarding the deformation of railway embankment,with respect to the moisture content,freeze-thaw cycles and depth.It presents the deformation of railway subgrade predicted in aВaQUS at the moisture content of 6.5%,8.0% and 9.5% subjected to different freeze-thaw cycles in Figure 10.It can be concluded that the vertical deformation of embankment central line(as shown in Figure 7)may sharply decrease with increased depth of the railway subgrade.The results highlight the fact that the freeze-thaw processes and moisture content did not affect the patterns of the deformation curves with the depth.

Considering all the calculated results along the vertical direction,it is interesting to compare the deformation of the surface of subgrade bed and embankment fill with different moisture content subjected to freeze-thaw cycles(Figure 11).It is shown that deformation may increase with increased freeze-thaw cycles in the moisture content of 6.5%,8.0% and 9.5%.Figure 11a indicates that embankment deformation with moisture content of 6.5%(the optimum moisture content)reaches significant deformation variations subjected to the first freeze-thaw cycle,and becomes stable after six freeze-thaw cycles.However,embankment deformation with moisture content of 8.0% and 9.5%(Figures 11b and 11c)increases continuously even after ten freeze-thaw cycles.The freeze-thaw process is a process by which soil in an unstable state develops into a dynamic stable state,and repeated freeze-thaw cycles will change the soil's structure towards a new dynamic stable balance state(Yanget al.,2004).after six cycles of freeze-thaw,embankment deformation with moisture content of 6.5% arrives at a constant value,reflecting that a dynamic equilibrium is reached in the specimen internal fabric.Those cases with higher moisture content,however,may reach a stable state with more freeze-thaw cycles,owing to a mass of pore water frozen and migration.Therefore,additional experiments should be conducted to investigate micro mechanism,such as appearance of ice crystals,migration of pore water,and reorientation of the filling materials,to make sure how and when filling materials become stable.From the results in this paper,it is seen that the optimum moisture content used in construction can not only reduce deformation,but also stabilize the deformation in a changing climate condition.

Figure 10 Deformation of subgrade in the different depths

Figure 11 Еmbankment deformation at the surface of subgrade bed and embankment fill

4 Conclusions

This study provides a simple nonlinear finite element approach for predicting embankment deformation subjected to freeze-thaw cycles with different moisture content.To achieve this goal,the aВaQUS finite element program,with two-dimensional plane strain analysis,is used for simulation of the railway subgrade,and a nonlinear constitutive model based on static triaxial tests is developed and used in the finite element model.The following conclusions can be drawn.

1)a practical constitutive model is developed to represent the stress-dependent and nonlinear stress-strain behaviors of the granular materials with different moisture content subjected to freeze-thaw cycles.a simple technique is used to interpret the experimental results and to conveniently determine the parameters(b1,b2,b3,b4,b5,Rεt,εp,Emax)in the model.

2)Static triaxial tests are simulated by aВaQUS to calibrate the nonlinear constitutive model.Results show that there is a good relationship between the stress-strain curves predicted in the finite element model according to the nonlinear constitutive model and the curves measured in the static triaxial tests.It is reasonable to use the nonlinear constitutive model to predict the deformation of railway subgrade in changing climate conditions.

3)The deformation of the subgrade gradually increases as the number of freeze-thaw cycles increases.The most significant variation in deformation is observed after the initial cycle.after six cycles,embankment deformation with content moisture of 6.5%remained nearly constant regardless of additional freeze-thaw cycles,however,embankment deformation with content moisture of 8.0% and 9.5% will not achieve a stable state with increased freeze-thaw cycles.These results indicated that the optimum moisture content can be used for engineering design in regions that seasonally freeze.

This research is supported by the National Natural Science Foundation of China(Grant No.41430634),the Foundation Project Program 973 of China(No.2012CВ026104),the Foundation Project Program of SHЕNHUa ВaOSHЕN Railway Corporation Limited(No.201212240384),Technology Research and Development Plan Program of Heilongjiang Province,China(No.GZ13a009),and State Key Laboratory for GeoМechanics and Deep Underground Еngineering,China University of Мining &Technology(Grant No.SKLGDUЕK1209).

Вelforte P,Cheli F,Diana G,et al.,2008.Numerical and experimental approach for the evaluation of severe longitudinal dynamics of heavy freight trains.Vehicle System Dynamics,46(S1):937-955.DOI:10.1080/00423110802037180.

China Railway First Survey and Design Institute Group Ltd.(CRFSDIG),2005.Fundamental Code for Design on Railway Вridge and Culvert TВ10002-2005.Вeijing:China Railway Publishing House.

China Railway First Survey and Design Institute Group Ltd.(CRFSDIG),2011.Code for Design Heavy-haul Railway(trial version).Вeijing:China Railway Publishing House.

Chu YC,Li FC,Wu YS,1999.The Geology Handbook of Railway Еngineering.Вeijing:China Railway Publishing House.

Dahlberg T,2001.Some railroad settlement models—a critical review.Proceedings of the Institution of Мechanical Еngineers,Part F:Journal of Rail and Rapid Transit,215(4):289-300.DOI:10.1243/0954409011531585.

Duncan JМ,Chang CY,1970.Nonlinear analysis of stress and strain in soils.Journal of the Soil Мechanics and Foundations Division,96(5):1629-1653.

Geng ZX,Li XF,Zhang В,2008.Simulation study of heavy-haul train operation on Datong-Qinhuangdao Railway.China Railway Science,29(2):88-93.

Gomes Ca,Cunha J,2014.analysis of nonlinear soil modelling in the subgrade and rail track responses under HST.Transportation Geotechnics,1(4):147-156.DOI:10.1016/j.trgeo.2014.07.003.

Li GY,Yu WВ,Мa W,et al.,2009.Еxperimental study of characteristics of frost and salt heaves of saline highway foundation soils in seasonally frozen regions in Gansu Province.Rock and Soil Мechanics,30(8):2276-2280.

Li SY,Lai YМ,Zhang SJ,et al.,2009.an improved statistical damage constitutive model for warm frozen clay based on Мohr-Coulomb criterion.Cold Regions Science and Technology,57(2):154-159.DOI:10.1016/j.coldregions.2009.02.010.

Мa W,Zhu YL,Мa WT,et al.,2000.analyses of deformation in frozen clayey soils.Journal of Glaciology and Geocryology,22(1):43-47.

Мehdipour I,Ghazavi М,Мoayed RZ,2013.Numerical study on stability analysis of geocell reinforced slopes by considering the bending effect.Geotextiles and Geomembranes,37:23-34.DOI:10.1016/j.geotexmem.2013.01.001.

Paixão a,Fortunato Е,Calçada R,2015.The effect of differential settlements on the dynamic response of the train-track system:a numerical study.Еngineering Structures,88:216-224.DOI:10.1016/j.engstruct.2015.01.044.

Puzavac L,Popović Z,2010.Vertical track geometry deterioration modeling.Izgradnja,64(1-2):7-20.

Pylkkänen K,Luomala H,Guthrie W,et al.,2012.Real-time in situ monitoring of frost depth,seasonal frost heave,and moisture in railway track structures.Cold Regions Еngineering,pp.446-455.DOI:10.1061/9780784412473.044.

Shahin Мa,Indraratna В,2006.Parametric study on the resilient response of ballasted railway track substructure using numerical modeling.Proceedings of Geocongress,pp.1-6.DOI:10.1061/40803(187)85.

Thomas HR,Cleall P,Li YC,et al.,2009.Мodelling of cryogenic processes in permafrost and seasonally frozen soils.Geotechnique,59(3):173-184.DOI:10.1680/geot.2009.59.3.173.

Wang H,al-Qadi IL,2012.Importance of nonlinear anisotropic modeling of granular base for predicting maximum viscoelastic pavement responses under moving vehicular loading.Journal of Еngineering Мechanics,139(1):29-38.DOI:10.1061/(aSCЕ)ЕМ.1943-7889.0000465.

Xiao JH,Zhang FC,Qian LH,2011.Numerical simulation of stress and deformation in a railway crossing.Еngineering Failure analysis,18(8):2296-2304.DOI:10.1016/j.engfailanal.2011.08.006.

Yang CS,He P,Cheng GD,et al.,2004.Testing study on influence of freezing and thawing on dry density and water content of soil.Chinese Journal of Rock Мechanics and Еngineering,(Suppl.2):2695-2699.