On using mesh-based and mesh-free methods in problems defined by Eringen’s non-local integral model: issues and remedies

With the recent success of nonlocal theories in modeling of engineering problems involving small intrinsic length scales, such as modeling of crack propagation, this paper addresses issues pertaining to cost-ineffectiveness of Eringen’s integral model. The cost effectiveness of the computation may be considered as a twofold issue; one pertaining to the non-local model and another pertaining to the numerical tool. First of all, we shall show that during the solution of problems with Eringen’s non-local integral model, there is no need to consider the integral model for the whole computational domain. In fact, the problems may be solved by just using the integral model close to the boundaries, i.e. a boundary layer effect, or around the points with singularities. In this paper we propose a partitioning strategy to remarkably reduce the computational cost. This may be considered as a gateway for solving some types of two-scale problems, e.g. those with macro/micro and nano scales, in which the small scale effects are localized just at parts of the domain. To demonstrate the efficiency of the numerical tools, we examine the performance of the finite element method (FEM), the element free Galerkin method (EFG) and the finite point method (FPM). This paves the way for using mesh-free methods in the solution of problems with non-local integral models. Examples with smooth and non-smooth solutions are considered for examining the efficiency of the methods. It will be shown that, by considering the boundary layer effect, the FEM and FPM will be efficient enough for being used in problems defined by Eringen’s non-local integral model.

     

Wave kinematic fields from the boundary integral method

The predictive potential of interior domain solutions from the boundary integral method for 2D extreme wave kinematics is explored. Comparisons with analytical solutions for near‐limit waves confirms the susceptibility of the boundary integral method to poor precision at near‐boundary locations. Additionally, these comparisons identify a domain‐wide precision challenge that is associated with the relatively rapid changes in water surface geometry and kinematics that are typical of extreme waves. A numerical evaluation of Green's integral around the boundary addresses these precision issues through formulation of the integration as a simultaneous system of ordinary differential equations at a cubic level of approximation. Careful attention is given to consistent interpolation of all contributions to the Green's integral. Copyright © 2005 John Wiley & Sons, Ltd.     

基于概率统计方法的微弱GPS信号捕获算法

微弱GPS信号的信噪比较低,需采用更长的相干和非相干积分时间增强信噪比,提出一种从概率统计角度分析相干和非相干积分运算过程的方法。首先给出了复数下变频信号输入的相关运算表达式;然后推导了相干积分的相关输出数学期望和方差,计算出单元捕获失败的概率密度和概率表达式;最后,用Matlab对GPS信号的相干和非相干积分进行仿真。仿真结果验证了捕获失败概率表达式的结论,即由于高斯近似和非相干平均产生的性能损失,非相干积分的捕获性能比相干积分弱3.5 dB。     

New Green’s function for stress field and a note of its application in quantum-wire structures

The elastic field caused by the lattice mismatch between the quantum wires and the host matrix can be modeled by a corresponding two-dimensional hydrostatic inclusion subjected to plane strain conditions. The stresses in such a hydrostatic inclusion can be effectively calculated by employing the Green’s functions developed by Downes and Faux, which tend to be more efficient than the conventional method based on the Green’s function for the displacement field. In this study, Downes and Faux’s paper is extended to plane inclusions subjected to arbitrarily distributed eigenstrains: an explicit Green’s function solution, which evaluates the stress field due to the excitation of a point eigenstrain source in an infinite plane directly, is obtained in a closed-form. Here it is demonstrated that both the interior and exterior stress fields to an inclusion of any shape and with arbitrarily distributed eigenstrains are represented in a unified area integral form by employing the derived Green’s functions. In the case of uniform eigenstrain, the formulae may be simplified to contour integrals by Green’s theorem. However, special care is required when Green’s theorem is applied for the interior field. The proposed Green’s function is particularly advantageous in dealing numerically or analytically with the exterior stress field and the non-uniform eigenstrain. Two examples concerning circular inclusions are investigated. A linearly distributed eigenstrain is attempted in the first example, resulting in a linear interior stress field. The second example solves a circular thermal inclusion, where both the interior and exterior stress fields are obtained simultaneously.     

Stroh formalism in evaluation of 3D Green’s function in thermomagnetoelectroelastic anisotropic medium

