[pt] Este trabalho tem por objetivo apresentar um sistema integrado para otimização de forma de estruturas planas que tenham comportamento elasto-plástico. A metodologia implementada propõe uma alternativa à forma conservadora com que tradicionalmente as estruturas têm sido otimizadas, ou seja, admitindo-se que as mesmas possuam comportamento linear elástico. O sistema computacional é denominado integrado pois reúne diversos módulos distintos para o tratamento do problema, como modelagem geométrica, geração de malhas de elementos finitos, análise não-linear da resposta da estrutura, análise de sensibilidade,programação matemática e otimização de estruturas. A geometria do contorno da estrutura plana é definida por meio de curvas (paramétricas)B-splines cúbicas. Estas, por sua vez, são determinadas em função de um conjunto de pontos de interpolação (pontos-chave) e condições de contorno em seus vértices extremos. A correta definição da geometria da estrutura é responsável pelo sucesso do processo de otimização. A resposta da estrutura às solicitações do carregamento externo é avaliada pelo método dos elementos finitos. Para isso, é necessário que o domínio da estrutura seja discretizado. No presente trabalho foi empregado um gerador automático de malhas não estruturadas de elementos finitos isoparamétricos. A configuração de equilíbrio da estrutura é obtida através de um procedimento iterativo/incremental envolvendo o método de Newton-Raphson. Localmente, o equilíbrio é satisfeito pela aplicação de um algoritmo implícito de integração de tensões nos pontos que violarem o critério de plastificação do material. A matriz tangente de rigidez é atualizada a cada iteração da análise e é obtida de forma consistente com o algoritmo de integração das tensões, preservando as características de convergência quadrática assintótica inerentes ao método de Newton- Raphson. No procedimento iterativo de otimização é empregado um algoritmo de programac¸ ão quadrática recursiva que requer a avaliação dos gradientes da função-objetivo e restrições. Para tal, foi implementado um método semi-analítico para a determinação das sensibilidades da resposta estrutural envolvidas nas expressôes dos gradientes citados. O método leva em consideração os efeitos da plastificação ocorrida durante o carregamento da estrutura e é dito -exato- por apresentar imprecisões apenas nos casos em que a magnitude da perturbação da variável é muito pequena, não podendo ser representada corretamente pelo hardware. Os exemplos analisados mostram que a consideração do comportamento elastoplástico da estrutura na otimização de sua forma leva a configurações mais eficientes do que aquelas obtidas admitindo-se a relação linear elástica entre deformações e tensões. / [en] The main goal of this work is to present an integrated system for the optimization of plane structures with elastoplastic behavior. The methodology proposes an alternative for the conservative way in which structures traditionally have been optimized, i.e., that they present linear elastic behavior. The computational system is said to be integrated because it congregates distinct modules for the solution of the problem, such as geometric modelling, finite element mesh generation, non-linear structural response analysis, sensitivity analysis, mathematical programming and optimization of structures. The geometry of the plane structure`s boundary is defined by cubic (parametric) B-splines curves. Those, in turn, are determined by a set of interpolation points (key points) and boundary constraints at their ends. The correct definition of the structure`s geometry is responsible for the success of the optimization process.The structural response to the applied loading is evaluated by the finite element method. For that, the domain of the structure must be discretized. In the present work, an automatic unstructured mesh generator of isoparametric finite elements has been used. The equilibrium layout of the structure is obtained by an iterative/incremental procedure using the standard Newton-Raphson method. Locally, the equilibrium is satisfied by applying an implicit stress return mapping algorithm at points which violate the yield criterion of the material. The tangent stiffness matrix is updated at each analysis iteration and it is obtained in a way which is consistent with the return mapping algorithm, so that the asymptotic quadratic rate of convergence of the Newton-Raphson method is preserved. The use of a quadratic recursive programming algorithm in the optimization procedure involves the gradient evaluation of the objective function and constraints. For that, a semi-analytical method for the calculation of the response sensitivities, which appear in the gradient expressions, has been implemented. The technique takes into account the plastic effects which take place during the loading of the structure and is considered - exact- up to round-off errors, which occurs when the magnitude of the perturbation is so small that the hardware cannot accurately represent it.The examples presented demonstrate that the consideration of the elastoplastic behavior of the material during the optimization process leads to structural layouts which are more efficient than of those obtained under the assumption of linear elastic relationship between strains and stresses. / [es] Este trabajo tiene por objetivo presentar un sistema integrado para otimización de forma de extructuras planas que tengan comportamiento elástico-plástico. LA metodología implementada propone una alternativa a la forma conservadora con que tradicionalmente las extructuras han sido optimizadas, o sea, admitiendo que las poseen um comportamiento lineal-elástico. EL sistema computacional se denomina integrado pues reúne diversos módulos para el tratamiento del problema, como modelage geométrica, generación de mallas de elementos finitos, análisis no lineal de la respuesta de la extructura, análisis de sensibilidad,programación matemática y otimización de extructuras. LA geometría del contorno de la extructura plana es definida por medio de curvas (paramétricas)B splines cúbicas. Estas, por su vez, son determinadas en función de un conjunto de puntos de interpolación (puntos claves) y condiciones de contorno en sus vértices extremos. La definición correta de la geometría de la extructura es responsable por el éxito del proceso de otimización. La respuesta de la extructura a las solicitudes de carga externa se evalúa por el método de los elementos finitos. Para esto, es necesario que el dominio de la extructura sea discretizado. En este trabajo se utiliza un generador automático de mallas no extructuradas de elementos finitos isoparamétricos. La configuración de equilíbrio de la extructura se obtiene a través de un procedimiento iterativo/incremental que envuelve el método de Newton Raphson. Localmente, el equilíbrio es satisfecho por la aplicación de un algoritmo implícito de integración de tensiones en los puntos que violen el critério de plastificación del material. La matriz tangente de rigidez se actualiza a cada iteración del análisis y se obtiene de forma consistente con el algoritmo de integración de las tensiones, preservando las características de convergencia cuadrática asintótica inherentes al método de Newton Raphson. En el procedimiento iterativo de otimización se utiliza un algoritmo de programación cuadrática recursiva que requiere la evaluación de los gradientes de la función objetivo y restricciones. Para tal, se implementó un método semi analítico para la determinación de las sensibilidades de la respuesta extructural envolvidas en las expresóes de los gradientes citados. EL método lleva en consideración el hecho de que la plastificación que ocurre durante la carga de la extructura y se dice exacta por presentar imprecisiones apenas en los casos en que la magnitud de la perturbación de la variable es muy pequeña, no puede ser representada correctamente por el hardware. Los ejemplos analizados muestran que la consideración del comportamiento elástico-plástico de la extructura en la otimización de su forma lleva la configuraciones más eficientes de que aquellas obtenidas admitiendo la relación lineal elástica entre deformaciones y tensiones.

