混凝土受压应力-应变全曲线方程(描述)

混凝土受压应力－应变

全曲线方程

混凝土受压应力－应变全曲线方程

混凝土的应力－应变关系是钢筋混凝土构件强度计算、超静定结构内力分析、结构延性计算和钢筋混凝土有限元分析的根底，几十年来，人们作了广泛的努力，研究混凝土受压应力－应变关系的非线性性质，探讨应力与应变之间合理的数学表达式，1942年，Whitney通过混凝土圆柱体轴压试验，提出了混凝土受压完整的应力应变全曲线数学表达式，得出了混凝土脆性破坏主要是由于试验机刚度缺乏造成的重要结论，这一结论于1948年由Ramaley和Mchenry的试验研究再次证实，1962年，Barnard在专门设计的具有较好刚性且能控制应变速度的试验机上，试验了一批棱柱体试件以及试件两靖被放大的圆柱体试件，试验再次证明，混凝土的突然破坏并非混凝土固有特性，而是试验条件的结果，即混凝土的脆性破坏可用刚性试验机予以防止，后来由很多学者〔如M.Sagin，P.T.Wang，过镇海等〕所进行的试验，都证明混凝土受压应力－应变曲线确实有下降段存在，那么混凝土受压应力与应变间的数学关系在下降段也必然存在，研究这一数学关系的工作一刻也没有停止。

钢筋混凝土结构是目前使用最为广泛的一种结构形式。但是，对钢筋混凝土的力学性能还不能说已经有了全面的掌握。近年来，随着有限元数值方法的开展和计算机技术的进步，人们已经可以利用钢筋混凝土有限元分析方法对混凝土结构作比拟精确的分析了。由于混凝土材料性质的复杂性，对混凝土结构进行有限元分析还存在不少困难，其中符合实际的混凝土应力应变全曲线确实定就是一个重要的方面。

1、混凝土单轴受压全曲线的几何特点

经过对混凝土单轴受压变形的大量试验大家一致公认混凝土单轴受压变过程的应力应变全曲线的形状有一定的特征。典型的曲线如图1所示，图中采用无量纲坐标。

s

c c E E N f y x 0,,===

σ

εε 式中，c f 为混凝土抗压强度；c ε为与c f 对应的峰值应变；0E 为混凝土的初始弹性模量；s E 为峰值应力处的割线模量。 此典型曲线的几何特性可用数学条件描述如下： ①x=0，y=0； ②0≤x<1，

2

2

x y

d d <0，即上升段曲线

x

y d d 单调减小，无拐点；

③C 点x=1处，

x

y d d =0和y=1.0，曲线单峰；

④D 点

2

2x y

d d =0处坐标x D >1.0，即下降段曲线上有一拐点；

⑤E 点

3

3x y

d d =0处坐标x E 〔≥x D 〕为下降段曲线上曲率最大点；

⑥当x →∞，y →0时，x

y d d →0；

⑦全部曲线x ≥0，0≤y ≤1.0。 这些几何特征与混凝土的受压变形和破坏过程完全对应，具有明确的物理意义。

2、混凝土单轴受压曲线

方程的比拟和分析

对于混凝土在单轴受压下的应力应变关系，已经做了大量的试验研究工作，在此根底上不少学者提出了多种混凝土受压应力应变曲线方程。 (1) Hongnestad 的模型

0.002

0.0038

f

c

f c

模型的上升段为二次抛物线，下降段为斜直线。

上升段：22

00022,x x y f c -=⇒⎥⎥⎦

⎤

⎢⎢⎣

⎡⎪⎪⎭

⎫

⎝⎛-=≤εε

εεσεε 下降段：115.085.015.01,0---=⇒⎥

⎦⎤⎢⎣⎡---=≤≤u u o u o c u x x

x y f εεεεσεεε 式中，c f ——峰值应力〔棱柱体抗压强度〕；

0ε——相应于峰值应力时的应变，取0ε＝0.002； u ε——极限压应变，取u ε＝0.0038。 混凝土受压应力应变曲线上升段，对x 求一阶导数：

x y 22-='

当x ＝1时，y '＝0；当x ＝0时，y '＝2。很容易得出曲线满足典型曲线的条件③。在Hongnestad 公式中y '＝2是一个固定值，所以Hongnestad 公式只能在工程上作为一个近似公式使用。

对x 求二阶导数，得：

2-=''y