The paper presents studies on the Green’s function for thermomagnetoelectroelastic medium and its reduction to the contour integral. Based on the previous studies the thermomagnetoelectroelastic Green’s function is presented as a surface integral over a half-sphere. The latter is then reduced to the double integral, which inner integral is evaluated explicitly using the complex variable calculus and the Stroh formalism. Thus, the Green’s function is reduced to the contour integral. Since the latter is evaluated over the period of the integrand, the paper proposes to use trapezoid rule for its numerical evaluation with exponential convergence. Several numerical examples are presented, which shows efficiency of the proposed approach for evaluation of Green’s function in thermomagnetoelectroelastic anisotropic solids.     

SCATTERING OF ANTI-PLANE SHEAR WAVES IN A FUNCTIONALLY GRADED MATERIAL STRIP WITH AN OFF-CENTER VERTICAL CRACK

The scattering problem of anti-plane shear waves in a functionally graded material strip with an off-center crack is investigated by use of Schmidt method. The crack is vertically to the edge of the strip. By using the Fourier transform, the problem can be solved with the help of a pair of dual integral equations that the unknown variable is the jump of the displacement across the crack surfaces. To solve the dual integral equations, the jump of the displacement across the crack surfaces was expanded in a series of Jacobi polynomials. Numerical examples were provided to show the effects of the parameter describing the functionally graded materials, the position of the crack and the frequency of the incident waves upon the stress intensity factors of the crack.     

An efficient approximate numerical method for modeling contact of materials with distributed inhomogeneities

This research explores the influence of distributed non-interpenetrating inhomogeneities on the contact of inhomogeneous materials via a new efficient numerical model based on Eshelby’s Equivalent Inclusion Method. The half-space contact of a sphere with an inhomogeneous material is considered, and the solutions take into account interactions between all inhomogeneities. The efficiency and solution accuracy of the proposed method are demonstrated through comparative studies with those of an existing numerical method and the finite element method. The influence of spatial inhomogeneity orientations on the contact elastic field is investigated and parametric studies are conducted for the effect of arbitrarily distributed inhomogeneities on the stress field of the materials. The significance of the influences of inhomogeneity distribution parameters on the inverse volumetric stress integral is quantified and the corresponding data are fitted into selected several formulas as a step towards understanding the rolling contact fatigue life of the materials.     

Regular wave integral approach to the prediction of hydrodynamic performance of submerged spheroid

A free surface Green function method is employed in numerical simulations of hydrodynamic performance of a submerged spheroid in a fluid of infinite depth. The free surface Green function consists of the Rankine source potential and a singular wave integral. The singularity of the wave integral is removed with the use of the Havelock regular wave integral. The finite boundary element method is applied in the discretisation of the fluid motion problem so that the panel integral of the Rankine source potential is evaluated by the Hess–Smith formula and the panel integral of the regular wave integral is evaluated in a straightforward way due to the regularity nature. Present method’s results are in good agreement with earlier numerical results.     

Capillary finite-amplitude waves with Tolman’s nonlinearity

Capillary finite-amplitude waves on a fluid surface are analyzed numerically by means of the conformal mapping method with allowance for Tolman’s nonlinearity which reflects the dependence of the surface tension coefficient on the surface curvature. The dependence of the wave profile and the wave velocity on the Tolman length is analyzed. Keywords: capillary waves, Tolman’s nonlinearity.     

A Comparison of Nonlinear Water Wave Models

We compare the numerical evolution of one-dimensional gravity waves in response to a traveling surface pressure pulse using a highly accurate boundary integral method and two relatively efficient approximate models (West et al. and Benney–Luke). In both water of finite-depth and in the deep-water limit, the steady state effect of the decaying pressure ramp is to create a profile which approximates a Stokes wave. Moreover, the transient surface profile appears to evolve through a series of Stokes waves of time-varying amplitude. Results show all three models yield similar predictions for lower amplitude waves, while the West et al. and boundary integral predictions differ from the Benney–Luke model at higher amplitudes.     

Analysis of damped guided waves using the method of multiple scales

