Hi all,
I wanted to apply Robin boundary conditions, in the nearlyincompressible hyperelastic model (Mooney Rivlin), using the 3D cylinder mesh. The Robin boundary condition should be like:
S·n + a·u = 0,
where S is the stress, n is the the normal vector to the surface, u is the displacement, and a is a parameter representing the elastic response due to the presence of surrounding tissue (for example 1·10^2),. Robin BC are imposed on Gamma_R = Gamma, the surface of the cylinder.
How can I implement this?
Thanks a lot!
Hi Patricia,
On 12/05/2014 05:10 PM, patricia...@upf.edu wrote:
Hi all,
I wanted to apply Robin boundary conditions, in the nearlyincompressible hyperelastic model (Mooney Rivlin), using the 3D cylinder mesh. The Robin boundary condition should be like:
S·n + a·u = 0,
where S is the stress, n is the the normal vector to the surface, u is the displacement, and a is a parameter representing the elastic response due to the presence of surrounding tissue (for example 1·10^2),. Robin BC are imposed on Gamma_R = Gamma, the surface of the cylinder.
How can I implement this?
The term that would be required for the Robin BC is not implemented yet. You would need something like \int_\Gamma_R a v·u dS, but transformed to the undeformed domain. (For small deformations, dw_surface_dot would do.)
Could you google a bit to see how people implement that? I did not find much of use.
r.
participants (2)

patricia...＠upf.edu

Robert Cimrman