共查询到20条相似文献，搜索用时 41 毫秒
1.
Keh‐Jian Shou 《国际地质力学数值与分析法杂志》2000,24(10):795-814
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. 相似文献
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.
A finite element approach based on an advanced multi‐surface kinematic constitutive model is used to evaluate the bearing capacity of footings resting on granular soils. Unlike simple elastic‐perfectly plastic models, often applied to granular foundation problems, the present model realistically accounts for stress dependency of the friction angle, strain softening–hardening and non‐associativity. After the model and its implementation into a finite element code are briefly discussed, the numerical difficulty due to the singularity at the footing edge is addressed. The bearing capacity factor Nγ is then calculated for different granular materials. The effect of footing size, shape, relative density and roughness on the ultimate bearing capacity are studied and the computed results compare very favourably with the general experimental trends. In addition, it is shown that the finite element solution can clearly represent counteracting mechanisms of progressive failure which have an important effect on the bearing capacity of granular foundations. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
4.
This paper presents a non‐linear soil–structure interaction (SSI) macro‐element for shallow foundation on cohesive soil. The element describes the behaviour in the near field of the foundation under cyclic loading, reproducing the material non‐linearities of the soil under the foundation (yielding) as well as the geometrical non‐linearities (uplift) at the soil–structure interface. The overall behaviour in the soil and at the interface is reduced to its action on the foundation. The macro‐element consists of a non‐linear joint element, expressed in generalised variables, i.e. in forces applied to the foundation and in the corresponding displacements. Failure is described by the interaction diagram of the ultimate bearing capacity of the foundation under combined loads. Mechanisms of yielding and uplift are modelled through a global, coupled plasticity–uplift model. The cyclic model is dedicated to modelling the dynamic response of structures subjected to seismic action. Thus, it is especially suited to combined loading developed during this kind of motion. Comparisons of cyclic results obtained from the macro‐element and from a FE modelization are shown in order to demonstrate the relevance of the proposed model and its predictive ability. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
5.
将地基视为多层各向同性弹性体,对Euler-Bernoulli梁进行有限单元离散分析,对地基-梁接触面采用边界积分法求解,根据地基-梁接触面的竖向位移协调和光滑接触条件,应用有限元与边界元耦合的方法推导出各向同性成层弹性地基上的Euler-Bernoulli梁的半数值半解析解。基于文中理论,编制了相应的程序,通过与现有文献对比,验证了文中理论的正确性,对比分析了分层地基与等效均匀地基两种模型上的弹性地基梁。研究结果表明,分层地基与等效均匀地基两种模型上弹性地基梁性状差异较大,实际弹性地基梁计算中应采用分层地基模型。 相似文献
6.
采用双参数地基模型和基于Mindlin板理论的8节点等参单元,根据Hamilton变分原理,推导了弹性地基上弹性板动力分析的有限元公式,将地基效应和板面质量分布归并到板的刚度矩阵和质量矩阵中,采用边界元法考虑板域外地基变形的影响。分析了四边自由矩形板的自由振动和强迫振动,计算的自由振动固有频率与实测结果相近,发生共振时板的强迫振动动力响应达到峰值,在给定的阻尼比条件下,无阻尼的稳态响应与有阻尼的稳态响应趋于相同,阻尼能够显著降低板的共振响应。 相似文献
7.
Sarah E. M. Gain Yoann Grau Hadrien Henry Elena Belousova Ivars Dainis William L. Griffin Suzanne Y. O'Reilly 《Geostandards and Geoanalytical Research》2019,43(3):339-354
Zircon megacrysts from the Mud Tank carbonatite, Australia, are being used in many laboratories as a reference material for LA‐ICP‐MS U‐Pb dating and trace element measurement, and LA‐MC‐ICP‐MS determination of Hf isotopes. We summarise a database of > 10000 analyses of Mud Tank zircon (MTZ), collected from 2000 to 2018 during its use as a secondary reference material for simultaneous U‐Pb and trace element analysis, and for Hf‐isotope analysis. Trace element mass fractions are highest in dark red‐brown stones and lowest in colourless and gem‐quality ones. Individual unzoned grains can be chemically homogeneous, while significant variations in trace element mass fraction are associated with oscillatory zoning. Chondrite‐normalised trace element patterns are essentially parallel over large mass fraction ranges. A Concordia age of 731.0 ± 0.2 Ma (2s, n = 2272) is taken as the age of crystallisation. Some grains show lower concordant to mildly discordant ages, probably reflecting minor Pb loss associated with cooling and the Alice Springs Orogeny (450–300 Ma). Our weighted mean 176Hf/177Hf is 0.282523 ± 10 (2s, n = 9350); the uncertainties on this ratio reflect some heterogeneity, mainly between grains. A few analyses suggest that colourless grains have generally lower 176Hf/177Hf. MTZ is a useful secondary reference material for U‐Pb and Hf‐isotope analysis, but individual grains need to be carefully selected using CL imaging and tested for homogeneity, and ideally should be standardised by solution analysis. 相似文献
8.
In recent years, the authors have proposed a new double‐node zero‐thickness interface element for diffusion analysis via the finite element method (FEM) (Int. J. Numer. Anal. Meth. Geomech. 2004; 28 (9): 947–962). In the present paper, that formulation is combined with an existing mechanical formulation in order to obtain a fully coupled hydro‐mechanical (or HM) model applicable to fractured/fracturing geomaterials. Each element (continuum or interface) is formulated in terms of the displacements (u) and the fluid pressure (p) at the nodes. After assembly, a particular expression of the traditional ‘u–p’ system of coupled equations is obtained, which is highly non‐linear due to the strong dependence between the permeability and the aperture of discontinuities. The formulation is valid for both pre‐existing and developing discontinuities by using the appropriate constitutive model that relates effective stresses to relative displacements in the interface. The system of coupled equations is solved following two different numerical approaches: staggered and fully coupled. In the latter, the Newton–Raphson method is used, and it is shown that the Jacobian matrix becomes non‐symmetric due to the dependence of the discontinuity permeability on the aperture. In the part II companion paper (Int. J. Numer. Anal. Meth. Geomech. 2008; DOI: 10.1002/nag.730 ), the formulation proposed is verified and illustrated with some application examples. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
9.
Wei Dong Guo 《国际地质力学数值与分析法杂志》2000,24(2):135-163
Viscoelastic or creep behaviour can have a significant influence on the load transfer (t–z) response at the pile–soil interface, and thus on the pile load settlement relationship. Many experimental and theoretical models for pile load transfer behaviour have been presented. However, none of these has led to a closed‐form expression which captures both non‐linearity and viscoelastic behaviour of the soil. In this paper, non‐linear viscoelastic shaft and base load transfer (t–z) models are presented, based on integration of a generalized viscoelastic stress–strain model for the soil. The resulting shaft model is verified through published field and laboratory test data. With these models, the previous closed‐form solutions evolved for a pile in a non‐homogeneous media have been readily extended to account for visco‐elastic response. For 1‐step loading case, the closed‐form predictions have been verified extensively with previous more rigorous numerical analysis, and with the new GASPILE program analysis. Parametric studies on two kinds of commonly encountered loading: step loading, ramp (linear increase followed by sustained) loading have been performed. Two examples of the prediction of the effects of creep on the load settlement relationship by the solutions and the program GASPILE, have been presented. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
10.
A. P. S. SELVADURAI 《国际地质力学数值与分析法杂志》1996,20(5):351-364
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. 相似文献
11.
针对多体相互作用体系的非连续变形分析问题和接触问题,采用Mohr-Coulomb屈服准则和关联流动法则以及接触界面上的非线性应力分布模式,考虑接触界面特性提出了非线性接触力元模型,以结点位移和界面相互接触应力同时作为独立未知变量,建立了离散系统的总体控制方程.进而,通过数值求解能够直接确定变形体内的应力与变形、界面上的接触应力与离散体的位移与运动.将这种以接触力元为基础的多体系统分析方法具体应用于基础与地基相互作用分析,通过数值计算与分析探讨了地基与基础的相对刚度、荷载大小及其偏心距、地基与基础间界面力学参数对接触界面的应力分布和地基变形的影响,所得结果为工程中考虑基础与地基相互作用影响的设计与分析提供了参考依据. 相似文献
12.
The mathematical structure and numerical analysis of classical small deformation elasto–plasticity is generally well established. However, development of large deformation elastic–plastic numerical formulation for dilatant, pressure sensitive material models is still a research area. In this paper we present development of the finite element formulation and implementation for large deformation, elastic–plastic analysis of geomaterials. Our developments are based on the multiplicative decomposition of the deformation gradient into elastic and plastic parts. A consistent linearization of the right deformation tensor together with the Newton method at the constitutive and global levels leads toward an efficient and robust numerical algorithm. The presented numerical formulation is capable of accurately modelling dilatant, pressure sensitive isotropic and anisotropic geomaterials subjected to large deformations. In particular, the formulation is capable of simulating the behaviour of geomaterials in which eigentriads of stress and strain do not coincide during the loading process. The algorithm is tested in conjunction with the novel hyperelasto–plastic model termed the B material model, which is a single surface (single yield surface, affine single ultimate surface and affine single potential surface) model for dilatant, pressure sensitive, hardening and softening geomaterials. It is specifically developed to model large deformation hyperelasto–plastic problems in geomechanics. We present an application of this formulation to numerical analysis of low confinement tests on cohesionless granular soil specimens recently performed in a SPACEHAB module aboard the Space Shuttle during the STS‐89 mission. We compare numerical modelling with test results and show the significance of added confinement by the thin hyperelastic latex membrane undergoing large stretching. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
13.
In this paper, a coupling method between finite element and analytical layer‐elements is utilized to analyze the time‐dependent behavior of a plate of any shape and finite rigidity resting on layered saturated soils. Based on the integral transform techniques together with the aid of an order reduction method, an analytical layer‐element solution is derived from the governing equations for three‐dimensional Biot consolidation with respect to a Cartesian coordinate system and then extended to be the fundamental solution for the layered saturated soil under a point load. The Mindlin plate is modeled by eight‐noded isoparametric elements. The governing equations of the interaction between soil and plate in the Laplace‐Fourier transformed domain are deduced by referring to the coupling theory of FEM/BEM, and the final solution is obtained by applying numerical inversion. Numerical examples concerned with the time‐dependent response of a plate are performed to demonstrate the influence of soil and plate properties on the interaction process. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
14.
The reflection and transmission phenomena of elastic waves incident at a corrugated interface between two dissimilar fibre‐reinforced elastic half‐spaces have been analysed. Using Rayleigh method of approximation, the expressions of the reflection and transmission coefficients are obtained in closed form for the plane interface as well as for the first order approximation of the periodic interface ζ = d cos px. All these reflection and transmission coefficients of regular and irregular waves are found to be the functions of angle of incidence and elastic parameters of the media. Moreover, the coefficients of irregularly reflected and transmitted waves are found to be proportional to the amplitude of the corrugated interface and are functions of the frequency of the incident wave. Numerical computations have been performed for a specific model to compute these coefficients and results obtained are shown graphically. The results of Singh and Singh (Sadhana 2004; 29 :249–257) and Ben‐Menahem and Singh (Seismic Waves and Sources. Springer: New York) have been derived from our analysis as particular cases. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
15.
弹性地基板的分析简化模型 总被引：1，自引：0，他引：1
将弹性半空间地基受任意竖向荷载作用下的静力位移积分变换解与弹性半空间地基上四边自由矩形板受任意竖向荷载作用的弯曲解析解相结合,建立了求解板下地基位移的方法。对一些算例进行大量数值计算分析,得出弹性半空间地基上四边自由矩形板板下地基水平位移和竖向位移的分布规律;并基于该位移分布规律,提出地基位移沿深度按一定的函数关系变化的假设,考虑板下地基水平位移,利用板地基系统的总势能最小原理,通过复杂的变分运算,得出弹性地基板的简化模型。在不考虑板下地基水平位移时,该模型退化成双参数地基模型,并给出了求解其上四边自由矩形板的近似边界条件。 相似文献
16.
We derive the governing equations for the dynamic response of unsaturated poroelastic solids at finite strain. We obtain simplified governing equations from the complete coupled formulation by neglecting the material time derivative of the relative velocities and the advection terms of the pore fluids relative to the solid skeleton, leading to a so‐called us ? pw ? pa formulation. We impose the weak forms of the momentum and mass balance equations at the current configuration and implement the framework numerically using a mixed finite element formulation. We verify the proposed method through comparison with analytical solutions and experiments of quasi‐static processes. We use a neo‐Hookean hyperelastic constitutive model for the solid matrix and demonstrate, through numerical examples, the impact of large deformation on the dynamic response of unsaturated poroelastic solids under a variety of loading conditions. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
17.
Bending stresses produced in pavement concrete slabs are influenced by the contact condition along the slab–foundation interface. A finite element model for the analysis of frictional contact between unbonded pavement layers is presented. A 2D plate element is proposed for the idealization of composite slabs. Interface elements are used to connect the unbonded layers within the composite slab and to model the interface between the bottom slab layer and the foundation layer. Interface elements have the ability to capture the separation and sliding between pavement layers, due to thermal loads, and to calculate the frictional traction at their interface. 相似文献
18.
土与结构间接触面的非线性弹性-理想塑性模型及其应用 总被引：19，自引：3，他引：16
将非线性弹性理论与弹塑性理论相结合,对于土与结构间接触面提出了一种非线性弹性-理想塑性模型,用于模拟土与结构相互作用体系的变形与破坏机理。推导建立了这种接触面单元的应力-应变关系和弹塑性系数矩阵,并且讨论了这种模型在基坑开挖与支护分析应用中所面临的数值计算问题。最后针对某一基坑工程实例,应用这种接触面模型进行了数值计算与分析,结果表明,该模型能够较合理地模拟接触面上的变形机理与受力状态。 相似文献
19.
This paper presents a relatively simple method for three‐dimensional liquefaction analysis of granular soil under offshore foundations. In this method, the Mohr–Coulomb model, which defines the elasto–plastic stress–strain relationship under monotonic loading, is modified to accommodate the plastic strains generated by cyclic loading. The effects of cyclic loading, evaluated from the results of laboratory tests on saturated samples of soil, are incorporated into the model. The method is implemented in an efficient finite element program for analyses of three‐dimensional consolidating soil. The practicability of the model is demonstrated by analysis of a typical offshore foundation, and the predictions of the numerical analysis are compared with the observed behaviour of the foundation. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
20.
Fast multipole displacement discontinuity method (FM‐DDM) for geomechanics reservoir simulations 
下载免费PDF全文
下载免费PDF全文 The displacement discontinuity method (DDM) is frequently used in geothermal and petroleum applications for modeling the behavior of fractures in linear‐elastic rocks. The DDM requires O(N2) memory and O(N3) floating point operations (where N is the number of unknowns) to construct the coefficient matrix and solve the linear system of equations by direct methods. Therefore, the conventional implementation of the DDM is not computationally efficient for very large systems of cracks, often limiting its application to small‐scale problems. This work presents an approach for solving large‐scale fracture problems using the fast multipole method (FMM). The approach uses both the DDM and a kernel‐independent version of the FMM along with a preconditioned generalized minimal residual algorithm to accelerate the solution of linear systems of equations using desktop computers. Using the fundamental solutions for constant displacement discontinuity in a two‐dimensional elastic medium, several numerical examples involving fracture networks representing fractured reservoirs are treated. Numerical results show good agreement with analytical solutions and demonstrate the efficiency of the FMM implementation of the DDM for large‐scale simulations. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献