Hertz - msolids/musen GitHub Wiki

Contact vector

$CP = O_1+ \frac{O_2-O_1}{r_1+r_2}\cdot r_1$

$\bar r_c = CP - O_1$

$\bar r_n= \frac{\bar r_c}{|\bar r_c|}$

$\bar{v}_{rel} = \bar{v}_2 - \bar{v}_1 + \bar{\omega}_1 \times \bar{r}_c - \bar{\omega}_2 \times \bar{r}_{c_2}$

$\bar{v}_{rel,n}=\bar{r}_n\cdot(\bar{r}_n\cdot\bar{v}_{rel})$

$\bar{v}_{rel,t}=\bar{v}_{rel}-\bar{v}_{rel,n}$

$R^* = \frac{r_1 \cdot r_2}{r_1 + r_2}$

$E^*=\left(\frac{1-\nu^2_1}{E_1}+\frac{1-\nu^2_2}{E_2}\right)^{-1}$

$\xi_n = r_1 + r_2 - |O_2 - O_1|$

$k_n=2E^*\sqrt{\xi_{n}\cdot R^*}$

$\bar F_n = - \bar r_n \cdot \frac{2}{3} \cdot \xi_n \cdot k_n$

$\Delta\bar\xi_t = \bar v_{rel,t} \cdot \Delta t$

$k_t=8\cdot G^*\cdot\sqrt{R^*\cdot\xi_n}$

$\Delta \bar{F}_t = [ k_t \cdot \Delta \bar{\xi}_t ]$

$\bar{F}_{t,pr}^{cor} = \bar{F}_{t,pr}-\bar{r}_n \cdot (\bar{r}_n \cdot \bar{F}_{t,pr} )$

$\bar{F}_{t,pr}^{cor} = \bar{F}_{t,pr}^{cor} \cdot |\bar F_{t,pr}|\, / \,|\bar F_{t,pr}^{cor} |$

$\bar F_t = \bar F_{t,pr}^{cor} + \Delta \bar F_t$

if $|\bar F_t| > \mu_{sl}\cdot |\bar F_n |$

then $\bar F_t = \mu_{sl} \cdot |\bar F_n | \cdot \frac{\bar F_t}{|\bar F_t|}$

$\bar M_{ro,1} = -\mu_{ro} \cdot |\bar F_n| \cdot r_1 \cdot \frac{\bar\omega_1}{|\bar\omega_1|}$

$\bar M_{ro,2} = -\mu_{ro} \cdot |\bar F_n| \cdot r_2 \cdot \frac{\bar\omega_2}{|\bar\omega_2|}$

$\bar F_{tot} = \bar F_n + \bar F_t$

$\bar F_{1} = \bar F_n + \bar F_t$

$\bar F_{2} = -\bar F_n + \bar F_t$

$\bar M_{tot,1} = \bar r_n \times \bar F_t \cdot r_1 + \bar M_{ro,1}$

$\bar M_{tot,2} = -\bar r_n \times \bar F_t \cdot r_2 + \bar M_{ro,2}$

Hertz H. (1882). Über die Berührung fester elastischer Körper. Journal die reine und angewandte Mathematik, 92, 156-171.

Symbol	Description
CP	Contact point [m]
$E^*$	Equivalent Young’s modulus [Pa]
$E_1,\,E_2$	Young’s moduli of contact partners [Pa]
$\bar F_n,\, \bar F_t$	Force in normal and tangential directions [N]
$\bar F_{t,pr}$	Tangential force on previous iteration [N]
$\bar F_{ro}$	Force due to the rolling friction [N]
$k_t$	Spring constant [N/m]
$\bar M_{ro}$	Moment due to the rolling friction [N]
$O_1,\,O_2$	Centers of contact partners [m]
$R^*$	Equivalent radius [m]
$r_1,\,r_2$	Particle radii [m]
$\bar r_c$	Contact vector [m]
$\bar r_n$	Normalized contact vector [-]
$\bar v_{rel}$	Relative velocity [m/s]
$\bar v_1,\, \bar v_2$	Velocities of contact partners [m/s]
$\Delta\bar\xi_t$	Increment of tangential displacement in the current step [m]
$\mu_{ro},\,\mu_{sl}$	Coefficient of rolling friction and sliding friction [-]
$\xi_n$	Normal overlap [m]
$\omega_1,\,\omega_2$	Rotation velocities of particles [rad/s]