Density-Based Isogeometric Topology Optimization of Shell Structures

Shell structures with high stiffness-to-weight ratios are desirable in various engineering applications. Topology optimization serves as a popular and effective tool for generating optimal shell structures. The solid isotropic material with penalization (SIMP) method is often chosen because of its simplicity and convenience. However, SIMP method is typically integrated with conventional Finite Element Analysis (FEA) which has limitations in computational accuracy. Achieving high accuracy with FEA necessitates a substantial number of elements, leading to computational burdens. In addition, the discrete representation of the material distribution function may result in rough boundaries. Owing to these limitations, this paper proposes an Isogeometric Analysis (IGA) based SIMP method for optimizing the topology of shell structures based on Reissner–Mindlin theory. This method uses Non-Uniform Rational B-Splines (NURBS) to represent both the shell structure and the material distribution function with the same basis functions, allowing for higher accuracy and smoother boundaries. The optimization model takes compliance as the objective function with a volume fraction constraint and the coefficients of the density function as design variables, resulting in an optimized shell structure defined by the material distribution function. To obtain fairing boundaries of the holes in the optimized shell structure, further process is conducted by fitting the boundaries with fair B-spline curves automatically. Furthermore, the proposed IGA-SIMP framework is applied to generate porous shell structures by imposing different local volume fraction constraints. Numerical examples are provided to demonstrate the feasibility and efficiency of the IGA-SIMP method, showing that it outperforms the FEA-SIMP method and produces smoother boundaries.