White box / Physical formulations / Contact and friction

Regularized Coulomb friction

Physical formulation · assumptions, essential equations, full derivation, algorithm, code and verification references.

Browse another sheet Physical formulations / Contact and friction
Physical formulationsContact and frictionSheet 4 / 5
Physical formulationPRIOR NARROW ACCEPTANCE · current-source requalification not established

The tangential trial force is projected onto a Coulomb bound set by compressive normal pressure.

01 Model and scope

The tangential trial force is compared with the Coulomb disk whose radius is μp. Inside the disk the point sticks; outside it the force is projected to the boundary and the point slips. The implemented regularization and frozen projection/normal constrain where this idealized law applies.

Route assumptions and calculation checks

  1. The documented route uses small transformations, fixed initial projection/normal, direct static solution, and stabilized translational DOFs.
  2. Stick/slip state, tangential force, friction bound, and dissipation diagnostics are exposed.
  3. Mesh sensitivity remains material; this is not general finite sliding, impact, self-contact, or frictional dynamics.

02 Essential equations and symbols

qtrial=kts;qmax=μpqtrial = kt s; qmax = \mu{} p∣∣qtrial∣∣≤qmax→stick,q=qtrial||qtrial|| \le{} qmax \to{} stick, q=qtrial∣∣qtrial∣∣>qmax→slip,q=qmaxqtrial/∣∣qtrial∣∣||qtrial|| > qmax \to{} slip, q=qmax qtrial/||qtrial||

Definitions and conventions

ss
Relative tangential slip in the contact frame.
kₜ
Tangential regularization stiffness.
μp\mu{}p
Coulomb force bound from friction coefficient and compressive pressure.
qq
Tangential contact force.
Regularized Coulomb friction · method mapNormal pressure sets the admissible tangential force; the trial traction is classified as stick or slip.
4 detailed steps

03 Step by step derivation

Follow the assumptions, equations and discrete operations in order. The full explanation stays visible.

1 · Separate normal pressure and tangential slip

First determine the compressive normal pressure from the declared contact route. In the tangent plane, measure relative slip since the last accepted state and form an elastic tangential trial force. The Coulomb limit is proportional to normal pressure, so friction capacity vanishes when the pair separates.

p≥0,qmax=μpp\ge0,\qquad q_{max}=\mu pqtr=kt str,∥qtr∥≤qmax  ⇒  stick\mathbf q^{tr}=k_t\,\mathbf s^{tr},\qquad \|\mathbf q^{tr}\|\le q_{max}\;\Rightarrow\;\text{stick}

2 · Project the trial traction onto the friction bound

If the trial traction lies outside the admissible disk, return it to the boundary in the direction of the trial vector. This is the discrete stick/slip projection. The implementation's regularization and tolerances determine how transitions are smoothed near the switching surface.

q={qtr,∥qtr∥≤μp,μp qtr∥qtr∥,∥qtr∥>μp.\mathbf q=\begin{cases}\mathbf q^{tr},&\|\mathbf q^{tr}\|\le\mu p,\\\mu p\,\dfrac{\mathbf q^{tr}}{\|\mathbf q^{tr}\|},&\|\mathbf q^{tr}\|>\mu p.\end{cases}

3 · Assemble force, tangent, and state

Apply equal-and-opposite tangential forces to the two sides of the pair and add the route-specific tangential tangent to equilibrium. Persist the stick/slip state and accepted tangential history only after the global increment converges; otherwise restore the previous accepted state.

ftc=CtTq,R=fext−fint−fnc−ftc\mathbf f_t^c=\mathbf C_t^T\mathbf q,\qquad \mathbf R=\mathbf f_{ext}-\mathbf f_{int}-\mathbf f_n^c-\mathbf f_t^cΔWdiss=−q⋅Δs≥0for a dissipative slip convention\Delta W_{diss}=-\mathbf q\cdot\Delta\mathbf s\ge0\quad\text{for a dissipative slip convention}

4 · Read the bounded claim correctly

Inspect pressure, tangential force, friction bound, stick/slip flag, and dissipation together. The documented case uses small transformations, a fixed initial projection and normal, a direct static solve, and stabilized translational DOFs. Mesh and regularization sensitivity remain material to the outcome.

∥q∥≤μp+τf\|\mathbf q\|\le\mu p+\tau_f

This record does not establish general finite sliding, impact, self-contact, or frictional dynamics. Use the cited contract for the exact regularization parameter and acceptance thresholds.

04 Algorithm and implementation

Read the linked versioned source for the actual assembly and solver behavior. The derivation above explains the mathematics and is not a verbatim source excerpt.

05 Verification and limitations

PRIOR NARROW ACCEPTANCE · current-source requalification not established

Each reference supports its own model, observable and analysis route. Read the original contract and limitations before transferring a conclusion.

06 Sources and related sheets

Continue the reasoning

Back to sheet contents ↑