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 [email protected] 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
<<attachment: ronan_perrussel.vcf>>
_______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
