dllrun / foam-extend-general-contact

---

Solid Contact Analysis Tool With Foam-Extend

1. Introduction

This code repository is based on the Finite Volume Method and it is forked from " https://github.com/Unofficial-Extend-Project-Mirror/foam-extend-foam-extend-4.0 ". The improved code library " solidModel " is an extension of the tools to predict deformation due to contact in solids. This includes the new tool " solidGeneralContact " for contact analysis - to allow the simulation of multiple solid bodies (e.g. particles) in contact.

It also provides the coupling of thermal contact boundary condition. Please note that this is only limited to one-to-one contact and no multiple solid bodies are allowed currently. Furthermore, the " elasticThermalSolidFoam " solver was improved to make it functional with the " solidContact " library.

---

2. Documentation

The documention below provide guidelines on how to install this foam-extend code repository. Please follow the link to the instruction for installation given below:

• Installation

Also, there is a link to the tutorial cases that have been developed for the purpose of demonstrating mechanical and thermal contact behavior of solids. Please follow the link to the example cases given below:

• Example Cases

---

4. How to cite this work?

If you are using this repository for your research, especially the library " solidModel " and the solver " elasticThermalSolidFoam ", then please cite it as follows.

  author      = {Dhakal, Ranjan},
  title       = {Solid contact analysis tool with Foam-extend},
  year        = {2022},
  institution = {Institute of Process and Particle Engineering, Graz University of Technology},
  url         = {https://github.com/dllrun/foam-extend-general-contact.git}

---

5. License

This toolkit is released under the GNU General Public License (version 3). More details can be found in the "LICENSE" file.

---

6. Disclaimer

This offering is not approved or endorsed by OpenCFD Limited, producer and distributor of the OpenFOAM software via www.openfoam.com, and owner of the OPENFOAM® and OpenCFD® trade marks.

---

7. Acknowledgement

This contribution is part of the Innovative Training Network (MSCA-ITN) “MATHEGRAM” Project. The project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 812638

---

Copyright Notice

• Copyright 2022 - Graz University of Technology (R. Dhakal, S. Radl).

• Version: foam-extend-4.0

• Licence: GPL-3.0 licence

• Last commit: June 2022

foam-extend-general-contact's Issues

Do not understand BLABLA BLA

this line is not clear to me:
https://github.com/dllrun/foam-extend-general-contact/blob/master/src/solidModels/fvPatchFields/solidContact/solidGeneralContactFvPatchVectorField.C#L2020

Please check if the current version compiles correctly

The current version is (commit 25b9e5a)
2. Please check if the code compiles correctly after pulling this commited version, (or adding functions until:

foam-extend-general-contact/src/solidModels/fvPatchFields/solidContact/solidGeneralContactFvPatchVectorField.C

Line 375 in 25b9e5a

void Foam::solidGeneralContactFvPatchVectorField::calcShadowZoneNames() const

3. And if it compiles, what about when you compile the code from the installation directory (i.e. entirely), is it possible to run a test case (of "solidMechanics" even without using this boundary condition)? Do you see the error "/bin/tools/RunFunctions: line 52: elasticSolidFoam: command not found"

"activeContactPairs" has not been declared in the header file

undeclared in .H file

foam-extend-general-contact/src/solidModels/fvPatchFields/solidContact/solidORGGeneralContactFvPatchVectorField.C

Line 1600 in 101679c

boolList activeContactPairs(shadowPatchNames().size(), false);

primitivePatchInterpolation has no member slaveToMaster

newGgiStandAlonePatchInterpolation class cannot found, so we used the primitivePatchInterpolation instead, however the masterToSlave function is missing.
Is it possible to replace "newGgiStandAlonePatchInterpolation" by some other class definition (from our repository) for the "zoneToZone" object (Line 511 and 517 of solidGeneralContactFvPatchVectorField.H) ?

Error message while compilation

When you compile the code using the version:
commit 0d7ba4a (HEAD -> master)
Author: Stefan Radl (as root) [email protected]
Date: Fri Jul 17 12:28:28 2020 +0200

1. Do you see "undefined reference to `Foam::zeroGradientFvPatchField::typeName" message? What is the error message it is refering to?

how is this function called

not sure how this function is called from normalContactModel.C

foam-extend-general-contact/src/solidModels/contactModels/normalContactModels/standardPenalty/standardPenalty.C

Line 143 in 101679c

void standardPenalty::correct

pointer "localSlavePtr_" is not defined --> resolve this issue

the case

https://github.com/dllrun/matheGRAMRanjanCases/tree/master/slidingFrictionGeneralHexa_2ParticlesInVertical

now runs fine, however, it bumgs into this error check:
https://github.com/dllrun/foam-extend-general-contact/blob/master/src/solidModels/fvPatchFields/solidContact/solidGeneralContactFvPatchVectorField.C#L2624

please resolve this

Number of Argument in correct member function of normalContactModel and frictionContactModel

When I try to compile the code (line 1913 of solidGeneralContactFvPatchVectorField.C) with
frictionModel(shadowI).correct
(
normalModel(shadowI).slavePressure(),
shadowPatchFaceNormals,
normalModel(shadowI).areaInContact(),
shadowPatchDD,
patchDDInterpToShadowPatch
);
An error: patchDDInterpToShadowPatch was not declared in this scope appears.
How to ensure that the correct argument is passed?

how is masterPatch_ of newGGIInterpolate.C set & why does masterPatch_.size() calculates the incorrect size of the master?

According to this error message -
https://github.com/dllrun/matheGRAMRanjanCases/blob/6d9749c4b29b31bfd4daffb0b905b0792849724a/slidingFrictionGeneralHexa_2ParticlesInVertical/log.slidingFrictionFoam#L3366 ,
it was directed from

foam-extend-general-contact/src/solidModels/numerics/newGGIInterpolation/newGGIInterpolate.C

Line 175 in 9f6b09d

"newGGIInterpolation::masterToSlave(const Field<Type> ff)"

.
So I checked if the patches size are correct and it was alright. But I couldn't understand how does this "masterPatch_" pointer

foam-extend-general-contact/src/solidModels/numerics/newGGIInterpolation/newGGIInterpolationTemplate.H

Line 148 in 9f6b09d

const MasterPatch& masterPatch_;

calculate the correct master patch?