Wind field simulation with isogeometric analysis

[EN]For wind field simulation with isogeometric analysis, firstly it is necessary to generate a spline parameterization of the computational domain, which is an air layer above the terrain surface. This parameterization is created with the meccano method from a digital terrain model. The main steps of the meccano method for tetrahedral mesh generation were introduced in [1, 2]. Based on the volume parameterization obtained by the method, we can generate a mapping from the parametric T-mesh to the physical space [3, 4]. Then, this volumetric parameterization is used to generate a cubic spline representation of the physical domain for the application of isogeometric analysis. We consider a mass-consistent model [5] to compute the wind field simulation in the three-dimensional domain from wind measurements or a wind forecasted by a meteorological model (for example, WRF or HARMONIE). From these data, an interpolated wind field is constructed. The mass-consistent model obtains a new wind field approaching the interpolated one, but verifying the continuity equation (mass conservation) for constant density and the impermeabilitycondition on the terrain. This adjusting problem is solved by introducing a Lagrange multiplier, that is the solution of a Poisson problem. The resulting field is obtained from the interpolated one and the gradient of the Lagrange multiplier. It is well known that if we use classical Lagrange finite elements, the gradient of the numerical solution is discontinuous over the element boundary. The advantage of using isogeometric analysis with cubic polynomial basis functions [6, 7] is that we obtain a C2 continuity for the Lagrange multiplier in the whole domain. In consequence, the resulting wind field is better approximated. Applications of the proposed technique are presented.Ministerio de Economía y Competitividad del Gobierno de España; Fondos FEDER; CONACYT-SENE

The Meccano Method for Mesh Generation and Solid Parameterization

[EN]In this conference, we present the main advances in the meccano method, which was introduced by the authors in 2006 for tetrahedral mesh generation

Spline parameterization method for 2D and 3D geometries based on T-mesh optimization

[EN]We present a method to obtain high quality spline parameterization of 2D and 3D geometries for their use in isogeometric analysis. As input data, the proposed method demands a boundary representation of the domain, and it constructs automatically a spline transformation between the physical and parametric domains. Parameterization of the interior of the object is obtained by deforming isomorphically an adapted parametric T-mesh onto the physical domain by applying a T-mesh untangling and smoothing procedure, which is the key of the method. Mesh optimization is based on the mean ratio shape quality measure. The spline representation of the geometry is calculated by imposing interpolation conditions using the data provided by one-to-one correspondence between the meshes of the parametric and physical domains. We give a detailed description of the proposed technique and show some examples. Also, we present some examples of the application of isogeometric analysis in geometries parameterized with our method.Secretaría de Estado de Universidades e Investigación del Ministerio de Economía y Competitividad del Gobierno de España y fondos FEDER; Programa de FPU 12/00202 del Ministerio de Educación, Cultura y Deporte; Programa de FPI propio de la Universidad de Las Palmas de Gran Canari

Construction of polynomial spline spaces over quadtree and octree T-meshes for its application in isogeometric analysis

[EN]We present a new strategy for construction spline spaces over hierarchical T-meshes with quad- and octree subdivision scheme.MINECO y FEDER; CONACYT-SENER Project, Fondo Sectorial, contract: 16372

Structural Shape Optimization Based On The Use Of Cartesian Grids

