Vibration of cylindrical shells with clamped ends

A general method is outlined for the determination of natural frequencies of cylindrical shells with any boundary conditions when the effects of rotatory inertia and transverse shear deformation are included in the analysis. This is applied to cylindrical shells with both ends clamped. It is shown that the inclusion of these effects tends to have a greater effect upon frequencies of cylindrical shells with clamped ends than it does for corresponding shells with simply supported ends, for which numerical results are available. The authors suggest an empirical relation, which together with the latter results enables rapid estimates to be made of the effects of rotatory inertia and shear deformation on the frequencies of a wide range of cylindrical shells with clamped ends. An assessment of the accuracy of the theory with these effects included is made by comparing frequencies with values from a three-dimensional elasticity theory, but this comparison has to be restricted to cylindrical shells with simply supported ends.     

The algorithms for calculating synthetic seismograms from a dipole source using the derivatives of Green’s function

The problem of calculating complete synthetic seismograms from a point dipole with an arbitrary seismic moment tensor in a plane parallel medium composed of homogeneous elastic isotropic layers is considered. It is established that the solutions of the system of ordinary differential equations for the motion–stress vector have a reciprocity property, which allows obtaining a compact formula for the derivative of the motion vector with respect to the source depth. The reciprocity theorem for Green’s functions with respect to the interchange of the source and receiver is obtained for a medium with cylindrical boundary. The differentiation of Green’s functions with respect to the coordinates of the source leads to the same calculation formulas as the algorithm developed in the previous work (Pavlov, 2013). A new algorithm appears when the derivatives with respect to the horizontal coordinates of the source is replaced by the derivatives with respect to the horizontal coordinates of the receiver (with the minus sign). This algorithm is more transparent, compact, and economic than the previous one. It requires calculating the wavenumbers associated with Bessel function’s roots of order 0 and order 1, whereas the previous algorithm additionally requires the second order roots.     

Models of the flexure-controlled strength,stiffness and cyclic deformation capacity of rectangular RC columns with smooth bars,including lap-splicing and FRP jackets

Cyclic tests of single concrete columns with smooth (plain) bars are not representative of building columns with lap splices at floor levels and story-long starter bars. Column specimens with fixity at top and bottom resemble building columns best, but few of those tested so far had smooth bars and even then without bar lap-splicing at floor level or FRP jackets at column ends. Empirical models based on single-column tests, especially the numerous ones with cantilever-type specimens, cannot be readily extended to columns with smooth bars in real-life buildings. Physical models of the Strut-and-Tie type are developed and are validated or calibrated through comparisons with laboratory tests. Their scope includes anchorage and splicing of bars with either 180° hooks or straight ends. Once validated, they are adapted to real-life multistory rectangular RC columns with smooth bars, in order to obtain the column properties of interest: the chord rotation at yielding and the cyclic ultimate chord rotation, with or without FRP jacketing. Different expressions apply to the top and bottom end of a column in a story, but a single one is used to estimate the column’s effective stiffness. Empirical alternatives fitted to the single-element test results have slightly less scatter than physical models, but caution is needed for their application to columns of real buildings. Simulations of the 3D seismic response of a plan wise asymmetric full size building, tested pseudo-dynamically before or after retrofitting all columns with FRPs or just two of them with RC jackets, provide certain confidence in the extension of the physical models for the estimation of the stiffness and ultimate deformation of columns with smooth bars in real-life buildings.     

估计辐射参数的井间电磁波层析成像技术

为了避免使用不合理初始辐射场强和方向性因子带来的误差,研究了估计辐射参数的井间电磁波层析成像技术.通过时域有限差分法模拟表明,天线长度与波长的比值、钻孔充填情况、钻孔周围介质的物性均会影响偶极天线的初始辐射场强或方向性因子;为此结合已知的分层资料,将它们设为未知参数,并设定初始辐射场强与发射点位置相关,方向性因子随射线角度而变化;采用正则化反演方法,由钻孔资料建立了模型方差目标函数,使得反演结果与钻遇的地质特征保持一致.通过理论模型试验和实例应用分析表明,相对于传统射线层析成像方法,估计辐射参数的正则化层析成像技术有助于提高反演的准确性.     

