混凝土材料不仅应用于房屋建筑结构，还大量应用于重力坝、高拱坝、核电站预应力混凝土安全壳以及隧道衬砌等复杂结构.这些结构大都在动力和多向应力状态下工作，对其加-卸载响应的研究，对材料应力-应变关系的描述是更全面掌握此材料特性及进行定量分析的基础.混凝土材料具有的静水压力相关性、剪切压缩-膨胀现象、应变软化行为等特性导致其荷载响应复杂多样，建立合理的本构模型是进行混凝土结构非线性分析的基础.

本构模型的发展经历了增量理论^[1]，剪切应变和随动硬化理论^[2-3]，多面模型^[4-5]，边界面塑性^[6-8].其中边界面模型以其形式简单、应用方便、可描述复杂的三轴循环加载等优点被许多学者用来描述混凝土在循环荷载下的性能.Dafalia等^[6]和Krieg^[9]首次提出边界面的概念，并将其应用于金属及土体.Fardis等^[7]第一次用此模型来模拟混凝土的性能：模型不包含屈服面，方程为增量形式，仅可模拟塑性变形，不可模拟软化行为.Fardis等^[8]提出增量形式的边界面模型，没有屈服面，可描述混凝土的弹性和塑性应变及软化现象；Yang等^[10]提出增量形式的边界面模型，包含两个加载面，考虑等向损伤，可描述混凝土软化行为；Suaris等^[11]提出全量形式的边界面模型，考虑了异向损伤，但不能模拟混凝土软化行为.屠永清^[12]、王湛等^[13]对边界面模型进行了一定研究, 并用边界面模型模拟混凝土在单轴单调、单轴循环荷载作用下的应力-应变曲线、应力-体积应变曲线.

一般的边界面模型中，塑性模量的确定不依赖于硬化函数和一致性条件，仅根据边界面和屈服面的位置关系确定，相比数学逻辑关系严密的经典塑性力学模型是一缺陷.屈服面的形式选取是研究混凝土弹塑性本构模型的一个重要方向.通常屈服面在子午面内是开口的，或者盖帽组合式.前者模型表明屈服应力可随静水压力无限增长，这与实验结果不符；后者模型或存在角点问题，不便进行数值计算，或曲线在连接点处连续，但公式较多.任放等^[14]建议的蛋形屈服面能够很好弥补以上两类屈服面的缺点，此类模型参数相对较少，且其形式更适用于混凝土材料，但被应用于本构模型及数值化较少，仅有少量文献^[15-16]涉及.

借鉴边界面模型优点，建立一个适用于循环加载作用下的混凝土本构模型.子午面上的屈服面采用了更符合实际且不通过原点的梨形面，将有效塑性偏应变作为硬化参数.通过引入剪胀参数，将塑性偏应变增量与剪切引起的塑性体积应变增量联系起来，进而推导增量形式的边界面弹塑性本构关系，实现数值模拟.结合课题组有关偏平面上循环加载的试验数据^[17]，对模型进行验证.结果表明该模型能够较合理模拟混凝土在偏平面上的应力-应变关系.

1 弹性和塑性应变分解

本文模型中，假设材料是各向同性的，弹性和弹塑性响应都是率无关的.定义应力和应变以压缩为正，拉伸为负.

总应变增量dε_ij由弹性和塑性两部分组成

$ {\rm{d}}{\varepsilon _{ij}} = {\rm{d}}\varepsilon _{ij}^{\rm{e}} + {\rm{d}}\varepsilon _{ij}^{\rm{p}}, $

(1)

式中dε_ij^e和dε_ij^p分别表示弹性和塑性应变增量.

弹性部分由胡克定律表示为

$ {\rm{d}}\varepsilon _{ij}^{\rm{e}} = {\rm{d}}e_{ij}^{\rm{e}} + \frac{{{\rm{d}}\varepsilon _{\rm{v}}^{\rm{e}}}}{3}{\delta _{ij}} = \frac{{{\rm{d}}{s_{ij}}}}{{2G}} + \frac{{{\rm{d}}{\sigma _{kk}}}}{{9K}}{\delta _{ij}}. $

