The Boundary Element Method for Ordinary State-Based Peridynamics

The peridynamics(PD),as a promising nonlocal continuum mechanics theory,shines in solving discontinuous problems.Up to now,various numerical methods,such as the peridynamic mesh-free particlemethod(PD-MPM),peridynamic finite element method(PD-FEM),and peridynamic boundary element method(PD-BEM),have been proposed.PD-BEM,in particular,outperforms other methods by eliminating spurious boundary softening,efficiently handling infinite problems,and ensuring high computational accuracy.However,the existing PD-BEM is constructed exclusively for bond-based peridynamics(BBPD)with fixed Poisson’s ratio,limiting its applicability to crack propagation problems and scenarios involving infinite or semi-infinite problems.In this paper,we address these limitations by introducing the boundary element method(BEM)for ordinary state-based peridynamics(OSPD-BEM).Additionally,we present a crack propagationmodel embeddedwithin the framework ofOSPD-BEM to simulate crack propagations.To validate the effectiveness of OSPD-BEM,we conduct four numerical examples:deformation under uniaxial loading,crack initiation in a double-notched specimen,wedge-splitting test,and threepoint bending test.The results demonstrate the accuracy and efficiency of OSPD-BEM,highlighting its capability to successfully eliminate spurious boundary softening phenomena under varying Poisson’s ratios.Moreover,OSPDBEMsignificantly reduces computational time and exhibits greater consistencywith experimental results compared to PD-MPM.

Topology Optimization of Sound-Absorbing Materials for Two-Dimensional Acoustic Problems Using Isogeometric Boundary Element Method

In this work,an acoustic topology optimizationmethod for structural surface design covered by porous materials is proposed.The analysis of acoustic problems is performed using the isogeometric boundary elementmethod.Taking the element density of porousmaterials as the design variable,the volume of porousmaterials as the constraint,and the minimum sound pressure or maximum scattered sound power as the design goal,the topology optimization is carried out by solid isotropic material with penalization(SIMP)method.To get a limpid 0–1 distribution,a smoothing Heaviside-like function is proposed.To obtain the gradient value of the objective function,a sensitivity analysis method based on the adjoint variable method(AVM)is proposed.To find the optimal solution,the optimization problems are solved by the method of moving asymptotes(MMA)based on gradient information.Numerical examples verify the effectiveness of the proposed topology optimization method in the optimization process of two-dimensional acoustic problems.Furthermore,the optimal distribution of sound-absorbingmaterials is highly frequency-dependent and usually needs to be performed within a frequency band.

Error Analysis of A New Higher Order Boundary Element Method for A Uniform Flow Passing Cylinders

A higher order boundary element method(HOBEM)is presented for inviscid flow passing cylinders in bounded or unbounded domain.The traditional boundary integral equation is established with respect to the velocity potential and its normal derivative.In present work,a new integral equation is derived for the tangential velocity.The boundary is discretized into higher order elements to ensure the continuity of slope at the element nodes.The velocity potential is also expanded with higher order shape functions,in which the unknown coefficients involve the tangential velocity.The expansion then ensures the continuities of the velocity and the slope of the boundary at element nodes.Through extensive comparison of the results for the analytical solution of cylinders,it is shown that the present HOBEM is much more accurate than the conventional BEM.

Three-dimensional acoustic propagation model for shallow waters based on an indirect boundary element method

This work has a two-fold purpose.On the one hand,the theoretical formulation of a three-dimensional(3D)acoustic propagation model for shallow waters with a constant sound speed is presented,based on the boundary element method(BEM),which uses a half-space Green function instead of the more conventional free-space Green function.On the other hand,a numerical implementation is illustrated to explore the formulation in simple idealized cases,controlled by a few parameters,which provides necessary tests for the accuracy and performance of the model.The half-space Green's function,which has been previously used in scattering and diffraction,adds terms to the usual expressions of the integral operators without altering their continuity properties.Verifications against the wavenumber integration solution of the Pekeris waveguide suggest that the model allows an adequate prediction for the acoustic field.Likewise,numerical experiments in relation to the necessary mesh size for the description of the water-marine sediment interface lead to the conclusion that a transmission loss prediction with acceptable accuracy can be obtained with the use of a limited mesh around the desired evaluation region.

Three-Dimensional Boundary Element Method Applied to Nonlinear Wave Transformation

For higher accuracy in simulating the transformation of three dimensional waves, in consideration of the advantages of constant panels and linear elements, a combined boundary elements is applied in this research. The method can be used to remove the transverse vibration due to the accumulation of computational errors. A combined boundary condition of sponge layer and Sommerfeld radiation condition is used to remove the reflected waves from the computing domain. By following the water particle on the water surface, the third order Stokes wave transform is simulated by the numerical wave flume technique. The computed results are in good agreement with theoretical ones.

BOUNDARY ELEMENT METHOD FOR MOVING AND ROLLING CONTACT OF 2D ELASTIC BODIES WITH DEFECTS

