共查询到20条相似文献，搜索用时 31 毫秒
1.
Over the last decade, researchers in the discontinuous deformation analysis (DDA) community have dedicated a great deal of effort to document the accuracy of the method by performing validation studies. This paper contains a summary of more than 100 published and unpublished validation studies which comprise the body of DDA validation information to which the authors have access. The studies are grouped into three general categories: (a) validation with respect to analytical solutions, (b) validation with respect to results of other numerical techniques, and (c) validation with respect to laboratory and field data. Three general techniques for validation are described: qualitative assessment visually examining runtime behaviour of simulations, semi‐quantitative assessment comparing numerical results of simulations, and quantitative where numerical simulation results are evaluated in detail with respect to similar analytical, laboratory or field results. We find that for many of the problems addressed by the papers in this review, DDA performs more than adequately for engineering analysis. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
2.
An innovative approach is presented, in which the discontinuous deformation analysis (DDA) method is used to estimate historic ground motions by back analysis of unique structural failures in archaeological sites. Two archaeological sites in Israel are investigated using this new approach and results are presented in terms of displacement evolution of selected structural elements in the studied masonry structure. The response of the structure is studied up to the point of incipient failure, in a mechanism similar to the one observed in the field. Structural response is found to be very sensitive to dynamic parameters of the loading function such as amplitude and frequency. Prior to back analysis of case studies, two validations are presented. Both compare the performance of DDA with analytical solutions and present strong agreement between the two. Using comprehensive sensitivity analyses, the most likely peak ground acceleration (PGA) and frequency that must have driven the observed block displacements are found for the two case studies—the Nabatean city of Mamshit and the medieval fortress of Nimrod in southern and northern Israel, respectively. It is found that horizontal peak ground accelerations (HPGA) of 0.5g and 1g were required to generate the observed deformations in Mamshit and Nimrod, respectively. Although these might seem too high, considering structural and topographic amplifications it is concluded that the analyses suggest ground motions of 0.2g at a frequency of 1.5 Hz for Mamshit and up to 0.4g at a frequency of 1 Hz for Nimrod. These values provide constraints on the seismic risk associated with these regions as appears in the local building code using a completely independent approach. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
3.
In this paper, we present a numerical model for simulating two-phase (oil–water and air–water) incompressible and immiscible flow in porous media. The mathematical model which is based on a fractional flow formulation is formed of two nonlinear partial differential equations: a mean pressure equation and a water saturation equation. These two equations can be solved in a sequential manner. Two numerical methods are used to discretize the equations of the two-phase flow model: mixed hybrid finite elements are used to treat the pressure equation, h-based Richards' equation and the diffusion term in the saturation equation, the advection term in the saturation equation is treated with the discontinuous finite elements. We propose a better way to calculate the nonlinear coefficients contained in our equations on each element of the discretized domain. In heterogeneous porous media, the saturation becomes discontinuous at the interface between two porous media. We show in this paper how to use the capillary pressure–saturation relationship in order to handle the saturation jump in the mixed hybrid finite element method. The two-phase flow simulator is verified against analytical solutions for some flow problems treated by other authors. 相似文献
4.
Phase transformation and crystal growth in nanoparticles may happen via mechanisms distinct from those in bulk materials. We combine experimental studies of as-synthesized and hydrothermally coarsened
titania (TiO2) and zinc sulfide (ZnS) with thermodynamic analysis, kinetic modeling and molecular dynamics (MD) simulations. The samples
were characterized by transmission electron microscopy, X-ray diffraction, synchrotron X-ray absorption and scattering, and
UV-vis spectroscopy. At low temperatures, phase transformation in titania nanoparticles occurs predominantly via interface nucleation at particle–particle contacts. Coarsening and crystal growth of titania nanoparticles can be described
using the Smoluchowski equation. Oriented attachment-based crystal growth was common in both hydrothermal solutions and under
dry conditions. MD simulations predict large structural perturbations within very fine particles, and are consistent with
experimental results showing that ligand binding and change in aggregation state can cause phase transformation without particle
coarsening. Such phenomena affect surface reactivity, thus may have important roles in geochemical cycling. 相似文献
5.
6.
Modeling the connectivity and intersection of hydraulically loaded cracks with in situ fractures in rock 
下载免费PDF全文
下载免费PDF全文 The interaction of an advancing hydraulically loaded crack and in situ fracture network can yield highly complex patterns. We model the connectivity of cells in a finite element domain and in a fracture network by a simplicial complex data structure. The complete adjacency information between cells is determined by one level down facet and one level up cofacet neighborhood information. Combined with a disjoint set data structure, explicit algorithms are derived to efficiently track network connectivity and load transfer between independent fracture sets. We also propose an approach to regularize the application of hydraulic load to newly intersected in situ cracks to smoothen the transition of pressure on intersected cracks from ambient to hydraulic pressure and to avoid the sudden loading of the entire length of these cracks. Numerical results demonstrate the performance of crack connectivity and load transfer models and the effect of regularization model. The results show that as the angle between an incoming hydraulically loaded crack and an in situ crack increases, the effect of in situ crack shifts from slight realignment to diversion/offsetting of the loaded crack. As the angle difference approaches the normal angle, the loaded crack tends to penetrate through the in situ crack. The proposed schemes are also used for transient simulation of 2D reservoirs with multiple perforations surrounded by in situ cracks with and without a bias in the distribution of their orientation. It was shown that from 2 perforations with angles closer to in situ cracks at low loading rates to all perforations at higher loading rates can result in active hydraulic crack propagation. The h‐adaptive method of asynchronous space‐time discontinuous Galerkin method is used to exactly track complex fracture patterns in these dynamic fracture simulations. 相似文献
7.
Fast multipole displacement discontinuity method (FM‐DDM) for geomechanics reservoir simulations 下载免费PDF全文
下载免费PDF全文 The displacement discontinuity method (DDM) is frequently used in geothermal and petroleum applications for modeling the behavior of fractures in linear‐elastic rocks. The DDM requires O(N2) memory and O(N3) floating point operations (where N is the number of unknowns) to construct the coefficient matrix and solve the linear system of equations by direct methods. Therefore, the conventional implementation of the DDM is not computationally efficient for very large systems of cracks, often limiting its application to small‐scale problems. This work presents an approach for solving large‐scale fracture problems using the fast multipole method (FMM). The approach uses both the DDM and a kernel‐independent version of the FMM along with a preconditioned generalized minimal residual algorithm to accelerate the solution of linear systems of equations using desktop computers. Using the fundamental solutions for constant displacement discontinuity in a two‐dimensional elastic medium, several numerical examples involving fracture networks representing fractured reservoirs are treated. Numerical results show good agreement with analytical solutions and demonstrate the efficiency of the FMM implementation of the DDM for large‐scale simulations. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
8.
Reducing vertical acceleration in ballast layers of railways to mitigate geometrical disorders in high‐speed lines 下载免费PDF全文
下载免费PDF全文 According to field feedbacks from high‐speed lines (HSL), the increase of train operating speeds is responsible for unusual fast evolving geometrical disorders in ballasted tracks. This paper deals with the search of solutions applicable at the design stage to mitigate these disorders. The starting point of the present work relies on the assumption, comforted by the literature, of a strong correlation between disorders and vertical accelerations in the ballast layer induced by the train passages. This led us focus herein on the calculation and the analysis of accelerations in the railway structure. The vertical accelerations (γz) are computed using the in‐house developed numerical program ViscoRail and on the basis of a reference HSL. These are shown to increase strongly with the train speed attesting to the link between the train speed and the geometrical disorders in ballast. Then, other simulations are run varying some structural parameters to evaluate their impact on the acceleration field γz. In that way, we show that decreasing the stiffness of the mechanical connection between the rails and the ballast, increasing the moment of inertia of the rails or the Young modulus of the sub‐ballast layer, leads to a decrease of γz and could provide solutions for the design of future HSL. The solution consisting in the incorporation of an asphalt sub‐ballast layer, as already experimented on sites, is finally examined in more details. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
9.
Tilmann Gneiting 《Mathematical Geology》1998,30(4):379-390
The turning bands method generates realizations of isotropic Gaussian random fields by means of appropriately summed line processes. For two-dimensional simulations the relation between the isotropic correlation function of the random field and the correlation function to be simulated along the lines is given by an integral equation of Abel type. We present closed form solutions of this integral equation for almost all two-dimensional correlation models encountered in practice and discuss their numerical implementation. As an additional benefit, our tables and illustrations serve as a concise guide to correlation models useful in geostatistics. 相似文献
10.
A modified three‐dimensional discontinuous deformation analysis (3D‐DDA) method is derived using four‐noded tetrahedral elements to improve the accuracy of current 3D‐DDA algorithm in practical applications. The analysis program for the modified 3D‐DDA method is developed in a C++ environment and its accuracy is illustrated through comparisons with several analytical solutions that are available for selected problems. The predicted solutions for these problems using the modified 3D‐DDA approach all show satisfactory agreement with the corresponding analytical results. Results presented in this paper demonstrate that the modified 3D‐DDA method with discontinuous modeling capabilities offers a useful computational tool to determine stresses and deformations in practical problems involving fissured elastic media with reasonable accuracy. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
11.
卫生填埋场复合边坡地震稳定性和永久变形分析 总被引：11，自引：3，他引：8
采用非连续变形分析方法,对卫生填埋场复合型边坡的地震稳定性和永久变形进行了详细地计算分析,探讨了边坡倾角、土与土工膜之间的摩擦系数、边坡的长度、覆盖土层的厚度以及地震加速度对边坡的稳定性和永久变形的影响。由边坡的破坏形式,讨论了复合型边坡的破坏机理,指出了文献[1]给出的计算公式的适用条件以及局限性,为卫生填埋场复合型衬垫层的稳定设计提供参考。 相似文献
12.
An extensive examination of the discontinuous deformation analysis (DDA) in block dynamic sliding modeling is carried out in this paper. Theoretical solutions for a single block sliding on an arbitrarily inclined plane by applying the horizontal/vertical seismic loadings to the sliding block as acceleration time histories or to the base as constraint displacement time histories are derived. As compared with the theoretical solutions, for a single block sliding, the DDA predicts the sliding displacements and block interaction forces accurately under various base incline angles and friction angles under both the harmonic loadings and a real seismic loading. The vertical seismic component may influence the block sliding displacements to different extent, and the DDA can capture these phenomena successfully and give accurate results. For the calculation of the single block relative sliding, both the theoretical and the DDA solutions indicate that applying the seismic accelerations as constraint displacement time histories (derived by integrating the seismic accelerations twice) to the base is equivalent to applying the seismic accelerations as volume forces to the sliding block in the opposite directions. The DDA modeling also demonstrates that this conclusion still stands for the case of multi‐block sliding. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
13.
Similar to plane strain, axisymmetric stress problem is also highly kinematics constrained. Standard displacement‐based finite element exhibits volumetric locking issue in simulating nearly/fully incompressible material or isochoric plasticity under axisymmetric loading conditions, which severely underestimates the deformation and overestimates the bearing capacity for structural/geotechnical engineering problems. The aim of this paper is to apply variational multiscale method to produce a stabilized mixed displacement–pressure formulation, which can effectively alleviate the volumetric locking issue for axisymmetric stress problem. Both nearly incompressible elasticity and isochoric J2 elastoplasticity are investigated. First‐order 3‐node triangular and 4‐node quadrilateral elements are tested for locking issues. Several representative simulations are provided to demonstrate the performance of the linear elements, which include the convergence study and comparison with closed‐form solutions. A comparative study with pressure Laplacian stabilized formulation is also presented. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
14.
Hongbin Zhan 《Mathematical Geology》1999,31(2):175-193
Closed-form steady-state analytical solutions for capture times are derived considering both pumping–injecting and pumping–pumping double-well systems. Numerical simulations are used to obtain capture times and the numerical results agree with analytical solutions. The influences of regional flows upon the capture zones are analyzed numerically. In a heterogeneously stratified aquifer, the mean and the variance of the capture time are calculated. It is found that hydraulic conductivity distribution plays an important role in determining the mean and the variance of the capture time. 相似文献
15.
In a companion Part I of this paper (Int. J. Numer. Anal. Meth. Geomech. 2008; DOI: 10.1002/nag.735 ), a coupled hydro‐mechanical (HM) formulation for geomaterials with discontinuities based on the finite element method (FEM) with double‐node, zero‐thickness interface elements was developed and presented. This Part II paper includes the numerical solution of basic practical problems using both the staggered and the fully coupled approaches. A first group of simulations, based on the classical consolidation problem with an added vertical discontinuity, is used to compare both the approaches in terms of accuracy and convergence. The monolithic or fully coupled scheme is also used in an application example studying the influence of a horizontal joint in the performance of a reservoir subject to fluid extraction. Results include a comparison with other numerical solutions from the literature and a sensitivity analysis of the mechanical parameters of the discontinuity. Some simulations are also run using both a full non‐symmetric and a simplified symmetric Jacobian matrix. On top of verifying the model developed and its capability to reflect the conductivity changes of the interface with aperture changes, the results presented also lead to interesting observations of the numerical performance of the methods implemented. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
16.
Eva R. Myers 《Physics and Chemistry of Minerals》1998,25(6):465-468
A statistical-thermodynamic approximation for order-disorder phase transitions in aluminosilicate solid solutions is presented.
The approximation involves estimating the number of configurations with long range order parameter Q and short range order parameter σ, using an approximation to the probability that a configuration with the correct amount
of long range order will also have the correct amount of short range order. This estimate is then used to give a free energy
F(Q, σ, x, T), where x is the concentration of Al in the structure, and hence quantities such as T
c
(x) are estimated. The predictions of the model for T
c
(x) and the critical concentration x
c
at which T
c
falls to zero are shown to be in good agreement with the results of Monte Carlo simulations.
Received: 14 May 1997 / Revised, accepted: 2 June 1997 相似文献
17.
Amir M. Abdolahzadeh Benoit Lacroix Vachon Alexandre R. Cabral 《Geotechnical and Geological Engineering》2011,29(5):783-802
One important step in the design of inclined covers with capillary barrier effect (CCBE) is the determination of the water
diversion length (DL). Numerical simulations can predict the DL more precisely than steady-state analytical solutions. Nevertheless,
as simplified methods have always been part of engineering design, the application of analytical solutions with conservative
boundary conditions, may allow engineers to make reasonable predictions, particularly during the pre-feasibility stage of
a project. In this study, a CCBE was designed, constructed and instrumented at the Saint-Tite-des-Caps landfill, Quebec, Canada.
This CCBE included a seepage control layer superimposing a sand-gravel capillary barrier. The seepage control layer was made
up of deinking by-products (DBP), an industrial by-product that was previously disposed of as waste. The capillary barrier
was designed using an adaptation of the Ross analytical solution and the scenario considered was that of steady-state flow
during constant seepage flow applied uniformly at the top of the sand-gravel capillary barrier. Although these conditions
appear simplistic, they were deemed reasonable because placement of the seepage control layer on the top of the capillary
barrier led to very low suctions at the interface, thereby allowing uniform downward seepage rates, limited by the saturated
hydraulic conductivity of the DBP. In this paper, a discussion about the behaviour of the cover system based on 4 years of
field data from several instruments is presented. The challenge of using DBP, more precisely the settlement of the DBP layer
and its impact on k
sat
, is also assessed. The DL was reassessed considering the new k
sat
. A discussion on the validity of employing analytical solutions to determine DL is also presented. This paper illustrates
how certain variables affect the design of inclined CCBEs that include a highly compressible material as seepage control layer. 相似文献
18.
Analytical solutions are presented for fluid‐saturated linear poroelastic beams under pure bending. The stress‐free boundary condition at the lateral surfaces is satisfied in the St Venant's sense and the Beltrami–Michell compatibility conditions are resolved rigorously, rendering the flexure of the beams analytically tractable. Two sets of formulations are derived based on the coupled and uncoupled diffusion equations respectively. The analytical solutions are compared with three‐dimensional finite element simulations. Both sets of analytical formulations are capable of capturing exactly both the initial (undrained) and the steady‐state (fully drained) deflection of the beams. However, the analytical solutions are found to be deficient during the transient phase. The cause for the deficiency of the transient analytical solutions is discussed. The accuracy of the analytical solutions improves as Poisson's ratio and the compressibility of the constituents of the porous beam increase, where the St Venant's edge effect at the lateral surfaces is mitigated. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
19.
Numerical simulations are essential tools for studying tsunami generation and evolution and finite-element (FE) methods are widely used, especially because of their capability in modeling water waves in basins with complex bathymetry and irregular coastlines. This paper presents the numerical simulation of an historical Italian tsunami that affected the Tyrrhenian coasts of Calabria and Sicily on 5 February 1783 following a strong destructive earthquake that was the first of a terrible sequence of seismic shocks terrifying the Calabrian population for more than two months. The numerical model is an FE model based on the nonlinear nondispersive shallow-water approximation of the Navier-Stokes equations. Since FE discretization schemes may lead to solutions undesirably affected by noise over coarse grids, in this study numerical noise is controlled by suitably smoothing the FE solution at regular time steps t
s. The performance of our smoothing algorithm is tested for significant linear cases for which an analytical solution is available. 相似文献
20.
This work presents analytical solutions for displacements caused by three‐dimensional point loads in a transversely isotropic full space, in which transversely isotropic planes are inclined with respect to the horizontal loading surface. In the derivation, the triple Fourier transforms are employed toyield integral expressions of Green's displacement; then, the triple inverse Fourier transforms and residue calculus are performed to integrate the contours. The solutions herein indicate that the displacements are governed by (1) the rotation of the transversely isotropic planes (?), (2) the type and degree of material anisotropy (E/E′, ν/ν′, G/G′), (3) the geometric position (r, φ, ξ) and (4) the types of loading (Px, Py, Pz). The solutions are identical to those of Liao and Wang (Int. J. Numer. Anal. Methods Geomechanics 1998; 22 (6):425–447) if the full space is homogeneous and linearly elastic and the transversely isotropic planes are parallel to the horizontal surface. Additionally, a series of parametric study is conducted to demonstrate the presented solutions, and to elucidate the effect of the aforementioned factors on the displacements. The results demonstrate that the displacements in the infinite isotropic/transversely isotropic rocks, subjected to three‐dimensional point loads could be easily determined using the proposed solutions. Also, these solutions could realistically imitate the actual stratum of loading situations in numerous areas of engineering. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献