A STUDY ON THE WEIGHT FUNCTION OF THE MOVING LEAST SQUARE APPROXIMATION IN THE LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.

RESEARCH ON THE COMPANION SOLUTION FOR A THIN PLATE IN THE MESHLESS LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless local boundary integral equation method is a currently developed numerical method, which combines the advantageous features of Galerkin finite element method(GFEM), boundary element method(BEM) and element free Galerkin method(EFGM), and is a truly meshless method possessing wide prospects in engineering applications. The companion solution and all the other formulas required in the meshless local boundary integral equation for a thin plate were presented, in order to make this method apply to solve the thin plate problem.

Improved non-singular local boundary integral equation method

When the source nodes are on the global boundary in the implementation of local boundary integral equation method （LBIEM）, singularities in the local boundary integrals need to be treated specially. In the current paper, local integral equations are adopted for the nodes inside the domain and moving least square approximation （MLSA） for the nodes on the global boundary, thus singularities will not occur in the new al- gorithm. At the same time, approximation errors of boundary integrals are reduced significantly. As applications and numerical tests, Laplace equation and Helmholtz equation problems are considered and excellent numerical results are obtained. Furthermore, when solving the Helmholtz problems, the modified basis functions with wave solutions are adapted to replace the usually-used monomial basis functions. Numerical results show that this treatment is simple and effective and its application is promising in solutions for the wave propagation problem with high wave number.

Roll Flattening Analytical Model in Flat Rolling by Boundary Integral Equation Method

In order to improve rolled strip quality, precise plate shape control theory should be established. Roll flat- tening theory is an important part of the plate shape theory. To improve the accuracy of roll flattening calculation based on semi infinite body model, especially near the two roll barrel edges, a new and more accurate roll flattening model is proposed. Based on boundary integral equation method, an analytical model for solving a finite length semi infinite body is established. The lateral surface displacement field of the finite length semi-infinite body is simulated by finite element method （FEM） and lateral surface displacement decay functions are established. Based on the boundary integral equation method, the numerical solution of the finite length semi-infinite body under the distribu ted force is obtained and an accurate roll flattening model is established. Different from the traditional semi-infinite body model, the matrix form of the new roll flattening model is established through the mathematical derivation. The result from the new model is more consistent with that by FEM especially near the edges.

Fast Multipole Accelerated Boundary Integral Equation Method for Evaluating the Stress Field Associated with Dislocations in a Finite Medium

In this paper,we develop an efficient numericalmethod based on the boundary integral equation formulation and new version of fast multipole method to solve the boundary value problem for the stress field associated with dislocations in a finite medium.Numerical examples are presented to examine the influence from material boundaries on dislocations.

Interaction of Linear Waves with Infinitely Long Horizontal Cylinders Studied by Boundary Element Method

The two-dimensional problems concerning the interaction of linear water waves with cylinders of arbitrary shape in two-layer deep water are investigated by use of the Boundary Integral Equation Method （BIEM）. Simpler new expressions for the Green functions are derived, and verified by comparison of results obtained by BIEM with those by an analytical method. Examined are the radiation and scattering of linear waves by two typical configurations of cylinders in two- layer deep water. Hydrodynamic behaviors including hydrodynamic coefficients, wave forces, reflection and transmission coef- ficients and energies are analyzed in detail, and some imeresting physical phenomena are observed.

A moving Kriging interpolation-based boundary node method for two-dimensional potential problems

In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This study combines the DIE method with the moving Kriging interpolation to present a boundary-type meshfree method, and the corresponding formulae of the MKIBNM are derived. In the present method, the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker＇s delta property, then the boundary conditions can be imposed directly and easily. To verify the accuracy and stability of the present formulation, three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.

A MESHLESS METHOD OF A THIN PLATE

A meshless approach to analysis of arbitrary Kirechhoff plates bythe local boundary integral equation (LBIF) method is presented. Themethod combines the advantageous features of all the three meth- ods:the Galerkin finite element method (GFEM), the boundary elementmethod (BEM) and the element- free Galerkin method (EFGM). It is atruly meshless method, which means that the discretization is inde-pendent of geometric subdivision into elements or cells, but is onlyboundary integration, however over a local boundary cen- tered) overa domain in question.

Application of A Fast Multipole BIEM for Flow Diffraction from A 3D Body

