📄 node17.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.6 Computing the stresses (D.S. TAE)</TITLE><LINK REL=STYLESHEET TYPE="text/css" HREF="./Modulef.css" TITLE="Modulef CSS"><meta name="description" value="2.6 Computing the stresses (D.S. TAE)"><meta name="keywords" value="Guide6"><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=tex2html683 HREF="node16.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/previous_motif.gif" ALT="previous"></A><A NAME=tex2html689 HREF="node11.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/up_motif.gif" ALT="up"></A><A NAME=tex2html691 HREF="node18.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/next_motif.gif" ALT="next"></A><A NAME=tex2html693 HREF="node2.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/contents_motif.gif" ALT="contents"></A><A HREF="../Guide6-18/node17.html"><IMG BORDER=0 SRC="../icons/zoom18.gif" ALIGN=BOTTOM ALT="[BIG]"></A><A HREF="../Guide6-14/node17.html"><IMG BORDER=0 SRC="../icons/zoom14.gif" ALIGN=BOTTOM ALT="[Normal]"></A><A HREF="../Guide6-10/node17.html"><IMG BORDER=0 SRC="../icons/zoom10.gif" ALIGN=BOTTOM ALT="[small]"></A><BR><B> Next: </B> <A NAME=tex2html692 HREF="node18.html">2.7 Computing the flux and temperature </A><B>Up: </B> <A NAME=tex2html690 HREF="node11.html">2 Manipulation of solutions</A><B> Prev: </B> <A NAME=tex2html684 HREF="node16.html">2.5 Computation of interface and reaction </A><B><A HREF="node2.html" >Contents</A></B><HR SIZE=3 WIDTH="75%"><H1><A NAME=SECTION03360000000000000000>2.6 Computing the stresses (D.S. TAE)</A></H1><P><P><P>Preprocessor <b> STREXX</b> constructs a D.S. <b> TAE</b> containing the stresses bycombining the array containing the element stresses of D.S.<b> TAE</b>, set up at the same time as the element matrices and right-hand-sides, with D.S. <b> B</b> containing the displacements, constructed at the same time the solution to theelasticity problem was computed.<P>All the finite elements in elasticity are possible.<P>Module STRESS, called by <b> STREXX</b> calculates the stresses in each element at the barycentre or at thenumerical integration points (see technical files) and, in addition,can print the stresses of all or part of the elements in the mesh. For a two-dimensional problem(the visualization of the three-dimensional case has not been implemented), structure <b> TAE</b> resulting from STRESS is used to visualize the stresses(see <b> TRSTXX</b> in part II of this user guide).<P><P><P>Program STRESS has the following form:<P><PRE> SUBROUTINE STRESS(M,XM,ISTRES,LINDIC,INDIC,NFBE,NIBE,NFBT,NIBT, + NFTAE,NITAE,NFNDL1,NINDL1,NFNDL2,NINDL2, + NFMAIL,NIMAIL,NFTAES,NITAES,NTTAES)C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++C MODULE STRESS : COMPUTE THE STRESSES DURING THE SOLUTION OF ANC -------------- ELASTICITY PROBLEM FROM THE ELEMENT STRESSESC AND THE SOLUTIONC ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++</PRE><P>where<P><UL><LI> M and XM represent the work array,<LI> ISTRES is the stress print option. ISTRES is set to IJK where:<DL COMPACT><DT>I <b> = 1</b><DD>: print all the elements,<DT>I <b> = 2</b><DD>: print the elements lying between numbers N1 and N2 given in INDIC,<DT>I <b> = 3</b><DD>: print the LINDIC elements whose numbers are given in INDIC,<DT>I <b> = 4</b><DD>: print the elements whose sub-domain numbers (LINDIC numbers) are given in INDIC,<P> </DL><DL COMPACT><DT>J <b> = 1</b><DD>, if the principal stresses should be printed, 0 if not,<P> </DL><DL COMPACT><DT>K <b> = 1</b><DD>, if the principal directions should be printed (I forced to equal 1), 0 if not,<P> </DL><LI> LINDIC, a function of I, is the length of array INDIC,<LI> INDIC, a function of I, is the following array:<DL COMPACT><DT>I <b> = 1</b><DD>, the array is not used,<DT>I <b> = 2</b><DD>, length 2, the array contains the number of the first and last element to print,<DT>I <b> = 3</b><DD>, the number of the LINDIC elements to print,<DT>I <b> = 4</b><DD>, the number of the LINDIC sub-domains to treat. Stresses of all the elements withthese sub-domain numbers will be selected.<P> </DL><LI> NFBE, NIBE are the file number and the level of D.S. <b> B</b> containing the displacements (NFBE <b>></b> 0: thefile is sequential, NFBE = 0, the D.S. is in main memory, NFBE <b><</b> 0, the file is in direct access),<LI> NFBT, NIBT, the same for the D.S. <b> B</b> containing the temperatures (thermoelastic case),<LI> NFTAE, NITAE are the file number and the level of D.S. <b> TAE</b> containing the elementarrays (including the element stresses) NFTAE <b>></b> 0,<LI> NFNDL1, NINDL1, used if ND = 0, are the file number and the level of the elasticity D.S. <b> NDL1</b> (ND = 0 indicates that the number of degrees of freedom per node is not constant),<LI> NFNDL2, NINDL2, used if ND = 0, are the file number and the level of the thermal D.S. <b> NDL1</b>, <LI> NFMAIL, NIMAIL are the file number and the level of D.S. <b> MAIL</b>,<LI> NFTAES, NITAES, NTTAES are the file number, the level and the number of arrays to beassociated with the O.D.S <b> TAE</b> which contains the stresses element by element (NFTAES <b>></b> 0).</UL><P>The input files necessary are NFBE, NFTAE, NFMAIL and, in certain cases, NFBT, NFNDL1 and NFNDL2.The output file is NFTAES.<P>There is no input data required specifically for this module except in the case where NTTAES is non-zero. In this case, the following data is required:<P><P><P>Loop I from 1 to NTTAES<UL><LI> name of array I in 4 characters,<LI> array type (1 = integer, 2 = real single precision, 4 = character, 5 = real double precision, ...),<LI> number of words in this array,<LI> variables in this array,<LI> comments pertaining to the contents of array I in 10 words of 4 characters.</UL> End of loop.<P><P><P><P><P><P>Note that depending on ISTRES, the computation of eigenvalues and eigenvectors is either performed or not.Subroutine VPM2DD computes the eigenvalues and eigenvectors for the two-dimensional case.VAM3DD computes the eigenvalues whereas VPM3DD computes the eigenvectors forthe three-dimensional case.<P><P><P>The <b> TAE</b> structure resulting from module STRESS contains the informations corresponding tothe stresses of all the elements in the mesh. Module COMTAE is employedto select certain elements only and to construct the corresponding D.S. <b> TAE</b>.<P><PRE> SUBROUTINE COMTAE(M,ISTRES,INDIC,LINDIC,NFTAE,NITAE,NFTAES,NITAES)C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++C MODULE COMTAE : COMPRESS D.S. TAES ( RESULTING FROM STRESS )C ------------- ONLY THE ELEMENTS OF INDIC AND ISTRES ARE RETAINEDC +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++</PRE><P>where:<P><UL><LI> M represent the work array,<LI> ISTRES is the stress print option. ISTRES is set to I00 where:<DL COMPACT><DT>I <b> = 1</b><DD>: select all the elements,<DT>I <b> = 2</b><DD>: select the elements lying between numbers N1 and N2 given in INDIC,<DT>I <b> = 3</b><DD>: select the LINDIC elements whose numbers are given in INDIC,<DT>I <b> = 4</b><DD>: select the elements whose sub-domain numbers (LINDIC numbers) are given in INDIC,<P> </DL><LI> INDIC, a function of I, is the following array:<DL COMPACT><DT>I <b> = 1</b><DD>, the array is not used,<DT>I <b> = 2</b><DD>, length 2, the array contains the number of the first and last element to select,<DT>I <b> = 3</b><DD>, the number of the LINDIC elements to select,<DT>I <b> = 4</b><DD>, the number of the LINDIC sub-domains to treat. Stresses of all the elements withthese sub-domain numbers will be selected.<P> </DL><LI> LINDIC, a function of I, is the length of array INDIC,<LI> NFTAE, NITAE are the file number and the level of D.S. <b> TAE</b> resulting from STRESS,<LI> NFTAES, NITAES are the file number and the level of O.D.S. <b> TAE</b> containing the stresses of the selected elements <BR> (NFTAES <b>></b> 0 and NFTAES <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img23.gif"> NFTAE).</UL><P>The input file required is NFTAE.The output file is NFTAES.<P>There is no input data required specifically for this module.<P><P><P><HR SIZE=3 WIDTH="75%"><IMG SRC="../icons/smallmod.gif" WIDTH=211 HEIGHT=50 ALIGN=BOTTOM ALT="Modulef"><A NAME=tex2html683 HREF="node16.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/previous_motif.gif" ALT="previous"></A><A NAME=tex2html689 HREF="node11.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/up_motif.gif" ALT="up"></A><A NAME=tex2html691 HREF="node18.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/next_motif.gif" ALT="next"></A><A NAME=tex2html693 HREF="node2.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/contents_motif.gif" ALT="contents"></A><A HREF="../Guide6-18/node17.html"><IMG BORDER=0 SRC="../icons/zoom18.gif" ALIGN=BOTTOM ALT="[BIG]"></A><A HREF="../Guide6-14/node17.html"><IMG BORDER=0 SRC="../icons/zoom14.gif" ALIGN=BOTTOM ALT="[Normal]"></A><A HREF="../Guide6-10/node17.html"><IMG BORDER=0 SRC="../icons/zoom10.gif" ALIGN=BOTTOM ALT="[small]"></A><BR><B> Next: </B> <A NAME=tex2html692 HREF="node18.html">2.7 Computing the flux and temperature </A><B>Up: </B> <A NAME=tex2html690 HREF="node11.html">2 Manipulation of solutions</A><B> Prev: </B> <A NAME=tex2html684 HREF="node16.html">2.5 Computation of interface and reaction </A><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 + -