---
title: Four-Parameter Fractional Maxwell Model
url: https://www.emergentmind.com/topics/four-parameter-fractional-maxwell-model-fmm
type: topic
---

# Four-Parameter Fractional Maxwell Model

The four-parameter Fractional Maxwell Model (FMM) is a class of Maxwell-type linear viscoelastic constitutive models in which one replaces the integer-order operators of classical spring–dashpot rheology by fractional operators, thereby obtaining relaxation, creep, and oscillatory responses governed by power laws or Mittag–Leffler functions rather than single exponentials. In the cited literature, the label refers to several closely related parametrizations: a series combination of two Scott–Blair elements with parameters \((\mathbb{E}_1,\mathbb{E}_2,\beta_1,\beta_2)\) [2210.01308], the equivalent two-springpot form \((\mathbb{V},\mathbb{G},\alpha,\beta)\) [2202.05878], a fractional Zener realization that is rheologically a fractional Maxwell branch in parallel with a spring [1212.4024], and modified Maxwell models in which the Caputo derivative is replaced by Prabhakar or Hadamard-type operators [1705.09246], [2204.02783].

## 1. Canonical meanings of “four-parameter FMM”

A standard and explicit four-parameter FMM is the series combination of two Scott–Blair elements. In the formulation used for fractional visco-elasto-plasticity, the constitutive equation is
\[
\sigma(t) + \frac{\mathbb{E}_2}{\mathbb{E}_1}\, \prescript{C}{0}{}\mathcal{D}_t^{\beta_2-\beta_1} \sigma(t)
=
\mathbb{E}_2\,\prescript{C}{0}{}\mathcal{D}_t^{\beta_2}\varepsilon(t),
\qquad
0<\beta_1<\beta_2<1,
\]
so the four parameters are \(\mathbb{E}_1,\mathbb{E}_2,\beta_1,\beta_2\) [2210.01308]. In the distributed-order Maxwell literature, the same structure is written as
\[
\sigma(t)+\frac{\mathbb{V}}{\mathbb{G}}\frac{d^{\alpha-\beta}\sigma(t)}{dt^{\alpha-\beta}}
=
\mathbb{V}\,\frac{d^{\alpha}\gamma(t)}{dt^{\alpha}},
\qquad
0<\beta\le \alpha<1,
\]
with four parameters \((\mathbb{V},\mathbb{G},\alpha,\beta)\); the classical Maxwell model is recovered at \(\alpha=1,\beta=0\) [2202.05878].

A second widespread meaning is the physically admissible four-parameter fractional Zener model,
\[
\sigma(t) + \tau_{\epsilon}^{\alpha}\frac{\partial^{\alpha}\sigma(t)}{\partial t^{\alpha}}
=
E_0\left[\epsilon(t)+\tau_{\sigma}^{\alpha}\frac{\partial^{\alpha}\epsilon(t)}{\partial t^{\alpha}}\right],
\]
with independent parameters \(E_0,\tau_\sigma,\tau_\epsilon,\alpha\) when \(\alpha=\beta\) [1212.4024]. In that representation the model is a fractional Maxwell element in parallel with an additional spring, and it is therefore a four-parameter fractional Maxwell-type structure in the spring–dashpot sense [1212.4024].

A third four-parameter interpretation appears in modified fractional Maxwell models with nonstandard kernels. The Hadamard-type model has constitutive law
\[
\sigma(t)+\widehat{O}_\nu^t \sigma(t)=\widehat{O}_\nu^t \epsilon(t), \qquad 0<\nu<1,
\]
and the physically meaningful parameter set is \((E,\eta_0,\theta,\nu)\), where \(E\) is the elastic modulus, \(\eta_0\) the initial viscosity, \(\theta\) the strain-hardening coefficient, and \(\nu\) the fractional order [2204.02783].

