弹性层状体系考虑层间接触非线性的有限单元法_张起森资料.pdfVIP

第2 卷第 3 期 土 木 工 程 学 报 1 9 5 9年 8 月 弹性层状体系考虑层间接触非 线性的有限单元分析法 张起森 郑健龙 ( 长沙交通学院 ) 【提耍】 本文根据弹性层状体系层间接触面的实际状态 , 提出了一种考虑层间非线性接 触的有限元增量子结构分析法 。 该方法把层间界面的接触状态当作应力的函数 , 根据剪应力 与抗剪力的对比关系 , 决定各接触点的具体接触性质 : 分离 、 无靡擦接触 ( 完全滑动 ) 、 有 摩擦接触 (完全连续或有摩擦滑动 ) 等 。 计算表明 , 这种分析方法较弹性层状体系理论假定 接触界面完全滑动或完全连续的分析结果 , 是更为合理的 , 一 、 前 言 在弹性层状体系的应力分析中 , 通常假定层间接触界面处于两种极端的情况 , 即完全连 续或完全滑动 。 按不同接触状态求得的应力 , 在某些情况下某些点处 , 有时相差几倍 , 甚致 十几倍 。 工程实际中 , 例如公路的机场道面的路面结构 , 建筑物的地基等 , 它们的层间接触 状态可能既非完全连续 , 亦非完全滑动 。 其层间接触状态是受应力所制约的 。 为了解决这一 问题 , 近年来在弹性层状体系理论分析中 , 又提出了一种所谓 “ 半接触状态 ” 的假定 , 即假 设层间界面是处于完全连续和完全滑动之间的一种状态 , 它有一定的传递剪应力的能力 。 显 然 , 这种假定在某些情况下可以得到比较符合实际的计算结果 。 但是 , 实际上 , 在弹性层状 体系中 , 层间接触界面不可能处于同一的接触状态 。 根据荷载作用应力和界面抗力的对比关 系 , 在同一接触界面上可能出现部分连续 , 部分滑动的混合情况 。 这种接触状态随应力变化 的性质 , 通常称之为接触非线性 。 对于这一类问题 , 用解析法求解一般是 难于实现的 。 为 此 , 我们用有限单元法研究了这类问题 。 该方法的实质是采用了增量 一 子结构法 。 即把荷载分 成若干增量级 , 对每一荷载增量级用子结构法进行求解 : 先求出界面的接触应力 , 然后据此 分别求出上 、 下两接触体的应力 。 方法是 , 把求得的界面接触应力 , 当作已知的节点力作用 在上 、 下两接触体上 , 用普通的有限元方法就可求得各节点的位移 、 应变和应力 。 这就是我 们所需的解答 。 由于界面的接触点相对于接触体的内点来讲 , 数厦是很少的 , 所以求解起来 比较方便 , 甚致使用较小的计算机也可解决较大的问题 。 这就是本方法的优点所在 。 根据计算结果对比看出 , 在弹性层状体系应力分析中 , 考虑层间界面接触非线性是很有 意义的 , 因为它可使层状体系理论的应力分析结果更符合实际 。 本文首先介绍弹性层状体系考虑层间接触非线性分析的增量 一 子结构法 , 然后以多层路面 ( 地基 ) 结构作为例子 , 给出了数值计算结果 , 并把它与弹性层状体系理论考虑层间完全连 续或完全滑动的解进行了对比 。 二 、 接触非线性分析的增盆一子结构法原理 ( 一 ) 在接触非线性分析中主要采用如下荟本假设 : 1 . 物体的组成材料是线性弹性的 , 变形很小 . 2 . 接触表面是连续 、 平顺的 , 3 . 作用在接触面上的摩擦力遵守库仑定律 , 4 . 接触面的力学边界和几何边界条件均用节点力和节点位移表示 。 (二 )续触表面的连续条件 设接触体系由两个接触体口 . 和D , 组成 , 它们之 一 间的接触面是个平面 。 接触 或分离只在 各可能接触点对范围内发生 , 即曾经是分离的点对只可能在本点对的邻域内重新接触 , 反之 亦然 。 取可能接触段非端部节点内角的等分角线为该点的祛方向 , 而端部节点的切方向则取 与该段的端部直线段重合 。 口 . 的可能 接触段各点单位切向量 t 和单位外法向量 n 组成右手 坐标系 , 用作该点的局部坐标 , 并和口 : 的对应点共用 ( 图 1 ) 。 二 粉Pt’ 又叭 图 1 F ig . 图 2 P ig 。 设某一时刻 t 的解已求得 。 记此时第 ‘可能接触点的 向径 、 位移 向量和 节点力分别为 心 、 川 、 川 、 衅及尸 : 、 尸 : 。 下一时段刁 t的增量解 , 根据 t + 刁 t时刻的不同接触状态 , 应满 足下述不同的连接条件 : 1 . 分离状态 尸 ; ‘ + J P I‘= P 二 ‘ + 才P : 。二 尸贾; + J 尸 : ‘ = P : ; + 刁P : 。二 0 ‘ 2 、 接触状态 ( 记拉应力匆正 , 压应力为负 ) 尸二 ‘十 刁尸二 . 二 一 (尸 : ; + 刁尸 : . ) 0 ( 1 ) 无摩擦情形 (摩擦系数f 二 0 ) , 即完全滑动 : 犷二 . + 才 “二 . = 犷