Rheological Interconversions
|
|
- Caitlin Gibbs
- 5 years ago
- Views:
Transcription
1 Technical White Paper Rheological Interconversions Fundamental Properties of Asphalts and Modified Asphalts III Product: FP 23 March 215 Prepared for Federal Highway Administration Contract No. DTFH61-7-D-5 By Ronald R. Glaser, Thomas F. Turner, Changping Sui, Stephen Salmans and Jean-Pascal Planche Western Research Institute 3474 North 3rd Street Laramie, WY
2 TABLE OF CONTENTS INTRODUCTION...1 Shifting...2 Prony Series Shear Relaxation from Complex Modulus...4 Shear Relaxation to Extensional Relaxation...1 Creep Compliance Deconvolution and Stiffness...12 Prony Series Distribution Improvement...12 EXPERIMENTAL...13 RESULTS AND DISCUSSION...13 Nifty-Shifty Results...14 Prony Herd Calculation Results...21 CONCLUSIONS...37 RECOMMENDATIONS...37 ACKNOWLEDGMENTS...38 DISCLAIMER...38 REFERENCES...38 i
3 LIST OF FIGURES Figure 1. Graph. Least squares and fractional exponent residual fits compared...3 Figure 2. Diagram. Alignment shifting method minimizing connecting line length...4 Figure 3. Graph. AZ1-3 complex modulus master curve reference temperature -12 C...14 Figure 4. Graph. AZ1-3 Arrhenius shift function...14 Figure 5. Graph. AZ1-3 WLF shift function...15 Figure 6. Graph. AZ1-3 polynomial Arrhenius shift function...15 Figure 7. Graph. AZ1-3 black space...16 Figure 8. Graph. MN1-5 complex modulus master curve reference temperature -18 C...16 Figure 9. Graph. MN1-5 Arrhenius shift function...17 Figure 1. Graph. MN1-5 WLF shift function...17 Figure 11. Graph. MN1-5 polynomial Arrhenius shift function...18 Figure 12. Graph. MN1-5 black space...18 Figure 13. Graph. YNP complex modulus master curve reference temperature -24 C...19 Figure 14. Graph. YNP Arrhenius shift function...19 Figure 15. Graph. YNP WLF shift function...2 Figure 16. Graph. YNP polynomial Arrhenius shift function...2 Figure 17. Graph. YNP black space...21 Figure 18. Graph. AZ1-3. Prony spectrum...22 Figure 19. Graph. AZ1-3. Measured and calculated phase angle...22 Figure 2. Graph AZ1-3. Complex modulus measured compared to complex modulus calculated from relaxation function convolution...23 Figure 21. Graph. AZ1-3. Shear relaxation from Prony series...23 ii
4 LIST OF FIGURES (continued) Figure 22. Graph. AZ1-3. Extensional relaxation from shear relaxation assuming.5 Poisson s ratio and time dependant Poisson s ratio...24 Figure 23. Graph. AZ1-3. Creep compliance from deconvolution calculation and reciprocal extensional relaxation approximation...24 Figure 24. Graph. AZ1-3. Stiffness from calculated creep compliance...25 Figure 25. Graph. MN1-5. Prony spectrum...25 Figure 26. Graph MN1-5. Measured and calculated phase angle...26 Figure 27. Graph MN1-5. Complex modulus measured compared to complex modulus calculated from relaxation function convolution...26 Figure 28. Graph. MN1-5. Shear relaxation from Prony series...27 Figure 29. Graph. MN1-5. Extensional relaxation from shear relaxation assuming.5 Poisson s ratio and time dependant Poisson s ratio...27 Figure 3. Graph. MN1-5. Creep compliance from deconvolution calculation and reciprocal extensional relaxation approximation...28 Figure 31. Graph. MN1-5. Stiffness from calculated creep compliance...28 Figure 32. Graph. YNP. Prony spectrum...29 Figure 33. Graph YNP. Measured and calculated phase angle...29 Figure 34. Graph. YNP. Complex modulus measured compared to complex modulus calculated from relaxation function convolution...3 Figure 35. Graph. YNP. Shear relaxation from Prony series...3 Figure 36. Graph. YNP. Extensional relaxation from shear relaxation assuming.5 Poisson s ratio and time dependant Poisson s ratio...31 Figure 37. Graph. YNP. Creep compliance from deconvolution calculation and reciprocal extensional relaxation approximation...31 Figure 38. Graph. YNP. Stiffness from calculated creep compliance...32 iii
5 LIST OF FIGURES (continued) Figure 39. Graph. DSR and BBR stiffness compared using a constant Poisson s ratio assumption in the interconversions...33 Figure 4. Graph. The stiffness curves generated using a range of temperatures in the DSR...34 Figure 41. Graph. The stiffness curves generated using a range of constant Poisson ratio value assumptions...35 Figure 42. Graph. The stiffness curves generated using a range of constant Poisson ratio value assumptions and time dependent Poisson s ratio...36 Figure 43. Graph. Poisson s ratio decay times as a function of temperature...37 iv
6 LIST OF TABLES Table 1. Poisson s ratio exponential decay parameters for three binders...36 v
7 RHEOLOGICAL INTERCONVERSIONS INTRODUCTION The topic of interconversion of the various rheological moduli functions has been of interest to rheologists for a considerable time. The utility of being able to measure a modulus function and then generate a more useful one provides a great deal of latitude in selecting measurement methods. In the area of asphalt binder rheology, measurements are taken at high and modest temperatures using the dynamic shear rheometer (DSR). Until recently, low temperature data was unobtainable with this method, and creep measurements of a beam using the bending beam rheometer (BBR) were employed to obtain low temperature measurements needed for performance grading (PG). This combination has been in routine use since it was developed under the Strategic Highway Research Program (SHRP), but is limited for some applications due to the large amount of material required to fabricate the binder beams. WRI developed the 4-mm DSR plate method with machine compliance corrections that allows the use of DSR on very small sample down to -4 C (Sui et al. 21). There are other advantages to using DSR in lieu of BBR. Input data for mechanical performance simulations generally require the time domain relaxation function, often expressed as a Prony series in material files, instead of dynamic modulus. Conversion of dynamic modulus data to relaxation is a one step process, while conversion from bending beam creep data involves several steps, and has been shown to be inconsistent with DSR testing and conversion on the same materials (Marasteanu and Anderson 1999). The mathematical manipulations required to convert beam creep modulus data to relaxation are quite involved, and approximate methods are often applied instead of the mathematically rigorous, but tedious methods. These are covered quite well in the Marasteanu (1999) thesis. We have developed a series of algorithms to make the conversions while adhering to fundamental solutions based on linear viscoelastic theory as much as possible. This effort has attempted, in part, to understand the discrepancies between relaxation as determined from BBR measurements and relaxation as determined from DSR measurements. The objectives were to produce more general and fundamental interconversions that would apply to highly modified asphalt binders, and to determine if the use of approximate methods may have played a role in discrepancies noticed in the literature. The computing power available in the average personal computer is such that more rigorous numerical methods are indeed practical. There is, however, always a concern over computation error propagation with some of these numerical methods. To compare DSR and BBR measurements, DSR data are converted to stiffness, and then compared to stiffness from BBR measurements. The computational sequence is: 1) Shift DSR isotherm data to a time-temperature superposed (TTS) master curve. 2) Obtain a Prony series representation directly from the complex modulus using an algorithm developed in this study that removes any error caused by erroneous phase angle measurements. 1
8 3) Convert shear relaxation to extensional relaxation (assuming flexural and extensional moduli are equivalent). These computations require the use of Poisson s ratio, and we employed a constant Poisson s ratio and two approaches to time dependant Poisson s ratio methods in the calculations. 4) Numerical time domain deconvolution was employed to convert the extensional modulus to creep compliance. 5) Creep compliance was inverted to get stiffness. 6) Compare the interconverted DSR stiffness to the stiffness obtained from BBR creep measurements. All of these computations are performed on two separate Microsoft Excel spreadsheets. The first of these is Nifty Shifty, which performs a wild shift according to a modulus sorted connecting line minimization routine The Prony Herd spreadsheet accepts the shifted master curve data and performs all of the interconversion steps described in this study. The Prony series computation was later modified by changing the form of the Maxwell element coefficient distribution. Shifting The SHIFT routine first sorts the data according to complex modulus. The data are then classified according to the isotherm temperature and then shifted using an alignment algorithm. This spreadsheet requires that all isotherms overlap to some extent. If you imagine lines connecting the data points, the unshifted data has many zigzag patterns in it. The alignment is accomplished by minimizing the relative dimensionless length of the connecting lines using a sum of fractional exponent residuals. This is a much better fitting criterion than the sum of least squares. The square of the residuals weights the aberrant points more heavily than the good points. The traditional equation for a dimensionless sum of squares residue is: n 1 i= 2 (( datai fiti) / datai) (1) (Note this form evaluates local relative error and works over any range of data, including the many decades encountered in rheological data.) While a better approach is the fractional exponent sum of residuals: n 1 i=.25 ( ABS(( datai fiti) / datai)) (2) The two methods are compared side by side in the following figure: 2
9 3 Dependent variable data relative fractional exponent (.25) standard least squares Independent variable Figure 1. Graph. Least squares and fractional exponent residual fits compared. The shift factors are adjusted for each point according to its isotherm temperature until the sum of the residuals (overlap sums) is minimized using Excel Solver. Good results are obtained with an exponent of.25. In usual math notation, the code line above is: Residual at a data point = ( (1 ω / ω ) ) 2 + (1 G* / G* ) 2 + (1 ω / ω ) 2 + (1 G* / G* ) n (3) where: ω = shifted frequency G* = complex modulus n = weighting exponent < 1 The shifted frequency is then computed. The following diagram (figure 2) geometrically illustrates the line length minimization principle: 3
10 Figure 2. Diagram. Alignment shifting method minimizing connecting line length. Minimization is accomplished using Solver programmatically by minimizing the residuals sums by adjusting the frequency. The shift is made to the warmest isotherm one isotherm at a time until all isotherms have been shifted. Once all the shift factors are computed, the shift functions are fit (the linear Arrhenius using standard least squares, and the WLF and quadratic Arrhenius using solver). Additional lines of code in the SHIFT subroutine manage the plots on the spreadsheet. Once the shift is accomplished and the various shift functions determined, the user can generate a new master curve at any desired temperature by entering the values on the spreadsheet. The curve is shifted according to the chosen shift function and temperature. The new temperature only changes the frequency column, and new reference temperature columns are provided on the spreadsheet. The glass transition temperature in this fit is not fixed, but is fit along with the other parameters. The master curves are displayed on a plot that also contains the original isotherms so the quality of the shift functions can be easily evaluated visually. A good shift function should place the shifted back data through the isotherm generated at the same temperature. Prony Series Shear Relaxation from Complex Modulus We start with the defining equation, the hereditary integral for visco-elastic materials: for a Maxwell Element the relaxation modulus is t g() t = ge β (4) r r where t = time (s) and β = decay constant (s) 4
11 The response of this element in a DSR is given by the hereditary integral: t ( t τ ) β dε (sin( ωτ ) σω () t = ge r dτ (5) dτ where: σ ω = stress (Pa) ω = angular frequency (rad/s) τ = variable of integration (s) ε = strain (m/m) The integral can now be evaluated: () t β 1 e ( ωsin( ωt) cos( ωt)) σ + ω β β () t = grω ε 2 1 ω + 2 β (6) Notice the exponential transient term. From a practical stand point, rheometric data should not be collected until the response stabilizes. Some preliminary plots indicate this happens quite quickly, usually in less than 1 cycle for a single Maxwell element. A rapid approach to steady state is not necessarily true for a series with long decay times, which can be misinterpreted as the equilibrium modulus for a visco-elastic solid. For our current inquiry, we are interested in steady state response, so we take the limit as t and luckily find that the exponential vanishes. 1 ( ωsin( ωt) + cos( ωt)) σ ω β () t = grω ε 2 1 ω + 2 β (7) Define A and B as 1 β B = 2 1 ω + 2 β (8) 5
12 ω A = 2 1 ω + 2 β (9) After a little algebra (surprisingly simple) one obtains, 1 1 ϕ = tan ωβ (1) So, the steady state response function is ( t) = grω A + B sin ωt+ tan ε σ ω ωβ (11) The full transient form is ( t) g A B sin ε t tan σ ω = rω + ω + ωβ () t β e β 1 ω 2 β (12) Considering the series form with several Maxwell elements, the integrals simply add. Now that we have the response for a single Maxwell element, most complex materials can be treated by describing the relaxation function as series of Maxwell Elements: n i g () t = g e β (13) r i= 1 ri, t The complex vector form of the response is most convenient to use, as we shall see when we calculate the phase angle of the series: 6
13 1 ( ωsin( ωt) + cos( ωt)) σ ω β () t = grω ε 2 1 ω + 2 β (14) g ω r 2 grω β = (sin( ωt)) + (cos( ωt)) ω + ω β β (15) where: [ A(sin( ωt)) ] [ B(cos( ωt)) ] = + (16) g r ω 2 g A r ω β = and B = ω + 2 β ω + 2 β So the series can be written as: n n σ ω = i ω + ε i= 1 i= 1 ( t) A(sin( t)) Bi(cos( ωt)) (17) Let A Σ n = A and i= 1 i B Σ n = Bi i= 1 then ϕ B 1 Σ = tan AΣ (18) In the complex plane A Σ represents the real part, and B Σ represents the imaginary part, and one can work this out with phasors. The next concern in the computation is getting the modulus, which occurs when the response is at its maximum. The maximum strain is at π/2. The stress peaks before the strain, so the proper angle at the peak is π/2 phase angle. Divide this by the frequency, ω, to obtain the time at peak stress: max σ ω G * ( ω) = = AΣ(sin( π / 2 ϕ)) + BΣ(cos( π / 2 ϕ)) ε (19) 7
14 So, the algorithm is fairly straight forward. For each frequency ω we calculate complex modulus G* and phase angle, ϕ, like this: Given the prony series constants, first sum the A and B terms. A Σ n = A, B i= 1 where: 2 giω Ai = 2 1 ω + 2 βi giω βi Bi = 2 1 ω + 2 βi Now get the phase angle: i Σ n = B (2, 21) i= 1 i ϕ B 1 Σ = tan AΣ (22) Then get the modulus: or max * σ ω ( ) G ω = = AΣ(sin( π / 2 ϕ)) + BΣ(cos( π / 2 ϕ)) ε σ 2 2 * = ω = Σ + Σ sin / 2 ε max G A B ( π ) (23) (24) σ ω G* = = A + B ε max 2 2 Σ Σ (25) σ ω 2 2 G* = = G' + G" ε max (26) Use the working equation above to fit the Prony coefficients using generalized gradient reduction (Excel Solver ). The Prony series to G(t) conversion is a straight-forward sum of exponentials. 8
15 The Prony Herd sheet constrains the coefficients to fit up to 3 log normal Gaussian distributions (Jongepier and Kuilman 1969). where: g j (27) t 3 6 βk, j Ecalc = gk, je prony series 3 distributions k= 1 j= 1 2 log1( βk) log1( centerk) log1( σ k ) 1 Gglassy ( relative H k ) e 2 log1( σ k ) = normalized distributed coefficients 3 6 gk, j k= 1 j= 1 The program solves for: relative H,log1( center ),log1( σ ), G k k k glassy by minimizing: n * * G data ( tk ) G calc ( tk ) errorfunc = G ( t ) * k = 1 data k 3 Using 3 log-normal distributions to describe the Prony series exponential decay coefficients works well for creating a relaxation function that re-convolves to the master-curve, and this approach was used for the comparison with BBR stiffness. Further studies with alternative distributions revealed that a generalized Gaussian distribution (the exponent now a variable instead of 2) would provide high quality relaxation computations for unmodified binders with only one distribution (3 adjustable parameters). The distribution is further constrained by setting the peak of the distribution at the decay time for the most rapidly decaying Maxwell element. The reduction in adjustable parameters results in a unique solution using non-linear fitting methods. g j 1 Gglassy ( relative H k ) e 2 log1( σ k ) = 3 6 gk, j k= 1 j= 1 β log1( βk) log1( centerk) log1( σ k ) (28) It is important to note that using non-linear regression methods, the phase angle is not required to obtain the shear relaxation modulus, and, can in fact be derived from the complex modulus alone. This fact is quite helpful when working with data containing poorly measured phase angles, a frequent occurrence with DSR measurements under low torque conditions. Once the Prony series distributions are fitted, and the Prony series spectrum known, the shear relaxation function is known and an x,y pair list can be generated to describe the material relaxation in time. 9
16 Shear Relaxation to Extensional Relaxation An approximate formula is presented by Lakes (1992) for the time dependency of Poisson s ratio, assuming the bulk modulus is constant: 1 Et () υ () t = (29) 2 6B where: υ(t) = time dependent Poisson s ratio E(t) = extensional modulus (relaxation) B = bulk modulus (assumed time invariant) So, the equation above combined with Et () = Gt ()2(1 + υ) (3) can be rearranged with a bit of algebra to produce: 3 Gt () Et () = 3 Gt () 1+ 6B (31) The above formulation does not work when comparing BBR data. In the BBR flexural setting, it appears that a Poisson s ratio model based upon the flow of a mobile phase does indeed reconcile the DSR and BBR data sets. An approximate formulation follows assuming Darcy flow: dv dt = KA dp η dh (32) where: V = volume t = time K = permeability A = area normal to flow η = viscosity dp / dh = pore pressure (stress) gradient and dlwh KLw dp dt = η dh (33) 1
17 where: w = width h = height of some small slice L = length Assuming L, h, K, η, and dp/dh don t change much or dw K dp = dt (34) w hη dh hη dh λ = (35) K dp η dh λ = (36) ' K dp where ' K = relative permeability, so w w dw = w t 1 dt λ (37) t wt () = we λ (38) ν = w/ L= Poisson's ratio (39) so, approximately or wt () we = L L t e λ t λ ν() t = ν (41) (4) which is about right if v() t = as t. More generally ( ν ) λ ν() t = v e + v (42) t 11
18 The source code can be easily modified to use other Poisson ratio models, including the incorrect constant assumption frequently employed by others. Creep Compliance Deconvolution and Stiffness From Ferry (198) the relationship between creep and extensional relaxation is given by: t Et ( τ ) D( τ ) d τ = t (43) A number of methods exist for deconvolving this equation to D(creep) based upon the idea that for small steps D can be considered a constant (or ramp) and taken outside the integral in small Euler or trapezoidal rule steps. E(t) is the extensional relaxation. A variety of combinations of trapezoidal or Euler rule to either kernel or both have been employed, but we have found that the double Euler is least likely to oscillate (Park and Kim 1999). The discrete form reduces to the following recursive formula. Note that the double trapezoidal rule form for t 1 is employed as the equal spaced form is very sensitive to this first step. Dt ( ) = 1/ Et ( ) (44) 3 Et ( 1)/ Et ( ) Dt ( 1) = Et ( ) + Et ( ) 1 (45) n 1 tn Dt ( i 1) Et ( n ti)( ti ti 1) 1 Dt ( n) = ( Et ( ) + Et ( t ))( t t ) n n 1 n n 1 (46) for equal spacing n 1 tn Dt ( i 1) Et ( n ti) t 1 Dt ( n) = ( Et ( ) + Et ( )) t 1 (47) The recursion starts at index 2, so the last equation is looped n=2 to n maximum to get all the point values for D. Stiffness is the inverse of creep compliance. Prony Series Distribution Improvement A second series of studies was performed in an effort to improve the relaxation computation. We replaced the Gaussian distribution of the Prony spectrum with a stretched exponential, improving fit quality while reducing the adjustable parameter count. The ill-posed nature of this 12
19 computation appears to be eliminated. Additional investigations of single isotherm relaxation computations were done to test the robustness of the method and to investigate the deviation of asphalt binders from ideal thermally simple behavior. EXPERIMENTAL This report consists of three separate studies. The first investigates the interconversion of low temperature DSR complex modulus data to stiffness, and then compares that to stiffness results from the BBR. The second study focuses on the relaxation computation alone using and improved Prony spectrum distribution shape. The third study investigates the limitations of timetemperature superposition and the use of relaxation modulus conversions of individual complex modulus isotherms. The data from the examination of 3 asphalts using the bending beam rheometer and the 4 mm plate dynamic shear rheometer were inter-converted using the methods described in the introduction and compared in this study. A detailed description of these measurements can be found in Sui et al. (21). These binders are MN1-5 from the Minnesota comparative pavement performance site, AZ1-1 from the Arizona comparative pavement performance site, and YNP from the Yellowstone comparative pavement performance site. These data were generated in support of the 4 mm plate DSR method, described in a topical report (Farrar et al. 215). RESULTS AND DISCUSSION To examine inter-conversion results for unmodified asphalts, three binders were examined where data already existed from 4 mm plate DSR testing and BBR runs on the same materials. These binders are all unmodified and were used in field performance test sites as well. The Minnesota test site binder, MN1-5, the Arizona test site binder, AZ1-3, and the Yellowstone test site binder, YNP were examined. The BBR testing on the MN1-5 binder was done at -18 C, the AZ1-3 binder was run at -12 C, and the YNP binder at -24 C. The first step in the interconversion process is to convert the DSR isotherms to a master curve shifted to the BBR test temperature. Three shift functions were then fit to the alignment based master curve. The standard Arrhenius plot was used to obtain a simple Arrhenius shift function, the WLF equation was fit without specifying the glass transition temperature. While the use of a known glass transition temperature is often employed, for asphalts this produces poor fits at low temperatures, and a curve fitted reference temperature (usually lower than the true glass transition temperature) provides excellent shift function fits. The free volume theory for the WLF does not appear to work very well for producing a precise shift function. An alternative to WLF is also used. This simple alternative is a quadratic fit on an Arrhenius plot, which essentially describes temperature dependant activation energy for viscous flow. All three methods work reasonably well. Figures 3-17 show the shifting results for these three binders. 13
20 Nifty-Shifty Results 1.E+1 shift G* Complex Modulus, Pa 1.E+9 1.E+8 1.E+7 1.E+6 1.E+5 1.E-9 1.E-4 1.E+1 1.E+6 1.E+11 Reduced time, s shift Tmax Data Shift New Arrhenius Shift New Poly Shift New WLF Figure 3. Graph. AZ1-3 complex modulus master curve reference temperature -12 C. Arrhenius Log at y = E+4x E-1 R² = E Series1 Linear (Series1) 1/(T )-1/(Tref ) Figure 4. Graph. AZ1-3 Arrhenius shift function. 14
21 WLF Log at data WLF T-Tref Figure 5. Graph. AZ1-3 WLF shift function. 3 POLY Log at data fit Poly. (data) y = 1.241E+6x E+4x E-2 R² = 9.996E-1-1.5E-3-1.E-3-5.E-4.E+ 5.E-4 1.E-3 1.5E-3 2.E-3 1/(T )-1/(Tref ) Figure 6. Graph. AZ1-3 polynomial Arrhenius shift function. 15
22 Black Space 6 5 Phase angle, degrees E+5 1.E+6 1.E+7 1.E+8 1.E+9 1.E+1 Complex Modulus, Pa Figure 7. Graph. AZ1-3 black space. 1.E+1 shift G* Complex Modulus, Pa 1.E+9 1.E+8 1.E+7 1.E+6 1.E+5 1.E+4 1.E-9 1.E-5 1.E-1 1.E+3 1.E+7 1.E+11 Reduced time,s shift Tmax Data Shift New Arrhenius Shift New Poly Shift New WLF Figure 8. Graph. MN1-5 complex modulus master curve reference temperature -18 C. 16
23 Log at Arrhenius y = 1.31E+4x E-1 R² = E Series1 Linear (Series1) 1/(T )-1/(Tref ) Figure 9. Graph. MN1-5 Arrhenius shift function. Log at WLF data WLF T-Tref Figure 1. Graph. MN1-5 WLF shift function. 17
24 POLY Log at data fit Poly. (data) -5-1 y = E+6x E+4x E-4 R² = E-1-1.5E-3-1.E-3-5.E-4 1.E-17 5.E-4 1.E-3 1.5E-3 2.E-3 1/(T )-1/(Tref ) Figure 11. Graph. MN1-5 polynomial Arrhenius shift function. Phase Angle, degrees Black Space E+4 1.E+5 1.E+6 1.E+7 1.E+8 1.E+9 1.E+1 Complex Modulus, Pa Figure 12. Graph. MN1-5 black space. 18
25 1.E+9 shift G* Complex Modulus, Pa 1.E+8 1.E+7 1.E+6 1.E+5 1.E+4 1.E+3 1.E+2 1.E-12 1.E-6 1.E+ 1.E+6 Reduced time, s 1.E+12 shift Tmax Data Shift New Arrhenius Shift New Poly Shift New WLF Figure 13. Graph. YNP complex modulus master curve reference temperature -24 C. 3 Arrhenius Log at Series y = E+4x E-1 R² = E-1 Linear (Series1) /(T )-1/(Tref ) Figure 14. Graph. YNP Arrhenius shift function. 19
26 WLF Log at data WLF T-Tref Figure 15. Graph. YNP WLF shift function POLY Log at data fit Poly. (data) -5-1 y = E+6x E+3x E-2 R² = E-1-1.5E-3-1.E-3-5.E-4 1.E-17 5.E-4 1.E-3 1.5E-3 2.E-3 1/(T )-1/(Tref ) Figure 16. Graph. YNP polynomial Arrhenius shift function. 2
27 8 Black Space 7 6 Phase Angle, degrees E+2 1.E+3 1.E+4 1.E+5 1.E+6 1.E+7 1.E+8 1.E+9 Complex Modulus, Pa Figure 17. Graph. YNP black space. Prony Herd Calculation Results The Prony Herd spreadsheet does all of the interconversions of the DSR to finally arrive at stiffness. The spreadsheet contains graphical displays of these computations and Graphs of these for three binders are shown in figures for the Darcy flow model case of time dependent Poisson s ratio: 21
28 Prony Spectrum 9.E+7 8.E+7 7.E+7 6.E+7 g i, Pa 5.E+7 4.E+7 3.E+7 2.E+7 1.E+7.E+ 1.E-19 1.E-15 1.E-11 1.E-7 1.E-3 1.E+1 1.E+5 1.E+9 decay constant, s Figure 18. Graph. AZ1-3. Prony spectrum. 6 Phase angle 5 Phase angle degrees calc data 1.E-9 1.E-7 1.E-5 1.E-3 1.E-1 1.E+1 1.E+3 angular freq radians Figure 19. Graph. AZ1-3. Measured and calculated phase angle. 22
29 Complex Modulus G(w) 1.E+1 1.E+9 1.E+8 1.E+7 1.E+6 1.E+5 1.E+4 1.E+3 1.E+2 1.E+1 1.E+ 1.E-9 1.E-7 1.E-5 1.E-3 1.E-1 1.E+1 1.E+3 Calc Data angular freq radians Figure 2. Graph. AZ1-3. Complex modulus measured compared to complex modulus calculated from relaxation function convolution. Shear Relaxation Modulus G(t) 1.E+9 G(t) Calc Calc 1.E+8.E+ 5.E+1 1.E+2 1.5E+2 2.E+2 2.5E+2 3.E+2 time, s Figure 21. Graph. AZ1-3. Shear relaxation from Prony series. 23
30 Relaxation Modulus E(t) 1.E+1 G(t) 1.E+9 E(t) v=.5 E(t) using v(t) 1.E+8-1.E+1 4.E+1 9.E+1 1.4E+2 1.9E+2 2.4E+2 2.9E+2 time, s Figure 22. Graph. AZ1-3. Extensional relaxation from shear relaxation assuming.5 Poisson s ratio and time dependant Poisson s ratio. 6.E-9 Creep Compliance, D(t) 5.E-9 4.E-9 D(t) 3.E-9 2.E-9 Calc 1/E(t) 1.E-9.E+ 1.2E+ 5.1E+1 1.E+2 1.5E+2 2.E+2 2.5E+2 3.E+2 time, s Figure 23. Graph. AZ1-3. Creep compliance from deconvolution calculation and reciprocal extensional relaxation approximation. 24
31 Stiffness, S(t) 1.E+1 G(t) 1.E+9 Calc 1.E+8.E+ 5.E+1 1.E+2 1.5E+2 2.E+2 2.5E+2 3.E+2 time, s Figure 24. Graph. AZ1-3. Stiffness from calculated creep compliance. Prony Spectrum 1.E+8 1.E+8 8.E+7 g i, Pa 6.E+7 4.E+7 2.E+7.E+ 1.E-2 1.E-16 1.E-12 1.E-8 1.E-4 1.E+ 1.E+4 1.E+8 decay constant s Figure 25. Graph. MN1-5. Prony spectrum. 25
32 Phase angle Phase angle degrees calc data -1 1.E-9 1.E-7 1.E-5 1.E-3 1.E-1 1.E+1 1.E+3 1.E+5 angular freq radians Figure 26. Graph. MN1-5. Measured and calculated phase angle. Complex Modulus G(w) 1.E+1 1.E+9 1.E+8 1.E+7 1.E+6 1.E+5 1.E+4 1.E+3 1.E+2 1.E+1 1.E+ 1.E-9 1.E-6 1.E-3 1.E+ 1.E+3 1.E+6 Calc Data angular freq radians Figure 27. Graph. MN1-5. Complex modulus measured compared to complex modulus calculated from relaxation function convolution. 26
33 Relaxation Modulus, E(t) 1.E+1 G(t) 1.E+9 E(t) v=.5 E(t) using v(t) 1.E+8-1.E+1 9.E+1 1.9E+2 2.9E+2 time, s Figure 28. Graph. MN1-5. Shear relaxation from Prony series. 1.E+1 Relaxation Modulus, E(t) G(t) 1.E+9 E(t) v=.5 E(t) using v(t) 1.E+8-1.E+1 4.E+1 9.E+1 1.4E+2 1.9E+2 2.4E+2 2.9E+2 time, s Figure 29. Graph. MN1-5. Extensional relaxation from shear relaxation assuming.5 Poisson s ratio and time dependant Poisson s ratio. 27
34 Creep Compliance, D(t) 1.2E-8 1.E-8 8.E-9 D(t) 6.E-9 4.E-9 Calc 1/E(t) 2.E-9.E+ 1.2E+ 5.1E+1 1.E+2 1.5E+2 2.E+2 2.5E+2 3.E+2 time, s Figure 3. Graph. MN1-5. Creep compliance from deconvolution calculation and reciprocal extensional relaxation approximation. 1.E+1 Stiffness, S(t) G(t) 1.E+9 Calc 1.E+8.E+ 5.E+1 1.E+2 1.5E+2 2.E+2 2.5E+2 3.E+2 time, s Figure 31. Graph. MN1-5. Stiffness from calculated creep compliance. 28
35 g i, Pa Prony Spectrum 1.E+8 9.E+7 8.E+7 7.E+7 6.E+7 5.E+7 4.E+7 3.E+7 2.E+7 1.E+7.E+ 1.E-16 1.E-12 1.E-8 1.E-4 1.E+ 1.E+4 1.E+8 1.E+12 decay constant, s Figure 32. Graph. YNP. Prony spectrum. Phase angle 8 7 Phase angle degrees calc data 1.E-12 1.E-1 1.E-8 1.E-6 1.E-4 1.E-2 1.E+ angular freq radians Figure 33. Graph YNP. Measured and calculated phase angle. 29
36 G(w) Complex Modulus 1.E+9 1.E+8 1.E+7 1.E+6 1.E+5 1.E+4 1.E+3 1.E+2 1.E+1 Calc Data 1.E+ 1.E-12 1.E-1 1.E-8 1.E-6 1.E-4 1.E-2 1.E+ angular freq radians Figure 34. Graph. YNP. Complex modulus measured compared to complex modulus calculated from relaxation function convolution. Shear Relaxation Modulus, G(t) G(t) 1.E+1 1.E+9 1.E+8 1.E+7 1.E+6 1.E+5 1.E+4 1.E+3 1.E+2 1.E+1 1.E+.E+ 5.E+1 1.E+2 1.5E+2 2.E+2 2.5E+2 3.E+2 time, s Calc Calc Figure 35. Graph. YNP. Shear relaxation from Prony series. 3
37 Relaxation Modulus, E(t) 1.E+1 G(t) 1.E+9 E(t) v=.5 E(t) using v(t) 1.E+8-1.E+1 4.E+1 9.E+1 1.4E+2 1.9E+2 2.4E+2 2.9E+2 time, s Figure 36. Graph. YNP. Extensional relaxation from shear relaxation assuming.5 Poisson s ratio and time dependant Poisson s ratio. 8.E-9 7.E-9 6.E-9 5.E-9 Creep Compliance, D(t) D(t) 4.E-9 3.E-9 2.E-9 1.E-9.E+ 1.2E+ 5.1E+1 1.E+2 1.5E+2 2.E+2 2.5E+2 3.E+2 time, s Calc 1/E(t) Figure 37. Graph. YNP. Creep compliance from deconvolution calculation and reciprocal extensional relaxation approximation. 31
38 Stiffness, S(t) 1.E+1 G(t) 1.E+9 Series1 1.E+8.E+ 5.E+1 1.E+2 1.5E+2 2.E+2 2.5E+2 3.E+2 time, s Figure 38. Graph. YNP. Stiffness from calculated creep compliance. The interconverted DSR stiffnesses were compared to the BBR flexural stiffnesses for the three binders studied. In general, the slope of the stiffness curve was higher for the BBR testing than from the DSR. An example of the DSR and BBR stiffness curves over a time frame relevant to the determination of m-value is shown in figure
39 S(t), Pa 1.E+9 9.E+8 8.E+8 7.E+8 6.E+8 5.E+8 4.E+8 3.E+8 2.E+8 BBR Stiffness -18C MN-5 Mn-5 DSR stiffness - 18C Poisson's=.5 1.E Time, s Figure 39. Graph. DSR and BBR stiffness compared using a constant Poisson s ratio assumption in the interconversion. The question of whether small errors in temperature measurements in either instrument affected results was investigated by repeating the stiffness calculations at different shift temperatures. Figure 4 shows that some of the error in absolute terms could be attributed to a few degrees measurement error in the temperature, but a difference in slope persists. A drift to warmer temperatures would correct this, but seems much more unlikely than a drift to colder temperatures. Notice, however, that a.5 C change over 2 minutes time would cause this much error. 33
40 1.E+9 9.E+8 8.E+8 S(t), Pa 7.E+8 6.E+8 5.E+8 4.E+8 BBR -18C DSR -18C DSR -17C DSR -16C DSR -15C 3.E time, s Figure 4. Graph. The stiffness curves generated using a range of temperatures in the DSR. A range of stiffnesses were computed based upon a range of assumed constant Poisson s ratio in the interconversion calculation (figure 41). 34
41 S(t) Pa 1.E+9 9.E+8 8.E+8 7.E+8 6.E+8 5.E+8 4.E+8 3.E+8 2.E+8 1.E+8 BBR -18C Poisson's=.5 Poisson's=.4 Poisson's=.3 Poisson's=.2 Poisson's=.1 Poisson's=. Poisson's=-.1.E Time, s Figure 41. Graph. The stiffness curves generated using a range of constant Poisson ratio value assumptions. Poisson s ratio for a visco-elastic material is time dependent, but generally is considered to increase as the material relaxes and flows (Tschoegl et al. 22). Certainly this is true if the material is nearly incompressible. However, the data from the BBR suggest that the Poisson s ratio is actually decreasing. This is only possible if the material consists of two phases, a mobile phase and a stationary phase. Compression of the top of the beam would drive flow toward the bottom, and tension at the bottom would pull flow toward it. An approximate Darcy s law was derived earlier in the introduction to describe this kind of flow resulting in a decaying exponential description of the Poisson s ratio in time. The decaying Poisson s ratio interconversion from DSR data is shown in figure
42 S(t) Pa 1.E+9 9.E+8 8.E+8 7.E+8 6.E+8 5.E+8 4.E+8 3.E+8 2.E+8 1.E+8 BBR -18C Poisson's=.5 Poisson's=.4 Poisson's=.3 Poisson's=.2 Poisson's=.1.E Time, s Poisson's=. Figure 42. Graph. The stiffness curves generated using a range of constant Poisson ratio value assumptions and time dependent Poisson s ratio. The time dependent Poisson ratio parameters for the three binders studied, ν, time zero value, ν, infinite time value, and λ, the decay constant, are tabulated below in table 1 for the three binders studied. Table 1. Poisson s ratio exponential decay parameters for three binders. Binder MN1-5 AZ1-3 YNP ν ν λ (s) The BBR data were collected at three temperatures depending on lower PG grade temperature. These correlate nicely (figure 43) with the decay times, suggesting that the decay function, no matter its physical meaning, might be used to reconcile BBR and DSR data, since at the same temperatures, the Poisson decay function would essentially be the same. It may be possible to match DSR derived stiffness to BBR stiffness consistently to obtain m values. 36
43 BBR Temperature, Celsius y = 1.2x - 3 R² = Figure 43. Graph. Poisson s ratio decay times as a function of temperature. CONCLUSIONS The rigorous approach to interconversion of DSR complex modulus to BBR stiffness does not improve agreement between the two methods over approximate approaches unless a time dependant Poisson s ratio is employed in the conversion from shear relaxation to flexural relaxation. The physical explanation for the Poisson s ratio time dependent form is not clearly understood, but might be due to flow of a fluid component in the asphalt. Other explanations, such as softening in the BBR solvent bath may be the cause for the differences. A stretched exponential to describe the Prony series coefficients in the computation of the relaxation from DSR complex modulis removes the uniqueness issues with this computational method. Non-linear regression can be used to get relaxation directly from complex modulis master curves, or individual isotherms, without phase angle data. The computational method can be used to compute phase angles that compare well with those reported from the instrument. RECOMMENDATIONS The current embodiment of this calculation exists as a Visual Basic for Applications (VBA) program in Microsoft Excel. A more user friendly stand alone application for routine use would be better suited for use by the asphalt industry and research community. 37
44 ACKNOWLEDGMENTS The authors gratefully acknowledge the Federal Highway Administration, U.S. Department of Transportation, for financial support of this project under contract no. DTFH61-7D-5. DISCLAIMER This document is disseminated under the sponsorship of the Department of Transportation in the interest of information exchange. The United States Government assumes no liability for its contents or use thereof. The contents of this report reflect the views of Western Research Institute which is responsible for the facts and the accuracy of the data presented herein. The contents do not necessarily reflect the official views of the policy of the Department of Transportation. REFERENCES Christensen, D. W., and D. A. Anderson, 1992, Interpretation of dynamic mechanical test data for paving grade asphalt cements. J. Assoc. Asphalt Paving. Technol., 61: Farrar, M., C. Sui, S. Salmans, and Q. Qin, 215, Technical white paper FP 8 Determining the Low Temperature Rheological Properties of Asphalt Binder Using a Dynamic Shear Rheometer (DSR). Prepared for Federal Highway Administration, Contract No. DTFH61-7-D-5, Fundamental Properties of Asphalts and Modified Asphalts, III, March 215. Ferry, J. D., 198, Viscoelastic Properties of Polymers, 3rd ed., John Wiley and Sons, Inc., New York. Jongepier, R., and B. Kuilman, 1969, Characteristics of the rheology of bitumens. Proc., Association of Asphalt Paving Technologists, 38: Lakes, R. S., 1992, The time-dependent Poisson's ratio of viscoelastic materials can increase or decrease. Cellular Polymers, 11: Marasteanu, M. O., and D. A. Anderson, 1999, Improved Model for Bitumen Rheological Characterization. Eurobitume Workshop on Performance Related Properties for Bituminous Binders, Luxembourg, Park, S. W., and Y. R. Kim, 1999, Interconversion between Relaxation Modulus and Creep Compliance for Viscoelastic Solids. J. Mat. Civ. Eng., 11 (1): Sui, C., M. J. Farrar, W. H. Tuminello, and T. F. Turner, 21, New Technique for Measuring Low-Temperature Properties of Asphalt Binders with Small Amounts of Material. Transportation Research Record, 2179:
45 Tschoegl, N. W., W. G. Knauss, and I. Emri, 22, Poisson s Ratio in Linear Viscoelasticity A Critical Review. Mechanics of Time-Dependent Materials, 6:
Malaysian Journal of Civil Engineering 22(1) : (2010)
MODELLING THE LINEAR VISCOELASTIC RHEOLOGICAL PROPERTIES OF BASE BITUMENS Nur Izzi Md Yusoff 1, Marc-Stéphane Ginoux 1, Gordon Dan Airey 1, Mohd Rosli Hainin 2 1 Nottingham Transportation Engineering Centre
More informationThe 2S2P1D: An Excellent Linear Viscoelastic Model
The 2S2P1D: An Excellent Linear Viscoelastic Model Md. Yusoff 1, N. I., Monieur, D. 2, Airey, G. D. 1 Abstract An experimental campaign has been carried out on five different unaged and five aged penetration
More informationContinuous Models for Characterizing Linear Viscoelastic Behavior of Asphalt Binders
Continuous Models for Characterizing Linear Viscoelastic Behavior of Asphalt Binders Dr. David A. Anderson Professor Emeritus Penn State University -and - Dr. Mihai Marasteanu Associate Professor University
More informationCookbook for Rheological Models Asphalt Binders
CAIT-UTC-062 Cookbook for Rheological Models Asphalt Binders FINAL REPORT May 2016 Submitted by: Offei A. Adarkwa PhD in Civil Engineering Nii Attoh-Okine Professor Pamela Cook Unidel Professor University
More informationRELIABILITY OF RHEOMETRIC MEASUREMENTS IN BITUMENS BY MEANS OF DYNAMIC SHEAR RHEOMETERS
RELIABILITY OF RHEOMETRIC MEASUREMENTS IN BITUMENS BY MEANS OF DYNAMIC SHEAR RHEOMETERS Antonio MONTEPARA, University of Parma, Parma, Italy Felice GIULIANI, University of Parma, Parma, Italy 1. INTRODUCTION
More informationPredicting Oxidative Age Hardening in Asphalt Pavement
Predicting Oxidative Age Hardening in Asphalt Pavement with Time and Depth: A First Approximation Michael J. Farrar (mfarrar@uwyo.edu) Will Grimes Fred Turner P. Michael Harnsberger Shin-Che Huang P3 Symposium
More informationSPECTRAL ANALYSIS AND THE INTERCONVERSION OF LINEAR VISCOELASTIC FUNCTIONS
SPECTRAL ANALYSIS AND THE INTERCONVERSION OF LINEAR VISCOELASTIC FUNCTIONS 1 INTRODUCTION TA Instruments Rheology Advantage spectral analysis and interconversion software was developed in collaboration
More informationALTERNATE SHIFT FACTOR RELATIONSHIP FOR DESCRIBING THE TEMPERATURE DEPENDENCY OF THE VISCO-ELASTIC BEHAVIOUR OF ASPHALT MATERIALS
ALTERNATE SHIFT FACTOR RELATIONSHIP FOR DESCRIBING THE TEMPERATURE DEPENDENCY OF THE VISCO-ELASTIC BEHAVIOUR OF ASPHALT MATERIALS Dr. Geoffrey M. Rowe* Abatech Inc. PO Box 356 Blooming Glen, PA 18911 +1
More informationWhy Characterize Binder Rheology With Mathematical Models?
Historical and Current Rheological Binder Characterization vs. Binder Performance Dr. David A. Anderson Consultant Dr. Geoffrey M. Rowe ABATECH INC. Dr. Donald Christensen Advanced Asphalt Technologies,
More informationRheology of Asphalt Binders and Implications for Performance
Rheology of Asphalt Binders and Implications for Performance November 10, 2009 Dr. David Anderson Consultant and Professor Emeritus Penn State University DA.SC@COMCAST.NET Dr. Geoffrey Rowe President Abatech
More informationInterconversion of Dynamic Modulus to Creep Compliance and Relaxation Modulus: Numerical Modeling and Laboratory Validation
Final Report September 2016 Interconversion of Dynamic Modulus to Creep Compliance and Relaxation Modulus: Numerical Modeling and Laboratory Validation SOLARIS Consortium, Tier 1 University Transportation
More informationPrediction of Complex Shear Modulus and Fracture Properties of Asphalt Binders with Oxidative Aging
Prediction of Complex Shear Modulus and Fracture Properties of Asphalt Binders with Oxidative Aging By Haifang Wen, PhD, PE Colf Distinguished Professor Director of Washington Center for Asphalt Technology
More informationGeoffrey Rowe, Abatech. 51st Petersen Asphalt Research Conference July 14-16, 2014 Laramie, Wyoming
Geoffrey Rowe, Abatech 51st Petersen Asphalt Research Conference July 14-16, 2014 Laramie, Wyoming 1 Work collected data, BBR DSR Review of Interconversions analysis DSR, use of G* and δ instead of S and
More informationLecture 7 Constitutive Behavior of Asphalt Concrete
Lecture 7 Constitutive Behavior of Asphalt Concrete What is a Constitutive Model? A constitutive model or constitutive equation is a relation between two physical quantities that is specific to a material
More informationDYNAMIC AND TRANSIENT TESTING OF ASPHALT BINDER AND PAVING MIX
66 6th RILEM Symposium PTEBM'03, Zurich, 2003 DYNAMIC AND TRANSIENT TESTING OF ASPHALT BINDER AND PAVING MIX L. Zanzotto, O.J. Vacin and J. Stastna University of Calgary, Canada Abstract: A commercially
More informationRheological Properties and Fatigue Resistance of Crumb Rubber Modified Bitumen
Rheological Properties and Fatigue Resistance of Crumb Rubber Modified Bitumen F. Khodary Department of Civil Engineering, Institute of traffic and transport, section of road and pavement engineering,
More informationAbvanced Lab Course. Dynamical-Mechanical Analysis (DMA) of Polymers
Abvanced Lab Course Dynamical-Mechanical Analysis (DMA) of Polymers M211 As od: 9.4.213 Aim: Determination of the mechanical properties of a typical polymer under alternating load in the elastic range
More informationGM Rowe and DA Anderson. Binder ETG Meeting, April 2014
GM Rowe and DA Anderson Binder ETG Meeting, April 2014 What is relaxation spectra Historical perspective What it means Other simplifications Field correlations Black space Discrete Model the asphalt by
More informationRheology considerations for recycled and recovered binder
Rheology considerations for recycled and recovered binder Geoffrey M. Rowe PAVEMENT PERFORMANCE PREDICTION SYMPOSIUM - Binder Issues Affecting the Durability of RAP and WMA Pavements Sponsored by Federal
More informationMSE 383, Unit 3-3. Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept.
Dynamic Mechanical Behavior MSE 383, Unit 3-3 Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept. Scope Why DMA & TTS? DMA Dynamic Mechanical Behavior (DMA) Superposition Principles
More informationRole of Binders in Pavement Performance
Role of Binders in Pavement Performance Presented by H. Bahia Research conducted by The University of Wisconsin-Asphalt Group The Pavement Performance Prediction Symposium Western Research Institute, Laramie,
More informationPhenomenological models for binder rutting and fatigue. University of Wisconsin Research Team
Phenomenological models for binder rutting and fatigue University of Wisconsin Research Team Outline of Talk The need for measuring damage behavior Current efforts, and evolution of ideas Gaps and what
More informationMechanical Models for Asphalt Behavior and Performance
Mechanical Models for Asphalt Behavior and Performance Introduction and Review of Linear Viscoelastic Behaviors About the webinar series Past, current, and future plan for webinar series Introduction to
More informationUSE OF BBR TEST DATA TO ENHANCE THE ACCURACY OF G* -BASED WITCZAK MODEL PREDICTIONS
USE OF BBR TEST DATA TO ENHANCE THE ACCURACY OF G* -BASED WITCZAK MODEL PREDICTIONS Mekdim T. Weldegiorgis * PhD. Candidate, Department of Civil Engineering, University of New Mexico, USA * 1 University
More informationRheology Based Specifications for Asphalt Binders used in Roofing
Rheology Based Specifications for Asphalt Binders used in Roofing By: Gaylon L. Baumgardner, Geoffrey M. Rowe, and Isaac L. Howard Introduction Specification requirements for roofing and industrial asphalts
More informationG. R. Strobl, Chapter 5 "The Physics of Polymers, 2'nd Ed." Springer, NY, (1997). J. Ferry, "Viscoelastic Behavior of Polymers"
G. R. Strobl, Chapter 5 "The Physics of Polymers, 2'nd Ed." Springer, NY, (1997). J. Ferry, "Viscoelastic Behavior of Polymers" Chapter 3: Specific Relaxations There are many types of relaxation processes
More informationNonlinear oscillatory shear tests applied to bitumen: analysis and application of large amplitude oscillatory shear (LAOS)
Nonlinear oscillatory shear tests applied to bitumen: analysis and application of large amplitude oscillatory shear (LAOS) Alec Otto Cookman 1, a, Michael J. Farrar 1, Jean-Pascal Planche 1 1 Western Research
More informationMechanical characterization of visco termo elastic properties of a polymer interlayer by dynamic tests
Mechanical characterization of visco termo elastic properties of a polymer interlayer by dynamic tests Laura Andreozzi 1, Silvia Briccoli Bati 2, Mario Fagone 2 Giovanna Ranocchiai 2 Fabio Zulli 1 1 Department
More informationAdvances in performance evaluation of asphalt binders
Advances in performance evaluation of asphalt binders Hussain U Bahia Professor University of Wisconsin, USA March 25, Lisbon, Potugal s:, design performance a, LNEC, 25 March 2010 Road Outline >Performance
More informationMATERIAL MECHANICS, SE2126 COMPUTER LAB 3 VISCOELASTICITY. k a. N t
MATERIAL MECHANICS, SE2126 COMPUTER LAB 3 VISCOELASTICITY N t i Gt () G0 1 i ( 1 e τ = α ) i= 1 k a k b τ PART A RELAXING PLASTIC PAPERCLIP Consider an ordinary paperclip made of plastic, as they more
More informationThe Effect of Aging on Binder Properties of Porous Asphalt Concrete
The Effect of Aging on Binder Properties of Porous Asphalt Concrete Eyassu T. Hagos Raveling Supervisors: Prof. Dr. Ir. A.A.A. Molenaar Assoc. Prof. Ir. M.F.C. van de Ven Porous Asphalt layer March 27,
More informationA Discussion of Peltier Cooling Compared to Chilled Air Cooling and 4 mm DSR Results (and related issues)
A Discussion of Peltier Cooling Compared to Chilled Air Cooling and 4 mm DSR Results (and related issues) Gerald Reinke MTE Services Binder ETG, May 3, 2017 Ames, IA STATEMENT OF PROBLEM In the course
More informationCONSISTENCY OF RHEOLOGICAL EXPERIMENTS FOR PSA CHARACTERIZATION
CONSISTENCY OF RHEOLOGICAL EXPERIMENTS FOR PSA CHARACTERIZATION Dr. Laura Yao, Senior Research Chemist, Scapa North America, Windsor, CT Robert Braiewa, Research Chemist, Scapa North America, Windsor,
More informationFlow and Transport. c(s, t)s ds,
Flow and Transport 1. The Transport Equation We shall describe the transport of a dissolved chemical by water that is traveling with uniform velocity ν through a long thin tube G with uniform cross section
More informationA Comparative Study of the Influence of Shift Factor Equations on Master Curve Construction
Technical Paper ISSN 997- Int. J. Pavement Res. Technol. (6): Copyright @ Chinese Society of Pavement Engineering A Comparative Study of the Influence of Shift Factor Equations on Master Curve Construction
More informationLINEAR AND NON-LINEAR VISCOELASTIC BEHAVIOUR OF BINDERS AND ASPHALTS
UNIVERSITY OF NOTTINGHAM SCHOOL OF CIVIL ENGINEERING LINEAR AND NON-LINEAR VISCOELASTIC BEHAVIOUR OF BINDERS AND ASPHALTS By Behzad Rahimzadeh, BSc, MSc, MIAT, MIHT Thesis submitted to the University of
More informationAblaos Asphalt Binder Large Amplitude Oscillatory Shear
Michael J. Farrar mfarrar@uwyo.edu Jean-Pascal Planche Alec Otto Cookman Steve Salmans Ablaos Asphalt Binder Large Amplitude Oscillatory Shear Binder ETG Baton Rouge, LA Sept. 16, 2014 Outline Linear Viscoelasticity
More informationAnalysis of Non-Linear Dynamic Behaviours in Asphalt Concrete Pavements Under Temperature Variations
ENOC 2017, June 25 30, 2017, Budapest, Hungary Analysis of Non-Linear Dynamic Behaviours in Asphalt Concrete Pavements Under Temperature Variations Amal Abdelaziz *, Chun-Hsing Ho *, and Junyi Shan * *
More informationLinear viscoelastic behavior
Harvard-MIT Division of Health Sciences and Technology HST.523J: Cell-Matrix Mechanics Prof. Ioannis Yannas Linear viscoelastic behavior 1. The constitutive equation depends on load history. 2. Diagnostic
More informationThe Rheology Handbook
Thomas G. Mezger The Rheology Handbook For users of rotational and oscillatory rheometers 2nd revised edition 10 Contents Contents 1 Introduction 16 1.1 Rheology, rheometry and viscoelasticity 16 1.2 Deformation
More informationRheology of Soft Materials. Rheology
Τ Thomas G. Mason Department of Chemistry and Biochemistry Department of Physics and Astronomy California NanoSystems Institute Τ γ 26 by Thomas G. Mason All rights reserved. γ (t) τ (t) γ τ Δt 2π t γ
More informationA.C. Pronk & A.A.A. Molenaar Road & Railway Engineering, Delft University of Technology, Delft, The Netherlands
The Modified Partial Healing Model used as a Prediction Tool for the Complex Stiffness Modulus Evolutions in Four Point Bending Fatigue Tests based on the Evolutions in Uni-Axial Push-Pull Tests. A.C.
More informationA Phenomenological Model for Linear Viscoelasticity of Monodisperse Linear Polymers
Macromolecular Research, Vol. 10, No. 5, pp 266-272 (2002) A Phenomenological Model for Linear Viscoelasticity of Monodisperse Linear Polymers Kwang Soo Cho*, Woo Sik Kim, Dong-ho Lee, Lee Soon Park, Kyung
More informationUsing Rheological Properties to Evaluate Storage Stability and Setting Behaviors of Emulsified asphalts
Using Rheological Properties to Evaluate Storage Stability and Setting Behaviors of Emulsified asphalts Huachun Zhai, Delmar Salomon * and Eric Milliron Idaho Asphalt Supply, Inc. 2627 Brandt Ave., Nampa,
More informationChapter 6 Molten State
Chapter 6 Molten State Rheology ( 流變學 ) study of flow and deformation of (liquid) fluids constitutive (stress-strain) relation of fluids shear flow shear rate ~ dγ/dt ~ velocity gradient dv 1 = dx 1 /dt
More informationComments on Use of Reference Fluid to Verify DSR
Comments on Use of Reference Fluid to Verify DSR David Anderson Professor Emeritus Penn State FHWA Asphalt Binder Expert Task Group Baton Rouge, LA September 16-17, 2014 Reference fluid how and why? Used
More informationMechanical Models for Asphalt Behavior and Performance
Mechanical Models for Asphalt Behavior and Performance All Attendees Are Muted Questions and Answers Please type your questions into your webinar control panel We will read your questions out loud, and
More informationEXPERIMENTALLY DETERMINING THE VISCOELASTIC BEHAVIOR OF A CURING THERMOSET EPOXY R. Thorpe 1, A. Poursartip 1*
19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS EXPERIMENTALLY DETERMINING THE VISCOELASTIC BEHAVIOR OF A CURING THERMOSET EPOXY R. Thorpe 1, A. Poursartip 1* 1 Composites Group, Dept. of Materials
More informationViscoelasticity, Creep and Oscillation Experiment. Basic Seminar Applied Rheology
Viscoelasticity, Creep and Oscillation Experiment Basic Seminar Applied Rheology Overview Repetition of some basic terms Viscoelastic behavior Experimental approach to viscoelasticity Creep- and recovery
More informationNEW APPROACH FOR THE FATIGUE CHARACTERISATION OF BITUMINOUS BINDERS
424 6th RILEM Symposium PTEBM'3, Zurich, 23 NEW APPROACH FOR THE FATIGUE CHARACTERISATION OF BITUMINOUS BINDERS Sándor Tóth Technical and Information Services on National Roads, Hungary Róbert Perlaki
More informationDetermination of the Relaxation Modulus of a Linearly Viscoelastic Material
Determination of the Relaxation Modulus of a Linearly Viscoelastic Material Joonas Sorvari and Matti Malinen Department of Physics, University of Kuopio, PO Box 1627, FIN-7211 Kuopio, Finland; E-mail joonas.sorvari@uku.fi
More informationRock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth
Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth References: Turcotte and Schubert, Geodynamics, Sections 2.1,-2.4, 2.7, 3.1-3.8, 6.1, 6.2, 6.8, 7.1-7.4. Jaeger and Cook, Fundamentals of
More informationEffect of polyphosphoric acid on bituminous binder oxidation
Effect of polyphosphoric acid on bituminous binder oxidation Ron Glaser 1, a, Fred Turner 1, b, Giovanni Onnembo 2, c, Jean-Pascal Planche 1, d 1 Western Research Institute, Laramie, WY, United States
More informationInvestigation on the cooling medium effect in the characterization of asphalt binder with the Bending Beam Rheometer (BBR)
Investigation on the cooling medium effect in the characterization of asphalt binder with the Bending Beam Rheometer (BBR) Journal: Manuscript ID cjce-2017-0586.r1 Manuscript Type: Article Date Submitted
More informationElements of Polymer Structure and Viscoelasticity. David M. Parks Mechanics and Materials II February 18, 2004
Elements of Polymer Structure and Viscoelasticity David M. Parks Mechanics and Materials II 2.002 February 18, 2004 Outline Elements of polymer structure Linear vs. branched; Vinyl polymers and substitutions
More informationDynamic Mechanical Analysis of Solid Polymers and Polymer Melts
Polymer Physics 2015 Matilda Larsson Dynamic Mechanical Analysis of Solid Polymers and Polymer Melts Polymer & Materials Chemistry Introduction Two common instruments for dynamic mechanical thermal analysis
More informationDynamic Modulus of Asphalt Mixtures for Development of
JOURNAL PAPER TYPE Journal of Testing and Evaluation, Sept 2002 Vol XX, No. X Paper ID: Available online at: www.astm.org Published XXXX Lee, K. 1, Kim, H. 2, Kim, N. 3, Kim, Y 4 Dynamic Modulus of Asphalt
More informationCM4655 Polymer Rheology Lab. Torsional Shear Flow: Parallel-plate and Cone-and-plate
CM4655 Polymer heology Lab Torsional Shear Flow: Parallel-plate and Cone-and-plate (Steady and SAOS) Professor Faith A. Morrison Department of Chemical Engineering Michigan Technological University r (-plane
More informationEstimating Damage Tolerance of Asphalt Binders Using the Linear Amplitude Sweep
Standard Method of Test for Estimating Damage Tolerance of Asphalt Binders Using the Linear Amplitude Sweep AASHTO Designation: TP 101-14 American Association of State Highway and Transportation Officials
More informationFINITE ELEMENT APPROACHES TO MESOSCOPIC MATERIALS MODELING
FINITE ELEMENT APPROACHES TO MESOSCOPIC MATERIALS MODELING Andrei A. Gusev Institute of Polymers, Department of Materials, ETH-Zürich, Switzerland Outlook Generic finite element approach (PALMYRA) Random
More informationThe determination of creep and relaxation functions from a single experiment
The determination of creep and relaxation functions from a single experiment A. Nikonov Center for Experimental Mechanics, University of Ljubljana, Ljubljana, Slovenia A. R. Davies Institute of Mathematical
More informationLOW TEMPERATURE RHEOLOGICAL AND FRACTURE PROPERTIES OF POLYMER-MODIFIED BITUMENS
LOW TEMPERATURE RHEOLOGICAL AND FRACTURE PROPERTIES OF POLYMER-MODIFIED BITUMENS L.Lapalu 1, J.P. Planche 1, D. Martin 1, D. Anderson 2, J.F. Gérard 3 1 ELF Research Center (Solaize) 2 The Pennsylvania
More information5 The Oldroyd-B fluid
5 The Oldroyd-B fluid Last time we started from a microscopic dumbbell with a linear entropic spring, and derived the Oldroyd-B equations: A u = u ρ + u u = σ 2 pi + η u + u 3 + u A A u u A = τ Note that
More informationTemperature Effects on LIGO Damped Coil Springs
LIGO-T97241--D HYTEC-TN-LIGO-18 Temperature Effects on LIGO Damped Coil Springs Franz Biehl July 28, 1997 Abstract A simple method for estimating a damping material stiffness and loss factor as a function
More informationViscoelasticity. Basic Notions & Examples. Formalism for Linear Viscoelasticity. Simple Models & Mechanical Analogies. Non-linear behavior
Viscoelasticity Basic Notions & Examples Formalism for Linear Viscoelasticity Simple Models & Mechanical Analogies Non-linear behavior Viscoelastic Behavior Generic Viscoelasticity: exhibition of both
More informationBitumen fatigue performance evaluation - with or without RAP - a real challenge
Bitumen fatigue performance evaluation - with or without RAP - a real challenge Ronald van Rooijen 1, a, Luc Planque 2, b, Marcel Wegmann 1, c, Hervé Di Benedetto 3, d, Salvatore Mangiafico 3, e, Cedric
More informationEshan V. Dave, Secretary of M&FGM2006 (Hawaii) Research Assistant and Ph.D. Candidate. Glaucio H. Paulino, Chairman of M&FGM2006 (Hawaii)
Asphalt Pavement Aging and Temperature Dependent Properties through a Functionally Graded Viscoelastic Model I: Development, Implementation and Verification Eshan V. Dave, Secretary of M&FGM2006 (Hawaii)
More informationFundamentals of Polymer Rheology. Sarah Cotts TA Instruments Rubber Testing Seminar CUICAR, Greenville SC
Fundamentals of Polymer Rheology Sarah Cotts TA Instruments Rubber Testing Seminar CUICAR, Greenville SC Rheology: An Introduction Rheology: The study of stress-deformation relationships =Viscosity =Modulus
More informationHow to measure the shear viscosity properly?
testxpo Fachmesse für Prüftechnik 10.-13.10.2016 How to measure the shear viscosity properly? M p v Rotation Capillary Torsten Remmler, Malvern Instruments Outline How is the Shear Viscosity defined? Principle
More informationSTUDY OF RHEOLOGICAL PROPERTIES OF BITUMINOUS BINDERS IN MIDDLE AND HIGH TEMPERATURES
Vol. 12, Issue 1/2016, 13-20 DOI: 10.1515/cee-2016-0002 STUDY OF RHEOLOGICAL PROPERTIES OF BITUMINOUS BINDERS IN MIDDLE AND HIGH TEMPERATURES Eva REMIŠOVÁ 1,*, Viera ZATKALIKOVÁ 2, František SCHLOSSER
More informationEstimating Fatigue Resistance Damage Tolerance of Asphalt Binders Using the Linear Amplitude Sweep
Standard Method of Test for Estimating Fatigue Resistance Damage Tolerance of Asphalt Binders Using the Linear Amplitude Sweep AASHTO Designation: TP 2b xx (LAST)101-1214 American Association of State
More informationPrediction of the complex modulus of asphalt mixture containing RAP materials from the rheological properties of mortars
1 1 1 1 1 1 1 1 0 1 0 1 0 1 Prediction of the complex modulus of asphalt mixture containing RAP materials from the rheological properties of mortars Chiara Riccardi 1, Augusto Cannone Falchetto *, and
More informationProposed Nomenclature for Steady Shear Flow and Linear Viscoelastic Behavior
Proposed Nomenclature for Steady Shear Flow and Linear Viscoelastic Behavior C. L. Sieglaff Citation: Transactions of the Society of Rheology 20, 311 (1976); doi: 10.1122/1.549415 View online: https://doi.org/10.1122/1.549415
More informationBohlin. Rheological instruments backed with rheological experience. Rheological properties
Rheological properties ṙ Bohlin Rheological instruments backed with rheological experience detailed specification sheets from /bohlingemini The Bohlin Gemini rheometers are an advanced range of compact
More informationRILEM inter-laboratory test on binder rheology test methods DSR and BBR.
RILEM inter-laboratory test on binder rheology test methods DSR and BBR. Dariusz Sybilski, IBDiM Poland Ann Vanelstrate, CRR Belgium RILEM Technical Committee 182 PEB Performance Testing and Evaluation
More informationViscoelastic Mechanical Analysis for High Temperature Process of a Soda-Lime Glass Using COMSOL Multiphysics
Viscoelastic Mechanical Analysis for High Temperature Process of a Soda-Lime Glass Using COMSOL Multiphysics R. Carbone 1* 1 Dipartimento di Ingegneria dei Materiali e della Produzione sez. Tecnologie
More informationThe Large Amplitude Oscillatory Strain Response of Aqueous Foam: Strain Localization and Full Stress Fourier Spectrum
The Large Amplitude Oscillatory Strain Response of Aqueous Foam: Strain Localization and Full Stress Fourier Spectrum By F. Rouyer, S. Cohen-Addad, R. Höhler, P. Sollich, and S.M. Fielding The European
More informationMeasuring structure of low viscosity fluids in oscillation using rheometers with and without a separate torque transducer
ANNUAL RANSACONS OF HE NORDC RHEOLOGY SOCEY, VOL., 003 Measuring structure of low viscosity fluids in oscillation using rheometers with and without a separate torque transducer Aly Franck A nstruments,
More informationMeasurement Engineering Group, Paderborn University, Warburger Straße 100, Paderborn, Germany
Nadine Feldmann 1, Fabian Bause 1, Bernd Henning 1 1 Measurement Engineering Group, Paderborn University, Warburger Straße 100, 33098 Paderborn, Germany feldmann@emt.uni-paderborn.de Abstract The present
More informationAnalysis of Hot-Mix Asphalt (HMA) Linear Viscoelastic and Bimodular Properties Using Uniaxial Compression and Indirect Tension (IDT) Tests
Analysis of Hot-Mix Asphalt (HMA) Linear Viscoelastic and Bimodular Properties Using Uniaxial Compression and Indirect Tension (IDT) Tests by Samer W. Katicha Dissertation submitted to the Faculty of the
More informationINVESTIGATING THE RHEOLOGICAL CHARACTERISTICS OF SOUTH AFRICAN ROAD BITUMENS
INVESTIGATING THE RHEOLOGICAL CHARACTERISTICS OF SOUTH AFRICAN ROAD BITUMENS G. MTURI 1, J. O CONNELL 2 and S.E. ZOOROB 3 1 Researcher (gmturi@csir.co.za), 2 Researcher, 3 Senior Researcher CSIR Built
More informationLecture 7: Rheology and milli microfluidic
1 and milli microfluidic Introduction In this chapter, we come back to the notion of viscosity, introduced in its simplest form in the chapter 2. We saw that the deformation of a Newtonian fluid under
More informationRheology of cellulose solutions. Puu Cellulose Chemistry Michael Hummel
Rheology of cellulose solutions Puu-23.6080 - Cellulose Chemistry Michael Hummel Contents Steady shear tests Viscous flow behavior and viscosity Newton s law Shear thinning (and critical concentration)
More informationViscoelastic Flows in Abrupt Contraction-Expansions
Viscoelastic Flows in Abrupt Contraction-Expansions I. Fluid Rheology extension. In this note (I of IV) we summarize the rheological properties of the test fluid in shear and The viscoelastic fluid consists
More informationRHEOLOGICAL COEFFICIENTS FOR MEDIA WITH MECHANICAL RELAXATION PHENOMENA
Communications to SIMAI Congress, ISSN 187-915, Vol. (7) DOI: 1.1685/CSC6157 RHEOLOGICAL COEFFICIENTS FOR MEDIA WITH MECHANICAL RELAXATION PHENOMENA A. CIANCIO, V. CIANCIO Department of Mathematics, University
More informationDrilling Fluid Thixotropy & Relevance
ANNUAL TRANSACTIONS OF THE NORDIC RHEOLOGY SOCIETY, VOL. 13, 2005 Drilling Fluid Thixotropy & Relevance Richard Jachnik1, 1Baker Hughes INTEQ, Stoneywood Park North, Dyce, Aberdeen, Scotland, UK ABSTRACT
More informationModule III - Macro-mechanics of Lamina. Lecture 23. Macro-Mechanics of Lamina
Module III - Macro-mechanics of Lamina Lecture 23 Macro-Mechanics of Lamina For better understanding of the macromechanics of lamina, the knowledge of the material properties in essential. Therefore, the
More informationVECD (Visco-ElasticContinuum Damage): State-of-the-art technique to evaluate fatigue damage in asphalt pavements
VECD (Visco-ElasticContinuum Damage): tate-of-the-art technique to evaluate fatigue damage in asphalt pavements M. Emin Kutay, Ph.D., P.E. Assistant Professor Michigan tate University History of the Viscoelastic
More informationFigure 1. Dimension of PSA in face paper laminate
EVALUATION OF DYNAMIC MECHANICAL PROPERTIES OF PRESSURE SENSITIVE ADHESIVE IN PAPER LANIMATES FOR POSTAGE STAMP APPLICATION: METHOD DEVELOPMENT AND ADHESIVE CHARACTERIZATION Hailing Yang, Graduate Student,
More informationInterfacial dynamics
Interfacial dynamics Interfacial dynamics = dynamic processes at fluid interfaces upon their deformation Interfacial rheological properties: elasticity, viscosity, yield stress, Relation between macroscopic
More informationSound Propagation through Media. Nachiketa Tiwari Indian Institute of Technology Kanpur
Sound Propagation through Media Nachiketa Tiwari Indian Institute of Technology Kanpur LECTURE-13 WAVE PROPAGATION IN SOLIDS Longitudinal Vibrations In Thin Plates Unlike 3-D solids, thin plates have surfaces
More informationBending of Simply Supported Isotropic and Composite Laminate Plates
Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,
More informationMulti-mode revisited
Multi-mode revisited Testing the application of shift factors S.J.M Hellenbrand 515217 MT 7.29 Coaches: Ir. L.C.A. van Breemen Dr. Ir. L.E. Govaert 2-7- 7 Contents Contents 1 Introduction 2 I Polymers
More informationLecture 4: viscoelasticity and cell mechanics
Teaser movie: flexible robots! R. Shepherd, Whitesides group, Harvard 1 Lecture 4: viscoelasticity and cell mechanics S-RSI Physics Lectures: Soft Condensed Matter Physics Jacinta C. Conrad University
More informationSEISMOLOGY. Master Degree Programme in Physics - UNITS Physics of the Earth and of the Environment ANELASTICITY FABIO ROMANELLI
SEISMOLOGY Master Degree Programme in Physics - UNITS Physics of the Earth and of the Environment ANELASTICITY FABIO ROMANELLI Department of Mathematics & Geosciences University of Trieste romanel@units.it
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationThermal Analysis of Polysaccharides Mechanical Methods
Biopolymer Solutions Thermal Analysis of Polysaccharides Mechanical Methods John Mitchell John.Mitchell@biopolymersolutions.co.uk 1 Topics Covered Introduction to polymer viscoelasticity Examples Thermal
More informationPERFORMANCE BASED EVALUATION OF CRACKING IN ASPHALT CONCRETE USING VISCOELASTIC AND FRACTURE PROPERTIES
University of New Hampshire University of New Hampshire Scholars' Repository Doctoral Dissertations Student Scholarship Spring 2018 PERFORMANCE BASED EVALUATION OF CRACKING IN ASPHALT CONCRETE USING VISCOELASTIC
More informationVISCOELASTIC PROPERTIES OF POLYMERS
VISCOELASTIC PROPERTIES OF POLYMERS John D. Ferry Professor of Chemistry University of Wisconsin THIRD EDITION JOHN WILEY & SONS New York Chichester Brisbane Toronto Singapore Contents 1. The Nature of
More informationYield Energy of Asphalt Binders Using the Dynamic Shear Rheometer
Standard Method of Test for Yield Energy of Asphalt Binders Using the Dynamic Shear Rheometer AASHTO Designation: T XXX-09 1. SCOPE 1.1. This test method covers the indication of asphalt binders resistance
More information