首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 609 毫秒
1.
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.  相似文献   

2.
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.  相似文献   

3.
This paper presents a three‐dimensional energy‐based solution for the time‐dependent response of a deeply embedded and unsupported semi‐infinite tunnel of circular cross‐section. The tunnel is taken to be excavated quasi‐instantaneously from an infinite rock body that initially exhibits an isotropic stress state and that is made up of a homogeneous, isotropic and viscoelastic material. The viscoelastic behaviour is modelled by means of Burger's model, and the rock is taken to behave volumetrically linear elastic and to exhibit exclusively deviatoric creep. This viscoelastic problem is transformed into the Laplace domain, where it represents a quasi‐elastic problem. The displacement fields in the new solution are taken to be the products of independent functions that vary in the radial and longitudinal directions. The differential equations governing the displacements of the system and appropriate boundary conditions are obtained using the principle of minimum potential energy. The solutions for these governing equations in the Laplace domain are then obtained analytically and numerically using a one‐dimensional finite difference technique. The results are then transformed back into the time domain using an efficient numerical scheme. The accuracy of the new solution is comparable with that of a finite element analysis but requires much less computation effort. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

4.
An analytical investigation of a half‐space containing transversely isotropic material under forced vertical and horizontal displacements applied on a rectangular rigid foundation is presented in this paper. With the goal of a rigorous solution to the shape‐ and rigidity‐ induced singular mixed boundary value problem, the formulation employs scalar potential representation, the Fourier expansion and the Hankel integral transforms method to obtain the surface arbitrary point‐load solution in cylindrical coordinate system. The obtained Green's functions are rewritten in rectangular coordinate system, allowing the response of the half‐space because of an arbitrary distributed load on a rectangular surface area be given in terms of a double integral. The numerical evaluations of stresses are done with the use of an element, which is singular at the edge and the corner of the rectangle. Upon the imposition of the rigidity displacement boundary condition for a rigid foundation and the use of a set of two‐dimensional adaptive‐gradient elements, which can capture the singular behavior in the contact stress effectively, a set of new numerical results are presented to illustrate the effect of transverse isotropy on the foundation response. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

5.
The scaled boundary finite‐element method, a semi‐analytical computational scheme primarily developed for dynamic stiffness of unbounded domains, is applied to the analysis of unsteady seepage flow problems. This method is based on the finite‐element technology and gains the advantages of the boundary element method as well. Only boundary of the domain is discretized, no fundamental solution is required and singularity problems can be modeled rigorously. Anisotropic and non‐homogeneous materials satisfying similarity are modeled with no additional efforts. In this study, firstly, formulation of the method for the transient seepage flow problems is derived followed by its solution procedures. The accuracy, simplicity and applicability of the method are demonstrated via four numerical examples of transient seepage flow – three of them are available in the literature. Homogenous, non‐homogenous, isotropic and anisotropic material properties are considered to show the versatility of the technique. Excellent agreement with the finite‐element method is observed. The method out‐performs the finite‐element method in modeling singularity points. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

6.
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.  相似文献   

7.
An analytical solution of the plane strain problem of the deformation of a homogeneous, isotropic, poroelastic layer of uniform thickness overlying a homogeneous, isotropic, elastic half‐space due to two‐dimensional seismic sources buried in the elastic half‐space has been obtained. The integral expressions for the displacements, stresses and pore pressure have been obtained using the stress function approach by applying suitable boundary conditions at the free surface and the interface. The solution obtained is in the Laplace–Fourier transform domain. The case of a vertical dip‐slip line dislocation for the oceanic crust model of Earth is studied in detail. Schapery's formula is used for the Laplace inversion and the extended Simpson's formula for the Fourier inversion. Diffusion of pore pressure in the layer is studied numerically. Contour maps showing the pore pressure in the poroelastic layer have been plotted. The effect of the compressibility of the solid and fluid constituents on pore pressure has also been studied. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

8.
An analytical solution is presented in this paper to study the time‐dependent settlement behaviour of a rigid foundation resting on a transversely isotropic saturated soil layer. The governing equations for a transversely isotropic saturated soil, within Biot's poroelasticity framework, are solved by means of Laplace and Hankel transforms. The problem is subsequently formulated in the Laplace transform domain in terms of a set of dual integral equations that are further reduced to a Fredholm integral equation of the second kind and solved numerically. The developed analytical solution is validated via comparison with the existing analytical solution for an isotropic saturated soil case, and adopted as a benchmark to examine the sensitivities of the mesh refinement and the locations of truncation boundaries in the finite element simulations using ABAQUS. Particular attention is paid to the influences of the degree of soil anisotropy, boundary drainage condition, and the soil layer thickness on the consolidation settlement and contact stress of the rigid foundation. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