| Realization | Four parameters | Representative constitutive form |
|---|---|---|
| Two Scott–Blair elements in series | \(\mathbb{E}_1,\mathbb{E}_2,\beta_1,\beta_2\) | \(\sigma+\frac{\mathbb{E}_2}{\mathbb{E}_1}D_t^{\beta_2-\beta_1}\sigma=\mathbb{E}_2D_t^{\beta_2}\varepsilon\) |
| Two-springpot Maxwell form | \(\mathbb{V},\mathbb{G},\alpha,\beta\) | \(\sigma+\frac{\mathbb{V}}{\mathbb{G}}d_t^{\alpha-\beta}\sigma=\mathbb{V}d_t^\alpha\gamma\) |
| Fractional Zener / Maxwell + spring | \(E_0,\tau_\sigma,\tau_\epsilon,\alpha\) | \(\sigma+\tau_\epsilon^\alpha \partial_t^\alpha\sigma=E_0(\epsilon+\tau_\sigma^\alpha\partial_t^\alpha\epsilon)\) |
| Modified Hadamard-type Maxwell | \(E,\eta_0,\theta,\nu\) | \(\sigma+\widehat{O}_\nu^t\sigma=\widehat{O}_\nu^t\epsilon\) |

The coexistence of these formulations shows that “four-parameter FMM” is not a single fixed equation but a family of closely related Maxwell-type fractional models distinguished by how the four parameters are assigned to elastic scales, viscous scales, and fractional orders.

## 2. Fractional elements, series composition, and operator structure

The elementary building block of many FMMs is the Scott–Blair or spring-pot element. In the McKinley-type notation used for dual cross-linked hydrogels, a single element satisfies
\[
\sigma(t)=G\,\tau^{\nu}\,{}^{C}D_t^{\nu}\gamma(t),
\]
with \(0\le \nu \le 1\); \(\nu=0\) gives a Hookean spring and \(\nu=1\) gives a dashpot with viscosity \(\eta=G\tau\) [2509.02219]. The hydrogels paper builds the FMM by connecting two such elements in series, imposing common stress and additive strains,
\[
\sigma_1(t)=\sigma_2(t)=\sigma(t), \qquad \gamma(t)=\gamma_1(t)+\gamma_2(t),
\]
with exponents ordered as \(0\le \beta<\alpha\le 1\) [2509.02219].

The same series rule underlies the four-parameter FM model in the return-mapping framework. Each element obeys \(\sigma_i(t)=\mathbb{E}_i\,\prescript{C}{0}{}\mathcal{D}_t^{\beta_i}\varepsilon_i(t)\), and elimination of the internal strains produces a single fractional differential equation for the total strain and common stress [2210.01308]. This construction is the direct fractional analogue of the classical Maxwell series connection, except that the discrete spring/dashpot dichotomy is replaced by a continuum of intermediate rheological behaviors indexed by the orders \(\beta_i\).

In the Hadamard-type modification, the operator itself carries time dependence. The integer-order precursor is
\[
\sigma + \frac{\eta_0+\theta t}{E}\frac{d\sigma}{dt}
=
(\eta_0+\theta t)\frac{d\epsilon}{dt},
\]
with linearly varying viscosity \(\eta(t)=\eta_0+\theta t\) [2204.02783]. Its fractional generalization replaces \(\frac{\eta_0+\theta t}{E}\frac{d}{dt}\) by a Hadamard-type fractional power \(\widehat{O}_\nu^t\), which is a Caputo-type derivative with respect to the function
\[
\psi(t)=\frac{E}{\theta}\ln(\eta_0+\theta t)
\]
and therefore implements memory in logarithmic time [2204.02783]. This changes the model from a fixed-kernel Caputo FMM to one in which memory and time-dependent viscosity are intertwined.

## 3. Relaxation, creep, and complex modulus

For the two-springpot FMM, the characteristic material functions are explicit. In normalized form, the relaxation modulus is
\[
G(t)=t^{-\beta_2}E_{\beta_1-\beta_2,\,1-\beta_2}\!\big(-t^{\beta_1-\beta_2}\big),
\]
and the creep compliance is
\[
J(t)=\frac{t^{\beta_1}}{\Gamma(1+\beta_1)}+\frac{t^{\beta_2}}{\Gamma(1+\beta_2)}.
\]
In the frequency domain, the normalized complex compliance is
\[
J^*(\omega)=(i\omega)^{-\beta_1}+(i\omega)^{-\beta_2},
\]
so the model exhibits two power-law regimes in both time and frequency [2202.05878].

