一种模拟非均匀介质中弹性波传播新的k-space方法          下载免费PDF全文

黄继伟  刘洪 《地球物理学报》2020,63(8):3091-3104

传统的伪谱(PS)方法,采用傅里叶变换(FT)计算空间导数具有很高的精度,每个波长仅需要两个采样点,而时间导数采用有限差分(FD)近似因而精度较低.当采用大时间步长时,由于时空精度不平衡,PS法存在不稳定性问题.原始的k-space方法可以有效地克服这些问题但是却无法适用于非均匀介质.为了提高原始k-space方法模拟非均匀介质波动方程的精度,我们提出了一种新的k-space算子族.它是用非均匀介质的变速度代替原k-space算子中的常数补偿速度构造得到,引入低秩近似可以高效求解.我们将构造的新的k-space算子应用于耦合的二阶位移波动方程,而不是交错网格一阶速度应力波动方程,使模拟弹性波的计算存储量减少.我们从数学上证明了基于二阶波动方程的k-space方法与基于一阶波动方程的k-space方法是等价的.数值模拟实验表明,与传统的PS、交错网格PS和原始的k-space方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本.  相似文献   

A nearly analytic symplectic partitioned Runge–Kutta method based on a locally one-dimensional technique for solving two-dimensional acoustic wave equations     

Chol Sim  Chunyou Sun  Nam Yun 《Geophysical Prospecting》2020,68(4):1253-1269

In this paper, we develop a new nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique for numerically solving two-dimensional acoustic wave equations. We first split two-dimensional acoustic wave equation into the local one-dimensional equations and transform each of the split equations into a Hamiltonian system. Then, we use both a nearly analytic discrete operator and a central difference operator to approximate the high-order spatial differential operators, which implies the symmetry of the discretized spatial differential operators, and we employ the partitioned second-order symplectic Runge–Kutta method to numerically solve the resulted semi-discrete Hamiltonian ordinary differential equations, which results in fully discretized scheme is symplectic unlike conventional nearly analytic symplectic partitioned Runge–Kutta methods. Theoretical analyses show that the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique exhibits great higher stability limits and less numerical dispersion than the nearly analytic symplectic partitioned Runge–Kutta method. Numerical experiments are conducted to verify advantages of the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique, such as their computational efficiency, stability, numerical dispersion and long-term calculation capability.  相似文献   

Bidirectional shaking table tests of a low-cost friction sliding system with flat-inclined surfaces     

Miguel B. Brito  Mitsuyoshi Akiyama  Yoshitaka Ichikawa  Hiroki Yamaguchi  Riki Honda  Naomitsu Ishigaki 《地震工程与结构动力学》2020,49(8):817-837

A novel low-cost friction sliding system for bidirectional excitation is developed to improve the seismic performance of reinforced concrete (RC) bridge piers. The sliding system is a spherical prototype developed by combining a central flat surface with an inclined spherical segment, characterized by stable oscillation and a large reduction in response accelerations on the flat surface. The inclined part provides a restoring force that limits the residual displacements of the system. Conventional steel and concrete are employed to construct a flat-inclined spherical surface atop an RC pier. The seismic forces are dissipated through the frictions generated during the sliding movements; hence, the seismic resilience of bridges can be ensured with a low-cost design solution. The proposed system is fabricated utilizing a mold created by a three-dimensional printer, which facilitates the use of conventional concrete to construct spherical shapes. The concrete surface is lubricated with a resin material to prevent abrasion from multiple input ground motions. To demonstrate the effectiveness of the system, bidirectional shaking table tests are conducted in the longitudinal and transverse directions of a scaled bridge model. The effect of the inclination angle and the flat surface size is investigated. The results demonstrate a large decrease in response acceleration when the system exhibits circular sliding displacement. Furthermore, the inclination angle that generates the smallest residual displacement is identified experimentally.  相似文献   

深厚软土场地基岩地震动反演——以天津波为例          下载免费PDF全文

梁建文  张萍  何颖 《世界地震工程》2020,(3):001-11

根据天津波地表加速度记录和相应实测场地钻孔资料，采用等效线性化方法模拟土体非线性，进行了深厚软土场地基岩地震动反演，研究了频率截断对深厚软土场地基岩地震动反演的影响，为天津波地表加速度记录提供了相应的基岩地震动输入。研究表明:不同截断频率情况下，反演得到的基岩地震动峰值和反应谱曲线形状均相差悬殊，但反演后再正演所得到的地面运动与原地表记录相比，加速度时程曲线形状和峰值大小、反应谱曲线形状和峰值大小均非常接近；可根据工程需要对地表加速度记录进行不同频率的截断，反演得到所需强度的基岩地震动。  相似文献   

