\documentclass[11pt,twoside]{article}\makeatletter

\IfFileExists{xcolor.sty}%
  {\RequirePackage{xcolor}}%
  {\RequirePackage{color}}
\usepackage{colortbl}
\usepackage{wrapfig}
\usepackage{ifxetex}
\ifxetex
  \usepackage{fontspec}
  \usepackage{xunicode}
  \catcode`⃥=\active \def⃥{\textbackslash}
  \catcode`❴=\active \def❴{\{}
  \catcode`❵=\active \def❵{\}}
  \def\textJapanese{\fontspec{Noto Sans CJK JP}}
  \def\textChinese{\fontspec{Noto Sans CJK SC}}
  \def\textKorean{\fontspec{Noto Sans CJK KR}}
  \setmonofont{DejaVu Sans Mono}
  
\else
  \IfFileExists{utf8x.def}%
   {\usepackage[utf8x]{inputenc}
      \PrerenderUnicode{–}
    }%
   {\usepackage[utf8]{inputenc}}
  \usepackage[english]{babel}
  \usepackage[T1]{fontenc}
  \usepackage{float}
  \usepackage[]{ucs}
  \uc@dclc{8421}{default}{\textbackslash }
  \uc@dclc{10100}{default}{\{}
  \uc@dclc{10101}{default}{\}}
  \uc@dclc{8491}{default}{\AA{}}
  \uc@dclc{8239}{default}{\,}
  \uc@dclc{20154}{default}{ }
  \uc@dclc{10148}{default}{>}
  \def\textschwa{\rotatebox{-90}{e}}
  \def\textJapanese{}
  \def\textChinese{}
  \IfFileExists{tipa.sty}{\usepackage{tipa}}{}
\fi
\def\exampleFont{\ttfamily\small}
\DeclareTextSymbol{\textpi}{OML}{25}
\usepackage{relsize}
\RequirePackage{array}
\def\@testpach{\@chclass
 \ifnum \@lastchclass=6 \@ne \@chnum \@ne \else
  \ifnum \@lastchclass=7 5 \else
   \ifnum \@lastchclass=8 \tw@ \else
    \ifnum \@lastchclass=9 \thr@@
   \else \z@
   \ifnum \@lastchclass = 10 \else
   \edef\@nextchar{\expandafter\string\@nextchar}%
   \@chnum
   \if \@nextchar c\z@ \else
    \if \@nextchar l\@ne \else
     \if \@nextchar r\tw@ \else
   \z@ \@chclass
   \if\@nextchar |\@ne \else
    \if \@nextchar !6 \else
     \if \@nextchar @7 \else
      \if \@nextchar (8 \else
       \if \@nextchar )9 \else
  10
  \@chnum
  \if \@nextchar m\thr@@\else
   \if \@nextchar p4 \else
    \if \@nextchar b5 \else
   \z@ \@chclass \z@ \@preamerr \z@ \fi \fi \fi \fi
   \fi \fi  \fi  \fi  \fi  \fi  \fi \fi \fi \fi \fi \fi}
\gdef\arraybackslash{\let\\=\@arraycr}
\def\@textsubscript#1{{\m@th\ensuremath{_{\mbox{\fontsize\sf@size\z@#1}}}}}
\def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}}
\def\abbr{}
\def\corr{}
\def\expan{}
\def\gap{}
\def\orig{}
\def\reg{}
\def\ref{}
\def\sic{}
\def\persName{}\def\name{}
\def\placeName{}
\def\orgName{}
\def\textcal#1{{\fontspec{Lucida Calligraphy}#1}}
\def\textgothic#1{{\fontspec{Lucida Blackletter}#1}}
\def\textlarge#1{{\large #1}}
\def\textoverbar#1{\ensuremath{\overline{#1}}}
\def\textquoted#1{‘#1’}
\def\textsmall#1{{\small #1}}
\def\textsubscript#1{\@textsubscript{\selectfont#1}}
\def\textxi{\ensuremath{\xi}}
\def\titlem{\itshape}
\newenvironment{biblfree}{}{\ifvmode\par\fi }
\newenvironment{bibl}{}{}
\newenvironment{byline}{\vskip6pt\itshape\fontsize{16pt}{18pt}\selectfont}{\par }
\newenvironment{citbibl}{}{\ifvmode\par\fi }
\newenvironment{docAuthor}{\ifvmode\vskip4pt\fontsize{16pt}{18pt}\selectfont\fi\itshape}{\ifvmode\par\fi }
\newenvironment{docDate}{}{\ifvmode\par\fi }
\newenvironment{docImprint}{\vskip 6pt}{\ifvmode\par\fi }
\newenvironment{docTitle}{\vskip6pt\bfseries\fontsize{22pt}{25pt}\selectfont}{\par }
\newenvironment{msHead}{\vskip 6pt}{\par}
\newenvironment{msItem}{\vskip 6pt}{\par}
\newenvironment{rubric}{}{}
\newenvironment{titlePart}{}{\par }

\newcolumntype{L}[1]{){\raggedright\arraybackslash}p{#1}}
\newcolumntype{C}[1]{){\centering\arraybackslash}p{#1}}
\newcolumntype{R}[1]{){\raggedleft\arraybackslash}p{#1}}
\newcolumntype{P}[1]{){\arraybackslash}p{#1}}
\newcolumntype{B}[1]{){\arraybackslash}b{#1}}
\newcolumntype{M}[1]{){\arraybackslash}m{#1}}
\definecolor{label}{gray}{0.75}
\def\unusedattribute#1{\sout{\textcolor{label}{#1}}}
\DeclareRobustCommand*{\xref}{\hyper@normalise\xref@}
\def\xref@#1#2{\hyper@linkurl{#2}{#1}}
\begingroup
\catcode`\_=\active
\gdef_#1{\ensuremath{\sb{\mathrm{#1}}}}
\endgroup
\mathcode`\_=\string"8000
\catcode`\_=12\relax

\usepackage[a4paper,twoside,lmargin=1in,rmargin=1in,tmargin=1in,bmargin=1in,marginparwidth=0.75in]{geometry}
\usepackage{framed}

\definecolor{shadecolor}{gray}{0.95}
\usepackage{longtable}
\usepackage[normalem]{ulem}
\usepackage{fancyvrb}
\usepackage{fancyhdr}
\usepackage{graphicx}
\usepackage{marginnote}

\renewcommand{\@cite}[1]{#1}


\renewcommand*{\marginfont}{\itshape\footnotesize}

\def\Gin@extensions{.pdf,.png,.jpg,.mps,.tif}

  \pagestyle{fancy}

\usepackage[pdftitle={Topology Optimization: Applications of VFLSM and SESO in the Generation of Three-Dimensional Strut-and-Tie Models},
 pdfauthor={}]{hyperref}
\hyperbaseurl{}

	 \paperwidth210mm
	 \paperheight297mm
              
\def\@pnumwidth{1.55em}
\def\@tocrmarg {2.55em}
\def\@dotsep{4.5}
\setcounter{tocdepth}{3}
\clubpenalty=8000
\emergencystretch 3em
\hbadness=4000
\hyphenpenalty=400
\pretolerance=750
\tolerance=2000
\vbadness=4000
\widowpenalty=10000

\renewcommand\section{\@startsection {section}{1}{\z@}%
     {-1.75ex \@plus -0.5ex \@minus -.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\Large\bfseries}}
\renewcommand\subsection{\@startsection{subsection}{2}{\z@}%
     {-1.75ex\@plus -0.5ex \@minus- .2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\Large}}
\renewcommand\subsubsection{\@startsection{subsubsection}{3}{\z@}%
     {-1.5ex\@plus -0.35ex \@minus -.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\large}}
\renewcommand\paragraph{\@startsection{paragraph}{4}{\z@}%
     {-1ex \@plus-0.35ex \@minus -0.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\normalsize}}
\renewcommand\subparagraph{\@startsection{subparagraph}{5}{\parindent}%
     {1.5ex \@plus1ex \@minus .2ex}%
     {-1em}%
     {\reset@font\normalsize\bfseries}}


\def\l@section#1#2{\addpenalty{\@secpenalty} \addvspace{1.0em plus 1pt}
 \@tempdima 1.5em \begingroup
 \parindent \z@ \rightskip \@pnumwidth 
 \parfillskip -\@pnumwidth 
 \bfseries \leavevmode #1\hfil \hbox to\@pnumwidth{\hss #2}\par
 \endgroup}
\def\l@subsection{\@dottedtocline{2}{1.5em}{2.3em}}
\def\l@subsubsection{\@dottedtocline{3}{3.8em}{3.2em}}
\def\l@paragraph{\@dottedtocline{4}{7.0em}{4.1em}}
\def\l@subparagraph{\@dottedtocline{5}{10em}{5em}}
\@ifundefined{c@section}{\newcounter{section}}{}
\@ifundefined{c@chapter}{\newcounter{chapter}}{}
\newif\if@mainmatter 
\@mainmattertrue
\def\chaptername{Chapter}
\def\frontmatter{%
  \pagenumbering{roman}
  \def\thechapter{\@roman\c@chapter}
  \def\theHchapter{\roman{chapter}}
  \def\thesection{\@roman\c@section}
  \def\theHsection{\roman{section}}
  \def\@chapapp{}%
}
\def\mainmatter{%
  \cleardoublepage
  \def\thechapter{\@arabic\c@chapter}
  \setcounter{chapter}{0}
  \setcounter{section}{0}
  \pagenumbering{arabic}
  \setcounter{secnumdepth}{6}
  \def\@chapapp{\chaptername}%
  \def\theHchapter{\arabic{chapter}}
  \def\thesection{\@arabic\c@section}
  \def\theHsection{\arabic{section}}
}
\def\backmatter{%
  \cleardoublepage
  \setcounter{chapter}{0}
  \setcounter{section}{0}
  \setcounter{secnumdepth}{2}
  \def\@chapapp{\appendixname}%
  \def\thechapter{\@Alph\c@chapter}
  \def\theHchapter{\Alph{chapter}}
  \appendix
}
\newenvironment{bibitemlist}[1]{%
   \list{\@biblabel{\@arabic\c@enumiv}}%
       {\settowidth\labelwidth{\@biblabel{#1}}%
        \leftmargin\labelwidth
        \advance\leftmargin\labelsep
        \@openbib@code
        \usecounter{enumiv}%
        \let\p@enumiv\@empty
        \renewcommand\theenumiv{\@arabic\c@enumiv}%
	}%
  \sloppy
  \clubpenalty4000
  \@clubpenalty \clubpenalty
  \widowpenalty4000%
  \sfcode`\.\@m}%
  {\def\@noitemerr
    {\@latex@warning{Empty `bibitemlist' environment}}%
    \endlist}

\def\tableofcontents{\section*{\contentsname}\@starttoc{toc}}
\parskip0pt
\parindent1em
\def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}}
\newenvironment{reflist}{%
  \begin{raggedright}\begin{list}{}
  {%
   \setlength{\topsep}{0pt}%
   \setlength{\rightmargin}{0.25in}%
   \setlength{\itemsep}{0pt}%
   \setlength{\itemindent}{0pt}%
   \setlength{\parskip}{0pt}%
   \setlength{\parsep}{2pt}%
   \def\makelabel##1{\itshape ##1}}%
  }
  {\end{list}\end{raggedright}}
\newenvironment{sansreflist}{%
  \begin{raggedright}\begin{list}{}
  {%
   \setlength{\topsep}{0pt}%
   \setlength{\rightmargin}{0.25in}%
   \setlength{\itemindent}{0pt}%
   \setlength{\parskip}{0pt}%
   \setlength{\itemsep}{0pt}%
   \setlength{\parsep}{2pt}%
   \def\makelabel##1{\upshape ##1}}%
  }
  {\end{list}\end{raggedright}}
\newenvironment{specHead}[2]%
 {\vspace{20pt}\hrule\vspace{10pt}%
  \phantomsection\label{#1}\markright{#2}%

  \pdfbookmark[2]{#2}{#1}%
  \hspace{-0.75in}{\bfseries\fontsize{16pt}{18pt}\selectfont#2}%
  }{}
      \def\TheFullDate{1970-01-01 (revised: 01 January 1970)}
\def\TheID{\makeatother }
\def\TheDate{1970-01-01}
\title{Topology Optimization: Applications of VFLSM and SESO in the Generation of Three-Dimensional Strut-and-Tie Models}
\author{}\makeatletter 
\makeatletter
\newcommand*{\cleartoleftpage}{%
  \clearpage
    \if@twoside
    \ifodd\c@page
      \hbox{}\newpage
      \if@twocolumn
        \hbox{}\newpage
      \fi
    \fi
  \fi
}
\makeatother
\makeatletter
\thispagestyle{empty}
\markright{\@title}\markboth{\@title}{\@author}
\renewcommand\small{\@setfontsize\small{9pt}{11pt}\abovedisplayskip 8.5\p@ plus3\p@ minus4\p@
\belowdisplayskip \abovedisplayskip
\abovedisplayshortskip \z@ plus2\p@
\belowdisplayshortskip 4\p@ plus2\p@ minus2\p@
\def\@listi{\leftmargin\leftmargini
               \topsep 2\p@ plus1\p@ minus1\p@
               \parsep 2\p@ plus\p@ minus\p@
               \itemsep 1pt}
}
\makeatother
\fvset{frame=single,numberblanklines=false,xleftmargin=5mm,xrightmargin=5mm}
\fancyhf{} 
\setlength{\headheight}{14pt}
\fancyhead[LE]{\bfseries\leftmark} 
\fancyhead[RO]{\bfseries\rightmark} 
\fancyfoot[RO]{}
\fancyfoot[CO]{\thepage}
\fancyfoot[LO]{\TheID}
\fancyfoot[LE]{}
\fancyfoot[CE]{\thepage}
\fancyfoot[RE]{\TheID}
\hypersetup{citebordercolor=0.75 0.75 0.75,linkbordercolor=0.75 0.75 0.75,urlbordercolor=0.75 0.75 0.75,bookmarksnumbered=true}
\fancypagestyle{plain}{\fancyhead{}\renewcommand{\headrulewidth}{0pt}}

\date{}
\usepackage{authblk}

\providecommand{\keywords}[1]
{
\footnotesize
  \textbf{\textit{Index terms---}} #1
}

\usepackage{graphicx,xcolor}
\definecolor{GJBlue}{HTML}{273B81}
\definecolor{GJLightBlue}{HTML}{0A9DD9}
\definecolor{GJMediumGrey}{HTML}{6D6E70}
\definecolor{GJLightGrey}{HTML}{929497} 

\renewenvironment{abstract}{%
   \setlength{\parindent}{0pt}\raggedright
   \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}}
   \vskip16pt
   \textcolor{GJBlue}{\large\bfseries\abstractname\space}
}{%   
   \vskip8pt
   \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}}
   \vskip16pt
}

\usepackage[absolute,overlay]{textpos}

\makeatother 
      \usepackage{lineno}
      \linenumbers
      
\begin{document}

             \author[1]{Hélio Luiz  Simonetti}

             \author[2]{Valério S.  Almeida}

             \author[3]{Virgil Del Duca  Almeida}

             \author[4]{Luttgardes de Oliveira  Neto}

             \author[5]{Marlan D. S.  Cutrim}

             \author[6]{Vitor Manuel A.  Leitão}

\renewcommand\Authands{ and }

\date{\small \em Received: 1 January 1970 Accepted: 1 January 1970 Published: 1 January 1970}

\maketitle


\begin{abstract}
        


This article presents the analysis of Strut-and-Tie Model (STM) in reinforced concrete 3D structures based on the study of topological optimization, so that the problem is formulated with the Smooth-ESO (SESO) discrete method, whose removal heuristic is bidirectional with discrete optimization procedure, and the Velocity Field Level Set Method (VFLSM), which is an inheritance of the classical continuum Level Set Method (LSM), but advances the design limits with a velocity field constructed from the rate of the design variables and base functions. The proposed approach is to couple both methods in conjunction with the Method of Moving Asymptotes (MMA), used to control the various design constraints that are the minimization of compliance and the Von Mises stress that has demonstrated more rational STM results. Additionally, it has been formulated a methodology for the automatic generation of optimal of 3D STM by using sensitivity analysis obtaining via derivatives of the Von Mises stress fields, finding the force paths prevailing compression in the directions of the strut and the tensile in the directions of the ties for the reinforcement insertion. All the codes are implemented with Matlab software and several comparison examples: Deep beam with opening, a pile cap, a bridge pier, and a single corbel, are presented to validate the present formulations and the results are compared with the literature.

\end{abstract}


\keywords{reinforced concrete, topology optimization, strut-and-tie model, SESO, VFLSM.}

\begin{textblock*}{18cm}(1cm,1cm) % {block width} (coords) 
\textcolor{GJBlue}{\LARGE Global Journals \LaTeX\ JournalKaleidoscope\texttrademark}
\end{textblock*}

\begin{textblock*}{18cm}(1.4cm,1.5cm) % {block width} (coords) 
\textcolor{GJBlue}{\footnotesize \\ Artificial Intelligence formulated this projection for compatibility purposes from the original article published at Global Journals. However, this technology is currently in beta. \emph{Therefore, kindly ignore odd layouts, missed formulae, text, tables, or figures.}}
\end{textblock*}


\begin{textblock*}{10cm}(1.05cm,3cm)
{{\textit{CrossRef DOI of original article:}} \underline{}}
\end{textblock*}\let\tabcellsep& 	 	 		 
\section[{I. Introduction}]{I. Introduction}\par
n the field of structural engineering, most concrete linear elements are designed by the classical theory of Bernoulli hypothesis. For a real physical analysis about behavior of these bending elements it is common to use the Strut-and-Tie Model (STM) that is a generalization of the classical analogy of the truss beam model. This analogy is shown by Ritter and Morsch at the beginning of the twentieth century, associated with the Reinforced Concrete (RC) beam in an equivalent truss structure (regions B, Fig.  {\ref 1}). The bar elements represent the fields of tensile and the compressed struts emerged inside the structural element as bending effects. The analogy has been improved and is still used by the technical standards in the design of reinforced concrete beams in flexural and shear force and laying down various criteria for determining safe limits in its procedures. However, the application of this hypothesis for any structural element can lead to over or under sizing of certain parts of the structure.\par
The Bernoulli hypothesis is valid for parts of the frame that there is no interference from other regions, such as sections near the columns, changing in geometry or other areas where the influence of strain due to shear efforts is not negligible. In this line, there are regions where the assumptions of Bernoulli do not adequately represent the bending structural behavior and the stress distribution. Structural elements such as beams, walls and pile caps and special areas such as beam-column connection, openings in beams and geometric discontinuities are examples. These regions, denominated "discontinuity regions D", are limited to distances of the dimension order of structural adjacent elements (Saint Venant's principle), that the shear stresses are applicable and the distribution of strains in the cross section is not linear. From the 80's, a Professor at the University of Stuttgart and other collaborators presented several researches that evaluated these regions more adequately, as \hyperref[b0]{[1]}, \hyperref[b1]{[2]}, \hyperref[b2]{[3]}, and other researchers as \hyperref[b3]{[4]}, \hyperref[b4]{[5]} and \hyperref[b5]{[6]}. The pioneering work by \hyperref[b0]{[1]} describes the STM more generally, covering the equivalent truss models and including these regions and special structural elements. The analogy used in the STM uses the same idea of the classical theory in order to define bars representing the flow of stress trying to create the shortest and more logical path loads. Several experimental evaluations have been studied to validate the STM applied to the RC design, as \hyperref[b6]{[7]}, \hyperref[b7]{[8]}, \hyperref[b8]{[9]} and \hyperref[b9]{[10]}. 
\section[{I}]{I}\par
The STM is recognized as a rational approach to the design of discontinuity regions and is incorporated in several current codes, such as ASCE-ACI 445 on Shear and Torsion \hyperref[b10]{[11]}, \hyperref[b11]{[12]}, \hyperref[b12]{[13]} and \hyperref[b13]{[14]}. These code provisions still require improvement due to uncertainties in the selection of optimal struts-and-ties, especially in the case of complex geometry or general load application conditions. Because of its simple model and needs the experience of the designer to select and distribute the elements of the model in order to represent the stresses path in a better way, it becomes evident the use of more reliable and automatic tools for defining its geometric and structural configuration. 
\section[{Fig. 1: D and B regions}]{Fig. 1: D and B regions}\par
To overcome these difficulties and improve the efficiency in building the optimal STM in RC structures, the theory of Topology Optimization (TO) has been used for two decades as an alternative and systematic approach consolidating itself as a fruitful path of design related research, once facilitates the shaping of materials under certain conditions. Many methods have been proposed for the solution of TO applied to STM, highlighting the use of the classical SIMP: \hyperref[b14]{[15]}, \hyperref[b15]{[16]}, \hyperref[b16]{[17]}, or ESO (Evolutionary Structural Optimization): \hyperref[b17]{[18]}, \hyperref[b18]{[19]}, \hyperref[b19]{[20]}, Liang et al. \hyperref[b20]{[21,}\hyperref[b21]{22,}\hyperref[b22]{23]}, Chen et al. \hyperref[b23]{[24]}, Zhong et al. \hyperref[b24]{[25]}, or variants, like BESO, Shobeiri et al. \hyperref[b25]{[26]}, RESO (Refined ESO), Leu et al. \hyperref[b26]{[27]} or SESO proposed by the present authors, Almeida et al. \hyperref[b27]{[28]}. SESO is based on the philosophy that if an element is not really necessary for the structure, its contribution to the structural stiffness is gradually diminished until it does not have any influence in the structure; that is, its removal is done smoothly, not radically as in the ESO method, that have been showed more efficient and robust and less sensitive to the discretization than ESO and faster than BESO, causing a decrease of the checkerboard formation.\par
In the last decade, the Level Set Method (LSM) has been highlighted in the field of TO, different from the conventional element wise density-based methods. LSM has clearer and smoother results and are flexible for complex topological changing, citing the pioneer's works of \hyperref[b28]{[29]}, \hyperref[b29]{[30]} and \hyperref[b30]{[31]}. The method describes the topological path by an implicit shape evolutive sequence by using a higher dimensional function to the design space for achieving the minimum energy under design constraints. Several other schemes have been included in the standard LSM to improve performance and achieve better results for general applications, like \hyperref[b31]{[32]}, \hyperref[b32]{[33]}. Wang and Kang \hyperref[b33]{[34,}\hyperref[b34]{35]} proposed the Velocity Field Level Set Method (VFLSM) which has been proved to be more efficient to deal with multiple constraints and design variables than LSM, but few works have been applied to STM by using VFLSM.\par
OT in solving problems in the field of 3D STM is not much explored for general D-regions, discouraged by the instabilities (checkerboard problem) inherent to SIMP, ESO/BESO or the complex formulation and high processing time of LSM/VFLSM. Thus, for stabilizing and accelerating the TO solution, several mathematical optimization methods have been proposed, such as Optimality Criteria, by Huang et al. \hyperref[b35]{[36]} with BESO, Augmented Lagrangian \hyperref[b36]{[37]} or \hyperref[b37]{[38]} with Level-Set, Lagrangian multiplier by \hyperref[b38]{[39]} and \hyperref[b39]{[40]} with LSM, or the Method of Moving Asymptotes (MMA), by \hyperref[b40]{[41]} with SIMP.\par
In the present work, aiming at the solution of 3D STM in general reinforced concrete problems, the SESO methods whose advantages are easy implementation and decrease of the checkerboard effect and the VFLSM, which deals well with shape and topological optimizations, are formulated together with the MMA optimization method to accelerate and stabilize 3D STM. It is also noteworthy new approach of sensitivity analysis is incorporated in these formulations for the automatic generation of struts-and-ties based on partial derivatives 
\section[{Global Journal of Researches in Engineering}]{Global Journal of Researches in Engineering}\par
© 2023 Global Journ als  {\ref ( )} with respect to Von Mises stresses. The volume constraints are considered in the analyses, as the implementation of a spatial filter and the conjugate gradient method with the incomplete Cholesky preconditioner to speed up the solution of the linear system of each step of the search. 
\section[{a) Problem Formulation}]{a) Problem Formulation}\par
Considering the classical topology problem for the maximum stiffness of statically loaded linear elastic structures, a TO mathematical formulation for continuum structure can be discussed. Considering the TO problem as minimizing the deformation energy of a given structure considering the equilibrium, it follows that W=2U. The problem can then be defined as: \hyperref[b0]{(1)} with ?? ?? being the element's elasticity matrix, ?? ?? is the element's strain vector, ?? ?? is the volume of an element, NE is the number of finite elements of the mesh, ?? is the stiffness matrix, ???? = ?? is the equilibrium equation, ?? is the vector of loads applied to the structure, ?? ?? is the design variable of the i-th element, ?? is the vector of design variables. 
\section[{b) Smooth Evolutionary Structural Optimization (SESO)}]{b) Smooth Evolutionary Structural Optimization (SESO)}\par
The ESO method, which heuristic is based on the gradual and systematic removal of elements whose contribution to the stiffness of the structure are insignificant, was proposed by Xie and Steven \hyperref[b41]{[42]}. The SESO method proposed by Simonetti et al. \hyperref[b42]{[43]} is based on the ESO philosophy and applies a weighting to the constitutive matrix so that the element that would be eliminated is maintained and receives a smoothing characteristic. This treatment procedure applies a degradation in the value of its initial stiffness in such, during the removal process, its influence can contribute and determine its permanence or its definitive withdrawal from the design domain. Thus, the elements located near the limit to the left of this maximum strain energy are kept in the structure, defining a smoother heuristic removal. In Fig.  {\ref 2}, ??(??) is the constitutive matrix of element j, ?? = ?? ???? + ?? ???? is the domain of elements that can be withdrawn, ?? ???? is the domain of elements that must be effectively removed, ?? ???? is the domain of elements that are returned to the structure, 0 ? ?????? ? 1 is a weighted function. Subject to:???? = ?? ??(??) = ? ?? ?? ?? ?? -?? ? 0 ???? ??=1 ?? = \{?? 1 ?? 2 ?? 3 ? ?? ?? \}, ?? ?? = 1 ?????? ?? = 0 ??(??, ??) = \{ ??(??, ??) > 0 ? ?? ? ??\textbackslash ???? ; ??(??, ??) = 0 ? ?? ? ???? ; ??(??, ??) < 0 ? ?? ? ??\textbackslash ?? \}\par
with ?????? ? \{(??, ??)???? 2 \} is any point in the design domain D and ?? is the solid domain boundary as shown in Fig. \hyperref[fig_2]{3} for a 2D case.\par
()\textbf{2}\par
In classical LSM for TO, such as \hyperref[b29]{[30]} and \hyperref[b30]{[31]}, the design evolution is based on the solution of the Hamilton-Jacobi partial differential Eq. (  {\ref 4}). Thus, it needs an appropriate choice of finite difference methods on a fixed cartesian mesh. In general, the design update involves differentiation, resetting and velocity extension.  Recently, Wang and Kang \hyperref[b33]{[34,}\hyperref[b34]{35]} proposed a 100-line Velocity Field Level Set (VFLS), implemented in Matlab code. The structural shape and topology are updated by a velocity field constructed with the base function and velocity design variables defined throughout the domain. Then, the velocity field determines the search direction of the shape and the topological evolution can be obtained by a generic mathematical programming algorithm, which makes it more convenient and efficient to deal with multiple constraints and types of design variables. For VFLS, we have:(5) ?? ?? (??) = ? ?? ?? ?? ?? (??) ?? ?? =1\par
with ?? ?? (?? = 1,2, ? , ??) are the velocity design variables at N velocity points distributed throughout the main design, and ?? ?? (??) are the basic functions. It is observed that when ?? ?? satisfies the properties of the Kronecker delta it has ?? ?? = ?? ?? of Eq.(3). 
\section[{II. Optimization Algorithm -Moving Asymptotes Method}]{II. Optimization Algorithm -Moving Asymptotes Method}\par
To accelerate and stabilize the present 3D STM in this paper, MMA is employed, which is a mathematical programming algorithm suitable for TO, capable of handling optimization of many constraints and design variables. At each step of the algorithm's iterative process, a convex approximation subproblem is generated and solved. The generation of these subproblems is controlled by the moving asymptotes, which can both stabilize and accelerate the convergence of the overall process, \hyperref[b44]{[44]}.\par
The optimal solution of the subproblem may or may not be accepted: if so, the outer iteration is completed; if not, a new inner iteration is performed, in which a new subproblem is generated and solved. The iterations are repeated until the values of the approximations of the objective function and the constraints become greater than or equal to the values of the original function when evaluated in the optimal solution of the subproblem, that is, until the conservative condition is satisfied for the functions involved. The approximations that characterize the MMA are rational functions whose asymptotes are updated at each iteration. It is noteworthy that the use of rational approximations is justified by the fact that in several structural engineering problems where reciprocal variables arise, that is, interaction and mutual effort, given the objective function or a constraint ??(??) , the approximation functions are given by:Global Journal of Researches in Engineering © 2023 Global Journ als ( ) E Volume Xx XIII Issue II V ersion I ??(??) ? ??(?? ?? ) + ? ? ?? ?? ?? ?? (?? ) -?? ?? + ?? ?? ?? ?? -?? ?? (?? ) ? ?? 1\par
where ?? ?? e ?? ?? are defined as:???? ????(??) ???? ?? > 0 ??????? ?? ?? = ??? ?? (?? ) -?? ?? (?? ) ? 2 ????(??) ???? ?? ?????? ?? ?? = 0 ???? ????(??) ???? ?? < 0 ??????? ?? ?? = -??? ?? (?? ) -?? ?? (?? ) ? 2 ????(??) ???? ?? ?????? ?? ?? = 0\par
For the optimization problem in compliance Eq. (  {\ref 8}), it is known that it is satisfied because???? (??) ???? ?? < 0.\par
Then the MMA provides the current design with an approximation of a linear programming problem of the type:\par
with?? ?? = \{??????|0.9?? ?? (?? ) + 0.1?? ?? (?? ) ? ?? ?? ? 0.9?? ?? (?? ) + 0.1?? ?? (?? ) \} ? ?? = 1,2, ? , ??\par
with ?? ?? (?? ) and ?? ?? (?? ) being lower and upper asymptotes, respectively, k is the current iteration, n the number of design variables, ?? ?? the design variable and ?? the prescribed volume. The following heuristic rule is used by \hyperref[b44]{[44]} for updating the asymptotes, for the first two outer iterations, when k =1 and k = 2 are adopted:?? ?? (?? ) + ?? ?? (?? ) = 2?? ?? (?? )\textbf{(10)}\par
?? ?? (?? ) -?? ?? (?? ) = 1\par
For ?? ? 3?? ?? (?? ) + ?? ?? (?? ) = 2?? ?? (?? ) ?? ?? (?? ) -?? ?? (?? ) = ?? ?? (?? ) with ?? ?? (?? ) = ? ????? ?? (?? ) -?? ?? (?? -1) ???? ?? (??-1) -?? ?? (?? -2) ? < 0 ????? ?? (?? ) -?? ?? (?? -1) ???? ?? (??-1) -?? ?? (?? -2) ? > 0 ????? ?? (?? ) -?? ?? (?? -1) ???? ?? (??-1) -?? ?? (?? -2) ? = 0\textbf{(11)}\par
where the values of ??, ?? and ?? were fitted in the respective numerical ranges 0.65 ? ?? ? 0.75 , 1.15 ? ?? ? 1.25 and 0.9 ? ?? ? 1.\par
It can be seen in Eq. ( \hyperref[formula_9]{11}) that the MMA saves the signal of three consecutive iterations. Thus, when the signals alternate, the MMA detects that the values of the design variables are oscillating, i.e., ??? ?? (?? ) -approximate the design point ?? ?? (?? ) . If the values of the design variables do not oscillate, i.e., ??? ?? (?? ) -?? ?? (?? -1) ???? ?? (??-1) -?? ?? (??-2) ? ? 0, then the MMA moves the asymptotes away from the design point in order to accelerate up convergence. There are two approaches to solving subproblems in MMA, the "dual approach" and the "primal-dual interior point approach". The dual approach is based on the dual Lagrangian relaxation corresponding to the subproblem, which seeks the maximization of a concave objective function without other constraints and the non-negativity condition on the variables. This dual problem can be solved by a modified Newton method, and then the dual optimal solution can be translated into a corresponding optimal solution of the primal subproblem, which is used in this paper. 
\section[{III. Methodology for Generation 3D Strut-and-Tie Models and the Final Flowchart}]{III. Methodology for Generation 3D Strut-and-Tie Models and the Final Flowchart}\par
To determine the path load of the 3D bars of the STM from the TO analysis, this paper employs a new procedure to evaluate the struts and ties by the signs of the derivatives of the Von Mises stress components. It is known that for 3D problems they can be written as (?? ?? ???? )² = 1 2 [(???? ?? (?? -1) ???? ?? (??-1) -?? ?? (??-2) ? , < 0 the asymptotes ???????????????? -? ??? ?? (??) -?? ?? (??) ? 2 ?? ?? -?? ?? (??) ???? ???? ?? ??? ? (??) ? ?? 1 ?????????????? ???? ?? ? ?? ?? -?? ? 0 ?? ? ?? ??\par
Taking the local calculation of the derivative of the von Mises stress of the element with respect to the components of the stress vector described respectively as:??(?? ?? ???? ) ???? ??11 = 1 2?? ?? ???? (2?? ??11 -?? ??22 -?? ??33 ) ??(?? ?? ???? ) ???? ??22 = 1 2?? ?? ???? (2?? ??22 -?? ??11 -?? ??33 ) ??(?? ?? ???? ) ???? ??33 = 1 2?? ?? ???? (2?? ??33 -?? ??11 -?? ??22 ) ??(?? ?? ???? ) ???? ??12 = 3?? ??12 2?? ?? ???? ??(?? ?? ???? ) ???? ??23 = 3?? ??23 2?? ?? ???? ??(?? ?? ???? ) ???? ??31 = 3?? ??31 2?? ?? ????\par
Considering Eq.(  {\ref 13}) and making??(?? ?? ???? (??)) ???? ??11\par
> 0 then the elements are preponderantly tensioned (blue color -ties) while??(?? ?? ???? (??)) ???? ??33\par
< 0 are preponderantly compressed (green color -strut). The flowchart presented in Fig.  {\ref 4} shows the original methodology presented in this section with the approach of using element sensitivity for automatic generation of STMs via stress derivatives, when a target volume is reached, the stopping criterion is reached. A set of techniques has not yet been presented in scientific articles on 3D models, so the results obtained in item 4 are compared with those proposed by \hyperref[b15]{[16]}, \hyperref[b25]{[26]} and \hyperref[b43]{[45]}. Highlights that the VFLSM method required a neighborhood filter to define the tensile (blue) and compression (green) regions. This filter is due to intermediate values that occur in continuous TO methods such as the intermediate densities that occur in the SIMP methodology. 
\section[{IV. Numerical Examples}]{IV. Numerical Examples}\par
The following examples of structures engineering focus on TO base on minimizing compliance for STMs. The geometry and boundary conditions for numerical applications are represented for each case. All numerical examples were processed on a Core i7-2370, 8th Gen notebook, 2.8 GHz CPU with 20.0 GB (RAM). 
\section[{a) Example 1 -Deep Beam with Opening}]{a) Example 1 -Deep Beam with Opening}\par
The example presents a simply supported deep beam with an opening at the bottom left corner. The beam has its span three times its height and it is defined in \hyperref[b46]{[46]}, where the simple bending structural behavior is no longer considered. A vertical downward force F=3000 kN is applied eccentrically at the top edge as shown in Fig.  {\ref 5}. The structure is discretized with a total of 65,420 hexahedral elements (SESO) and 65,420 tetrahedral elements (LSM) (Fig.  {\ref 5} shows the design domain and its boundary conditions). In this configuration, the force in off-center position and the opening positioned near the left low end create a situation that changes the internal stress flow in the structure, between the load and the supports. The tie elements, resulting from tensile stresses, are positioned at the extremities of the strut elements, resulting from compressive stresses, geometrically defining the final model. Fig. \hyperref[fig_4]{6} provides the optimal topologies of the optimization procedures for the SESO (Fig. \hyperref[fig_4]{6a} and Fig. \hyperref[fig_4]{6b}) and VFSLM (Fig. \hyperref[fig_4]{6c}) methods, with a final volume fraction equal to 32\%. The optimal configurations have similarity to the classical STM presented by \hyperref[b0]{[1]} and later by \hyperref[b19]{[20]}. The computational cost presented by SESO using Optimality Criteria \hyperref[b47]{[47]} is approximately 40\% lower than the SESO and VFLSM methods using the MMA. It can be also noticed in Fig. \hyperref[fig_4]{6} that the optimal settings obtained by the VFLSM formulation clearly defines distinct elements (strut or tie) near the lateral faces of the deep beam, resulting in a more discrete STM, compared to the optimal settings presented by the SESO method. The classic model, Fig. \hyperref[fig_5]{7}, denotes three diagonal struts starting from the region of load application, one of them external directed to the closest support, another contouring the opening and directed to the support, and a third internal one. The ends of the struts are connected by tie composing the final structure of the STM.   \hyperref[b43]{[45]} In Fig. \hyperref[fig_6]{8} it can be seen that the SESO formulation (Fig. \hyperref[fig_6]{8a}) results in a setting similar to the classical model, but the VFLSM formulation (Fig. \hyperref[fig_6]{8b}) presents a model with a discretely simpler setting, with the internal strut in the vertical direction, unifying at a lower point the two ties. This setting simplifies the design procedure and the reinforcement detailing of the reinforced concrete structure, in the practical and executive sense, although the classic model makes it possible to calculate the complementary reinforcement around the opening.  
\section[{b) Example 2 -Pile Cap}]{b) Example 2 -Pile Cap}\par
In this example, a building foundation structure is dimensioned as a pile cap according to the dimensions shown in Fig. \hyperref[fig_7]{9}, for consideration as a rigid block and to enable the analysis by the STM concept. The pile cap is subjected to a vertical force of 4,000 kN located at the center of the upper face. The material properties used are the compressive strength of the concrete cylinder is 32 MPa. The Young's modulus of the concrete E c = 25,000 MPa and Poisson's coefficient ?? = 0.15. The filter radius ?? ?????? = 1.5 ???? mm and the volume fraction of 22.5\%, a rejection ratio, RR = 1\% and the evolution ratio ER = 2\% were specified in the optimization process. 
\section[{Global Journal of Researches in Engineering}]{Global Journal of Researches in Engineering}\par
© 2023 Global Journ als In the numerical simulations, to discretize the domain of the structure, a refined mesh of 40x20x40 was used, totaling 32,000 hexahedral elements (SESO) with 1mm reference side was used and a mesh of 32,000 tetrahedral elements (VFLSM). The results obtained as final optimal topologies of this problem for these meshes are represented in Fig.  {\ref 10} and can be compared with the results with those presented by \hyperref[b15]{[16]} and \hyperref[b25]{[26]}, see Fig.  {\ref 11}.\par
The optimal topology is basically composed of discrete elements represented in the principal stress flows. These optimal settings are adequate to perform the detailing and dimensioning of the required reinforcement, as well as strength checks. In this structure, the vertical load is distributed in four struts inclined toward the supports represented by vertical piles. The models highlight elements at the base of the pile cap, representing the tensile stresses, where a plane frame of ties balances the strut ends generated by the 3D structure in both horizontal directions, Fig.  {\ref 10}, where it can be seen the optimum topologies for the two methods, SESO and VFLSM. In the automatic generation of the strut models, it was considered the main flows of distinct stresses by colors, where the region of compression struts is green color and the region of tensile ties is blue. Although the models result quite similar, when approaching this problem, one must consider the increased computational burden associated with a 3D structure; a solid mesh usually requires that many elements be investigated at an adequate level of detail, with notable consequences on the number of equations and variables. Seeking to minimize this aspect of the processing, the system of equations received the implementation of a sparse approximation preconditioner for the inverse matrix. With this routine active, the computational cost of SESO-3D for this problem was decreased from 8,000 sec to 1,854 sec (4.3 times less) while VFLSM had a decrease from 8,000 sec to 3,851sec (2.1 times less).\par
The dimensioning of the reinforcement of this model is performed, as already presented in \hyperref[b27]{[28]}; in the calculations of the dimensions of the model elements, namely, inclined compression strut -column-pile and horizontal tension tie -pile-pile, the geometry of the problem presented \hyperref[b25]{[26]} is used: In the inclined strut, the verification of the compressive stresses is performed according to and the area of the strut required for the design strength of the concrete not to be exceeded:?????????? ???? ?????????? -?? 1 =?? ???? ,?????????? = ?? ?? * ?? 1 0.8 * ?? ???? = 1.4 * 2,031 0.8 * (3.2/1.4) = 1.555 ????²\par
By way of comparison, in \hyperref[b25]{[26]}, the results of this sizing are ?? ????,?????? = 41.66 ???? 2 and ?? ????,?????????? = 1,659 ????². The differences in values (3.5\% and 6.3, respectively) are due to different calculation criteria between the technical standards used, but values of the same order of magnitude can be considered. Fig. \hyperref[fig_10]{12} shows the reinforcement arrangement for the pile cap. The SESO and VFLSM methods using the MMA as accelerator are applied to a structure representing a column receiving loading from the bridge superstructure, represented by four vertical forces, as shown in Fig. \hyperref[fig_11]{13}. The concrete material properties, rejection ratio (RR), evolutionary ratio (ER) and filter radius are the same as in the previous example. For the numerical simulations, in the SESO method the bridge support is discretized using a fine mesh of 85x55x20 hexahedral elements of eight nodes, with reference side of 1 mm, while in the VFLSM method the mesh used has 85x20x55, totaling 93,500 tetrahedral finite elements. 
\section[{Global Journal of Researches in Engineering}]{Global Journal of Researches in Engineering}\par
The compliance history and the performance of the methods during the optimization procedure are plotted in Fig.  {\ref 14}. It can be seen in Fig.  {\ref 14b} that the performance index perfectly captures the changes in compliance and increases from unity to a maximum value of 2.3, stabilizing quickly around 2.1, the value at the optimal iteration.\par
The history of the optimization procedure via SESO and VFLSM for the bridge pier are shown in Fig.  {\ref 1} and Fig. \hyperref[fig_4]{16}. The optimal topologies were achieved at iterations 82 and 100 with final volumes equal to 20\% of the initial volume and a computational cost of 4,315.8 sec for SESO while VFLSM showed a computational cost of 5,486.5 sec.  The optimal settings with highlights of the distinct regions by colors are presented in Fig. \hyperref[fig_5]{17} for the two methods proposed in this paper. In these representations, a vertical axial force is expected to balance the symmetric external loads in the region of the base constraint. The applied vertical forces, in fact, are transferred to the column axis by means of two inclined struts and two vertical struts that merge into two in the proximity of the top region of the vertical element, driving the load distribution to the lower region where are the base supports. Note that the SESO method creates a unified region at the base while the VFLSM method sets up two parallel vertical paths. In addition, a horizontal tensile tie is arranged at the top of the body receiving the applied forces, which ensures the "T" geometry of the structure and configures the struts equilibrium in the load application zones. From a numerical point of view, the result obtained is optimal and configures the symmetry defined by the position of the design load. For automatic generation of STM models in the VFLSM method, it was necessary to implement the derivatives of von Mises stresses in the code proposed by \hyperref[b33]{[34]}. stress flows (green) similar to those of the SESO method, Fig. \hyperref[fig_5]{17a}, highlighting the robustness of both methods for creating strut-and-tie models. With the objective of investigate the effects of D-regions, three holes were inserted in the horizontal element of the bridge pier structure, and the number of finite elements of the mesh was reduced to 88,700, as shown in Fig. \hyperref[fig_6]{18}. The optimal topologies of the SESO and VFLSM models are represented in Fig. \hyperref[fig_7]{19}, where the struts are represented by green color and the ties by blue color.\par
The optimum results obtained demonstrate that the presence of geometric discontinuities produces changes in the stress flows, that seek to contour the discontinuities, describing practically vertical struts in the horizontal body of the bridge pier from the points of load application. These struts bend below the openings to meet at the top of the vertical element, creating points of deviation that need to be equilibrated by tensile ties. In Fig. \hyperref[fig_7]{19}, it can be seen the representations of STM elements created as described.\par
This modification with the presence of the openings affects the STMs models significantly, and the real load transfer mechanism can change with the dimensions of the openings. The optimization histories are shown in Fig.  {\ref 20} and Fig.  {\ref 21}, by the SESO and VFLSM formulations, respectively. The SESO and VFLSM methods were also experimented with for modeling struts-and-ties in a single corbel attached to a column. A simple structure can eventually result in an intricate STM as the dimensions and load arrangements can be defined. The geometry and dimensions of the structure are shown in Fig.  {\ref 22}. This single corbel is subjected to a concentrated load of 1 kN. The compressive strength of the concrete used in this example is 32 MPa. Young's modulus of the concrete E = 28,567 MPa and the Poisson's ratio ?? = 0.15 were defined in the analysis. A prescribed fraction volume V = 0.22 m 3 and an evolution ratio of ER = 2\% was specified in the optimization process.\par
In the SESO method the structure was discretized with a mesh of 44x12x108 unit hexahedral finite elements. The performance of the structure was monitored throughout the optimization procedure and, despite the breaks in the load transfer mechanisms due to element removal, the structure did not present failure modes and the performance index remained higher than 1, stabilizing at 1.6. In the VFLSM method, the same mesh was used, totaling 57,024 tetrahedral elements. Figures \hyperref[fig_2]{23a, b, c} and {\ref d} show that the optimal topologies obtained by the two models are different and checkerboard patterns were not detected. It is noted in observation made in the deep beam example that both formulations, SESO and VFLSM, define settings differently for elements of strut-and-tie models. Discrete elements are configured on the side faces of the models 
\section[{Global Journal of Researches in Engineering}]{Global Journal of Researches in Engineering}\par
© 2023 Global Journ als ( ) in different regions, while complete planes are shaped in other regions, with no common convention between the two formulations. The presented results show that both SESO-3D and VFLSM-3D are able to provide the prediction of the load transfer mechanism in reinforced concrete structures, even considering the structural domain thickness in the configuration of the component elements of the models. The STMs are presented in Fig. \hyperref[fig_16]{24}, it can be seen that these models are different and capable of clearly representing the location of the struts, ties and nodal zones. These results can be compared with those presented by \hyperref[b15]{[16]} and \hyperref[b25]{[26]}. It is also highlighted that the parameters of the MMA optimizer were changed to ?? ? 0.98, ?? ? 1.25 e ?? ? 0.75 proportion a more feasible topology for design. Fig. \hyperref[fig_16]{24b} shows the optimal setting of the VFLSM used for automatic creation of the STM models; both formulations exhibit distinct tensile (blue) and compressed (green) regions, even in the width of the structural domain. Table \hyperref[tab_1]{1} highlights the computational cost of SESO and VFSLM in all the examples presented in this paper evidencing the better performance for SESO-OC and SESO-MMA compared to VFLSM-MMA.  This paper aimed to extend the application of TO in 3D elasticity to obtain the best solution to STM problems. It brought some processes as innovation, such as the use of the SESO method and the VFLSM employed in conjunction with the OC and MMA methods to accelerate and stabilize the analyses; so that, the first method demonstrated to be more efficient when employed with the SESO, about 2 to 3 times faster in all the examples evaluated. It is highlighted that in these processes the incorporation of the linear solution by the conjugate gradient method with the incomplete Cholesky preconditioner further enhanced the computational cost. In the automated generation of the final designs of the STM, the procedure of obtaining struts and ties computed by the partial derivatives of the stresses of each element was applied highlighting that this novelty is easy to implement and the use of a spatial modal filter in the stress field was enough to completely eliminate the checkerboard. From the automatic generations performed, it was possible to design an example according to the recurring norm in an expeditious manner, in which the required reinforcement areas were evaluated and compared, demonstrating a good similarity. All codes were implemented in the high level language Matlab, which is easily accessible and extensible for future incorporation of other more realistic models, such as a rheological model more suitable for concrete. The study of STM using optimization applied to both materials (steel and concrete), leading to dimensioning and detailing of RC structural elements under the reliability-based topology optimization (RBTO) paradigm, taking advantage of the efficiency and stability of the procedures, are the highlights in the formulations developed in this paper.\begin{figure}[htbp]
\noindent\textbf{211}\includegraphics[]{image-2.png}
\caption{\label{fig_0}Fig. 2 : 1 2 1 2}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-3.png}
\caption{\label{fig_1}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{3}\includegraphics[]{image-4.png}
\caption{\label{fig_2}Fig. 3 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{45}\includegraphics[]{image-5.png}
\caption{\label{fig_3}Fig. 4 :Fig. 5 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{6}\includegraphics[]{image-6.png}
\caption{\label{fig_4}Fig. 6 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{7}\includegraphics[]{image-7.png}
\caption{\label{fig_5}Fig. 7 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{8}\includegraphics[]{image-8.png}
\caption{\label{fig_6}Fig. 8 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{9}\includegraphics[]{image-9.png}
\caption{\label{fig_7}Fig. 9 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{1011}\includegraphics[]{image-10.png}
\caption{\label{fig_8}Fig. 10 :Fig. 11 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-11.png}
\caption{\label{fig_9}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{12}\includegraphics[]{image-12.png}
\caption{\label{fig_10}Fig. 12 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{13}\includegraphics[]{image-13.png}
\caption{\label{fig_11}Fig. 13 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{141516}\includegraphics[]{image-14.png}
\caption{\label{fig_12}Fig. 14 :Fig. 15 :Fig. 16 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-15.png}
\caption{\label{fig_13}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{17182021}\includegraphics[]{image-16.png}
\caption{\label{fig_14}Fig. 17 :Fig. 18 :Fig. 20 :Fig. 21 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{2223}\includegraphics[]{image-17.png}
\caption{\label{fig_15}Fig. 22 :Fig. 23 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{24}\includegraphics[]{image-18.png}
\caption{\label{fig_16}Fig. 24 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-19.png}
\caption{\label{figure19}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-20.png}
\caption{\label{figure20}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-21.png}
\caption{\label{figure21}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{} \par 
\begin{longtable}{P{0.1367132867132867\textwidth}P{0.7132867132867132\textwidth}}
\tabcellsep Year 2023\\
\tabcellsep 27\\
\tabcellsep Volume Xx XIII Issue II V ersion I\\
\tabcellsep ( ) E\\
\tabcellsep Global Journal of Researches in Engineering\\
2 + ?? 23 2 + ?? 31 2 )\tabcellsep (12)\\
\tabcellsep © 2023 Global Journ als\end{longtable} \par
 
\caption{\label{tab_0}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{1} \par 
\begin{longtable}{P{0.26857638888888885\textwidth}P{0.17118055555555556\textwidth}P{0.17413194444444444\textwidth}P{0.17708333333333334\textwidth}P{0.059027777777777776\textwidth}}
Structure\tabcellsep Computation al Cost (SESO-OC) (minute)\tabcellsep Computational (SESO-MMA) Cost (minute)\tabcellsep Computational (VFLSM) Cost (minute)\tabcellsep Strut-and-tie Models\\
Deep beam with opening\tabcellsep 90.7\tabcellsep 109.8\tabcellsep 229.81\tabcellsep \\
Pile Cap\tabcellsep 21.4\tabcellsep 38.0\tabcellsep 64.2\tabcellsep \\
Bridge pier\tabcellsep 57.9\tabcellsep 95.0\tabcellsep 159.3\tabcellsep \\
Bridge pier with three holes\tabcellsep 56.5\tabcellsep 95.9\tabcellsep 158.6\tabcellsep \\
Single Corbel\tabcellsep 45.4\tabcellsep 71.9\tabcellsep 103.7\tabcellsep \end{longtable} \par
 
\begin{quote}
\par
Global Journal of Researches in Engineering© 2023 Global Journ als ( ) E Volume Xx XIII Issue II V ersion I 38 Year 2023\end{quote}

\caption{\label{tab_1}Table 1 :}\end{figure}
 			\label{foot_0}\footnote{\label{foot_0} Topology Optimization: Applications of VFLSM and SESO in the Generation of Three-Dimensional Strut-and-Tie Models} 			\label{foot_1}\footnote{\label{foot_1} Year 2023 Topology Optimization: Applications of VFLSM and SESO in the Generation of Three-Dimensional Strut-and-Tie Models} 		 		\backmatter   			 
\subsection[{Acknowledgment}]{Acknowledgment}\par
The authors acknowledge the Sao Paulo State Research Foundation (FAPESP) under Grant Number 2016/02327-5 for their financial support and CNPq (National Council of Scientific and Technological Development) under Grant Numbers 305093/2018-5 and Federal Institute of Education, Science and Technology of Minas Gerais under Grant Numbers 23792.001327/2022-49. 			  			  				\begin{bibitemlist}{1}
\bibitem[Yamasaki et al. ()]{b37}\label{b37} 	 		‘A level set based topology optimization method using the discretized signed distance function as the design variables’.  		 			S Yamasaki 		,  		 			T Nomura 		,  		 			A Kawamoto 		,  		 			K Sato 		,  		 			K Izui 		,  		 			S Nishiwaki 		.  	 	 		\textit{Struct Multidisc Optim}  		2009. 41  (5)  p. .  	 
\bibitem[Luo et al. ()]{b31}\label{b31} 	 		‘A level set method for structural shape and topology optimization using radial basis functions’.  		 			Z Luo 		,  		 			L Tong 		,  		 			Z Kang 		.  	 	 		\textit{Comput Struct}  		2009. 87 p. .  	 
\bibitem[Wang et al. ()]{b29}\label{b29} 	 		‘A level set method for structural topology optimization’.  		 			M Y Wang 		,  		 			X Wang 		,  		 			D Guo 		.  	 	 		\textit{Comput Methods Appl Mech Eng}  		2003. 192 p. .  	 
\bibitem[Zhong et al. ()]{b24}\label{b24} 	 		‘A new evaluation procedure for the strut-and-tie models of the disturbed regions of reinforced concrete structures’.  		 			J T Zhong 		,  		 			L Wang 		,  		 			P Deng 		,  		 			M Zhou 		.  	 	 		\textit{Engineering Structures}  		2017. 148 p. .  	 
\bibitem[Xie and Steven ()]{b41}\label{b41} 	 		‘A simple evolutionary procedure for structural optimization’.  		 			Y M Xie 		,  		 			G P Steven 		.  	 	 		\textit{Computers \& Structures}  		1993. 49  (5)  p. .  	 
\bibitem[Simonetti et al. ()]{b42}\label{b42} 	 		‘A smooth evolutionary structural optimization procedure applied to plane stress problem’.  		 			H L Simonetti 		,  		 			V S Almeida 		,  		 			L O Neto 		.  	 	 		\textit{Eng Struct}  		2014. 75 p. .  	 
\bibitem[AASHTO LRFD Bridge Design Specifications; American Association of State Highway and Transportation Officials ()]{b11}\label{b11} 	 		\textit{AASHTO LRFD Bridge Design Specifications; American Association of State Highway and Transportation Officials},  		2007. Washington, DC, USA: American Association of State Highway and Transportation Officials.  	 
\bibitem[Vaquero and Bertero (2020)]{b17}\label{b17} 	 		‘Automatic Generation of Proper Strut-and-Tie Model’.  		 			Sebastian F Vaquero 		,  		 			Raul D Bertero 		.  	 	 		\textit{ACI Structural Journal}  		Nov. 2020. 117  (6)  p. 81.  	 
\bibitem[Ali and White (2001)]{b18}\label{b18} 	 		‘Automatic Generation of Truss Model for Optimal Design of Reinforced Concrete Structures’.  		 			M A Ali 		,  		 			R N White 		.  	 	 		\textit{ACI Structural Journal}  		July-Aug. 2001. 98  (4)  p. .  	 
\bibitem[Shobeiri and Ahmadi-Nedushan ()]{b25}\label{b25} 	 		‘Bi-directional evolutionary structural optimization for strut-and-tie modelling of three-dimensional structural concrete’.  		 			V Shobeiri 		,  		 			B Ahmadi-Nedushan 		.  	 	 		\textit{Engineering Optimization}  		2017. 49  (12)  p. .  	 
\bibitem[Almeida et al. ()]{b27}\label{b27} 	 		‘Comparative analysis of strut-and-tie models using Smooth Evolutionary Structural Optimization’.  		 			V S Almeida 		,  		 			H L Simonetti 		,  		 			L O Neto 		.  	 	 		\textit{Engineering Structures}  		2013. 56 p. .  	 
\bibitem[Bogomolny and Amir ()]{b16}\label{b16} 	 		‘Conceptual design of reinforced concrete structures using topology V. Conclusions © 2023 Global Journ als optimization with elastoplastic material modeling’.  		 			M Bogomolny 		,  		 			O Amir 		.  	 	 		\textit{Int. J. Numer. Meth. Engng}  		2012. 90 p. .  	 
\bibitem[Dunning et al. ()]{b39}\label{b39} 	 		‘Coupled aerostructural topology optimization using a level set method for 3D aircraft wings’.  		 			P D Dunning 		,  		 			B K Stanford 		,  		 			H A Kim 		.  	 	 		\textit{Struct Multidiscip Optim}  		2014. 51  (5)  p. .  	 
\bibitem[Schlaich and Schäfer ()]{b1}\label{b1} 	 		‘Design and detailing of structural concrete using strut-and-tie models’.  		 			J Schlaich 		,  		 			K Schäfer 		.  	 	 		\textit{Struct Eng}  		1991. 69  (6)  p. .  	 
\bibitem[Design of Concrete Structures for Buildings ()]{b13}\label{b13} 	 		\textit{Design of Concrete Structures for Buildings},  		1984. Toronto, ON, Canada: Canadian Standards Association (CSA). CSA.  	 	 (CAN3-A23.3-M84) 
\bibitem[Muttoni et al. ()]{b5}\label{b5} 	 		\textit{Design of Concrete Structures with Stress Fields, 1st ed},  		 			A Muttoni 		,  		 			J Schwartz 		,  		 			B Thürlimann 		.  		1996. Basel: Birkhäuser Verlag.  	 
\bibitem[Kwak and Noh ()]{b19}\label{b19} 	 		‘Determination of strut-and-tie models using evolutionary structural optimization’.  		 			H G Kwak 		,  		 			S H Noh 		.  	 	 		\textit{Engineering Structures}  		2006. 28  (10)  p. .  	 
\bibitem[Deutsche Norm, Concrete, Reinforced and Prestressed Concrete Structures-Part 1:3 Design and Construction, Corrigenda to DIN]{b12}\label{b12} 	 		\textit{Deutsche Norm, Concrete, Reinforced and Prestressed Concrete Structures-Part 1:3 Design and Construction, Corrigenda to DIN},  		 1045-1:2001- 07.  		 	 
\bibitem[different material properties in tension and compression Struct Multidisc Optim ()]{b45}\label{b45} 	 		‘different material properties in tension and compression’.  		 \xref{http://dx.doi.org/10.1007/s00158-011-0633-z}{10.1007/s00158-011-0633-z}.  		 \url{https://doi.org/10.1007/s00158-011-0633-z}  	 	 		\textit{Struct Multidisc Optim}  		2011. 44 p. .  	 
\bibitem[Tuchscherer et al. ()]{b7}\label{b7} 	 		‘Evaluation of existing strut-And-tie methods and recommended improvements’.  		 			Robin G Tuchscherer 		,  		 			David B Birrcher 		,  		 			Christopher S Williams 		,  		 			Dean J Deschenes 		,  		 			Oguzhan Bayrak 		.  	 	 		\textit{ACI structural journal}  		2014. 111  (6)  p. .  	 
\bibitem[Chen et al. ()]{b23}\label{b23} 	 		‘Evaluation of strut-and-tie model applied to deep beam with opening’.  		 			B S Chen 		,  		 			M J Hagenberger 		,  		 			J E Breen 		.  	 	 		\textit{ACI Structural Journal}  		2002. 99  (2)  p. .  	 
\bibitem[Chen et al. ()]{b8}\label{b8} 	 		‘Evaluation of Strut-and-Tie Modeling Applied to Dapped Beam with Opening’.  		 			B S Chen 		,  		 			M J Hagenberger 		,  		 			J E Breen 		.  	 	 		\textit{ACI Structural Journal}  		2002. 99 p. .  	 
\bibitem[Huang et al. ()]{b35}\label{b35} 	 		‘Evolutionary topological optimization of vibrating continuum structures for natural frequencies’.  		 			X Huang 		,  		 			Z H Zuo 		,  		 			Y M Xie 		.  	 	 		\textit{Computer Structure}  		2010. 88  (5-6)  p. .  	 
\bibitem[Maxwell and Breen ()]{b9}\label{b9} 	 		‘Experimental Evaluation of Strut-and-Tie Model Applied to Deep Beam with Opening’.  		 			B S Maxwell 		,  		 			J E Breen 		.  	 	 		\textit{ACI Structural Journal}  		2000. 97 p. .  	 
\bibitem[Liang et al. ()]{b21}\label{b21} 	 		‘Generating Optimal Strut-and-Tie Models in Prestressed Concrete Beams by Performance-Based Optimization’.  		 			Q Q Liang 		,  		 			Y M Xie 		,  		 			G P Steven 		.  	 	 		\textit{ACI Structural Journal}  		2001. 98  (2)  p. .  	 
\bibitem[Bruggi ()]{b15}\label{b15} 	 		‘Generating strut-and-tie patterns for reinforced concrete structures using topology optimization’.  		 			M Bruggi 		.  	 	 		\textit{Computers and Structures}  		2009. 87 p. .  	 
\bibitem[Victoria et al.]{b43}\label{b43} 	 		\textit{Generation of strut-and-tie models by topology design using},  		 			M Victoria 		,  		 			O M Querin 		,  		 			P Martí 		.  		 	 
\bibitem[Coelho and Rodrigues ()]{b40}\label{b40} 	 		‘Hierarchical topology optimization addressing material design constraints and application to sandwich-type structures’.  		 			P G Coelho 		,  		 			H C Rodrigues 		.  	 	 		\textit{Struct Multidiscip Optim}  		2015. 52  (1)  p. .  	 
\bibitem[Van Dijk et al. ()]{b32}\label{b32} 	 		‘Level-set methods for structural topology optimization: a review’.  		 			N P Van Dijk 		,  		 			K Maute 		,  		 			M Langelaar 		,  		 			F Van Keulen 		.  	 	 		\textit{Struct Multidiscip Optim}  		2013. 48 p. .  	 
\bibitem[Wang and Kang ()]{b33}\label{b33} 	 		\textit{MATLAB implementations of velocity field level set method for topology optimization: an 80-line code for 2D and a 100-line code for 3D problems. Structural and Multidisciplinary Optimization},  		 			Y Wang 		,  		 			Z Kang 		.  		2021. 64 p. .  	 
\bibitem[Svanberg ()]{b44}\label{b44} 	 		‘Method of moving asymptotes -a new method for structural optimization’.  		 			K Svanberg 		.  	 	 		\textit{Int. J. Num. Meth. Eng}  		1987. 24 p. .  	 
\bibitem[Marti ()]{b3}\label{b3} 	 		\textit{On plastic analysis of reinforced concrete. Institute of Structural Engineers},  		 			P Marti 		.  		1980. 104.  	 
\bibitem[Xia et al. ()]{b14}\label{b14} 	 		‘Optimization-based three-dimensional strut-and-tie model generation for reinforced concrete’.  		 			Yi Xia 		,  		 			; Matthijs Langelaar 		,  		 			; Max 		,  		 			A N Hendriks 		.  	 	 		\textit{Computer-Aided Civil and Infrastructure Engineering}  		2021. 36 p. .  	 
\bibitem[Liang et al. ()]{b22}\label{b22} 	 		‘Performance-Based Optimization for Strut-Tie Modeling of Structural Concrete’.  		 			Q Q Liang 		,  		 			B Uy 		,  		 			G P Steven 		.  	 	 		\textit{Journal of Structural Engineering}  		2002. 128  (5)  p. .  	 
\bibitem[Recent Approaches to Shear Design of Structural Concrete Journal of Structural Engineering (1998)]{b10}\label{b10} 	 		‘Recent Approaches to Shear Design of Structural Concrete’.  		 ASCE-ACI Committee 445.  	 	 		\textit{Journal of Structural Engineering}  		Dec. 1998. 124  (12)  p. .  	 
\bibitem[Simonetti et al.]{b47}\label{b47} 	 		‘Reliability-Based Topology Optimization: An Extension of the SESO and SERA Methods for Three-Dimensional Structures’.  		 			H L Simonetti 		,  		 			V S Almeida 		,  		 			F De Assis Das Neves 		,  		 			V Del Duca Almeida 		,  		 			L De Oliveira Neto 		.  		 \xref{http://dx.doi.org/10.3390/app12094220.©2023GlobalJournals}{10.3390/app12094220.©2023GlobalJournals}.  		 \url{https://doi.org/10.3390/app12094220.©2023GlobalJournals}  	 	 		\textit{Appl. Sci. 2022}  		12 p. 4220.  	 
\bibitem[Geevar and Menon ()]{b6}\label{b6} 	 		‘Strength of Reinforced Concrete Pier Caps-Experimental Validation of Strut-and-Tie Method’.  		 			Indu Geevar 		,  		 			; Devdas Menon 		.  	 	 		\textit{ACI Structural Journal}  		2019. 116  (1)  p. .  	 
\bibitem[Schlaich and Anagnostou ()]{b2}\label{b2} 	 		‘Stress fields for nodes of strut-and-tie models’.  		 			M Schlaich 		,  		 			G Anagnostou 		.  	 	 		\textit{Struct Eng}  		1990. 116  (1)  p. .  	 
\bibitem[Sethian and Wiegmann ()]{b28}\label{b28} 	 		‘Structural boundary design via level set and immersed interface methods’.  		 			J A Sethian 		,  		 			A Wiegmann 		.  		 \xref{http://dx.doi.org/10.1006/jcph.2000.6581}{10.1006/jcph.2000.6581}.  		 \url{http://dx.doi.org/10.1006/jcph.2000.6581}  	 	 		\textit{Journal of Computational Physics}  		2000. 163 p. .  	 
\bibitem[Kong and Sharp ()]{b46}\label{b46} 	 		‘Structural idealization for deep beams with web openings’.  		 			F K Kong 		,  		 			G R Sharp 		.  	 	 		\textit{Magazine Concrete Res}  		1977. 29 p. .  	 
\bibitem[Allaire et al. ()]{b38}\label{b38} 	 		‘Structural optimization by the level-set method’.  		 			G Allaire 		,  		 			F Jouve 		,  		 			A-M Toader 		.  	 	 		\textit{Free boundary problems},  		 			P Colli 		,  		 			C Verdi 		,  		 			A Visintin 		 (eds.)  		2003. Birkhäuser Basel. p. .  	 
\bibitem[Allaire et al. ()]{b30}\label{b30} 	 		‘Structural optimization using sensitivity analysis and a level-set method’.  		 			G Allaire 		,  		 			F Jouve 		,  		 			A M Toader 		.  	 	 		\textit{J Comput Phys}  		2004. 194  (1)  p. .  	 
\bibitem[Leu et al. ()]{b26}\label{b26} 	 		‘Strut-and-tie design methodology for threedimensional reinforced concrete structures’.  		 			L J Leu 		,  		 			C W Huang 		,  		 			C S Chen 		,  		 			Y P Liao 		.  	 	 		\textit{Journal of structural Engineering}  		2006. New York. 132  (6)  p. 929.  	 
\bibitem[Strömberg ()]{b36}\label{b36} 	 		‘Topology optimization of structures with manufacturing and unilateral contact constraints by minimizing an adjustable compliance-volume product’.  		 			N Strömberg 		.  	 	 		\textit{Struct Multidiscip Optim}  		2010. 42  (3)  p. .  	 
\bibitem[Liang et al. ()]{b20}\label{b20} 	 		‘Topology Optimization of Strut-and-Tie Models in Reinforced Concrete Structures Using an Evolutionary Procedure’.  		 			Q Q Liang 		,  		 			Y M Xie 		,  		 			G P Steven 		.  	 	 		\textit{ACI Structural Journal}  		2000. 97  (2)  p. .  	 
\bibitem[Schlaich et al. (1987)]{b0}\label{b0} 	 		‘Toward a consistent design of structural concrete’.  		 			J Schlaich 		,  		 			K Schafer 		,  		 			M Jennewein 		.  	 	 		\textit{PCI-Journal}  		May/June, 1987. 32  (3)  p. .  	 
\bibitem[Marti ()]{b4}\label{b4} 	 		‘Truss models in detailing’.  		 			P Marti 		.  	 	 		\textit{Concr. Int}  		1985. 7  (12)  p. .  	 
\bibitem[Wang et al. ()]{b34}\label{b34} 	 		‘Velocity field level-set method for topological shape optimization using freely distributed design variables’.  		 			Y Wang 		,  		 			Z Kang 		,  		 			P Liu 		.  	 	 		\textit{Int J Numer Methods Eng}  		2019. 120 p. .  	 
\end{bibitemlist}
 			 		 	 
\end{document}
