共查询到17条相似文献，搜索用时 70 毫秒
1.
数值流形方法是包含流形元、有限元及DDA在内的数值方法体系,建立流形元与DDA块体的接触方程,则可实现流形方法框架下的连续介质和散体系统共同作用模拟。针对填石路堤工程,编制了大型数值计算程序,采用块体随机生成、块体粒径控制及块体自然堆积的方法建立散体系统的DDA模型,对路堤的分层铺设、碾压及工后沉降变形等进行模拟分析。通过算例表明,在数值流形方法框架下,采用流形元与DDA共同作用的方法,可以很好地对同时存在连续变形和散体大变形的体系进行计算分析,其对该类问题的模拟更接近分析对象的实际情况,有助于从根本上揭示分析对象变形的细观机制和规律,并能考察更多因素对工程问题的影响。 相似文献
2.
非连续变形分析（DDA）是一种针对块体系统变形和位移求解的数值计算方法。引入Newmark方法于结构动力学微分方程中,考虑惯性力和阻尼力作用,改进时间步长自动调节,并实现DDA求解程序;比较研究Newmark方法中的线性加速法、常加速法和平均加速法在DDA程序中计算的收敛速度,讨论块体系统动力学计算过程中DDA方法对惯性力和阻尼力的添加和删除,并提出根据计算精度要求的误差控制实现方案。将改进的DDA方法模拟一个典型的煤与瓦斯突出过程,取得了满意的计算结果,该改进算法为DDA方法处理动力学问题提供新的途径 相似文献
3.
《岩土力学》2019,(12):4838-4847
动力时程分析中,结构地震反应的计算精度很大程度上取决于瑞利阻尼系数的确定。通过对某一成层场地在不同地震动作用下的动力反应分析,比较了目前工程上常用瑞利阻尼构建方法的时域计算结果和频域等效线性化计算结果,探讨了由频域分析中的滞回阻尼比构建时域分析中瑞利阻尼系数的新方法。在此基础上,对比研究了不同瑞利阻尼构造方法对某单层双跨地下车站结构地震反应的影响。结果表明:对于场地地震反应而言,提出的瑞利阻尼系数的时域计算结果与频域解的平均误差最小;对于地下结构地震反应而言,改进完整形式瑞利阻尼系数的时域计算结果与该方法计算结果基本一致。提出的瑞利阻尼系数计算方法精度较高,操作简单,可在场地和地下结构地震反应分析中推广使用。 相似文献
4.
非连续变形分析（DDA）是一种隐式求解的动力学计算方法,且采用在块体界面加减刚硬弹簧的方式来满足块体界面无张拉和无嵌入的接触准则,其中时间步长和弹簧刚度两个物理量的取值直接影响DDA的计算结果。基于对DDA时间步和弹簧刚度在程序运行过程中的调整策略和块体接触的简化力学概念模型,研究了惯性力在DDA收敛求解中的作用过程。采用数值模拟试验对自由落体和斜面单滑块模型在3种力学状态下的相关力学问题进行了数值模拟研究,通过对自由落体运动的模拟,研究了时间步长单一因素对计算结果的影响规律,并初步确定了时间步长的合理取值区间。在此基础上,采用斜面单滑块模型,研究了时间步长和弹簧刚度对计算结果的共同影响,确定了不同时间步长条件下弹簧刚度的合理取值区间。研究成果表明,合适的时间步长和弹簧刚度的取值组合构成一个单连通参数取值域,当时间步和弹簧刚度的取值组合位于此“域”范围内时,DDA的计算结果是合理的。 相似文献
5.
6.
为探究精密仪器厂房防微振结构设计方法,本文分别通过现场实测和数值模拟方法对苏州某高科技电子工业厂区一单层轻钢结构厂房微振动水平进行了评价和分析。对该厂区进行场地环境振动实测,对比了厂房建成前后地坪振动特性;利用ANSYS软件建立了厂房有限元模型,分析了Rayleigh阻尼和底板厚度对厂房振动特性的影响。研究结果表明:厂房建成后底板振动满足振动限值要求;厂房的桩筏基础有明显的减振效果,但对不同方向和频段地面的振动影响有所差异;阻尼比和底板厚度对结构主频影响较小,而对振幅影响较大,阻尼比在1~3 Hz和10~50 Hz的频段内对振幅影响较大,底板厚度在50 Hz以下对振幅影响较大。因此,在厂房防微振设计时,适当增加底板厚度和阻尼可减小中低频振动,从而达到较好的减振效果。 相似文献
7.
通过全面考虑土-结构动力相互作用体系中上部结构阻尼、地基土黏滞材料阻尼以及地基基础平动与转动时的辐射阻尼,推导了土-结构相互作用体系的耦合阻尼比公式。结合国内外关于阻尼调整系数的研究成果,提出了基于加速度反应谱的阻尼影响因子概念,建立阻尼影响因子 与《规范》中地震影响系数 的关系式。通过耦合阻尼比公式以及提出的 与 的关系式,对土-结构相互作用体系地震作用加以研究。结果表明考虑土与上部结构的相互作用后能够增大体系的阻尼性能,并大大降低上部结构的地震作用,提高体系的整体抗震性能。其成果为未来抗震防灾研究与工程的抗震设计提供了理论依据。 相似文献
8.
《岩土力学》2017,(7):2103-2112
针对非连续变形分析中开合迭代难以收敛的难题,基于块体接触约束状态和块体位移之间的关系,提出了基于逼近阶跃函数和拉格朗日插值的改进DDA方法。采用双曲正切函数来逼近阶跃函数,利用阶跃函数将块体接触约束状态用块体位移来表达,以此来替代开合迭代,避免了开合迭代难以收敛的难题。利用拉格朗日插值原理,推导得到只含有块体位移为未知量的块体系统势能函数,并利用变尺度法来求解总体势能函数的极值以得到块体位移。分别结合滑块模型和地下洞室模型,分析了改进DDA方法的计算精度和计算速度,验证了文中提出的改进DDA方法的正确性和稳定性。研究表明:基于逼近阶跃函数和拉格朗日插值的改进DDA方法具有较高的精度,且相比较传统DDA方法而言,具有更为稳定的和更为强健的计算收敛性。因此,基于逼近阶跃函数和拉格朗日插值的改进DDA方法是一种稳定有效的数值计算方法,为解决非连续变形中开合迭代难以收敛的问题提供了新思路。 相似文献
9.
10.
传统块体理论仅针对块体的静力稳定性分析,而未涉及到动力稳定分析。因此,借用Newmark法累积滑移量的概念,基于块体理论对定位块体在地震历时作用下的动力稳定性进行了研究,提出了以块体累积滑动位移和动安全系数两个指标共同评价块体动力稳定性。提出一种接触面识别及面积计算算法,它考虑了接触面积变化对块体稳定的不利影响,也可用于块体滑动的大变形分析,考虑了块体滑动过程中结构面抗剪强度的演化,可对复杂的凸形块体及凹形块体进行静、动力稳定性分析。将其应用于三峡升船机乘船箱左侧边坡的块体稳定分析中,结果显示块体在强震作用下发生了一段距离的滑动,但震后仍可保持稳定 相似文献
11.
Qinghui Jiang Yifeng Chen Chuangbing Zhou Man-chu Ronald Yeung 《Rock Mechanics and Rock Engineering》2013,46(6):1443-1460
The discontinuous deformation analysis (DDA) is a numerical method for modeling discontinuous deformation behaviour of jointed rocks. In this paper, two basic problems are discussed related to kinetic energy dissipation and the convergence criterion for the DDA method when it is applied to geotechnical engineering. In view of the fact that the deformation and progressive failure can be treated as a quasi-static process with low kinetic energy dissipation rates, this paper introduces a viscous damping component to absorb discrete blocks’ kinetic energy, establishes the global equations of motion of the discrete block system that take damping effects into account, investigates the energy dissipation mechanism when solving a static or quasi-static problem, and defines the convergence criteria of displacement, kinetic energy and unbalanced force for DDA solutions when the system arrives at a stable state. 相似文献
12.
13.
S. Amir Reza Beyabanaki Ahmad Jafari S. Omid Reza Biabanaki M. Ronald Yeung 《Computers and Geotechnics》2009
This paper presents a new numerical model that can add a finite element mesh into each block of the three-dimensional discontinuous deformation analysis (3-D DDA), originally developed by Gen-hua Shi. The main objectives of this research are to enhance DDA block’s deformability. Formulations of stiffness and force matrices in 3-D DDA with conventional Trilinear (8-node) and Serendipity (20-node) hexahedral isoparametric finite elements meshed block system due to elastic stress, initial stress, point load, body force, displacement constraints, inertia force, normal and shear contact forces are derived in detail for program coding. The program code for the Trilinear and Serendipity hexahedron elements have been developed, and it has been applied to some examples to show the advantages achieved when finite element is associated with 3-D DDA to handle problems under large displacements and deformations. Results calculated for the same models by use of the original 3-D DDA are far from the theoretical solutions while the results of new numerical model are quite good in agreement with theoretical solutions; however, for the Trilinear elements, more number of elements are needed. 相似文献
14.
A New Viscous Boundary Condition in the Two-Dimensional Discontinuous Deformation Analysis Method for Wave Propagation Problems 总被引：1，自引：1，他引：0
Viscous boundaries are widely used in numerical simulations of wave propagation problems in rock mechanics and rock engineering. By using such boundaries, reflected waves from artificial boundaries can be eliminated; therefore, an infinite domain can be modeled as a finite domain more effectively and with a much greater accuracy. Little progress has been made, thus far, with the implementation and verification of a viscous boundary in the numerical, discrete element, discontinuous deformation analysis (DDA) method. We present in this paper a new viscous boundary condition for DDA with a higher absorbing efficiency in comparison to previously published solutions. The theoretical derivation of the new viscous boundary condition for DDA is presented in detail, starting from first principles. The accuracy of the new boundary condition is verified using a series of numerical benchmark tests. We show that the new viscous boundary condition works well with both P waves as well as S waves. 相似文献
15.
在动三轴试验过程中,加载频率对动模量和阻尼比的影响目前还没有定论。基于室内动三轴试验结果,采用二维颗粒流程序PFC2D模拟了筑坝反滤料在不同加载频率（f = 1.0、0.5、0.1、0.05 Hz）条件下的室内动三轴试验。数值模拟结果与试验结果基本吻合,试验结果表明,（1）当 f ≥ 0.1 Hz时,频率对动模量和阻尼比的影响可以忽略;当 f < 0.1 Hz时,模量在小应变情况下会有所降低,但随着应变的增大会趋于一致。（2）阻尼比在f = 0.05 Hz时波动比较大,但与其他频率条件下的阻尼比相比都在一个较小的区域内,可以认为频率对阻尼比基本没有影响 相似文献
16.
《Geomechanics and Geoengineering》2013,8(3):221-236
This paper presents a new point-to-face contact algorithm for contacts between two polyhedrons with planar boundaries. A new discrete numerical method called three-dimensional discontinuous deformation analysis (3-D DDA) is used and formulations of normal contact submatrices based on the proposed algorithm are derived. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. This approach does not need to use an iterative algorithm in each time step to obtain the contact plane, unlike the ‘Common-Plane’ method applied in the existing 3-D DDA. In the present 3-D DDA method, block contact constraints are enforced using the penalty method. This approach is quite simple, but may lead to inaccuracies that may be large for small values of the penalty number. The penalty method also creates block contact overlap, which violates the physical constraints of the problem. These limitations are overcome by using the augmented Lagrangian method that is used for normal contacts in this research. This point-to-face contact model has been programmed and some illustrative examples are provided to demonstrate the new contact rule between two blocks. A comparison between results obtained by using the augmented Lagrangian method and the penalty method is presented as well. 相似文献
17.
The complexity of determining strain associated with shear modulus and damping ratio in torsional tests has been resolved by means of several approaches. The stress integration approach is adequate when generating the plots of equivalent radius ratio versus strain more effectively over any range of strains in resonant column and torsional shear (RC/TS) tests. The stress integration approach was applied for hyperbolic, modified hyperbolic, and Ramberg–Osgood models in evaluating damping ratio. This study showed that using a single value of equivalent radius ratio in evaluating damping ratio is not appropriate. The combined hysteretic‐nonviscous damping model was developed and employed to consider the increased damping behavior at small strains using the stress integration approach. The results suggest that adding viscous behavior has no significant effect strain calculations in RC/TS testing. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献