node44.html

来自「htmdoc for html coding」· HTML 代码 · 共 747 行 · 第 1/2 页

<P>
<DT>Array MAI5:
<DD> Number of unknowns per type of nodes. <BR>
<P>
This integer array contains  NTYNOE words:
<P>
<UL><LI> MAI5(i), the number of unknowns of node type i
</UL>
<P>
<DT>Array MAI6:
<DD> Unknown numbers. <BR>
<P>
This integer array contains <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img127.gif"> MAI5(i)  words:
<P>
<UL><LI> MAI6(j), the number of unknowns j
</UL>
<P>
<DT>Array MAI7:
<DD> Number of degrees of freedom per unknown. <BR>
<P>
This integer array contains <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img127.gif"> MAI5(i)  words:
<P>
<UL><LI> MAI7(j), the number of degrees of freedom of unknown j
</UL>
<P>
<DT>Array MAI8:
<DD> The mnemonic code corresponding to the degree of freedom. <BR>
<P>
This integer array contains <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img128.gif"> MAI7(i)   words:
<P>
<UL><LI> MAI8(j), the mnemonic code corresponding to degree of freedom j
</UL>
<P>
The codes are the following;
<UL><LI> 1: VN: nodal value
<LI> 2: DX: X derivative at a node
<LI> 3: DY: Y derivative at a node
<LI> 4: DZ: Z derivative at a node
<LI> 5: DTG: derivative tangential to an edge or face at a node
<LI> 6: DDXX: second derivative with respect to X at a node
<LI> 7: DDYY: second derivative with respect to Y at a node
<LI> ... (see subroutine <b> MNEMO</b> below)
</UL>
<P>
<DT>Array MAI9:
<DD> Number of mnemonics per point type. <BR>
<P>
This integer array contains  NTYPOI words:
<P>
<UL><LI> MAI9(i), the number of mnemonics of point type i
</UL>
<P>
<DT>Array MAIA:
<DD> The mnemonic number per point type. <BR>
<P>
This integer array contains  <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img129.gif"> MAI9(i)  words:
<P>
<UL><LI> MAIA(i), the mnemonic code of point type i (see <b> MAI8</b> or <b> MNEMO</b>),
</UL>
<P>
<DT>Array MAIB:
<DD> Libraries and finite elements.  <BR>
<P>
This integer and character array contains  4 NTYELM  words:
<P>
<UL><LI> MAIB(1,i), the library name of a  type i finite element
<LI> MAIB(2,i), the first word of the name of a type i finite element 
<LI> MAIB(3,i), the second words of the name of a  type i  finite element 
<LI> MAIB(4,i), the code number of a type i finite element.
</UL>
<P>
<DT>Array MAIC:
<DD> Geometry of the finite element types. <BR>
<P>
This integer array contains  NTYELM words:
<P>
<UL><LI> MAIC(i), the geometric code of a  type i finite element
</UL>
<P>
This code is described in DS <b> NOPO</b> (3: triangle, 4: quadrangle, etc.).
<P>
<DT>Array MAID:
<DD> Number of nodes per finite element type. <BR>
<P>
This integer array contains  NTYELM words:
<P>
<UL><LI> MAID(i), the number of nodes of a finite element of type i
</UL>
<P>
<DT>Array MAIE:
<DD> Number of points per finite element type. <BR>
<P>
This integer array contains  NTYELM words:
<P>
<UL><LI> MAID(i), the number of points of type i  finite element
</UL>
<P>
<DT>Array MAIF:
<DD> The number of associated arrays per finite element type. <BR>
<P>
This integer array contains  NTYELM words:
<P>
<UL><LI> MAIF(i), the number of arrays for finite element type i
</UL>
<P>
<DT>Array MAIG:
<DD> Pointer in <b> MAIJ</b>.<BR>
<P>
This integer array contains  NTYELM  words:
<P>
<UL><LI> MAIG(i), address in <b> MAIJ</b> of the words preceding the nodal type numbers of a type i
finite element
</UL>
<P>
<DT>Array MAIH:
<DD> Pointer in  <b> MAIJ</b>. <BR>
<P>
this integer array contains  NTYELM words:
<P>
<UL><LI> MAIG(i), address in  <b> MAIJ</b> of the word preceding the point type numbers of a type i finite element
</UL>
<P>
<DT>Array MAII:
<DD> Pointer in  <b> MAIJ</b>. <BR>
<P>
this integer array contains NTYELM words:
<P>
<UL><LI> MAII(i), address in <b> MAIJ</b> of the first word describing the arrays associated with a type i 
finite element
</UL>
<P>
<DT>remark:
<DD> These addresses are relative to the beginning of array <b> MAIJ</b>.
<P>
 <DT>Array  MAIJ:
<DD> Pointer in  <b> MAIK</b> and description of the associated arrays. <BR>
<P>
This integer array contains:
<P>

<UL><LI> Loop i from 1 to NTYELM
  <UL><LI> address in array <b> MAIK</b> of the first coordinate of the first node  of the i-th
     type element
     <LI>
     Loop j from 1 to  MAID(i)
     <UL><LI> node j type 
     </UL> 
     End of loop j
<P>
     <LI> address in array <b> MAIK</b> of the first coordinate of the first point of the i-th
     type element
     <LI> loop j from 1 to MAIE(i)
     <UL><LI> point j type 
     </UL> 
     End of loop j
     <LI> if MAIF(i) = 0, go to the end of element type loop
     <LI> Loop j from 1 to MAIF(i)
     <UL><LI> array type
        <LI> number of indices
        <LI> maximum value of each index
     </UL>
     End of loop j
  </UL> 
  End of loop i
</UL>
<P>
LMAIJ, the length of array <b> MAIJ</b>, equals the number of words required to
store all this information. In the file, this value can be recuperated as each record starts with
the number of words which  follow.
<P>
<DT>Array MAIK:
<DD> Point and node coordinates of the reference element. <BR>
<P>
This array, of type given in <b> MAIJ</b>,  contains:
<P>

<UL><LI> Loop i from 1 to  NTYELM
  <UL><LI> Loop j from 1 MAID(i)
     <UL><LI> the NDIM coordinates of node j in the reference element
     </UL> 
     End of loop j
     <LI> Loop j from 1 to MAIE(i)
     <UL><LI> the NDIM coordinates of point j in the reference element
     </UL> 
     End of loop j
  </UL> 
  End of loop i
</UL>
<P>
Nodes and points are arranged in the order defined by DS  <b> NOPO</b> (see this DS). If the reference  element  
is typically two-dimensional, but is considered as three-dimensional, it suffices to
set the z components to 0.
<P>
<DT>Array MAIL:
<DD> Element description. <BR>
<P>
This double precision complex array contains:
<P>

