git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@14620 f3b2605a-c512-4ea7-a41b-209d697bcdaa

This commit is contained in:
sjplimp
2016-02-15 18:37:06 +00:00
parent 531cd5fe2c
commit c43da51b44

View File

@ -1,123 +0,0 @@
.. index:: pair_style body/rounded/polygon
pair_style body/rounded/polygon command
=======================================
Syntax
""""""
.. parsed-literal::
pair_style body/rounded/polygon cutoff_global cutoff_inner delta_ua c_n
* cutoff_global = global cutoff for COM interactions (distance units)
* cutoff_inner = cutoff for vertex-vertex and vertex-edge interactions (distance units)
* delta_ua = characteristic contact length (distance units)
* c_n = normal friction coefficient (energy time/distance^2)
Examples
""""""""
.. parsed-literal::
pair_style body/rounded/polygon 6.0 0.5 0.5 0.1
pair_coeff 1 1 1.0 100.0 12.0
Description
"""""""""""
Style *body/rounded/polygon* implements the pairwise body/body
interactions between 2D convex polygons, as described
in :ref:`Fraige <Fraige>`. This pair style is designed for use with
the "rounded/polygon" body style, which is specified as
an argument to the "atom-style body/rounded/polygon" command.
See the :doc:`body <body>` doc page for more details about the body
styles LAMMPS supports.
The "rounded/polygon" style treats a body particle as a discrete
element model (DEM) polygon composed of N vertices.
The coordinates of a body particle are its center-of-mass (COM).
The initial position of the vertices are specified the data file.
The current implementation follows the multiple contact points model as
described in :ref:`Fraige <Fraige>`, but neglects and the contact history
the tangential contact forces for now.
The parameters *delta_ua* and *c_n* are the characteristic contact length
and the normal friction coefficient, respectively. See Eqs. (1) and (5)
in :ref:`Fraige <Fraige>` for more details.
For style *body/rounded/polygon*, the following coefficients
must be defined for each pair of atoms types via the :doc:`pair_coeff <pair_coeff>`
command as in the examples above:
* sigma = vertex rounded diameter (distance units)
* k_n = normal repulsion strength (energy/distance^2)
* k_na = normal attraction strength (energy/distance^2)
* cutoff = cutoff for COM interactions (distance units)
The last parameter is optional. If not specified, the global
cutoff specified in the pair_style command is used.
Note that the current implementation assumes that the edge length
of the polygons (L) should be at least twice greater than
the vertex rounded radius (sigma) so as to distinguish vertex-vertex from
vertex-edge contacts.
.. note::
Because this pair style requires the body velocities to compute
the friction component of the forces, users should specify
*comm_modify vel yes* in the input script with this pair style.
----------
**Mixing, shift, table, tail correction, restart, rRESPA info**:
For atom type pairs I,J and I != J, the epsilon and sigma coefficients
and cutoff distance for all of this pair style can be mixed. The
default mix value is *geometric*. See the "pair_modify" command for
details.
This pair style does not support the :doc:`pair_modify <pair_modify>`
shift, table, and tail options.
This pair style does not write its information to :doc:`binary restart files <restart>`.
This pair style can only be used via the *pair* keyword of the
:doc:`run_style respa <run_style>` command. It does not support the
*inner*, *middle*, *outer* keywords.
----------
Restrictions
""""""""""""
This style is part of the BODY package. It is only enabled if LAMMPS
was built with that package. See the :ref:`Making LAMMPS <2_3>` section for more info.
Defining particles to be bodies so they participate in body/body or
body/particle interactions requires the use of the :doc:`atom_style body <atom_style>` command.
Related commands
""""""""""""""""
:doc:`pair_coeff <pair_coeff>`
**Default:** none
.. _Fraige:
**(Fraige)** F. Y. Fraige, P. A. Langston, A. J. Matchett, J. Dodds,
Particuology, 6, 455 (2008).
.. _lws: http://lammps.sandia.gov
.. _ld: Manual.html
.. _lc: Section_commands.html#comm