[pt] Este trabalho consiste no desenvolvimento e implementação de um sistema computacional integrado para Modelagem Geo- étrica, Geração de Malhas, Análise Numérica, Auto- Adaptação do tipo h e Otimização de Forma e Espessura em Cascas. O Modelo Geométrico adotado consiste em representar a superfície por meio de B-Splines Racionais Recortadas, admitindo variação de espessura segundo as mesmas funções que descrevem a superfície. Para a utilização nos módulos de Auto-Adaptação e Otimização, Geradores de Malhas Não-Estruturadas sobre superfícies paramétricas foram empregados. Com base em um gerador de malhas triangulares que utiliza a técnica de avanço de Fronteira, dois geradores de malhas quadrilaterais foram desenvolvidos. Os elementos finitos empregados são baseados nas hipóteses de Reissner-Mindlin e no conceito de degeneração de sólidos. São considerados os elementos tradicionais, baseados puramente em deslocamentos, bem como elementos de formulação mais robusta, com base em campos assumidos de deformação. Um procedimento para Auto-Adaptação de Malhas do tipo foi desenvolvido para o Modelo Geométrico e para os elementos considerados. As malhas obtidas a partir de níveis de erro prescritos permitem aos modelos a obtenção de suas taxas de convergência mesmo em problemas com singularidades e efeitos de fronteira. A Análise de Sensibilidade, ferramenta essencial nos procedimentos de Otimização, é feita com a utilização do Método Semi-Analítico, considerando as características especiais dos elementos de casca. Foi desenvolvida para elementos degenerados de casca uma versão do Método Semi- Analítico que elimina sua imprecisão eventual, através da diferenciação exata das componentes de deslocamento de corpo rígido. Para os elementos baseados puramente em deslocamentos, o Método Analítico também foi desenvolvido. Os módulos de Otimização de Forma e Espessura trabalham sobre diversas possibilidades de definição das variáveis de projeto e com diferentes algoritmos de Programação Matemá tica. Permitem também a Otimização de Forma com consideração de Auto-Adaptatividade para obter as malhas durante o processo de mudança de forma. Devido à interconexão entre os diversos módulos e com o objetivo de facilitar a comunicação e a implementação dos mesmos, o sistema computacional foi completamente desenvolvido utilizando técnicas de Programação Orientada para Objetos. / [en] The purpose of this thesis is the development of an integrated computational system for Geometric Modelling, Unstructured Mesh Generation, Numerical Analysis, Adaptivity and Shape Optimization of Shells.The Geometric Model is composed of Non-Uniform Rational B-Spline Surfaces (NURBS), further modified by trimming loops described in the parametric plane. Smooth thickness variations may be modelled by means of the same functions that describe the surface geometry.For Unstructured Mesh Generation, two algorithms were implemented for quadrilateral elements and one for triangles. The triangular Mesh Generator is based on the Advancing Front Technique applied to parametric surfaces. The quadrilateral Mesh Generators employ the indirect approach for converting the triangular meshes to pure quadrilateral ones. The Finite Element formulation is based on the degenerated isoparametric approach. Pure displacement based elements, as well as assumed strain robust shell elements, are employed in the Analysis, Adaptivity and Optimization modules. A procedure for h-Adaptive Mesh Refinement was developed for the shell models. For this purpose an Error Estimator scheme, based upon a Zienkiewicz-Zhu Patch Recovery Technique, adapted to handle curved shell surfaces, was employed. The adaptive procedure allows the convergence rates of the Finite Element Model to be maintained even in the presence of singularities and boundary layers.For the Sensitivity Evaluation, the well- known Semi-Analytical Method is employed and adapted for the degenerated solid shell element approach. In order to solve the inaccuracy problems inherent to the application of the method for certain types of structures, the recent Refined Semi- Analytical Method, is extended for degenerated shell elements. For the pure displacement-based elements, the Analytical Method is also developed. The Shape and Thickness Optimization modules work with a wide variety of design variable descriptions, different mathematical Programming algorithms, Sensitivity schemes and Finite Element Models. The possibility of h-Adaptive Mesh refinement in conjunction with Shape Optimization is also considered in this stage. In order to ease up code expansion, communication and data exchange between the modules,the computational system was fully developed employing Object-Oriented Programming techniques.