<UL><LI> Loop i from 1 to NE
  <UL><LI> NEL, the element number
     <LI> NTYE, the element type (between 1 and NTYELM)
     <LI> NMAE, 0 if the element is internal, or the number of words necessary to store the reference numbers (see below)
     <LI> NDSDE, the sub-domain number  
     <LI> NTALIR, the number of arrays associated with the element (see below)
     <LI> if NCACAR (of <b> MAI2</b> is non-zero: NCCAR, 0 if all the element  points do not have characteristics
           other than the coordinates, or the address in DS <b> COOR</b> of the first characteristic other than the coordinates  
     <LI> NNO, the number of nodes in the element
     <LI> Loop j from 1 to NNO
     <UL><LI> NONO(j), the number of node j in the element
     </UL> 
     End of loop j
     <LI> if NCOPNP (see array <b> MAI2</b>) is zero:  
     <UL><LI> NPO, the number of points in the element
        <LI> Loop j from 1 to NPO  
        <blockquote> NOPO(j), the number of point j in the element
        </blockquote> 
        End of loop j
     </UL>  
     <LI> if NTALIR <b>&gt;</b> 1: loop over the associated arrays to read and for each one its coefficients.
     <LI> if NMAE is non-zero, NTALIR is at least 1 and we have:
        <UL><LI> INING, indicator of the level of information on the reference numbers (see array <b> NOP5</b> of
                      DS <b> NOPO</b>).
<P>
<P>
           Loop from 2 to NMAE
           <OL><LI> if INING = 3, vertex references,
              <LI> if INING = 2, edge and vertex references,
              <LI> if INING = 1, face, edge and vertex references.
           </OL>
<P>
           End of loop
        </UL></UL> 
  End of loop i
</UL>
<P>
 </DL>
<P>
<P>
<P>
<H2><A NAME=SECTION05192000000000000000>1.8.2 Corresponding tools</A></H2>
<P>
<P>
<P>
<H3><A NAME=SECTION05192100000000000000> Reading and writing a DS MAIL</A></H3>
<P>
<P>
<P>
If NOPTFI=0, the DS is category 1; consequently, in this case, a DS residing in main memory is
written wholly on sequential file via module <b> SDECRI</b>. 
Similarly, it is read from file and transferred wholly into main memory via
 module <b> SDLECT</b>.
<P>
If NOPTFI=1, the DS is category 2; array <b> MAIL</b> is read and written outside of these modules. Each 
element constitutes a record:
<P>
<blockquote> WRITE(NFMAIL) LE,(M(IAMAIL-1+I),I=1,LE)  
</blockquote>
<P>
<P>
<P>
<H3><A NAME=SECTION05192200000000000000> Printing a DS MAIL and plotting a mesh</A></H3>
<P>
<P>
<P>
The contents (total or partial)  of a DS <b> MAIL</b> can be printed by the following modules:
<P>
<DL COMPACT><DT>IMMAIL:
<DD> the whole DS (according to the print parameter IMPRE)
<P>
<PRE>      SUBROUTINE  IMMAIL(M,NFMAIL,NIMAIL,IMPRE)
C AIM : PRINT THE ARRAYS OF DS MAIL, INCREASING WITH IMPRE
C ---  FILE NFMAIL (OR MC IF 0), LEVEL NIMAIL
</PRE>
<P>
<DT>INFOMA:
<DD> print a particular value contained in the DS
<P>
<PRE>      SUBROUTINE INFOMA(M,IOPT,TYPESD,NIVEAU,TAB2,TAB3,TABC,TABI,TABJ,
     +                  IAMAIL,NUMREC)
C AIM : GIVE INFORMATION ON A DS MAIL
</PRE>
<P>
Preprocessors <b> IMAGXX</b> and <b> INFOXX</b> are used if a conversational call of <b> IMMAIL</b> and
 <b> INFOMA</b>, respectively, is desired.
<P>
 </DL>
<P>
The mesh contained in a file can be plotted via preprocessor <b> TRMCXX</b> (see [<A HREF="node57.html#mod_96">96</A>])
in two dimensions, or <b> TRC3XX</b> in three dimensions which also uses DS <b> COOR</b>.
<P>
<P>
<P>
<H3><A NAME=SECTION05192300000000000000> Creation and modification modules of DS MAIL</A></H3>
<P>
<P>
<P>
DS <b> MAIL</b> is constructed via modules <b> COMACO</b>, <b> COMAC2</b> or <b> COMAC3</b>.
<P>
<P>
<P>
<H3><A NAME=SECTION05192400000000000000> Manipulation utilities of DS MAIL</A></H3>
<P>
<P>
<P>
 Numerous utilities exist for manipulating the contents of a DS <b> MAIL</b>, amongst which
the utilities facilitating the exploration of array <b> MAIL</b>.
<P>
<PRE> 

      SUBROUTINE CLMAI5(MAI2,MAIC,MAIL,COO4,IA,NTYE,NDSDE,NNO,NONO,NPO
     +                  ,NOPOI,COORP,NREF,NCGE,IFLAG)
C AIM :  READ ARRAY MAILAND EXPLAIN IT (!)
</PRE>
<P>
<PRE> 

      SUBROUTINE AELESU(M,LONREE,NCOPNP,NTYE,NTALIR,IAMAII,IAMAIJ,IA,
     +                  NMA)
C AIM : CALCULATE THE ADDRESS OF THE WORD PRECEDING THAT OF THE NEXT ELEMENT
</PRE>
<P>
<PRE> 

       FUNCTION MNEMO(NUMERO,MNEMON,CODE)
C AIM : ENTER THE CORRESPONDENCE BETWEEN A NUMBER AND ITS MNEMONIC
C ---   MNEMO ALWAYS RETURNS THE MNEMONIC CODE IN AN INTEGER
      .
      .
      PARAMETER (NBRE = 24)
      INTEGER CODE
      CHARACTER*4 MNEM(NBRE),MNEMON
      DATA MNEM/'VN  ','DX  ','DY  ','DZ  ','DTG ','DN  ',
     +          'DDXX','DDYY','DDZZ','DDXY','DDYZ','DDZX','DDTG',
     +          'DDNN','DDTN','DLSP','IADN','TRAN','ITI1','ITI2',
     +          'ITI3','ITI4','ROTX','ROTY'/
C     
     .
     .
</PRE>
<P>
<PRE> 

      SUBROUTINE TYPGEL(N,NM1,NM2,NM3)
C AIM :  WRITE CLEARLY IN 3 WORDS NM1, NM2, NM3 THE
C ---    GEOMETRIC TYPE OF THE ELEMENT OF CODE N
</PRE>
<P>

<P>
 <HR SIZE=3 WIDTH="75&#37;"><IMG SRC="../icons/smallmod.gif" WIDTH=211 HEIGHT=50 ALIGN=BOTTOM
	ALT="Modulef"><A NAME=tex2html965 HREF="node43.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/previous_motif.gif"
	ALT="previous"></A><A NAME=tex2html971 HREF="node35.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/up_motif.gif"
	ALT="up"></A><A NAME=tex2html973 HREF="node45.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/next_motif.gif"
	ALT="next"></A><A NAME=tex2html975 HREF="node2.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/contents_motif.gif"
	ALT="contents"></A><A NAME=tex2html976 HREF="node58.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/index_motif.gif"
	ALT="index"></A><A HREF="../Guide2-18/node44.html"><IMG BORDER=0 SRC="../icons/zoom18.gif" ALIGN=BOTTOM
	ALT="[BIG]"></A><A HREF="../Guide2-14/node44.html"><IMG BORDER=0 SRC="../icons/zoom14.gif" ALIGN=BOTTOM
	ALT="[Normal]"></A><A HREF="../Guide2-10/node44.html"><IMG BORDER=0 SRC="../icons/zoom10.gif" ALIGN=BOTTOM
	ALT="[small]"></A><BR>
<B> Next: </B> <A NAME=tex2html974 HREF="node45.html">1.9 DS MILI</A>
<B>Up: </B> <A NAME=tex2html972 HREF="node35.html">1 Description of DS by type</A>
<B> Prev: </B> <A NAME=tex2html966 HREF="node43.html">1.7 DS GEOM</A>
<B><A HREF="node58.html"
	>Index</A></B>
<B><A HREF="node2.html"
	>Contents</A></B>
<BR> <HR>
<P><ADDRESS>
</ADDRESS>
</BODY>
</HTML>