I would like to be able to calculate and visualize the traction vector field for a given surface in a problem with hyperelasticity and Dirichlet boundary conditions. Below is sample code for an example problem in which a displacement is enforced on one side of a cube while the other is held in place…

You should be able to do something along the lines of: Iterate connected mesh entities in parallel

Calculating the traction vector field for a surface

dokken November 27, 2020, 2:27pm 6

Note that to visualize DG functions, you need to use XDMFFile.write_checkpoint, see: Loading xdmf data back in

Topic		Replies	Views
Storing traction and displacement in each iteration	2	1220	December 18, 2021
How to formulate Dirichlet boundary problem with traction free surfaces on the rest of the boundaries?	1	520	March 10, 2021
Retrieving imposed traction boundary in a linear elastic problem	4	581	October 7, 2021
Traction Boundary Conditions in 2D Mechanics Problems General	3	233	April 24, 2024
Non-trivial boundary conditions for linear elasticity General	3	55	November 5, 2024