首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 468 毫秒
1.
A method is presented for coupling cubic‐order quadrilateral finite elements with the finite side of a new coordinate ascent hierarchical infinite element. At a common side shared by a hierarchical infinite element and an arbitrary number of finite elements, the displacements are minimized in the least square sense with respect to the degrees‐of‐freedom of the finite elements. This leads to a set of equations that relate the degrees‐of‐freedom of the finite and hierarchical infinite elements on the shared side. The method is applied to a non‐homogeneous cross‐anisotropic half‐space subjected to a non‐uniform circular loading with Young's and shear moduli varying with depth according to the power law. A constant mesh constructed from coupled finite and hierarchical infinite elements is used and convergence is sought simply by increasing the degree of the interpolating polynomial. The displacements and stresses produced by conical and parabolic circular loads applied on the surface are obtained. The efficiency of the proposed method is demonstrated through convergence and comparison studies. New results produced by a frusto‐conical circular load applied on the surface of a half‐space made up of heavily consolidated London clay are provided. The non‐homogeneity parameter and degree of anisotropy are shown to influence the soil response. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

2.
Analysis of large deformation of geomaterials subjected to time‐varying load poses a very difficult problem for the geotechnical profession. Conventional finite element schemes using the updated Lagrangian formulation may suffer from serious numerical difficulties when the deformation of geomaterials is significantly large such that the discretized elements are severely distorted. In this paper, an operator‐split arbitrary Lagrangian–Eulerian (ALE) finite element model is proposed for large deformation analysis of a soil mass subjected to either static or dynamic loading, where the soil is modelled as a saturated porous material with solid–fluid coupling and strong material non‐linearity. Each time step of the operator‐split ALE algorithm consists of a Lagrangian step and an Eulerian step. In the Lagrangian step, the equilibrium equation and continuity equation of the saturated soil are solved by the updated Lagrangian method. In the Eulerian step, mesh smoothing is performed for the deformed body and the state variables obtained in the updated Lagrangian step are then transferred to the new mesh system. The accuracy and efficiency of the proposed ALE method are verified by comparison of its results with the results produced by an analytical solution for one‐dimensional finite elastic consolidation of a soil column and with the results from the small strain finite element analysis and the updated Lagrangian analysis. Its performance is further illustrated by simulation of a complex problem involving the transient response of an embankment subjected to earthquake loading. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

3.
A fully coupled formulation of a hydro‐thermo‐poro‐mechanical model for a three‐phase black oil reservoir model is presented. The model is based upon the approach proposed by one of the authors which fully couples geomechanical effects to multiphase flow. Their work is extended here to include non‐isothermal effects. The gas phase contribution to the energy equation has been neglected based on a set of assumptions. The coupled formulation given herein differs in several ways when compared to the earlier work and an attempt is made to link the flow based formulation and mixture theory. The Finite Element Method is employed for the numerical treatment and essential algorithmic implementation is discussed. Numerical examples are presented to provide further understanding of the current methodology. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

4.
Hydraulic fracturing (HF) of underground formations has widely been used in different fields of engineering. Despite the technological advances in techniques of in situ HF, the industry uses semi‐analytical tools to design HF treatment. This is due to the complex interaction among various mechanisms involved in this process, so that for thorough simulations of HF operations a fully coupled numerical model is required. In this study, using element‐free Galerkin (EFG) mesh‐less method, a new formulation for numerical modeling of hydraulic fracture propagation in porous media is developed. This numerical approach, which is based on the simultaneous solution of equilibrium and continuity equations, considers the hydro‐mechanical coupling between the crack and its surrounding porous medium. Therefore, the developed EFG model is capable of simulating fluid leak‐off and fluid lag phenomena. To create the discrete equation system, the Galerkin technique is applied, and the essential boundary conditions are imposed via penalty method. Then, the resultant constrained integral equations are discretized in space using EFG shape functions. For temporal discretization, a fully implicit scheme is employed. The final set of algebraic equations that forms a non‐linear equation system is solved using the direct iterative procedure. Modeling of cracks is performed on the basis of linear elastic fracture mechanics, and for this purpose, the so‐called diffraction method is employed. For verification of the model, a number of problems are solved. According to the obtained results, the developed EFG computer program can successfully be applied for simulating the complex process of hydraulic fracture propagation in porous media. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

