📄 node17.html
字号:
------------------------------ &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE APNOPO &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ? - - ---END</PRE><P><P><P>The resulting data file, created by following the steps above, looks as follows:<P><P><P><PRE>'TEST EXAMPLE 'COURBES 1 $ IMPRE COURBE02(X,Y)= X**2+Y**2-0.25; FIN'POIN ' 1 6 $ IMPRE NPOINT $$ NOP NOREF(NOP) X(NOP). Y(NOP). $ 1 1 0.000000E+00 -.100000E+01 2 1 0.100000E+01 -.100000E+01 3 1 0.100000E+01 0.100000E+01 4 1 0.000000E+00 0.100000E+01 5 1 0.000000E+00 0.500000E+00 6 1 0.000000E+00 -.500000E+00'LIGN ' 1 6 $ IMPRE NDLM $$ NOLIG NOELIG NEXTR1 NEXTR2 NOREFL NFFRON RAISON $ 1 5 1 2 1 0 0.100000E+01 2 10 2 3 1 0 0.100000E+01 3 5 3 4 1 0 0.100000E+01 4 3 4 5 1 0 0.100000E+01 5 6 6 5 2 10 0.100000E+01 6 3 6 1 1 0 0.100000E+01'QUAC ' 1 1 1 6 1 $ IMPRE NIVEAU NUDSD NBRELI NS1L $ LIST OF THE CONTOUR LINES : 1 2 3 4 5 6 5 -1 $ IMAX NQUAD'ADPO ' 1 1 2 1 0 $ IMPRE NIVEA1 NIVEA2 N1 ISET 0 1 $ ISEQ NOESOM'RENC ' 1 2 1 $ IMPRE NIVEA1 NIVEA2 'SAUV ' 1 1 0 $ IMPRE NINOPO NTNOPO TEST.NOPO $ FILE NAME'FIN '</PRE><P><P><P>The resulting mesh is stored in a data structure, under the name you specified when creating the data file, and can be plotted by calling the preprocessor <b> TRNOXX</b>, and specifying the D.S. <b> NOPO</b> file name, forexample:<P><TABLE COLS=3><COL ALIGN=LEFT><COL ALIGN=CENTER><COL ALIGN=LEFT><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP>Call preprocessor: </TD><TD VALIGN=BASELINE ALIGN=CENTER NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> TRNOXX</b> </TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=CENTER NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Specify data structure file name: </TD><TD VALIGN=BASELINE ALIGN=CENTER NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> TEST.NOPO</b> </TD></TR></TABLE><P>The resulting mesh is shown in figure <A HREF="node17.html#figex2">2.2</A>.<P><P><A NAME=1663> </A><IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img142.gif"><BR><STRONG>Figure 2.2:</STRONG> <i> The mesh (plot generated by <b> TRNOXX</b>)</i><A NAME=figex2> </A><BR><P><P>Now that you have successfully generated a mesh, stored in a D.S. <b> NOPO</b>, we may proceed to the next step.<P><P><P><DT>Step 2:<DD> <b> Choice of finite elements</b> <BR><P>Once the D.S. <b> NOPO</b> is created, we are ready to generate the interpolation functions. As discussed earlier, the example under consideration is a thermal type problem. We therefore choose finiteelement from the <b> THER</b>MAL library. We are concerned with a two-dimensional domain, and decide to use two-dimensional (2) straight (D) and curved (C)triangular (TRIA) finite elements with polynomial basis functions of degree 2 (P2 Lagrange), denoted <b> TRIA 2P2D</b> and <b> TRIA 2P2C</b> in the MODULEF library.<P>The preprocessor <b> COMAXX</b> is utilized to perform this step. As for step 1, the data file is created by choosing option CREATE from the main menu, and responding to the various prompts.<P>The input is described in detail below:<P><UL><LI> The name of the data file must be entered (ex: <b> COMA.DATA</b>).<P><LI> The problem dimension must be specified (ex: 2).<P><LI> The use of curved finite elements necessitates the entering of the curve on which the projections will be performed.<P> The curve data is input in exactly the same manner as in step 1.<P><LI> The following information must be input:<P> <UL><LI> number of sub-domains: (ex: 1), <LI> number of reference numbers: (ex: 2), <LI> number of curved lines: (ex: 1), and <LI> curved line number: the reference number designated to each curve, in step 1 (ex: 2). </UL><P><LI> The finite elements used are specified as follows:<P> <UL><LI> The library containing the desired finite elements is specified (ex: <b> THER</b>).<P> A list of the available libraries can be obtained by typing a " ? " in response to the prompt " LIBRARY NAME ? ". <LI> The finite element name is entered for each type of element used <BR> (ex: <b> TRIA 2P2D</b> and <b> TRIA 2P2C</b>) <BR><P> Type " ? " when in doubt. </UL><P><LI> The following file names must be input:<P> <UL><LI> the name of the I.D.S. <b> NOPO</b> file, created in step 1, storing the mesh (ex: <b> TEST.NOPO</b>), <LI> the name of the file containing the resulting O.D.S. <b> MAIL</b> describing the mesh including all information regarding the interpolation (ex: <b> TEST.MAIL</b>), and <LI> the name of the file containing the resulting O.D.S. <b> COOR</b> storing the coordinates of points (ex: <b> TEST.COOR</b>). </UL><P><LI> There are no associated data arrays in the present example (ex: N).<P></UL><P>Once the data file has been created we may executethis module by choosing option EXECUTE from the main menu.<P>The input is described below:<P><UL><LI> The name of the data file must be entered (ex: <b> COMA.DATA</b>).<P><LI> The output parameter ( 0 <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img137.gif"> <b> IMPRE</b> <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img137.gif"> 10) must be specified (ex: 5).<P></UL><P>The data creation procedure and the output corresponding to this step is printed below:<P><P><P><PRE>B$ comaxx **************** EXECUTING COMACO : ( INTERPOLATION ) **************** -- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ? - - ---C************************************************************************** CONVERSATIONAL DATA CREATION OF COMACO************************************************************************** -- NAME OF FILE TO CREATE ?COMA.DATA -- SPACE DIMENSION ?2 -- ARE YOU USING A FUNCTIONS PARSER (YES-NO-?) ?Y PRINT PARAMETER (I) IMPRE NUMBER OF CURVES TO DEFINE (I) NFONC -- IMPRE NFONC ?1 1 THE 1 CURVES ARE NUMBERED FROM 1 TO 1 (YES-NO) ?N ** ENTER THE EXPRESSION F(X,Y) ( = 0. ) ENDING WITH ; ** MAXIMUM LENGHT A64 ( 64 SPACES ) ** EXAMPLE :' X**2 + Y**2 - 4;' OR : X*X+Y*Y-4; CURVE NUMBER 1 ?2 -- CURVE 1 WITH NUMBER 2 :X**2+Y**2-0.25; -- NUMBER OF SUB-DOMAINS ?1 -- NUMBER OF REFERENCE NUMBERS ?2 -- NUMBER OF CURVED LINES ?1 -- NUMBER OF THE CURVED LINE : 1 ?2 ** CALL THE ELEMENT LIBRARY ( IDTECA=0 ) ** -- LIBRARY NAME ( NAME IN A4 OR ? ) ?THER -- NUMBER OF STRAIGHT ELEMENT TYPES IN THE SUB-DOMAIN: 1 ?1 -- NAME OF THE STRAIGHT ELEMENT ( IN A4,1X,A4 OR ? ) ?TRIA 2P2D -- NUMBER OF CURVED ELEMENT TYPES IN THE SUB-DOMAIN : 1 ?1 -- NAME OF THE CURVED ELEMENT ( IN A4,1X,A4 OR ? ) ?TRIA 2P2C -- FILE NAME OF THE INPUT D.S. NOPO ?TEST.NOPO -- FILE NAME OF THE OUTPUT D.S. MAIL ?TEST.MAIL -- FILE NAME OF THE OUTPUT D.S. COOR ?TEST.COOR -- ARE THERE TABLES TO ASSOCIATE (YES-NO) ?N -- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ? - - ---E -- NAME OF THE DATA FILE ? COMA.DATA -- PRINT PARAMETER OF THE D.S. ? 1 ++ OPEN(11,FILE='test.nopo',SPEC='OLD,UNFORMATTED',RECL=0) ++ OPEN(12,FILE='test.mail',SPEC='UNFORMATTED',RECL=0) ++ OPEN(13,FILE='test.coor',SPEC='UNFORMATTED',RECL=0) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE COMACO : END OF MODULE COMACO &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& TABLE M A I 2 -------------- SPACE DIMENSION (NDIM) : 2 NUMBER OF UNKNOWNS OF THE VARIATIONAL FORMULATION (NINCFV) : 1 NUMBER OF TYPES OF ELEMENTS (NTYELM) : 2 NUMBER OF REFERENCES (NNR) : 2 NUMBER OF TYPES OF NODES (NTYNOE) : 1 NUMBER OF TYPES OF POINTS (NTYPOI) : 1 NUMBER OF WORDS OF TABLE MAI3 (LMAI3) : 36 POINTS ARE ONLY DEFINED BY COORDINATES NODES AND POINTS ARE NOT IDENTICAL (NCOPNP) : 0 MAX NUMBER OF WORDS FOR AN ELEMENT (NMMAEL) : 28 NUMBER OF WORDS OF THE TABLE MAIL (LMAIL) : 1421 STORAGE OPTION (NOPFI) : 0 TABLE C O O 2 ------------- TYPE OF THE TABLE COO4 (NTYT) : 2 NUMBER OF ITS INDEX (NINDI) : 2 SPACE DIMENSION (NDIM) : 2 MAX VALUE OF THE SECOND INDEX (M2) : 65 SEGMENTATION CODE (NCODS) : 1 NUMBER OF BLOCS (NBLOC) : 1 TYPE OF THE COORDINATE AXIS (NTACOO) : 1 -- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ? - - ---END</PRE><P><P><P>The resulting data file is shown below:<P><P><P><PRE> 1 $ USING THE FUNCTION PARSER 1 $ IMPRE COURBE02(X,Y)= X**2+Y**2-0.25; FIN 2 1 0 $ NDIM NDSD NBSDC 2 1 $ NNR NBLC 2 THER $ LIBRARY NAME 1 $ NTYED OF SUB-DOMAIN 1 TRIA 2P2D $ THE NAME OF THE STRAIGHT ELEMENTS 1 $ NTYEC OF SUB-DOMAIN 1 TRIA 2P2C $ THE NAME OF THE CURVED ELEMENTS TEST.NOPO $ FILE NAME 0 $ AND LEVEL OF THE D.S. NOPO TEST.MAIL $ FILE NAME 0 $ AND LEVEL OF THE D.S. MAIL TEST.COOR $ FILE NAME 0 $ AND LEVEL OF THE D.S. COOR 0 0 $ NTMAIL NTCOOR</PRE><P><P><P><DT>Step 3:<DD> <b> Entering of physical data, and creation of element arrays</b> <BR><P>Once the mesh is constructed and the interpolation is defined, it is necessary to enter the physical data of the problem under consideration. As discussed earlier in section <A HREF="node15.html#secint_fea">1.1.2</A>, the present example can be solved in terms of the classical thermal problem by specifying some parameters in order to obtain a simplified system, which describes the problem under consideration. Thus, by specifying the following parameters<DL COMPACT><DT><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img143.gif"><DD> = 1,<DT><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img144.gif">⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -