Re: [Getfem-users] Assign different material constants

[Ronan Perrussel]

Le 25/08/2010 15:23, Artur A. Poźniak a écrit :
> Dear Developers and Users,
>
> I created a microstructure for an elastic analysis (2D). So I have a
> box with some regions inside it. My question is how can I assign
> different material constants for these regions? I expect that while
> generating mesh in GMSH the ''physical surface`` would let me
> distinguish different surfaces, but do I need different brick for each
> region?
>
> Thanks in advance!
> Artur Pozniak
>
> _______________________________________________
> Getfem-users mailing list
> address@hidden
> https://mail.gna.org/listinfo/getfem-users
>
>    
Dear Artur,

you have to describe your material by a finite element representation. 
For instance (it is for the matlab interface, quickly extracted from one 
script):
""
mfd = gfMeshFem(m, 1);
mfd.set('fem', gfFem('FEM_PK_DISCONTINUOUS(2,0)'));

% Definition of the medium parameter.
A = ones(mfd.get('nbdof'), 1);
for k = my_physical_region_indices_in_gmsh
ind = mfd.get('dof on region', k);
A(ind) = my_material_paramater(k);
end;
K = gf_asm('laplacian', mi, mf, mfd, A);
""
I hope it helps.

Best regards,
Ronan

ronan_perrussel.vcf
Description: Vcard