5.
A hierarchical mathematical model for analyses of coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes at high temperature is presented. The concretes are modelled as unsaturated deforming reactive porous media filled with two immiscible pore fluids, i.e. the gas mixture and the liquid mixture, in immiscible–miscible levels. The thermo‐induced desalination process is particularly integrated into the model. The chemical effects of both the desalination and the dehydration processes on the material damage and the degradation of the material strength are taken into account. The mathematical model consists of a set of coupled, partial differential equations governing the mass balance of the dry air, the mass balance of the water species, the mass balance of the matrix components dissolved in the liquid phases, the enthalpy (energy) balance and momentum balance of the whole medium mixture. The governing equations, the state equations for the model and the constitutive laws used in the model are given. A mixed weak form for the finite element solution procedure is formulated for the numerical simulation of chemo‐thermo‐hygro‐mechanical behaviours. Special considerations are given to spatial discretization of hyperbolic equation with non‐self‐adjoint operator nature. Numerical results demonstrate the performance and the effectiveness of the proposed model and its numerical procedure in reproducing coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes subjected to fire and thermal radiation. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

6.
A new formulation of the element‐free Galerkin (EFG) method is developed for solving coupled hydro‐mechanical problems. The numerical approach is based on solving the two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Spatial variables in the weak form, i.e. displacement increment and pore water pressure increment, are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on a penalty method. Numerical stability of the developed formulation is examined in order to achieve appropriate accuracy of the EFG solution for coupled hydro‐mechanical problems. Examples are studied and compared with closed‐form or finite element method solutions to demonstrate the validity of the developed model and its capabilities. The results indicate that the EFG method is capable of handling coupled problems in saturated porous media and can predict well both the soil deformation and variation of pore water pressure over time. Some guidelines are proposed to guarantee the accuracy of the EFG solution for coupled hydro‐mechanical problems. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

7.
  相似文献   

8.
The diffusion equation of suspended sediment concentration in a wide sediment‐laden stream flow is dependent on the vertical gradient of streamwise velocity and the sediment diffusivity. This study aims at investigating the influence of the streamwise velocity laws on the suspended sediment concentration distributions, resulting from the solution of the diffusion equation. Firstly, the sediment concentration distributions are obtained numerically from the solution of the diffusion equation using different velocity laws and compared with the experimental data. It is found that the power‐law approximation produces good computational results for the concentration distributions. The accuracy of using a power‐law velocity model is comparable with the results obtained from other classical velocity laws, namely log‐law, log wake‐law and stratified log‐law. Secondly, a novel analytical solution is proposed for the determination of sediment concentration distribution, where a power‐law, wall‐concentration profile is coupled with a concentration wake function. The power‐law model (for velocity and concentration) is calibrated using the experimental data, and then a generalized wake function is obtained by choosing a suitable law. The developed power‐law model involving the wake function adjusted by an exponent predicts the sediment concentration distributions quite satisfactorily. Finally, a new explicit formula for the suspended‐load transport rate is derived from the proposed theory, where numerical computation of integrals, as needed in the Einstein theory, is avoided.  相似文献   

9.
The flow of water–kaolinite mixtures exhibits a non‐Newtonian nature that differs from the flow of Newtonian fluid. The varying viscosities and shear history of non‐Newtonian fluid flows necessitate the use of a rheology model in moving particle semi‐implicit (MPS) for the numerical studies. On the other hand, the Lagrangian method has the advantage of handling free surface flows with large deformation and fragmentation. This study proposes a mesh‐free Lagrangian method, namely, the MPS method, together with a simple rheology model to investigate the non‐Newtonian free surface flows. The rheological parameters required in the rheology model are determined based upon experiments. The proposed method is applied to a water–kaolinite mixture collapse problem and is proved to be capable of reproducing the significant flow features observed in laboratory experiments. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

10.
The propagation characteristic of Rayleigh waves in a fluid‐saturated non‐homogeneous poroelastic half‐plane is addressed. Based on Biot's theory for fluid‐saturated media, which takes the inertia, fluid viscosity, mechanical coupling, compressibility of solid grains, and fluid into account, the dispersion equations of Rayleigh waves in fluid‐saturated non‐homogeneous soils/rocks are established. By considering the shear modulus of solid skeleton variation with depth exponentially, a small parameter, which reflects the relative change of shear modulus, is introduced. The asymptotic solution of the dispersion equation expressing the relationship between the phase velocity and wave number is obtained by using the perturbation method. In order to analyze the effects of non‐homogeneity on the propagation characteristic of Rayleigh waves, the variation of the phase velocity with the wave number is presented graphically and discussed through numerical examples. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