9.
An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

10.
In this paper, a single-region boundary element method (BEM) is presented for the analysis of transversely isotropic bi-material halfspaces with arbitrarily inclined planes of isotropy and material interfaces. The proposed BEM uses the fundamental solution of a transversely isotropic bi-material fullspace and five boundary element techniques. Infinite boundary elements are introduced to consider the far-fields of a transversely isotropic bi-material halfspace. The effective integration methods are proposed for dealing with various integrals in the discretized boundary integral equation. The stresses at internal points are obtained using the coordinate transformation of kernel functions, and the stresses on the boundary surface are calculated using an improved traction recovered method. Numerical verifications of displacements and stresses for a benchmark problem are conducted, and excellent agreement with previously published results is obtained. Numerical examples are presented to illustrate the influence of non-horizontal or horizontal planes of isotropy in bi-material halfspaces on the displacements and stresses induced by the tractions on the horizontal boundary surface. Results reveal that the elastic fields vary clearly with the dip angle of the isotropic plane and the stresses across the bi-material interface are closely related to the ratios of the elastic parameters of the bi-material.  相似文献   

11.
In this paper, a closed-form solution is presented for the stress and displacement distributions throughout a hollow cylinder subjected to uniform pressures acting on its internal and external boundary surfaces under plane strain conditions. The material is assumed to be elastoplastic, obeying a Mohr–Coulomb failure criterion, and exhibiting dilatant plastic deformation according to a non-associated flow rule. The newly developed analytical solution is verified through comparison with the solutions obtained from an infinite boundary problem (for which a closed-form solution exists), and numerical analyses using the program FLAC. The solution is also compared with the results of a borehole collapse test on a thick-walled hollow cylinder of synthetic shale. The analytical solution can be used to calculate the stress and displacement distributions around boreholes and other cylindrical cavities under both infinite and finite boundary conditions under both drained and undrained conditions. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

12.
A theoretical formulation is presented for the determination of the interaction of a vertically loaded disc embedded in a transversely isotropic half‐space. By means of a complete representation using a displacement potential, it is shown that the governing equations of motion for this class of problems can be uncoupled into a fourth‐order partial differential equation. With the aid of Hankel transforms, a relaxed treatment of the mixed‐boundary value problem is formulated as dual integral equations, which, in turn, are reduced to a Fredholm equation of the second kind. In addition to furnishing a unified view of existing solutions for zero and infinite embedments, the present treatment reveals a severe boundary‐layer phenomenon, which is apt to be of interest to this class of problems in general. The present solutions are analytically in exact agreement with the existing solutions for a half‐space with isotropic material properties. To confirm the accuracy of the numerical evaluation of the integrals involved, numerical results are included for cases of different degrees of the material anisotropy and compared with existing solutions. Further numerical examples are also presented to elucidate the influence of the degree of the material anisotropy on the response. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

13.
This paper investigates the two‐dimensional flow problem through an anisotropic porous medium containing several intersecting curved fractures. First, the governing equations of steady‐state fluid flow in a fractured porous body are summarized. The flow follows Darcy's law in matrix and Poiseuille's law in fractures. An infinite transversal permeability is considered for the fractures. A multi‐region boundary element method is used to derive a general pressure solution as a function of discharge through the fractures and the pressure and the normal flux on the domain boundary. The obtained solution fully accounts for the interaction and the intersection between fractures. A numerical procedure based on collocation method is presented to compute the unknowns on the boundaries and on the fractures. The numerical solution is validated by comparing with finite element solution or the results obtained for an infinite matrix. Pressure fields in the matrix are illustrated for domains containing several interconnected fractures, and mass balance at the intersection points is also checked. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

14.
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.  相似文献   

15.
Similar to plane strain, axisymmetric stress problem is also highly kinematics constrained. Standard displacement‐based finite element exhibits volumetric locking issue in simulating nearly/fully incompressible material or isochoric plasticity under axisymmetric loading conditions, which severely underestimates the deformation and overestimates the bearing capacity for structural/geotechnical engineering problems. The aim of this paper is to apply variational multiscale method to produce a stabilized mixed displacement–pressure formulation, which can effectively alleviate the volumetric locking issue for axisymmetric stress problem. Both nearly incompressible elasticity and isochoric J2 elastoplasticity are investigated. First‐order 3‐node triangular and 4‐node quadrilateral elements are tested for locking issues. Several representative simulations are provided to demonstrate the performance of the linear elements, which include the convergence study and comparison with closed‐form solutions. A comparative study with pressure Laplacian stabilized formulation is also presented. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

