1

Within Fenics what ordering is used for storing vectors as upper triangular matrices? Th:e question really is I define a green Lagrange strain tensor in the following way:

V = VectorFunctionSpace(mesh, "Lagrange", 1)
u  = Function(V)
I = Identity(V.cell().d)    # Identity tensor
F = I + grad(u)             # Deformation gradient
C = F.T*F                   # Right Cauchy-Green tensor
E = (C-I)/2                 #green strain tensor

But when I accesses the elements later,

psi = stress1*E[0,0]+stress2*E[1,1]+stress3*(E[2,2])+stress4*E[1,2]+stress5*E[0,2]+stress6*E[0,1]

I assumed Voigt notation but it occurs to me that this may not be the case. Using this method to initialize the strain tensor is this the notation I should be using or is there a different type or initialization that would be more appropriate?

Thanks for any help.

Nick
  • 158
  • 1
  • 7
  • Hopefully you have also tried http://fenicsproject.org/support/. The subject is not very familiar to me. –  Feb 03 '14 at 06:06
  • No I hadn't. Thanks for the tip, its been a while since I have used fenics and it appears they have started their own Q & A section since I was last there. – Nick Feb 03 '14 at 06:31
  • "Pure" FEniCS question (i.e., about implementation details) have indeed moved to http://fenicsproject.org/qa. (Since it's a different platform, your question cannot be migrated unfortunately, so you'll have to reask there.) Questions with a mathematical or modeling background are still welcome here, though! (Not your question is not welcome; it just makes sense to collect all support-type questions in a single place.) – Christian Clason Feb 03 '14 at 07:52

1 Answers1

0

Adding the answer here in case others stumble on it. Answer is from Anders Logg, a fenics dev

Your definition of the tensors F, C, E is correct. FEniCS does not know anything about Voigt notation, so you need to access the components directly as E[0, 0] etc.

Nick
  • 158
  • 1
  • 7