[pt] Este trabalho propõe uma metodologia para a otimização de forma de estruturas geometricamente não-lineares. O objetivo desta metodologia é evitar os problemas de instabilidade apresentados por estruturas otimizadas de acordo com a formulação clássica. Ela foi implementada para problemas bidimensionais e os resultados obtidos na otimização de diferentes estruturas demonstraram o seu sucesso. Utilizando-se conceitos de modelagem geométrica, a forma da estrutura é defini-da através das curvas de seu contorno. Assim, a representação paramétrica de curvas e a definição destas em função de um conjunto de pontos de interpolação (pontos-chave) são discutidas detalhadamente. A ênfase é dada à interpolação através de B-splines,devido a sua grande flexibilidade. O problema de otimização é definido com base no modelo geométrico e as variáveis de projeto são as coordenadas dos pontos-chave. A simetria da estrutura é garantida através da ligação de variáveis. A estrutura é analisada através de elementos isoparametricos planos. Assim, antes de realizar a análise, é necessário discretizar a estrutura em um conjunto de elementos finitos. Para realizar esta tarefa foram implementados diferentes algoritmos de geração de malhas, tanto estruturadas quanto não-estruturadas. O método de Newton-Raphson é utilizado pa- ra determinar a configuração de equilíbrio e diferentes métodos podem ser aplicados para determinar os pontos críticos. Devido aos problemas de convergência apresentados pelos métodos diretos para a determinação dos pontos crticos, um método semi-direto foi desenvolvido neste trabalho. Os resultados obtidos na análise de diferentes exemplos mostraram a adequação dos elementos finitos e dos métodos numéricos implementados. Os algoritmos de programação matemática utilizados neste trabalho precisam dos gradientes da função objetivo e das restrições, que são calculadas com base nos gradientes das respostas da estrutura. Partindo-se de equações gerais válidas para quaisquer elementos,foram desenvolvidas expressões analíticas que permitem o cálculo exato das sensibilidades de elementos finitos isoparamétricos formulados através do procedimento Lagrangiano Total. O desenvolvimento e a implementação de expressões semelhantes para elementos mais complexos é uma tarefa bastante árdua. Por outro lado, o método das diferenças fi- nitas é simples e genérico, mas muito caro computacionalmente. O método semi-analítico mantémm as vantagens da utilização de diferenças finitas e possui um custo computacional baixo, porém pode apresentar sérios problemas de preciso. Devido a estes motivos, foi desenvolvido neste trabalho um procedimento para melhorar a qualidade das sensibilidades semi-analíticas de estruturas geometricamente não-lineares. O procedimento é baseado na diferenciação exata dos movimentos de corpo rígido do elemento utilizado. Os resultados numéricos obtidos demonstraram a sua eficácia. / [en] This work presents a methodology for shape optimization of geometrically nonlinear structures. The main purpose is to avoid the stability problems generated by optimization based on linear behavior. The methodology was implemented for two-dimensional problems, and several structures were successfully optimized. Using geometrical modeling concepts, the shape of the structure is defined by its boundary curves. Therefore, parametric representation and curve definition by a set of key points are discussed in detail. Due to its flexibility in shape definition, particular attention is given to interpolation using B- splines. The optimization problem is defined based on the geometrical model and the design variables are the positions of key points. Design variable linking can be applied to enforce symmetry.The structure it is analyzed using plane isoparametric elements. Thus, is necessary to perform the discretization of the structure in a set of finite elements. Different algorithms were implemented to generate structured or unstructured finite element meshes. The standard Newton-Raphson method is applied to find the equilibrium configuration, and different methods can be used to evaluate critical points. Due to the convergence problems presented by direct methods, a new semi-direct method was developed. The numerical results show the suitability of the finite elements and numerical methods implemented in the present work.The mathematical programming algorithms used in this work require the evaluation of design sensitivities in order to compute the search direction of the optimization process.Using basic sensitivity equations, which are independent from the particular element, analytical expressions were developed for the sensitivity computation of isoparametric elements formulated according to the Total Lagrangian approach. Applying the analytical method for more complex elements is very cumbersome and error prone. On the other hand, the finite difference method is simple and generic, but its computational cost is prohibitive. The semi-analytical method preserves the advantages of the use of finite differences and has a low computational cost, but presents severe accuracy problems. Hence, a method based on the exact differentiation of the rigid body motions was developed in this work to improve the accuracy of the semi- analytical sensitivities of geometrically nonlinear structures. The numerical examples show that this method eliminates the abnormal errors presented by the semi- analytical sensitivities. / [es] Este trabajo propone una metodología para la optimización de forma de extructuras geométricamente no lineares. EL objetivo de esta metodología es evitar los problemas de inestabilidad que presentan las extructuras optimizadas de acuerdo con la formulación clásica. Ella fue implementada para problemas bidimensionales y los resultados obtenidos en la optimización de diferentes extructuras demuestran su éxito. Utilizando conceptos de modelaje geométrica, la forma de la extructura se define a través de las curvas de contorno. Así, la representación paramétrica de curvas y su definición en función de un conjunto de puntos de interpolación (puntos clave) son discutidas detalladamente. Se le da especial a la interpolación a través de B splines,debido a su gran flexibilidad. El problema de optimización se define con base en el modelo geométrico y las variables de proyecto son las coordenadas de los puntos clave. La simetría de la extructura se garante a través de la llamada de variables. La extructura se analiza a través de elementos isoparamétricos planos. Antes de realizar el análisis, es necesario discretizar la extructura en un conjunto de elementos finitos. Para realizar esta tarea fueron implementados diferentes algoritmos de generación de mallas, tanto extructuradas como no extructuradas. EL método de Newton Raphson es utilizado para determinar la configuración de equilibrio y pueden ser aplicados diferentes métodos para determinar los puntos críticos. Debido a los problemas de convergencia presentados por los métodos directos para la determinación de los puntos crticos, se desarrolló un método semidirecto. Los resultados obtenidos en el análisis de diferentes ejemplos muestran la adque los elementos finitos los métodos numéricos implementados son adequados Los algoritmos de programación matemática utilizados en este trabajo necesitan los gradientes de la función objetivo y de las restricciones, que son calculadas con base en los gradientes de las respuestas de la extructura. Partiendo de las ecuaciones generales válidas para cualesquiera elementos, fueron desarrolladas expresiones analíticas que permiten el cálculo exacto de las sensibilidades de elementos finitos isoparamétricos formulados a través del procedimiento Lagrangiano Total. EL desarrollo y la implementación de expresiones semejantes para elementos más complexos es una tarea bastante árdua. Por otro lado, el método de las diferencias finitas es simple y genérico, pero muy caro computacionalmente. EL método semianalítico mantiene las ventajas de la utilización de diferencias finitas y posee un costo computacional bajo, pero puede presentar serios problemas de precisión. Es por ello que se desarrolló en este trabajo un procedimiento para mejorar la calidad de las sensibilidades semianalíticas de extructuras geométricamente no lineares. EL procedimiento tiene como base la diferenciación exacta de los movimentos de cuerpo rígido del elemento utilizado. Los resultados numéricos obtenidos demuestran su eficacia.

