1、钢筋混凝土破坏准则及本构关系钢筋混凝土破坏准则及本构关系 在主应力空间中，在主应力空间中，与各坐标轴保持等距的各点连结成为静水与各坐标轴保持等距的各点连结成为静水压力轴（即各点应力状态均满足：压力轴（即各点应力状态均满足：1=2=3）。）。此轴必通过坐标原点，且与各坐标轴的夹角相等，此轴必通过坐标原点，且与各坐标轴的夹角相等，均为均为)3/1cos(arc 静水压力轴上一点与坐静水压力轴上一点与坐标原点的距离称为标原点的距离称为静水压静水压力（力（）；）；其值为其值为3个主应力在静水个主应力在静水压力轴上的投影之和，故：压力轴上的投影之和，故：cot132133313/)(mI-1-3-231
2、2+(1,2)-(1,2)静水压力轴静水压力轴垂直于静水压力轴的平面为偏平面。垂直于静水压力轴的平面为偏平面。3个主应力轴在偏平面上的投影各成个主应力轴在偏平面上的投影各成120o角。同一偏平面上的每一点的角。同一偏平面上的每一点的3个主个主应力之和为一常数：应力之和为一常数：I1为应力张量为应力张量ij的第一不变量的第一不变量1321Iconst 偏平面与破坏包络曲面的交线成为偏平面包络线。不同静水压力下的偏平面偏平面与破坏包络曲面的交线成为偏平面包络线。不同静水压力下的偏平面包络线构成一族封闭曲线。包络线构成一族封闭曲线。偏平面包络线为三折对称，有夹角偏平面包络线为三折对称，有夹角60o范
3、围内的曲线段，和直范围内的曲线段，和直线段一起共同构成全包络线。取主应力轴正方向处为线段一起共同构成全包络线。取主应力轴正方向处为=0o，负，负方向处为方向处为=60o，其余各处为，其余各处为0o60o。在偏平面上，包络线上一点至静水压力轴的距离称为偏应力在偏平面上，包络线上一点至静水压力轴的距离称为偏应力 r。偏应力在偏应力在=0o处最小处最小(rt），随），随角逐渐增大，至角逐渐增大，至=60o处为最大处为最大(rc），故），故rt rc。一些特殊应力状态的混凝土强度点，在破坏包络面上占有特定的一些特殊应力状态的混凝土强度点，在破坏包络面上占有特定的位置。从工程观点，混凝土沿各个方向的力学
4、性能可看作相同，即位置。从工程观点，混凝土沿各个方向的力学性能可看作相同，即立方体试件的多轴强度只取决于应力比例立方体试件的多轴强度只取决于应力比例 1：2：3，而与各应力，而与各应力的作用方向的作用方向X、Y、Z无关。例如：无关。例如：混凝土的单轴抗压强度混凝土的单轴抗压强度 fc 和抗拉强度和抗拉强度 ft 不论作用在哪一个方向，不论作用在哪一个方向，都有相等的强度值都有相等的强度值。在包络面各有在包络面各有3个点，分别位于个点，分别位于3个坐标轴的负、个坐标轴的负、正方向；正方向；同理，混凝土的二轴等压（同理，混凝土的二轴等压（1=0，f2=f3=fcc）和等拉（）和等拉（3=0，f1=
5、f2=ftt）强度）强度位于坐标平面内的两个坐标轴的等分线上，位于坐标平面内的两个坐标轴的等分线上，3个坐标面内各有一点；个坐标面内各有一点；混凝土的三轴等拉强度（混凝土的三轴等拉强度（fl=f2=f3=fttt)只有一点且落在静水压力轴的正方向。只有一点且落在静水压力轴的正方向。对于任意应力比对于任意应力比(flf2f3)的三轴受压、受拉或拉压应力状态，从工程观点的三轴受压、受拉或拉压应力状态，从工程观点考考虑混凝土的各向同性，虑混凝土的各向同性，可由坐标或主应力可由坐标或主应力(fl，f2，f3)值的轮换（破坏横截面三重值的轮换（破坏横截面三重对称），在应力空间中各画出对称），在应力空间中
6、各画出6个点，位于同一偏平面上，且夹角个点，位于同一偏平面上，且夹角值相等。值相等。破坏包络曲面的三维立体图既不便绘制，又不适于理解和应用，常改用拉破坏包络曲面的三维立体图既不便绘制，又不适于理解和应用，常改用拉压子午面和偏平面上的平面图形来表示。压子午面和偏平面上的平面图形来表示。拉压子午面拉压子午面为为静水压力轴与任一主应力轴（如图中的静水压力轴与任一主应力轴（如图中的3轴）组成的平面，轴）组成的平面，同时通过另两个主应力轴（同时通过另两个主应力轴（1，2）的等分线）的等分线。此平面与破坏包络面的交。此平面与破坏包络面的交线，分别称为拉、压子午线。线，分别称为拉、压子午线。1、拉子午线的应
7、力条件为、拉子午线的应力条件为1 2=3，线上特征强度点有单轴受拉，线上特征强度点有单轴受拉(ft,0,0)和二轴等压和二轴等压(0,-fcc,-fcc）在偏平）在偏平面上的夹角为面上的夹角为=0o;2、压子午线的应力条件则为、压子午线的应力条件则为1=2 3，线上有单轴受压，线上有单轴受压(0,0,-fc)和二和二轴等拉轴等拉(ftt,ftt,0)，在偏平面上的夹角，在偏平面上的夹角=60o。3、拉、压子午线与静水压力轴同交、拉、压子午线与静水压力轴同交于一点，即三轴等拉于一点，即三轴等拉(fttt,fttt,fttt)。拉、。拉、压子午线至静水压力轴的垂直距离压子午线至静水压力轴的垂直距离
8、即为偏应力即为偏应力 rt 和和 rc。=0o=60o 拉压子午线的命名，并非指应力状态的拉或压，而是相应于拉压子午线的命名，并非指应力状态的拉或压，而是相应于三轴试验过程。三轴试验过程。若试件先施加静水应力若试件先施加静水应力1=2=3，后在一轴，后在一轴1上施加拉力，上施加拉力，得得1 2=3，称拉子午线；，称拉子午线；若试件先施加静水应力若试件先施加静水应力1=2=3，后在另一轴，后在另一轴3上施加压力，上施加压力，得得1=2 3，称压子午线。，称压子午线。另外也可以理解为以单轴拉、另外也可以理解为以单轴拉、压条件定义拉、压子午线，即单压条件定义拉、压子午线，即单轴拉状态所在的子午线成为
9、拉子轴拉状态所在的子午线成为拉子午线，而单轴压状态所在的子午午线，而单轴压状态所在的子午线成为压子午线。线成为压子午线。试验研究指出，混凝土的三维试验研究指出，混凝土的三维破坏面也可用三维主应力空间破破坏面也可用三维主应力空间破坏曲面的圆柱坐标坏曲面的圆柱坐标，r，来描述，来描述，其本身也是应力不变量其本身也是应力不变量。=0o=60o12oNr31=2=3oct3oct3圆柱坐标系及主应圆柱坐标系及主应力空间应力分解力空间应力分解，r，的几何表示的几何表示12oNP(1,2,3)r3e=60o=0orcrt拉子午线拉子午线压子午线压子午线偏平面-3+3-(1,2)等应力轴和一个主应力轴组成的
10、平等应力轴和一个主应力轴组成的平面通过另两个主应力轴的等分线面通过另两个主应力轴的等分线转换过转换过程归纳程归纳偏平面偏平面1-12-2-33rN静水应力偏斜应力偏斜应力平面中矢量的方向P与时间有关的混凝土受力性能，如定应力或变应力作用下的徐变（松弛）、收缩、；E0 混凝土单轴受压的初始切线弹性模量；另一方面，若按混凝土规范三轴抗压强度设计值进行验算，相同应力比例下的三轴抗压强度仅为：也可满足承载能力要求，混凝土的单轴受压和受拉应力-应变关系；此二模型的主要计算原则简要介绍如下，详细的计算公式、推导过程和参数值等请见有关文献。将这些钢筋混凝土的特殊本构关系引入结构的非线性分析，完全有理由称之为
11、钢筋混凝土力学。结构是连续体，在荷载作用下会发生变形和位移，但仍应为连续体。由于混凝土大量地采用地方性材料，施工制作工艺和质量控制水平出入较大，使混凝土的实际力学性能有较大的变异性和离散度。28 fc，S0=-0.由上式可知，将上图的坐标缩小当按照荷载步每一步进行计算时，由当前应力（变）状态确定E1、E2和值，代入式（15）或式（16）的本构方程计算应力（或应变）的增量，和前此的增量和累加作为下一步计算的依据，直至完成全部荷载步的分析，获得最终的应力（变）状态和受力的非线性全过程。计算破坏准则的各项指标和参数值：P(1,2,3)二轴等压(fcc=1.当按照荷载步长逐步地进行计算时，由当前的混凝
12、土主应力（1、2、3）值确定非线性指数，再由式（10）和式（12）计算Es和s代入各向同性的线弹性本构关系式(1)或式(3)，即可行结构的有限元分析运算。将以上图形绕坐标原点逆时针方向旋转一角度将以上图形绕坐标原点逆时针方向旋转一角度(90o)，得到以，得到以静水压力轴静水压力轴()为横坐标、偏应力为横坐标、偏应力(r)为纵坐标的拉、压子午线。为纵坐标的拉、压子午线。于是，空间的破坏包络面于是，空间的破坏包络面改为由子午面和偏平面上的包络曲线改为由子午面和偏平面上的包络曲线来表达。破坏面上任一点的直角坐标来表达。破坏面上任一点的直角坐标(fl，f2，f3)改为由圆柱坐改为由圆柱坐标标(,r,)
13、来表示来表示，换算关系为：，换算关系为：)6/()2(cos33/)()()(33/)(321213232221321rfffffffffrfffoctoct 由上式可知，将上图的坐标缩由上式可知，将上图的坐标缩小小 可以用八面体正应力（可以用八面体正应力（oct）和剪应力（和剪应力（oct）坐标代替静水）坐标代替静水压力和偏应力坐标，得到相应的压力和偏应力坐标，得到相应的拉、压子午线和破坏包络线。拉、压子午线和破坏包络线。3 根据试验结果绘制的拉、压子午线和偏平面包络线。根据试验结果绘制的拉、压子午线和偏平面包络线。子午线按照偏平面夹角划分，试验点的子午线按照偏平面夹角划分，试验点的=306
14、0o 分别列在横分别列在横坐标轴的上、下。坐标轴的上、下。试验时测试试验时测试=0o60o的扇形的扇形（其他的扇形是对称的）（其他的扇形是对称的）偏平面包络线则以八面体应力值分段给出。图中曲线为混凝偏平面包络线则以八面体应力值分段给出。图中曲线为混凝土破坏准则的理论值。土破坏准则的理论值。根据国内外混凝土多轴强根据国内外混凝土多轴强度的大量试验资料分析，破度的大量试验资料分析，破坏包络曲面的几何形状具有坏包络曲面的几何形状具有如下特征：如下特征：曲面连续、光滑、外凸；曲面连续、光滑、外凸；对静水压力轴三折对称，对静水压力轴三折对称，当应力状态为静水应力与单当应力状态为静水应力与单向拉应力叠加时
15、，向拉应力叠加时，=0o，故，故=0o的子午线称为受拉子午的子午线称为受拉子午线。如将单向拉应力换为压线。如将单向拉应力换为压应力，则相应于受压子午线，应力，则相应于受压子午线，=60o。破坏曲线与等应力轴破坏曲线与等应力轴有关。在有关。在轴的正向，静水压力轴的拉端轴的正向，静水压力轴的拉端封闭，顶点为三轴等拉应力状态；在封闭，顶点为三轴等拉应力状态；在轴的负向，压端开口，不与轴的负向，压端开口，不与静水压力轴相交，破坏曲线的开口随静水压力轴相交，破坏曲线的开口随轴绝对值的增大而增大；轴绝对值的增大而增大；子午线上各点的偏应力或子午线上各点的偏应力或八面体剪应力值，八面体剪应力值，随静水压随静
16、水压力或八面体正应力的力或八面体正应力的代数值代数值的减小而单调增大，但斜率的减小而单调增大，但斜率渐减，有极限值；渐减，有极限值；偏平面上的封闭曲线三折偏平面上的封闭曲线三折对称，其形状对称，其形状随静水压力或随静水压力或八面体正应力值的减小，由八面体正应力值的减小，由近似三角形近似三角形(rtrc0.5)逐渐逐渐外凸饱满，过渡为一圆外凸饱满，过渡为一圆(rtrc=1）。）。4.7.2破坏准则破坏准则 将混凝土的破坏包络曲面用数学函数加以描述，作将混凝土的破坏包络曲面用数学函数加以描述，作为判定混凝土是否达到破坏状态或极限强度的条件，为判定混凝土是否达到破坏状态或极限强度的条件，称为破坏准则
17、或强度准则。称为破坏准则或强度准则。虽然它不属基于机理分析、虽然它不属基于机理分析、具有明确物理概念的强度理论，但它是大量试验结果具有明确物理概念的强度理论，但它是大量试验结果的总结，具有足够的计算准确性，对实际工程有重要的总结，具有足够的计算准确性，对实际工程有重要的指导意义。的指导意义。1、分类：、分类：借用古典强度理论的观点和计算式借用古典强度理论的观点和计算式;以混凝土多轴强度试验资料为基础的经验回归式；以混凝土多轴强度试验资料为基础的经验回归式；以包络曲面的几何形状特征为依据的纯数学推导式，以包络曲面的几何形状特征为依据的纯数学推导式，参数值由若干特征强度值标定。各个准则的表达方式参
18、数值由若干特征强度值标定。各个准则的表达方式和简繁程度各异，适用范围和计算精度差别大，使用和简繁程度各异，适用范围和计算精度差别大，使用时应认真选择。时应认真选择。因此，建立或选择本构关系是结构非线性分析的关键问题，成为近20年混凝土结构的一个重要研究方向。采用上述应力量致使准则的数学形式差别很大，不便作深入对比分析。即使如此，仍然因为混凝土材性的离散、变形成分的多样和影响因素的众多等而在一定范围内变动。一些常用的、有代表性的混凝土破坏准则列于下表,同时给出了原始表达式和统一表达式，可看到两者中参数的互换关系。以包络曲面的几何形状特征为依据的纯数学推导式，参数值由若干特征强度值标定。至今国内外
19、在混凝土结构的非线性有限元分析中常用的非线弹性本构模型有 Ottosen 模型和 Darwin Pecknold 模型等，有些国际规范也明确建议采用这两种模型。式（5）中的刚度矩阵对称，独立的弹性常数只有6个，加上式（6）中的3个常数，故正交异性材料的独立弹性常数共为9个。确定了合适的本构关系、进行非线性的全过程分析，有可能改变目前的钢筋混凝土结构的内力弹性分析和截面承载力经验性计算等不尽理想的景况，走向更完善、准确的理论解方向。的破坏形态和应力-应变曲线的分类；按此公式可计算各种应力状态下的混凝土多轴强度理论值，并绘制子午线和偏平面包络线，以及二轴和三轴包络线。其中包含了3个弹性常数：当=0
20、o增加到60o时，ct逐渐增加至cc，符合光滑、外凸的特性；混凝土破坏准则中包含的5个参数，可以用全部试验数据进行回归分析拟定，也可在破坏包络面上，或拉、压子午线上选定任意5个特征强度值加以标定。混凝土当然不是线弹性材料，上述线弹性本构关系用于分析混凝土结构时，其适用范围和计算精度显然都受限制，因而建立和发展了非线（性）弹性类本构关系。较好的准则：过王、Ottosen和Podgorski准则；1fc，解得的4个参数为：a=1.可利用线弹性本构关系的已有分析和计算程序；此强度值大于按下图所给的混凝土多轴抗压强度设计值。混凝土的三轴抗压强度为：2、著名的古典强度理论包括：、著名的古典强度理论包括：
21、最大主拉应力理论（最大主拉应力理论（Rankine)；最大主拉应变理论（最大主拉应变理论（Mariotto）；）；最大剪应力理论最大剪应力理论(Tresca)；统计平均剪应力理论（统计平均剪应力理论（Von Mises)；Mohr-Coulomb理论；理论；Drucker-Prager理论。理论。共同特点：共同特点：针对某种特定材料而提出，对于解释材料破坏的内在原因和针对某种特定材料而提出，对于解释材料破坏的内在原因和规律有明确的理论（物理）观点，有相应的试验验证，破坏包规律有明确的理论（物理）观点，有相应的试验验证，破坏包络面的几何形状简单，计算式简明，只含络面的几何形状简单，计算式简明，只
22、含1个或个或2个参数，其值个参数，其值易于标定。因而，它们应用于相适应的材料时，可在工程实践易于标定。因而，它们应用于相适应的材料时，可在工程实践中取得良好的效果。例如中取得良好的效果。例如.Von Mises准则适用于塑性材料（如软准则适用于塑性材料（如软钢），在金属的塑性力学中应用最广；钢），在金属的塑性力学中应用最广；Mohr-Coulomb准则反准则反映了材料抗拉和抗压强度不等（映了材料抗拉和抗压强度不等（ftfc）的特点，适用于脆性的）的特点，适用于脆性的土壤、岩石类材料，在岩土力学中广为应用。土壤、岩石类材料，在岩土力学中广为应用。3、以混凝土多轴强度试验资料为基础的经验回归式、以
23、混凝土多轴强度试验资料为基础的经验回归式 随试验数据的积累，许多研究人员提出了若干基于试验结果、随试验数据的积累，许多研究人员提出了若干基于试验结果、较为准确、但数学形式复杂的混凝土破坏准则。准则中一般需较为准确、但数学形式复杂的混凝土破坏准则。准则中一般需要包含要包含45个参数。个参数。这些破坏准则的原始表达式中采用了不同的应力量作这些破坏准则的原始表达式中采用了不同的应力量作为变量，分为变量，分5种：种：主应力主应力fl,f2,f3;应力不变量应力不变量Il,J2,J3;静水压力和偏应力静水压力和偏应力,r,;八面体应力八面体应力 oct,oct;平均应力平均应力m,m。采用上述应力量致使
24、准则的数学形式差别很大，不采用上述应力量致使准则的数学形式差别很大，不便作深入对比分析。但这些应力量借助下列基本公式便作深入对比分析。但这些应力量借助下列基本公式可以很方便地互相变换：可以很方便地互相变换：采用上述应力量致使准则的数学形式差别很大，不便作深人采用上述应力量致使准则的数学形式差别很大，不便作深人对比分析。但这些应力量借助下列基本公式可以很方便地互相对比分析。但这些应力量借助下列基本公式可以很方便地互相变换：变换：3oct31.52332132123213212213232221013210J22JJ33cos3 30262322232cos353323)()()(333或moct
25、moctcmoctcfffrfffJffffffrJfffffffIffff 最终可统一用相对八面体强度（最终可统一用相对八面体强度（0=oct/fc和和0=oct/fc）表达，）表达，经归纳得子午线方程的经归纳得子午线方程的3种基本形式：种基本形式：最终可统一用相对八面体强度（最终可统一用相对八面体强度（0=oct/fc和和0=oct/fc）表达，经归纳得子午线方程的）表达，经归纳得子午线方程的3种基本形式：种基本形式：HGFEDCBA)(0020002000 一些常用的、有代表性的混凝土破坏准则列于下表一些常用的、有代表性的混凝土破坏准则列于下表,同时给出了原始表达式和统一表达式，可看到两
26、者中同时给出了原始表达式和统一表达式，可看到两者中参数的互换关系。参数的互换关系。过镇海、王传志、张秀琴等搜集了国内外大量的混过镇海、王传志、张秀琴等搜集了国内外大量的混凝士多轴强度试验数据，与按上述准则计算的理论值凝士多轴强度试验数据，与按上述准则计算的理论值进行全面比较，根据三项标准：进行全面比较，根据三项标准：计算值与试验强度的相符程度；计算值与试验强度的相符程度；适用的应力范围宽窄；适用的应力范围宽窄；理论破坏包络面几何特征的合理性等加以评定。理论破坏包络面几何特征的合理性等加以评定。所得结论为：所得结论为：较好的准则：过较好的准则：过王、王、Ottosen和和Podgorski准则；
27、准则；一般的准则：一般的准则：Hsieh-Ting-Chen，Kotsovos，Willam-Warnke准则；准则；较差准则：较差准则：Bresler-Pister准则。准则。在结构的有限元分析中，可根据结构的应力范围和在结构的有限元分析中，可根据结构的应力范围和准确度要求选用合理的混凝土破坏准则。准确度要求选用合理的混凝土破坏准则。混凝土的二轴强度fci按确定的破坏准则进行计算，二轴峰值应变ci和等效泊松比，区别二轴压-压、拉-压和拉-拉应力状态，各按经验公式进行计算。将这些钢筋混凝土的特殊本构关系引入结构的非线性分析，完全有理由称之为钢筋混凝土力学。不同静水压力下的偏平面包络线构成一族封
28、闭曲线。确定这5个参数采用的混凝土特征强度值为：式中 Kii,ii正应力ii和正应变ii间的刚度系数，即弹性模量；横向变形系数、即泊松比；根据国内外混凝土多轴强度的大量试验资料分析，破坏包络曲面的几何形状具有如下特征：因此，在设计和分析结构时应选择合理和适用的本构模型。试验研究指出，混凝土的三维破坏面也可用三维主应力空间破坏曲面的圆柱坐标，r，来描述，其本身也是应力不变量。(非弹性模型)八面体应力 oct,oct;上两个本构模型都是以混凝土单轴受压应力-应变曲线方程为基础建立的，所得多轴应力-应变理论曲线与单轴和二轴受压类的试验曲线相符，但对三轴受压、受拉应力-应变曲线有不同程度的误差。这是最
29、简单、最基本的材料本构关系。较好的准则：过王、Ottosen和Podgorski准则；例如，本构关系有弹性的、塑性的，还有与时间相关的黏弹性、黏塑性的，与温度相关的热弹性、热塑性等。式中 E1、E2、E3 3个相互垂直方向的弹性模量；1fc，解得的4个参数为：a=1.任何类型、体系和受力状况的结构或其局部都可依靠非线性分析方法求解。分别代入上式，用迭代法计算的参数值：为了保证本构关系的可靠性，上述方法按优选次序排列。因而，它们应用于相适应的材料时，可在工程实践中取得良好的效果。4、以包络曲面的几何形状特征为依据的纯数学推导公式以包络曲面的几何形状特征为依据的纯数学推导公式 模式规范模式规范CE
30、B FIP MC90C采纳了采纳了Ottosen准则。它根据偏平面准则。它根据偏平面包络线由三角形过渡为圆形的特点、应用薄膜比拟法：即在等边包络线由三角形过渡为圆形的特点、应用薄膜比拟法：即在等边三角形边框上蒙上一薄膜，承受均匀压力后薄膜鼓起，等高线的三角形边框上蒙上一薄膜，承受均匀压力后薄膜鼓起，等高线的形状由外向内的变化恰好相同据此建立了二阶偏微分方程，求形状由外向内的变化恰好相同据此建立了二阶偏微分方程，求解后转换得到以应力不变量表达的破坏准则式：解后转换得到以应力不变量表达的破坏准则式：011222cccfIbfJfJa)3cos(cos313cos 0cos3 ,30 )3cos(c
31、os31cosr1 0cos3 ,30 211o211okkkk时即当时即当011222cccfIbfJfJa)3cos(cos313cos 0cos3 ,30 )3cos(cos31cosr1 0cos3 ,30 211o211okkkk时即当时即当其中：其中：a和和b决定子午线的形状，决定子午线的形状，k1和和k2分别决定偏平面包分别决定偏平面包络线的大小和形状。络线的大小和形状。标定参数值的标定参数值的4个特征强度值取为：个特征强度值取为：单轴抗压单轴抗压(-fc)、单轴抗拉、单轴抗拉(ft）、二轴等压）、二轴等压(fcc=1.16 fc)三轴抗压强度三轴抗压强度22/J ,5/,602
32、10ccffI三轴抗压强度三轴抗压强度按下式计算各特征强度的按下式计算各特征强度的22/J ,5/,60210ccffI代入代入值 ,J,21I3oct31.52332132123213212213232221013210J22JJ33cos3 30262322232cos353323)()()(333或moctmoctcmoctcfffrfffJffffffrJfffffffIffff011222cccfIbfJfJa)3cos(cos313cos 0cos3 ,30 )3cos(cos31cosr1 0cos3 ,30 211o211okkkk时即当时即当 得得4阶联立方程，解得各参数值。
33、若取阶联立方程，解得各参数值。若取ft=0.1fc，解得，解得的的4个参数为：个参数为：a=1.2759,b=3.1962 k111.7365，k2=0.9801Hsieh-Ting-Chen和和Podgorski准则是对准则是对Ottosen准则的简准则的简化和修正。化和修正。我国的我国的混凝土结构设计规范混凝土结构设计规范附录附录C.4中采纳了过中采纳了过王准则，王准则，其与试验结果相符较好、以八面体应力无量纲量表达、应用幕其与试验结果相符较好、以八面体应力无量纲量表达、应用幕函数拟合混凝土的破坏包络面，一般计算式为函数拟合混凝土的破坏包络面，一般计算式为:octoctoctctdcoct
34、coctdcoctfffarcrfffffffffcccfcfbacbaf232cos 31)()()(31313)23(sin)23(cos/32121323222132125.10004.7.3、规范中的破坏准则、规范中的破坏准则 破坏准则的计算公式破坏准则的计算公式式中式中5个参数都有明确的几何（物理）意义：个参数都有明确的几何（物理）意义：当当 a=0,max时，时，0时时0有极限值（高压应力状态），即有极限值（高压应力状态），即25.1000)23(sin)23(cos/ctdcoctcoctdcoctcccfcfbacbaf参数参数b，当，当oct/fc=0时，时，b=oct/fc
35、即包络面或子午线与静水即包络面或子午线与静水压力轴交点的坐标；故压力轴交点的坐标；故b值为混凝土三轴等拉强度（值为混凝土三轴等拉强度（f1=f2=f3=fttt）与单轴抗压强度的比值与单轴抗压强度的比值 符合破坏曲面包络线随符合破坏曲面包络线随oct的增大由近似三角形趋向圆柱面过渡的增大由近似三角形趋向圆柱面过渡的特性；即，此时，拉、压子午线与静水压力轴平行切等距的特性；即，此时，拉、压子午线与静水压力轴平行切等距（rc=rt），偏平面上包络线为一半径），偏平面上包络线为一半径a的圆，破坏包络面趋于圆的圆，破坏包络面趋于圆柱形。柱形。maxcoctfactttffb 0d1.0时，时，=0o时
36、时c=ct，=60o时。时。c=cc，代人上式分别得拉、，代人上式分别得拉、压子午线，即为拉、压子午线对应的剪切强度。压子午线，即为拉、压子午线对应的剪切强度。当当=0o增加到增加到60o时，时，ct逐渐增加至逐渐增加至cc，符合光滑、外凸，符合光滑、外凸的特性；的特性；dcoctcoctcoctfcfbaf/其导数在其导数在 oct/fc=b处的数值为处的数值为，即切线垂直于横，即切线垂直于横坐标，拉、压子午线在此处连续，破坏包络面顶点处坐标，拉、压子午线在此处连续，破坏包络面顶点处连续、光滑；连续、光滑；25.1000)23(sin)23(cos/ctdcoctcoctdcoctcccfc
37、fbacbaf 另外，由于该破坏准则是根据包括整个应力空间另外，由于该破坏准则是根据包括整个应力空间8个个象限的各种应力状态的上千个试验点建立起来的，所象限的各种应力状态的上千个试验点建立起来的，所以它不仅在中、高静水压力区域实验值符合较好，而以它不仅在中、高静水压力区域实验值符合较好，而且在拉区乃至三向等拉状态也能较好地反映实际受力且在拉区乃至三向等拉状态也能较好地反映实际受力情况。情况。该准则适用于平面应力、平面应变、三向受压、三该准则适用于平面应力、平面应变、三向受压、三向受拉、乃至三向拉压等多种应力状态，且计算简单，向受拉、乃至三向拉压等多种应力状态，且计算简单，便于工程设计和非线性分
38、析应用。便于工程设计和非线性分析应用。25.1000)23(sin)23(cos/ctdcoctcoctdcoctcccfcfbacbaf 计算参数值的确定计算参数值的确定 混凝土破坏准则中包含的混凝土破坏准则中包含的5个参数，可以用全部试验个参数，可以用全部试验数据进行回归分析拟定，也可在破坏包络面上，或拉、数据进行回归分析拟定，也可在破坏包络面上，或拉、压子午线上选定任意压子午线上选定任意5个特征强度值加以标定。前者计个特征强度值加以标定。前者计算工作量大，一般取用后者。算工作量大，一般取用后者。单轴抗压和抗拉强度是混凝土的基本强度指标，应单轴抗压和抗拉强度是混凝土的基本强度指标，应作为首
39、选的二个特征强度值。其余作为首选的二个特征强度值。其余3个特征强度可以个特征强度可以选用：包络面顶端，即拉压子午线交点处的三轴等拉选用：包络面顶端，即拉压子午线交点处的三轴等拉强度；试验数量较多的二轴等压强度；和一个强度较强度；试验数量较多的二轴等压强度；和一个强度较高的常规三轴抗压强度高的常规三轴抗压强度(0 f1=f2 f3，=60o）。这）。这样使拉、压子午线上各有样使拉、压子午线上各有3个控制点，可以较好地拟个控制点，可以较好地拟合试验结果。合试验结果。将这将这5个特征值的应力状态分别代入式个特征值的应力状态分别代入式计算计算octoctoctfffarcrfffffffff232co
40、s 31)()()(31313321213232221321coctfcoctf并代人破坏准则计算式，并代人破坏准则计算式，可得可得5个联立方程如下：个联立方程如下：25.1000)23(sin)23(cos/ctdcoctcoctdcoctcccfcfbacbafdddddScSba TFcFbaF.FcFbaScSbaTcba0t00tc60c6060c 3/3/47140 0 3/13/10.4714 从这些方程求解从这些方程求解5个参数值，个参数值，难有显式解，可采用迭代难有显式解，可采用迭代法进行数值计算：法进行数值计算：由式直接得：由式直接得：dddddScSba TFcFbaF.
41、FcFbaScSbaTcba0t00tc60c6060c 3/3/47140 0 3/13/10.4714 ctttffFbc ffFffttttt、其中：其中：由其余由其余4式消去参数式消去参数a，有：，有：F.TScFcFbSbFcFcFbbTSccbSb471403/3/13/13/3/3/1 4714.03/13/1 0d0tt0dct60d60cc60由式由式F.TScFcFbSb471403/3/0d0tt0得参数得参数 d 的计算式：的计算式：3/log47140log3/3/log47140log0000tt0FbSbnF.TFbSbScFcF.Td13/0ttScFcn其中取
42、4714.03/13/1 60d60cc60TSccbSb由式由式1/1606060cc4714.03/13/1KTSbbSccd160113/1KSKcc取取得得由式由式FcFcFbb13/13/3/3/1 dct取取2/1ct13/13/3/13/KFbFbcFcd得得3312FKccct最后由式中任意一式计算参数最后由式中任意一式计算参数a，取式得：，取式得：dtFcFbFa3/3/4714.0 在设定了在设定了5个特征强度值后、即个特征强度值后、即S60、T60、S0、T0等值已知，等值已知，可应用这些方程进行迭代计算，以确定混凝土破坏准则的可应用这些方程进行迭代计算，以确定混凝土破坏
43、准则的5个参个参数值。其步骤如下：数值。其步骤如下：计算参数计算参数b；ctttffFb 设定设定n（1）的初始值，如）的初始值，如n0=0.98；代入代入计算参数计算参数d；3/log47140log3/3/log47140log0000tt0FbSbnF.TFbSbScFcF.Td由式由式代入代入1/1606060cc4714.03/13/1KTSbbSccd2/1ct13/13/3/13/KFbFbcFcd计算计算K1和和K2；由式由式160113/1KSKcc3312FKccct计算参数计算参数cc和和ct；代入代入13/0ttScFcn得得n的第一次近似值的第一次近似值n1，计算误，
44、计算误差，差，01nn 若不满足精度要求（取若不满足精度要求（取0.0001），），则按步骤继续迭代计算；则按步骤继续迭代计算；代入代入dtFcFbFa3/3/4714.0计算参数计算参数a。确定这确定这5个参数采用的混凝土特征强度值为：个参数采用的混凝土特征强度值为：单轴抗压（单轴抗压（-fc)；单轴抗拉单轴抗拉(ft=0.1 fc，F=0.1）；）；二轴等压二轴等压(fcc=1.28 fc，S0=-0.8533，T0=0.6034）；）；三轴等拉三轴等拉(fttt=0.9 ft，=0.9）；）；三轴抗压强度三轴抗压强度 （=60o，S60=oct/fc=4，T60=oct/fc=2.7）。
45、）。分别代入上式，用迭代法计算的参数值：分别代入上式，用迭代法计算的参数值：a6.9638 b=0.09 d=0.9297 ct12.2445 cc7.3319 按此公式可计算各种应力状态下的混凝土多轴强度按此公式可计算各种应力状态下的混凝土多轴强度理论值，并绘制子午线和偏平面包络线，以及二轴和理论值，并绘制子午线和偏平面包络线，以及二轴和三轴包络线。按此准则计算的混凝土多轴强度值与国三轴包络线。按此准则计算的混凝土多轴强度值与国内外的试验结果比较吻合。内外的试验结果比较吻合。将所得参数值代入基本方程，即得混凝土的破坏准将所得参数值代入基本方程，即得混凝土的破坏准则公式：则公式：25.1929
46、7.0)23(sin3319.7)23(cos2445.12/09.09638.6cfcffcoctcoctcoct 需要说明，选用的上述需要说明，选用的上述5个特征强度值，是分析了国个特征强度值，是分析了国内外众多研究者的试验结果而确定的，与此相应的混内外众多研究者的试验结果而确定的，与此相应的混凝土破坏准则（上两式）可适用于各种试验条件和全凝土破坏准则（上两式）可适用于各种试验条件和全部多轴应力范围，总体计算准确度较高。如果针对某部多轴应力范围，总体计算准确度较高。如果针对某一种特定的混凝土材料，或者在有限的应力比或静水一种特定的混凝土材料，或者在有限的应力比或静水压力范围（如二轴应力状态
47、）内，为了得到更准确的压力范围（如二轴应力状态）内，为了得到更准确的破坏准则，可以通过试验测定，或参照已有试脸资料破坏准则，可以通过试验测定，或参照已有试脸资料另行设定另行设定5个特征强度值，用上述迭代法计算参数值，个特征强度值，用上述迭代法计算参数值，得相应的破坏准则计算式。得相应的破坏准则计算式。25.19297.0)23(sin3319.7)23(cos2445.12/09.09638.6cfcffcoctcoctcoct4.7.4多轴强度验算举例多轴强度验算举例 二维和三维结构在线弹性或非线性分析后获得了混凝土的多二维和三维结构在线弹性或非线性分析后获得了混凝土的多轴应力状态，轴应力状
48、态，可按多轴强度设计值进行验算可按多轴强度设计值进行验算（如（如4.5所述），所述），也可采用破坏准则进行验算也可采用破坏准则进行验算，通常将混凝土的破坏准则编成程，通常将混凝土的破坏准则编成程序，附在结构分析之后，由计算机完成混凝土的应力分析和多序，附在结构分析之后，由计算机完成混凝土的应力分析和多轴强度验算。轴强度验算。下面列举几个手算例题，说明具体的计算方法和步骤，有助下面列举几个手算例题，说明具体的计算方法和步骤，有助于对混凝土破坏准则的理解。于对混凝土破坏准则的理解。例例4-7 混凝土三向受压，应力比为混凝土三向受压，应力比为1:2:3-0.15:-0.3:-1，用上述破坏准则计算相
49、应的多轴强度值。用上述破坏准则计算相应的多轴强度值。解：设三轴抗压强度为：解：设三轴抗压强度为：cxff3另二个方向分别为：另二个方向分别为：3.0 15.021ccxffxff其中其中 x 为待定值。为待定值。计算破坏准则的各项指标和参数值：任何类型、体系和受力状况的结构或其局部都可依靠非线性分析方法求解。应力不变量Il,J2,J3;式中 E1、E2、E3 3个相互垂直方向的弹性模量；Y和H合称为耦合刚度系数（模量）G 剪切模量（N/mm2）。静水压力轴上一点与坐标原点的距离称为静水压力（）；解此超越方程得：x=1.任何类型、体系和受力状况的结构或其局部都可依靠非线性分析方法求解。混凝土的单
50、轴受压和受拉应力-应变关系；至今国内外在混凝土结构的非线性有限元分析中常用的非线弹性本构模型有 Ottosen 模型和 Darwin Pecknold 模型等，有些国际规范也明确建议采用这两种模型。I1为应力张量ij的第一不变量12、应力22对方向1的横向变形系数（泊松比），其余类推。一些特殊应力状态的混凝土强度点，在破坏包络面上占有特定的位置。上两个本构模型都是以混凝土单轴受压应力-应变曲线方程为基础建立的，所得多轴应力-应变理论曲线与单轴和二轴受压类的试验曲线相符，但对三轴受压、受拉应力-应变曲线有不同程度的误差。取主应力轴正方向处为=0o，负方向处为=60o，其余各处为0o60o。最后由
