Study on Characteristics of 3-D Translating-Pulsating Source Green Function of Deep-Water Havelock Form and Its Fast Integration Method

The singularities, oscillatory performances and the contributing factors to the 3-＇D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.

Continuum Skyrme Hartree-Fock-Bogoliubov theory with Green’s function method for neutron-rich Ca,Ni,Zr,and Sn isotopes

The possible exotic nuclear properties in the neutron-rich Ca,Ni,Zr,and Sn isotopes are examined with the continuum Skyrme Hartree-Fock-Bogoliubov theory in the framework of the Green’s function method.The pairing correlation,the couplings with the continuum,and the blocking effects for the unpaired nucleon in odd-A nuclei are properly treated.The Skyrme interaction SLy4 is adopted for the ph channel and the density-dependentinteraction is adopted for the pp chan-nel,which well reproduce the experimental two-neutron separation energies S_(2n)and one-neutron separation energies Sn.It is found that the criterion S_(n)>0 predicts a neutron drip line with neutron numbers much smaller than those for S_(2n)>0.Owing to the unpaired odd neutron,the neutron pairing energies−E_(pair)in odd-A nuclei are much lower than those in the neighbor-ing even-even nuclei.By investigating the single-particle structures,the possible halo structures in the neutron-rich Ca,Ni,and Sn isotopes are predicted,where sharp increases in the root-mean-square(rms)radii with significant deviations from the traditional rA^(1∕3)rule and diffuse spatial density distributions are observed.Analyzing the contributions of various partial waves to the total neutron densityρlj(r)∕ρ(r)reveals that the orbitals located around the Fermi surface-particularly those with small angular momenta-significantly affect the extended nuclear density and large rms radii.The number of neutrons Nλ(N_(0))occupying above the Fermi surfacen(continuum threshold)is discussed,whose evolution as a function of the mass number A in each isotope is consistent with that of the pairing energy,supporting the key role of the pairing correlation in halo phenomena.

Adomian Decomposition Method with Green’s Function for Solving Tenth-Order Boundary Value Problems

In this paper, the Adomian decomposition method with Green’s function (Standard Adomian and Modified Technique) is applied to solve linear and nonlinear tenth-order boundary value problems with boundary conditions defined at any order derivatives. The numerical results obtained with a small amount of computation are compared with the exact solutions to show the efficiency of the method. The results show that the decomposition method is of high accuracy, more convenient and efficient for solving high-order boundary value problems.

In-plane dynamic Green's functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic half-space

The dynamic stiffness method combined with the Fourier transform is utilized to derive the in-plane Green’s functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic（TI）half-space.The loaded layer is fixed to obtain solutions restricted in it and the corresponding reactions forces,which are then applied to the total system with the opposite sign.By adding solutions restricted in the loaded layer to solutions from the reaction forces,the global solutions in the wavenumber domain are obtained,and the dynamic Green’s functions in the space domain are recovered by the inverse Fourier transform.The presented formulations can be reduced to the isotropic case developed by Wolf（1985）,and are further verified by comparisons with existing solutions in a uniform isotropic as well as a layered TI halfspace subjected to horizontally distributed loads which are special cases of the more general problem addressed.The deduced Green’s functions,in conjunction with boundary element methods,will lead to significant advances in the investigation of a variety of wave scattering,wave radiation and soil-structure interaction problems in a layered TI site.Selected numerical results are given to investigate the influence of material anisotropy,frequency of excitation,inclination angle and layered on the responses of displacement and stress,and some conclusions are drawn.

Green Functions with Pulsating Sources in A Two-Layer Fluid of Finite Depth

The derivation of Green function in a two-layer fluid model has been treated in different ways. In a two-layer fluid with the upper layer having a free surface, there exist two modes of waves propagating due to the free surface and the interface. This paper is concerned with the derivation of Green functions in the three dimensional case of a stationary source oscillating. The source point is located either in the upper or lower part of a two-layer fluid of finite depth. The derivation is carried out by the method of singularities. This method has an advantage in that it involves representing the potential as a sum of singularities or multipoles placed within any structures being present. Furthermore, experience shows that the systems of equations resulted from using a singularity method possess excellent convergence characteristics and only a few equations are needed to obtain accurate numerical results. Validation is done by showing that the derived two-layer Green function can be reduced to that of a single layer of finite depth or that the upper Green function coincides with that of the lower, for each case. The effect of the density on the internal waves is demonstrated. Also, it is shown how the surface and internal wave amplitudes are compared for both the wave modes. The fluid in this case is considered to be inviscid and incompressible and the flow is irrotational.

Side Wall Effects on the Hydrodynamics of a Floating Body by ImageGreen Function Based on TEBEM