[pt] Esse trabalho propõe uma implementação numérica para otimização de forma em problemas bi-dimensionais de elasticidade. O objetivo principal é propor uma metodologia eficiente e robusta para solução de problemas de otimização de forma considerando a minimização de concentração de tensões. Na implementação proposta, a análise estrutural é realizada pelo Método dos Elementos de Contorno Simétrico de Galerkin (MECSG), evitando-se assim a dispendiosa etapa de geração da malha. A avaliação das tensões no contorno é obtida por meio de um método preciso, ideal para problemas com concentrações de tensões. Outro aspecto relevante na implementação é a adequada partição das equações do MECSG de forma a reduzir, consideravelmente, o esforço computacional associado à etapa da análise estrutural. O problema de otimização é resolvido utilizando-se um método de otimização moderno, conhecido como Programação Cônica de Segunda Orderm (PCSO). Especificamente, busca-se a reposta do problema de otimização não linear por meio da solução de uma sequência de subproblemas de PCSO. / [en] In this work a numerical implementation of shape optimization in two-dimensional linear elasticity problems is proposed. The main goal is to propose a robust and efficient methodology for the solution of shape optimization problems regarding the minimization of stress concentration effects. In the proposed implementation, the structural analysis is performed by the Symmetric Galerkin Boundary Element Method (SGBEM), thus disposing of the mesh generation burden. The boundary stress evaluation is carried out by an accurate approach which is ideally suited for problems with stress concentrations. Another relevant feature of the proposed implementation is a suitable partition of the SGBEM equations which aims at reducing the computational effort associated with the structural analysis stage. The solution for the optimization problem is obtained by means of a modern numerical optimization method, the so-called Second Order Conic Programming (SOCP). Specifically, the solution for the non-linear optimization is sought by solving a sequence of SOCP subproblems.

