如何在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 ），开裂分析状态变化剧烈，往往是接近收敛但出现振荡，放松收敛可以保证在较松的准则下收敛，但可得到足以满足要求的结果。
- 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。