Can FEniCS do some part elastic and other part of the body elastoplastic?
Thanks for suggestions.
FEniCS is a generalised tool for approximation of the solution of PDEs by the finite element method, So yes, why not?
1 Like
I think as long as you can translate your weak form into UFL the existing machinery should allow you to do that. See, maybe this and the replies in that discussion to maybe get started.
1 Like
Hi,
you can adapt the following demo https://comet-fenics.readthedocs.io/en/latest/demo/2D_plasticity/vonMises_plasticity.py.html to your need by arbitrarily setting a very large value of the yield limit in your region which should remain elastic.
Alternatively, you can also split the form into two parts, one elastic and one elastoplastic on the corresponding subdomains. Note however, that without using the mixed-dimensional branch, you will still have to define the various Quadrature functions throughout the whole domain.
The alternative way is what @kamensky suggested in this topic
Yes,
to be more precise, the only difference for plasticity compared to hyperelasticity of @kamensky’s answer is that internal state variables are required. Hence, they must be defined either throughout the whole domain and set to zero in the elastic part or defined only on the corresponding subdomain but this requires the mixed-dimensional branch.
2 Likes