Mesh deformation resulting in intersecting boundaries

Numerical methods and mathematical models of Elmer
Post Reply
skaal
Posts: 6
Joined: 20 Nov 2023, 18:41
Antispam: Yes

Mesh deformation resulting in intersecting boundaries

Post by skaal »

Hi all,

I have the following rather weird behavior:
I set up a simulation with a rather thin and soft material. When solving for the deformation, I see the upper wall protruding through the lower one, which seems a bit weird, physically ;) :
protruding_walls.png
(132.55 KiB) Not downloaded yet
I cast this into a very simple and reduced test case, that shows exactly this behavior in the second time step. Solver input and mesh file are attached.

Am I am doing something wrong here? Did I find a bug? Is this a genuine feature?
Attachments
testcase.sif
(2.88 KiB) Downloaded 20 times
plate.unv
(62.45 KiB) Downloaded 18 times
Rich_B
Posts: 423
Joined: 24 Aug 2009, 20:18

Re: Mesh deformation resulting in intersecting boundaries

Post by Rich_B »

Hello,

Attached is a modified sif file, adding a slowly increasing load over time. This improves the behavior somewhat, up until the solver stops early.

Rich.
Attachments
testcase2.sif
(2.92 KiB) Downloaded 18 times
kevinarden
Posts: 2315
Joined: 25 Jan 2019, 01:28
Antispam: Yes

Re: Mesh deformation resulting in intersecting boundaries

Post by kevinarden »

It is possible in the FEA simulation to shove one side of the plate through it's thickness and out the other side. This is not possible in real word physics. However, you only gave the FEA model linear elastic material properties, and therefore it doesn't know that it should be shearing and failing. You basically defined F=Kx, Hooke's law. Force = stiffness times displacement. If you want any thing else to happen during the simulation than you would have to define many more physical terms in the FEA simulation. It is possible to simulate geometric and material non-linear behavior. However the code solved exactly what you modeled.

The result looks exactly like it should given the boundary conditions, load, and solver used.
skaal
Posts: 6
Joined: 20 Nov 2023, 18:41
Antispam: Yes

Re: Mesh deformation resulting in intersecting boundaries

Post by skaal »

Thanks!

It seems that this works quite a bit better when running with a nonlinear / plastic material model.
@kevinarden: your UMAT library works pretty well, here :)

However, also revisiting the time step helped quite a lot -- the case is significantly more stable with a larger time step, even when using the linear material model as in the provided sif file.
Post Reply