An nth high order perturbation-based stochastic isogeometric method and implementation for quantifying geometric uncertainty in shell structures

Abstract This paper presents an n-th high order perturbation-based stochastic isogeometric Kirchhoff–Love shell method, formulation and implementation for modeling and quantifying geometric (thickness) uncertainty in thin shell structures. Firstly, the Non-Uniform Rational B-Splines (NURBS) is used to describe the geometry and interpolate the variables in a deterministic aspect. Then, the shell structures with geometric (thickness) uncertainty are investigated by developing an nth order perturbation-based stochastic isogeometric method. Here, we develop the shell stochastic formulations in detail (particularly, expand the random input (thickness) and IGA Kirchhoff-Love shell element based state functions analytically around their expectations via n-th order Taylor series using a small perturbation parametere), whilst freshly providing the Matlab core codes helpful for implementation. This work includes three key novelties: 1) by increasing/utilizing the high order of NURBS basis functions, we can exactly represent shell geometries and alleviate shear locking, as well as providing more accurate deterministic solution hence enhancing stochastic response accuracy. 2) Via increasing the nth order perturbation, we overcome the inherent drawbacks of first and second-order perturbation approaches, and hence can handle uncertainty problems with some large coefficients of variation. 3) The numerical examples, including two benchmarks and one engineering application (B-pillar in automobile), simulated by the proposed formulations and direct Monte Carlo simulations (MCS) verify that thickness randomness does strongly affect the response of shell structures, such as the displacement caused by uncertainty can increase up to 35%; Moreover, the proposed formulation is effective and significantly efficient. For example, compared to MCS, only 0.014% computational time is needed to obtain the stochastic response.

[1]  M. Ortiz,et al.  Subdivision surfaces: a new paradigm for thin‐shell finite‐element analysis , 2000 .

[2]  Guangyao Li,et al.  Isogeometric generalized n th order perturbation-based stochastic method for exact geometric modeling of (composite) structures: Static and dynamic analysis with random material parameters , 2019, Computer Methods in Applied Mechanics and Engineering.

[3]  George Stefanou,et al.  Stochastic finite element analysis of shells , 2002 .

[4]  Marcin Kamiński,et al.  Generalized perturbation-based stochastic finite element method in elastostatics , 2007 .

[5]  Timon Rabczuk,et al.  An alternative alpha finite element method with discrete shear gap technique for analysis of isotropic Mindlin-Reissner plates , 2011 .

[6]  Guangyao Li,et al.  Analysis of plates and shells using an edge-based smoothed finite element method , 2009 .

[7]  Stéphane Bordas,et al.  Isogeometric analysis of functionally graded plates using a refined plate theory , 2014 .

[8]  Manolis Papadrakakis,et al.  The effect of material and thickness variability on the buckling load of shells with random initial imperfections , 2005 .

[9]  Marcin Marek Kaminski,et al.  Generalized stochastic perturbation technique in engineering computations , 2010, Math. Comput. Model..

[10]  Stéphane Bordas,et al.  Isogeometric analysis of laminated composite and sandwich plates using a new inverse trigonometric shear deformation theory , 2014 .

[11]  H. Nguyen-Xuan,et al.  Isogeometric analysis of large-deformation thin shells using RHT-splines for multiple-patch coupling , 2017 .

[12]  Manolis Papadrakakis,et al.  A computationally efficient method for the buckling analysis of shells with stochastic imperfections , 2009 .

[13]  Kumar K. Tamma,et al.  Proper orthogonal decomposition and Monte Carlo based isogeometric stochastic method for material, geometric and force multi-dimensional uncertainties , 2018, Computational Mechanics.

[14]  J. C. Simo,et al.  A CLASS OF MIXED ASSUMED STRAIN METHODS AND THE METHOD OF INCOMPATIBLE MODES , 1990 .

[15]  G. Stefanou.  The stochastic finite element method: Past, present and future , 2009 .

[16]  L. Heltai,et al.  Nonsingular isogeometric boundary element method for Stokes flows in 3D , 2014 .

[17]  Guangyao Li,et al.  A fast uncertainty analysis for material nonlinear problems based on reanalysis and semi-analytic strategy , 2016, Adv. Eng. Softw..

[18]  John A. Evans,et al.  Isogeometric analysis using T-splines , 2010 .

[19]  R. Ghanem,et al.  Stochastic Finite Elements: A Spectral Approach , 1990 .

[20]  Elena Atroshchenko,et al.  Weakening the tight coupling between geometry and simulation in isogeometric analysis: From sub‐ and super‐geometric analysis to Geometry‐Independent Field approximaTion (GIFT) , 2016, 1706.06371.

[21]  R. Ghanem,et al.  Iterative solution of systems of linear equations arising in the context of stochastic finite elements , 2000 .