(2)

式中：de_ij^e为弹性偏应变增量；dε_v^e为弹性体积应变增量；ds_ij为偏应力增量；dσ_ij为应力增量；δ_ij为Kronecker张量；G为弹性剪切模量；K为弹性体积模量.G和K可由杨氏模量E和泊松比υ求出.

塑性部分表示为

$ {\rm{d}}\varepsilon _{ij}^{\rm{p}} = {\rm{d}}e_{ij}^{\rm{p}} + \frac{{{\rm{d}}\varepsilon _{\rm{v}}^{\rm{p}}}}{3}{\delta _{ij}}, $

(3)

式中de_ij^p为塑性偏应变增量，dε_v^p为塑性体积应变增量.

2 塑性部分

边界面塑性本构关系主要由屈服面方程、边界面方程、映射法则和内变量演化方程4个要素组成.

2.1 屈服面和边界面的形式

对于岩土、混凝土类工程材料，边界面在子午面内的形状有以下几种：直线^[19-20]、单个椭圆^[21-22]、由封闭逐渐开口的非均匀硬化加载面^[23-24]、椭圆-双曲线^[18]、椭圆-抛物线^[25]的多段曲线组合、水滴形状或蛋形^{[14-16, 26-27]}.

本文边界面模型包含一个由f=0决定的屈服面和一个由F=0决定的边界面.借鉴任放等^[14]研究成果，提出以下梨形函数作为屈服面函数：

$ f = {\left( {\frac{{{\sigma _{\rm{m}}} - c}}{a}} \right)^2} + {\left( {\frac{{1 - p}}{{1 + p\frac{{{\sigma _{\rm{m}}} - c}}{a}}}} \right)^2}{\left( {\frac{{\sqrt {{J_2}} }}{{bg\left( \theta \right)}}} \right)^2} - 1, $

(4)

而

$ \begin{array}{l} {\sigma _{\rm{m}}} = \frac{1}{3}{\sigma _{kk}} = \frac{1}{3}{I_1},\\ {J_2} = \frac{1}{2}\left[ {2G\left( {{e_{ij}} - e_{ij}^{\rm{p}}} \right) - {\alpha _{ij}}} \right]\left[ {2G\left( {{e_{ij}} - e_{ij}^{\rm{p}}} \right) - {\alpha _{ij}}} \right]. \end{array} $

式中：σ_m为平均应力，I₁为应力σ_ij的第一不变量；e_ij和e_ij^p分别为偏应变张量和塑性偏应变张量；J₂为偏应力s_ij的第二不变量；α_ij为屈服面在π平面上的中心；g(θ)为屈服面在π平面上的形状函数，θ为罗德角；p为子午面内屈服面形状的修正系数，当-1＜p＜0时，为剑桥弹头形曲面，当0＜p＜1时，为梨形曲面，当p=0时，为椭圆，此时a、b及c分别表示椭圆的长短半轴及在σ_m轴上的中心.

a、b、c和α_ij为内变量，在加载过程中会发生变化，其变化规律用微分方程描述，这些方程称为演化方程.其中a和b控制屈服面的大小，c和α_ij控制屈服面的位置.假设在变化的过程中a与b始终保持固定比例，记b/a=t₁.

通过改进摩尔-库伦破坏面函数，消除棱角影响，提出g(θ)的函数形式如下

$ g\left( \theta \right) = \frac{1}{{\left( {1 + \kappa } \right) + {{\left( {1 - \kappa } \right)}^2}{{\sin }^2}3\theta - \left( {1 - \kappa } \right)\cos 3\theta }}, $

(5)

其中κ为π平面上屈服面形状的修正系数.

边界面函数采用与屈服面函数相似的形式，取为

