Energy Method Solution for the Vertical Deformation of Longitudinally Coupled Prefabricated Slab Track

Upwarping on the longitudinally coupled prefabricated slab track system caused by the rising temperature is a common distress, which deteriorates the mechanical properties of the coupled slabs and the vertical stability of slabs. The objective of this paper is to quantify the upwarping deformation on the slab subjected to temperature force and to find out the influence of different factors on the upwarping phenomenon of the slabs. An analytical expression is deduced using energy method, and a finite element model is also established to verify the analytical solution’s adequacy. The following main findings are drawn: (1) when the amplitude of the initial elastic misalignmentwe is equal to a half of the amplitude of the initial plastic misalignmentwp and the half-wavelength lmin takes the most unfavorable value, the maximum relative error between the analytical solution and the result in FEM is only 2.64%, which demonstrates that the analytical solution correlates well with the FEM results. (2) lmin is closely related withwp. With the increase of wp, lmin becomes longer, and the maximum length of the half-wavelength is 7.769m. (3) When the total amplitude we + wp exaggerates, the slab will be much likely to suffer upwarping.


Introduction
Since slab track is characterized by its high stability, availability, reliability, and low maintenance associated costs, it has been a developing trend for high-speed railways [1].By September in 2016, China has completed the construction of 20,000 km high-speed railways among which the slab track technology accounts for about 78%.The longitudinally coupled prefabricated slab track (LCPST) which is reinnovated from the prototype of Bögl slab of Germany by China has been widely used on high-speed railway or Passenger Dedicated Lines (PDLs) in China with the total length of about 4637 km (double track).In 2008, this kind of track was firstly put into operation on the Beijing-Tianjin intercity PDL and gradually popularized on Beijing-Shanghai High-Speed Railway, Shanghai-Hangzhou Railway, Shijiazhuang-Wuhan PDL, and so forth in the following years [1].So far, the first LCPST has been served for almost 9 years in China.
The LCPST design uses prefabricated concrete slab made of reinforced concrete which is 0.2 m thick, 2.55 m wide, and 6.45 m long.The main components of the LCPST system include rail CN60 (Chinese standard rail with rail mass per meter 60 kg/m), WJ-8 type of elastic fastening system invented in China, prefabricated slabs, cement asphalt mortar (CA mortar), and hydraulic material layer, which are displayed in Figure 1.In this system, breaking points are notched every 650 mm at the face of the slab to form dummy joints.The threaded steels of the prefabricated slab are coupled in a longitudinal direction by tumbuckles [2]; after that sealing concrete is cast into the space between the slabs to form broad joint and narrow joint as shown in Figure 2. The width of a broad joint is 0.2 m.The weight of a slab is about 9 t (without rail fastening and rails).
As can been seen, the LCPST stands out its advantages not only in mechanical and structural aspects but also in technology's.For instance, preferable continuity or integrity is insured due to tumbuckles.Rail ditches are processed with the help of Computer Numeric Control (CNC) grinder with high precision to guarantee the accurate adjustment [1].Exact level and alignment are measured efficiently owing to an advanced adjustment device and there is a same structural type of LCPST at subgrade, bridge, and tunnel.Nevertheless, external factors such as differential drying shrinkage and creep of the material and temperature contribute a straight slab to initial upward curl.The slab and CA mortar of LCPST are bonded together due to the bonding behavior of CA mortar.However, this kind of bonding is deteriorated from the long-term operation; therefore the longitudinal forces in the coupled slabs caused by an increase of temperature may result in irrevocable loss of bonding strength.This loss leads to diminished mechanical performance between the slab and CA mortar, and sometimes forms debonding, which is displayed in Figure 3.Because the concrete at broad joint and the coupled slab are not cast in synchrony, a weak location can be found, and then cracking at this place is most likely to occur.If slabs with initial misalignments experience higher temperature, the top surface of slabs will elongate, while the bottom will compress, causing coupled slabs to a potential upwarping, which is demonstrated as a practical problem in LCPST.In the summer of 2013, surveys about upwarping at Beijing-Shanghai PDL and Ningbo-Hangzhou High-Speed Railway have been done by the research team of the author from Southwest Jiaotong University; both  showed that upwarping of slab occurred preferably around broad joint.Upwarping of coupled slab at bridge at Beijing-Shanghai PDL was measured to be approximately 5 mm in total from its initial configuration, whereas it was about 20 mm at K170 km + 545 m (the mileage of this position) of Ningbo-Hangzhou High-Speed Railway (see Figure 4).Slab upwarping not only may influence the geometric alignment of railway and the ride comfort but also brings about potential danger to the bearing capacity of track structure.Thus, research on this particular topic is required.Some analytical models have been developed to study upwarping of the coupled slab.Ren et al. preliminary analyzed the causes of upwarping of continuously reinforced concrete slabs for a twin-block slab track on the subgradebridge transition section and sorted out a repair scheme for upwarping of track slab through data analysis [3]; Liu and Feng put forward the repair measure for upwarping defect of LCPST [4] longitudinal force of Bögl longitudinal connected ballastless track on high-speed railway bridge [8].
Previous mentioned research work mostly covered FEM and experiment; the analytical solution has been rarely applied to study on the upwarping deformation of LCPST, but similar problem on the buckling of Continuous Welded Rail track has been investigated a lot using the analytical solution.The main difference between models of upwarping of slab track and buckling of Continuous Welded Rail (CWR) consists in different component of potential energy; nevertheless the analysis methods hold many similarities.A rise of temperature can generate compressive stresses in the rails in CWR tracks.Since the longitudinal expansion of the rails is constrained, generated compressive stresses may cause track buckling in the horizontal plane.Early in 1978, Kerr has presented a beam model to calculate the buckling load of CWR [9].Martínez et al. present an analytical model to obtain the buckling load of a CWR track [10].Luo et al. studied the relationship between the rail's longitudinal temperature force and the natural frequency [11].Manap has explored the changes of longitudinal stress distribution in the rails due to mechanical loading applied by a travelling train [12].The stability equation of CWR track (unified formula for short) in China has been derived since the 1970s based on the principle of stationary potential energy [13].And the equation which has been compiled in the Code for Design of Railway Continuous Welded Rail [14] works well in practical applications.
Unlike the previously developed studies on the upwarping of LCPST, which mainly rely on FEM or experiment, this paper is based on the minimization of potential energy and takes coupled slabs with initial misalignments into consideration.The analytical expression is obtained based on the energy method.