Engineering ground motion record selection in the ITalian ACcelerometric Archive

This paper introduces REXELite, an internet version of REXEL, a software for automatic selection of ground motion suites for nonlinear dynamic analysis of structures. REXELite was developed with the aim of integrating an advanced earthquake records’ repository, such as the ITalian ACcelerometric Archive (ITACA), with a tool to define seismic input for engineering seismic analysis according to international standards (with priority to Europe). In fact, REXELite allows to define target design spectra according either to Eurocode 8 or to the Italian building code, and to search ITACA for suitable sets of seven records (comprised of one or two horizontal ground motion components) matching such target spectra: on average, in a user-specified period range, and with the desired tolerance. The records in the set also have, individually and according to some criteria, the most similar spectral shape with respect to that of the code. Selection options include magnitude, source-to-site distance, soil conditions and, if desired, linear scaling of records to reduce further record-to-record variability of the selected suite.     

1997～1998年伽师强震群小地震震源机制解系统聚类分析

用P、S波最大振幅比法计算1997~1998年伽师强震群MS2.5~3.7小地震震源机制解,对小地震震源机制解进行系统聚类分析。结果表明,伽师小地震主要以走滑和正断层为主,与该时段发生的伽师6级强震具有相似的断错性质。1997年伽师小地震主压应力P轴方位为北东或近南北向,与1997~1998年6级强震主压应力P轴方位一致。1997年4月6日6.4级和1998年8月2日6.1级两次6级正断层地震前出现地震空段现象。     

芦山7.0级地震近断层地震动的方向性

基于芦山7.0级地震中断层距小于100 km自由场台站的强震动记录观测数据,研究此次地震近断层地震动的方向性特性,并探讨方向性特性与震源破裂机制、断层距离和空间方位的关系.研究结果表明：（1）与距断层较远记录不同,近断层地震动在不同的观测方向上表现出显著的强度差异,存在明显的极大和极小作用方向.在不同的方向上,最大加速度反应可达最小值的4倍以上;（2）这种方向性差异在T=1.0 s以上的长周期段更为明显,在T=0.1 s以下的短周期段,地震动随方向变化的差异较小.地震动强度随方向变化的差异随周期增大而增大,不同方向上的加速度反应谱值的最大值与最小值之比从周期T=0.01 s时的约1.7增大到周期T=10 s时的约2.4;（3）在距离断层约35 km以内,地震动具有明显方向性,地震动卓越方向具有垂直断层走向的特征,随断层距的增大,这种方向性不明显.从不同方向上地震动强度的差异来看,随断层距增大,地震动强度在不同方向上的差异在减小,表现为各个周期的最大值/中值和最大值/最小值比值均随断层距离增大缓慢减小;（4）近断层地震动的方向性特性主要受断层上、下盘的相对运动所控制,其在长周期的卓越方向与水平同震位移方向一致,且该卓越方向上的地震动强度绝对大小与地震破裂造成的静态位移明显相关,表现为地震动强度随水平同震位移的增大而增大.     

武汉市区土体剪切波速与深度的经验关系

通过对武汉市区3个主要地质单元共8 305个剪切波速数据的分析整理,分别运用线性函数、一元二次多项式函数、指数函数对武汉市区不同地貌单元不同土类的剪切波速与深度的关系进行统计回归,得到其经验关系。结果表明,武汉市区土体剪切波速与埋深相关性比较明显;一元二次多项式函数的拟合效果最好。将实测数据与利用经验关系得到的预测值进行对比检验,两者基本吻合,可供武汉市区场地剪切波速数据缺乏时参考使用。     

MODAL ANALYSIS OF SOFT-SOIL SITES INCLUDING RADIATION DAMPING

