静力弹塑性损伤实用本构模型及其稳定性分析.docxVIP

静力弹塑性损伤实用本构模型及其稳定性分析 混凝土结构可承受影响、地震和爆炸等动态载荷。基于混凝土静态力学性能的分析结果与实际情况之间存在较大差异。因此，考虑到适应性率的结构模式在研究中发挥着重要作用。在动态负荷的作用下，混凝土材料表现出适应性，即材料的损伤和非线性程度随着适应性的增加而降低，抗张力和抗压强度随着适应性的增加而增加。混凝土材料的适应性在影响负荷（如碰撞和爆炸）和地震负荷下的作用下是明显的。 研究人员很早就注意到上述现象,提出了多种动力作用下的混凝土本构模型理论,然而,由于这些理论模型本身以及数值实现算法的复杂性,工程实际中应用最为广泛的仍然是一些经验模型,如根据经验提高材料的单轴抗拉强度(或单轴抗压强度)等简化处理方法.考虑应变加载速率的混凝土动力本构模型通常是在率不相关本构模型的基础上做如下扩展得到:① 将率不相关本构模型破坏面扩大,在实际工程中应用较多,但没有物理背景;② 将弹塑性应力—应变关系粘滞化,使得塑性应力—应变关系率相关,在动力加载过程中混凝土微粒塑性应变的发展存在滞后;③ 考虑损伤演化的率相关性,考虑了动力加载情况下混凝土材料损伤的滞后. 文献通过将塑性应变和损伤变量均粘滞化,将其提出的率不相关弹塑性损伤本构模型进行动力推广,建立了相应的率相关本构模型,但模型较为复杂,使用不方便.文献认为:高应变率对微裂纹发展,即损伤演化的迟滞作用导致了混凝土材料非线性程度降低及动力强度提高.文献通过对损伤变量进行粘滞化,提出了适合大型钢筋混凝土结构分析的混凝土率相关弹塑性损伤本构模型,文献采用和文献类似的方法将提出的静力弹塑性损伤本构模型粘滞化,并通过算例表明应变率效应对实际结构的动力反应有一定的影响. 损伤力学能够从本质上描述混凝土材料内部微裂纹的开裂、发展所引起的材料宏观非线性行为,其中损伤变量的发展演化就是为了模拟材料内部微裂纹的开裂发展,从而在损伤力学中可以将损伤变量粘滞化,通过延缓损伤变量的演化来模拟高应变率对微裂纹发展的迟滞作用.笔者在以前的研究中,建议了一个混凝土率不相关弹塑性损伤本构模型,该模型能较好地描述混凝土材料在静力或伪静力加载下的各种非线性行为.基于以上分析,本文将笔者建议的模型进行动力推广,使其能够描述动力加载下混凝土材料的非线性行为. 阻尼是结构振动过程中一种特有的能量耗散机制.目前应用最为广泛的是Rayleigh阻尼,Rayleigh阻尼假定阻尼矩阵C是质量矩阵M和刚度矩阵K的线性组合,数值计算表明:单纯采用质量比例阻尼也无法减小高频噪音的影响.考虑到质量比例阻尼本身并无物理依据,很多学者都仅仅考虑刚度比例阻尼的影响.本文将刚度阻尼应力引入到建立的弹塑性损伤本构模型中,使其能在材料层次考虑阻尼的影响.利用该模型对Koyna混凝土重力坝进行了地震动作用下的非线性时程分析,分析结果与文献符合较好,表明了建立模型的有效性. 本文首先简要介绍笔者建议的静力弹塑性损伤本构模型,然后对其进行扩展使之能够考虑应变率效应和刚度阻尼的影响,最后通过算例从结构位移反应和数值稳定性两个方面探讨应变率效应和刚度阻尼对结构动力反应的影响. 1 本构模型的开发 损伤力学能描述材料的各种非线性本构行为,如刚度退化应变软化等,目前用结合塑性力学和损伤力学的弹塑性损伤本构模型来模拟混凝土材料得到了较为广泛的认同.由于塑性力学本身的复杂性,目前用大多弹塑性损伤本构模型数值处理复杂计算,效率低稳定性不好,且涉及参数较多,一般只给出了特定情况下的材料参数取值,因此较难在实际工程中应用.笔者在以前的研究中提出了一个率不相关弹塑性实用本构模型,通过对弹性Helmholtz自由能进行修正,使得模型能较为准确地模拟混凝土材料在双轴、三轴加载下的非线性本构行为,其弹塑性部分使用经验弹塑性模型降低模型的数值复杂性,下面首先对笔者开发的静力弹塑性损伤本构模型做简要介绍. 1.1 应力球量拟合 损伤能量释放率采用如下公式计算: Y±=-?ψe?d±=12γ|ˉσ±||ˉσ|(ˉσ∶ˉD∶ˉσ-19((ˉσm)2δ∶ˉD∶δ))(1) 式中:标记“∶”为二阶缩并积;Y±表示正、负损伤能量释放率;ψe表示修正后的弹性Helmholtz自由能;d±表示受拉、受压损伤变量;ˉσ表示有效应力张量;ˉσ±为ˉσ的正、负分量;ˉD为有效柔度张量;ˉσm表示有效应力球量;‖X‖=X∶X;γ为本文引入的折减系数;δ表示向量(1,1,1,0,0,0). 考虑到双轴受压加载下(?σ3?σ2?σ1=0,其中:?σ1??σ2??σ3分别为第一、第二、第三主应力)材料强度的提高和加载应力比密切相关,另外不同应力状态下,应力偏量大小相等时(不考虑方向),应力比能在一定程度上反映应力球量的大小,本文将γ取为应力比的表达式: γ=1-0.265(1-Η(?σ2))√|