node62.html

htmdoc for html coding

 <UL><LI> the sub-domain numbers: option JOPTSD and arrays NUMAT and ITYP:
<P>
  <P><A NAME=5143>&#160;</A><IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img179.gif">
<BR><STRONG>Figure 3.28:</STRONG> <i> Assignment of sub-domain numbers</i><A NAME=5139iAssignmentofsubdomainnumbersi5139>&#160;</A><BR>
<P>
<P>
  <UL><LI> JOPTSD <b>&lt;</b> 0: each element has the value -JOPTSD as sub-domain number.
<P>
   <LI> JOPTSD = 1: the domain is constituted of identical bands in terms of materials. Array
   NUMAT(1:N1,1:N2) associates the value NUMAT(I,J), as sub-domain number, to all the elements 
   whose map on the <b>z=0</b> plane is the element in position I, J.
<P>
   <LI> JOPTSD = 2: the domain (hexaheral elements) is constituted of different bands in terms of materials.
<P>
    Array   NUMAT(1:N1,1:N2,1:max(NTYP)) and array ITYP(1:NB)  associates  the value NUMAT(I,J,NTYP), as sub-domain
    number, to all the elements  whose map on the <b>z=0</b> plane is the element in position I, J,  
    where NTYP = ITYP(K) (i.e. NTYP is the type of band K).
<P>
   <LI> JOPTSD = 3: the domain (pentahedral or tetrahedral elements) is constituted of identical bands 
    in terms of the materials. Array  <BR>[4]
    NUMAT(2(1:N1),1:N2) associates the value NUMAT(I,J), as sub-domain number, to all the elements whose map on the
     <b>z=0</b> plane is one of the 2 triangles in position I,J.
<P>
   <LI> JOPTSD = 4: the domain (pentahedral or tetrahedral elements) is constituted of different bands in terms 
    of materials. Array  <BR>[4]
    NUMAT(2(1:N1),1:N2,1:max(NTYP)) and array ITYP(1:NB) associates the value NUMAT(I,J,NTYP), as sub-domain
    number, to all the elements whose map on the <b>z=0</b> plane is one of the 2 triangles in position I, J, 
    where NTYP = ITYP(K) (i.e. NTYP is the type of band K).
<P>
  </UL>
<P>
  <LI> The reference numbers  (see  figure <A HREF="#figrefgel">3.29</A>):
<P>
  <P><A NAME=5153>&#160;</A><IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img180.gif">
<BR><STRONG>Figure 3.29:</STRONG> <i> Assignment of the reference numbers</i><A NAME=5149iAssignmentofthereferencenumbersi5149>&#160;</A><BR>
<P>
<P>
  <UL><LI> points:  the points of the final mesh  are, either the vertices of the topological cube, or points on its 
   edges, or points on its faces, or lastly, internal points. Consequently, they have the following 
  successive  reference numbers: 19 to 26, 7 to 18, 1 to 6, or 0.
   <LI> edges:  the edges of the final mesh are, either supported by the edges of the topological cube, or
   supported by its faces, or lastly, internal edges. Consequently, they have the following
   successive reference numbers:  7 to 18, 1 to 6, or 0.
   <LI> faces: the faces of the final mesh are, either supported by the faces of the topological cube, or
    are internal faces. Consequently, they have the following successive  reference numbers:
   1 to 6, or 0.
  </UL></UL></UL>
<P>
<P>
<P>
The point numbering is done according to the indications depicted in  figure <A HREF="#fignumgel">3.30</A>.  
The element numbering is done analogously, layer by layer.
<P>
<P><A NAME=5165>&#160;</A><IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img181.gif">
<BR><STRONG>Figure 3.30:</STRONG> <i> Point numbering  by module <b> GEL3D1</b></i><A NAME=5163iPointnumberingbymodule5159bGEL3D1b5159i5163>&#160;</A><BR>
<P>
<P>
<H3><A NAME=SECTION00643200000000000000> Calling the  module</A></H3>
<P>
To call module  <b> GEL3D1</b><A NAME=5106>&#160;</A>, it is necessary to write a main program as follows:
<P>
<P>
<P>
<UL><LI> dimension array M (in the blank common)
<LI> declare the integer array NUMAT(.), array with 2 or 3 indices according to the case
<LI> declare the integer array ITYP(NB) where NB is the number of layers
<LI> declare the real array X(.), array with 1, 2 or 3 indices according to the case
<LI> declare in the same way the real arrays  Y(.) and H(.) 
<LI> declare the externals F1, F2 and F3 (according to the case)
<LI> call <b> INITI</b>   <BR>
<P>
<LI> call the module   <BR>
<P>
 CALL GEL3D1 (M,IOPTM,JOPTSD,NOPTDE,N1,N2,NB,X,Y,H,<BR> 
.  F1,F2,F3,NUMAT,ITYP,NFNOPO,NINOPO,NFAUX)   <BR>
<P>
where
<P>
<UL><LI> M is the super array
<LI> IOPTM is the data description option
<LI> JOPTSD is the sub-domain numbering option
<LI> NOPTDE is the element splitting option
<LI> N1 is the number of points in the <b>x</b> direction
<LI> N2 is the number of points in the <b>y</b> direction
<LI> NB is the  number of layers (i.e. the number of points-1 in the <b>z</b> direction)
<LI> X,Y and H are the arrays with 1, 2 or 3 indices to define the coordinates of the points 
<LI> F1, F2 and F3 are the names of the 3 external functions to supply (according to the value of IOPTM)
  to define the point coordinates. These functions have the following form:
<UL><LI> F1(I,J,XR,YR,K) defines the abscissis of the points from I  and J, the position of the point in the
 plane of definition, 
XR and YR,  the coordinates of  point with indices I and J,  and K , the sectional index (i.e. K=1,NB+1) 
<LI> F2(I,J,XR,YR,K) has the same form and defines the ordinates of the  points
<LI> F3(I,J,XR,YR,K)  has the same form and defines the sides of the points
</UL>
<LI> NUMAT is an array with 2 or 3 indices to assign the sub-domain numbers
<LI> ITYP is the array determining the type of layers (array linked to NUMAT)
<LI> NF(NI)NOPO is the support number of the O.D.S. <b> NOPO</b> and its level
<LI> NFAUX has the value 0
<P>
</UL></UL>
<P>
<b> Remark:</b>   <BR>
<P>
In order to change the references which are automatically assigned, either module <b> MODNOP</b>, or
preprocessor <b> MODNXX</b>, or the keyword <b> NUME</b> of <b> APNOP3</b>,  can be used.
<P>
<HR SIZE=3 WIDTH="75&#37;"><IMG SRC="../icons/smallmod.gif" WIDTH=211 HEIGHT=50 ALIGN=BOTTOM
	ALT="Modulef"><A NAME=tex2html1615 HREF="node61.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/previous_motif.gif"
	ALT="previous"></A><A NAME=tex2html1621 HREF="node59.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/up_motif.gif"
	ALT="up"></A><A NAME=tex2html1623 HREF="node63.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/next_motif.gif"
	ALT="next"></A><A NAME=tex2html1625 HREF="node2.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/contents_motif.gif"
	ALT="contents"></A><A NAME=tex2html1626 HREF="node28.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/index_motif.gif"
	ALT="index"></A><A HREF="../Guide3-18/node62.html"><IMG BORDER=0 SRC="../icons/zoom18.gif" ALIGN=BOTTOM
	ALT="[BIG]"></A><A HREF="../Guide3-14/node62.html"><IMG BORDER=0 SRC="../icons/zoom14.gif" ALIGN=BOTTOM
	ALT="[Normal]"></A><A HREF="../Guide3-10/node62.html"><IMG BORDER=0 SRC="../icons/zoom10.gif" ALIGN=BOTTOM
	ALT="[small]"></A><BR>
<B> Next: </B> <A NAME=tex2html1624 HREF="node63.html">3.4.4 Examples of use</A>
<B>Up: </B> <A NAME=tex2html1622 HREF="node59.html">3.4 Method for a hexahedral topology</A>
<B> Prev: </B> <A NAME=tex2html1616 HREF="node61.html">3.4.2 Description of the input data</A>
<B><A HREF="node28.html"
	>Index</A></B>
<B><A HREF="node2.html"
	>Contents</A></B>
<BR> <HR>
<P><ADDRESS>
</ADDRESS>
</BODY>
</HTML>