双向粘弹性节理岩体时域等效分析

以双向粘弹性节理岩体为研究对象,基于一种简单等效假定和自适应算法,并结合有限元方法,得到递推格式的等效本构模型和等效粘弹性场计算模型.通过算例,将等效解与基于ANSYS的非均质解进行比较.从计算精度和效率平衡的角度而言,等效模型的结果是令人满意的.     

基于时域自适应算法的单向粘弹性节理岩体的等效分析

在粘弹性节理岩体的数值仿真中,必需合理地计及节理的影响.若将节理与岩体作为独立的材料组分考虑,在节理密集的情况下,时空两方面的计算开销可能是很难承受的.本文通过时域展开技术与一个简单的等效假定,得到了递推格式的粘弹性单向节理岩体的等效本构关系,及等效场的递推求解格式,并由此数值模拟了一个粘弹性节理岩体中的洞室,计算结果与考虑独立组分的ANSYS粘弹性分析进行了比较,二者的计算平均误差约为11.03%,但前者的计算效率为后者16.81倍.本文工作有可能为粘弹性节理岩体的数值模拟提供一条新的途径.     

秦山三期核电工程反应堆地段节理裂隙模拟

以秦山三期核电工程反应堆地段作为研究对象 ,针对节理裂隙在岩体中的分布具有随机性的特点 ,运用随机方法对其分布规律进行了研究。根据实测的节理裂隙产状、间距及迹线长度 ,推断节理裂隙概率分布特征。同时充分考虑到节理裂隙分布的不均匀性和方向性等特点 ,采用非平稳态随机过程模拟节理裂隙间距 ,通过 Monte- Carlo方法得到等效的节理裂隙网络。根据统计和拟合的结果 ,对反应堆地段岩体的节理裂隙进行了评价。     

Jennings边坡稳定性分析的概率方法(英文)

Jennings边坡 (简称J边坡 )是一种常见于格状与层状岩体中的典型岩质边坡模型 ,论文在岩质边坡中对J边坡运用蒙特卡罗方法进行概率分析 ,考虑陡、缓倾角两组节理随机组合 ,对计算范围内节理分布进行位置搜索 ,并考虑水压力及地震灾害条件下 ,计算相应坡坏概率及其可靠度指标。使用该方法对江西某铜矿高边坡工程进行了计算 ,取得令人满意的结果     

岩体裂隙扩展过程的数值模拟

本文根据断裂力学理论及节理岩体的等效连续模型,探讨了岩体裂隙扩展过程的数值模拟方法,文末算例说明了方法的有效性。     

深部节理岩体塑性损伤耦合微面模型

采用微面模型理论和损伤力学方法,建立了节理岩体的弹塑性损伤耦合微面模型. 在节理岩体的微面上,将岩体视为由节理面与岩石组成的二元介质,以节理连通率作为岩体沿该方向的面积损伤变量,考虑微面法向拉应力和压应力下的不同塑性变形和损伤耦合作用机制,基于塑性理论建立了节理岩体的微面塑性损伤增量本构关系. 采用微面物理量与宏观物理量的几何约束模型,根据微面方向积分导出了节理岩体的宏观弹塑性增量本构关系. 编制了节理岩体微面模型的MARC有限元子程序,对节理岩体的单轴拉伸、压缩试验和泥浆压力作用下的井壁稳定问题进行了数值模拟研究. 数值计算结果表明,该模型能很好地揭示载荷作用下节理岩体的各向异性非弹性变形和次生节理演化过程.      

隧道围岩强度不均地段塌方成因及其处理方法

采用微面模型理论和损伤力学方法,建立了节理岩体的弹塑性损伤耦合微面模型. 在节理岩体的微 面上,将岩体视为由节理面与岩石组成的二元介质,以节理连通率作为岩体沿该方向的面积 损伤变量,考虑微面法向拉应力和压应力下的不同塑性变形和损伤耦合作用机制,基于塑性 理论建立了节理岩体的微面塑性损伤增量本构关系. 采用微面物理量与宏观物理量的几何约 束模型,根据微面方向积分导出了节理岩体的宏观弹塑性增量本构关系. 编制了节理岩体微 面模型的MARC有限元子程序,对节理岩体的单轴拉伸、压缩试验和泥浆压力作用下的井壁稳 定问题进行了数值模拟研究. 数值计算结果表明,该模型能很好地揭示载荷作用下节理岩体 的各向异性非弹性变形和次生节理演化过程.     

裂隙岩体损伤局部化破坏分岔模型及其应用

