Stress analysis without meshing: isogeometric boundary-element method

Published on Apr 3, 2013
· DOI :10.1680/EACM.11.00024
Haojie Lian13
Estimated H-index: 13
(Cardiff University),
Robert Napier Simpson12
Estimated H-index: 12
(Cardiff University),
Stéphane Bordas65
Estimated H-index: 65
(University of Luxembourg)
Sources
Abstract
The focus of this paper is the description and numerical validation of a computational method where stress analysis can be performed directly from computer-aided design data without mesh generation. The clear benefit of the approach is that no mesh needs to be generated prior to running the analysis. This is achieved by utilising the isogeometric concept whereby computer-aided design data are used to construct not only the geometry discretisation but also the displacement and traction approximations. In this manner, significant savings can be made in the engineering design and analysis process. This paper also demonstrates that, compared with a standard boundary-element method implementation using quadratic Lagrangian shape functions, superior accuracy is achieved using the present approach for the same number of degrees of freedom. It further illustrates practical applications of the method, comparing against results obtained with a standard boundary-element method and finite-element method for verification. In addition, a propeller is analysed as a sample to show the ability of the present method to handle complex three-dimensional geometries.
📖 Papers frequently viewed together
283 Citations
332 Citations
3,558 Citations
References14
Newest
#1Robert Napier Simpson (Cardiff University)H-Index: 12
#2Stéphane Bordas (Cardiff University)H-Index: 65
Last. Timon Rabczuk (Bauhaus University, Weimar)H-Index: 101
view all 4 authors...
The concept of isogeometric analysis, where functions that are used to describe geometry in CAD software are used to approximate the unknown fields in numerical simulations, has received great attention in recent years. The method has the potential to have profound impact on engineering design, since the task of meshing, which in some cases can add significant overhead, has been circumvented. Much of the research effort has been focused on finite element implementations of the isogeometric conce...
263 CitationsSource
#1İlker Temizer (Leibniz University of Hanover)H-Index: 19
#2Peter Wriggers (Leibniz University of Hanover)H-Index: 77
Last. Thomas J. R. Hughes (University of Texas at Austin)H-Index: 152
view all 3 authors...
Abstract We study NURBS-based isogeometric analysis of contact problems and compare with standard C 0 -continuous Lagrange finite elements. A knot-to-surface (KTS) algorithm is developed to treat the contact constraints with NURBS contact surface discretizations. Qualitative studies deliver satisfactory results for various finite deformation frictionless thermoelastic contact problems. Quantitative studies based on the Hertz problem suggest the need for a relaxation of the mechanical contact con...
206 CitationsSource
#1Annalisa BuffaH-Index: 49
#2Giancarlo Sangalli (UNIPV: University of Pavia)H-Index: 35
Last. Rafael VázquezH-Index: 16
view all 3 authors...
Abstract We introduce a new discretization scheme for Maxwell equations in two space dimension. Inspired by the new paradigm of Isogeometric analysis introduced in Hughes et al. (2005) [16] , we propose an algorithm based on the use of bivariate B-splines spaces suitably adapted to electromagnetics. We construct B-splines spaces of variable interelement regularity on the parametric domain. These spaces (and their push-forwards on the physical domain) form a De Rham diagram and we use them to sol...
237 CitationsSource
#1David J. Benson (UCSD: University of California, San Diego)H-Index: 50
#2Yuri Bazilevs (UCSD: University of California, San Diego)H-Index: 77
Last. Thomas J. R. Hughes (University of Texas at Austin)H-Index: 152
view all 4 authors...
Abstract A Reissner–Mindlin shell formulation based on a degenerated solid is implemented for NURBS-based isogeometric analysis. The performance of the approach is examined on a set of linear elastic and nonlinear elasto-plastic benchmark examples. The analyses were performed with LS-DYNA, an industrial, general-purpose finite element code, for which a user-defined shell element capability was implemented. This new feature, to be reported on in subsequent work, allows for the use of NURBS and ot...
472 CitationsSource
#1Wolfgang A. Wall (TUM: Technische Universität München)H-Index: 62
#2Moritz A. Frenzel (TUM: Technische Universität München)H-Index: 2
Last. Christian J. Cyron (TUM: Technische Universität München)H-Index: 18
view all 3 authors...
An old dilemma in structural shape optimization is the needed tight link between design model or geometric description and analysis model. The intention of this paper is to show that isogeometric analysis offers a potential and promising way out of this dilemma. To this end we show a structural shape optimization framework based on the isogeometric analysis approach. With the discretization based on Non-Uniform Rational B-Splines (NURBS) the analysis model represents the structural geometry exac...
387 CitationsSource
#1J. A. Cottrell (University of Texas at Austin)H-Index: 7
#2Alessandro Reali (University of Texas at Austin)H-Index: 45
Last. Thomas J. R. Hughes (University of Texas at Austin)H-Index: 152
view all 4 authors...
This paper begins with personal recollections of John H. Argyris. The geometrical spirit embodied in Argyris’s work is revived in the sequel in applying the newly developed concept of isogeometric analysis to structural vibration problems. After reviewing some fundamentals of isogeometric analysis, application is made to several structural models, including rods, thin beams, membranes, and thin plates. Rotationless beam and plate models are utilized as well as three-dimensional solid models. The...
785 CitationsSource
#1Thomas J. R. HughesH-Index: 152
#2J. A. Cottrell (University of Texas at Austin)H-Index: 7
Last. Yuri Bazilevs (University of Texas at Austin)H-Index: 77
view all 3 authors...
Abstract The concept of isogeometric analysis is proposed. Basis functions generated from NURBS (Non-Uniform Rational B-Splines) are employed to construct an exact geometric model. For purposes of analysis, the basis is refined and/or its order elevated without changing the geometry or its parameterization. Analogues of finite element h - and p -refinement schemes are presented and a new, more efficient, higher-order concept, k -refinement, is introduced. Refinements are easily implemented and e...
3,558 CitationsSource
Jul 1, 2003 in SIGGRAPH (International Conference on Computer Graphics and Interactive Techniques)
#1Thomas W. Sederberg (BYU: Brigham Young University)H-Index: 49
#2Jianmin Zheng (BYU: Brigham Young University)H-Index: 27
Last. Ahmad Nasri (AUB: American University of Beirut)H-Index: 17
view all 4 authors...
This paper presents a generalization of non-uniform B-spline surfaces called T-splines. T-spline control grids permit T-junctions, so lines of control points need not traverse the entire control grid. T-splines support many valuable operations within a consistent framework, such as local refinement, and the merging of several B-spline surfaces that have different knot vectors into a single gap-free model. The paper focuses on T-splines of degree three, which are C2 (in the absence of multiple kn...
661 CitationsSource
#1Les A. Piegl (USF: University of South Florida)H-Index: 41
#2Wayne TillerH-Index: 24
Note: Includes references and index Reference Record created on 2004-09-07, modified on 2016-08-08
3,511 Citations
#2R. ButterfieldH-Index: 19
1,094 Citations
Cited By28
Newest
#1Haifeng Gao (TUT: Taiyuan University of Technology)H-Index: 1
#1Haifeng Gao (TUT: Taiyuan University of Technology)H-Index: 1
Last. Toshiro Matsumoto (Nagoya University)H-Index: 22
view all 5 authors...
Abstract The paper presents a level set based topology optimization method for unidirectional phononic structures with finite layers of lattice cells . Boundary element method(BEM) is employed as the numerical approach to solve the acoustic problems governed by Helmholtz equation . A sized reduced coefficient matrix is derived due to the iteration forms for input and output quantities on the periodic boundary of unit cells. Topological derivatives are formulated by boundary integral equation com...
Source
#1Baotong Li (Xi'an Jiaotong University)H-Index: 14
#2Senmao Ding (Xi'an Jiaotong University)H-Index: 1
Last. Jun Hong (Xi'an Jiaotong University)H-Index: 19
view all 6 authors...
Abstract In this article, we focus on a design problem of planar compliant mechanisms within the framework of isogeometric topology optimization. An integrated model is developed to identify the optimal deformation transferring path for precise motion output. The model comprises two coupled computational layers: the upper layer for geometry representation and the lower layer for sensitivity calculation. In the upper layer, the structural geometries are described explicitly by parameterized level...
Source
#1M. Aminzadeh (University of Shahrood)H-Index: 1
#2S. M. Tavakkoli (University of Shahrood)H-Index: 6
In this research, optimum shape of plate structures is sought to maximize the energy dissipation via structural shape optimization. To achieve this, isogeometric analysis (IGA) is utilized for structural analysis of plates considering elasto-plastic behavior of materials. The von Mises material model is employed for this purpose. Non-uniform rational B-splines basis functions are used for both geometry definition and approximating the unknown deformation field. The optimization problem is to max...
2 CitationsSource
#1Quansheng Zang (DUT: Dalian University of Technology)H-Index: 4
#2Jun Liu (DUT: Dalian University of Technology)H-Index: 15
Last. Gao Lin (DUT: Dalian University of Technology)H-Index: 14
view all 4 authors...
Abstract An isogeometric boundary element method (IGABEM) is proposed for solving the steady-state heat transfer problems with concentrated/surface internal heat sources. The isogeometric boundary element method (IGABEM) possesses the advantages of both the isogeometric analysis (IGA) and the boundary element method (BEM), the non-uniform rational B-spline (NURBS) basis functions used in the (computer-aided design) CAD system are flexible and stable in dealing with irregular boundaries, due to w...
1 CitationsSource
#1Thibault Jacquemin (University of Luxembourg)H-Index: 2
#2Satyendra Tomar (University of Luxembourg)H-Index: 15
Last. Stéphane Bordas (PRC: China Medical University (PRC))H-Index: 65
view all 5 authors...
We provide a primer to numerical methods based on Taylor series expansions such as generalized finite difference methods and collocation methods. We provide a detailed benchmarking strategy for these methods as well as all data files including input files, boundary conditions, point distribution and solution fields, so as to facilitate future benchmarking of new methods. We review traditional methods and recent ones which appeared in the last decade. We aim to help newcomers to the field underst...
3 CitationsSource
#1Leilei Chen (Xinyang Normal University)H-Index: 11
#2Y. Zhang (Xinyang Normal University)H-Index: 1
Last. Stéphane Bordas (Cardiff University)H-Index: 65
view all 6 authors...
Abstract Isogeometric boundary element method is a computer simulation algorithm that can directly utilize the data of a geometric model represented by its surface. Catmull-Clark subdivision surfaces are a widely used technique in 3D computer graphics to construct complicated geometries. In the present work, we combine Catmull-Clark subdivision surfaces with isogeometric boundary element method to simulate acoustic propagation in semi-infinite domains. The present method can produce fast and hig...
4 CitationsSource
#1Quansheng Zang (DUT: Dalian University of Technology)H-Index: 4
#2Jun Liu (DUT: Dalian University of Technology)H-Index: 15
Last. Gao Lin (DUT: Dalian University of Technology)H-Index: 14
view all 4 authors...
Abstract This paper proposes an isogeometric boundary element (IGABEM) with the aim to solve the 2D steady-state heat conduction problems under spatially varying conductivity and internal heat source. The IGABEM boundary-domain integral equation is derived underlying the divergence theorem of Gauss and the Laplace equation. The non-uniform rational B-spline (NURBS) basis used in the CAD/CAE industries is applied to approximate both the boundary of problem domain and the unknown physical quantiti...
1 CitationsSource
#1Yanpeng Gong (Beijing University of Technology)H-Index: 9
#2Chunying Dong (BIT: Beijing Institute of Technology)H-Index: 19
Last. Jon Trevelyan (Durham University)H-Index: 18
view all 5 authors...
Abstract The isogeometric boundary element method (IGABEM) has great potential for the simulation of elasticity problems because of its exact geometric representation and good approximation properties. These advantages can be extended to thin structures, including coatings, but the development of an accurate and efficient method to deal with the large number of nearly singular integrals existing in the IGABEM presents a great challenge for very thin sections. In this paper, we propose a new sinh...
8 CitationsSource
#1Zhijun Liu (Lanzhou University)H-Index: 4
#1Zhijun Liu (Lanzhou University)H-Index: 1
Last. Feng Liu (TJU: Tianjin University)H-Index: 134
view all 4 authors...
Abstract Numerous approaches have been proposed to enhance the accuracy and convergence of the numerical manifold method (NMM) in recent years, but most, if not all, of these approaches cannot ensure C 1 continuity. Hermitian interpolation is an effective approach for obtaining high-order approximations. However, the requirement of rectangular meshes hinders the application of this approach in the finite element method. Taking advantage of the freedom in meshing in NMM, Hermitian interpolation i...
8 CitationsSource
#1Yonghong Wu (BIT: Beijing Institute of Technology)H-Index: 63
#2Chunying Dong (BIT: Beijing Institute of Technology)H-Index: 19
Last. H. S. Yang (BIT: Beijing Institute of Technology)H-Index: 1
view all 3 authors...
Abstract This paper provides an Isogeometric Indirect Boundary Element Method (IGIBEM) based on NURBS (Non-Uniform Rational B-Splines) and PHT-splines (polynomial splines over hierarchical T-meshes) for analyzing the three-dimensional (3D) acoustic problems. In the classical procedure, the geometries are discretized by Lagrange polynomials elements, which leads to both substantial geometrical error and time-consuming meshing steps. However, these deficiencies can be eliminated by the isogeometri...
8 CitationsSource