首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 390 毫秒
1.
The present paper examines the elastostatic problem pertaining to the axisymmetric loading of a rigid circular foundation resting on the surface of a non-homogeneous elastic half-space. The non-homogeneity corresponds to a depth variation in the linear elastic shear modulus according to the exponential form G(z)=G1+G2e-ζz. The equations of elasticity governing this type of non-homogeneity are solved by employing a Hankel transform technique. The mixed boundary value problem associated with the indentation of the half-space by the rigid circular foundation is reduced to a Fredholm integral equation which is solved via a numerical technique. The numerical results presented in the paper illustrate the influence of the near-surface elastic non-homogeneity on the settlement of the foundation.  相似文献   

2.
A two-dimensional analytical solution for stress, strain rate, and velocity is obtained for parallel-sided and wedge-shaped blocks with generalized viscous rheology (linearly viscous and power-law) deforming in plane strain. The main assumptions used in the derivation of the solution are that the material is incompressible, the longitudinal gradient in shear stress is much less than the vertical gradient of vertical normal stress, and the longitudinal strain rate varies linearly in the horizontal direction. Velocity boundary conditions are specified at the top of the block, and shear stress boundary conditions at the base of the block. In the one-dimensional case (where stress and strain rate do not vary in the longitudinal direction), the solution reduces to a well-known solution for the deformation of parallel-sided ice sheets [Nye, J. F. (1957) The distribution of stress and velocity in glaciers and ice sheets. Proceedings of the Royal Society of London A-239, 113–133]. The stress equilibrium for tapered wedges [Platt, J. P. (1986) Dynamics of orogenic wedges and the uplift of high-pressure metamorphic rocks. Geological Society of America Bulletin 97, 1037–1053] is a special case of the present stress solution. Implementation of the solution requires the subdivision of the wedge into vertical segments, and yields the tectonic normal and shear stresses that must be applied to the rear of a block with specified rheology in order to maintain a given longitudinal strain rate. The solution makes it possible to model deformation patterns analytically with longitudinally varying strain rate (including coeval compression and extension) and with vertical components of velocity reflecting the effects of underplating.  相似文献   

3.
王云岗  林宏剑 《岩土力学》2007,28(Z1):259-262
竖向圆形荷载作用下弹性半空间问题的位移和应力解是桩基分析的基础。利用Hankel积分变换,首先导出了弹性地基半空间位移与应力的积分形式的通解。通过适当地引入边界条件和界面位移和应力的连续条件,求得了内部作用竖向圆形荷载时弹性地基半空间位移与应力的积分形式解。在此基础上,给出了不同深度处荷载作用投影范围内竖向位移和竖向正应力的平均值。数值结果验证了解析解的正确性。  相似文献   

4.
By using the lower bound limit analysis in conjunction with finite elements and linear programming, the bearing capacity factors due to cohesion, surcharge and unit weight, respectively, have been computed for a circular footing with different values of ?. The recent axisymmetric formulation proposed by the authors under ?=0 condition, which is based on the concept that the magnitude of the hoop stress (σθ) remains closer to the least compressive normal stress (σ3), is extended for a general c–? soil. The computational results are found to compare quite well with the available numerical results from literature. It is expected that the study will be useful for solving various axisymmetric geotechnical stability problems. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

5.
The fully coupled Biot quasi-static theory of linear poroelasticity is used to study the consolidation of a poroelastic half-space caused by axisymmetric surface loads. The fluid and solid constituents of the poroelastic medium are compressible and its permeability in the vertical direction is different from its permeability in the horizontal direction. An analytical solution of the governing equations is obtained by taking the displacements and the pore pressure as the basic state variables and using a combination of the Laplace and Hankel transforms. The problem of an axisymmetric normal load is discussed in detail. An explicit analytical solution is obtained for normal disc loading. Detailed numerical computations reveal that the anisotropy in permeability as well as the compressibilities of the fluid and solid constituents of the poroelastic medium have significant effects on the consolidation of the half-space. The anisotropy in permeability may accelerate the consolidation process and may lead to a dilution in the theoretical prediction of the Mandel-Cryer effect. The compressibility of the solid constituents may also accelerate the consolidation process.  相似文献   

