test8.di

matlab有限元分析工具,比经较全面的一个手册,请大家下载呀

17:26:59.5 27-Aug-1999
-- >> % QMG test 8: a cube in standard orientation versus

17:26:59.5 27-Aug-1999
-- >> % a tilted cube.

17:26:59.5 27-Aug-1999
-- >> 

17:26:59.5 27-Aug-1999
-- >> global QMG_DATA

17:26:59.5 27-Aug-1999
-- >> cube  = gmoffread(fullfile(QMG_DATA, 'cube.off'));

17:26:59.6 27-Aug-1999
-- >> 

17:26:59.6 27-Aug-1999
-- >> show = 0;

17:26:59.6 27-Aug-1999
-- >> global interactive

17:26:59.6 27-Aug-1999
-- >> if length(interactive) > 0
  show = 1;
end

17:26:59.6 27-Aug-1999
-- >> 

17:26:59.7 27-Aug-1999
-- >> mesh = gmmeshgen(cube, 'show', show);

17:26:60.0 27-Aug-1999
-- >> 

17:26:60.0 27-Aug-1999
-- >> asp = gmchecktri(cube,mesh)
Maximum aspect ratio =        2.70243 achieved in
simplex #12 of topological entity region (3:0) which has vertices 13 4 21 26 
Maximum global side length =  1.9105
Minimum global altitude =     0.597927
Number of nodes = 27 number of elements = 48

asp =

    2.7024


17:27:00.0 27-Aug-1999
-- >> [scrap,numvtx] = size(mesh{4});

17:27:00.0 27-Aug-1999
-- >> global aspprod

17:27:00.0 27-Aug-1999
-- >> global meshsizesum

17:27:00.0 27-Aug-1999
-- >> 

17:27:00.0 27-Aug-1999
-- >> if length(aspprod) > 0
  aspprod = aspprod * asp;
  meshsizesum = meshsizesum + numvtx;
end

17:27:00.1 27-Aug-1999
-- >> 

17:27:00.1 27-Aug-1999
-- >> q = [ -0.84910409704140  -0.40704628877252   0.33665345859509 0
   0.23026379507690  -0.85881459325243  -0.45762001605420 0
   0.47539543231325  -0.31104792751790   0.82295101295431 0];

17:27:00.1 27-Aug-1999
-- >> 

17:27:00.1 27-Aug-1999
-- >> tiltcube = gmapply(q,cube);

17:27:00.1 27-Aug-1999
-- >> 

17:27:00.1 27-Aug-1999
-- >> mesh2 = gmmeshgen(tiltcube, 'show', show);

17:27:01.1 27-Aug-1999
-- >> 

17:27:01.1 27-Aug-1999
-- >> % Display the properties of the second mesh

17:27:01.1 27-Aug-1999
-- >> 

17:27:01.1 27-Aug-1999
-- >> asp2  = gmchecktri(tiltcube, mesh2);
Maximum aspect ratio =        78.5844 achieved in
simplex #606 of topological entity region (3:0) which has vertices 108 6 110 195 
Maximum global side length =  1.55283
Minimum global altitude =     0.00634797
Number of nodes = 228 number of elements = 867

17:27:01.3 27-Aug-1999
-- >> [scrap,numvtx2]  = size(mesh2{4});

17:27:01.3 27-Aug-1999
-- >> 

17:27:01.3 27-Aug-1999
-- >> 

17:27:01.4 27-Aug-1999
-- >> if length(aspprod) > 0
  aspprod = aspprod * asp2;
  meshsizesum = meshsizesum + numvtx2;
end

17:27:01.4 27-Aug-1999
-- >> 

17:27:01.4 27-Aug-1999
-- >> 

17:27:01.4 27-Aug-1999
-- >> % ------------------------------------------------------------------

17:27:01.4 27-Aug-1999
-- >> % Copyright (c) 1999 by Cornell University.  All rights reserved

17:27:01.4 27-Aug-1999
-- >> % See the accompanying file 'Copyright' for authorship information,

17:27:01.4 27-Aug-1999
-- >> % the terms of the license governing this software, and disclaimers

17:27:01.4 27-Aug-1999
-- >> % concerning this software.

17:27:01.4 27-Aug-1999
-- >> % ------------------------------------------------------------------

17:27:01.4 27-Aug-1999
-- >> % This file is part of the QMG software.  

17:27:01.4 27-Aug-1999
-- >> % Version 2.0 of QMG, release date RELDATE.

17:27:01.4 27-Aug-1999
-- >> % ------------------------------------------------------------------

17:27:01.4 27-Aug-1999
-- >>