首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 672 毫秒
1.
王路君  艾智勇 《岩土力学》2018,39(9):3139-3146
利用解析层元法推导温度荷载作用下非稳态热传导时层状路面体系的温度响应解答。从热弹性理论平面应变问题的控制方程出发,借助于Laplace-Fourier积分变换,推导出单层介质及下卧半平面的精确刚度矩阵即解析层元,结合有限层法原理及边界条件,组装并求解总刚度矩阵,得到其在变换域内的解答,最后通过相应的积分逆变换得到物理域内的真实解。由于该法刚度矩阵元素中不含正指数项,计算时不会出现溢出或病态矩阵的现象。编译了相应的计算程序,所得结果与有限元模拟结果吻合较好。在此基础上,对有限深度和半平面两种假定条件下的解答进行对比分析,并分析层状路面体系中位移和温度随时间的变化趋势及沿深度的分布规律。分析表明:温度场具有一定的影响深度,超过此深度,有限深度与半平面理论解答基本一致;温度荷载的影响深度与其强度有关,强度越大,其影响深度越深。  相似文献   

2.
艾智勇  张逸帆 《岩土力学》2016,37(5):1243-1248
运用对偶积分方程来求解层状横观各向同性地基与墙下条形基础的共同作用问题。从直角坐标平面应变问题控制方程出发,通过傅里叶(Fourier)变换和层间连续性条件,可以得到层状横观各向同性地基的传递矩阵解。基于该传递矩阵解,并利用条形基础与地基接触的混合边值条件,推导出一组关于基础挠度和地基反力的对偶积分方程。考虑墙下条形基础受到竖向集中荷载的情况,利用弹性薄板理论先求解出条形基础挠度;随后应用雅可比(Jacobi)正交多项式和级数展开的方法,将对偶积分方程转化为线性代数方程组进行求解。编制了相应的计算程序,其计算结果与有限元软件ABAQUS的结果基本吻合,从而验证了所提理论的正确性。算例分析表明,板土相对刚度与地基成层性对地基反力、地表沉降和沿z轴竖向正应力有很大的影响。  相似文献   

3.
改进的双参数层状地基模型   总被引:8,自引:3,他引:5  
高广运  邱畅  王贻荪 《岩土力学》2003,24(2):159-163
文克尔地基模型广泛应用于弹性地基梁、板的设计计算,但由于该模型过于简化,不能正确反映地基土的工程性质。为此,提出了一种改进的双参数层状地基模型。该模型是Vlazov模型的延伸和发展,由一系列的弹性层组成,对每一层的应力应变分布做了一定的假设。通过积分变换的方法,可以求出该模型每一层表面位移与力的关系,进而形成层刚度矩阵。按照有限元法的原理,将每一层的刚度矩阵凝聚成总体刚度矩阵求解,与Vlazov模型计算结果的比较,证明该方法是正确的,还得到了荷载作用于层状地基模型内部的解,为双参数层状地基模型用于桩基分析打下了理论基础。  相似文献   

4.
艾智勇  慕金晶 《岩土力学》2018,39(7):2632-2638
基于Biot固结理论,运用解析层元方法求解竖向简谐荷载作用下二维层状饱和地基的动力响应问题。从直角坐标平面应变问题控制方程出发,通过Fourier-Laplace变换将偏微分方程组转化为常微分方程组,求解得到单层饱和地基的解析层元。结合层间连续条件和边界条件,组装得到多层饱和地基的总刚度矩阵方程,进而求得变换域内的解。借助Fourier-Laplace逆变换的数值积分方法,获得平面应变动力问题在物理域内的解,编制了相应的计算程序,其计算结果与已有文献结果吻合较好。通过算例分析了荷载圆频率、荷载作用深度及地基成层性对地基竖向位移的影响。计算结果表明:随荷载圆频率的增大,地基竖向位移先增加后减小;地基竖向位移在荷载作用点处呈现波峰,且受表层土性的影响较大。  相似文献   

5.
艾智勇  张逸帆  王路君 《岩土力学》2018,39(5):1885-1890
利用扩展精细积分法求解横观各向同性地基的平面应变问题。扩展精细积分法具有高精度和较高计算效率的特点,是求解微分方程的有效方法,相比于解析法可以节省大量理论推导工作量。从直角坐标下弹性力学控制方程出发,推导出Fourier变换域内地基的常微分矩阵方程;之后对地基微层元进行消元合并,进一步得到荷载作用在地基内部时层状地基的扩展精细积解。与已有文献的对比验证了方法的精确性,并分析了横观各向同性参数、层状性质和荷载作用点对计算结果的影响。结果表明:土体竖向位移随着横观各向同性参数n的增大而减小,而随着横观各向同性参数m的增大而增大;荷载作用点 的变化只对作用点以上的土体有影响,而上层土体的模量对竖向位移计算结果的影响更为显著,土体成层性对沉降的影响要比对竖向应力的影响更为显著。  相似文献   

6.
根据梁与地基之间的内力变形协调关系推导出考虑摩擦效应的文克尔弹性地基上杆系有单元法的单元刚度矩阵,计算分析表明,当地基比较坚 硬且人的接触较为粗糙是时,摩擦效应将会较显著地减小结的内力与变形。  相似文献   

7.
成层地基轴对称问题的基本解   总被引:1,自引:2,他引:1  
基于各向同性弹性体基本方程的转换及Hankel变换,得到了轴对称荷载作用下,无限弹性层轴向不同深度经Hankel变换的位移应力向量间的传递矩阵,运用该传递矩阵可求解成层的、层间完全接触情形和各向同性空间轴对称问题。  相似文献   

