[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Getfem-users] Assign different material constants
From: |
Ronan Perrussel |
Subject: |
Re: [Getfem-users] Assign different material constants |
Date: |
Wed, 25 Aug 2010 16:47:26 -0300 |
User-agent: |
Mozilla/5.0 (X11; U; Linux i686; en-US; rv:1.9.1.11) Gecko/20100713 Lightning/1.0b1 Thunderbird/3.0.6 |
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