Experimental investigation on a cabin-suspended catamaran in terms of motion reduction and wave energy harvesting by means of a semi-active motion control system     

《Applied Ocean Research》2019

A novel concept catamaran equipped with a suspended cabin, named Wave Harmonizer Type 4 (WHzer-4), is proposed and evaluated. The mass-spring-mass system is constructed by mounting four sets of suspensions in-between the cabin and the twin-hull. Two sets of dual motor/generators (M/Gs) are attached on the center beam of the cabin's deck fore and aft. Each shaft-end of the dual M/Gs is connected to the twin-hull through a rack-pinion gear unit. In this way the vertical relative motion between the cabin and the twin-hull can be transferred into the rotational motion of the M/Gs, and vice versa. A semi-active motion control system, which contains a proportional-integral (PI) controller, is designed and applied to each of the dual M/Gs for the aim of absorbing wave energy under the condition of suppressing the local vertical velocity of the cabin as much as possible. A 1/5 scale model ship with a length of 1.6 m is built, and a forced-oscillation bench test is implemented to validate the performance of the control system. Then, a series of towing tank tests is carried out in regular head waves. The heave and pitch responses of the cabin, those of the twin-hull and the corresponding wave energy capture width ratio (CWR) at five control scenarios and two reference scenarios are investigated. Discussion on the results of the tank test shows that the motion reduction of the cabin and the wave energy harvesting can be achieved simultaneously at a few wave conditions. However, at other conditions, although noticeable amount of wave energy is harvested, motion reduction of the heave and pitch of the cabin could not be obtained at the same time. It is suggested that varying the gain settings of the PI controllers according to the location of the controllers may improve the effectiveness of the proposed control system.  相似文献   

Hydrodynamic performance of a T-shaped floating breakwater     

《Applied Ocean Research》2019

The comprehensive utilization of floating breakwaters, specially acting as a supporting structure for offshore marine renewable energy explorations, has received more and more attention recently. Based on linear water-wave theory, the hydrodynamic performance of a T-shaped floating breakwater is semi-analytically investigated through the matched eigenfunction expansion method (MEEM). Auxiliary functions, to speed up the convergence and improve the accuracy in the numerical computations, are introduced to represent the singular behavior of fluid field near the lower salient corners of the structure. The effects of the height and installation position of the vertical screen on the reflection and transmission coefficients, dynamic response and wave forces are examined. It is found that the presence of the screen shifts the resonance frequency of RAO for both surge and pitch modes to the low-frequency area, while has no effect on heave mode. The identical added masses, damping and transmission coefficients can be obtained in the cases where the screen holds the same distance away from the longitudinal central axis of the upper box-type structure. Moreover, a relatively small pitch response can be achieved in a wide wave–frequency range, when the breakwater is Γ-shaped.  相似文献   

Modelling of viscoacoustic wave propagation in transversely isotropic media using decoupled fractional Laplacians     

Zhihao Qiao  Chengyu Sun  Jie Tang 《Geophysical Prospecting》2020,68(8):2400-2418

A new wave equation is derived for modelling viscoacoustic wave propagation in transversely isotropic media under acoustic transverse isotropy approximation. The formulas expressed by fractional Laplacian operators can well model the constant-Q (i.e. frequency-independent quality factor) attenuation, anisotropic attenuation, decoupled amplitude loss and velocity dispersion behaviours. The proposed viscoacoustic anisotropic equation can keep consistent velocity and attenuation anisotropy effects with that of qP-wave in the constant-Q viscoelastic anisotropic theory. For numerical simulations, the staggered-grid pseudo-spectral method is implemented to solve the velocity–stress formulation of wave equation in the time domain. The constant fractional-order Laplacian approximation method is used to cope with spatial variable-order fractional Laplacians for efficient modelling in heterogeneous velocity and Q media. Simulation results for a homogeneous model show the decoupling of velocity dispersion and amplitude loss effects of the constant-Q equation, and illustrate the influence of anisotropic attenuation on seismic wavefields. The modelling example of a layered model illustrates the accuracy of the constant fractional-order Laplacian approximation method. Finally, the Hess vertical transversely isotropic model is used to validate the applicability of the formulation and algorithm for heterogeneous media.  相似文献   

