---
title: Dynamic Resistive Force Theory (DRFT)
url: https://www.emergentmind.com/topics/dynamic-resistive-force-theory-drft
type: topic
---

# Dynamic Resistive Force Theory (DRFT)

Searching arXiv for recent and foundational papers on Dynamic Resistive Force Theory and closely related resistive-force frameworks.
Dynamic Resistive Force Theory (DRFT) is a reduced-order framework for predicting forces on bodies moving through granular media across quasistatic and strongly dynamic regimes. Its core construction inherits the local superposition structure of granular Resistive Force Theory (RFT): an intruder is decomposed into small surface elements, each element is assigned a local traction determined by depth, orientation, and motion direction, and the net force is obtained by integration over the intruder surface. DRFT extends that quasistatic picture by introducing dynamic corrections motivated by continuum frictional-plasticity with free-surface separation, especially an inertial term proportional to $\rho v^2$ and a structural correction associated with changes in the free surface around the intruder [2005.10976]. In the broader literature, the term is sometimes used narrowly for this continuum-informed dynamic granular model and sometimes more loosely for DRFT-like local superposition models embedded in robotic locomotion studies [1911.00079].

## 1. Concept and scope

In its most explicit formulation, DRFT is a granular intrusion model that preserves the locality and additivity of classical RFT while extending it into rate-dependent regimes [2005.10976]. The intended domain includes wheel locomotion, plate intrusion, and legged locomotion near free surfaces, where grains display solid-like and fluid-like behavior together with an ejected gas-like phase. The framework is therefore neither a grain-resolved simulation nor a full constitutive theory; it is a reduced-order force law designed to capture the dominant mechanisms that determine intrusion forces.

The quasistatic ancestor of DRFT is granular RFT. For a quasi-2D intruder surface subset $S$, that formulation writes the resistive force as
\[
(f_x,f_z)=\int_S \big(\alpha_x(\beta,\gamma),\ \alpha_z(\beta,\gamma)\big)\,H(z)\,z\,dS,
\]
where $\beta$ is local surface orientation, $\gamma$ is velocity direction, $H(z)$ suppresses resistance above the free surface, and the depth factor reflects gravity-scaled resistance [1510.02966]. DRFT retains this local decomposition but supplements it with dynamic corrections that arise when the momentum balance term $\rho \dot v_i$ becomes consequential even though the constitutive law remains rate-independent [2005.10976].

The literature summarized here also shows that the name is not used uniformly. A resistive force model for legged locomotion on granular media adopts the same segment-by-segment superposition logic and is described as being inspired by earlier resistive force models for sand-swimming, but it is not called DRFT in the paper [1911.00079]. A later MuJoCo implementation explicitly distinguishes its quasistatic 3D-RFT core from full Dynamic RFT because it omits the inertial term [2606.19504]. Outside granular media, related generalizations of RFT appear under different motivations, such as load-dependent coefficients for helical filaments and nonlocal operators converging to RFT in Stokes flow [2503.20520], [2604.09496].

## 2. Quasistatic foundations and continuum rationale

A central question in the DRFT literature is why a local force superposition should work at all in granular media. A continuum-mechanical answer is provided by a frictional-plasticity model with a Drucker–Prager/Mohr–Coulomb-type yield law, incompressible plastic flow during yielding, and an opening rule that sets stress to zero where density drops below a critical packing $\rho_c$ [1510.02966]. In dense regions,
\[
\sigma_{ij} = -P \delta_{ij} + 2\mu_c P \, D'_{ij}/\dot{\gamma},
\]
with
\[
D_{ij}=\frac{\partial v_i/\partial x_j+\partial v_j/\partial x_i}{2},\qquad
D'_{ij}=D_{ij}-\delta_{ij}D_{kk}/3,\qquad
\dot{\gamma}=\sqrt{2D'_{ij}D'_{ij}},
\]
while the open-state rule is
\[
\sigma_{ij}=0 \quad \text{if } \rho<\rho_c.
\]
The momentum balance is
\[
\frac{\partial \sigma_{ij}}{\partial x_j}+\rho g_i=\rho \dot{v}_i.
\]

