MICROMESH_ELEMENTS Computes the micro sampling domain triangulation
MICROELEMENTS = MICROMESH_ELEMENTS(NMicro)
Computes the triangulation of the micro sampling domain using
square micro elements and returns the node numbers of the
triangulation.
N=NMICRO is the number of degrees of freedom per space dimension in the
micro domain
MICROELEMENTS has dimension (NMicro-1)^2 x 4 and contain the node
numbers of the triangulation of the micro sampling domain.
This function should not be modified.
The code is available at http://anmc.epfl.ch/ and described in
further detail in
A. Abdulle and A. Nonnenmacher
"A short and versatile finite element multiscale code for
homogenization problems"
Computer Methods in Applied Mechanics and Engineering,
http://dx.doi.org/10.1016/j.cma.2009.03.019
Please cite this article in any publication describing research
performed using the software.
Email : assyr.abdulle@epfl.ch and achim.nonnenmacher@epfl.ch
Last updated : 04/29/2009 with MATLAB 7.4
FE_HMM2D is Copyright (C) 2009 A. Abdulle and A. Nonnenmacher.
The software is provided free for non-commercial use unter the terms of
the GNU General Public License. See "copyright.m" for full details.0001 function [MicroElements]=micromesh_elements(N) 0002 %MICROMESH_ELEMENTS Computes the micro sampling domain triangulation 0003 % MICROELEMENTS = MICROMESH_ELEMENTS(NMicro) 0004 % Computes the triangulation of the micro sampling domain using 0005 % square micro elements and returns the node numbers of the 0006 % triangulation. 0007 % 0008 % N=NMICRO is the number of degrees of freedom per space dimension in the 0009 % micro domain 0010 % 0011 % MICROELEMENTS has dimension (NMicro-1)^2 x 4 and contain the node 0012 % numbers of the triangulation of the micro sampling domain. 0013 % 0014 % 0015 % This function should not be modified. 0016 % 0017 % 0018 % The code is available at http://anmc.epfl.ch/ and described in 0019 % further detail in 0020 % 0021 % A. Abdulle and A. Nonnenmacher 0022 % "A short and versatile finite element multiscale code for 0023 % homogenization problems" 0024 % Computer Methods in Applied Mechanics and Engineering, 0025 % http://dx.doi.org/10.1016/j.cma.2009.03.019 0026 % 0027 % Please cite this article in any publication describing research 0028 % performed using the software. 0029 % 0030 % 0031 % Email : assyr.abdulle@epfl.ch and achim.nonnenmacher@epfl.ch 0032 % Last updated : 04/29/2009 with MATLAB 7.4 0033 % 0034 % FE_HMM2D is Copyright (C) 2009 A. Abdulle and A. Nonnenmacher. 0035 % The software is provided free for non-commercial use unter the terms of 0036 % the GNU General Public License. See "copyright.m" for full details. 0037 0038 % 0039 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0040 0041 0042 MicroElements=zeros( (N-1)*(N-1), 4); 0043 for j=0:N-2 0044 MicroElements(1+j*(N-1):N-1+j*(N-1),:)=... 0045 [(1+j*N:N-1+j*N)', (N+1+j*N:N+N-1+j*N)',... 0046 (N+2+j*N:N+N+j*N)', (2+j*N:N+j*N)']; 0047 end 0048 0049 end