共查询到20条相似文献，搜索用时 31 毫秒
1.
在地震动数值模拟方法中,谱元法和有限元法是应用较广泛的两种方法。基于经典的Lamb问题模型,首先推导给出地表竖向位移的解析解答。然后分别利用常用的四阶谱元法和线性有限元法,模拟了地表脉冲力源作用下模型的位移响应。考虑有意义的最短波长内的采样点个数及单元高宽比的变化,对比了两种方法的模拟精度;结果表明:对于谱元法,观测点与波源之间需至少包含两个网格,在此条件下,最短波长内包含一个网格（最短波长内5个采样点）时,数值解与解析解的误差小于1%,已达很高的精度;对于有限元法,最短波长内需包含10个网格时才能达到这一精度。此外,在满足网格尺寸要求的前提下,单元水平向与垂直向尺寸的比值在1∶1到5∶1的范围内时,谱元法和有限元法的模拟精度均变化不大。因此,单位波长内采样点个数相同时,谱元法的模拟精度比有限元法高的多,同时,在一定范围内两种方法的模拟结果对于宽高比的变化不敏感。 相似文献
2.
《地震研究进展（英文）》2022,2(4):100118
The pushover method for underground structures is a seismic analysis method featured by high calculation accuracy and a simple implementation process. The method has been widely used in seismic design and other related scientific research; however, the influence of different soil-structure flexibility ratios on the accuracy of this method is still not well understood. In this study, we select the cross-section structures beneath the Daikai subway station as the research object and establish 12 finite element analysis models with different soil-structure flexibility ratios using ABAQUS. All models are computed by the dynamic time-history method or the pushover method. Furthermore, the dynamic time-history solution result is taken as the standard solution, and the precision and application of the pushover analysis method are discussed based on the parameters of peak interlayer displacement and peak internal force of the middle column section. The results show that the soil-structure flexibility ratio has a significant influence on the calculation accuracy of the pushover method, and the calculation accuracy of this method is the most ideal when the soil-structure flexibility is equal to 1. The research results can provide significant references for the seismic design of underground structures or the improvement of simplified seismic analysis methods. 相似文献
3.
The dynamic element method has been shown previously to provide a computational advantage over the ordinary finite element method for various beam elements. The Taylor expansions are computed here for the dynamic shape functions (two terms) and dynamic stiffness matrix (four terms) for the axisymmetric vibrations of an annular plate element. The complicated matrices which result are made more tractable by expressing them as power series in powers of the aspect ratio. The percentage error in the natural frequencies is then calculated using both the two- and the three-term dynamic stiffness matrix, demonstrating the increased accuracy for a given number of elements. 相似文献
4.
盆地效应作为场地效应的一种特殊类型,其对地震动的影响仍在持续研究当中。建立二维梯形沉积盆地模型,模型边界采用透射人工边界,利用显式有限元法作为求解方法。根据Buckingham提出的π定理定义了多个影响盆地效应的无量纲参数,包括盆地深宽比、介质波速比、盆地深度与入射波长之比、盆地开口宽度与入射波长之比。以Ricker子波作为输入波,研究SH波垂直入射盆地模型时,以上无量纲参数对盆地放大效应的影响。研究发现盆地深宽比越大盆地诱发的聚焦效应越明显,盆地深宽比越小盆地边缘效应越明显。在盆地的深宽比和入射波确定时,随着盆地开口宽度与入射波长之比不断增大盆地的最值放大系数逐渐从盆地中心移向盆地边缘,表现出越来越明显的边缘放大效应。其中波速比差异对放大系数的影响最大,波速比差异越大介质的波阻抗比越大,放大效应越显著。通过将影响盆地放大效应的参数无量纲化,使梯形盆地对地震动放大效应的研究的结论更具有普遍意义。 相似文献
5.
数字岩心微观孔隙结构十分复杂,有限元模拟物性参数与弹性参数之间关系是非线性的,直接反演其物性参数准确度低、稳定性差.本文发展了一种数字岩石物理逆建模方法,实现了基于数字岩心的储层参数有效预测.从数字岩心基函数的构建出发,基于有限元方法,计算了一系列具有等间距物性参数值(孔隙度、泥质含量和含水饱和度)的数字岩心弹性参数(体积模量、剪切模量和密度),通过插值算法建立了数字岩心弹性参数三维数据集,从而实现了弹性模量的有限元数值解的快速构建;然后搜索弹性参数的单值等值面,通过等值面的空间交会得到交点,完成储层参数预测.测试结果表明:基于数字岩心逆建模理论的储层参数预测结果与实际模型一致,具有可行性,并且可以通过增加插值点数目提高预测的准确性;孔隙度和泥质含量预测结果稳定性很好,而含水饱和度对噪声的加入较为敏感. 相似文献
6.
The scaled boundary finite‐element method has been developed for the dynamic analysis of unbounded domains. In this method only the boundary is discretized resulting in a reduction of the spatial dimension by one. Like the finite‐element method no fundamental solution is required. This paper extends the scaled boundary finite‐element method to simulate the transient response of non‐homogeneous unbounded domains with the elasticity modulus and mass density varying as power functions of spatial coordinates. To reduce the number of degrees of freedom and the computational cost, the technique of reduced set of base functions is applied. The scaled boundary finite‐element equation for an unbounded domain is reformulated in generalized coordinates. The resulting acceleration unit‐impulse response matrix is obtained and assembled with the equation of motion of standard finite elements. Numerical examples of non‐homogeneous isotropic and transversely isotropic unbounded domains demonstrate the accuracy of the scaled boundary finite‐element method. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
7.
对岷县地震台、庆阳黄土塬2~12m深度土层的粉土、粉质黏土样品,利用室内动三轴试验和扫描电镜试验,测定阻尼比和微观结构特征,研究了粉土、粉质黏土的阻尼比随深度的变化及应力应变关系,采用灰色关联度方法,分析了黄土竖向（沉积方向）和横向（垂直于沉积方向）的微观结构参数对阻尼比影响的主次关系及其权重。研究结果表明:小应变下,粉土和粉质黏土阻尼比随着深度的增加而减小;粉土和粉质黏土的微观参数与阻尼比有一定的关联性,粉土竖向的微观参数对阻尼比的影响较为敏感;粉质黏土横向的微观参数对阻尼比的影响更加敏感;影响权重方面,粉质黏土的平均形状系数影响权重最大,对阻尼比的影响最大;粉土的平均形状系数、分形维数和概率熵影响权重最大,对阻尼比影响最敏感。平均面积无论在关联度还是影响权重方面对阻尼比影响最小,该结论粉土与粉质黏土保持一致。 相似文献
8.
本文推导了声波方程频散函数, 分析了伪谱法的空间网格大小和采样周期对数值频散的影响, 通过数值模拟实验得到了最佳空间参数选择方法。 结果表明: 伪谱方法稳定数值模拟的最大空间采样间距选取原则是使中波长(奈奎斯特频率的一半)的采样点数为2个; 对于所有维度, 稳定性随空间采样间距的增加而增加, 但不易变化太大, 变化太大时需要适当减小震源子波的主频, 以满足空间合理采样; 空间采样间距的大小设置, 需要考虑满足采样定理和稳定性计算条件, 并且稳定性条件对空间采样间距的要求更加严格; 伪谱法数值模拟的最佳(数值频散最小)空间参数选择为中波长2个采样点, 对应主波长约6~7个采样点。 以上研究对于采用伪谱法进行声波方程数值模拟过程中, 如何合理选择模拟参数提供一些参考。 相似文献
9.
M.Th. Van Genuchten 《Advances in water resources》1983,6(2):106-111
This paper describes a Galerkin-type finite element solution of the two-dimensional saturated-unsaturated flow equation. The numerical solution uses an incomplete (reduced) set of Hermitian cubic basis functions and is formulated in terms of normal and tangential coordinates. The formulation leads to continuous pressure gradients across interelement boundaries for a number of well-defined element configurations, such as for rectangular and circular elements. Other elements generally lead to discontinuous gradients; however, the gradients remain uniquely defined at the nodes. The method avoids calculation of second-order derivatives, yet retains many of the advantages associated with Hermitian elements. A nine-point Lobatto-type integration scheme is used to evaluate all local element integrals. This alternative scheme produces about the same accuracy as the usual 9- or 16-point Gaussian quadrature schemes, but is computationally more efficient. 相似文献
10.
Local and global shape functions in a boundary element formulation for the calculation of traffic induced vibrations 总被引：1，自引：1，他引：1
This paper compares the use of local and global shape functions in a boundary element method that is used in a prediction model for traffic induced vibrations. The boundary element formulation describes the interaction problem between a linear elastic layered half-space and a longitudinally invariant structure representing a road or a railway track. The boundary element formulation in the frequency–wavenumber domain is obtained by means of a weighted residual method. Constant element shape functions, as well as Legendre and Chebyshev shape functions are considered. Their effect on both accuracy and computational effort is investigated. The presence of a singularity in the Chebyshev based shape functions allows to obtain a better approximation for the soil tractions. The theory is applied to road traffic induced vibrations where the response is calculated in a large number of output points. 相似文献
11.
This paper presents a step-by-step procedure using the three-dimensional boundary element approach to study the behavior of semi-circular canyons under seismic shear waves. The boundary element code TDASC allows utilization for various canyon geometries, evaluation of concurrent seismic waves and calculation of the ground motions on canyons due to an excitation at any arbitrary point of the incident field. Considering the widening ratio of the canyon(including prismatic, semi-prismatic and non-prismatic canyons), wave characteristics(wavelength, dimensionless period, direction) and maximum amplification pattern, the solution was applied to carry out a series of parametric studies. It was shown that canyon form can significantly affect the displacement amplification, especially at the points located on its edges. By increasing the wave dimensionless frequency(η 1), the amplification pattern becomes more complex. On the basis of the results from a variety of considered cases, a new expression has been presented for the limiting wavelength beyond which the widening of the canyon will not have a major effect on the displacement amplification. To verify the reliability of the proposed approach, the obtained results, expressed in terms of displacement amplitude, were compared with those from the available published literature and a reasonably good agreement was observed. 相似文献
12.
《Advances in water resources》2004,27(2):129-140
The discontinuous spectral Galerkin method uses a finite-element discretization of the groundwater flow domain with basis functions of arbitrary order in each element. The independent choice of the basis functions in each element permits discontinuities in transmissivity in the flow domain. This formulation is shown to be of high order accuracy and particularly suitable for accurately calculating the flow field in porous media. Simulations are presented in terms of streamlines in a bidimensional aquifer, and compared with the solution calculated with a standard finite-element method and a mixed finite-element method. Numerical simulations show that the discontinuous spectral Galerkin approximation is more efficient than the standard finite-element method (in computing fluxes and streamlines/pathlines) for a given accuracy, and it is more accurate on a given grid. On the other hand the mixed finite-element method ensures the continuity of the fluxes at the cell boundaries and it is particular efficient in representing complicated flow fields with few mesh points. Simulations show that the mixed finite-element method is superior to the discontinuous spectral Galerkin method producing accurate streamlines even if few computational nodes are used. The application of the discontinuous Galerkin method is thus of interest in groundwater problems only when high order and extremely accurate solutions are needed. 相似文献
13.
A direct time domain boundary element method is presented based on the Stokes fundamental solutions, discretized in both time and space, and an efficient time step-by-step solution that minimizes the accumulation of errors. A non-singular numerical integration procedure, in the Cauchy sense, is proposed for the generation of the associated influence matrices. This methodology is shown to be efficient for the solution of a number of computationally intensive problems in the area of soil–structure interaction. In addition, an algorithm for the direct calculation of the response of massive foundations to externally applied forces and/or obliquely incident seismic waves is introduced. The accuracy and computational efficiency of the proposed methodologies is established through a number of comparison studies. 相似文献
14.
ＪｏｉｎｔｓｏｌｕｔｉｏｎｏｆｌｅｖｅｌｉｎｇａｎｄｇｒａｖｉｔｙｄａｔａｉｎｃｒｕｓｔａｌｄｅｆｏｒｍａｔｉｏｎｍｏｎｉｔｏｒｉｎｇＺＵＳＨＥＮＧＺＨＡＮＧ１）（张祖胜）ＹＵＡＮＸＩＹＡＮＧ２）（杨元喜）ＨＡＮＲＯＮＧＳＵＮ１）（孙汉荣）Ｚ... 相似文献
15.
The analytical solution and algorithm for simulating the electric potential in an arbitrarily anisotropic multilayered medium produced by a point DC source is here proposed. The solution is presented as a combination of Hankel transforms of integer order and Fourier transforms based on the analytical recurrent equations obtained for the potential spectrum. For the conversion of the potential spectrum into the space domain, we have applied the algorithm of the Fast Fourier Transform for logarithmically spaced points. A comparison of the modelling results with the power‐series solution for two‐layered anisotropic structures demonstrated the high accuracy and computing‐time efficiency of the method proposed. The results of the apparent‐resistivity calculation for both traditional pole‐pole and tensor arrays above three‐layered sequence with an azimuthally anisotropic second layer are presented. The numerical simulations show that both arrays have the same sensitivity to the anisotropy parameters. This sensitivity depends significantly on the resistivity ratio between anisotropic and adjacent layers and increases for the models with a conductive second layer. 相似文献
16.
The problem of resolving or parameterising small-scale processes in oceanographic models and the extent to which small-scale
effects influence the large scale are briefly discussed and illustrated for a number of cases. For tides and surges in near-shore
regions, the advantages of using a graded mesh to resolve coastal and estuarine small-scale features are demonstrated in terms
of a west coast of Britain unstructured mesh model. The effect of mesh resolution upon the accuracy of the overall solution
is illustrated in terms of a finite element model of the Irish Sea and Mersey estuary. For baroclinic motion at high Froude
number, the effect of resolving small-scale topography within a non-hydrostatic model is illustrated in terms of tidally induced
mixing at a single sill, or two closely spaced sills. The question of how to parameterise small-scale non-linear interaction
processes that lead to significant mixing, in a form suitable for coarser grid hydrostatic models, is briefly considered.
In addition, the importance of topographically induced mixing that occurs in the oceanic lateral boundary layer, namely, the
shelf edge upon the large-scale ocean circulation is discussed together with the implications for coarse grid oceanic climate
models. The use of unstructured grids in these models to enhance resolution in shelf-edge regions in a similar manner to that
used in storm surge models to enhance near coastal resolution is suggested as a suitable “way forward” in large-scale ocean
circulation modelling. 相似文献
17.
This paper presents the concept of using an additional generator to prevent ground vibrations. A linear, transversally isotropic three dimensional half-space with the hysteretic damping model, acted upon by a harmonic vertical excitation is assumed. Equations of motion for the transversally isotropic ground model with the absorbing boundary conditions are presented and numerically integrated using FlexPDE software, based on the finite element method. The efficiency of the solution is analysed in terms of reducing the vertical and horizontal components of ground surface vibrations. Results in the form of a dimensionless amplitude reduction factor are presented for four different locations of a generator. The influence of the soil parameters and layers locations on the additional generator's efficiency is investigated. The vibration reduction efficiency in a four-story building is also presented. 相似文献
18.
19.
A novel finite element method has been described in this paper for the solution of the steady-state scalar wave equation in two dimensions. In this method the physical domain of the problem is mapped into an image domain; the governing equation and the prescribed boundary conditions are also appropriately transformed. Because logarithmic mapping functions are used, the physical domain is logarithmically condensed in the image plane. The method is therefore particularly suitable for the efficient and economical solution of large or very large aspect ratio problems. The high degree of accuracy which the method is capable of is demonstrated by means of two typical examples. Possible extension of the method to complex problems including non-linearity, multiplicity, etc. is also discussed. 相似文献
20.
B. Nath 《地震工程与结构动力学》1981,9(1):33-51
A novel method has been proposed for the finite element solution of the steady-state scalar wave equation in three-dimensions. In this the governing equation and the prescribed boundary conditions in the physical space are transformed into a spherical polar space in which the radial direction is logarithmically condensed; the physical problem domain is also mapped into the new space. The transformed equation is then solved in the mapped domain using conventional finite elements. Because physical dimensions of the problem are logarithmically condensed in the proposed spherical polar space, the method is particularly suitable for solving truly three-dimensional problems in which the aspect ratio(s) is large or very large. A number of illustrative examples considered show that the proposed method is capable of a high degree of accuracy, achieved efficiently and economically. A hybrid scheme has also been proposed for dealing with awkward-shaped domains. 相似文献