Within this continuum model, the local RFT force law emerges naturally rather than empirically. For a granular “garden hoe” geometry, dimensional analysis yields
\[
\mathbf{F}_{\text{gran}}=\rho_c g L^3\,\mathbf{\Psi}_{\text{gran}}(\beta,\gamma,\mu_c,\mu_w),
\]
or, after absorbing material constants,
\[
\mathbf{F}_{\text{gran}}=L^3\,\mathbf{\Psi}_{\text{gran}}(\beta,\gamma).
\]
For a small plate element this implies
\[
(\alpha_x(\beta,\gamma),\alpha_z(\beta,\gamma))=\frac{2\mathbf{\Psi}_{\text{gran}}(\beta,\gamma)}{\sin\beta}.
\]
Integrating those local laws reproduces the same scaling and structure as the full continuum solution for that class of geometries [1510.02966].

This continuum derivation also explains why granular RFT is stronger than its viscous-fluid analogue. In Stokes flow, the force scales as
\[
\mathbf{F}_{\text{visc}}=\eta L v\,\mathbf{\Psi}_{\text{visc}}(\beta,\gamma),
\]
whereas a naive RFT construction over a patch of size $\lambda$ gives
\[
\mathbf{F}^{\text{RFT}}_{\text{visc}}=\frac{L^2}{\lambda}v\,\mathbf{\Psi}_{\text{visc}}(\beta,\gamma).
\]
The unavoidable dependence on the arbitrary length $\lambda$ weakens strong superposition in viscous flow. By contrast, granular frictional plasticity, hydrostatic pressure dependence, and the inability to sustain tension create a local, depth-scaled response that is especially favorable to RFT-style reduction [1510.02966].

A plausible implication is that DRFT should be expected to work best when those quasistatic premises remain meaningful locally even as macro-inertial effects alter the global intrusion response. That is precisely the regime formalized in the dynamic extension.

## 3. DRFT formulation and force decomposition

The 2020 dynamic intrusion paper identifies three force contributions that together define DRFT: a static quasistatic resistive term, a dynamic inertial correction, and a dynamic structural correction [2005.10976]. The static term is the classical local traction
\[
\mathbf{t} = \boldsymbol{\alpha}(\beta,\gamma)\,H(-z)|z|,
\]
where $\boldsymbol{\alpha}(\beta,\gamma)=(\alpha_x,\alpha_z)$ is the empirical traction-per-depth function and $z$ is height relative to the free surface.

The dynamic inertial correction is a normal traction
\[
-\mathbf{n}\lambda \rho v_n^2,
\]
where $\mathbf{n}$ is the outward normal, $v_n$ is the normal component of the surface velocity, $\rho$ is the effective bulk density, and $\lambda$ is an $O(1)$ fitting constant. This term is the momentum-flux contribution that emerges from the inertial term in the momentum equation.

The dynamic structural correction modifies the depth entering the quasistatic law:
\[
\tilde z = z+\delta h.
\]
The full DRFT traction law is then
\[
\mathbf{t}=\boldsymbol{\alpha}(\beta,\gamma)\,H(-\tilde z)|\tilde z|-\mathbf{n}\lambda \rho v_n^2.
\]
For an intruder surface $S$, the net force is
\[
\mathbf{F}=\int_S \left[\boldsymbol{\alpha}(\beta,\gamma)\,H(-\tilde z)|\tilde z|-\mathbf{n}\lambda \rho v_n^2\right] dA.
\]

| Contribution | Expression | Physical role |
|---|---|---|
| Static quasistatic term | $\boldsymbol{\alpha}(\beta,\gamma)\,H(-z)|z|$ | Depth- and orientation-dependent resistance |
| Dynamic inertial correction | $-\mathbf{n}\lambda \rho v_n^2$ | Momentum transfer to grains |
| Dynamic structural correction | $\tilde z=z+\delta h$ | Free-surface modification around the intruder |

The limiting cases are explicit. In the quasistatic limit, $v_n\to 0$ and $\delta h\to 0$, so DRFT reduces to classical RFT. If $\delta h\approx 0$, one obtains a static-plus-$\rho v^2$ law. If the free surface changes strongly with speed, the structural correction can dominate the rate effect [2005.10976].

