G-S变换的快速算法

在电磁场瞬变响应的数值计算中 ,常采用G S变换法作逆拉氏变换 .它是纯实数运算 ,而且只需对较少的拉氏变换变量s值作计算 （通常对每一采样时间选用 1 2个s值 ） ,因而是一种计算速度较快的算法 .但是 ,要对大量采样时间作计算 ,其计算量仍太大 .本文基于拉氏变换的延迟定理 ,建立了一种新的G S变换算法 .数值检验结果表明 ,新算法可成级次地减少对大量采样时间作G S变换的计算量 ,显著提高电磁场瞬变响应的计算速度 .     

双重介质地层变井筒储集效应典型曲线特征分析

本文在文献[1]对均质地层分析结果的基础上,用拉氏变换及stehfest数值反演方法研究了针对双重介质地层井筒储集系数呈线形变化,指数变化等情形下有效井径数学模型的近似解,绘制相应的典型曲线,并对影响曲线形态的参数进行探讨,这对于用实际资料进行试井解释时合理使用其早期数据有很好的参考作用。     

多边形有限单元形函数的几何构造方法

根据非均质材料的细观结构,划分多边形的计算网格。采用多边形单元进行有限元分析,实现了基于材料真实结构的数值模拟。给出了多边形有限单元形函数的几何构造方法,构造了一个辅助多边形,采用散度定理推导出多边形单元形函数的表达式。分析总结了多边形Wachspress插值、Laplace插值和平均值插值的构造方法和性质。     

渗透破坏若干问题的解析解及机理分析

讨论了柱坐标下不同介质有限区域内用分离变量法求解Laplace 方程的问题。通过一系列计算试验,给出了柱坐标下Laplace 方程在不同介质有限区域内用分离变量方法求的解,并用该解对渗透破坏现象进行了讨论和实例验证。认为渗透系数分布不均匀是发生渗透破坏的根本原因,渗透破坏朝着渗透系数差异较大的方向发展,并对管涌的复杂机理进行了引申阐述。该工作对于流土、管涌机理的更深入研究具有重要的意义。     

循环荷载下黏弹性饱和土层的一维固结

针对单层黏弹性地基Merchant模型,运用Laplace变换,求得频域内单层黏弹性地基的一维固结解。通过Laplace逆变换,计算了单层黏弹性地基在任意循环荷载下的有效应力及平均固结度。此外,结合工程实例,研究了Merchant模型各参数对循环荷载下黏弹性地基固结的影响。结果表明,在黏弹性地基的固结过程中,有效应力和沉降的发展速率是不一致的,黏壶的存在使地基固结初期的有效应力增长加快,而使固结后期的有效应力增长减慢,同时使变形的发展滞后于有效应力的发展。研究结果亦表明,循环荷载下土体的固结对独立弹簧模量的变化要比Kelvin体中弹簧模量的变化敏感。     

A new analytical model to study the influence of weld on the vertical dynamic response of prestressed pipe pile

This investigation is concerned with the mathematical analysis of a viscoelastic prestressed pipe pile embedded in multilayered soil under vertical dynamic excitation. The pile surrounding soil is governed by the plane strain model, and the soil plug is assumed to be an additional mass connected to the pipe pile shaft by applying the distributed Voigt model. Meanwhile, the prestressed pipe pile is assumed to be a vertical, viscoelastic, and hollow cylinder governed by the one‐dimensional wave equation. Then, analytical solutions of the dynamic response of the pipe pile in the frequency domain are derived by means of the Laplace transform and impedance function transfer method. Subsequently, the corresponding quasi‐analytical solution in the time domain for the case of the prestressed pipe pile undergoing a vertical semi‐sinusoidal exciting force applied at the pile top is obtained by employing the inverse Fourier transform. Utilizing these solutions, selected results for the velocity admittance curve and the reflected wave curve are presented for different heights of the soil plug to examine the influence of weld properties on the vertical dynamic response of prestressed pipe pile. The reasonableness of the theoretical model is verified by comparing the calculated results based on the presented solutions with measured results. Copyright © 2017 John Wiley & Sons, Ltd.     

基于改进泰勒级数法的总强度磁异常稳健向下延拓

位场向下延拓隶属于经典不适定问题,观测数据中细微的误差在向下延拓过程中都被严重放大,甚至会掩盖真实信息。如何精确地求解总强度磁异常（Bm）在垂直方向的各阶导数,是利用泰勒级数实现稳健向下延拓的关键。为此,本文首先分析了调和函数的相关性质,从理论上证明了Bm为准调和函数的结论,在精确计算各阶垂向导数基础上,提出利用改进泰勒级数实现磁场稳健向下延拓。为降低边界效应对向下延拓计算结果的影响,提出采用半余弦函数对磁场在4个方向上进行平滑扩边处理。通过球体与长方体仿真试验以及航空、船载实测磁场数据对提出方法进行了验证。结论表明,提出的技术方法可实现磁场稳健向下延拓,当观测数据无噪声时,计算结果精度要明显优于现行的FFT法、常规泰勒级数法以及积分-迭代法;当观测数据含有噪声时,本文方法和积分-迭代法计算结果精度相当。     

Computational issues and applications of line-elements to model subsurface flow governed by the modified Helmholtz equation

Two new approaches are presented for the accurate computation of the potential due to line elements that satisfy the modified Helmholtz equation with complex parameters. The first approach is based on fundamental solutions in elliptical coordinates and results in products of Mathieu functions. The second approach is based on the integration of modified Bessel functions. Both approaches allow evaluation of the potential at any distance from the element. The computational approaches are applied to model transient flow with the Laplace transform analytic element method. The Laplace domain solution is computed using a combination of point elements and the presented line elements. The time domain solution is obtained through a numerical inversion. Two applications are presented to transient flow fields, which could not be modeled with the Laplace transform analytic element method prior to this work. The first application concerns transient single-aquifer flow to wells near impermeable walls modeled with line-doublets. The second application concerns transient two-aquifer flow to a well near a stream modeled with line-sinks.     

An improved numerical simulation mode of nonlinear wave with consideration of high order terms

In this paper the 0-1 combined BEM is adopted to subdivide the computational domain boundary,and to discretize the Green’s integral expression based on Laplace equation.The FEM is used to subdivide the wave surface and deduce the surface equation which satisfies the nonlinear boundary conditions on the surface.The equations with potential function and wave surface height as an unknown quantity by application of Taylor expansion approach can be solved by iteration within the time step.In m-time iteration within the computational process of time step(n-1)Δt to nΔt,the results of the previous iteration are taken as the initial value of the two-order unknown terms in the present iteration.Thus,an improved tracking mode of nonlinear wave surface is established,and numerical results of wave tank test indicate that this mode is improved obviously and is more precise than the previous numerical model which ignored the two-order unknown terms of wave surface location and velocity potential function in comparison with the theoretical values.     

Fractional Reaction-Diffusion Equations

In a series of papers, Saxena et al. (2002, 2004a, 2004b) derived solutions of a number of fractional kinetic equations in terms of generalized Mittag-Leffler functions which provide the extension of the work of Haubold and Mathai (1995, 2000). The subject of the present paper is to investigate the solution of a fractional reaction-diffusion equation. The results derived are of general nature and include the results reported earlier by many authors, notably by Jespersen et al. (1999) for anomalous diffusion and del-Castillo-Negrete et al. (2003) for reaction-diffusion systems with Lévy flights. The solution has been developed in terms of the H-function in a compact form with the help of Laplace and Fourier transforms. Most of the results obtained are in a form suitable for numerical computation.