6.
The paper presents a comparison between the behavior of slurry-consolidated Kaolin specimens tested under axisymmetric Conventional Triaxial Compression (CTC) and Plane Strain (PS) loading conditions. The PS experiments were conducted on an instrumented apparatus capable of capturing the onset of shear bands. The specimens were consolidated and then sheared under undrained condition. The PS specimens failed via a well-defined shear band that began to develop during the hardening stress regime whereas CTC specimens failed through a diffuse bulging mode. The undrained shear strength (S u) of CTC experiments is smaller than the S u of PS experiments. However, the normalized S u of CTC experiments is very close to PS1, which was consolidated under similar K o condition.  相似文献   

7.
Activation of major faults in Bulgaria and northern Greece presents significant seismic hazard because of their proximity to populated centers. The long recurrence intervals, of the order of several hundred years as suggested by previous investigations, imply that the twentieth century activation along the southern boundary of the sub-Balkan graben system, is probably associated with stress transfer among neighbouring faults or fault segments. Fault interaction is investigated through elastic stress transfer among strong main shocks (M ≥ 6.0), and in three cases their foreshocks, which ruptured distinct or adjacent normal fault segments. We compute stress perturbations caused by earthquake dislocations in a homogeneous half-space. The stress change calculations were performed for faults of strike, dip, and rake appropriate to the strong events. We explore the interaction between normal faults in the study area by resolving changes of Coulomb failure function (ΔCFF) since 1904 and hence the evolution of the stress field in the area during the last 100 years. Coulomb stress changes were calculated assuming that earthquakes can be modeled as static dislocations in an elastic half-space, and taking into account both the coseismic slip in strong earthquakes and the slow tectonic stress buildup associated with major fault segments. We evaluate if these stress changes brought a given strong earthquake closer to, or sent it farther from, failure. Our modeling results show that the generation of each strong event enhanced the Coulomb stress on along-strike neighbors and reduced the stress on parallel normal faults. We extend the stress calculations up to present and provide an assessment for future seismic hazard by identifying possible sites of impending strong earthquakes.  相似文献   

8.
Summary. Most currently used techniques for analysing the stability of near surface structures, such as rock slopes, are based on the application of the effective Coulomb shear strength parameters cohesion c′, and the angle of friction φ′ on some known or anticipated shear surface subjected to an effective normal stress σ′n. The most widely used of these techniques are the variants of the method of slices and related upper bound techniques. If the Hoek-Brown criterion is to be used to model the strength of near surface fractured rocks, it is necessary to determine equivalent Coulomb shear strength parameters for the specified level of effective normal stress. Calculation of the equivalent Coulomb parameters for the Hoek-Brown criterion for cases when a ≠ 0.5 is not a straightforward matter. A simple procedure for calculating instantaneous values of c′i and φ′i has been developed based on spreadsheet calculations and the application of a numerical optimisation routine. This procedure can also be applied to calculating the Hoek-Brown envelope plotted in shear stress/normal stress space. A simple closed form solution for c′i and tan φ′i has also been developed for the special case when a = 1. A three-dimensional version of the Hoek-Brown criterion has been developed by combining it with the Drucker-Prager criterion. This new yield criterion has been implemented by numerical solution of the governing equations. A simplification of this three-dimensional yield criterion has been developed by introducing an intermediate principal stress weighting factor. Comparison with published results demonstrates that this simplified criterion has the capacity to model the results of true triaxial tests for a range of different rock types over a wide range of stress levels. The new three-dimensional yield criterion has the advantage that its input parameters can be determined from routine uniaxial compression tests and mineralogical examination.  相似文献   