$ F = {\left( {\frac{{{{\bar \sigma }_{\rm{m}}} - \bar c}}{{\bar a}}} \right)^2} + {\left( {\frac{{1 - p}}{{1 + p\frac{{{{\bar \sigma }_{\rm{m}}} - \bar c}}{{\bar a}}}}} \right)^2}{\left( {\frac{{\sqrt {{{\bar J}_2}} }}{{\bar bg\left( {\bar \theta } \right)}}} \right)^2} - 1, $

(6)

式中σ_m、J₂、a、b、c和θ分别为边界面中对应的量，表达形式和含义均类似于以上屈服面中的定义.边界面在π平面上的中心用α_ij表示.

记$ {\mathit{h}_{\rm{0}}} = \frac{\mathit{b}}{{1 - \mathit{p}}}, {\mathit{\bar h}_{\rm{0}}} = \frac{{\mathit{\bar b}}}{{1 - \mathit{p}}}$，分别由式(4)和(6)得：当σ_m=c时，有$ \frac{{\sqrt {{\mathit{J}_{\rm{2}}}} }}{{\mathit{g}\left( \mathit{\theta } \right)}} = {\mathit{h}_{\rm{0}}}$；当σ_m=c时，有$ \frac{{\sqrt {{{\mathit{\bar J}}_{\rm{2}}}} }}{{\mathit{g}\left( {\mathit{\bar \theta }} \right)}} = {{\mathit{\bar h}}_{\rm{0}}}$.边界面和屈服面与子午面的交线见图 1.

2.2 映射法则

映射关系定义为：如图 1所示，A为屈服面上的当前应力点，A为边界面上一点，若两点有相同的单位法向向量，就认为A为A的映射点.

假设以下式子成立：

$ \theta = \bar \theta , $

(7)

$ \frac{{{\sigma _{\rm{m}}} - c}}{a} = \frac{{{{\bar \sigma }_{\rm{m}}} - \bar c}}{{\bar a}}, $

(8)

$ \frac{{\sqrt {{J_2}} }}{b} = \frac{{\sqrt {{{\bar J}_2}} }}{{\bar b}}. $

(9)

若屈服面上应力点对应的θ、σ_m、J₂及边界面上应力点对应的θ、σ_m、J₂满足式(7)~(9)，则两个点有相同的单位法向向量，即满足本文定义的映射关系.

定义距离d₁和d₂，其表达式分别如下：

$ {d_1} = {{\bar I}_1} - {I_1}, $

(10)

$ {d_2} = \sqrt {\left( {{{\bar s}_{ij}} - {s_{ij}}} \right)\left( {{{\bar s}_{ij}} - {s_{ij}}} \right)} , $

(11)

式中I₁为边界面上映射应力σ_ij的第一不变量.

为方便后续应用，引入与d₂相关的量化距离：

$ \delta = \frac{{{{\left( {{d_2}} \right)}_0} - {d_2}}}{{{{\left( {{d_2}} \right)}_0}}}, $

(12)

式中(d₂)₀为当前应力点由弹性区域刚到达屈服面时对应的d₂值.之后若持续加载，其保持不变；若卸载进入弹性区域，则下一次加载时重新计算(d₂)₀.

2.3 硬化准则

屈服面和边界面的大小分别由内变量a、b和a、b确定，位置分别由内变量c、α_ij和c、α_ij确定.

2.3.1 屈服面的内变量演化方程

假设在加载过程中，屈服面的大小和形状不发生变化，但位置能够变化.

构造屈服面内变量的演化方程分别为

$ {\rm{d}}a = 0, $

(13)

$ {\rm{d}}b = 0, $

(14)

$ {\rm{d}}c = {\rm{d}}\bar c + {\gamma _1}{d_1} \cdot {\rm{d}}\lambda + {\gamma _2}\left( {{\sigma _{\rm{m}}} - c} \right)\frac{{{\rm{d}}\left( {\bar a} \right) - {\rm{d}}\left( a \right)}}{a}, $

