ABAQUS混凝土弥散开裂模型Standard模块

1、.混凝土弥散开裂模型程序： Abaqus/StandardAbaqus/CAE参考1、“材料库：概述”章2、“非弹性行为”章3、*混凝土4、*拉伸硬化5、*剪力传递6、*破坏面比率7、“定义混凝弥散开裂”in“定义塑性”章概述Abaqus/Standard 中的混凝土弥散开裂模型：1、为所有结构中的混凝土提供了普遍建模功能,包括梁，桁架，壳体和固体；2、可以用于素混凝土，虽然其主要用于钢筋混凝土；3、可以与钢筋（ rebar ）一起使用来模拟混凝土加筋；4、被设计用于混凝土在低围压下的基本单调张拉应变的情况；5、由各向同性硬化屈服面（压应力为主导时激活）和一个独立的“裂缝检测面”（由某点的破坏

2、是否由于开裂导致所决定）组成。6、用非各向同性的损伤弹性概念（弥散开裂）来描述开裂破坏后材料反应的可Word 资料.逆部分7、要求使用线弹性材料模型来定义弹性性能；并不能和局部方向一起使用。 (see“Orientations, ”Abaqus 中可以使用的混凝土模型的讨论见“Inelastic behavior, ”加筋混凝土结构的增强材料以钢筋作为典型代表，它是单轴应变理论单元（杆件），它可以单独定义或者嵌入有方向的表面。钢筋主要与塑性模型一起使用来描述钢筋材料的行为，并叠加在用来模拟混凝土的标准单元类型网格上。使用这种建模方法，混凝土材料行为独立于钢筋材料。钢筋/ 混凝土的接触面连结效应

3、，如粘结滑移和销钉作用，是通过在混凝土模型中引入某种“拉伸硬化”来近似模拟裂缝处通过钢筋的荷载转移。在复杂问题中定义钢筋会是冗长的，但钢筋的准确定义是非常重要的，因为在模型的关键位置缺少钢筋会引起分析失败。See “Definingreinforcement, ”开裂此模型的本意是用于在相当低的围压下（小于混凝土单轴受压极限的45个量级）的相对单调荷载下的混凝土行为。裂缝检测Word 资料.开裂被认为是力学行为的最重要方面，它以开裂和开裂后行为为代表，并主导控制模型。开裂被认为当应力达到被称为“裂纹探测面”的破坏面时出现。这个破坏面是一个等效压应力p 和 Mises 等效偏应力 q 之间的线性

4、关系，如图所示。当检测到裂缝时， 裂缝的方向被保存用于随后的分析。由于开口裂缝的应力分量不被包含在用于检测新增裂缝破坏面的定义中，同一点的再次开裂被限制为垂直于原开裂方向。裂缝是不能恢复的：他们在分析的过程中一直存在(但是可以开口或闭合 )。任意一点只能产生不多于3 条裂缝（在平面应力状态下为2 条，单轴应力状态下为 1 条）。检测到裂缝后，裂缝使用损伤弹性模型来影响计算。非各向同性的损伤弹性材料的更多详细讨论见节。弥散裂缝混凝土模型是弥散开裂模型，这意味着它不能模拟单独存在的“巨大”裂缝。有限单元模型在每个积分点会各自进行基本计算。通过裂缝影响积分点处的应力和材料刚度的方法，裂缝的存在作用进

5、入了计算中。拉伸硬化贯穿裂缝直接张拉产生的后破坏行为使用拉伸硬化来模拟，它允许用户定义开裂混凝土的应变软化行为。 这一行为也可以被用于以简单的方法模拟混凝土和钢筋的相互作用效应。 混凝土弥散开裂模型中需要定义拉伸硬化。用户可以依据后继破坏应力应变关系的方法或应用断裂能开裂准则来指定拉伸硬化。Word 资料.破坏后的应力应变关系钢筋混凝土的应变软化行为准则通常意味着指定贯穿裂缝的破坏后的应力应变函数。当没有钢筋或者其数量很少时，这种准则往往引起分析结果的网格敏感性，这意味着当网格被细化后，由于网格细化会导致更窄的裂缝带，有限元预测不能收敛到唯一解。 这种问题特别发生在结构中仅有少数离散裂缝并且网

6、格细化没有引起额外的裂缝的情况。如果裂缝是均匀分布的 （由于钢筋效应或由于稳定弹性材料的存在导致的，如板弯曲情况），网格敏感性则不值得关注。在钢筋混凝土的实际计算中， 网格通常是每个单元都包含钢筋。钢筋和混凝土的相互作用有利于减小网格敏感性，只要在混凝土中提供合理的拉伸硬化参数来模拟这种相互作用（图）。图“拉伸硬化”模型必须估计拉伸硬化的影响； 这依赖于钢筋密度、 钢筋和混凝土的粘结、 混凝土骨料对钢筋直径的相对大小、网格等因素。 对于在较细网格下的大量混凝土材Word 资料.料，当总应变10 倍于破坏应变时，对于相当重的网格过细的钢筋混凝土模型，合理的起点被假定为破坏后应变软化使应力在破坏应

7、变的10 倍的范围内线性降低为 0。标准混凝土破坏时的应变为 410 量级，这意味着拉伸硬化使应力在总应 3变大约为 10的时候降为 0 是合理的。对特殊情况，必须校准这个参数。Abaqus/Standard 中拉伸硬化参数的选择是十分重要的，一般说来，越大的拉伸硬化越容易得到数值解。 太小的拉伸硬化将会导致混凝土局部裂缝破坏从而在模型的总体响应中引入临时性的不稳定行为。少数实际设计中存在以上行为，因此分析模型中存在这类响应通常表明拉伸硬化取值不合理的较低。InputFileUsage:Use both of the following options:*CONCRETE*TENSION ST

8、IFFENING,TYPE=STRAIN (default)Abaqus/CAEUsage:Property module: material editor:MechanicalPlasticityConcrete Smeared Cracking:SuboptionsTension Stiffening: Type: Strain断裂能开裂准则如之前所讨论的， 当混凝土模型的重要部分没有钢筋的时候，定义拉伸硬化的应变软化方法可能在结果中引入不合理的网格敏感性。Crisfield (1986)讨论了这一问题并且给出结论，他认为 Hillerborg (1976) 的提议对许多实际问题足够精确。

9、 Hillerborg 用脆性Word 资料.断裂概念定义单位面积开裂产生裂缝所需的能量作为一个材料参数。使用这种方法，混凝土的脆性行为以应力- 位移响应为特性， 而不是应力 - 应变响应。 受拉时混凝土试件会在某些截面出现贯穿裂缝。当构件被充分张拉后， 绝大多数应力消失（因此弹性应变很小），试件的长度主要由裂缝的开口大小决定。开口的大小不依赖于试件的长度 (图。图断裂能开裂模型应用在有限元模型中应用应力 位移相关的概念，需要定义与质点相关的特征长度。裂缝特征长度是基于单元的几何特征和公式表述来确定的：对于一阶单元，该长度是穿过一个单元的线段的典型长度，对于二阶单元是穿过单元线段的典型长度的

10、1/2 。对于梁单元和桁架单元，该特征长度是单元轴向的特征长度，对于膜单元和壳单元则是参考面的特征长度（可以理解为膜或壳单元二维尺寸乘积的Word 资料.平方根），对于轴对称单元可以理解为r z 平面上的特征长度，对于粘聚单元可视为其组成厚度。 之所以通过这种方法定义裂缝特征长度，是因为裂缝产生方向无法预知。然而，裂缝的产生方向对高宽比（二维尺寸比）较大的单元的力学行为有很大不同的影响， 有些效应会因为网格敏感性而仍然存在，因此推荐最好尽可能地将单元划分为正方形网格。这种模拟混凝土脆性响应的方法需要指定线性逼近失效后应变软化时应力为 0 的位移（见图）失效应力 ,出现在失效应变（由破坏应力除以

11、弹性模量来定义）处；然而，应力在最终位移处降为0，这不依赖于试件长度。这意味着只用当试件足够短，以至于失效应变小于以下位移对应的应变值时，位移- 荷载试件才可以在破坏后保持静力平衡。其中， L 是试件的长度。InputFileUsage:Use both of the following options:*CONCRETE*TENSION STIFFENING,TYPE=DISPLACEMENTAbaqus/CAEUsage:Property module: material editor:MechanicalPlasticityConcrete Smeared Cracking:Subopt

12、ionsTension Stiffening: Type:DisplacementWord 资料.获取极限位移极限位移。可以用单位面积断裂能来评价，为，其中是混凝土能承受的最大拉伸应力。的典型值是从普通混凝土的 0.05mm（2×10 3英寸）到高强混凝土的（ × 3英寸）。的典型值大约为104 ，因此要求mm (20 in)。0.08mm 310临界长度如果试件长于临界长度L，当试件在固定的位移下开裂，相比于开裂过程消散的能量，更多的应变能将被储存在试件中。所以，一些应变能必须转换成动能，甚至在指定位移荷载下产生动态破坏，这就意味着， 在有限元中用此方法时， 单元的特征尺

13、寸必须小于临界长度，或者必须考虑额外（动态）的因素。分析输入文件处理器检查使用这一混凝土模型的每一个单元的特征长度，且不允许任何单元的特征长度超过。 用户必须用更小的单元重新划分网格或者用拉伸硬化的应力- 应变定义。由于断裂能方法通常仅用于素混凝土中，网格划分很少设置限制条件。开裂剪力传递当混凝土开裂时，它的剪切强度被削弱。这种现象通过减小剪切模量来定义，它是贯穿裂缝的开口应变的函数。用户还可以为闭合裂缝指定减小的剪切模量。当通过裂缝的正应力变为压力时，这一减小后的剪切模量也会产生作用。由于裂缝的存在，新的剪切刚度将会退化。裂缝的剪切模量定义为，其中G 是未开裂混凝土的弹性剪切模量，是乘法Wo

14、rd 资料.因子。剪力传递模型假定当裂缝宽度增加时开口裂缝的剪切刚度线性降低到0：其中是通过裂缝的直接应变，是一个用户指定值。这个模型同时假定开裂后闭合的裂缝剪切模量也降低了：这里用户指定。和 可以选择性的依赖温度或预定义场变量定义。如果混凝土弥散开裂模型的材料定义中不包含剪力传递，Abaqus/Standard将自动使用剪力传递的默认值，这样剪力响应就不会受到开裂的影响（完全剪力传递）。这一假定通常是合理的：在很多情况下，结构的整体响应并不是非常依赖于剪力传递量。InputFileUsage:Use both of the following options:*CONCRETE*SHEAR

15、RETENTIONAbaqus/CAEUsage:Property module: material editor:Word 资料.MechanicalPlasticityConcrete Smeared Cracking:SuboptionsShear Retention压缩行为当主应力分量中压应力占主导地位时，混凝土响应由弹塑性理论模拟，它采用依据等效压应力p 和 Mises 等效偏应力 q 写成的一种简单的屈服面形式， 其具体介绍见图采用了关联流动法则和各向强化法则。这种模型显著简化了真实行为。关联流动法则假定通常过度预测了非弹性体积应变。由于忽略了第三应力不变量， 屈服面不能和三轴受拉

16、和受压试验数据准确吻合。当混凝土的应变超过了极限应力点，假定弹性响应不受非弹性变形的影响是不现实的。另外，混凝土在承受非常高的压应力时表现出非弹性响应：而模型中没有试图考虑将这种行为。引入受压行为的相关简化是为了提高计算效率。尤其是当关联流动法则的假定没有试验数据检验时， 只要问题中压应力的范围不大，它可以保证结果是可接受的接近测量值。 从计算角度看，关联流动假定导致整体本构模型的雅可比矩阵（“材料刚度矩阵”）充分对称，这样整体平衡方程的解通常就不需要非对称方程的解。在牺牲计算成本的情况下，所有的此类限制都可以被去除。用户可以定义素混凝土在单轴压缩下弹性范围为外的应力- 应变行为。受压应力可以

17、写成成塑性应变的表格函数，如果需要，也可以写成成温度和场变量的函数。压应力和压应变必须为正值（绝对值）。应力应变曲线可以定义超过极限应力进入应变软化范围。Word 资料.InputFileUsage:*CONCRETEAbaqus/CAEUsage:Property module: material editor:MechanicalPlasticityConcrete Smeared Cracking单轴和多轴行为混凝土模型中的开裂和压缩响应同时用构件的单轴响应来详细阐述，如图所示。图素混凝土的单轴行为混凝土受压时最初表现出弹性响应。随着应力增大， 一些不可恢复的 （非弹性）应变产生，材料行

18、为发生软化。达到极限压应力后，材料强度持续降低直到不能承受任何压力。 如果在非弹性应变出现后的某一点卸载，卸载响应比初始弹Word 资料.性响应弱：说明弹性被破坏了。 模型中忽略了此效应， 因为我们假定模型主要应用与单调应变， 仅偶尔出现反向卸载。 当单轴混凝土试件试件受拉时，在达到极限压应力的7%-10%前的表现为弹性，随后裂缝产生。裂缝的形成非常快，即使使用最硬的测试仪器也很难观测到其真实行为。模型假定开裂引起损伤， 这意味着开口裂缝可以由弹性刚度的损失来描述。同时假定开裂过程不存在永久变形。如果裂缝应力变为受压时，裂缝可以完全闭合。在多轴应力状态下，这些理论通过应力空间上的破坏面和流动法

19、则的概念推而广之。破坏面应与实验数据吻合。 使用的破坏面如图、所示。图平面应力下的屈服和破坏面图（ p q）平面下的屈服和破坏面Word 资料.破坏面用户可以指定失效比来定义破坏面的形状（可能是温度和预定义常变量的函数）。可以指定四种失效比：1、极限双轴压应力与极限单轴压应力的比值。2、破坏时的单轴拉应力值与到单轴极限压应力的比值的绝对值。3、双轴极限压应力中塑性应变主成分的大小与单轴极限压应力下的塑性应变的比值。4、平面应力下其余主应力为极限压应力值时的开裂拉伸主应力与单轴受拉开裂应力的比值。如果用户不指定，以上比值将使用默认值。InputFileUsage:*FAILURE RATIOSA

20、baqus/CAEUsage:Property module: material editor:Word 资料.MechanicalPlasticityConcrete Smeared Cracking:SuboptionsFailure Ratios应变反向响应由于模型是用来解决相对单调应变问题的，它并没有试图给出循环响应或由主压应力下由非弹性应变引起的弹性刚度降低的问题。尽管如此，甚至有可能在已经设计好的模型应用中， 应变轨迹也不完全是放射状的，因此模型应该能够用合理的方式预示出偶然的应变反向和应变路径方向的改变。当主应力由压应力控制时，“受压”屈服面的各向同性硬化形成了模型非弹性反应预测

21、的基础。校正混凝土的最简化模型 （使用所有默认值并假定温度和场变量独立）至少需要单轴压缩和单轴拉伸两组实验来校正。为了得到准确的破坏后行为， 可能需要另外的实验。单轴压缩和拉伸实验单轴受压实验包括在两片刚性压板间压缩试件。加载方向上的力以及位移被记录下来。从这个记录中，可以直接提取出所需要的混凝土模型应力- 应变曲线。单轴拉伸实验较难实现， 这意味着必须有刚性测试仪器能够记录下来破坏后的响应。这一实验常常实现不了， 可以假定混凝土拉伸破坏强度（通常大约是抗压强度的 7%10%）。拉伸开裂应力的选择是十分重要的，如果使用很小的开裂应力（小于抗压强度的1/100 或 1/1000 ）将会导致数值问题。Word 资料.开裂后拉伸行为后破坏效应的校核依赖于添加在混凝土中的钢筋率。对于素混凝土的模拟，应该用应力 - 位移拉伸硬化模型。的一般值为对普通混凝土0.05 mm (2 ×103in)到高强混凝土 0.08 mm (3 × 10 3 in)。对于钢筋混凝土的模拟，应该用应力 - 应变拉伸硬化模型。对于在较细网格下的大量混凝土材料，当总应变10 倍于破坏应变时，对于相当重的网格过细的钢筋混凝土模型，合理的起点被假定为