We analytically investigate the influence of damping on Lamb waves, which are a specific type of guided wave in two-dimensional plates. Considering material attenuation, we suppose that Lamé constants are complex numbers. This leads to the associated wavenumbers being complex, with the imaginary part of the wavenumber being associated with effect of attenuation of the guided wave. In this paper, we show how dispersion curves and attenuation coefficients can be obtained using the self-adjointness and the method of multiple scales (MMS), which is a type of perturbation method. Using the self-adjointness and the MMS, we can calculate the frequency- and wavenumber-dependent attenuation coefficients from the integral values and boundary values of a corresponding eigenfunction with respect to each propagation mode. This analytical method can yield not only dispersion curves but also mode-by-mode attenuation coefficients regardless of the numerical initial values, unlike numerical approaches using the Newton method. Thus, the proposed method can more easily calculate the attenuation coefficients with respect to a particular mode than conventional methods. Furthermore, the results obtained by proposed method were in good agreement with those obtained by the semi-analytical finite element (SAFE) method, which validates the proposed method.     

Three-dimensional desingularized boundary integral methods for potential problems

The concept of desingularization in three-dimensional boundary integral computations is re-examined. The boundary integral equation is desingularized by moving the singular points away from the boundary and outside the problem domain. We show that the desingularization gives better solutions to several problems. As a result of desingularization, the surface integrals can be evaluated by simpler techniques, speeding up the computation. The effects of the desingularization distance on the solution and the condition of the resulting system of algebraic equations are studied for both direct and indirect versions of the boundary integral method. Computations show that a broad range of desingularization distances gives accurate solutions with significant savings in the computation time. The desingularization distance must be carefully linked to the mesh size to avoid problems with uniqueness and ill-conditioning. As an example, the desingularized indirect approach is tested on unsteady non-linear three-dimensional gravity waves generated by a moving submerged disturbance; minimal computational difficulties are encountered at the truncated boundary.     

Compressional waves in fluid-saturated porous solid containing a penny-shaped crack

Diffraction of normal compression waves by a penny-shaped crack in a fluid-saturated porous medium is investigated. Two wave types are considered, namely, compressional wave of the first kind, and the second kind. The former, also known as fast wave, propagates primarily through the solid, whereas the latter or slow wave, propagates mainly in the fluid. Each wave propagates in the medium along with induced wave of the same type in the companion constituent of the material. Application of Biot’s theory in conjunction with integral transform technique reduces the problem to a mixed boundary-value problem whose solution is in turn governed by a Fredholm integral equation of the second kind. Near-field and far-field solutions are obtained in terms of the dynamic stress-intensity factor and the scattering cross section, respectively. They are of particular importance to the linear elastic fracture mechanics (LEFM) and in the scattering theory of elastic waves. The mode I stress-intensity factors are computed numerically for a set of selected material property values, and shown graphically for various mass density and viscosity-to-permeability ratios. The obtained results reveal significant impact of the presence of pore fluid upon the stress-intensity factors, both magnitudes and frequencies at their peak values. The influence of the fluid is also observed from the calculated scattering cross sections of the scattered far-field. Accuracy of the present solution procedure is verified by comparing the numerical results with existing results in the limiting case of dry elastic materials.     

Continuation of travelling‐wave solutions of the Navier–Stokes equations

An efficient way of obtaining travelling waves in a periodic fluid system is described and tested. We search for steady states in a reference frame travelling at the wave phase velocity using a first‐order pseudospectral semi‐implicit time scheme adapted to carry out the Newton's iterations. The method is compared to a standard Newton–Raphson solver and is shown to be highly efficient in performing this task, even when high‐resolution grids are used. This method is well suited to three‐dimensional calculations in cylindrical or spherical geometries. Copyright © 2006 John Wiley & Sons, Ltd.     

Water wave radiation by a submerged rough disc

A thin circular body is submerged below the free surface of deep water. The problem is reduced to a hypersingular integral equation over the boundary of the body. Using a perturbation method, the problem is then reformulated by a sequence of simpler hypersingular equations over a flat disc making it well suited for an efficient previously used solution method. The first order approximation is computed and the hydrodynamic force due to heaving radiation motion are presented in terms of the added mass and damping coefficients for a polynomial cap and for a rough disc, modelled by a superposition of sinusoidal surfaces defined by randomly generated parameters. The solution exhibits larger maxima associated with smaller volume of submergence of the body. A slight shift of the damping coefficient maxima to lower frequencies is noticed for the caps. Rough discs with similar statistical properties exhibit different behaviours. Thus, it is the exact specific form of the rough disc that dictates the hydrodynamic force.     