Methodology
The main difference between models of upwarping of coupled slabs and buckling of CWR consists in different components of the total potential.However, the analysis method and the initial state hold many similarities.The slab under study is slender; therefore, it can be modeled as an Euler-Bernoulli beam.Both models share the similar flexural behavior when both are subjected to temperature forces in the longitudinal direction.The flexural behavior of the coupled slabs model is shown in Figure 5.
2.1.Theoretical Background.Unified formula has been applied to analyze the stability of Continuous Welded Rail, and it works well in practice.It is assumed that the coupled slabs are considered as a beam with two initial misalignments (the initial elastic misalignment and the initial plastic misalignment), which share the same wavelength with the length of the coupled slabs, lying on the homogeneous medium.During the buckling process, the compression deformation energy of rail, deformation energy of track frame, and transverse deformation energy of ballast bed were calculated, respectively [15].Applying the principle of stationary potential energy, the unified formula was released by MOR (the former Ministry of Railway of China) in 1977 [13].As for the coupled slab, the hypotheses are the same as the rail's previously mentioned.The upwarping load of compressed slabs is the one for which the slabs find themselves in a neutral equilibrium state.For this reason, according to the energy method, the value of the upwarping load can be obtained from the system's total potential , which is a function of the generalized coordinates of the structural system.

