📄 node16.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.5 Computation of interface and reaction forces</TITLE><LINK REL=STYLESHEET TYPE="text/css" HREF="./Modulef.css" TITLE="Modulef CSS"><meta name="description" value="2.5 Computation of interface and reaction forces"><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=tex2html672 HREF="node15.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/previous_motif.gif" ALT="previous"></A><A NAME=tex2html678 HREF="node11.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/up_motif.gif" ALT="up"></A><A NAME=tex2html680 HREF="node17.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/next_motif.gif" ALT="next"></A><A NAME=tex2html682 HREF="node2.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/contents_motif.gif" ALT="contents"></A><A HREF="../Guide6-18/node16.html"><IMG BORDER=0 SRC="../icons/zoom18.gif" ALIGN=BOTTOM ALT="[BIG]"></A><A HREF="../Guide6-14/node16.html"><IMG BORDER=0 SRC="../icons/zoom14.gif" ALIGN=BOTTOM ALT="[Normal]"></A><A HREF="../Guide6-10/node16.html"><IMG BORDER=0 SRC="../icons/zoom10.gif" ALIGN=BOTTOM ALT="[small]"></A><BR><B> Next: </B> <A NAME=tex2html681 HREF="node17.html">2.6 Computing the stresses (D.S. TAE)</A><B>Up: </B> <A NAME=tex2html679 HREF="node11.html">2 Manipulation of solutions</A><B> Prev: </B> <A NAME=tex2html673 HREF="node15.html">2.4 Changing the type of a </A><B><A HREF="node2.html" >Contents</A></B><HR SIZE=3 WIDTH="75%"><H1><A NAME=SECTION03350000000000000000>2.5 Computation of interface and reaction forces</A></H1><P><P><P>In general, the stresses computed by the finite element method are not contained in one element or another.The stress vectors (normal force, tangential force) along an interface (a line or a surface)infers that the stresses are not in equilibrium on this interface (the actionis different from the reaction).<P>Preprocessor <b> FOINXX</b> constructs the D.S. <b> B</b> of interface and reaction forces obtained bycombining the element stress arraysof 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, obtained at the same time the solution of theelasticity problem was computed.<P>Module FOINRE, called by <b> FOINXX</b> computes the interface and reaction forces inall or part of the elements in the mesh. In the case of a two-dimensional problem, structure <b> B</b> resulting from FOINRE is used to visualize the interface andreaction forces via<b> TRMCXX</b> (see part II of this user guide). For the three-dimensional case, structure <b> B</b> resultingfrom FOINRE is visualized via <b> TRC3XX</b> (see part II as before).<P><P><P>Program FOINRE has the following form:<P><PRE> SUBROUTINE FOINRE(M,ISTRES,LINDIC,INDIC,NFTAE,NITAE,NFB,NIB, + NFMAIL,NIMAIL,NFNDL1,NINDL1,NFB2,NIB2,NTB2,NOT1,NOT2)C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++C AIM : COMPUTE THE INTERFACE AND REACTIONAL FORCESC --- RESTORE HE D.S. AND ADDRESS THE WORK ARRAYSC ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++</PRE><P>where<P><UL><LI> M is the work array,<LI> ISTRES is an option for selecting those elements to be treated, with:<DL COMPACT><DT>1<DD>: select all the elements,<DT>2<DD>: select the elements lying between numbers N1 and N2 given in INDIC,<DT>3<DD>: select the LINDIC elements whose numbers are given in INDIC,<DT>4<DD>: select the elements whose sub-domain numbers (LINDIC numbers) are given in INDIC,<P> </DL><LI> LINDIC, a function of ISTRES, is the length of array INDIC,<LI> INDIC, a function of ISTRES, is the following array:<DL COMPACT><DT>ISTRES <b> = 1</b><DD>, the array is not used,<DT>ISTRES <b> = 2</b><DD>, length 2, this array contains the number of the first and last element to treat,<DT>ISTRES <b> = 3</b><DD>, the numbers of the LINDIC elements to treat,<DT>ISTRES <b> = 4</b><DD>, the numbers of the LINDIC sub-domains to treat. The elements which have these sub-domain numbers will then be selected.<P> </DL><LI> NFTAE, NITAE are the file number and level of D.S. <b> TAE</b> containing the element arrays. NFTAE <b>></b> 0.<LI> NFB, NIB are the file number and level of D.S. <b> B</b> containing the displacements (NFB <b>></b> 0: thefile is sequential, NFBE = 0, the D.S. is in main memory),<LI> NFMAIL, NIMAIL are the file number and level of D.S. <b> MAIL</b>,<LI> NFNDL1, NINDL1, used if ND = 0, are the the file number and level of elasticity D.S. <b> NDL1</b> (ND = 0 indicates that the number of degrees of freedom per node is not constant),<LI> NFB2, NIB2, NTB2 are the file number, the level and number of arrays to be associated to the O.D.S. <b> B</b>which will contain the interface forces of the selected elements (NFB2 <IMG BORDER=0 ALIGN=MIDDLE ALT="" SRC="img22.gif"> 0),<LI> NOT1 is the rank in <b> TAE</b> of element stiffness array,<LI> NOT2 is the rank in <b> TAE</b> of the element RHS array,</UL><P>The input files necessary are NFB, NFTAE, NFMAIL and, in certain cases, NFNDL1.The output file is NFB2.<P>There is no data to input specifically for this module except if NTB2 is non-zero.In this case, the following input is required;<P><P><P>Loop I from 1 to NTB2<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>The output structure, <b> B</b>, is described in [<A HREF="node109.html#guide2"><A NAME=tex2html333 HREF="../Guide2/welcome.html">MODULEF User Guide - 2</A></A>].<P><P><P><HR SIZE=3 WIDTH="75%"><IMG SRC="../icons/smallmod.gif" WIDTH=211 HEIGHT=50 ALIGN=BOTTOM ALT="Modulef"><A NAME=tex2html672 HREF="node15.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/previous_motif.gif" ALT="previous"></A><A NAME=tex2html678 HREF="node11.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/up_motif.gif" ALT="up"></A><A NAME=tex2html680 HREF="node17.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/next_motif.gif" ALT="next"></A><A NAME=tex2html682 HREF="node2.html"><IMG BORDER=0 ALIGN=BOTTOM SRC="../icons/contents_motif.gif" ALT="contents"></A><A HREF="../Guide6-18/node16.html"><IMG BORDER=0 SRC="../icons/zoom18.gif" ALIGN=BOTTOM ALT="[BIG]"></A><A HREF="../Guide6-14/node16.html"><IMG BORDER=0 SRC="../icons/zoom14.gif" ALIGN=BOTTOM ALT="[Normal]"></A><A HREF="../Guide6-10/node16.html"><IMG BORDER=0 SRC="../icons/zoom10.gif" ALIGN=BOTTOM ALT="[small]"></A><BR><B> Next: </B> <A NAME=tex2html681 HREF="node17.html">2.6 Computing the stresses (D.S. TAE)</A><B>Up: </B> <A NAME=tex2html679 HREF="node11.html">2 Manipulation of solutions</A><B> Prev: </B> <A NAME=tex2html673 HREF="node15.html">2.4 Changing the type of a </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 + -