Static equilibrium of a thermoelastic half-plane: Green’s functions and solutions in integrals

This article presents in a closed form new influence functions of a unit point heat source on the displacements for three boundary value problems of thermoelasticity for a half-plane. We also obtain the corresponding new integral formulas of Green’s and Poisson’s types that directly determine the thermoelastic displacements and stresses in the form of integrals of the products of specified internal heat sources or prescribed boundary temperature and constructed already thermoelastic influence functions (kernels). All these results are presented in terms of elementary functions in the form of three theorems. Based on these theorems and on derived early by author the general Green-type integral formula, we obtain in elementary functions new solutions to two particular boundary value problems of thermoelasticity for half-plane. The graphical presentation of the temperature and thermal stresses of one concrete boundary value problems of thermoelasticity for half-plane also is included. The proposed method of constructing thermoelastic Green’s functions and integral formulas is applicable not only for a half-plane, but also for many other two- and three-dimensional canonical domains of different orthogonal coordinate systems.     

A 2.5-D dynamic model for a saturated porous medium. Part II: Boundary element method

The 3-D boundary integral equation is derived in terms of the reciprocal work theorem and used along with the 2.5-D Green’s function developed in Part I [Lu, J.F., Jeng, D.S., Williams, S., submitted for publication. A 2.5-D dynamic model for a saturated porous medium: Part I. Green’s function. Int. J. Solids Struct.] to develop the 2.5-D boundary integral equation for a saturated porous medium. The 2.5-D boundary integral equations for the wave scattering problem and the moving load problem are established. The Cauchy type singularity of the 2.5-D boundary integral equation is eliminated through introduction of an auxiliary problem and the treatment of the weakly singular kernel is also addressed. Discretisation of the 2.5-D boundary integral equation is achieved using boundary iso-parametric elements. The discrete wavenumber domain solution is obtained via the 2.5-D boundary element method, and the space domain solution is recovered using the inverse Fourier transform. To validate the new methodology, numerical results of this paper are compared with those obtained using an analytical approach; also, some numerical results and corresponding analysis are presented.     

用有限积分计算曲率复杂分布下的结构挠度

有限积分法是Brown和Trahair在求解微分方程时采用的数值解法, 其核心环节是已知函数z= z(x)的导函数z′ =z′(x)的某些值的情况下数值分析z的方法. 由曲率φ计算挠度, 实质意义上是由z″计算z的数学问题. 基于有限积分法给出的z与z″之间及z′与z″之间的数值关系, 通过矩阵运算推导得到了挠曲矩阵,通过引入转换式φ= -z″得到了曲率挠度关系式, 讨论了几种常见边界条件下的曲率挠度关系, 提出了曲率复杂分布情况下结构挠度计算的有限积分方法.     

An optimal method for seismic drift design of concrete buildings using gradient and Hessian matrix calculations

This paper describes a novel seismic optimal design method for the reinforced concrete frame. First, an optimal mathematical model with time-dependent constraints, i.e., inter-story drift constraints, is established for achieving minimum weight design. Second, the inequality constraint problem with time-dependent constraints is converted into a sequence of appropriately formed unconstrained problems using the integral interior point penalty function method. Third, an efficient algorithm of the first and second derivatives of the inter-story drift with respect to design variables is formulated based on Newmark-β method. Gradient and Hessian matrix of the integral interior penalty function are also computed. Fourth, Marquardt’s method is employed to solve a sequence of unconstrained problems. Finally, the minimum weight design of a three-story, two-bay planar frame is demonstrated using the new optimization method and the augmented Lagrange multiplier method. The comparative results show the seismic optimal design method presented in this paper is more efficient than the augmented Lagrange multiplier method in terms of computational time. The proposed new method is an effective and efficient approach for minimum weight design of the reinforced concrete frames subjected to earthquake excitation.     

关于线性热释电弹性介质的互等功定理及应用

根据线弹性理论以及电介质理论推出了类似于Ｂｅｔｔｉ互等功定理的热电弹性耦合体─—热释电弹性体的互等功方程；作为一种应用推导了在边界积分方法中被广泛采纳的热电弹性Ｓｏｍｉｇｌｉａｎａ方程．