[22]  Zhiguo Wang,et al.  3-D elasto-plastic large deformations: IGA simulation by Bézier extraction of NURBS , 2017, Adv. Eng. Softw..

[23]  Yang Zeng,et al.  A Copula-based perturbation stochastic method for fiber-reinforced composite structures with correlations , 2017 .

[24]  Xiangyang Cui,et al.  Stochastic analysis using the generalized perturbation stable node-based smoothed finite element method , 2016 .

[25]  T. Rabczuk,et al.  A two-dimensional Isogeometric Boundary Element Method for elastostatic analysis , 2012 .

[26]  Thomas J. R. Hughes,et al.  Isogeometric Analysis: Toward Integration of CAD and FEA , 2009 .

[27]  Anh-Vu Vuong,et al.  ISOGAT: A 2D tutorial MATLAB code for Isogeometric Analysis , 2010, Comput. Aided Geom. Des..

[28]  Roger A. Sauer,et al.  A NURBS-based Inverse Analysis for Reconstruction of Nonlinear Deformations of Thin Shell Structures , 2018, ArXiv.

[29]  Ming Li,et al.  Constructing IGA-suitable planar parameterization from complex CAD boundary by domain partition and global/local optimization , 2017, ArXiv.

[30]  Guirong Liu,et al.  Generalized stochastic cell-based smoothed finite element method (GS_CS-FEM) for solid mechanics , 2013 .

[31]  Chang-Koon Choi,et al.  Stochastic finite element analysis of plate structures by weighted integral method , 1996 .

[32]  Di Wu,et al.  Spectral stochastic isogeometric analysis for linear stability analysis of plate , 2019, Computer Methods in Applied Mechanics and Engineering.

[33]  Guangyao Li,et al.  Exact and efficient isogeometric reanalysis of accurate shape and boundary modifications , 2017 .

[34]  J. Trevelyan,et al.  An isogeometric boundary element method for elastostatic analysis: 2D implementation aspects , 2013, 1302.5305.

[35]  Alex H. Barbat,et al.  Monte Carlo techniques in computational stochastic mechanics , 1998 .

[36]  Michael B. Giles,et al.  Multilevel Monte Carlo methods , 2013, Acta Numerica.

[37]  Chao Wang,et al.  Implementation of isogeometric boundary element method for 2-D steady heat transfer analysis , 2018, Adv. Eng. Softw..

[38]  Daniel Rypl,et al.  Object oriented implementation of the T-spline based isogeometric analysis , 2012, Adv. Eng. Softw..

[39]  Guangyao Li,et al.  Modeling and simulation of steady heat transfer analysis with material uncertainty: Generalized nth order perturbation isogeometric stochastic method , 2018, Numerical Heat Transfer, Part A: Applications.

[40]  Di Wu,et al.  Spectral stochastic isogeometric analysis of linear elasticity , 2018 .

[41]  Roland Wüchner,et al.  Isogeometric shell analysis with Kirchhoff–Love elements , 2009 .

[42]  Marcin Marek Kaminski,et al.  The Stochastic Perturbation Method for Computational Mechanics , 2013 .

[43]  Timon Rabczuk,et al.  A software framework for probabilistic sensitivity analysis for computationally expensive models , 2016, Adv. Eng. Softw..

[44]  P. Kerfriden,et al.  Shape optimization directly from CAD: An isogeometric boundary element approach using T-splines , 2017 .

[45]  Hung Nguyen-Xuan,et al.  Isogeometric Analysis of Laminated Composite Plates Using the Higher-Order Shear Deformation Theory , 2015 .

[46]  G. I. Schuëller,et al.  Buckling analysis of cylindrical shells with random geometric imperfections , 2003 .

[47]  George Stefanou,et al.  Stochastic finite element analysis of shells with combined random material and geometric properties , 2004 .

[48]  P. Kerfriden,et al.  Isogeometric boundary element methods for three dimensional static fracture and fatigue crack growth , 2017 .

[49]  Manolis Papadrakakis,et al.  Improving the computational efficiency in finite element analysis of shells with uncertain properties , 2005 .

[50]  Yongjie Zhang,et al.  Adaptive and Quality Quadrilateral/Hexahedral Meshing from Volumetric Data. , 2006, Computer methods in applied mechanics and engineering.

[51]  M. Scott,et al.  Acoustic isogeometric boundary element analysis , 2014 .

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

[53]  Gerhart I. Schuëller,et al.  Efficient stochastic structural analysis using Guyan reduction , 2011, Adv. Eng. Softw..

[54]  Bing Wang,et al.  Stochastic stable node-based smoothed finite element method for uncertainty and reliability analysis of thermo-mechanical problems , 2020 .

