📄 node58.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.3.6 Examples of use</TITLE><LINK REL=STYLESHEET TYPE="text/css" HREF="./Modulef.css" TITLE="Modulef CSS"><meta name="description" value="3.3.6 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=tex2html1567 HREF="node57.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/previous_motif.gif" ALT="previous"></A><A NAME=tex2html1571 HREF="node52.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/up_motif.gif" ALT="up"></A><A NAME=tex2html1573 HREF="node59.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/next_motif.gif" ALT="next"></A><A NAME=tex2html1575 HREF="node2.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/contents_motif.gif" ALT="contents"></A><A NAME=tex2html1576 HREF="node28.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/index_motif.gif" ALT="index"></A><A HREF="../Guide3-18/node58.html"><IMG BORDER=0 SRC="../icons/zoom18.gif" ALIGN=BOTTOM ALT="[BIG]"></A><A HREF="../Guide3-14/node58.html"><IMG BORDER=0 SRC="../icons/zoom14.gif" ALIGN=BOTTOM ALT="[Normal]"></A><A HREF="../Guide3-10/node58.html"><IMG BORDER=0 SRC="../icons/zoom10.gif" ALIGN=BOTTOM ALT="[small]"></A><BR><B> Next: </B> <A NAME=tex2html1574 HREF="node59.html">3.4 Method for a hexahedral topology</A><B>Up: </B> <A NAME=tex2html1572 HREF="node52.html">3.3 Extrusion method (cylindrical topology)</A><B> Prev: </B> <A NAME=tex2html1568 HREF="node57.html">3.3.5 The preprocessor MA23XX</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=SECTION00636000000000000000>3.3.6 Examples of use</A></H2><P><P><P>Some examples of calling module <b> MA2D3E</b> in batch and calls via preprocessor <b> MA23XX</b> are given below.<P>First, an example with <b> TTYPE = -1</b>, and the same with <b> TTYPE = - 2 </b>, i.e., a step-by-step definitionof the section starting from the basis via the general matrix, MAT, or via function <b> XYZ23</b>.<P>The calling programs, the output obtained during execution of case 1, and a plot of the resulting mesh (figure <A HREF="node58.html#figma2312">3.14</A>) is shown below.<P><P><P>Definition via MAT (TTYPE = -1):<P><PRE>C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++C EXAMPLE OF CALLING MODULE MA2D3E : TEST WITH A MATRIXC +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ PARAMETER (LM = 100000 , MAXFON = 4 , MAXSEC = 9 , MAXDES = 10 ) COMMON M(LM) INTEGER NUCOU(3,MAXFON),DESREF(5,MAXDES) REAL MAT(4,4,MAXFON) REAL ZINT(0:MAXSEC) LOGICAL BFONC,COLLER,RAPIDE EXTERNAL XYZ23,XYZ33C IMPRE = 5 CALL INITIS(M,LM,IMPRE,0)C ------ THE I.D.S. NOPO AND THE O.D.S. NOPO NFNOPO = 10 NINOPO = 1 CALL OUVRIR(NFNOPO,'nopo2d.ma23.2', + 'OLD,UNFORMATTED',0,IOSTAT) NFNOPS = 11 NINOPS = 2 CALL OUVRIR(NFNOPS,'nopo3d.ma23.1','UNFORMATTED',0,IOSTAT) C -----------------------------------------------C ------ GEOMETRIC TRANSFORMATIONS ------ C ----------------------------------------------- C --- THE BASIS: BFONC = .FALSE. ZINT(0) = -1. NBCOU = 4C --- MATRICES BETWEENC THE BASIS AND THE ONE ABOVE DEFINED BY MATC TO CONSTRUCT SECTIONS 1 TO 4 NBFONC = 4 DO ISEC=1,4 NUCOU(1,ISEC) = ISEC - 1 NUCOU(2,ISEC) = ISEC NUCOU(3,ISEC) = -1 ENDDO C DO ISEC=1,4 MAT(1,1,ISEC) = 0.8 MAT(1,2,ISEC) = 0. MAT(1,3,ISEC) = 0. MAT(1,4,ISEC) = 0. MAT(2,1,ISEC) = 0. MAT(2,2,ISEC) = 0.8 MAT(2,3,ISEC) = 0. MAT(2,4,ISEC) = 0. MAT(3,1,ISEC) = 0. MAT(3,2,ISEC) = 0. MAT(3,3,ISEC) = 1. MAT(3,4,ISEC) = 1. + ISEC MAT(4,1,ISEC) = 0. MAT(4,2,ISEC) = 0. MAT(4,3,ISEC) = 0. MAT(4,4,ISEC) = 1. ENDDO C -----------------------------------------------C ------ DEFINITION OF THE REFERENCES ------ C ----------------------------------------------- NBDES = 2C --- BASIS: FACES EDGE AND POINTS SET TO REFERENCE 1 : IDES = 1 DESREF(1,IDES) = 0 DESREF(2,IDES) = 0 DESREF(3,IDES) = 0 DESREF(4,IDES) = -100 DESREF(5,IDES) = 1 C ---C --- ABOVE: FACES EDGES AND POINTS SET TO REFERENCE 2 : IDES = 2 DESREF(1,IDES) = 9 DESREF(2,IDES) = 9 DESREF(3,IDES) = 0 DESREF(4,IDES) = -100 DESREF(5,IDES) = 2 C -----------------------------------------------C ------ CALL THE MODULE ------ C ----------------------------------------------- EPS = 0.001 COLLER = .FALSE. RAPIDE = .TRUE.C CALL MA2D3E(M,M,NFNOPO,NINOPO,NFNOPS,NINOPS,NBCOU,EPS,NBFONC, + NUCOU,MAT,COLLER,RAPIDE,DESREF,NBDES,ZINT,BFONC, + XYZ23,XYZ33) END</PRE><P><P><A NAME=4953> </A><IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img153.gif"><BR><STRONG>Figure 3.14:</STRONG> <i> Example <b> MA2D3E</b> (step-by-step definition via MAT or XYZ23)</i><A NAME=figma2312> </A><BR><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 : 16/10/91 AUTHOR : dutoit ++ OPEN(10,FILE='nopo2d.ma23.2',SPEC='OLD,UNFORMATTED',RECL=0) ++ OPEN(11,FILE='nopo3d.ma23.1',SPEC='UNFORMATTED',RECL=0) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE MA2D3E : &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& HEIGHT OF THE BOTTOM BASIS : -0.1000000E+01 NUMBER OF GEOMETRIC TRANSFORMATIONS : 4 LIST OF TRANSFORMATIONS : SECTION 0 A 1 TYPE OF TRANSFORMATION -1 SECTION 1 A 2 TYPE OF TRANSFORMATION -1 SECTION 2 A 3 TYPE OF TRANSFORMATION -1 SECTION 3 A 4 TYPE OF TRANSFORMATION -1 NUMBER OF TRANSFERS OF NUMBERS : 2 LIST OF TRANSFERS : SECTION 0 A 0 TYPE OF TRANSFER 0 : -100 GIVES 1 SECTION 9 A 9 TYPE OF TRANSFER 0 : -100 GIVES 2 MESH OBTAINED : SAVE D.S. TABLE N O P 2 -------------- MESH CHARACTERISTICS : SPACE DIMENSION (NDIM ) : 3 MAXIMUM NUMBER OF REFERENCES (NDSR ) : 1 MAXIMUM NUMBER OF SUB-DOMAINS (NDSD ) : 1 NODES AND POINTS COINCIDE (NCOPNP) : 1 NUMBER OF ELEMENTS (NE ) : 880 NUMBER OF PENTAHEDRA (NPENT) : 880 NUMBER OF BOUNDARY ELEMENTS (NEF ) : 72 NUMBER OF NODES (NOE ) : 710 NUMBER OF NODES PER SEGMENT (NO EXTREMITIES) : 0 TYPE OF COORDINATE VALUES (NTYCOO) : REEL1MOT MAX DIFFERENCE +1 BETWEEN 2 NODES OF AN ELEMENT : 271 NUMBER OF COARSE ELEMENTS (NBEGM ) : 0 NUMBER OF WORDS IN ARRAY NOP5 (LNOP5 ) : 13420 REFERENCE AXIS X,Y,Z (NTACOO) : 1 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MA2D3E &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&</PRE><P><P><P>A definition via XYZ23 (TTYPE = -2):<P><PRE>C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++C EXAMPLE OF CALLING MODULE MA2D3E : TEST WITH XYZ23 THROUGHOUTC +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ PARAMETER (LM = 100000 , MAXFON = 4 , MAXSEC = 9 , MAXDES = 10 ) COMMON M(LM) INTEGER NUCOU(3,MAXFON),DESREF(5,MAXDES) REAL MAT(4,4,MAXFON) REAL ZINT(0:MAXSEC) LOGICAL BFONC,COLLER,RAPIDE EXTERNAL XYZ23,XYZ33C IMPRE = 5 CALL INITIS(M,LM,IMPRE,0)C ------ THE I.D.S. NOPO AND THE O.D.S. NOPO ------ NFNOPO = 10 NINOPO = 1 CALL OUVRIR(NFNOPO,'nopo2d.ma23.2', + 'OLD,UNFORMATTED',0,IOSTAT) NFNOPS = 11 NINOPS = 2 CALL OUVRIR(NFNOPS,'nopo3d.ma23.2','UNFORMATTED',0,IOSTAT) C -----------------------------------------------C ------ GEOMETRIC TRANSFORMATIONS ------ C ----------------------------------------------- C --- THE BASIS : BFONC = .FALSE. ZINT(0) = -1. NBCOU = 4C --- FUNCTIONS XYZ23C BETWEEN THE BASE THE ONE ABOVE DEFINED BY XYZ23C TO CONSTRUCT SECTIONS 1 TO 4 NBFONC = 4 DO ISEC=1,4 NUCOU(1,ISEC) = ISEC - 1 NUCOU(2,ISEC) = ISEC NUCOU(3,ISEC) = -2 ENDDO C C -----------------------------------------------C ------ DEFINITION OF THE REFERENCES ------ C ----------------------------------------------- NBDES = 2C --- BASIS : FACES EDGES AND POINTS SET TO REFERENCE 1 : IDES = 1 DESREF(1,IDES) = 0 DESREF(2,IDES) = 0 DESREF(3,IDES) = 0 DESREF(4,IDES) = -100 DESREF(5,IDES) = 1 C ---C --- ABOVE : FACES EDGES AND POINTS SET TO REFERENCE 2 : IDES = 2 DESREF(1,IDES) = 9 DESREF(2,IDES) = 9 DESREF(3,IDES) = 0 DESREF(4,IDES) = -100 DESREF(5,IDES) = 2 C -----------------------------------------------C ------ CALL THE MODULE ------ C ----------------------------------------------- EPS = 0.001 COLLER = .FALSE. RAPIDE = .TRUE.C CALL MA2D3E(M,M,NFNOPO,NINOPO,NFNOPS,NINOPS,NBCOU,EPS,NBFONC, + NUCOU,MAT,COLLER,RAPIDE,DESREF,NBDES,ZINT,BFONC, + XYZ23,XYZ33) END CC +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ C SUBROUTINE XYZ23C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ SUBROUTINE XYZ23(I,NUM,NREF,X,Y,Z,X3,Y3,Z3) REAL*4 DZ(4) DATA DZ/ 2., 5. , 9., 14./C X3 = 0.8 **I * X Y3 = 0.8 **I * Y Z3 = Z + DZ(I) END</PRE><P>An example with <b> TTYPE = -3</b> for sections 5 to 9 and <b> TTYPE = -6</b> for the first sections. Consequently, we have:<UL><LI> an interpolation, of sides given in array ZINT between the basis (z= ZINT(0)) and a virtual sectionof side 10 described by XYZ23, for sections 1 to 5<LI> the construction, section by section, of sections 6 to 9 via XYZ33.</UL> The references of the 3D mesh are attributed via array DESREF(.,.) (see the calling program forthe specification of the transfers and figure <A HREF="node58.html#figma2336">3.15</A> for the result obtained).<P>The calling program, the output obtained during execution and a plot of the resulting mesh is shown below.<P><PRE>C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++C EXAMPLE OF CALLING MODULE MA2D3E : TEST NUMBER 3C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ PARAMETER (LM = 100000 , MAXFON = 2 , MAXSEC = 9 , MAXDES = 10 ) COMMON M(LM) INTEGER NUCOU(3,MAXFON),DESREF(5,MAXDES) REAL MAT(4,4,MAXFON) REAL ZINT(0:MAXSEC) LOGICAL BFONC,COLLER,RAPIDE EXTERNAL XYZ23,XYZ33C IMPRE = 5 CALL INITIS(M,LM,IMPRE,0)C ------ THE I.D.S NOPO AND THE O.D.S NOPO ------ NFNOPO = 10 NINOPO = 1 CALL OUVRIR(NFNOPO,'nopo2d.ma23.2', + 'OLD,UNFORMATTED',0,IOSTAT) NFNOPS = 11 NINOPS = 2 CALL OUVRIR(NFNOPS,'nopo3d.ma23.36','UNFORMATTED',0,IOSTAT) C -----------------------------------------------C ------ GEOMETRICAL TRANSFORMATIONS ------ C ----------------------------------------------- C --- LA BASE : BFONC = .FALSE. ZINT(0) = -1.C --- THE INTERMEDIARY SIDES : NBCOU = 9 ZINT(1) = 0. ZINT(2) = 2. ZINT(3) = 4. ZINT(4) = 6. ZINT(5) = 7. C --- A FUNCTION : A MANUAL INTERPOLATION C BETWEEN THE BASIS AND THE ONE ABOVE DEFINED BY XYZ23 C FOR SIDE 10. ( SEE XYZ23 ) C TO CONSTRUCT SECTIONS 1 TO 5 NBFONC = 2 C IFONC = 1 NUCOU(1,IFONC) = 0 NUCOU(2,IFONC) = 5 NUCOU(3,IFONC) = -6 C --- A FUNCTION : LOCAL DEFINITION BY XYZ33 C TO CONSTRUCT SECTIONS 6 TO 9 IFONC = 2 NUCOU(1,IFONC) = 5 NUCOU(2,IFONC) = 9 NUCOU(3,IFONC) = -3C -----------------------------------------------C ------ DEFINITION OF THE REFERENCES ------ C ----------------------------------------------- ⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -