node15.html

htmdoc for html coding

of a generalized quadrangle. The input data consists of the contour of the domain in terms of lines (keyword
 <b> LIGN</b>).
<P>
<UL><LI> QUAC<A NAME=617>&#160;</A>  (A) : the keyword to activate
  <LI> IMPRE (I) : print parameter for execution  
  <LI> LEVEL (I) : the level of the resulting D.S. 
  <LI> NUDSD (I) : the sub-domain number to assign to the elements created
  <LI> NBRELI (I) : the number of lines constituting the contour of the domain
  <LI> NS1L (I) : the direction of the first line of the contour: if NS1L = 1 , the line
        runs from  NEXTR1 towards NEXTR2, if  NS1L = -1, the line runs in the opposite  direction 
        (i.e. from NEXTR2 to NEXTR1)
  <LI> Loop from  J=1 to  NBRELI, the number of lines constituting the contour
  <UL><LI>  NOLIG(J) (I) : the number of line J
  </UL>
  <LI> end of loop J 
  <LI> IMAX (I) : the number of points on the first side of the quadrangle to be meshed (endpoints included)
  <LI> NQUAD (I): element splitting option:
  <UL><LI> 1 : mesh consisting of quadrangles,
    <LI> 0 : mesh consisting  of  triangles (choice of the best angle to cut the quadrangles),
    <LI> -1 : same, with special treatment of the topological corners of the initial  quadrangle,
    <LI> -2 (-3) : regular mesh (choice of the same diagonal  for each splitting of quadrangles) with
         special treatment  at the corners,
    <LI> -4 (-5) : regular mesh without  special treatment at the corners.
  </UL></UL>
<P>
<P>
<P>
<b> Remarks:</b>
<P>
<UL><LI> Two consecutive lines must have a common endpoint.
<LI> The number of points on the contour must be even. In fact, the domain is considered as a quadrilateral
and the number of points on two opposite sides must be the same. Consequently, the knowledge of IMAX is
sufficient to determine the splitting completely.
<P>
<LI> IMAX <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img6.gif"> 2, same for JMAX (the equivalent of  IMAX for the 2 other sides, calculated implicitly).
</UL>
<P>
Figures <A HREF="#figquac1">1.2</A> to <A HREF="#figquac21">1.5</A> shows the effect of parameter NQUAD.
<P>
 fig5bd.psfig:quac1Regular splittingfig5be.psfig:quac11Other choice
<P>
 fig5bf.psfig:quac2Other splittingfig5bh.psfig:quac21Other choice
<P>
<b> TRIC </b> : Calls the mesh generator <b> TRICOO</b> which corresponds to an algebraic method<A NAME=642>&#160;</A>
for the case of a generalized triangle. The input data consists of the contour of the domain in terms of the lines
 (keyword <b> LIGN</b>).
<P>
<UL><LI> TRIC<A NAME=645>&#160;</A>  (A) : the keyword to activate 
  <LI> IMPRE (I) : print parameter for execution  
  <LI> LEVEL (I) : the level of the resulting D.S.
  <LI> NUDSD (I) : the sub-domain number to assign to the elements created
  <LI> NBRELI (I) : the number of lines constituting the contour of the domain
  <LI> NS1L (I) : the direction of the first contour line: if NS1L = 1 , the line runs from endpoint
    NEXTR1 towards endpoint NEXTR2, if NS1L = -1, the line
   runs in the opposite direction (i.e. from  NEXTR2 to NEXTR1)
  <LI> Loop from J=1 to NBRELI, the number of contour lines
  <UL><LI>  NOLIG(J) (I) : the number of line  J
  </UL>
  <LI> end of loop J 
  <LI> IMAX (I) : the number of points on the first side of the triangle to mesh (endpoints included)
</UL>
<P>
<P>
<P>
<b> Remarks:</b>
<P>
<UL><LI> Two consecutive lines must have a common endpoint.
<P>
<LI> The number of points on the contour must be a multiple  of 3. In fact, the domain corresponds to a 
triangle and the number of points on its topological sides must be equal. Consequently, the knowledge
of  IMAX  suffices to determine the splitting completely.
<LI> IMAX <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img6.gif"> 2.
</UL>
<P>
<b> TRIA  </b> : Calls the mesh generator <b> TRIFRO</b> which corresponds to a 
advancing-front  method<A NAME=654>&#160;</A>.
The input data consists of the contour of the domain in terms of lines 
(keyword <b> LIGN</b>).  Internal points
can also be specified. <BR> Furthermore, this mesh generator accepts contours
 with several components ; it allows to generate triangles or quadrangles 
meshes and/or adaptive meshes.
<P>
<UL><LI> TRIA<A NAME=657>&#160;</A>  (A) : the keyword to activate 
  <LI> IMPRE, NIVEAU, NUDSD, NBRELI, NPROPA as for keyword TRIH
  <LI> Loop from J=1 to NBRELI the number of contour lines
  <UL><LI>  NOLIG(J) (I) : the number of the line J of the contour.
  </UL>
  <LI> end of loop  J 
  <LI> NCOMP, NBRINT like the keyword  TRIH
  <LI> IOPT (I) : front management option (i.e.  selection of the
        departure segment)
  <UL><LI> 0 : mesh consisting of  triangles - slow option 
    <LI> 1 : mesh consisting of   triangles  - fast option  (this option is advised)
    <LI> 10 :  mesh consisting of  triangles - slow option - adaptive version
    <LI> 11 :  mesh consisting of  triangles - fast option - adaptive version
    <LI> 20 : mesh consisting of  quadrangles - slow option 
    <LI> 21 :  mesh consisting of  quadrangles - fast option
    <LI> 30 : mesh consisting of  quadrangles - slow option - adaptive version
    <LI> 31 : mesh consisting of  quadrangles - fast option - adaptive version
  </UL>
  <LI> if IOPT = 10, 11, 30 or 31
        <UL><LI> NOMFIC (C) : name of the file that contains the initial mesh
        <LI> NIVEAU (I) : level of this mesh
        <LI> NOMFIC (C) : name of the file that contains the
data structure B, this data structure B contains for each point of the initial
 mesh  the size of the cells in its neighborhood
        <LI> NIVEAU (I) : level of this data structure
        </UL>
<P>
  <LI> Loop from  K=1 to  NCOMP the number of contour components
  <UL><LI>  NTACO(K) (I) : the number of lines of component  K of the contour  (see TRIH).
  </UL>
  <LI> end of loop K 
  <LI> Loop from L=1 to NBRINT the number of fixed internal points 
 (if <b>&gt;</b> 0)
  <UL><LI>  NUMPIN(L) (I) : the internal point number L (see TRIH) ;
    <LI>  DIAM(L) (R) : the desired diameter for the triangles created around point L.
  </UL>
  <LI> end of loop L 
</UL>
<P>
<P>
<P>
<b> Remarks:</b>
 <P>
<P>
This keyword requires the same data than the following  keyword  (TRIH), except
 the IOPT parameter describe above. 
<P>
<P>
<b> TRIH </b> : Calls the mesh generator <b> TRIHER</b> which corresponds
to a Voronoi method<A NAME=672>&#160;</A>.
The input data consists of the contour of the domain in terms of lines (keyword <b> LIGN</b>). Internal points
can also be specified. <BR> 
Furthermore, this mesh generator accepts contours with several components. It allows to generate adaptive meshes. In its standard version (non adaptive) internal edges can also be specified.
<P>
<UL><LI> TRIH<A NAME=675>&#160;</A>  (A) : the keyword to activate 
  <LI> IMPRE (I) : print parameter for execution  
  <LI> LEVEL (I) : the level of the resulting D.S.
  <LI> NUDSD (I) : the sub-domain number to assign to the elements created
  <LI> NBRELI (I) : the number of contour lines of the domain
  <LI> NPROPA (I) : propagation type and coefficient  (e.g. 1, see remark)
  <LI> Loop from  J=1 to NBRELI, the number of contour lines
  <UL><LI>  NOLIG(J) (I) : the number of line  J of the contour described in the positive direction for
   the external part and in the inverse direction for the eventual internal parts.
  </UL>
  <LI> end of loop J 
  <LI> NCOMP (I) : the number of components connected to the boundary<FONT color="#B8860B" SIZE=-1><I>(<IMG ALIGN=middle SRC="../icons/foot_motif.gif"
	ALT="[note]"> i.e. the number of internal components
       +1, in the preprocessor we also ask the number of internal components (holes and internal edges).)</I></FONT>
  <LI> NBRINT (I) : the number of fixed internal points
  <LI> IOPT (I) : internal point creation option:
  <UL><LI> 0 : the mesh generator does not generate any internal points
    <LI> 1 : the  mesh generator generates all the internal points
    <LI> -1 : the mesh generator generates internal points in addition to those points already given (NBRINT)
    <LI> 11 :  mesh consisting of  triangles - adaptive version
  </UL>
<P>
<P>
  <LI> if IOPT =  11
        <UL><LI> NOMFIC (C) : name of the file that contains the initial mesh
        <LI> NIVEAU (I) : level of this mesh
        <LI> NOMFIC (C) :  name of the file that contains the
data structure B, this data structure B contains for each point of the initial
 mesh  the size of the cells in its neighborhood
        <LI> NIVEAU (I) : level of this data structure
        </UL>
  <LI> Loop from K=1 to NCOMP, the number of contour components
  <UL><LI>  NTACO(K) (I) : the number of lines of component K of the contour. The first line of this
 component runs in the direction NEXTR1 to NEXTR2 if NTACO(K) <b>&gt;</b> 0, and in the opposite direction if not.
  </UL>
  <LI> end of loop K 
  <LI> Loop from  L=1 to NBRINT, the number of fixed internal points  (if <b>&gt;</b> 0)
  <UL><LI>  NUMPIN(L) (I) : the internal point number L (taken from the list of points defined with keyword
     <b> POIN</b>);
    <LI>  DIAM(L) (R) : the desired diameter for the triangles created around point L.
  </UL>
  <LI> end of loop L 
</UL>
<P>
<P>
<P>
<b> Remarks:</b>
<P>
<UL><LI> Two consecutive lines (on the same contour component) must have a common endpoint.
<LI> A line cannot appear twice in the contour.
<LI> The fixed internal  points allow, on the one hand, to force the mesh to approach these points and, on
the other hand, to specify the density of the mesh in their neighborhood (according to the value of DIAM).
<LI> The parameter NPROPA is used when the internal points are created. They are distributed 
 in a geometric manner if  NPROPA is positive and in a arithmetic manner if it is negative.
<P>
<LI> The position of the internal points is a function of the point positions on the given boundary. 
A fine mesh at the level of the boundary introduces a fine mesh of the interior. To dilute the influence 
of the data on the interior of the domain, the value NPROPA can be changed (in fact, its absolute value <b>a</b>). 
From <b>a</b>, the propagation coefficient is calculated: <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img7.gif">.
<P>
Therefore,  <b>a=1</b> has a neutral effect; for <b>a=2,3,4,5,</b> etc., we obtain successively <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img8.gif"> etc.
<P>
<LI> In its standard version (non adaptative) internal edges can be specified,
these edges can intersect the components.
<P>
<LI> During the use of the  adaptative version the user has to give the following supplementary data
  <UL><LI> the description of the domain : points and lines that describe the components  (reference of the points must be negative)
  <LI> an initial mesh of the domain
  <LI> a data structure B  that contains for each point of the initial
 mesh  the size of the cells in its neighborhood
  </UL></UL>
<P>
<b> MANU </b> : Calls the mesh generator <b> CONOPO</b> which corresponds to a manual description<A NAME=705>&#160;</A>,
 i.e., from points (keyword <b> POIN</b>).
<P>
<UL><LI> MANU<A NAME=708>&#160;</A>  (A) : the keyword to activate 
  <LI> IMPRE (I) : print parameter for execution  
  <LI> LEVEL (I) : the level of the resulting D.S.
  <LI> NE (I) : the number of elements to create
  <LI> Loop from J=1 to NE, the number of elements
  <UL><LI> NCGE (I) : type of element: 1 <IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img5.gif"> point, 2 <IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img5.gif"> segment,
     3 <IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img5.gif"> triangle, <BR> 4 <IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img5.gif"> quadrangle
    <LI> NUDSD (I) : the sub-domain number of the element
    <LI> Loop from  K=1 to the number of element vertices
    <DL COMPACT><DT>.
<DD> NUMSOM(K) (I) : the number of element vertices
<P>
 </DL>
    <LI> end of loop K <BR> 
    <LI> Loop from L=1 to the number of element edges   (if NTYP <b>&gt;</b> 2)
    <DL COMPACT><DT>.