[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Getfem-users] Re : Generic elliptic Brick
From: |
Yves Renard |
Subject: |
Re: [Getfem-users] Re : Generic elliptic Brick |
Date: |
Thu, 10 Feb 2011 13:24:41 +0100 |
User-agent: |
KMail/1.13.5 (Linux/2.6.32-28-server; KDE/4.4.5; x86_64; ; ) |
Dear Jean-François,
You would better use the new bricks because the old ones won't be improved
anymore.
You are right, the correct syntax is
md.add_generic_elliptic_brick(mim, 'varname', 'C')
where you first add the tensor C with
add_initialized_data('A',C)
A priori (I am not completely familiar with arrays in python) C could be a
simple vector with the vlaues in the fortran order or an array with several
dimensions.
yves.
On mercredi 9 février 2011, jeff witz wrote:
> I try to use this method anyway, but I don't succeed setting the 'A'
> parameter
>
> I defined the good tensor for a simple composite case :
>
> C=gf.MdBrick.get(b0, 'param','A')
> C[:]=0
>
> # Definition of the isotropic transverse tensor in is orthotropic plan in
> plan stress
> El=120e6
> Et=8e6
> G=4
> Nult=0.3
> Nutl=0.01875
> C[0,0,0,0]=El/(1-Nult*Nutl)
> C[0,0,1,1]=Nult*Et/(1-Nult*Nutl)
> C[1,1,0,0]=C[0,0,1,1]
> C[0,1,0,1]=G
>
>
> But once there I don't find any method allowing me to set this tensor to
> 'A' og the brick.
>
> I also try to use the new way for bricks as to know :
> md.add_generic_elliptic_brick(....)
>
> but I don't find the good syntaxe.
> Could you give me the good syntaxe ?
>
>
> I think if I can use this last method I could set the model with :
> add_initialized_data('A',C)
--
Yves Renard (address@hidden) tel : (33) 04.72.43.87.08
Pole de Mathematiques, INSA-Lyon fax : (33) 04.72.43.85.29
20, rue Albert Einstein
69621 Villeurbanne Cedex, FRANCE
http://math.univ-lyon1.fr/~renard
---------