📄 node18.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>2.2 Test 1: Thermal Problem</TITLE><LINK REL=STYLESHEET TYPE="text/css" HREF="./Modulef.css" TITLE="Modulef CSS"><meta name="description" value="2.2 Test 1: Thermal Problem"><meta name="keywords" value="Guide1"><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=tex2html414 HREF="node17.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/previous_motif.gif" ALT="previous"></A><A NAME=tex2html420 HREF="node16.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/up_motif.gif" ALT="up"></A><A NAME=tex2html422 HREF="node19.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/next_motif.gif" ALT="next"></A><A NAME=tex2html424 HREF="node2.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/contents_motif.gif" ALT="contents"></A><A HREF="../Guide1-18/node18.html"><IMG BORDER=0 SRC="../icons/zoom18.gif" ALIGN=BOTTOM ALT="[BIG]"></A><A HREF="../Guide1-14/node18.html"><IMG BORDER=0 SRC="../icons/zoom14.gif" ALIGN=BOTTOM ALT="[Normal]"></A><A HREF="../Guide1-10/node18.html"><IMG BORDER=0 SRC="../icons/zoom10.gif" ALIGN=BOTTOM ALT="[small]"></A><BR><B> Next: </B> <A NAME=tex2html423 HREF="node19.html">2.3 Test 2: 2D Elasticity Problem </A><B>Up: </B> <A NAME=tex2html421 HREF="node16.html">2 Conversational tests</A><B> Prev: </B> <A NAME=tex2html415 HREF="node17.html">2.1 Introductory Example</A><B><A HREF="node2.html" >Contents</A></B><HR SIZE=3 WIDTH="75%"><H1><A NAME=SECTION04220000000000000000>2.2 Test 1: Thermal Problem</A></H1><P><P><P>This test has been proposed by NAFEMS (National Agency for Finite Element Methods & Standards) [<A HREF="node65.html#nafems1">Cameron et al.</A>].<P><P><P><H2><A NAME=SECTION04221000000000000000>2.2.1 Description</A></H2><P><P><P>Consider the thermal problem described below:<P><P><A NAME=2042> </A><IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img165.gif"><BR><STRONG>Figure 2.5:</STRONG> <i> Two-dimensional thermal test with heat transfer</i><A NAME=fig1> </A><BR><P><P>The internal heat production is zero.<P><TABLE COLS=2><COL ALIGN=RIGHT><COL ALIGN=LEFT><TR><TD VALIGN=BASELINE ALIGN=RIGHT NOWRAP>Boundary conditions:</TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD></TR><TR><TD VALIGN=BASELINE ALIGN=RIGHT NOWRAP> on boundary <b>AB</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> T = 100癈 . </TD></TR><TR><TD VALIGN=BASELINE ALIGN=RIGHT NOWRAP> on boundary <b>DA</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Zero flux. </TD></TR><TR><TD VALIGN=BASELINE ALIGN=RIGHT NOWRAP> on boundaries <b>BC</b>, <b>CD</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Heat transfer. </TD></TR><TR><TD VALIGN=BASELINE ALIGN=CENTER NOWRAP COLSPAN=2> Surrounding temperature 0癈</TD></TR></TABLE><P>Material properties:<P>Isotropic material with thermal conductivity:<DIV ALIGN=center><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img166.gif"></DIV> Coefficient of heat transfer on boundaries <b>BC</b>, <b>CD</b>:<DIV ALIGN=center><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img167.gif"></DIV><P><P><P><H2><A NAME=SECTION04222000000000000000>2.2.2 The boundary value problem</A></H2><P><P><P><DIV ALIGN=center><IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img168.gif"></DIV><P>The type of finite element utilized in order to solve this problem: <b> TRIA 2P2D</b>. <BR><P><P><A NAME=2051> </A><IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img169.gif"><BR><STRONG>Figure 2.6:</STRONG> <i> Suggested mesh</i><A NAME=fig2> </A><BR><P><P><b> Result:</b> The known temperature at point <b>E</b>, or at vertex 23 with the coordinates (0.6 , 0.2) is:<DIV ALIGN=center><IMG BORDER=0 ALIGN=BOTTOM ALT="" SRC="img170.gif"></DIV><P><P><P><H2><A NAME=SECTION04223000000000000000>2.2.3 Execution of the test</A></H2><P><P><P>As shown in the previous section, we execute the test by following various steps, indicated below. The data files have already been created for your convenience, i.e., all that remains to be done is to call the preprocessor and tochoose the EXECUTE option from the main menu.<P><P><P><b> Step 1:</b> Mesh<P><TABLE COLS=3><COL ALIGN=LEFT><COL ALIGN=LEFT><COL ALIGN=LEFT><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP>Call preprocessor: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> APNOXX</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Libraries used: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP COLSPAN=2> <b> CONV</b>, <b> NOP2</b>, <b> NOPO</b>, <b> UTII</b>, <b> UTSD</b></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Modules used: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> APNOPO</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> ... 2D mesh generator [<A HREF="node65.html#mod_104"><A NAME=tex2html37 HREF="../Guide3/welcome.html">MODULEF User Guide - 3</A></A>] </TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> QUACOO</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> ... generation of quadrangles </TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> MODNOP</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> ... refinement transformation </TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> ADPNOP</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> ... addition of nodes </TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> GIBBS</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> ... renumbering </TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Specify input data file: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> D1NOPO.D</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Output data structure: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> T1NOPO</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR></TABLE><P><P><P><b> Step 2:</b> Interpolation [<A HREF="node65.html#mod_13">13</A>]<P><TABLE COLS=3><COL ALIGN=LEFT><COL ALIGN=LEFT><COL ALIGN=LEFT><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP>Call preprocessor: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> COMAXX</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Libraries used: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP COLSPAN=2> <b> CONV</b>, <b> COSD</b>, <b> THER</b>, <b> UTII</b></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Module used: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> COMACO</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Element type: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> TRIA 2P2D</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Specify input data file: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> D1COMA.D</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Output data structures: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> T1MAIL</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> T1COOR</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR></TABLE><P><P><P><b> Step 3:</b> Calculation of the element arrays (constant case) [<A HREF="node65.html#mod_14">14</A>]<P><TABLE COLS=3><COL ALIGN=LEFT><COL ALIGN=LEFT><COL ALIGN=LEFT><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP>Call preprocessor: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> THECXX</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Libraries used: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP COLSPAN=2><b> CONV</b>, <b> COSD</b>, <b> THER</b>, <b> UTII</b></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Module used: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> THERCT</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Specify input data file: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> D1TAE.D</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Output data structure: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> T1TAE</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR></TABLE><P><P><P><b> Step 4:</b> Description of boundary conditions [<A HREF="node65.html#mod_18">18</A>]<P><TABLE COLS=3><COL ALIGN=LEFT><COL ALIGN=LEFT><COL ALIGN=LEFT><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP>Call preprocessor: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> COBDXX</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Libraries used: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP COLSPAN=2> <b> CONV</b>, <b> COSD</b>, <b> UTII</b></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Modules used: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> COBDC1</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> COTYNO</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> CALMAI</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Specify input data file: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> D1BDCL.D</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Output data structure: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> T1BDCL</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR></TABLE><P><P><P><b> Step 5:</b> Assembly and solution by a direct method (Cholesky) <BR>[4][<A HREF="node65.html#mod_5"><A NAME=tex2html39 HREF="../Guide5/welcome.html">MODULEF User Guide - 5</A></A>]<P><TABLE COLS=3><COL ALIGN=LEFT><COL ALIGN=LEFT><COL ALIGN=LEFT><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP>Call preprocessor: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> CHOLXX</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Libraries used: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP COLSPAN=2> <b> CONV</b>, <b> COSD</b>, <b> RESO</b>, <b> RESB</b>, <b> UTII</b>, <b> UTSD</b></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Modules used: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> PREPAC</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> ... allocation of pointers </TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> ASSMUA</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> ... assembly of matrix in main memory (m.m.) </TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> ASEMBV</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> ... assembly of RHS vector in main memory </TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> CLIMPC</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> ... factorization by Cholesky </TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> DRCHPC</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> ... forward- and backward substitution </TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Specify input data file: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> D1B.D</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR><TR><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> Output data structure: </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP> <b> T1B</b> </TD><TD VALIGN=BASELINE ALIGN=LEFT NOWRAP></TD></TR></TABLE><P>The number of words necessary to store the matrix : 922<P>The number of words necessary to store the corresponding right-hand-side (RHS) vector : 217<P><P><P><H2><A NAME=SECTION04224000000000000000>2.2.4 Memory usage and execution time</A></H2>
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -