共查询到20条相似文献，搜索用时 436 毫秒
1.
J. MullikenD.C. Rizos 《Soil Dynamics and Earthquake Engineering》2012,34(1):78-88
This work presents an efficient and stable methodology for the coupling of Finite Element Methods (FEM) and Boundary Element Methods (BEM) that is independent of the particular solver and allows for independent temporal discretizations among solvers. The approach satisfies explicitly compatibility conditions and equilibrium of forces at the contact interfaces. Although the proposed approach has been developed in view of the soil-rail-vehicle dynamic interaction problem in High Speed Rail applications, it is expressed in a general form applicable to any multi-domain, multi-phase transient problem. The method development and formulations are presented in detail. Verification and application studies demonstrate the accuracy, efficiency and versatility of the method for the direct time domain solution of dynamic problems including structure-structure interaction and soil-structure interaction. The proposed approach demonstrates high accuracy and efficiency to that of direct coupling solutions and more rigorous methods. 相似文献
2.
A study on the dynamic response of a railway track is presented via a 3-D formulation based on the frequency domain Boundary Element Method (BEM) and the Finite Element Method (FEM). The railway track consists of a group of surface, massive, rigid footings resting on a viscoelastic half-space and connected by an overlaying rail structure. The BEM, employing the full-space fundamental solutions and quadrilateral elements, is used for the simulation of the elastic half-space while the FEM is used to model the rigid footings and the rail superstructure. The loading function consists of a set of externally applied, harmonic or transient loads. Frequency as well as transient, by way of FFT, results are presented for various modes of vibration. Various numerical studies assess the through-the-soil interaction of the adjacent footings, the influence of soil damping, the effect of the overlaying structure on the frequency content of the system, and the effective simulation of an infinitely long railway track by a truncated one. 相似文献
3.
热演化历史的研究对于盆地分析和油气勘探具有重要意义.北黄海盆地是中国近海海域油气勘探程度较低的盆地之一,迄今未见专门的热演化研究报道.笔者利用北黄海盆地中生代砂岩的磷灰石裂变径迹分析结果,结合地质条件约束,模拟获得了中生代以来盆地的热演化史.结果表明,北黄海盆地经历了两次增温和两次冷却的热演化过程,并在100-80 Ma时盆地的热历史出现明显变化,表明在晚白垩世北黄海盆地发生过一次较大的构造-热事件.磷灰石裂变径迹分析所表明的北黄海盆地的热历史与盆地原型演化阶段相对应,而这种热历史和盆地的演化过程与区域构造背景相关.磷灰石裂变径迹所揭示的热演化史对于深化认识北黄海盆地的地质演化过程和油气勘探潜力具有重要意义. 相似文献
4.
This paper introduces a new method for simulating large-scale subsurface contaminant transport that combines an Analytic Element Method (AEM) groundwater flow solution with a split-operator Streamline Method for modeling reactive transport. The key feature of the method is the manner in which the vertically integrated AEM flow solution is used to construct three-dimensional particle tracks that define the geometry of the Streamline Method. The inherently parallel nature of the algorithm supports the development of reactive transport models for spatial domains much larger than current grid-based methods. The applicability of the new approach is verified for cases with negligible transverse dispersion through comparisons to analytic solutions and existing numerical solutions, and parallel performance is demonstrated through a realistic test problem based on the regional-scale transport of agricultural contaminants from spatially distributed sources. 相似文献
5.
We developed a new semi-analytical source zone depletion model (SZDM) for multicomponent light nonaqueous phase liquids (LNAPLs) and incorporated this into an existing screening model for estimating cleanup times for chemical spills from railroad tank cars that previously considered only single-component LNAPLs. Results from the SZDM compare favorably to those from a three-dimensional numerical model, and from another semi-analytical model that does not consider source zone depletion. The model was used to evaluate groundwater contamination and cleanup times for four complex mixtures of concern in the railroad industry. Among the petroleum hydrocarbon mixtures considered, the cleanup time of diesel fuel was much longer than E95, gasoline, and crude oil. This is mainly due to the high fraction of low solubility components in diesel fuel. The results demonstrate that the updated screening model with the newly developed SZDM is computationally efficient, and provides valuable comparisons of cleanup times that can be used in assessing the health and financial risk associated with chemical mixture spills from railroad-tank-car accidents. 相似文献
6.
A frequency‐dependent and intensity‐dependent macroelement for reduced order seismic analysis of soil‐structure interacting systems 
下载免费PDF全文
下载免费PDF全文 The computational demand of the soil‐structure interaction analysis for the design and assessment of structures, as well as for the evaluation of their life‐cycle cost and risk exposure, has led the civil engineering community to the development of a variety of methods toward the model order reduction of the coupled soil‐structure dynamic system in earthquake regions. Different approaches have been proposed in the past as computationally efficient alternatives to the conventional finite element model simulation of the complete soil‐structure domain, such as the nonlinear lumped spring, the macroelement method, and the substructure partition method. Yet no approach was capable of capturing simultaneously the frequency‐dependent dynamic properties along with the nonlinear behavior of the condensed segment of the overall soil‐structure system under strong earthquake ground motion, thus generating an imbalance between the modeling refinement achieved for the soil and the structure. To this end, a dual frequency‐dependent and intensity‐dependent expansion of the lumped parameter modeling method is proposed in the current paper, materialized through a multiobjective algorithm, capable of closely approximating the behavior of the nonlinear dynamic system of the condensed segment. This is essentially the extension of an established methodology, also developed by the authors, in the inelastic domain. The efficiency of the proposed methodology is validated for the case of a bridge foundation system, wherein the seismic response is comparatively assessed for both the proposed method and the detailed finite element model. The above expansion is deemed a computationally efficient and reliable method for simultaneously considering the frequency and amplitude dependence of soil‐foundation systems in the framework of nonlinear seismic analysis of soil‐structure interaction systems. 相似文献
7.
瞬变电磁法的三维有限元正演通常采用齐次边界条件,为满足该边界条件,需要构建较大尺寸的模型,这降低了正演问题的求解速度.针对该问题,本文采用吸收边界条件代替齐次边界条件,以缩小模型体积,加快正演速度:首先,从时间域麦克斯韦方程组出发,推导了基于库伦规范的矢量势的微分控制方程,结合一阶吸收边界条件推导了相应的的弱形式方程;在此基础上采用一阶四面体矢量单元进行单元分析、Newmark法进行时间离散,实现了瞬变电磁法的快速三维正演.通过均匀半空间模型的解析解,H型地电断面的CR1Dmod解和相应模型有限元解的对比,验证了本文算法的正确性.均匀半空间模型分别采用吸收边界条件和齐次边界条件的正演结果对比表明:吸收边界条件确实可以提高三维正演的精度或者缩小模型尺寸、加快计算速度. 相似文献
8.
We present advances in compositional modeling of two-phase multi-component flow through highly complex porous media. Higher-order methods are used to approximate both mass transport and the velocity and pressure fields. We employ the Mixed Hybrid Finite Element (MHFE) method to simultaneously solve, to the same order, the pressure equation and Darcy's law for the velocity. The species balance equation is approximated by the discontinuous Galerkin (DG) approach, combined with a slope limiter. In this work we present an improved DG scheme where phase splitting is analyzed at all element vertices in the two-phase regions, rather than only as element averages. This approximation is higher-order than the commonly employed finite volume method and earlier DG approximations. The method reduces numerical dispersion, allowing for an accurate capture of shock fronts and lower dependence on mesh quality and orientation. Further new features are the extension to unstructured grids and support for arbitrary permeability tensors (allowing for both scalar heterogeneity, and shear anisotropy). The most important advancement in this work is the self-consistent modeling of two-phase multi-component Fickian diffusion. We present several numerical examples to illustrate the powerful features of our combined MHFE–dg method with respect to lower-order calculations, ranging from simple two component fluids to more challenging real problems regarding CO2 injection into a vertical domain saturated with a multi-component petroleum fluid. 相似文献
9.
The use of ballastless slab track is significantly increasing in HST line technology. This development is due to some structural and operational advantages over ballasted track. In addition, floating slab tracks can be used to control ground-borne vibrations generated by surface and underground rail transportation systems. In this paper, a general and fully three dimensional multi-body-finite element-boundary element model is used to study vibrations due to train passage on ballast and non-ballast tracks. The vehicle is modelled as a multi-body system, the track, in both cases, using finite elements and the soil is represented using boundary elements. The three components of the load are considered; the quasi-static excitation (force generated by moving axle loads), the parametric excitation due to discrete supports of the rails and the excitation due to wheel and rail roughness and track unevenness. Track receptances are computed for both track systems and vibrations induced by high-speed train passage at the track and the free-field are evaluated for different train speeds. Soil behaviour changes significantly with the track system. Finally, a floating slab track is studied to show how this type of solution leads to a significant vibration reduction for surface tracks. 相似文献
10.
《Advances in water resources》2001,24(6):621-630
One of the more advanced approaches for simulating groundwater flow in fractured porous media is the discrete-fracture approach. This approach is limited by the large computational overheads associated with traditional modeling methods. In this work, we apply the Lanczos reduction method to the modeling of groundwater flow in fractured porous media using the discrete-fracture approach. The Lanczos reduction method reduces a finite element equation system to a much smaller tridiagonal system of first-order differential equations. The reduced system can be solved by a standard tridiagonal algorithm with little computational effort. Because solving the reduced system is more efficient compared to solving the original system, the simulation of groundwater flow in discretely fractured media using the reduction method is very efficient. The proposed method is especially suitable for the problem of large-scale and long-term simulation. In this paper, we develop an iterative version of Lanczos algorithm, in which the preconditioned conjugate gradient solver based on ORTHOMIN acceleration is employed within the Lanczos reduction process. Additional efficiency for the Lanczos method is achieved by applying an eigenvalue shift technique. The “shift” method can improve the Lanczos system convergence, by requiring fewer modes to achieve the same level of accuracy over the unshifted case. The developed model is verified by comparison with dual-porosity approach. The efficiency and accuracy of the method are demonstrated on a field-scale problem and compared to the performance of classic time marching method using an iterative solver on the original system. In spite of the advances, more theoretical work needs to be carried out to determine the optimal value of the shift before computations are actually carried out. 相似文献
11.
A study on the dynamic response of three-dimensional flexible foundations of arbitrary shape, embedded in a homogenous, isotropic and linear elastic half-space is presented. Both massive and massless foundations are considered. The soil-foundation system is subjected to externally applied forces, and/or to obliquely incident seismic waves. The numerical method employed is a combination of the frequency domain Boundary Element Method, which is used to simulate the elastic soil medium, and the Finite Element Method, on the basis of which the stiffness matrix of the foundation is obtained. The foundation and soil media are combined by enforcing compatibility and equilibrium conditions at their common interface. Both relaxed and completely bonded boundary conditions are considered. The accuracy of the proposed methodology is partially verified through comparison studies with results reported in the literature for rigid embedded foundations. 相似文献
12.
Ilya Prutkin Peter VajdaRobert Tenzer Miroslav Bielik 《Journal of Applied Geophysics》2011,75(3):472-478
We present a novel methodology for 3D gravity/magnetic data inversion. It combines two algorithms for preliminary separation of sources and an original approach to 3D inverse problem solution. The first algorithm is designed to separate sources in depth and to remove the shallow ones. It is based on subsequent upward and downward data continuation. For separation in the lateral sense, we approximate the given observed data by the field of several 3D line segments. For potential field data inversion we apply a new method of local corrections. The method is efficient and does not require trial-and-error forward modeling. It allows retrieving unknown 3D geometry of anomalous objects in terms of restricted bodies of arbitrary shape and contact surfaces. For restricted objects, we apply new integral equations of gravity and magnetic inverse problems. All steps of our methodology are demonstrated on the Kolarovo gravity anomaly in the Danube Basin of Slovakia. 相似文献
13.
本文实现了一种面向目标自适应海洋可控源电磁三维矢量有限元方法.为满足三维复杂电性结构模拟的需求,网格剖分采用非结构化六面体.在组装刚度矩阵之后,形成的大型复数线性方程组分解为等价的实数形式,利用带预条件的广义最小残差法进行求解.在获得微分方程的解之后,为提高解的准确性,通过面向目标的自适应误差估计来指示网格细化,重点加密能使观测点数值模拟精度提高的网格.对于大规模三维数据,为了使模型空间的并行计算达到均衡负载的效果,我们使用METIS函数库来进行网格计算任务量的划分.最后,通过对比一维解析解与三维自适应矢量有限元计算结果,验证了程序的正确性;通过自适应过程中误差指示子的分布,验证了面向目标自适应的有效性;通过对三维复杂模型进行均衡负载下的并行计算,测试了程序的可扩展性. 相似文献
14.
In this paper, we develop a methodology for early detection of potential CO2 leakage from geological storage formations using pressure and surface-deformation anomalies. The basic idea is based on the fact that leakage-induced pressure signals travel much faster than the migrating CO2; thus such anomalies may be detected early enough for risk management measures taking effect in avoiding substantial CO2 leaks. The early detection methodology involves automatic inversion of anomalous brine leakage signals with efficient forward pressure and surface-deformation modeling tools to estimate the location and permeability of leaky features in the caprock. We conduct a global sensitivity analysis to better understand under which conditions pressure anomalies can be clearly identified as leakage signals, and evaluate signal detectability for a broad parameter range considering different detection limits and levels of data noise. The inverse methodology is then applied to two synthetic examples of idealized two-aquifer-and-one aquitard storage systems, with an injection well and a leaky well, for different monitoring scenarios. In Example 1, only pressure data at the monitoring and injection wells are used for leakage detection. Our results show that the accuracy of leakage detection greatly depends on the level of pressure data noise. In Example 2, joint inversion of pressure and surface-deformation measurements significantly improves the speed of convergence toward the true solution of the leakage parameters and enables early leakage detection. In both examples, successful detection is achieved when two monitoring wells are appropriately placed within up to 4 km from the leaky well. 相似文献
15.
在前人工作的基础上,本文推导了电导率任意各向异性介质的海洋可控源电磁三维谱元法正演方程.采用一次场/二次场分离算法结合混合阶矢量基函数,可以有效避免源点的奇异性的影响,从而提高数值解的精度.采用任意六面体单元离散研究区域,有利于模拟复杂地形和地电结构.利用不完全LU分解的Induced Dimension Reduction(IDR(s))迭代算法求解线性方程组,有效地提高了求解的效率.设计典型的地电模型进行正演计算,并将计算结果与有限元解进行对比,对比结果表明本文提出的基于混合阶矢量基函数的海洋可控源电磁三维谱元数值模拟算法是正确的、有效的.本文算法具有良好的通用性,可推广用于电导率呈任意各向异性的陆地电磁、井中电磁等数值模拟研究. 相似文献
16.
本文提出了一种计算不规则起伏地形中SH波散射的有效方法——局域边界元法.本方法基于传统边界元法,为计算复杂地表散射问题提供了一种更加高效的解决方案.根据地震波满足的边界积分方程中牵引力格林函数的特性,我们将自由边界分解成水平部分和起伏部分.通过公式推导,可将水平部分的位移由起伏部分的位移通过格林函数线性叠加表示,因此只需对起伏部分的位移进行直接求解,从而极大地减少了待求解的未知数个数,显著提高了计算效率.通过与半圆形山谷SH波平面波入射的解析解比较,验证了方法的正确性.数值模型比较显示,局域边界元模拟结果与传统边界元数值解完全吻合,但是大幅提高了计算效率.因此,局域边界元法可以作为模拟不规则地形中地震波散射的有效工具. 相似文献
17.
Classical least‐squares techniques (Moore–Penrose pseudoinverse) are covariance based and are therefore unsuitable for the solution of very large‐scale linear systems in geophysical inversion due to the need of diagonalisation. In this paper, we present a methodology to perform the geophysical inversion of large‐scale linear systems via the discrete wavelet transform. The methodology consists of compressing the linear system matrix using the interesting properties of covariance‐free orthogonal transformations, to design an approximation of the Moore–Penrose pseudoinverse. We show the application of the discrete wavelet transform pseudoinverse to well‐conditioned and ill‐conditioned linear systems. We applied the methodology to a general‐purpose linear problem where the system matrix has been generated using geostatistical simulation techniques and also to a synthetic 2D gravimetric problem with two different geological set‐ups, in the noise‐free and noisy cases. In both cases, the discrete wavelet transform pseudoinverse can be applied to the original linear system and also to the linear systems of normal equations and minimum norm. The results are compared with those obtained via the Moore–Penrose and the discrete cosine transform pseudoinverses. The discrete wavelet transform and the discrete cosine transform pseudoinverses provide similar results and outperform the Moore–Penrose pseudoinverse, mainly in the presence of noise. In the case of well‐conditioned linear systems, this methodology is more efficient when applied to the least‐squares system and minimum norm system due to their higher condition number that allows for a more efficient compression of the system matrix. Also, in the case of ill‐conditioned systems with very high underdetermined character, the application of the discrete cosine transform to the minimum norm solution provides very good results. Both solutions might differ on their regularity, depending on the wavelet family that is adopted. These methods have a general character and can be applied to solve any linear inverse problem arising in technology, particularly in geophysics, and also to non‐linear inversion by linearisation of the forward operator. 相似文献
18.
Manuel Carbajal-Romero Esteban Flores-Mendez Norberto Flores-Guzmán Jaime Núñez-Farfán Enrique Olivera-Villaseñor Francisco José Sánchez-Sesma 《Geofísica Internacional》2013,52(1):21-30
The present work shows the propagation of Scholte interface waves at the boundary of a fluid in contact with an elastic solid, for a broad range of solid materials. It has been demonstrated that by an analysis of diffracted waves in a fluid it is possible to infer the mechanical properties of the elastic solid medium, specifically, its propagation velocities. For this purpose, the diffracted wave field of pressures and displacements, due to an initial wave of pressure in the fluid, are expressed using boundary integral representations, which satisfy the equation of motion. The source in the fluid is represented by a Hankel’s function of second kind and zero order. The solution to this wave propagation problem is obtained by means of the Indirect Boundary Element Method, which is equivalent to the well-known Somigliana representation theorem. The validation of the results is carried out by using the Discrete Wave Number Method and the Spectral Element Method. Firstly, we show spectra of pressures that illustrate the behavior of the fluid for each solid material considered, then, we apply the Fast Fourier Transform to show results in time domain. Snapshots to exemplify the emergence of Scholte’s waves are also included. 相似文献
19.
The level set methodology for time-optimal path planning is employed to predict collision-free and fastest-time trajectories for swarms of underwater vehicles deployed in the Philippine Archipelago region. To simulate the multiscale ocean flows in this complex region, a data-assimilative primitive-equation ocean modeling system is employed with telescoping domains that are interconnected by implicit two-way nesting. These data-driven multiresolution simulations provide a realistic flow environment, including variable large-scale currents, strong jets, eddies, wind-driven currents, and tides. The properties and capabilities of the rigorous level set methodology are illustrated and assessed quantitatively for several vehicle types and mission scenarios. Feasibility studies of all-to-all broadcast missions, leading to minimal time transmission between source and receiver locations, are performed using a large number of vehicles. The results with gliders and faster propelled vehicles are compared. Reachability studies, i.e., determining the boundaries of regions that can be reached by vehicles for exploratory missions, are then exemplified and analyzed. Finally, the methodology is used to determine the optimal strategies for fastest-time pick up of deployed gliders by means of underway surface vessels or stationary platforms. The results highlight the complex effects of multiscale flows on the optimal paths, the need to utilize the ocean environment for more efficient autonomous missions, and the benefits of including ocean forecasts in the planning of time-optimal paths. 相似文献