基于比例边界有限元法计算应力强度因子的不确定量化分析

PDF(2907 KB)

振动与冲击 ›› 2024, Vol. 43 ›› Issue (5) : 250-259.

论文

胡昊文1，2，陈灯红1，2，王乾峰1，2，胡记磊1，2，骆欢1，2

作者信息 +

UQ analysis in calculating SIFs based on SBFEM

HU Haowen1,2,CHEN Denghong1,2,WANG Qianfeng1,2,HU Jilei1,2,LUO Huan1,2

Author information +

文章历史 +

摘要

应力强度因子（SIFs）是预测荷载作用下结构中裂纹产生和扩展的重要参数。半解析的比例边界有限元法（SBFEM）结合了有限元和边界元法的优势，在裂纹尖端或存在奇异应力的区域不需要局部网格细化，可以直接提取应力强度因子。本文在比例边界有限元法计算应力强度因子的框架下，引入随机参数进行蒙特卡罗模拟(MCs)，并提出一种新颖的基于蒙特卡罗模拟的不确定量化（UQ）分析。与直接的蒙特卡罗模拟不同，采用奇异值分解（SVD）构造低阶的子空间，降低系统的自由度，并使用径向基函数（RBF）对子空间进行近似，通过子空间的线性组合获得新的结构响应，实现基于MCs的快速不确定量化分析。考虑不同荷载状况下，结构形状参数和材料属性参数对应力强度因子的影响，使用改进的蒙特卡罗模拟计算应力强度因子的统计特征，量化不确定参数对结构的影响。最后通过若干算例验证了本文算法的准确性和有效性。

Abstract

Stress intensity factors (SIFs) are important parameters in predicting the initiation and propagation of cracks in structures subjected to loads. The semi-analytical scaled boundary finite element method (SBFEM) combines the advantages of finite element and boundary element method, and without local mesh refinement around crack tips or regions existing stress singularities. The stress intensity factors can be directly extracted. Within the framework of the scaled boundary finite element method in evaluating stress intensity factors, random parameters are introduced for Monte Carlo simulation (MCs), a novel uncertainty quantification (UQ) analysis based on Monte Carlo simulation is proposed in this manuscript. Unlike direct Monte Carlo simulation, singular value decomposition (SVD) is used to construct lower order subspaces to reduce the degree of freedom of the system. Radial basis function (RBF) is used to approximate the subspaces and obtain new structural responses through linear combinations of subspaces, achieving fast uncertainty quantification analysis based on MCs. Considering the effects of structural shape parameters and material property parameters on stress intensity factors under different load conditions, an improved Monte Carlo simulation is used to calculate the statistical characteristics of stress intensity factors and quantify the impact of uncertain parameters on the structure. Finally, the accuracy and effectiveness of the algorithm proposed in this paper are verified by several examples.

导出引用

胡昊文1，2，陈灯红1，2，王乾峰1，2，胡记磊1，2，骆欢1，2. 基于比例边界有限元法计算应力强度因子的不确定量化分析[J]. 振动与冲击, 2024, 43(5): 250-259

HU Haowen1,2,CHEN Denghong1,2,WANG Qianfeng1,2,HU Jilei1,2,LUO Huan1,2. UQ analysis in calculating SIFs based on SBFEM[J]. Journal of Vibration and Shock, 2024, 43(5): 250-259

参考文献

[1] Li C, Niu Z, Hu Z, et al. Effectiveness of the stress solutions in notch/crack tip regions by using extended boundary element method[J]. Engineering Analysis with Boundary Elements, 2019, 108(1): 1-13. [2] Peng X, Atroshchenko E, Kerfriden P, et al. Isogeometric boundary element methods for three dimensional static fracture and fatigue crack growth[J]. Computer Methods in Applied Mechanics & Engineering, 2017, 316(1): 151-185. [3] Song C, Wolf J P. The scaled boundary finite-element method—alias consistent infinitesimal finite-element cell method—for elastodynamics[J]. Computer Methods in Applied Mechanics & Engineering, 1997, 147(3-4): 329-355. [4] Song C. The Scaled Boundary Finite Element Method:Introduction to Theory and Implementation[M]. Sydney:John Wiley & Sons, 2018. [5] Wolf J P, Song C. Finite element modelling of unbounded media[M]. Chichester: John Wiley & Sons, 1996. [6] Bazyar M H, Song C. A continued-fraction-based high-order transmitting boundary for wave propagation in unbounded domains of arbitrary geometry[J]. International Journal for Numerical Methods in Engineering, 2008, 74(2): 209-237. [7] Qu Y, Chen D, Liu L, et al. A direct time-domain procedure for the seismic analysis of dam–foundation–reservoir systems using the scaled boundary finite element method[J]. Computers and Geotechnics, 2021, 138: 104364 [8] 高毅超，梅真. 基于ABAQUS的动水压力波双渐近透射边界单元及应用[J]. 振动与冲击, 2019, 38(10): 38-42. GAO Yi-chao, MEI Zhen. A doubly asymptotic open boundary element for hydrodynamic pressure wave analysis based on the ABAQUS and its application[J]. Journal of Vibration and Shock, 2019, 38(10): 38-42. [9] 钟红, 暴艳利, 林皋. 基于多边形比例边界有限元的重力坝裂缝扩展过程模拟[J]. 水利学报, 2014, 45(S1): 30-37. ZHONG Hong, BAO Yan-li, LIN Gao. Modelling of crack propagation of gravity dams based on scaled boundary polygons[J]. Journal of Hydraulic Engineering, 2014, 45(S1): 30-37. [10] 施明光, 徐艳杰, 张楚汉,等. 比例边界有限元模拟裂纹和夹杂动力相互作用[J]. 振动与冲击, 2016, 35(4): 15-21. SHI Ming-guang, XU Yan-jie, ZHANG Chu-han, et al. Simulation of dynamic interactions between a crack and inclusions with scaled boundary finite element method[J]. Journal of Vibration and Shock, 2016, 35(4): 15-21. [11] 章鹏,杜成斌,江守燕. 比例边界有限元法求解裂纹面接触问题[J].力学学报, 2017, 49(6): 1335-1347. ZHANG Peng, DU Cheng-bin, JIANG Shou-yan. Crack face contact problem analysis using scaled boundary finite element method[J]. Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(6): 1335-1347. [12] Yu B, Sun W, Wei P, et al. The pixel-based quadtree SBFEM with the parameter level set method for identifying cracks and voids[J]. Computational Mechanics, 2022, 70(5): 911–929. [13] Chen D, Birk C, Song C, et al. A high-order approach for modelling transient wave propagation problems using the scaled boundary finite element method[J]. International Journal for Numerical Methods in Engineering, 2014, 97(13): 937–959. [14] Chen D, Wang M, Dai S. A high-order doubly asymptotic continued-fraction solution for frequency-domain analysis of vector wave propagation in unbounded domains[J]. Soil Dynamics and Earthquake Engineering, 2018, 113: 230–240. [15] Chen C, Lin G, Hu Z. An innovative and efficient solution for axisymmetric contact problem between structure and half-space[J]. Engineering analysis with boundary elements, 2022, 142: 10–27. [16] 张俊奇，余波，宋崇民. 比例边界有限元法的最新研究进展[J]. 应用力学学报, 2022, 39(6): 1038-1054． ZHANG Jun-qi, YU Bo, SONG Chong-min．The latest progress in research on the scaled boundary finite element method[J]．Chinese journal of Applied Mechanics, 2022, 39(6): 1038-1054． [17] Song C, Tin-Loi F, Gao W. A definition and evaluation procedure of generalized stress intensity factors at cracks and multi-material wedges[J]. Engineering Fracture Mechanics, 2010, 77(12): 2316–2336. [18] Song C, Wolf J P. Semi-analytical representation of stress singularities as occurring in cracks in anisotropic multi-materials with the scaled boundary finite-element method[J]. Computers& Structures, 2002, 80(2): 183–197. [19] Song C, Vrcelj Z. Evaluation of dynamic stress intensity factors and t-stress using the scaled boundary finite-element method[J]. Engineering Fracture Mechanics, 2008, 75(8): 1960–1980. [20] Song C. Evaluation of power-logarithmic singularities, t-stresses and higher order terms of in-plane singular stress fields at cracks and multi-material corners[J]. Engineering Fracture Mechanics, 2005, 72(10): 1498–1530. [21] 汤涛, 周涛. 不确定性量化的高精度数值方法和理论[J].中国科学：数学, 2015, 45(7): 891-928. TANG Tao, ZHOU Tao. Recent developments in high order numerical methods for uncertainty quantification[J]. Scientia Sinica(mathematica), 2015, 45(7): 891-928. [22] Hurtado J E, Barbat A H. Monte carlo techniques in computational stochastic mechanics[J]. Archives of Computational Methods in Engineering, 1998, 5(1): 3–30. [23] Chowdhury M S, Song C, Gao W. Probabilistic fracture mechanics by using monte carlo simulation and the scaled boundary finite element method[J]. Engineering Fracture Mechanics, 2011, 78(12): 2369–2389. [24] 谷晓忱. 并行蒙特卡罗计算硬件加速器的关键技术研究[D].长沙:国防科技大学, 2012. GU Xiao-chen. Research on hardware accelerator for parallel Monte Carlo simulation[D]. Changsha: National University of Defense Technology, 2012. [25] Fishman G S. Monte Carlo: Concepts, Algorithms, and Applications[M]. New York: Springer-Verlag, 1996. [26] Ding C, Deokar R R, Ding Y, et al. Model order reduction accelerated monte carlo stochastic isogeometric method: Analyzing structures considering high-dimensional and independent material uncertainties[J]. Computer Methods in Applied Mechanics & Engineering, 2019, 349(1): 266–284. [27] Song C. The scaled boundary finite element method in structural dynamics[J]. International Journal for Numerical Methods in Engineering, 2009, 77(8): 1139–1171. [28] Li S, Trevelyan J, Wu Z, et al. An adaptive SVD-Krylov reduced order model for surrogate based structural shape optimization through isogeometric boundary element method[J]. Computer Methods in Applied Mechanics & Engineering, 2019, 349(1): 312–338. [29] Wu Z, Wang D, Okolo N P, et al. Unified estimate of gaussian kernel width for surrogate models[J]. Neurocomputing, 2016, 203(26): 41–51. [30] Smith D J, Ayatollahi M R, Pavier M J. The role of t-stress in brittle fracture for linear elastic materials under mixed-mode loading [J]. Fatigue Fracture of Engineering Materials Structures, 2001, 24(2): 137–150. [31] Ding C, Hu X, Cui X, et al. Isogeometric generalized nth order perturbation-based stochastic method for exact geometric modeling of (composite) structures: Static and dynamic analysis with random material parameters[J]. Computer Methods in Applied Mechanics &Engineering, 2019, 346(4): 1002–1024. [32] Stazi F L , Budyn E , Chessa J ,et al. An extended finite element method with higher-order elements for curved cracks[J].Computational Mechanics, 2003, 31(1-2): 38-48. [33] Chen D, Dai S. Dynamic fracture analysis of the soil-structure interaction system using the scaled boundary finite element method[J]. Engineering Analysis with Boundary Elements,2017,77(4):26–35.