接触面单元翻译flac3d.docVIP

3.2 简述 在FLAC3D中，接触面是由一系列三角形单元的组合而形成的，其中每一个三角形单元都是由3个节点组成。接触面单元可以在空间的任何位置创建。一般来说，接触面单元是和实体单元的表面粘合在一起的；每一个单元的四边形平面由两个三角形接触面单元组成，这样，接触面的节点就自动在每一个接触面单元的顶点创建而成了。但其他网格的表面和一个接触面单元发生接触时，那么，接触面的节点就会发现这种接触，并且由法向，切向刚度及滑动的相关性质所定义。 每一个单元是按照权重的方式分布面积到节点上的。每一个接触面的节点都有一个相关联的特征面积。整个接触面面积就被这些节点分开，如下图3.1所示，阐述了接触面单元，接触面节点和单独节点的特征面积之间的相互关系。 值的注意的是，FLAC3D中的接触面是单面的，这个不同于二维FLAC中接触面为双面的格式规定。我们可以把FLAC3D中的接触面想象成’退缩的包装’，把它拉伸延长覆盖所需表面，这使得和其他面相互渗透作用变的敏感起来。 这种基本的接触关系是由接触面节点和实体单元外表面（目标面）建立起来的，接触力的法向方向也是由目标面的法向方向决定的。 在每一时步中，对每一个接触面节点和接触的目标面来说，他们的绝对法向和相对剪切速度都要计算的，他们的值被接触面的本构模型用来计算法向力和剪切力。其中，本构模型是由线性库仑剪切屈服强度准则定义的，它限制了剪切力只能作用在接触面的节点上，法向刚度，切向刚度，抗拉强度，剪切黏结强度和剪胀角等，在剪切强度的限值到达后，将会在接触的目标作用面上形成一个有效法向应力增量。　默认情况下，在接触面有效应力计算中使用了孔隙水压力，这可以由　interface I effective on/off 命令来生成或解除。下图3.2阐述了作用在接触面的节点上的本构模型的各组成部分: 其中，描述弹性接触面响应的法向和剪切力是由在（t+?t）时步中，使用下列关系计算得到的： 这里，是在（t+?t）时步中的法向力 是在（t+?t）时步中的剪切力 　是接触面节点插入到目标面的法向的渗透距离（量纲：位移） 是相对剪切位移增量 是接触面初始化