LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Commands

### compute smd/tlsph_shape command

Syntax:

    compute ID group-ID smd/tlsph_shape

- ID, group-ID are documented in compute command
- smd/tlsph_shape = style name of this compute command

Examples:

    compute 1 all smd/tlsph_shape

Description:

Define a computation that outputs the current shape of the volume associated with a particle as a rotated ellipsoid. It is only meaningful for particles which interact according to the Total-Lagrangian SPH pair style.

See this PDF guide to use Smooth Mach Dynamics in LAMMPS.

Output info:

This compute calculates a per-particle vector of vectors, which can be accessed by any command that uses per-particle values from a compute as input. See How-to discussions, section 6.15 for an overview of LAMMPS output options.

The per-particle vector has 7 entries. The first three entries correspond to the lengths of the ellipsoid's axes and have units of length. These axis valus are computed as the contact radius times the xx, yy, or zz components of the Green-Lagrange strain tensor associated with the particle. The next 4 values are quaternions (order: q, x, y, z) which describe the spatial rotation of the particle relative to its initial state.

Restrictions:

This compute is part of the USER-SMD package. It is only enabled if LAMMPS was built with that package. See the Making LAMMPS section for more info.

This quantity will be computed only for particles which interact with the Total-Lagrangian SPH pair style.

Related commands:

smd/contact_radius

Default: none