📄 node17.html
字号:
TABLE M I L 2 -------------- INFORMATION ( NUMBER OF SETS ) OR NOT ( 0 ) -------------------------------------------------- | | | | | | SUB-DOMAINS | ELEMENTS | | | | | -------------------------------------------------- | VOLUMES | 0 | 0 | | SURFACES | 1 | 0 | | LINES | 0 | 0 | | POINTS | 0 | 0 | -------------------------------------------------- -- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ? - - ---END</PRE><P><P><P>The resulting data file, generated by executing preprocessor <b> FOMIXX</b>, is shown below:<P><P><P><PRE> TEST.FORC $ FILE NAME OF THE D.S. FORC 1 0 $ ITS LEVEL AND NUMBER OF TABLES TEST.MILI $ FILE NAME OF THE D.S. MILI 1 1 $ ITS LEVEL AND NUMBER OF TABLES $ DATA RELATIVE TO THE D.S. FORC $ 1 $ NDSM 0 1 0 0 0 0 0 0 1 0 1 2 $ NOSD NFRO NOPT ITRAIT 1 1 $ LOPT IADR $ DATA RELATIVE TO THE D.S. MILI $ /CAR 2 4 $ NAME TYPE NBREMOT 0.100000E+01 $ /CAR( 1) 0.100000E+01 $ /CAR( 2) 0.000000E+00 $ /CAR( 3) 0.100000E+01 $ /CAR( 4) ' RO AND THE THERMAL CONDUCTIVITY ' $ CONTENTS 0 1 0 0 0 0 0 0 1 0 1 1 $ NOSD NFRO NOPT ITRAIT /CAR 1 $ NTABL IADR</PRE><P><P><P>The user defined subroutine <b> FORCE</b> corresponding to the present example is printed below:<P><P><P><PRE> SUBROUTINE FORCEV(M,LOPT,X,NDIM,NPO,TAR,LTAR,LVECT,IADR,I1, 1 NARE,IA)C.....................................................................C AIM:C INPUT THE HEAT QUANTITIES IN ORDER TO CALCULATE THEC ELEMENT RHS VECTORS FOR A THERMAL PROBLEM.C.....................................................................C INPUT PARAMETERS:C M SUPER ARRAYC LOPT COFORC PARAMETERC X ELEMENT COORDINATE TABLE C NDIM DIMENSIONC NPO NUMBER OF POINTS FOR THE ELEMENTC LTAR NUMBER OF WORDS OF THE ARRAY TAR FOR EVRY ELEMENTC LVECT NUMBER OF SAME COLOR ELEMENTSC IADR COFORC PARAMETERC I1 0 USUALLY (SEE DOC 14)C NARE FACE NUMBER (3D), EDGE NUMBER (2D)C IA ADRESS -1 IN MAIL OF THE 1ST WORD OF THE CURRENT ELEMENTCC OUTPUT PARAMETERS:C TAR ARRAY TAR(LTAR,LVECT)C......................................................................C INTEGER M(*),IP(3) REAL X(NPO,NDIM),TAR(*) DATA IP/2,3,1/ IF (NPO .EQ. 3) THENCC THE ELEMENT IS STRAIGHTC DO 1 I = 1, 3 XM = 0.5*(X(I,1)+X(IP(I),1)) YM = 0.5*(X(I,2)+X(IP(I),2)) TAR(I) = 6*XM*YM*(3-2*(XM**2)-2*(YM**2)+(XM**2)*(YM**2)) 1 CONTINUE ELSE IF (NPO.EQ. 6) THENCC THE ELEMENT IS CURVEDC DO 2 I = 1, 6 TAR(I) = 6*X(I,1)*X(I,2)*(3-2*X(I,1)**2-2*X(I,2)**2 + +X(I,1)**2*X(I,2)**2) 2 CONTINUE ELSECC ERROR MESSAGEC IMPRIM = IINFO('I') WRITE(IMPRIM,1000) 'ERROR FORCE' STOP END IF 1000 FORMAT(1X,A) END</PRE><P><P><P><b> Comments:</b> The two types of finite elements do not need the same amount of data. It is possible to distinguish between a straight element and a curved element by the number of points (3 for a straight element, 6 for a curved element).<P>The element arrays are created by preprocessor <b> THELXX</b>. The input data file is generated by choosingoption CREATE from the main menu and supplying the relevant responses to the remaining options orquestions.<P>The input is described in detail below:<P><UL><LI> The name of the data file must be entered (ex: <b> TAE.DATA</b>).<P><LI> The file names of the files containing the following data structures must be entered:<P> <UL><LI> I.D.S. <b> MAIL</b> containing the mesh and interpolation data (ex: <BR>[4] <b> TEST.MAIL</b>), <LI> I.D.S. <b> COOR</b> containing the coordinates of points (ex: <b> TEST.COOR</b>), and <LI> O.D.S. <b> TAE</b> containing the element arrays: mass, stiffness, loads, and stresses (ex: <b> TEST.TAE</b>). </UL><P><LI> There are no associated data arrays associated (ex: N).<P><LI> We do not make use of the <b> POBA</b> file (ex: N).<P><LI> We utilize the I.D.S. <b> MILI</b> (ex: Y) and input:<P> <UL><LI> the name of the I.D.S. <b> MILI</b> file containing the material characteristics (ex: <b> TEST.MILI</b>). </UL><P><LI> We utilize the I.D.S. <b> FORC</b> (ex: Y) and input:<P> <UL><LI> the name of the file containing the I.D.S. <b> FORC</b> describing the RHS loads (ex: <b> TEST.FORC</b>). </UL><P><LI> Specify the type of the problem under consideration:<P> <DL COMPACT><DT>NPROV:<DD> (ex: 1) <BR> 1: thermal, <BR> 2: elastic, <BR> 3: streamlines.<BR><P> </DL><P><LI> The mass matrix (not diagonal), stiffness matrix and RHS vector must be calculated, the flux is not calculated for this example (ex: Y, N, Y, Y, and N).<P> There is a linear combination between mass/stiffness (ex: Y and 1). The coefficients are stored in a array, /LIN in this example, and contains the values (6, 1) (keeping in mind the variational formulation of our problem).<P></UL><P>Having created the data file, we proceed by choosing option EXECUTE from themain menu. However, before being able to execute <b> THELXX</b>, subroutine <b> FORCE</b> must be compiled and linked to preprocessor <b> THELXX</b>.<P>The input is described below:<P><UL><LI> The name of the data file must be entered (ex: <b> TAE.DATA</b>).<P><LI> The output parameter ( 0 <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img137.gif"> <b> IMPRE</b> <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img137.gif"> 10) must be specified (ex: 5).<P></UL><P>The data creation procedure and screen output of preprocessor <b> THELXX</b> is given below:<P><P><P><PRE>B$ thelxx **************** EXECUTING THELAS (NEW) : ( ELEMENT COMPUTATIONS ) **************** -- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ? - - ---C************************************************************************** CONVERSATIONAL DATA CREATION OF THENEW************************************************************************** -- NAME OF FILE TO CREATE ?TAE.DATA -- FILE NAME OF THE INPUT D.S. MAIL ?TEST.MAIL -- FILE NAME OF THE INPUT D.S. COOR ?TEST.COOR -- FILE NAME OF OUTPUT D.S. TAE ?TEST.TAE -- ARE THERE ASSOCIATED TABLES (YES-NO) ?N -- DO YOU USE THE POBA FILE (YES-NO) ? ( BASIS POLYNOMIALS ... IN DIRECT ACCESS )N -- DO YOU USE THE D.S. MILI (YES-NO) ?Y -- FILE NAME OF THE INPUT D.S. MILI ?TEST.MILI -- DO YOU USE THE D.S. FORC (YES-NO) ?Y -- FILE NAME OF THE INPUT D.S. FORC ?TEST.FORC -- TYPE OF PROBLEM (NPROV) ? 1 : THERMAL , 2 : ELASTIC , 3 : STREAMLINES1 -- COMPUTE THE MASS MATRIX (YES-NO) ?Y -- IT MAY BE DIAGONAL (YES-NO) ?N -- COMPUTE THE STIFFNESS MATRIX (YES-NO) ?Y -- COMPUTE THE RHS VECTORS (YES-NO) ?Y -- COMPUTE THE FLUX (YES-NO) ?N -- ARE THERE LINEAR MASS-STIFFNESS COMBINATIONS (YES-NO) ?Y -- NUMBER OF LINEAR MASS-STIFFNESS COMBINATIONS ?1 -- NAME (4 CHARACTERS ) OF THE TABLE WITH THESE VALUES ?/LIN -- COEFFICIENTS OF COMBINATION NUMBER : 1 ?6 1 -- FILE NAME OF COLOURS (NOMCOU) OR EMPTY AND MAX LENGTH OF THE VECTORS (LVECT)VIDE 1 -- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ? - - ---E !! YOU USE THE LAST VERSION OF THELAS YOUR DATA FILE MUST END WITH : -- NAME OF THE DATA FILE ? TAE.DATA -- PRINT PARAMETER ?1 ++ OPEN(11,FILE='test.mail',SPEC='OLD,UNFORMATTED',RECL=0) ++ OPEN(12,FILE='test.coor',SPEC='OLD,UNFORMATTED',RECL=0) ++ OPEN(13,FILE='test.tae',SPEC='UNFORMATTED',RECL=0) ++ OPEN(14,FILE='test.mili',SPEC='OLD,UNFORMATTED',RECL=0) ++ OPEN(15,FILE='test.forc',SPEC='OLD,UNFORMATTED',RECL=0) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE THENEW : ** WARNING READRE : ARRAY TRAV HAS ZERO LENGTH WITH ADRES SET TO 1 TABLE T A E 2 -------------- NUMBER OF ELEMENTS (NE) : 72 NUMBER OF NODES (NOE) : 171 NUMBER OF TABLES ASSOCIATED AT EACH ELEMENT (NTACE) : 2 MAX NUMBER OF NODES BY ELEMENT (NNOMAX) : 6 CONSTANT NUMBER OF D.O.F. BY NODES ( 0 IF NOT ) (ND) : 1 MAXIMUM NUMBER OF D.O.F. PER NODE (NDLMAX) : 1 NUMBER OF TYPES OF ELEMENTS (NTYELM) : 2 SUM OF THE NUMBER OF NODES OF THESE TYPES (LNOET) : 12 PROBLEM : 1 THERMAL, 2 ELASTIC, 3 OTHER (NPROV) : 1 REQUESTED OPTIONS (NOPTNT) : 0 NUMBER OF THE ASSOCIATED TABLE COLLECTING : THE MASS MATRIX : 0 THE STIFFNESS MATRIX : 0 THE FIRST LINEAR COMBINATION : 1 THE LAST ONE : 1 THE RIGHT HAND SIDE(S) : 2 THE STRESSES OR THE FLOW : 0 END OF MODULE THENEW &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ? - - ---END</PRE><P><P><P>The resulting data file is given below:<P><P><P><PRE> TEST.MAIL $ FILE NAME 1 $ AND LEVEL OF THE D.S. MAIL TEST⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -