📄 node63.html
字号:
<!DOCTYPE HTML PUBLIC "-//IETF//DTD HTML 3.2 Final//FR"><!-- Converted with LaTeX2HTML 95.1 (Fri Jan 20 1995) --><!-- by Nikos Drakos (nikos@cbl.leeds.ac.uk), CBLU, University of Leeds --><!-- Modified Simulog 03/97 --><HTML><HEAD><TITLE>3.4.4 Examples of use</TITLE><LINK REL=STYLESHEET TYPE="text/css" HREF="./Modulef.css" TITLE="Modulef CSS"><meta name="description" value="3.4.4 Examples of use"><meta name="keywords" value="Guide3"><meta name="resource-type" value="document"><meta name="distribution" value="global"></HEAD><BODY BGCOLOR="#FFFFFF"><P> <IMG SRC="../icons/smallmod.gif" WIDTH=211 HEIGHT=50 ALIGN=BOTTOM ALT="Modulef"><A NAME=tex2html1627 HREF="node62.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/previous_motif.gif" ALT="previous"></A><A NAME=tex2html1631 HREF="node59.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/up_motif.gif" ALT="up"></A><A NAME=tex2html1633 HREF="node64.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/next_motif.gif" ALT="next"></A><A NAME=tex2html1635 HREF="node2.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/contents_motif.gif" ALT="contents"></A><A NAME=tex2html1636 HREF="node28.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/index_motif.gif" ALT="index"></A><A HREF="../Guide3-18/node63.html"><IMG BORDER=0 SRC="../icons/zoom18.gif" ALIGN=BOTTOM ALT="[BIG]"></A><A HREF="../Guide3-14/node63.html"><IMG BORDER=0 SRC="../icons/zoom14.gif" ALIGN=BOTTOM ALT="[Normal]"></A><A HREF="../Guide3-10/node63.html"><IMG BORDER=0 SRC="../icons/zoom10.gif" ALIGN=BOTTOM ALT="[small]"></A><BR><B> Next: </B> <A NAME=tex2html1634 HREF="node64.html">4 List of the modules</A><B>Up: </B> <A NAME=tex2html1632 HREF="node59.html">3.4 Method for a hexahedral topology</A><B> Prev: </B> <A NAME=tex2html1628 HREF="node62.html">3.4.3 Calling of module GEL3D1</A><B><A HREF="node28.html" >Index</A></B><B><A HREF="node2.html" >Contents</A></B><HR SIZE=3 WIDTH="75%"><H2><A NAME=SECTION00644000000000000000>3.4.4 Examples of use</A></H2><P>As an illustration, we now give three examples, each of which is accompanied by a listing of the calling program, the output obtained during execution (with print parameter: IMPRE = 5) and a plot of thefinal mesh.<P>The first example is the mesh of a cube into tetrahedra. The data is equidistant (step <b>h=0.4</b>) in the 3 directions, with a single sub-domain number (- JOPTSD ) assigned. The final mesh is shown in figure <A HREF="node63.html#figgel1">3.31</A>.<P>The calling program is the following:<P><P><P><PRE>C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++C EXAMPLE 1 : MODULE GEL3D1C --------- ( PARALLELEPEDICAL DOMAIN ==> TETRAHEDRA )C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ PARAMETER ( LM = 18000 ) COMMON M(LM) REAL H(1) INTEGER NUMAT(1),ITYP(1) EXTERNAL F1,F2,F3C IMPRE = 5 CALL INITIS(M,LM,IMPRE,0) C ------ THE RESULTING FILE NOPO : NFNOPO = 10 NINOPO = 0 CALL OUVRIR(NFNOPO,'GEL1.NOPO','UNFORMATTED',0,IOSTAT) C ------ THE GEOMETRY : N1 = 8 N2 = 7 NB = 4 H(1) = .4 IOPTM = 0 NOPTDE = 5 C ------ THE NUMBERS : JOPTSD = -1C NFAUX = 0 C ------ CALL THE MODULE ------ CALL GEL3D1(M,IOPTM,JOPTSD,NOPTDE,N1,N2,NB,H,H,H,F1,F2,F3, + NUMAT,ITYP,NFNOPO,NINOPO,NFAUX) END</PRE><P><P><P>The listing of the output obtained during execution is given below:<P><P><P><PRE> M M OOO DDDD U U L EEEEE FFFFF MM MM O O D D U U L E F M M M O O D D U U L EEEE FFFF M M O O D D U U L E F M M OOO DDDD UUU LLLLL EEEEE F VERSION 91 DATE : 18/10/91 AUTHOR : dutoit ++ OPEN(10,FILE='gel1.nopo',SPEC='UNFORMATTED',RECL=0) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE GEL3D1 : &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& GENERATION OF 840 TETRAHEDRA AND 280 NODES LENGTH OF M (OVER-ESTIMATION) TO GENERATE NOPO : 17238 NUMBER OF POINTS IN X (NX) : 8 NUMBER OF POINTS IN Y (NY) : 7 NUMBER OF BANDS (NB) : 4 GENERATION OPTION (IOPT) : 0 NUMBERING OPTION OF THE D.S. (JOPT) : -1 SPLITTING OPTION (NOPT) : 5 STEP IN X (X(1)) : 0.4000E+00 STEP IN Y (Y(1)) : 0.4000E+00 STEP IN Z (H(1)) : 0.4000E+00 LENGTH CALCULATED OF M FOR D.S. NOPO 14726 END OF GENERATION OF NOPO ON FILE : 10 TABLE N O P 2 -------------- MESH CHARACTERISTICS : SPACE DIMENSION (NDIM ) : 3 MAXIMUM NUMBER OF REFERENCES (NDSR ) : 26 MAXIMUM NUMBER OF SUB-DOMAINS (NDSD ) : 1 NODES AND POINTS COINCIDE (NCOPNP) : 1 NUMBER OF ELEMENTS IN THE MESH (NE ) : 840 NUMBER OF TETRAHEDRA (NTET ) : 840 NUMBER OF BOUNDARY ELEMENTS (NEF ) : 640 NUMBER OF NODES (NOE ) : 280 NUMBER OF NODES PER SEGMENT (NO EXTREMITIES) : 0 TYPE OF COORDINATE VALUES (NTYCOO) : REEL1MOT MAX DIFFERENCE +1 BETWEEN 2 NODES OF AN ELEMENT : 64 NUMBER OF COARSE ELEMENTS (NBEGM ) : 0 NUMBER OF WORDS IN ARRAY NOP5 (LNOP5 ) : 13808 REFERENCE AXIS X,Y,Z (NTACOO) : 1 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE GEL3D1 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&</PRE><P><P><A NAME=5182> </A><IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img182.gif"><BR><STRONG>Figure 3.31:</STRONG> <i> First example of a mesh generated by <b> GEL3D1</b></i><A NAME=figgel1> </A><BR><P><P>The second example is a mesh of a cube into hexahedra. The data consists of several bands of varying heights (array H(.)).Several sub-domain numbers are assigned (arrays ITYP and NUMAT). The final mesh is shown infigure <A HREF="node63.html#figgel2">3.32</A>.<P>The calling program is the following:<P><P><P><PRE>C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++C EXAMPLE 2 : MODULE GEL3D1C --------- ( PARALLELEPEDICAL DOMAIN ==> HEXAHEDRA ) C ( DIFFERENT HEIGHT BANDS, SEVERAL MEDIA )C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ PARAMETER ( LM = 18000 ) COMMON M(LM) REAL X(7),Y(7),H(6) INTEGER NUMAT(6,6,2),ITYP(5) EXTERNAL F1,F2,F3C IMPRE = 5 CALL INITIS(M,LM,IMPRE,0) C ------ THE RESULTING FILE NOPO : NFNOPO = 10 NINOPO = 0 CALL OUVRIR(NFNOPO,'GEL2.NOPO','UNFORMATTED',0,IOSTAT) C ------ THE GEOMETRY : N1 = 7 X(1) = 0. X(2) = 1. X(3) = 1.5 X(4) = 2. X(5) = 3.5 X(6) = 5.5 X(7) = 7.5 N2 = 7 Y(1) = 0. Y(2) = 1. Y(3) = 1.5 Y(4) = 2. Y(5) = 3. Y(6) = 4. Y(7) = 5. NB = 5 H(1) = 0. H(2) = 1. H(3) = 1.5 H(4) = 2. H(5) = 3. H(6) = 4. IOPTM = 2 NOPTDE = 0 C ------ THE NUMBERS : JOPTSD = 2C ------ THE LAYER TYPES : DO 1 I=1,5 ITYP(I) = 1 1 CONTINUE ITYP(2) = 2 ITYP(3) = 2 C ------ THE SUB-DOMAINS : DO 2 I=1,6 DO 2 J=1,6 DO 2 K=1,2 NUMAT(I,J,K) = 1 2 CONTINUE NUMAT(4,4,1) = 2 NUMAT(4,4,2) = 3C NFAUX = 0 C ------ CALL THE MODULE ------ CALL GEL3D1(M,IOPTM,JOPTSD,NOPTDE,N1,N2,NB,X,Y,H,F1,F2,F3, + NUMAT,ITYP,NFNOPO,NINOPO,NFAUX) END</PRE><P><P><P>The listing of the output during execution is given below:<P><P><P><PRE> M M OOO DDDD U U L EEEEE FFFFF MM MM O O D D U U L E F M M M O O D D U U L EEEE FFFF M M O O D D U U L E F M M OOO DDDD UUU LLLLL EEEEE F VERSION 91 DATE : 21/10/91 AUTHOR : dutoit ++ OPEN(10,FILE='gel2.nopo',SPEC='UNFORMATTED',RECL=0) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE GEL3D1 : &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& GENERATION OF 180 HEXAHEDRA AND 294 NODES LENGTH OF M (OVER-ESTIMATION) TO GENERATE NOPO : 6684 NUMBER OF POINTS IN X (NX) : 7 NUMBER OF POINTS IN Y (NY) : 7 NUMBER OF BANDS (NB) : 5 GENERATION OPTION (IOPT) : 2 NUMBERING OPTION OF THE D.S. (JOPT) : 2 SPLITTING OPTION (NOPT) : 0 ABSCISSIS OF POINTS (X(*)) : 0.0000E+00 0.1000E+01 0.1500E+01 0.2000E+01 0.3500E+01 0.5500E+01 0.7500E+01 ORDINATES OF POINTS (Y(*)) : 0.0000E+00 0.1000E+01 0.1500E+01 0.2000E+01 0.3000E+01 0.4000E+01 0.5000E+01 HEIGHTS OF POINTS (H(*)) : ⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -