横观各向同性饱和土体三维粘弹性动力分析

采用针对横观各向同性饱和土体u-w形式三维粘弹性动力方程,考虑土骨架的粘弹性性质且基于粘弹性理论,通过运用Fourier 展开、Laplace 和Hankel 积分变换方法和引入中间变量,将含有粘弹性参数的六元二阶偏微分运动控制方程组,化为2组各含4个未知变量的常微分方程组,从而给出了柱坐标系下粘弹性横观各向同性饱和土体在非轴对称动力荷载作用下的瞬态反应的土骨架位移分量、孔隙流体相对于土骨架的位移分量瞬态反应一般解。在此基础上,引入初始条件和边界条件,对垂直向和水平向动力荷载作用下半空间边值问题进行了求解。根据动力时域解答的一般解,利用Laplace和Hankel 数值逆变换技术,编制了相应的数值计算程序。并进行了实例验证和弹性、粘弹性解的对比分析。结果表明,在进行横观各向同性饱和土体动力分析时,考虑土骨架的粘弹性是必要的。     

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.     

基于LGCP的城市管理事件空间点模式分析

数字化城管系统累积了大量城管事件历史数据,充分挖掘事件背后的空间分布模式和事件成因机制能够为城管部门的管控工作提供决策支持。论文利用Log Gaussian Cox Processes(LGCP)模型分析了西北某地H市P区的街面秩序类、市容环境类和宣传广告类城管事件之间的空间分布差异和事件成因影响差异。研究发现：① 3类城管事件都呈现出明显的空间聚集,其空间聚集尺度最远不超过924 m;② 各类事件聚集的特征各不相同,街面秩序类贴近城区主要干道,呈路网状。市容环境类表现出在城区中心块状聚集,周边地区零星分散的特征。宣传广告类靠近交通干线呈长条状,靠近商业中区域呈块状分布;③ 城区内不同类别的POI对城管事件的影响大小不同。购物服务类、医疗保健类和居民住宅类表现出最显著的影响,说明特定区域内人群的流量和密度是影响城管事件分布的重要因素,人群的流动和聚集会加剧城管事件数量的增加。研究结果能够满足城管部门的城管事件空间分布热点识别以及事件成因分析的需求。     

关于固体表面上液体球冠的平衡条件问题——兼评"冷凝器壁面滴状冷凝的热力学机理及最佳接触角"等文章

对于第一类液滴(尺度远大于界面层的厚度),无论是远离固体壁面的液体球或附着在壁面上的球冠,其内外压力差(简称"附加压力")均适用经典Laplace公式,并且特别对球冠情况给出了一种新的整体性证明.还澄清有关争论:指出[曹治党、郭愚1999物理学报481823]一文对附壁面第一类液体球冠所推导出的附加压力与接触角有关的公式是错误的,而[闵敬春2002物理学报512730]是正确的.     

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.     

Complex spatial feedbacks of tephra redistribution,ice melt and surface roughness modulate ablation on tephra covered glaciers

Tephra fallout from the 2011 Grímsvötn eruption onto Svínafellsjökull, Iceland, created an ice‐ash landscape of a type that is rarely studied but is nevertheless common in glacio‐volcanic regions. We used terrestrial laser scanning (TLS) to measure ice surface topography and absorption at high spatial resolution, confirming ablation rates either reduce or increase under thick (insulating) and thin (reduced albedo) ash deposits, respectively. Fourier transform analysis of the TLS data identified that a three‐fold increase in aerodynamic roughness was attributable to an increase in larger (> 0·2 m) surface features. Moreover, TLS measurements revealed the importance of ash redistribution by meltwater in generating differential melting which modifies roughness and ash patchiness, such that the net effect of these spatial ash–ice feedbacks was to reduce ablation rates by up to 59%. The modulating effects of these previously undocumented feedbacks on ablation rates are, therefore, significant and must be correctly parameterized if ash‐covered glacier mass balances are to be predicted correctly. Copyright © 2012 John Wiley & Sons, Ltd.     

Axisymmetric wave propagation in multilayered poroelastic grounds due to a transient acoustic point source

This paper deals with the study of axisymmetric wave propagation in various acoustic/porous stratified media coupling configurations. It presents the theoretical developments of a semi-analytical method, its validation for a limit test-case half-space ground, and an extension to a realistic multilayered seabed, when spherical waves are emitted from a transient point source in water.