Isogeometric shape optimization of nonlinear, curved 3D beams and beam structures

• Published in: Computer methods in applied mechanics and engineering Vol. 345; pp. 26 - 51
• Main Authors: Weeger, Oliver, Narayanan, Bharath, Dunn, Martin L.
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier B.V 01.03.2019; Elsevier BV
• Subjects: 3D beams; Algorithms; Beams (structural); Collocation method; Collocation methods; Curved beams; Curves; Deformation; Formability; Freeform fabrication; Geometrically exact beam theory; Isogeometric analysis; Lattices (mathematics); Mathematical analysis; Mathematical models; Mechanical engineering; Mechanical properties; Metamaterials; Nonlinear equations; Nonlinear optimization; Numerical methods; Sensitivity analysis; Shape optimization; Three dimensional models; Trusses; Isogeometric analysis; Collocation method; Geometrically exact beam theory; 3D beams; Shape optimization; Nonlinear optimization
• ISSN: 0045-7825; 1879-2138
• DOI: 10.1016/j.cma.2018.10.038

Straight beams, rods and trusses are common elements in structural and mechanical engineering, but recent advances in additive manufacturing now also enable efficient freeform fabrication of curved, deformable beams and beam structures, such as microstructures, metamaterials and conformal lattices. To exploit this new design freedom for applications with nonlinear mechanical behavior, we introduce an isogeometric method for shape optimization of curved 3D beams and beam structures. The geometrically exact Cosserat rod theory is used to model nonlinear 3D beams subject to large deformations and rotations. The initial and current geometry are parameterized in terms of NURBS curves describing the beam centerline and an isogeometric collocation approach is used to discretize the strong form of the balance equations. Then, a nonlinear optimization problem is formulated in order to optimize the positions of the control points of the NURBS curve that describes the beam centerline, i.e., the geometry or shape of the beam. To solve the design problem using gradient-based algorithms, we introduce semi-analytical, inconsistent analytical and fully analytical approaches for calculation of design sensitivities. The methods are numerically validated and their performance is investigated, before the applicability and versatility of our 3D beam shape optimization method is illustrated in various numerical applications, including optimization of an auxetic 3D metamaterial. •Shape of curved, freeform, nonlinear 3D beams is optimized in terms of centerline curve.•NURBS-based integration of isogeometric collocation for nonlinear analysis and optimization.•Control points of centerline curves serve directly as design variables for optimization.•Incomplete analytical and fully analytical sensitivities are derived for gradient-based optimization.•Numerical validation and application to optimization of an auxetic metamaterial.