预应力混凝土组合单元的结构多场效应数值的分析.doc

文章编号：0559-9350(2010)00-0000-00 预应力混凝土组合单元的结构多场效应数值分析 陈松，柯敏勇，刘海祥，洪晓林 (南京水利科学研究院，江苏 南京，210024) 摘 要：为开展外界温度、荷载和预应力作用的混凝土结构数值分析，提出了一种针对预应力混凝土结构分析的单元模型。采用单独的空间杆单元模拟预应力钢筋，通过自由度变换把空间杆单元的节点自由度用实体单元节点自由度来表示，在单元网格划分时可不必考虑预应力钢筋方位，且无须固定在实体单元节点上，能直接求解预应力对结构产生的变形和内力。应用本文提出的单元信息生成技术，可方便地生成组合单元中的信息，编制了预应力钢筋混凝土结构分析的温度场和徐变应力场三维有限元程序。通过算例的对比分析表明，该程序能够有效模拟预应力混凝土结构分析在温度、荷载和预应力作用条件下的结构响应。 关键词：组合单元；温度场；徐变应力场；预应力混凝土 中图分类号：TU375;O242.21 文献标识码：A 1 前言 有关三维有限元分析结构温度场或者应力场的文献比较常见[1-4]，而对于工程中像渡槽、桥梁等结构工作环境复杂，其混凝土结构既有外界温度、荷载和预应力作用同时又受内部钢筋约束，用三维有限元进行结构分析就相当困难了。对这些结构的分析通常忽略了温度的影响，且研究多是利用板、壳单元或杆单元进行一维或平面分析[5,6]，或者在温度场和温度徐变应力的计算中不考虑钢筋和预应力的影响[7,8]。 在预应力混凝土的计算中，对预应力钢筋的处理也很复杂，相关学者在单元模型上做了一些研究，但有限元网格划分和单元信息生成仍然是一个难题[5]，有的单元模型存在缺陷[9]，可能会出现某些区域内为空，即不存在实体结构或认为弹模为零。若要精确计算，则应采用空间实体单元模型计算出混凝土中的应力分布，并采用空间杆单元模拟预应力钢筋，进行预应力效应分析。本文提出了一种新的预应力组合单元模型，该模型最大优点是：①划分混凝土单元网格时可不受预应力钢筋方位的限制，大大地增强了结构在单元网格划分的自由性；②采用了八节点六面体单元，能够描述复杂结构的几何特征（如变截面结构和有中空的结构）；③输入信息量较少，能广泛应用于各种复杂结构形式。在组合单元的基础上编制了预应力钢筋混凝土结构温度场、徐变应力场的三维有限元程序，为预应力结构的多场效应分析提供了数值模拟方法。 2 有限元控制方程 2.1 基本假定 本文着重研究预应力及外部因素对结构下挠的影响。为了减少单元量和简化计算，在结构分析前对混凝土结构与材料特性作出如下假定：⑴普通钢筋、预应力筋与混凝土充分粘结，相互间不发生相对位移；⑵只计温度的改变对结构变形（应力）影响，略去结构变形对温度场改变的作用。 2.2 预应力混凝土的单元模型及其刚度矩阵 预应力混凝土实体单元中包含混凝土、普通钢筋和预应力筋，采用组合式模型建立单元。其中，混凝土采用八节点六面体实体单元，为了简化模型普通钢筋采用“整体式”模型，预应力钢筋则采用杆单元。本文提出的预应力混凝土实体单元模型是在相关研究[10,11] 基础上通过自由度变换，把空间杆单元的节点自由度用实体单元的节点自由度来表示，再利用 收稿日期：2009-7-2 基金项目：国家863计划项目()；西部交通科技项目() 作者简介：陈 松（1980-），男，四川达州人，工学博士研究生，（E-mail）chensong621002@163.com 虚功原理推导得到等效刚度阵和等效节点荷载的一种有限元计算方法，单元刚度分别考虑混凝土、普通钢筋和预应力钢筋的贡献，即 （1） 式中：下标“e”、“c”、“s”、“p”分别表示组合实体单元刚度阵、混凝土单元刚度阵、普通钢筋单元刚度阵及预应力钢筋单元刚度阵。 2.2.1 混凝土及普通钢筋对组合单元刚度的贡献 混凝土及普通钢筋的刚度矩阵表达式为： （2） （3） 式中： 式中：为钢筋的弹模；，和分别为沿X，Y和Z方向的配筋率。 2.2.2 预应力钢筋对组合单元刚度的贡献 当组合单元划分的足够小时，其内部预应力钢筋可由直线代替，可将预应力钢筋作为杆单元处理。 预应力钢筋的杆单元映射到局部坐标系成为一曲线，杆单元的整体坐标系S映射成为局部坐标系r。设杆单元在整体坐标系中的节点坐标为、、，引入杆单元插值函数，可得到杆单元上任意一点的坐标为：