Dynamic loading tests and seismic response analysis of a stud-type damper composed of multiple friction units with disc springs     

Takeshi Sano  Kazutaka Shirai  Yasumasa Suzui  Yoshikazu Utsumi 《地震工程与结构动力学》2020,49(13):1259-1280

This paper presents an experimental study on the performance of a shear-sliding stud-type damper composed of multiple friction units with high-tension bolts and disc springs. A numerical evaluation of the response reduction effects achieved by the stud-type damper is also presented. In dynamic loading tests, the behavior of stud-type multiunit friction damper specimens was investigated. Three different full-scale damper specimens, which were composed of five, six, or seven friction units with two or four sliding surfaces, were incorporated into loading devices for testing. The stud-type friction dampers demonstrated stable rigid-plastic hysteresis loops without any remarkable decrease in the sliding force even when subjected to repetitive loading, in addition to showing no unstable behavior such as lateral buckling. The damper produced a total sliding force approximately proportional to the number of sliding surfaces and friction units. The total sliding force of the stud-type damper can thus be estimated by summing the contributions of each friction unit. In an earthquake response simulation, the control effects achieved by stud-type dampers incorporated into an analytical high-rise building model under various input waves, including long-period, long-duration and pulse-like ground motions, were evaluated. A satisfactory response reduction was obtained by installing the developed stud-type dampers into the main frame without negatively impacting usability and convenience in terms of building planning.  相似文献   

Impact of fractional probability distributions on statistics of hydraulic conductivity,dynamics of groundwater flow and solute transport at a low-permeability site     

Chengpeng Lu  Wei Qin  Ching-Sheng Huang  Yong Zhang  Tema Koketso Ealotswe 《水文研究》2020,34(21):4112-4127

The fractional Brownian motion (fBm) and fractional Lévy motion (fLm) can easily describe the geometry and the statistical structure of hydraulic conductivity (K) for real-world. However, the fBm and fLm models have not been systematically evaluated when building the K field for a low-permeability site. In this study, both the fBm and fLm are used to simulate the low-K field at NingCheGu (NCG), Tianjin, China. Groundwater flow and solute transport are then computed using MODFLOW and MT3DMS, respectively, and the influence of the fBm/fLm models for K on groundwater flow and solute transport is discussed. Results show that the fLm fits better the statistics of the low-K medium than fBm, and the random logarithmic K (LnK) field generated by fLm is more stable because the resultant LnK field captures more of the measured properties at the field site than that generated by fBm. In contrast, the LnK generated by fBm is more likely to form both high-K channels and low-K barriers. The fBm therefore predicts more extreme behaviours in flow and transport, including the preferential flow, low-concentration blocks and solute retention. The overall groundwater renewal period and solute travel time for the fLm simulation are slightly shorter than those for fBm. The impacts of the fLm and fBm models on the statistics of the resultant LnK fields and the dynamics of groundwater flow and solute transport revealed by this study shed light on the selection and evaluation of the fractional probability distribution models in capturing the K fields for low-K media.  相似文献   

Efficient strong motion duration of pulse-like records for nonlinear structural analyses     

Constantinos C. Repapis  Petros P. Mimoglou  Vicky V. Dimakopoulou  Ioannis N. Psycharis  Ioannis M. Taflampas 《地震工程与结构动力学》2020,49(5):479-497

The special interest produced by near-field directivity records and their effect on structural response has given a new significance in the velocity time history, its pulse-like content, and relevant parameters and indices. Recent research has shown that directivity pulses inherent in these records govern the linear and the nonlinear response of a wide range of structures. Based on this observation, it is suggested in this paper that a truncated ground motion, limited to the duration of the predominant velocity pulse, can be efficiently used to predict the structural response, instead of the base motion with the total duration, reducing significantly the required runtimes. The proposed methodology is verified for a series of medium to high rise reinforced concrete buildings, for which nonlinear time-history analyses are performed for a vast suite of pulse-like near-field records applied as base excitations with their total duration and the proposed truncated one. Comparison of the results for the response displacements and forces shows very good agreement, permitting the acceptance of the pulse duration as the efficient strong motion time interval of the original record, which determines the response and, thus, it can be used for nonlinear structural analyses.  相似文献