← Back to team overview

dolfin team mailing list archive

Re: [Question #142110]: defining a varying 2nd order tensor

 

Question #142110 on DOLFIN changed:
https://answers.launchpad.net/dolfin/+question/142110

Sungick Kim posted a new comment:
https://answers.launchpad.net/dolfin/+question/102501

I can write the tensor that way. FEniCS doesn't complain about it.
But I think that it's not solving problem.
Here are code and the result

<my code>
f00 = Expression("(10.45-9.96)/10.45 *u1c")
f11 = Expression("(5.58-6.05)/5.58 *u1c")
f22 = Expression("(4.86-4.74)/4.86 *u1c")
z = Constant(0.0)
F2= as_matrix( (f00,z,z),(z,f11,z),(z,z,f22) )

F1 = I + grad(u2c)          
F = F1*F2                 # Deformation gradient


<result>

(1) computing Transport problem
  Solving nonlinear variational problem
    Starting Newton solve.
      Applying boundary conditions to linear system.
      Applying boundary conditions to linear system.
      Solving linear system of size 7381 x 7381 (PETSc LU solver, umfpack).
      Applying boundary conditions to linear system.
      Newton iteration 1: r (abs) = 0.000e+00 (tol = 1.000e-10) r (rel) = -nan (tol = 1.000e-14)

-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.