(15)

$ \begin{array}{l} {\rm{d}}{\alpha _{ij}} = {\rm{d}}{{\bar \alpha }_{ij}} + \frac{{{m_1}{{\left( {{d_2}} \right)}_0}\left( {{{\bar s}_{ij}} - {s_{ij}}} \right){\rm{d}}{e_{\rm{p}}}}}{{\left[ {1 + {m_2}{{\left( {{{\left( {{d_2}} \right)}_0}} \right)}^2}} \right]{{\left( {{d_2} + \mathit{\Delta} } \right)}^{0.5}}}} + \\ \;\;\;\;\;\;\;\;\left[ {2G\left( {{e_{ij}} - e_{ij}^{\rm{p}}} \right) - {\alpha _{ij}}} \right]\frac{{{\rm{d}}\bar b - {\rm{d}}b}}{b}. \end{array} $

(16)

显然式(13)和(14)能够保证屈服面的大小不变.式(15)和(16)是通过分析试验数据和研究已有的演化方程得到的.式(15)中，dλ为非负比例因子；γ₁和γ₂为正值材料常数.式(16)中，Δ为一个较小正值，用来避免分母为零，本文取Δ=10^-8；m₁和m₂为正值材料常数；de_p= $ \sqrt {{\rm{d}}\mathit{e}_{\mathit{ij}}^{\rm{p}}{\rm{d}}\mathit{e}_{\mathit{ij}}^{\rm{p}}} $为等效塑性偏应变增量.

2.3.2 边界面的内变量演化方程

假设在加载过程中，边界面的大小发生变化，边界面的形状不变，边界面在静水压力轴方向的位置发生变化，但是边界面在偏平面内的位置不发生变化.在加载过程中，b和a都发生变化，但是由于边界面形状不变，二者的比值为定值.这里取b/a=b/a=t₁.

假设a的变化一部分是由静水压力引起，另一部分由偏应变引起，取

$ {\rm{d}}\bar a = {\rm{d}}\left( {\bar \mu \cdot \bar r} \right), $

(17)

式中:μ为与静水压力相关的函数，r为与偏应变相关的函数.μ与r的演变公式分别为

$ {\rm{d}}\bar \mu = {l_1}\frac{{{\rm{d}}{I_1}}}{{\left| {{\rm{d}}{I_1}} \right|}}{\rm{d}}\lambda , $

(17a)

$ {\rm{d}}\bar r = \frac{{{l_2}\chi }}{{\chi {e_{\rm{p}}} + 0.1}}{\rm{d}}{e_{\rm{p}}}. $

(17b)

式(17a)中，l₁为正值材料常数.式(17b)中，l₂为负值材料常数；χ=χ₁e^χ₂I'₁, χ₁为正值材料常数，χ₂为负值材料常数，I'₁=I₁/f_c，f_c为混凝土的单轴抗压强度；e_p=∫de_p为等效塑性偏应变.

边界面的其余内变量的演化方程为

$ {\rm{d}}\bar b = {t_1}{\rm{d}}\bar a, $

(18)