采用概率统计方法分析节理裂隙岩体的几何特征，定义了反映裂隙岩体几何特征的组构张量．基于不可逆热力学理论，通过裂纹扩展的细观分析，得出了损伤的发展机理和演化方程，把损伤演化和裂隙的几何特征的变化联系起来，建立了弹塑性损伤本构关系．为分析含有节理裂隙岩体在发生局部化破坏时的特征，通过对发生局部化时的裂隙岩体的分析。构造了适用于节理裂隙岩体局部化分析的不连续分岔模型．利用非线性规划数值解法，可以得出局部化破坏的方向特征．在有限元方法中，根据该模型给出了节理裂隙岩体相关的算例，分析表明该模型用于分析裂隙岩体的局部化破坏是有效的．     

层状岩体力学特性模拟实验研究

对层状岩体变形特性和强度特性的正确认识,合理选取其变形和强度参数,是水电站等大型岩体工程论证与设计需要解决的关键问题.为此,利用模拟实验对层状岩体的变形特性和强度特性进行了研究,得到了模拟试样的变形模量和强度随试样所含节理数、轴向应力与节理面间的夹角和围压之间的变化关系.结果表明:当最大主应力与节理面间的夹角由0°增加...     

层状盐岩力学特性研究进展

盐岩具有孔隙度低,渗透率小,损伤自恢复和塑性变形能力大等特性.本文着眼于深部盐岩能源地下储备,在简要介绍国内外盐岩力学特性研究现状基础上,着重综述了近期关于层状沉积型盐岩工程力学特性研究的进展.层状盐岩压缩试验表明泥质硬石膏夹层对盐岩体的变形和破坏特性有明显的影响,揭示了复合岩体的变形和破坏机制;层状盐岩流变试验表明层状盐岩的蠕变主要由盐岩层控制,硬夹层的存在对溶腔围岩的长期蠕变起到抑制作用;为了深入了解界面的力学特性,开展了层状盐岩直接剪切试验、间接拉伸试验和界面微观分析试验,试验表明层状盐岩体中盐岩夹层界面具有较强的黏结力,这十分有利于盐岩溶腔的密闭性和稳定性;针对特殊的层状岩体提出了复合岩体Cosserat介质扩展本构模型,并应用于层状盐岩构造中能源储库稳定性分析中.最后指出了层状盐岩工程力学特性研究中值得关注的几个重点研究方向.      

用均匀化理论分析蜂窝结构的等效弹性参数

在线弹性范围内,根据均匀化理论,并结合有限元方法推导出适用于二维周期性结构的均匀化的有限元格式（Homo FEM）,计算出不同相对密度下的规则蜂窝结构的等效弹性模量Ee和泊松比νe.同时,利用蜂窝结构的代表胞元模型,用常规的有限元方法（FEM）模拟计算出相应的等效弹性参数.最后将两种数值计算结果与己有的理论公式进行了比较和分析讨论.结果表明:在考察的相对密度全场范围内（0～0.4）,HOmO FEM得到的蜂窝结构的 Ee和νe 与 FEM使用平面实体单元模拟计算得到的结果一致吻合,反映出 Homo FEM数值方法的客观准确性和可行性.而 Gibson公式和 W-K得到的等效弹性模量值 Ee只是在较小相对密度的情况下（小于0.15）与数值计算结果吻合.当结构的相对密度较大时,必须考虑胞棱附近区域由应力集中导致的复杂的应力和应变分布的影响.     

基于迭代法的非线性弹性均质化研究

微观结构对复合材料的宏观力学性能具有至关重要的影响, 通过合理设计复合材料微观结构可以得到期望的宏观性能. 均质化方法作为一种有效的设计方法, 它从微观结构的角度出发, 利用均匀化的概念, 实现了对复合材料宏观力学性能的预测和设计. 而当考虑非线性因素, 均质化的实现就非常困难. 本文利用双渐近展开方法, 将位移按照宏观位移和微观位移展开, 推导了非线性弹性均质化方程. 通过直接迭代法, 对非线性弹性均质化方程进行了求解, 并给出了具体的迭代方法和实现步骤. 本文基于迭代步骤和非线性弹性均质化方程编写MATLAB 程序, 对3种典型本构关系的周期性多孔材料平面问题进行了计算, 对比细致模型的应变能、最大位移和等效泊松比, 对程序及迭代方法的准确性进行了验证. 之后对一种三元橡胶基复合材料进行多尺度均质化, 将其分为芯丝尺度和层间尺度. 用线弹性的均质化方法得到了芯丝尺度的等效弹性参数, 并将其作为层间尺度的材料参数. 在层间尺度应用非线性弹性均质化方法对结构进行计算, 得到材料的宏观等效性能, 并以实验结果为基准进行评价.      

Modeling Free-Surface Seepage Flow in Complicated Fractured Rock Mass Using a Coupled RPIM-FEM Method