Tesis por compendioAs ever more challenging designs are required in present-day industries, the traditional trial-and-error procedure frequently used for designing mechanical parts slows down the design process and yields suboptimal designs, so that new approaches are needed to obtain a competitive advantage. With the ascent of the Finite Element Method (FEM) in the engineering community in the 1970s, structural shape optimization arose as a promising area of application. However, due to the iterative nature of shape optimization processes, the handling of large quantities of numerical models along with the approximated character of numerical methods may even dissuade the use of these techniques (or fail to exploit their full potential) because the development time of new products is becoming ever shorter. This Thesis is concerned with the formulation of a 3D methodology based on the Cartesian-grid Finite Element Method (cgFEM) as a tool for efficient and robust numerical analysis. This methodology belongs to the category of embedded (or fictitious) domain discretization techniques in which the key concept is to extend the structural analysis problem to an easy-to-mesh approximation domain that encloses the physical domain boundary. The use of Cartesian grids provides a natural platform for structural shape optimization because the numerical domain is separated from a physical model, which can easily be changed during the optimization procedure without altering the background discretization. Another advantage is the fact that mesh generation becomes a trivial task since the discretization of the numerical domain and its manipulation, in combination with an efficient hierarchical data structure, can be exploited to save computational effort. However, these advantages are challenged by several numerical issues. Basically, the computational effort has moved from the use of expensive meshing algorithms towards the use of, for example, elaborate numerical integration schemes designed to capture the mismatch between the geometrical domain boundary and the embedding finite element mesh. To do this we used a stabilized formulation to impose boundary conditions and developed novel techniques to be able to capture the exact boundary representation of the models. To complete the implementation of a structural shape optimization method an adjunct formulation is used for the differentiation of the design sensitivities required for gradient-based algorithms. The derivatives are not only the variables required for the process, but also compose a powerful tool for projecting information between different designs, or even projecting the information to create h-adapted meshes without going through a full h-adaptive refinement process. The proposed improvements are reflected in the numerical examples included in this Thesis. These analyses clearly show the improved behavior of the cgFEM technology as regards numerical accuracy and computational efficiency, and consequently the suitability of the cgFEM approach for shape optimization or contact problems.La competitividad en la industria actual impone la necesidad de generar nuevos y mejores diseños. El tradicional procedimiento de prueba y error, usado a menudo para el diseño de componentes mecánicos, ralentiza el proceso de diseño y produce diseños subóptimos, por lo que se necesitan nuevos enfoques para obtener una ventaja competitiva. Con el desarrollo del Método de los Elementos Finitos (MEF) en el campo de la ingeniería en la década de 1970, la optimización de forma estructural surgió como un área de aplicación prometedora. El entorno industrial cada vez más exigente implica ciclos cada vez más cortos de desarrollo de nuevos productos. Por tanto, la naturaleza iterativa de los procesos de optimización de forma, que supone el análisis de gran cantidad de geometrías (para las se han de usar modelos numéricos de gran tamaño a fin de limitar el efecto de los errores intrínsecamente asociados a las técnicas numéricas), puede incluso disuadir del uso de estas técnicas. Esta Tesis se centra en la formulación de una metodología 3D basada en el Cartesian-grid Finite Element Method (cgFEM) como herramienta para un análisis numérico eficiente y robusto. Esta metodología pertenece a la categoría de técnicas de discretización Immersed Boundary donde el concepto clave es extender el problema de análisis estructural a un dominio de aproximación, que contiene la frontera del dominio físico, cuya discretización (mallado) resulte sencilla. El uso de mallados cartesianos proporciona una plataforma natural para la optimización de forma estructural porque el dominio numérico está separado del modelo físico, que podrá cambiar libremente durante el procedimiento de optimización sin alterar la discretización subyacente. Otro argumento positivo reside en el hecho de que la generación de malla se convierte en una tarea trivial. La discretización del dominio numérico y su manipulación, en coalición con la eficiencia de una estructura jerárquica de datos, pueden ser explotados para ahorrar coste computacional. Sin embargo, estas ventajas pueden ser cuestionadas por varios problemas numéricos. Básicamente, el esfuerzo computacional se ha desplazado. Del uso de costosos algoritmos de mallado nos movemos hacia el uso de, por ejemplo, esquemas de integración numérica elaborados para poder capturar la discrepancia entre la frontera del dominio geométrico y la malla de elementos finitos que lo embebe. Para ello, utilizamos, por un lado, una formulación de estabilización para imponer condiciones de contorno y, por otro lado, hemos desarrollado nuevas técnicas para poder captar la representación exacta de los modelos geométricos. Para completar la implementación de un método de optimización de forma estructural se usa una formulación adjunta para derivar las sensibilidades de diseño requeridas por los algoritmos basados en gradiente. Las derivadas no son sólo variables requeridas para el proceso, sino una poderosa herramienta para poder proyectar información entre diferentes diseños o, incluso, proyectar la información para crear mallas h-adaptadas sin pasar por un proceso completo de refinamiento h-adaptativo. Las mejoras propuestas se reflejan en los ejemplos numéricos presentados en esta Tesis. Estos análisis muestran claramente el comportamiento superior de la tecnología cgFEM en cuanto a precisión numérica y eficiencia computacional. En consecuencia, el enfoque cgFEM se postula como una herramienta adecuada para la optimización de forma.Actualment, amb la competència existent en la industria, s'imposa la necessitat de generar nous i millors dissenys . El tradicional procediment de prova i error, que amb freqüència es fa servir pel disseny de components mecànics, endarrereix el procés de disseny i produeix dissenys subòptims, pel que es necessiten nous enfocaments per obtindre avantatge competitiu. Amb el desenvolupament del Mètode dels Elements Finits (MEF) en el camp de l'enginyeria en la dècada de 1970, l'optimització de forma estructural va sorgir com un àrea d'aplicació prometedora. No obstant això, a causa de la natura iterativa dels processos d'optimització de forma, la manipulació dels models numèrics en grans quantitats, junt amb l'error de discretització dels mètodes numèrics, pot fins i tot dissuadir de l'ús d'aquestes tècniques (o d'explotar tot el seu potencial), perquè al mateix temps els cicles de desenvolupament de nous productes s'estan acurtant. Esta Tesi se centra en la formulació d'una metodologia 3D basada en el Cartesian-grid Finite Element Method (cgFEM) com a ferramenta per una anàlisi numèrica eficient i sòlida. Esta metodologia pertany a la categoria de tècniques de discretització Immersed Boundary on el concepte clau és expandir el problema d'anàlisi estructural a un domini d'aproximació fàcil de mallar que conté la frontera del domini físic. L'utilització de mallats cartesians proporciona una plataforma natural per l'optimització de forma estructural perquè el domini numèric està separat del model físic, que podria canviar lliurement durant el procediment d'optimització sense alterar la discretització subjacent. A més, un altre argument positiu el trobem en què la generació de malla es converteix en una tasca trivial, ja que la discretització del domini numèric i la seua manipulació, en coalició amb l'eficiència d'una estructura jeràrquica de dades, poden ser explotats per estalviar cost computacional. Tot i això, estos avantatges poden ser qüestionats per diversos problemes numèrics. Bàsicament, l'esforç computacional s'ha desplaçat. De l'ús de costosos algoritmes de mallat ens movem cap a l'ús de, per exemple, esquemes d'integració numèrica elaborats per poder capturar la discrepància entre la frontera del domini geomètric i la malla d'elements finits que ho embeu. Per això, fem ús, d'una banda, d'una formulació d'estabilització per imposar condicions de contorn i, d'un altra, desevolupem noves tècniques per poder captar la representació exacta dels models geomètrics Per completar la implementació d'un mètode d'optimització de forma estructural es fa ús d'una formulació adjunta per derivar les sensibilitats de disseny requerides pels algoritmes basats en gradient. Les derivades no són únicament variables requerides pel procés, sinó una poderosa ferramenta per poder projectar informació entre diferents dissenys o, fins i tot, projectar la informació per crear malles h-adaptades sense passar per un procés complet de refinament h-adaptatiu. Les millores proposades s'evidencien en els exemples numèrics presentats en esta Tesi. Estes anàlisis mostren clarament el comportament superior de la tecnologia cgFEM en tant a precisió numèrica i eficiència computacional. Així, l'enfocament cgFEM es postula com una ferramenta adient per l'optimització de forma.Marco Alacid, O. (2017). Structural Shape Optimization Based On The Use Of Cartesian Grids [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/86195TESISCompendi

An extension of shape sensitivity analysis to an immersed boundary method based on Cartesian grids

[EN] Gradient-based shape optimization processes of mechanical components require the gradients (sensitivity) of the magnitudes of interest to be calculated with sufficient accuracy. The aim of this study was to develop algorithms for the calculation of shape sensitivities considering geometric representation by parametric surfaces (i.e. NURBS or T-splines) using 3D Cartesian h-adapted meshes independent of geometry. A formulation of shape sensitivities was developed for an environment based on Cartesian meshes independent of geometry, which implies, for instance, the need to take into account the special treatment of boundary conditions imposed in non body-fitted meshes. The immersed boundary framework required to implement new methods of velocity field generation, which have a primary role in the integration of both the theoretical concepts and the discretization tools in shape design optimization. Examples of elastic problems with three-dimensional components are given to demonstrate the efficiency of the algorithms.The authors wish to thank the Spanish Ministerio de Economia y Competitividad for the financial support received through the project DPI2013-46317-R and the FPI program (BES-2011-044080), and the Generalitat Valenciana through the Project PROMETEO/2016/007.Marco, O.; Ródenas, JJ.; Fuenmayor Fernández, F.; Tur Valiente, M. (2018). An extension of shape sensitivity analysis to an immersed boundary method based on Cartesian grids. Computational Mechanics. 62(4):701-723. https://doi.org/10.1007/s00466-017-1522-0S701723624Abel JF, Shephard MS (1979) An algorithm for multipoint constraints in finite element analysis. Int J Numer Methods Eng 14(3):464–467Akgün MA, Garcelon GH, Haftka RT (2001) Fast exact linear and nonlinear structural reanalysis and the sherman-morrison-woodbury formulas. Int J Numer Methods Eng 50(7):1587–1606Arora JS (1993) An exposition of the material derivative approach for structural shape sensitivity analysis. Comput Methods Appl Mech Eng 105(1):41–62Arora JS, Lee TH, Cardoso JB (1992) Structural shape sensitivity analysis: relationship between material derivative and control volume approaches. AIAA J 30(6):1638–1648Belegundu D, Zhang S, Manicka Y, Salagame R (1991) The natural approach for shape optimization with mesh distortion control. Technical report. Penn State University, State CollegeBennett JA, Botkin ME (1985) Structural shape optimization with geometric problem description and adaptive mesh refinement. AIAA J 23(3):459–464Bischof C, Carle A, Khademi P, Mauer A (1996) The adifor 2.0 system for the automatic differentiation of fortran 77 programs. IEEE Comput Sci Eng 3(3):18–32Braibant V, Fleury C (1984) Shape optimal design using b-splines. Comput Methods Appl Mech Eng 44(3):247–267Bugeda G, Oliver J (1993) A general methodology for structural shape optimization problems using automatic adaptive remeshing. Int J Numer Methods Eng 36(18):3161–3185Cho S, Ha SH (2009) Isogeometric shape design optimization: exact geometry and enhanced sensitivity. Struct Multidiscip Optim 38(1):53–70Choi K, Kim N (2005) Structural sensitivity analysis and optimization, mechanical engineering series, vol 1. Springer, Berlin, HeidelbergChoi KK, Chang KH (1994) A study of design velocity field computation for shape optimal design. Finite Elem Anal Des 15(4):317–341Choi KK, Duan W (2000) Design sensitivity analysis and shape optimization of structural components with hyperelastic material. Comput Methods Appl Mech Eng 187(1–2):219–243Choi KK, Twu SL (1989) Equivalence of continuum and discrete methods of shape design sensitivity analysis. AIAA J 27(10):1419–1424Choi MJ, Cho S (2014) Isogeometric shape design sensitivity analysis of stress intensity factors for curved crack problems. Comput Methods Appl Mech Eng 279:469–496Chowdhury MS, Song C, Gao W (2014) Shape sensitivity analysis of stress intensity factors by the scaled boundary finite element method. Eng Fract Mech 116:13–30Doctor LJ, Torborg JG (1981) Display techniques for octree-encoded objects. IEEE Comput Graph Appl 1(3):29–38El-Sayed MEM, Zumwalt KW (1991) Efficient design sensitivity derivatives for multi-load case structures as an integrated part of finite element analysis. Comput Struct 40(6):1461–1467Escobar JM, Montenegro R, Rodríguez E, Cascón JM (2014) The meccano method for isogeometric solid modeling and applications. Eng Comput 30(3):331–343Farhat C, Lacour C, Rixen D (1998) Incorporation of linear multipoint constraints in substructure based iterative solvers. Part 1: a numerically scalable algorithm. Int J Numer Methods Eng 43(6):997–1016Fuenmayor FJ, Domínguez J, Giner E, Oliver JL (1997) Calculation of the stress intensity factor and estimation of its error by a shape sensitivity analysis. Fatigue Fract Eng Mater Struct 20(5):813–828Fuenmayor FJ, Oliver JL, Ródenas JJ (1997) Extension of the Zienkiewicz-Zhu error estimator to shape sensitivity analysis. Int J Numer Methods Eng 40(8):1413–1433Gil AJ, Arranz-Carreño A, Bonet J, Hassan O (2010) The immersed structural potential method for haemodynamic applications. J Comput Phys 229(22):8613–8641Giner E, Fuenmayor FJ, Besa AJ, Tur M (2002) An implementation of the stiffness derivative method as a discrete analytical sensitivity analysis and its application to mixed mode in LEFM. Eng Fract Mech 69(18):2051–2071Griewank A, Juedes D, Utke J (1996) ADOL-C, a package for the automatic differentiation of algorithms written in C/C++. ACM Trans Math Softw (TOMS) 22(2):131–167Ha SH, Choi K, Cho S (2010) Numerical method for shape optimization using T-spline based isogeometric method. Struct Multidiscip Optim 42(3):417–428Haftka RT (1993) Semi-analytical static nonlinear structural sensitivity analysis. AIAA J 31(7):1307–1312Haftka RT, Adelman H (1989) Recent developments in structural sensitivity analysis. Struct Optim 1(3):137–151Haftka RT, Barthelemy B (1989) Discretization methods and structural optimization—procedures and applications. In: Proceedings of a GAMM-Seminar October 5–7, 1988, Siegen, FRG, Springer, Berlin, Heidelberg, chap On the Accuracy of shape sensitivity derivatives, pp 136–144Haslinger J, Jedelsky D (1996) Genetic algorithms and fictitious domain based approaches in shape optimization. Struc Optim 12:257–264Hughes TJR, Cottrell JA, Bazilevs Y (2005) Isogeometric analysis: CAD, finite elements, NURBS, exact geometry, and mesh refinement. Comput Methods Appl Mech Eng 194:4135–4195Iman MH (1982) Three-dimensional shape optimization. Int J Numer Methods Eng 18(5):661–673Jackins CL, Tanimoto SL (1980) Oct-tree and their use in representing three-dimensional objects. Comput Graph Image Process 14(3):249–270Kanninen MF, Popelar CH (1985) Advanced fracture mechanics. Oxford Engineering Science Series. Oxford University Press, Oxfordvan Keulen F, Haftka R, Kim N (2005) Review of options for structural design sensitivity analysis. Part I: linear systems. Comput Methods Appl Mech Eng 194(30–33):3213–3243Kibsgaard S (1992) Sensitivity analysis-the basis for optimization. Int J Numer Methods Eng 34(3):901–932Kiendl J, Schmidt R, Wüchner R, Bletzinger KU (2014) Isogeometric shape optimization of shells using semi-analytical sensitivity analysis and sensitivity weighting. Comput Methods Appl Mech Eng 274:148–167Kirsch U (2002) Design-oriented analysis: a unified approach. Springer Netherlands, DordrechtKunisch K, Peichl G (1996) Numerical gradients for shape optimization based on embedding domain techniques. Comput Optim 18:95–114Lee BY (1996) Consideration of body forces in axisymmetric design sensitivity analysis using the bem. Comput Struct 61(4):587–596Lee BY (1997) Direct differentiation formulation for boundary element shape sensitivity analysis of axisymmetric elastic solids. Int J Solids Struct 34(1):99–112Li FZ, Shih CF, Needleman A (1985) A comparison of methods for calculating energy release rates. Eng Fract Mech 21(2):405–421Li K, Qian X (2011) Isogeometric analysis and shape optimization via boundary integral. Comput Aided Des 43(11):1427–1437Lian H, Bordas SPA, Sevilla R, Simpson RN (2012) Recent developments in the integration of computer aided design and analysis. Comput Technol Rev 6:1–36Lian H, Kerfriden P, Bordas SPA (2016) Implementation of regularized isogeometric boundary element methods for gradient-based shape optimization in two-dimensional linear elasticity. Int J Numer Methods Eng 106(12):972–1017Liu L, Zhang Y, Hughes TJR, Scott MA, Sederberg TW (2014) Volumetric T-spline construction using boolean operations. Eng Comput 30(4):425–439Liu WK, Tang S (2007) Mathematical foundations of the immersed finite element method. Comput Mech 39(3):211–222Liu WK, Liu Y, Darell D, Zhang L, Wang XS, Fukui Y, Patankar N, Zhang Y, Bajaj C, Lee J, Hong J, Chen X, Hsu H (2006) Immersed finite element method and its applications to biological systems. Comput Methods Appl Mech Eng 195(13):1722–1749Marco O, Sevilla R, Zhang Y, Ródenas JJ, Tur M (2015) Exact 3D boundary representation in finite element analysis based on Cartesian grids independent of the geometry. Int J Numer Methods Eng 103:445–468Marco O, Ródenas JJ, Navarro-Jiménez JM, Tur M (2017) Robust h-adaptive meshing strategy considering exact arbitrary CAD geometries in a Cartesian grid framework. Comput Struct 193:87–109Meagher D (1980) Octree encoding: a new technique for the representation, manipulation and display of arbitrary 3-D objects by computer. Tech. Rep. IPL-TR-80-11 I, Rensselaer Polytechnic InstituteMoita JS, Infanta J, Mota CM (2000) Sensitivity analysis and optimal design of geometrically non-linear laminated plates and shells. Comput Struct 76(1–3):407–420Nadal E (2014) Cartesian grid FEM (cgFEM): high performance h-adaptive FE analysis with efficient error control. Application to structural shape optimization. Ph.D. Thesis. Universitat Politècnica de ValènciaNadal E, Ródenas JJ, Albelda J, Tur M, Tarancón JE, Fuenmayor FJ (2013) Efficient finite element methodology based on cartesian grids: application to structural shape optimization. Abstract and Applied Analysis 2013Navarrina F, López-Fontán S, Colominas I, Bendito E, Casteleiro M (2000) High order shape design sensitivities: a unified approach. Comput Methods Appl Mech Eng 188(4):681–696Neittaanmäki P, Salmenjoki K (1989) Comparison of various techniques for shape design sensitivity analysis. In: Computer aided optimum design of structures: recent advances. Springer, Berlin, Heidelberg, pp 367–377Nguyen VP, Anitescu C, Bordas SPA, Rabczuk T (2015) Isogeometric analysis: an overview and computer implementation aspects. Math Comput Simul 117:89–116Ozaki I, Kimura F, Berz M (1995) Higher-order sensitivity analysis of finite element method by automatic differentiation. Comput Mech 16(4):223–234Pandey PC, Bakshi P (1999) Analytical response sensitivity computation using hybrid finite elements. Comput Struct 71(5):525–534Peskin CS (1977) Numerical analysis of blood flow in the heart. J Comput Phys 25:220–252Phelan DG, Haber RB (1989) Sensitivity analysis of linear elastic systems using domain parametrization and a mixed mutual energy principle. Comput Methods Appl Mech Eng 77(1–2):31–59Piegl L, Tiller W (1995) The NURBS book. Springer, New YorkPoldneff MJ, Rai IS, Arora JS (1993) Implementation of design sensitivity analysis for nonlinear structures. AIAA J 31(11):2137–2142Qian X (2010) Full analytical sensitivities in NURBS based isogeometric shape optimization. Comput Methods Appl Mech Eng 199(29–32):2059–2071Rice JR (1968) A path independent integral and the approximate analysis of strain concentration by notches and cracks. J Appl Mech 35(2):379–386Ródenas JJ (2001) Error de Discretización en el Cálculo de Sensibilidades mediante el Método de los Elementos Finitos. Ph.D. ThesisRódenas JJ, Fuenmayor FJ, Tarancón JE (2003) A numerical methodology to assess the quality of the design velocity field computation methods in shape sensitivity analysis. Int J Numer Methods Eng 59(13):1725–1747Rodenas JJ, Tur M, Fuenmayor FJ, Vercher A (2007) Improvement of the superconvergent patch recovery technique by the use of constraint equations: the SPR-C technique. Int J Numer Methods Eng 70(6):705–727Ródenas JJ, Bugeda G, Albelda J, Oñate E (2011) On the need for the use of error-controlled finite element analyses in structural shape optimization processes. Int J Numer Methods Eng 87(11):1105–1126Rogers DF (2001) An introduction to NURBS: with historical perspective. Elsevier, LondonSchramm U, Pilkey WW (1993) The coupling of geometric descriptions and finite elements using NURBs a study in shape optimization. Finite Elem Anal Des 15(1):11–34Sevilla R, Fernández-Méndez S, Huerta A (2011) 3D-NURBS-enhanced finite element method (NEFEM). Int J Numer Methods Eng 88(2):103–125Sevilla R, Fernández-Méndez S, Huerta A (2011) NURBS-enhanced finite element method (NEFEM): a seamless bridge between CAD and FEM. Arch Comput Methods Eng 18(4):441–484Shiriaev D, Griewank A (1996) ADOL-F: automatic differentiation of fortran codes. Comput Differ Tech Appl Tools (SIAM) 1:375–384Shivakumar KN, Raju IS (1992) An equivalent domain integral method for three-dimensional mixed-mode fracture problems. Eng Fract Mech 42(6):935–959Silva CAC, Bittencourt ML (2007) Velocity fields using NURBS with distortion control for structural shape optimization. Struct Multidiscip Optim 33(2):147–159Tur M, Albelda J, Nadal E, Ródenas JJ (2014) Imposing dirichlet boundary conditions in hierarchical cartesian meshes by means of stabilized lagrange multipliers. Int J Numer Methods Eng 98(6):399–417Tur M, Albelda J, Marco O, Ródenas JJ (2015) Stabilized method to impose dirichlet boundary conditions using a smooth stress field. Comput Methods Appl Mech Eng 296:352–375Wujek BA, Renaud JE (1998) Automatic differentiation for more efficient system analysis and optimization. Eng Optim 31(1):101–139Yang RJ, Fiedler MJ (1987) Design modelling for large-scale three-dimensional shape optimization problems. ASME Comput Eng 15:177–182Yoon BG, Belegundu AD (1988) Iterative methods for design sensitivity analysis. AIAA J 26(11):1413–1415Zhang L, Gerstenberger A, Wang X, Liu WK (2004) Immersed finite element method. Comput Methods Appl Mech Eng 293(21):2051–2067Zhang W, Beckers P (1989) Comparison of different sensitivity analysis approaches for structural shape optimization. In: Computer aided optimum design of structures: recent advances. Springer, Berlin, Heidelberg, pp 347–356Zhang Y, Wang W, Hughes TJR (2013) Conformal solid T-spline construction from boundary T-spline representations. Comput Mech 6(51):1051–1059Zienkiewicz OC, Zhu JZ (1987) A simple error estimator and adaptive procedure for practical engineering analysis. Int J Numer Methods Eng 24(2):337–357Zienkiewicz OC, Taylor RL, Zhu J (2013) The finite element method: its basis and fundamentals, 7th edn. Butterworth-Heinemann, Oxfor

Structural shape optimization using Cartesian grids and automatic h-adaptive mesh projection

[EN] We present a novel approach to 3D structural shape optimization that leans on an Immersed Boundary Method. A boundary tracking strategy based on evaluating the intersections between a fixed Cartesian grid and the evolving geometry sorts elements as internal, external and intersected. The integration procedure used by the NURBS-Enhanced Finite Element Method accurately accounts for the nonconformity between the fixed embedding discretization and the evolving structural shape, avoiding the creation of a boundary-fitted mesh for each design iteration, yielding in very efficient mesh generation process. A Cartesian hierarchical data structure improves the efficiency of the analyzes, allowing for trivial data sharing between similar entities or for an optimal reordering of thematrices for the solution of the system of equations, among other benefits. Shape optimization requires the sufficiently accurate structural analysis of a large number of different designs, presenting the computational cost for each design as a critical issue. The information required to create 3D Cartesian h- adapted mesh for new geometries is projected from previously analyzed geometries using shape sensitivity results. Then, the refinement criterion permits one to directly build h-adapted mesh on the new designs with a specified and controlled error level. Several examples are presented to show how the techniques here proposed considerably improve the computational efficiency of the optimization process.The authors wish to thank the Spanish Ministerio de Economia y Competitividad for the financial support received through the project DPI2013-46317-R and the FPI program (BES-2011-044080), and the Generalitat Valenciana through the project PROMETEO/2016/007.Marco, O.; Ródenas, J.; Albelda Vitoria, J.; Nadal, E.; Tur Valiente, M. (2017). Structural shape optimization using Cartesian grids and automatic h-adaptive mesh projection. Structural and Multidisciplinary Optimization. 1-21. https://doi.org/10.1007/s00158-017-1875-1S121MATLAB version 8.3.0.532 (R2014a) (2014) Documentation. The Mathworks, Inc., Natick, MassachusettsAbel JF, Shephard MS (1979) An algorithm for multipoint constraints in finite element analysis. Int J Numer Methods Eng 14(3):464–467Amestoy P, Davis T, Duff I (1996) An approximate minimum degree ordering algorithm. SIAM J Matrix Anal Appl 17(4):886–905Barth W, Stürzlinger W (1993) Efficient ray tracing for Bezier and B-spline surfaces. Comput Graph 17 (4):423–430Bennett J A, Botkin M E (1985) Structural shape optimization with geometric problem description and adaptive mesh refinement. AIAA J 23(3):459–464Braibant V, Fleury C (1984) Shape optimal design using b-splines. Comput Methods Appl Mech Eng 44 (3):247–267Bugeda G, Oliver J (1993) A general methodology for structural shape optimization problems using automatic adaptive remeshing. Int J Numer Methods Eng 36(18):3161–3185Bugeda G, Ródenas J J, Oñate E (2008) An integration of a low cost adaptive remeshing strategy in the solution of structural shape optimization problems using evolutionary methods. Comput Struct 86(13–14):1563–1578Chang K, Choi K K (1992) A geometry-based parameterization method for shape design of elastic solids. Mech Struct Mach 20(2):215–252Cho S, Ha S H (2009) Isogeometric shape design optimization: exact geometry and enhanced sensitivity. Struct Multidiscip Optim 38(1):53–70Belegundu D, Zhang YMS, Salagame R (1991) The natural approach for shape optimization with mesh distortion control. Tech. rep., Penn State UniversityDavis T A, Gilbert J R, Larimore S, Ng E (2004) An approximate column minimum degree ordering algorithm. ACM Trans Math Softw 30(3):353–376Doctor L J, Torborg J G (1981) Display techniques for octree-encoded objects. IEEE Comput Graph Appl 1(3):29–38Dunning P D, Kim H A, Mullineux G (2011) Investigation and improvement of sensitivity computation using the area-fraction weighted fixed grid FEM and structural optimization. Finite Elem Anal Des 47(8):933–941Düster A, Parvizian J, Yang Z, Rank E (2008) The finite cell method for three-dimensional problems of solid mechanics. Comput Methods Appl Mech Eng 197(45-48):3768–3782Escobar J M, Montenegro R, Rodríguez E, Cascón J M (2014) The meccano method for isogeometric solid modeling and applications. Eng Comput 30(3):331–343Farhat C, Lacour C, Rixen D (1998) Incorporation of linear multipoint constraints in substructure based iterative solvers. Part 1: a numerically scalable algorithm. Int J Numer Methods Eng 43(6):997–1016Fries T P, Omerović S (2016) Higher-order accurate integration of implicit geometries. Int J Numer Methods Eng 106(5):323–371Fuenmayor F J, Oliver J L (1996) Criteria to achieve nearly optimal meshes in the h-adaptive finite element mehod. Int J Numer Methods Eng 39(23):4039–4061Fuenmayor F J, Oliver J L, Ródenas J J (1997) Extension of the Zienkiewicz-Zhu error estimator to shape sensitivity analysis. Int J Numer Methods Eng 40(8):1413–1433García-Ruíz M J, Steven G P (1999) Fixed grid finite elements in elasticity problems. Eng Comput 16 (2):145–164Gill P, Murray W, Saunders M, Wright M (1984) Procedures for optimization problems with a mixture of bounds and general linear constraints. ACM Trans Math Software 10:282–298González-Estrada O A, Nadal E, Ródenas J J, Kerfriden P, Bordas S P A, Fuenmayor F J (2014) Mesh adaptivity driven by goal-oriented locally equilibrated superconvergent patch recovery. Comput Mech 53(5):957–976Ha S H, Choi K K, Cho S (2010) Numerical method for shape optimization using T-spline based isogeometric method. Struct Multidiscip Optim 42(3):417–428Haftka R T, Grandhi R V (1986) Structural shape optimization: A survey. Comput Methods Appl Mech Eng 57(1):91–106Haslinger J, Jedelsky D (1996) Genetic algorithms and fictitious domain based approaches in shape optimization. Struc Optim 12:257–264Hughes T J R, Cottrell J A, Bazilevs Y (2005) Isogeometric Analysis: CAD, Finite Elements, NURBS, Exact Geometry, and Mesh Refinement. Comput Methods Appl Mech Eng 194:4135–4195Jackins C L, Tanimoto S L (1980) Oct-tree and their use in representing three-dimensional objects. Comput Graphics Image Process 14(3):249–270Kajiya J T (1982) Ray Tracing Parametric Patches. SIGGRAPH Comput Graph 16(3):245–254van Keulen F, Haftka R T, Kim N (2005) Review of options for structural design sensitivity analysis. Part I: linear systems. Comput Methods Appl Mech Eng 194(30-33):3213–3243Kibsgaard S (1992) Sensitivity analysis-the basis for optimization. Int J Numer Methods Eng 34(3):901–932Kikuchi N, Chung K Y, Torigaki T, Taylor J E (1986) Adaptive finite element methods for shape optimization of linearly elastic structures. Comput Methods Appl Mech Eng 57(1):67–89Kim N H, Chang Y (2005) Eulerian shape design sensitivity analysis and optimization with a fixed grid. Comput Methods Appl Mech Eng 194(30–33):3291–3314Kudela L, Zander N, Kollmannsberger S, Rank E (2016) Smart octrees: Accurately integrating discontinuous functions in 3d. Comput Methods Appl Mech Eng 306(1):406–426Kunisch K, Peichl G (1996) Numerical gradients for shape optimization based on embedding domain techniques. Comput Optim 18:95–114Li K, Qian X (2011) Isogeometric analysis and shape optimization via boundary integral. Computer-Aided Design 43(11):1427–1437Lian H, Kerfriden P, Bordas S P A (2016) Implementation of regularized isogeometric boundary element methods for gradient-based shape optimization in two-dimensional linear elasticity. Int J Numer Methods Eng 106 (12):972–1017Liu L, Zhang Y, Hughes T J R, Scott M A, Sederberg T W (2014) Volumetric T-spline Construction using Boolean Operations. Eng Comput 30(4):425–439Marco O, Sevilla R, Zhang Y, Ródenas J J, Tur M (2015) Exact 3D boundary representation in finite element analysis based on Cartesian grids independent of the geometry. Int J Numer Methods Eng 103:445–468Marco O, Ródenas J J, Fuenmayor FJ, Tur M (2017a) An extension of shape sensitivity analysis to an immersed boundary method based on cartesian grids. Computational Mechanics SubmittedMarco O, Ródenas J J, Navarro-Jiménez JM, Tur M (2017b) Robust h-adaptive meshing strategy for arbitrary cad geometries in a cartesian grid framework. Computers & Structures SubmittedMeagher D (1980) Octree Encoding: A New Technique for the Representation, Manipulation and Display of Arbitrary 3-D Objects by Computer. Tech. Rep. IPL-TR-80-11 I, Rensselaer Polytechnic InstituteMoita J S, Infante J, Mota C M, Mota C A (2000) Sensitivity analysis and optimal design of geometrically non-linear laminated plates and shells. Comput Struct 76(1–3):407–420Nadal E (2014) Cartesian Grid FEM (cgFEM): High Performance h-adaptive FE Analysis with Efficient Error Control. Application to Structural Shape Optimization. PhD Thesis. Universitat Politècnica de ValènciaNadal E, Ródenas J J, Albelda J, Tur M, Tarancón J E, Fuenmayor F J (2013) Efficient finite element methodology based on cartesian grids: application to structural shape optimization. Abstr Appl Anal 2013:1–19Najafi A R, Safdari M, Tortorelli D A, Geubelle P H (2015) A gradient-based shape optimization scheme using an interface-enriched generalized FEM. Comput Methods Appl Mech Eng 296:1–17Nguyen V P, Anitescu C, Bordas S P A, Rabczuk T (2015) Isogeometric analysis: An overview and computer implementation aspects. Math Comput Simul 117:89–116Nishita T, Sederberg TW, Kakimoto M (1990) Ray Tracing Trimmed Rational Surface Patches. SIGGRAPH Comput Graph 24(4):337–345Nocedal J, Wright SJ (2006) Numerical optimization, 2nd edn. Springer-Verlag, New YorkPandey P C, Bakshi P (1999) Analytical response sensitivity computation using hybrid finite elements. Comput Struct 71(5):525–534Parvizian J, Düster A, Rank E (2007) Finite Cell Method: h- and p- Extension for Embedded Domain Methods in Solid Mechanics. Comput Mech 41(1):121–133Peskin C S (1977) Numerical Analysis of Blood Flow in the Heart. J Comput Phys 25:220–252Poldneff M J, Rai I S, Arora J S (1993) Implementation of design sensitivity analysis for nonlinear structures. AIAA J 31(11):2137–2142Powell M (1983) Variable metric methods for constrained optimization. In: Bachem A, Grotschel M, Korte B (eds) Mathematical Programming: The State of the Art, Springer, Berlin, Heidelberg, pp 288–311Qian X (2010) Full analytical sensitivities in NURBS based isogeometric shape optimization. Comput Methods Appl Mech Eng 199(29–32):2059–2071Riehl S, Steinmann P (2014) An integrated approach to shape optimization and mesh adaptivity based on material residual forces. Comput Methods Appl Mech Eng 278:640–663Riehl S, Steinmann P (2016) On structural shape optimization using an embedding domain discretization technique. Int J Numer Methods Eng 109(9):1315–1343Ródenas J J, Tarancón J E, Albelda J, Roda A, Fuenmayor F J (2005) Hierarchical Properties in Elements Obtained by Subdivision: a Hierarquical h-adaptivity Program. In: Díez P, Wiberg N E (eds) Adaptive Modeling and Simulation, p 2005Ródenas J J, Corral C, Albelda J, Mas J, Adam C (2007a) Nested domain decomposition direct and iterative solvers based on a hierarchical h-adaptive finite element code. In: Runesson K, Díez P (eds) Adaptive Modeling and Simulation 2007, Internacional Center for Numerical Methods in Engineering (CIMNE), pp 206–209Ródenas J J, Tur M, Fuenmayor F J, Vercher A (2007b) Improvement of the superconvergent patch recovery technique by the use of constraint equations: the SPR-C technique. Int J Numer Methods Eng 70(6):705–727Ródenas J J, Bugeda G, Albelda J, Oñate E (2011) On the need for the use of error-controlled finite element analyses in structural shape optimization processes. Int J Numer Methods Eng 87(11):1105–1126Schillinger D, Ruess M (2015) The finite cell method: A review in the context of higher-order structural analysis of cad and image-based geometric models. Arch Comput Meth Eng 22(3):391– 455Sevilla R, Fernández-Méndez S, Huerta A (2011a) 3D-NURBS-enhanced Finite Element Method (NEFEM). Int J Numer Methods Eng 88(2):103–125Sevilla R, Fernández-Méndez S, Huerta A (2011b) Comparison of High-order Curved Finite Elements. Int J Numer Methods Eng 87(8):719–734Sevilla R, Fernández-Méndez S, Huerta A (2011c) NURBS-enhanced Finite Element Method (NEFEM): A Seamless Bridge Between CAD and FEM. Arch Comput Meth Eng 18(4):441–484Sweeney M, Bartels R (1986) Ray tracing free-form b-spline surfaces. IEEE Comput Graph Appl 6(2):41–49Toth D L (1985) On Ray Tracing Parametric Surfaces. SIGGRAPH Comput Graph 19(3):171–179Tur M, Albelda J, Nadal E, Ródenas J J (2014) Imposing dirichlet boundary conditions in hierarchical cartesian meshes by means of stabilized lagrange multipliers. Int J Numer Methods Eng 98(6):399–417Tur M, Albelda J, Marco O, Ródenas J J (2015) Stabilized Method to Impose Dirichlet Boundary Conditions using a Smooth Stress Field. Comput Methods Appl Mech Eng 296:352–375Yao T, Choi KK (1989) 3-d shape optimal design and automatic finite element regridding. Int J Numer Methods Eng 28(2):369–384Zhang L, Gerstenberger A, Wang X, Liu W K (2004) Immersed Finite Element Method. Comput Methods Appl Mech Eng 293(21):2051–2067Zhang Y, Wang W, Hughes T J R (2013) Conformal Solid T-spline Construction from Boundary T-spline Representations. Comput Mech 6(51):1051–1059Zienkiewicz O C, Zhu J Z (1987) A Simple Error Estimator and Adaptive Procedure for Practical Engineering Analysis. Int J Numer Methods Eng 24(2):337–35

Generalized averaged Gaussian quadrature and applications

A simple numerical method for constructing the optimal generalized averaged Gaussian quadrature formulas will be presented. These formulas exist in many cases in which real positive GaussKronrod formulas do not exist, and can be used as an adequate alternative in order to estimate the error of a Gaussian rule. We also investigate the conditions under which the optimal averaged Gaussian quadrature formulas and their truncated variants are internal

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described