A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes.

Monte Carlo Simulation of Fractures Using Isogeometric Boundary Element Methods Based on POD-RBF

This paper presents a novel framework for stochastic analysis of linear elastic fracture problems.Monte Carlo simulation(MCs)is adopted to address the multi-dimensional uncertainties,whose computation cost is reduced by combination of Proper Orthogonal Decomposition(POD)and the Radial Basis Function(RBF).In order to avoid re-meshing and retain the geometric exactness,isogeometric boundary element method(IGABEM)is employed for simulation,in which the Non-Uniform Rational B-splines(NURBS)are employed for representing the crack surfaces and discretizing dual boundary integral equations.The stress intensity factors(SIFs)are extracted by M integral method.The numerical examples simulate several cracked structures with various uncertain parameters such as load effects,materials,geometric dimensions,and the results are verified by comparison with the analytical solutions.

Introducing an effective coherence function to generate non-uniform ground motion on topographic site using time-domain boundary element method

In this study,a comprehensive parametric analysis was performed on non-uniform excitation of V-shaped topography using the boundary element method in time domain.For this purpose,wave scattering analysis was carried out on a topography subjected to the SV-wave for different predominant frequencies and shape ratios.Based on the numerical results,new coherence and time delay functions are proposed to generate non-uniform ground motion for topographic irregularities.The efficiency and accuracy of the proposed functions for real engineering problems are indicated by comparison with observations reported in previous literature.

A Combined Shape and Topology Optimization Based on Isogeometric Boundary Element Method for 3D Acoustics

A combined shape and topology optimization algorithm based on isogeometric boundary element method for 3D acoustics is developed in this study.The key treatment involves using adjoint variable method in shape sensitivity analysis with respect to non-uniform rational basis splines control points,and in topology sensitivity analysis with respect to the artificial densities of sound absorption material.OpenMP tool in Fortran code is adopted to improve the efficiency of analysis.To consider the features and efficiencies of the two types of optimization methods,this study adopts a combined iteration scheme for the optimization process to investigate the simultaneous change of geometry shape and distribution of material to achieve better noise control.Numerical examples,such as sound barrier,simple tank,and BeTSSi submarine,are performed to validate the advantage of combined optimization in noise reduction,and to demonstrate the potential of the proposed method for engineering problems.

Machine Learning Enhanced Boundary Element Method:Prediction of Gaussian Quadrature Points

This paper applies a machine learning technique to find a general and efficient numerical integration scheme for boundary element methods.A model based on the neural network multi-classification algorithmis constructed to find the minimum number of Gaussian quadrature points satisfying the given accuracy.The constructed model is trained by using a large amount of data calculated in the traditional boundary element method and the optimal network architecture is selected.The two-dimensional potential problem of a circular structure is tested and analyzed based on the determined model,and the accuracy of the model is about 90%.Finally,by incorporating the predicted Gaussian quadrature points into the boundary element analysis,we find that the numerical solution and the analytical solution are in good agreement,which verifies the robustness of the proposed method.

The Boundary Element Method for Pile-Soil Interaction

In this paper, several mathmatical models for the pile- soil interaction are outlined. The Boundary Element Method is one of the very effective methods for the reasonable models of elasticity and elastoplasticity. The major of this paper is concerned with the Boundary Element Method for the pile-soil interaction, including general methods and calculating formulation of static and dynamic analysis of the pile and pile groups. Some results of analysis are also given.

STOCHASTIC BOUNDARY ELEMENT METHOD IN ELASTICITY

The stochastic boundary element method is developed to analyze elasticity problems with random material and/or geometrical parameters and ran- domly perturbed boundaries. Based on the first-order Taylor series expansion, the boundary integration equations concerning the mean and deviation of the displace- ments are derived, respectively. It is found that the randomness of material param- eters is equivalent to a random body force, so the mean and covariance matrices of unknown boundary displacements and tractions can be obtained. Furthermore, the mean and covariance of displacements and stresses at inner points can also be obtained. Numerical examples show that the proposed stochastic boundary element method gives satisfactory solutions, as compared with those obtained by theoretical analysis or other numerical methods.

Ship Hull Slamming Analysis with Nonlinear Boundary Element Method

A boundary element method is developed for calculating the flare ship hull slamming problem. The nonlinear free surface elevation and the linear element assumption are employed. The method has been verified by comparisons with results for the water entry of wedges with various deadrise angles. Numerical results show that the pressure distribution varies greatly with the ship hull with different curvilinear equations, and the slamming features are also different; From the numerical simulation, the authors found that the structural damage of the flare hull might be caused by the increasing hydrodynamic pressure over an extensive area on the flare when the upper part of the flare comes into contact with water.

Resolving Domain Integral Issues in Isogeometric Boundary Element Methods via Radial Integration:A Study of Thermoelastic Analysis

The paper applied the isogeometric boundary element method(IGABEM)to thermoelastic problems.The Non-Uniform Rational B-splines(NURBS)used to construct geometric models are employed to discretize the boundary integral formulation of the governing equation.Due to the existence of thermal stress,the domain integral term appears in the boundary integral equation.We resolve this problem by incorporating radial integration method into IGABEM which converts the domain integral to the boundary integral.In this way,IGABEM can maintain its advantages in dimensionality reduction and more importantly,seamless integration of CAD and numerical analysis based on boundary representation.The algorithm is verified by numerical examples.

A new magneto-cardiogram study using a vector model with a virtual heart and the boundary element method

A cardiac vector model is presented and verified, and then the forward problem for cardiac magnetic fields and electric potential are discussed based on this model and the realistic human torso volume conductor model, including lungs. A torso-cardiac vector model is used for a 12-lead electrocardiographic （ECG） and magneto-cardiogram （MCG） simulation study by using the boundary element method （BEM）. Also, we obtain the MCG wave picture using a compound four-channel HTc.SQUID system in a magnetically shielded room. By comparing the simulated results and experimental results, we verify the cardiac vector model and then do a preliminary study of the forward problem of MCG and ECG. Therefore, the results show that the vector model is reasonable in cardiac electrophysiology.

Forward and inverse problem for cardiac magnetic field and electric potential using two boundary element methods

This paper discusses the forward and inverse problem for cardiac magnetic fields and electric potentials. A torso-heart model established by boundary element method （BEM） is used for studying the distributions of cardiac magnetic fields and electric potentials. Because node-to-node and triangle-to-triangle BEM can lead to discrepant field distributions, their properties and influences are compared. Then based on constructed torso-heart model and supposed current source functional model-current dipole array, the magnetic and electric imaging by optimal constrained linear inverse method are applied at the same time. Through figure and reconstructing parameter comparison, though the magnetic current dipole array imaging possesses better reconstructing effect, however node-to-node BEM and triangleto-triangle BEM make little difference to magnetic and electric imaging.

Weakly Singular Symmetric Galerkin Boundary Element Method for Fracture Analysis of Three-Dimensional Structures Considering Rotational Inertia and Gravitational Forces

The Symmetric Galerkin Boundary Element Method is advantageous for the linear elastic fracture and crackgrowth analysis of solid structures,because only boundary and crack-surface elements are needed.However,for engineering structures subjected to body forces such as rotational inertia and gravitational loads,additional domain integral terms in the Galerkin boundary integral equation will necessitate meshing of the interior of the domain.In this study,weakly-singular SGBEM for fracture analysis of three-dimensional structures considering rotational inertia and gravitational forces are developed.By using divergence theorem or alternatively the radial integration method,the domain integral terms caused by body forces are transformed into boundary integrals.And due to the weak singularity of the formulated boundary integral equations,a simple Gauss-Legendre quadrature with a few integral points is sufficient for numerically evaluating the SGBEM equations.Some numerical examples are presented to verify this approach and results are compared with benchmark solutions.

Diagonal form fast multipole boundary element method for 2D acoustic problems based on Burton-Miller boundary integral equation formulation and its applications

This paper describes formulation and implementation of the fast multipole boundary element method （FMBEM） for 2D acoustic problems. The kernel function expansion theory is summarized, and four building blocks of the FMBEM are described in details. They are moment calculation, moment to moment translation, moment to local translation, and local to local translation. A data structure for the quad-tree construction is proposed which can facilitate implementation. An analytical moment expression is derived, which is more accurate, stable, and efficient than direct numerical computation. Numerical examples are presented to demonstrate the accuracy and efficiency of the FMBEM, and radiation of a 2D vibration rail mode is simulated using the FMBEM.

Subdivision Surface-Based Isogeometric Boundary Element Method for Steady Heat Conduction Problems with Variable Coefficient

The present work couples isogeometric analysis(IGA)and boundary element methods(BEM)for three dimensional steady heat conduction problems with variable coefficients.The Computer-Aided Design(CAD)geometries are built by subdivision surfaces,and meantime the basis functions of subdivision surfaces are employed to discretize the boundary integral equations for heat conduction analysis.Moreover,the radial integration method is adopted to transform the additional domain integrals caused by variable coefficients to the boundary integrals.Several numerical examples are provided to demonstrate the correctness and advantages of the proposed algorithm in the integration of CAD and numerical analysis.

ELECTRO-MECHANICAL COUPLING ANALYSIS OF MEMS STRUCTURES BY BOUNDARY ELEMENT METHOD

In this paper,we present the applications of Boundary Element Method(BEM) to simulate the electro-mechanical coupling responses of Micro-Electro-Mechanical systems(MEMS). The algorithm is programmed in our research group based on BEM modeling for electrostatics and elastostatics.Good agreement is shown while the simulation results of the pull-in voltages are compared with the theoretical/experimental ones for some examples.