共查询到20条相似文献，搜索用时 125 毫秒
1.
提出了砂土中水平受荷单桩的修正应变楔方法。从3个方面对应变楔模型进行了修正,首先假设三维被动土楔后桩身的水平位移呈非线性变化,从而使得楔形体内的水平应变不再是一常量,而是沿深度变化。其次引入了两个双曲线型的模型分别用以模拟楔形体中土体的水平应力增量-应变关系和桩-土界面处桩侧剪应力-位移关系。然后通过算例验证了修正方法的有效性,模拟结果与实测结果吻合较好。最后分析了各个修正因素对计算结果影响,结果显示,非线性位移假设的引入对原应变楔模型计算的水平承载力偏高有明显地改善,新的水平应力增量-应变关系和桩侧剪应力-位移关系的引入使得应变楔方法更加简便有效。 相似文献
2.
非连续变形分析（DDA）是一种针对块体系统变形和位移求解的数值计算方法。引入Newmark方法于结构动力学微分方程中,考虑惯性力和阻尼力作用,改进时间步长自动调节,并实现DDA求解程序;比较研究Newmark方法中的线性加速法、常加速法和平均加速法在DDA程序中计算的收敛速度,讨论块体系统动力学计算过程中DDA方法对惯性力和阻尼力的添加和删除,并提出根据计算精度要求的误差控制实现方案。将改进的DDA方法模拟一个典型的煤与瓦斯突出过程,取得了满意的计算结果,该改进算法为DDA方法处理动力学问题提供新的途径 相似文献
3.
4.
任意形状块体的接触模拟是发展三维不连续变形分析方法（3D-DDA）的关键问题。提出了三维离散块体单元的 边-边接触模型,并给出了边-边接触进入面的判定方法和嵌入准则。基于矢量分析和罚函数法,推导了接触子矩阵的表达式,算例分析结果表明了该模型的有效性。 相似文献
5.
真空预压加固软土地基变形与固结计算研究 总被引：1，自引：0，他引：1
考虑真空度的衰减情况,对真空预压加固软基的变形和固结度计算方法进行了研究。先根据真空预压时砂井的真空压力状态建立了真空预压的空间轴对称变形模型,用位移法推导该模型在等应变条件下的变形及体积应变计算公式。在此基础之上结合Hansbo砂井地基固结理论和真空预压的边界条件,推导了忽略竖向渗流情况下的真空预压加固软土地基的固结解析解。比较文中计算方法与已有计算理论和现场试验实测资料结果表明,体积应变的计算对真空预压的孔压和固结度的计算有较大的影响,而直接引用传统堆载预压的体积应变计算公式计算会导致较大误差。在固结度计算中,采用Hansbo的近似方法能满足计算精度的要求,所得计算结果与实测结果吻合较好 相似文献
6.
《岩土力学》2020,(2)
详细地介绍了基于修正对称和反对称分解(MSAD)的三维数值流形元法,并提出一个针对三维数值流形元法(3D NMM)中应用罚函数法施加位移约束和材料边界条件时罚系数的选取公式。在基于MSAD的三维数值流形元法中,引入了Bathe隐式时间积分方案,编写了基于统一强度理论和非关联流动法则的理想塑性本构模型,实现了三维弹塑性开挖问题模拟。将基于MSAD的三维数值流形元法应用到非线性动力学研究中,案例研究结果表明:Bathe隐式时间积分方案和基于MSAD的三维数值流形元法在处理大转动和长持续时间的非线性动力学问题时能够很好地保障模拟结果的稳定性,同时保证守恒体系动能和角动量的守恒。再次验证了MSAD理论,在模拟大转动问题时,MSAD具有很好的稳定性和较高的计算精度,能够合理地从变形梯度增量中分离出转动和应变,精确地更新转动应力,而不会产生错误体积膨胀问题。 相似文献
7.
《岩土力学》2021,(2)
详细地介绍了基于修正对称和反对称分解(MSAD)的三维数值流形元法,并提出一个针对三维数值流形元法中应用罚函数法施加位移约束和材料边界条件时罚系数的选取公式。在基于MSAD的三维数值流形元法中,引入了Bathe隐式时间积分方案,编写了基于统一强度理论和非关联流动法则的理想塑性本构模型,实现了三维弹塑性开挖问题模拟。将基于MSAD的三维数值流形元法应用到非线性动力学研究中,案例研究结果表明:Bathe隐式时间积分方案和基于MSAD的三维数值流形元法在处理大转动和长持续时间的非线性动力学问题时能够很好地保障模拟结果的稳定性,同时保证守恒体系动能和角动量的守恒。再次验证了MSAD理论,在模拟大转动问题时,MSAD具有很好的稳定性和较高的计算精度,能够合理地从变形梯度增量中分离出转动和应变,精确地更新转动应力,而不会产生错误体积膨胀问题。 相似文献
8.
地震永久变形是评价土石坝抗震稳定性的一个重要依据,其合理计算需要同时考虑材料的残余体积应变和残余剪切应变的影响。因此,在过去提出的永久变形等效节点力计算方法的基础上,通过考虑剪切压缩耦合作用来反映残余体积应变对地震永久变形的影响,并推导了相应的模量矩阵,明确了参数选择和计算步骤。对典型单元体的残余变形模式进行了分析,计算得到的剪切应变和体积应变大小均与理论值相符,验证了模量矩阵的可靠性。在此基础上,对心墙堆石坝和面板堆石坝的算例进行了计算,重点分析了竖直向、顺河向和横河向3个方向永久变形分布,并与实测规律进行比较。计算结果表明,考虑残余体应变之后,所提方法能够反映土石坝地震永久变形以沉降为主、水平变形相对较小的特点,可更为合理地反映地震永久变形的规律。 相似文献
9.
在基坑设计计算中,就深基坑桩锚支护结构的受力与变形研究而言,目前运用较为广泛的是改进增量法理论,荷载增量的准确计算是该理论是否应用于实践的关键点,荷载增量的计算中是否应考虑开挖面以下土体因抗力系数变动引起弹性抗力变化释放的反力增量成为研究焦点,而对土体弹性抗力的影响因素却少有研究。鉴于此,本文在计算荷载增量时,综合考虑位移土压力及开挖面以下土体因抗力系数变动引起弹性抗力变化的影响,提出了更符合工程实际的增量计算方法。结合工程实例,编制了MATLAB计算程序,将其所得结果与未考虑前进行了对比,并进一步采用有限元软件PLAXIS对工程实例进行模拟,验证了此方法的合理性,结果表明本文提出的基于改进增量法的桩锚支护结构位移与内力计算方法是合理、可行的。最后对影响土体弹性抗力的两大因素——土体黏聚力c与内摩擦角φ,通过控制变量法分别进行了讨论分析,结果显示随着土体黏聚力c或内摩擦角φ的减小,桩身位移呈现逐渐增大的趋势,但最大位移都在基坑变形允许的范围内。 相似文献
10.
《岩土力学》2017,(7):2103-2112
针对非连续变形分析中开合迭代难以收敛的难题,基于块体接触约束状态和块体位移之间的关系,提出了基于逼近阶跃函数和拉格朗日插值的改进DDA方法。采用双曲正切函数来逼近阶跃函数,利用阶跃函数将块体接触约束状态用块体位移来表达,以此来替代开合迭代,避免了开合迭代难以收敛的难题。利用拉格朗日插值原理,推导得到只含有块体位移为未知量的块体系统势能函数,并利用变尺度法来求解总体势能函数的极值以得到块体位移。分别结合滑块模型和地下洞室模型,分析了改进DDA方法的计算精度和计算速度,验证了文中提出的改进DDA方法的正确性和稳定性。研究表明:基于逼近阶跃函数和拉格朗日插值的改进DDA方法具有较高的精度,且相比较传统DDA方法而言,具有更为稳定的和更为强健的计算收敛性。因此,基于逼近阶跃函数和拉格朗日插值的改进DDA方法是一种稳定有效的数值计算方法,为解决非连续变形中开合迭代难以收敛的问题提供了新思路。 相似文献
11.
Updating the block configuration on the basis of additive decomposition and its linearized expression of the displacement increment leads to the low calculation accuracy of the original discontinuous deformation analysis (DDA) and false volume expansion. In this study, the displacement expressions of a small deformation, a large rotation, and the corresponding velocity and acceleration terms on the basis of the initial configuration are presented using multiplicative decomposition. With the use of the principle of virtual work, the stiffness matrix, mass matrix, and force vector of blocks are obtained. Compared with the original DDA, each of the block deformation parameters has obvious physical meaning as a parameter of mechanics, which can be obtained by adding the incremental deformation components of each time step directly without co-ordinate transformation. Moreover, the proposed modification automatically considers the block deformation produced by centrifugal and Coriolis forces. The analysis of some typical numerical examples have verified the accuracy of the strain and stress calculated by the proposed method, and the current configuration is updated by the total displacements, which completely overcomes the false volume expansion and provides reasonable linear strains. 相似文献
12.
Conventional discontinuous deformation analysis (DDA) results in a change of block volume, which is known as free expansion, during rotation calculations because of the use of a linear displacement function to simulate the behavior of a block with a rigid body and elastic behaviors. This study demonstrates that the linear displacement function also generates unsolved elastic distortion, especially when the block undergoes large rotation in each calculation step. The distortion disturbs the contact judgment in the open–close iteration and update calculations of vertex coordinates, stresses, velocities, etc. at the end of each calculation step. A new procedure follows the flow chart of the original DDA, but it adopts additional codes for the coordinate-transformation calculations in vertex coordinate, stress, and velocity updates. When the vertex coordinates are updated, vertex displacements caused by strains are calculated before involving the block-rotation term in the displacement function to mitigate the elastic distortion. In addition, new codes compile formulas to transform stresses and velocities with block rotation. The new DDA ensures the correctness of rotating elastic calculations to solve practical falling rock problems with a large rotational angle in each calculation step. 相似文献
13.
An extensive examination of the discontinuous deformation analysis (DDA) in block dynamic sliding modeling is carried out in this paper. Theoretical solutions for a single block sliding on an arbitrarily inclined plane by applying the horizontal/vertical seismic loadings to the sliding block as acceleration time histories or to the base as constraint displacement time histories are derived. As compared with the theoretical solutions, for a single block sliding, the DDA predicts the sliding displacements and block interaction forces accurately under various base incline angles and friction angles under both the harmonic loadings and a real seismic loading. The vertical seismic component may influence the block sliding displacements to different extent, and the DDA can capture these phenomena successfully and give accurate results. For the calculation of the single block relative sliding, both the theoretical and the DDA solutions indicate that applying the seismic accelerations as constraint displacement time histories (derived by integrating the seismic accelerations twice) to the base is equivalent to applying the seismic accelerations as volume forces to the sliding block in the opposite directions. The DDA modeling also demonstrates that this conclusion still stands for the case of multi‐block sliding. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
14.
A Boundary Element based Discontinuous Deformation Analysis (BE‐DDA) method is developed by implementing the improved dual reciprocity boundary element method into the open close iterations based DDA. This newly developed BE‐DDA is capable of simulating both the deformation and movement of blocks in a blocky system. Based on geometry updating, it adopts an incremental dynamic formulation taking into consideration initial stresses and dealing with external concentrated and contact forces conveniently. The boundaries of each block in the discrete blocky system are discretized with boundary elements while the domain of each block is divided into internal cells only for the integration of the domain integral of the initial stress term. The contact forces among blocks are treated as concentrated forces and the open–close iterations are applied to ensure the computational accuracy of block interactions. In the current method, an implicit time integration scheme is adopted for numerical stability. Three examples are used to show the effectiveness of the algorithm in simulating block movement, sliding, deformation and interaction of blocks. At last, block toppling and tunnel stability examples are conducted to demonstrate that the BE‐DDA is applicable for simulation of blocky systems. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
15.
Discontinuous deformation analysis (DDA), a discrete numerical analysis method, is used to simulate the behaviour of falling rock by applying a linear displacement function in the computations. However, when a block rotates, this linear function causes a change in block size called the free expansion phenomenon. In addition, this free expansion results in contact identification problems when the rotating blocks are close to each other. To solve this problem of misjudgment and to obtain a more precise simulation of the falling rock, a new method called Post‐Contact Adjustment Method has been developed and applied to the program. The basic procedure of this new method can be divided into three stages: using the linear displacement function to generate the global matrix, introducing the non‐linear displacement function to the contact identification, and applying it to update the co‐ordinates of block vertices. This new method can be easily applied to the original DDA program, demonstrating better contact identification and size conservation results for falling rock problems than the original program. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
16.
锚索抗滑桩的设计计算包括锚索拉力计算和桩体内力计算,而锚索拉力又直接影响桩体内力,因此,锚索拉力计算对桩-锚结构设计十分关键。在常规和修改变形协调条件的基础上,考虑锚索预应力作用在滑动面处桩身产生的位移和转角对第i排锚索作用点处桩身位移的贡献,除去其他排锚索预应力Rj0对桩身作用产生的位移,同时考虑锚索拉力增量的水平分量作用对桩体变形的影响,对滑动面处桩身位移、转角以及其他排锚索拉力对某一排锚索作用点处桩身位移的贡献进行优化,建立优化的桩-锚变形协调方程。结合结构力学中的虚功原理和图乘法理论,编写Visual C++计算程序,通过滑坡工程实例对锚索预应力和总拉力进行计算。3种变形协调条件的计算结果表明：优化的变形协调条件不仅能使抗滑桩在预应力施加阶段处于有利的受力和变形状态,而且适用于锚索与水平面成不同角度的情况,弥补了常规和修改变形协调条件的不足。 相似文献
17.
传统的非连续变形分析法（DDA）法采用简单的线性位移模式计算效率高,描述大块体的高阶多项式位移模式在一定程度保留了该特点,并提高了计算精度。近年来流行的耦合有限元、自然单元的DDA法实质上是引入相应的插值形函数构成块体位移函数,计算相对低效,但具有计算更精细、更容易施加边界条件等优点。为结合传统DDA法与DDA耦合法各自的优点,建立了一种同时利用传统DDA法线性位移模式与耦合型DDA法非线性位移模式的混合法。该方法非线性模式主要针对大块体,采用了自然单元插值,缘于其具有一定无网格特征,且效率比有限元高。建立了混合模式下的整体矩阵并推导出接触等因素刚度子矩阵和荷载子向量的具体表达式。该方法建模更加方便合理,计算精度、效率介于线性模式的传统DDA法和非线性位移模式的耦合法之间。通过基本算例验证了混合法的有效性,并给出了节理围岩-隧道衬砌整体分析模型的计算结果,体现了新方法的优越性。 相似文献
18.
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. 相似文献
19.
The high computational costs associated with the implicit formulation of discontinuous deformation analysis (DDA) have been one of the major obstacles for its implementation to engineering problems involving jointed rock masses with large numbers of blocks. In this paper, the Newmark-based predictor-corrector solution (NPC) approach was modified to improve the performance of the original DDA solution module in modeling discontinuous problems. The equation of motion for a discrete block system is first established with emphasis on the consideration of contact constraints. A family of modified Newmark-based predictor-corrector integration (MNPC) scheme is then proposed and implemented into a unified analysis framework. Comparisons are made between the proposed approach and the widely used constant acceleration (CA) integration approach and central difference (CD) approach, regarding the stability and numerical damping features for a single-degree-of-freedom model, where the implications of the proposed approach on open-close iteration are also discussed. The validity of the proposed approach is verified by several benchmarking examples, and it is then applied to two typical problems with different numbers of blocks. The results show that the original CA approach in DDA is efficient for the simulation of quasi-static deformation of jointed rock masses, while the proposed MNPC approach leads to improved computational efficiency for dynamic analysis of large-scale jointed rock masses. The MNPC approach therefore provides an additional option for efficient DDA of jointed rock masses. 相似文献
20.
David M. Doolin 《国际地质力学数值与分析法杂志》2005,29(12):1199-1207
Displacement boundary constraints in discontinuous deformation analysis (DDA) are applied using stiff penalty springs. A co‐ordinate‐free formulation for displacement boundary constraints is presented here for DDA, which unifies previous derivations for points of fixity, and for points constrained to induce or prohibit block motion in specified directions as a function of location or time. Examples for each type of constraint are used to illustrate the behaviour of the algorithm and provide a link with previous formulations for each case. The new, unified formulation has five benefits: (1) simple to express algorithmically; (2) easy to program and verify; (3) penalty values in different directions may be chosen to allow fixed points, lines, curves or planes; (4) formulation works for 2D and 3D; (5) displacement constraint may be a function of time or location or both. Feedback in the algorithm may induce internal resonance in homogeneously deformable discrete elements used in DDA, and resonance in block‐to‐block contact interactions. Consequently, high mass problems with insufficient damping may suffer from excessive ‘vibrational hammering’, inducing physically implausible behaviour such as elastic rebound. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献