Ingeniero Civil Mecánico / Los fenómenos físicos, presentes en las ciencias y en las diferentes áreas de la ingeniería, a menudo son modelados por Ecuaciones Diferenciales Parciales (EDP). Los problemas de valor de frontera resultantes en muchos casos carecen de soluciones analíticas. Para resolver tales problemas, uno puede hacer suposiciones que simplifiquen el problema, o usar métodos numéricos para aproximar la solución. Dentro de los métodos numéricos actualmente existentes, el más popular es el Método de Elementos Finitos (FEM), que es la base de diferentes programas comerciales, como ADINA o ANSYS, entre muchos otros. La desventaja de este método es la gran cantidad de recursos computacionales y los tiempos de iteración requeridos para obtener una solución precisa del problema. Dada esta desventaja, Hughes desarrolló el Análisis IsoGeométrico (IGA). Este método permite integrar el modelo CAD con el Análisis de Elementos Finitos (FEA), por lo tanto, reduce los tiempos y los recursos necesarios para obtener una solución precisa. Pero a su vez, el IGA no tiene flexibilidad para obtener soluciones de ciertos problemas, ya que usa las mismas funciones bases para parametrizar tanto la geometría como el campo de solución. Debido a esto último, surge el Análisis IsoGeométrico Generalizado (GIFT) como una generalización del IGA, este método utiliza diferentes funciones bases para parametrizar la geometría del objeto y el campo de solución, permitiendo la selección de funciones que se adapten mejor al problema estudiado. En trabajos anteriores, el GIFT ha sido aplicado a problemas de la Ecuación de Laplace y de Elasticidad Lineal. El objetivo principal de este trabajo es estudiar el rendimiento del GIFT para problemas de flexión y de vibraciones de placas delgadas. El estudio consiste en implementar el GIFT para 3 placas diferentes y comparar los resultados numéricos con lo predicho por la Teoría de Placas de Kirchhoff-Love (KLPT). Se consideran una placa de geometría circular simple, una placa de geometría circular de dos parches y una placa cuadrada con un agujero de forma compleja, modelada por 8 parches. Las placas están parametrizadas por NURBS, mientras que las soluciones se aproximan por un parche usando NURBS o B-Splines. Los resultados se muestran en términos de curvas de convergencia, modos de vibración y frecuencias naturales. Los resultados numéricos se comparan con las soluciones analíticas para problemas con geometría simple y con la solución FEM para el problema de una placa más compleja. El análisis realizado indica que, para la misma parametrización de geometría (uniforme), (a) la solución se puede aproximar mediante un parche NURBS o B-Splines, manteniendo inalterada la geometría original, (b) los resultados obtenidos con las aproximaciones de campo NURBS y B-Splines son idénticas, (c) la tasa de convergencia depende del grado de aproximación de la solución. Para parametrizaciones geométricas no uniformes, el método no produce una tasa de convergencia óptima o resultados suficientemente precisos, al igual que el IGA tradicional.

Tesis para optar al grado de Magíster en Ciencias de la Ingeniería, Mención Mecánica / Both the CAD software and FEM software have a significant impact on engineering nowadays. Even though both are powerful tools for design and analysis, the main drawback is that CAD geometries and Finite Element models do not entirely match, which results in the necessity to re-parameterize the geometry many times during the solution cycle in FEM. Isogeometric Analysis (IGA) was proposed to fulfill this gap and create the direct link between the CAD design and FEM analysis. The main idea of IGA is to substitute the shape functions used in FEM by the shape functions used in the CAD software. In particular, one of the main drawbacks of NURBS basis functions, and therefore of IGA, is the lack of local refinement, which makes them computationally highly expensive in applications that demands a non-uniform refinement of the geometry. Polynomial splines over Hierarchical T-meshes (PHT-splines) were introduced by Deng et al. as a type of spline that allows local refinement and adaptability by means of a polynomial basis capable of parameterizing the geometry. In this work, we demonstrate the application of PHT-splines for two type of problems: time-harmonic acoustic problems, modeled by the Helmholtz equation, and fracture mechanics of thin plate problems, modeled by the Kirchhoff-Love theory. Solutions of the Helmholtz equation have two features: global oscillations associated with the wave number and local gradients caused by geometrical irregularities. The results show that after a sufficient number of degrees of freedom is used to approximate global oscillations, adaptive refinement can capture local features of the solution. The residual-based and recovery-based error estimators are compared and the performance of $p$-refinement is investigated. Moreover, an eXtended Geometry Independent Field approximaTion (XGIFT) formulation based on Polynomials Splines Over Hierarchical T-meshes (PHT-splines) for modeling both static and dynamic fracture mechanic problems for plates described by the Kirchhoff-Love theory is presented. Adaptive refinement is employed using a recovery-based error estimator. Results show that adaptive refinement can capture local features of the solution around the crack tip, improving results in both static and dynamic examples. In both cases, the simulations are done in the context of recently introduced Geometry Independent Field approximaTion (GIFT), where PHT-splines are only used to approximate the solution, while the computational domain is parameterized with NURBS. This approach builds on the natural adaptation ability of PHT-splines and avoids the re-parameterization of the NURBS geometry during the solution refinement process.