采用ANSYS对混凝土徐变收缩的加载研究.docxVIP

-1-

采用ANSYS对混凝土徐变收缩的加载研究

胡伟伟

河海大学土木工程系，南京（210098）

摘要：概述了采用ANSYS进行混凝土徐变收缩的分析方法，采用自身材料库中所提供金属蠕变本构关系并利用ansys软件自带的参数化设计语言APDL，结合有关规范中对徐变的定义，从而进行对混凝土收缩、徐变的加载实现，最后通过实例进行比较分析，得到了一些有益的结论。

关键词：ANSYS；收缩徐变；APDL；有限元分析

1引言

收缩与徐变是混凝土材料本身固有的时变特性[1]，随着时间的推移而不断变化。其中收缩只与混凝土的构成成分、养护条件和外界环境有关；而徐变则不仅与上述因素有关，还与实际结构所承受的荷载大小和加载龄期有着密切的关系。因客观因素的复杂性，混凝土徐变引起的结构徐变变形及次内力计算的精确分析十分困难，如何在大型软件中引入可考虑徐变的材料本构关系是进行徐变分析的关键。

2徐变理论的基本假定

影响混凝土徐变的因素具有强烈的随机性[2]，而且就许多因素本身来说，也是随机的。因而在实际工程中常采用如下假定：

（1）徐变线性理论。徐变应变ξc与初始弹性应变ξe之比即徐变系数?=ξc/ξe。与持续应力的大小无关，即徐变应变与初始弹性应变成线性关系。

（2）不考虑结果配筋的影响，把结构当作素混凝土。这对预应力混凝土结构（因含筋率较小）是适应的，但对不同材料或相同材料（弹性模量相差较大）组成的复合结构是不适用的。

（3）混凝土的弹性模量假定为常值。试验证明，混凝土弹性模量随时间变化而变化，一般可增加10%—15%。但考虑到徐变系数的计算值中已部分包含了这一因素，可取常值计算。

3ansys中混凝土徐变分析方法

3.1ansys中徐变理论

徐变是指在不变荷载作用下，材料持续变形的特性。如果位移固定，反力或应力将随时间而变小，这种特性有时也称为应力松弛。见图1

图1应力松弛

-2-

徐变的三个阶段如图2所示，在初始徐变阶段，应变率随时间而减小，这个阶段一般发生在一个相当短的时期。在第二期徐变阶段，有一个常应变率，所以应变以常速率发展，在第三期徐变阶段，应变率迅速增加直到材料失效。

图2常应力作用下的徐变

通过一个方程来模拟徐变行为，此方程描述了在实验中观测到的主要特征（特别是在一维的拉伸实验中）。这个方程以徐变应变率的方式表示出来，其形式如下：

上式中，A,B,C,D是从实验中得到的材料常数，常数本身也可能是应力，应变，时间或温度的函数，这种形式的方程被称为状态方程。

以求解算法划分，徐变分析有“隐式”和“显式”2种方法，在此仅介绍显式徐变方法。

显式徐变方法求解徐变应变使用了欧拉朝前法，以时间步开始时的应力、应变为基础计算徐变应变率，在每个时间步长内，徐变应变率被假定是常数，因此有：

使用这种方法所对应的徐变应变曲线如下图3所示，为了减小误差，需要较小的时间步长，特别是在徐变应变率变化很大的区域。

图3显式蠕变对应的蠕应变曲线

-3-

ANSYS程序提供了很多徐变准则，在此仅介绍研究比较多应用比较广泛的C6=0初始徐变方程：

Δξcr=C1σC2ξC3e?C4/TΔt

其中：ξ等效应变（以修正的总应变为基础）

σ——等效应力

T——绝对温度（所有给定温度加上TOFFST）t——在子步结束的时间

C1、C2、C3、C4,为通过TBDATA命令所输入的材料常数

3.2方程中材料常数的确定

由于采用的是线性徐变理论，即徐变的变化率与应力无关，而与应变呈线性关系，所以C2=0，C3=1。又因为不考虑温度对徐变的影响，所以C4=0，最后徐变准则可简化为：Δξcr=C1ξΔt

因在应力不变的条件下，总应变：

ξ(t)=ξcr(t)+ξe=?(t,τ)ξe+ξe

式中，ξcr为徐变应变，ξe为初始弹性应变，t为时间，τ为加载龄期，?为加载龄期为τ时刻计算龄期t时刻的徐变系数。在ansys中

Δξcr=C1ξΔt=C1(ξe+ξcr)Δt=C1[ξe+?(t,τ)ξe]Δt所以

3.3徐变系数的计算

在确定了徐变系数与ansys计算参数C1之间的关系后，就可以用apdl语言编写计算徐变的程序，程序中徐变系数的计算可采用有关规范，从而可以计算出相应时间的徐变系数和creep准则中的C1值，再将该相应时间内的徐变系数和C1赋值给数组，从而实现ansys对徐变的计算。

3.4算例

计算模型