The reduction from continuum mechanics to DRFT is physics-guided rather than a strict asymptotic derivation. The continuum analysis shows that if a slow quasistatic solution is scaled to a faster one with similar flow topology, then the inertial term implies an added interface stress proportional to $\rho V^2$. The structural correction is introduced because the quasistatic term depends on depth below the free surface; once the free surface itself is altered by intrusion, the local resistive force changes even if the constitutive law remains rate-independent [2005.10976].

## 4. Robotic locomotion and intrusion applications

DRFT and DRFT-like models have been used across wheels, articulated legs, bipedal feet, and full multibody robot simulators. The most direct dynamic application in the summarized literature is grousered wheel locomotion. In that setting, a purely quasistatic RFT captures low-speed behavior but fails at higher spin rates, where translation saturates, slip increases, and sinkage increases. Adding only the inertial $\rho v^2$ term improves the model little and can predict the wrong sinkage trend; the key missing mechanism is the structural correction, because high spin expels material behind the wheel, lowers the rear free surface, and reduces support in the trailing shear zone. For the wheel case, the free-surface drop is approximated as
\[
\delta h \approx r\left(\frac{r\omega^2}{g}\right),
\]
applied to the rear zone, and with $\lambda=1$ the resulting DRFT reproduces the observed velocity and sinkage trends well [2005.10976].

A closely related but explicitly quasistatic application is rigid-wheel interaction with dry granular media under forced-slip conditions. There, granular RFT is implemented as a local surface-stress superposition model in which each wheel surface element contributes according to depth, orientation, and movement direction, subject to a leading-edge hypothesis. The study concludes that, for the range of inputs considered, RFT can be reliably used to predict rigid wheel granular media interactions with accuracy exceeding that of traditional terramechanics methodology in several circumstances, while still remaining fundamentally quasistatic rather than a full DRFT law [1901.10667].

For legged locomotion on granular media, a resistive force model in the vertical plane applies the same local-superposition principle to complex limbs. The leg is divided into small segments, each segment is assigned a local force determined by depth $|z|$, angle of attack $\beta$, and angle of intrusion $\gamma$, and the net force is obtained by integration:
\[
F_{z,x} = \int \alpha_{z,x}(\beta_s,\gamma_s)\,|z_s|\, dA_s.
\]
In loosely packed poppy seeds of packing fraction about $0.58$, with a model aluminum plate of area $3.8 \times 2.5\ \text{cm}^2$ moved at $1\ \text{cm/s}$, both $f_z$ and $f_x$ were found to be approximately linear in depth. Applied to an L-leg and a reversed L-leg, the model predicts $F_z(\theta)$ and $F_x(\theta)$ without any fitting parameters, matching measured forces to within about $15\%$. Embedded into a multi-body dynamics simulation using MBDyn for a robot about $15\ \text{cm}$ long and $150\ \text{g}$ in mass, the simulated average forward speed matches experiments within $15\%$ for L-legs across the tested frequencies and within $30\%$ for reversed L-legs up to about $3\ \text{Hz}$ [1911.00079].

A more explicit DRFT extension appears in the foot-shape-dependent resistive force model for bipedal walkers on sand. That model adds both inertial drag and an effective intrusion depth correction. After discretizing the foot into plates, the dynamic contribution is
\[
d{F}_{i,v} = -\lambda_{v}\,\rho\, \|\langle \boldsymbol{v}, \boldsymbol{n}\rangle\|^2 \hat{\boldsymbol{v}}\, dS_i,
\]
while the corrected depth is
\[
\tilde{z}= z_0 + \lambda_{h}\sqrt{\|\boldsymbol{v}\|\, z_0\, g(\beta,\phi_s)},
\qquad
g(\beta,\phi_s)=\left(\frac{1}{\tan\beta}+\frac{1}{\tan\phi_s}\right)^{-1}.
\]
The full plate-level force is
\[
d {F}_i = f_{i,1}(\boldsymbol{v},{e}_1)\, d{F}_{i,1} + f_{i,23}(\boldsymbol{v},{e}_1)\, d\tilde{F}_{i,23} + d{F}_{i,v}.
\]
On fine sand with particle size $0.2$–$0.5\ \text{mm}$, calibrated values include $\zeta=2.06$, $\lambda_h=1.93$, and $\lambda_v=1.1$. Across slow, medium, and high gait conditions, the corrected model reduces RMSE relative to conventional 3D-RFT for $F_x$, $F_y$, and $F_z$, with especially large improvements at medium and high gait speeds. The same study reports that the elliptical foot is better for high-speed walking because it required less work than flat foot in fast gait conditions, whereas the flat foot is better for slow walking because it saved energy under slow motion [2403.03460].

