Hello everyone, I’m trying to solve a very simple problem from linear elasticity, where a square mesh of 2x2 elements has a different shear modulus at each element (the shear modulus is piecewise constant, DG degree 0). Then apply a shear to the mesh, and compute the resulting strain. My problem is…

There are two things that you are ignoring. The first one is that even if you would like to have a linear stress, you do not want a continuous lagrange representation. Therefore you need to use write_checkpoint. Using a function space that is linear discontinuous you obtain the attached image. T =…

When using CG degree 1, the result corresponds to degree 2

dokken March 30, 2020, 6:15am 2

To visualize other elements than CG 1 in paraview, you Need to use XDMFFile.write_checkpoint, see: Loading xdmf data back in

Topic		Replies	Views
Expression for higher order element	1	679	October 18, 2019
How to apply quadratic element type	8	1679	March 25, 2020
Inaccurate solution in linear elasticity simulation with volume forces and quadrature elements(FEniCS 2019.1.0) Legacy DOLFIN mesh , boundary-conditions	0	29	March 31, 2025
Why use DG space to project stress/strain？ mathematics	6	1516	July 14, 2020
Varying material properties based on data points General	0	445	February 8, 2023

When using CG degree 1, the result corresponds to degree 2

Related topics