|
|
| Line 1: |
Line 1: |
| − | <!-- metadata commented in wiki content
| |
| | | | |
| | + | == Abstract == |
| | | | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| |
| − | '''NUMERICAL MODELING OF SOLID ELEMENTS OF STRUCTURAL CARBON STEEL, STAINLESS STEEL AND INCONEL ALLOY 718 VIA FINITE ELEMENT METHOD'''</div>
| |
| | | | |
| − | <div id="_GoBack" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| |
| − | L. D. de Almeida <sup>1</sup>, M. A. A. Lima <sup>2</sup> e P. A. S. Rocha <sup>3</sup></div>
| |
| | | | |
| − | Federal University of Ouro Preto, Department of Civil Engineering, Campus Morro do Cruzeiro, CEP35400–000, Ouro Preto – MG, Brazil, 55–31–35591546.
| + | == Full document == |
| − | | + | <pdf>Media:Draft_Almeida_708527476-4349-document.pdf</pdf> |
| − | Electronic mails: [mailto:lukermandouglas@yahoo.com.br lukermandouglas@yahoo.com.br] <sup>1</sup>, [mailto:miguel_angelo05@yahoo.com.br miguel_angelo05@yahoo.com.br]<sup>2</sup>, [mailto:paulorocha@em.ufop.br paulorocha@em.ufop.br] <sup>3</sup>.
| + | |
| − | -->
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | '''ABSTRACT'''</div>
| + | |
| − | | + | |
| − | The objective of this work is to carry out numerical analyzes of solid elements related to structural carbon steel, both for medium strength steel (MS250) and for high strength steel (HS350); stainless steel 304 (SS) and the inconel 718 superalloy (SI). For this, a mathematical formulation based on the Finite Element Method (FEM) was computationally implemented, which used the FORTRAN programming language, whose purpose was to obtain the values of stresses, deformations and nodal displacements of the steels and the alloy, taking into account that the materials present elastic-linear behavior. As for the stress analyses, the von Mises rupture criterion method was used. For the discretization of solid elements, 4-node tetrahedral (T4), 8-node hexahedral (H8) and 20-node hexahedral (H20) finite elements were used, with each node having three degrees of freedom. The Gauss-Legendre method (Gauss quadrature) was also used for the numerical resolution of the integrals. In order to verify and validate the responses obtained from the implemented computational program, comparisons were made with results found in the literature, and then numerical analyzes of the materials were carried out. This work intends to efficiently contribute to the calculation of stresses, strains and displacements in solids.
| + | |
| − | | + | |
| − | Keywords: Finite Element Method, Numerical Analysis, Fortran, Gauss-Legendre Quadrature.
| + | |
| − | | + | |
| − | =1. INTRODUCTION= | + | |
| − | | + | |
| − | The use of computational tools for solving structural analysis problems has become increasingly effective. Since the 1950s, the use of computational mechanics to solve existing phenomena in engineering has been increasingly used.
| + | |
| − | | + | |
| − | Throughout history, engineers and mathematicians have developed various methods for finding approximate solutions. The first methods were based on assumptions and simplifications that facilitated calculations, but often mischaracterized the problem and made the solution imprecise. With the invention of computers, the time and cost to perform a large number of operations was drastically reduced, which made numerical methods increasingly used and popular.
| + | |
| − | | + | |
| − | One of the analysis methods used in several researches in engineering areas to solve physical problems whose mathematical model is represented by partial differential equations is the Finite Element Method (FEM).
| + | |
| − | | + | |
| − | Since 1967, many books have been written about the Finite Element Method, with emphasis on Professor [1], and also [2], [3], [4] and [5]. During the same period, many journals presented papers on the method.
| + | |
| − | | + | |
| − | Therefore, the present research aims to carry out numerical (mechanical) analyzes of solid elements referring to structural carbon steel (MS250 and HS350) and stainless steel (SS304), in addition to the inconel superalloy (SI718). Such analyzes will be carried out from the numerical modeling of solid elements based on a computational program implemented in FORTRAN language [6], based on the Finite Element Method, in which three finite elements will be used for the modeling of the problems, namely: o 4-node tetrahedral finite element (T4), the 8-node hexahedral finite element (H8) and the 20-node hexahedral finite element (H20).
| + | |
| − | | + | |
| − | Thus, the aim of this research is to obtain the values of stresses, deformations and displacements along the solids under study, in which these will be subjected to certain types of loads. In order to validate the numerical results obtained with the aid of the developed computer program, the answers will be compared with results found in the literature.
| + | |
| − | | + | |
| − | =2. FORMULATIONS=
| + | |
| − | | + | |
| − | | + | |
| − | [[Image:Draft_Almeida_708527476-picture-Caixa de texto 7.svg|center|28px]]
| + | |
| − | Figure 1, shows the master elements of the 4-node tetrahedral finite element (T4) and 8-node hexahedral finite element (H8), respectively.
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | + | |
| − | | + | |
| − | {|
| + | |
| − | |-
| + | |
| − | | [[Image:Draft_Almeida_708527476-picture-Caixa de texto 13.svg|center|33px]]
| + | |
| − | | [[Image:Draft_Almeida_708527476-picture-Caixa de texto 11.svg|center|38px]]
| + | |
| − | |}
| + | |
| − | | + | |
| − | {|
| + | |
| − | |-
| + | |
| − | | [[Image:Draft_Almeida_708527476-picture-Caixa de texto 697.svg|center|33px]]
| + | |
| − | | [[Image:Draft_Almeida_708527476-picture-Caixa de texto 676.svg|center|28px]]
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | {|
| + | |
| − | |-
| + | |
| − | | [[Image:Draft_Almeida_708527476-picture-Caixa de texto 698.svg|center|39px]]
| + | |
| − | | [[Image:Draft_Almeida_708527476-image1.png|center|600px]]
| + | |
| − | |}
| + | |
| − | </div> | + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">'''Figure 1 – Master elements'''</span></div>
| + | |
| − | | + | |
| − | Next, the formulations of solid finite elements studied in this research are presented.
| + | |
| − | | + | |
| − | ===2.1 4-node Tetrahedral Finite Element ===
| + | |
| − | | + | |
| − | The shape functions corresponding to this element are:
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | <math display="inline">{N}_{1}=\, \xi</math> '' '' <math display="inline">{N}_{2}=</math><math>\, \eta</math> '' '' <math display="inline">{N}_{3}=\varphi</math>
| + | |
| − | | colspan='2' style="text-align: right;vertical-align: top;"|(1)
| + | |
| − | |-
| + | |
| − | | colspan='2' style="text-align: center;vertical-align: top;"|<math>{N}_{1}+{N}_{2}+{N}_{3}+{N}_{4}=1</math>
| + | |
| − | | style="text-align: right;vertical-align: top;"|(2)
| + | |
| − | |-
| + | |
| − | | style="vertical-align: top;"|
| + | |
| − | | colspan='2' style="text-align: right;vertical-align: top;"|
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | The displacement vector will be described as:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>q=\left[ {q}_{1}\quad {q}_{2}\quad {q}_{3}\quad {q}_{4}\quad {q}_{5}\quad {q}_{6}\quad {q}_{7}\quad {q}_{8}\quad {q}_{9}\quad {q}_{10}\quad {q}_{11}\quad {q}_{12}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(3)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | Being the relation between the vector of the field of displacements and the vector of nodal displacements given by:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{u=N\, \, q}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref114239928'></span>(4)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | Where N is the matrix representing the shape functions, given by:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{N=\, }\left[ \begin{matrix}{N}_{1}&0\\0&{N}_{1}\\0&0\end{matrix}\quad \begin{matrix}0&{N}_{2}\\0&0\\{N}_{1}&0\end{matrix}\quad \begin{matrix}0&0\\{N}_{2}&0\\0&{N}_{2}\end{matrix}\quad \begin{matrix}{N}_{3}&0\\0&{N}_{3}\\0&0\end{matrix}\quad \begin{matrix}0&{N}_{4}\\0&0\\{N}_{3}&0\end{matrix}\quad \begin{matrix}0&0\\{N}_{4}&0\\0&{N}_{4}\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref114239940'></span>(5)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | Then, with the help of Eq. (4) and Eq. (5), it is possible to conclude that:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>u={N}_{1}{q}_{1}+{N}_{2}{q}_{4}+{N}_{3}{q}_{7}+{N}_{4}{q}_{10}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(6)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>v={N}_{1}{q}_{2}+{N}_{2}{q}_{5}+{N}_{3}{q}_{8}+{N}_{4}{q}_{11}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(7)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>w={N}_{1}{q}_{3}+{N}_{2}{q}_{6}+{N}_{3}{q}_{9}+{N}_{4}{q}_{12}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(8)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | Since the function u depends on x, y and z, and that these depend on the natural coordinates ξ, η and φ, then the function u is also dependent on ξ, η and φ. However, there is:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\left\{ \begin{matrix}\frac{{\partial }_{u}}{{\partial }_{\xi }}\\\frac{{\partial }_{u}}{{\partial }_{\eta }}\\\frac{{\partial }_{u}}{{\partial }_{\varphi }}\end{matrix}\right\} =</math><math>\left[ \begin{matrix}\frac{{\partial }_{x}}{{\partial }_{\xi }}&\frac{{\partial }_{y}}{{\partial }_{\xi }}&\frac{{\partial }_{z}}{{\partial }_{\xi }}\\\frac{{\partial }_{x}}{{\partial }_{\eta }}&\frac{{\partial }_{y}}{{\partial }_{\eta }}&\frac{{\partial }_{z}}{{\partial }_{\eta }}\\\frac{{\partial }_{x}}{{\partial }_{\varphi }}&\frac{{\partial }_{y}}{{\partial }_{\varphi }}&\frac{{\partial }_{z}}{{\partial }_{\varphi }}\end{matrix}\right] \left\{ \begin{matrix}\frac{{\partial }_{u}}{{\partial }_{x}}\\\frac{{\partial }_{u}}{{\partial }_{y}}\\\frac{{\partial }_{u}}{{\partial }_{z}}\end{matrix}\right\}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(9)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | Since the Jacobian matrix is given by:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{J=}\left[ \begin{matrix}\frac{{\partial }_{x}}{{\partial }_{\xi }}&\frac{{\partial }_{y}}{{\partial }_{\xi }}&\frac{{\partial }_{z}}{{\partial }_{\xi }}\\\frac{{\partial }_{x}}{{\partial }_{\eta }}&\frac{{\partial }_{y}}{{\partial }_{\eta }}&\frac{{\partial }_{z}}{{\partial }_{\eta }}\\\frac{{\partial }_{x}}{{\partial }_{\varphi }}&\frac{{\partial }_{y}}{{\partial }_{\varphi }}&\frac{{\partial }_{z}}{{\partial }_{\varphi }}\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref114242322'></span>(10)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | Considering that the matrix '''A''' is the inverse matrix of the Jacobian matrix, it comes:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{A=}{\boldsymbol{J}}^{-1}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(11)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | Getting to:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\left\{ \begin{matrix}\frac{{\partial }_{u}}{{\partial }_{x}}\\\frac{{\partial }_{u}}{{\partial }_{y}}\\\frac{{\partial }_{u}}{{\partial }_{z}}\end{matrix}\right\} =</math><math>\boldsymbol{A\, }\left\{ \begin{matrix}\frac{{\partial }_{u}}{{\partial }_{\xi }}\\\frac{{\partial }_{u}}{{\partial }_{\eta }}\\\frac{{\partial }_{u}}{{\partial }_{\varphi }}\end{matrix}\right\}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(12)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | We have that the relation between the strain vector and the displacement vector is:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{\epsilon =B\, q}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref114251660'></span>(13)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | Knowing that the strain vector is defined by:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{\epsilon =\, }{\left[ \begin{matrix}{\epsilon }_{x}&{\epsilon }_{y}&{\epsilon }_{z}\end{matrix}\quad \begin{matrix}{\gamma }_{zy}&{\gamma }_{zx}&{\gamma }_{yx}\end{matrix}\right] }^{T}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(14)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | After some mathematical manipulations it is possible to conclude that the matrix '''B''' is equal to:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{B}=\left[ \begin{matrix}\begin{matrix}{A}_{11}&0&0\\0&{A}_{21}&0\\0&0&{A}_{31}\end{matrix}\quad \begin{matrix}{A}_{12}&0&0\\0&{A}_{22}&0\\0&0&{A}_{32}\end{matrix}\quad \begin{matrix}{A}_{13}&0&0\\0&{A}_{23}&0\\0&0&{A}_{33}\end{matrix}\quad \begin{matrix}-{\tilde{A}}_{1}&0&0\\0&-{\tilde{A}}_{2}&0\\0&0&-{\tilde{A}}_{3}\end{matrix}\\\begin{matrix}0&{A}_{31}&{A}_{21}\\{A}_{31}&0&{A}_{11}\\{A}_{21}&{A}_{11}&0\end{matrix}\quad \begin{matrix}0&{A}_{12}&{A}_{22}\\{A}_{32}&0&{A}_{12}\\{A}_{22}&{A}_{12}&0\end{matrix}\quad \begin{matrix}0&{A}_{33}&{A}_{23}\\{A}_{33}&0&{A}_{13}\\{A}_{23}&{A}_{13}&0\end{matrix}\quad \begin{matrix}0&-{\tilde{A}}_{3}&-{\tilde{A}}_{2}\\-{\tilde{A}}_{3}&0&-{\tilde{A}}_{1}\\-{\tilde{A}}_{2}&-{\tilde{A}}_{1}&0\end{matrix}\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(15)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | Given that:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>-{\tilde{A}}_{1}=\left[ {A}_{11}+{A}_{12}+{A}_{13}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(16)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>-{\tilde{A}}_{2}=\left[ {A}_{21}+{A}_{22}+{A}_{23}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(17)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>-{\tilde{A}}_{3}=\left[ {A}_{31}+{A}_{32}+{A}_{33}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(18)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | The stiffness of the element can be obtained based on the internal strain energy equation, given by:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{U}_{e}=\frac{1}{2}\, {\boldsymbol{q}}^{T}\, {\boldsymbol{B}}^{T}\, \boldsymbol{D}\, \boldsymbol{B}\, \boldsymbol{q\, }\int_{e}^{}dV</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(19)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | where the element stiffness matrix will be defined by:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{\boldsymbol{k}}_{e}={V}_{e}\, {\boldsymbol{B}}^{T}\boldsymbol{\, D\, B}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(20)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | In turn, the body force (corresponding to its own weight) will be given by:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{\boldsymbol{f}}_{e}=\int_{0}^{1}\int_{0}^{1-\xi }\int_{0}^{1-\xi -\eta }{\boldsymbol{N}}^{T}\, \boldsymbol{f\, }\left| \mathrm{det}\,J\right| \boldsymbol{\, }d\varphi \, d\eta \, d\xi</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref114249134'></span>(21)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | ===2.2 8-node Hexahedral Finite Element ===
| + | |
| − | | + | |
| − | The Lagrange shape functions are represented as:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{i}=\frac{1}{8}\, \, (1+{\xi }_{i}\, \xi )(1+{\eta }_{i}\, \eta )(1+{\varphi }_{i}\, \varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(22)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | The element stiffness matrix corresponding to the hexahedral finite element with 8 nodes is defined as:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{\boldsymbol{k}}^{e}=\int_{-1}^{1}\int_{-1}^{1}\int_{-1}^{1}{\boldsymbol{B}}^{T}\boldsymbol{\, D\, B\, }\left| det\, J\right| \boldsymbol{\, }d\varphi \, d\eta \, d\xi \,</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref114781062'></span>(23)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | Remembering that the integrals will be solved numerically with the aid of the Gauss-Legendre Method (Gauss Quadrature).
| + | |
| − | | + | |
| − | Considering that the gamma matrix is the inverse matrix of the Jacobian matrix, we have:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{\Gamma =}{\boldsymbol{J}}^{-1}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref114301265'></span>(24)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | Hence, it comes to:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\left\{ \begin{matrix}\frac{{\partial }_{u}}{{\partial }_{x}}\\\frac{{\partial }_{u}}{{\partial }_{y}}\\\frac{{\partial }_{u}}{{\partial }_{z}}\end{matrix}\right\} =</math><math>\boldsymbol{\Gamma \, }\left\{ \begin{matrix}\frac{{\partial }_{u}}{{\partial }_{\xi }}\\\frac{{\partial }_{u}}{{\partial }_{\eta }}\\\frac{{\partial }_{u}}{{\partial }_{\varphi }}\end{matrix}\right\}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(25)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | The matrix B corresponding to the hexahedral finite element is represented as:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{B}=\boldsymbol{H}\, \, {\boldsymbol{\Gamma }}_{\boldsymbol{u}}\boldsymbol{\, DN}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(26)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | From Eq. (27) it is possible to conclude that:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\left[ \begin{matrix}\begin{matrix}\frac{{\partial }_{u}}{{\partial }_{\xi }}\\\frac{{\partial }_{u}}{{\partial }_{\eta }}\\\frac{{\partial }_{u}}{{\partial }_{\varphi }}\end{matrix}\\\begin{matrix}\frac{{\partial }_{v}}{{\partial }_{\xi }}\\\frac{{\partial }_{v}}{{\partial }_{\eta }}\\\frac{{\partial }_{v}}{{\partial }_{\varphi }}\end{matrix}\\\begin{matrix}\frac{{\partial }_{w}}{{\partial }_{\xi }}\\\frac{{\partial }_{w}}{{\partial }_{\eta }}\\\frac{{\partial }_{w}}{{\partial }_{\varphi }}\end{matrix}\end{matrix}\right] =</math><math>\boldsymbol{DN\, \, q}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref114303759'></span>(27)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | since the DN matrix has 9 rows and 24 columns and will be organized using the sub-matrices proposed in this work and presented as:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN1=}\left[ \begin{matrix}\frac{{\partial }_{N1}}{{\partial }_{\xi }}&0&0\\\frac{{\partial }_{N1}}{{\partial }_{\eta }}&0&0\\\frac{{\partial }_{N1}}{{\partial }_{\varphi }}&0&0\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N2}}{{\partial }_{\xi }}&0&0\\\frac{{\partial }_{N2}}{{\partial }_{\eta }}&0&0\\\frac{{\partial }_{N2}}{{\partial }_{\varphi }}&0&0\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(28)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN2=}\left[ \, \begin{matrix}\frac{{\partial }_{N3}}{{\partial }_{\xi }}&0&0\\\frac{{\partial }_{N3}}{{\partial }_{\eta }}&0&0\\\frac{{\partial }_{N3}}{{\partial }_{\varphi }}&0&0\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N4}}{{\partial }_{\xi }}&0&0\\\frac{{\partial }_{N4}}{{\partial }_{\eta }}&0&0\\\frac{{\partial }_{N4}}{{\partial }_{\varphi }}&0&0\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(29)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN3=}\left[ \, \, \begin{matrix}\frac{{\partial }_{N5}}{{\partial }_{\xi }}&0&0\\\frac{{\partial }_{N5}}{{\partial }_{\eta }}&0&0\\\frac{{\partial }_{N5}}{{\partial }_{\varphi }}&0&0\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N6}}{{\partial }_{\xi }}&0&0\\\frac{{\partial }_{N6}}{{\partial }_{\eta }}&0&0\\\frac{{\partial }_{N6}}{{\partial }_{\varphi }}&0&0\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(30)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN4=}\left[ \, \, \begin{matrix}\frac{{\partial }_{N7}}{{\partial }_{\xi }}&0&0\\\frac{{\partial }_{N7}}{{\partial }_{\eta }}&0&0\\\frac{{\partial }_{N7}}{{\partial }_{\varphi }}&0&0\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N8}}{{\partial }_{\xi }}&0&0\\\frac{{\partial }_{N8}}{{\partial }_{\eta }}&0&0\\\frac{{\partial }_{N8}}{{\partial }_{\varphi }}&0&0\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(31)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN5=}\left[ \, \, \begin{matrix}0&\frac{{\partial }_{N1}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N1}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N1}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N2}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N2}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N2}}{{\partial }_{\varphi }}&0\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(32)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN6=}\left[ \, \, \begin{matrix}0&\frac{{\partial }_{N3}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N3}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N3}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N4}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N4}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N4}}{{\partial }_{\varphi }}&0\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(33)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN7=}\left[ \, \, \begin{matrix}0&\frac{{\partial }_{N5}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N5}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N5}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N6}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N6}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N6}}{{\partial }_{\varphi }}&0\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(34)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN8=}\left[ \, \, \begin{matrix}0&\frac{{\partial }_{N7}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N7}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N7}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N8}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N8}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N1}}{{\partial }_{\varphi }}&0\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(35)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN9=}\left[ \, \, \begin{matrix}0&0&\frac{{\partial }_{N1}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N1}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N1}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N2}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N2}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N2}}{{\partial }_{\varphi }}\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(36)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN10=}\left[ \, \, \begin{matrix}0&0&\frac{{\partial }_{N3}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N3}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N3}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N4}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N4}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N4}}{{\partial }_{\varphi }}\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(37)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN11=}\left[ \, \, \begin{matrix}0&0&\frac{{\partial }_{N5}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N5}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N5}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N6}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N6}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N6}}{{\partial }_{\varphi }}\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(38)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN12=}\left[ \, \, \begin{matrix}0&0&\frac{{\partial }_{N7}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N7}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N7}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N8}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N8}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N8}}{{\partial }_{\varphi }}\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(39)
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | or in a compact form as:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN=\, }\left[ \begin{matrix}\boldsymbol{DN1}&\boldsymbol{\vdots }&\boldsymbol{DN2}\\\boldsymbol{\cdots }&\boldsymbol{\cdots }&\boldsymbol{\cdots }\\\begin{matrix}\boldsymbol{DN5}\\\boldsymbol{\cdots }\\\boldsymbol{DN9}\end{matrix}&\begin{matrix}\boldsymbol{\vdots }\\\boldsymbol{\cdots }\\\boldsymbol{\vdots }\end{matrix}&\begin{matrix}\boldsymbol{DN6}\\\boldsymbol{\cdots }\\\boldsymbol{DN10}\end{matrix}\end{matrix}\boldsymbol{\quad }\begin{matrix}\boldsymbol{\vdots }&\boldsymbol{DN3}&\boldsymbol{\vdots }\\\boldsymbol{\cdots }&\boldsymbol{\cdots }&\boldsymbol{\cdots }\\\begin{matrix}\boldsymbol{\vdots }\\\boldsymbol{\cdots }\\\boldsymbol{\vdots }\end{matrix}&\begin{matrix}\boldsymbol{DN7}\\\boldsymbol{\cdots }\\\boldsymbol{DN11}\end{matrix}&\begin{matrix}\boldsymbol{\vdots }\\\boldsymbol{\cdots }\\\boldsymbol{\vdots }\end{matrix}\end{matrix}\boldsymbol{\quad }\begin{matrix}\boldsymbol{DN4}\\\boldsymbol{\cdots }\\\begin{matrix}\boldsymbol{DN8}\\\boldsymbol{\cdots }\\\boldsymbol{DN12}\end{matrix}\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(40)
| + | |
| − | |-
| + | |
| − | | style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | Based on the equation referring to the deformation vector, it is possible to conclude that:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{\epsilon =H\, }\left[ \begin{matrix}\begin{matrix}\frac{{\partial }_{u}}{{\partial }_{x}}\\\frac{{\partial }_{u}}{{\partial }_{y}}\\\frac{{\partial }_{u}}{{\partial }_{z}}\end{matrix}\\\begin{matrix}\frac{{\partial }_{v}}{{\partial }_{x}}\\\frac{{\partial }_{v}}{{\partial }_{y}}\\\frac{{\partial }_{v}}{{\partial }_{z}}\end{matrix}\\\begin{matrix}\frac{{\partial }_{w}}{{\partial }_{x}}\\\frac{{\partial }_{w}}{{\partial }_{y}}\\\frac{{\partial }_{w}}{{\partial }_{z}}\end{matrix}\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(41)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | where the matrix H will be expressed by:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{H=}\left[ \begin{matrix}\begin{matrix}1&0&0\\0&0&0\\0&0&0\end{matrix}\quad \begin{matrix}0&0&0\\0&1&0\\0&0&0\end{matrix}\quad \begin{matrix}0&0&0\\0&0&0\\0&0&1\end{matrix}\\\begin{matrix}0&1&0\\0&0&0\\0&0&1\end{matrix}\quad \begin{matrix}1&0&0\\0&0&1\\0&0&0\end{matrix}\quad \begin{matrix}0&0&0\\0&1&0\\1&0&0\end{matrix}\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(42)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | Hence, from Eq. (43):
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\left[ \begin{matrix}\begin{matrix}\frac{{\partial }_{u}}{{\partial }_{x}}\\\frac{{\partial }_{u}}{{\partial }_{y}}\\\frac{{\partial }_{u}}{{\partial }_{z}}\end{matrix}\\\begin{matrix}\frac{{\partial }_{v}}{{\partial }_{x}}\\\frac{{\partial }_{v}}{{\partial }_{y}}\\\frac{{\partial }_{v}}{{\partial }_{z}}\end{matrix}\\\begin{matrix}\frac{{\partial }_{w}}{{\partial }_{x}}\\\frac{{\partial }_{w}}{{\partial }_{y}}\\\frac{{\partial }_{w}}{{\partial }_{z}}\end{matrix}\end{matrix}\right] \boldsymbol{=\, }{\boldsymbol{\Gamma }}_{u}\boldsymbol{\, \, }\left[ \begin{matrix}\begin{matrix}\frac{{\partial }_{u}}{{\partial }_{\xi }}\\\frac{{\partial }_{u}}{{\partial }_{\eta }}\\\frac{{\partial }_{u}}{{\partial }_{\varphi }}\end{matrix}\\\begin{matrix}\frac{{\partial }_{v}}{{\partial }_{\xi }}\\\frac{{\partial }_{v}}{{\partial }_{\eta }}\\\frac{{\partial }_{v}}{{\partial }_{\varphi }}\end{matrix}\\\begin{matrix}\frac{{\partial }_{w}}{{\partial }_{\xi }}\\\frac{{\partial }_{w}}{{\partial }_{\eta }}\\\frac{{\partial }_{w}}{{\partial }_{\varphi }}\end{matrix}\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref114781238'></span>(43)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | It is possible to arrive at '''Γ<sub>u</sub>''', defined by the following relation:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{\boldsymbol{\Gamma }}_{u}=\left[ \begin{matrix}\Gamma (\xi ,\eta ,\varphi )&\mathit{\boldsymbol{\vdots }}&\mathit{\boldsymbol{0}}\\\mathit{\boldsymbol{\cdots }}&\mathit{\boldsymbol{\cdots }}&\mathit{\boldsymbol{\cdots }}\\\begin{matrix}\mathit{\boldsymbol{0}}\\\mathit{\boldsymbol{\cdots }}\\\mathit{\boldsymbol{0}}\end{matrix}&\begin{matrix}\mathit{\boldsymbol{\vdots }}\\\mathit{\boldsymbol{\cdots }}\\\mathit{\boldsymbol{\vdots }}\end{matrix}&\begin{matrix}\Gamma (\xi ,\eta ,\varphi )\\\mathit{\boldsymbol{\cdots }}\\\mathit{\boldsymbol{0}}\end{matrix}\end{matrix}\mathit{\boldsymbol{\quad }}\begin{matrix}\mathit{\boldsymbol{\vdots }}&\mathit{\boldsymbol{0}}\\\mathit{\boldsymbol{\cdots }}&\mathit{\boldsymbol{\cdots }}\\\begin{matrix}\mathit{\boldsymbol{\vdots }}\\\mathit{\boldsymbol{\cdots }}\\\mathit{\boldsymbol{\vdots }}\end{matrix}&\begin{matrix}\mathit{\boldsymbol{0}}\\\mathit{\boldsymbol{\cdots }}\\\Gamma (\xi ,\eta ,\varphi )\end{matrix}\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(44)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | Since the matrix Γ(ξ,η,φ) has 3 rows and 3 columns.
| + | |
| − | | + | |
| − | ===2.3 20-node Hexahedral Finite Element ===
| + | |
| − | | + | |
| − | Figure 2 shows the master cube corresponding to the 20-node hexahedral finite element.
| + | |
| − | | + | |
| − | <div id="_Toc115359752" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | '''Figure 2 - 20-node Hexahedral Master Element'''</div>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | [[Image:Draft_Almeida_708527476-image2.png|402px]] </div>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Source: Author</span></div>
| + | |
| − | | + | |
| − | Analogously to the 8-node hexahedral finite element, the shape functions of the 20-node hexahedral finite element will be defined as:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{1}=\frac{1}{8}\, \, \left( 1-\, \xi \right) \left( 1-\, \eta \right) \left( 1-\right. </math><math>\left. \, \varphi \right) (-2-\xi -\eta -\varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: bottom;width: 5px;text-align: right;white-space: nowrap;"|(45)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{2}=\frac{1}{4}\, \, (1-{\xi }^{2})(1-\, \eta )(1-\, \varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: bottom;width: 5px;text-align: right;white-space: nowrap;"|(46)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{3}=\frac{1}{8}\, \, \left( 1+\, \xi \right) \left( 1-\, \eta \right) \left( 1-\right. </math><math>\left. \, \varphi \right) (-2+\xi -\eta -\varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: bottom;width: 5px;text-align: right;white-space: nowrap;"|(47)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{4}=\frac{1}{4}\, \, (1+\xi )(1-{\eta }^{2})(1-\varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: bottom;width: 5px;text-align: right;white-space: nowrap;"|(48)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{5}=\frac{1}{8}\, \, \left( 1+\, \xi \right) \left( 1+\, \eta \right) \left( 1-\right. </math><math>\left. \, \varphi \right) (-2+\xi +\eta -\varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: bottom;width: 5px;text-align: right;white-space: nowrap;"|(49)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{6}=\frac{1}{4}\, \, (1-{\xi }^{2})(1+\, \eta )(1-\, \varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: bottom;width: 5px;text-align: right;white-space: nowrap;"|(50)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{7}=\frac{1}{8}\, \, \left( 1-\, \xi \right) \left( 1+\, \eta \right) \left( 1-\right. </math><math>\left. \, \varphi \right) (-2-\xi +\eta -\varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: bottom;width: 5px;text-align: right;white-space: nowrap;"|(51)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{8}=\frac{1}{4}\, \, (1-\xi )(1-{\eta }^{2})(1-\varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: bottom;width: 5px;text-align: right;white-space: nowrap;"|(52)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{9}=\frac{1}{4}\, \, (1-\xi )(1-\eta )(1-{\varphi }^{2})</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: bottom;width: 5px;text-align: right;white-space: nowrap;"|(53)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{10}=\frac{1}{4}\, \, (1+\xi )(1-\eta )(1-{\varphi }^{2})</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: bottom;width: 5px;text-align: right;white-space: nowrap;"|(54)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{11}=\frac{1}{4}\, \, (1+\xi )(1+\eta )(1-{\varphi }^{2})</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(55)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{12}=\frac{1}{4}\, \, (1-\xi )(1+\eta )(1-{\varphi }^{2})</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: bottom;width: 5px;text-align: right;white-space: nowrap;"|(56)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{13}=\frac{1}{8}\, \, \left( 1-\, \xi \right) \left( 1-\, \eta \right) \left( 1+\right. </math><math>\left. \, \varphi \right) (-2-\xi -\eta +\varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: bottom;width: 5px;text-align: right;white-space: nowrap;"|(57)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{14}=\frac{1}{4}\, \, (1-{\xi }^{2})(1-\, \eta )(1+\, \varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: bottom;width: 5px;text-align: right;white-space: nowrap;"|(58)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{15}=\frac{1}{8}\, \, \left( 1+\, \xi \right) \left( 1-\, \eta \right) \left( 1+\right. </math><math>\left. \, \varphi \right) (-2+\xi -\eta +\varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(59)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{16}=\frac{1}{4}\, \, (1+\xi )(1-{\eta }^{2})(1+\varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(60)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{17}=\frac{1}{8}\, \, \left( 1+\, \xi \right) \left( 1+\, \eta \right) \left( 1+\right. </math><math>\left. \, \varphi \right) (-2+\xi +\eta +\varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(61)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{18}=\frac{1}{4}\, \, (1-{\xi }^{2})(1+\, \eta )(1+\, \varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(62)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{19}=\frac{1}{8}\, \, \left( 1-\, \xi \right) \left( 1+\, \eta \right) \left( 1+\right. </math><math>\left. \, \varphi \right) (-2-\xi +\eta +\varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(63)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{N}_{20}=\frac{1}{4}\, \, (1-\xi )(1-{\eta }^{2})(1+\varphi )</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(64)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | The matrix '''B''' corresponding to the hexahedral finite element is represented as:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{B}=\boldsymbol{H}\, \, {\boldsymbol{\Gamma }}_{\boldsymbol{u}}\boldsymbol{\, DN}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(65)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | From Equation (66) it is possible to conclude that:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\left[ \begin{matrix}\begin{matrix}\frac{{\partial }_{u}}{{\partial }_{\xi }}\\\frac{{\partial }_{u}}{{\partial }_{\eta }}\\\frac{{\partial }_{u}}{{\partial }_{\varphi }}\end{matrix}\\\begin{matrix}\frac{{\partial }_{v}}{{\partial }_{\xi }}\\\frac{{\partial }_{v}}{{\partial }_{\eta }}\\\frac{{\partial }_{v}}{{\partial }_{\varphi }}\end{matrix}\\\begin{matrix}\frac{{\partial }_{w}}{{\partial }_{\xi }}\\\frac{{\partial }_{w}}{{\partial }_{\eta }}\\\frac{{\partial }_{w}}{{\partial }_{\varphi }}\end{matrix}\end{matrix}\right] =</math><math>\boldsymbol{DN\, \, q}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref114309277'></span>(66)
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | Hence, following the same valid considerations for the 8-node hexahedral finite element, it is possible to arrive at the '''DN''' matrix, which has 9 rows and 60 columns, being formed by the sub-matrices described below, and subsequently the matrix '''B.'''
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | | colspan='6' style="text-align: center;vertical-align: top;text-align: center;white-space: nowrap;"|<math>\boldsymbol{DN1=}\left[ \begin{matrix}\frac{{\partial }_{N1}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N1}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N1}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N2}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N2}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N2}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N3}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N3}}{{\partial }_{\eta }}&0\\\, 0&0&\frac{{\partial }_{N3}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N4}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N4}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N4}}{{\partial }_{\varphi }}\end{matrix}\right]</math>
| + | |
| − | | colspan='2' style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(67)
| + | |
| − | |-
| + | |
| − | | colspan='5' style="text-align: center;vertical-align: top;text-align: center;white-space: nowrap;"|<math>\boldsymbol{DN2=}\left[ \, \begin{matrix}\frac{{\partial }_{N5}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N5}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N5}}{{\partial }_{\varphi }}\end{matrix}\, \, \begin{matrix}\frac{{\partial }_{N6}}{{\partial }_{\xi }}&0&\begin{matrix}0&\quad \begin{matrix}\frac{{\partial }_{N7}}{{\partial }_{\xi }}&\, 0\end{matrix}\end{matrix}\\0&\frac{{\partial }_{N6}}{{\partial }_{\eta }}&\, \, \begin{matrix}0&\quad \begin{matrix}0&\quad \frac{{\partial }_{N7}}{{\partial }_{\eta }}\end{matrix}\end{matrix}\\0&0&\begin{matrix}\frac{{\partial }_{N6}}{{\partial }_{\varphi }}&\begin{matrix}0&\, \, 0\end{matrix}\end{matrix}\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N8}}{{\partial }_{\xi }}&\begin{matrix}0\quad &\, 0\end{matrix}\\0&0&\begin{matrix}\frac{{\partial }_{N8}}{{\partial }_{\eta }}\quad &0\, \end{matrix}\\\frac{{\partial }_{N7}}{{\partial }_{\varphi }}&0&\quad \, \begin{matrix}0&\quad \frac{{\partial }_{N8}}{{\partial }_{\varphi }}\end{matrix}\end{matrix}\right]</math>
| + | |
| − | | colspan='2' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(68)
| + | |
| − | |-
| + | |
| − | | colspan='2' style="text-align: center;vertical-align: top;text-align: center;white-space: nowrap;"|<math>\boldsymbol{DN3=}\left[ \, \begin{matrix}\frac{{\partial }_{N9}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N9}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N9}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N10}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N10}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N10}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N11}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N11}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N11}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N12}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N12}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N12}}{{\partial }_{\varphi }}\end{matrix}\, \right]</math>
| + | |
| − | | colspan='5' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(69)
| + | |
| − | |-
| + | |
| − | | colspan='4' style="text-align: center;vertical-align: top;text-align: center;white-space: nowrap;"|<math>\boldsymbol{DN4=}\left[ \, \begin{matrix}\frac{{\partial }_{N13}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N13}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N13}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N14}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N14}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N14}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N15}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N15}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N15}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N16}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N16}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N16}}{{\partial }_{\varphi }}\end{matrix}\, \right]</math>
| + | |
| − | | colspan='3' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(70)
| + | |
| − | |-
| + | |
| − | | colspan='3' style="text-align: center;vertical-align: top;text-align: center;white-space: nowrap;"|<math>\boldsymbol{DN5=}\left[ \, \begin{matrix}\frac{{\partial }_{N17}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N17}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N17}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N18}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N18}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N18}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N19}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N19}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N19}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}\frac{{\partial }_{N20}}{{\partial }_{\xi }}&0&0\\0&\frac{{\partial }_{N20}}{{\partial }_{\eta }}&0\\0&0&\frac{{\partial }_{N20}}{{\partial }_{\varphi }}\end{matrix}\, \right]</math>
| + | |
| − | | colspan='4' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(71)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN6=}\left[ \, \begin{matrix}0&\frac{{\partial }_{N1}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N1}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N1}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N2}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N2}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N2}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N3}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N3}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N3}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N4}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N4}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N4}}{{\partial }_{\varphi }}&0\end{matrix}\, \right]</math>
| + | |
| − | |}
| + | |
| − | | colspan='6' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(72)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN7=}\left[ \, \begin{matrix}0&\frac{{\partial }_{N5}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N5}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N5}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N6}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N6}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N6}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N7}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N7}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N7}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N8}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N8}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N8}}{{\partial }_{\varphi }}&0\end{matrix}\, \right]</math>
| + | |
| − | |}
| + | |
| − | | colspan='6' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(73)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN8=}\left[ \, \begin{matrix}0&\frac{{\partial }_{N9}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N9}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N9}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N10}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N10}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N10}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N11}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N11}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N11}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N12}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N12}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N12}}{{\partial }_{\varphi }}&0\end{matrix}\, \right]</math>
| + | |
| − | |}
| + | |
| − | | colspan='6' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(74)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN9=}\left[ \, \begin{matrix}0&\frac{{\partial }_{N13}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N13}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N13}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N14}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N14}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N14}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N15}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N15}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N15}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N16}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N16}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N16}}{{\partial }_{\varphi }}&0\end{matrix}\, \right]</math>
| + | |
| − | |}
| + | |
| − | | colspan='6' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(75)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN10=}\left[ \, \begin{matrix}0&\frac{{\partial }_{N17}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N17}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N17}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N18}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N18}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N18}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N19}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N19}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N19}}{{\partial }_{\varphi }}&0\end{matrix}\quad \begin{matrix}0&\frac{{\partial }_{N20}}{{\partial }_{\xi }}&0\\0&\frac{{\partial }_{N20}}{{\partial }_{\eta }}&0\\0&\frac{{\partial }_{N20}}{{\partial }_{\varphi }}&0\end{matrix}\, \right]</math>
| + | |
| − | |}
| + | |
| − | | colspan='6' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(76)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN11=}\left[ \, \begin{matrix}0&0&\frac{{\partial }_{N1}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N1}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N1}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N2}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N2}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N12}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N3}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N3}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N3}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N4}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N4}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N4}}{{\partial }_{\varphi }}\end{matrix}\, \right]</math>
| + | |
| − | |}
| + | |
| − | | colspan='6' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(77)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN12=}\left[ \, \begin{matrix}0&0&\frac{{\partial }_{N5}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N5}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N6}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N6}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N6}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N7}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N7}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N7}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N8}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N8}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N8}}{{\partial }_{\varphi }}\end{matrix}\, \right]</math>
| + | |
| − | |}
| + | |
| − | | colspan='6' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(78)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN13=}\left[ \, \begin{matrix}0&0&\frac{{\partial }_{N9}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N9}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N9}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N10}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N10}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N10}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N11}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N11}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N11}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N12}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N12}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N12}}{{\partial }_{\varphi }}\end{matrix}\, \right]</math>
| + | |
| − | |}
| + | |
| − | | colspan='6' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(79)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN14=}\left[ \, \begin{matrix}0&0&\frac{{\partial }_{N13}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N13}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N13}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N14}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N14}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N14}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N15}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N15}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N15}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N16}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N16}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N16}}{{\partial }_{\varphi }}\end{matrix}\, \right]</math>
| + | |
| − | |}
| + | |
| − | | colspan='6' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(80)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN15=}\left[ \, \begin{matrix}0&0&\frac{{\partial }_{N17}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N17}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N17}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N18}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N18}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N18}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N19}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N19}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N19}}{{\partial }_{\varphi }}\end{matrix}\quad \begin{matrix}0&0&\frac{{\partial }_{N20}}{{\partial }_{\xi }}\\0&0&\frac{{\partial }_{N20}}{{\partial }_{\eta }}\\0&0&\frac{{\partial }_{N20}}{{\partial }_{\varphi }}\end{matrix}\, \right]</math>
| + | |
| − | |}
| + | |
| − | | colspan='6' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(81)
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\boldsymbol{DN}\boldsymbol{=\, }\left[ \begin{matrix}\boldsymbol{DN1}&\boldsymbol{\vdots }&\boldsymbol{DN2}\\\boldsymbol{\cdots }&\boldsymbol{\cdots }&\boldsymbol{\cdots }\\\begin{matrix}\boldsymbol{DN6}\\\boldsymbol{\cdots }\\\boldsymbol{DN11}\end{matrix}&\begin{matrix}\boldsymbol{\vdots }\\\boldsymbol{\cdots }\\\boldsymbol{\vdots }\end{matrix}&\begin{matrix}\boldsymbol{DN7}\\\boldsymbol{\cdots }\\\boldsymbol{DN12}\end{matrix}\end{matrix}\boldsymbol{\quad }\begin{matrix}\boldsymbol{\vdots }&\boldsymbol{DN3}&\boldsymbol{\vdots }\\\boldsymbol{\cdots }&\boldsymbol{\cdots }&\boldsymbol{\cdots }\\\begin{matrix}\boldsymbol{\vdots }\\\boldsymbol{\cdots }\\\boldsymbol{\vdots }\end{matrix}&\begin{matrix}\boldsymbol{DN8}\\\boldsymbol{\cdots }\\\boldsymbol{DN13}\end{matrix}&\begin{matrix}\boldsymbol{\vdots }\\\boldsymbol{\cdots }\\\boldsymbol{\vdots }\end{matrix}\end{matrix}\boldsymbol{\quad }\begin{matrix}\begin{matrix}\boldsymbol{DN4}&\mathit{\boldsymbol{\vdots \, }}&\boldsymbol{DN5}\end{matrix}\\\begin{matrix}\mathit{\boldsymbol{\cdots }}&\mathit{\boldsymbol{\cdots }}&\mathit{\boldsymbol{\cdots }}\end{matrix}\\\begin{matrix}\begin{matrix}\boldsymbol{DN9}&\mathit{\boldsymbol{\, \, \vdots }}&\boldsymbol{DN10}\end{matrix}\\\begin{matrix}\mathit{\boldsymbol{\cdots }}&\mathit{\boldsymbol{\cdots }}&\mathit{\boldsymbol{\cdots }}\end{matrix}\\\begin{matrix}\boldsymbol{DN14}&\mathit{\boldsymbol{\vdots }}&\boldsymbol{DN15}\end{matrix}\end{matrix}\end{matrix}\right]</math>
| + | |
| − | |}
| + | |
| − | | colspan='6' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(82)
| + | |
| − | |-
| + | |
| − | | style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|
| + | |
| − | | colspan='6' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | ===2.4 Calculation of von Mises Stress===
| + | |
| − | | + | |
| − | For elements that are subject to plane stresses, it is possible to inform that the von Mises stress is represented by the following relation, since the yield stress must be greater than the calculated von Mises stress:
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{\sigma }_{VM}=\sqrt{{{\sigma }_{x}}^{2}+{{\sigma }_{y}}^{2}-{\sigma }_{x}{\sigma }_{y}+{{3\tau }_{xy}}^{2}}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(83)
| + | |
| − | |}
| + | |
| − | | + | |
| − | =3. APPLICATION=
| + | |
| − | | + | |
| − | Three examples will be presented using the implemented computational code. The mechanical properties of the following steels and alloys were used: for the structural carbon steel of medium and high mechanical strength (MS250 and HS350, respectively), the Modulus of Elasticity (E) of 200 GPa and the Poisson coefficient ( u) of 0.3; for stainless steel (SS304), the Modulus of Elasticity (E) of 193 GPa and the coefficient of Poisson (υ) of 0.27 and for the superalloy inconel (SI718), the Modulus of Elasticity (E) of 206 GPa and the Poisson coefficient (υ) of 0.28.
| + | |
| − | | + | |
| − | For the yield stress of the materials, the following values were used: 250 MPa and 350 MPa for structural carbon steel of medium and high mechanical strength, respectively; 215 MPa for stainless steel and 820 MPa for inconel alloy.
| + | |
| − | | + | |
| − | ===3.1 Example 1: 4-node tetrahedral element ===
| + | |
| − | | + | |
| − | The example shown in Fig. 3 refers to the modeling of the solid that was discretized with only a 4-node tetrahedral finite element.
| + | |
| − | | + | |
| − | <span id='_Ref114663382'></span><span id='_Ref114663365'></span><div id="_Toc115359754" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | '''Figure 3 - Example T4'''</div>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">''' [[Image:Draft_Almeida_708527476-image3.png|252px]] '''</span></div>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Source: Author</span></div>
| + | |
| − | | + | |
| − | For the analysis of the nodal displacement in the solid above, the coordinates of the nodes, given in millimeters, were defined as:
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 1 (0.25.25) 2 (0.0.25) 3 (25.0.25) 4 (0.0.0)</div>
| + | |
| − | | + | |
| − | In order to validate the developed computational program, Tab. 1 the comparison between the result obtained by the research and that found in the literature.
| + | |
| − | | + | |
| − | <span id='_Ref114664053'></span><span id='_Ref114664049'></span><div id="_Toc115359717" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | '''Table 1 – Comparison of T4 Nodal Displacent in z'''</div>
| + | |
| − | | + | |
| − | {| style="width: 71%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''Nodes'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''Present work'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''[7]'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Present work /'''
| + | |
| − | | + | |
| − | '''Literatura'''
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''1'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|-0.01341
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|-0.01340
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|1.00075
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''2'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|0.00000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|0.00000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|1.00000
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''3'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|0.00000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|0.00000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|1.00000
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''4'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|0.00000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|0.00000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|1.00000
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Source: Author</span></div>
| + | |
| − | | + | |
| − | Table 2 presents the results obtained by the implemented code for the materials used in this research.
| + | |
| − | | + | |
| − | <span id='_Ref114664366'></span><div id="_Toc115359718" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | '''Table 2 - T4 Nodal Displacent in z'''</div>
| + | |
| − | | + | |
| − | {| style="width: 67%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''Nodes'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''MS250'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''HS350'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''SS304'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''SI718'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''1'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|-0.0139
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|-0.0139
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|-0.0141
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|-0.0133
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''Nodes'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''MS250'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''HS350'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''SS304'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''SI718'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''2'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|0.0000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.0000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|0.0000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.0000
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''3'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|0.0000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.0000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|0.0000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.0000
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''4'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|0.0000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.0000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|0.0000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.0000
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <div id="_Toc115361972" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Source: Author</span></div>
| + | |
| − | | + | |
| − | ===3.2 Example 2: Corbel with straight voute (H8)===
| + | |
| − | | + | |
| − | A corbel with straight voute discretized with only one hexahedral finite element of 8 nodes (H8) was statically analyzed. In Fig 3, it is possible to verify the dimensions of the structure, points of application of loads and fixed nodes, in addition to the discretization to be used.
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|'''Figure 4 - Corbel with straight voute'''
| + | |
| − | | + | |
| − | [[Image:Draft_Almeida_708527476-image4.png|306px]]
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Source: Author</span></div>
| + | |
| − | | + | |
| − | Table 3 shows the displacements of the most requested nodes obtained by this work on the z axis.
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | '''Table 3 - Corbel with straight voute – Nodal displacements in z'''</div>
| + | |
| − | | + | |
| − | {| style="width: 77%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''Nodes'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''MS250'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''HS350'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''SS304'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''SI718'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''1'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00000
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''4'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.04364
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.04364
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.04472
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.04207
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''7'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.04725
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.04725
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.04857
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.04564
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''8'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.04725
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.04725
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.04857
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.04564
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Source: Author</span></div>
| + | |
| − | | + | |
| − | The table shows the results obtained for the highest values of von Mises stresses.
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | '''Table 4 - Corbel with straight voute - Von Mises stresses'''</div>
| + | |
| − | | + | |
| − | {| style="width: 79%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|'''Strain'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''MS250'''
| + | |
| − | | + | |
| − | '''(N/mm²)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|'''HS350'''
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | '''(N/mm²)'''</div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | '''SS304'''</div>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | '''(N/mm²)'''</div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | '''SI718'''</div>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | '''(N/mm²)'''</div>
| + | |
| − | | + | |
| − | | + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <math display="inline">{\sigma }_\mbox{1}</math> </div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 16.45000</div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 16.45000</div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 16.45000</div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 16.54000</div>
| + | |
| − | | + | |
| − | | + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <math>{\sigma }_\mbox{2}</math> </div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 18.01000</div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 18.01000</div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 18.01000</div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 18.17000</div>
| + | |
| − | | + | |
| − | | + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <math display="inline">{\sigma }_\mbox{7}</math> </div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 16.45000</div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 16.45000</div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 16.45000</div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 16.54000</div>
| + | |
| − | | + | |
| − | | + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <math display="inline">{\sigma }_\mbox{8}</math> </div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 18.01000</div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 18.01000</div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 18.01000</div>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | 18.17000</div>
| + | |
| − | | + | |
| − | | + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Source: Author</span></div>
| + | |
| − | | + | |
| − | ===3.3 Example 3: Corbel with parabolic voute (H20)===
| + | |
| − | | + | |
| − | Finally, a parabolic voute corbel discretized with only a 20-node hexahedral finite element (H20) was statically analyzed. In Fig. 5, it is possible to verify the dimensions of the structure, points of application of the loads, crimped nodes, in addition to the numbering used for the nodes to carry out the discretization of the structure via the Finite Element Method.
| + | |
| − | | + | |
| − | <span id='_Ref115200172'></span><span id='_Ref115200168'></span><div id="_Toc115359759" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | '''Figure 5 – Corbel with parabolic voute'''</div>
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;width: 100%;"|[[Image:Draft_Almeida_708527476-image9.png|324px]]
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Source: Author</span></div>
| + | |
| − | | + | |
| − | Table 5 shows the displacements, in the most requested nodes, obtained by the present work in the z axis.
| + | |
| − | | + | |
| − | <span id='_Ref115287106'></span><span id='_Toc115359745'></span>'''Table 5 – Corbel with parabolic voute - Nodal displacements in z'''
| + | |
| − | | + | |
| − | {| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''Nodes'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''MS250'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''HS350'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''SS304'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''SI718'''
| + | |
| − | | + | |
| − | '''(mm)'''
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|'''1'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00000
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|'''3'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00222
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00222
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00239
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00222
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|'''15'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.00191
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.00191
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.00195
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.00184
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|'''16'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.00392
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.00392
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.00413
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.00385
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|'''18'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.00148
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.00148
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.00158
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-0.00147
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;"|'''20'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00432
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00432
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00453
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.00424
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Source: Author</span></div>
| + | |
| − | | + | |
| − | Finally, Tab. 6 and Tab. 7 show the results obtained for the von Mises stresses and the elementary stresses, respectively.
| + | |
| − | | + | |
| − | <span id='_Ref115287219'></span><div id="_Toc115359746" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | '''Table 6 - Corbel with parabolic voute - von Mises stresses'''</div>
| + | |
| − | | + | |
| − | {| style="width: 79%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''Strain'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''MS250'''
| + | |
| − | | + | |
| − | '''(N/mm²)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''HS350'''
| + | |
| − | | + | |
| − | '''(N/mm²)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|'''SS304'''
| + | |
| − | | + | |
| − | '''(N/mm²)'''
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''SI718'''
| + | |
| − | | + | |
| − | '''(N/mm²)'''
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;width: 15%;"|<math display="inline">{\sigma }_\mbox{1}</math>
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.76680
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.76680
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.75990
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.76640
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;width: 15%;"|<math>{\sigma }_\mbox{2}</math>
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.84610
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.84610
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.91860
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.89950
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;width: 15%;"|<math display="inline">{\sigma }_\mbox{3}</math>
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2.49200
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2.49200
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2.59500
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2.58600
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;width: 15%;"|<math display="inline">{\sigma }_\mbox{4}</math>
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|4.00300
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|4.00300
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|4.18400
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|4.15100
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;width: 15%;"|<math display="inline">{\sigma }_\mbox{5}</math>
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2.94500
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2.94500
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2.84600
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2.88800
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;width: 15%;"|<math display="inline">{\sigma }_\mbox{6}</math>
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2.81200
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2.81200
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|3.07700
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|3.01100
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;width: 15%;"|<math display="inline">{\sigma }_\mbox{7}</math>
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|13.99000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|13.99000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|15.30000
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|14.86000
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: bottom;width: 15%;"|<math display="inline">{\sigma }_\mbox{8}</math>
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|1.86800
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|1.86800
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|1.99200
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|1.96700
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Source: Author</span></div>
| + | |
| − | | + | |
| − | =4. CONCLUSIONS=
| + | |
| − | | + | |
| − | It is noteworthy that the computational program implemented for the 4-node tetrahedral 4 element, when demonstrated in example 1, presented the values of the nodal displacements with a percentage difference of 0.00075% when compared with the values found in the literature. As for the solids analyzed in this work, the lowest nodal displacement occurred for the superalloy inconel 718, with stainless steel having the highest nodal displacement.
| + | |
| − | | + | |
| − | For the examples in which the solids were solved with the 8-node and 20-node hexahedral finite elements, the integrals were solved numerically with the aid of the Gauss-Legendre Method (Gauss Quadrature). When analyzing the solids with the materials under study, it was found that for the analysis of nodal displacements, the superalloy inconel 718 obtained the lowest value of nodal displacement and stainless steel, in general, had the highest value. of displacement. In addition, it is added that in the case of structures with curvilinear geometries, H<sub>20</sub> would be the most recommended finite element to carry out the analysis.
| + | |
| − | | + | |
| − | For such analyses, it is possible to highlight how much the Modulus of Elasticity and the Poisson coefficient influence the resistance capacity of a solid when it is requested by a load, in addition to showing that the Modulus of Elasticity is inversely proportional to its deformability, thus showing that the greater the Modulus of Elasticity, the less the structure becomes deformed.
| + | |
| − | | + | |
| − | Since in the present research, for the analysis of stresses, the von Mises method was used for all examples, it is concluded that none of the analyzed materials came to fail, due to this, the fact that the stress values of von Mises obtained did not exceed the yield stress values of the analyzed materials.
| + | |
| − | | + | |
| − | It is concluded, therefore, that the developed implementation was satisfactory, contributing with precise values of tension, deformation and displacement in solids subjected to a certain type of loading.
| + | |
| − | | + | |
| − | <span id='_Toc115361979'></span>
| + | |
| − | | + | |
| − | ==REFERENCES ==
| + | |
| − | | + | |
| − | [1] ZIENKIEWICZ, O. C., La Methods' des elements fínits (translated from the English), Pluralis, France, 1976.
| + | |
| − | | + | |
| − | [2] GALLAGHER, R. H., Introduction aux elements finis(translated from the English by J .L. Claudon), Pluralis, France, 1976.
| + | |
| − | | + | |
| − | [3] ABSI, E., Methode de calcul numérique en elaticité, Eyrolles, France, 1978.
| + | |
| − | | + | |
| − | [4] ROCKEY, K.C., Evans,H.R., Griffiths, D.W. and Nethercot, D.A., Elements finis, (translated from the English by C. Gomes), Eyrolles, France, 1979.
| + | |
| − | | + | |
| − | [5] IMBERT, J. F., Analyse des structures par elémentsfinis, Cepadues Ed., France, 1979.
| + | |
| − | | + | |
| − | [6] CHAPMAN, S. J., 2007, “Fortran 90/95 for Scientists and Engineers”, McGraw-Hill, 2nd ed.
| + | |
| − | | + | |
| − | [7] Chandrupatla, T.R., Belegundu, A.D., 2012, “Introduction to Finite Elements in Engineering”, 4rd ed. Pearson Education Limited, Edinburgh Gate, Reino Unido.
| + | |