[55]  Hendrik Speleers,et al.  THB-splines: The truncated basis for hierarchical splines , 2012, Comput. Aided Geom. Des..

[56]  M. Kamiński,et al.  STOCHASTIC PERTURBATION APPROACH TO ENGINEERING STRUCTURE VIBRATIONS BY THE FINITE DIFFERENCE METHOD , 2002 .

[57]  G. Y. Li,et al.  Accurate analysis and thickness optimization of tailor rolled blanks based on isogeometric analysis , 2016 .

[58]  Alessandro Reali,et al.  GeoPDEs: A research tool for Isogeometric Analysis of PDEs , 2011, Adv. Eng. Softw..

[59]  Pierre Kerfriden,et al.  Adaptive Isogeometric analysis for plate vibrations: An efficient approach of local refinement based on hierarchical a posteriori error estimation , 2018, Computer Methods in Applied Mechanics and Engineering.

[60]  Di Wu,et al.  Spectral stochastic isogeometric analysis of free vibration , 2019, Computer Methods in Applied Mechanics and Engineering.

[61]  Hyuk-Chun Noh,et al.  Stochastic isogeometric analysis of free vibration of functionally graded plates considering material randomness , 2017 .

[62]  Xiangyang Cui,et al.  A Smoothed Finite Element Method (SFEM) for Linear and Geometrically Nonlinear Analysis of Plates and Shells , 2008 .

[63]  Marcin Kamiński.  Stochastic second-order BEM perturbation formulation , 1999 .

[64]  Yuri Bazilevs,et al.  Rotation free isogeometric thin shell analysis using PHT-splines , 2011 .

[65]  Roger Ghanem,et al.  Ingredients for a general purpose stochastic finite elements implementation , 1999 .

[66]  Guangyao Li,et al.  An isogeometric independent coefficients (IGA-IC) reduced order method for accurate and efficient transient nonlinear heat conduction analysis , 2018 .

[67]  Vinh Phu Nguyen,et al.  Isogeometric analysis: An overview and computer implementation aspects , 2012, Math. Comput. Simul..

[68]  Werner Wagner,et al.  Influence of stochastic geometric imperfections on the load-carrying behaviour of thin-walled structures using constrained random fields , 2018, Computational Mechanics.

[69]  Stéphane Bordas,et al.  Accelerating Monte Carlo estimation with derivatives of high-level finite element models , 2017 .

[70]  Hyuk-Chun Noh,et al.  A formulation for stochastic finite element analysis of plate structures with uncertain Poisson's ratio , 2004 .

[71]  T. Hughes,et al.  Finite Elements Based Upon Mindlin Plate Theory With Particular Reference to the Four-Node Bilinear Isoparametric Element , 1981 .

[72]  C. Y. Dong,et al.  An adaptive isogeometric boundary element method for predicting the effective thermal conductivity of steady state heterogeneity , 2018, Advances in Engineering Software.

[73]  O. C. Zienkiewicz,et al.  Reduced integration technique in general analysis of plates and shells , 1971 .

[74]  C. Soares,et al.  Spectral stochastic finite element analysis for laminated composite plates , 2008 .

[75]  H. Nguyen-Xuan,et al.  An extended isogeometric thin shell analysis based on Kirchhoff-Love theory , 2015 .

[76]  Thomas-Peter Fries,et al.  Fast Isogeometric Boundary Element Method based on Independent Field Approximation , 2014, ArXiv.

[77]  Sharif Rahman,et al.  A perturbation method for stochastic meshless analysis in elastostatics , 2001 .

[78]  H. Nguyen-Xuan,et al.  Isogeometric analysis using polynomial splines over hierarchical T-meshes for two-dimensional elastic solids , 2011 .

[79]  Zhenyu Liu,et al.  Isogeometric analysis of large thin shell structures based on weak coupling of substructures with unstructured T-splines patches , 2019, Adv. Eng. Softw..

[80]  Stéphane Bordas,et al.  Implementation of regularized isogeometric boundary element methods for gradient‐based shape optimization in two‐dimensional linear elasticity , 2016 .

[81]  Kumar K. Tamma,et al.  Model order reduction accelerated Monte Carlo stochastic isogeometric method for the analysis of structures with high-dimensional and independent material uncertainties , 2019, Computer Methods in Applied Mechanics and Engineering.

[82]  Guirong Liu,et al.  A generalized probabilistic edge-based smoothed finite element method for elastostatic analysis of Reissner–Mindlin plates , 2018 .

[83]  Benjamin Marussig,et al.  A Review of Trimming in Isogeometric Analysis: Challenges, Data Exchange and Simulation Aspects , 2017, Archives of Computational Methods in Engineering.

[84]  R. Echter,et al.  A hierarchic family of isogeometric shell finite elements , 2013 .