Spelling suggestions: "subject:"goaloriented adaptive"" "subject:"goalorientated adaptive""
1 |
Cartesian grid FEM (cgFEM): High performance h-adaptive FE analysis with efficient error control. Application to structural shape optimizationNadal Soriano, Enrique 14 February 2014 (has links)
More and more challenging designs are required everyday in today¿s industries.
The traditional trial and error procedure commonly used for mechanical
parts design is not valid any more since it slows down the design process and
yields suboptimal designs. For structural components, one alternative consists
in using shape optimization processes which provide optimal solutions.
However, these techniques require a high computational effort and require
extremely efficient and robust Finite Element (FE) programs. FE software
companies are aware that their current commercial products must improve in
this sense and devote considerable resources to improve their codes. In this
work we propose to use the Cartesian Grid Finite Element Method, cgFEM
as a tool for efficient and robust numerical analysis. The cgFEM methodology
developed in this thesis uses the synergy of a variety of techniques to achieve
this purpose, but the two main ingredients are the use of Cartesian FE grids
independent of the geometry of the component to be analyzed and an efficient
hierarchical data structure. These two features provide to the cgFEM
technology the necessary requirements to increase the efficiency of the cgFEM
code with respect to commercial FE codes. As indicated in [1, 2], in order to
guarantee the convergence of a structural shape optimization process we need
to control the error of each geometry analyzed. In this sense the cgFEM code
also incorporates the appropriate error estimators. These error estimators are
specifically adapted to the cgFEM framework to further increase its efficiency.
This work introduces a solution recovery technique, denoted as SPR-CD, that in combination with the Zienkiewicz and Zhu error estimator [3] provides very
accurate error measures of the FE solution. Additionally, we have also developed
error estimators and numerical bounds in Quantities of Interest based
on the SPR-CD technique to allow for an efficient control of the quality of
the numerical solution. Regarding error estimation, we also present three new
upper error bounding techniques for the error in energy norm of the FE solution,
based on recovery processes. Furthermore, this work also presents an
error estimation procedure to control the quality of the recovered solution in
stresses provided by the SPR-CD technique. Since the recovered stress field
is commonly more accurate and has a higher convergence rate than the FE
solution, we propose to substitute the raw FE solution by the recovered solution
to decrease the computational cost of the numerical analysis. All these
improvements are reflected by the numerical examples of structural shape optimization
problems presented in this thesis. These numerical analysis clearly
show the improved behavior of the cgFEM technology over the classical FE
implementations commonly used in industry. / Cada d'¿a dise¿nos m'as complejos son requeridos por las industrias actuales.
Para el dise¿no de nuevos componentes, los procesos tradicionales de prueba y
error usados com'unmente ya no son v'alidos ya que ralentizan el proceso y dan
lugar a dise¿nos sub-'optimos. Para componentes estructurales, una alternativa
consiste en usar procesos de optimizaci'on de forma estructural los cuales
dan como resultado dise¿nos 'optimos. Sin embargo, estas t'ecnicas requieren
un alto coste computacional y tambi'en programas de Elementos Finitos (EF)
extremadamente eficientes y robustos. Las compa¿n'¿as de programas de EF
son conocedoras de que sus programas comerciales necesitan ser mejorados
en este sentido y destinan importantes cantidades de recursos para mejorar
sus c'odigos. En este trabajo proponemos usar el M'etodo de Elementos Finitos
basado en mallados Cartesianos (cgFEM) como una herramienta eficiente
y robusta para el an'alisis num'erico. La metodolog'¿a cgFEM desarrollada en
esta tesis usa la sinergia entre varias t'ecnicas para lograr este prop'osito, cuyos
dos ingredientes principales son el uso de los mallados Cartesianos de EF independientes
de la geometr'¿a del componente que va a ser analizado y una
eficiente estructura jer'arquica de datos. Estas dos caracter'¿sticas confieren
a la tecnolog'¿a cgFEM de los requisitos necesarios para aumentar la eficiencia
del c'odigo cgFEM con respecto a c'odigos comerciales. Como se indica en
[1, 2], para garantizar la convergencia del proceso de optimizaci'on de forma
estructural se necesita controlar el error en cada geometr'¿a analizada. En
este sentido el c'odigo cgFEM tambi'en incorpora los apropiados estimadores de error. Estos estimadores de error han sido espec'¿ficamente adaptados al
entorno cgFEM para aumentar su eficiencia. En esta tesis se introduce un
proceso de recuperaci'on de la soluci'on, llamado SPR-CD, que en combinaci'on
con el estimador de error de Zienkiewicz y Zhu [3], da como resultado medidas
muy precisas del error de la soluci'on de EF. Adicionalmente, tambi'en se han
desarrollado estimadores de error y cotas num'ericas en Magnitudes de Inter'es
basadas en la t'ecnica SPR-CD para permitir un eficiente control de la calidad
de la soluci'on num'erica. Respecto a la estimaci'on de error, tambi'en se presenta
un proceso de estimaci'on de error para controlar la calidad del campo
de tensiones recuperado obtenido mediante la t'ecnica SPR-CD. Ya que el
campo recuperado es por lo general m'as preciso y tiene un mayor orden de
convergencia que la soluci'on de EF, se propone sustituir la soluci'on de EF por
la soluci'on recuperada para disminuir as'¿ el coste computacional del an'alisis
num'erico. Todas estas mejoras se han reflejado en esta tesis mediante ejemplos
num'ericos de problemas de optimizaci'on de forma estructural. Los resultados
num'ericos muestran claramente un mejor comportamiento de la tecnolog'¿a
cgFEM con respecto a implementaciones cl'asicas de EF com'unmente usadas
en la industria. / Nadal Soriano, E. (2014). Cartesian grid FEM (cgFEM): High performance h-adaptive FE analysis with efficient error control. Application to structural shape optimization [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/35620 / TESIS
|
2 |
GOAL-ORIENTED ERROR ESTIMATION AND ADAPTIVITY FOR HIERARCHICAL MODELS OF THIN ELASTIC STRUCTURESBILLADE, NILESH S. 01 July 2004 (has links)
No description available.
|
3 |
Goal-Oriented Adaptivity using Unconventional Error Representations / Adaptabilité ciblée basée sur des représentations d'erreur non classiquesDarrigrand, Vincent 01 September 2017 (has links)
Dans un contexte d'adaptabilité ciblée, l'erreur commise sur une quantité d'intérêt peut être représentée grâce aux erreurs globales des problèmes direct et adjoint. Cette représentation de l'erreur est majorée par la somme des indicateurs d'erreurs élémentaires. Ces derniers sont alors utilisés pour produire des raffinements de maillage optimaux. Dans ces travaux, nous proposons de représenter l’erreur du problème adjoint via un opérateur alternatif. L’avantage principal de notre approche est que lorsque l'on choisit correctement l'opérateur alternatif, la majoration correspondante de l'erreur à la quantité d'intérêt devient plus précise, pour autant l'adaptabilité issue de l'utilisation de ces nouveaux indicateurs s'en trouve améliorée. Ces représentations peuvent être employées pour concevoir des algorithmes adaptatifs en espace (h), en ordre d’approximation (p) ou les deux (hp), basés sur la norme d’énergie ou bien ciblés sur une quantité d'intérêt. Bien que la méthode puisse être appliquée à une large gamme de problèmes, nous nous concentrons tout d’abord sur des problèmes unidimensionnels (1D), comme le problème d’Helmholtz et le problème de convection-diffusion stationnaire à convection dominante. Les résultats numériques en 1D montrent que, pour les problèmes de propagation d'ondes, les avantages de notre méthode sont notoires lorsque l'on considère l'opérateur de Laplace pour la représentation de l'erreur. Plus précisément, les majorations issues de la nouvelle représentation sont plus précises que celles provenant de la méthode classique et ce si l'on considère l'énergie globale ou bien une quantité d'intérêt particulière. Le phénomène est d’autant plus notable lorsque l'erreur de dispersion (pollution) est significative. Le problème 1D de convection-diffusion stationnaire à convection dominante avec des conditions limites de Dirichlet homogènes présente une couche limite qui produit une perte de stabilité numérique. La nouvelle représentation d'erreur délivre des majorations plus précises. Lorsqu’appliquée à une p-adaptabilité ciblée, la représentation d'erreur alternative permet une capture plus efficace la couche limite, malgré les oscillations numériques parasites existantes. Devant ces résultats encourageants, nous nous penchons sur l'équation d'Helmholtz à deux et trois dimensions (2D et 3D). Nous montrons, au travers de multiples simulations numériques, que les majorations fournies par les représentations d'erreur alternatives sont plus précises que celle de la représentation classique. Lorsque l'on utilise les indicateurs d'erreur alternatifs, un processus naïf de p-adaptabilité ciblée converge, tandis que dans les mêmes conditions, la méthode classique échoue et requiert l'utilisation d'un opérateur de projection ou d'autre techniques pour récupérer la convergence. Dans ce travail, nous fournissons également des directives pour déterminer les opérateurs qui fournissent des représentations d’erreur induisant de majorations précises. Des résultats similaires sont aussi établis tant pour un problème 2D de convection-diffusion stationnaire à convection dominante que pour des problèmes 2D ayant des coefficients de matériaux discontinus. Nous considérons un problème de diagraphie ultra-sonique en cours de forage pour illustrer l'applicabilité de la méthode proposée. / In Goal-Oriented Adaptivity (GOA), the error in a Quantity of Interest (QoI) is represented using global error functions of the direct and adjoint problems. This error representation is subsequently bounded above by element-wise error indicators that are used to drive optimal refinements. In this work, we propose to replace, in the error representation, the adjoint problem by an alternative operator. The main advantage of the proposed approach is that, when judiciously selecting such alternative operator, the corresponding upper bound of the error representation becomes sharper, leading to a more efficient GOA. These representations can be employed to design novel h, p, and hp energy-norm and goal-oriented adaptive algorithms. While the method can be applied to a variety of problems, in this Dissertation we first focus on one-dimensional (1D) problems, including Helmholtz and steady state convection-dominated diffusion problems. Numerical results in 1D show that for the Helmholtz problem, it is advantageous to select the Laplace operator for the alternative error representation. Specifically, the upper bounds of the new error representation are sharper than the classical ones used in both energy-norm and goal-oriented adaptive methods, especially when the dispersion (pollution) error is significant. The 1D steady state convection-dominated diffusion problem with homogeneous Dirichlet boundary conditions exhibits a boundary layer that produces a loss of numerical stability. The new error representation based on the Laplace operator delivers sharper error upper bounds. When applied to a p-GOA, the alternative error representation captures earlier the boundary layer, despite the existing spurious numerical oscillations. We then focus on the two- and three-dimensional (2D and 3D) Helmholtz equation. We show via extensive numerical experimentation that the upper bounds provided by the alternative error representations are sharper than the classical ones. When using the alternative error indicators, a naive p-adaptive process converges, whereas under the same conditions, the classical method fails and requires the use of the so-called Projection Based Interpolation (PBI) operator or some other technique to regain convergence. We also provide guidelines for finding operators delivering sharp error representation upper bounds. / En un contexto de adaptatividad orientada a un objetivo, el error en una cantidad de interés está representado a través de los errores globales de los problemas directo y adjunto. Esta representación del error se acota superiormente por una suma de indicadores de error de cada elemento. Estos se utilizan para producir refinamientos óptimos. En este trabajo, proponemos representar el error del problema adjunto utilizando un operador alternativo. La principal ventaja de nuestro enfoque es que cuando se elige correctamente dicho operador alternativo, la correspondiente cota superior se vuelve más cercana al error en la cantidad de interés, lo que permite una adaptatividad más eficiente. Estas representaciones pueden ser utilizadas para diseñar algoritmos adaptativos en h, p o hp, basados en la norma de la energía o para aproximar una cantidad de interés específica. Aunque el método propuesto se puede aplicar a una amplia gama de problemas, en esta tesis doctoral nos centramos primero en problemas unidimensionales (1D), tales como el problema de Helmholtz y el problema estacionario de convección-difusión con convección dominante. Los resultados numéricos en 1D muestran que, para los problemas de propagación de ondas, las ventajas de este método son notorias cuando se considera el operador de Laplace para la representación del error. Específicamente, las cotas superiores derivadas de la nueva representación son más cercanas a la cantidad de interés que las del método convencional. Esto es cierto tanto para la norma de la energía global como para una cantidad de interés particular, especialmente cuando el error de dispersión es significativo. El problema estacionario 1D de convección-difusión con convección dominante y con condiciones de Dirichlet homogéneas tiene una capa límite que produce una pérdida de estabilidad numérica. La nueva representación del error proporciona cotas superiores más cercanas a la cantidad de interés. Cuando se aplica a un algoritmo adaptativo en p orientado a un objetivo, la representación alternativa del error captura antes la capa límite, a pesar de las existentes oscilaciones numéricas no físicas. En esta tesis doctoral, también nos centramos en la ecuación de Helmholtz en dos y tres dimensiones (2D y 3D). Mostramos a través de múltiples experimentos numéricos que las cotas superiores proporcionadas por las representaciones alternativas del error son más cercanas a la cantidad de interés que cuando uno considera la representación clásica. Al utilizar los indicadores alternativos del error, un algoritmo adaptativo en p sencillo converge, mientras que en las mismas condiciones, el método convencional falla y requiere el uso de operadores de proyección o de otras técnicas para recuperar la convergencia. En este trabajo, también determinamos operadores que proporcionan representaciones del error que inducen cotas superiores más ajustadas. Establecemos resultados similares tanto para el problema estacionario de convección-difusión con convección dominante en 2D como para problemas 2D con materiales discontinuos. Finalmente, se considera un problema sónico en pozos petrolíferos para ilustrar la aplicabilidad del método propuesto.
|
Page generated in 0.11 seconds