Engineering Analysis with Boundary Elements最新文献

{"title":"An alternative dual reciprocity BEM for P-SV wave propagation problems: A comparative study","authors":"Pouya Kavandi ,&nbsp;Mehdi Panji ,&nbsp;Navid Ganjian ,&nbsp;Jafar Asgari Marnani","doi":"10.1016/j.enganabound.2025.106238","DOIUrl":"10.1016/j.enganabound.2025.106238","url":null,"abstract":"<div><div>This research introduces a dual reciprocity boundary element method (BEM) designed to analyze the transient scattering of vertically travelling incident <em>P</em>-<em>SV</em> waves. By using static fundamental solutions and appropriate predictor operations, the domain inertia integrals from the equilibrium equation were transformed into boundary integral equations. The computable format of the integral equations was achieved by incorporating the effects of free-field displacements into the equations. After coding the formulation, the proposed technique's accuracy and efficiency were assessed through various wave scattering problems, such as an encased round hole, a half-circle depression, and a double-peaked mound, all exposed to vertically incoming <em>P</em>-<em>SV</em> waves. An extensive comparison was conducted with the full-plane time-domain boundary element technique available in the literature, focusing on accuracy and computation time. The findings indicated that despite the complexity and necessity for interior points in dual reciprocity models, it is more advantageous than the traditional full-plane time-domain approach, as it substantially reduces analysis time.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106238"},"PeriodicalIF":4.2,"publicationDate":"2025-03-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143735323","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"High-precision meshless method for 3D radiation diffusion problem in sphere and cylinder","authors":"Nan Ma,&nbsp;Qiuyan Xu,&nbsp;Zhiyong Liu,&nbsp;Jiye Yang","doi":"10.1016/j.enganabound.2025.106206","DOIUrl":"10.1016/j.enganabound.2025.106206","url":null,"abstract":"<div><div>The problem of radiation diffusion is extremely challenging due to the complex physical processes and nonlinear characteristics of the equation involved. In this paper, we propose a class of high-precision meshless methods for 3D nonlinear radiation diffusion equations applicable to spherical and cylindrical walls. Firstly, when the energy density is linearly related to temperature, we use a full-implicit difference scheme to discretize the time term, and then approximate the spatial term using radial basis functions to construct a new solution scheme for solving the 3D linear radiation diffusion equation. Secondly, when dealing with the nonlinear relationship between energy density and temperature, we successfully reduced the complexity of problem to be by linearizing <span><math><msup><mrow><mi>T</mi></mrow><mrow><mn>4</mn></mrow></msup></math></span>. Then, we use radial basis functions to approximate unknown functions and established a large class of solving schemes, which solved by the Kansa’s method. Finally, we validate the efficiency and high accuracy of the proposed methods through a series of numerical examples on spherical and cylindrical walls. In summary, the meshless numerical solution methods proposed in this paper not only avoids the complexity of meshing in irregular areas, but also provides a new and high-precision numerical solution method for the 3D radiation diffusion equation.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106206"},"PeriodicalIF":4.2,"publicationDate":"2025-03-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143735324","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"On the improvement of the local boundary conditions in GFEMgl","authors":"Túlio R.E. Marques ,&nbsp;Gabriela M. Fonseca ,&nbsp;Rafael M. Lins ,&nbsp;Felício B. Barros","doi":"10.1016/j.enganabound.2025.106219","DOIUrl":"10.1016/j.enganabound.2025.106219","url":null,"abstract":"<div><div>In this work, the ZZ-BD recovered stress field is first used to enhance the data transferred from the global to the local scale models in the Generalized Finite Element Method with Global–Local enrichments (GFEM^gl). The recovered stress field is constructed by solving a block-diagonal system of equations resulting from an <span><math><msub><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msub></math></span> approximate function projection associated with the singular stress field in the crack tip neighboring. In GFEM^gl analysis, the global solution is imposed as Dirichlet or Cauchy-type boundary conditions in the local domain. In the former case, only displacements are considered. The main contribution of this work lies in the definition of the Cauchy boundary conditions, where the stress field is combined with the displacements. A two-dimensional plate problem with an edge crack under mixed opening mode is solved using GFEM^gl. Stress intensity factors are extracted from global and local problems using the Interaction Integral strategy. Numerical results indicate that the Cauchy boundary conditions with the ZZ-BD recovered stress field provide a more accurate solution than raw or average stress fields, as well as regular Dirichlet boundary conditions. The effects of using a buffer zone in the local problem are also examined. Finally, the Interaction Integral performance strategy is investigated, with the key parameter being the circumference radius that intersects the elements where the stress intensity factors are extracted. An investigation is performed into the local and global problems, and a range of these parameters is identified to minimize errors in the stress intensity factors.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106219"},"PeriodicalIF":4.2,"publicationDate":"2025-03-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143726064","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Thermo-magneto-mechanical analysis of curved laminated structures with arbitrary variation of the material properties and novel recovery procedure","authors":"Francesco Tornabene ,&nbsp;Matteo Viscoti ,&nbsp;Rossana Dimitri ,&nbsp;Timon Rabczuk","doi":"10.1016/j.enganabound.2025.106232","DOIUrl":"10.1016/j.enganabound.2025.106232","url":null,"abstract":"<div><div>The paper introduces a novel methodology based on a generalized formulation and higher-order-theories for the fully-coupled multifield analysis of laminated curved structures subjected to thermal, magnetic, and mechanical loads. The formulation follows the Equivalent Single Layer approach, taking into account a generalized through-the-thickness expansion of displacement field components, scalar magnetic potential, and temperature variation with respect to the reference configuration. In addition, specific thickness functions are selected according to the Equivalent Layer Wise methodology, allowing the imposition of particular values of configuration variables in specific regions of the structure. The lamination scheme includes smart materials derived from an analytical homogenization technique, with material properties varying arbitrarily along the thickness direction within each layer. The fundamental relations are derived under thermodynamic equilibrium using curvilinear principal coordinates, and a semi-analytical Navier solution is derived for specific geometric, material, and loading conditions. A recovery procedure using Generalized Differential Quadrature is presented for reconstructing three-dimensional primary and secondary variables. In addition, a novel recovery procedure is presented for the first time, based on a Generalized Integral Quadrature. The model is validated through numerical examples involving straight and curved panels with various multifield load distributions, showing consistency and the computational efficiency when compared to three-dimensional reference solutions. New coupling effects between physical problems are explored, and parametric investigations highlight the influence of key governing parameters. Unlike the existing literature, this paper presents an efficient and accurate methodology for analyzing laminated smart structures of various curvatures with multifield couplings, not usually addressed by commercial software. This theory allows for arbitrary variations in multifield properties without using three-dimensional models that can be computationally expensive. In this way, novel possible design applications of smart materials and structures are offered in many engineering fields.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106232"},"PeriodicalIF":4.2,"publicationDate":"2025-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143726062","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A virtual material point peridynamic model for failure investigation of anisotropic laminated composites","authors":"Xiongwu Yang ,&nbsp;Dongsheng Mao ,&nbsp;Zhanhui Liu","doi":"10.1016/j.enganabound.2025.106236","DOIUrl":"10.1016/j.enganabound.2025.106236","url":null,"abstract":"<div><div>In this study, a new virtual material point peridynamic model (abbreviated as VMPPD) is proposed to capture the fracture behavior of composite laminates with arbitrary fiber orientation. The unique feature is that virtual material points serve as intermediate variables to achieve load transfer in a regularized discrete grid. As a result, a PD model for describing the reinforcement characteristics of composite materials was developed, and the limitation for these conventional PD models in characterizing fiber orientation was overcome. In the damage stage, the damage state and mechanical characteristics of composite materials are described by stiffness reduction technique to four material points rather than permanent termination to one material point. Furthermore, it is very convenient to introduce the concept of single-layer algorithms for multi-layer laminated structures under the VMPPD framework to achieve high efficiency. The tensile example of laminates has demonstrated that the VMPPD model can accurately predict the anisotropic characteristics of composite materials with a reliable numerical accuracy. It can also be observed from the matrix-dominated composite example that the fracture behavior of the laminate can be adaptively captured without any other numerical guidance techniques.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106236"},"PeriodicalIF":4.2,"publicationDate":"2025-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143726063","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Artificial neural network-based approach for simulating influenza dynamics: A nonlinear SVEIR model with spatial diffusion","authors":"Rahat Zarin","doi":"10.1016/j.enganabound.2025.106230","DOIUrl":"10.1016/j.enganabound.2025.106230","url":null,"abstract":"<div><div>Artificial Neural Networks (ANNs) have revolutionized machine learning by enabling systems to learn from data and generalize to new, unseen examples. As biologically inspired models, ANNs consist of interconnected neurons organized in layers, mimicking the human brain’s functioning. Their ability to model complex, nonlinear processes makes them powerful tools in various domains. In this study, the author apply ANNs to simulate the dynamics of a nonlinear Influenza transmission model with spatial diffusion. The model comprises five compartments: Susceptible (<span><math><mrow><mi>S</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow></math></span>), Vaccinated (<span><math><mrow><mi>V</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow></math></span>), Exposed (<span><math><mrow><mi>E</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow></math></span>), Infected (<span><math><mrow><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow></math></span>), and Recovered (<span><math><mrow><mi>R</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow></math></span>), governed by a system of partial differential equations (PDEs). We employ the Levenberg–Marquardt backpropagation algorithm to train the ANN, utilizing reference datasets generated through meshless and finite difference methods in MATLAB. The performance of the ANN is validated through mean square error (MSE) metrics, achieving a mean square error as low as <span><math><mrow><mn>1</mn><msup><mrow><mn>0</mn></mrow><mrow><mo>−</mo><mn>12</mn></mrow></msup></mrow></math></span>. Regression and state transition plots illustrate the training, testing, and validation processes. Furthermore, absolute error analyses across various components of the system confirm the robustness and accuracy of the proposed approach. The data were split into 81% for training, with 9% each for testing and validation.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106230"},"PeriodicalIF":4.2,"publicationDate":"2025-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143726065","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A novel fracture model for composite laminates based on bond-based peridynamics","authors":"Guanghui Zhang,&nbsp;Zili Dai","doi":"10.1016/j.enganabound.2025.106229","DOIUrl":"10.1016/j.enganabound.2025.106229","url":null,"abstract":"<div><div>Traditional model for composite laminae based on bond-based peridynamics (BB-PD) involves only two material parameters, which is insufficient to fully describe the complicated engineering properties of composite laminae. This limitation results in constrained Poisson's ratio and shear modulus in the PD model. In this study, a novel fracture model for composite laminae is proposed based on BB-PD, which includes the fiber bond, matrix bond, transverse bond, and tangential stiffness between particles, thereby overcoming the limitations of Poisson's ratio and shear modulus in the traditional BB-PD model. By employing the equivalence of strain energy density, expressions for the four micro-modulus parameters of the model are derived. Additionally, this study proposes a new surface correction method, based on the energy method, to correct surface effects in the model and reduce numerical errors. By stacking laminae, the fracture model is further extended into a 3D one for composite laminates. The energy criterion is then used to derive the critical value of micro elastic strain energy density, which can be utilized to evaluate the damage condition of composite materials. To validate the proposed model, this study simulates relevant experimental cases and analyzes the displacement results and fracture behavior.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106229"},"PeriodicalIF":4.2,"publicationDate":"2025-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143714396","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"RBF based backward differentiation methods for stiff differential equations","authors":"A. Sreedhar,&nbsp;Manoj Kumar Yadav,&nbsp;Chirala Satyanarayana","doi":"10.1016/j.enganabound.2025.106215","DOIUrl":"10.1016/j.enganabound.2025.106215","url":null,"abstract":"<div><div>Numerical solutions of initial value problems (IVPs) for stiff differential equations via explicit methods such as Euler’s method, trapezoidal method and Runge–Kutta methods suffer from stability issues and demand unacceptably small time steps. Backward differentiation formulas (BDF), a class of implicit methods, have been successfully used for resolving stiff IVPs. Classical BDF methods are derived using polynomial basis functions. In this paper, we develop radial basis function based finite difference (RBF-FD) type BDF methods for solving stiff problems. Therefore, we obtain analytical expressions for Gaussian and Multiquadric based RBF-BDF schemes along with their local truncation errors. Then we discuss the stability, order, consistency and convergence of RBF-BDF methods, which also depend on the free shape parameter. Finally, we validate the proposed methods by solving some benchmark problems. In order to gain enhanced accuracy, we adaptively choose the shape parameter such that local truncation error is minimized at each time-step. RBF-BDF methods of order two to six achieve at least one order greater accuracy and order of convergence than corresponding classical BDF schemes.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106215"},"PeriodicalIF":4.2,"publicationDate":"2025-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143704887","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Uncertainty quantification for the 3D half-space sound scattering problem of IGABEM based on the Catmull–Clark subdivision surfaces","authors":"Xiaohui Yuan ,&nbsp;Ruijin Huo ,&nbsp;Qingxiang Pei ,&nbsp;Gaochao Zhao ,&nbsp;Yongsong Li","doi":"10.1016/j.enganabound.2025.106222","DOIUrl":"10.1016/j.enganabound.2025.106222","url":null,"abstract":"<div><div>The generalized <span><math><mi>n</mi></math></span>th-order perturbation method for the quantitative uncertainty analysis in half-space acoustic problems proposed in this study is based on the isogeometric boundary element method, where the acoustic wave frequency is defined as a stochastic variable. We derive the Taylor series expansion and the kernel function formulation of the acoustic boundary integral equation for the half-space acoustic problem, and obtain the sound pressure’s <span><math><mi>n</mi></math></span>th-order derivative with respect to the acoustic wave frequency. In addition, we employ Burton–Miller method to deal with the fictitious frequency problem of external sound field and apply fast multipole method to accelerate the matrix–vector product computation. The statistical characterization of the acoustic state function is obtained based on the <span><math><mi>n</mi></math></span>th-order perturbation theory. Finally, the accuracy and efficacy of the uncertainty quantization algorithm is confirmed by three numerical examples.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106222"},"PeriodicalIF":4.2,"publicationDate":"2025-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143704888","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Correlation between acoustic emission characteristics and shear behavior of rock fracture","authors":"Yang Wu ,&nbsp;Zhihong Zhao ,&nbsp;Jinfan Chen ,&nbsp;Jintong Zhang ,&nbsp;Xingguang Zhao","doi":"10.1016/j.enganabound.2025.106234","DOIUrl":"10.1016/j.enganabound.2025.106234","url":null,"abstract":"<div><div>Rock fractures significantly diminish shear strength and stability of rock masses. Understanding the shear behavior of fractured rock and associated energy release is essential for disaster prediction in rock engineering. This study investigates the shear behavior and damage evolution of intact and fractured rock samples through the analysis of acoustic emission (AE) characteristics. We conduct a series of direct shear tests and discrete element numerical simulations calibrated by experimental results. AE signals are simultaneously monitored during shear, and AE simulations are performed based on moment tensor inversion theory. We propose an effective method to determine crack damage stress, defined as the maximum gradient point of cumulative AE energy-shear displacement curve in the pre-peak stage. The results demonstrate that the distribution of explosive AE events is relatively concentrated along the fracture profile or failure zone, while implosive and shear events are principally located in the contacting asperities or fractures newly induced by shear. Meanwhile, a quantified relationship between shear parameters and AE energy has been established to assess shear properties and predict energy release. Its feasibility is validated by experimental results. The study contributes to providing a basis for analyzing rock failure in engineering and early warnings for rock mass disasters.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106234"},"PeriodicalIF":4.2,"publicationDate":"2025-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143714395","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}