The FM model of the return-mapping paper gives the same structural result in a non-normalized notation:
\[
G^{FM}(t)
=
\mathbb{E}_1\,t^{-\beta_1}
E_{\beta_2-\beta_1,\,1-\beta_1}\!\left(
-\frac{\mathbb{E}_1}{\mathbb{E}_2}t^{\beta_2-\beta_1}
\right).
\]
The paper notes that this yields stretched-exponential behavior at short times and power-law behavior at long times [2210.01308].

In the hydrogel parametrization with shared modulus and time scale, the creep compliance is
\[
J(t)=\frac{1}{G}\left[
1+\frac{1}{\Gamma(1+\alpha)}\left(\frac{t}{\tau}\right)^\alpha
+\frac{1}{\Gamma(1+\beta)}\left(\frac{t}{\tau}\right)^\beta
\right],
\]
while the relaxation modulus is
\[
G(t)=G(\alpha-\beta)\,
E_{1-\alpha,\;1-\alpha+\beta}
\!\left[-\left(\frac{t}{\tau}\right)^{1-\alpha+\beta}\right].
\]
This representation is explicitly used to fit creep, stress relaxation, and oscillatory shear with the four parameters \(G,\tau,\alpha,\beta\) [2509.02219].

A one-order Maxwell form, often used as a target or limiting case, is
\[
\sigma(t)+\tau_0^\alpha D_t^\alpha \sigma(t)=b\,D_t^\alpha\varepsilon(t),
\qquad \alpha\in(0,1),
\]
which under constant strain gives
\[
\sigma(t)=\sigma_0\,E_\alpha\big[-(t/\tau_0)^\alpha\big],\qquad
G(t)=G_0\,E_\alpha\big[-(t/\tau_0)^\alpha\big].
\]
Its asymptotics are
\[
G(t)\approx
G_0\left[1-\frac{1}{\Gamma(\alpha+1)}\left(\frac{t}{\tau_0}\right)^\alpha\right]
\quad (t\ll \tau_0),
\]
and
\[
G(t)\sim
\frac{G_0}{\Gamma(1-\alpha)}
\left(\frac{t}{\tau_0}\right)^{-\alpha}
\quad (t\gg \tau_0),
\]
so the decay is algebraic and the model has no finite equilibrium modulus in this pure Maxwell form [2402.04969].

The modified Hadamard-type model changes the asymptotics qualitatively. Under constant strain, in normalized units \((E=\eta_0=\theta=1)\),
\[
\sigma(t)=E_{\nu,1}\big(-\ln^\nu(1+t)\big),
\]
and
\[
E_{\nu,1}\big(-\ln^\nu(1+t)\big)
\sim
\frac{\ln^{-\nu}(1+t)}{\Gamma(1-\nu)},
\qquad t\to\infty.
\]
For \(\nu=1\), the relaxation becomes
\[
\sigma(t)=e^{-\ln(1+t)}=(1+t)^{-1}.
\]
The paper explicitly analyzes relaxation and states that creep compliance is not derived [2204.02783].

## 4. Generalizations and unifying formulations

The four-parameter FMM sits inside several larger fractional-calculus frameworks. A particularly broad one is the Maxwell–Prabhakar model,
\[
\sigma(t)+a\,{}^{C}\!D_{\alpha,\beta,\omega}^{\gamma}\sigma(t)
=
b\,{}^{C}\!D_{\alpha,\beta,\omega}^{\gamma}\varepsilon(t),
\]
with kernel parameters \(\alpha,\beta,\gamma,\omega\) and mechanical parameters \(a,b\) [1705.09246]. In this formulation, the classical fractional Maxwell model is recovered either by setting \(\gamma=0\) or by setting \(\omega=0\), since both reductions collapse the Prabhakar kernel to the Caputo one [1705.09246]. The same paper also shows formal equivalence, in Laplace-domain creep compliance, to fractional Voigt and fractional Zener models for specific parameter sets, making the Maxwell–Prabhakar equation a unifying super-structure rather than a single named FMM [1705.09246].