11.
This paper presents a u‐p (displacement‐pressure) semi‐Lagrangian reproducing kernel (RK) formulation to effectively analyze landslide processes. The semi‐Lagrangian RK approximation is constructed based on Lagrangian discretization points with fixed kernel supports in the current configuration. As a result, it tracks state variables at discretization points while allowing extreme deformation and material separation that is beyond the capability of Lagrangian formulations. The u‐p formulation following Biot theory is incorporated into the formulation to describe poromechanics of saturated geomaterials. In addition, a stabilized nodal integration method to ensure stability of the domain integration and kernel contact algorithms to model contact between bodies are introduced in the u‐p semi‐Lagrangian RK formulation. The proposed method is verified with several numerical examples and validated with an experimental result and the field data of an actual landslide.  相似文献   

12.
The early Pleistocene clastic succession of the Peri‐Adriatic basin, eastern central Italy, records the filling of a series of piggyback sub‐basins that formed in response to the development of the eastward‐verging Apennine fold‐thrust belt. During the Gelasian (2·588 to 1·806 Ma), large volumes of Apennine‐derived sediments were routed to these basins through a number of slope turbidite systems. Using a comprehensive outcrop‐based dataset, the current study documents the depositional processes, stratigraphic organization, foraminiferal age and palaeodepth, and stratigraphic evolution of one of these systems exposed in the surroundings of the Castignano village. Analysis of foraminiferal assemblages consistently indicates Gelasian deposition in upper bathyal water depths. Sediments exposed in the study area can be broken into seven main lithofacies, reflecting specific gravity‐induced depositional elements and slope background deposition: (i) clast‐supported conglomerates (conglomerate channel‐fill); (ii) amalgamated sandstones (late stage sandstone channel‐fill); (iii) medium to thick‐bedded tabular sandstones (frontal splay sandstones); (iv) thin to thick‐bedded channelized sandstones (sandy channel‐fill); (v) medium to very thin‐bedded sandstones and mudstones (levée‐overbank deposits); (vi) pebbly mudstones and chaotic beds (mudstone‐rich mass‐transport deposits); and (vii) massive mudstones (hemipelagic deposits). Individual lithofacies combine vertically and laterally to form decametre‐scale, disconformably bounded, fining‐upward lithofacies successions that, in turn, stack to form slope valley fills bounded by deeply incised erosion surfaces. A hierarchical approach to the physical stratigraphy of the slope system indicates that it has evolved through multiple cycles of waxing then waning flow energy at multiple scales and that its packaging can be described in terms of a six‐fold hierarchy of architectural elements and bounding surfaces. In this scheme, the whole system (sixth‐order element) is comprised of three distinct fifth‐order stratigraphic cycles (valley fills), which define sixth‐order initiation, growth and retreat phases of slope deposition, respectively; they are separated by discrete periods of entrenchment that generated erosional valleys interpreted to record fifth‐order initiation phases. Backfilling of individual valleys progressed through deposition of two vertically stacked lithofacies successions (fourth‐order elements), which record fifth‐order growth and retreat phases. Fourth‐order initiation phases are represented by erosional surfaces bounding lithofacies successions. The component lithofacies (third‐order element) record fourth‐order growth and retreat phases. Map trends of erosional valleys and palaeocurrent indicators converge to indicate that the sea floor bathymetric expression of a developing thrust‐related anticline markedly influenced the downslope transport direction of gravity currents and was sufficient to cause a major diversion of the turbidite system around the growing structure. This field‐based study permits the development of a sedimentological model that predicts the evolutionary style of mixed coarse‐grained and fine‐grained turbidite slope systems, the internal distribution of reservoir and non‐reservoir lithofacies within them, and has the potential to serve as an analogue for seismic or outcrop‐based studies of slope valley fills developed in actively deforming structural settings and under severe icehouse regimes.  相似文献   

13.
A Lagrangian particle‐based method, smooth particle hydrodynamics (SPH), is used in this paper to model the flow of self‐compacting concretes (SCC) with or without short steel fibres. An incompressible SPH method is presented to simulate the flow of such non‐Newtonian fluids whose behaviour is described by a Bingham‐type model, in which the kink in the shear stress vs shear strain rate diagram is first appropriately smoothed out. The viscosity of the SCC is predicted from the measured viscosity of the paste using micromechanical models in which the second phase aggregates are treated as rigid spheres and the short steel fibres as slender rigid bodies. The basic equations solved in the SPH are the incompressible mass conservation and Navier–Stokes equations. The solution procedure uses prediction–correction fractional steps with the temporal velocity field integrated forward in time without enforcing incompressibility in the prediction step. The resulting temporal velocity field is then implicitly projected on to a divergence‐free space to satisfy incompressibility through a pressure Poisson equation derived from an approximate pressure projection. The results of the numerical simulation are benchmarked against actual slump tests carried out in the laboratory. The numerical results are in excellent agreement with test results, thus demonstrating the capability of SPH and a proper rheological model to predict SCC flow and mould‐filling behaviour. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

14.
This paper presents a coupled hydro‐mechanical formulation for the simulation of non‐planar three‐dimensional hydraulic fractures. Deformation in the rock is modeled using linear elasticity, and the lubrication theory is adopted for the fluid flow in the fracture. The governing equations of the fluid flow and elasticity and the subsequent discretization are fully coupled. A Generalized/eXtended Finite Element Method (G/XFEM) is adopted for the discretization of the coupled system of equations. A Newton–Raphson method is used to solve the resulting system of nonlinear equations. A discretization strategy for the fluid flow problem on non‐planar three‐dimensional surfaces and a computationally efficient strategy for handling time integration combined with mesh adaptivity are also presented. Several three‐dimensional numerical verification examples are solved. The examples illustrate the generality and accuracy of the proposed coupled formulation and discretization strategies. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

15.
A characteristic finite element (CFE) algorithm for modelling contaminant transport problems coupled with non‐linear competitive adsorption is presented. An alternative algorithm, termed as the transport‐equilibrium Petrov–Galerkin (TEPG) methods in this paper, is also presented for comparison. Through analyses of a number of examples with Peclet number ranging from zero to infinity, it is shown that the CFE algorithm is very competitive with the middle–point TEPG method in terms of accuracy, stability and efficiency. The fully explicit and fully implicit TEPG methods are found to be less appropriate for transport problems coupled with non‐linear equilibrium equations. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

16.
This paper presents a non‐linear soil–structure interaction (SSI) macro‐element for shallow foundation on cohesive soil. The element describes the behaviour in the near field of the foundation under cyclic loading, reproducing the material non‐linearities of the soil under the foundation (yielding) as well as the geometrical non‐linearities (uplift) at the soil–structure interface. The overall behaviour in the soil and at the interface is reduced to its action on the foundation. The macro‐element consists of a non‐linear joint element, expressed in generalised variables, i.e. in forces applied to the foundation and in the corresponding displacements. Failure is described by the interaction diagram of the ultimate bearing capacity of the foundation under combined loads. Mechanisms of yielding and uplift are modelled through a global, coupled plasticity–uplift model. The cyclic model is dedicated to modelling the dynamic response of structures subjected to seismic action. Thus, it is especially suited to combined loading developed during this kind of motion. Comparisons of cyclic results obtained from the macro‐element and from a FE modelization are shown in order to demonstrate the relevance of the proposed model and its predictive ability. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

17.
A mathematical treatment is presented for the forced vertical vibration of a padded annular footing on a layered viscoelastic half‐space. On assuming a depth‐independent stress distribution for the interfacial buffer, the set of triple integral equations stemming from the problem is reduced to a Fredholm integral equation of the second kind. The solution method, which is tailored to capture the stress concentrations beneath footing edges, is highlighted. To cater to small‐scale geophysical applications, the model is used to investigate the near‐field effects of ground‐loading system interaction in dynamic geotechnical and pavement testing. Numerical results indicate that the uniform‐pressure assumption for the contact load between the composite disc and the ground which is customary in dynamic plate load testing may lead to significant errors in the diagnosis of subsurface soil and pavement conditions. Beyond its direct application to non‐intrusive site characterization, the proposed solution can be used in the seismic analysis of a variety of structures involving annular foundation geometries. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

