混凝土的开裂有限元分析

裂纹尖端应力场

应力强度因子

弹性理论得到的裂纹尖端应力
σx = σy = τ xy =
KI 2πr K II
θ θ 3θ cos (1 − sin sin ) 2 2 2
θ θ 3θ cos (1 + sin sin ) 2 2 2 2πr θ θ 3θ cos (sin cos ) 2 2 2
清华大学研究生课程——《钢筋混凝土有限元》
开裂单元
多裂缝模型



当主应力方向和裂缝方向差距较大时，由于裂 面剪力锁死往往导致出现错误结果 多裂缝模型认为当主应力方向和裂缝方向夹角 超过一定范围θ后，原有的裂缝闭合，重新在 新的主应力方向生成新的裂缝 θ ＝30度（6裂缝模型）45度（4裂缝模型）
清华大学研究生课程——《钢筋混凝土有限元》
清华大学研究生课程——《钢筋混凝土有限元》
清华大学研究生课程——《钢筋混凝土有限元》
固定裂缝模型与转动裂缝模型

固定裂缝模型


分布裂缝模型中，一般认为，当单元内 部的最大拉应力达到开裂应力时，混凝 土即开裂 混凝土开裂后，改变混凝土材料为各向 异性材料，主应力和主应变方向可能不 再一致。同时，初始裂缝方向和主应力 方向也不再一致，裂缝表面将出现剪应 力
分离裂缝模型的具体步骤

开裂标准和裂缝发展方向
主拉应力 虚拟裂缝模型 断裂力学方法
开裂标准和裂缝发展方向 裂缝发展与模型网格调整 裂面行为
清华大学研究生课程——《钢筋混凝土有限元》
清华大学研究生课程——《钢筋混凝土有限元》

实际裂缝
虚拟裂缝
清华大学研究生课程——《钢筋混凝土有限元》
混凝土的开裂与裂缝处理

脆性材料与半脆性材料



混凝土的一个重要特点是它在较低的应力水平 下就会开裂，且很多混凝土结构都是带裂缝工 作的 开裂后的混凝土其力学行为与未开裂混凝土有 很大不同，能否正确模拟开裂后的混凝土是混 凝土有限元分析中的关键问题 混凝土中大量的裂缝是对基于连续体力学的有 限元方法的一个重要挑战
裂纹尖端的应力趋向于无穷大 应力强度理论已经不再适合 引入断裂强度因子等概念来描述裂纹尖 端附近的应力场
K III 2πr
r->0，应力->无穷
清华大学研究生课程——《钢筋混凝土有限元》
清华大学研究生课程——《钢筋混凝土有限元》
应力强度因子计算

裂缝发展判据

无限大平板
K I = lim[ 2πr (σ y )θ =0 ] = σ πa
1 2 ∆ = (1 − ν 2 ) − ν 2 (1 + ν ) Ec Et
清华大学研究生课程——《钢筋混凝土有限元》
清华大学研究生课程——《钢筋混凝土有限元》
忽略开裂条件下泊松比
Et Dcr ' = 0 0 1 Ec νEc 0 Ec 1 −ν 2 0 0 0 0 0 0 Ec 2(1 + ν ) 0 0 0 0 ηEc 2(1 + ν ) 0
清华大学研究生课程——《钢筋混凝土有限元》
清华大学研究生课程——《钢筋混凝土有限元》
断裂力学的基本理论 混凝土的开裂有限元分析
江见鲸 陆新征 清华大学土木工程系 2004

基本开裂形式

张开型(I) 滑开型(II) 撕开型(III)

复合裂缝
清华大学研究生课程——《钢筋混凝土有限元》
清华大学研究生课程——《钢筋混凝土有限元》
ηE c 2(1 + ν )
0

为了防止过大的负刚度导致计算不收敛，一般 整体迭代计算不使用切线刚度矩阵，一般加以 适当修正
清华大学研究生课程——《钢筋混凝土有限元》
清华大学研究生课程——《钢筋混凝土有限元》
坐标转换矩阵

评述


混凝土开裂计算时需要进行坐标转换，在裂缝 方向建立应力应变关系。在求解单元刚度矩阵 的时候还要转回到整体坐标系中去。 整体坐标系下刚度矩阵 [ Dcr ] = [ R]T [ Dcr ' ][ R]
评述

分布裂缝模型

分离裂缝模型的虽然历史悠久、概念清晰，但 是由于混凝土结构中开裂问题的复杂性，以及 网格重划分技术的限制，目前主要用于分析有 少量裂缝的素混凝土结构，分析的问题也大部 分是平面问题。在大坝、岩石等领域有着较多 的应用。在通用有限元程序DIANA中，也集成 了这种裂缝模型。
分布裂缝模型的实质之将实际的混凝土裂缝“弥散”到 整个单元中，将混凝土材料处理为各向异性材料，利 用混凝土的材料本构模型来模拟裂缝的影响。这样一 来，当混凝土某一单元的应力超过了开裂应力，则只 需将材料本构矩阵加以调整，无需改变单元形式或重 新划分单元网格，易于有限元程序实现，因此得到了 非常广泛的应用。
应力应变关系和迭代矩阵

裂缝坐标系下的应力增量为：
∆ε 11 ∆σ 11 ∆ε ∆σ 22 22 ∆ε 33 ∆σ 33 = Dcr ' ∆γ 12 ∆σ 12 ∆γ 23 ∆σ 23 ∆γ 31 ∆σ 31
m12 m m
2 2 2 3
l12 2 l2 l2 [ R] = 3 2l1l 2 2l l 23 2l3l1
n12
2 n2 2 n3
l1m1 l 2 m2 l3 m3 l1m2 + l 2 m1 l 2 m3 + l3 m2 l3 m1 + l1m3
m1n1 m2 n 2 m3 n3 m1n2 + m2 n1 m 2 n3 + m 3 n 2 m3 n1 + m1n3
清华大学研究生课程——《钢筋混凝土有限元》
清华大学研究生课程——《钢筋混凝土有限元》
初始开裂判断

裂缝坐标下的应力应变矩阵
(1 −ν 2 ) ∆ Dcr ' = (1 + ν ) ∆ 1 Ec ( −ν 2 ) ∆ Et
开裂前，认为混凝土为弹性材料


认为混凝土开裂后，裂缝方向将永远保 持不变，即R矩阵始终不变 计算过程中首先将应力应变通过R矩阵转 换到裂缝坐标系下，计算此时的正应力、 正应变和剪应力、剪应变，以及裂缝坐 标系下的本构矩阵。迭代求解得到新的 荷载步的应力应变关系。 一个积分点最多出现3条彼此垂直的裂缝
清华大学研究生课程——《钢筋混凝土有限元》
该模型被称为共轴转动裂缝模型
清华大学研究生课程——《钢筋混凝土有限元》
清华大学研究生课程——《钢筋混凝土有限元》
开裂混凝土受剪行为
3.5 3.0 Shear stress (MPa) 2.5 2.0 1.5 1.0 0.5 0.0 0 1000 2000 3000 Shear strain (µε) 转动裂缝模型 多裂缝模型 固定裂缝模型
清华大学研究生课程——《钢筋混凝土有限元》
裂缝发展与模型网格调整

预设裂缝


由于分离裂缝模型是使用单元边界来模拟裂缝， 因此随着裂缝的发生和发展，需要不断调整单 元网格划分。 对于裂缝路径明确的问题，可以采用预设裂缝 的方法来模拟。 分离裂缝模型多用于分析只有一条或几条关键 裂缝的素混凝土或少筋混凝土结构，例如水坝 等
ν
(1 + ν ) ∆ ν Ec ( +ν ) ∆ Et 1 Ec ( −ν 2 ) ∆ Et
ν
0 0 0
0 0 0 0 Ec 2(1 + ν )
sym
ηEc 2(1 + ν )

当主拉应力达到开裂应力时，在垂直主 拉应力的方向出现第一条裂缝
0 0 0 0 ηE c 2(1 + ν ) 0
清华大学研究生课程——《钢筋混凝土有限元》
处理裂缝的主要方式

分离裂缝模型

分离裂缝模型 分布裂缝模型 内嵌裂缝单元模型
将裂缝处理为单元边界，一旦出现裂缝就调整 节点位置或增加新的节点，并重新划分单元网 格，使裂缝处于单元边界与边界之间。
裂缝尖端
清华大学研究生课程——《钢筋混凝土有限元》
清华大学研究生课程——《钢筋混凝土有限元》
2m1m2 2m2 m3 2m3 m1
2n1n2 2 n 2 n3 2n3 n1
n1l 2 + n2l1 n2 l3 + n3l 2 n3l1 + n1l3 n1l1 n2l 2 n3 l 3

分布裂缝模型因为将单个的裂缝连续化，不需 要改变有限元网格划分，特别适用于有限元分 析，因而在各个有限元软件中都得到了广泛的 应用，各种非线性有限元软件都或多或少采用 了分布裂缝模型 对于正常配筋的大型混凝土构件，分布裂缝模 型在一定程度上更接近工程实际情况（比分离 裂缝模型），分析结果也较好
网格依赖性

在长期应用与研究中发现，分布裂缝模型裂缝 比较容易沿着网格划分方向发展，这种现象被 称为网格依赖型。
清华大学研究生课程——《钢筋混凝土有限元》
清华大学研究生课程——《钢筋混凝土有限元》
单元尺寸效应
F
L
N 个单元
F
N=10
N=1
清华大学研究生课程——《钢筋混凝土有限元》
清华大学研究生课程——《钢筋混凝土有限元》
r →0
K II = τ πa

材料的断裂韧度 当应力强度因子大于材料的断裂韧度时， 裂缝将扩展 用途

K III = τ l πa