Another extension is the generalized distributed-order Maxwell model. Its constitutive equation in Laplace form is
\[
\tilde{\sigma}(s)=\frac{b_1(s)b_2(s)}{b_1(s)+b_2(s)}\,\tilde{\gamma}(s),
\qquad
b_i(s)=\int_0^1 c_i(\alpha)s^\alpha\,d\alpha.
\]
If the weighting functions are delta distributions,
\[
c_1(\alpha)=\delta(\alpha-\beta_1),\qquad
c_2(\alpha)=\delta(\alpha-\beta_2),
\]
then \(b_1(s)=s^{\beta_1}\), \(b_2(s)=s^{\beta_2}\), and the generalized model reduces exactly to the two-springpot FMM [2202.05878]. The distributed-order formulation therefore replaces the two discrete fractional orders of the FMM by two continuous order spectra.

The Hadamard-type model provides a different generalization strategy. Instead of broadening the order distribution, it changes the time variable itself through
\[
s=\frac{E}{\theta}\ln(\eta_0+\theta t),
\]
so that in the transformed variable the model behaves like a standard Caputo-type fractional Maxwell equation, while in physical time it produces ultra-slow, logarithmic relaxation [2204.02783]. This suggests that some nonstandard four-parameter FMMs are best understood not as additional spring–dashpot branches but as Maxwell equations written in a non-Euclidean time scale.

## 5. Admissibility, network interpretation, and wave propagation

Physical admissibility is a central issue for four-parameter fractional Maxwell-type models. For the fractional Zener form,
\[
\sigma + \tau_\epsilon^\alpha \partial_t^\alpha \sigma
=
E_0\left(\epsilon+\tau_\sigma^\alpha \partial_t^\alpha \epsilon\right),
\]
the admissible regime is
\[
\alpha=\beta,\qquad E_0\ge 0,\qquad \tau_\epsilon^\alpha\ge \tau_\sigma^\alpha>0,
\]
which ensures thermodynamic constraints and monotonic relaxation [1212.4024]. In the frequency domain its complex modulus is
\[
G^*(\omega)
=
E_0\,\frac{1+(\tau_\sigma i\omega)^\alpha}{1+(\tau_\epsilon i\omega)^\alpha}
=
E_\infty + \frac{E_1(i\omega\tau_M)^\alpha}{1+(i\omega\tau_M)^\alpha},
\]
with
\[
E_\infty = E_0\left(\frac{\tau_\sigma}{\tau_\epsilon}\right)^\alpha,\quad
E_1 = E_0-E_\infty,\quad
\tau_M=\tau_\epsilon.
\]
This is precisely the modulus of a fractional Maxwell branch in parallel with a spring [1212.4024].

The same paper also connects the model to a Maxwell–Wiechert continuum. For the admissible case \(\alpha=\beta\), the generalized compressibility is equivalent to a continuous distribution of Maxwell elements with spectrum
\[
\kappa_{\nu\text{ML}}(\Omega)
=
\frac{1}{\pi}\,
\frac{\big(\tau_\sigma^\alpha-\tau_\epsilon^\alpha\big)\Omega^{\alpha-1}\sin(\alpha\pi)}
{\Omega^{2\alpha}+2\Omega^\alpha\cos(\alpha\pi)+1},
\]
which has three asymptotic power-law regions [1212.4024]. In wave propagation this yields three attenuation regimes,
\[
\alpha_k(\omega)\propto
\omega^{1+\alpha},\qquad
\omega^{1-\alpha/2},\qquad
\omega^{1-\alpha},
\]
at low, intermediate, and high frequencies, respectively [1212.4024]. This is one reason the four-parameter Maxwell/Zener family is used in acoustics, elastography, and medical applications [1212.4024].