A novel numerical model based on the image Green function and first-order Taylor expansion boundary element method(TEBEM), which can improve the accuracy of the hydrodynamic simulation for the non-smooth body, was developed to calculate the side wall effects on first-order motion responses and second-order drift loads upon offshore structures in the wave tank. This model was confirmed by comparing it to the results from experiments on hydrodynamic coefficients, namely the first-order motion response and second-order drift load upon a hemisphere, prolate spheroid, and box-shaped barge in the wave tank. Then,the hydrodynamics of the KVLCC2 model were also calculated in two wave tanks with different widths. It was concluded that this model can predict the hydrodynamics for offshore structures effectively, and the side wall has a significant impact on the firstorder quantities and second-order drift loads, which satisfied the resonant rule.

Time Domain Rankine-Green Panel Method for Offshore Structures

To solve the numerical divergence problem of the direct time domain Green function method for the motion simulation of floating bodies with large flare, a time domain hybrid Rankine-Green boundary element method is proposed. In this numerical method, the fluid domain is decomposed by an imaginary control surface, at which the continuous condition should be satisfied. Then the Rankine Green function is adopted in the inner domain. The transient free surface Green function is applied in the outer domain, which is used to find the relationship between the velocity potential and its normal derivative for the inner domain. Besides, the velocity potential at the mean free surface between body surface and control surface is directly solved by the integration scheme. The wave exciting force is computed through the convolution integration with wave elevation, by introducing the impulse response function. Additionally, the nonlinear Froude-Krylov force and hydrostatic force, which is computed under the instantaneous incident wave free surface, are taken into account by the direct pressure integration scheme. The corresponding numerical computer code is developed and first used to compute the hydrodynamic coefficients of the hemisphere, as well as the time history of a ship with large flare; good agreement is obtained with the analytical solutions as well as the available numerical results. Then the hydrodynamic properties of a FPSO are studied. The hydrodynamic coefficients agree well with the results computed by the frequency method; the influence of the time interval and the truncated time is investigated in detail.

NEW TYPE OF DUAL BEM AND GREEN'S-FUNCTION-LIBRARY STRATEGY FOR FRACTURE ANALYSIS IN COMPLEX STRUCTURES

A new type of dual boundary integral equations(DBIE)is presented first,through which,a smaller system of equations needs to be solved in fracture analysis.Then a non-conforming crack tip element in two-dimensional problems is proposed.The exact formula for the hypersingular integral over the non-con- forming crack tip element is given next.By virtue of Green's-function-library strategy,a series of stress in- tensity factors(SIF)of different crack orientations,locations and/or sizes in a complicated structure can be obtained easily and efficiently.Finally,several examples of fracture analysis in two dimensions are given to demonstrate the accuracy and efficiency of the method proposed.

用磁流源并矢Green函数的递推矩阵方法计算层状各向异性地层中多分量感应测井响应

采用递推矩阵方法计算各向异性介质的磁流源并矢Green函数,并利用上述并矢Green函数对层状各向异性倾斜地层中多分量感应测井的响应进行数值模拟,分析了线圈距、层厚、倾角和围岩对多分量感应测井响应的影响.计算发现,共面视电导率比共轴视电导率变化规律复杂,不能反映地层电导率的真实情况.各向异性地层共轴视电导率随井眼相对倾角的增加而减小,而共面视电导率随井眼相对倾角的增加而增加.仪器的垂向分辨率、围岩各向异性对目的层响应的影响程度取决于线圈距,围岩各向异性对高电导率目的层中测井响应的影响大于对高电阻率目的层中测井响应的影响.

无航速具有自由液面的三维脉动源Green函数的数值积分方法

提出了一种计算无航速具有自由液面的三维脉动源Green函数的方法.采用零阶Bessel函数的Laplace变换,将无航速具有自由液面的三维脉动源Green函数变形为易于截断处理的形式,实现无穷积分向有限积分的转变.采用样条插值技术和自适应方法解决了高振荡函数的数值积分的速度和精度问题.对半球浮体进行了数值计算,结果令人满意,表明了该方法的有效性和实用性.

三维层状饱和地基动力Green函数的薄层法求解

建立了三维层状饱和地基模型并在模型底部施加无反射边界条件,引入Fourier变换和正交变换将饱和地基Biot动力方程化为Love和Rayleigh模态方程,将地基沿深度方向划分为数个薄层单元并对单元位移场引入线性插值函数,从而将模态方程离散化,通过求解本征值问题最终给出了三维层状饱和地基的频域动力Green函数。分析了均质饱和地基的Love和Rayleigh模态的弥散特性,求解了均质和成层饱和半空间的Lamb问题,并与Philippacopoulos的数值积分解对比,说明了所求Green函数的正确性。所得Green函数无需关于Hankel函数进行数值逆变换且能直接考虑地基的成层性,提高了计算效率。