Hongnestad 公式满足条件②。受压应力应变曲线下降段的形状，更敏感地反映混凝土的延性和破坏过程的缓急，以往的曲线公式都不能很好的反映混凝土受压应力应变曲线的下降段，Hongnestad 公式不满足典型曲线下降段的要求。

Hongnestad 的模型一般可以作为钢筋混凝土简支梁的实例分析，采用三维模型，对矩形截面钢筋混凝土简支梁进行模拟分析。梁单元类型采用ANSYS 中的6面体8节点单元。在ANSYS 中需要输入的物理参数有弹性模量E 和泊松比μ，参考?混凝土结构设计标准?(GB50010--2002)规定的材料力学指标的标准值，查得相应的取值，对混凝土简支梁进行数值分析。Hongnestad 的模型已经纳入CEB-FIP MC90等混凝土结构设计标准。

(2) Saenz 的模型

表达式：2

)2(1x

x N Nx

y +-+=

在混凝土应力应变曲线上升段需要满足条件①②③⑦，显然Saenz 公式

满足条件①⑦。下面看是否满足②③。 上升段曲线对x 求一阶导数得：

[]

2

2)2(1x x N Nx

N y +-+-=

'

容易得：N y y x x ='='==01,

0，满足条件③。

Saenz 公式的s

x E E N N y 0

0,=='=其值对于不同强度的混凝土是变化的。 曲线对x 求二阶导数：

[][]

3

223)2(1)

2(26)2(8)84(22x x N N N Nx x N N N N Nx y +-+------+=

''

那么：21

01840)2(2N N y N N y x x --=''--=''==

因为s

E E N 0

=

显然N>1且N 的值是变化的，对于Saenz 公式只有2≥N 时条件②才满足，所以只有当2≥N ，即混凝土的初始弹性模量和峰值割线模量的比值大于等于2时，采用Saenz 公式才是适宜的。当N 小于2时，Saenz 公式那么不能成立。实际应用中，当遇到这种情况时，总是强令N=2，这样处理显然是不合理的。同时Saenz 公式不能反映强度等级低的混凝土峰值局部 比强度等级高的混凝土峰值局部更为扁平这一事实。即不能满足特征⑧。

在工程应用中，Saenz 公式就可以作为FRP 约束混凝土应力应变的曲线模型，进行建模分析。Saenz 基于Pantazopoulou 的研究成果，引入体积应变v ε。

l c l c v εεεεεεθ2+=++=

式中，θεεε,,l c 分别为轴向应变、横向应变和环向应变，对于圆柱体，

θεε=l ，规定压应变为正，拉应变或者膨胀应变为负。

将FRP 约束混凝土应力应变曲线分成3段：

)0005.00(≤≤=c c c c E εεσ )00206.00005.0(sec ≤≤=c c

c E εεσ

)00206.0()00206.0(0,cc c c ct v c c E f εεεσ≤≤-+=

式中，co

c f E '=5700；0,v c f 为体积应变为0时的轴向应力； )00206.0()(0,--'=cc v c cc ct f f E ε。

在第二阶段约束混凝土轴向应变与横向应变的关系为

()00206.00005.000156.00005.0000618.02.02

≤≤⎪

⎭

⎫

⎝⎛---=c c c l εεεε

割线模量为β

εA c

E E +=11

sec ，式中，A ε为面积应变，对于圆柱体，

l A εε2=；β为割线模量软化率，310)44.114.3(-⨯+-=p μβ，其中p μ为极限 横向应变与轴向应变比值绝对值。

2

1C

le

cc

lu p K C -=-

=εεμ 式中，21,C C 为参数，分别取6.21和0.63；rup f f lu ,εξε=；le K 为FRP ，侧向有效刚度，0

2c f f le f D E t K '=

。本模型先通过式2

1C

le

cc

lu p K C -=-

=εεμ计算cc ε，再根据式c c E εσsec =计算cc

f '。 上述模型是在FRP 约束混凝土应力应变关系双直线特征的根底上建立的分段式模型，它回避了FRP 约束力的变化过程，极大简化了计算过程，适用范围较广，但它的精度受峰值点或极限点应力、应变的计算影响较大，且

没有明确的物理含义。

(3) 清华大学过镇海教授提出的模型

过镇海教授提出的应力应变全曲线模型为两段式模型。

)1(1

)1(10)2()23(232⎪⎩⎪

⎨

⎧≥+-≤≤-+-+=x x

x x x x a x a ax y α

式中，α,a 分别为使用年限t 的函数。

由公式中参数α,a 的物理意义可知：a 值小和α值大，那么曲线陡，曲线下的面积小，说明此混凝土的塑性变形小，剩余强度低，破坏过程急速， 材质较脆，接近于使用年限长的混凝土；反之，a 值大和α值小，那么混凝土变形大，剩余强度较高，破坏缓慢延性较好，适用于使用年限短的混凝土。本着这样原那么，将公式的混凝土应力应变曲线上升段、下降段与试验所测的不同使用年限的既有混凝土的应力应变全曲线上升段、下降段分别相比拟，选取一个吻合程度最好的值，具体数值见表l 。

根据a 、α值与使用年限t 的关系，对其进行非线性拟合，可由以下公式确定：

)2(6.193.0)46

.261

(-

+=e

a

)3(75.007.11

)76

.131(-+=

e

α

图1 参数a 与使用年限的关系 图2 参数α与使用年限的关系

图1、图2表示参数a 、α试验值和理论计算值的比拟，吻合程度较好。 这样，将a 值和α值直接代入式〔1〕，就可以得到不同使用年限既有未碳化混凝土的应力应变曲线方程，如式〔4〕所示：

)76

.131()46

.261

(23

275.007.11

)

4(6.193.01

)1(10)2()23(-

-+=

+=⎪⎩⎪

⎨

⎧≥+-≤≤-+-+=e

e a x x x x x x a x a ax y αα

式中，t 为混凝土的使用年限。

图3 不同使用年限混凝土应力应变曲线试验值与计算值的比拟

同强度养护28d 的新混凝土和既有混凝土的试验平均应力应变曲线与

按式〔4〕计算的理论曲线的比拟如图3所示，试验曲线与理论曲线吻合得较好。

清华大学叶知满对掺F 矿粉或粉煤灰高强混凝土应力应变全曲线试验研究时，对下降段曲线采取的就是过镇海教授的模型。

x

x d x

y +-=

2)1(

式中，d x f y c c

,,εεσ

==

——下降段参数，经统计可得d 与c f 关系式为

〔参图4〕 24)1026.7(c f d -⨯=

图4 下降段参数d 随c f 变化关系

图5给出了理论方程与实测曲线的比拟，可知理论方程与实测曲线吻合较好。

图5 理论曲线与实测曲线的比拟

3、结束语

建立混凝土轴压应力应变全曲线的数学模型，首先要弄清楚应力应变全 曲线的几何特点，观察和分析实测应力应变全曲线，通过与典型试验曲线的比拟，分析Hongnestad 公式、Saenz 公式和过镇海提出的公式在混凝土受压应力应变曲线上升段、下降段的适用范围，以及各自的拟合情况。Hongnestad 公式在上升段拟合较好，但下降段不分强度大小，斜率一律为－0.15，与典型试验曲线相差太大，不切实际，不适用；Saenz 公式在上升段曲线存在明显的问题，在s E E 0偏离2较多时，不能很好反映混凝土单向

受压试验所表现出的全部特征；清华大学提出的公式在上升段和下降段的拟合都较好，建议采用。较好的应力应变表达式，首先应该符合实测的应力应变曲线，由于它的影响因素较多，且有相当的分散性，所以公式应尽可能有较大范围的适用性及灵活性。

参考文献

[1] 叶列平. 混凝土结构[M]. 清华大学出版社，2005

[2] 过镇海，时旭东. 钢筋混凝土原理和分析[M]. 清华大学出版社，2003

[3] 曹居易. 混凝土的应力应变关系[J].

[4] 李义强,王新敏,陈士通. 混凝土单轴受压应力-应变曲线比拟[J]. 公路交通科技,2005

[5] 邓宗才，李建辉. FRP 约束混凝土应力应变曲线模型研究[J].

应用根底与工程科学学报，2021

[6] 叶知满. 掺F 矿粉或粉煤灰高强混凝土应力应变全曲线试验研究[J].

青岛建筑工程学院学报，1996

[7] 袁发顺，胡夏闽. 混凝土应力应变上升段曲线的一个最正确模式[J]. 南京建筑工程学院学报，1995

[8] 孟丽岩，王凤来，王涛. 既有混凝土应力应变曲线方程的研究[J].

[9] 刘世华. 混凝土全应力应变曲线分析[J]. 沈阳大学学报，2000