$ {\rm{d}}\bar c = \left( {1 + {t_2}\ln \left( {1 + 35\left| {{{I'}_1}} \right|} \right)} \right){\rm{d}}\bar \mu , $

(19)

$ {\rm{d}}{{\bar \alpha }_{ij}} = 0, $

(20)

式中t₁和t₂为正值材料常数.

式(19)中，内变量c将沿着I₁变化的方向移动.通过选取合适的参数，同时结合内变量a的演化，可以描述混凝土材料在经历三向等压荷载历史后，单轴抗压强度会有所降低的特性.当进行三向等压加载时，a增大，c沿着I₁增大方向移动，边界面在高静水压力区一侧不断膨胀，而在另一侧收缩，这种变化形式比较类似于文献[27]提出的混凝土不均匀强(软)化.

2.4 塑性应变增量

对于混凝土材料，静水压力仅引起体积变化，由静水压力所产生的塑性体积应变，通常考虑一个塑性体积模量即可近似计算.而剪应力的变化既会产生剪应变，又会引起体积应变的发生.

对于剪切引起的塑性变形，采用非相关流动法则.把剪切引起的塑性应变增量的方向n_ij^p分解为两项

$ n_{ij}^{\rm{p}} = n_{ij}^{\left( 1 \right)} + n_{ij}^{\left( 2 \right)}. $

(21)

其中

$ n_{ij}^{\left( 1 \right)} = \frac{{\partial f}}{{\partial {\sigma _{ij}}}} - \frac{1}{3}\frac{{\partial f}}{{\partial {I_1}}}{\delta _{ij}}, $

(21a)

$ n_{ij}^{\left( 2 \right)} = \frac{{\omega \zeta }}{{{{\left( {100qe + 1} \right)}^{2.5}}}}{\delta _{ij}}. $

(21b)

在式(21a)中

$ \frac{{\partial f}}{{\partial {\sigma _{ij}}}} = \frac{{\partial f}}{{\partial {I_1}}}\frac{{\partial {I_1}}}{{\partial {\sigma _{ij}}}} + \frac{{\partial f}}{{\partial {I_2}}}\frac{{\partial {I_2}}}{{\partial {\sigma _{ij}}}} + \frac{{\partial f}}{{\partial g\left( \theta \right)}}\frac{{\partial g\left( \theta \right)}}{{\partial {\sigma _{ij}}}}. $

(22)

在式(21b)中

$ \begin{array}{l} \omega = \left( {{e_{11}} - {e_{22}}} \right)\left( {n_{11}^{\left( 1 \right)} - n_{22}^{\left( 1 \right)}} \right) + \\ \;\;\;\;\;\;\left( {{e_{22}} - {e_{33}}} \right)\left( {n_{22}^{\left( 1 \right)} - n_{33}^{\left( 1 \right)}} \right) + \\ \;\;\;\;\;\;\left( {{e_{33}} - {e_{11}}} \right)\left( {n_{33}^{\left( 1 \right)} - n_{11}^{\left( 1 \right)}} \right), \end{array} $

$ qe = \sqrt {\frac{1}{2}\left[ {{{\left( {{e_{11}} - {e_{22}}} \right)}^2} + {{\left( {{e_{22}} - {e_{33}}} \right)}^2} + {{\left( {{e_{33}} - {e_{11}}} \right)}^2}} \right]} . $

式(21)中，第一项n_ij⁽¹⁾描述的是剪切引起的塑性应变中的关联流动部分；第二项n_ij⁽²⁾为修正项，修正剪切引起的塑性应变与其关联流动部分的偏差.

式(21b)是通过分析体积应变ε_v与广义剪应变qe之间的关系曲线得到的.用ζ反映剪切应变对体积应变的影响，其表达式为

$ \zeta = {\zeta _1} + {\zeta _2}, $

(23)

其中

$ {\zeta _1} = {z_1}\left( {{k_1}{{I'}_1} + {k_2} - \delta } \right){{\rm{e}}^{ - {z_2}{{\left( {\delta - 0.9} \right)}^2}}}, $

$ {\zeta _2} = {z_3}\left( {qe - q{e_{\max }}} \right){{\rm{e}}^{ - {z_4}{{\left( {qe - q{e_{\max }}} \right)}^2}}}. $

式中：ζ₁反映剪切应变引起的体积压缩；ζ₂反映剪切应变引起的体积膨胀；z₁、z₂、z₃和z₄均为正值材料常数；qe_max为qe经历的最大值.

因此，剪切引起的塑性应变增量d(ε_ij^p)_s可表示为

$ {\rm{d}}{\left( {\varepsilon _{ij}^{\rm{p}}} \right)_{\rm{s}}} = n_{ij}^{\rm{p}}{\rm{d}}\lambda = {\rm{d}}e_{ij}^{\rm{p}} + \frac{{{\delta _{ij}}}}{3}{\rm{d}}{\left( {\varepsilon _{\rm{v}}^{\rm{p}}} \right)_{\rm{s}}}. $

(24)

其中

$ {\rm{d}}e_{ij}^{\rm{p}} = n_{ij}^{\left( 1 \right)}{\rm{d}}\lambda , $

(24a)

$ {\rm{d}}{\left( {\varepsilon _{\rm{v}}^{\rm{p}}} \right)_{\rm{s}}} = n_{kl}^{\left( 2 \right)}{\delta _{kl}}{\rm{d}}\lambda . $

(24b)

式中d(ε_v^p)_s表示剪切引起的塑性体积应变增量.

由静水压力引起的塑性体积应变增量可表示为

$ {\rm{d}}{\left( {\varepsilon _{\rm{v}}^{\rm{p}}} \right)_{\rm{h}}} = \frac{{{\rm{d}}{\sigma _{kk}}}}{{9{K^{\rm{p}}}}}, $

(25)

式中:K^p为塑性体积模量.

则总塑性体积应变增量表示为

$ {\rm{d}}\varepsilon _{\rm{v}}^{\rm{p}} = {\rm{d}}{\left( {\varepsilon _{\rm{v}}^{\rm{p}}} \right)_{\rm{s}}} + {\rm{d}}{\left( {\varepsilon _{\rm{v}}^{\rm{p}}} \right)_{\rm{h}}}. $

(26)

2.5 一致性条件

屈服面函数f的形式有3种：f=f₁(σ_ij, Z)，f=f₂(ε_ij, Z)，f=f₃(I₁, e_ij, Z).式中Z表示内变量的集合，本文中Z≡(e_ij^p, α_ij, a, b, c).

在整个塑性加载过程中，一致性条件成立.根据3种形式的屈服函数，相应有3种形式的微分方程df=0.这里只对第3种进行详细讨论.

结合式(13)~(20)和式(24a)，dλ可通过一致性条件得到.在整个塑性加载过程中，df₃=0，进一步有

$ \begin{array}{l} \frac{{\partial {f_3}}}{{\partial {I_1}}}{\rm{d}}{I_1} + \frac{{\partial {f_3}}}{{\partial {e_{ij}}}}{\rm{d}}{e_{ij}} + \frac{{\partial {f_3}}}{{\partial e_{ij}^{\rm{p}}}}{\rm{d}}e_{ij}^{\rm{p}} + \frac{{\partial {f_3}}}{{\partial {\alpha _{ij}}}}{\rm{d}}{\alpha _{ij}} + \\ \frac{{\partial {f_3}}}{{\partial a}}{\rm{d}}a + \frac{{\partial {f_3}}}{{\partial b}}{\rm{d}}b + \frac{{\partial {f_3}}}{{\partial c}}{\rm{d}}c = 0. \end{array} $

(27)

由式(27)可求得dλ，然后带入式(24a)、(24b)求出de_ij^p，d(ε_v^p)_s，进而得到增量应力-应变关系.

2.6 加卸载准则

为对塑性变形对应的加卸载进行判断，需引入与屈服函数f相关的加载函数L.

每一种形式的屈服面函数对应有一组加卸载判别公式，共有3组：

当f₁(σ_ij, Z)=0, L₁=$ \frac{{\partial {\mathit{f}_{\rm{1}}}}}{{\partial {\mathit{\sigma }_{\mathit{ij}}}}}{\rm{d}}{\mathit{\sigma }_{\mathit{ij}}}$;

当f₂(ε_ij, Z)=0, L₂=$ \frac{{\partial {\mathit{f}_{\rm{2}}}}}{{\partial {\mathit{\varepsilon }_{\mathit{ij}}}}}{\rm{d}}{\mathit{\varepsilon }_{\mathit{ij}}}$;

当f₃(I₁, e_ij, Z)=0, L₃=$ \frac{{\partial {\mathit{f}_{\rm{3}}}}}{{\partial {\mathit{I}_{\rm{1}}}}}{\rm{d}}{\mathit{I}_{\rm{1}}} + \frac{{\partial {\mathit{f}_{\rm{3}}}}}{{\partial {\mathit{e}_{\mathit{ij}}}}}{\rm{d}}{\mathit{e}_{\mathit{ij}}}$.

仅对f₃ 对应的加卸载准则进行讨论：

1) 当f₃＜0时，弹性状态，内变量不变；

2) 当f₃=0时，塑性状态，材料状态点在屈服面上，应用加卸载判别公式，有a)若L₃＞0，则为应力球量和应变偏量联合加载，使用式(13)~(20)来计算内变量的增量；b)若L₃≤0，则为卸载或中性变载，内变量增量为零.f₁、f₂对应的加卸载准则的判别方法与f₃的相似.

3 模型计算与试验数据比较 3.1 材料参数的选取

模型验证采用的是课题组在大连理工大学进行的混凝土三轴试验的部分数据^[17].经过对数据处理和试算，所有模拟采用的同一组参数见表 1.

除以上参数，模型中的参变量a、a、b、b、c和c需要赋予初值，分别为a₀=2 MPa，a₀=38 MPa，b₀=t₁a₀，b₀=t₁a₀，c₀=1.2 MPa和c₀=21.6 MPa.另K^p=2 620 MPa，f_c=35 MPa.

3.2 计算流程图

计算程序采用fortran语言编写.对静水压力加载部分，以平均应力驱动开始计算；对于偏应力加载部分，以偏应变驱动开始计算.前者计算过程较简单，仅给出偏应变驱动的计算过程流程图(见图 2).

3.3 模型验证

利用上述计算程序，结合课题组有关偏平面加载的部分试验数据与本文的计算结果进行比较验证.

试验是在若干个偏平面上沿压缩路径分别进行了单调加载和循环加载.加载分为两个阶段，先施加应力球量，然后再施加应力偏量，按照应力增量dσ_y=dσ_z=-dσ_x/2的约束关系改变荷载.本文采用的试验数据中，σ_m选取了3档：σ_m=30、45、60 MPa；单轴压缩时θ=0°，按照每档σ_m与θ的组合可分为3组加载路径.

模型计算结果与试验结果的比较见图 3~6.由于篇幅限制且试块的离散性较低^[17]，对于单调路径，选取每组路径的2块试验数据参与比较分析；对于循环路径，选取每组路径的1块试验数据参与比较分析.

从计算曲线来看，1)单调加载下，模拟得到的应力-应变曲线在两个阶段的表现都与试验曲线比较吻合(见图 3(a)~(b))；2)多轴加载下，无论单调或者循环应力状态，由本文模型计算得到的偏应力-偏应变曲线，其变化规律及趋势与试验结果基本一致(见图 3(b)、4(a)~(b)、5(a)~(b)和6(a)~(b))，循环加载曲线的滞回环宽度较小，应是某些方程形式的确定不合适；3)体积应变的变化与试验曲线在低围压时偏差较多，高围压时更吻合，是由于选取的参数ζ的表达式不合适，但是基本趋势是正确的(见图 3(c)、4(c)、5(c)和6(c)).总体上模型预测结果和试验结果是比较接近的.

4 结论

提出一种可描述材料循环加载的边界面塑性本构模型，其中的屈服面和边界面是封闭且不过原点的梨形曲面.引入的映射法则，可以确保当前应力点和映射应力点的法线方向始终相同.模型的加卸载判据中的微分变量是应力的第一不变量和应变偏量，根据它们的增量来判断加卸载，模型能够模拟材料的应变-软化行为.在考虑剪切引起塑性变形时，采用了非相关流动法则.

把模型的计算结果与试验结果进行了比较，模型能够模拟混凝土在多轴加载后的多种特性.