A preliminary empirical model for frequency-dependent attenuation of Fourier amplitude spectra In Serbia from the Vrancea earthquakes

We present the frequency-dependent attenuation model for Fourier amplitude spectra of strong earthquake ground motion in Serbia from intermediate depth earthquakes in the Vrancea source zone in Romania. The development of this type of scaling is the essential first step toward developing the corresponding attenuation and scaling equations for pseudo relative velocity spectra (PSV), which are necessary for seismic macro- and microzoning in the territory of Serbia. Such scaling is necessary because the Vrancea source zone produces large earthquakes with shaking that attenuates differently from the local earthquakes in Serbia. Development of such a scaling model is associated with several difficulties, the principal one being the lack of recorded strong motion accelerograms at epicentral distances exceeding 300 km. To reduce uncertainties with such scaling, we require our preliminary scaling equations to be consistent with independent estimates of seismic moment, stress drop, and radiated wave energy. In the future, when the recorded strong motion data from Vrancea earthquakes increases several-fold of what it is today, it will become possible to perform this analysis again, thus leading to more reliable and permanent scaling estimates.     

非连续变形分析（DDA）线性方程组的高效求解算法

非连续变形分析(DDA)方法对大规模工程问题的数值模拟耗时太长,其中线性方程组求解耗时可占总计算时间的70%以上,因此,高效的线性方程组解法是重要研究课题。首先,阐述了适用于DDA方法的基于块的行压缩法和基于试验-误差迭代格式的非0位置记录;然后,针对DDA的子矩阵技术,将块雅可比迭代法 (BJ)、预处理的块共轭梯度法 (PCG,包括Jacobi-PCG、SSOR-PCG) 引入DDA方法,重点研究了线性方程组求解过程中的关键运算;最后,通过两个洞室开挖算例,分析了各线性方程组求解算法在DDA中的计算效率。研究表明：与迭代法相比,直解法无法满足大规模工程计算需要;BJ迭代法与块超松弛迭代法(BSOR)的效率差别不大,但明显不如PCG迭代法。因此,建议采用PCG迭代法求解DDA线性方程组,特别是SSOR-PCG值得推广;如果开展并行计算研究,Jacobi-PCG是较好的选择,当刚度矩阵惯性优势明显时,BJ迭代法同样有效。     

墙下条形基础与层状横观各向同性地基共同作用

运用对偶积分方程来求解层状横观各向同性地基与墙下条形基础的共同作用问题。从直角坐标平面应变问题控制方程出发,通过傅里叶（Fourier）变换和层间连续性条件,可以得到层状横观各向同性地基的传递矩阵解。基于该传递矩阵解,并利用条形基础与地基接触的混合边值条件,推导出一组关于基础挠度和地基反力的对偶积分方程。考虑墙下条形基础受到竖向集中荷载的情况,利用弹性薄板理论先求解出条形基础挠度;随后应用雅可比（Jacobi）正交多项式和级数展开的方法,将对偶积分方程转化为线性代数方程组进行求解。编制了相应的计算程序,其计算结果与有限元软件ABAQUS的结果基本吻合,从而验证了所提理论的正确性。算例分析表明,板土相对刚度与地基成层性对地基反力、地表沉降和沿z轴竖向正应力有很大的影响。     

隧道围岩松动圈的现场测试与理论分析

围岩松动圈是制定隧道支护措施和整体稳定性评价的重要参数和依据,如何快速准确地确定松动圈范围已备受关注。以张-石（张家口至石家庄）高速公路上的多条隧道工程为背景,采用声波探测和理论计算相结合的方法,对不同埋深、围岩等级条件下的隧道松动圈进行了研究。基于D-P准则推导了围岩应力松动圈的弹塑性理论和损伤理论计算方法,其计算结果均小于单孔声波法测试结果,但两者变化规律类似。由于考虑了岩体的峰后软化特性,损伤理论分析法计算结果比弹塑性理论更加接近于实测值,运用损伤理论对围岩（特别是低级别围岩）松动圈进行预测是可行的,研究结果对其他类似工程具有借鉴作用。     

日本东北M_W9.0地震海啸在港池及邻近区域诱发的涡流危险性计算与评估分析

