Hi Patricia,

On 12/05/2014 05:10 PM, patricia...@upf.edu wrote:

Hi all,

I wanted to apply Robin boundary conditions, in the nearly-incompressible 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.