Dear Colin and David, I have been trying things and looking all over the place but I cannot get to do what I would like to. So please allow me to come back to you in order to make my question clearer. I have a Gmsh 2D mesh where I have defined Physical Surface corresponding to regions having different values of the Lamé's coefficients. Then I would like to set the values of these coefficients with respect to the names or numbers of the regions in my bilinear form. Said that way, it seems obvious but I couldn't find similar tests or examples. Thanks for your help. Nicolas Le mer. 27 mai 2020 à 09:37, Karin&NiKo <niko.karin@gmail.com> a écrit :
Dear Colin and David, Thank you very much. I am trying this.
Le mer. 27 mai 2020 à 09:26, Ham, David A <david.ham@imperial.ac.uk> a écrit :
Dear Nicolas,
You can simply use a DG0 Function for each Lame coefficient, and interpolate a spatially varying expression into it. For a case like this where you want a step function in the coefficient, you probably want to use a ufl.conditional (effectively an “if” expression) in the expression you interpolate.
Regards,
David
*From: *<firedrake-bounces@imperial.ac.uk> on behalf of Karin&NiKo < niko.karin@gmail.com> *Date: *Wednesday, 27 May 2020 at 07:31 *To: *firedrake <firedrake@imperial.ac.uk> *Subject: *[firedrake] Bi-material example
Dear Firedrakers,
I would like to model a bi-material medium. One can think of an elasticity problem where the domain is composed of 2 parts, one being made of steel and the other one of rubber. Then I would like to assign different Lamé coefficients to these parts.
I do not find such an example in Firedrake tests.
Could you please give me some piece of advice?
Best regards,
Nicolas