[Getfem-users] Generic elliptic Brick

[jeff witz]

Bonjour,

I'm trying to use the generic elliptic brick under the python binding, and I get some troubles to undertand how to set the tensor for a classical anisotropic elasticity.

When I set  : 

b0 = gf.MdBrick('generic elliptic',mim,mfu,'matrix')
# print the parameter list of the generic elliptic brick
 gf.MdBrick.get(b0, 'param list') # give ('A',) 
# print the parameter A
print(gf.MdBrick.get(b0, 'param','A')) # by default the A tensor is set to identity


I get :
[[ 1.  1.]
 [ 1.  1.]]

This can't define a 4rth order tensor.

I set this way : 

b0 = gf.MdBrick('generic elliptic',mim,mfu,'tensor')
# print the parameter list of the generic elliptic brick
 gf.MdBrick.get(b0, 'param list') # give ('A',) 
# print the parameter A
print(gf.MdBrick.get(b0, 'param','A')) # by default the A tensor is set to identity

I get an object that seems to have the good shape : (2,2,2,2)

Is there a way to handle 4rth order tensor in the normailized Voigt basis. I want use this code in a course so i would like ti get the same formalism.

Best regards.


-- 
Jean-François WITZ