A Fast Multipole Method (FMM) is developed as a numerical approach to the reduction of the computational cost and requirement memory capacity for a large in solving large-scale problems. In this paper it is applied to the boundary integral equation method (BIEM) for current diffraction from arbitrary 3D bodies. The boundary integral equation is discretized by higher order elements, the FMM is applied to avoid the matrix/vector product, and the resulting algebraic equation is solved by the Generalized Conjugate Residual method (GCR). Numerical examination shows that the FMM is more efficient than the direct evaluation method in computational cost and storage of computers.

Dynamic rupture process of the 1999 Chi-Chi,Taiwan,earthquake

In this study, we preliminarily investigated the dynamic rupture process of the 1999 Chi-Chi, Taiwan, earthquake by using an extended boundary integral equation method, in which the effect of ground surface can be exactly included. Parameters for numerical modeling were carefully assigned based on previous studies. Numerical results indicated that, although many simplifications are assumed, such as the fault plane is planar and all heterogeneities are neglected, distribution of slip is still consistent roughly with the results of kinematic inversion, implying that for earthquakes in which ruptures run up directly to the ground surface, the dynamic processes are controlled by geometry of the fault to a great extent. By taking the common feature inferred by various kinematic inversion studies as a restriction, we found that the critical slip-weakening distance Dc should locate in a narrow region [60 cm, 70 cm], and supershear rupture might occur during this earthquake, if the initial shear stress before the mainshock is close to the local shear strength.

Nonlinear Capillary-Gravity Waves Produced by a Vertically Oscillating Plate

The nonlinear capillary-gravity wave produced by a vertically oscillating plate, in which the contact-angle model is considered, is studied by use of the Boundary Integral Equation Method (BIEM). The present numerical experiment shows that the code is robust and efficient for modeling the generation and propagation of capillary-gravity waves. It is found that the wave heights of stationary periodic nonlinear waves radiated away from the plate are dependent on the parameters involved in the contact-angle model. The effect of the contact-angle hysteresis and the nonlinearity of capillary-gravity waves on the wave profile is discussed in the paper.

Focusing phenomenon and near-trapped modes of SH waves

In this study, the null-field boundary integral equation method （BIEM） and the image method are used to solve the SH wave scattering problem containing semi-circular canyons and circular tunnels. To fully utilize the analytical property of Circular geometry, the polar coordinates are used to expand the closed-form fundamental solution to the degenerate kernel, and the Fourier series is also introduced to represent the boundary density. By collocating boundary points to match boundary condition on the boundary, a linear algebraic system is constructed. The unknown coefficients in the algebraic system can be easily determined. In this way, a semi-analytical approach is developed. Following the experience of near-trapped modes in water wave problems of the full plane, the focusing phenomenon and near-trapped modes for the SH wave problem of the half-plane are solved, since the two problems obey the same mathematical model. In this study, it is found that the SH wave problem containing two semi-circular canyons and a circular tunnel has the near-trapped mode and the focusing phenomenon for a special incident angle and wavenumber. In this situation, the amplification factor for the amplitude of displacement is over 300.

Analysis of the role of branching angle in the dynamic rupture process on a 3-D branching fault system

The fault branching phenomenon,which may heavily influence the patterns of rupture propagation in fault systems,is one of the geometric complexities of fault systems that is widely observed in nature.In this study,we investigate the effect of the branching angle on the rupture inclination and the interaction between branch planes in two-fork branching fault systems by numerical simulation and theoretical analysis based on Mohr’s circle.A friction law dependent on normal stress is used,and special attention is paid to studying how ruptures on the upper and lower branch planes affect the stress and rupture on each other separately.The results show that the two branch planes affect each other in different patterns and that the intensity of the effect changes with the branching angle.The rupture of the lower branch plane has a negative effect on the rupture of the upper branch plane in the case of a small branching angle but has almost no negative effect in the case of a large branching angle.The rupture of the upper branch plane,however,suppresses the rupture of the lower branch plane regardless of whether the branching angle is large or small.

PERTURBATION SOLUTION TO 3-D NONLINEAR SUPERCAVITATING FLOW PROBLEMS

A 3-D nonlinear problem of supercavitating flow past an axisymmetric body at a small angle of attack is investigated by means of the perturbation method and Fourier-cosine-expansion method. The first three order perturbation equations are derived in detail and solved numerically using the boundary integral equation method and iterative techniques. Computational results of the hydrodynamic characteristics and cavity shapes of each order are presented for nonaxisymmetric supercavitating flow past cones with various apex-angles at differ- ent cavitation numbers. The numerical results are found in good agreement with experimental data.

A High Frequency Boundary Element Method for Scattering by Convex Polygons with Impedance Boundary Conditions

We consider scattering of a time harmonic incident plane wave by a convex polygon with piecewise constant impedance boundary conditions.Standard finite or boundary element methods require the number of degrees of freedom to grow at least linearly with respect to the frequency of the incident wave in order to maintain accuracy.Extending earlier work by Chandler-Wilde and Langdon for the sound soft problem,we propose a novel Galerkin boundary element method,with the approximation space consisting of the products of plane waves with piecewise polynomials supported on a graded mesh with smaller elements closer to the corners of the polygon.Theoretical analysis and numerical results suggest that the number of degrees of freedom required to achieve a prescribed level of accuracy grows only logarithmically with respect to the frequency of the incident wave.

Hydroelastic Behaviour and Analysis of Marine Structures

Hydroelasticity of marine structures with and without forward speed is studied directly using time dependent Boundary Integral Equation Method with Neumann-Kelvin linearisation where the potential is considered as the impulsive velocity potential.The exciting and radiation hydrodynamic parameters are predicted in time with transient wave Green function whilst the structural analysis is solved with Euler-Bernoulli beam method at which modeshapes are defined analytically.The modal analysis is used to approximate the hydroelastic behaviour of the floating systems through fully coupling of the structural and hydrodynamic analyses.As it is expected,it is found with numerical experience that the effects of the rigid body modes are greater than elastic modes in the case of stiff structures.The predicted numerical results of the present in-house computational tool ITU-WAVE are compared with experimental results for validation purposes and show the acceptable agreements.

INTERACTION OF ELASTIC WAVES WITH A PERIODIC ARRAY OF COLLINEAR INPLANE CRACKS

A boundary integral equation method is applied to the study of the interaction of plane elastic waves with a periodic array of collinear inplane cracks.Numerical results are presented for the dynamic stress in- tensity factors.The effects of the wave type,wave frequency,wave incidence angle,and crack spacing on the dynamic stress intensity factors are analyzed in detail.

Dynamic interaction of twin vertically overlapping lined tunnels in an elastic half space subjected to incident plane waves

The scattering of plane harmonic P and SV waves by a pair of vertically overlapping lined tunnels buried in an elastic half space is solved using a semi-analytic indirect boundary integration equation method. Then the effect of the distance between the two tunnels, the stiffness and density of the lining material, and the incident frequency on the seismic response of the tunnels is investigated. Numerical results demonstrate that the dynamic interaction between the twin tunnels cannot be ignored and the lower tunnel has a significant shielding effect on the upper tunnel for high-frequency incident waves, resulting in great decrease of the dynamic hoop stress in the upper tunnel; for the low-frequency incident waves, in contrast, the lower tunnel can lead to amplification effect on the upper tunnel. It also reveals that the frequency-spectrum characteristics of dynamic stress of the lower tunnel are significantly different from those of the upper tunnel. In addition, for incident P waves in low-frequency region, the soft lining tunnels have significant amplification effect on the surface displacement amplitude, which is slightly larger than that of the corresponding single tunnel.

Diffraction of plane P waves around an alluvial valley in poroelastic half-space

Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green＇s fimctions of line source in poroelastic half-space, the scattered waves are constructed using the fictitious wave sources close to the interface of the valley and the density of ficti- tious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, and the comparison between the degenerated solutions and available results in single-phase case. Finally, the nature of diffraction of plane P waves around an alluvial valley in poroelastic half-space is investigated in detail through nu- merical examples.

Nearly Singular Integrals in 3D Stokes Flow

A straightforward method is presented for computing three-dimensional Stokes flow,due to forces on a surface,with high accuracy at points near the surface.The flowquantities arewritten as boundary integrals using the free-spaceGreen’s function.To evaluate the integrals near the boundary,the singular kernels are regularized and a simple quadrature is applied in coordinate charts.High order accuracy is obtained by adding special corrections for the regularization and discretization errors,derived here using local asymptotic analysis.Numerical tests demonstrate the uniform convergence rates of the method.