[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Getfem-users] Adding an initialized fem data in a hyperelastic prob

From: Konstantinos Poulios
Subject: Re: [Getfem-users] Adding an initialized fem data in a hyperelastic problem
Date: Fri, 6 Oct 2017 18:05:25 +0200

Dear David,

Could you also send us a result of your simulation without the active contraction term, only with internal pressure? This is to estimate how large is the contribution of T_a compared to the contribution of the applied pressure p.

Best regards

On Fri, Oct 6, 2017 at 3:07 PM, Yves Renard <address@hidden> wrote:

Dear David,

I do not see any mistake. You can of course simplify the _expression_ in


which should be faster.


Le 04/10/2017 à 18:04, David Danan a écrit :
Dear Getfem users,

i would like to take into account the fiber direction in my model;
to do so i computed the vector associated to each point of the mesh in the order given by mf.point_of_basic_dof, stored it in a vector, used an initialized fem data


and included it in the model (basically, add the component T_a*fibre*fibre^T to the second Piola-Kirchhoff tensor in the reference configuration)
      getfem::add_nonlinear_generic_assembly_brick(model, mim,"T_a*(Id(3)+Grad_u)*[fibre(1)*fibre(1),fibre(1)*fibre(2),fibre(1)*fibre(3);fibre(2)*fibre(1),fibre(2)*fibre(2),fibre(2)*fibre(3);fibre(3)*fibre(1),fibre(3)*fibre(2),fibre(3)*fibre(3)]:Grad_Test_u");
(this is where i have some doubts)

The results i obtained are qualitatively correct but a visible difference still remains between the reference solution and my own solution.

It corresponds to the last problem described in this article
the first two were already validated, so such a difference should (a priori) come from the fiber i guess.

I have tried naively to interpolate the fiber data computed on a P3 mesh on the original P2 mesh but it didn't change anything.

Some help will be really appreciated, i may have overlooked something.

Thanks in advance,

PS: Please, find enclosed the reference solution (deformed.vtk) and the solution obtained (Fibre.vtk)


  Yves Renard (address@hidden)       tel : (33)
  Pole de Mathematiques, INSA-Lyon             fax : (33)
  20, rue Albert Einstein
  69621 Villeurbanne Cedex, FRANCE


reply via email to

[Prev in Thread] Current Thread [Next in Thread]