linear discrete element method

Community: MathModDB

mathematical model

Available identifiers

MaRDI QIDQ6675389

variant of the discrete element method that employs linear contact models for inter-particle forces

List of contained entities

Coulomb friction condition between two particles	$\begin{matrix} ‖ 𝑭_{i j}^{T, c o n s} ‖ & \leq μ ‖ 𝑭_{i j}^{N} ‖ with \\ 𝑭_{i j}^{T, c o n s} & = - k_{i j}^{T} 𝝃_{i j} = - k_{i j}^{T} (𝒙_{C_{j i}} - 𝒙_{C_{i j}} - ⟨ 𝒙_{C_{j i}} - 𝒙_{C_{i j}}, 𝒏_{i j} ⟩ 𝒏_{i j}) \end{matrix}$
	$𝑭_{i j}^{N}$ represents normal interaction force of two particles
	$𝒏_{i j} \in ℝ^{3}$ represents unit normal vector
	$μ$ represents coefficient of friction
	$k_{i j}^{T}$ represents tangential stiffness
	$x_{C_{i j}}$ represents contact point of particles

normal interaction force of two particles	$\begin{matrix} 𝑭_{i j}^{N} & = (k_{i j}^{N} δ_{i j} + d_{i j}^{N} {\dot{δ}}_{i j}) 𝒏_{i j} \\ δ_{i j} & = ⟨ 𝒙_{i} - 𝒙_{j}, 𝒏_{i j} ⟩ \\ {\dot{δ}}_{i j} & = ⟨ 𝒗_{i} - 𝒗_{j}, 𝒏_{i j} ⟩ \\ r_{i j} & = (r_{i} + r_{j}) / 2 \\ k_{i j}^{N} & = E_{N} π r_{i j} / 2 \\ d_{i j}^{N} & = D_{N} 2 \sqrt{k_{i j}^{N} m_{i j}} \\ m_{i j} & = \frac{m_{i} m_{j}}{m_{i} + m_{j}} \\ 𝒏_{i j} & = \frac{𝒙_{i} - 𝒙_{j}}{‖ 𝒙_{i} - 𝒙_{j} ‖} \end{matrix}$
	$D_{N}$ represents damping coefficient
	$E_{N}$ represents Young's modulus
	$𝑭_{i j}^{N}$ represents interaction force
	$𝒗_{i} \in ℝ^{3}$ represents particle velocity
	$r_{i j}$ represents particle radius
	$𝒏_{i j} \in ℝ^{3}$ represents unit normal vector
	$𝒙_{i} \in ℝ^{3}$ represents particle position
	$k_{i j}^{N}$ represents normal stiffness
	$m_{i}$ represents particle mass

tangential interaction force of two particles	$\begin{matrix} 𝑭_{i j}^{T, p r e} & = & - k_{i j}^{T} 𝝃_{i j} - d_{i j}^{T} {\dot{𝝃}}_{i j} \\ 𝐢𝐟 ‖ 𝑭_{i j}^{T, c o n s} ‖ & \leq & μ ‖ 𝑭_{i j}^{N} ‖ 𝐭𝐡𝐞𝐧 𝑭_{i j}^{T} = 𝑭_{i j}^{T, p r e}, \\ 𝐞𝐥𝐬𝐞 𝑭_{i j}^{T} & = & μ ‖ 𝑭_{i j}^{N} ‖ 𝒕, \\ 𝐰𝐡𝐞𝐫𝐞: \\ 𝑭_{i j}^{T, c o n s} & = & - k_{i j}^{T} 𝝃_{i j}, \\ 𝝃_{i j} & = & 𝒙_{C_{j i}} - 𝒙_{C_{i j}}, \\ {\dot{𝝃}}_{i j} & = & ⟨ 𝒗_{i} - 𝒗_{j}, 𝒕 ⟩ 𝒕, \\ 𝒕 & = & \frac{𝝃_{i j}}{‖ 𝝃_{i j} ‖}, \\ k_{i j}^{T} & = & k_{i j}^{N} / 1.2, \\ k_{i j}^{N} & = & \frac{π}{2} \frac{r_{i} + r_{j}}{2} E_{N}, \\ d_{i j}^{T} & = & 2 D_{T} \sqrt{k_{i j}^{N} \frac{m_{i} m_{j}}{m_{i} + m_{j}}} . \end{matrix}$
	$𝑭_{i j}^{T} \in ℝ^{3}$ represents tangential interaction force of two particles
	$𝒕$ represents unit tangent vector
	$𝒗_{i} \in ℝ^{3}$ represents particle velocity
	$𝒙_{C_{i j}} \in ℝ^{3}$ represents contact point of particles
	$D_{T}$ represents damping coefficient
	$E_{N}$ represents Young's modulus
	$k_{i j}^{T}$ represents tangential stiffness
	$k_{i j}^{N}$ represents normal stiffness
	$m_{i}$ represents particle mass
	$r_{i}$ represents particle radius

torque of particle	$\begin{matrix} 𝑻_{i} & = (𝒙_{a_{i j}} - 𝒙_{i}) \times 𝑭_{i j}^{T} \\ 𝒙_{a_{i j}} & = 𝒙_{i} + \frac{r_{i}}{r_{i} + r_{j}} (𝒙_{i} - 𝒙_{j}) \end{matrix}$
	$𝒙_{i} \in ℝ^{3}$ represents particle position
	$𝐅_{i j}^{T}$ represents tangential interaction force of two particles
	$r_{i}$ represents particle radius
	$T$ represents torque

Further items linking to linear discrete element method

Item	Property
discrete element method	specialized by
recurrent neural network surrogate for discrete element method	contains

This page was built for model: linear discrete element method