16.
王少杰  吕爱钟  张晓莉 《岩土力学》2018,39(12):4437-4447
将围岩和衬砌分别视作均质、连续的线弹性正交各向异性和各向同性体,并充分考虑衬砌的支护滞后效应和隧洞运行时的内水压力作用,运用复变函数方法中的幂级数解法,提出了正交各向异性岩体中任意形状水工隧洞的力学解析方法。以直墙半圆拱形水工隧洞为例,所获得的解析解可精确满足衬砌内边界的应力边界条件以及围岩与衬砌接触面的应力、位移连续条件,同时还将解析结果与ANSYS数值结果对比分析,吻合良好。利用获得的解析解,讨论了围岩开挖面上不同的各向异性程度、不同的弹性对称面角度以及隧洞内不同的水压荷载对衬砌以及围岩上应力和位移分布的影响。  相似文献   

17.
A superposition scheme is proposed to obtain a fundamental solution for boundary elements in multi‐layered elastic media. A three‐layered elastic region is obtained by superposing two sets of bonded half‐planes and subtracting one infinite plane. Therefore, the solution for an element in a layered media can be expressed in terms of bonded half‐plane solutions and an infinite‐plane solution. The major advantages of this superposition scheme are: (1) it is unnecessary to introduce elements at the interface, (2) it can be extended to higher‐order element, and (3) it may be applicable to three dimensions easily. The accuracy and performance of the developed model is illustrated by two examples. For the problem of a pressurized two‐dimensional crack within a three‐layered system, the comparison with other numerical results shows the model is quite accurate and efficient. The model is also used for a study of a practical two‐dimensional mining problem in South Africa, i.e. stoping through a dyke with material properties different from the host rock. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

18.
A transversely isotropic linear elastic half‐space, z?0, with the isotropy axis parallel to the z‐axis is considered. The purpose of the paper is to determine displacements and stresses fields in the interior of the half‐space when a rigid circular disk of radius a completely bonded to the surface of the half‐space is rotated through a constant angle θ0. The region of the surface lying out with the circle r?a, is free from stresses. This problem is a type of Reissner–Sagoci mixed boundary value problems. Using cylindrical co‐ordinate system and applying Hankel integral transform in the radial direction, the problem may be changed to a system of dual integral equations. The solution of the dual integral equations is obtained by an approach analogous to Sneddon's (J. Appl. Phys. 1947; 18 :130–132), so that the circumferential displacement and stress fields inside the medium are obtained analytically. The same problem has already been approached by Hanson and Puja (J. Appl. Mech. 1997; 64 :692–694) by the use of integrating the point force potential functions. It is analytically proved that the present solution, although of a quite different form, is equivalent to that given by Hanson and Puja. To illustrate the solution, a few plots are provided. The displacements and the stresses in a soil deposit due to a rotationally symmetric force or boundary displacement may be obtained using the results of this paper. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

19.
In spite of the gradual development of three-dimensional analysis packages utilizing finite element models or finite difference algorithms for stress–strain calculations, two-dimensional (2D) analysis is still used as the primary engineering tool for practical analysis of tunnel behavior and tunnel support performance for design—particularly at the preliminary stage of a project. The applicability of 2D finite element analysis or analytical convergence confinement solutions to staged support installation depend on the application of an assumed or validated longitudinal displacement profile. Convergence in commonly applied 2D staged models is controlled by boundary displacement or internal pressure relaxation. While there have been developments to improve this methodology, this often assumes independence between the ground reaction curve and the support resistance, independence between the longitudinal displacement profile to support application, and the assumption that non-isotropic stresses and non-circular geometries can be handled in the same way as circular tunnels in isotropic conditions. This paper examines the validity of these assumptions and the error inherent in these extensions to 2D tunnel analysis. Anisotropic stresses and lagged (staged) excavation present a particular problem. Practical solutions are proposed for support longitudinal displacement LDPs in simplified conditions.  相似文献   

20.
In this paper a micro‐polar continuum approach is proposed to model the essential properties of cohesionless granular materials like sand. The model takes into account the influence of particle rotations, the mean grain size, the void ratio, the stresses and couple stresses. The constitutive equations for the stresses and couple stresses are incrementally non‐linear and based on the concept of hypoplasticity. For plane strain problems the implementation of the model in a finite element program is described. Numerical studies of the evolution of micro‐polar effects within a granular strip under plane shearing are presented. It is shown that the location and evolution of shear localization is strongly influenced by the initial state and the micro‐polar boundary conditions. For large shearing the state quantities tend towards a stationary state for which a certain coupling between the norm of the stress deviator and the norm of the couple stress tensor can be derived. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

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

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