Structural shape optimization by IGABEM and particle swarm optimization algorithm

Published on Mar 1, 2018in Engineering Analysis With Boundary Elements2.884
· DOI :10.1016/J.ENGANABOUND.2017.12.007
S. H. Sun1
Estimated H-index: 1
(Hohai University),
Tiantang Yu32
Estimated H-index: 32
(Hohai University)
+ 2 AuthorsTrinh Bui4
Estimated H-index: 4
(TITech: Tokyo Institute of Technology)
Abstract In this paper, a new approach is developed for structural shape optimization, which consists in coupling the particle swarm optimization (PSO) algorithm and the isogeometric boundary element method (IGA-BEM). The IGA-BEM is based on the combination of the isogeometric analysis (IGA) and the boundary element method (BEM), where Non-Uniform Rational B-Splines (NURBS) are employed as shape functions for geometry parameterization and approximation of the field variables. The method inherits the main advantage of the IGA-based shape optimization, i.e., the control points are used as design variables, thus the design model, analysis model and optimization model are uniformly described with the NURBS, providing easy communication between the three models and resulting in a smooth optimized boundary. However, the main feature of the proposed method is the use of PSO, which provides an attractive gradient-free alternative to complicated sensitivity analysis. The efficiency and accuracy of the proposed approach are demonstrated through four 2D shape optimization examples.
📖 Papers frequently viewed together
3,633 Citations
3,547 Citations
76 Citations
#1Haojie Lian (Cardiff University)H-Index: 13
#2Pierre Kerfriden (Cardiff University)H-Index: 14
Last. Stéphane BordasH-Index: 65
view all 3 authors...
Abstract We develop a T-spline isogeometric boundary element method (IGABEM) (Simpson et al., 2012; Scott et al., 2013; Simpson et al., 2014) to shape sensitivity analysis and gradient-based shape optimization in three dimensional linear elasticity. Contrary to finite element based isogeometric analysis (IGA) approaches, no parametrization of the volume is required. Hence, the iterative optimization algorithm can be implemented directly from CAD without any mesh generation or postprocessing step...
75 CitationsSource
#1S.K. Singh (IITR: Indian Institute of Technology Roorkee)H-Index: 6
#2Indra Vir Singh (IITR: Indian Institute of Technology Roorkee)H-Index: 31
Last. Tinh Quoc Bui (TITech: Tokyo Institute of Technology)H-Index: 54
view all 5 authors...
Abstract In the present work, a simple, efficient and accurate Bezier extraction based T-spline XIGA (BEBT-XIGA) has been proposed for the crack simulations. Bezier extraction maps T-spline basis functions over the element in the forms of Bernstein polynomial basis defined over the Bezier element. T-spline basis functions over the element can be written as linear combination of Bernstein polynomials basis. Bezier extraction approach generates an element extraction operator and Bezier element for...
57 CitationsSource
#1Xuan Peng (Cardiff University)H-Index: 11
#2Elena Atroshchenko (University of Chile)H-Index: 11
Last. Stéphane Bordas (University of Luxembourg)H-Index: 65
view all 4 authors...
Abstract We present a novel numerical method to simulate crack growth in 3D, directly from the Computer-Aided Design (CAD) geometry of the component, without any mesh generation. The method is an isogeometric boundary element method (IGABEM) based on non-uniform rational B-splines (NURBS). NURBS basis functions are used for the domain and crack representation as well as to approximate the physical quantities involved in the simulations. A stable quadrature scheme for singular integration is prop...
120 CitationsSource
#1Xuan Peng (Cardiff University)H-Index: 11
#2Elena Atroshchenko (University of Chile)H-Index: 11
Last. Stéphane Bordas (Cardiff University)H-Index: 65
view all 4 authors...
We propose a method for simulating linear elastic crack growth through an isogeometric boundary element method directly from a CAD model and without any mesh generation. To capture the stress singularity around the crack tip, two methods are compared: (1) a graded knot insertion near crack tip; (2) partition of unity enrichment. A well-established CAD algorithm is adopted to generate smooth crack surfaces as the crack grows. The M integral and \(J_k\) integral methods are used for the extraction...
58 CitationsSource
#1Tiantang Yu (Hohai University)H-Index: 32
#2Shuohui Yin (XTU: Xiangtan University)H-Index: 18
Last. Nuttawit Wattanasakulpong (Mahanakorn University of Technology)H-Index: 16
view all 5 authors...
Abstract Practical applications such as airplane wings are usually subjected to combined thermal and mechanical loads, and they hence are prone to buckling failure. Preceding works on the buckling of advanced materials, e.g., functionally graded materials, under combined thermal and mechanical loads are rather rare in literature. In this paper, we report new numerical results of thermal-mechanical buckling of functionally graded rectangular and skew plates (FGPs) under combined thermal and mecha...
68 CitationsSource
#1Shuohui Yin (XTU: Xiangtan University)H-Index: 18
#2Tiantang Yu (Hohai University)H-Index: 32
Last. Satoyuki Tanaka (Hiroshima University)H-Index: 22
view all 5 authors...
Abstract Numerical investigation of buckling and free vibration of functionally graded plates considering in-plane material inhomogeneity is presented in this paper. A novel and effective approach based on isogeometric analysis (IGA) and higher-order shear deformation theory (HSDT) is developed for fulfilling that purpose, analyzing the critical buckling parameter and natural eigenvalues of functionally graded plates involving the in-plane material inhomogeneity. The HSDT allows us to account fo...
69 CitationsSource
#1Chensen Ding (Hunan University)H-Index: 8
#2Xiangyang Cui (Hunan University)H-Index: 35
Last. G.Y. Li (Hunan University)H-Index: 11
view all 3 authors...
Combined isogeometric analysis (IGA) and shape optimization, this paper provides a method that can accurately analyze Tailor Rolled Blanks (TRB) and optimize its thickness profile continuously. While existed traditional methods can hardly solve these problems. TRB has a continuous transition between the thick and thin and a better surface quality. We found in IGA, two-dimension NURBS could present TRB geometry model exactly, but Kirchhoff shell element whose control points are allocated differen...
22 CitationsSource
#1Gernot Beer (Graz University of Technology)H-Index: 20
#2Benjamin Marussig (Graz University of Technology)H-Index: 10
Last. Thomas Peter Fries (Graz University of Technology)H-Index: 23
view all 5 authors...
Abstract In this work a novel approach is presented for the isogeometric Boundary Element analysis of domains that contain inclusions with different elastic properties than the ones used for computing the fundamental solutions. In addition the inclusion may exhibit inelastic material behavior. In this paper only plane stress/strain problems are considered. In our approach the geometry of the inclusion is described using NURBS basis functions. The advantage over currently used methods is that no ...
32 CitationsSource
#1B. H. Nguyen (Vietnamese-German University)H-Index: 7
#2Han D. Tran (Vietnamese-German University)H-Index: 5
Last. Timon Rabczuk (Bauhaus University, Weimar)H-Index: 101
view all 5 authors...
Abstract The isogeometric analysis is applied for the weakly singular symmetric Galerkin boundary element method (SGBEM) to analyze quasi-static elastic problems including crack problems in two-dimensional domains. This method takes the advantages from the common boundary representation of the isogeometric analysis and the boundary element method. The background of the developed method is to use non-uniform rational B-splines (NURBS) for the Galerkin approximation of both geometry and field vari...
74 CitationsSource
#1Haojie Lian (Cardiff University)H-Index: 13
#2Pierre Kerfriden (Cardiff University)H-Index: 21
Last. Stéphane Bordas (Cardiff University)H-Index: 65
view all 3 authors...
The present work addresses shape sensitivity analysis and optimization in two-dimensional elasticity with a regularized isogeometric boundary element method (IGABEM). Non-uniform rational B-splines are used both for the geometry and the basis functions to discretize the regularized boundary integral equations. With the advantage of tight integration of design and analysis, the application of IGABEM in shape optimization reduces the mesh generation/regeneration burden greatly. The work is distinc...
76 CitationsSource
Cited By54
#2Cosmin AnitescuH-Index: 15
Last. Timon RabczukH-Index: 101
view all 4 authors...
#1Zhicheng Wang (HIT: Harbin Institute of Technology)H-Index: 1
#2Zhenggang Cao (HIT: Harbin Institute of Technology)H-Index: 6
Last. Ying Sun (HIT: Harbin Institute of Technology)H-Index: 3
view all 4 authors...
Abstract In this paper, a novel shape optimization method is proposed to improve the mechanical performance of the free-form grid structures while satisfying the architectural requirements. The height of control points are considered as the optimization variables, the structural strain energy and the comprehensive quantitative index are considered as optimization objectives. The sensitivity hybrid multi-objective evolutionary algorithm (referred to as SH-MOEA) is developed and the shape optimiza...
#1Yunluo Yu (University of Birmingham)H-Index: 2
#2Guang Pu (University of Birmingham)H-Index: 2
Last. Kyle Jiang (University of Birmingham)H-Index: 26
view all 4 authors...
#1Rui Dai (BIT: Beijing Institute of Technology)H-Index: 1
#2Chunying Dong (BIT: Beijing Institute of Technology)H-Index: 19
Last. Deyong Sun (BIT: Beijing Institute of Technology)H-Index: 1
view all 4 authors...
Abstract null null Randomly distributed liquid inclusions in elastic matrix exist in nature, biology, material science and other fields. In this paper, the isogeometic analysis boundary element method (IGABEM) is applied to study the mechanical properties of the elastic matrix with liquid inclusions, in which the liquid inclusion is assumed to be linearly compressible and the interface tension is neglected. Singularity subtraction technique (SST) and Tells transformation method are respectively ...
Abstract null null null The capability of boundary element methods (BEM) for solving three-dimensional time harmonic Helmholtz acoustic scattering problems is presented in the framework of the isogeometric analysis (IGA). Both the CAD geometry and the physical boundary variables are approximated using Non-uniform Rational B-splines basis functions (NURBS) in an isogeometric setting. A detailed comparison between two null BEM null methods: the conventional boundary integral equation ( null null n...
2 CitationsSource
#1Antonio Rodrigues Neto (USP: University of São Paulo)H-Index: 1
#2Edson Denner Leonel (USP: University of São Paulo)H-Index: 12
Abstract null null The adequate coupling of dissimilar materials and structural element types provides effective engineering designs. Particularly, it enables high stiffness and low weight mechanical components, which are largely desired in numerous applications. Besides, Isogeometric Analysis (IGA) has recently stood out as a robust and accurate approach for the mechanical modelling of complex 3D engineering designs directly from Computer-Aided Design (CAD) tools. In this context, this study pr...
This paper presents a shape optimisation system to design the shape of an acoustically-hard object in the three-dimensional open space. Boundary element method (BEM) is suitable to analyse such an exterior field. However, the conventional BEM, which is based on piecewise polynomial shape and interpolation functions, can require many design variables because they are usually chosen as a part of the nodes of the underlying boundary element mesh. In addition, it is not easy for the conventional met...
#1Liu Wenfei (Taizhou University)
#2Yongju Zhang (Taizhou University)H-Index: 1
view all 5 authors...
Fatigue cracks usually occur at welded joints of heavy haul wagon body, in order to effectively curb the occurrence of such problems from the design stage, it is researched from structural design s...
#1Daniele M. S. Paulino (PSU: Pennsylvania State University)
#2Edson Denner Leonel (USP: University of São Paulo)H-Index: 12
Topology optimization is an effective approach for the efficient layout design of structures and their components. This approach is well-established in the solid mechanics’ domain for linear elastic and small-displacement conditions. However, the topology optimization of elastic structures under large-displacement conditions has been marginally addressed in the literature. This study proposes a numerical formulation for the topology optimization analysis of plane structures subjected to geometri...
#1Zhao Jing (NPU: Northwestern Polytechnical University)
#2Qin Sun (NPU: Northwestern Polytechnical University)H-Index: 5
Last. Ke Liang (NPU: Northwestern Polytechnical University)H-Index: 5
view all 4 authors...
Due to the classical thin plate theory based governing equations of composite rectangular plates for dynamic problems are linear homogeneous equations on bending stiffness coefficients, a general s...