海啸造成的灾害与损失并非都与淹没有关,特别是港口中海啸诱导的强流会对船只及海事设施产生重要的影响及损害.由于海啸流观测数据稀缺及海啸诱导涡流机制的不确定性,过去60年海啸科学主要集中于对海啸波特征及淹没过程的研究与分析,海啸流模拟及验证工作开展较少,导致对海啸流基本特征及其造成灾害现象的曲解.开展海啸诱导的涡流研究及预警服务显得尤为重要及紧迫.考虑快速海啸预警需要,综合对比海啸诱导涡流的物理框架及模型方法,探索兼顾效率与计算精度的海啸流模拟方法是本文的核心工作及出发点.通过分析浅层湍流相干结构(TCS)产生的主要物理耗散机制,确定了考虑2D水平耗散机制的非线性浅水方程可用于海啸涡流的模拟分析.基于高精、高分辨率有限体积模型Geoclaw建立了三个精细化的港口海啸流模型,模型分辨率为5m.利用基于海啸浮标反演的海啸源模型作为初始条件,模拟分析了日本东北地震海啸在远场的海啸波流特征.海啸波流特征模拟结果与观测吻合较好,结果可信.对比发现:波驱动的自由表面流,小的位相或波幅误差就会导致大的流速误差,流的模拟和预报相对波幅来说更具挑战性.研究了海啸波流能量在港池中的分布特征,得到:港池入口及防波堤两端常被强流控制,具有极高的危险性;相对于波幅的空间变化,海啸流具有更强的空间敏感性;所建立的高分辨率海啸模型模拟再现了日本海啸在近场的涡旋结构,给出了与观测基本一致的涡流特征.最后,引入海啸流危险等级标准,分析了港口海啸流危险性等级分布、船只疏散的安全深度及回港的时间周期.针对港口、海湾同时考虑海啸波流特征的海啸预警与评估对于港口应急管理者科学决策具有重要意义.     

一类三次超越多项式零点的分布及其在时滞生物系统的应用

本文讨论了一类三次超越多项式零点的分布,给出了零点分布在左半复平面、穿过虚轴和进入右半复平面的充分条件.作为应用,考虑了带时滞的三种群食物链模型、带时滞的睾酮分泌模型、带时滞的HIV体内感染动力学模型、带时滞的葡萄糖-胰岛素模型、带时滞的肿瘤-免疫系统作用模型的稳定性和分支问题.     

Analytical modelling of a large‐scale dynamic testing facility

An analytical model, which aims at reproducing the response of a large‐scale dynamic testing facility, that is a system composed of the specimen/shaking table/reaction‐mass/airbags/dampers/soil is developed. The Lagrangian of the system is derived, under the assumption of large displacements and rotations. A set of four nonlinear differential equations is obtained and solved with numerical methods. Preliminary verifications of the derived model are carried out by reproducing both well‐known results in the literature as well as those of a lumped model employed in the design of an existing dynamic testing facility. The case‐study for validating the nonlinear equations of motion is the shaking table of the EUCENTRE Laboratory. Copyright © 2011 John Wiley & Sons, Ltd.     

岩石的临界孔隙度反演及横波预测(英文)

临界孔隙度模型是利用岩石的临界孔隙度来计算岩石骨架的弹性模量,岩石的临界孔隙度值受到很多因素的影响,而实际应用中通常无法获得准确的临界孔隙度值,只能选取经验临界孔隙度值,就会给岩石物理建模带来误差。本文提出了一种利用纵波速度反演岩石的临界孔隙度的方法,并且把它应用于横波预测中。实验室和测井数据应用结果表明本文提出的方法可以降低以往选取经验值带来的不确定性,并且能够为横波预测提供准确的临界孔隙度值,提高了横波预测的精度。     

Numerical simulation of irregular wave overtopping against a smooth sea dike

Based on the filtered Navier-Stokes equations and Smagorinsky turbulence model,a numerical wave flume is developed to investigate the overtopping process of irregular waves over smooth sea dikes.Simulations of fully nonlinear standing wave and regular wave’s run-up on a sea dike are carried out to validate the implementation of the numerical wave flume with wave generation and absorbing modules.To model stationary ergodic stochastic processes,several cases with different random seeds are computed for each specified irregular wave spectrum.It turns out that the statistical mean overtopping discharge shows good agreement with empirical formulas,other numerical results and experimental data.