Dear all, I am trying to implement plasticity in FEniCS. I find an extremely valuable tutorial here (Courtesy @bleyerj ). Thanks a ton @bleyerj . I tried to match the result of the FEniCS code with the Abaqus so that I can build upon that. But unfortunately, the displacement at the end of the simul…

Hi, thanks for the comment on the plasticity example. I haven’t compared to Abaqus but I get identical results between the FeniCS implementation and another completely independent implementation of FE plasticity. I am not an experienced Abaqus user but your time goes up to 1., what is your imposed…

Hi, try this for plotting via matplotlib import matplotlib.pyplot as plt from matplotlib import cm p = plot(p_avg, cmap=cm.jet) plt.colorbar(p) plt.show() this works for me [pavg]

Hi @bhaveshshrimali Thanks a ton ! It is working now. How to change the limits of the colorbar in such plots ? I haven’t used matplotlib much before and very new to it. Thanks again.

Matplotlib has an extensive documentation with plenty of examples, see https://matplotlib.org/3.2.1/api/_as_gen/matplotlib.pyplot.colorbar.html

Thank you for the direction ! @dokken

Hi, I modified the Elasto-plastic demo to solve a thermomechanical problem in a staggered manner. I need to reset the equivalent plastic strain, p=0, once the temperature exceeds a certain value(e.g. 1514 Kelvin ) I have tried using an Expression that depends on the temperature field T # therma…

Thank you for your reply. because the Function spaces of the thermal and mechanical problems are different, I thought of defining a new parameter key using the quadrature space function W0 and then locally project the Expression using the local_project function. At the key should represent T_ann …

The code does not Need to solve the whole problem, only Give the error. Thus, if the product consist of two functions from different spaces, just define the two functions and do the assign operation .

Hi, If p and key are both defined on the quadrature function space W0 you can perform the multiplication on the array of dofs directly: p.vector().set_local(p.vector().get_local()*key.vector().get_local())

Thank you so much. This works now.

Plasticity in fenics

dokken May 24, 2020, 3:27pm 25

Please supply a minimal example as described at Read before posting: How do I get my question answered? as it for instance is not clear what kind of variable key is.

Topic		Replies	Views
Newton solver converge initially, then diverge and blow up Legacy DOLFIN	8	564	January 15, 2024
Problems encountered in solving thermoelastic problems General	4	311	May 24, 2024
Incremental implementation of hyperelasticity variational formulation	13	272	October 8, 2024
Implementing dynamic analysis for von-Mises plasticity General	5	207	October 30, 2024
Problem implementing effective stress function method with FEniCS variational formulation	4	1402	August 27, 2021

Plasticity in fenics

Related topics