📄 node63.html
字号:
0.0000E+00 0.1000E+01 0.1500E+01 0.2000E+01 0.3000E+01 0.4000E+01 TYPE OF BANDS : 1 2 2 1 1 ARRAY NUMAT OF TYPE 1 : 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ARRAY NUMAT OF TYPE 2 : 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 LENGTH CALCULATED OF M FOR D.S. NOPO 6684 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 ) : 3 NODES AND POINTS COINCIDE (NCOPNP) : 1 NUMBER OF ELEMENTS IN THE MESH (NE ) : 180 NUMBER OF HEXAHEDRA (NHEX ) : 180 NUMBER OF BOUNDARY ELEMENTS (NEF ) : 132 NUMBER OF NODES (NOE ) : 294 NUMBER OF NODES PER SEGMENT (NO EXTREMITIES) : 0 TYPE OF COORDINATE VALUES (NTYCOO) : REEL1MOT MAX DIFFERENCE +1 BETWEEN 2 NODES OF AN ELEMENT : 57 NUMBER OF COARSE ELEMENTS (NBEGM ) : 0 NUMBER OF WORDS IN ARRAY NOP5 (LNOP5 ) : 5724 REFERENCE AXIS X,Y,Z (NTACOO) : 1 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE GEL3D1 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&</PRE><P><P><A NAME=5199> </A><IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img183.gif"><BR><STRONG>Figure 3.32:</STRONG> <i> A second example of a mesh generated by <b> GEL3D1</b></i><A NAME=figgel2> </A><BR><P><P>The third example is a mesh of a cube into pentahedra. The data isinput via functions F1, F2 and F3, and a different sub-domain number is assigned to the layers (array NUMAT). The mesh obtained is shown in figure <A HREF="node63.html#figgel3">3.33</A>.<P>The calling program is the following:<P><P><P><PRE>C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++C EXAMPLE 3 : MODULE GEL3D1C --------- ( PARALLELEPEDICAL DOMAIN ==> PENTAHEDRA ) C ( UTILISATION OF THE 3 FUNCTIONS F1, F2 AND F3 )C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ PARAMETER ( LM = 30000 ) COMMON M(LM) COMMON /CENTR/XC,YC,ANG REAL X(8,5),Y(8,5),H(1) INTEGER NUMAT(14,4),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,'GEL3.NOPO','UNFORMATTED',0,IOSTAT) C ------ THE GEOMETRY : XC = -7. YC = 0. XPI = 3.141592654 ANG = XPI / 20. C N1 = 8 N2 = 5 XH1 = 2. XH2 = 1. DO 1 I=1,8 X(I,1) = ( I - 1 ) * XH1 X(I,5) = ( I - 1 ) * XH2 + 3.5 X(I,2) = .25 * ( 3.* X(I,1) + X(I,5) ) X(I,3) = .50 * ( X(I,1) + X(I,5) ) X(I,4) = .25 * ( X(I,1) + 3.*X(I,5) ) Y(I,1) = 0. Y(I,2) = 1.5 Y(I,3) = 3. Y(I,4) = 4.5 Y(I,5) = 6. 1 CONTINUE NB = 20 IOPTM = 7 NOPTDE = 1 C ------ THE NUMBERS : JOPTSD = 3C ------ THE SUB-DOMAINS : DO 2 I=1,14 DO 2 J=1,4 NUMAT(I,J) = 1 2 CONTINUE NUMAT(5,2) = 2 NUMAT(5,3) = 2 NUMAT(6,3) = 2 NUMAT(7,3) = 2C 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 C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++C THE THREE FUNCTIONS F1, F2 AND F3C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ FUNCTION F1(I,J,X,Y,K) COMMON /CENTR/ XC,YC,ANG F1 = XC + ( X - XC ) * COS( ( K - 1 ) * ANG ) END FUNCTION F2(I,J,X,Y,K) COMMON /CENTR/ XC,YC,ANG F2 = Y END FUNCTION F3(I,J,X,Y,K) COMMON /CENTR/ XC,YC,ANG F3 = ( X - XC ) * SIN( ( K - 1 ) * ANG ) END</PRE><P><P><P>The listing of the output generated 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='gel3.nopo',SPEC='UNFORMATTED',RECL=0) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE GEL3D1 : &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& GENERATION OF 1120 PENTAHEDRA AND 840 NODES LENGTH OF M (OVER-ESTIMATION) TO GENERATE NOPO : 29758 NUMBER OF POINTS IN X (NX) : 8 NUMBER OF POINTS IN Y (NY) : 5 NUMBER OF BANDS (NB) : 20 GENERATION OPTION (IOPT) : 7 NUMBERING OPTION OF THE D.S. (JOPT) : 3 SPLITTING OPTION (NOPT) : 1 INITIAL COORDINATES GIVEN : X Y DE I,J 1 1 : 0.0000E+00 0.0000E+00 2 1 : 0.2000E+01 0.0000E+00 3 1 : 0.4000E+01 0.0000E+00 4 1 : 0.6000E+01 0.0000E+00 5 1 : 0.8000E+01 0.0000E+00 6 1 : 0.1000E+02 0.0000E+00 7 1 : 0.1200E+02 0.0000E+00 8 1 : 0.1400E+02 0.0000E+00 1 2 : 0.8750E+00 0.1500E+01 2 2 : 0.2625E+01 0.1500E+01 3 2 : 0.4375E+01 0.1500E+01 4 2 : 0.6125E+01 0.1500E+01 5 2 : 0.7875E+01 0.1500E+01 6 2 : 0.9625E+01 0.1500E+01 7 2 : 0.1138E+02 0.1500E+01 8 2 : 0.1313E+02 0.1500E+01 1 3 : 0.1750E+01 0.3000E+01 2 3 : 0.3250E+01 0.3000E+01 3 3 : 0.4750E+01 0.3000E+01 4 3 : 0.6250E+01 0.3000E+01 5 3 : 0.7750E+01 0.3000E+01 6 3 : 0.9250E+01 0.3000E+01 7 3 : 0.1075E+02 0.3000E+01 8 3 : 0.1225E+02 0.3000E+01 1 4 : 0.2625E+01 0.4500E+01 2 4 : 0.3875E+01 0.4500E+01 3 4 : 0.5125E+01 0.4500E+01 4 4 : 0.6375E+01 0.4500E+01 5 4 : 0.7625E+01 0.4500E+01 6 4 : 0.8875E+01 0.4500E+01 7 4 : 0.1013E+02 0.4500E+01 8 4 : 0.1138E+02 0.4500E+01 1 5 : 0.3500E+01 0.6000E+01 2 5 : 0.4500E+01 0.6000E+01 3 5 : 0.5500E+01 0.6000E+01 4 5 : 0.6500E+01 0.6000E+01 5 5 : 0.7500E+01 0.6000E+01 6 5 : 0.8500E+01 0.6000E+01 7 5 : 0.9500E+01 0.6000E+01 8 5 : 0.1050E+02 0.6000E+01 ARRAY NUMAT OF TYPE 1 : 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 2 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 LENGTH CALCULATED OF M FOR D.S. NOPO 28318 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 ) : 2 NODES AND POINTS COINCIDE (NCOPNP) : 1 NUMBER OF ELEMENTS IN THE MESH (NE ) : 1120 NUMBER OF PENTAHEDRA (NPENT) : 1120 NUMBER OF BOUNDARY ELEMENTS (NEF ) : 760 NUMBER OF NODES (NOE ) : 840 NUMBER OF NODES PER SEGMENT (NO EXTREMITIES) : 0 TYPE OF COORDINATE VALUES (NTYCOO) : REEL1MOT MAX DIFFERENCE +1 BETWEEN 2 NODES OF AN ELEMENT : 49 NUMBER OF COARSE ELEMENTS (NBEGM ) : 0 NUMBER OF WORDS IN ARRAY NOP5 (LNOP5 ) : 25720 REFERENCE AXIS X,Y,Z (NTACOO) : 1 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE GEL3D1 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&</PRE><P><P><P><P><A NAME=5215> </A><IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img184.gif"><BR><STRONG>Figure 3.33:</STRONG> <i> Third example of a mesh generated by <b> GEL3D1</b></i><A NAME=figgel3> </A><BR><P><HR SIZE=3 WIDTH="75%"><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><BR> <HR><P><ADDRESS></ADDRESS></BODY></HTML>⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -