结构变形缝部位应力场的有限元分析.pdfVIP

结构变形缝部位应力场的有限元分析 1 2 张兴其 ，王大伟 （河海大学土木工程学院南京210098 ） E-mail ：jinlingwind@126.com 摘 要：本文主要采用空间六面体八节点等参单元[1]，编制程序数值模拟结构及其变形缝周 围的应力场；同时针对该问题也使用ANSYS模型在分缝处施加不抗拉的杆单元作类似分析。 关键词：混凝土；变形缝；应力场；有限单元法；数值分析 1. 引 言 在各类结构、大体积混凝土的设计施工中，常常根据规范要求和实际需要设置变形缝， 主要包括伸缩缝、沉降缝、防震缝等，而由于地基变形或不均匀沉降等原因，在变形缝周围 两侧的结构会向中间挤压产生压应力。因此在结构整体分析中,变形缝周围的应力分布状况 往往也是设计者比较关心的问题，如何正确地通过数值分析模拟变形缝周边的应力场也就显 得至关重要。目前各种构筑物数值分析的主要方法是有限单元法，本文将结合计算实例采用 空间六面体八结点等参单元[1]，并在分缝处定义不抗拉单元来计算分缝部位的应力场；同 时还将运用ANSYS 对同一问题进行计算，使两种结果相互印证检验。 2. 基本理论和数值计算的实现方法 2.1 不抗拉单元数值计算的基本原理与实现 对于一些不抗拉材料（如破碎岩体、土体等），在非线性分析中常采用不抗拉单元进行 数值模拟计算。对于本文中所讨论的变形缝就可作为能够抗压，但抗拉强度为零的不抗拉单 元处理。参照文献[2]中对复合地基中类似单元的处理，本文将整体结构剖分为空间六面体 单元，在分缝处耦合空间六面体不抗拉单元，先求出全部单元的高斯点应力，与材料内的初 应力迭加，再计算出主应力 （σ1、σ2 、σ3 ），然后检验不抗拉单元内的高斯点是否出现主 拉应力，当高斯点主应力σ10 （或σ20 或σ30 ）时，可类似于不抗拉单元受拉破坏，而 单元破坏以后它们所不能承担的应力就要向其他单元或高斯点转移，应力重分布。这时需要 对各单元高斯点的应力进行调整，进行应力迁移计算。应力迁移的计算步骤[2]如下： σ 1）首先计算出全部单元的高斯点应力{ }和主应力（σ1、σ2 、σ3 ）； 2 ）对于不抗拉单元，如果出现主拉应力，则对高斯点应力进行修正（假定修正前后主 应力的方向不变），即若σ10，则取σ1=0；若σ20 ，则取σ2=0 ；若σ30，则取σ3=0 ； σ 3 ）然后将修正后的主应力转化为应力分量形式{ }，求出修正前后的应力差 （即需要 迁移的应力值）； （10） ∆σ = σ - σ { } { } { } 4 ）根据变分原理，由上一步中所求出的需迁移的应力值构成的不平衡节点荷载增量； T （11） ∆R B ∆σ dx dy d z { } ∑∫∫∫v [ ] { } e 5 ）将（11）式所求出的不平衡节点荷载增量作用在原来结构上，求出由此产生的应力 增量，并叠加到原来的应力中去，并计算主应力； 6 ）重复步骤 2 ）—5 ），直到所有单元需要迁移的应力值或不平衡节点荷载增量小于容 许应力为止。 - 1 -