共查询到20条相似文献，搜索用时 31 毫秒
1.
2.
3.
《Advances in water resources》2002,25(2):159-177
Time integration methods that adapt in both the order of approximation and time step have been shown to provide efficient solutions to Richards' equation. In this work, we extend the same method of lines approach to solve a set of two-phase flow formulations and address some mass conservation issues from the previous work. We analyze these formulations and the nonlinear systems that result from applying the integration methods, placing particular emphasis on their index, range of applicability, and mass conservation characteristics. We conduct numerical experiments to study the behavior of the numerical models for three test problems. We demonstrate that higher order integration in time is more efficient than standard low-order methods for a variety of practical grids and integration tolerances, that the adaptive scheme successfully varies the step size in response to changing conditions, and that mass balance can be maintained efficiently using variable-order integration and an appropriately chosen numerical model formulation. 相似文献
4.
5.
层土动力Green函数数值计算中若干问题的探讨 总被引：2，自引：0，他引：2
动力Ｇｒｅｅｎ子数数值计算技术是动力边界元法得以实施的关键技术之一。本文 无对两类层状介质动力Ｇｒｅｅｎ函数的频域解进行了比较,在此基础上,文中着重探讨了在土－结构相互作用问题中应用显式动力Ｇｒｅｅｎ函数时的若干数值计算问题,包括土层的子层离散和波数的计算以及半空间的模拟等。 相似文献
6.
The numerical evaluation of certain integral transforms is required for the interpretation of some geophysical exploration data. Digital linear filter operators are widely used for carrying out such numerical integration. It is known that the method of Wiener–Hopf minimization of the error can be used to design very efficient, short digital linear filter operators for this purpose. We have found that, with appropriate modifications, this method can also be used to design longer filters. Two filters for the Hankel J0 transform (61-point and 120-point operators), and two for the Hankel J1 transform (47-point and 140-point operators) have been designed. For these transforms, the new filters give much lower errors compared to all other known filters of comparable, or somewhat longer, size. The new filter operators and some results of comparative performance tests with known integral transforms are presented. These filters would find widespread application in many numerical evaluation problems in geophysics. 相似文献
7.
《Advances in water resources》2007,30(4):1027-1045
Streamline methods have shown to be effective for reservoir simulation. For a regular grid, it is common to use the semi-analytical Pollock’s method to obtain streamlines and time-of-flight coordinates (TOF). The usual way of handling irregular grids is by trilinear transformation of each grid cell to a unit cube together with a linear flux interpolation scaled by the Jacobian. The flux interpolation allows for fast integration of streamlines, but is inaccurate even for uniform flow. To improve the tracing accuracy, we introduce a new interpolation method, which we call corner-velocity interpolation. Instead of interpolating the velocity field based on discrete fluxes at cell edges, the new method interpolates directly from reconstructed point velocities given at the corner points in the grid. This allows for reproduction of uniform flow, and eliminates the influence of cell geometries on the velocity field. Using several numerical examples, we demonstrate that the new method is more accurate than the standard tracing methods. 相似文献
8.
Iterative implicit integration procedure for hybrid simulation of large nonlinear structures 总被引：1，自引：0，他引：1
A fully implicit iterative integration procedure is presented for local and geographically distributed hybrid simulation of the seismic response of complex structural systems with distributed nonlinear behavior. The purpose of this procedure is to seamlessly incorporate experimental elements in simulations using existing fully implicit integration algorithms designed for pure numerical simulations. The difficulties of implementing implicit integrators in a hybrid simulation are addressed at the element level by introducing a safe iteration strategy and using an efficient procedure for online estimation of the experimental tangent stiffness matrix. In order to avoid physical application of iterative displacements, the required experimental restoring force at each iteration is estimated from polynomial curve fitting of recent experimental measurements. The experimental tangent stiffness matrix is estimated by using readily available experimental measurements and by a classical diagonalization approach that reduces the number of unknowns in the matrix. Numerical and hybrid simulations are used to demonstrate that the proposed procedure provides an efficient method for implementation of fully implicit numerical integration in hybrid simulations of complex nonlinear structures. The hybrid simulations presented include distributed nonlinear behavior in both the numerical and experimental substructures. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
9.
The existing on‐line numerical integration algorithms are derived from the Newmark method, which is based on an approximation of derivatives in the differential equation. The state–space procedure (SSP), based on an interpolation of the discrete excitation signals for piecewise convolution integral, has been confirmed as more reliable than the Newmark method in terms of numerical accuracy and stability. In an attempt to enhance the pseudodynamic test, this study presents an on‐line integration algorithm (referred to as the OS–SSP method) via an integration of the state–space procedure with Nakashima's operator‐splitting concept. Numerical stability and accuracy assessment of the proposed algorithm in addition to the explicit Newmark method and the OS method were investigated via an eigenvalue, frequency‐domain and time‐domain analysis. Of the on‐line integration algorithms investigated, the OS–SSP method is demonstrated as the most accurate method with an acceptable stability (although not unconditionally stable) characteristic. Therefore, the OS–SSP method is the most desirable method for pseudodynamic testing if the numerical stability criterion (Δt/T⩽0.5) is ensured for every vibration mode involved. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
10.
Developing robust and efficient numerical solution methods for Richards' equation (RE) continues to be a challenge for certain problems. We consider such a problem here: infiltration into unsaturated porous media initially at static conditions for uniform and non-uniform pore size media. For ponded boundary conditions, a sharp infiltration front results, which propagates through the media. We evaluate the resultant solution method for robustness and efficiency using combinations of variable transformation and adaptive time-stepping methods. Transformation methods introduce a change of variable that results in a smoother solution, which is more amenable to efficient numerical solution. We use adaptive time-stepping methods to adjust the time-step size, and in some cases the order of the solution method, to meet a constraint on nonlinear solution convergence properties or a solution error criterion. Results for three test problems showed that adaptive time-stepping methods provided robust solutions; in most cases transforming the dependent variable led to more efficient solutions than untransformed approaches, especially as the pore-size uniformity increased; and the higher-order adaptive time integration method was robust and the most efficient method evaluated. 相似文献
11.
A stream tube integration method is introduced to solve transient subsurface fluid flow problems. The method combines a geometry-embedded form of Darcy's Law and the notion of location of average. Two types of problems, transient radial flow to a well of finite radius in an areally infinite aquifer and in a double porosity system, are solved by the stream tube integration method and the integral finite difference method. Results of the solutions show that the stream tube integration method, with fixed coarse mesh, are more accurate and better behaved than the integral finite difference method, with fine mesh. The fixed mesh stream tube integration method is readily extended to the moving mesh method. With much coarse mesh, the moving mesh technique can obtain the same accurate results as the fixed mesh stream tube integration method. It is suggested that the stream tube integration method is a viable way to state, solve, interpret and verify numerical solutions. The method provides efficient computation and improved accuracy for analysing subsurface fluid flow. © 1998 John Wiley & Sons, Ltd. 相似文献
12.
A semi-analytical approach for analyzing ground vibrations caused by trains moving over elevated bridges 总被引：3，自引：0，他引：3
A semi-analytic approach is presented for the three-dimensional analysis of ground vibrations induced by trains moving over elevated bridges. The train is modeled as two sets of moving loads, with one for the front wheels and the other for the rear ones, the elevated bridge as a series of elastically supported beams, and the ground as a viscoelastic half space. Three key elements are considered in the solution: (1) the analytic solution for the vibration of an elastically supported beam under a series of moving loads, (2) the impedance of the foundation–soil system, and (3) Green's function for an elastic half space under a harmonic point load. Such an approach allows us to consider the structural dynamics of the elevated bridge, the foundation–soil interaction, and the wave propagation characteristics in the half space. From the numerical examples studied, the proposed approach was demonstrated to be accurate and efficient. The framework of analysis described herein can be generalized to solve problems with complex foundations and layered soils. 相似文献
13.
Implementation and application of the unconditionally stable explicit parametrically dissipative KR‐α method for real‐time hybrid simulation 
下载免费PDF全文
下载免费PDF全文 Chinmoy Kolay James M. Ricles Thomas M. Marullo Akbar Mahvashmohammadi Richard Sause 《地震工程与结构动力学》2015,44(5):735-755
In real‐time hybrid simulations (RTHS) that utilize explicit integration algorithms, the inherent damping in the analytical substructure is generally defined using mass and initial stiffness proportional damping. This type of damping model is known to produce inaccurate results when the structure undergoes significant inelastic deformations. To alleviate the problem, a form of a nonproportional damping model often used in numerical simulations involving implicit integration algorithms can be considered. This type of damping model, however, when used with explicit integration algorithms can require a small time step to achieve the desired accuracy in an RTHS involving a structure with a large number of degrees of freedom. Restrictions on the minimum time step exist in an RTHS that are associated with the computational demand. Integrating the equations of motion for an RTHS with too large of a time step can result in spurious high‐frequency oscillations in the member forces for elements of the structural model that undergo inelastic deformations. The problem is circumvented by introducing the parametrically controllable numerical energy dissipation available in the recently developed unconditionally stable explicit KR‐α method. This paper reviews the formulation of the KR‐α method and presents an efficient implementation for RTHS. Using the method, RTHS of a three‐story 0.6‐scale prototype steel building with nonlinear elastomeric dampers are conducted with a ground motion scaled to the design basis and maximum considered earthquake hazard levels. The results show that controllable numerical energy dissipation can significantly eliminate spurious participation of higher modes and produce exceptional RTHS results. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
14.
15.
各向异性介质模型电性结构复杂,如何进行合理的网格剖分成为获得高精度正演结果的关键,为此本文开展时间域航空电磁各向异性大地三维自适应有限元正演算法研究.通过结合非结构时间域有限元算法和自适应网格优化技术,实现各向异性介质条件下三维时间域航空电磁自适应正演.考虑到时间域航空电磁响应随时间的衰减特性,为了综合评价不同时刻的后验误差,本文将时间作为加权因子,调整各个时刻后验误差的相对权重,进而实现对浅部和深部网格的同步优化.通过与一维解析结果进行对比验证了本文算法的可靠性.数值实验结果显示电导率各向异性对自适应网格影响严重,其最大主轴电导率的数值及其分布特征直接决定了网格加密效果.此外,各向异性对时间域航空电磁三分量响应的分布形态和异常幅值也会产生严重影响,利用全域视电阻率极性图,可以很好地识别各向异性主轴方向. 相似文献
16.
E. Bertrand A.-M. Duval J. R��gnier R. M. Azzara F. Bergamashi P. Bordoni F. Cara G. Cultrera G. Di Giulio G. Milana J. Salichon 《Bulletin of Earthquake Engineering》2011,9(3):809-823
During the microzonation studies of the April 6th, 2009 L??Aquila earthquake, we observed local seismic amplifications in the Roio area??a plane separated from L??Aquila city center by mount Luco. Six portable, digital instruments were deployed across the plain from 15 April to mid-May 2009. This array recorded 152 aftershocks. We analyzed the ground motion from these events to determine relative site amplification within the plain and on surrounding ridges. Horizontal over vertical spectral ratio on noise data (HVSRN), aftershock recordings (HVEQ) and standard spectral ratio (SSR) showed amplifications at 1.3 and 4.0?Hz on quaternary deposits. Seismic amplifications in the frequency range of 4 and 6?Hz were also observed on a carbonate ridge of Colle di Roio, on the northwestern border of the plateau. A small amplification was noticed near the top of mount Luco, another rocky site. Large discrepancies in the amplification levels between methods have been observed for these sites, but the HVSRN, HVEQ and SSR gave similar results at the stations located in the Roio plain. On the rocky sites, the SSR was more reliable than the HVSRN at estimating the transfer function of the site, even if the resonance frequency seemed to be well detected by the latter method. 相似文献
17.
《Advances in water resources》1998,22(2):185-195
We present an efficient numerical method for solving stochastic porous media flow problems. Single-phase flow with a random conductivity field is considered in a standard first-order perturbation expansion framework. The numerical scheme, based on finite element techniques, is computationally more efficient than traditional approaches because one can work with a much coarser finite element mesh. This is achieved by avoiding the common finite element representation of the conductivity field. Computations with the random conductivity field only arise in integrals of the log conductivity covariance function. The method is demonstrated in several two- and three-dimensional flow situations and compared to analytical solutions and Monte Carlo simulations. Provided that the integrals involving the covariance of the log conductivity are computed by higher-order Gaussian quadrature rules, excellent results can be obtained with characteristic element sizes equal to about five correlation lengths of the log conductivity field. Investigations of the validity of the proposed first-order method are performed by comparing nonlinear Monte Carlo results with linear solutions. In box-shaped domains the log conductivity standard deviation σY may be as large as 1.5, while the head variance is considerably influenced by nonlinear effects as σY approaches unity in more general domains. 相似文献
18.
Minxian Cheng 《地震工程与结构动力学》1988,16(2):227-244
The convergence and stability of step-by-step integration schemes used in the inelastic dynamic analysis of structures and their corresponding criteria were studied for a restoring force model with negative-stiffness. Convergence conditions and stability conditions 1, 2 or 3 were established. The numerical stability of the integration under negative-stiffness belongs to the category of relative stability; consequently, the concepts and the conclusions concerning numerical stability in the case of positive-stiffness (which belongs to absolute stability) cannot be used. Research into several step-by-step integration methods usually employed in inelastic dynamic analysis has shown great differences in numerical stability for models with negative-stiffness as compared with positive-stiffness models. The central difference method is convergent and unconditionally stable in the case of negative-stiffness, though it is only conditionally stable in the case of positive-stiffness. The Houbolt method satisfies the requirement for convergence; its stability, however, depends not only on the integration step size Δt but also on the stiffness ratio β for the model with negative-stiffness, unlike the unconditional stability for the model with positive-stiffness. The Newmark constant acceleration method is convergent and unconditionally stable in the case of negative-stiffness just like it is in the case of positive-stiffness. 相似文献
19.
本文将应用2.5维数值模式匹配算法研究建立柱状横向同性地层中偏心条件下磁流源并矢Green函数的高效算法,并利用该算法高效计算多分量阵列感应井眼校正库.首先,通过Fourier展开技术将偏心条件下磁流源并矢Green函数的数值模拟转化为一系列的轴对称问题;然后利用模式匹配算法求解轴对称问题,得到柱状介质中偏心条件下磁流源并矢Green函数的半解析解;最后,结合中国石油集团测井有限公司开发的三维感应测井仪器(TDIT)的结构参数,推导出复合线圈系视电导率张量的计算公式.在此基础上,根据井场实际地层资料归纳出理论地层模型相应参数的变化范围,高效高精度建立一套直接应用于工业生产的多分量阵列感应井眼校正库,并利用井眼校正库详细地分析考察各模型参数对TDIT响应的影响. 相似文献
20.
重力勘探中复杂条件下的三维正演计算量大存储要求高,使得这种条件下重力勘探高效、精细正反演变得困难.针对这一问题,提出一种空间-波数混合域数值模拟方法,该方法将空间域引力位积分进行水平方向二维傅里叶变换,将三维空间域卷积问题转换为多个不同波数之间相互独立的空间垂向一维积分问题,一维积分垂向可离散为多个单元积分之和,每个单元采用二次形函数表征密度变化,可得出单元积分的解析表达式.该方法计算量和存储需求少,算法高度并行;保留垂向为空间域,优势之一在于可根据实际情况合理调整单元疏密程度,准确模拟任意复杂地形和密度异常体的重力异常,兼顾计算精度与计算效率;优势之二在于用形函数拟合求得积分的解析解,计算精度和效率高;充分利用一维形函数积分的高效和高精度,不同波数之间一维积分高度并行性及快速傅里叶变换的高效性,实现重力异常场三维数值模拟.设计棱柱体模型,通过数值解和解析解对比验证了该方法的正确性、适用性和高效性.针对任意复杂地形条件下的重力场及其张量的模拟问题,提出一种快速算法,对其有效性进行了验证.探究标准FFT法的截断效应对计算精度的影响,对比分析Gauss-FFT法和标准FFT扩边法两种方法的计算精度和效率,总结了二者的选取策略,结果表明选用标准FFT扩边法计算效率更高.实际地形的数值模拟表明本文算法适用于任意复杂地形的高效计算. 相似文献