9.
A boundary element approach for time harmonic axisymmetric problems using the complete space point load fundamental solution is presented. The fundamental solution is integrated numerically along the azimuthal co-ordinate of each axisymmetric element. To increase the accuracy of the numerical integration a simple co-ordinate transformation is proposed. The approach is applied to the computation of the dynamic stiffness functions of rigid circular foundations on layered viscoelastic soils. Three different sites are considered: a uniform half-space, a soil layer on a half-space, and a soil consisting of four horizontal layers and a compliant half-space. The numerical results obtained by the proposed approach for surface circular foundations are very close to corresponding published results obtained by different procedures.  相似文献   

10.
We consider the multigrid solution of a discretized singular integral equation whose solution represents the displacement discontinuity distribution across a pressurized crack or the slip on a fault subjected to a prescribed shear stress. The multigrid technique reduces the operation count for a crack model having N degrees of freedom from O(N3) operations for standard stationary iterative methods to O(N2) operations. In the numerical simulations performed the multigrid approach proves to be extremely efficient even for small values of N. We use Fourier analysis to determine the spectral properties of the coarse grid correction process and the effect of a number of different interpolation operators on the multigrid algorithm. We show that the multigrid technique can be combined with the process of lumping remote influences to yield an algorithm that involves O(N) operations. The performance of multigrid iteration in a nonlinear environment is explored by considering seams filled with nonlinear material. For this purpose a segmented multigrid algorithm is developed that allows for different seam constitutive relations to be used along different line segments. The O(N2) operation count characteristic of linear multigrid iteration is shown to persist in this nonlinear environment and the segmented multigrid approach is shown to provide significant computational savings for values of N well within the range for typical problems that occur in practice.  相似文献   

11.
Quantitative Parameters for Rock Joint Surface Roughness   总被引：17，自引：5，他引：12  
Summary The morphologies of two artificial granite joints (sanded and hammered surfaces), one artificial regularly undulated joint and one natural schist joint, were studied. The sanded and hammered granite joints underwent 5 cycles of direct shear under 3 normal stress levels ranging between 0.3–4 MPa. The regularly undulated joint underwent 10 cycles of shear under 6 normal stress levels ranging between 0.5–5 MPa and the natural schist replicas underwent a monotonous shear under 5 normal stress levels ranging between 0.4–2.4 MPa. In order to characterize the morphology of the sheared joints, a laser sensor profilometer was used to perform surface data measurements prior to and after each shear test. Rather than describing the morphology of the joints from the single profiles, our characterization is based on a simultaneous analysis of all the surface profiles. Roughness was viewed as a combination of a primary roughness and a secondary roughness. The surface angularity was quantified by defining its three-dimensional mean angle, θs, and the parameter Z2s. The surface anisotropy and the secondary roughness were respectively quantified by the degree of apparent anisotropy, k a, and the surface relative roughness coefficient, R s. The surface sinuosity was quantified by the surface tortuosity coefficient, T s.  Comparison between the means of the classical linear parameters and those proposed shows that linear parameters underestimate the morphological characteristics of the joint surfaces. As a result, the proposed bi-dimensional and tri-dimensional parameters better describe the evolution of the joints initial roughness during the course of shearing.  相似文献   

12.
H.J. Melosh 《Tectonophysics》1976,35(4):363-390
This paper investigates the effect of shear heating in the asthenosphere on the thermal structure of the upper mantle. Equations describing the motion of the lithosphere over the asthenosphere in the presence of a strongly temperature-dependent stress-strain rate relation are derived and solved with the help of several approximations. These approximations are shown to be valid under conditions appropriate for the earth.Two sets of solutions are found. For one set (the “subcritical” solutions) a normal shear stress—velocity relation is found for small stresses. The velocity increases as the stress increases, reaching a maximum velocity σc for a critical stress σc. The subcritical solutions have a negligible effect on the thermal structure of the earth, even at the critical stress. The other set of solutions (the “supercritical” solutions) has the bizarre property that a decrease of applied shear stress leads to an increase of velocity. Thus, as the shear stress goes to zero, the velocity becomes infinite. At larger shear stresses the velocity decreases until it reaches σc at a stress σc (the two sets of solutions share this point in common). There are no steady solutions of any kind for shear stresses in excess of σc. We discard the supercritical solutions as candidates for the thermal structure of the earth on the basis of their instability to small perturbations of applied stress and temperature.The realm of subcritical solutions (stress less than σc, velocity less than σc) thus defines a regime of plate motion in which the thermal effects of shear heating are negligible. If the shear stresses acting on plates exceed σc, however, new physical processes must come into play to dissipate the excess heat generated. Assuming that the velocities of plates on the earth today are less than σc, relative to the deep mantle, a strict upper limit of a few tens of bars can be derived for σc, corresponding to effective viscosities of ca. 1019 poise in the asthenosphere.  相似文献   

13.
A simplified analysis of the problem of horizontal soil stress changes around circular displacement piles, caused by pile placement, is presented. Pile installation is assumed to cause soil displacements in the horizontal direction only, thus yielding an axisymmetric problem. The soil surrounding the pile is modelled as a weightless non-linear elastic material. Material non-linearity is handled in a simplified manner by adopting secant shear moduli defined in terms of a proportionality coefficient and a softening factor. The resulting equilibrium equation is solved analytically and an expression is obtained which is also conveniently presented in graphical form. The derived expression can be used to estimate horizontal soil stress changes and is incorporated into a simple procedure to estimate the ultimate load or the efficiency of pile groups. Comparisons are made between efficiencies calculated according to this procedure and efficiencies measured in full-scale group load tests in sand.  相似文献   

14.
The paper deals with the propagation of shock waves at the surface of soils. Heterogeneity and damping are introduced into analytical half-space solutions. The suggested model explains two phenomena, often observed with shock propagation in actual soils, that differ from the behaviour of the homogeneous half-space: the pronounced decay of the disturbances with distance and the elongation of the disturbance into a train of waves. The effects of heterogeneity and damping are discussed quantitatively. The response of footings on heterogeneous soils has been investigated by several authors. Awojobi4 considered the Gibson soil in which the shear modulus increases linearly with depth. Luco5 and Gazetas and Roesset6 investigated a multi-layered soil, the shear modulus being constant within each layer. Gazetas7, using a technique suggested by Gupta8 extended this method to layers with linearly varying shear modulus. Little work is available on the propagation of waves in heterogeneous bodies. Some results concerning the modes and the mode shapes in heterogeneous soils were reported by Ewing, Jardetzky and Press9 and Bath.10 The modes have to be superposed in an appropriate way to obtain the displacement field at the surface. This has been approximately achieved by the finite element formulations of Lysmer, 11, Lysmer and Waas12 and Waas.13 Auersch14 applied this method to a homogeneous layer on a rigid base. He found some dispersion of the Rayleigh wave within a narrow frequency range. Finite elements combined with discrete Laplace transforms, however, consume much computer time. Rao and Goda15 and Rao16 calculated surface vibrations of a half-space with exponentially varying shear modulus and density. Their method is similar to Lamb'S1 procedure for the homogeneous half-space. Only one mode–the Rayleigh wave–occurs in their heterogeneous half-space. The examples show the considerable effect of heterogeneity on wave propagation. In the present paper, more general variations of the shear modulus are considered.  相似文献   

15.
Vertical and horizontal deformations of surface footings have been studied for an inhomogeneous elastic half-space in which the shear modulus increases with an arbitrary power of depth, n, and Poisson's ratio is constant. A general solution for displacements has been obtained first for point loads applied in vertical and horizontal directions. These are then used in obtaining closed-form solutions for displacements of uniformly loaded circular and rectangular footings. Finally, a numerical method is described that can be used to analyse a rigid footing of an arbitrary shape, and results for rigid rectangular footings are given.  相似文献   

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

