📄 node17.html
字号:
<DD> = <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img145.gif">, a <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img146.gif"> unit matrix,<DT><b>g</b><DD> = 0 ,<DT><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img147.gif"><DD> = <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img148.gif">, and<DT><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img149.gif"><DD> = 0,<P> </DL>we obtain the simplified form of the mass and stiffness matrices and force vector, which enables us to solve the present problem as a thermal type problem (see [<A HREF="node65.html#mod_100">100</A>]).<P>The data description is performed with preprocessor <b> FOMIXX</b>. We choose to define the function<IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img147.gif"> via the subroutine <b> FORCEV</b> and the rest of the data <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img150.gif"> in array <b> /CAR</b>.Subroutine <b> FORCEV</b> is printed below to aid you.<P>The input is described in detail below:<P><UL><LI> The name of the data file must be entered (ex: <b> FOMI.DATA</b>).<P><LI> If more information is required regarding the input data, enter " Y " to the prompt " HELP ".<P><LI> The file names of the files containing the following data structures must be entered, along with their level numbers and number of associated arrays:<P> <UL><LI> O.D.S. <b> FORC</b> describing the loads (ex: <b> TEST.FORC</b> 1 0), and <LI> O.D.S. <b> MILI</b> describing the characteristics of materials (ex: <b> TEST.MILI</b> 1 1). </UL><P><LI> Firstly, we will describe the manner of furnishing the data necessary to calculate the RHS vector(s) (for the computation of <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img147.gif">). <BR><P> The number of load cases, or RHS vectors, is given by:<P> <DL COMPACT><DT>NDSM:<DD> (ex: 1)<P> </DL><P><LI> The processing is done globally (by reference numbers or sub-domain numbers) (ex: Y).<P><LI> The following values must then be input:<P> <DL COMPACT><DT>NDIM:<DD> dimension (ex: 2), <DT>NDSD:<DD> total number of sub-domains (ex: 1), and <DT>NBFR:<DD> total number of reference numbers (ex: 0).<P> </DL><P><LI> There is surface data applicable to the present sub-domain (ex: Y), and nothing more (ex: N). <BR><P><LI> A description of the surfaces loads must be input:<P> <DL COMPACT><DT>NOPT:<DD> (ex: 1) <BR> <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img151.gif"> for a isotropic material, and <BR> <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img152.gif"> otherwise. <DT>ITRAIT:<DD> (ex: 2) <BR> The data is furnished by: <BR> 1: array, <BR> 2: subroutine <b> FORCE</b>, or <BR> 3: user defined subroutine.<P> </DL><P><LI> Information concerning subroutine <b> FORCE</b> must be entered (ITRAIT=2):<P> <DL COMPACT><DT>LOPT:<DD> subroutine option for this load case (ex: 1), and <DT>IADR:<DD> start address in the subroutine (ex: 1).<P> </DL><P><LI> We will now proceed to describe the data corresponding to the matrices (ex: <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img143.gif"> and <b>[k]</b> corresponding to the mass and stiffness matrices respectively).<P><LI> For each array associated (see above) with the D.S. <b> MILI</b> (ex: 1), we must specify:<P> <DL COMPACT><DT>NOMTAB:<DD> the name of the array (ex: /CAR), <DT>NTYTAB:<DD> its type (ex: 2), <BR> 1: integer, <BR> 2: single precision, <BR> 5: double precision, <BR> 3: logical, and <BR> 4: character, <BR> generally, for a thermal problem : single precision <IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img153.gif"> type 2, and <BR> in elasticity : double precision <IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img153.gif"> type 5 and <DT>NBVAR:<DD> the number of variables (ex: 4).<P> </DL><P><LI> The data to be stored in the array specified above must then be entered:<P> In the present example, to calculate (recall equations (<A HREF="node15.html#eqclass">1.7</A>) and (<A HREF="node15.html#eqsimp">1.8</A>)):<P> <DL COMPACT><DT><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img154.gif">:<DD> set <b>[k] = I</b> and <b>g = 0</b>, <DT><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img155.gif">:<DD> set <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img60.gif">, and <DT><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img156.gif">:<DD> set <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img157.gif"> and <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img158.gif">,<P> </DL><P> If the data values do not enter the problem under consideration, for example <b>g</b> and <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img149.gif"> are zero, we may concatenate the data.<P> Thus, for the present example, we input only:<P> <DL COMPACT><DT>/CAR(1):<DD> <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img143.gif"> (ex: 1), <DT>/CAR(2):<DD> <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img159.gif"> (ex: 1), <DT>/CAR(3):<DD> <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img160.gif"> (ex: 0), and <DT>/CAR(4):<DD> <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img161.gif"> (ex: 1).<P> </DL><P> A comment, written between single quotes, must be input to describe the above array of data <DIV ALIGN=center>(ex: ' RO AND THERMAL CONDUCTIVITY (UNIT MATRIX) '). </DIV><P> You also have the opportunity to re-enter the data in the array under consideration.<P><LI> As before, the processing is done globally (of elements) (by reference numbers or sub-domain numbers) (ex: Y),<P><LI> after which the following values must then be input:<P> <DL COMPACT><DT>NDIM:<DD> dimension (ex: 2), <DT>NDSD:<DD> total number of sub-domains (ex: 1), and <DT>NBFR:<DD> total number of reference numbers (ex: 0).<P> </DL><P><LI> The description of the surfaces is now given by:<P> <DL COMPACT><DT>NOPT:<DD> (ex: 1, as before) <BR> <b> NOPT</b> <b>></b> 0 for a isotropic material, and <BR> <b> NOPT</b> <b><</b> 0 otherwise. <DT>ITRAIT:<DD> (ex: 1, instead of 2 previously) <BR> The data is furnished by: <BR> 1: array, <BR> 2: subroutine <b> MILIEU</b>, or <BR> 3: user defined subroutine.<P> </DL><P><LI> Lastly, the following entities must be specified:<P> <DL COMPACT><DT>NTABL:<DD> the name of the data array for the present case (ex: /CAR), and <DT>IADR:<DD> the return address (ex: 1).<P> </DL></UL><P>After completing the data file you will be confronted with the main menu, from which you can choose to either create another data file (C), execute the module (E) once the data file is created, or terminate the process (END). Suppose that we choose to execute <b> FOMIXX</b> be typing a " E ", then theinput data are the following:<P><UL><LI> The name of the data file must be entered (ex: <b> FOMI.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 is printed below, followed by the output generated by executing the preprocessor:<P><P><P><PRE>B$ fomixx *************************** EXECUTING COFORC AND COMILI : ( DATA ) *************************** -- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ? - - ---C ************************************************** CONVERSATIONAL DATA CREATION FOR COFORC AND COMILI ************************************************** -- NAME OF FILE TO CREATE (A) ? FOMI.DATA -- HELP (YES-NO) ?N -- NOMFIC (FORC) NIFORC NTFORC ?TEST.FORC 1 0 -- NOMFIC (MILI) NIMILI NTMILI ?TEST.MILI 1 1 ** WE WILL DESCRIBE THE MANNER IN WHICH TO ENTER THE DATA NECESSARY TO CALCULATE THE RHS VECTOR(S) -- NDSM ?1 -- PROCESSING BY GROUP (SUB-DOMAIN OR REFERENCE NUMBERS) (YES-NO) ?Y -- NDIM NDSD NBFR ?2 1 0 SURFACE DATA S-D : 1 (YES-NO) ?Y OTHER INFORMATION (OYES-NO) ?N -- NOPT ITRAIT ( SURFACE REFERENCE : 1 ) ?1 2 -- LOPT IADR ?1 1 ** DATA DESCRIPTION FOR THE MATRI(X)(CES) NUMBER OF ASSOCIATED TABLES FOR THE D.S. MILI : 1 FOR EACH GIVEN TABLE ENTER: ITS NAME (A) NOMTAB ITS TYPE ( 1 TO 6 ) (I) NTYTAB ** INTEGER 1 , SP 2 , DP 5 , LOGICAL 3 , CHARACTER 4... ** HEAT PROBLEM : SINGLE PRECISION => TYPE 2 ** ELASTICITY PROBLEM : DOUBLE PRECISION => TYPE 5 ITS NUMBER OF VARIABLES (I) NBVAR ** TABLE 1 ASSOCIATED WITH THE D.S. MILI -- NOMTAB NTYTAB NBVAR ? /CAR 2 4 ** THE VARIABLES OF THE TABLE /CAR -- /CAR( 1) = ? 1 -- /CAR( 2) = ? 1 -- /CAR( 3) = ? 0 -- /CAR( 4) = ? 1 ** COMMENTS ABOUT THE TABLE /CAR (60 CHARACTERS MAXI BETWEEN ')'RO AND THE THERMAL CONDUCTIVITY' -- RE-ENTER THE CHARACTERISTICS OF THE TABLE :/CAR (YES-NO) ?N -- PROCESSING BY GROUP (SUB-DOMAIN OR REFERENCE NUMBERS) (YES-NO) ?Y -- NDIM NDSD NBFR ?2 1 0 SURFACE DATA S-D : 1 (YES-NO) ?Y OTHER INFORMATION (YES-NO) ?N -- NOPT ITRAIT ( SURFACE REFERENCE : 1 ) ?1 1 -- NTABL IADR ?/CAR 1 -- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ? - - ---E -- NAME OF DATA FILE (A) ?FOMI.DATA -- PRINT PARAMETER OF THE D.S. ?1 ++ OPEN(11,FILE='test.forc',SPEC='UNFORMATTED',RECL=0) ++ OPEN(12,FILE='test.mili',SPEC='UNFORMATTED',RECL=0) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE COFORC : &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE COFORC &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE COMILI : &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE COMILI &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& TABLE F O R 2 -------------- INFORMATION ( NUMBER OF SETS ) OR NOT ( 0 ) -------------------------------------------------- | | | | | | SUB-DOMAINS | ELEMENTS | | | | | -------------------------------------------------- | VOLUMES | 0 | 0 | | SURFACES | 1 | 0 | | LINES | 0 | 0 | | POINTS | 0 | 0 | --------------------------------------------------⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -