[juanblanco]

Hello,

I ran the tutorial plateHole with solidDisplacementFoam but I modified the initial conditions and set the traction boundary condition for displacement to zero, i.e. no forces acting on the plate. Also, I turned on the thermal stresses calculations with a constant initial temperature with adiabatic boundary conditions.

From displacement equations I would expect that with a zero gradient temperature field the displacement will drop to zero, it did not (see figures attached). I analyzed the thermal gradient (I attach the figure) and, even though it is really small (~1e-13), it seems that affects the result from the displacement calculation. To test this, I turned off the thermal stress calculation and without the thermal term in the diplacement equation I do obtain a zero displacement field.

I need to run a transient coupled to this which will add a source term to the energy equation that eventually fades away so I will eventually arrive to the problem state before.

To sum up, the problem is that when approaching the stationary an incredibly small gradient modifies appreciably the displacement field.

Observations:
*Increasing outer loops and/or increasing number of energy equation solving steps did not work out.
*A steady-state ddt scheme has problems with the energy equation, floating point sends an error.

Thank you for your time. If anybody knows a solution it will be welcomed.