Isogeometric configuration design sensitivity analysis of finite deformation curved beam structures using Jaumann strain formulation

Abstract Using an isogeometric approach, a continuum-based configuration design sensitivity analysis (DSA) method is developed for curved Kirchhoff beams with multi-patch junctions. Under the total Lagrangian formulation, large deformations considering the initial curvature of curved beams are described by geometrically exact beam theory (GEBT) and Jaumann strain formulation. In the isogeometric approach, the higher order continuity and the exact description of initial geometry are naturally embedded using NURBS basis functions. In multi-patch models, C 0 -continuity of physical displacement or C 1 -continuity of displacement component at junction is weakly imposed using the Lagrange multiplier method. The superior accuracy of isogeometric analysis (IGA) is verified through the comparison with the results of finite element analysis (FEA) using cubic Hermite interpolation. In the DSA, a material derivative is utilized and the kinematical description of GEBT is consistently employed to express orientation design variations. Contrary to the IGA-based DSA, the Hermite basis function explicitly depends on design in the FEA-based DSA due to its element length parameter. Moreover, since the design velocity field is approximated using the nodal velocity imposed at nodal tangential vector, the amount of design perturbations should be very small to obtain precise design sensitivity.

[1]  Wolfgang A. Wall,et al.  An objective 3D large deformation finite element formulation for geometrically exact curved Kirchhoff rods , 2014 .

[2]  P. Frank Pai,et al.  Highly Flexible Structures : Modeling, Computation, and Experimentation , 2007 .

[3]  Kyung K. Choi,et al.  Meshfree analysis and design sensitivity analysis for shell structures , 2002 .

[4]  Leopoldo Greco,et al.  An implicit G1 multi patch B-spline interpolation for Kirchhoff–Love space rod , 2014 .

[5]  Seung-Hyun Ha,et al.  Isogeometric Shape Design Optimization of Geometrically Nonlinear Structures# , 2013 .

[6]  Alessandro Spadoni,et al.  Isogeometric rotation-free analysis of planar extensible-elastica for static and dynamic applications , 2015 .

[7]  Kyung K. Choi,et al.  Configuration design sensitivity analysis of built-up structures , 1992 .

[8]  Enzo Marino,et al.  Isogeometric collocation for three-dimensional geometrically exact shear-deformable beams , 2016 .

[9]  Habibou Maitournam,et al.  Improved numerical integration for locking treatment in isogeometric structural elements, Part I: Beams , 2014 .

[10]  Les A. Piegl,et al.  The NURBS Book , 1995, Monographs in Visual Communication.

[11]  Seonho Cho,et al.  Isogeometric shape design sensitivity analysis of stress intensity factors for curved crack problems , 2014 .

[12]  M. Crisfield,et al.  Objectivity of strain measures in the geometrically exact three-dimensional beam theory and its finite-element implementation , 1999, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[13]  Yuri Bazilevs,et al.  The bending strip method for isogeometric analysis of Kirchhoff–Love shell structures comprised of multiple patches , 2010 .

[14]  P. Frank Pai,et al.  Three kinematic representations for modeling of highly flexible beams and their applications , 2011 .

[15]  Alessandro Reali,et al.  Avoiding shear locking for the Timoshenko beam problem via isogeometric collocation methods , 2012 .

[16]  P. Frank Pai,et al.  Large-deformation tests and total-Lagrangian finite-element analyses of flexible beams , 2000 .

[17]  F. Filippou,et al.  Response Gradients for Nonlinear Beam-Column Elements under Large Displacements , 2007 .

[18]  Roland Wüchner,et al.  Nonlinear isogeometric spatial Bernoulli Beam , 2016 .

[19]  Alain Combescure,et al.  Locking free isogeometric formulations of curved thick beams , 2012 .

[20]  Yuri Bazilevs,et al.  Isogeometric rotation-free bending-stabilized cables: Statics, dynamics, bending strips and coupling with shells , 2013 .

[21]  Seung-Hyun Ha,et al.  Isogeometric shape design optimization: exact geometry and enhanced sensitivity , 2009 .

[22]  Johannes Gerstmayr,et al.  A continuum mechanics based derivation of Reissner’s large-displacement finite-strain beam theory: the case of plane deformations of originally straight Bernoulli–Euler beams , 2009 .

[23]  Anthony N. Palazotto,et al.  Polar Decomposition Theory in Nonlinear Analyses of Solids and Structures , 1995 .

[24]  Zafer Gürdal,et al.  Isogeometric sizing and shape optimisation of beam structures , 2009 .

[25]  Alessandro Reali,et al.  Locking-free isogeometric collocation methods for spatial Timoshenko rods , 2013 .

[26]  T. Hughes,et al.  Isogeometric analysis : CAD, finite elements, NURBS, exact geometry and mesh refinement , 2005 .

[27]  Kyung K. Choi,et al.  Design sensitivity analysis and optimization of non‐linear transient dynamics. Part II—configuration design , 2000 .

[28]  Luca Dedè,et al.  B-spline goal-oriented error estimators for geometrically nonlinear rods , 2012 .

[29]  Erik Lund,et al.  A Method of “Exact” Numerical Differentiation for Error Elimination in Finite-Element-Based Semi-Analytical Shape Sensitivity Analyses* , 1993 .

[30]  D. F. Rogers,et al.  An Introduction to NURBS: With Historical Perspective , 2011 .

[31]  J. C. Simo,et al.  A finite strain beam formulation. The three-dimensional dynamic problem. Part I , 1985 .

[32]  Seonho Cho,et al.  Isogeometric configuration design optimization of built-up structures , 2015 .

[33]  Jose Manuel Valverde,et al.  Invariant Hermitian finite elements for thin Kirchhoff rods. I: The linear plane case ☆ , 2012 .

[34]  Thomas J. R. Hughes,et al.  Blended isogeometric shells , 2013 .