A different admissibility result appears in the nonlinear Rational Extended Thermodynamics embedding of the fractional Maxwell law. There, exact coincidence with the fractional Maxwell relaxation
\[
\sigma(t)=k_0\,E_\alpha\big[-(t/\tau_0)^\alpha\big]
\]
is obtained for a special viscous energy \(e^{(V)}\), constant strain \(\varepsilon(t)=\varepsilon_0\), and initial condition \(\sigma(0)=k_0\) [2402.04969]. The same work shows that bounded viscous energy requires \(\alpha\in(1/2,1)\), while the broader compatibility range is \(\alpha\in(1/2,1]\), with \(\alpha=1\) giving the classical Maxwell model [2402.04969]. This locates the FMM inside a nonlinear, thermodynamically structured theory rather than treating it solely as an empirical kernel.

## 6. Computation, fitting, applications, and unresolved points

For numerical implementation, the return-mapping framework for fractional visco-elasto-plasticity discretizes the Caputo operators with the implicit L1 scheme and derives a discrete constitutive update for the FM model,
\[
\sigma_{n+1}
=
\frac{
C_1^M\left[\varepsilon_{n+1}-\varepsilon_n+\mathcal{H}^{\beta_2}\varepsilon\right]
+
C_2^M\left[\sigma_n-\mathcal{H}^{\beta_2-\beta_1}\sigma\right]
}{
1+C_2^M
},
\]
with model-dependent projection term
\[
C_{RM}^{ve}=\frac{C_1^M}{1+C_2^M}
\]
in the plastic correction phase [2210.01308]. The framework is fully implicit for linear viscoelastic models, semi-implicit for the quasi-linear case, and the reported numerical behavior is “at least first-order accurate for general loading conditions,” with a “reduction of \(50\%\) in CPU time” relative to existing approaches in the visco-plastic range [2210.01308].

The hydrogel study provides a direct materials application of a four-parameter FMM. The model is fitted simultaneously to stress relaxation, creep, and oscillatory shear using the parameters \(G,\tau,\alpha,\beta\), with creep and stress relaxation fitted over \(1\) to \(100\) s and oscillatory shear over \(0.1\) to \(100\) rad/s [2509.02219]. The cost function is a weighted sum of logarithmic least-squares residuals with
\[
w_{SR}=w_C=w_{FS}=0.333,
\]
and optimization is performed in Excel using the Generalized Reduced Gradient method [2509.02219]. For the PMA–Fe\(^{3+}\) dual-crosslinked hydrogels, \(\log G\) increases from \(4.58\) to \(5.33\) as salt concentration rises from \(0\) to \(1\) M, \(\log \tau\) decreases from about \(5.3\) to about \(3.6\), \(\alpha\) is approximately \(0.31\)–\(0.38\) at low/intermediate salt and about \(0.68\) at \(1\) M, and \(\beta\) remains in the range \(0.11\)–\(0.24\) [2509.02219].

Several unresolved points recur across the literature. In the hydrogel application, “The physical meaning of the two Scott Blair elements is, at present, somewhat obscure,” and the suggested mapping to covalent and ionic subnetworks is explicitly described as tentative [2509.02219]. In the Maxwell–Prabhakar model, a complete characterization of physically acceptable parameter sets is “left for future work” [1705.09246]. In the Hadamard-type model, only relaxation is worked out explicitly, and no creep formula is written [2204.02783]. In the wave-propagation literature, the microscopic physical origins of fractional behavior are described as still debated, even though the Maxwell–Wiechert equivalence gives a concrete interpretation in terms of many relaxation processes with a power-law distribution of relaxation times [1212.4024].

Taken together, these results place the four-parameter FMM at the center of modern fractional rheology: it is simultaneously a compact constitutive model, a special case of broader operator families, a reduced representation of distributed Maxwell networks, and a practical fitting model for materials with broad relaxation spectra. Its precise mathematical form varies across subfields, but the persistent features are a Maxwell backbone, four independent constitutive parameters, and non-exponential memory encoded by fractional calculus.

Source: https://www.emergentmind.com/topics/four-parameter-fractional-maxwell-model-fmm