共查询到20条相似文献，搜索用时 31 毫秒
1.
An exact stiffness matrix method is presented to evaluate the dynamic response of a multi-layered poroelastic medium due to time-harmonic loads and fluid sources applied in the interior of the layered medium. The system under consideration consists of N layers of different properties and thickness overlying a homogeneous half-plane or a rigid base. Fourier integral transform is used with respect to the x-co-ordinate and the formulation is presented in the frequency domain. Fourier transforms of average displacements of the solid matrix and pore pressure at layer interfaces are considered as the basic unknowns. Exact stiffness (impedance) matrices describing the relationship between generalized displacement and force vectors of a layer of finite thickness and a half-plane are derived explicitly in the Fourier-frequency space by using rigorous analytical solutions for Biot's elastodynamic theory for porous media. The global stiffness matrix and the force vector of a layered system is assembled by considering the continuity of tractions and fluid flow at layer interfaces. The numerical solution of the global equation system for discrete values of Fourier transform parameter together with the application of numerical quadrature to evaluate inverse Fourier transform integrals yield the solutions for poroelastic fields. Numerical results for displacements and stresses of a few layered systems and vertical impedance of a rigid strip bonded to layered poroelastic media are presented. The advantages of the present method when compared to existing approximate stiffness methods and other methods based on the determination of layer arbitrary coefficients are discussed. 相似文献
2.
Three-dimensional time-harmonic Green's functions of saturated soil under buried loading 总被引：3，自引：0，他引：3
Three-dimensional time-harmonic response of a poroelastic half space subjected to an arbitrary buried loading is investigated. The analysis starts with the field equations in cylindrical coordinates based on Biot's general theory of poroelasticity. General solutions for the displacements are first derived using the Fourier expansions and Hankel integral transform with respect to the circumferential and radial coordinates, respectively. The transformed-domain solutions are obtained in explicit form. The physical-domain displacements and stress components are then obtained numerically by inverse integral transform. Comparisons illustrating the accuracy of the developed approach are made with existing solutions for an elastic half space, which is reduced directly from the general solution developed in the paper. Numerical results are presented for the displacements of a saturated soil subjected to a horizontal internal excitation. 相似文献
3.
A. K. Sarker M. Sc. 《Pure and Applied Geophysics》1963,55(1):42-52
Summary In this paper, surface displacements of theSH type of motion due to time-dependent body forces and due to sudden introduction of a discontinuity in the shearing stress in a semi-infinite, homogeneous, isotropic, viscoelastic medium have been worked out by the use of Laplace and Fourier transforms. The displacements due to a point force and due to a stress discontinuity over a fixed region have been calculated numerically. 相似文献
4.
The dynamic response of a double infinite beam system connected by a viscoelastic foundation under the harmonic line load is studied. The double infinite beam system consists of two identical and parallel beams, and the two beams are infinite elastic homogeneous and isotropic. A viscoelastic layer connects the two beams continuously. To decouple the two coupled equations governing the response of the double infinite beam system, a variable substitution method is introduced. The frequency domain solutions of the decoupled equations are obtained by using Fourier transforms as well as Laplace transforms successively. The time domain solution in the generalized integral form are then obtained by employing the corresponding inverse transforms, i.e. Fourier transform and inverse Laplace transform. The solution is verified by numerical examples, and the effects of parameters on the response are also investigated. 相似文献
5.
A time domain boundary element in a cylindrical coordinate system is developed for the analysis of wave propagation in a half space. The integral formulation is based on Graffi's dynamic reciprocal theorem and Stokes' fundamental solutions. The field quantities (displacements and tractions) are expressed as products of Fourier series in the tangential direction and linear polynomials in the other spatial directions. Gaussian integration is used to integrate the non-singular parts of the integral equations, whereas the integration of the singular components, which are either of order 1/r or 1/r2, is handled by special numerical schemes. In the time marching aspect, the field quantities are assumed to vary linearly in the temporal direction as well. Examples for wave propagation due to various forms of surface excitations are reported to demonstrate the accuracy of the method. 相似文献
6.
Philip Underwood 《地震工程与结构动力学》1973,2(3):219-233
The frequency versus wave number characteristics of four O(ΔS2) finite difference formulations for one-dimensional linear shell (ring and axisymmetric) equations are investigated and compared with the exact continuum characteristics. It is found that three of the formulations give virtually identical results. These are half-spacing techniques with equilibrium in terms of displacements or resultants and whole-spacing with equilibrium in terms of displacements. The formulation based on whole-spacing with equilibrium in terms of resultants produces some dramatically different results. These discrepancies partially explain some late time instability problems and critical time step behaviour that have been reported by other investigators. 相似文献
7.
Three-dimensional simulation of track on poroelastic half-space vibrations due to a moving point load 总被引：1，自引：0，他引：1
An analytical approach is used to investigate dynamic responses of a track system and the poroelastic half-space soil medium subjected to a moving point load under three-dimensional condition. The whole system is divided into two separately formulated substructures, the track sub-system and the ground. The ballast supporting rails and sleepers is placed on the surface of the ground. The rail is modeled by introducing the Green function for an infinitely long Euler beam subjected to the action of the moving point load and the reaction of sleepers represented by a continuous mass. Using the double Fourier transform, the governing equations of motion are then solved analytically in the frequency–wave-number domain. The time domain responses are evaluated by the inverse Fourier transform computation for a certain load velocities. Computed results show that dynamic responses of the soil medium are considerably affected by the fluid phase as well as the load velocity. 相似文献
8.
9.
Irene M. Moroz 《地球物理与天体物理流体动力学》2013,107(3-4):313-314
Abstract The model equations describing two-dimensional thermohaline convection of a Boussinesq fluid in a rotating horizontal layer are known to support multiple instabilities, depending on the values of certain control parameters (Arneodo et al., 1985). Most of these multiple instabilities have already been studied for double or triple diffusive convection, where behaviours ranging from simple steady to irregular motions have been found. Here we consider the one remaining bifurcation mentioned by Arneodo et al. (1985): the interaction between a steady and an oscillatory convection roll when the linear spectrum for a single wavenumber comprises one zero and one pair of purely imaginary eigenvalues. The method of centre manifolds and normal forms is used to derive evolution equations for the amplitudes of the convection rolls close to bifurcation and the behaviours associated with the equations is discussed. 相似文献
10.
This paper proposes a new approach for forecasting continuous indoor air quality time series and in particular the concentration of a common air pollutant in offices like formaldehyde. Forecasting is achieved through the combination of the spectral band decomposition using fast Fourier transform and nonlinear time series modeling. Two nonlinear models have been tested: a threshold autoregressive (TAR) model and a Chaos dynamics-based modeling. This study shows the benefit of the Fourier decomposition coupled with nonlinear modeling of each extracted component, compared to forecasting applied directly on the raw data. Both TAR and Chaos dynamics models are able to reproduce nonlinearities, with slightly better performance in the case of the second model. These hybrid models provide good performance on forecast time horizon up to 12 h ahead. 相似文献
11.
Selection of the wavenumbers k using an optimization method for the inverse Fourier transform in 2.5D electrical modelling 总被引：1，自引：0，他引：1
An optimization method is used to select the wavenumbers k for the inverse Fourier transform in 2.5D electrical modelling. The model tests show that with the wavenumbers k selected in this way the inverse Fourier transform performs with satisfactory accuracy. 相似文献
12.
A. D. Sneyd 《地球物理与天体物理流体动力学》2013,107(1-3):141-151
Abstract Formation of electric current sheets in the corona is thought to play an important role in solar flares, prominences and coronal heating. It is therefore of great interest to identify magnetic field geometries whose evolution leads to variations in B over small length-scales. This paper considers a uniform field B 0[zcirc], line-tied to rigid plates z = ±l, which are then subject to in-plane displacements modeling the effect of photospheric motion. The force-free field equations are formulated in terms of field-line displacements, and when the imposed plate motion is a linear function of position, these reduce to a 4 × 4 system of nonlinear, second-order ordinary differential equations. Simple analytic solutions are derived for the cases of plate rotation and shear, which both tend to form singularities in certain parameter limits. In the case of plate shear there are two solution branches—a simple example of non-uniqueness. 相似文献
13.
The pseudodynamic (PSD) test method imposes command displacements to a test structure for a given time step. The measured restoring forces and displaced position achieved in the test structure are then used to integrate the equations of motion to determine the command displacements for the next time step. Multi‐directional displacements of the test structure can introduce error in the measured restoring forces and displaced position. The subsequently determined command displacements will not be correct unless the effects of the multi‐directional displacements are considered. This paper presents two approaches for correcting kinematic errors in planar multi‐directional PSD testing, where the test structure is loaded through a rigid loading block. The first approach, referred to as the incremental kinematic transformation method, employs linear displacement transformations within each time step. The second method, referred to as the total kinematic transformation method, is based on accurate nonlinear displacement transformations. Using three displacement sensors and the trigonometric law of cosines, this second method enables the simultaneous nonlinear equations that express the motion of the loading block to be solved without using iteration. The formulation and example applications for each method are given. Results from numerical simulations and laboratory experiments show that the total transformation method maintains accuracy, while the incremental transformation method may accumulate error if the incremental rotation of the loading block is not small over the time step. A procedure for estimating the incremental error in the incremental kinematic transformation method is presented as a means to predict and possibly control the error. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
14.
Fourier transform techniques have been used to calculate the theoretical filter (amplitude) response function of Nth order vertical derivative continuation operation. The amplitude response functions of the vertical gradient and its continuation follow from the same. These response functions are subsequently used to calculate the weighting coefficients suitable for two dimensional equispaced data. A shortening operator has been incorporated to limit the extent of the operator. For comparative study, some of the developed coefficient sets and the one presented in this paper are analysed in the frequency domain and their merits and demerits are discussed. 相似文献
15.
基于孔隙介质的Biot理论,首先利用Laplace变换,给出圆柱坐标系下横观各向同性饱和弹性多孔介质在变换域上的波动方程;将波动方程解耦后,根据方位角的Fourier展开和径向Hankel变换,求解了Biot波动方程,得到以土骨架位移、孔隙水压力和土介质总应力分量的积分形式的一般解;借助一般解,建立了有限厚度饱和土层和饱和半空间的精确动力刚度矩阵,并由土层的层间界面连续条件建立三维非轴对称层状饱和地基的总刚度方程;在此基础上,系统研究了横观各向同性饱和半空间体在内部集中荷载激励下的动力响应,并给出了问题的瞬态解答.该研究为运用边界元法求解饱和地基动力响应奠定了理论基础. 相似文献
16.
Andrei M. Ionov 《Geophysical Prospecting》2013,61(1):104-119
The dynamic response of a semi‐infinite fluid‐filled borehole embedded in an elastic half‐space under a concentrated normal surface load is analysed in the long‐wavelength limit. The solution of the problem is obtained with integral transforms in the form of a double integral with respect to the slowness and frequency. The partial P‐ and SV‐wave responses are further transformed to path integrals along Cagniard paths in the complex slowness plane. Unlike the traditional Cagniard‐de Hoop technique based on the Laplace transform of time dependence, this paper is based on the Fourier transform. The tube‐wave response is presented as a causal integral over a slowness range. The resultant representation in the time‐domain is suitable for the numerical evaluation of the complete response in the fluid‐filled borehole, especially at large distances. Asymptotic analysis of seismic phases arising in the borehole is performed on the basis of the obtained solution. The complete asymptotic wavefield consists in P‐ and SV‐waves, the Rayleigh wave and the low‐frequency Stoneley (tube) wave. Pressure synthetics obtained by the use of the asymptotic formulas are shown to be in good agreement with straightforward calculations. 相似文献
17.
Dynamic analysis of slab track on multi-layered transversely isotropic saturated soils subjected to train loads 总被引：1，自引：0，他引：1
The dynamic responses of a slab track on transversely isotropic saturated soils subjected to moving train loads are investigated by a semi-analytical approach. The track model is described as an upper Euler beam to simulate the rails and a lower Euler beam to model the slab. Rail pads between the rails and slab are represented by a continuous layer of springs and dashpots. A series of point loads are formulated to describe the moving train loads. The governing equations of track-ground systems are solved using the double Fourier transform, and the dynamic responses in the time domain are obtained by the inverse Fourier transform. The results show that a train load with high velocity will generate a larger response in transversely isotropic saturated soil than the lower velocity load, and special attention should be paid on the pore pressure in the vicinity of the ground surface. The anisotropic parameters of a surface soil layer will have greater influence on the displacement and excess pore water pressure than those of the subsoil layer. The traditional design method taking ground soil as homogeneous isotropic soil is unsafe for the case of RE 1 and RG 1, so a transversely isotropic foundation model is of great significance to the design for high train velocities. 相似文献
18.
W. G. CLEMENT 《Geophysical Prospecting》1973,21(1):125-145
The geophysicist involved in the analysis of two-dimensional data should have an understanding of the two-dimensional finite Fourier transform and the mechanics of two-dimensional filtering. Frequency aliasing must be considered when working with sampled data. In two dimensions it is advantageous to consider aliasing in terms of the overlap of the repeating spectra inherent in the finite Fourier transform. Two-dimensional filtering can be performed as a transient convolution in the space domain, as cyclic convolution utilizing the frequency domain or as the multiplication of polynomials using the z-transform. If the “edge” effects are removed, the results of the three methods are identical. 相似文献
19.
借助Walsh变换实现引力位球谐函数的快速Fourier变换导出了球谐函数的Walsh-Fourier变换、转换矩阵的快速Walsh-Hadamard变换算法及其数据压缩方法还讨论了Walsh-Fouriede换的特性及其在球谐分析中的应用研究表明：当序率和频率等同时．Walsh．Fourier变换和Fourier变换的结果完全一致,两者曲线形态相同;按双精度运算,两种方法的计算准确度均可达到±（10－15-10－14）;Walsh－Fourler变换可以用实数变换取代Fourier变换的复数变换;快速Walsh－Hadamard变换速度提高的幅度将随着阶数的增加而递增：Walsh-Fourier变换可以用于序率和频率等同或不等同的情形Walsh-Fourler变换可在计算精度、数据压缩和位场谱表示方面好于Fourier变换 相似文献
20.
S. J. William Stauder 《Pure and Applied Geophysics》1959,44(1):135-143
Summary The fault plane solution for the Alaska earthquake of October 24, 1927, is obtained and is found to be similar to that of the more recent earthquake of July 10, 1958, which occurred in the same general region. Using the data from the plane of polarization of theS waves, the mechanism at the focus is seen to conform to a single couple or faulting model rather than to a double couple. TheS waves, further, select one of the nodal planes forP as the fault plane. 相似文献