Generalized nth-Order Perturbation Method Based on Loop Subdivision Surface Boundary Element Method for Three-Dimensional Broadband Structural Acoustic Uncertainty Analysis

In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Miller method is employed to solve the problem of non-unique solutions that may be encountered in the external acoustic field,and the nth-order discretization formulation of the boundary integral equation is derived.In addition,the computation of loop subdivision surfaces and the subdivision rules are introduced.In order to confirm the effectiveness of the algorithm,the computed results are contrasted and analyzed with the results under Monte Carlo simulations(MCs)through several numerical examples.

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.

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.

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.

Elastic one-return boundary element method and hybrid elastic thin-slab propagator for strong-contrast media with sharp boundaries

In this paper, we developed the theory and algorithm of an elastic one-way boundary element method（BEM） and a corresponding hybrid elastic thin-slab propagator for earth media with sharp boundaries between strong contrast media. This approach can takes the advantage of accurate boundary condition of BEM and completely overcomes the weak contrast limitation of the perturbationtheory based one-way operator approach. The one-way BEM is a smooth boundary approximation, which avoids huge matrix operations in exact full BEM. In addition, the one-way BEM can model the primary-only transmitted and reflected waves and therefore is a valuable tool in elastic imaging and inversion. Through numerical tests for some simple models,we proved the validity and efficiency of the proposed method.

SEMI-ANALYTICAL FINITE ELEMENT METHOD FOR FICTITIOUS CRACK MODEL IN FRACTURE MECHANICS OF CONCRETE

Based on the Hamiltonian governing equations of plane elasticity for sectorial domain, the variable separation and eigenfunction expansion techniques were employed to develop a novel analytical finite element for the fictitious crack model in fracture mechanics of concrete. The new analytical element can be implemented into FEM program systems to solve fictitious crack propagation problems for concrete cracked plates with arbitrary shapes and loads. Numerical results indicate that the method is more efficient and accurate than ordinary finite element method.

Substructure Method of Three-Dimensional Semi-Analytic Boundary Element and Its Application to Seismic Scattering of Hill Topography

In this paper,a substructure method of three-dimensional semi-analytic boundary element is established.The seismic scattering by three-dimensional topography of a hill can be analyzed by the method in frequency domain.Using this method,the computational effort and storage space are reduced considerably.Finally,analytic results are given.

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.

HIGH ACCURACY FINITE VOLUME ELEMENT METHOD FOR TWO-POINT BOUNDARY VALUE PROBLEM OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS

In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.

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.

An efficient source wavefield reconstruction scheme using single boundary layer values for the spectral element method

In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.

Increment-Dimensional Scaled Boundary Finite Element Method for Solving Transient Heat Conduction Problem

An increment-dimensional scaled boundary finite element method (ID-SBFEM) is developed to solve the transient temperature field.To improve the accuracy of SBFEM,the effect of high frequency factor on dynamic stiffness is considered,and the first-order continued fraction technique is used.After the derivation,the SBFE equations are obtained,and the dimensions of thermal conduction,the thermal capacity matrix and the vector of the right side term in the equations are doubled.An example is presented to illustrate the feasibility and good accuracy of the proposed method.

Isogeometric Boundary Element Analysis for 2D Transient Heat Conduction Problem with Radial Integration Method

This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structures,are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations.Bezier extraction technique is employed to accelerate the evaluation of NURBS basis functions.We adopt a radial integration method to address the additional domain integrals.The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis.

CALCULATION OF MILL RIGIDITY BY THREE DIMENSION CONTACT BOUNDARY ELEMENT METHOD

Vertical rigidity of the space self adaptive 530 high rigidity mill is calculated by applying the boundary element method (BEM) of three dimension elastic contact problem,which can update the existed deforming separation calculating theory and corresponding methods of material mechanics,elastic mechanics and finite element method.The method has less hypotheses and stronger synthesis in contact type calculating model.The advantages of the method are high calculating rate,high calculating accuracy,etc..

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.

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.

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.

Modeling of Fully Nonlinear Wave Radiation by Submerged Moving Structures Using the Higher Order Boundary Element Method

The higher-order boundary element method is applied to the numerical simulation of nonlinear waves radiated by a forced oscillating fully submerged vertical circular cylinder. In this time-domain approach, the mixed boundary value problem based on an Eulerian description at each time step is solved using the higher order boundary element method. The 4th-order Runge–Kutta scheme is adopted to update the free water surface boundary conditions expressed in a Lagrangian formulation. Following completion of the numerical model, the problems of radiation(heave) of water waves by a submerged sphere in finite depth are simulated and the computed results are verified against the published numerical results in order to ensure the effectiveness of the model. The validated numerical model is then applied to simulate the nonlinear wave radiation by a fully submerged vertical circular cylinder undergoing various forced sinusoidal motion in otherwise still conditions. The numerical results are obtained for a series of wave radiation problems; the completely submerged cylinder is placed in surging, heaving and combined heave-pitching motions with different drafts, amplitudes and frequencies. The corresponding numerical results of the cylinder motions, wave profiles, and hydrodynamic forces are then compared and explained for all the cases.