17.
Earlier solutions (Bull. Seismol. Soc. Amer. 1985; 75 :1135–1154; Bull. Seismol. Soc. Amer. 1992; 82 :1018–1040) of deformations caused by the movements of a shear or tensile fault in an isotropic half-space for finite rectangular sources of strain nucleus have been extended for a transversely isotropic half-space. Results of integrating previous solutions (Int. J. Numer. Anal. Meth. Geomech. 2001; 25 (10): 1175–1193) of deformations due to a shear or tensile fault in a transversely isotropic half-space for point sources of strain nucleus over the fault plane are presented. In addition, a boundary element (BEM) model (POLY3D:A three-dimensional, polygonal element, displacement discontinuity boundary element computer program with applications to fractures, faults, and cavities in the Earth's crust. M.S. Thesis, Stanford University, Department of Geology, 1993; 62) is given. Different from similar researches (e.g. Thomas), the Akaike's view on Bayesian statistics (Akaike Information Criterion Statistics. D. Reidel Publication: Dordrecht, 1986) is applied for inverting deformations due to a fault to obtain displacement discontinuities on the fault plane. An example is given for checking displacements predicted by proposed analytical expressions. Another example is generated for the use of proposed BEM model. It demonstrates the effectiveness of this model in exploring displacement behaviours of a fault. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

18.
层状地基中单桩负摩擦问题积分方程解法   总被引：9，自引：1，他引：8  
高绍武  王建华  毛娜 《岩土力学》2005,26(9):1456-1460
利用Biot固结理论和积分方程方法研究了表面有堆载的层状地基中单桩负摩擦问题。根据层状饱和土的圆形载荷基本解得出了单桩在圆形均布载荷作用下在时间域内的第二类Fredholm积分方程组。运用Laplace变换对上述积分方程组进行简化。再结合传递和刚度矩阵传递到各个层中去,对变换域内的积分方程采用Schapery 逆变换方法得到时域内单桩的近似积分方程。求解积分方程组并进行相应的数值逆变换,就可得出层状地基中的单桩在表面圆形均布载荷作用下的位移、轴力、孔压和桩侧摩阻力随时间的变化情况。计算结果表明,桩侧剪力和孔压分层明显。  相似文献   

19.
李镜培  方睿  李林  唐剑华 《岩土力学》2016,37(3):679-696
以静压沉桩后桩周土体的应力状态为初始条件,根据桩周土体孔隙比、渗透系数和有效应力之间的相关性,在考虑固结系数随固结时间变化的条件下改进了轴对称条件下的太沙基固结控制方程。随后,采用分离变量法和离散化分析推导得出了桩周超孔隙水压力消散的半解析半数值解,并将解答与实测数据进行对比验证。在此基础上,采用空间滑动面理论改进的修正剑桥模型（SMP-MCC）定义土体三维不排水抗剪强度,研究了静压桩周土体强度、剪切模量随固结时间的变化规律。研究结果表明：由于解答考虑了固结系数随固结时间的变化,因而与实测结果吻合良好;土体压缩指数与渗透指数之比对土体固结系数和孔压消散速率影响较大;当土体压缩指数与渗透指数之比为1时,土体固结系数保持不变,解答退化为经典的太沙基轴对称固结方程;土体强度和剪切模量随固结时间的增长而逐步增加,固结完成后其值超越了土体原位强度和原位剪切模量。  相似文献   

20.
The Biot linearized quasi-static theory of fluid-infiltrated porous materials is used to formulate the problem of the two-dimensional plane strain deformation of a multi-layered poroelastic half-space by surface loads. The Fourier-Laplace transforms of the stresses, displacements, pore pressure and fluid flux in each homogeneous layer of the multi-layered half-space are expressed in terms of six arbitrary constants. Generalized Thomson-Haskell matrix method is used to obtain the deformation field. Simplified explicit expressions for the elements of the 6 × 6 propagator matrix for the poroelastic medium are obtained. As an example of the possible applications of the analytical formulation developed, formal solution is given for normal strip loading, normal line loading and shear line loading.  相似文献   

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

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