8.
艾智勇  曹国军 《岩土力学》2011,32(Z2):58-63
利用弹性矩形板与多层地基表面的竖向位移协调条件与光滑接触条件,由横观各向同性多层地基应力与位移非耦合的传递矩阵解,推导出弹性矩形板下竖向应力和位移的解析解。在此基础上,编制了相应的程序,并进行了数值计算。计算结果表明:矩形板刚度对板底竖向位移及板中心下的竖向应力有着较为显著的影响;板底竖向位移及板中心下的竖向应力随着板刚度的增加而减小,相同荷载作用下横观各向同性地基与均匀各向同性地基模型的计算结果差异较大,实际工程中很有必要采用更符合土体性质的横观各向同性地基模型  相似文献   

9.
对弹性地基上不同刚度的基础板,利用弹性地基上梁的计算方法进行了不同集中荷载工况下的变形计算。根据计算结果,分析了基础刚度对于弹性地基上基础的最大沉降的影响以及刚度对于基础整体变形的影响。  相似文献   

10.
赵强  焦玉勇  张秀丽  谢壁婷  王龙  黄刚海 《岩土力学》2019,40(11):4515-4522
非连续变形分析方法(DDA)是一种平行于有限元法的新型数值计算方法,该方法基于最小势能原理,把每个离散块体的变形、运动和块体之间的接触统一到平衡方程中进行隐式求解。然而,传统DDA方法在计算过程中需组装整体刚度矩阵并联立求解方程组,在用于大型岩土工程问题的三维数值模拟时占用内存较大、耗时较长、计算效率极低。因此,提出一种基于显式时间积分的三维球颗粒DDA方法。该方法在求解过程中不需要组装整体刚度矩阵,在求解加速度时,由于质量矩阵为对角矩阵,可存储为一维向量占用内存较少,且可分块逐自由度求解,效率较高,在接触判断上采用最大位移准则简化了接触算法,采用较小的时步,保证了计算的精确性;通过几个典型算例验证了该方法的准确性及计算效率。  相似文献   

11.
A boundary integral equation approach is presented for the two-dimensional plane strain analysis of horizontally layered elastic systems resting on a rough and rigid base. The validity of Somigliana identity for unbounded layers is discussed and an Infinite Boundary Element is proposed in order to take into account the lateral unboundedness of the physical problem in the discretized equations. This element is implemented in a computer code based on a successive stiffness solution procedure. A series of numerical tests concerning single- and multi-layered problems is illustrated, and the performance of the proposed solution method is compared with that of more traditional boundary element and finite element techniques.  相似文献   

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

13.
A method for the analysis of the consolidation of a horizontally layered soil under plane conditions is developed. The method depends upon the transformation of the governing equations by a Fourier trasform. This transformation has the effect of reducing the partial differential equations of consolidation to ordinary differential equations. The ordinary differential equations are then solved using a finite layer or finite difference approach. Once the solution in the transformed plane has been found, the actual solution is synthesized by Fourier inversion. The method leads to a considerable reduction in the amount of core storage necessary for solution and enables the solution of quite significant problems to be obtained on a mini-computer.  相似文献   

14.
This paper presents a model for the analysis of plane waves diffraction at a cavity in an infinite homogeneous poroelastic saturated medium, lined by a lining composed of four equal segments. An elastic boundary layer is placed between the cavity lining and the infinite porous medium. The boundary layer is simulated by ‘elastic boundary conditions’ in which the bulk matrix stress is proportional to the relative displacement between the lining and the surrounding medium matrix boundary. In addition, fluid impermeability through the intermediate layer is assumed. For the frequencies, that differ from the pseudoresonanse frequencies, the problem was reduced to the problem of an ideal elastic medium. A closed‐form analytical solution of the problem was obtained using Fourier–Bessel series, the convergence of which was proven. It was shown that the number of series terms required to obtain a desired level of accuracy can be determined in advance. The influence of the medium porosity on the medium dynamic stress concentration was studied. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

15.
裂隙型单斜介质中弹性系数的计算及波的传播特性研究   总被引:3,自引:1,他引:3  
根据Hudson等关于裂隙介质弹性系数计算的扰动理论及Bond变换矩阵原理,给出了各向同性介质中含多组垂直裂隙时等效弹性系数的计算方法。计算了含两组斜交的垂直裂隙形成的单斜各向异性介质中的等效弹性系数,并根据Christoffel方程推导、得出水平面内平面波传播的相速度和群速度随方位变化的特性。  相似文献   

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

17.
用双曲线数学模型描述被动土压力变化和挡土结构位移变化的关系。基于Mindlin弹性解,推导了土体某深度矩形区域作用均布水平荷载时该矩形区域内任意点水平位移的表达式,并用该表达式确定双曲线模型初始刚度。应用对数螺旋线被动土压力理论求解被动土压力,把该法求解的做图过程转化为数学迭代过程。算例表明,这样的被动土压力理论概念明确,结果准确。   相似文献   

18.
Measurements are conducted with small samples in the laboratory and thus for all practical purposes the medium is macroscopically homogeneous. On the other hand, the uncertainties and the irregular changes in situ are macroscopic inhomogeneities. This work is an attempt to account for these stochastic changes in the elastic properties and density in a rational manner. The method used is that of Karal and Keller which is based on the use of the Green's function and neglect of third-order correlations. The resulting integral equations are solved by Laplace transform. The analysis indicates that the energy decay in the mean motion through random mode coupling introduces damping into even a purley'elastic medium and enhances the damping in a significant manner in a hysteretic viscoelastic medium. This consideration is important in relating the damping and dispersion characteristics of wave in situ to those measured in the laboratory. The formulation is extended to multilayer systems through transfer matrices and to arbitrary inputs by Fourier transform. Sample calculations are presented for single and multilayer systems to obtain response spectra and for the response to Gaussian and actual earthquake input motions.  相似文献   

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

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