Model Assumptions.
While the debonding of CA mortar occurs, CA mortar does not work well and has terrible bonding behavior, so that the coupled slab and CA mortar are separated.In this paper, coupled slabs are solely focused on, and the effect of CA mortar is ignored.For simplicity, there are some assumptions which are given as follows: (a) The effect of rail, fastening system is all neglected.
(b) During the upwarping process, it is assumed that the ends of the slabs remain where they stay, with hinge support at both ends.
(c) It is assumed that coupled slabs and CA mortar are completely separated in the debonding zone.
(d) The chord length  is assumed to be unchanged during the process of upwarping.
(e) The influence of adjacent broad joints is not taken into account in the analysis.(f) The stiffness of a broad joint is considered as a reduction of the stiffness of the coupled slab.
(g) The concept of the initial misalignments is similar to that in the model of buckling; both include an initial elastic misalignment and an initial plastic misalignment and share the same wavelength with the beam.
Just as previously mentioned, coupled slabs can be simulated as an equivalent beam, with the same cross section and inertia as a real slab.Surveys at Beijing-Shanghai PDL and Ningbo-Hangzhou High-Speed Railway revealed that upwarping of the coupled slabs occurred preferably around broad joint, where there is a stiffness difference between the prefabricated slab concrete and the cast in-site concrete, so the upwarping zone ranges from A to B in Figure 6.It is hypothesized that, during the upwarping process, the coupled slabs will deform by amplifying the initial misalignments (see Figure 6).Upwarping of slabs holds similar behaviors to buckling of CWR, especially in the regularity of deformation.For the buckling model, the beam is loaded by axial compressive forces and is elastically supported by laterally distributed resistance which is proportional to the lateral deformation.While the compressive forces develop, the beam will deform by amplifying the initial misalignments.For the upwarping model, the beam is loaded by axial compressive forces and is subject to a vertical distributed resistance (dead weight) which is proportional to the vertical deformation, and the deformation form is similar to that of buckling.In [16], Liu and Guang compared different types of deformation curve in the unified formula, and the result of calculation shows that sinusoidal curve is the most similar to the conclusions from German railway, so sinusoidal curve is used to calculate the stability of CWR.Also, the vertical deformation of the coupled slabs   under the action of the compressive temperature force is assumed to be a sinusoidal curve with the following equation: where  denotes the amplitude of the upwarping deformation and  denotes the half-wavelength of the slabs, just as illustrated in Figure 6.
Under the effect of temperature, the misalignment shown in Figure 6 will occur on the slabs.The initial elastic misalignment   is taken as where   denotes the amplitude of the initial elastic misalignment.
The coupled effect of the manufacture error, creep, and thermal expansion of concrete will cause unrecovered deformation of the slab, which results in an initial plastic misalignment.With the amplitude   , the expression of the initial plastic misalignment   is The function of vertical curve is to provide a gradual transition from one grade to another for smooth riding.Vertical curve represents one kind of initial misalignments.We can find in [17] that the arc is assigned to use in high-speed railway design in China.Since the sagitta of vertical curve is an irrational expression and is difficult for integration, an approximate equation is taken instead.With the radius of vertical curve , the initial plastic misalignment (  ) could be expressed as Thus, taking the initial elastic and plastic misalignments of the coupled slabs into consideration, the total deformation curves of LCPST with temperature force ( total ) or without temperature force ( 0total ) are, respectively, (

Model Description.
Considering all the hypotheses previously mentioned, the model can be accomplished using a compressive temperature load  which is applied at both ends of the misalignment and causes additional bending moment  to resist the upwarping.The beam (the blue one in Figure 7) with consideration of its gravity in the calculated model, which represents the coupled slab, has a half-wavelength , a uniform cross section area , a moment of inertia , and a modulus of elasticity .Besides, we simulate the broad joint using the second beam (the yellow one in Figure 7), which has a length Δ (Δ = 0.2 m), a same cross section area as the coupled slab, a moment of inertia  1 , and a modulus of elasticity  1 .And  denotes the force per unit length generated by the gravity of coupled slabs.
The total potential energy  of the model is the sum of the deformation energy associated with the compressing the slabs  1 , the energy associated with bending the slabs  2 , and gravitational potential energy of slabs  3 .According to the previous assumptions, the parameters that describe the deformed shape of the structural system are the halfwavelength of the slabs  and the amplitude of the upwarping deformation .The total potential energy of the model based on principle of minimum stationary potential energy [13] is obtained by Since assumption (d) is given, that is,  = 0, hence the equilibrium equation can be written as

The Deformation Energy Associated with Compressing the Slabs W 1 .
Once the model is subject to axial compression, the compression deformation energy of the slabs can be expressed as The parameter Δ is the change in the arc length (not the chord length) of the coupled slabs, Δ = Δ  − Δ 0 , in which Δ  represents arc-chord difference after deformation and Δ 0 represents arc-chord difference before deformation.According to Figure 8, we have in which   is arc-differential after deformation.Using a binomial series expansion for the quadratic terms and neglecting higher-order terms, we have Similarly, where  0 is arc-differential before deformation.
Integrating the difference   −  0 over the beam's span, the resulting deformation Δ is The energy becomes

The Energy Associated with Bending the Slabs 𝑊 2 .
During the upwarping process of coupled slabs, deformation energy could be stored inside and can be written as (1/2) ∫ , where  is the angle of rotation induced by the upwarping of the slabs [15].The energy associated with bending the slabs is composed of two parts, first of which is the energy of deformation ∫  0   , stored in the beam because of the deformation caused by   (the internal bending moment of the initial elastic misalignment), while the other is the energy of deformation (1/2) ∫  0   , stored due to the increase of the internal bending moment   during the process of upwarping.Hence, where  =    ,   =    , and EI represents flexural rigidity of the slab (EI will be replaced by  1  1 at the zone of the broad joint) and  represents the reduction factor of flexural rigidity of track slab due to the influence of dummy joints, and the formula is obtained: 2.3.3.Gravitational Potential Energy of Slab  3 .As mentioned in Section 2.3, the gravity of the slabs is assumed to be the distribution force along the longitudinal direction of slabs.The gravitational potential energy of track slab could be expressed by where  stands for the density of the slab:  = 2500 kg/m 3 . is the area of the cross section,  = 0.51 m 2 , and  is gravitational acceleration:  = 9.8 m/s 2 .

Analytical Solution.
Substituting ( 14), (16), and ( 18) into (7), the calculation formula of temperature force is  =  ( 1  1 − ) (Δ +  sin (Δ/))  4 ( +   ) /4 5 +  5 ( +   ) /4 4 + According to Code for Design of High Speed Railway (TB 10621-2014) [18], 20000 m or 25000 m is selected as the radius of vertical curve depending on the speed of passenger train.Considering the most unfavorable condition, we choose 20000 m.We take  as the length of a slab.Because the slab is made of concrete C55 (concrete grade), the elastic modulus of the slab is 35500 MPa.The cross section of the slab is a 2.55 m × 0.02 m rectangle, so the inertia of the slab is 1.7 × 10 −6 m 4 . = 35500 × 10 6 × 1.7 × 10 −6 = 60350 kN⋅m 2 .Based on assumption (f) in the second section of the second chapter, Model Assumptions, we take half stiffness of the slab as the stiffness of broad joint ( 1  1 = 0.5).
Table 1 shows that the damage of the broad joint (0.5EI) influences just 1.7% of the total temperature force compared with the nondamage of the broad joint.Also the radius of the vertical curve ( = 20000 m) influences just 1.78% of the total temperature force compared to the slab without the vertical curve ( = ∞).
It could be determined from tentative calculation that the damage of the broad joint and the vertical curve provides less contribution to the temperature force.So they could be ignored in the following verification procedure, and the formula could be simplified as follows: It could be obtained from (20) that  is closely related to half-wavelength (), amplitude of the initial elastic misalignment (  ) and amplitude of the initial plastic misalignment (  ), amplitude of the deformation (), bending rigidity of the track slab EI, and dead weight of track slab within the half-wavelength range.However, parameters such as   and   are hard to be ascertained.To solve this problem, we take the stationary value of the compression into account.
Since  is a function of the compression  and  is a parameter representing the half-wavelength, (20) represents a functional relationship among the axial compression, the upwarping half-wavelength, and the deformation of the slabs.From (20), it can be found that when the amplitude of the deformation remains unchanged, the longer the halfwavelength develops, the less the energy associated with bending the slab is and also the more the gravitational potential energy of track slab is, and vice versa.Therefore, it can be inferred that there must be a most unfavorable half-wavelength of deformation  min which corresponds to a minimum of the temperature force.According to the definition of functional derivation, the half-wavelength could be obtained by and hence The upwarping deformation has a maximum value with respect to the rising of temperature force.The following formula is obtained by taking derivative of temperature force  as a stationary value: We can derive When   =  min , the relationship among the amplitude of the deflection of the coupled slabs, the amplitude of the initial elastic misalignment, and the amplitude of the initial plastic misalignment is   =  +   .Since the parameters   and   are variable, we assume that  =   to a certain condition.Once  is given, the other two parameters   and   can be determined.Substituting the above conditions ( =  min ,  =   , and   =  +   ) in (20), if the temperature force  is taken as a stationary value, then we have

Formula Verification and Parametric Study
The primary objective of this section is to quantify the upwarping in coupled slabs and determine its effect on the bending behavior of the coupled slab.Finite element method was utilized to verify whether the analytical solution correlates with the results of the FEM.As mentioned previously from the survey at Ningbo-Hangzhou High-Speed Railway in 2013, the maximum value of the amplitude of the total deformation   +  +   may reach 20 mm.Since we have such a relationship,   =  +   , we can find that 2  = 20 mm, so   = 10mm.On the other hand, the initial plastic misalignment is caused by the coupled effect of the manufacture error, creep, and thermal expansion of concrete, and such a misalignment in the coupled slabs holds the same reasons as the CWR.Since these adverse effects have a greater influence on the coupled slabs and the amplitude of the initial plastic misalignment in CWR is 3 mm [15], we calculate the temperature force with the amplitude of the initial plastic misalignment in LCPST starting from 4 mm.Therefore, the amplitude of the initial plastic misalignment   in this section ranges from 4 mm to 10 mm.

Model and Parameters in ANSYS.
The slab with liner density of 2500 kg/m 3 is applied and simulated by the continuously supported Euler beam.The element type of beam 188 in the finite element model is adopted in ANSYS software, which is suitable for analyzing slender to moderately stubby/thick beam structures.The element is based on Timoshenko beam theory which is a first-order shear-deformation theory: transverse-shear strain is constant through the cross section.This element uses the slenderness ratio of a beam structure ( 2 /()) to judge the applicability of the element, where  is shear modulus; it can take a value between E/3 and E/2, and other parameters are mentioned previously.A slenderness ratio greater than 30 is recommended; our model's slenderness ratio is 416.025,far more than 30, so it is suitable for simulating the slab.As shown in Figure 9, the mesh density of the slab in length is 10 times of the slab length, and those in the height and width are 18 and 36, respectively.Thus, the total number of elements used in the FE model is 10 ×  (length of the model) × 18 × 36.The nonlinear spring is used to simulate the supporting between the slab and CA mortar.COMBIN39 is a unidirectional element with nonlinear generalized forcedeflection capability that can be used in any analysis.So the supporting of CA mortar can be simulated with COMBIN39 spring element, with the thickness of 0.03 m.The elastic modulus of the supporting of CA mortar with 8,000 MPa [19] is applied in the calculation model.The supporting area of CA mortar is equal to the bottom of ×2.55 (slab width) m 2 .With the number of divided elements of track slab , the stiffness of the spring  yields Hence, the relationship between the supporting force of the spring and displacement  is where   is the interaction force between the supporting of the CA mortar and the coupled slab and  is the relative displacement between the supporting of the CA mortar and the coupled slab, which is negative when the spring is compressed.The lower end of spring is applied with full restraint, as displayed in Figure 10.The elastic modulus of prefabricated slab  is taken as 3,5500 MPa, and the thickness and width of slab are taken as 0.2 m and 2.55 m, respectively, and the coefficient of linear expansion  is 10 −5 .The calculation model is shown in Figure 10, in which the slabs are subject to both the effect of temperature and the dead weight.

Formula Verification.
With the initial conditions of   = 2 mm,   = 4 mm, and  = 2 mm, the most unfavorable half-wavelength of the deformation  min of 6.18 m is calculated out according to (22) or (24).And we can also get the temperature force  (15,383.377kN) by using (25); the value of the temperature force is converted to temperature and then applied in the finite element model with the initial misalignments of   +   .According to the comparison between the obtained vertical maximum displacement ( = 1.995 mm) and the given parameter  in the analytical solution, the relative error is only 0.25%.From the compared results displayed in Table 2, we can obtain the fact that the calculation result of derivation formula is in accordance with the calculation result of ANSYS. Figure 11 shows the contour plots of ANSYS calculation results.

Parametric Study.
The presented model was used to perform a parametric study and to evaluate the influence of the considered parameters on the upwarping vertical deformation of LCPST.The study in this section analyzed the effect of the amplitude of the initial elastic misalignment   , the amplitude of the initial plastic misalignment   , the amplitude of the deformation of the slabs , and the most unfavorable half-wavelength  min .It could be seen from Table 2 that the relative error between parameter  for formula calculation and deformation output of finite element method is relatively small.In Table 2, the changes of the parameters (  ,   ) do not vary in random, they affect each other, and  min is the value with respect to   .In general, they represent a series of peak values, which consist in a series of discrete points, which indicates the deformation of track slab, and other parameters,   ,   , and  min , in the formula are not independent, and the deformation of the upwarping is getting larger along with the increase of initial misalignments.
Figure 12 and (24) show that the most unfavorable half-wavelength  min is closely related to the amplitude of the initial plastic misalignment.The most unfavorable halfwavelength of deformation increases with the increase of initial misalignments.Since the upwarping generally occurs at broad joints, the length of one longitudinal coupled slab is 6.45 m and the calculated longest half-wavelength is 7.769 m, which is no more than the length of 2 slabs; hence, the influence of adjacent broad joints could be ignored.It could be seen from Table 3 that when the amplitude of the initial elastic misalignment   is 2 mm and initial plastic misalignment's   is 4 mm, the most unfavorable halfwavelength  equals 6.18 m; the temperature force necessary for upwarping of the track slab up to 2 mm is 15383.377kN; when   is 5 mm and   is 10 mm,  equals 7.77 m; the temperature force necessary for upwarping of the track slab up to 5 mm is 9729.302kN, which illustrates that the more terrible the initial misalignments defect of track slab is, the longer the most unfavorable half-wavelength is and also the more unstable the track slab is; that is, unwarping is more likely to occur.Therefore, it is recommended to reduce the initial misalignment caused by the coupled effects of creep, thermal expansion, and manufacture error of track slabs to The amplitude of the initial plastic mislignment w p The most unfavorable half-wavelength 0.004 0.005 0.006 0.007 0.008 0.009 0.010 6.0  prevent the coupled slabs from large upwarping within a smaller temperature force.

Conclusions
In this paper, a study on the effect of upwarping on the longitudinally coupled prefabricated slab track was conducted by establishing an analytical expression using energy method which applies energy method to quantify the vertical upwarping deformation of the coupled slabs.The main contribution of the expression is that it can quantify the vertical upwarping of the coupled slabs with both of the gravity of the slabs and the initial misalignments of the slabs taken into account.
Then a finite element model of the coupled slabs for the upwarping was built to demonstrate the analytical solution's adequacy.The analytical solution may accurately predict the observed behavior of the presented problem.Based on the limited results in this study, several conclusions can be drawn as follows: (1) Based on the results between the analytical solution and the solution in FEM, it can be concluded that when the amplitude of the initial elastic misalignment   is equal to a half of the amplitude of the initial plastic misalignment   , the relative error between the deformation parameters from formula calculation and deformation output from FEM is relatively small (the maximum value of the relative error is 2.64%).The response predictions from the analytical solution correlate well with the solution in FEM.(2) The most unfavorable half-wavelength of the deformation  min is closely related with the amplitude of the initial plastic misalignment   .With the increase of the amplitude   , the most unfavorable half-wavelength  min becomes longer.The calculated maximum half-wavelength is 7.769 m, thus less than the length of 2 slabs (the length of LCPST is 6.45 m), so the influence of adjacent broad joints could be ignored, which is in accordance with assumption (e) in Section 2.2.(3) When the amplitude of the initial misalignments   +   exaggerates, the vertical deformation of the slabs will be amplified with a smaller temperature force than that of the slab with smaller initial misalignments; namely, the slab with larger initial misalignments will be much likely to suffer upwarping.It is thus recommended to reduce the initial plastic misalignment caused by the coupled effects of creep, manufacture error, and thermal expansion of the coupled slabs to prevent the coupled slabs from large upwarping within a smaller temperature force.Compressive temperature load :

Nomenclature
Total potential energy of the slabs  1 : Deformation energy associated with compressing the slabs  2 : Energy associated with bending the slabs  3 : Gravitationalpotentialenergyoftheslabs : Cross section area of the slabs Δ: Change in the arc length (not the chord length) of the coupled slabs Δ  : Arc-chord difference after deformation Δ 0 : Arc-chord difference before deformation  0 : Arc-differential before deformation   : Arc-differential after deformation : Internal bending moment during the process of deformation   : Internal bending moment of the initial elastic misalignment : Angle of rotation induced by the upwarping of the slabs : Flexural rigidity of the slab  1  1 : Flexural rigidity of broad joint : Force per unit length generated by the gravity of coupled slabs Δ: Length of broad joint : Reduction factor of flexural rigidity of track slab due to the influence of dummy joints   : The interaction force between the supporting of the CA mortar and the coupled slab.

Figure 5 :
Figure 5: Flexural behavior of the coupled slabs model.

Figure 8 :
Figure 8: Arch-chord difference of the coupled slabs before and after deformation.

Figure 9 :
Figure 9: Mesh density of the slab.

Figure 10 :
Figure 10: Finite element model for the coupled slabs with initial misalignments in ANSYS.

Figure 12 :
Figure 12: Relationship between the amplitude of the initial plastic misalignment and the most unfavorable half-wavelength of deformation.

Table 1 :
The tentative calculation about the damage of the broad joint/radius of the slab.

Table 2 :
Comparison of vertical deformation  between analytical solution and the result of the FEM.

Table 3 :
Relationship among the amplitude of the initial misalignments, most unfavorable half-wavelength, and temperature force.