Nonlinear Analysis of Orthotropic Membrane and Shell Structures Including Fluid-Structure Interaction.

Abstract

Problemas de interacciónn fluido-estructura representan hoy en día un gran desafío en diferentes áreas de ingeniería y ciencias aplicadas. Dentro de las aplicaciones en ingeniería civil, el flujo del viento y los movimientos estructurales pueden ocasionar inestabilidades aeroelásticas en construcciones tales como puentes de gran luz, rascacielos y cubiertas estructurales ligeras. Por otro lado, aplicaciones en biomecánica están interesadas en el estudio de hemodinámica, por ejemplo: flujo sanguíneo en arterias, donde grandes deformaciones de las venas interactúan con fluidos.En la parte estructural de este trabajo, una nueva metodología para el análisis geométricamente no-lineal ortótropo de membranas y láminas sin grados de libertad de rotación es desarrollada basándose en la orientación de la fibra principal del material. Una consecuencia directa de la estrategia de orientación de fibras es la posibilidad de analizar membranas y láminas pretensadas cuya configuración inicial está fuera del plano. Por otra parte, ya que la teoría convencional de membranas permite que existan tensiones de compresión, un modelo de arrugado basado en la modificación de la ecuación constitutiva se presenta. El desarrollo estructural es modelado con elementos finitos provenientes de las ecuaciones de la elastodinámica.<br/>La parte de fluidos de este trabajo está gobernada por las ecuaciones de Navier-<br/>Stokes para flujos incompresibles, las cuales son modeladas por interpolaciones estabilizadas de elementos finitos. Ya que la solución monolítica de dichas ecuaciones tiene la desventaja que consumen mucho tiempo en la solución de grandes sistemas de ecuaciones, el método de pasos fraccionados se usa para aprovechar las ventajas computacionales que brinda gracias al desacoplamiento de la presión del campo de las velocidades. Además, el esquema α-generalizado para integración en el tiempo para fluidos es adaptado para que se use con la t´ecnica de los pasos fraccionados.<br/>El problema de interacción fluido-estructura es formulado como un sistema de tres campos: la estructura, el fluido y el movimiento de la malla. El movimiento del dominio del fluido es tomado en cuenta mediante la formulación arbitraria Lagrangiana-Euleriana, para la cual se usan dos estrategias de movimiento de malla.<br/>Para el acoplamiento entre el fluido y la estructura se usa un acoplamiento fuerte por bloques usando la técnica de Gauss-Seidel. Debido a que la interacción entre el fluido y la estructura es altamente no-lineal, se implementa el método de relajación basado en la técnica de Aitken, la cual acelera la convergencia del problema.<br/>Finalmente varios problemas se presentan en los diferentes campos, los cuales verifican la eficiencia de los algoritmos implementados. Nowadays, fluid-structure interaction problems are a great challenge of different fields in engineering and applied sciences. In civil engineering applications, wind flow and structural motion may lead to aeroelastic instabilities on constructions such as long-span bridges, high-rise buildings and light-weight roof structures. On the other hand, biomechanical applications are interested in the study of hemodynamics, i.e. blood flow through large arteries, where large structural membrane deformations interact with incompressible fluids.<br/>In the structural part of this work, a new methodology for the analysis of geometrically nonlinear orthotropic membrane and rotation-free shell elements is developed based on the principal fiber orientation of the material. A direct consequence of the fiber orientation strategy is the possibility to analyze initially out-ofplane prestressed membrane and shell structures. Additionally, since conventional membrane theory allows compression stresses, a wrinkling algorithm based on modifying the constitutive equation is presented. The structure is modeled with finite elements emerging from the governing equations of elastodynamics.<br/>The fluid portion of this work is governed by the incompressible Navier-Stokes equations, which are modeled by stabilized equal-order interpolation finite elements.<br/>Since the monolithic solution for these equations has the disadvantage that take great computer effort to solve large algebraic system of equations, the fractional step methodology is used to take advantage of the computational efficiency given by the uncoupling of the pressure from the velocity field. In addition, the generalized-α time integration scheme for fluids is adapted to be used with the fractional step technique.<br/>The fluid-structure interaction problem is formulated as a three-field system: the structure, the fluid and the moving fluid mesh solver. Motion of the fluid domain is accounted for with the arbitrary Lagrangian-Eulerian formulation with two different mesh update strategies. The coupling between the fluid and the structure is performed with the strong coupling block Gauss-Seidel partitioned technique.<br/>Since the fluid-structure interaction problem is highly nonlinear, a relaxation technique based on Aitken's method is implemented in the strong coupling formulation to accelerate the convergence.<br/>Finally several example problems are presented in each field to verify the robustness and efficiency of the overall algorithm, many of them validated with reference solutions.