For the one-dimensional analysis of soft-soil layers on an elastic half-space, a general form of analytical solution is developed for converting radiation damping due to energy leaking back to the half-space into equivalent modal damping, allowing the modal analysis technique to be extended to a site where radiation damping has to be accounted for. Closed-form solutions for equivalent modal damping ratios and effective modal participation factors are developed for a single layer with a shear wave velocity distribution varying from constant to linearly increasing with depth. Compact and recursive forms of solutions for equivalent modal damping ratios are developed for a system with an arbitrary number of homogeneous layers on an elastic half-space. Comparisons with numerical solutions show that the modal solutions are accurate. The nominal frequency of a site, i.e. the inverse of four times the total shear wave travel time through the layers, is an important parameter for estimating the high mode frequencies. A parameter study shows that for the same impedance ratio of the bottom layer to the elastic half-space, a system of soil layers with an increasing soil rigidity with depth has, in general, larger peak modal amplifications at the ground surface than does a single homogeneous layer on an elastic half-space, while a system with a decreasing soil rigidity with depth has smaller modal peak amplifications. © 1997 by John Wiley & Sons, Ltd.     

多分辨率地震信号反褶积

基于二进小波变换提出了一种新的反褶积方法─-多分辨率地震信号反褶积．在地震信号二进小波变换域中的各尺度上分别进行其分辨率随小波尺度变化的反褶积,利用不同分辨率反褶积结果之间的相关性,以及测量噪声随尺度的衰减特性,从低分辨率反褶积结果逼近高分辨率反褶积结果．理论分析和实验表明,该方法有较高的精度,并且在较低信噪比情况下有好的效果．     

基于有限差分正演的带地形三维大地电磁反演方法

本研究实现了一套基于有限差分（FD）方法的大地电磁测深数据带地形三维反演算法及代码.其中,在大地电磁场正演数值模拟方面,开发了起伏地形条件下基于交错网格剖分、有限差分方法的大地电磁测深三维正演代码;在满足平面波场假设的前提下,使用长方体网格剖分模拟三维起伏地形,实现了带地形三维正演计算;并设计理论模型进行试算,经试算结果与前人的有限元法计算结果对比,验证了所研发的带地形三维正演计算的正确性与可靠性.在反演方面,本研究基于非线性共轭梯度方法编写了大地电磁测深带地形三维反演代码,试验了不同的共轭梯度搜索因子β,避免了目标函数对海森矩阵（参数二次导数矩阵）的显式计算和存储,初步实现了大地电磁资料的带地形三维反演.最后,对一系列理论模型进行正演计算,利用其生成的合成数据模拟实测数据进行反演,并与现有的不带地形大地电磁测深三维反演结果比较,检验了所研发的带地形三维反演计算的可靠性与稳定性.     

Heating-induced changes in the anisotropy of magnetic susceptibility of impact breccias, Chicxulub Crater (Mexico)

Initial results of a thermal treatment study on the anisotropy of magnetic susceptibility (AMS) of impact breccias from Chicxulub crater are used to investigate the nature of the magnetic fabrics. Chicxulub impact breccias are heterogeneous materials, with carbonate, basement and melt clasts within carbonate-rich or melt-rich matrix. Samples studied come from the carbonate-rich basal unit Lower Suevite in the Yaxcopoil-1 borehole impactite sequence (core depth interval: 885?C895 m). The Lower Suevite is characterized by mixed prolate and oblate ellipsoids with shallow to steep principal susceptibility axes, which had been related to emplacement as an excavation flow with ground-surge components during the early cratering stages. Thermal treatment results in changes in the fabrics with a tendency to oblate fabrics. Stepwise thermal treatment up to 700°C reveals different behaviors for the oblate, neutral and prolate fabrics marked by changes in AMS parameters and principal susceptibility axis orientations. A sample with oblate fabrics and vertical minimum axes showed an increase of magnetic susceptibility at high temperatures, indicating formation of secondary magnetite and fabric enhancement. A sample with neutral ellipsoid showed heating-induced changes towards oblate fabrics and vertical minimum susceptibility axes. Samples characterized by prolate ellipsoids with horizontal maximum axes showed no directional changes. In a sample with apparent intermediate or inverse fabrics, vertical maximum axes showed changes to horizontal inclinations, with the intermediate and maximum axes switching positions. Changes induced by stepwise thermal treatment appear useful to characterize the fabrics of impact lithologies. Further investigation of heating-induced effects in mineralogy, grain size and textural changes is, however, required to relate the different behaviors observed after stepwise thermal treatment with the magnetic mineralogy and emplacement mode of the breccias.     

Cyanobacterial crust and soil particle detachment: a rain‐chamber experiment

Cyanobacteria are known to form a crust on soil surfaces holding soil particles together and thereby offering resistance to erosion. A controlled experiment was carried out to throw light on this issue. The experiment consisted of subjecting erosion cups filled with soil to artificial rainfall in the laboratory. Three sets of erosion cups, each set consisting of six, were used. One set consisted of soil with inoculated cyanobacteria and the second set consisted of soil with naturally colonized cyanobacteria, both over a period of about 8 months. The third set consisted of soil with no bacterial growth. The results indicate that the soil erosion cups with the inoculated cyanobacterial crust had at least one order of magnitude less erosion for coarse soils, and about two orders of magnitude less erosion for fine soils, compared with erosion cups with no bacterial crust. Copyright © 2002 John Wiley & Sons, Ltd.     

Theory of the drift of asymmetrically-heated lithospheric plates

An approximate solution to the problem of the drift of a raft formed of two parallel wires, differentially heated, rigidly coupled and floating in a fluid of finite thickness and with linear viscosity, has been obtained and is shown to agree well with experiment up to a multiplicative constant. In addition, the solution to the problem of the drift of a solid raft with a suspended single-wire heater also shows good agreement with experiment up to a multiplicative constant. In both cases, the rafts drift with constant velocity. For small amounts of heat, the drift velocity is proportional to the first power of the heat input; for large amounts of heat, the drift velocity is proportional to the square root of the heat input. Within imponderable factors of an order of magnitude, the drift velocities are appropriate for drift of lithospheric plates containing both an oceanic and a continental part.     

四川地区强震前弱震空区的变化特征

本文使用四川台网地震目录全部资料,适当参考邻省地震目录,对于四川省从1970年下半年以来发生的15组5级以上地震(另有5组缺乏完整资料无法讨论),作 R-t 图、D-t 图寻找时间空区,并作相应的平面空区图.通过对比分析,发现6组6级以上和4组5级地震前震中附近出现空区现象.根据形成空区的变化过程可分为先扩(大)后缩(小)型和逐渐扩大的开放型两种类型.平面空区的位置和形状往往与构造展布有关.有5组5级震例,或因震前震中区小震不断,或因震中及其附近地区小震太少,看不出明显的空区.本文利用出现空区的10个震例,总结出空区参数与未来地震的震级的经验关系:M=5.50 lg Rmax-2.331.32 M=2.98 lgT+2.941.06 M=2.52lg(RmaxT)-0.510.96文中将上述结果与大竹政和的结果作了对比.空区不同变化形态可用金森博雄概括的突起点模式加以解释.      

Analysis of the seismicity in the central part of the Pyrenees (France), and tectonic implications

Two temporary seismological networks have been set up in 2000 and 2002 in the central part of the Pyrenees, in a region, which appears as a transition between two domains where both the seismic activity and the tectonic regime are different. Together with the permanent networks, they allowed us to obtain precise hypocenter locations for more than 400 events with local magnitudes ranging from 1.5 to 4.6, as well as 30 new focal mechanisms. The seismicity is distributed in several clusters, which are not located along the North Pyrenean Fault, considered as the major tectonic accident resulting from the suture of the Iberian and Eurasian convergent plates when the range formed. Several small fault segments dipping to the north are identified. The maximum focal depth varies from 10 to 20 km, with variations which are roughly parallel to those of the Moho, indicating a thickening of the seismogenic layer to the east of the studied area. The obtained focal solutions reveal a predominance of normal faulting to the West and reverse faulting to the East, with strike-slip motions in between. The largest fault segment to the East, with a length of about 20 km, could possibly be related to a large historical event which occurred in 1660, with intensity IX, close to cities which have become since then important touristic centres.     

基于OpenSees的混凝土结构连续倒塌分析中的从属节点建模技术