18.
Research in 2006 and 2007 dealing with laser ablation‐(multicollector)‐inductively coupled plasma‐mass spectrometry, LA‐(MC)‐ICP‐MS, involved studies concerned with optimising the technique itself, as well as applying the method to a variety of problems in the Earth sciences. The causes of elemental and isotopic fractionation produced during laser ablation continues to be of considerable interest, with evidence mounting that processes occurring both at the ablation site and in the argon plasma of the ICP are culpable. There is growing excitement in the use of femtosecond lasers for LA‐(MC)‐ICP‐MS, with the hope that they reduce or eliminate melting and non‐congruent volatilisation at the ablation site and thus approach stoichiometric sampling. Ablation chamber design emerged as a serious concern, particularly with respect to achieving the rapid washout needed for fine‐scale compositional mapping of geological objects. LA‐MC‐ICP‐MS provided data for a wide range of isotopic systems, especially hafnium, but also B, S, Mg, Cu, Fe, Sr, Nd, Pb and U. Measurement uncertainties in LA‐ICP‐MS were discussed by several researchers, and are critically reviewed here ‐ total uncertainties for trace element concentration measurements of silicates including errors on the calibration values of common reference materials are ～10% (95% confidence limits), though the precision of individual spot measurements (50 to 100 μm) is much better, ～1% RSD, using a 193 nm laser and a sector field‐ICP‐MS. LA‐ICP‐MS U‐Pb ages for zircon and other U‐rich accessory phases are claimed by most geoanalysts to have 2s uncertainties of ～0.7 and 1.3% respectively but the actual accuracy of the method is probably only as good as ～2% (2s), when uncertainties associated with laser‐induced Pb/U fractionation are included.  相似文献   

19.
The study of surface wave in a layered medium has their possible application in geophysical prospecting. In the present work, dispersion equation for torsional wave in an inhomogeneous isotropic layer between inhomogeneous isotropic half‐spaces has been derived. Two cases are discussed separately for torsional wave propagation in inhomogeneous layer between homogeneous and non‐homogeneous half‐spaces, respectively. Further, two possible modes for torsional wave propagation are obtained in case of inhomogeneous layer sandwiched between non‐homogeneous half‐spaces. Closed form solutions for displacement in the layer and half‐spaces are obtained in each case. The study reveals that the layer width, layer inhomogeneity, frequency of inhomogeneity, as well as inhomogeneity in the half‐space has significant effect on the propagation of torsional surface waves. Displacement and implicit dispersion equation for torsional wave velocities are expressed in terms of Heun functions and their derivatives. Effects of inhomogeneity on torsional wave velocity are also discussed graphically by plotting the dimensionless phase velocity against dimensionless and scaled wave number for different values of inhomogeneity parameter. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

20.
In this article, we evaluate geomechanics of fluid injection from a fully penetrating vertical well into an unconsolidated formation confined with stiff seal rocks. The coupled behavior of an isotropic, homogeneous sand layer is studied under injection pressures that are high enough to induce plasticity yet not fracturing. Propagation of the significant influence zone surrounding the injection borehole, quantified by the extent of the plastic domain in the elasto‐plastic model, is examined for the first time. First, a new fully coupled axisymmetric numerical model is developed. A comprehensive assessment is performed on pore pressures, stresses/strains, and failure planes during the entire transient period of an injection cycle. Results anticipate existence of five distinctive zones in terms of plasticity state: liquefaction at wellbore; two inner plastic domains surrounding the wellbore, where failure occurs along two planes and major principal stress is in vertical direction; remaining of the plastic domain, where formation fails along one plane and major principal stress is in radial direction; and a non‐plastic region. Failure mechanism at the wellbore is found to be shear followed by liquefaction. Next, a novel methodology is proposed based on which new weakly coupled poro‐elasto‐plastic analytical solutions are derived for all three stress/strain components. Unlike previous studies, extension of the plastic zone is obtained as a function of injection pressure, incorporating plasticity effects on the subsequent elastic domain. Solutions, proven to be a good approximation of numerical simulations, offer a huge advantage as the run time of coupled numerical simulations is considerably long. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号