Page 172 - MODELING OF ASPHALT CONCRETE
P. 172

150    Cha pte r  S i x


                    {r , E (i = 1,…, m) and E } or {t , D  (j = 1,…, n), D and h } are known and the target time
                      i  i              e    j  j           j     0
                    constants specified, the unknown set of constants can be determined through a system
                    of linear algebraic equations. For example, creep compliance in its Prony series form, D
                                                                                            j
                    (j   1,…, n), can be determined from relaxation modulus E(t) as follows:
                                                   []{} =  { }                          (6-23)
                                                           B
                                                    AD
                    or A D = B (summed on j; j = 1,…, n; k = 1,…, p) where
                        kj  j  k
                                                        E
                                                    m
                                                         i
                                       E
                                 A = { (1 −  e −(/τ j ) ) + ∑  ρ ρ i − τ  e (  −(/ρ i ) i  −  e −  k t ( / )ρ i  )  ρ i  ≠ τ j  (6-24)
                                              t k
                                                             t k
                                        e
                                   kj
                                                   i=1  i  j
                    or
                                                     m
                                       E (1 −  e − (/τ j ) ) + ∑  tE i  e (  − (/ρ i ) )  ρ i  = τ j
                                                       k
                                              t k
                                                       τ
                                                            t k
                                        e
                                                     =
                                                    i 1  j
                    and
                                              ⎛     m       ⎞  ⎛    m   ⎞ ⎞
                                                                 e ∑
                                                 e ∑
                                        B =−1  ⎜ ⎝ E +  i=1  E e −(/ρ i ) ⎟ ⎠  ⎜ ⎝  E +  i=1 E i ⎟ ⎠  (6-25)
                                                         t k
                                                       i
                                         k
                       The symbol t  (k   1,…, p) denotes a discrete time corresponding to the upper limit
                                  k
                    of integration in Eq. (6-21). Once the model constants D , D , and t are determined, then
                                                                  g  j     j
                    function D(t) can be obtained in its Prony series form. Similarly, E(t) can be determined
                    from D(t) by determining the unknown constants of the Prony series representation of
                    E(t): E (i   1,…, m). As for D and E , they are related as follows:
                         i                  g     e
                                                         1
                                                  D ≡
                                                   g      m
                                                       e ∑
                                                      E +   E i
                                                          i=1
                                                                                        (6-26)
                                                  E ≡    1
                                                   e      n
                                                      D + ∑ D
                                                       g     j
                                                          j=1
                       Equations (6-23) through (6-26) were used in converting E(t) data of an AC mixture
                    to D(t). The source and target Prony series coefficients are presented in Table 6-1, with
                    corresponding fitted data plotted in Fig. 6-5.
                       For additional details and derivations of the above equations, the reader is referred
                    to Park and Schapery (1999) where other approximate interconversion techniques
                    available in literature are documented and compared. Only a limited number of common
                    techniques are presented in the next section.
                    Quasi-Elastic Interrelationship
                    The most intuitive and crude interrelationship between the E(t) and D(t) is the one
                    based on quasi-elastic approximation:
                                               Et D t() ≅ 1  (for > 0)                  (6-27)
                                                ()
                                                             t
   167   168   169   170   171   172   173   174   175   176   177