Mechanical Systems: Computing library for simulations in continuum and discrete mechanics

MechSys is a programming library for the implementation of simulation tools in mechanics. Its source code is mainly written in C++ with easier to use templates for further custumization.

Potential applications include, for instance, solid mechanics using the finite element method, mechanics of granular materials using the discrete element method, and both computational fluid dynamics and computational electromagnetics using the lattice Boltzmann method.

See: References below for papers produced with MechSys.

See: Mailing list for a forum for discussions and Q&A: MechSys' mailing list.

The Discrete Element Method in MechSys can handle particles of many different shapes.

MechSys can be used to develop simulation tools for continuum mechanical systems as well as for systems based on the concept of discrete/molecular dynamics.

Table of Contents


This library includes (among others):


The following list presents some Conference and Journal papers that have employed MechSys for calculations.

Acknowledgement of these papers is highly appreciated especially if you use MechSys in your work.

Behraftar S, Galindo-Torres SA, Scheuermann A, Williams DJ, Marques EAG, Janjani H (2015) Validation of a novel discrete-based model for fracturing of brittle materials. International Journal of Rock Mechanics and Mining Sciences (submitted) ,
Galindo-Torres SA, Scheuermann A, Puscasu R (2015) A Lattice Boltzmann solver for Maxwell equations in dielectric media. Journal of Geophysics and Engineering (submitted) ,
To H.D., Scheuermann A., Galindo-Torres SA (2015) Sequential sphere packing by trilateration equations. Granular Matter (submitted)
To H.D., Scheuermann A., Galindo-Torres SA (2014) Primary fabric fraction analysis of granular soils. Acta Geotechnica Special Issue: Hydraulic instabilities in granular materials Theory and Experiments , http://dx.doi.org/10.1007/s11440-014-0353-9
Galindo-Torres SA, Pedroso DM, Williams DJ, Muhlhaus HB (2014) An Analysis of the Strength of Anisotropic Granular Assemblies via Discrete Methods. Applied Mechanics and Materials , http://www.scientific.net/AMM.553.525
Galindo-Torres SA, Scheuermann A, Williams DJ, Muhlhaus HB (2014) Micro-Mechanics of Contact Erosion. Applied Mechanics and Materials , http://www.scientific.net/AMM.553.513
To H.D., Scheuermann A., Galindo-Torres SA (2014) A Numerical Approach for the Determination of the Primary Fabric of Granular Soils. Applied Mechanics and Materials , http://www.scientific.net/AMM.553.489
Galindo-Torres SA, Scheuermann A, Williams DJ, Muhlhaus HB (2013) A micro-mechanical approach for the study of contact erosion. Acta Geotechnica Special Issue: Hydraulic instabilities in granular materials Theory and Experiments , dx.doi.org/10.1007/s11440-013-0282-z
Galindo-Torres SA (2013) A coupled Discrete Element Lattice Boltzmann Method for the simulation of fluid-solid interaction with particles of general shapes Computer Methods in Applied Mechanics and Engineering , dx.doi.org/10.1016/j.cma.2013.06.004
Galindo-Torres SA, Pedroso DM, Williams DJ, Muhlhaus HB (2013) Strength of non-spherical particle assemblies under triaxial and shearing loading configurations. Granular Matter, dx.doi.org/10.1007/s10035-013-0428-6
Galindo-Torres SA, Scheuermann A, Li L, Pedroso DM, Williams DJ (2013) A lattice Boltzmann model for studying transient effects during imbibition-drainage cycles in unsaturated soils. Computer Physics Communications, 184(4), 1086-1093. dx.doi.org/10.1016/j.cpc.2012.11.015
Galindo-Torres SA, Scheuermann A, Li L (2012) A numerical study of the permeability in a tensorial form for laminar flow in anisotropic media. Physical Review E, 86, 046306. dx.doi.org/10.1103/PhysRevE.86.046306
Galindo-Torres SA, Pedroso DM, Williams DJ, Li L (2012) Breaking processes in three-dimensional bonded granular materials with general shapes. Computer Physics Communications, 183(2), 266-277. dx.doi.org/10.1016/j.cpc.2011.10.001
Pedroso DM, Galindo-Torres SA, Williams DJ (2010) On the mechanical behaviour of dry cohesionless soils by DEM simulations. Proc 5th International Workshop on New Frontiers in Computational Geotechnics (NFCOMGEO-V), July 26-27 (2010), Brisbane, Australia, ISBN:9784877971779
Galindo-Torres SA, Pedroso DM (2010) Molecular dynamics simulations of complex shaped particles using Voronoi-based spheropolyhedra. Physical Review E, 81, 061303. dx.doi.org/10.1103/PhysRevE.81.061303
Galindo-Torres SA, Marroquin F, Wang Y, Pedroso DM, Castano J (2009) Molecular dynamics simulation of complex particles in three dimensions and the study of friction due to nonconvexity. Physical Review E, 79, 060301. dx.doi.org/10.1103/PhysRevE.79.060301
Buzzi O, Pedroso DM, Giacomini A (2008) Caveats on the implementation of the Generalized Material Point Method, Computer Modeling in Engineering & Sciences, 31(2), 85-106. dx.doi.org/10.3970/cmes.2008.031.085

Copyright of the papers rests with a variety of different parties. They are made available solely on the same basis that the researcher would supply a single copy of a paper to an individual. It is up to you to ensure that your use of a paper does not infringe copyright law.


Some animations produced by MechSys' simulations are shown below. Additional animations can be found in the MechSysLib YouTube channel or in Sergio Torres' Youtube channel for simulations created from the DEM and the LBM modules.

LBM Simulations

An imbibition-drainage simulation within a porous medium

A cavitation simulation where a fluid forms a liquid bubble

DEM Simulations

A triaxial test with general shaped particles

A brazilian test to determine the tensile strenght of rocks

DEM+LBM Simulations

A DEM dolphin spinning and diving into an LBM pool

A DEM duck with lower density of liquid LBM fluid floating on it

EMLBM Simulations

A cross section area of a dipole antenna, both the electric (parallel to the plane) and the magnetic (perpendicular to the plane) are shown

An electric pulse (the electric field vector is going out from the screen) propagates between two media of different dielectric permittivity changing its speed and wavelength

Mailing List

Forum for discussions and Q&A: MechSys' mailing list.