A recent open-source simulator integrates 3D granular RFT into MuJoCo by replacing rigid-contact ground forces with a distributed granular contact model. Each body is discretized into plate elements, only leading surfaces contribute, and forces are applied at plate centroids using MuJoCo sites and `mj.applyFT`. The implementation uses
\[
\mathbf{F}_j=\mathbf{F}_1+\mathbf{F}_2+\mathbf{F}_3
\]
with depth-, orientation-, and motion-direction-dependent empirical terms, and it introduces single exponential moving average smoothing,
\[
\mathbf{F}_{sn} = \alpha_F \mathbf{F}_0 + (1-\alpha_F)\mathbf{F}_{s,n-1},
\qquad \alpha_F=0.1,
\]
to stabilize low-speed dynamics. The paper explicitly states that it did not implement the inertial term of Dynamic RFT because simulated speeds were below about $12.6\ \text{cm/s}$ tip speed and below the cited threshold of about $20\ \text{cm/s}$. Even so, the resulting RFT-SiM predicts walking distance and foot sinkage of a 12-degree-of-freedom hexapod robot within $20\%$ of experiments in sand in the abstract, while the detailed validation reports walking-distance underestimation of about $9.6\%$ to $34\%$ and sinkage overestimation by about $20\%$ [2606.19504].

## 5. Internal torque, body loading, and generalized resistive-force models

Although DRFT is most often used for net ground-reaction prediction, the same local-force integration logic also determines internal loading. In steady undulatory locomotion through resistive-force-dominated media, the body is prescribed to undergo a traveling curvature wave
\[
\kappa(s,t)=A\xi \sin\!\big[2\pi(\xi s+t)\big],
\]
with wave speed
\[
v_\kappa=\frac{1}{\xi}.
\]
Local forces satisfy
\[
F_\perp = C_\perp v_\perp,\qquad F_\parallel = C_\parallel v_\parallel,
\]
with $C_\perp=2$ and $C_\parallel=1$ for a thin cylinder in viscous drag, and the internal torque at body position $s$ is obtained by integrating external resistive loads over the anterior body segment:
\[
T(s)=-\int_s^1 [\mathbf r(l)-\mathbf r(s)]\times \mathbf F(l)\,dl.
\]
In the small-amplitude approximation,
\[
T(s)=-\int_s^1 (l-s)\,F_y\,dl.
\]
The main result is that torque typically forms a traveling wave for $0.5\lesssim \xi <1.8$, but the speed ratio $v_T/v_\kappa$ decreases from about $5.2$ to $1.3$ as $\xi$ increases from $0.5$ to $1.8$. Near $\xi\approx 2$, the pattern transitions to a two-wave-like structure, and for $\xi>2$ more complex patterns appear because phase-shifted force contributions cancel in the distance-weighted integral. The paper explicitly frames this as a resistive-force-theory study of internal torque rather than a new force law, but it shows how distributed resistive loads reorganize internal actuation demands [1707.00492].

In low-Reynolds-number hydrodynamics, generalized resistive-force models introduce different kinds of corrections. For rigid helical filaments with an attached spherical load of radius $R$, the local RFT law is written as
\[
\mathbf w(t_0) = \alpha_\perp \mathbf f_\perp(t_0) + \alpha_\parallel \mathbf f_\parallel(t_0),
\]
but the optimal coefficients are shown to depend on the load through
\[
r=\frac{3R}{4\pi n b}.
\]
The resulting coefficients $\alpha_\perp(R)$ and $\alpha_\parallel(R)$ are load dependent rather than universal. For a normal polymorph, using Lighthill’s free-helix coefficients yields only about $3.1\%$ error in the in-plane flow components $(u,v)$ but about $111.3\%$ error in the axial component $w$, and for larger loads the swimming-speed error reaches nearly $50\%$. The parallel resistance coefficient changes by up to about $20\%$ as the load grows [2503.20520]. This is a generalized RFT, but the paper itself states that it is not “dynamic” in the time-dependent-memory sense.

A distinct mathematical generalization appears in the nonlocal curve evolution for an immersed elastic filament in a 3D Stokes fluid. There, the RFT operator
\[
\mathcal L_{\epsilon,\mathrm{RFT}}(X)=\frac{\log\epsilon}{4\pi}(I+X_s\otimes X_s)
\]
is compared with a pseudodifferential operator $\overline{\mathcal L}_\epsilon$ built from the exact Fourier symbol of the straight-cylinder slender-body Neumann-to-Dirichlet map. The resulting evolution
\[
\frac{\partial X}{\partial t} = -\overline{\mathcal L}_\epsilon(X)\big[(X_{sss}-\tau X_s)_s\big], \qquad X_s^2=1,
\]
is globally well posed in the natural energy space, and after the rescaled time $\underline t=(\log\epsilon)t$ the solutions converge to RFT dynamics as $\epsilon\to 0$ with
\[
\|X-Y\|_{L^\infty_t H^2_s} + (\log\epsilon)^{1/2}\|X-Y\|_{L^2_t\dot H^{7/2}_s} \le (\log\epsilon)^{-1/2} C_{\rm RFT}.
\]
This establishes, in a rigorous Stokes-flow setting, how local resistive-force dynamics emerge from a more detailed nonlocal model [2604.09496].

## 6. Validity regime, limitations, and nomenclature

DRFT is best understood as a reduced-order model whose reliability depends on the physical regime and on what is meant by “dynamic.” In the granular continuum foundation, the quasistatic RFT reduction is strongest when the material response is local, frictional, depth-scaled, and effectively decoupled from global geometry. When rate dependence or size effects become important, the local resistive plots acquire explicit dependence on additional groups such as
\[
I_G=\frac{v^2 d^2}{gL^2}, \qquad \frac{d}{L},
\]
and the simple universal RFT law no longer survives unchanged [1510.02966].

The 2020 DRFT formulation is also explicit about scope. It assumes additivity of static and dynamic terms, treats the structural correction $\delta h$ as an effective representation of free-surface change, and is intended as a bridge between continuum simulation, empirical quasistatic RFT, and higher-speed $\rho v^2$ intrusion laws. It is not a replacement for continuum simulation, is best suited to shallow intrusions where RFT assumptions still hold, and has been validated primarily for dry, non-cohesive granular media [2005.10976].

A frequent misconception is that any granular resistive-force model embedded in a dynamic simulator is automatically DRFT. Several papers draw a sharper distinction. The MuJoCo implementation states that it uses the quasistatic 3D-RFT core only and explicitly did not implement the inertial term because the motions were below the regime where inertial sand effects become significant. The same paper also lists limitations including no grain-scale effects, no surface deformation, no inertial correction, sensitive calibration, and computational cost about an order of magnitude slower than default MuJoCo contacts due to many force application points [2606.19504]. Similarly, the legged-locomotion model on granular media is DRFT-like in construction but is presented as a resistive force model rather than a named DRFT formalism [1911.00079].

Another nomenclature issue is that related generalized RFT models in fluids need not be “dynamic” in the granular DRFT sense. The load-dependent helical-filament theory modifies local coefficients through global force balance with an attached load rather than through inertial or free-surface corrections [2503.20520]. Conversely, the undulatory-locomotion torque study belongs to the RFT family because torque emerges from integrating local external loads, but its key contribution is an explanation of internal torque-wave formation and transition rather than a new resistive-force law [1707.00492].

Taken together, these studies define DRFT not as a single universally fixed formula, but as a family of local-additive reduced models whose most developed form in granular media consists of quasistatic RFT plus a momentum-flux correction and a free-surface structural correction. The framework’s central promise is computational economy without abandoning the geometry-aware surface integration that makes resistive-force methods useful across locomotion and intrusion problems [2005.10976].

Source: https://www.emergentmind.com/topics/dynamic-resistive-force-theory-drft