共查询到20条相似文献，搜索用时 31 毫秒
1.
Non-equilibrium thermodynamics and Biot poro-elasticity have been combined to give a coupled hydro-mechanical formulation for unsaturated rock. Darcy’s law for unsaturated flow has been derived from the dissipation process by using standard arguments of non-equilibrium thermodynamics, whereas Helmholtz free energy has been used to derive the relationship between stress and pore pressure changes. The resulting general framework accommodates both large and small deformation theories. When small deformations are assumed, the formulation is comparable with coupled equations derived using an alternative approach. For illustrative purposes, the formulation has been used to analyse a seasonally affected tunnel model. Numerical results for the desaturation in winter and resaturation in summer, of the zone near the tunnel wall, have been evaluated and compared with the findings of other researchers. 相似文献
2.
This paper examines difficulties encountered when using the Green-Nagdi rate and presents an updated configuration formulation for problems dealing with large deformations. It is shown that a natural link can be established between the updated reference configuration and linearized updated Lagrangian formulations, if Truesdell's objective stress rate is used. Owing to this connection, any limitations associated with the described formulation would be expected to also apply to an updated Lagrangian formulation and vice versa. It is pointed out that, if second-order effects are induced, the Truesdell rate appears in the rate equation for equilibrium even when strains are infinitesimal. By examining a simple-shear problem, where hypoelastic material behaviour with ‘Hooke's law’ is assumed, it is shown that the spherical state of stress is coupled with the shear stress when using the Truesdell or Green-Nagdi stress rate. It is suggested that objective stress rates, introduced in large-deformation formulations, may contribute to erroneous predictions if a constitutive law is developed independent of the adopted stress rate. 相似文献
3.
Parameterization of sea surface drag under varying sea state and its dependence on wave age 总被引：1，自引：0，他引：1
R. Rajesh Kumar B. Prasad Kumar A. N. V. Satyanarayana D. Bala Subrahamanyam A. D. Rao S. K. Dube 《Natural Hazards》2009,49(2):187-197
Momentum and energy exchange at air–sea interface through wind stress is very important for air–sea interaction studies, ocean
modeling, and climate studies. The accurate representation of wind stress, in terms of drag coefficient, is a key factor in
estimating the momentum transfer at the interface. The drag coefficient, in general, estimated using bulk formulae does not
take into account the influence of wave age. This study examines the dependence of wave age on computed surface drag coefficient
obtained by combining the Toba 3/2-power law with Froude number scaling, resulting in a new drag formulation (hereafter referred
as RP formulation). We demonstrate that our proposed formulation is in good conjunction with established theories for both
young and mature waves. Our investigation shows the theoretical formulation advocated earlier by Guan and Xie (hereafter referred
as GX) overestimated the surface drag for mature waves as wind speed tends to increase. In addition, the formulation by GX
was not verified by observational data. In the present work, for validation purpose, we use time series measurement of meteorological
and oceanographic data from a deep water location in the Indian Ocean which was tested with both RP and GX formulations. We
find that the proposed RP formulation, which embeds the 3/2 power of wave-age, shows a better match for both young and mature
waves with the results of Janssen compared to the hypothesis of conventional wave age used by GX. 相似文献
4.
This paper sets forth the theoretical background and basic numerical expressions for the incorporation of elastic-plastic constitutive equations for ductile rock into a finite element computer code. The derivation of an expression for the total strain rate is performed both for a total stress formulation and for a formulation that employs the concept of effective stress for inelastic behaviour. Specific expressions for the incremental strain rate are presented for the case of a porous material having a quadratic initial yield surface and observing the associated flow rule with a special hardening law for subsequent plastic deformation. A final section of the paper summarizes the expressions required to insert the quadratic yield surface model into a finite element code. 相似文献
5.
A new constitutive formulation for simulating the behaviour of nearly saturated sands under seismic loads is presented. The formulation is based on combining the Henry's law for dissolution of gas in water, the ideal or perfect gas law and the law of conservation of mass. The effects of transient air dissolution in water on the compressibility of partially saturated soils are also taken into account. The model was calibrated based on numerical simulations of isotropically consolidated cyclic triaxial tests conducted on partially saturated samples of Toyoura sand. A multi‐yield plasticity soil constitutive model implemented in the finite element code DYNAFLOW was used for these numerical simulations. It is shown that the formulation proposed here is able to reasonably predict the soil cyclic undrained behaviour at various degrees of saturation (95% and higher). Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
6.
This paper presents a novel formulation for defining soil failure. It plots in the principal stress space as a surface with the shape ranging between an approximation of the Matsuoka–Nakai and of the Mohr–Coulomb criteria depending on the value of a single parameter. The new function can be used as a replacement of the original equations of these well‐established criteria for implementing in a program for numerical analyses, and it is particularly effective for approximating the Matsuoka–Nakai criterion. Both the Mohr–Coulomb and the Matsuoka–Nakai failure criteria present numerical difficulties during implementation and also at run‐time. In the case of the Matsuoka–Nakai, the new formulation plots in the first octant only, whereas the original criterion plots in all octants, which causes severe convergence problems particularly for those Gauss points with low stress state, such as those on the side of a shallow footing. When the shape parameter is set to reproduce the Mohr–Coulomb failure criterion, on the other hand, the new formulation plots as a pyramid with rounded edges. Moreover, as the new function is at least of class C2, the second derivatives are continuous, thus ensuring quadratic convergence of the Newton's method used within the integration scheme of the constitutive law. The proposed formulation can also provide both sharp and rounded apex of the surface at the origin of the stress space by setting accordingly one additional parameter. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
7.
Coupled hydromechanical‐fracture simulations of nonplanar three‐dimensional hydraulic fracture propagation 
下载免费PDF全文
下载免费PDF全文 This paper presents an algorithm and a fully coupled hydromechanical‐fracture formulation for the simulation of three‐dimensional nonplanar hydraulic fracture propagation. The propagation algorithm automatically estimates the magnitude of time steps such that a regularized form of Irwin's criterion is satisfied along the predicted 3‐D fracture front at every fracture propagation step. A generalized finite element method is used for the discretization of elasticity equations governing the deformation of the rock, and a finite element method is adopted for the solution of the fluid flow equation on the basis of Poiseuille's cubic law. Adaptive mesh refinement is used for discretization error control, leading to significantly fewer degrees of freedom than available nonadaptive methods. An efficient computational scheme to handle nonlinear time‐dependent problems with adaptive mesh refinement is presented. Explicit fracture surface representations are used to avoid mapping of 3‐D solutions between generalized finite element method meshes. Examples demonstrating the accuracy, robustness, and computational efficiency of the proposed formulation, regularized Irwin's criterion, and propagation algorithm are presented. 相似文献
8.
A five-step procedure involving mathematical formulation, identification and determination of parameters and verification is presented for development and selection of appropriate and reliable constitutive law(s) for geologic media. Comprehensive analyses are performed toward determination of an appropriate law for a (artificial) soil. The most suitable model is obtained by critical evaluation of four different plasticity models; here verification and comparisons of predictions with observations from laboratory tests, and with those from two boundary value problems are used as the basis of the selection. The model thus selected is found to be appropriate for applications to relevant practical problems. 相似文献
9.
A theoretical formulation and a numerical solution method are proposed for the problem of the time dependent consolidation of an elasto-plastic soil subject to finite deformations. The soil is assumed to be a two-phase material with a skeleton which may yield according to a general yield criterion with plastic flow governed by a general flow law, and whose pore fluid flows according to Darcy's Law. Governing equations are cast in a rate form and constitutive laws are expressed in a frame indifferent manner. The method of analysis is illustrated by several examples of practical interest for both a soil with an elastic skeleton and a soil with an elasto-plastic skeleton which obeys a Morh–Coulomb yield criterion and a non-associated flow law. 相似文献
10.
In this paper, the problem of propagation of localized deformation associated with formation of macrocracks/shear bands is studied in both tensile and compressive regimes. The main focus here is on enhancement of the constitutive law with embedded discontinuity to provide a discrete representation of the localization phenomenon. This has been accomplished by revising the formulation and coupling it with the level‐set method for tracing the propagation path. Extensive numerical studies are conducted involving various fracture modes, ranging from brittle to frictional, and the results are compared with the experimental data as well as those obtained using XFEM methodology. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
11.
This paper presents a new, fully-coupled, hydro-mechanical (HM) formulation for a finite-discrete element method computer code. In the newly-developed, hydraulic solver, fluid flow is assumed to occur through the same triangular mesh used for the mechanical calculations. The flow of a viscous, compressible fluid is explicitly solved based on a cubic law approximation. The implementation is verified against closed-form solutions for several flow problems. The approach is then applied to a field-scale simulation of fluid injection in a jointed, porous rock mass. Results show that the proposed method can be used to obtain unique geomechanical insights into coupled HM phenomena. 相似文献
12.
A general thermo-hydro-mechanical framework for the modelling of internal erosion is proposed based on the theory of mixtures applied to two-phase porous media. The erodible soil is partitioned in two phases: one solid phase and one fluid phase. The solid phase is composed of nonerodible grains and erodible particles. The fluid phase is composed of water and fluidized particles. Within the fluid phase, species diffuse. Across phases, species transfer. The modelling of internal erosion is contributed directly by mass transfer from the solid phase towards the fluid phase. The constitutive relations governing the thermomechanical behaviour, generalised diffusion, and transfer are structured by the dissipation inequality. The particular case of soil suffusion is investigated with a focus on constitutive laws. A new constitutive law for suffusion is constructed based on thermodynamic conditions and experimental investigations. This erosion law is linearly related to the power of seepage flow and to the erosion resistance index. Owing to its simplicity, this law tackles the overall trend of the suffusion process and permits the formulation of an analytical solution. This new model is then applied to simulate laboratory experiments, by both analytical and numerical methods. The comparison shows that the newly developed model, which is theoretically consistent, can reproduce correctly the overall trend of the cumulated eroded mass when the permeability evolution is small. In addition, the results are provided for four different materials, two different specimen sizes, and various hydraulic loading paths to demonstrate the applicability of the new proposed law. 相似文献
13.
14.
V. Navarro A. Yustres M. Candel J. López E. Castillo 《Computers and Geotechnics》2010,37(7-8):837-845
This article discusses how sensitivity analysis is a sound assessment tool for selecting the most efficient stabilization method of slopes at failure. A discretized form of the variational approach is used not only for performing sensitivity analysis but to locate the critical slip surface, i.e., the sensitivity analysis is carried out in the same way as it is done in optimization problems. This method supplies a robust formulation and methodology for obtaining the sensitivities of the safety factor with respect to both the soil parameters and the slope profile, stating the slope stabilization design as a relatively simple minimization problem. Two well known examples, as the Selset landslide and the Sudbury Hill slip are used to illustrate the application of the method and to highlight both its capabilities and limitations. 相似文献
15.
In this paper, a simple bounding surface plasticity model is used to reproduce the yielding and stress–strain behavior of the structured soft clay found at Shanghai of China. A series of undrained triaxial tests and drained stress probe tests under isotropic and anisotropic consolidation modes were performed on undisturbed samples of Shanghai soft clay to study the yielding characteristics. The degradation of the clay structure is modeled with an internal variable that allows the size of the bounding surface to decay with accumulated plastic strain. An anisotropic tensor and rotational hardening law are introduced to reflect the initial anisotropy and the evolution of anisotropy. Combined with the isotropic hardening rule, the rotational hardening rule and the degradation law are incorporated into the bounding surface formulation with an associated flow rule. Validity of the model is verified by the undrained isotropic and anisotropic triaxial test and drained stress probe test results for Shanghai soft clay. The effects of stress anisotropy and loss of structure are well captured by the model. 相似文献
16.
Cumulative Benioff strain-release, modified Omori's law and transient behaviour of rocks 总被引：1，自引：1，他引：1
Irreversible thermodynamic theories with internal state variables can be used to derive a general constitutive law for both transient and steady-state behaviours of rocks. This constitutive law can represent the concepts of damage and damage evolution in either the fibre-bundle model or continuum damage mechanics. We have previously proposed an empirically based constitutive law for both the transient and steady-state behaviours of rocks ultimately derived from laboratory experimental data. We show here that this law is concordant with the general constitutive law derived from irreversible thermodynamic theories, and that the relaxation modulus has a temporal power–law that depends on a structural fractal property of rocks. Our constitutive law predicts forms for the cumulative Benioff strain-release for precursory seismic activations and the modified Omori's laws of aftershocks, both aspects of the temporal fractal properties of seismicity. These seismic properties can also be derived by the fibre-bundle model or continuum damage mechanics. Our model suggests that these time-scale invariant processes of seismicity may be regulated by the fractal structures of crustal rocks. 相似文献
17.
18.
We present a stabilized extended finite element formulation to simulate the hydraulic fracturing process in an elasto‐plastic medium. The fracture propagation process is governed by a cohesive fracture model, where a trilinear traction‐separation law is used to describe normal contact, cohesion and strength softening on the fracture face. Fluid flow inside the fracture channel is governed by the lubrication equation, and the flow rate is related to the fluid pressure gradient by the ‘cubic’ law. Fluid leak off happens only in the normal direction and is assumed to be governed by the Carter's leak‐off model. We propose a ‘local’ U‐P (displacement‐pressure) formulation to discretize the fluid‐solid coupled system, where volume shape functions are used to interpolate the fluid pressure field on the fracture face. The ‘local’ U‐P approach is compatible with the extended finite element framework, and a separate mesh is not required to describe the fluid flow. The coupled system of equations is solved iteratively by the standard Newton‐Raphson method. We identify instability issues associated with the fluid flow inside the fracture channel, and use the polynomial pressure projection method to reduce the pressure oscillations resulting from the instability. Numerical examples demonstrate that the proposed framework is effective in modeling 3D hydraulic fracture propagation. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
19.
20.
In this paper, a coupled constitutive model is proposed for anisotropic damage and permeability variation in brittle rocks under deviatoric compressive stresses. The formulation of the model is based on experimental evidences and main physical mechanisms involved in the scale of microcracks are taken into account. The proposed model is expressed in the macroscopic framework and can be easily implemented for engineering application. The macroscopic free enthalpy of cracked solid is first determined by approximating crack distribution by a second‐order damage tensor. The effective elastic properties of damaged material are then derived from the free enthalpy function. The damage evolution is related to the crack growth in multiple orientations. A pragmatic approach inspired from fracture mechanics is used for the formulation of the crack propagation criterion. Compressive stress induced crack opening is taken into account and leads to macroscopic volumetric dilatancy and permeability variation. The overall permeability tensor of cracked material is determined using a micro–macro averaging procedure. Darcy's law is used for fluid flow at the macroscopic scale whereas laminar flow is assumed at the microcrack scale. Hydraulic connectivity of cracks increases with crack growth. The proposed model is applied to the Lac du Bonnet granite. Generally, good agreement is observed between numerical simulations and experimental data. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献