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