共查询到20条相似文献，搜索用时 750 毫秒
1.
Yuki Yamakawa Kiyohiro Ikeda Isao Saiki Jacques Desrues Reiko J. Tanaka 《国际地质力学数值与分析法杂志》2018,42(1):3-33
Shear bands with characteristic spatial patterns observed in an experiment for a cubic or parallelepiped specimen of dry dense sand were simulated by numerical bifurcation analysis using the Cam‐clay plasticity model. By incorporating the subloading surface concept into the plasticity model, the model became capable of reproducing hardening/softening and contractive/dilative behavior observed in the experiment. The model was reformulated to be compatible with the multiplicative hyperelasto‐plasticity for finite strains. This enhanced constitutive model was implemented into a finite‐element code reinforced by a stress updating algorithm based on the return‐mapping scheme, and by an efficient numerical procedure to compute critical eigenvectors of elastoplastic tangent stiffness matrix at bifurcation points. The emergence of diamond‐ and column‐like diffuse bifurcation modes breaking uniformity of the materials, followed by the evolution of shear bands through strain localization, was observed in the analysis. In the bifurcation analysis of plane strain compression test, unexpected bifurcation modes, which broke out‐of‐plane uniformity and led to 3‐dimensional diamond‐like patterns, were detected. Diffuse bifurcations, which were difficult to observe by experiments, have thus been found as a catalyst creating diverse shear band patterns. 相似文献
2.
This paper examines the performance of the Jacobi preconditioner when used with two Krylov subspace iterative methods. The number of iterations needed for convergence was shown to be different for drained, undrained and consolidation problems, even for similar condition number. The differences were due to differences in the eigenvalue distribution, which cannot be completely described by the condition number alone. For drained problems involving large stiffness ratios between different material zones, ill‐conditioning is caused by these large stiffness ratios. Since Jacobi preconditioning operates on degrees‐of‐freedom, it effectively homogenizes the different spatial sub‐domains. The undrained problem, modelled as a nearly incompressible problem, is much more resistant to Jacobi preconditioning, because its ill‐conditioning arises from the large stiffness ratios between volumetric and distortional deformational modes, many of which involve the similar spatial domains or sub‐domains. The consolidation problem has two sets of degrees‐of‐freedom, namely displacement and pore pressure. Some of the eigenvalues are displacement dominated whereas others are excess pore pressure dominated. Jacobi preconditioning compresses the displacement‐dominated eigenvalues in a similar manner as the drained problem, but pore‐pressure‐dominated eigenvalues are often over‐scaled. Convergence can be accelerated if this over‐scaling is recognized and corrected for. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
3.
It is shown that the property of the scale invariance of the eigenvalues and eigenmodes of a finite element can be used as a basis to calculate good approximations to the analytical magnitudes of eigenvalues. This requires the subdivision of the element into a mesh of small elements with the same shape as the large element, the enforcement of the modal boundary displacements of the large element to the mesh of small elements and finally the application of the conditions of both the nodal equilibrium and the equality of the nodal work at both scales. Due to the self‐similarity of the elements at all scales the authors propose to call this method the fractal approach. The method is applied to calculate the hour‐glass eigenvalue of a plane square 4‐node quad for isotropic linear elastic material. The resulting hour‐glass eigenvalue is shown to be a good approximation of the analytical magnitude as derived in a companion paper. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
4.
Stress and strain fields for overdamped slug tests in aquifer materials,and resulting conservation equation 
下载免费PDF全文
下载免费PDF全文 Robert P. Chapuis 《国际地质力学数值与分析法杂志》2017,41(18):1908-1921
Two theories may be used to analyze overdamped variable‐head (slug) tests in aquifer materials. The first theory assumes that the solid matrix strain has a negligible influence. The second theory takes into account some elastic and immediate strain. Something is wrong with these theories because they yield different hydraulic conductivity values. This paper explains what is wrong after establishing the strain‐stress elastic equations for a slug test in 2 ideal conditions: plane strain and spherical symmetry. The equations show that the radial contraction (or elongation) and tangential elongation (or contraction) yield a null volumetric strain. As a result, the conservation is described by the Laplace equation, which is used by the first theory. This first theory is the only one to yield correct solutions. The diffusion equation, with storativity, which is used by the second theory, is physically ill founded for slug tests in aquifers. This new proof scientifically confirms previously raised doubts and experimental proofs that the second theory is ill founded. 相似文献
5.
The present paper investigates bifurcation analysis based on the second‐order work criterion, in the framework of rate‐independent constitutive models and rate‐independent boundary‐value problems. The approach applies mainly to nonassociated materials such as soils, rocks, and concretes. The bifurcation analysis usually performed at the material point level is extended to quasi‐static boundary‐value problems, by considering the stiffness matrix arising from finite element discretization. Lyapunov's definition of stability (Annales de la faculté des sciences de Toulouse 1907; 9 :203–274), as well as definitions of bifurcation criteria (Rice's localization criterion (Theoretical and Applied Mechanics. Fourteenth IUTAM Congress, Amsterdam, 1976; 207–220) and the plasticity limit criterion are revived in order to clarify the application field of the second‐order work criterion and to contrast these criteria. The first part of this paper analyses the second‐order work criterion at the material point level. The bifurcation domain is presented in the 3D stress space as well as 3D cones of unstable loading directions for an incrementally nonlinear constitutive model. The relevance of this criterion, when the nonlinear constitutive model is expressed in the classical form (dσ = Mdε) or in the dual form (dε = Ndσ), is discussed. In the second part, the analysis is extended to the boundary‐value problems in quasi‐static conditions. Nonlinear finite element computations are performed and the global tangent stiffness matrix is analyzed. For several examples, the eigenvector associated with the first vanishing eigenvalue of the symmetrical part of the stiffness matrix gives an accurate estimation of the failure mode in the homogeneous and nonhomogeneous boundary‐value problem. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
6.
Softening solids are analysed under impact loading using a new numerical method which allows displacement discontinuities to propagate arbitrarily through a finite element mesh. The Dirac‐delta distributions that arise in the strain field of classical continuum theory in the presence of strain softening are interpreted as discontinuities in the displacement field. A new finite element procedure with Heaviside jumps added to the underlying displacement interpolation basis is able to capture displacement jumps independent of the spatial discretisation. The amplitudes of displacement jumps are represented by extra degrees of freedom at existing nodes. Numerical results for mode‐I and mode‐II failure due to impact loading are presented. The numerical results highlight the objectivity of the approach with respect to spatial discretisation under dynamic loading conditions. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
7.
This paper presents a novel mass-conservative mixed multiscale method for solving flow equations in heterogeneous porous media. The media properties (the permeability) contain multiple scales and high contrast. The proposed method solves the flow equation in a mixed formulation on a coarse grid by constructing multiscale basis functions. The resulting velocity field is mass-conservative on the fine grid. Our main goal is to obtain first-order convergence in terms of the mesh size which is independent of local contrast. This is achieved, first, by constructing some auxiliary spaces, which contain global information that cannot be localized, in general. This is built on our previous work on the generalized multiscale finite element method (GMsFEM). In the auxiliary space, multiscale basis functions corresponding to small (contrast-dependent) eigenvalues are selected. These basis functions represent the high-conductivity channels (which connect the boundaries of a coarse block). Next, we solve local problems to construct multiscale basis functions for the velocity field. These local problems are formulated in the oversampled domain, taking into account some constraints with respect to auxiliary spaces. The latter allows fast spatial decay of local solutions and, thus, allows taking smaller oversampled regions. The number of basis functions depends on small eigenvalues of the local spectral problems. Moreover, multiscale pressure basis functions are needed in constructing the velocity space. Our multiscale spaces have a minimal dimension, which is needed to avoid contrast dependence in the convergence. The method’s convergence requires an oversampling of several layers. We present an analysis of our approach. Our numerical results confirm that the convergence rate is first order with respect to the mesh size and independent of the contrast. 相似文献
8.
Scaling variables and interpretation of eigenvalues in principal component analysis of geologic data
A. T. Miesch 《Mathematical Geology》1980,12(6):523-538
The dominant feature distinguishing one method of principal components analysis from another is the manner in which the original data are transformed prior to the other computations. The only other distinguishing feature of any importance is whether the eigenvectors of the inner product-moment of the transformed data matrix are taken directly as the Q-mode scores or scaled by the square roots of their associated eigenvalues and called the R-mode loadings. If the eigenvectors are extracted from the product-moment correlation matrix, the variables, in effect, were transformed by column standardization (zero means and unit variances), and the sum of the p-largest eigenvalues divided by the sum of all the eigenvalues indicates the degree to which a model containing pcomponents will account for the total variance in the original data. However, if the data were transformed in any manner other than column standardization, the eigenvalues cannot be used in this manner, but can only be used to determine the degree to which the model will account for the transformed data. Regardless of the type of principal components analysis that is performed—even whether it is Ror Q-mode—the goodness-of-fit of the model to the original data is given better by the eigenvalues of the correlation matrix than by those of the matrix that was actually factored. 相似文献
9.
Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix Kep (ILU(0)‐Kep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐Kep. To further reduce computational cost, the global stiffness matrix Kep is divided into two parts. The first part is the linear elastic stiffness matrix Ke, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐Ke preconditioner is more time effective than IDR(s = 6) with ILU(0)‐Kep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
10.
11.
A challenging computational problem arises when a discrete structure (e.g. foundation) interacts with an unbounded medium (e.g. deep soil deposit), particularly if general loading conditions and non‐linear material behaviour is assumed. In this paper, a novel method for dealing with such a problem is formulated by combining conventional three‐dimensional finite‐elements with the recently developed scaled boundary finite‐element method. The scaled boundary finite‐element method is a semi‐analytical technique based on finite‐elements that obtains a symmetric stiffness matrix with respect to degrees of freedom on a discretized boundary. The method is particularly well suited to modelling unbounded domains as analytical solutions are found in a radial co‐ordinate direction, but, unlike the boundary‐element method, no complex fundamental solution is required. A technique for coupling the stiffness matrix of bounded three‐dimensional finite‐element domain with the stiffness matrix of the unbounded scaled boundary finite‐element domain, which uses a Fourier series to model the variation of displacement in the circumferential direction of the cylindrical co‐ordinate system, is described. The accuracy and computational efficiency of the new formulation is demonstrated through the linear elastic analysis of rigid circular and square footings. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
12.
Unbounded plane stress and plane strain domains subjected to static loading undergo infinite displacements, even when the zero displacement boundary condition at infinity is enforced. However, the stress and strain fields are well behaved, and are of practical interest. This causes significant difficulty when analysis is attempted using displacement‐based numerical methods, such as the finite‐element method. To circumvent this difficulty problems of this nature are often changed subtly before analysis to limit the displacements to finite values. Such a process is unsatisfactory, as it distorts the solution in some way, and may lead to a stiffness matrix that is nearly singular. In this paper, the semi‐analytical scaled boundary finite‐element method is extended to permit the analysis of such problems without requiring any modification of the problem itself. This is possible because the governing differential equations are solved analytically in the radial direction. The displacement solutions so obtained include an infinite component, but relative motion between any two points in the unbounded domain can be computed accurately. No small arbitrary constants are introduced, no arbitrary truncation of the domain is performed, and no ill‐conditioned matrices are inverted. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
13.
Effects of recoverable deformation induced anisotropy in the elastic stiffness of isotropic materials are described. In isotropic materials, thermodynamics predicts coupling of hydrostatic and deviatoric responses. It is shown that the coupling of the two responses is more significant than previously recognized in the literature. Properly accounting for the coupling of hydrostatic and deviatoric responses requires re‐evaluating elastic materials characterization data, allowing for the coupled response. The result is an apparent decrease in the pressure sensitivity of the elastic shear modulus. The decrease in the pressure sensitivity of the shear modulus leads to stress paths that are more tangential to the yield surface in stress space, resulting in an increase in predicted elastic strain at each step of an elastic–plastic stress update. Consequently, predicted plastic strains and, in particular, volumetric plastic strains, are smaller than if recoverable deformation induced anisotropy had been neglected. The result is an associated plasticity model, which appears to be non‐associated. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
14.
Poro‐elasto‐plastic response of an unconsolidated formation confined with stiff seal rocks under radial injection 下载免费PDF全文
下载免费PDF全文 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. 相似文献
15.
This paper presents an elasto‐plastic model for non‐linear analyses of cement‐treated sand. Various laboratory tests were systematically carried out to investigate the pre‐peak and post‐peak behaviours of a cement‐treated sand. On the basis of these experimental results, the new model was built within the framework of a relatively simple elasto‐plastic theory. Two failure criteria are employed to express tensile and shear failure characteristics observed in the experimental results of the cement‐treated sand. The proposed model can describe strain‐hardening and strain‐softening responses under both failure modes. In the strain‐softening rules, the smeared crack concept is used, and a characteristic length is considered to avoid the issue of mesh‐size dependency. Since the failure criterion and strain‐hardening/softening rules are based on the experimental evidences, the model is relatively easy to understand and the parameters used in the model have clear physical meaning. The proposed model was applied to simulate the behaviour of cement‐treated sand in various laboratory tests, allowing for a reasonable comprehensive evaluation. It was demonstrated that the proposed model is suitable for describing both the tensile and shear failure behaviours of cement‐treated sand. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
16.
Shear‐lag analysis is used to obtain closed‐form solutions for the problem of a stiff inclusion embedded in a poroelastic soil matrix. The following assumptions are made: the soil matrix and the inclusion are elastic; plane strain conditions apply; and shear stresses at the soil‐inclusion interface follow Coulomb's friction law. Two solutions are obtained, the first one for drained conditions where no excess pore pressures are generated, and the second one for undrained conditions where excess pore pressures are produced and the soil does not change volume during pullout. The solutions are verified by comparing analytical predictions with numerical results obtained using a finite element method. Predictions from the analytical solutions are also compared with results from experiments conducted in a large‐scale pullout box. Both comparisons show good agreement. The analytical solution shows that the pullout capacity in drained and undrained conditions is overall independent of the relative stiffness of the soil and the inclusion. The most important factor controlling the pullout capacity is the coefficient of friction between the soil and the inclusion. Both drained and undrained pullout capacities increase with the coefficient of friction; although the drained capacity shows a proportional increase, it is not so for the undrained capacity. The ratio of undrained to drained pullout capacity is about 0.9 for friction coefficients smaller than 0.2, but can be as small as 0.6 for a coefficient of friction of 1.0. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
17.
Timo Saksala 《国际地质力学数值与分析法杂志》2010,34(13):1362-1386
This paper presents a damage–viscoplastic cap model for rocks with brittle and ductile behavior under low‐velocity impact loading, which occurs, e.g. in percussive drilling. The model is based on a combination of the recent viscoplastic consistency model by Wang and the isotropic damage concept. This approach does not suffer from ill posedness—caused by strain softening—of the underlying boundary/initial‐value problem since viscoplasticity provides a regularization under dynamic loading by introducing an internal length scale. The model uses the Drucker–Prager (DP) yield function with the modified Rankine criterion as a tension cut‐off and a parabolic cap surface as a compression cut‐off. The parabolic cap is smoothly fitted to the DP cone. The strain softening law in compression is calibrated with the degradation index concept of Fang and Harrison. Thereby, the model is able to capture the brittle‐to‐ductile transition and hardening behavior of geomaterials under highly confined compression, which is the prevailing stress state under a bit‐button in percussive drilling. Rock strength heterogeneity is characterized statistically at the structural level using the Weibull distribution. An explicit time integrator is chosen for solving the FE‐discretized equations of motion. The contact constraints due to the impact of an indenter are imposed with the forward increment Lagrange multiplier method that is compatible with explicit time integrators. The model is tested at the material point level with various uniaxial and triaxial tests. At the structural level confined compression, uniaxial tension tests and a rock sample under low‐velocity impact are simulated. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
18.
在岩土工程中,裂隙岩体经常处于峰后变形状态,研究裂隙岩体的峰后应力-应变关系对预测结构的稳定性有重要意义。基于峰后软化阶段强度参数的逐渐演化行为,首先提出一个求岩石峰后应力-应变关系和裂隙峰后应力-切向位移关系的一般方法。然后采用摩尔-库仑强度准则,以岩石的最大主应变和裂隙的切向位移作为软化参数,假设强度参数为软化参数的分段线性函数,分岩体沿裂隙滑移破坏和沿岩石剪切破坏两种情况,提出多组贯穿裂隙岩体峰后应力-应变关系式的求法。最后,在算例中,分岩体沿裂隙滑移破坏和沿岩石剪切破坏两种情况给出了裂隙岩体的峰后应力-应变曲线,讨论了裂隙的平均间距、法向刚度和剪切刚度对峰后应变的影响。结果表明,裂隙的平均间距、法向刚度和剪切刚度越小,裂隙岩体的轴向应变越大。 相似文献
19.
A model of a multivariate covariance function with an ellipsoidal directional correlation scale has been developed. The axes of the ellipsoidal scale are related to the eigenvalues and eigenvectors of a matrix B which characterizes the ellipsoid of the range of influence. The matrix B is found to be related to a matrix T which can be estimated directly from sparse sampling data and can be used to determine estimates of the matrix B. The method has been applied to both two-dimensional and three-dimensional cases. The numerical results show that the satisfactory accuracy is obtained with sparse sampling data from an anisotropic random function. 相似文献
20.
Xianda Shen Jihui Ding Chlo Arson Judith S. Chester Frederick M. Chester 《国际地质力学数值与分析法杂志》2021,45(1):28-44
The dependence of rock behavior on the deformation rate is still not well understood. In salt rock, the fundamental mechanisms that drive the accumulation of irreversible deformation, the reduction of stiffness, and the development of hysteresis during cyclic loading are usually attributed to intracrystalline plasticity and diffusion. We hypothesize that at low pressure and low temperature, the rate‐dependent behavior of salt rock is governed by water‐assisted diffusion along grain boundaries. Accordingly, a chemo‐mechanical homogenization framework is proposed in which the representative elementary volume (REV) is viewed as a homogeneous polycrystalline matrix that contains sliding grain‐boundary cracks. The slip is related to the mass of salt ions that diffuse along the crack surface. The relationship between fluid inclusion‐scale and REV‐scale stresses and strains is established by using the Mori–Tanaka homogenization scheme. It is noted from the model that a lower strain rate and a larger number of sliding cracks enhance stiffness reduction and hysteresis. Thinner sliding cracks (i.e., thinner brine films) promote stiffness reduction and accelerate stress redistributions. The larger the volume fraction of the crack inclusions, the larger the REV deformation and the larger the hysteresis. Results presented in this study shed light on the mechanical behavior of salt rock that is pertinent to the design of geological storage facilities that undergo cyclic unloading, which could help optimize the energy production cycle with low carbon emissions. 相似文献