共查询到20条相似文献，搜索用时 31 毫秒
1.
This paper presents an analytical solution for cavity expansion in thermoplastic soil considering non‐isothermal conditions. The constitutive relationship of thermoplasticity is described by Laloui's advanced and unified constitutive model for environmental geomechanical thermal effect (ACMEG‐T), which is based on multi‐mechanism plasticity and bounding surface theory. The problem is formulated by incorporating ACMEG‐T into the theoretical framework of cavity expansion, yielding a series of partial differential equations (PDEs). Subsequently, the PDEs are transformed into a system of first‐order ordinary differential equations (ODEs) using a similarity solution technique. Solutions to the response parameters of cavity expansion (stress, excess pore pressure, and displacement) can then be obtained by solving the ODEs numerically using mathematical software. The results suggest that soil temperature has a significant influence on the pressure‐expansion relationships and distributions of stress and excess pore pressure around the cavity wall. The proposed solution quantifies the influence of temperature on cavity expansion for the first time and provides a theoretical framework for predicting thermoplastic soil behavior around the cavity wall. The solution found in this paper can be used as a theoretical tool that can potentially be employed in geotechnical engineering problems, such as thermal cone penetration tests, and nuclear waste disposal problems. 相似文献
2.
A novel procedure associated with the precise integration method (PIM) and the technique of dual vector is proposed to effectively calculate the magnitude and distribution of deformations in a homogeneous multilayered transversely isotropic medium. The planes of transverse isotropy are assumed to be parallel to the horizontal surface of the soil system. The linearly elastic medium is subjected to four types of vertically acting axisymmetric loads prescribed either at the external surface or in the interior of the soil medium. There are no limits for the thicknesses and number of soil layers to be considered. By virtue of the governing equations of motion and the constitutive equations of the transversely isotropic elastic body, and based on the Hankel integral transform and a dual vector formulation in a cylindrical coordinate system, the partial differential motion equations can be converted into first‐order ordinary differential matrix equations. Applying the approach of PIM, it is convenient to obtain the solutions of ordinary differential matrix equations for the continuously homogeneous multilayered transversely isotropic elastic soil in the transformed domain. The PIM is a highly accurate algorithm to solve the sets of first‐order ordinary differential equations, which can ensure to achieve any desired accuracy of the solutions. What is more, all calculations are based on the standard method with the corresponding algebraic operations. Computational efforts can be reduced to a great extent. Finally, numerical examples are provided to illustrate the accuracy and effectiveness of the proposed approach. Some more cases are analyzed to evaluate the influences of the elastic parameters of the transversely isotropic media on the load‐displacement responses. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
3.
In this article, an approach for the efficient numerical solution of multi-species reactive transport problems in porous media
is described. The objective of this approach is to reformulate the given system of partial and ordinary differential equations
(PDEs, ODEs) and algebraic equations (AEs), describing local equilibrium, in such a way that the couplings and nonlinearities
are concentrated in a rather small number of equations, leading to the decoupling of some linear partial differential equations
from the nonlinear system. Thus, the system is handled in the spirit of a global implicit approach (one step method) avoiding
operator splitting techniques, solved by Newton’s method as the basic algorithmic ingredient. The reduction of the problem
size helps to limit the large computational costs of numerical simulations of such problems. If the model contains equilibrium
precipitation-dissolution reactions of minerals, then these are considered as complementarity conditions and rewritten as
semismooth equations, and the whole nonlinear system is solved by the semismooth Newton method. 相似文献
4.
A method for the analysis of the consolidation of a horizontally layered soil under plane conditions is developed. The method depends upon the transformation of the governing equations by a Fourier trasform. This transformation has the effect of reducing the partial differential equations of consolidation to ordinary differential equations. The ordinary differential equations are then solved using a finite layer or finite difference approach. Once the solution in the transformed plane has been found, the actual solution is synthesized by Fourier inversion. The method leads to a considerable reduction in the amount of core storage necessary for solution and enables the solution of quite significant problems to be obtained on a mini-computer. 相似文献
5.
E. Alonso L. R. Alejano F. Varas G. Fdez‐Manin C. Carranza‐Torres 《国际地质力学数值与分析法杂志》2003,27(13):1153-1185
A literature review has shown that there exist adequate techniques to obtain ground reaction curves for tunnels excavated in elastic‐brittle and perfectly plastic materials. However, for strain‐softening materials it seems that the problem has not been sufficiently analysed. In this paper, a one‐dimensional numerical solution to obtain the ground reaction curve (GRC) for circular tunnels excavated in strain‐softening materials is presented. The problem is formulated in a very general form and leads to a system of ordinary differential equations. By adequately defining a fictitious ‘time’ variable and re‐scaling some variables the problem is converted into an initial value one, which can be solved numerically by a Runge–Kutta–Fehlberg method, which is implemented in MATLAB environment. The method has been developed for various common particular behaviour models including Tresca, Mohr–Coulomb and Hoek–Brown failure criteria, in all cases with non‐associative flow rules and two‐segment piecewise linear functions related to a principal strain‐dependent plastic parameter to model the transition between peak and residual failure criteria. Some particular examples for the different failure criteria have been run, which agree well with closed‐form solutions—if existing—or with FDM‐based code results. Parametric studies and specific charts are created to highlight the influence of different parameters. The proposed methodology intends to be a wider and general numerical basis where standard and newly featured behaviour modes focusing on obtaining GRC for tunnels excavated in strain‐softening materials can be implemented. This way of solving such problems has proved to be more efficient and less time consuming than using FEM‐ or FDM‐based numerical 2D codes. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
6.
7.
This paper presents a complete finite‐element treatment for unsaturated soil problems. A new formulation of general constitutive equations for unsaturated soils is first presented. In the incremental stress–strain equations, the suction or the pore water pressure is treated as a strain variable instead of a stress variable. The global governing equations are derived in terms of displacement and pore water pressure. The discretized governing equations are then solved using an adaptive time‐stepping scheme which automatically adjusts the time‐step size so that the integration error in the displacements and pore pressures lies close to a specified tolerance. The non‐linearity caused by suction‐dependent plastic yielding, suction‐dependent degree of saturation, and saturation‐dependent permeability is treated in a similar way to the elastoplasticity. An explicit stress integration scheme is used to solve the constitutive stress–strain equations at the Gauss point level. The elastoplastic stiffness matrix in the Euler solution is evaluated using the suction as well as the stresses and hardening parameters at the start of the subincrement, while the elastoplastic matrix in the modified Euler solution is evaluated using the suction at the end of the subincrement. In addition, when applying subincrementation, the same rate is applied to all strain components including the suction. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
8.
An alternative approach for quasi‐static large deformation analysis of saturated porous media using meshfree method 
下载免费PDF全文
下载免费PDF全文 An alternative coupled large deformation formulation combined with a meshfree approach is proposed for flow–deformation analysis of saturated porous media. The formulation proposed is based on the Updated Lagrangian (UL) approach, except that the spatial derivatives are defined with respect to the configuration of the medium at the last time step rather than the configuration at the last iteration. In this way, the Cauchy stresses are calculated directly, rendering the second Piola–Kirchhoff stress tensor not necessary for the numerical solution of the equilibrium equations. Moreover, in contrast with the UL approach, the nodal shape function derivatives are calculated once in each time step and stored for use in subsequent iterations, which reduces the computational cost of the algorithm. Stress objectivity is satisfied using the Jaumann stress rate, and the spatial discretisation of the governing equations is achieved using the standard Galerkin method. The equations of equilibrium are satisfied directly, and the nonlinear parts of the system matrix are derived independent of the stresses of the medium resulting in a stable numerical algorithm. Temporal discretisation is effected based on a three‐point approximation technique that avoids spurious ripple effects and has second‐order accuracy. The radial point interpolation method is used to construct the shape functions. The application of the formulation and the significance of large deformation effects on the numerical results are demonstrated through several numerical examples. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
9.
In order to capture the influence of the cavity expansion velocity, this paper presents a semianalytical solution for dynamic spherical cavity expansion in modified Cam Clay (MCC) soil. The key problem is solving the six coupled partial differential equations (PDEs) of cavity expansion, in which the dynamic term is considered in the stress equilibrium equation. The similarity transformation technique is used to transform the PDEs into ordinary differential equations (ODEs). Subsequently, the numerical method using the function “ODE45” in MATLAB is selected to solve the ODEs, which allows the stress and excess pore pressure around the expanding spherical cavity wall to be obtained. The proposed semianalytical solution for dynamic spherical cavity expansion was validated by comparting the degenerate solution with the published quasistatic solution for the MCC model. Parametric study was then conducted to capture the influence of the cavity wall velocity on the cavity expansion response. The proposed solution has potential application to geotechnical problems such as dynamic pile driving, the dynamic cone penetration test, and so forth. 相似文献
10.
3D dynamic response of a transversely isotropic multilayered medium subjected to a moving load 下载免费PDF全文
下载免费PDF全文 The dynamic problem of a transversely isotropic multilayered medium is reducible to quasi‐static problem by introducing a moving system that travels synchronously with the load. Based on the governing equations in the moving system, the ordinary differential equations in the Fourier transformed domain are deduced. An extended precise integration method is adopted to solve the ordinary differential equations, and the solution in the physical domain is recovered by the inverse Fourier transform. Numerical examples are performed to verify the accuracy of the presented method and to analyze the influence of material properties and the load characteristic. 相似文献
11.
A large strain analysis of undrained expansion of a spherical/cylindrical cavity in a soil modelled as non‐linear elastic modified Cam clay material is presented. The stress–strain response of the soil is assumed to obey non‐linear elasticity until yielding. A power‐law characteristic or a hyperbolic stress–strain curve is used to describe the gradual reduction of soil stiffness with shear strain. It is assumed that, after yielding, the elasto‐plastic behaviour of the soil can be described by the modified Cam clay model. Based on a closed‐form stress–strain response in undrained condition, a numerical solution is obtained with the aid of simple numerical integration technique. The results show that the stresses and the pore pressure in the soil around an expanded cavity are significantly affected by the non‐linear elasticity, especially if the soil is overconsolidated. The difference between large strain and small strain solutions in the elastic zone is not significant. The stresses and the pore pressure at the cavity wall can be expressed as an approximate closed‐form solution. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
12.
The Crank–Nicolson scheme has second‐order accuracy, but often leads to oscillations affecting numerical stability. On the other hand, the implicit scheme is free from oscillation, but it has only first‐order accuracy. In this work, a three‐point discretization scheme with variable time step is presented for the time marching of parabolic partial differential equations. The method proposed has second‐order accuracy, is unconditionally stable and dampens spurious oscillations of the numerical results. The application and effectiveness of the new method are demonstrated through several numerical examples. It is shown that, unlike the Crank–Nicolson method, the approach proposed produces no oscillatory response irrespective of the time step adopted. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
13.
Extended precise integration method for axisymmetric thermo‐elastic problem in transversely isotropic material 下载免费PDF全文
下载免费PDF全文 This paper presents a numerical solution for the analysis of the axisymmetric thermo‐elastic problem in transversely isotropic material due to a buried heat source by means of extended precise integral method. By virtue of the Laplace–Hankel transform applied into the basic governing equations, an ordinary differential matrix equation is achieved, which describes the relationship between the generalized stresses and displacements in transformed domain. An extended precise integration method is introduced to solve the aforementioned matrix equation, and the actual solution in the physical domain is acquired by inverting the Laplace–Hankel transform. Numerical examples are carried out to demonstrate the accuracy of the proposed method and elucidate the influence of the character of transverse isotropy, the anisotropy of linear expansion coefficient, the anisotropy of thermal diffusivity, and medium's stratification on the thermo‐elastic response. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
14.
The formulation of macroscopic poroelastic behavior of a jointed rock is investigated within the framework of a micro–macro approach. The joints are modeled as interfaces, and their behavior is modeled by means of generalized poroelastic state equations. Starting from Hill's lemma extended for a jointed medium and extending the concept of strain concentration to relate the joint displacement jump to macroscopic strain, the overall poroelastic constitutive equations for the jointed rock are formulated. The analysis emphasizes the main differences and similarities of the resulting behavior with respect to that characterizing ordinary porous media. It is shown that, unlike ordinary porous media, conditions on the poroelastic parameters of joints are required for the macroscopic drained stiffness to entirely define the poroelastic behavior. This is achieved, for instance, if the joint network is characterized by a unique Biot coefficient. Extension of the analysis to non‐linear poroelasticity is also outlined. Finally, the theoretical formulation is applied to two particular cases of jointed rock for which explicit expressions of the overall poroelastic parameters are derived. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
15.
A numericl method for solving consolidation problems of layered soils is developed. Starting from the governing differential equations for the coupled poro-elastic medium, the governing partial differential equations are reduced to ordinary differential equations by means of the appropriate displacement functions and Laplace-Fourier transformation. Once the fundamental solution in the transformed domain has been found, the solution in the physical domain is obtained by numerically inverting the transformations. A series of soil consolidation problems have been solved and validated against existing solutions in order to compare the feasibility and the accuracy of the present technique. 相似文献
16.
Based on the Fredlund consolidation theory of unsaturated soil, exact solutions of the governing equations for one‐dimensional consolidation of single‐layer unsaturated soil are presented, in which the water permeability and air transmission are assumed to be constants. The general solution of two coupled homogeneous governing equations is first obtained. This general solution is expressed in terms of two functions psi1 and ψ2, where ψ1 and ψ2, respectively, satisfy two second‐order partial differential equations, which are in the same form. Using the method of separation of variables, the two partial differential equations are solved and exact solutions for three typical homogeneous boundary conditions are obtained. To obtain exact solutions of nonhomogeneous governing equations with three typical nonhomogeneous boundary conditions, the nonhomogeneous boundary conditions are first transformed into homogeneous boundary conditions. Then according to the method of undetermined coefficients and exact solutions of homogenous governing equations, the series form exact solutions are put forward. The validity of the proposed exact solutions is verified against other analytical solutions in the literature. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
17.
Mathematical simulation of non‐isothermal multiphase flow in deformable unsaturated porous media is a complicated issue because of the need to employ multiple partial differential equations, the need to take into account mass and energy transfer between phases and because of the non‐linear nature of the governing partial differential equations. In this paper, an analytical solution for analyzing a fully coupled problem is presented for the one‐dimensional case where the coefficients of the system of equations are assumed to be constant for the entire domain. A major issue is the non‐linearity of the governing equations, which is not considered in the analytical solution. In order to introduce the non‐linearity of the equations, an iterative discretized procedure is used. The domain of the problem is divided into identical time–space elements that cover the time–space domain. A separate system of equations is defined for each element in the local coordinate system, the initial and boundary conditions for each element are obtained from the adjacent elements and the coefficients of the system of equations are considered to be constant in each step. There are seven governing differential equations that should be solved simultaneously: the equilibrium of the solid skeleton, mass conservation of fluids (water, water vapor and gas) and energy conservation of phases (solid, liquid and gas). The water vapor is not in equilibrium with water and different phases do not have the same temperature. The governing equations that have been solved seem to be the most comprehensive in this field. Three examples are presented for analyzing heat and mass transfer in a semi‐infinite column of unsaturated soil. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
18.
Grégoire Allaire Robert Brizzi Jean-François Dufrêche Andro Mikelić Andrey Piatnitski 《Computational Geosciences》2013,17(3):479-495
In this work, we undertake a numerical study of the effective coefficients arising in the upscaling of a system of partial differential equations describing transport of a dilute N-component electrolyte in a Newtonian solvent through a rigid porous medium. The motion is governed by a small static electric field and a small hydrodynamic force, around a nonlinear Poisson–Boltzmann equilibrium with given surface charges of arbitrary size. This approach allows us to calculate the linear response regime in a way initially proposed by O’Brien. The O’Brien linearization requires a fast and accurate solution of the underlying Poisson–Boltzmann equation. We present an analysis of it, with the discussion of the boundary layer appearing as the Debye–Hückel parameter becomes large. Next, we briefly discuss the corresponding two-scale asymptotic expansion and reduce the obtained two-scale equations to a coarse scale model. Our previous rigorous study proves that the homogenized coefficients satisfy Onsager properties, namely they are symmetric positive definite tensors. We illustrate with numerical simulations several characteristic situations and discuss the behavior of the effective coefficients when the Debye–Hückel parameter is large. Simulated qualitative behavior differs significantly from the situation when the surface potential is given (instead of the surface charges). In particular, we observe the Donnan effect (exclusion of co-ions for small pores). 相似文献
19.
Most analytical or semi‐analytical solutions of the problem of load‐settlement response of axially loaded piles are based on the assumption of zero radial displacement. These solutions also are only applicable to piles embedded in either a homogeneous or a Gibson soil deposit. In reality, soil deposits consist of multiple soil layers with different properties, and displacements in the radial direction within the soil deposit are not zero when the pile is loaded axially. In this paper, we present a load‐settlement analysis applicable to a pile with circular cross section installed in multilayered elastic soil that accounts for both vertical and radial soil displacements. The analysis follows from the solution of the differential equations governing the displacements of the pile–soil system obtained using variational principles. The input parameters needed for the analysis are the pile geometry and the elastic constants of the soil and pile. We compare the results from the present analysis with those of an analytical solution that considers only vertical soil displacements. The analysis presented in this paper also provides useful insights into the displacement and strain fields around axially loaded piles. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
20.
A new approach on numerical modeling of wave propagation is introduced and is used to analyze the effect of earthquake magnitudes (ground motion amplitudes) on wave propagation. In this method, the sum of the maximum amplitudes of the first output model at time 0 s and rest of the output models at different times are normalized to unity. Considering this as a constraint, the sum of the weighted‐squared Fourier amplitudes is minimized by using the Lagrange multiplier method. The proposed method can reveal the relationship of actual time histories by showing simple clear peaks. This method is used to analyze the time histories of various earthquake events at different vertical array sites of the Kashiwazaki–Kariwa nuclear power plant of Tokyo electric power company (TEPCO). The wave arrival times obtained from this method and down‐hole measurements are compared. The results show increase in the arrival times at surface layer when the magnitude of earthquake is large. The results reveal that the amplitudes of small magnitude earthquakes at depths are small and are largely amplified at surface, whereas in case of large magnitude earthquakes, the amplitudes are large at depths and are deamplified at surface reflecting the effects of the strain‐dependent soil properties that result in non‐linear site response to strong shaking. The results also show that the reflected peak amplitudes are higher for small magnitude earthquakes than for large magnitude earthquakes. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献