混凝土本构关系总结

作业1：总结典型的混凝土本构模型类型，并就每种类型给出有代表性的几个模型按照力学理论基础的不同，已有的本构模型大致分为以下几种类型：以弹性理论为基础的线弹性和非线性弹性本构模型；以经典塑性理论为基础的弹全塑性和弹塑性硬化本构模型；用内时理论描述的混凝土本构模型等。1、混凝土单轴受力应力—应变关系1.1混凝土单向受压应力—应变关系1、saenz等人的表达式saenz等人（1964年）所提出的应力—应变关系为0230000=1(2)(21)()()SEEE?001ESuu1E0图1混凝土单轴受压应力--应变关系2、Hognestad的表达式Hognestad建议的模型，其应力—应变曲线的上升段为二次抛物线，下降段为斜直线，如图2所示，表达式为2000=[2()]0000=[1-0.15()]cu0cu00uu图2Hognestand建议的应力--应变关系3、GB50010—2002建议公式我国《混凝土结构设计规范》所推荐的混凝土轴心受压应力—应变关系为01（上升段）3000[(32)(2)()]aaa01（下降段）00200/(-+c1）式中，a表示应力—应变曲线的上升段参数；c为下降段参数。4、CEB—FIP建议公式CEB—FIP模式规范建议的单轴受压应力—应变关系为20000(/)(/)1(2)(/)kk式中，k为系数，00(1.1)(/)CkE,CE为混凝土纵向弹性模量。2、混凝土非线性弹性本构模型1、混凝土非线性弹性全量型本构模型当材料刚度矩阵[]D用材料弹性模量E和泊松比表达，则为全量E-型；如果材料的刚度矩阵[]D用材料模量K和剪变模量G表达，则为全量K—G型。在全量本构模型中，关键是要合理确定材料参数E和随应力状态变化的规律。Ottosen本构模型的建立过程可分为四个步骤：建立强度和开裂准则；定义非线性指标；建议采用的割线模量SE;建议采用的泊松比s。引入一非线性指数β，表示当前应力(321,,)距破坏(包络面)的远近，以反映塑性变形的发展程度。假定主应力1和2值保持不变，3(压应力)增大至3f时混凝土破坏，则33f混凝土的多轴应力—应变关系仍采用单轴受压的Sargin方程22(1)()1(2)()cccccADfAD但用多轴应力状态的相应值代替033,//icfpffscffsEEAEEEEEE式中各符号的意义见图3。将以上两式代入后，得一元二次方程，解之即得混凝土的多轴割线模量：11222222DEEEEEEEEffiifiis0cuu(a)E0Epffuu(b)EiEfEB00(c)00.81.0it图3Ottosen本构模型（a）单轴受压应力-应变曲线（b）多轴应力-应变曲线（c）泊松比其中，与均以受压为正值；cf为混凝土单轴抗压强度；cEEA/0；0E为混凝土初始弹性模量；cE为混凝土应力达cf时的割线模量；c为应力达峰值时的应变；D为系数，对曲线上升段影响不大，而对下降段影响很大，如图所示。限制0≤D≤1.0，D愈大，则曲线下降愈平缓。这一曲线基本上可以反映混凝土应力应变关系全曲线的主要特征，因而在混凝土有限元分析中应用很广。由上述分析可见，Ottosen模型是以非线性弹性理论为基础，仅仅适当地改变了割线模量SE和泊松比s；所采用的参数仅采用单轴试验数据便可确定；给出的与单轴受压应力—应变全曲线特征相同的一般三轴受压应力—应变曲线，以及峰值应力点和软化段，使计算简单。2.2非线性弹性本构关系—增量型采用全量形式对按比例一次加载的条件是合适的，它与加载路径无关。在逐级加载以及非比例加载情况下采用全量形式会感到困难，这时采用增量形式比较合理。因为采用非线性弹性理论，所以仍假定应力状态与应变状态有一一对应关系，材料参数是应力状态(或应变状态)的函数。但这时不采取全量形式，而采用应力增量与应变增量的形式，材料本构矩阵将应力增量与应变增量联系起来。各向同性增量本构模型（1）含一个可变模量tE的各向同性模量假定泊松比为不随应力状态变化的常数，而用随应力状态变化的变切线模量tE取代弹性常数E，并采用应力和应变增量，则可得下列增量应力—应变关系：1(1)(12)ttijijkkijEEddd式中，ij，ij分别表示应力和应变增量张量；tE为应力—应变曲线上任一点的切线模量，对多轴应力状态下混凝土，可根据等效单轴应力—应变关系确定。（2）含两个可变模量tK和tG的各向同性模量增量应力—应变关系2(32)ijtijttoctijdGdKGd1、双轴正交各向异性增量本构模型混凝土在开裂，尤其是接近破坏时，不再表现出各向同性的性质，而呈现出明显的各向异性的性质。因此，用各向异性描述混凝土开裂后的性能更为合理。（1）正交各向异性本构关系假定混凝土为正交各向异性材料，并且在各级荷载增量内力—应变呈线弹性关系，其应力增量与应变增量关系为321ddd=11211222121212010100(1)EEEdEEEdGd式中，1E、2E为施加一级荷载后在主应力方向的等效切线模量；1d、2d、3d为由荷载增量引起的应力增量；21,为在方向1，2的应力对方向2，l所引起的影响(泊松比)。关于泊松比的值，一般由试验确定，一般在0.15~0.20之间，的值的大小对计算结果影响相对较小，Darwin和Pecknold建议：双向受压时2.0一向受压、一向受拉和双向受拉时41224.06.02.0ccff这一模型被广泛采用。（2）等效单轴应力—应变关系1E和2E值可由双轴受力的应力--应变关系求得，但由于这种应力—应变关系包含了泊松比和微裂缝的影响，因此必须在双轴受力的应力—应变关系中消除泊松比的影响。线弹性材料在单轴受压和双轴受压时的应力—应变曲线如图4所示，可见，由于泊松比的影响，在双轴受压状态下的材料刚度增大了。?i0图4线弹性材料的等效单轴应变双向受压单向受压E00对于混凝土这样的非弹性材料，在双轴受力时刚度增大，除了泊松比的影响外，还有其内部微裂缝开展的影响。消除泊松比影响并考虑微裂缝影响的单轴应力—应变关系，称为等效单轴应力—应变关系，如图5所示，其曲线形状随应力12(/)变化。这样切线模量可由等效单轴应力—应变关系确定。?ic0图5等效单轴应变--应变曲线0Darwin—Pecknold模型在双向应力状态下，Darwin和Pecknold于1977年提出了一个正交异性的应力应变关系矩阵，他们认为在消除了泊松比的影响后，等效的应力应变关系仍可用Saenz公式：20021iciuiciuciuiEEE其中，i为主应力方向（i=1,2）；ic为相当于最大压应力ic时的轴向应变值；iu为等效的单向应变，即;ciiuEiciccE/是最大压应力时的割线模量。任意应力时的切线模量iE:20220[1()][1(2)()()]iiiiiciiuiuisicicEdEEdE混凝土弹塑性本构模型目前采用的弹塑性本构关系可分为两种：全量理论和增量理论。全量理论是塑性小变形理论的简称，适用于简单加载情况，在按比例加载的情况下一般可获得比较满意的结果。增量理论又称流动理论，是描述材料在塑性状态时应力与应变速度或应变增量之间关系的理论。1、混凝土弹塑性增量理论弹塑性增量理论需要对屈服准则、流动法则和硬化法则做出假定。TTffDDdDdffAD上式即为增量理论的弹塑性本构矩阵的一般表达式。式中的A表示硬化参数，其值由材料试验确定。对于“做功硬化”材料，参数A等于在产生塑性变形过程中所做的功。2、混凝土弹塑性全量理论1、全量理论的基本假定（1）假定体积的改变是弹性的，且与平均应力成正比，而塑性变形体积不变不可压缩。（2）假定应变偏量和应力偏量相似且同轴。（3）单一曲线假定：对于同一种材料，无论应力状态如何，其等效应力与等效应变之间有确切的关系。2、应力—应变关系弹塑性应力—应变关系采用如下形式弹性阶段2ijijseG塑性阶段'2ijijseG式中，G为剪变模量；'G为等效应变i的函数。3、全量理论的弹塑性矩阵应力—应变关系的矩阵表达式epD其中epD为弹塑性矩阵，其表达式为121100012100012000300231230232epED对称其中，2=-1-3（12）3、混凝土损伤本构模型在外荷载和环境作用下，由于细观结构缺陷（如微裂缝、微空洞等）引起的材料和结构的劣化过程，称为损伤。从上世纪80年代开始，国内外学者开展了对混凝土材料损伤模型的研究。从混凝土材料本构关系入手，在本构模型中加入损伤变量来反映材料的损伤演变过程。比较典型的损伤模型有：Mazars损伤模型，Loland损伤模型以及分线段模型和分段曲线模型等。下面介绍下Mazars损伤模型：单轴受力状态下混凝土损伤本构模型（1）单轴受拉损伤模型Mazars将脆性材料的应力—应变曲线分为上升段和下降段分别予以描述，如图6所示。?cc图6Mazars损伤模型的相关关系曲线0?cc0u~?c0uD以峰值点（c，c）作为分界点，当≤c时，认为--曲线为线性关系，材料无损伤（D=0）。当＞c时，--曲线按指数规律下降，材料产生的损伤（D＞0）。其受拉应力—应变关系如下：0001expctcctcEEb式中，0E是线弹性阶段的弹性模量；ta，tb是材料常数，其中下标t表示受拉。（2）单轴受压损伤模型设1，2和3为主应变，当沿1方向单向受压时，1=＜0，2=3=，则单向受压时的等效应变为222123=++=-2式中，角括号定义为/2xxx。单向受压时，也是将应力—应变曲线分为上升段和下降段：当c时，认为曲线为线性关系，材料无损伤（D=0）;当＞c时，材料产生损伤（D＞0）。其受压应力—应变关系如下001exp2cccccccEEb