Feasible span-to-depth ratios for many modern bridge systems have been identified and documented in literature. No such parameters have been adequately identified in terms of proportioning bridge piers. This thesis includes a study of 22 existing reinforced concrete highway bridges and their respective pier systems to determine the state-of-the-art in design. The effect of different geometric and material parameters such as concrete strength, reinforcement ratio and slenderness ratio on the structural behavior of individual piers and multiple pier systems was examined. Approximate methods, which may be used for the purposes of preliminary design are discussed and reviewed. Serviceability and ultimate limit states design aids that can be used to identify appropriate preliminary cross-sectional pier dimensions and reinforcement ratios for individual piers given various slenderness ratios were developed. The structural behavior as well as an approach to the preliminary design of multiple pier bridge systems is presented.

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第18585号 / 工博第3946号 / 新制||工||1606(附属図書館) / 31485 / 京都大学大学院工学研究科機械理工学専攻 / (主査)教授 西脇 眞二, 教授 田畑 修, 教授 松原 厚 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

[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.

This study examines the use of embedded shape memory alloy (SMA)actuators for adaptive control of the themomechanical response of composite structures. Control of static and dynamic responses are demonstrated including thermal buckling, thermal post-buckling, vibration, sonic fatigue, and acoustic transmission. A thermomechanical model is presented for analyzing such shape memory alloy hybrid composite (SMAHC) structures exposed to thermal and mechanical loads. Also presented are (1) fabrication procedures for SMAHC specimens, (2) characterization of the constituent materials for model quantification, (3) development of the test apparatus for conducting static and dynamic experiments on specimens with and without SMA, (4) discussion of the experimental results, and (5) validation of the analytical and numerical tools developed in the study. The constitutive model developed to describe the mechanics of a SMAHC lamina captures the material nonlinearity with temperature of the SMA and matrix material if necessary. It is in a form that is amenable to commercial finite element (FE) code implementation. The model is valid for constrained, restrained, or free recovery configurations with appropriate measurements of fundamental engineering properties. This constitutive model is used along with classical lamination theory and the FE method to formulate the equations of motion for panel-type structures subjected to steady-state thermal and dynamic mechanical loads. Mechanical loads that are considered include acoustic pressure, inertial (base acceleration), and concentrated forces. Four solution types are developed from the governing equations including thermal buckling, thermal post-buckling, dynamic response, and acoustic transmission/radiation. These solution procedures are compared with closed-form and/or other known solutions to benchmark the numerical tools developed in this study. Practical solutions for overcoming fabrication issues and obtaining repeatable specimens are demonstrated. Results from characterization of the SMA constituent are highlighted with regard to their impact on thermomechanical modeling. Results from static and dynamic tests on a SMAHC beam specimen are presented, which demonstrate the enormous control authority of the SMA actuators. Excellent agreement is achieved between the predicted and measured responses including thermal buckling, thermal post-buckling, and dynamic response due to inertial loading. The validated model and thermomechanical analysis tools are used to demonstrate a variety of static and dynamic response behaviors associated with SMAHC structures. Topics of discussion include the fundamental mechanics of SMAHC structures, control of static (thermal buckling and post-buckling) and dynamic responses (vibration, sonic fatigue, and acoustic transmission), and SMAHC design considerations for these applications. The dynamic response performance of a SMAHC panel specimen is compared to conventional response abatement approaches. SMAHCs are shown to have significant advantages for vibration, sonic fatigue, and noise control. / Ph. D.

A análise de compósitos laminados apresenta grandes desafios, pois, diferentemente dos materiais isotrópicos homogêneos, os compósitos laminados são constituídos de materiais heterogêneos e anisotrópicos. Além disso, as distribuições de tensões interlaminares obtidas com as formulações convencionais são descontínuas e imprecisas. Sua melhoria, portanto, é imprescindível para buscar e modelar critérios de falha relacionados às estruturas formadas por compósitos laminados. Diante disso, este trabalho se concentrou no desenvolvimento e implementação computacional de um elemento finito posicional de pórtico plano laminado cuja cinemática é descrita ao longo da espessura do laminado de acordo com a teoria Layerwise. A formulação do elemento considera a não linearidade geométrica, originada pela ocorrência de grandes deslocamentos e rotações, e admite deformações moderadas, em função da lei constitutiva de Saint-Venant-Kirchhoff. O desenvolvimento deste trabalho se iniciou com uma preparação teórica sobre mecânica dos sólidos deformáveis e métodos numéricos para que fossem adquiridos os subsídios teóricos necessários ao desenvolvimento de códigos computacionais, à interpretação dos resultados e à tomada de decisões quando das análises numéricas. A formulação desenvolvida é Lagrangiana total com emprego do método dos elementos finitos baseado em posições. Inicialmente o elemento finito posicional de pórtico plano homogêneo é proposto, uma vez que sua cinemática possibilita uma expansão natural para o caso laminado. Os graus de liberdade são compostos por posições nodais e por vetores generalizados que representam o giro e a variação na altura da seção transversal. A eficiência do elemento é constatada através de análises realizadas em problemas de pórtico sujeitos a grandes deslocamentos e rotações. Os resultados obtidos apresentaram excelente concordância com soluções numéricas e analíticas disponíveis na literatura. Uma expansão natural da cinemática é empregada na formulação do elemento laminado. Os graus de liberdade do elemento são as posições nodais e as componentes de vetores generalizados associados às seções transversais de cada lâmina. Dessa forma, as lâminas têm liberdade para variação de espessura e giro independente das demais, mas com as posições compatibilizadas nas interfaces. Os resultados de análises numéricas realizadas em vários exemplos demonstram a eficiência da formulação proposta, pois as distribuições de deslocamentos e tensões ao longo da espessura do laminado apresentaram excelente concordância com as obtidas a partir de análises numéricas utilizando um elemento finito bidimensional em uma discretização bastante refinada. Os exemplos analisados contemplam problemas com seção laminada fina ou espessa. / The analysis of laminated composites presents challenges because, unlike homogeneous isotropic materials, the laminated composites are made up of heterogeneous and anisotropic materials. Moreover, the distribution of interlaminar stresses obtained with conventional formulations are discontinuous and inaccurate. His improvement is therefore essential to check and modeling failure criteria related to structures formed by laminates. Thus, this work focused on developing and computational implementation of a positional finite element of laminated plane frame whose kinematics is described throughout the thickness of the laminate according to Layerwise theory. The formulation element considers the geometric nonlinearity, caused by the occurrence of large displacements and rotations, and admits moderate deformation, in the constitutive law function of Saint-Venant-Kirchhoff. The development of this work began with a theoretical preparation on mechanics of deformable solids and numerical methods for the acquired of the theoretical support needed for the development of computational codes, interpretation of results and decision-making when of the numerical analyzes. The developed formulation is total Lagrangian with use of the finite element method based on positions. Initially the positional finite element of homogeneous plane frame is proposed, since their kinematic enables a natural expansion for the laminate case. The degrees of freedom are composed of nodal positions and generalized vectors representing the spin and the variation in the height of the cross section. The efficiency of the element is verified through analyzes performed in frame problems subject to large displacements and rotations. The results showed excellent agreement with numerical and analytical solutions available in the literature. A natural expansion of the kinematics is used in the formulation of the laminate element. The degrees of freedom of the element are the nodal positions and components of the generalized vectors associated to cross-sections of each lamina. Thus, the laminas are free for the thickness variation and for independent spin, but with the positions matched in the interfaces. The results of numerical analysis performed in various examples show the effectiveness of the proposed formulation, since the distributions of displacements and stresses through the thickness of the laminate agreed well with those obtained from numerical analysis using a discretization with two-dimensional finite elements in a very refined. The examples discussed include problems with thin or thick laminated section.

No presente trabalho são avaliados os principais parâmetros de projeto na análise da estabilidade global de estruturas em concreto pré-moldado de múltiplos pavimentos, considerando a não-linearidade física (NLF), não-linearidade geométrica (NLG) e ligações semi-rígidas. Os sistemas estruturais analisados são constituídos por pórticos com ligação viga-pilar semi-rígida e pilares engastados na fundação. As ligações viga-pilar têm sua tipologia definida por chumbadores retos e capa de concreto moldado in loco com armadura de continuidade. A NLF é avaliada segundo a construção dos diagramas M x N x 1/r, onde foi considerado de forma consistente o efeito da força normal, da armadura passiva, da armadura ativa e a reologia do concreto. São propostas funções e definidos coeficientes redutores de rigidez que foram comparados com valores apontados em expressões normativas que contemplam de forma aproximada a NLF. A NLG é avaliada com o auxílio do programa ANSYS e na forma aproximada segundo o coeficiente z e o método da carga lateral fictícia (P-). Apresentam-se modelos analíticos de caracterização de rigidez e resistência de ligações viga-pilar ao momento fletor negativo e positivo. Na análise numérica é feito um estudo de caso completo de uma estrutura típica de múltiplos pavimentos em concreto pré-moldado com o auxílio do programa ANSYS, avaliando-se diferentes formas de consideração da NLF e da NLG. É analisada a distribuição de esforços na estrutura frente às combinações de ações utilizadas e os modelos de comportamento das ligações, para algumas variações de geometria e carregamento. Dentre as conclusões, pode ser citado que os coeficientes redutores obtidos segundo o diagrama M x N x 1/r divergem das indicações normativas para consideração simplificada de NLF. Quanto à consideração simplificada de análise NLG, o coeficiente z apresentou resultados melhores na previsão dos esforços de segunda ordem em relação aos obtidos por 0,95.z. / In this research it was evaluated the main design parameters in the analysis of global stability of multi-storey precast concrete framed structures considering the physical nonlinearity (PNL), geometric nonlinearity (GNL) and semi-rigid connections behavior. The structural systems considered are those consisting of frames with beam-to-column semi-rigid connection and columns inset in the foundation. The beam-to-column connections are defined by two dowels and cast-in-place concrete cap with longitudinal reinforced bars. The PNL is evaluated according the M x N x 1/r diagrams, where are considered the effect of normal force, the reinforcing steel, prestressing steel and creep. Functions are proposed and reducing stiffness coefficients are defined, with were compared with codes that include the approximated PNL. The GNL is evaluated according to ANSYS software by no approximate analisys and its approximate according to the coefficient z and (P-) method. This research presents an analytical model characterization of stiffness and strength of beam-to-column connections to the negative and positive bending moment. In the numerical analysis is made a full case study of a typical multi-storey precast concrete framed structure with the aid of the ANSYS software, evaluating different forms of consideration of the PNL and the GNL. It analysed the bending moment distribution according the combinations of actions used and the behavioral models of connections, for some variations at geometry and loading. The reduction coefficients of flexural rigidity obtained according the M x N x 1/r diagrams differ from normative indications for simplified account of PNL. In GNL simplified analysis, the coefficient z showed better results in predicting the second-order effects with respect to those obtained by the coefficient 0,95. z.

Sheet metal forming is used in a wide range of industrial processes ranging from tube manufacturing to automobile and aviation industry. It includes processes like stamping, bending, stretching, drawing and wheeling. In the past few years materials for sheet metal forming and, technology have improved a lot. The improved materials have higher strength and more ductility than conventional sheet steel and therefore they have to be worked differently. For such steels conventional methods can not be applied totally. So there is a need for constant improvement in technology. Trial and error method currently in use increases lead time and is not economic also. To overcome the problems, use of simulation software in metal forming processes has increased significantly. The rapid development of software technology accompanied with lower cost computer hardware have enabled many manufacturing operations to be modeled cost-effectively that only a few years ago would have been considered impractical. However there are some difficulties in simulation of sheet metal forming process. For example it is never an easy task to select the correct software for a particular process. Various authors ascribe different causes for the difficulties. Among them the prominent ones are lacunae in elasto-plastic modeling, material behaviour, involved complexities and a lack of good elements. Apart from that the demands of sheet metal processes are increasing both with respect to the tolerance requirements of the finished part and with regard to geometric complexity of the part to be formed. A few years ago finite elements have been developed using Papcovitch-Neuber solutions of the Navier equation for the displacement function. Among these elements PN5X1 has the abilities to predict both displacements and stresses accurately. And recently the element is extended to include material nonlinearity and is working well for the small deformation range. To use this element for sheet metal forming it is necessary that the element should predict correct results for large deformations. In the present work we have further extended this element for large displacements and large rotation. In the literature there are various algorithms recommended for use with large deformation. Among them we have selected a suitable algorithm and verified its usefulness. First we have taken a simple truss and applied loads to cause large deflection. We observe adequate convergence with the chosen algorithm and then we extend it to PN5X1. in large deformation analysis, equilibrium is computed about the deformed shape. In the chosen algorithm we apply incremental loading and within each load step loop we iterate for equilibrium. We ensure error free solution (equilibrium) before additional loading is introduced. With the help of flowchart these processes have been depicted. A computer program in C, based on the above incremental method and equilibrium check has been written. For the purpose of verification of the program, we have solved some benchmark tests. We start with linear cases and then attempt a number of geometric nonlinear problems like- cantilever subjected to end shear, pinched cylinder with open end etc. We have also included the classical benchmark problem of the cantilever subjected to end moment. The present algorithm gives solutions which are in excellent agreement with those reported in the literature. Finally, we look at some aspects of the problem which require further investigation.

La modélisation numérique des structures composites à renfort fibreux de géométrie complexe constitue un axe de recherche majeur afin de prédire correctement leur comportement mécanique. Dans ce contexte, l’étude menée dans ce travail de thèse porte sur le développement de nouveaux éléments finis basés sur une approche numérique multi-échelle, appelée Approche de la Fibre Projetée (AFP). Cette approche a l’avantage de tenir compte de la présence des fibres au sein d’un espace matrice sans les discrétiser, ce qui limite considérablement la taille du système à résoudre. Pour analyser le comportement des structures composites, plusieurs éléments finis 2D et 3D ont été développés et implémentés dans le code ABAQUS via la routine UEL. Plusieurs cas tests de validation sont considérés pour tester la précision et l’efficacité des éléments finis proposés et les résultats obtenus sont globalement en bon accord avec les solutions de référence. De plus, l’intérêt de la nouvelle approche (AFP) est particulièrement mis en exergue en étudiant des structures composites complexe à renfort 3D : une plaque sandwich cousue et une plaque sandwich à âme creuse renforcée par des fibres en forme de « 8 ». / Numerical modeling of composite materials and structures with complex geometry of fiber reinforcement, such as stitched composites, constitutes a major research axis in order to correctly predict their mechanical behavior. Within this context, this study focuses on the development of new linear and nonlinear specific finite elements based on a multiscale numerical approach, called the Projected Fiber Approach (PFA). This numerical approach has the advantage of taking into account the presence of fi bers, long or short and distributed randomly or specifically, within a matrix space without discretizing them. Consequently, the obtained system of equations size is equivalent to that without reinforcement (matrix), which considerably reduces the computational cost. To analyze the linear and geometrically nonlinear behaviors of composite structures, two membrane finite elements, named PFT3 and PFQ4, and a 3D solid finite element, named PFH8, were developed and implemented into the commercial finite element code ABAQUS via the user element subroutine (UEL). Several numerical linear and nonlinear tests are considered to assess the accuracy and efficiency of the proposed composite finite elements, and the obtained results are globally in good agreement with the reference solutions. Moreover, the major interest of the PFA approach is particularly emphasized by studying two 3D complex reinforced composite structures: a stitched sandwich plate and a hollow core sandwich plate reinforced by "8" shape fibers.

[pt] O objetivo deste trabalho é complementar a ferramenta de dimensionamento de pórticos planos de concreto armado já existente no Ftool, programa educacional amplamente difundido no meio acadêmico. Para tanto, foi introduzido o cálculo e dimensionamento de pilares de concreto armado à flexão composta reta. Foi adicionada uma nova seção transversal, referente aos pilares retangulares com armaduras simétricas. Além disso, foi necessária a inclusão de um método de análise não linear geométrica simplificado que fosse compatível com a filosofia do Ftool, aliando simplicidade e eficiência: o método dos Dois Ciclos Iterativos. A ferramenta para análise não linear geométrica pode ou não ser utilizada em conjunto com a ferramenta para o dimensionamento de estruturas de concreto armado, sendo possível a análise não linear geométrica de pórticos planos constituídos de outros materiais. A metodologia utilizada para o dimensionamento dos pilares em concreto armado é a que se baseia nas zonas de solicitação e foi adequada para estar de acordo com a norma brasileira vigente, a ABNT NBR 6118:2014, assim como o dimensionamento de vigas existente em uma versão anterior dessa ferramenta. Dessa forma, é possível exibir resultados para pórticos planos compostos por vigas e pilares em concreto armado, nos mesmos moldes da versão anterior, com diagramas para as armaduras longitudinal e transversal disponíveis nos modos necessária e adotada. / [en] The main objective of this work is to complement the reinforced concrete plane frames design tool already existing in Ftool, an educational tool widely known in academia. Therefore, the design of reinforced concrete columns was introduced. A new cross section for rectangular columns with symmetrical steel reinforcement was added. In addition to that, the inclusion of a simplified nonlinear geometric analysis that was in accordance to the philosophy of Ftool, combining simplicity and efficiency, was needed: the Two cycles iterative method. The nonlinear geometric analysis tool may or may not be used together with the reinforced concrete plane frames design tool, thus allowing for geometric nonlinear analyses of plane frames of other materials. The methodology used for the reinforced concrete frames design was based on solicitation zones and was adapted to be in accordance with the Brazilian code, the ABNT NBR 6118:2014, as was the existing reinforced concrete beams design tool. It is now possible to show results for plane frames composed of reinforced concrete columns and beams just as in the previous version of the reinforced concrete design tool, showing diagrams for the necessary and adopted longitudinal and transversal steel reinforcement.

Esse trabalho trata da análise computacional de problemas aeroelásticos. Os casos de interesse envolvem escoamentos compressíveis em torno de estruturas de cascas de materiais compósitos laminados. A solução do problema de interação fluido- estrutura é conduzida por meio do esquema particionado denominado improved serial staggered procedure (ISS), o que permite a utilização de algoritmos independentes para a análise de cada componente do sistema. Um elemento triangular plano com três nós destinado à análise de cascas de materiais compósitos laminados é obtido com base na formulação de dois elementos de alta performance desenvolvidos para a análise de membranas e de placas. A flexibilidade ao corte transversal, a ausência de travamento por razão de aspecto em problemas com flexão no plano da casca são características do elemento. A análise dinâmica não linear é realizada através do procedimento corrotacional de conservação aproximada de energia (AECCP) que foi implementado com base na formulação corrotacional independente de elemento (element independent CR formulation – EICR). Essa abordagem permite que problemas transientes com grandes rotações e forte não linearidade geométrica sejam integrados por longos períodos de tempo, mantendo a estabilidade, a precisão e a eficiência da solução. A simulação do escoamento é realizada através do método explícito de Taylor-Galerkin de dois passos utilizando elementos finitos tetraédricos em malhas não estruturadas para a discretização espacial. As equações governantes do escoamento são dadas segundo a formulação lagrangeana-euleriana arbitrária (ALE), permitindo que malhas móveis sejam empregadas na análise dos problemas aeroelásticos. Um método simples e eficiente é adotado para a movimentação da malha do fluido com base na distância do nó às superfícies de contorno do escoamento. A utilização de malhas não coincidentes na discretização do fluido e da estrutura é possibilitada pelo emprego do esquema de projeção nodal para a transferência de informações entre os meios na região de interface. A apresentação, implementação e verificação de cada componente da ferramenta de análise é realizada independentemente graças à natureza do esquema de solução particionado. Por fim são analisados problemas de interação fluido-estrutura, onde as respostas obtidas são comparadas com resultados experimentais e numéricos de outros autores. Também são propostos e analisados novos problemas envolvendo estruturas de compósitos laminados. / This work deals with computational analysis of aeroelastic problems, specifically those related to compressible flows around laminated composite shell structures. The fluid- structure problem is solved by a partitioned scheme called improved serial staggered procedure (ISS), allowing to use independent algorithms for the analysis of each system component. A 3-node triangular flat element for laminated composite shells is obtained based on the formulation of two high performance plate and membrane elements. Shear flexibility and no aspect ratio locking on in-plane bending problems are some of the element characteristics. The nonlinear dynamic analysis is performed using the approximately energy conserving corotational procedure (AECCP), which is implemented based on the element independent corotational formulation (EICR). This approach enables long time transient problems with strong nonlinearities and large rotations to be efficiently solved while keeping stability and accuracy. The flow simulation is performed by the explicit two step Taylor-Galerkin method employing tetrahedral finite elements in unstructured meshes for the space discrezation. The flow governing equations are given in arbitrary lagrangian-eulerian formulation, enabling aeroelastic simulations with moving meshes. A simple and efficient method is adopted to move the fluid mesh based on the distance of each node to the boundary surfaces. Transference of information between non-matching fluid and structure meshes are performed by the node-projection scheme. Each component of the aeroelastic analysis tool is independently presented, implemented and verified due to the partitioned nature of the scheme adopted for the solution of the coupled system. Some aeroelastic problems are analyzed and results are compared to other experimental and numerical works. New problems of aeroelastic analysis of laminated composite structures are proposed and solved.