The prediction of the free-surface seepage flow behavior in fractured rock mass is of significance in geotechnical engineering. There are two major issues in solving the seepage flow in complicated fractured rock mass based on the fractured porous medium (FPM) flow model, in which groundwater is assumed to flow simultaneously in both rock matrix and embedded fractures: One is the mesh generation of rock mass in the presence of the fracture network, especially when there exist a large number of stochastic fractures; the other is that a robust iteration algorithm is required since the free surface is unknown at the beginning of solution. Aiming at these two issues, this paper proposes a novel numerical method by coupling radial point interpolation method (RPIM) and finite element method (FEM), in which RPIM is utilized to model the rock matrix and FEM is utilized to model the fractures. On the basis of the variational inequality (VI) theory for free-surface seepage analysis, the computation formulations of the numerical method are derived and the corresponding computation program is developed. Three examples are solved with the present method. It is found that the VI theory can be extended to solve the free-surface seepage problem based on the FPM flow model. A crucial advantage of the present method is that the mesh generation can be greatly simplified. The present method has been verified to be a robust, efficient and reliable method for modeling the groundwater flow in complicated fractured rock mass.     

分数阶流变岩体中深埋圆形支护隧道力学分析

为用更少的参数描述不同岩体的实际流变特性,引入分数阶Kelvin粘弹模型并对深埋圆形支护隧道的开挖过程进行力学分析．利用Laplace变换原理、分数阶微积分得到了衬砌支护力和洞口位移的时效解答,并对解答的正确性进行了验证．讨论了不同分数阶阶次下支护力、洞口位移及等效应力随时间变化的规律．分析结果表明,阶次较小的模型,洞口位移在短期内迅速增大,随后缓慢增长,在经历很长时间后仍在持续增长,且等效应力在相同时刻下更接近于0,更容易发生破坏,因此在工程中对阶次较小模型对应的岩体要关注其长时间后的蠕变位移和受力状态．     

Cosserat model for periodic masonry deduced by nonlinear homogenization

The paper deals with the problem of the determination of the in-plane behavior of periodic masonry material. The macromechanical equivalent Cosserat medium, which naturally accounts for the absolute size of the constituents, is derived by a rational homogenization procedure based on the Transformation Field Analysis. The micromechanical analysis is developed considering a Cauchy model for masonry components. In particular, a linear elastic constitutive relationship is considered for the blocks, while a nonlinear constitutive law is adopted for the mortar joints, accounting for the damage and friction phenomena occurring during the loading history. Some numerical applications are performed on a Representative Volume Element characterized by a selected commonly used texture, without performing at this stage structural analyses. A comparison between the results obtained adopting the proposed procedure and a nonlinear micromechanical Finite Element Analysis is presented. Moreover, the substantial differences in the nonlinear behavior of the homogenized Cosserat material model with respect to the classical Cauchy one, are illustrated.     

Effective Diffusion Coefficient: From Homogenization to Experiment

In this paper, an example of the application of the homogenization approach (asymptotic expansion technique) to predict the effective diffusion coefficient for an equivalent continuum, together with the experimental verification of the theoretical results is presented. The experimental setup was constructed for the measurements of diffusion in a model periodic porous medium made of Plexiglas. The computer program using the FEM was elaborated to solve the local boundary value problem for a period and to calculate the effective diffusion coefficient. The comparison between the theory and the experiment indicates good agreement between the numerical and experimental values of the effective diffusion coefficient. Interpretation of the test data from the point of view of the homogenization theory is also incorporated.     

三维Cosserat连续体模型与微结构尺寸相关效应的有限元分析

分析了三维Cosserat连续体理论中的应力应变特征,推导了三维Cosserat连续体的有限元方程,基于ABAQUS计算软件提供的用户单元子程序(UEL)接口编写了弹性Cosserat连续体三维20节点有限元程序,并分析了微悬臂梁自由端的挠度问题和微杆扭转问题。通过与基于经典连续体理论的解析解及有限元数值计算结果进行比较,表明所发展的三维Cosserat连续体有限元能有效地模拟微结构尺寸相关效应问题,即随着微结构尺寸与材料内部长度参数的接近,基于Cosserat连续体有限元分析得到的微梁的挠度以及微杆的转角与经典连续体的解析解及有限元解相比越来越小;反之,Cosserat连续体有限元的计算结果与经典连续体的解析解及有限元数值解较为一致。     

AutoCAD图形交换文件DXF在面向对象程序设计中的应用

利用FORTRAN90语言并结合AutoCAD图形交换文件DXF编制出一个面向对象 的有限元程序, 执行该程序, 实现了包括网格自动生成、单元节点信息、等效节点载荷、约 束条件及变形、等应力线图等前后处理的功能. 它克服了以往在有限元教学中面对大量数据 需要输入和整理一大堆输出的结果数据, 做到了有限元前后处理的可视化, 让学生初步掌握 有限元的面向对象编程的技术.