fix mesh/surface command

fix mesh/surface/planar command

Syntax:

fix ID group-ID mesh/surface file filename premesh_keywords premesh_values mesh_keywords mesh_values surface_keywords surface_values
fix ID group-ID mesh/surface/planar file filename premesh_keywords premesh_values mesh_keywords mesh_values surface_keywords surface_values 

• ID, is documented in fix command, the group-ID is ignored for this command.
• mesh/surface or mesh/surface/planar = style name of this fix command
• file = obligatory keyword
• filename = name of STL or VTK file containing the triangle mesh data
• zero or more premesh_keywords/premesh_value pairs may be appended
• premesh_keyword = type or precision or heal or verbose
• type value = atom type (material type) of the wall imported from the STL file precision value = length mesh nodes this far away at maximum will be recognized as identical (length units) heal value = auto_remove_duplicates or no verbose value = yes or no zero or more mesh_keywords/mesh_value pairs may be appended
• mesh_keyword = scale or move or rotate or temperature

    scale value = factor
      factor = factor to scale the mesh in x-, y-, and z-direction (double value)
    move values = mx my mz
      mx my mz = move the mesh by this extent in x-, y-, and z-direction (length units)
    rotate values = axis ax ay az angle ang
      axis = obligatory keyword 
      ax, ay, az = axis vector for rotation (length units) 
      angle = obligatory keyword 
      ang = angle to rotate the geometry around the specified axis (in degrees) 
    temperature value = T0
      T0 = Temperature of the wall (temperature units) 
  

• zero or more surface_keywords/surface_value pairs may be appended
• surface_keyword = surface_vel or surface_ang_vel or curvature

    surface_vel values = vx vy vz
      vx vy vz = conveyor belt surface velocity  (velocity units)
    surface_ang_vel values = origin ox oy oz axis ax ay az omega om
      origin = mandatory keyword
      ox oy oz = origin of rotation (length units)
      axis = mandatory keyword
      ax ay az = axis vector for rotation (length units) 
      omega = mandatory keyword
      om = rotaional velocity around specifyied axis (rad/time units)
    curvature value = c
      c = maximum angle between mesh faces belonging to the same surface (in degree) 
  

Examples:

fix cad all mesh/surface file mesh.stl type 1 

LIGGGHTS vs. LAMMPS Info:

This command is not available in LAMMPS.

Description:

This fix allows the import of triangual surfeace mesh wall geometry for granular simulations from ASCII STL files or legacy ASCII VTK files. Style mesh/surface is a general surface mesh, and mesh/surface/planar represents a planar mesh. mesh/surface/planar requires the mesh to consist of only 1 planar face.

Generall, you can apply scaling, offset and rotation to the imported mesh via keywords scale, move, rotate. Operations are applied in the order as they are specified. The group-ID is ignored for this command. For periodic boundaries, the mesh is mapped.

One fix represents one wall with a specific material, where the material is identified via keyword type. If multiple meshes with different materials are desired, the respective meshes must be imported with different fix mesh/surface commands.

With the temperature keyword, you can define a constant temperature for a mesh in conjunction with heat conduction via fix heat/gran. Note that the actual calculation of the heat transfer happens only if you use the mesh in conjunction with a granular wall, see fix wall/gran.

With the optional surface_vel keyword, you can specify the imported mesh as conveyor belt. The velocity direction for each mesh face is given by the projection of the conveyor belt velocity parallel to the mesh face, the velocity magnitude for each mesh face is equal to the conveyor belt velocity. This ensures the model makes sense also in case the mesh is curved. Likewise, the optional rotation model activated via keyword surface_ang_vel mimics rotational motion of the mesh (e.g. for modeling a shear cell)

The precision keyword specifies how far away mesh nodes can be at maximum to be recognized as identical.

If LIGGGHTS stalls because of duplicate elements, you can try setting heal to auto_remove_duplicates. LIGGGHTS will then try to heal the geometry by removing duplicate elements.

IMPORTANT NOTE: You should check the changes to the geometry, e.g. by using a dump mesh/stl command.

The curvature keyword lets you specify up to which angle between two triangles the triangles should be treated as belonging to the same surface (e.g. useful for bends). This angle is used to decide if (a) contact history is copied from one triangle to the other as the contact point proceeds and (b) if edge and corner interaction is calculated.

Quality checks / error and warning messages:

LIGGGHTS checks a couple of quality criteria upon loading a mesh. LIGGGHTS tries to give you as much information about the issue as possible.

Warning messages: