弹塑性本构关系的认识及其在钢筋混凝土中的应用浅谈_塑

1弹塑性本构关系的认识及其在钢筋混凝土结构中的应用浅谈摘要：本文首先对弹塑性本构关系和钢筋混凝土材料的本构模型作了简要概述，然后结合上课所学知识和自己阅读的几篇文章，从材料的屈服准则、流动准则、硬化准则和加载卸载准则等四个方面详细阐述了弹塑性本构关系。最后，结合上述准则简要论述了混凝土这一常用材料在地震作用下的弹塑性本构关系。关键词：弹塑性本构关系，钢筋混凝土，地震UnderstandingofElastoplasticConstitutiveRelationandaBrifeTalkofItsAapplicationtoReinforcedConcreteStructureAbstract：Thispaperfirstlymakesabriefoverviewaboutelastoplasticconstitutiverelationandreinforcedconcreteconstitutivemodel.Then，elaboratingtheelastoplasticconstitutiverelationfromthefouraspectsofmaterialyieldcriterion，flowrule，hardeningrule，loadingandunloadingcriterionbasedonwhatIhavelearnedinclassandreadingfromafewarticles.Lastly，asimplyintroductionontheelastoplasticconstitutiveofreinforcedconcreteunderearthquakeisdemonstrated.Keywords：elastoplasticconstitutiverelation;reinforcedconcretestructure;earthquake1引言钢筋混凝土结构材料的本构关系对钢筋混凝土结构有限元分析结果有重大的影响，如果选用的本构关系不能很好地反映材料的各项力学性能，那么其它计算再精确也无法反映结构的实际受力特征。所谓材料的本构关系，主要是指描述材料力学性质的数学表达式。用什么样的表达式来描述材料受力后的变化规律呢？不同的学者根据材料的性质、受力条件和大小、试验方法以及不同的理论模型等因素综合考虑，建立了许多种钢筋混凝土材料的本构关系表达式。材料的本构关系所基于的理论模型主要有：弹性理论、非线性弹性理论、弹塑性理论、粘弹性理论、粘弹塑性理论、断裂力学理论、损伤力学理论、内时理论等。迄今为止，由于钢筋混凝土材料的复杂因素，还没有一种理论模型被公认为可以完全描述钢筋混凝土材料的2本构关系。有些本构关系虽然能比较好地反映材料的应力应变关系，但是由于试验条件的不同，使得各表达式的结果具有很大的离散性。本文基于钢筋混凝土材料的弹塑性本构模型，确定了地震作用下钢筋混凝土结构的本构模型。这里的弹塑性问题主要是不依赖于时间的弹塑性问题，弹塑性材料进入塑性的特征是当荷载卸去后存在不可恢复的永久变形，因而在涉及卸载的情况下，应力应变之间不再存在唯一的对应关系，这是区别于非线性弹性材料的基本属性。为了更有效地利用通用有限元软件，有必要了解塑性力学的基本法则、弹塑性有限元分析的基本原理[2]。2结构的弹塑性本构模型分析2.1弹塑性力学的基本准则弹塑性理论提供了描述材料弹塑性发展的数学关系。在弹塑性理论中有三个重要的准则：屈服准则、流动准则、硬化准则和加载卸载准则[7]。2.1.1初始屈服准则初始屈服准则决定了材料由弹性变形进入塑性变形的初始应力状态。对于初始各向同性材料，在一般应力状态下开始进入塑性变形的条件是：),(000kFF（2.1）式中：表示应力张量，0k是给定的材料参数。),(000kFF的几何意义可以理解为应力空间的一个超曲面，此曲面称之为初始屈服面。通常采用的屈服条件有：（1）Von-Mises屈服准则Von-Mises屈服准则数学表达式可简化成：0)(),(0000kfkFF（2.2）其中：)(3131:21)(321200mmyIskssf式中：yo为由单向拉伸试验所得的材料屈服强度；s为偏斜应力张量；m为静水应力；I为单位张量。且s和等效应力有以下关系：2231:21Jss（2.3）2J为第二应力不变量，将上式代入（2.2）式，则有0y，所以此Von-Mises屈服准则的力学意义是：当等效应力等于材料的初始屈服应力yo时，材料开始进入塑性变形，几何意义是：在偏斜应力空间内，它代表一个以yo3/2为半径的超球面，即材料用偏斜应力张量表示的应力状态在超球面以内，材料是弹性的；当应力状态到达球面时，材料开始3进入塑性变形。在三维主应力空间，Von-Mises屈服准则表示为：031])()()[(6120213232221y（2.4）式中：1，2和3是三个主应力。该式的几何意义是：在三维主应力空间内，初始屈服面是以321为轴线的圆柱面。此面和过原点并垂直于直线321的平面的交线，即屈服函数0F在平面上的轨迹是以yo为半径的圆周，如图2.1（a）所示。而在03的平面上屈服函数的轨迹是一椭圆，该椭圆的长半轴为yo2，短半轴为yo3/2，如图2.1（b）所示。（2）Tresca屈服准则对受有三个主应力1，2和3作用的固体，其屈服条件应满足下列方程：0])][()][()[(202132023220221yyy（2.5）（a）平面上的屈服轨迹（b）03平面上的屈服轨迹图2.1屈服轨迹此式的力学意义是：当最大剪应力等于初始剪切屈服应力时，材料开始进入塑性变形。则可推出屈服准则：3,2,1,0jiyji（2.6）几何上，（2.5）式表示一个在主应力空间内，以321为轴线并内接于Von-MISes屈服轨迹的正六边形，如图2.1（a）所示。同样，在03的平面内，Tresea屈服轨迹内接于Von-Mises屈服轨迹的六边形，如图2.1（b）。比较以上屈服条件，为计算方便，本文中钢筋的有限元分析采用Von-Mises屈服条件。2.1.2流动准则4流动法则用来规定材料进入塑性应变后的塑性应变增量在各个方向上的分量和应力分量以及应力增量之间的关系。Von-Mises流动法则假设塑性应变增量可从塑性势导出，即Qddp（2.7）其中：pd是塑性应变增量；d是正的待定有限量，它的具体数值与材料硬化法则有关；Q是塑性势函数，一般说它是应力状态和塑性应变的函数。对于稳定的应变硬化材料（随着载荷增大，如果材料的应力增量击和应变增量d所做的功为正功，即0:dddW，此类材料称为稳定材料），Q通常取和后继屈服函数F相同的形式，称之为和屈服函数相关联的塑性势。对于关联塑性情况，流动法则表示为：fdFddp（2.8）从微分学得知，/FF定义的向量正是沿着应力空间后继屈服面0F的法线方向，所以Von-Mises流动法则又称为法向流动法则。2.1.3硬化法则硬化法则是用来规定材料进入塑性变形后的后继屈服面函数（又称加载函数或加载曲面）在应力空间中变化的规则。一般来说，后继屈服函数可以采用以下形式：0),(kF（2.9）其中：k是硬化参数，它依赖于变形的历史，通常是等效塑性应变了的函数。对于理想塑性材料，因无硬化效应，显然后继屈服函数和初始屈服函数相同，即0),(),(00kFkF（2.10）对于硬化材料，与图2.2所示的不同硬化特征相对应，通常采用的硬化法则有：（1）各向同性硬化法则各向同性硬化法则规定，当材料进入塑性变形以后，加载曲面在各方向均匀的向外扩张，但其形状、中心及其在应力空间中的方位均保持不变。例如对于03的情形，初始屈服轨迹和后继屈服轨迹如图2.3（a）所示。如采用Von-Mises屈服条件，则各向同性硬化的后继屈服函数可以表示为：0),(kfkF（2.11）其中：)(31:212pykssf式中的y是现时的弹塑性应力，它是等效塑性应变p的函数，p的表达式为：2/1):32(ppppddd（2.12）5)(py可从材料的单轴拉伸试验的~曲线得到。定义pypddE（2.13）为材料的塑性模量，又称之为硬化系数（见图2.4）。它与弹性模量E及切向模量)/(ddEEtt的关系为ttpEEEEE（2.14）需要指出，各向同性法则主要适用于单调加载情形。如果用于卸载情形，它只适用于反向屈服应力1y数值上等于应力反转点1r的材料。而通常材料是不具有这种性质的，因此在塑性力学中还发展了其他的硬化法则。图2.2各种硬化塑性的特征（a）各向同性硬化（b）Prager运动硬化（c）Zeigler运动硬化6图2.3各种硬化法则示意图图2.4单轴拉伸情况下的弹塑性硬化系数（2）运动硬化法则此法则规定材料在进入塑性以后，加载曲面在应力空间作一刚体移动，但其形状、大小和方位均保持不变。后继屈服函数可表示为0),,(0kF（2.15）其中：0k即初始屈服条件；是加载曲面的中心在应力空间内的移动张量，与材料的硬化特性以及变形的历史有关。根据的具体规定的不同，运动硬化法则可以分为Prager运动硬化法则和Zeigler修正运动硬化法则。2.1.4加载、卸载准则该准则用以判别从一塑性状态出发是继续塑性加载还是弹性卸载，这是计算中判定是否继续塑性变形以及决定采用弹塑性本构关系还是弹性本构关系所必需的。该准则可表述如下：（1）若0F且0:df，则继续进行塑性加载；（2）若0F且0:df，则由塑性变为弹性卸载；（3）若0F且0:df，则应区分下面两种情况，即①对于理想弹塑性材料，此情况是塑性加载，因为在此条件下可以继续塑性流动；②对于硬化材料，此情况是中性变载，即仍保持在塑性状态，但不发生新的塑性流动（0pd）。以上各式中的/ff，按不同材料特性而采用不同的屈服函数形式确定。对于理想塑性材料以及采用各向同性硬化法则的材料，则sf（2.16）对于采用运动硬化法则和混合硬化法则的材料，有7sf（2.17）3钢筋混凝土结构材料弹塑性本构模型3.1基本假设当材料的应力点己处于屈服面上继续弹塑性加载时，需要应用弹一塑性增量的应力应变关系进行弹塑性行为的分析。建立本构模型时应遵循以下假设[9]：（1）采用相关联的流动法则，即在塑性状态下塑性势面与材料的屈服面重合；（2）采用Drucker-Prager模型来描述混凝土的弹塑性本构模型；（3）考虑混凝土和钢筋的屈服强度是塑性应变的函数；（4）钢筋服从Von-Mises屈服准则。3.2结构的弹塑性本构模型首先建立增量形式的虚位移原理[2]。如果tt时刻的应力和体积载荷FF及边界载荷TT满足平衡条件，则在满足几何协调条件的虚位移)(的总虚功等于零，即0)()()(VSttVtdSTTudVFF（3.1）将本构方程表达式代入上式，则可得StVtVtSVVepdSTudVFddSTudVFdD)():(0（3.2）采用与Aitken加速收敛方法相结合的常刚度迭代方案相应的胡克定律增量形式为etetettpetetetetetetettpetDDDDDDDDD):(::)(:):()(:00（3.3）其中：ettetetDDD,,0分别是其材料常数GE,,在tttt,,0时刻数值时的弹性张量；et是t时刻的弹性应变张量。将上式代入有限元离散后的虚位移公式，则可得到用初始弹性刚度矩阵表示的有限元方程，其矩阵形式如下dVDDDD