[pt] Este trabalho apresenta um estudo do cálculo de sensibilidades necessário para a análise de problemas inversos e de otimização, usando o método híbrido dos elementos de contorno. Com esta finalidade, é desenvolvida uma formulação que permite obter as sensibilidades à mudança de forma, por diferenciação implícita das integrais de contorno, de uma estrutura já discretizada. Demonstra-se que as sensibilidades das matrizes obtidas desta formulação apresentam propriedades espectrais definidas, que são derivadas da formulação básica do método híbrido dos elementos de contorno. Todo o desenvolvimento é feito para um problema da elastostática tridimensional, embora sejam apresentadas apenas aplicações de problemas bidimensionais e de potencial, como casos particulares. As singularidades que surgem na integração no cálculo das sensibilidades são facilmente solucionáveis a partir das integrais da formulação básica do método híbrido dos elementos de contorno. As implementações numéricas são feitas utilizando a linguagem de programação Maple V release 3. Para ambos os casos, de potencial e elasticidade bidimensional, são usados elementos lineares para a representação do contorno. São apresentadas comparações entre os resultados analíticos obtidos através desta formulação com os resultados obtidos usando a técnica de diferenças finitas (centradas), com o objetivo de demonstrar a eficiência e precisão da metodologia aqui desenvolvida. / [en] The present work describes a formulation for computing design sensitivities required in inverse problems and shape optimization of solid objects, in the frame of the hybrid boundary element method. The so-called direct differentiation method is applied in order to calculate the gradients, i.e. the implicit diferentiation of the discretized boundary is performed, resulting in a general and efficient analysis technique for shape design sensitivity analysis of all structural quantities. It is demonstrated that the resulting sensitivities matrices present some useful spectral properties, which are related to the matrix spectral properties of the basic hybrid formulation. This formulation is valid for tridimensional solids, although only potential and bidimensional applications are considered as particular cases. The singularities that appear in the resulting boundary integrals are exactly the same which have already been dealt with in the basic formulation. The analytical and numerical procedures were performed by using the mathematical package Maple V release 3. Linear boundary elements were used for both potential and elasticity problems. Numerical results obtained by the present procedure are compared to finite differences results to demonstrate the effectiveness of the present formulation. / [es] Este trabajo presenta un estudio del cálculo de sensibilidades, que tiene gran importancia en el análisis de problemas inversos y de optimización, usando el método híbrido de los elementos de contorno. Con esta finalidad, se desarrolla una formulación que permite obtener las sensibilidades al cambio de forma de una extructura ya discretizada, por diferenciación implícita de las integrales de contorno. Se demuestra que las sensibilidades de las matrices obtenidas por esta formulación presentan propriedades espectrales definidas, que son derivadas de la formulación básica del método híbrido de los elementos de contorno. El desarrollo de la formulación se realiza para un problema de elastostática tridimensional, aunque se presentan apenas las aplicaciones de problemas bidimensionales y de potencial, como casos particulares. Las singularidades que surgen en la integración en el cálculo de las sensibilidades pueden ser fácilmente resueltas a partir de las integrales de la formulación básica del método híbrido de los elementos de contorno. La implementación numérica utiliza el lenguaje de programación Maple V release 3. Para los casos de potencial y elasticidad bidimensional, se utilizan elementos lineales para la representación del contorno. Se comparan los resultados analíticos obtenidos a través de esta formulación con los resultados obtenidos usando la técnica de diferencias finitas (centradas), con el objetivo de demostrar la eficiencia y precisión de la metodología aqui desarrollada.

