Dear François Bouchut and FF++ users
I’d like to apply the code developed on the basis of the above discussion to a hexagon.
Assuming a linear elastic body, we consider a boundary condition where the right half of the regular hexagon is under tension with a strain of 1 (see attached image 1).
However, the result is as shown in the attached image 2, where a mesh inversion occurs near the corners of the hexagon.
It may be possible to suppress this if a boundary condition is given at this corner point where the vertical displacement is zero…
If you have any good ideas, please let me know.
250204_Hexagon.edp (4.2 KB)