直角平面区域内固定圆形刚性夹杂问题的Green函数解

利用复变函数法、多极坐标移动技术研究了直角平面区域内含有固定圆形夹杂时的反平面问题Green函数解.首先构造出不含夹杂的完整直角平面区域内满足边界应力条件的入射位移场;其次,建立直角平面区域内固定圆形夹杂对该入射场产生的满足直角边界应力自由条件的散射波解,并由叠加原理得到介质内的总波场.最后利用夹杂边界处的位移条件确定出散射波解中的未知系数,最终得到问题的Green函数解,还通过算例讨论了夹杂边界处的径向应力和环向应力随不同波数、角度和不同夹杂位置及不同点源位置的变化情况.算例结果表明了该文方法的有效实用性.

圆形区域多圆孔对SH波散射Green函数解

利用多极坐标变换、Graf加法公式及波函数展开法研究圆形区域内多圆孔缺陷存在时对圆形域边界作用反平面稳态脉冲SH波散射Green函数解。给出各圆孔对SH波散射产生含未知系数的散射波函数,用叠加原理写出圆形域介质内总波场,据边界条件、Graf加法公式及Fourier积分变换获得确定未知系数的无穷线代数方程组;由计算精度进行有限项截断求解,确定问题的Green函数解。用算例结果验证该方法的可行性。

变刚度梁挠度曲线的Green函数法求解

在变刚度梁的线弹性问题中,求解梁受静力荷载的挠度曲线常用解法有积分法与单位荷载法.本文从变刚度梁挠度曲线的微分方程出发,给出了变刚度梁挠度曲线的Green函数法解答,并分析了该解法的优点.从推导结果可以看到,本文提出的公式具有统一、精确、简洁、适合电算的特点,在编制杆系结构计算软件中将具有重要应用价值.

平面分层媒质空域Green函数快速计算(Ⅰ)

格林函数法是电磁场数值分析中的一种基本方法,平面分层媒质的空域格林函数通常表达成Sommer-feld积分形式。文中针对该问题快速求解的几种基本方法进行了系统分析,对各种分析方法进行了分类,根据研究进展情况,概述了传统离散复镜像法(DCIM)、二阶DCIM、增强型DCIM(EDCIM)、柱面波逼近和支割线积分以及其他方法,阐明了各种方法的优缺点;对计算得到的闭式Green函数在具体电磁散射与辐射计算中产生的误差控制、定阶、维数、媒质有耗问题进行了分析,给出了一种简便的有耗媒质DICM方法,数值计算结果显示文中方法的正确性。

高频电磁场问题中的并矢Green函数方法

高频电磁辐射、散射问题中激励与响应之间的线性关系成立。这一关系可用一并矢变换函数给出。因此高频电磁场问题关键在于寻求并矢变换函数。并矢格林函数理论能够给出典型电磁场问题激励与响应之间的线性关系，借助高频电磁场的局域性原理、传播与绕射遵循的费马原理，极易求得高频电磁场问题的并矢变换函数。本文讨论了高频电磁场问题中的并矢格林函数方法，研究了并矢格林函数方法在寻求并矢变换函数中的应用，得到了部分工程问题的理论结果。

Poisson方程边值问题的Green函数法

归纳出在求解Poisson方程边值问题时引入Green函数的物理意义,讨论了电像法、级数法和保角变换法等三种方法求解Green函数的适用范围.

任意边界条件的连续板的Green函数法计算

本方法将板在特定域中的 Green 函数作为影响函数,首先根据连续板的外边界条件以及内支承条件建立方程,求出虚拟域中的 Green 函数“源”以及连续板的内支承反力,继而由求得的“源值”、内支承反力和板上的已知荷载确定板内任意点的挠度和内力.方法简单,易于编程序,且未知量比边界元分域法进一步减少,适应性强,不受板的形状以及边界条件的限制。文中附若干算例,并将其数值结果与有限元法比较,表明本方法精度良好.

准Green函数Galerkin逼近的权范数估计的探讨

用权范数法探讨并得到了准Green函数的一些性质及其Galerkin逼近的新估计.

饱和多孔介质的二维Green函数及地震反应的BEM计算

按照De Hoop的方法,通过三维解答的积分得到集中力作用下饱和多孔介质U-P形式的二维Green函数,利用耦合纵波的解耦及分波类型的奇异点处理,完成了集中力作用下半空间多孔介质动力响应的BEM计算,通过Duhamel公式积分分别给出了在输入EI Centro和2008年汶川地震加速度记录后的饱和土隧道的地震响应计算。研究结果表明:隧道的地震放大作用明显;EI Centro或汶川2008地震记录作用下,隧道底部的响应均要大于中间和顶部的响应。