[pt] Estruturas treliçadas espaciais são arranjos geométricos de barras amplamente utilizados em coberturas de edificações. Diversos fatores favorecem o seu uso, tais como a capacidade de vencer grandes vãos e a facilidade em assumir diversas formas. A busca pela geometria ótima é um objetivo importante no projeto de estruturas, onde o interesse principal é minimizar o custo da estrutura. O objetivo deste trabalho é apresentar um sistema computacional capaz de minimizar o peso de estruturas treliçadas cuja geometria é definida por curvas de Bézier. Portanto, os pontos de controle das curvas de Bézier são utilizados como variáveis de projeto. As áreas das seções transversais das barras e a altura da treliça também são consideradas como variáveis de projeto e restrições sobre a tensão de escoamento e a tensão crítica de Euler são impostas no problema de otimização. A estrutura é analisada por meio do método dos elementos finitos considerando a hipótese do comportamento linear físico e geométrico. Os algoritmos de otimização usados neste trabalho utilizam o gradiente da função objetivo e das restrições em relação às variáveis de projeto. O sistema computacional desenvolvido neste trabalho foi escrito em linguagem MATLAB e conta com uma integração com o SAP2000 por meio da OAPI (Open Application Programming Interface). Os resultados numéricos obtidos demonstram a eficiência e a aplicabilidade deste sistema. / [en] Spatial truss structures are geometrical arrangements of bars widely used in building roofs. Several factors favor their use, such as the ability to overcome large spans and the capability of assuming a variety of configurations. The search for optimal geometry is an important goal in the design of structures, where the main interest is to minimize the cost of the structure. The objective of this work is to present a computational system capable of minimizing the weight of truss structures whose geometry is defined by Bézier curves. Therefore, the control points of the Bézier curves are used as design variables. The cross-sectional areas of the bars and the truss height are also considered as design variables and constraints on the yield stress and Euler critical stress are imposed on the optimization problem. The structure is analyzed using truss elements considering the physical and geometric linear behavior. The optimization algorithms used in this work require the gradient of the objective function and constraints with respect to the design variables. The computational system developed in this work was written in MATLAB and has an integration with SAP2000 through the OAPI (Open Application Programming Interface). The obtained numerical results demonstrate the efficiency and applicability of the developed system.