对重大工程结构进行强地震作用下的连续倒塌全过程分析并建立相应的设计与控制方法,已成为当前地震工程领域的发展趋势。目前,由于数值求解方面的困难,绝大多数针对极端作用下的结构连续倒塌的研究止步于数值临界状态的界定,在分析过程中不能实时地对结构构件的损伤状态进行监测并根据构件的损伤状态对分析模型进行修改。为了实现连续倒塌过程中构件的逐步失效,在OpenSees程序中,基于Beam with Hinges Element构建了端部带附属节点的Beam with Hinges Element,并根据构件失效情况对附属节点的多点约束进行控制。采用三次静力凝聚方法和Newmark-beta法,给出了对此改进建模技术可信的理论背景,方法的准确性通过一个简单框架算例得到了验证。     

流体对石灰岩断层摩擦滑动影响的实验研究

在气体介质三轴高温岩石力学实验仪器上,采用意大利Scaglia Bianca石灰岩,在温度50～300℃、围压150MPa,含50MPa孔隙压、无孔隙压含饱和水和完全干燥三种条件下,开展摩擦滑动实验.实验力学数据和显微结构表明,完全干燥样品在120℃时出现慢滑移,实验样品中没有出现溶解与沉淀.无孔隙压含饱和水条件下,100℃、120℃、150℃条件下出现典型的慢滑移,实验样品中含有微弱的溶解与沉淀;300℃条件下出现黏滑,实验样品中出现沉淀.在含50MPa孔隙压条件下,50℃时的实验表现为典型的稳滑,实验样品中以溶解为主;在100～150℃时,出现慢滑移,实验样品中以溶解为主,沉淀为辅;在200～300℃时,出现典型黏滑,实验样品中以沉淀为主.实验结果表明,石灰岩断层泥摩擦滑动稳定性随温度变化,受流体中碳酸钙的溶解和沉淀作用控制,因此,流体中矿物质的饱和度这一化学性质对断层带的摩擦强度和摩擦滑动稳定性具有显著影响.     

Non‐uniqueness with refraction inversion – the Mt Bulga shear zone

The tau‐p inversion algorithm is widely employed to generate starting models with many computer programs that implement refraction tomography. However, this algorithm can frequently fail to detect even major lateral variations in seismic velocities, such as a 50 m wide shear zone, which is the subject of this study. By contrast, the shear zone is successfully defined with the inversion algorithms of the generalized reciprocal method. The shear zone is confirmed with a 2D analysis of the head wave amplitudes, a spectral analysis of the refraction convolution section and with numerous closely spaced orthogonal seismic profiles recorded for a later 3D refraction investigation. Further improvements in resolution, which facilitate the recognition of additional zones with moderate reductions in seismic velocity, are achieved with a novel application of the Hilbert transform to the refractor velocity analysis algorithm. However, the improved resolution also requires the use of a lower average vertical seismic velocity, which accommodates a velocity reversal in the weathering. The lower seismic velocity is derived with the generalized reciprocal method, whereas most refraction tomography programs assume vertical velocity gradients as the default. Although all of the tomograms are consistent with the traveltime data, the resolution of each tomogram is comparable only with that of the starting model. Therefore, it is essential to employ inversion algorithms that can generate detailed starting models, where detailed lateral resolution is the objective. Non‐uniqueness can often be readily resolved with head wave amplitudes, attribute processing of the refraction convolution section and additional seismic traverses, prior to the acquisition of any borehole data. It is concluded that, unless specific measures are taken to address non‐uniqueness, the production of a single refraction tomogram that fits the traveltime data to sufficient accuracy does not necessarily demonstrate that the result is either correct, or even the most probable.     

一个中国大陆地区7级以上地震判定指标被再次验证

笔者曾利用统计显著性检验方法给出了中国大陆地区5级以上地震年频次与7级以上地震发生关系的判定指标,他们是：①当5级以上地震年频次N5＜10时,其后第6年我国大陆地区将发生7级以上地震;②当5级以上地震年频次N5＞31时,其后第6年我国大陆地区发生7级以上地震的可能性不大.2007年,中国大陆地区5级以上地震年频次为6次,预示2013年中国大陆地区可能会发生7级以上地震.2013年4月20日四川芦山7.0级地震再次使前一个判定指标得到验证.