fix multicontact/halfspace command¶
Syntax¶
fix ID group multicontact/halfspace geometric_prefactor gp_value
- ID, group are documented in fix command
- multicontact/halfspace = style name of this fix command
- geometric_prefactor/gp_value = optional keyword-value pair
geometric_prefactor gp_value = gamma gamma is an empirical factor accounting for the geometry
Examples¶
fix mc all multicontact/halfspace
fix mc all multicontact/halfspace geometric_prefactor 1.8
Description¶
Implements the multicontact model by (Brodu et al.) . This model computes a per-contact deformation for each particle based on the other contacts this particles has. A particle i with contact ij has the following new radius (when computing the contact laws with particle j):
r_i + sum_k delta_{ij->ik}
where r_i is the default radius of particle i, sum_k is the sum over all particles in contact with i (k!=*j*). The delta value is given by
delta_{ij->ik} = - gamma (1 + nu) F_ik / (2 pi Y d_{ik->ij}) * [(**n**_ik . **u**_{ik->ij})(**n**_ij . **u**_{ik->ij}) + (3 - 4 nu) **n**_ik . **n**_ij - (1 - 2 nu) (**n**_ik + **u**_{ik->ij}) . **n**_ij) / (1 + **n**_ik **u**_{ik->ij}) ]
where
- gamma = geometric prefactor
- nu = Poisson ratio
- F_ik = absolute value of normal force acting at contact ik
- Y = Youngs modulus
- d_{ik->ij} = distance from contact ik to contact ij
- **u**_{ik->ij} = unit vector pointing from contact ik to contact ij
- **n**_ik = normal vector of contact ik
- **n**_ij = normal vector of contact ij
Restart, fix_modify, output, run start/stop, minimize info¶
No information about this fix is written to binary restart files.
fix_modify cannot be used on the parameter of this fix.
Restrictions¶
Requires the use of the gran surface model multicontact