MICROMESH_COORDS Computes node coordinates for a micro sampling domain MICROCOORDINATES = MICROMESH_COORDS(NMICRO, CENTER, DELTA) Computes the coordinates of the nodes for the micro sampling domain meshed using rectangles. NMICRO is the number of degrees of freedom per space dimension in the micro domain CENTER contains the coordinates of the center of the sampling domain. These usually are the barycenters for macro triangles or the quadrature nodes for macro parallelograms. DELTA is the width or height of the square 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 [MicroCoordinates]=micromesh_coords(NMicro, center, delta) 0002 %MICROMESH_COORDS Computes node coordinates for a micro sampling domain 0003 % MICROCOORDINATES = MICROMESH_COORDS(NMICRO, CENTER, DELTA) 0004 % Computes the coordinates of the nodes for the micro sampling domain 0005 % meshed using rectangles. 0006 % 0007 % NMICRO is the number of degrees of freedom per space dimension in the 0008 % micro domain 0009 % 0010 % CENTER contains the coordinates of the center of the sampling domain. 0011 % These usually are the barycenters for macro triangles or the quadrature nodes 0012 % for macro parallelograms. 0013 % 0014 % DELTA is the width or height of the square sampling domain. 0015 % 0016 % 0017 % This function should not be modified. 0018 % 0019 % 0020 % The code is available at http://anmc.epfl.ch/ and described in 0021 % further detail in 0022 % 0023 % A. Abdulle and A. Nonnenmacher 0024 % "A short and versatile finite element multiscale code for 0025 % homogenization problems" 0026 % Computer Methods in Applied Mechanics and Engineering, 0027 % http://dx.doi.org/10.1016/j.cma.2009.03.019 0028 % 0029 % Please cite this article in any publication describing research 0030 % performed using the software. 0031 % 0032 % 0033 % Email : assyr.abdulle@epfl.ch and achim.nonnenmacher@epfl.ch 0034 % Last updated : 04/29/2009 with MATLAB 7.4 0035 % 0036 % FE_HMM2D is Copyright (C) 2009 A. Abdulle and A. Nonnenmacher. 0037 % The software is provided free for non-commercial use unter the terms of 0038 % the GNU General Public License. See "copyright.m" for full details. 0039 0040 % 0041 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0042 0043 0044 % divide each dimension into N equidistant nodes 0045 x=linspace(center(1)-.5*delta, center(1)+.5*delta, NMicro); 0046 y=linspace(center(2)-.5*delta, center(2)+.5*delta, NMicro); 0047 0048 % build 2d-grid 0049 [xx,yy] = meshgrid(x,y); 0050 % and reshape so we get one matrix with the coordinates 0051 MicroCoordinates = [xx(:) yy(:)]; 0052 0053 end