7.0 44 (2) Ministerie van Verkeer en Waterstaat Dienst Weg- en Waterbouwkunde Directoraat-Generaal Rijkswaterstaat EQUIVALENT LAYER THEORIES "State of the art report" W-DWW-94-904 NIET UITLEENBAAR BUITEN DWW BIDOC DWW - 015-2518363 B I D O C (bibliotheek en documentatie) Dienst Weg- en Waterbouwkunde Postbus 5044, 2600 GA DELFT Tel. 015 -2518 363/364 1 .Ü-UH I I I 1 1 I I EQUIVALENT LAYER THEORIES "State of the art report" A.C. Pronk Road and Hydraulic Engineering Division Rijkswaterstaat The Netherlands (bibliotheek en documentatie) I I Dienst Wfig- en Waterbouwkunde Postbus 5CM4,2S0O GA OL1.FI Te!. 015-2518 363/564 I i i f f l i i I I i Disclaimer (Dutch) Dit werkdocument wordt uitgegeven om geïnteresseerden de gelegenheid te bieden om van de voortgang van het desbetreffend onderzoek, e.d. kennis te nemen. Benadrukt wordt dat de gezichtspunten in dit werkdocument niet noodzakelijk overeen behoeven te komen met de officiële gezichtspunten of het beleid van de directeur-generaal van de Rijkswaterstaat Met de in dit werkdocument gegeven informatie dient derhalve met de nodige voorzichtigheid te worden omgegaan, aangezien de hierin vermelde conclusies in de loop van het verder onderzoek of anderszins mogelijk herzien dienen te worden. I I I I I I I i t i i § i i Abstract An overview is given on the most known equivalent layer theories. In the equivalent layer theory, overlaying (pavement construction) layers with different thicknesses and moduli are combined into one layer with an equivalent thickness Heq. Special attention is given to the applicability in case of calculations using the Hogg and Westergaard models for concrete pavements. It is shown, that in principle the equivalent layer thickness theories can not be used in the Westergaard model because this model (Winkler foundation) is not comparable with the linear elastic isotropic half space or multi layer model (Burmister model) on which all the equivalent layer theories are based. Only in the case of the "Supplementary Theory" model of Westergaard in which the reaction of the foundation is 'corrected' to that of a linear elastic half space (HoggJeuffroy model) a (unique) relationship can be established between the modulus of subgrade reaction k and the Youngs modulus Eo for both the deflection and the stress. Therefore the relationship given by Eisenmann : k = Eo / Heq is not correct. Contents 1. Equivalent Layerthickness Theory of Palmer and Barber 2. Odemark's Equivalent Layerthickness Theory 2.1 Deflections 2.2 Stresses and Strains 3. Ullidtz's elaboration of Odemark's Theory 3.1 General 3.2 Special application for a stress dependent subgrade modulus 3.2.1 Introduction 3.2.2 Strain calculations in a two layer system for a circular load 3.2.3 Calculations of deflections 4. Other variants on the Equivalent Layer Thickness Theory 4.1 General 4.2 Nijboer's Equivalent Layer Thickness Theory 4.3 Pronk's Equivalent Layer Thickness Theory 4.4 Comparison of the several approaches 4.5 Thenn de Barros Equivalent Layer Thickness Theory 5. Concrete Pavements 5.1 General 5.2 Two layer model 5.3 Three layer model References Appendix I ; Overview of formulas 1 i I 1. Equivalent Layerthickness Theory of Palm er and Barber In this not so well-known theory [1] the equivalency of the layers is based on a "bending" stiffness, that is to say, each horizontal infinite layer is characterized by: 1/3 L = H± (1) (1-v?) The behaviour of a layer denoted by: (Hir Eif v) is said to be equal of that of another layer with the parameters (Hkr Ekf vk) if the following relation is valid: HiQi = HkQk or 1/3 (2) 1/3 (1-v?) If the layers are numbered starting from above (surface : 1 ; subgrade : n) the procedure is as follows: The equivalent layer thickness (Heq,) of all layers above layer (i+1) with characteristics (H kv E M ,v k1 ) is calculated using equation (3). Heqd = H.1 Q, X x- gj-1 •H, x•i+1 (3) J=l The deflection W is now written as the summation of the compressions in all the layers above the subgrade and the deflection of the subgrade: j=n-l W = +Wn = Wn The compression in a layer with number i (i=1,2,... n-1) is calculated using the deflection formulas for a half space W{zrE,v} (see Appendix I) according to: (5) The deflection of the subgrade is calculated with: Wn = According to [1] this approximation of the deflection in a linear elastic multi layer system is quite sufficient if the ratio of the Youngs moduli of layer i and layer i+1 is bigger than 3 (E/El+1 > 3). 2. Ödemark's Equivalent Layerthickness Theory 2.1 Deflections In the equivalent layerthickness theory of Ödemark [2, 3] two new equivalent layer thicknesses are calculated for each pavement layer (except of course for the infinite subgrade). 1/3 (7) Heq1 = Geqi = f2Hi with f± and f2 are constants and E*+1 = F i EUk, Huk, \>i+k } (8) k = 1, 2, 3, . . up to the subgrade modulus Eo The second equivalent layerthickness Geq, is only used in the calculation of the (vertical) compression of.a construction layer. The first equivalent layerthickness Heq, is known as the equivalent layerthickness. Calculation of E-^ M and v^ ui Odemark did not give the actual calculation 'for the replacing Poisson ratio v'M because in his original work all the Poisson ratio's v, were taken equal. In view of the small influence of the Poisson ratio on the calculations it can be justified to put vM = vM. However, the calculation of the replacing modulus E*^ is quite complex. The calculation is based on the assumption that the deflection at the top of layer i can be calculated on one hand as the deflection in a half space with modulus E*, at a depth z,* = I[Heq k ] from k = 1 to k = (i-1) and on the other hand as the sum of the compression of layer i and the deflection at the bottom of layer i. I I I The compression of layer i will be calculated in this case as the difference between the deflections in a half space with modulus £•, at the depths z(* = §[Heqk] from k = 1 to k = (i-1) and z," = zs* + The deflection at the bottom of layer i is calculated as the deflection in a half space with modulus E*M at a depth z,'** = §[Heqk] from k = 1 to k = i. By performing this procedure for all layers a relation will be established between E„ E*, and E*^, for all values of i. Because the, for the layer j just above the subgrade, modulus E*H equals the subgrade modulus Eo all the replacing moduli E*, can be calculated. This procedure will be explained for a three layer model with the parameters E1r vv H v E2, v2, H2, Eo, vo. The deflection in a half space with modulus E and Poisson ratio v at a depth z is defined as W(z,E,v). For the first pavement layer the two equivalent layerthicknesses are defined by: 1/3 j i' Heqx = E v?) (9) Geq± = The two equivalent layerthicknesses for the second pavement layer (just above the subgrade) are given by: Heq2 = f±H2x E2 ( 1 - vg) 1/3 (10) Geq2 = f22H. "2 The compression of the two pavement layers are calculated as follows: AWX = W{0,E1,v1) - W {Geqx,E1,v1) (ïi) AW2 = W {Heq1, E2,v2) - W iHeq1 + Geq2, E2, v2) The deflection at the top of the subgrade is determined by: Wo = (12) The yet still unknown replacing modulus E2* is determined by the requirement that the deflection at the top of the second pavement layer (AW2+W0) equals the deflection in a half space with the modulus E 2 at a depth z* = Heq, : W {Heq1rE2\v2} AW2 + Wo = (13) The deflection W in a halfspace (£,v) at a depth z due a homogeneous divided circular load P with radius a is given by: Wiz,E, v) = with (1 + v) P cos(x) x 1 + - 2v naE 1 + s i n (x) tg{x) = — d. (14) The term in equation (14) which contaïns the factor (1-2v) disappears for v = 0,5. In that case (1^=1^=1^3=0,5) the requirement \N(Heqv£*2,vJ = AW2+WO will lead to equation (15): f{a2 + Heql) . E2 a (15) (Heg1 + Heq2)2) .Eo a (Heq1 + Geq2)2) y/( .E2 1/3 with Heq1 = f1H1x 1/3 Heq2 = fxH2 (16) 12. 'Jo Geq2 = f2H2 Usually a value of 0,9 is taken for the constants f, and f2. In table 1 the ratio E2'/Eo is giyen for several combinations of the ratio's E,/Eo, E2/Eo, H / a and H2/a. The ratio E2*/Eo can be seen as the factor by which the bearing capacity of the subgrade (E*2 is a kind of effective subgrade modulus) has been increased by the application of the foundation (H2,E2,v2). As shown in table 1 the ratio E/E,, has limited influence on the increase for low values of the ratio E2/Eo. The increase becomes even smaller if the stiffness modulus of the first layer E, increases. It also appears that the influences of the ratio's H/a and H2/a are mainly determined by the ratio H,/H2 which implies that the increase factor E27EO nearly depends on the radius a of the load plate. 2.2 Stresses and Strains For the calculations of the stresses and strains in the interface between two layers the equivalent layer thickness Heqi is used. The formulas for the stresses and strains in a halfspace (E,v) on the center line at a depth z due to a homogenous divided circular load P are denoted by : ar(z,v) : Radial stress = Tangential stress (center line) Vertical stress Radial strain = Tangential strain (center line) Vertical strain Complete friction is assumed between the pavement layers -> the horizontal strain at the bottom of layer i (erl) equals the horizontal strain at the top of layer i+1 (eM+1). Table 1. Combinations Increase factor E2*/Eo H/a H2/H, E,/Eo 1 2 0,5 40 1,10 1,19 1,28 1 2 0,5 20 1,12 1,24 1,35 1 1 1,0 80 1,15 1,29 1,44 1 1 1,0 40 1,18 1,35 1,53 1 1 1.0 20 1,21 1,42 1,64 2 2 1,0 80 1,16 1,31 1,46 2 2 1,0 40 1,19 1,38 1,58 2 2 1,0 20 1,23 1,47 1,73 2 1 2,0 80 1,27 1,57 1,90 2 1 2,0 40 1,32 1,67 2,09 2 1 2,0 20 1,36 1,79 2,30 H2/a E2/Eo = 2 E2/Eo = 8 E2/Eo = 4 The horizontal strain at the bottom of layer i is calculated as : k=i z= k=l The vertical strain at the top of layer i+1 (evi+1) is calculated as: k=i e v , i + i = e v ( z , £•;+!, v i + 1 } with z = £ Heqk (18) k=l The horizontal and vertical stresses at the top of layer i+1 (a-hM-, and <7vM ) are calculated as: r, i + 1 = oriz, v i+1 ) and k=i z (19) = k=i The horizontal strain at the bottom of layer i (sri) can not be calculated by: k=i-l ezi * £ r ( z , ^ i / v i } with z = J^ [Heqk] + Geq£ Jc=l (20 > The horizontal stress at the bottom of layer i (<ro) can be calculated with the aid of the following relations: for the lower layer i : and for the upper layer: (21) interface : : frietion continuity For a two layer system (HirE1f 1^=0,5; Eo=100 MPa,vo=0,5) the horizontal strain has been calculated at thé bottom of the first layer due to a circular load P = 50 kN. The En value was varied from 2000 to 8000 MPa and the ratio H/a was varied from 0,5 to 2,0. The results of the calculations are given in table 2 together with BISAR calculations. I I Table 2. cr>1 = £r {Heq1r 100,0,5} H/a E1=2000 [MPa] ^=4000 [Mpa] £r.i E1=8000 [Mpa] BISAR E1=2000 [MPa] E1=4000 [MPa] E1=8000 [MPa] I 0,5 1647 1320 990 780 531 343 1,0 705 482 320 421 257 152 I 2,0 209 135 86 161 93 53 1 I i 1 i i i 3. Ullidtz's Elaboration of Odemark's Theory 3.1 General Ullidtz [4] extended the theory of Odemark in one way by a more accurate determination of the correction factor f, in the equivalent layerthickness calculation of Odemark. Ullidtz denoted this factor as f. Instead of a constant value (= 0,9) as in the original Odemark's theory, the correction factor f depends on the ratio of the layerthickness H, the radius a and the obtained quantity. An overview for a two layer systeem (HvEvvvEo,vj is given in table 3. Table 3 Correction cofficients f All Hf/a values H,/a = 2/3 H/a = 4/3 W [um] [Pa] [m/m] [m/m] [m/m] [m/m] 0,9 0,9 1,25 1,00 1,00 0,90 H/a > 2 Gr [m/m] 0,90 [m/m] 0,85 A value for f of 0,8 is recommended in the case of three and more layer systems with one exception for the calculations of the strains and stresses in the first interface where a value of f = 1,0 is advised. However, it should be remarked that Ullidtz did not follow the original concept of Odemark for the determination of an equivalent cq replacing modulus E,' in the case of three and more layer systems. The implication is that the equivalent layerthickness of layer i depends only on the material properties of layer i+1 just below layer i and not on the layers below layer i+1. 1/3 Heqi = fxHix (22) Furthermore, in many cases Ullidtz neglects the influences of the Poisson ratios. Besides the variation in the correction factor f the elaboration of Ullidtz resembles more the equivalent layertheory of Palmer and Barber. However, the merit of the elaboration of Odemark's theory by Ullidtz is the fact that Ullidtz also developed an application of the equivalent layerthickness theory for multilayer systems for cases in which the moduli of the subgrade and foundation are stress dependent [4]. Next to this important application, Ullidtz developed a methodology for the calcultion of stresses and strains in the interfaces of multi-layer systems, due to several circular loadings (twin wheel). In this methodology, the pointload formulas of Boussinesq are used, but the equivalent layerthickness Heq (or actual depth z) is replaced by a fictive depth z \ in order to 'simulate' the influence of the actual circular load. The second term at the right hand side of the equation above, resembles the influence of the circular loaqV -Aa. iiffpjQyepieiQt'fer 3his transformation is Pronk [5]. Heq * = Heq x Heq+X2 a 3.2 Special application for a stress dependent subgrade modulus 3.2.1 Introduction The equivalent layer methodology of Ullidtz can also be used if the subgrade modulus Eo (= material modulus) is stress dependent i.c. the modulus can be described by: Eo = C0*x (25) in which <x, is the major principal stress, a0 is a reference stress of 0,1 MPa and C'o is a constant. The procedure for a two layer system is as follows: I I I I 1 I 1 I 1 I I I I I I 1 I I I I Surface modulus E*Oheq at the interface: n O. (26) eg EZ* 'eg = [ 1 - 2 2 3 ] xC0*x Equivalent layer thickness: (27) E,o.h,eg Major principal stress at the interface: 3 2 2 °1 h ~ ö ' = oT, n x 1 a - (28) • h ' e«r 1 + 2 a The calculation is as follows: 1) Take an estimate for heq. 2) Calculate cr,Mer 3) Calculate E'oheq. 4) Calculate a new estimate for heq. The correction coëfficiënt f is taken equal to 1 for the calculation of stresses and deflections. For the calculations of strains the following values depending on heq [m] should be used: ehorizontal ^vertical < 0.1 1,25 1,15 0,1<-<0,3 1,10 1,05 >0,3 1,00 1,00 If the construction consists out of more than one layer above the subgrade, the equivalent layer thicknesses should be calculated as: *^-E ii i • •EQ, 1 3 (29) 3.2.2 Strain calculations in a two layer model for a circular load The required formulas are: ov(0,z) = ov(0,0) (z/a) (z/a) x 1 - 1+v xo v (0,0) x (z/a) E th(0,z) x (z/a)2) 1 - ( l - 2 v ) x (1+ ( 1 + U / a ) 22)\ (30) (31) 2 = (32) For the calculation of the vertical subgrade strain the depth z is replaced by the equivalent layer thichness heq (calculated with the 'strain' values for the correction coëfficiënt f) and the (material) modulus Eo should be used. In case of the horizontal asphalt strain two procedures are possible if complete friction exists between the layers. In both methods the equivalent layer thichkness heq should be calculated with the f values for the horizontal strain. method 1: - Calculate a "vertical stress" av(0rheq). - Calculate a "vertical asphalt strain" ev(0,heq) using the asphalt modulus Ev - Calculate the horizontal asphalt strain eh(0rheq) from these figures. method 2: - Calculate a "vertical subgrade strain" ev(0,heq) using the subgrade material modulus Eo. - Calculate a "horizontal subgrade strain (= horizontal asphalt strain) eh(0,heq) using the subgrade material modulus EO. I I I I I I I I I I I I I I I I I I I I 3.2.3 Calculations of deflections For the calculation of the deflection the value for the correction coëfficiënt f is taken equal to 1. The deflection at the interface with the subgrade can be calculated with the following formula: xaxo7(0,0) ö(0,z) = (33) 1+(l-2v ) x \ For the modulus E the surface modulus E*Oheq should be used instéad of the material modulus Eo. For relatively large heq values (>a) it is possible to use the simple Boussinesq formulas for a point load P on top of a half space. To simulate the circular load the depth z should be replaced by the formula: 0, 6 a: (34) eq Using the point load formulas it is possible to calculate the strains and stress due to a twin wheel like in the Shell Pavement Design Manual. This procedure is used in the program "I LEAP DEEP". 4. Other Variants on the Equivalent Layerthickness Theory 4.1 General For the reduction or simplification of a multi layer system to a two layer system or even to a half space several other theories than given here are known. An extended overview is given by Poulos [6]. In this paragraph only the equivalent layerthickness theories of Nijboer [7], Pronk [8] and Thenn de Barros [9] are dealed with. 4.2 Nijboer's Equivalent Layerthickness Theory A pavement construction can almost be simulated by or modelled into a three layer system (H1,Evv1,H2rE2,v2,Eo,vo). When Nijboer developed his theory, the multi-layer model of Burmister was already known, but no computer programs like BISTRO or BISAR were evailable at that time. For the calculations of strains, stresses and deflections graphs were used. These graphs were developed by Odemark and were only valid for two layer systems. Nijboer developed a methodology for the combination of the two upper pavement layers to one layer in order to be able to use the Odemark's graphs. Nijboer used his theory (in which v1~v2=ve,} especially for the stiffness measurements which were carried out with ther Road Vibration Machine (RVM). The point of departure in all most each equivalent layerthickness theorie is the requirement that the stress distribution in the interface will be the same. Odemark tried to fulfill this major requirement by the requirement that the 'bending moments' (expressed by the quantity E.h3) of the original layers and the equivalent layer are the same. Nijboer used the same approach in the replacement of two layers by one equivalent layer but he used also a second requirement for the place of the neutral axis. His second requirement was that the actual depth of the neutral axis should also be the same. However, at this point it should be noticed that in contrast with ordinary beams, it is not possible to talk about a neutral axis in the case of a multilayer system. In that case several neutral "surfaces" can occur. The properties of the new equivalent layer (H',E\v'=v,=v2) are determined by the following two equations: = (A + N)2x A* + 4A3N + 6A2N + AAN + N2 x (A2 + 2A + N)3 x (A + N) . (N + 1 ) with A = —- and 2 N = H Hl E 2 Nijboer (I) Besides this replacement of the two pavement layers by one different layer (Nijboer I: H*,E*)r Nijboer also gave another equivalent one layer substitution (Nijboer II). In his second method, Nijboer dropped the second requirement for the place of the neutral axis and replaced it by the requirement, that the equivalent layerthickness H* should equal the sum of the layerthicknesses of the two original layers. His second approach is given by the following two equations: E* (II) = A ' + 6A N+ 4 M *A3N+ " 3 (A+l) (A + N) H* (JJ) + j y 2 xK O7) - K + H9 (38) Nijboer(II) 4.3 Pronk's Equivalent Layer Thickness Theory Pronk [8] gives a supplement on the methodology of Nijboer. For the second requirement, Pronk suggested that the modulus of the equivalent layer, should equal the modulus of the lower original layer. This requirement leads to the following two equations for the calculation of the equivalent layer (H*,E*): E* H* = (39) = + 6A2N + 4AN + N2 (A + l ) 3 x (A + N) 1/3 x (Hr + H2) (40) Pronk The last equation resembles more the equivalency as used in the theory of Odemark. Although these three methods (Nijboer (I), Nijboer (II) and Pronk) have been developed for the simplification of a three layer model into a equivalent two layer model, the methodology can easily be extended to more than three layers. The procedure is as follows: 1) The first (H^E,) and second (H2,E2) pavement layers are transformed into one new equivalent layer (H2\E2*). 2) This new layer and the next layer new equivalent layer (H3*,E3*). 3) The steps above are repeated until all layers above the subgrade (E,,) are transformed into one new equivalent layer (H*,, £'). are transformed into the following After this procedure the solution techniques for a two layer model can be used. 4.4 Compar/son of the several approaches A comparison has been made for the several approaches, using a four layer model of which the material properties and dimensions are given in Table 4a and 4b. Table 4a: Construction A "Full depth asphalt" Load P = 2 x 25 kN ; SPDM Complete friction between layers Layer number E, [MPa] Asphalt 1 2000 0,35 0,04 Asphalt 2 4000 0,35 0,04 Asphalt 3 6000 0,35 0,10 Subgrade o 100 0,35 v, [m/m] H, [m] For these two four layer models the horizontal strain (er) at the bottom of the lowest pavement layer (i=3) and the vertical strain (ev) at the top of the subgrade have been calculated for the several equivalent layerthickness approaches. The results are given in Table 5a and 5b Table 4b: Construction B "Asphalt on top of an unbound base" Complete friction between layers Load P = 2 x 25 kN ; SPDM H, [m] Layer number E, [MPa] Asphalt 1 2000 0,35 0,04 Asphalt 2 4000 0,35 0,04 Asphalt 3 6000 0,35 0,10 400 0,35 0,20 100 0,35 Base 4 Subgrade o Vi [m/m] oo Table 5a : Results for the full depth asphalt construction (A) e BISAR Nijboer (I) Nijboer (II) Pronk Ullidtz er 168 127 196 171 221 (185) ev 530 320 561 525 535 (455) Table 5b : Results for the asphalt construction on top of an unbound base (B) e BISAR er 122 85 138 126 112 (96) ev 385 212 377 373 280 (233) Nijboer (I) Nijboer (II) Pronk Ullidtz For the correction coëfficiënt f in the approach of Ullidtz a value of 0,9 has been used. The strain values (between brackets) are calculated if a value of f = 1,0 is adopted. Better comparison is reached if for construction A the values 1,065 (er) and 0,905 (ev) are used and the values 0,85 (er) and 0,745 (ev) for construction B. base. I I I I I I I I I I 1 I I I I I I I I I The replacing properties for the new equivalent layer in the approaches of Nijboer, Pronk and Ullidtz are given below: BlSAR Calculations for two layer system with Eo = 100 MPa and v0 = 0,35 Nijboer (I) Nijboer (II) Pronk E* = 1374 Mpa; E* = 3951 Mpa; E* = 6000 Mpa (=E3); H* = 0,405 m H* = 0,180 m (= H1+H2+H3) H* = 0,1591 m Calculations with "I LEAP DEEP" : Boussinesq formulas (point load) with a correction for the equivalent depth (circular load). Ullidtz 4a Heq= 0,6369 m (er asphalt and ev subgrade) Eeq= 100 MPa Ullidtz 4b Heq= 0,4012 m (er asphalt) Eeq= 100 MPa Ullidtz 4b Heq= 0,9543 m (ev subgrade) Eeq= 100 MPa In case of the "Nijboer" and "Pronk" formulas, the first step was to take the upper two asphalt layers together into an equivalent layer. The second step was to combine this equivalent upper layer with the lowest asphalt layer into a new equivalent layer. In the case of the "Pronk" formulas, this approach gives satisfying results for a general modelling of different asphalt layers (differences in E moduli due to differencesin temperature) in a common pavement type (asphalt, unbound foundation and subgrade) into one new equivalent asphalt layer [8]. 4.5 Thenn de Barros's Equivalent Layer Thickness Theory Next to the approaches developed by Nijboer and Pronk, another equivalent approach is developed by Thenn de Barros [9]. In his approach Thenn de Barros uses the requirement for the equivalent layerthickness H" as Nijboer (II), but the formula for the replacing stiffness modulus E* is different. Nevertheless, the answers for E* do not differ much, as can be seen in Table 6. n.1/3 v2 E* = H2 N1/3 + A 1 +A H* = 13 - (41) 3 v <E2 H2 Thenn de Barros (42) Table 6 Ratios E*/E2 for the Nijboer II and the Thenn de Barros model. Ratios E/E2 ; Y\2/H, E/E2 Hj/H, Ratio E*/E2 Thenn de Barros Nijboer (II) 0,5 0,5 0,64 0,69 0,5 1,0 0,72 0,69 0,5 2,0 0,81 0,71 2,0 0,5 1,62 1,43 2,0 1,0 1,44 1,38 2,0 2,0 1,28 1,37 5. Application of the Equivalent Layer Theory in Concrete Pavement Design 5.1 General The equivalent layerthickness approaches, as described in the former paragraphs, are sometimes used in calculations for concrete pavements. Often, a relation is made with regard to the Herz/Westergaard model, in which the behaviour of the subgrade is decribed by one modulus of subgrade reaction (k) [10, 11] and the Hogg/Jeuffroy model [12, 13]. However, one forgets that the equivalent layer approaches are fundamentally based on a linear elastic multi-layer model (Boussinesq/Burmister), in which complete friction between the layers is assumed. In the Westergaard and Hogg models, the pavement is modelled as a thin (or thick) plate frictionless resting on the subgrade. Hogg also gives the formulas for a two layer model, in which complete friction between plate and subgrade is assumed, but this model is only used by Wiseman [14]. The subgrade in the (ordinary) Westergaard models, is modelled by a spring constant k (modulus of subgrade reaction). The deflection W is linear proportional to the vertical stress av. In the case of a point load, both the stress and the deflection will go to infinity at the surface in a halfspace. So, the modellation of the subgrade by a system of uncoupled/linked springs is mainly used for pavement layers (plates) resting on the subgrade. The subgrade is considered to be a dense liquid (Winkler foundation) and there is no horizontal stress transfer in the subgrade, in contrast with a linear elastic subgrade, in which horizontal load transfer occurs (Poisson ratio). In a linear elastic layer a deflection/deformation can occur, without the presence of a vertical stress. In spite of these differences the linear elastic theory of Burmister has been misused for the 'determination' of relationships between the modulus of subgrade reaction k and the Youngs' modulus Eo as will be shown. The starting point will be a halfspace loaded by a homogenous divided load P (constant contact stress cr0) on a circle with radius a. I I I I I I I I I I I I I I I I I I I I I • The definition for k is given by: 0 Ir = JS (43) I W Therefore k can be calculated from the equation belovv (assuming that the deflection below the loading plate is constant). k = Z— 01 by na2N0 . (44) I • I M P = 2 7t fovir}rdi o = 2nk(w{r)idr o Outside the area r=a the vertical stress at the surface is zero. So, in case of a Winkler foundation, no deflection may occur at the surface for r > a (in the case of a halfspace). However, in a linear elastic halfspace the surface deflection for r > a is not zero. In spite of the non-compatibility of these two types of foundations, relationships between k and Eo are often made. The k values based on these (wrong) relationships will be denoted by k*. Homogeneous Isotropic Linear Elastic Halfspace : I - Flexible loading (constant contact stress cr0 over the plate) -> The deflection under the loading plate is not constant! na2wc 8 • - Rigid loading (constant deflection Wo under the plate) " •j P = na2Wo 1 _ • • (45) In this formula the deflection Wc is the central deflection at r=0. ic# = • - v2) 2a(1 2E na (1 - v2) Wiseman uses a variant on the integral definition for the calculation of k given by: .ie» ^ P a r , ^ 27i / w{i)rdx l » r , , 2% (47 > \wii)xdr (outside the plate (r>a) the deflection is zero for a halfspace). The two first relationships depend on the radius a and so the calculated parameter kwill not be a constant, but will be inversely proportional to the applied radius of the I I I loading plate. Therefore, the two mechanical models (Burmister-Westergaard) are principally not compatible. The third definition has its advantage when it is used for two or more layer systems (see paragraph 5.2), in which the compressions of the pavement layers are negligible (plates). This definition gives a right k value, if the subgrade reacts as a Winkler foundation and the volume of the whole deflection bowl is measured. However, it should be noticed that for large distances from the center, the surface deflection will become negative [15] in the Westergaardmodel. I ™ Westergaard also presented a model in which a correction term was introduced in order to get a better match between theory and practice. In this so called Supplementary Theory, he presented an extreme case in which the subgrade reaction resembles that of a linear elastic subgrade (Hogg model; no friction). Westergaard himself stated that in that case the quantity k = k x L became a constant which implied that the original concept of k = constant has been leaved. — I Therefore the use of the linear elastic multi-layer model for the establishment of relationships k - Eo, can only be justified if the extreme case is used of the Supplementary Theory. It is recommended to use in that case a relationship between k x L = K - Eo rather than a relationship between k and Eo and also to calculate the stresses and strains using the Westergaard formulas, which are expressed in this quantity K. Besides the Westergaard model (Winkler foundation), the Pastemak model [16] can also be used. This model can be considered as a general Winkler foundation with a kind of horizontal coupling between the vertical springs k, denoted by a damping constant k*. The Westergaard model is a special case of the Pastemak model (k' = 0). | • * • « J m | • 8 • I 5.2 Two Layer Systems I In order to get still a relationship between k and Eo, one might calculate the vertical stress and deflection at the interface between the pavement layer and the subgrade in a linear elastic two layer system. The parameter k will be defined as the ratio of the calculated stress and deflection. However, the deflection in a point is not solely determined by the vertical stress, but this fact is completely disregarded. In a linear elastic two layer model (HvEvvvEo,v() the deflection W can be written as the sum of the deflection in the subgrade (Wo) and the compression of the pavement layer (AW,) -> W = Wo + AW, . If the equivalent layer theory of Odemark will be used, these two terms can be calculated with the aid of the formulas for a halfspace. AW1 = wiO,E1,v1} - W{Geq1,E1,v1) Wo = wiHeqlfEo,vJ wi th Geqx = 0 , 9 / ^ I " m • I (48) 1/3 and I * • • Heq1 = 0, 9 Hx x I I I I • • • Because on one hand there will be p.d. no compression of the pavement layer in the plate theory and on the other hand the compression for stiff pavement layers like concrete will be very small in the linear elastic theory it is justified to neglect the compression term AW V Therefore a definition for the analog parameter k* might be: # = | I f . I 1 1 1 • ° [1 + s i n ( P ) 3 x [l - sin 3 ((3) 3 a [1 + vo3 xcos(P) x [2 - 2v o + s i n ( p ) 3 (5°> tan<P, - The layerthickness H, will be nearly always bigger than the radius a :: (H, > a) and the same is even more true for the ratio of the moduli E, and Eo :: (E, » Ë^. These conditions imply that : COS(P) ° . .._. sin(P) = a HeQl (51) Heq, + a2 Before these approximations are used in the relationship above, the term (1 - sin3(/?)) has to be rewritten: [ l - s i n 3 ( P ) 3 = l - s i n ( p ) x [ l - c o s 2 ( P ) ] *cos 2 (P) • M9) Using the formulas for a halfspace with a circular load P (constant contact stress a as given in the appendix, thè following relationship will be found: ,:# _ E I oviHeqi,Eo,vo) M H E J (52) Using the approximations above the following relationship will be obtained: (1 + v o ) x (3 - 2v o ) xHeqx (53) '2 (1 - v2o) xHeqi » (3 - 2 v o ) If the Poisson ratio vo for the subgrade is taken equal to 0.5, the relationship established in this way will be: I I A wrong relationship between fc*(denoted as £**) and Eor is given by Eisenmann [17J: E = ° (55) Heq1 I Eisenmann deduced this relationship by combining the equivalent layerthickness theory of Odemark, together with the Boussinesq formulas for a point load P on a halfspace [18]. For positions on the centerline the vertical stress <7v(z) and deflection W(z) are given by: for v ° = 0 f 5 ifc follows _ o iz) • _ E —ïj—r- = —WXZJ z For the depth z Eisenmann uses the equivalent layerthickness Heq, for which in reference [18] the Poisson ratio v, is taken equal to 0,5. In later reference [17] Eisenmann uses a more comon value of 1/6 « 0,16 for the Poisson ratio v, of the (concrete) pavement layer, which leads to: 1/3 *1 Heqx = 0 , 825 (57) The large difference between the two relationships (^-E o and fc*#-E„) can be traced back to the use of a pointload P instead of a circular load P. Pointload P Circular 3P _ 2 load P (z>a) P_ 71 Z2 2 71 Z vo) x ( 3 - 2 v o ) 2 P ( 1 - v2o) (58) 2nzEo 2 n zE o Afterwards calculated fc*# value is used in the Westergaard model. In the Westergaard model a circular load is applied. So, only for this reason alone the procedure as suggested by Eisenmann ought to be dissuaded. Another objection against the procedures above is the fact that, in the Westergaard model, the pavement layer (plate) rests frictionless on the foundation i'.c. the subgrade which is not the case in the linear elastic models used in the equivalent layer theory. | • • I I I If as a starting point the Hogg model without friction between plate and the linear elastic subgrade (E^vJ is chosen, an impression of the influence of this ignoration can be obtained. According to Eisenmann [18] at the interface due to a point load the vertical stress av {H,} is given by: x 6EO(1 - y\) Ex(l 1/3 (59) -" v|) Using Odemark's expression for the equivalent layerthickness Heq, this equation can be rewritten as: ov = 1,078 x (60) 2 it Heql In the same way the deflection W {0} = W {Hl} can be obtained: W{0} = ^& x SD = 3,955x 2 D ( 1 - vl) 2/3 PU - v|) (61) 2 7i EoHeqx + v o ) x (3 - 2v o ) 2nEoHeqx 3 , 9 5 5 x (1 - v o ) (3 - 2v o ) The relationship for /c*#f is given by: 3,955 'o x = Cx (62) The value of the coëfficiënt C in the relationship above is given in Table 7 as a function of the Poisson ratio vn. Table 7 C = Coeffiient C as a function of the Poisson ratio 0,15 0,25 0,35 0,45 0,50 I 0,836 0,872 0,932 1,025 1,090 f It should be remarked that both the stress av as well as the deflection W in the Hogg model (no-friction), do not differ much from the stress av and deflection W in a linear elastic two layer model (Burmister) with complete friction between the two layers and using a point load. The ratio WHogg/WBurmister = [3,955 x(1-iO3'/K3-2«O] varies from 1,245 at vo = 0,15 to 0,989 at vo = 0,50. The influence of friction or no-friction between the layers in the Hogg model can be obtained from the two following expressions for the deflection: Hogg;no-friction *Hogg; friction o = O, 1 9 2 5 X (63) 2 (1 - Vl) W(no-friction) W {friction) and ( 1 + v j x (3 - 4 v J _ (1 + vo) x (3 - 4 v 4X (1 - vl) The ratio of the defiections Kv,) are presented in table 8 as a function of the Poisson ratio vo. Table 8 1/ "0 The ratio of the defiections in the two Hogg models (friction versus nofriction) as a function of the Poisson ratio = f(O = 0,15 0,25 0,35 0,45 0,50 0,793 0,763 0,723 0,668 0,630 As already stated above, it is quite curious that a remarkable resemblance exists between the stresses, strains and defiections in the Hogg model, with no-friction between the layers and a linear elastic two layer model with friction between the layers, when in both models a pointload P is applied. From the exercise above, it can be concluded that the (original) Westergaard model is not compatible with a Hogg model or a linear elastic Burmister model. Only in the extreme case of the Supplementary Theory of Westergaard [11], which resembles the Hogg model (no-friction) [12], a relationship between k and Eo might be established. However, in that case, the parameter k has no real physical meaning and may only be considered as a mathematical quantity. By equalization of the defiections in both models the following relationship can be obtained: ie = 0 , 1 2 8 x Heq1 for v o = 0 , 5 i t follows: (64) k = 0,17x Heg1 If the relationship above is used the extreme case of the Supplementary Theory of Westergaard has to be used. I I I I I I I I I I I I I I I I I I I I 5.3 Th ree Layer Models A three layer model can be simplified into a two layer model by either: A: Combination of the two pavement layers into one equivalent layer (Concrete layer on a foundation of lean concrete) or by: To consider the second pavement layer as an increase of the bearing B: capacity of the subgrade (Concrete on a foundation of unbound material) Possibility A: If the two pavement layers are taken together in order to calculate the stresses and strains at the bottom of the second pavement layer the following procedures can be used: A1) Equivalent layer theory of Barber and Palmer (chapter 1) A2) Equivalent layer theory of Odemark (chapter 2) A3) Equivalent layer theory of Ullidtz (chapter 3) A4) Equivalent layer theory of Nijboer and Pronk (chapter 4) Possibility B: If the second pavement layer is considered as an increase in the bearing capacity of the subgrade the increase can be obtained using: B1) I I I I I I I I I Equivalent layer theory of Odemark (chapter 2) Needless to say that, in all the theories mentioned above complete frictions between all layers is assumed. For systems where no-friction between the first and second pavement layer is assumed (Jeuffroy/Bachelez model), no equivalent layer theories are found in the literature. If complete friction between all layers is assumed, calculations at the top of the subgrade can be performed, using the (total) equivalent layerthickness for layer 1 and layer 2 which is given in the case yo-v^=v2 by: Heq2 = 0 , 9 x H, x + H2x\ (65) 1/3 1/3 = 0,9x X H2 Si • E o E2 • E2 The term H1.{(E1.EO)/(E2.E*2)}1/3 can be considered as an effective increase of the second pavement layer. The expression above can also be written as: .1/3 Heq2 = 0,9x Er (66) 1/3 H2 xE2 + H±xE- Besides the factor Eo/E*2 this procedure can be seen as the summation of the terms H,E1/3 If there is no friction (complete slip) between layer 1 and 2 in analogy with the expression above one might define an equivalent layer as: 1/3 (67) Also in the case of three layer systems, the equivalent layer theory is misused to obtain a /f* - Eo relationship based on the ratio <7V/W. Eisenmann [17] gives the following relationship, if the compression of the second pavement layer, may not be neglected. Heq1 x Heqt Heq2 x( ' Heqt Heq± (68) with : Heqt = Heq1 + Heq2 1/3 Heq1 = 1/3 ; Heq2 = 0,9 xH2x E, The starting point for this relationship was a pointload and the second pavement layer was considered as an increase of the subgrade properties. Furthermore, the influences of the Poissonratios were ignored by taking vo - v, = v2 = 0,5. In view of all these neglections the following more 'fundamentaP approach is advised if one would establish a relationship for /f*: v Wo For the calculations of HeqvHeq2, and Geq2 the Odemark's formulas are used. (69) References 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. Palmer, L.A. and Barber, E.S., "Soil displacement under a circular loaded area" Proceedings Highway Research Board, Vol, 20, 1940 Odemark, N., "Investigations to the Elastic Properties of the Soils and the Design of Pavements to the Theory of Elasticity", Stockholm, 1949 Ullidtz, P., "En Studie af to Dybdeasfaltbeasestelser", Ph.D. dissertion, The Technical University of Denmark, Lyngby, 1973 Ullidtz, P., "Some Simple Methods Determining the Critical Strains in Road Structures", Technical University of Denmark, Department for Road construction, Transportation and Town Planning, Lyngby, 1974 l Pronk, A.C., "Berekening van Deflecties, Rekken en Spanningen in een Half Oneindig Medium en in 2- of Meerlagen Systemen", Report DIM-R-80-75, Road and Hydraulic Engineering Division, The Netherlands, 1980. Poulos, H.G. and Davis, E.H., "Elastic Solutions for Soil and Rock Mechanics", John Wiley & Sons, Inc., New York-London-Sydney-Toronto, 1974 Nijboer, L.W., "Dynamic Investigations of Road Constructions", Shell Bitumen Monograph no. 2. Pronk, A.C., "Equivalente laagdikte theorieën", Internal Report, TW-N-86-42, Road and Hydraulic Engineering Division, The Netherlands, 1986. Thenn de Barros, S.r "Deflection Factor Charts for Two- and Three-layer Elastic Systems", Highway Research Record, No.145, 1966 Westergaard, H.M., "Stresses in Concrete Pavements Computed by Theoretical Analysis", Public Roads7, no. 2, 1926 Westergaard, H.M.H., "Analytical Tools for Judging Results of Structural Tests of Concrete Pavements", Public Roads 14, no. 10, 1933 Hogg, A.H.A., "Equilibrium of a Thin Plate, symetrically loaded, resting on an Elastic Foundation of Infinite Depth", London, Edinberg and Dublin, Philosophical Magazine and Journal of Science Vol. 25, (168), 1938. Jeuffroy, G., "Note sur Ie Comportement des Chaussées", Annales des Ponts et Chaussées, no. 3, 1959 14. 15. 16. 17. 18. Wiseman, G., "Flexible Pavement Evaluation Using Herz Theory", Transportation Engineering Journal, ASCE, Vol. 99, No.TE3, Proceeding Paper 9921, 1973. Pronk, A.C., "A Further Investigation into the Westergaard Formulas for Interior Loading and the (Mis)use of the Modulus of the Subgrade Reaction k", 2nd International Workshop on Theoretical Design of Concrete Pavements, 1990, Siguenzq, Spain. Pronk, A.C., "The Pastemak Foundation - An Attractive Altemative for the Winkler Foundation", Proceedings of the 5th International Conference on Concrete Pavement Design and Rehabilitation, Purdue University, West Lafayette, Indiana, 1993. Eisenmann, J., "Betonfahrbanen", Verlag von Wilhelm Ernst & Sohn, Berlin-MunchenDusseldorf, 1979. Meier, H., Eisenmann, J. and Koroneos, E., "Beanspruchung der Strasse unter Verkehrslast", Forschunsarbeiten aus dem Strassenwesen, Neue Folge, Heft 76, Kirschbaum Verlag, Bad Godesberg. APPENDIX: 1) The The The The Overview of Formulas Formulas for the deflections, stresses and strains in a halfspace due to a point load P (Boussinesq). horizontal distance to the centerline is denoted by: r [m] depth below the point load is denoted by: z [m] distance to the point load is denoted by: R [m] angle between the centerline and the position r,z: a [°] (70) R = Jz2 + sin(a) = cos(a) = R (71) _z ~R 2 W = P ( 1 + v) x 2 ( 1 - v ) + c o s ( « ) 2nRE (72) I I I I I I I I I I I I I I I I I I I I I a = 3xPxcos3(g) 2uR2 I ( 3xcos(a) xsin2(a) | I px (73) 1 2 ] + C O S (Ot) 2 2nR ev = ^ x ( l + v) x 3cos 3 (g) - 2vxcos(g) „ 6) I e • r - - x ( l + v ) x ' • • 3 - 3 x c o s ( a ) + (3 - 2 v ) x c o s ( a ) _ • (77) 1 2v ~ , .. 1 + c o s (a) - cos(a) + - — For point on the centerline these formulas can be simplified to: (79) I I E 27i z or{z) = oc(z) = - Px ° I I (1 ' 2 v ) 4 11 Z 2 (H) (1 - v) xor(z} - er{zl = zt\ 2) \xoviz} (83) Formulas for the deflection in a halfspace due to a homogenous loaded circular plate. xPx ƒ _1_ 2x(l-v)x[T 2 -r i ]-C 2 x T dy X \ \ 1- [ a i y] 2 _ i-[^y] 2 \ \ a i (84) v/l-- y 2 • r2 z2 2 y + + -JT2 r On the surface (z=0) this equation reduces to a Hypergeometric serie. (85) (2n) ! T a r (n\)2 3) n +1 Formulas for the deflections, stresses and strains in a halfspace on the centerline under a circular homogenous divided load P. tan(p) = — a (86) - 2v) 1 + Wiz) =P x ( l + v)xcos(P)x n aE w{0) - 2P (88) - sin3(p)] n (89) I I I I I I I I I I I I I I I I I I I I P ariz} = o t{z) = na2 x (90) ±-2ï. - (1 + v) xsin(p) + -ixsin 3 (p)l u £i J _ -PU + V) w na2E [ l - s i n 3 ( P ) oi) +2xvx[sin(P) - 1 ] ] {zY= x na2E (92) I I f z{z) B AMz) = 4) • I I = — L 2^ x s i n ( p ) x c o s 2 ( p ) 2na E 3 P 2 4na E i f v = 0 , 5 os) xsin(p) xcos2(p) ; (v-0.5)(94) The requirement that there will be no volume change (ev + 2 er = 0 ) for v = 0,5 is fulfilled. The shear stress r2 equals in that case - e r x E If the depth z is smaller than two times the radius of the loading plate (z < 2.a) the compression of the material between 0 and z can be very well estimated by: " I 3 - (1 - v) xsin(P) + - | x s i n ( p ) | er{z) = et{z}= • — "2 2XV lt should be noticed explicitely that in these equations compression (negative strain) has a positive sign (because the vertical loading is taken positive). If the Poisson ratio v equals 0,5 the formulas for the strains can even more simplified: I I 1 2 2na E x[2x(l -v 2 ) - 0 , 7 ] os) Formulas for the deflections, stresses and strains in a halfspace on the centerline under a circular rigid load P. If the loading plate is stiff (rigid) instead of a flexible plate (constant contact stress cr=P//7.a2) the following formulas have to be applied. p Contactstress oo = — — • for r < a (96) 2 2 2naja -r 1 + 3x (97) = Px W{Z) . (98) (1 - v) x7C - 2 x a r c t a n ( —) + ax 1 + *M = 1 - 2 (1 - v) 1 + 5) 1 + z2 a2 z2' a2 2 (99) Formulas for the deflection in a stress dependent halfspace due to a homogenous loaded plate If the modulus of the halfspace is stress dependent according to: Material modulus E = Cx with o± = major principal n stress (100) os = referencestress n = coëfficiënt C = constant than the surface modulus calculated from deflection measurements at r = 0 is given by: n = (1-272) xCx (101) I I I I I I I I I I I I I I I I I I I I For r = O the following equation yields: n 2 a = (1-2x2) xCx ofl-D2) (102) The major principal stress can be approximated by: aAzY= 1 3 a2 * r n 2 The surface modulus for r o 0 is given by: n I I I I I I I I I I I It should be remarked that by application of these equations in the equivalent layer theory the surface modulus Es must be used for the calculation of the defiections and that the modulus-E-must be used for the calculation of strain and stresses.
© Copyright 2024 ExpyDoc