📄 node51.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.2.6 Some utilization examples</TITLE><LINK REL=STYLESHEET TYPE="text/css" HREF="./Modulef.css" TITLE="Modulef CSS"><meta name="description" value="3.2.6 Some utilization examples"><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=tex2html1479 HREF="node50.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/previous_motif.gif" ALT="previous"></A><A NAME=tex2html1483 HREF="node45.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/up_motif.gif" ALT="up"></A><A NAME=tex2html1485 HREF="node52.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/next_motif.gif" ALT="next"></A><A NAME=tex2html1487 HREF="node2.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/contents_motif.gif" ALT="contents"></A><A NAME=tex2html1488 HREF="node28.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/index_motif.gif" ALT="index"></A><A HREF="../Guide3-18/node51.html"><IMG BORDER=0 SRC="../icons/zoom18.gif" ALIGN=BOTTOM ALT="[BIG]"></A><A HREF="../Guide3-14/node51.html"><IMG BORDER=0 SRC="../icons/zoom14.gif" ALIGN=BOTTOM ALT="[Normal]"></A><A HREF="../Guide3-10/node51.html"><IMG BORDER=0 SRC="../icons/zoom10.gif" ALIGN=BOTTOM ALT="[small]"></A><BR><B> Next: </B> <A NAME=tex2html1486 HREF="node52.html">3.3 Extrusion method (cylindrical topology)</A><B>Up: </B> <A NAME=tex2html1484 HREF="node45.html">3.2 Multi-block method</A><B> Prev: </B> <A NAME=tex2html1480 HREF="node50.html">3.2.5 Preprocessor COLIXX</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=SECTION00626000000000000000>3.2.6 Some utilization examples</A></H2><P><P><P>The following examples are given:<P><UL><LI> an example of a call to <b> COLIB2</b> with the calling program and a plot of the resulting mesh;<LI> an example of a call to <b> COLIBR</b> with the calling program, the data file and a plot of the mesh created in the 2D and 3D case;<LI> an example of a call to <b> COLIBH</b> with the calling program and a plot of the mesh created;<LI> another example of a mesh created via several calls to <b> COLIBH </b> and the gluing together ofmeshes thus obtained.</UL><P>The calling program of <b> COLIB2</b>, for a pentahedral element, is given below. Figure <A HREF="node51.html#figpentacol2">3.3</A> showsthe mesh obtained for <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img129.gif"> and <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img130.gif">, where <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img118.gif"> and <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img119.gif"> designate the number of points on the edges of the triangular faces and the number of these points on the other edges, respectively.<P><PRE>C ++++++++++++++++++++++++++++++++++++++++++++++C TEST COLIBR ( COLIB2 WITH DATA FILE ) C ( A 2 DIMENSIONAL EXAMPLE )C ++++++++++++++++++++++++++++++++++++++++++++++ PARAMETER ( LM = 300000 ) COMMON M(LM)C IMPRE = -2 NNN = 0 CALL INITIS(M,LM,IMPRE,NNN) NFNOPO = 10 CALL TRUNIT (NFNOPO) CALL OUVRIS(NFNOPO,'tunnel.nopo','UNKNOWN,UNFORMATTED',2000000) C NFDATA = 11 CALL TRUNIT (NFDATA) CALL OUVRIS(NFDATA,'tunnel.data','UNKNOWN,FORMATTED',2000000) C CALL LECTEU(NFDATA) CALL COLIBR(M,NFNOPO,NINOPO) STOP END</PRE><P><P><A NAME=4572> </A><IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img131.gif"><BR><STRONG>Figure 3.3:</STRONG> <i> Mesh of a pentahedron via <b> COLIB2</b></i><A NAME=figpentacol2> </A><BR><P><P>The calling program for <b> COLIBR</b> is shown below for the case corresponding to a 2D mesh constructionfrom 17 triangular and quadrilateral blocks. The resulting mesh is shown in figure <A HREF="node51.html#figtunnel">3.4</A>.<P><PRE>C ++++++++++++++++++++++++++++++++++++++++++++++C TEST COLIBR ( COLIB2 WITH DATA FILE ) C ( A 2 DIMENSIONAL EXAMPLE )C ++++++++++++++++++++++++++++++++++++++++++++++ PARAMETER ( LM = 300000 ) COMMON M(LM)C IMPRE = -2 NNN = 0 CALL INITIS(M,LM,IMPRE,NNN) NFNOPO = 10 CALL TRUNIT (NFNOPO) CALL OUVRIS(NFNOPO,'tunnel.nopo','UNKNOWN,UNFORMATTED',2000000) C NFDATA = 11 CALL TRUNIT (NFDATA) CALL OUVRIS(NFDATA,'tunnel.data','UNKNOWN,FORMATTED',2000000) C CALL LECTEU(NFDATA) CALL COLIBR(M,NFNOPO,NINOPO) STOP END</PRE><P><P><P>The data file required is the following:<P><P><P><PRE>$ NBS NA NBGRO IMPRE NBFR NDIM $ 19 35 17 2 0 2$ NOP NOREF X Y Z ? $ 1 1 0.000000E+00 0.000000E+00 0.000000E+00 2 2 0.110000E+02 0.000000E+00 0.000000E+00 3 3 0.110000E+02 0.900000E+01 0.000000E+00 4 4 0.110000E+02 0.180000E+02 0.000000E+00 5 5 0.200000E+01 0.180000E+02 0.000000E+00 6 6 0.000000E+00 0.180000E+02 0.000000E+00 7 7 0.000000E+00 0.400000E+01 0.000000E+00 8 8 0.800000E+01 0.400000E+01 0.000000E+00 9 9 0.800000E+01 0.800000E+01 0.000000E+00 10 10 0.600000E+01 0.110000E+02 0.000000E+00 11 11 0.200000E+01 0.110000E+02 0.000000E+00 12 12 0.000000E+00 0.800000E+01 0.000000E+00 13 13 0.200000E+01 0.600000E+01 0.000000E+00 14 14 0.550000E+01 0.600000E+01 0.000000E+00 15 15 0.550000E+01 0.750000E+01 0.000000E+00 16 16 0.490000E+01 0.900000E+01 0.000000E+00 17 17 0.290000E+01 0.910000E+01 0.000000E+00 18 18 0.200000E+01 0.750000E+01 0.000000E+00 19 19 0.400000E+01 0.750000E+01 0.000000E+00$ NA IEX1 IEX2 IDEC NP NOREF $ 1 1 2 0 5 1 2 2 3 0 5 1 3 3 4 0 5 1 4 4 5 0 5 1 5 5 6 0 5 1 6 6 12 0 4 1 7 12 7 0 5 2 8 7 1 0 4 1 9 7 8 0 5 2 10 8 2 0 4 0 11 8 9 0 5 2 12 9 10 0 5 2 13 9 3 0 4 0 14 10 4 0 4 0 15 10 11 0 5 2 16 12 11 0 5 2 17 11 5 0 4 0 18 13 7 0 5 0 19 14 8 0 5 0 20 15 9 0 5 0 21 16 10 0 5 0 22 17 11 0 5 0 23 18 12 0 5 0 24 13 14 0 5 3 25 14 15 0 5 3 26 15 16 1 5 3 5.45 7.75 0. 5.4 8. 0. 5.3 8.25 0. 5.2 8.50 0. 5.1 8.75 0. 27 16 17 1 5 3 4.5 9.3 0. 4.25 9.45 0. 4.0 9.5 0. 3.5 9.45 0. 3.25 9.3 0. 28 17 18 1 5 3 2.60 8.75 0. 2.40 8.50 0. 2.25 8.25 0. 2.10 8.0 0. 2.05 7.75 0. 29 18 13 0 5 3 30 19 13 0 5 0 31 19 14 0 5 0 32 19 15 0 5 0 33 19 16 0 5 0 34 19 17 0 5 0 35 19 18 0 5 0$ COARSE ELEMENTS $ $ I NCGE VERTICES NCGE NDSD $ 1 4 1 2 8 7 4 1 2 4 2 3 9 8 4 1 3 4 3 4 10 9 4 1 4 4 10 4 5 11 4 1 5 4 11 5 6 12 4 1 6 4 7 13 18 12 4 2 7 4 7 8 14 13 4 2 8 4 8 14 15 9 4 2 9 4 15 9 10 16 4 2 10 4 16 10 11 17 4 2 11 4 18 17 11 12 4 2 12 3 13 14 19 3 3 13 3 14 15 19 3 3 14 3 15 16 19 3 3 15 3 16 17 19 3 3 16 3 19 18 17 3 3 17 3 13 19 18 3 3</PRE><P><P><A NAME=4590> </A><IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img132.gif"><BR><STRONG>Figure 3.4:</STRONG> <i> 2D mesh via <b> COLIBR</b></i><A NAME=figtunnel> </A><BR><P><P>The calling program for <b> COLIBR</b> is given below for the case corresponding tothe construction of a 3D mesh from 12 pentahedral and hexahedral blocks.The resulting mesh is shown in figure <A HREF="node51.html#figparttun">3.5</A>.<P><PRE>C ++++++++++++++++++++++++++++++++++++++++++++++C TEST COLIBR ( COLIB2 WITH DATA FILE ) C ( A 3 DIMENSIONAL EXAMPLE )C ++++++++++++++++++++++++++++++++++++++++++++++ PARAMETER ( LM = 300000 ) COMMON M(LM)C IMPRE = -2 NNN = 0 CALL INITIS(M,LM,IMPRE,NNN) NFNOPO = 10 CALL TRUNIT (NFNOPO) CALL OUVRIS(NFNOPO,'parttun.nopo','UNKNOWN,UNFORMATTED',2000000) C NFDATA = 11 CALL TRUNIT (NFDATA) CALL OUVRIS(NFDATA,'parttun.data','UNKNOWN,FORMATTED',2000000) C CALL LECTEU(NFDATA) CALL COLIBR(M,NFNOPO,NINOPO) STOP END</PRE><P><P><P>The data file required is the following:<P><P><P><PRE>$ NBS NA NBGRO IMPRE NBFR NDIM $ 26 61 12 2 0 3$ NOP NOREF X Y Z ? $ 1 7 0.000000E+00 0.400000E+01 0.000000E+00 2 8 0.800000E+01 0.400000E+01 0.000000E+00 3 9 0.800000E+01 0.800000E+01 0.000000E+00 4 10 0.600000E+01 0.110000E+02 0.000000E+00 5 11 0.200000E+01 0.110000E+02 0.000000E+00 6 12 0.000000E+00 0.800000E+01 0.000000E+00 7 13 0.200000E+01 0.600000E+01 0.000000E+00 8 14 0.550000E+01 0.600000E+01 0.000000E+00 9 15 0.550000E+01 0.750000E+01 0.000000E+00 10 16 0.490000E+01 0.900000E+01 0.000000E+00 11 17 0.290000E+01 0.910000E+01 0.000000E+00 12 18 0.200000E+01 0.750000E+01 0.000000E+00 13 19 0.400000E+01 0.750000E+01 0.000000E+00 14 7 0.000000E+00 0.400000E+01 1.000000E+01 15 8 0.800000E+01 0.400000E+01 1.000000E+01 16 9 0.800000E+01 0.800000E+01 1.000000E+01 17 10 0.600000E+01 0.110000E+02 1.000000E+01 18 11 0.200000E+01 0.110000E+02 1.000000E+01 19 12 0.000000E+00 0.800000E+01 1.000000E+01 20 13 0.200000E+01 0.600000E+01 1.000000E+01 21 14 0.550000E+01 0.600000E+01 1.000000E+01 22 15 0.550000E+01 0.750000E+01 1.000000E+01 23 16 0.490000E+01 0.900000E+01 1.000000E+01 24 17 0.290000E+01 0.910000E+01 1.000000E+01 25 18 0.200000E+01 0.750000E+01 1.000000E+01 26 19 0.400000E+01 0.750000E+01 1.000000E+01$ NA IEX1 IEX2 IDEC NP NOREF $ 1 6 1 0 5 2 2 1 2 0 5 2 3 2 3 0 5 2⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -