📄 node58.html
字号:
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) ENDC +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ C SUBROUTINE XYZ33C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ SUBROUTINE XYZ33(I,NUM,NREF,X,Y,Z,X3,Y3,Z3)C X3 = 0.10 * I * X Y3 = 0.10 * I * Y Z3 = Z + 0.450 ENDC +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ C SUBROUTINE XYZ23C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ SUBROUTINE XYZ23(I,NUM,NREF,X,Y,Z,X3,Y3,Z3)C X3 = X / 2.5 Y3 = Y / 2.5 Z3 = 10. END</PRE><P><P><A NAME=4974> </A><IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img154.gif"><BR><STRONG>Figure 3.15:</STRONG> <i> Example <b> MA2D3E</b> (an interpolation followed by step-by-step definitions)</i><A NAME=figma2336> </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 TO 5 TYPE OF TRANSFORMATION -6 SECTION 5 TO 9 TYPE OF TRANSFORMATION -3 NUMBER OF TRANSFERS OF NUMBERS : 2 LIST OF TRANSFERS : SECTION 0 TO 0 TYPE OF TRANSFER 0 : -100 GIVES 1 SECTION 9 TO 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 ) : 1980 NUMBER OF PENTAHEDRA (NPENT) : 1980 NUMBER OF BOUNDARY ELEMENTS (NEF ) : 72 NUMBER OF NODES (NOE ) : 1420 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 ) : 29040 REFERENCE AXIS X,Y,Z (NTACOO) : 1 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MA2D3E &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&</PRE><P>An example with <b> TTYPE = -4</b>, i.e., starting from the 2D mesh data (a ring with reference 1 forthe items of the internal circle and reference 2 for the external circle), we define the basis via ZINT(0) = 10. and construct sections 1 to 4 of the sides given in array ZINT(1:4) by translation.<P>The references of the 3D mesh are attributed via array DESREF(.,.) (see the calling program for thespecification of the transfers and figure <A HREF="node58.html#figma234">3.16</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 4C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ PARAMETER (LM = 100000 , MAXFON = 1 , 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 CALL TRUNIT (NFNOPO) NINOPO = 1 CALL OUVRIR(NFNOPO,'NOPO2D.MA23.1','OLD,UNFORMATTED',0,IOSTAT) NFNOPS = 11 CALL TRUNIT (NFNOPS) NINOPS = 2 CALL OUVRIR(NFNOPS,'NOPO3D.MA23.4','UNFORMATTED',0,IOSTAT) C -----------------------------------------------C ------ GEOMETRICAL TRANSFORMATIONS ------ C ----------------------------------------------- C --- The BASIS : BFONC = .FALSE. ZINT(0) = 10. C --- THE INTERMEDIARY SIDES : NBCOU = 4 ZINT(1) = 11. ZINT(2) = 12. ZINT(3) = 13. ZINT(4) = 14.5 C --- A FUNCTION : A MANUAL TRANSLATION NBFONC = 1 C IFONC = 1 NUCOU(1,IFONC) = 0 NUCOU(2,IFONC) = 4 NUCOU(3,IFONC) = -4C -----------------------------------------------C ------ DEFINITION OF THE REFERENCES ------ C ----------------------------------------------- NBDES = 6C --- BASIS : FACES EDGES AND POINTS SET TO REFERENCE 3 : IDES = 1 DESREF(1,IDES) = 0 DESREF(2,IDES) = 0 DESREF(3,IDES) = 0 DESREF(4,IDES) = -100 DESREF(5,IDES) = 3 C ---C --- ABOVE : FACES EDGES AND POINTS SET TO REFERENCE 4 : C ABOVE : THEN EDGE 1 (2D) ==> EDGE 7 (3D)C ABOVE : THEN EDGE 2 (2D) ==> EDGE 8 (3D) IDES = 2 DESREF(1,IDES) = 4 DESREF(2,IDES) = 4 DESREF(3,IDES) = 0 DESREF(4,IDES) = -100 DESREF(5,IDES) = 4 IDES = 3 DESREF(1,IDES) = 4 DESREF(2,IDES) = 4 DESREF(3,IDES) = 4 DESREF(4,IDES) = 1 DESREF(5,IDES) = 7 IDES = 4 DESREF(1,IDES) = 4 DESREF(2,IDES) = 4 DESREF(3,IDES) = 4 DESREF(4,IDES) = 2 DESREF(5,IDES) = 8 C --- VERTICAL FACES : EDGE 1 (2D) ==> FACE 5 (3D)C EDGE 2 (2D) ==> FACE 6 (3D) IDES = 5 DESREF(1,IDES) = 0 DESREF(2,IDES) = 4 DESREF(3,IDES) = 3 DESREF(4,IDES) = 1 DESREF(5,IDES) = 5 IDES = 6 DESREF(1,IDES) = 0 DESREF(2,IDES) = 4 DESREF(3,IDES) = 3 DESREF(4,IDES) = 2 DESREF(5,IDES) = 6 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=4992> </A><IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img155.gif"><BR><STRONG>Figure 3.16:</STRONG> <i> Example <b> MA2D3E</b> (translation)</i><A NAME=figma234> </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.1',SPEC='OLD,UNFORMATTED',RECL=0) ++ OPEN(11,FILE='nopo3d.ma23.4',SPEC='UNFORMATTED',RECL=0) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE MA2D3E : &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& HEIGHT OF THE BOTTOM BASIS : 0.1000000E+02 NUMBER OF GEOMETRIC TRANSFORMATIONS : 1 LIST OF TRANSFORMATIONS : SECTION 0 TO 4 TYPE OF TRANSFORMATION -4 NUMBER OF NUMBER TRANSFERS : 6 LIST OF TRANSFERS : SECTION 0 TO 0 TYPE OF TRANSFER 0 : -100 GIVES 3 SECTION 4 TO 4 TYPE OF TRANSFER 0 : -100 GIVES 4 SECTION 4 TO 4 TYPE OF TRANSFER 4 : 1 GIVES 7 SECTION 4 TO 4 TYPE OF TRANSFER 4 : 2 GIVES 8 SECTION 0 TO 4 TYPE OF TRANSFER 3 : 1 GIVES 5 SECTION 0 TO 4 TYPE OF TRANSFER 3 : 2 GIVES 6 MESH OBTAINED : SAVE D.S. TABLE N O P 2 -------------- MESH CHARACTERISTICS : SPACE DIMENSION (NDIM ) : 3 MAXIMUM NUMBER OF REFERENCES (NDSR ) : 8 MAXIMUM NUMBER OF SUB-DOMAINS (NDSD ) : 1 NODES AND POINTS COINCIDE (NCOPNP) : 1 NUMBER OF ELEMENTS IN THE MESH (NE ) : 144 NUMBER OF HEXAHEDRA (NHEX ) : 144 NUMBER OF BOUNDARY ELEMENTS (NEF ) : 28 NUMBER OF NODES (NOE ) : 240 NUMBER OF NODES PER SEGMENT (NO EXTREMITIES) : 0 TYPE OF COORDINATE VALUES (NTYCOO) : REEL1MOT MAX DIFFERENCE +1 BETWEEN 2 NODES OF AN ELEMENT : 78 NUMBER OF COARSE ELEMENTS (NBEGM ) : 0 NUMBER OF WORDS IN ARRAY NOP5 (LNOP5 ) : 4968 REFERENCE AXIS X,Y,Z (NTACOO) : 1 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MA2D3E &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&</PRE><P>An example with <b> TTYPE = -5</b>, i.e., definition of the upper basis (via XYZ23 at side Z=10.) and interpolation between the basis (ZINT(0) = -1.) and this section to construct 5 sections withequidistant sides between that of the basis and that of the upper section.<P>The references of the 3D mesh are attributed via via array DESREF(.,.) (see the calling program forthe specification of the transfers and figure <A HREF="node58.html#figma236">3.18</A> for the result obtained).<P>Extracts of the calling program and a plot of the resulting mesh is shown below.<P><PRE>C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++C EXAMPLE OF CALLING MODULE MA2D3E : TEST NUMBER 5C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ PARAMETER (LM = 100000 , MAXFON = 1 , 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.5','UNFORMATTED',0,IOSTAT) C -----------------------------------------------C ------ GEOMETRICAL TRANSFORMATIONS ------ C ----------------------------------------------- C --- LA BASE : BFONC = .FALSE. ZINT(0) = -1.C --- THE INTERMEDIARY SIDES : NBCOU = 5C --- A FUNCTION : A MANUAL INTERPOLATION C BETWEEN THE BASIS AND THE ONE ABOVE DEFINED BY XYZ23 C FOR SIDE 10. ( SEE XYZ23 ) NBFONC = 1 C IFONC = 1 NUCOU(1,IFONC) = 0 NUCOU(2,IFONC) = 5 NUCOU(3,IFONC) = -5C -----------------------------------------------C ------ DEFINITION OF THE REFERENCES ------ C ----------------------------------------------- NBDES = 5C --- 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 --- .... .... (SEE EXAMPLE BELOW) ....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)⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -