test1.di

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

17:19:56.2 27-Aug-1999
-- >> %% QMG test 1

17:19:56.3 27-Aug-1999
-- >> %% This is an approximation to the mesh on the cover of Claes 

17:19:56.3 27-Aug-1999
-- >> %% Johnson's book ``Numerical solutions of partial differential

17:19:56.3 27-Aug-1999
-- >> %% equations by the finite element method'' , Cambridge U. Press 1987.

17:19:56.3 27-Aug-1999
-- >> %% It has greater refinement near the origin.

17:19:56.3 27-Aug-1999
-- >> global interactive

17:19:56.3 27-Aug-1999
-- >> wedgeangle = pi / 6.0;

17:19:56.3 27-Aug-1999
-- >> %% Make a circular approximation with 4 Bezier curves from pi/6 to 2pi.

17:19:56.3 27-Aug-1999
-- >> [verts,codes] = gm_circ_approx(wedgeangle, 2*pi, 4);

17:19:56.4 27-Aug-1999
-- >> newverts = [0,0;verts];

17:19:56.4 27-Aug-1999
-- >> newcodes = [0;codes];

17:19:56.4 27-Aug-1999
-- >> domain = gm_cpoly(newverts, newcodes);

17:19:56.5 27-Aug-1999
-- >> % The size function is of the form h0*r^p+f, where p is determined

17:19:56.5 27-Aug-1999
-- >> % by (4.31) in Johnson, and f is determined so that the H1 error on

17:19:56.5 27-Aug-1999
-- >> % the innermost triangles matches the H1 error everywhere else.

17:19:56.5 27-Aug-1999
-- >> h0 =  0.2;

17:19:56.5 27-Aug-1999
-- >> gamma = pi / (2*pi - wedgeangle);

17:19:56.5 27-Aug-1999
-- >> f  = h0 ^ (1.0/gamma);

17:19:56.5 27-Aug-1999
-- >> p = 1.0 - gamma;

17:19:56.5 27-Aug-1999
-- >> sizectl = sprintf('(formula (%d * (%%0^2+%%1^2)^(%d/2.0) + %d))', h0, p, f);

17:19:56.5 27-Aug-1999
-- >> show =  0;

17:19:56.5 27-Aug-1999
-- >> if length(interactive) > 0
  disp(sprintf('size control string = %s', sizectl));
  show = 1;
end

17:19:56.5 27-Aug-1999
-- >> mesh = gmmeshgen(domain, 'size', sizectl, 'show', show);

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

17:20:01.4 27-Aug-1999
-- >> asp  = gmchecktri(domain,mesh);
Maximum aspect ratio =        10.2397 achieved in
simplex #359 of topological entity mregion (2:0) which has vertices 13 11 230 
Maximum global side length =  0.200806
Minimum global altitude =     0.00727931
Number of nodes = 259 number of elements = 447

17:20:02.1 27-Aug-1999
-- >> global aspprod

17:20:02.1 27-Aug-1999
-- >> global meshsizesum

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

17:20:02.1 27-Aug-1999
-- >> if length(interactive) > 0
  gmviz(mesh)
end

17:20:02.1 27-Aug-1999
-- >> % ------------------------------------------------------------------

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

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

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

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

17:20:02.2 27-Aug-1999
-- >> % ------------------------------------------------------------------

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

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

17:20:02.2 27-Aug-1999
-- >> % ------------------------------------------------------------------

17:20:02.2 27-Aug-1999
-- >>