工程力学专题.pptVIP

河海大学; 一 前言 二 不连续介质力学问题研究进展 三 界面应力元基本理论与方法 四 相关专题与工程应用 五 结语与展望 ;一、前 言;二、不连续介质力学问题研究进展; 结合单元法、节理单元法和薄层单元法都是在连续介质的有限元法框架内，引入特殊的界面单元或简单的力学元件以试图反映不连续性。 结合单元法与节理单元法是在两个连续体 单元之间设置法向与切向弹簧，一旦法向弹簧力为拉力则认为两个连续体之间发生脱离，若切向弹簧力超过其极限强度，则认为两个连续体之间发生相对滑移。 薄层单元法假定块体间界面是由特殊介质组成的厚度很薄的实体单元，这种单元的本构模型取决于周围介质的力学性态由试验直接获得，上述三种方法作为解决不连续变形问题的第一步突破，且便于实施，在工程中得到了较为广泛的应用。 ; 需要指出的是，这类方法是在连续模型的框架内强行引入不连续特性，往往会导致数值方程组一定程度的病态，而且界面上单元数目不能过多，否则界面接触状态会发生混乱而致使数值计算不稳定或发散。 在应用上则表现为难以精细模拟实际岩体中存在的大量走向各异的软弱结构面，一般是采用地质概化的手段加以处理，此外，模型中涉及到的法向和切向刚度系数Kn和Ks的确定也具有很大随意性，为此类模型或多或少蒙上了一些阴影。  ; 用于模拟多节理岩体的等效连续体模型和损伤模型是从另一个角度建立的试图反映材料不连续性的连续介质模型，其关键在于如何建立反映不连续性的参数或变量的演化方程。 损伤模型的研究方法是选择表征损伤(即不连续性)的合适的状态变量——损伤变量，通过实验的途径或连续介质热力学的途径，确定含损伤变量的损伤材料本构方程和损伤演化方程，并对上述方程作出可能的简化，再与连续介质力学的其它场方程一起，形成连续损伤力学的初边值问题的微分提法或变分提法，以求解物体的应力应变场和损伤场。 ; 但问题在于损伤变量是一种含义广泛的内部状态变量，它表征的是材料内部微结构(如位错密度、微裂纹??孔隙等)的不连续性,而难以描述宏观的结构面带来的不连续性，且在具体应用中仍然面临着如何确定有关参数的困难，损伤的扩展规律以及致损的原因也有待深入研究。 第一类模型还有近年来迅速发展的扩展有限元法。 ;模拟裂缝扩展的扩展有限元法;附加非连续函数 ;距离函数 ;含裂纹尖端的单元处理 ; 控制方程 ;区域离散化 ;B矩阵 ;单元刚度矩阵组装;积分策略 ;算例;算例; 为了突破连续介质计算模型的框架，第二类模型近年来获得了快速发展，并已在固体力学、 岩土力学、散体运输、材料、矿冶等多学科中引起广泛关注，其中比较有代表性的方法有离散单元法、刚体—弹簧元模型、无单元法、非连续变形分析模型、数值流形元法、自然单元法以及近年来兴起的PD方法等。 ; 2 离散单元法 由Cundall 提出用以模拟岩石边坡的渐进破坏。该方法以结构面切割而成的离散块体为基本单元，其几何形状取决于岩土结构中不连续面的空间位置及其产状，应用牛顿运动定律描述各块体的运动过程，块体可以发生有限移动与回转(体现了变形的不连续性)，进一步又考虑了岩块自身的变形，Maini与 Cundall于1978年开始研制的离散元通用分析软件UDEC可以模拟岩块破碎与爆炸的运动过程，近年来这一方法又与有限元法或边界元法耦合应用，有效地求解了远场岩体作为连续介质，近场为不连续介质的耦合问题，并逐步应用于各种离散体结构的分析和模拟中。 ; 离散元中一个基本假设是块体运动时动能将转化为热能而耗散掉，因此在计算中，即使是静力问题也必须人为地引入粘性阻尼器以使系统达到平衡，运动趋于稳态，反映在其求解方法上有动态松驰法和静态松驰法。为了保证迭代能进行下去，这二种解法的计算时间步长均需取得足够小，耗时较长，且块体之间阻尼系数的确定带有较大的随意性，另外，离散元对运动学的考虑尚不够充分。上述原因常使数值计算不收敛，难以保证获得各时刻的真实平衡状态。 国内的主要工作者：张楚汉、王咏嘉等 ; 3 刚体—弹簧元模型 最早由日本东京大学Kawai教授在1976年提出，该方法的基本思想在于以离散块体单元的形心位移为基本未知量，用分片的刚体位移模式逼近真实的结构位移场，以单元之间的联接弹簧来反映结构内部的弹性变形，并用弹簧应力表征结构内部的应力，根据计算力学原理建立按位移求解的支配方程。由于该法位移模式为刚体位移模式，故也有刚体有限元法之称。 ;