如何在ANSYS中模拟钢筋混凝土的计算模型
如何在ANSYS中模拟钢筋混凝土的计算模型
如何在ANSYS中模拟钢筋混凝土的计算模型最近做了点计算分析，结合各论坛关于这方面的讨论，就一些问题探讨如下，不当之处敬请指正。
一、关于模型钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种，即分离式、分布式和组合式模型。
考虑钢筋和混凝土之间的粘结和滑移，则采用引入粘结单元的分离式模型；假定混凝土和钢筋粘结很好，不考虑二者之间的滑移，则三种模型都可以；分离式和分布式模型适用于二维和三维结构分析，后者对杆系结构分析比较适用。
裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型，后者目前尚处研究之中，主要应用的是前两种。
离散裂缝模型和分布裂缝模型各有特点，可根据不同的分析目的选择使用。
随着计算速度和网格自动划分的快速实现，离散裂缝模型又有被推广使用的趋势。
就ANSYS而言，她可以考虑分离式模型(solid65+link8，认为混凝土和钢筋粘结很好，如要考虑粘结和滑移，则可引入弹簧单元进行模拟，比较困难！)，也可采用分布式模型(带筋的solid65)。
而其裂缝的处理方式则为分布裂缝模型。
二、关于本构关系混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类，其中研究最多的是非线性弹性和弹塑性本构关系，其中不乏实用者。
混凝土破坏准则从单参数到五参数模型达数十个模型，或借用古典强度理论或基于试验结果等，各个破坏准则的表达方式和繁简程度各异，适用范围和计算精度差别也比较大，给使用带来了一定的困难。
就ANSYS而言，其问题比较复杂些。
1 ANSYS混凝土的破坏准则与屈服准则是如何定义的？采用tb,concr,matnum则定义了W-W破坏准则(failure criterion)，而非屈服准则(yield criterion)。
W-W破坏准则是用于检查混凝土开裂和压碎用的，而混凝土的塑性可以另外考虑（当然是在开裂和压碎之前）。
理论上破坏准则(failure criterion)和屈服准则(yield criterion)是不同的，例如在高静水压力下会发生相当的塑性变形，表现为屈服，但没有破坏。
钢筋混凝土建模步骤
钢筋混凝土建模步骤在土木工程结构中，最为常用的一种结构形式就是钢筋混凝土结构，在各类房屋、水坝、桥梁、道路中都有广泛应用。
ANSYS软件提供了专门的钢筋混凝土单元和材料模型。
本算例将介绍ANSYS软件分析混凝土一些基本应用。
(1) 首先建立有限元模型，这里我们选用ANSYS软件自带的专门针对混凝土的单元类型Solid 65，进入ANSYS主菜单Preprocessor->ElementType->Add/Edit/Delete，选择添加Solid 65号混凝土单元。
(2) 点击Element types窗口中的Options，设定Stress relax after cracking 为Include，即考虑混凝土开裂后的应力软化行为，这样在很多时候都可以提高计算的收敛效率。
(3) 下面我们要通过实参数来设置Solid 65单元中的配筋情况。
进入ANSYS主菜单Preprocessor-> Real Constants->Add/Edit/Delete，添加实参数类型1与Solid 65单元相关，输入钢筋的材料属性为2号材料，但不输入钢筋面积，即这类实参数是素混凝土的配筋情况。
(4) 再添加第二个实参数，输入X方向配筋为0.05，即X方向的体积配筋率为5%。
(5) 下面输入混凝土的材料属性。
混凝土的材料属性比较复杂，其力学属性部分一般由以下3部分组成：基本属性，包括弹性模量和泊松比；本构关系，定义等效应力应变行为；破坏准则，定义开裂强度和压碎强度。
下面分别介绍如下。
(6) 首先进入ANSYS主菜单Preprocessor-> Material Props-> Material Models，在Define Material Model Behavior 窗口中选择Structural-> Linear -> Elastic-> Isotropic，输入弹性模量和泊松比分别为30e9和0.2(7) 下面输入混凝土的等效应力应变关系，这里我们选择von Mises屈服面，该屈服面对于二维受力的混凝土而言精度还是可以接受的。
手把手教你学ansys--钢筋混凝土梁
大家好，我是水哥。
水哥ansys使用经验三年多，既做过重大科研项目，也做过许多实际项目，对ansys的使用有一定的心得体会，本着分享经验的精神，今日以一个钢筋混凝土梁的建模求解过程来简单说明ansys的基本操作步骤。
（我的ansys14.0）总的说来，无论小项目还是大项目，总体过程无非定义单元、定义材料、物理模型、有限元模型、加载、约束、求解、查看这几个过程，和我们工程类的设计软件例如PKPM、SAP2000等基本过程都差不多，只不过最大的区别在于ansys的建模实在是蛋疼了一些。
废话不多说，以下面的一根悬臂钢筋混凝土梁来教新手如何快速进入ansys 的大堂，每一步都有GUI操作，完了之后会有相应部分的命令流，这里多说一句，一个ansys的使用高手必然是一个精通apdl编程的能手，所以我建议新手在学习的时候最好以apdl入手，GUI操作辅助，这样在学习的时候能节省大量时间，而不会浪费在GUI毫无意义的重复操作上。
此题如下：悬臂梁如下，梁宽200mm，梁顶有两根直径为16的钢筋，钢筋中心距梁边的距离为40mm，在梁端附近受集中力P=100KN的作用.要求对此悬臂梁进行完全线弹性分析，结果要显示主应力迹线。
材料参数：混凝土弹性模量为3000MPa,钢筋的弹性模量取200GPa，不考虑材料自重。
（建模时注意单位的协调性）一、题目解读与材料单元定义注意此题要求进行完全线弹性分析，此话的意义在于我们可以用除solid65 以外的其他实体单元。
在ansys单元中，solid65是专门用于模拟钢筋混凝土构件的单元，但用此单元模拟时，一般是需要考虑材料的非线性，也即是多用于构件的非线性分析，并且需要材料的本构模型。
此题要求我们做弹性分析，我们可以用link8单元模拟钢筋，用solid45单元模拟混凝土，注意高版本的ansys已经将许多低阶单元合并掉了，以ansys14为例，在link单元中只有180，而低阶的link8、link10等已被合并。
(完整版)ansys钢筋混凝土梁的建模方法约束方程法
用约束方程法模拟钢筋混凝土梁结构问题描述建立钢筋线对钢筋线划分网格后形成钢筋单元bhP 位移载荷建立混凝土单元对钢筋线节点以及混凝土节点之间建立约束方程后施加约束以及位移载荷进入求解器进行求解；钢筋单元的受力云图混凝土的应力云图混凝土开裂fini/clear,nostart/config,nres,5000/filname,yue su fang cheng 5 jia mi hun nin tu /prep7/title,rc-beamb=150h=300a=30l=2000displacement=5!定义单元类型et,1,solid65et,2,beam188et,3,plane42!定义截面类型sectype,1,beam,csolid,,0secoffset,centsecdata,8,0,0,0,0,0,0,0,0,0sectype,2,beam,csolid,,0secoffset,centsecdata,4,0,0,0,0,0,0,0,0,0!定义材料属性，混凝土材料属性mp,ex,1,24000mp,prxy,1,0.2tb,conc,1,1,9tbdata,,0.4,1,3,-1!纵向受拉钢筋mp,ex,2,2e5mp,prxy,2,0.3tb,bkin,2,1,2,1tbdata,,350!横向箍筋，受压钢筋材料属性mp,ex,3,2e5mp,prxy,3,0.25tb,bkin,3,1,2,1tbdata,,200!生成钢筋线k,,k,,bkgen,2,1,2,,,hk,,a,ak,,b-a,akgen,2,5,6,,,h-2*akgen,21,5,8,,,,-100 *do,i,5,84,1l,i,i+4*enddo*do,i,5,85,4l,i,i+1l,i,i+2*enddo*do,i,8,88,4l,i,i-1l,i,i-2*enddo!受拉钢筋lsel,s,loc,y,alsel,r,loc,x,alsel,a,loc,x,b-a lsel,r,loc,y,acm,longitudinal,line type,2mat,2secnum,1 lesize,all,50lmesh,allallscmsel,u,longitudinalcm,hooping reinforcement,line!箍筋，受压钢筋type,2mat,2secnum,2lesize,all,50lmesh,all/eshape,1!将钢筋节点建为一个集合cm,steel,node!生成面单元，以便拉伸成体单元a,1,2,4,3lsel,s,loc,y,0lsel,a,loc,y,hlesize,all,,,10lsel,alllsel,s,loc,x,0lsel,a,loc,x,blesize,all,,,20type,3amesh,all!拉伸成混凝土单元type,1real,3mat,1extopt,esize,30extopt,aclear,1vext,all,,,,,-lalls!建立约束方程cmsel,s,hooping reinforcement cmsel,a,longitudinalnsll,s,1ceintf,,ux,uy,uzallsel,all!边界条件约束nsel,s,loc,y,0nsel,r,loc,z,0d,all,uyd,all,uxnsel,s,loc,y,0nsel,r,loc,z,-ld,all,uyd,all,ux!施加外部荷载/solunsel,allnsel,s,loc,y,hnsel,r,loc,z,-1000d,all,uy,-displacement alls!求解nlgeom,on nsubst,200 outres,all,all neqit,100pred,oncnvtol,f,,0.05,2,0.5 allselsolvefinish/post1allselplcrack,0,1plcrack,0,2!时间历程后处理/post26nsel,s,loc,z,-l/2*get,Nmin,node,0,num,min nsol,2,nmin,u,yprod,3,2,,,,,,-1nsel,s,loc,y,0nsel,r,loc,z,0*get,Nnum,node,0,count *get,Nmin,node,0,num,min n0=Nminrforce,5,Nmin,f,y*do,i,2,ndinqr(1,13)ni=ndnext(n0)rforce,6,ni,f,yadd,5,5,6n0=ni*enddoprod,7,5,,,,,,1/1000/axlab,x,uy/axlab,y,p(kn) xvar,3 plvar,7。
ANSYS分析钢筋混凝土结构技巧及实例详解
0 前言利用ANSYS分析钢筋混凝土结构时，其有限元模型主要有分离式和整体式两种模型。
这里结合钢筋混凝土材料的工作特性，从模型建立到非线性计算再到结果分析的全过程讲述了利用ANSYS进行钢筋混凝土结构分析的方法与技巧，并以钢筋混凝土简支梁为例，采用分离式有限元模型，说明其具体应用。
1 单元选取与材料性质1. 1 混凝土单元ANSYS中提供了上百种计算单元类型,其中Solid65单元是专门用于模拟混凝土材料的三维实体单元。
该单元是八节点六面体单元,每个节点具有三个方向的自由度( UX , UY , UZ) 。
在普通八节点线弹性单元Solid45 的基础上,该单元增加了针对于混凝土的材性参数和组合式钢筋模型,可以综合考虑包括塑性和徐变引起的材料非线性、大位移引起的几何非线性、混凝土开裂和压碎引起的非线性等多种混凝土的材料特性。
使用Solid65 单元时,一般需要为其提供如下数据:1)、实常数(Real Constants) :定义弥散在混凝土中的最多三种钢筋的材料属性,配筋率和配筋角度。
对于墙板等配筋较密集且均匀的构件,一般使用这种整体式钢筋混凝土模型。
如果采用分离式配筋,那么此处则不需要填写钢筋实常数。
2)、材料模型(Material Model) :在输入钢筋和混凝土的非线性材料属性之前,首先必须定义钢筋和混凝土材料在线弹性阶段分析所需的基本材料信息,如:弹性模量,泊松比和密度。
3)、数据表(Data Table) :利用数据表进一步定义钢筋和混凝土的本构关系。
对于钢筋材料,一般只需要给定一个应力应变关系的数据表就可以了,譬如双折线等强硬化(bilinear isotropic hardening)或随动硬化模型( kinematic hardening plasticity)等。
而对于混凝土模型,除需要定义混凝土的本构关系外,还需要定义混凝土材料的破坏准则。
在ANSYS中,常用于定义混凝土本构关系的模型有:1)多线性等效强化模型(Multilinear isotropic hardening plas2ticity ,MISO模型)，MISO模型可包括20条不同温度曲线，每条曲线可以有最多100个不同的应力-应变点;2)多线性随动强化模型(Multilinear kinematic hardening plas2ticity ,MKIN 模型)，MKIN 模型最多允许5个应力-应变数据点;3)Drucker2Prager plasticity(DP)模型。
ANSYS应用实例：钢筋混凝土简支梁数值模拟
（ii ）纵向钢筋：PIPE20 （iii ）横向箍筋：PIPE202.2 材料性质（i ）、混凝土材料表5-4 混凝土材料的输入参数一览表[16~19]·单轴受压应力-应变曲线（εσ-曲线）在ANSYS ○R程序分析中，需要给出混凝土单轴受压下的应力应变曲线。
在本算例中，混凝土单轴受压下的应力应变采用Sargin 和Saenz 模型[17,18]：221⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛-+=c c s c c E E E εεεεεσ （5-30）式中取4'4')108.0028.1(c c c f f -=ε；断面图配筋图断面图配筋图断面图配筋图RCBEAM-01 RCBEAM-02 RCBEAM-03图5-12 各梁FEM模型断面图（a）单元网格图（b）钢筋单元划分图图5-13 算例（一）的FEM模型图2.4 模型求解在ANSYS○R程序中，对于非线性分析，求解步的设置很关键，对计算是否收敛关系很大，对于混凝土非线性有限元分析，在计算时间容许的情况下，较多的求解子步（Substeps）或较小的荷载步和一个非常大的最大子步数更容易导致收敛[2]。
在本算例中，设置了100个子步。
最终本算例收敛成功，在CPU为P41.6G、内存为256MB的微机上计算，耗时约为8小时。
2.5 计算结果及分析2.5.1 荷载—位移曲线图5-14为ANSYS○R程序所得到的各梁的荷载-跨中挠度曲线，从图中可以看出：（i）、梁RCBEAM-01：曲线形状能基本反映钢筋混凝土适筋梁剪切破坏的受力特点，而且荷载-跨中挠度曲线与钢筋混凝土梁的弯剪破坏形态非常类似，即当跨中弯矩最大截面的纵筋屈服后，由于裂缝的开展，压区混凝土的面积逐渐减小，在荷载几乎不增加的情况下，压区混凝土所受的正应力和剪应力还在不断增加，当应力达到混凝土强度极限时，剪切破坏发生，荷载突然降低。
（ii）、梁RCBEAM-02：荷载-跨中挠度曲线与超筋梁的试验荷载-跨中挠度曲线很相似，在荷载达到极限情况下，没有出现屈服平台，而是突然跌落。
钢筋混凝土简支梁模拟计算 ansys命令流
lsel,r,loc,y,0
dl,all,,uy
dl,all,,uz
lsel,all
lsel,s,loc,z,3பைடு நூலகம்00
lsel,r,loc,y,0
dl,all,,uy
lsel,all
ksel,s,loc,x,0
ksel,r,loc,y,0
!方法：分离式；solid65和link8
!材料：混凝土采用concr和钢筋为弹性材料，但不考虑压碎
!---------------------------------------------------
!为方便，假定钢筋置于梁底两侧.
!===================================================
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
!施加约束
/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
/post1
pldisp,1
钢筋混凝 ansys算例
钢筋混凝土结构开裂计算方案1．概述开裂计算是工程中比较关心的问题，但一直是有限元分析的一个难点，涉及到材料本构、计算收敛性等诸多问题。
ANSYS ＋CivilFEM 提供了钢筋混凝土结构开裂计算功能，其中土木专用模块CivilFEM 提供的非线性混凝土计算适用于混凝土梁结构的非线性计算（包括开裂），可以直接通过截面定义钢筋，从而模拟钢筋混凝土梁。
但对于更一般的结构，用梁单元来模拟不一定合适，需要采用更一般的单元，ANSYS 提供了专用的钢筋混凝土实体单元SOLID65来模拟钢筋混凝土结构，该单元材料采用混凝土材料模型，可定义混凝土的开裂、压碎准则。
另外可以定义钢筋方向和体积率，可用来模拟钢筋混凝土的破坏。
本文将通过算例对ANSYS ＋CivilFEM 开裂计算的效果进行探讨，并针对一些计算难点提出初步的解决方案。
2．CivilFEM 开裂计算CivilFEM 适合于梁结构开裂分析，另外为了与后面SOLID65单元开裂计算结果进行比较，先探讨了CivilFEM 的开裂计算。
CivilFEM 开裂计算需要考虑的要点：1、激活CivilFEM 非线性模块（~CFACTIV ,NLC,Y ），这是CivilFEM 非线性计算的前提。
2、即使事实上为小变形，也必须打开几何非线性效应（NLGEOM,ON ），否则无法激活非线性迭代。
3、通常应该关闭求解控制（SOLCONTROL,OFF ），由于CivilFEM 非线性计算通过修改实常数的等效方法，自动求解控制反而可能导致发散。
4、在收敛不好的情况下，可以增加子步数、打开自动步长（AUTOTS,ON ）或可以给定一个比较大的迭代数（NEQIT,NUM ）,以改善收敛，线性搜索有时也可以改善收敛（LNSRCH,ON ）。
5、有些情况下上述调整可能仍然无法保证收敛，这通常发生在一些开裂、受压区状态转换的临界点，尤其在动力分析中更易出现，可以结合两个办法克服，一是放松收敛准则（CNVTOL ），开裂分析状态变化剧烈，往往是接近收敛但出现振荡，放松收敛可以保证在较松的准则下收敛，但可得到足以满足要求的结果。
用ANSYS建立钢筋混凝土梁模型
用ANSYS 建立钢筋混凝土梁模型问题描述：钢筋混凝土梁在受到中间位移荷载的条件下的变形以及个组成部分的应力情况。
一、用合并节点的方法模拟钢筋混凝土梁1．用solid65号单元以及beam188单元时材料特性钢材的应力应变关系混凝土的弹性模量采用线弹性B=150mmH=300mm图1钢筋混凝土结构尺寸图建立钢筋线对钢筋线划分网格后形成钢筋单元建立混凝土单元合并单元节点后施加约束以及位移载荷进入求解器进行求解钢筋单元的受力云图混凝土的应力云图混凝土开裂2 使用单元solid45号单元与beam188 钢筋的应力应变关系不变，而混凝土应力应变关系为：混凝土单元钢筋单元力与位移曲线合并节点时的命令流:fini/clear,nostart/config,nres,5000/prep7/title,rc-beamb=150h=300a=30l=2000fcu=40ec=2.85e4displacement=10!定义单元类型et,1,solid45et,2,beam188et,3,plane42!定义截面类型sectype,1,beam,csolid,,0 secoffset,centsecdata,8,0,0,0,0,0,0,0,0,0 sectype,2,beam,csolid,,0 secoffset,centsecdata,4,0,0,0,0,0,0,0,0,0!定义材料属性，混凝土材料属性mp,ex,1,ecmp,prxy,1,0.2tb,kinh,1,,16tbpt,,0.000179067,5.10tbpt,,0.000358133,9.67tbpt,,0.0005372,1.37e1tbpt,,0.000716267,1.72e1 tbpt,,0.000895333,2.01e1 tbpt,,0.0010744,2.26e1tbpt,,0.001253467,2.44e1 tbpt,,0.001432533,2.58e1 tbpt,,0.0016116,2.66e1tbpt,,0.001790667,2.69e1 tbpt,,0.0019916,2.65e1tbpt,,0.002393467,2.57e1 tbpt,,0.002795333,2.48e1 tbpt,,0.0031972,2.40e1 tbpt,,0.003599067,2.32e1 tbpt,,0.0038,2.28e1tb,conc,1,1,9tbdata,,0.4,1,3,-1!纵向受拉钢筋mp,ex,2,2e5mp,prxy,2,0.3tb,bkin,2,1,2,1tbdata,,350!横向箍筋，受压钢筋材料属性mp,ex,3,2e5mp,prxy,3,0.25tb,bkin,3,1,2,1tbdata,,200!生成钢筋线k,,k,,bkgen,2,1,2,,,hk,,a,ak,,b-a,akgen,2,5,6,,,h-2*akgen,21,5,8,,,,-100*do,i,5,84,1l,i,i+4*enddo*do,i,5,85,4l,i,i+1l,i,i+2*enddo*do,i,8,88,4l,i,i-1l,i,i-2*enddo!受拉钢筋lsel,s,loc,y,alsel,r,loc,x,alsel,a,loc,x,b-alsel,r,loc,y,acm,longitudinal,linetype,2mat,2secnum,1lesize,all,50lmesh,allallscmsel,u,longitudinalcm,hooping reinforcement,line!箍筋，受压钢筋type,2mat,2secnum,2lesize,all,50lmesh,all/eshape,1!将钢筋节点建为一个集合cm,steel,node!生成面单元，以便拉伸成体单元a,1,2,4,3lsel,s,loc,y,0lsel,a,loc,y,hlesize,all,,,8lsel,alllsel,s,loc,x,0lsel,a,loc,x,blesize,all,,,10type,3amesh,all!拉伸成混凝土单元type,1real,3mat,1extopt,esize,20extopt,aclear,1vext,all,,,,,-lalls!合并节点nummrg,allnumcmp,all!边界条件约束nsel,s,loc,y,0nsel,r,loc,z,0d,all,uyd,all,uxnsel,s,loc,y,0nsel,r,loc,z,-l d,all,uyd,all,ux!施加外部荷载/solunsel,allnsel,s,loc,y,hnsel,r,loc,z,-1000d,all,uy,-displacementalls!求解nlgeom,onnsubst,50outres,all,allneqit,50pred,oncnvtol,f,,0.05,2,0.5allselsolvefinish/post1allsel/device,vector,1!时间历程后处理/post26nsel,s,loc,z,-l/2*get,Nmin,node,0,num,min nsol,2,nmin,u,yprod,3,2,,,,,,-1nsel,s,loc,y,0nsel,r,loc,z,0*get,Nnum,node,0,count *get,Nmin,node,0,num,min n0=Nminrforce,5,Nmin,f,y*do,i,2,ndinqr(1,13)ni=ndnext(n0)rforce,6,ni,f,yadd,5,5,6n0=ni*enddoprod,7,5,,,,,,1/1000/axlab,x,uy/axlab,y,p(kn)xvar,3plvar,7二、用约束方程法模拟钢筋混凝土梁1．用solid65号单元以及beam188单元时混凝土以及钢筋采用线弹性关系：建立钢筋线对钢筋线划分网格后形成钢筋单元建立混凝土单元对钢筋线节点以及混凝土节点之间建立约束方程后施加约束以及位移载荷进入求解器进行求解；钢筋单元的受力云图混凝土的应力云图混凝土开裂2 使用单元solid45号单元与beam188 使用混凝土的本构关系曲线钢材的本构关系曲线钢筋的von mises 应力混凝土的应力用在solid45号单元下，用合并节点法、约束方程法建立模中钢筋与混凝土之间的关系的时候的一个力与位移全程曲线的比较。
ansys 钢筋混凝土建模
ansys 钢筋混凝土建模Ansys 钢筋混凝土建模在现代工程领域中，钢筋混凝土结构的应用极为广泛，从高楼大厦到桥梁隧道，从水利设施到工业厂房，无一不见其身影。
为了确保这些结构的安全性、可靠性和经济性，对其进行准确的力学分析至关重要。
Ansys 作为一款功能强大的有限元分析软件，为钢筋混凝土建模提供了高效且精确的解决方案。
钢筋混凝土是一种由钢筋和混凝土两种材料共同作用的复合材料。
混凝土具有较高的抗压强度，但抗拉强度较低；而钢筋则具有良好的抗拉性能。
在实际结构中，两者协同工作，共同承受外力。
因此，在Ansys 中进行钢筋混凝土建模时，需要准确地模拟这两种材料的特性以及它们之间的相互作用。
首先，我们来谈谈混凝土的建模。
在 Ansys 中，混凝土通常可以采用实体单元进行模拟。
对于混凝土的本构关系，我们可以选择合适的模型，如经典的混凝土损伤塑性模型（Concrete Damaged Plasticity Model）。
这个模型能够较好地考虑混凝土在受压和受拉时的非线性行为，包括混凝土的开裂、压碎等现象。
在定义混凝土的材料参数时，需要输入诸如弹性模量、泊松比、抗压强度、抗拉强度等参数。
这些参数的准确取值对于模型的准确性至关重要。
一般来说，可以通过实验测试或者参考相关的规范和标准来获取这些参数。
接下来是钢筋的建模。
钢筋在 Ansys 中有多种建模方法，常见的有两种：一种是使用杆单元（Link Element）来模拟钢筋，另一种是将钢筋嵌入到混凝土实体单元中（Embedded Element）。
使用杆单元模拟钢筋时，需要定义钢筋的截面积、弹性模量、屈服强度等参数。
这种方法计算效率较高，但对于钢筋与混凝土之间的粘结滑移行为模拟不够精确。
将钢筋嵌入到混凝土实体单元中的方法能够更准确地考虑钢筋与混凝土之间的相互作用，但计算量相对较大。
在这种方法中，需要确保钢筋单元与混凝土单元之间的节点协调。
在钢筋混凝土建模中，还需要考虑钢筋与混凝土之间的粘结滑移。
ANSYS-钢筋混凝土柱在集中力作用下ANSYS建模求解分析步骤与过程
柱承受集中力下的应力集中概念理解第二次作业规格：300mmX400mmX2000mm；端部荷载：240吨，作用在长边中心方向偏离中心100mm处；材料属性：箍筋直径d=10mm，E=2.0E5，v=0.3纵筋直径d=32mm，E=2.1E5，v=0.3混凝土E=3.0E4，v=0.167；钢板厚60mm，E=2.1E5，v=0.3钢筋保护层厚度：25mm；默认单位：长度：mm，应力：MPa，力：N单元类型：混凝土：SOLID45；钢筋：BEAM4；钢板：SHELL631 混凝土SOLID452 钢筋BEAM43 钢板SHELL63编号材料 E v1 混凝土 3.0E4 0.1672 纵筋/钢板 2.1E5 0.33 箍筋 2.0E5 0.31 纵筋截面积：804.252 箍筋截面积：78.543 钢板厚度：60按上述三表完成单元类型，材料属性和实常数的定义后，使用Creat命令创建一300mmX400mmX2000mm的长方体模拟混凝土；至于纵筋和箍筋，采用旋转工作平面进行布尔切割生成；加载时的钢板在网格划分完成后，定义加载面的属性生成。
纵筋和箍筋的生成方法如下：1：移动工作平面至相应位置：2：切割（注意切割时选择Pick All，如果没有完全切透，在单元划分时角部单元可能不是规则的六面体单元，甚至不能划分）切割命令切割完成3：划分前属性赋予在定义材料、单元及实常数时默认1号均是混凝土的属性，故在此只需指定纵筋和箍筋的属性，钢板等划分完成后再指定加载面属性，并进行划分。
选择箍筋箍筋属性赋予类似的，定义纵筋的属性：材料：2；实常数：1；单元类型：2BEAM4 4：划分单元，由于钢筋保护层厚度25mm，并且要保证后面施加力时长边偏离中心100mm处有节点，采用全局单元边长25mm。
先进行体划分，再分别选出纵筋和箍筋进行线划分，最后选择加载面，赋予属性后，进行面划分。
全局单元边长设为25mm进行体、面、线划分时要注意选择划分钢筋钢板属性赋予划分好的钢板5：使用Merge命令把相同的节点重合（划分时已控制好节点是重合的，故不用Merge也可以）6：施加约束和作用力由于在建立钢筋模型时把混凝土切成很多块，其中在加载面上共分为九块，在施加约束时要选中这九个面一起施加固接约束；划分网格时单元边长为25mm，保证了加载位置有节点。
有限元ansya分析钢筋混凝土步骤
分析过程（1）首先建立有限元模型，这里我们选用ANSYS软件自带的专门针对混凝土的单元类型Solid 65，进入ANSYS主菜单Preprocessor->Element Type->Add/Edit/Delete，选择添加Solid 65号混凝土单元。
(2) 点击Element types窗口中的Options，设定Stress relax after cracking为Include，即考虑混凝土开裂后的应力软化行为，这样在很多时候都可以提高计算的收敛效率。
(3) 下面我们要通过实参数来设置Solid 65单元中的配筋情况。
进入ANSYS主菜单Preprocessor-> Real Constants->Add/Edit/Delete，添加实参数类型1与Solid 65单元相关，输入钢筋的材料属性为2号材料，但不输入钢筋面积，即这类实参数是素混凝土的配筋情况。
(4) 再添加第二个实参数，输入X方向配筋为0.05，即X方向的体积配筋率为5%。
(5) 下面输入混凝土的材料属性。
混凝土的材料属性比较复杂，其力学属性部分一般由以下3部分组成：基本属性，包括弹性模量和泊松比；本构关系，定义等效应力应变行为；破坏准则，定义开裂强度和压碎强度。
下面分别介绍如下。
(6) 首先进入ANSYS主菜单Preprocessor-> Material Props-> Material Models，在Define Material Model Behavior 窗口中选择Structural-> Linear -> Elastic-> Isotropic，输入弹性模量和泊松比分别为30e9和0.2(7) 下面输入混凝土的等效应力应变关系，这里我们选择von Mises屈服面，该屈服面对于二维受力的混凝土而言精度还是可以接受的。
在Define Material Model Behavior 窗口中选择Structural-> Nonlinear-> Inelastic-> Rate Independent-> Isotropic Hardening Plasticity-> Mises Plasticity-> Multilinear，输入混凝土的等效应力应变曲线如下图所示。
用ANSYS对钢筋混凝土梁进行计算模拟
一、用钢筋混凝土简支梁的数值模拟为实例，对ANSYS的使用方法进行说明钢筋混凝土简支梁，尺寸为长2000mm，宽150mm，高300mm。
混凝土采用C30，钢筋全部采用HRB335，跨中集中荷载P作用于一刚性垫板上，垫板尺寸为长150mm，宽100mm。
建立分离式有限元模型，混凝土采用SOLID65单元，钢筋采用LINK8单元，不考虑钢筋和混凝土之间的粘结滑移。
创建分离式模型时，将几何实体以钢筋位置切开，划分网格时将实体的边线定义为钢筋即可。
加载点以均布荷载近似代替钢垫板，支座处则采用线约束和点约束相结合。
单元尺寸以50mm左右为宜。
二、命令流!钢筋混凝土简支梁数值分析!分离式模型FINISH/CLEAR/PREP7!1.定义单元与材料属性ET,1,SOLID65,,,,,,,1ET,2,LINK8MP,EX,1,13585 !混凝土材料的初始弹模以及泊松比MP,PRXY,1,0.2FC=14.3 !混凝土单轴抗压强度和单轴抗拉强度FT=1.43TB,CONCR,1TBDA TA,,0.5,0.95,FT,-1 !定义混凝土材料及相关参数，关闭压碎TB,MISO,1,,11 !定义混凝土应力应变曲线，用MISO模型TBPT,,0.0002,FC*0.19TBPT,,0.0004,FC*0.36TBPT,,0.0006,FC*0.51TBPT,,0.0008,FC*0.64TBPT,,0.0010,FC*0.75TBPT,,0.0012,FC*0.84TBPT,,0.0014,FC*0.91TBPT,,0.0016,FC*0.96TBPT,,0.0018,FC*0.99TBPT,,0.002,FCTBPT,,0.0033,FC*0.85MP,EX,2,2.0E5 !钢筋材料的初始弹模以及泊松比MP,PRXY,2,0.3TB,BISO,2TBDA TA,,300,0 !钢筋的应力应变关系，用BISO模型PI=ACOS(-1)R,1R,2,0.25*PI*22*22R,3,0.25*PI*10*10TBPLOT,MISO,1 !混凝土材料的数据表绘图TBPLOT,BISO,2 !钢筋材料的数据表绘图!2.创建几何模型BLC4,,,150,300,2000*DO,I,1,19 !切出箍筋位置WPOFF,,,100VSBW,ALL*ENDDOWPCSYS,-1WPOFF,,,950 !切出拟加载面VSBW,ALLWPOFF,,,100VSBW,ALLWPCSYS,-1WPROTA,,-90WPOFF,,,30VSBW,ALLWPOFF,,,240VSBW,ALLWPCSYS,-1WPOFF,30WPROTA,,,90VSBW,ALLWPOFF,,,45VSBW,ALLWPOFF,,,45VSBW,ALLWPCSYS,-1!3.划分钢筋网格ELEMSIZ=50 !网格尺寸变量，设置为50mm LSEL,S,LOC,X,30LSEL,R,LOC,Y,30LA TT,2,2,2LESIZE,ALL,ELEMSIZLMESH,ALLLSEL,S,LOC,X,75 LSEL,R,LOC,Y,30LA TT,2,2,2LESIZE,ALL,ELEMSIZ LMESH,ALLLSEL,S,LOC,X,120 LSEL,R,LOC,Y,30LA TT,2,2,2LESIZE,ALL,ELEMSIZ LMESH,ALLLSEL,S,LOC,X,30 LSEL,R,LOC,Y,270LA TT,2,3,2LESIZE,ALL,ELEMSIZ LMESH,ALLLSEL,S,LOC,X,120 LSEL,R,LOC,Y,270LA TT,2,3,2LESIZE,ALL,ELEMSIZ LMESH,ALLLSEL,S,TAN1,Z LSEL,R,LOC,Y,30,270 LSEL,R,LOC,X,30,120 LSEL,U,LOC,X,75 LSEL,U,LOC,Z,0 LSEL,U,LOC,Z,2000 LSEL,U,LOC,Z,950 LSEL,U,LOC,Z,1050 LA TT,2,3,2LESIZE,ALL,ELEMSIZ LMESH,ALLLSEL,ALL!4.划分混凝土网格V A TT,1,1,1 MSHKEY,1ESIZE,ELEMSIZ VMESH,ALLALLSEL,ALL!5.施加荷载和约束LSEL,S,LOC,Y,0LSEL,R,LOC,Z,100DL,ALL,,UYLSEL,S,LOC,Y,0LSEL,R,LOC,Z,1900DL,ALL,,UYDK,KP(0,0,100),UX,,,,UZDK,KP(0,0,1900),UXP0=180000Q0=P0/150/100ASEL,S,LOC,Z,950,1050ASEL,R,LOC,Y,300SFA,ALL,1,PRES,Q0ALLSEL,ALLFINISH!6.求解控制设置与求解/SOLUANTYPE,0NSUBST,60OUTRES,ALL,ALLAUTOS,ONNEQIT,20CNVTOL,U,,0.015ALLSEL,ALLSOLVEFINISH!7.进入POST1查看结果/POST1SET,LASTPLDISP,1 !设定最后荷载步，查看变形ESEL,S,TYPE,,2ETABLE,SAXL,LS,1PLLS,SAXL,SAXL !绘制钢筋应力图ESEL,S,TYPE,,1/DEVICE,VECTOR,ONPLCRACK,1,1 !绘制裂缝和压碎图三、计算结果图1 混凝土材料的数据表绘图图2 钢筋材料的数据表绘图图3 钢筋的模拟图4 混凝土梁的模拟图5 梁在荷载作用下Y方向上的位移图图5 梁在荷载作用下Z方向上的应力图。
ANSYS 钢筋混凝土模型
TBTEMP,0
TBDATA,,310e6,2e9,,,,
!建立模型
K,1,,,,
K,2,3,,,
KGEN,2,1,2 ,1 , ,.05, , ,0
KGEN,2,1,2 ,1 , ,0.5, , ,0
A,1,2,4,3
A,3,4,6,5
!输出每步结果
OUTRES,ALL,1
!每步至少迭代25次
CUTCONTROL,NOITERPRED,1
SOLVE
有两种方法：一种是整体式建模；另一种是分离式建模。
分离式通过网格划分，得到钢筋位置；一般采用link单元模拟钢筋。
举例说明：
整体式命令流如下：
/PREP7
! 单元属性
ET,1,SOLID65
KEYOPT,1,1,0
KEYOPT,1,5,0
KEYOPT,1,6,0
KEYOPT,1,7,1
!破坏准则
TB,CONC,1,1,9,
TBTEMP,0
TBDATA,,.5,.9,3e6,30e6,,
TBDATA,,,,1,,,
MPTEMP,,,,,,,,
MPTEMP,1,0
!钢材基本属性
MPDATA,EX,2,,200e9
MPDATA,PRXY,2,,.27
!屈服准则
D,35, , , , , ,ALL, , , , ,
D,48, , , , , ,UZ, , , , ,
!加竖向位移荷载
D,75, ,-0.01, , , ,UY, , , , ,
D,1230步
NSUBST,20,1000000,20
关于ANSYSY中钢筋混凝土模型建立的探讨
关于ANSYSY中钢筋混凝土模型建立的探讨摘要：本文主要以ansys ls-dyna进行钢筋混凝土数值模拟过程中，模型的建立与选择中时常遇到的各种问题进行探讨。
包括SolidWorks软件及hypermesh软件的联合应用及仿真模拟。
尽可能的减少繁琐的步骤，更快更精准的模拟出结果关键词：ansys ls-dyna；SolidWorks；hypermesh；仿真模拟1.钢筋混凝土模型建立的概况在我们日常用ANSYS进行数值模拟时，经常会遇到钢筋混凝土结构。
要满足工程实际情况，更加真实地反映真实地工效果。
就需要我们正确的建立钢筋混凝土模型，这是关键的第一步。
2.SolidWorksSolidWorks是建立钢筋混凝土模型的最简单的三维模型建立的工具。
当我们进行复杂模型的建立时，ansys自带的建模界面并不符合国人的建模习惯，也不简便，而SolidWorks的中文界面和强大的建模程序，可以大大降低工作量节省很多的时间。
并且，导入ansys中后可以完美结合。
并不会出现什么错误。
3.Hypermesh当我们用ANSYS软件进行网格划分时，会出现各种小问题，对于复杂模型特别是钢筋混凝土模型划分是十分不利的。
Hypermesh软件的强大的划分网格的功能就可以弥补这个缺陷。
4.三个软件的联合应用我们首先用SolidWorks软件进行钢筋混凝土模型的建立，然后倒入hypermesh中进行网格划分及材料模型定义，边界条件的设立，荷载的施加等等。
最后导入ansys ls-dyna求解器进行求解。
这三者的联合使用是目前最方便快捷的方法。
ANSYS分析流程图5.钢筋混凝土结构有限元模型根据钢筋的处理方法，钢筋混凝土有限元模型主要分为三种类型，即分离，整体和组合模型。
5.1分离模型分离的模型将混凝土和钢筋视为不同的单位，即将混凝土和钢筋分为足够小的单位，并分别求解两者的刚度矩阵。
作为细长材料，通常可以忽略横向剪切强度并将其视为线性元素。
ANSYS钢筋混凝土
ANSYS钢筋混凝土概述ANSYS是一种强大的工程模拟软件，可用于钢筋混凝土的分析和设计。
本文档将介绍如何使用ANSYS进行钢筋混凝土建模、分析和评估。
建模在ANSYS中建立钢筋混凝土模型的第一步是创建几何模型。
可以使用ANSYS中的几何建模工具，或者从其他CAD软件导入一个现有的几何模型。
确定好模型的尺寸和形状后，可以开始定义材料和断面。
材料定义在ANSYS中，钢筋混凝土可以由两种材料组成：混凝土和钢筋。
可以通过输入混凝土和钢筋的弹性模量、抗拉强度、抗压强度等材料属性来定义它们。
还可以定义混凝土的增量模型，以考虑非线性行为。
断面定义在建立钢筋混凝土模型时，需要定义构件的断面属性，包括形状和尺寸。
可以选择矩形、圆形、T型等断面形状，并输入各个参数。
此外，还可以定义钢筋的位置和数量。
分析建立了钢筋混凝土模型后，可以进行各种分析，包括静力分析、动力分析和热力分析。
静力分析静力分析用于评估结构在外部载荷作用下的行为。
可以施加不同类型的载荷，比如集中载荷、分布载荷、温度载荷等。
通过静力分析可以计算出结构的位移、应力和应变，并评估结构的安全性。
动力分析动力分析用于评估结构在地震、风载等动力荷载作用下的行为。
可以施加模拟实际工况的时间历程载荷，分析结构的振动特性、动态应力等。
动力分析可以帮助工程师设计出更稳定和抗震的结构。
热力分析热力分析用于评估结构在温度变化下的行为。
可以考虑温度梯度引起的热应力，以及温度载荷引起的结构位移和变形。
热力分析可以用于设计具有温度变化的结构，如核电站、高温窑炉等。
评估通过ANSYS进行钢筋混凝土分析后，可以评估结构的安全性和性能。
可以根据应力和应变结果，进行疲劳分析、损伤评估，评估结构的寿命和性能。
此外，还可以进行参数化分析，改变不同的模型参数，比如尺寸、材料属性等，评估对结构的影响。
通过评估不同方案的结果，工程师可以选择最优的设计方案。
结论ANSYS是一种强大而灵活的工程模拟软件，为钢筋混凝土的建模、分析和评估提供了有效的工具。
ANSYS应用实例：钢筋混凝土简支梁数值模拟
（ii ）纵向钢筋：PIPE20 （iii ）横向箍筋：PIPE202.2 材料性质（i ）、混凝土材料表5-4 混凝土材料的输入参数一览表[16~19]·单轴受压应力-应变曲线（εσ-曲线）在ANSYS ○R程序分析中，需要给出混凝土单轴受压下的应力应变曲线。
在本算例中，混凝土单轴受压下的应力应变采用Sargin 和Saenz 模型[17,18]：221⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛-+=c c s c c E E E εεεεεσ （5-30）式中取4'4')108.0028.1(c c c f f -=ε；断面图配筋图断面图配筋图断面图配筋图RCBEAM-01 RCBEAM-02 RCBEAM-03图5-12 各梁FEM模型断面图（a）单元网格图（b）钢筋单元划分图图5-13 算例（一）的FEM模型图2.4 模型求解在ANSYS○R程序中，对于非线性分析，求解步的设置很关键，对计算是否收敛关系很大，对于混凝土非线性有限元分析，在计算时间容许的情况下，较多的求解子步（Substeps）或较小的荷载步和一个非常大的最大子步数更容易导致收敛[2]。
在本算例中，设置了100个子步。
最终本算例收敛成功，在CPU为P41.6G、内存为256MB的微机上计算，耗时约为8小时。
2.5 计算结果及分析2.5.1 荷载—位移曲线图5-14为ANSYS○R程序所得到的各梁的荷载-跨中挠度曲线，从图中可以看出：（i）、梁RCBEAM-01：曲线形状能基本反映钢筋混凝土适筋梁剪切破坏的受力特点，而且荷载-跨中挠度曲线与钢筋混凝土梁的弯剪破坏形态非常类似，即当跨中弯矩最大截面的纵筋屈服后，由于裂缝的开展，压区混凝土的面积逐渐减小，在荷载几乎不增加的情况下，压区混凝土所受的正应力和剪应力还在不断增加，当应力达到混凝土强度极限时，剪切破坏发生，荷载突然降低。
（ii）、梁RCBEAM-02：荷载-跨中挠度曲线与超筋梁的试验荷载-跨中挠度曲线很相似，在荷载达到极限情况下，没有出现屈服平台，而是突然跌落。
ANSYS 钢筋混凝土建模
ANSYS 钢筋混凝土建模一、简介钢筋混凝土有限元建模的方法与结果评价（前后处理），是对钢筋混凝土结构进行数值模拟的重要步骤，能否把握模型的可行性、合理性，如何从计算结果中寻找规律，是有限元理论应用于实际工程的关键一环。
Blackeage以自己做过的一组钢筋混凝土暗支撑剪力墙的数值模拟为例，从若干方面提出一些经验与建议。
希望大家一起讨论、批评指正(******************.cn)。
程序：ANSYS单元：SOLID65、BEAM188建模方式：分离暗支撑剪力墙结构由北京工业大学曹万林所提出，简言之就是一种在普通钢筋配筋情况下，加配斜向钢筋的剪力墙结构。
二、单元选择以前经常采用的钢筋混凝土建模方法是通过SOLID65模拟混凝土，通过SOLID65的实常数指定钢筋配筋率，后来发现这种整体式的模型并不理想，而且将钢筋周围的SOLID65单元选择出来，再换算一个等效的配筋率，工作量也并不小。
最关键的是采用整体式模型之后，得不出什么有意义的结论，弄一个荷载-位移曲线出来又和实验值差距比较大。
只有计算的开裂荷载与实验还算是比较接近，但这个手算也算得出来的东西费劲去装模作样的建个模型又有什么意义？所以，这次我尝试采用分离式的模型，钢筋与混凝土单元分别建模，采用节点共享的方式。
建模时发现，只要充分、灵活地运用APDL的技巧，处理好钢筋与混凝土单元节点的位置，效率还是很高的。
暗支撑剪力墙数值模型看过很多的资料，分离式模型是用LINK8与SOLID65的组合方式，这样做到是非常直观，因为LINK8是spar类型的单元，每个节点有3个自由度，这与SOLID65单元单节点自由度数量是一致的。
但是问题也就由此产生，当周围的混凝土开裂或是压碎时，SOLID65将不能对LINK8的节点提供足够地约束（如下图箭头方向），从而导致总刚矩阵小主元地出现影响计算精度，或者干脆形成瞬变体系导致计算提前发散。
LINK8+SOLID65的问题如果采用梁单元模拟暗钢筋，就算包裹钢筋的混凝土破坏了，钢筋单元本身仍可对连接点提供一定的侧向刚度（其实钢筋本身就是有一定抗弯刚度的），保证计算进行下去。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
如何在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,则定义的屈服准则是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。
支持此说法的还有现行铁路桥规的抗剪计算理论，以及原公路桥规的容许应力法的抗计剪计算。
2 定义混凝土的应力应变曲线单向应力应变曲线很多，常用的可参考国标混凝土结构规范，其中给出的应力应变曲线是二次曲线＋直线的下降段，其参数的设置按规范确定即可。
当然如有实测的应力应变曲线更好了。
四、关于收敛的问题ANSYS混凝土计算收敛(数值)是比较困难的，主要影响因素是网格密度、子步数、收敛准则等，这里讨论如下。
1 网格密度：网格密度适当能够收敛。
不是网格越密越好，当然太稀也不行，这仅仅是就收敛而言的，不考虑计算费用问题。
但是究竟多少合适，没有找到规律，只能靠自己针对情况慢慢试算。
2 子步数：NSUBST的设置很重要，设置太大或太小都不能达到正常收敛。
这点可以从收敛过程图看出，如果F范数曲线在[F]曲线上面走形的很长，可考虑增大nsubst。
或者根据经验慢慢调正试算。
3 收敛精度：实际上收敛精度的调正并不能彻底解决收敛的问题，但可以放宽收敛条件以加速吧。
一般不超过5%(缺省是0.5%)，且使用力收敛条件即可。
4 混凝土压碎的设置：不考虑压碎时，计算相对容易收敛；而考虑压碎则比较难收敛，即便是没有达到压碎应力时。
如果是正常使用情况下的计算，建议关掉压碎选项；如果是极限计算，建议使用concr＋MISO 且关闭压碎检查；如果必设压碎检查，则要通过大量的试算（设置不同的网格密度、NSUBST）以达到目的，但也很困难。
5 其他选项：如线性搜索、预测等项也可以打开，以加速收敛，但不能根本解决问题。
6 计算结果：仅设置concr，不管是否设置压碎，其一般P-F曲线接近二折线；采用concr+miso则P-F曲线与二折线有差别，其曲线形状明显是曲线的。
*******************************************************************************例题1!----------------------------------------------------!题目：钢筋混凝土简支梁模拟计算!方法：分离式；solid65和link8!材料：混凝土采用concr和钢筋为弹性材料，但不考虑压碎!---------------------------------------------------!为方便，假定钢筋置于梁底两侧.!===================================================/config,nres,2000/prep7!定义单元及其材料特性等rd0=20.0 !钢筋直径et,1,solid65et,2,link8mp,ex,1,33e3mp,prxy,1,0.20r,1hntra=28hntrl=2.6tb,concr,1tbdata,,0.7,1.0,hntrl,-1mp,ex,2,2.1e5mp,prxy,2,0.30r,2,acos(-1)*0.25*rd0*rd0!定义梁体即单元划分blc4, , ,100,200,3000/view,1,1,1,1/ang,1gplot!定义网分时边长控制lsel,s,loc,z,1,2999lsel,r,loc,y,0latt,2,2,2lesize,all,,,20 !钢筋网格数目lmesh,alllsel,s,loc,z,0lesize,all,,,4 !截面上的网格数目4x4 vsel,allvatt,1,1,1mshape,0,3dmshkey,1vmesh,allallsel,allfinish/solu!施加约束lsel,s,loc,z,0lsel,r,loc,y,0dl,all,,uydl,all,,uzlsel,alllsel,s,loc,z,3000lsel,r,loc,y,0dl,all,,uylsel,allksel,s,loc,x,0ksel,r,loc,y,0dk,all,uxallsel,all!施加荷载qmz=0.3asel,s,loc,y,200sfa,all,1,pres,qmzallsel,allnsubst,40outres,all,alltime,qmz*10neqit,40solvefinish/post1pldisp,1etable,zxyl,ls,1plls,zxyl,zxyl,1finish/post26nsol,2,33,u,yprod,3,1,,,,,,1/100prod,4,2,,,,,,-1xvar,4plvar,3******************************************************************************* 例题2!----------------------------------------------------!题目：钢筋混凝土简支梁模拟计算!方法：分离式；solid65和link8!材料：混凝土采用concr+Miso和钢筋为弹性材料，但不考虑压碎!增加网格密度!---------------------------------------------------!为方便，假定钢筋置于梁底两侧.!===================================================/config,nres,2000/prep7!定义单元及其材料特性等rd0=20.0 !钢筋直径et,1,solid65et,2,link8mp,ex,1,33e3mp,prxy,1,0.20r,1hntra=28hntrl=2.6tb,concr,1tbdata,,0.7,1.0,hntrl,-1tb,miso,1,,14tbpt,,0.0002,hntra*0.19tbpt,,0.0004,hntra*0.36tbpt,,0.0006,hntra*0.51tbpt,,0.0008,hntra*0.64tbpt,,0.0010,hntra*0.75tbpt,,0.0012,hntra*0.84tbpt,,0.0014,hntra*0.91tbpt,,0.0016,hntra*0.96tbpt,,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.1e5mp,prxy,2,0.30r,2,acos(-1)*0.25*rd0*rd0 !定义梁体即单元划分blc4, , ,100,200,3000/view,1,1,1,1/ang,1gplot!定义网分时边长控制lsel,s,loc,z,1,2999lsel,r,loc,y,0latt,2,2,2lesize,all,,,20lmesh,alllsel,s,loc,z,0lesize,all,,,4vsel,allvatt,1,1,1mshape,0,3dmshkey,1vmesh,allallsel,allfinish/solu!施加约束lsel,s,loc,z,0lsel,r,loc,y,0dl,all,,uydl,all,,uzlsel,alllsel,s,loc,z,3000lsel,r,loc,y,0dl,all,,uylsel,allksel,s,loc,x,0ksel,r,loc,y,0dk,all,uxallsel,all!施加荷载qmz=0.3asel,s,loc,y,200 sfa,all,1,pres,qmz allsel,alloutres,all,all time,qmz*10 nsubst,40neqit,40solvefinish/post1pldisp,1etable,zxyl,ls,1 plls,zxyl,zxyl,1 finish/post26nsol,2,33,u,y prod,3,1,,,,,,1/100 prod,4,2,,,,,,-1 xvar,4plvar,3。