三向应力状态下高性能混凝土内时损伤本构模型分析研究.docVIP

三向应力状态下高性能混凝土内时损伤本构模型分析研究 　　【摘 要】本文采用了内时理论与损伤力学相结合的混凝土内时损伤本构模型，该模型的优点在于混凝土的弹塑性特性由内时理论描述，而微裂缝由损伤力学来描述，前者使本构模型摆脱了一般弹塑性模型中的屈服面的概念，从而使模型的参数大大减少，应用到实际工程中更加简便。 　　【关键词】高性能混凝土；三向应力；内时损伤 　　随着混凝土的广泛应用，对混凝土本构关系的研究也取得了许多成果，各种混凝土本构模型不断被提出，但是这些模型彼此之间差异较大、各有优缺点，而且适用于普通混凝土的本构理论对高性能混凝土不一定适用，因此很有必要对高性能混凝土的本构关系展开研究。本文指出了研究三向应力状态下高性能混凝土本构关系的意义，为高性能混凝土的应用提供了一定的理论基础，有助于提高工程技术人员对高性能混凝土内部组份之间作用机理的认识，同时可为进一步研究混凝土在多轴应力作用下的力学性能提供参考。 　　1、内时损伤本构模型的建立 　　1.1引入损伤变量的内时本构方程 　　（1）内蕴时间 　　Valanis[1]应变空间的度量张量采用的是一个和材料性质相关的四阶正定张量 。即使应变状态完全相同，但对于不同材料应变路径上的距离变化仍是不同的。这样Valanis就在应变空间中定义了一条与材料性质及其变化密切有关的“记忆路径”，因而塑性变形下应力响应的路径依赖性质或历史依赖能够通过内蕴时间标度 来加以描述。 　　在理论和试验数据分析的基础上，Bazant [4]对Valanis的内蕴时间关系式进行了简化，Bazant认为 = ≈0 因此取 = =0；由于混凝土材料的塑性体积应变主要受内部微裂缝的影响，当把微裂缝影响看作损伤问题分离出来后，混凝土的塑性体积应变就可忽略不计，且塑性偏应变与塑性体积应变在形成广义内摩擦力时的相互影响也可不考虑，所以可认为 = = =0。 　　（2）内时本构方程 　　耗散型材料本构方程的形式不变性定律[5]：如果能够恰当地定义某些内蕴时间标度Z以使得广义内摩擦力正比于相应的内变量对Z的变化率， 那么所研究的耗散材料的本构方程的形式将与内摩擦力正比于相应内变量速率的粘弹性材料的本构方程形式完全一样。 　　（3）混凝土损伤变量的定义 　　将混凝土的裂缝看作是损伤问题，并将损伤和非损伤从混凝土中分离出来，那么非损伤部分实质上就相当于混凝土的完全“紧密”状态。对于损伤部分则通过采用损伤变量引入损伤本构理论，而非损伤部分则依然服从弹塑性内时本构关系。 　　（4）引入损伤变量的混凝土内时弹塑性关系式 　　采用等效应变假设，认为在结构内任一点，对应于名义应力 与有效应力 的应变相等，都为 。将内蕴时间模型中 及 用损伤定义中的有效应力 、 代替，即得含有效应力 的内时弹塑性关系式。 　　1.2混凝土损伤变量的确定 　　损伤变量 、 的确定采用大连理工宋玉普[6-7]给出的表达式： 　　当 0时， ；式中： 的微分 。 　　1.3增量形式的内时损伤本构模型将式 及 写成增量形对于 0的情况，由 得 对于 ≤0的情况，由 将 的表达式 分别代入上两式，得 当 0时， = 　　当 ≤0的情况， 　　2、内时损伤本构模型的有限元实施 　　本文采用损伤理论和内蕴时间理论相结合的内时损伤理论，基于虚功原理运用计算机编制了非线性有限元程序，下面是计算的具体步骤。 　　矩阵形式的增量方程 　　（1） 　　式中： 为应变矩阵； 为矩阵转置符号， 为微元体积， 为等效节点力增量列阵，用单元节点位移 表示式 中的应变增量，并代入式（1）得到 （2 （3式中 为虚拟塑性力矩阵 （4）式中 的计算分两种情况：当 ≥0时，按式 计算；当 0时按式 计算。 　　由弹塑性分析框图1可知，在第i步的应力状态（ ， ）和第（i+1）步的总应变 已知时，则可以从弹塑性分析中求出第（i+1）步有效应力 ；该步中的塑性应变偏量和球量可以根据 和式 求出；然后根据损伤量的演化方程及结合有效应力表达式，就得到了第（i+1）步的总应力 。通过迭代方法以保证在荷载增量和施加在结构上产生的内部节点力达到平衡，迭代收敛准则采用位移范数准则。 　　上述计算方法与一般计算方法的区别在于不需要判断单元是否进入屈服的附加准则，这正体现了内时理论不以屈服面的概念作为基础，不必遵守流动法则的优越之处。 　　3、本构模型的验证 　　高性能混凝土内时损伤本构模型的计算曲线与高性能混凝土在多轴压应力状态下的实际应力—应变关系曲线相比较，见图2。 　　4、结论 　　由于混凝土内时损伤本构模型的特点在于混凝土的弹塑性特性由内时理论描述，而微裂缝由损伤力学来描述，前者使本构模型摆脱了一般弹塑性模型中的屈服面的概念，从而使模型的参数