[Xulia]

Hello !
I am running a solid contact problem using solids4foam . It is a cylinder inserted on a box and I apply a load on the top lid of the cylinder. The box is fixed by two patches as shown in the picture (I also highlight the load in the picture)

setup_cylinder.png


I am using the unsLinGeometry model and the problem is time independent. After running for 5k iterations, I obtain a residual of 6.e-05. If I visualize the results and then use the paraview "warp by vector" filter, the displacement of the box is as expected: shifting downwards. However, the cylinder pops up from the box, as if I were pulling from it rather than pushing:

warped.jpg



My D file is as follows:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object D;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 0 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
box_0_main_0
{
type solidTraction;
traction uniform ( 0 0 0 );
pressure uniform 0;
value uniform (0 0 0);
}
box_0_f_0
{
type fixedDisplacement;
value uniform (0 0 0);
}
box_0_f_1
{
type fixedDisplacement;
value uniform (0 0 0);
}
box_0_cylinder_0
{
type solidContact;
shadowPatch cylinder_0_box_0;
master no;
value uniform (0 0 0);
}
cylinder_0_main_0
{
type solidTraction;
traction uniform ( 0 0 0 );
pressure uniform 0;
value uniform (0 0 0);
}
cylinder_0_load_0
{
type solidTraction;
traction uniform ( 0 0 0 );
pressure uniform 319883;
value uniform (0 0 0);
}
cylinder_0_box_0
{
type solidContact;
master yes;
rigidMaster no;
shadowPatch box_0_cylinder_0;
normalContactModel standardPenalty;
standardPenaltyNormalModelDict
{
relaxationFactor 0.1;
penaltyScale 1;
}
frictionContactModel standardPenalty;
standardPenaltyFrictionModelDict
{
relaxationFactor 0.1;
penaltyScale 1;
frictionLaw coulomb;
frictionLawDict
{
frictionCoeff 0.5;
}
}
value uniform (0 0 0);
}
}


// ************************************************** *********************** //


I am happy to share the full case if it helps.

Is this the expected behavior ? I was hoping to see the full geometry going downwards.

any suggestions/ ideas are very welcome.