ANSYS中混凝土的计算问题
最近做了点计算分析，结合各论坛关于这方面的讨论，就一些问题探讨如下，不当之处敬请指正。
一、关于模型
钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种，即分离式、分布式和组合式模型。考虑钢筋和混凝土之间的粘结和滑移，则采用引入粘结单元的分离式模型；假定混凝土和钢筋粘结很好，不考虑二者之间的滑移，则三种模型都可以；分离式和分布式模型适用于二维和三维结构分析，后者对杆系结构分析比较适用。裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型，后者目前尚处研究之中，主要应用的是前两种。离散裂缝模型和分布裂缝模型各有特点，可根据不同的分析目的选择使用。随着计算速度和网格自动划分的快速实现，离散裂缝模型又有被推广使用的趋势。
就ANSYS而言，她可以考虑分离式模型(solid65+link8，认为混凝土和钢筋粘结很好，如要考虑粘结和滑移，则可引入弹簧单元进行模拟，比较困难！)，也可采用分布式模型(带筋的solid65)。而其裂缝的处理方式则为分布裂缝模型。
二、关于本构关系
混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类，其中研究最多的是非线性弹性和弹塑性本构关系，其中不乏实用者。混凝土破坏准则从单参数到五参数模型达数十个模型，或借用古典强度理论或基于试验结果等，各个破坏准则的表达方式和繁简程度各异，适用范围和计算精度差别也比较大，给使用带来了一定的困难。
就ANSYS而言，其问题比较复杂些。
1 ANSYS混凝土的破坏准则与屈服准则是如何定义的？
采用tb,concr,matnum则定义了W-W破坏准则(failure criterion)，而非屈服准则(yi eld criterion)。W-W破坏准则是用于检查混凝土开裂和压碎用的，而混凝土的塑性可以另外考虑（当然是在开裂和压碎之前）。理论上破坏准则(failure criterion)和屈服准则(y ield criterion)是不同的，例如在高静水压力下会发生相当的塑性变形，表现为屈服，但没有破坏。而工程上又常将二者等同，其原因是工程结构不容许有很大的塑性变形，且混凝土等材料的屈服点不够明确，但破坏点非常明确。
定义tb,concr matnum后仅仅是定义了混凝土的破坏准则和缺省的本构关系，即W—W破坏准则、混凝土开裂和压碎前均为线性的应力应变关系，而开裂和压碎后采用其给出的本构关系。但屈服准则尚可另外定义(随材料的应力应变关系，如tb,MKIN,则定义的屈服准则是Von Mises,流动法则、硬化法则也就确定了)。
2 定义tb,concr后可否定义其它的应力应变关系
当然是可以的，并且只有在定义tb,concr后，有些问题才好解决。例如可以定义tb,mi so,输入混凝土的应力应变关系曲线(多折线实现)，这样也就将屈服准则、流动法则、硬化法则等确定了。
这里可能存在一点疑问，即ANSYS中的应力应变关系是拉压相等的，而混凝土材料显然不是这样的。是的，因为混凝土受拉段非常短，认为拉压相同影响很小，且由于定义的tb, concr中确定了开裂强度，所以尽管定义的是一条大曲线，但应用于受拉部分的很小。
三、具体的系数及公式
1 定义tb,concr时候的两个系数如何确定？
一般的参考书中，其值建议先取为0.3~0.5(江见鲸)，原话是“在没有更仔细的数据时，不妨先取0.3~0.5进行计算”，足见此0.3~0.5值的可用程度。根据我的经验和理由，建议此值取大些，即开裂的剪力传递系数取0.5,（定要>0.2）闭合的剪力传递系数取1.0。支持此说法的还有现行铁路桥规的抗剪计算理论，以及原公路桥规的容许应力法的抗计剪计算。
2 定义混凝土的应力应变曲线
单向应力应变曲线很多，常用的可参考国标混凝土结构规范，其中给出的应力应变曲线是二次曲线＋直线的下降段，其参数的设置按规范确定即可。当然如有实测的应力应变曲线更好了。
四、关于收敛的问题
ANSYS混凝土计算收敛(数值)是比较困难的，主要影响因素是网格密度、子步数、收敛准则等，这里讨论如下。
1 网格密度：网格密度适当能够收敛。不是网格越密越好，当然太稀也不行，这仅仅是就收敛而言的，不考虑计算费用问题。但是究竟多少合适，没有找到规律，只能靠自己针对情况慢慢试算。
2 子步数：NSUBST的设置很重要，设置太大或太小都不能达到正常收敛。这点可以从收敛过程图看出，如果F范数曲线在[F]曲线上面走形的很长，可考虑增大nsubst。或者根据经验慢慢调正试算。
3 收敛精度：实际上收敛精度的调正并不能彻底解决收敛的问题，但可以放宽收敛条件以加速吧。一般不超过5%(缺省是0.5%)，且使用力收敛条件即可。
4 混凝土压碎的设置：不考虑压碎时，计算相对容易收敛；而考虑压碎则比较难收敛，即便是没有达到压碎应力时。如果是正常使用情况下的计算，建议关掉压碎选项；如果是极限计算，建议使用concr＋MISO且关闭压碎检查；如果必设压碎检查，则要通过大量的试算（设置不同的网格密度、NSUBST）以达到目的，但也很困难。
5 其他选项：如线性搜索、预测等项也可以打开，以加速收敛，但不能根本解决问题。
6 计算结果：仅设置concr，不管是否设置压碎，其一般P-F曲线接近二折线；采用conc r+miso则P-F曲线与二折线有差别，其曲线形状明显是曲线的。
*********************************************************** ********************
例题1
!----------------------------------------------------
!题目：钢筋混凝土简支梁模拟计算
!方法：分离式；solid65和link8
!材料：混凝土采用concr和钢筋为弹性材料，但不考虑压碎
!---------------------------------------------------
!为方便，假定钢筋置于梁底两侧.
!=============================================== ====
/config,nres,2000 !子步数最大为2000
/prep7
!定义单元及其材料特性等
rd0=20.0 !钢筋直径
et,1,solid65
et,2,link8
mp,ex,1,33e3
mp,prxy,1,0.20
r,1
hntra=28
hntrl=2.6
tb,concr,1 !混凝土材料1
tbdata,,0.7,1.0,hntrl,-1
mp,ex,2,2.1e5
mp,prxy,2,0.30
r,2,acos(-1)*0.25*rd0*rd0
!定义梁体即单元划分
blc4, , ,100,200,3000
/view,1,1,1,1
/ang,1
gplot
!定义网分时边长控制
lsel,s,loc,z,1,2999
lsel,r,loc,y,0
latt,2,2,2
lesize,all,,,20 !钢筋网格数目
lmesh,all
lsel,s,loc,z,0
lesize,all,,,4 !截面上的网格数目4x4
vsel,all
vatt,1,1,1
mshape,0,3d
mshkey,1
vmesh,all
allsel,all
finish
/solu
!施加约束
lsel,s,loc,z,0 lsel,r,loc,y,0
dl,all,,uy
dl,all,,uz
lsel,all
lsel,s,loc,z,3000 lsel,r,loc,y,0
dl,all,,uy
lsel,all
ksel,s,loc,x,0 ksel,r,loc,y,0
dk,all,ux
allsel,all
!施加荷载
qmz=0.3
asel,s,loc,y,200 sfa,all,1,pres,qmz allsel,all
nsubst,40 outres,all,all time,qmz*10 neqit,40
solve
finish
/post1
pldisp,1
etable,zxyl,ls,1 plls,zxyl,zxyl,1 finish
/post26
nsol,2,33,u,y prod,3,1,,,,,,1/100 prod,4,2,,,,,,-1 xvar,4
plvar,3
*********************************************************** ********************
例题2
!----------------------------------------------------
!题目：钢筋混凝土简支梁模拟计算
!方法：分离式；solid65和link8
!材料：混凝土采用concr+Miso和钢筋为弹性材料，但不考虑压碎
!增加网格密度
!---------------------------------------------------
!为方便，假定钢筋置于梁底两侧.
!=============================================== ====
/config,nres,2000
/prep7
!定义单元及其材料特性等
rd0=20.0 !钢筋直径
et,1,solid65
et,2,link8
mp,ex,1,33e3
mp,prxy,1,0.20
r,1
hntra=28
hntrl=2.6
tb,concr,1
tbdata,,0.7,1.0,hntrl,-1
tb,miso,1,,14
tbpt,,0.0002,hntra*0.19
tbpt,,0.0004,hntra*0.36
tbpt,,0.0006,hntra*0.51
tbpt,,0.0008,hntra*0.64
tbpt,,0.0010,hntra*0.75
tbpt,,0.0012,hntra*0.84
tbpt,,0.0014,hntra*0.91
tbpt,,0.0016,hntra*0.96
tbpt,,0.0018,hntra*0.99
tbpt,,0.0020,hntra*1.00
tbpt,,0.0025,hntra*0.95
tbpt,,0.0030,hntra*0.90
tbpt,,0.0035,hntra*0.85 tbpt,,0.0040,hntra*0.80 mp,ex,2,2.1e5
mp,prxy,2,0.30
r,2,acos(-1)*0.25*rd0*rd0 !定义梁体即单元划分
blc4, , ,100,200,3000
/view,1,1,1,1
/ang,1
gplot
!定义网分时边长控制
lsel,s,loc,z,1,2999
lsel,r,loc,y,0
latt,2,2,2
lesize,all,,,20
lmesh,all
lsel,s,loc,z,0
lesize,all,,,4
vsel,all
vatt,1,1,1
mshape,0,3d
mshkey,1
vmesh,all
allsel,all
finish
/solu
!施加约束
lsel,s,loc,z,0
lsel,r,loc,y,0
dl,all,,uy
dl,all,,uz
lsel,all
lsel,s,loc,z,3000
lsel,r,loc,y,0
dl,all,,uy
lsel,all
ksel,s,loc,x,0
ksel,r,loc,y,0
dk,all,ux
allsel,all
!施加荷载
qmz=0.3
asel,s,loc,y,200 sfa,all,1,pres,qmz allsel,all
outres,all,all time,qmz*10 nsubst,40
neqit,40
solve
finish
/post1
pldisp,1
etable,zxyl,ls,1 plls,zxyl,zxyl,1 finish
/post26
nsol,2,33,u,y prod,3,1,,,,,,1/100 prod,4,2,,,,,,-1 xvar,4
plvar,3
1. 讨论两种Ansys求极限荷载的方法 （1）力加载 可以通过对应的方法（比如说特征值屈曲）估计结构的极限荷载的大致范围，然后给结构施加一个稍大的荷载，打开自动荷载步二分法进行非线性静力分析，最后计算会因不收敛终止，则倒数第二个子步对应的就是结构的极限荷载；另外，也可以选择弧长法，采用足够的子步（弧长法可以一直分析到极限承载力之后的过程）同样可以从绘制的荷载位移曲线或计算结果中找出结构的极限荷载。 （2）位移加载 给结构施加一个比较大的位移，打开自动荷载步二分法进行非线性分析，保证足够的子步数，这样也可以分析到极限荷载以后，通过绘制荷载位移曲线或查看相应结果文件也可知道结构的极限荷载。 希望众高手讨论一下 （1）弧长法求极限荷载的收敛性问题，如何画到荷载位移曲线的下降段？ （2）位移法求极限荷载的具体步骤？ 2. 需要注意的问题 1. 由于SOLID 65单元本身是基于弥散裂缝模型和最大拉应力开裂判据，因此在很多情况下会因为应力集中而使混凝土提前破坏，从而和试验结果不相吻合，因此，在实际应用过程中应该对单元分划进行有效控制，根据作者经验，当最小单元尺寸大于5cm 时，就可以有效避免应力集中带来的问题； 2. 支座是另一个需要注意的问题。在有限元分析中，很多时候约束是直接加在混凝土节点上，这样很可能在支座位置产生很大的应力集中，从而使支座附近的混凝土突然破坏，造成求解失败。因此，在实际应用过程中，应该适当加大支座附近单元的尺寸或者在支座上加一些弹性垫块，避免支座的应力集中； 3. 六面体的SOLID 65 单元一般比四面体的单元计算要稳定且收敛性好，因此，只要条件允许，应该尽量使用六面体单元； 4. 正确选择收敛标准，一般位移控制加载最好用位移的无穷范数控制收敛，而用力控制加载时可以用残余力的二范数控制收敛。在裂缝刚刚出现和接近破坏的阶段，可以适当放松收敛标准，保证计算的连续性； 3. 关于下降段的问题 1）在实际混凝土中都有下降段，但是在计算的时候要特别小心下降段的问题。 2）下降段很容易导致计算不收敛，有时为了计算的收敛要避免设置下降段，采用rush模型。 3）利用最大压应变准则来判断混凝土是否破坏。 4. Solid65单元中的破坏准则 1）采用Willam&Warnke五参数破坏准则 2）需要参数： 单轴抗拉强度，单轴，双轴抗压强度，围压压力，在围压作用下双轴，单轴抗压强度
ABAQUS用连续介质的方法建立描述混凝土模型不采用宏观离散裂纹的方法描述裂纹的水平的在每一个积分点上单独计算其中。 低压力混凝土的本构关系包括： Con crete Smeared cracki ng model (ABAQUS/Sta ndard) Concrete Brittle cracki ng model (ABAQUS/Explicit) Con crete Damage plasticity model 高压力混凝土的本构关系： Cap model 1、ABAQUS/Standard 中的弥散裂缝模型Concrete Smeared cracking model (ABAQUS/Standard): 只能用于ABAQUS/Standard 中 裂纹是影响材料行为的最关键因素，它将导致开裂以及开裂后的材料的各向异性 用于描述：单调应变、在材料中表现出拉伸裂纹或者压缩时破碎的行为 在进行参数定义式的Keywords: *CONCRETE *TENSION STIFFENING *SHEAR RETENTION *FAILURE RATIOS 2、ABAQUS/Explicit 中脆性破裂模型Concrete Brittle cracking model (ABAQUS/Explicit): 适用于拉伸裂纹控制材料行为的应用或压缩失效不重要，此模型考虑了由于裂纹引起的材料 各向异性性质，材料压缩的行为假定为线弹性，脆性断裂准则可以使得材料在拉伸应力过大 时失效。 在进行参数定义式的Keywords *BRITTLE CRACKING, *BRITTLE FAILURE, *BRITTLE SHEAR 3、塑性损伤模型Concrete Damage plasticity model : 适用于混凝土的各种荷载分析，单调应变，循环荷载，动力载荷，包含拉伸开裂(cracking)和压缩破碎(crushing)，此模型可以模拟硬度退化机制以及反向加载刚度恢复的混凝土力学特性 在进行参数定义式的Keywords: *CONCRETE DAMAGED PLASTICITY *CONCRETE TENSION STIFFENING *CONCRETE COMPRESSION HARDENING *CONCRETE TENSION DAMAGE
一、关于模型 钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种，即分离式、分布式和组合式模型。考虑钢筋和混凝土之间的粘结和滑移，则采用引入粘结单元的分离式模型；假定混凝土和钢筋粘结很好，不考虑二者之间的滑移，则三种模型都可以；分离式和分布式模型适用于二维和三维结构分析，后者对杆系结构分析比较适用。裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型，后者目前尚处研究之中，主要应用的是前两种。离散裂缝模型和分布裂缝模型各有特点，可根据不同的分析目的选择使用。随着计算速度和网格自动划分的快速实现，离散裂缝模型又有被推广使用的趋势。 就ANSYS而言，她可以考虑分离式模型(solid65+link8，认为混凝土和钢筋粘结很好，如要考虑粘结和滑移，则可引入弹簧单元进行模拟，比较困难！)，也可采用分布式模型(带筋的solid65)。而其裂缝的处理方式则为分布裂缝模型。 二、关于本构关系 混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类，其中研究最多的是非线性弹性和弹塑性本构关系，其中不乏实用者。混凝土破坏准则从单参数到五参数模型达数十个模型，或借用古典强度理论或基于试验结果等，各个破坏准则的表达方式和繁简程度各异，适用范围和计算精度差别也比较大，给使用带来了一定的困难。 就ANSYS而言，其问题比较复杂些。 1 ANSYS混凝土的破坏准则与屈服准则是如何定义的？ 采用tb,concr,matnum则定义了W-W破坏准则(failure criterion)，而非屈服准则(yield criterion)。W-W破坏准则是用于检查混凝土开裂和压碎用的，而混凝土的塑性可以另外考虑（当然是在开裂和压碎之前）。理论上破坏准则(failure criterion)和屈服准则(yield criterion)是不同的，例如在高静水压力下会发生相当的塑性变形，表现为屈服，但没有破坏。而工程上又常将二者等同，其原因是工程结构不容许有很大的塑性变形，且混凝土等材料的屈服点不够明确，但破坏点非常明确。 定义tb,concr matnum后仅仅是定义了混凝土的破坏准则和缺省的本构关系，即W—W破坏准则、混凝土开裂和压碎前均为线性的应力应变关系，而开裂和压碎后采用其给出的本构关系。但屈服准则尚可另外定义(随材料的应力应变关系，如tb,MKIN,则定义的屈服准则是Von Mises,流动法则、硬化法则也就确定了)。 2 定义tb,concr后可否定义其它的应力应变关系 当然是可以的，并且只有在定义tb,concr后，有些问题才好解决。例如可以定义tb,miso,输入混凝土的应力应变关系曲线(多折线实现)，这样也就将屈服准则、流动法则、硬化法则等确定了。 这里可能存在一点疑问，即ANSYS中的应力应变关系是拉压相等的，而混凝土材料显然不是这样的。是的，因为混凝土受拉段非常短，认为拉压相同影响很小，且由于定义的tb,concr 中确定了开裂强度，所以尽管定义的是一条大曲线，但应用于受拉部分的很小。 三、具体的系数及公式 1 定义tb,concr时候的两个系数如何确定？ 一般的参考书中，其值建议先取为0.3~0.5(江见鲸)，原话是“在没有更仔细的数据时，不妨先取0.3~0.5进行计算”，足见此0.3~0.5值的可用程度。根据我的经验和理由，建议此值取大些，即开裂的剪力传递系数取0.5,（定要>0.2）闭合的剪力传递系数取1.0。支持此说法的还有 现行铁路桥规的抗剪计算理论，以及原公路桥规的容许应力法的抗计剪计算。
第４３卷第１期２０１５年１月同济大学学报（自然科学版） ＪＯＵＲＮＡＬ ＯＦ ＴＯＮＧＪ Ｉ ＵＮＩＶＥＲＳＩＴＹ（ＮＡＴＵＲＡＬ ＳＣＩＥＮＣＥ）Ｖｏｌ．４３ Ｎｏ．１ Ｊ ａｎ．２０１５文章编号：０２５３－３７４Ｘ（２０１５）０１－０００１－０７ ＤＯＩ：１０．１１９０８／ｊ ．ｉｓｓｎ．０２５３－３７４ｘ．２０１５．０１．００１收稿日期：２０１４－０３－ ０４基金项目：国家自然科学基金重大国际合作项目（５１２６１１２０３７４）；国家自然科学基金集成项目（９１３１５３０１ ）第一作者：冯德成（１９８７—），男，博士生，主要研究方向为结构非线性分析．Ｅ－ｍａｉｌ：ａｕｆｄｃ＠１６３．ｃｏｍ通讯作者：李 杰（１９５７—） ，男，教授，博士生导师，工学博士，主要研究方向为混凝土随机损伤力学、随机动力系统分析与生命线工程抗灾．Ｅ－ｍａｉｌ：ｌｉｊｉｅ＠ｔｏｎｇｊ ｉ．ｅｄｕ．ｃｎ非均匀受压下的箍筋约束混凝土本构模型 冯德成１，万增勇１，李 杰１，２ （１．同济大学土木工程学院，上海２０００９２；２．同济大学土木工程防灾国家重点实验室，上海２０００９２ ）摘要：以Ｍａｎｄｅｒ提出的箍筋约束混凝土模型为基础，考虑构件非均匀受压下截面应变梯度对箍筋约束效应的影响，引入偏心率系数反映非均匀受压下偏心率对箍筋有效约束力的影响，建立了一类新的箍筋约束混凝土模型．将这一模型与柔度法纤维梁柱单元相结合，实现了在计算过程中动态更新构件不同位置、 不同受力状态下的截面偏心率以及相应的约束混凝土应力－应变关系． 对钢筋混凝土柱的分析结果表明建立的模型物理意义明确、计算精度较高． 关键词：箍筋约束混凝土；非均匀受压；偏心率；柔度法梁柱单元 中图分类号：ＴＵ５２８．１ 文献标志码：Ａ Ｈｏｏｐ Ｒｅｉｎｆｏｒｃｅｍｅｎｔ Ｃｏｎｆｉｎｅｄ ＣｏｎｃｒｅｔｅＣｏｎｓｔｉｔｕｔｉｖｅ Ｍｏｄｅｌ ｆｏｒ Ｎｏｎ－ｕｎｉｆｏｒｍｌｙＣｏｍｐｒｅｓｓｉｏｎ ＦＥＮＧ Ｄｅｃｈｅｎｇ１ ，ＷＡＮ Ｚｅｎｇｙｏｎｇ１ ，ＬＩ Ｊ ｉｅ１，２ （１．Ｃｏｌｌｅｇｅ ｏｆ Ｃｉｖｉｌ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｔｏｎｇｊｉ Ｕｎｉｖｅｒｓｉｔｙ，Ｓｈａｎｇｈａｉ２０００９２，Ｃｈｉｎａ；２．Ｓｔａｔｅ Ｋｅｙ Ｌａｂｏｒａｔｏｒｙ ｏｆ Ｄｉｓａｓｔｅｒ Ｒｅｄｕｃｔｉｏｎ ｉｎ ＣｉｖｉｌＥｎｇｉｎｅｅｒｉｎｇ ，Ｔｏｎｇｊｉ Ｕｎｉｖｅｒｓｉｔｙ，Ｓｈａｎｇｈａｉ ２０００９２，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｂａｓｅｄ ｏｎ ｔｈｅ ｃｏｎｆｉｎｅｄ ｃｏｎｃｒｅｔｅ ｍｏｄｅｌ ｐｒｏｐｏｓｅｄ ｂｙＭａｎｄｅｒ，ｔｏ ｃｏｎｓｉｄｅｒ ｔｈｅ ｅｆｆｅｃｔ ｏｆ ｔｈｅ ｓｅｃｔｉｏｎａｌ ｓｔｒａｉｎ ｇｒａｄｉｅｎｔｏｎ ｔｈｅ ｃｏｎｆｉｎｅｍｅｎｔ ｅｆｆｅｃｔ ｕｎｄｅｒ ｎｏｎ－ｕｎｉｆｏｒｍｌｙ ｃｏｍｐｒｅｓｓｉｏｎ，ａｎｅｗ ｃｏｎｆｉｎｅｄ ｃｏｎｃｒｅｔｅ ｍｏｄｅｌ ｉｓ ｄｅｖｅｌｏｐｅｄ ｉｎ ｔｈｉｓ ｐａｐｅｒ ｂｙｉｎｔｒｏｄｕｃｉｎｇ ｔｈｅ ｅｃｃｅｎｔｒｉｃｉｔｙ ｒａｔｉｏ ｆａｃｔｏｒ ｔｏ ｒｅｆｌｅｃｔ ｔｈｅ ｉｎｆｌｕｅｎｃｅｏｆ ｔｈｅ ｅｃｃｅｎｔｒｉｃｉｔｙ ｏｎ ｔｈｅ ｃｏｎｆｉｎｉｎｇ ｆｏｒｃｅ．Ｍｅａｎｗｈｉｌｅ，ｂｙｃｏｍｂｉｎｉｎｇ ｔｈｅ ｍｏｄｅｌ ｗｉｔｈ ｔｈｅ ｆｉｂｅｒ ｆｏｒｃｅ－ｂａｓｅｄ ｂｅａｍ－ｃｏｌｕｍｎｅｌｅｍｅｎｔ，ｉｔ ｃａｎ ａｄｊｕｓｔ ｔｈｅ ｅｃｃｅｎｔｒｉｃｉｔｙ ｒａｔｉｏ ａｎｄ ｔｈｅｃｏｒｒｅｓｐｏｎｄｉｎｇ ｓｔｒｅｓｓ－ｓｔｒａｉｎ ｒｅｌａｔｉｏｎｓｈｉｐ ｏｆ ｔｈｅ ｓｅｃｔｉｏｎ ａｔｄｉｆｆｅｒｅｎｔ ｌｏｃａｔｉｏｎｓ ａｎｄ ｄｉｆｆｅｒｅｎｔ ｌｏａｄｉｎｇ ｓｔａｔｅｓ ｄｕｒｉｎｇｃａｌｃｕｌａｔｉｏｎ．Ｔｈｅ ａｎａｌｙｓｉｓ ｏｆ ｒｅｉｎｆｏｒｃｅｄ ｃｏｎｃｒｅｔｅ ｃｏｌｕｍｎｓｉｌｌｕｓｔｒａｔｅｓ ｔｈａｔ ｔｈｅ ｍｏｄｅｌ ｈａｓ ａ ｃｌｅａｒ ｐｈｙｓｉｃａｌ ｍｅａｎｉｎｇ ａｎｄｓｈｏｗｎ ｔｏ ｂｅ ｅｆｆｅｃｔｉｖｅ．Ｋｅｙ ｗｏｒｄｓ：ｈｏｏｐ ｒｅｉｎｆｏｒｃｅｍｅｎｔ ｃｏｎｆｉｎｅｄ ｃｏｎｃｒｅｔｅ；ｎｏｎ－ｕｎｉｆｏｒｍｌｙ ｃｏｍｐｒｅｓｓｉｏｎ；ｅｃｃｅｎｔｒｉｃｉｔｙ ｒａｔｉｏ；ｆｏｒｃｅ－ｂａｓｅｄ ｂｅａｍ－ｃｏｌｕｍｎ ｅｌｅｍｅｎｔ 有关约束混凝土的研究已有近百年的历史． 一般认为，这一历史最早可以追溯到１９０３年 Ｃｏｎｓｉｄｅｒｅ［１］发现利用螺旋箍筋能有效提高轴心受压 柱的承载力．１９２８年，Ｒｉｃｈａｒｔ［２］ 首次定量地研究了 液体围压对混凝土圆柱体轴压性能的影响， 并提出了相应的约束混凝土抗压强度以及峰值应变的计算 公式；１９５５年，Ｃｈａｎ［３］ 在试验的基础上提出了箍筋 约束混凝土的应力－应变关系模型，并认为，箍筋的 约束作用仅仅体现在对峰值应变的提高方面，而对强度影响甚微．此后的发展，多沿着试验研究—理论解释的基本路线，试图根据试验结果提出相应的约束混凝土的应力－应变关系模型．１９７１年，Ｋｅｎｔ和 Ｐａｒｋ［４］ 总结了前人的研究结果， 提出了一个上升段为二次抛物线、下降段为直线且斜率由体积配箍率、 混凝土强度和箍筋间距等因素决定的应力－应变关系模型． 该模型是这一时期的集大成之作，应用最为广泛，其表达形式也多为后来的研究者所采纳． ２０世纪７ ０年代之前的研究也具明显的时代局限性．由于当时的结构设计思想主要停留在承载能力设计阶段， 因此，对于材料本构关系下降段的关注不多；并且，由于试验设备的限制，难以准确测定混凝土应力－应变曲线的下降段．这些因素使得基于试 验提出的本构关系模型的下降段十分粗糙［５］ ． 尽管如此， 这一时期对于箍筋约束效应的认识以及其基本影响因素的辨识仍然为后来的研究提供了框架和基础． １９８２年，Ｓｃｏｔｔ等［６］ 在Ｋｅｎｔ－ Ｐａｒｋ模型的基础上考虑了应变率的影响；同年，Ｓｈｅｉｋｈ和Ｕｚｕｍｅｒｉ［７］ 发 现了矩形截面中的约束“ 拱效应”，并提出了有效约
摘要 本文介绍ANSYS 模拟钢筋混凝土梁的过程，讨论了有限元模型的建立以及在 ANSYS 中的实现，给出了用分离式配筋方法对混凝土梁的分析的一般过程。并给出了详细的命令流过程。并在此基础上对混凝土梁进行了分析，讨论了在力的作用下混凝土梁的塑形变形和裂缝的发展过程。 关键词 Ansys 混凝土梁 分离式配筋 The analysis of mechanics of a reinforced concrete based on ANSYS Abstract This paper introduces ANSYS simulation of the reinforced concrete beam process, discusses the establishment of the finite element model and the realization, and gives the ANSYS reinforcement method with separate the analysis of concrete beams of the general process. And gives the detailed command flow process. Based on the analysis of concrete beams, and discussed the concrete beam under the action of forces of the body deformation and fracture process. Keywords Ansys concrete beams reinforced separated 1 引言 由于钢筋混凝上材料性质复杂，使其表现出明显的非线性行为[1]。长期以来采用线弹性理论的设计方法来研究钢筋混凝上结构的应力或内力，显然不太合理，尽管有此理论是基于人量试验数据上的经验公式，还是不能准确反映混凝上的力学性能，特别是受力复杂的重要结构，必须采用三维钢筋混凝上非线性有限元方法才能很好地掌握其力学性能。利用ANSYS 对钢筋混凝上结构弹塑性的仿真分析，可以对结构自开始受荷载直到破坏的全过程进行分析，获得不同阶段的受力性能。本文将以混凝土梁的弹塑性分析为例，介绍在Ansys 中分析材料非线性问题的具体实现方法。 2 问题介绍 如图所示的钢筋混凝土梁[2]，横截面尺寸为200400b h mm mm ?=?,梁的跨度为3.0L m =，支座宽度为250mm 采用C20混凝土，梁内受拉纵筋3φ20，架立筋采用2φ12， 箍筋采用φ6@150，钢筋保护层厚度为25mm 。如图一。 图一 对于梁中所采用的所有钢筋，弹性模量为5 2.110MPa ?，抗拉强度设计值210MPa ， 密度33 7.810/kg m ?，泊松比为0.3。
一、混凝土本构关系模型 1.混凝土单轴受压应力-应变关系 （1）Saenz 等人的表达式 Saenz 等人（1964年）所提出的应力-应变关系为： ])()()( /[30 200εεεεεεεσd c b a E +++= （2）Hognestad 的表达式 Hognestad 建议模型，其上升段为二次抛物线，下降段为斜直线。所提出的应力-应变关系为： cu cu εεεσσεεσσεεεεεεεε≤≤-=≤-=--000 02,)]( 15.01[,])(2[0 00 （3）我国《混凝土结构设计规范》（GB50010-2010）中的混凝土受压应力-应变曲线，其表达式为： 1,)1(1 ,)1(2>+-=≤+-= x x x x y x x n nx y c n α r c x ,εε= ，r c f y ,σ= ，r c r c c r c c f E E n ,,,-=εε c α是混凝土单轴受压时的应力应变曲线在下降段的参数值，r c f ,是混凝土单轴抗压的 强度代表值，r c ,ε是与单轴抗压强度r c f ,相对应的混凝土峰值压应变。 2.混凝土单轴受拉应力-应变关系 清华大学过镇海等根据实验结果得出混凝土轴心受拉应力-应变曲线： 1 ],)1(/[)/(1 ,])(2.0)(2.1[7 .16≥+-?=≤-=t t t t t t t t t t εε εεεεεεεεεεασεεσσσ 3.混凝土线弹性应力-应变关系 张量表达式，对于未开裂混凝土，其线弹性应力应变关系可用不同材料常数表达，其中用材料弹性模量E 和泊松比v 表达的应力应变关系为： ij kk E ij E ij ij kk E ij E ij δσσεδεεσν ν νννν-=+=+-++1)21)(1(1
1 横向配筋的作用 混凝土结构中的配筋有两种：直接钢筋和间接钢筋。直接配筋即沿构件轴力或主应力方向设置的纵向钢筋，直接承担拉力或者压力，钢筋的应力与轴力方向一致；间接配筋又称横向配筋，沿与压应力与最大主压应力垂直的方向设置，通过约束混凝土的横向变形，提高轴向抗压承载力。 横向配筋有多种，比如螺旋（圆形）箍筋、矩形箍筋、钢管、焊接网片等。其主要作用是约束其内部混凝土的横向变形，使之处于三轴受压应力状态，从而提高了其强度和变形能力。 下面就箍筋对混凝土的约束作用做以简单分析。 箍筋的作用有许多种， ?抗剪。除了直接承受剪力外，还间接限制了斜裂缝的开展宽度，增强了腹部混凝土的骨料咬合力；还约束了纵筋对混凝土保护层的撕脱，增大了 钢筋的销栓力；同时，纵筋与腹筋形成的骨架使内部混凝土受到约束， 这也有利于抗剪； ?通过减小纵筋的自由长度，防止纵筋受力后压屈，充分发挥其抗压强度，同时也起到固定纵筋位置的作用； ?对于密排箍筋，通过约束核心区混凝土，提高了混凝土的抗压强度及延性（极限变形能力）； ?长期荷载作用下，可以承受因混凝土收缩和环境湿度变化等产生的横向应力，以防止或减少纵向裂缝； 其中，通过约束核心区混凝土，提高受压混凝土的抗压强度及延性，对于地震区的混凝土结构尤为重要。适当地增加箍筋和改进构造形式成为提高结构抗震性能的最简单、经济和有效的措施之一。 2 影响箍筋约束作用的因素 箍筋对约束混凝土的增强作用，除了受被约束混凝土自身强度的影响外，主要取决于它能够施加在核心区混凝土表面的约束力的大小。约束力越大，对混凝土的增强就越多。约束力主要受以下几个因素影响： ?体积配箍率。体积配箍率隐含反应了四个因素：箍筋强度、直径、间距及（计算配箍方向的）核心区宽度（对于螺旋或圆形配箍的圆形截面，指 核心区直径）。箍筋的强度和直径直接决定了箍筋所能提供的约束力的 大小，箍筋间距及核心区宽度则影响约束力在相邻箍筋间的分布。对于 矩形截面，通常两个方向上的尺寸和配箍形式不一样，因此提供的约束 力也不一样，所以应分别计算两个方向的配箍率。
ANSYS 钢筋混凝土建模 一、简介 钢筋混凝土有限元建模的方法与结果评价（前后处理），是对钢筋混凝土结构进行数值模拟的重要步骤，能否把握模型的可行性、合理性，如何从计算结果中寻找规律，是有限元理论应用于实际工程的关键一环。Blackeage以自己做过的一组钢筋混凝土暗支撑剪力墙的数值模拟为例，从若干方面提出一些经验与建议。希望大家一起讨论、批评指正(wang.jian@https://www.doczj.com/doc/c06030398.html,)。 程序：ANSYS 单元：SOLID65、BEAM188 建模方式：分离 暗支撑剪力墙结构由北京工业大学曹万林所提出，简言之就是一种在普通钢筋配筋情况下，加配斜向钢筋的剪力墙结构。 二、单元选择 以前经常采用的钢筋混凝土建模方法是通过SOLID65模拟混凝土，通过SOLID65的实常数指定钢筋配筋率，后来发现这种整体式的模型并不理想，而且将钢筋周围的SOLID65单元选择出来，再换算一个等效的配筋率，工作量也并不小。最关键的是采用整体式模型之后，得不出什么有意义的结论，弄一个荷载-位移曲线出来又和实验值差距比较大。只有计算的开裂荷载与实验还算是比较接近，但这个手算也算得出来的东西费劲去装模作样的建个模型又有什么意义？ 所以，这次我尝试采用分离式的模型，钢筋与混凝土单元分别建模，采用节点共享的方式。建模时发现，只要充分、灵活地运用APDL的技巧，处理好钢筋与混凝土单元节点的位置，效率还是很高的。 暗支撑剪力墙数值模型 看过很多的资料，分离式模型是用LINK8与SOLID65的组合方式，这样做到是非常直观，因为LINK8是spar类型的单元，每个节点有3个自由度，这与SOLID65单元单节点自由度数量是一致的。但是问题也就由此产生，当周围的混凝土开裂或是压碎时，SOLID65将不能对LINK8的节点提供足够地约束（如
ANSYS中混凝土的计算问题【精华】 最近做了点计算分析，结合各论坛关于这方面的讨论，就一些问题探讨如下，不当之处敬请指正。 一、关于模型 钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种，即分离式、分布式和组合式模型。考虑钢筋和混凝土之间的粘结和滑移，则采用引入粘结单元的分离式模型；假定混凝土和钢筋粘结很好，不考虑二者之间的滑移，则三种模型都可以；分离式和分布式模型适用于二维和三维结构分析，后者对杆系结构分析比较适用。裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型，后者目前尚处研究之中，主要应用的是前两种。离散裂缝模型和分布裂缝模型各有特点，可根据不同的分析目的选择使用。随着计算速度和网格自动划分的快速实现，离散裂缝模型又有被推广使用的趋势。 就ANSYS而言，她可以考虑分离式模型(solid65+link8，认为混凝土和钢筋粘结很好，如要考虑粘结和滑移，则可引入弹簧单元进行模拟，比较困难！)，也可采用分布式模型(带筋的solid65)。而其裂缝的处理方式则为分布裂缝模型。 二、关于本构关系 混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类，其中研究最多的是非线性弹性和弹塑性本构关系，其中不乏实用者。混凝土破坏准则从单参数到五参数模型达数十个模型，或借用古典强度理论或基于试验结果等，各个破坏准则的表达方式和繁简程度各异，适用范围和计算精度差别也比较大，给使用带来了一定的困难。 就ANSYS而言，其问题比较复杂些。 1 ANSYS混凝土的破坏准则与屈服准则是如何定义的？ 采用tb,concr,matnum则定义了W-W破坏准则(failure criterion)，而非屈服准则(yield criterion)。W-W破坏准则是用于检查混凝土开裂和压碎用的，而混凝土的塑性可以另外考虑（当然是在开裂和压碎之前）。理论上破坏准则(failure criterion)和屈服准则(yield criterion)是不同的，例如在高静水压力下会发生相当的塑性变形，表现为屈服，但没有破坏。而工程上又常将二者等同，其原因是工程结构不容许有很大的塑性变形，且混凝土等材料的屈服点不够明确，但破坏点非常明确。 定义tb,concr matnum后仅仅是定义了混凝土的破坏准则和缺省的本构关系，即W—W破坏准则、混凝土开裂和压碎前均为线性的应力应变关系，而开裂和压碎后采用其给出的本构关系。但屈服准则尚可另外定义(随材料的应力应变关系，如tb,MKIN,则定义的屈服准则是V on Mises,流动法则、硬化法则也就确定了)。
. ABQUS中的三种混凝土本构模型 ABAQUS 用连续介质的方法建立描述混凝土模型不采用宏观离散裂纹的方法描述裂纹的水平的在每一个积分点上单独计算其中。 低压力混凝土的本构关系包括： Concrete Smeared cracking model (ABAQUS/Standard) Concrete Brittle cracking model (ABAQUS/Explicit) Concrete Damage plasticity model 高压力混凝土的本构关系： Cap model 1、ABAQUS/Standard中的弥散裂缝模型Concrete Smeared cracking model (ABAQUS/Standard)：——只能用于ABAQUS/Standard中 裂纹是影响材料行为的最关键因素，它将导致开裂以及开裂后的材料的各向异性 用于描述：单调应变、在材料中表现出拉伸裂纹或者压缩时破碎的行为 在进行参数定义式的Keywords： *CONCRETE *TENSION STIFFENING *SHEAR RETENTION *FAILURE RATIOS 2、ABAQUS/Explicit中脆性破裂模型Concrete Brittle cracking model (ABAQUS/Explicit) ： 适用于拉伸裂纹控制材料行为的应用或压缩失效不重要，此模型考虑了由于裂纹引起的材料各向异性性质，材料压缩的行为假定为线弹性，脆性断裂准则可以使得材料在拉伸应力过大时失效。 在进行参数定义式的Keywords *BRITTLE CRACKING, *BRITTLE FAILURE, *BRITTLE SHEAR 3、塑性损伤模型Concrete Damage plasticity model： 适用于混凝土的各种荷载分析，单调应变，循环荷载，动力载荷，包含拉伸开裂(cracking)和压缩破碎(crushing)，此模型可以模拟硬度退化机制以及反向加载刚度恢复的混凝土力学特性 在进行参数定义式的Keywords： *CONCRETE DAMAGED PLASTICITY *CONCRETE TENSION STIFFENING *CONCRETE COMPRESSION HARDENING *CONCRETE TENSION DAMAGE *CONCRETE COMPRESSION DAMAGE 1 / 1'.
一、用钢筋混凝土简支梁的数值模拟为实例，对ANSYS的使用方法进行说明 钢筋混凝土简支梁，尺寸为长2000mm，宽150mm，高300mm。混凝土采用C30，钢筋全部采用HRB335，跨中集中荷载P作用于一刚性垫板上，垫板尺寸为长150mm，宽100mm。 建立分离式有限元模型，混凝土采用SOLID65单元，钢筋采用LINK8单元，不考虑钢筋和混凝土之间的粘结滑移。创建分离式模型时，将几何实体以钢筋位置切开，划分网格时将实体的边线定义为钢筋即可。加载点以均布荷载近似代替钢垫板，支座处则采用线约束和点约束相结合。单元尺寸以50mm左右为宜。 二、命令流 !钢筋混凝土简支梁数值分析 !分离式模型 FINISH /CLEAR /PREP7 !1.定义单元与材料属性 ET,1,SOLID65,,,,,,,1 ET,2,LINK8 MP,EX,1,13585 !混凝土材料的初始弹模以及泊松比 MP,PRXY,1,0.2 FC=14.3 !混凝土单轴抗压强度和单轴抗拉强度 FT=1.43 TB,CONCR,1 TBDA TA,,0.5,0.95,FT,-1 !定义混凝土材料及相关参数，关闭压碎 TB,MISO,1,,11 !定义混凝土应力应变曲线，用MISO模型 TBPT,,0.0002,FC*0.19 TBPT,,0.0004,FC*0.36 TBPT,,0.0006,FC*0.51 TBPT,,0.0008,FC*0.64 TBPT,,0.0010,FC*0.75 TBPT,,0.0012,FC*0.84 TBPT,,0.0014,FC*0.91 TBPT,,0.0016,FC*0.96 TBPT,,0.0018,FC*0.99 TBPT,,0.002,FC TBPT,,0.0033,FC*0.85 MP,EX,2,2.0E5 !钢筋材料的初始弹模以及泊松比 MP,PRXY,2,0.3 TB,BISO,2 TBDA TA,,300,0 !钢筋的应力应变关系，用BISO模型
ABQUS中的三種混凝土本構模型 ABAQUS?用連續介質的方法建立描述混凝土模型不采用宏觀離散裂紋的方法描述裂紋的水平的在每一個積分點上單獨計算其中。 低壓力混凝土的本構關系包括： Concrete Smeared cracking model (ABAQUS/Standard) Concrete Brittle cracking model (ABAQUS/Explicit) Concrete Damage plasticity model 高壓力混凝土的本構關系： Cap model 1、ABAQUS/Standard中的彌散裂縫模型Concrete Smeared cracking model (ABAQUS/Standard)： ——只能用于ABAQUS/Standard中 裂紋是影響材料行為的最關鍵因素，它將導致開裂以及開裂后的材料的各向異性 GAGGAGAGGAFFFFAFAF
用于描述?：單調應變?、在材料中表現出拉伸裂紋或者壓縮時破碎的行為 在進行參數定義式的Keywords： *CONCRETE *TENSION STIFFENING *SHEAR RETENTION *FAILURE RATIOS 2、ABAQUS/Explicit中脆性破裂模型Concrete Brittle cracking model (ABAQUS/Explicit) ： 适用于拉伸裂纹控制材料行为的应用或压缩失效不重要，此模型考虑了由于裂纹引起的材料各向异性性质，材料压缩的行为假定为线弹性，脆性断裂准则可以使得材料在拉伸应力过大时失效。 在进行参数定义式的Keywords *BRITTLE CRACKING, *BRITTLE FAILURE, *BRITTLE SHEAR GAGGAGAGGAFFFFAFAF
附录一 动力弹塑性分析的材料非线性参数取值 一 混凝土材料: 混凝土材料采用塑性损伤模型(Plastic-Damaged Model)(1). 根据GB 50010-2002 混凝土强度分类 如下: C25, C30, C35, C40, C45, C50, C55, C60, C65, C70, C75, C80 (1) 弹性模量: 按(2)表4.1.5, 单位kN/m 2 (2) 泊松比, 统一取 0.2 (参阅(2)的4.1.8) (3) 剪切模量: 按(2)表4.1.5中的0.4 倍采用(参阅(2)的4.1.8). (4) 密度(2): 2.5 T/m 3 (5) 单轴应力-应变关系 混凝土材料轴心抗压和轴心抗拉强度标准值按(2)表4.1.3采用. A: 单轴受压, 其应力-应变关系方程如下(参阅(2)C.2.1, P206): 当1≤x 时 32)2()23(x αx ααy a a a -+-+= 当1≥x 时 x x αx y d +-=2)1( c εεx = *= c f σy
在 0 – 0.7f c 的应力范围为线弹性, 其弹性模量按表1. 大于0.7f c 为塑性范围, 应力-塑性应变关系如下: E σεεc c in c -= B: 单轴受拉, 其应力-应变关系方程如下(参阅(2)C.2.2, P208): 当1≤x 时 62.02.1x x y -= 当1≥x 时 x x αx y t +-=7.1)1( t εεx = *= t f σy 在 0 – f t 的应力范围为线弹性, 其弹性模量按表1. 大于f t 为塑性 范围, 应力-塑性应变关系如下: E σεεt t ck t -= 据此得到下列各等级混凝土材料在拉和压屈服后的应力(kN/m 2)-塑性应变关系: *Material, Name=C25 *Concrete compression hardening 应力(kN/m 2) 塑性应变 11690., 0 16700., 0.000808693 13239.8, 0.00233739 9841.27, 0.00386389 7674.36, 0.0053464 6248.49, 0.00680245 5255.01, 0.00824305 4527.98, 0.00967414 3974.73, 0.011099 3540.4, 0.0125197 *Concrete tension stiffening 1797.8, 0 1780., 0.000025515 1191.06, 0.000135635
【原创】钢筋混凝土分离式建模方法（含ANSYS命令流） 钢筋混凝土, 分离式, 建模, ANSYS, 命令钢筋混凝土, 分离式, 建模, ANSYS, 命令 一、简介 钢筋混凝土有限元建模的方法与结果评价（前后处理），是对钢筋混凝土结构进行数值模拟的重要步骤，能否把握模型的可行性、合理性，如何从计算结果中寻找规律，是有限元理论应用于实际工程的关键一环。Blackeage以自己做过的一组钢筋混凝土暗支撑剪力墙的数值模拟为例，从若干方面提出一些经验与建议。希望大家一起讨论、批评指正(wang.jian@https://www.doczj.com/doc/c06030398.html,)。 程序：ANSYS 单元：SOLID65、BEAM188 建模方式：分离 暗支撑剪力墙结构由北京工业大学曹万林所提出，简言之就是一种在普通钢筋配筋情况下，加配斜向钢筋的剪力墙结构。 二、单元选择 以前经常采用的钢筋混凝土建模方法是通过SOLID65模拟混凝土，通过SOLID65的实常数指定钢筋配筋率，后来发现这种整体式的模型并不理想，而且将钢筋周围的SOLID65单元选择出来，再换算一个等效的配筋率，工作量也并不小。最关键的是采用整体式模型之后，得不出什么有意义的结论，弄一个荷载-位移曲线出来又和实验值差距比较大。只有计算的开裂荷载与实验还算是比较接近，但这个 手算也算得出来的东西费劲去装模作样的建个模型又有什么意义？ 所以，这次我尝试采用分离式的模型，钢筋与混凝土单元分别建模，采用节点共享的方式。建模时发现，只要充分、灵活地运用APDL 的技巧，处理好钢筋与混凝土单元节点的位置，效率还是很高的。 [center] 暗支撑剪力墙数值模型[/center] 看过很多的资料，分离式模型是用LINK8与SOLID65的组合方式，这样做到是非常直观，因为LINK8是spar类型的单元，每个节点有3个自由度，这与SOLID65单元单节点自由度数量是一致的。但是问题也就由此产生，当周围的混凝土开裂或是压碎时，SOLID65将不能对LINK8的节点提供足够地约束（如下图箭头方向），从而导致总刚矩阵小主元地出现影响计算精度，或者干脆形成瞬变体系导致计 算提前发散。 [center] LINK8+SOLID65的问题[/center] 如果采用梁单元模拟暗钢筋，就算包裹钢筋的混凝土破坏了，钢筋单元本身仍可对连接点提供一定的侧向刚度（其实钢筋本身就是有一定抗弯刚度的），保证计算进行下去。ANSYS中的梁单元比较多，建议选取beam188单元。beam188支持弹塑性分析、自定义截面。可以用内力计算结果按截面插值得出应力结果，这样，SOLID65+beam188不仅解决了SOLID65+beam188的小主元问题，而且可以方便地控 制钢筋单元的划分密度，也扩充了钢筋单元输出信息。 三、单元组合方式 将剪力墙中所有钢筋单元（包括暗柱、梁的纵、箍筋、暗支撑钢筋、暗支撑箍筋、暗分布筋）单独建模，为了能够与混凝土单元节点共享，将混凝土单元细化，单元高度设为暗柱箍筋间距与墙片分布筋间距的最大公约数。 钢筋与混凝土单元节点共享。不考虑粘接-滑移影响。其实由于混凝土单元已经细化过了，钢筋周围的混凝土由于钢筋作用而开裂之
ABAQUS钢筋混凝土损饬塑性模型有限元分析 发表时间：2009-10-12 刘劲松刘红军来源：万方数据 钢筋混凝土材料，是一种非匀质的力学性能复杂的建筑材料。随着计算机和有限元方法的发展，有限元法已经成为研究混凝土结构的一个重要的手段。由于数值计算具有快速、代价低和易于实现等诸多优点，这种分析方法已经广泛用于实际工程中。然而，要在有限元软件中尽可能准确地模拟混凝土这种材料，是不容易的，国内外学者提出了基于各种理论的混凝土本构模型。但是迄今为止，还没有一种理论被公认为可以完全描述混凝土的本构关系。 ABAQUS是大型通用的有限元分析软件，其在非线性分析方面的巨大优势，获得了广大用户的认可，在结构分析领域的应用趋于广泛。本文把规范建议的混凝土本构关系，应用到损伤塑性模型，对一悬臂梁进行了精细的有限元建模计算和探讨。 1 混凝土损伤塑性模型 ABAQUS在钢筋混凝土分析上有很强的能力。它提供了三种混凝土本构模型：混凝土损伤塑性模型，混凝土弥散裂缝模型和ABAQUS/Explicit中的混凝土开裂模型。其中混凝土损伤塑性模型可以用于单向加载、循环加载以及动态加载等场合，它使用非关联多硬化塑性和各向同性损伤弹性相结合的方式描述了混凝土破碎过程中发生的不可恢复的损伤。这一特性使得损伤塑性模型具有更好的收敛性。 2 模型材料的定义 2.1 混凝土的单轴拉压应力-应变曲线 本模型中选用的混凝土本构关系是《混凝土结构设计规范》所建议的曲线，其应力应变关系可由函数表达式定义。 2.2 钢筋的本构关系 钢筋采用本构关系为强化的二折线模型，无刚度退化。折线第一上升段的斜率，为钢筋本身的弹性模量，第二上升段为钢筋强化段，此时的斜率大致可取为第一段的1/100。 2.3 损伤的定义 损伤是指在单调加载或重复加载下，材料性质所产生的一种劣化现象，损伤在宏观方面的表现就是(微)裂纹的产生。材料的损伤状态，可以用损伤因子来描述。根据前面确定的混凝土非弹性阶段的应力一应变关系。可求得损伤因子的数值。 2.4混凝土塑性数值的计算 混凝土在单向拉伸，压缩试验中得到的数据，通常是以名义应变和名义应力表示的，为了准确地描述大变形过程中截面积的改变，需要使用真实应变和真实应力，可通过它们之间的换算公式计算。真实应变是由塑性应变和弹性应变两部分构成的。在ABAQUS中定义塑性材料参数时，需要使用塑性应变。 3 钢筋混凝土悬臂梁实例分析 3.1 模型设计 该悬臂梁的具体情况如图1所示，梁截面尺寸为200mm×300mm，梁长1500mm；纵筋为HRB335钢筋，箍筋为HPB235钢筋，混凝土强度等级为C30。混凝土和钢筋的各力学参数均取自《混凝土结构设计规范》的标准值。
ABQUS中的三种混凝土本构模型 ABAQUS 用连续介质的方法建立描述混凝土模型不采用宏观离散裂纹的方法描述裂纹的水平的在每一个积分点上单独计算其中。 低压力混凝土的本构关系包括： Concrete Smeared cracking model (ABAQUS/Standard) Concrete Brittle cracking model (ABAQUS/Explicit) Concrete Damage plasticity model 高压力混凝土的本构关系： Cap model 1、ABAQUS/Standard中的弥散裂缝模型Concrete Smeared cracking model (ABAQUS/Standard)：——只能用于ABAQUS/Standard中 裂纹是影响材料行为的最关键因素，它将导致开裂以及开裂后的材料的各向异性 用于描述：单调应变、在材料中表现出拉伸裂纹或者压缩时破碎的行为 在进行参数定义式的Keywords： *CONCRETE *TENSION STIFFENING *SHEAR RETENTION *FAILURE RATIOS 2、ABAQUS/Explicit中脆性破裂模型Concrete Brittle cracking model (ABAQUS/Explicit) ： 适用于拉伸裂纹控制材料行为的应用或压缩失效不重要，此模型考虑了由于裂纹引起的材料各向异性性质，材料压缩的行为假定为线弹性，脆性断裂准则可以使得材料在拉伸应力过大时失效。 在进行参数定义式的Keywords *BRITTLE CRACKING, *BRITTLE FAILURE, *BRITTLE SHEAR 3、塑性损伤模型Concrete Damage plasticity model： 适用于混凝土的各种荷载分析，单调应变，循环荷载，动力载荷，包含拉伸开裂(cracking)和压缩破碎(crushing)，此模型可以模拟硬度退化机制以及反向加载刚度恢复的混凝土力学特性 在进行参数定义式的Keywords： *CONCRETE DAMAGED PLASTICITY *CONCRETE TENSION STIFFENING *CONCRETE COMPRESSION HARDENING *CONCRETE TENSION DAMAGE *CONCRETE COMPRESSION DAMAGE