SOLVERS

UNDER DEVELOPMENT

We have developed a specialized, arc-length based Lagrangian FE solver with physically accurate surface interaction modeling accelerated using k-d Trees (Roy and Sundaram, IJNME 2023). An implicit dynamic version of this has also been developed and is currently being MPI accelerated.

STABLE

VISCOCONFCON is a fast Singular Integral / Integro-Differential Equation (SIE / SIDE) based solver for geometrically conforming partial slip frictional contacts in linear viscoelasticity.

It supports arbitrary pin-loads, including cyclic fretting-type loads, bulk plate-stresses, linear viscoelastic solid materials with arbitrarily complex viscous networks (multiple Prony terms) and a feature-rich postprocessor, including rapid edge-of-contact and subsurface stress calculations.

It has the ability to perform dissipation calculations and has been extensively validated with other numerical prodedures (FEA). The front-end and high-level functions in Viscoconfcon are written in MATLAB, with MEX/C++ optimized routines for core and kernel functions.

Download: VISCOCONFCON is now available for download as a zip file. It also contains my old elastic / elastic and rigid / elastic conforming contact codes. Please note that this code is provided as-is. When you download it, you do it with the understanding that it comes with NO guarantees of support or any assumption of responsibility on my part. Please DO NOT modify, publish, or redistribute it in any form without my consent.

Any version of MATLAB 2014a or later should be adequate should execute the code. Please note that the precompiled MEX .mexw32 files should work in both 32- and 64- bit Windows envirnoments, but will require recompilation with an appropriate C++ compiler on non-Windows systems. Please contact me if you want the core files for linux.

Input specifications are in u_input_viscoconfcon.m. Call visco_confcon_new to run the solver.

For background, please see the paper: Satish K. Dayalan, Narayan K. Sundaram, "Partial slip contact of a rigid pin and a linear viscoelastic plate", International Journal of Solids and Structures, 2016, 100-101, 319-331