Hyper Elastic demo code, output remains the same cube

sachin_goyal · March 26, 2019, 9:33am

I ran the demo_hyperelastic code on https://fenicsproject.org/docs/dolfin/1.4.0/python/demo/documented/hyperelasticity/python/documentation.html
The code runs successfully, but when I use ParaView to open the output pvd file, It is just a normal standard cube without any desired deformation. I simply ran the given demo code with degree as 10. Would be great if someone can help me out.

Rudy · March 26, 2019, 1:39pm

Depending on the specific version of Paraview, go to Filters > Warp By Vector and click apply.

marcomusy · March 26, 2019, 1:47pm

… or try out this:

github.com

marcomusy/vtkplotter/blob/master/examples/other/dolfin/ex06_elasticity2.py

"""
Show fenics mesh and displacement solution.

Refer to original script for the details:
https://fenicsproject.org/docs/dolfin/2018.1.0/python/
       demos/hyperelasticity/demo_hyperelasticity.py.html
"""
print(__doc__)

########################################################### dolfin
from dolfin import *

# Create mesh and define function space
mesh = UnitCubeMesh(12, 12, 12)
V = VectorFunctionSpace(mesh, "Lagrange", 1)

# Mark boundary subdomains
left  = CompiledSubDomain("near(x[0], side) && on_boundary", side=0.0)
right = CompiledSubDomain("near(x[0], side) && on_boundary", side=1.0)

This file has been truncated. show original

Topic		Replies	Views
Using the output instead of saving it in .pvd mesh	3	1384	May 10, 2020
Hyper elasticity tutorial in 2D and 3D comparison General	10	199	January 30, 2024
Different results for the same mesh created differently (Hyperelasticity) mesh	6	1353	August 17, 2020
Hyperelastic model with two different materials General	1	127	February 4, 2024
Visualize deformations of a hyperbolic elasticity system in ParaView I/O	0	278	July 7, 2022

Hyper Elastic demo code, output remains the same cube

Related Topics