基于abaqus的钢筋混凝土结构损伤塑性模型.docxVIP

基于abaqus的钢筋混凝土结构损伤塑性模型 1 近似法即局部法 钢筋混凝土结构广泛应用于土木工程中。一般来说，采用线性理论计算和分析结构的力和电压，并通过设计方法确定组件的负荷能力。确定钢筋混凝土梁的有限负荷能力时，根据大量实验数据，确定截面破坏中的负荷分布情况和计算的公式。该方法对于常规设计是适用的,但一般适用于梁、板等简单构件,当结构截面形式复杂或者结构复杂时缺少相应的计算公式,结构设计时的创新需建立在对结构极限承载力准确分析的基础上,需要考虑结构中材料尤其是混凝土的损伤、塑性、断裂对结构进行非线性分析。 随着电子计算机的飞速发展,有限元的出现和广泛应用大大地促进了钢筋混凝土结构的实体分析和非线性分析,在土木工程领域出现了很多著名的通用有限元软件,其中ABAQUS由于其具有的较大非线性计算能力而得到广泛应用。 2 abaqus损伤塑性模型 混凝土材料单轴拉伸和压缩应力—应变曲线(图1)显示:应力—应变曲线存在下降段,应力超过峰值以后卸载时存在刚度退化现象,这与混凝土材料内部存在微小裂缝(损伤)有关。断裂力学和损伤力学是研究带缺陷材料的常用方法,混凝土中的大量微小损伤及分布的裂缝的情况不是很适合简化为一条或者几条裂缝采用断裂力学的方法来研究。损伤力学研究带微裂缝和微空洞的有损材料,研究其在受力时损伤的扩展和演化,比较适合应用于混凝土这种材料。ABAQUS混凝土本构中提供了损伤塑性模型,该模型将损伤指标引入混凝土模型,能较好地模拟混凝土的卸载刚度因损伤增加而降低的特点,可用于钢筋混凝土结构的非线性分析。 ABAQUS提供的混凝土损伤塑性模型是依据Lubliner,Lee和Fenves(1998)提出的损伤塑性模型确定的,其目的是为分析在循环加载和动态加载条件下混凝土结构的力学响应提供普适的材料模型,它考虑了材料拉压性能的差异,主要用于模拟低静水压力下由损伤引起的不可恢复的材料退化。 ABAQUS损伤塑性模型将损伤指标引入混凝土模型,对混凝土的弹性刚度矩阵加以折减,以模拟混凝土的卸载刚度随损伤增加而降低的特点;能够模拟各种结构(梁、桁架、壳和实体)中混凝土和其他准脆性材料;采用各向同性弹性损伤结合各向同性拉伸和压缩塑性理论来表征混凝土的非弹性行为;将非关联硬化引入混凝土弹塑性本构模型中,以期更好地模拟混凝土的受压弹塑性行为;可以人为控制裂缝闭合前后的行为,更好模拟反复荷载下混凝土的反应。 ABAQUS损伤塑性模型中,考虑损伤时的有效应力可以表示为: 其中,d为损伤因子,其值在0(无损)到1(完全失效)之间变化,为有效应力。 有效应力和弹性应力之间的关系可以表示为: 其中,D0el为材料的初始(无损)刚度,εpl为塑性应变。 应力—应变关系为弹性标量损伤关系: 模型中采用损伤指标dc、dt分别来表示混凝土在受压和受拉时损伤引起的刚度退化,以单轴受力为例可以表示为: 由式(4)和式(5)及图1可知,混凝土在受压和受拉损伤后,刚度不是初始的E0,而是分别退化为(1-dc)E0、(1-dt)E0。 3 损伤模型中损伤因子的值方法 3.1 应力—损伤因子计算公式推导 采用ABAQUS损伤塑性模型进行计算时,需要用户分别输入材料的受压、受拉应力—应变曲线,以及材料受压、受拉时损伤因子—非弹性应变曲线。当缺少实验数据尤其是损伤因子—非弹性应变曲线实验数据时,可以结合《混凝土结构设计规范》(GB 50010—2002)中提供的混凝土应力—应变曲线并根据能量等效原理计算得出所需参数。规范中所提供的应力—应变曲线如图2和图3所示。 图中,fc、ft分别为混凝土单轴抗压、抗拉强度;εc为与fc对应的混凝土峰值压应变;εt为与ft对应的混凝土峰值拉应变。 混凝土单轴受压应力—应变曲线方程如下。当x≤1时: 当x1时: 混凝土单轴受拉应力—应变曲线方程如下。当x≤1时: 当x1时: 其中,参数αa、αt、αd的计算见规范中表C.2.1。 根据Sidiroff的能量等价原理,应力作用在受损材料产生的弹性余能与作用在无损材料产生的弹性余能在形式上相同,只要将应力改为等效应力,或将弹性模量改为损伤时的等效弹性模量即可。 无损伤材料弹性余能: 等效有损伤材料弹性余能: 为有效应力: 于是得Ed=E0(1-d)2,则进一步可以得到: 将式(6)至式(9)分别代入可得出以下各式。 单轴受压损伤因子计算公式如下。 当x1时: 单轴受拉损伤因子计算公式如下。 当x1时: 其中,kc=fc*/(εcE0),kt=ft*/(εtE0)。具体计算取值参见《混凝土结构设计规范》(GB 50010—2002)中C.2.2。 由公式(16)及公式(17)可以通过规范中提供的应力—应变曲线计算ABAQUS损伤塑性模型所需的损伤因子参数。 3.2 预应力—损伤因子计算示例 文本提出