如何在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使用经验三年多，既做过重大科研项目，也做过许多实际项目，对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算例
钢筋混凝土结构开裂计算方案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应用实例：钢筋混凝土简支梁数值模拟-推荐下载
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术通关，1系电过，力管根保线据护敷生高设产中技工资术艺料0不高试仅中卷可资配以料置解试技决卷术吊要是顶求指层，机配对组置电在不气进规设行范备继高进电中行保资空护料载高试与中卷带资问负料题荷试2下卷2，高总而中体且资配可料置保试时障卷，各调需类控要管试在路验最习；大题对限到设度位备内。进来在行确管调保路整机敷使组设其高过在中程正资1常料中工试，况卷要下安加与全强过，看度并22工且22作尽22下可22都能22可地护以缩1关正小于常故管工障路作高高；中中对资资于料料继试试电卷卷保破连护坏接进范管行围口整，处核或理对者高定对中值某资，些料审异试核常卷与高弯校中扁对资度图料固纸试定，卷盒编工位写况置复进.杂行保设自护备动层与处防装理腐置，跨高尤接中其地资要线料避弯试免曲卷错半调误径试高标方中高案资等，料，编试要5写、卷求重电保技要气护术设设装交备备置底4高调、动。中试电作管资高气，线料中课并敷3试资件且、设卷料中拒管技试试调绝路术验卷试动敷中方技作设包案术，技含以来术线及避槽系免、统不管启必架动要等方高多案中项；资方对料式整试，套卷为启突解动然决过停高程机中中。语高因文中此电资，气料电课试力件卷高中电中管气资壁设料薄备试、进卷接行保口调护不试装严工置等作调问并试题且技，进术合行，理过要利关求用运电管行力线高保敷中护设资装技料置术试做。卷到线技准缆术确敷指灵设导活原。。则对对：于于在调差分试动线过保盒程护处中装，高置当中高不资中同料资电试料压卷试回技卷路术调交问试叉题技时，术，作是应为指采调发用试电金人机属员一隔，变板需压进要器行在组隔事在开前发处掌生理握内；图部同纸故一资障线料时槽、，内设需，备要强制进电造行回厂外路家部须出电同具源时高高切中中断资资习料料题试试电卷卷源试切，验除线报从缆告而敷与采设相用完关高毕技中，术资要资料进料试行，卷检并主查且要和了保检解护测现装处场置理设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
有限元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 整体式钢筋混凝土建模算例
ANSYS整体式钢筋混凝土模型算例分析在土木工程结构中，最为常用的一种结构形式就是钢筋混凝土结构，在各类房屋、水坝、桥梁、道路中都有广泛应用。
ANSYS软件提供了专门的钢筋混凝土单元和材料模型。
本算例将介绍ANSYS软件分析混凝土一些基本应用。
(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屈服面，该屈服面对于二维受力的混凝土而言精度还是可以接受的。
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，保证了加载位置有节点。
ANSYS钢筋混凝土
ANSYS钢筋混凝土概述ANSYS是一种强大的工程模拟软件，可用于钢筋混凝土的分析和设计。
本文档将介绍如何使用ANSYS进行钢筋混凝土建模、分析和评估。
建模在ANSYS中建立钢筋混凝土模型的第一步是创建几何模型。
可以使用ANSYS中的几何建模工具，或者从其他CAD软件导入一个现有的几何模型。
确定好模型的尺寸和形状后，可以开始定义材料和断面。
材料定义在ANSYS中，钢筋混凝土可以由两种材料组成：混凝土和钢筋。
可以通过输入混凝土和钢筋的弹性模量、抗拉强度、抗压强度等材料属性来定义它们。
还可以定义混凝土的增量模型，以考虑非线性行为。
断面定义在建立钢筋混凝土模型时，需要定义构件的断面属性，包括形状和尺寸。
可以选择矩形、圆形、T型等断面形状，并输入各个参数。
此外，还可以定义钢筋的位置和数量。
分析建立了钢筋混凝土模型后，可以进行各种分析，包括静力分析、动力分析和热力分析。
静力分析静力分析用于评估结构在外部载荷作用下的行为。
可以施加不同类型的载荷，比如集中载荷、分布载荷、温度载荷等。
通过静力分析可以计算出结构的位移、应力和应变，并评估结构的安全性。
动力分析动力分析用于评估结构在地震、风载等动力荷载作用下的行为。
可以施加模拟实际工况的时间历程载荷，分析结构的振动特性、动态应力等。
动力分析可以帮助工程师设计出更稳定和抗震的结构。
热力分析热力分析用于评估结构在温度变化下的行为。
可以考虑温度梯度引起的热应力，以及温度载荷引起的结构位移和变形。
热力分析可以用于设计具有温度变化的结构，如核电站、高温窑炉等。
评估通过ANSYS进行钢筋混凝土分析后，可以评估结构的安全性和性能。
可以根据应力和应变结果，进行疲劳分析、损伤评估，评估结构的寿命和性能。
此外，还可以进行参数化分析，改变不同的模型参数，比如尺寸、材料属性等，评估对结构的影响。
通过评估不同方案的结果，工程师可以选择最优的设计方案。
结论ANSYS是一种强大而灵活的工程模拟软件，为钢筋混凝土的建模、分析和评估提供了有效的工具。
ANSYS中混凝土的计算问题 本构模型
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)是不同的，例如在高静水压力下会发生相当的塑性变形，表现为屈服，但没有破坏。
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钢筋混凝土分离式建模!跨中施加110KN的集中力FINISH $/CLEAR$/PREP7!AS0=380.1 $AS1=50.3 $A=30 $B=150!H=300 $L=2650 $L0=125!ET,1,SOLID65!KEYOPT,1,1,0!KEYOPT,1,5,1!KEYOPT,1,6,3!KEYOPT,1,7,1!ET,2,LINK180!ET,3,SOLID185,,3R,1,AS0$R,2,AS1$R,3MP,EX, 1,2.4E4 $MP,PRXY,1,0.2$FC=25!TB,CONCR,1,1,9!TBDATA,,0.35,0.75,3.1125,-1!TB,MISO,1,,15!TBPT,,0.0002,4.8$TBPT,,0.0004,9.375$TBPT,,0.0006,13.51! TBPT,,0.0008,17.02$TBPT,,0.001,19.83,$TBPT,,0.0012,21.95! TBPT,,0.0014,23.43$TBPT,,0.0016,24.365$TBPT,,0.0018,24.856! TBPT,,0.002,FC $TBPT,,0.0038,FC !TBPLOT!MP,EX,2,2E5$MP,PRXY,2,0.25!TB,BKIN,2$TBDATA,,360!MP,EX,3,2E5$MP,PRXY,3,0.25!TB,BKIN,3$TBDATA,,210!N,1,,B $N,9$FILL,1,9!NGEN,11,9,1,9,1,,,A!NGEN,2,1000,1,99,1,75!NGEN,3,1000,1001,1099,1,50!NGEN,7,1000,3001,3099,1,75!NGEN,7,1000,12001,12099,1,75!NGEN,2,1000,18001,18099,1,50!/VIEW,1,-1,-1,1!TYPE,2 $REAL,2$MAT,3!*DO,II,11,16,1$E,II,II+1 $*ENDDO!*DO,II,83,88,1$E,II,II+1 $*ENDDO!*DO,II,11,74,9$E,II,II+9 $*ENDDO!*DO,II,17,80,9$E,II,II+9 $*ENDDO! EGEN,20,1000,1,28,1!*DO,II,83,18083,1000$E,II,II+1000$*ENDDO!*DO,II,89,18089,1000$E,II,II+1000$*ENDDO!TYPE,2$REAL,1$MAT,2*DO,II,11,18011,1000$E,II,II+1000$*ENDDO!*DO,II,17,18017,1000$E,II,II+1000$*ENDDO!/ESHAPE,1$EP!BLC4,,,L/2,B,H $BLC4,75,,100,B,-40$WPOFFS,,,H$BLC4,625,,200,B,40 $WPCSYS,-1!WPOFFS,75$WPROTA,,,90$VSBW,ALL!WPOFFS,,,100$VSBW,ALL$WPOFFS,,,450$VSBW,ALL!WPOFFS,,,200$VSBW,ALL$WPOFFS,,,450$VSBW,ALL!WPCSYS,-1$ALLSEL!LSEL,S,LOC,Y,0$LSEL,A,LOC,Y,150$LSEL,R,LOC,X,0$LESIZE,ALL ,,,10!LSEL,S,LOC,Z,0$LSEL,A,LOC,Z,300$LSEL,R,LOC,X,0$LESIZE,AL L,,,8!LSEL,S,LOC,Z,300$LSEL,R,LOC,Y,0$LSEL,R,LOC,X,0,75$LESIZE, ALL,75!LSEL,S,LOC,Z,300$LSEL,R,LOC,Y,0$LSEL,R,LOC,X,75,175$LESIZ E,ALL,25!LSEL,S,LOC,Z,300$LSEL,R,LOC,Y,0$LSEL,R,LOC,X,175,625LSEL,S,LOC,Z,300$LSEL,R,LOC,Y,0$LSEL,R,LOC,X,625,825LSEL,S,LOC,Z,300$LSEL,R,LOC,Y,0$LSEL,R,LOC,X,825,1275LSEL,S,LOC,Z,300$LSEL,R,LOC,Y,0$LSEL,R,LOC,X,1275,1325$L ESIZE,ALL,50!LSEL,S,LOC,Z,340$LSEL,R,LOC,X,625$LESIZE,ALL,,,8!LSEL,S,LOC,Y,0$LSEL,R,LOC,X,625$LSEL,R,LOC,Z,300,340$LESI ZE,ALL,,,1!LSEL,S,LOC,Z,-40$LSEL,R,LOC,X,75$LESIZE,ALL,,,8!LSEL,S,LOC,Y,0$LSEL,R,LOC,Z,-40$LESIZE,ALL,,,4!LSEL,S,LOC,Y,0$LSEL,R,LOC,X,75$LSEL,R,LOC,Z,0,-40$LESIZE,ALL,,,1!VSEL,S,LOC,Z,0,H $VATT,1,3,1MSHAPE,0,3D $MSHKEY,1VMESH,ALL $ALLSELVSEL,S,LOC,Z,-40,0$VSEL,A,LOC,Z,H,H+40$VATT,2,3,3MSHAPE,0,3D $MSHKEY,1$VMESH,ALL/VIEW,1,-0.2,-1,1$EPLOT$ALLSELNUMMRG,ALL $NUMCMP,ALL$EPLOT/SOLU!NSEL,S,LOC,Z,-40$NSEL,R,LOC,X,L0$D,ALL,UY,,,,,UZ$ALLSEL!ASEL,S,LOC,X,L/2$DA,ALL,SYMM$ALLSEL!NSEL,S,LOC,Z,H+40$NSEL,R,LOC,X,725!*GET,NODE1,NODE,,COUNT $F,ALL,FZ,-110000/NODE1$ALLSEL!ANTYPE,STATIC$NLGEOM,ON$NSUBST,80$OUTRES,ALL,ALL$AUTOTS,1$LN SRCH,1!CNVTOL,F,,0.05,2, 0.5$ALLSEL!SOLVE $FINISH/POST1 $SET,LAST $PRRSOL,FZ!SET,LAST $PLDISP,1!ESEL,S,TYPE,,2$ETABLE,SAXL,LS,1$PLLS,SAXL,SAXL! ESEL,S,TYPE,,1$/DEVICE,VECTOR,ON$PLCRACK!/POST26$NSOL,2,NODE(L/2,B/2,0),U,Z! PROD,3,2,,,,,,-1$PROD,4,1,,,LOAD,,,110!/AXLAB,X,MID-UZ(MM)$/AXLAB,Y,P(KN)! XVAR,3$PLVAR,4。
ansys建钢管混凝土模型
ansys建钢管混凝土模型/CLEAR!(1)工作环境设置/FILNAME,CLO,1 ！指定文件名/TITLE,CLO ！指定图形标题/PREP7!(2)定义单元类型ET,1,LINK8 !定义钢筋的单元类型ET,2,SOLID65 !定义混凝土的单元类型KEYOPT,2,1,0KEYOPT,2,3,0KEYOPT,2,5,1KEYOPT,2,6,3KEYOPT,2,7,1KEYOPT,2,8,0ET,3,SOLID45 !定义钢管的单元类型KEYOPT,3,1,0KEYOPT,3,2,1KEYOPT,3,4,0KEYOPT,3,5,0KEYOPT,3,6,0ET,4,SHELL181KEYOPT,4,1,0KEYOPT,4,3,0KEYOPT,4,8,2KEYOPT,4,9,0KEYOPT,4,10,0!(3)定义实常数R,1,50.24E6, ，!定义钢筋的截面面积R,2, , , , , , ,R,3, ,R,4,0.06,0.06,0.06,0.06,0,0,!(4)定义钢材的材料模型及参数MP,EX,1,2.06E11MP,PRXY,1,0.25MP,DENS,1,7850TB,BISO,1,1,2, 1 !双线性各向同性强化模型TBDATA,,235E6,2.06E10 !理想弹塑性模型MP,EX,1,2.06E5MP,PRXY,1,0.25MP,DENS,1,7850 E-12TB,BISO,1,1,2, 1TBDATA,,235E-3,2.06E4!量纲问题，统一量纲为mm,tonne,s,oC,N,MPa !(5)定义混凝土材料MP,EX,2,32500E6MP,PRXY,2,0.173MP,DENS,2,2450TB,MISO,2,2,14,0 !多线性各向同性强化模型TBTEMP,0TBPT,,0.0000593,1.9286E6 !应力应变数据表TBPT,,0.0003,5.5782E6TBPT,,0.0006,10.5336E6TBPT,,0.0009,14.8662E6TBPT,,0.0012,18.576E6TBPT,,0.0016,22.5536E6TBPT,,0.002,25.4239E6TBPT,,0.0033,27.739E6TBPT,,0.005,22.267E6TBPT,,0.01,10.7514E6TBPT,,0.015,7.911E6TBPT,,0.02,6.658E6TBPT,,0.025,5.915E6TBPT,,0.03,5.40255E6TB,CONC,2,1,9,TBDATA,,0.3,0.9,1.71E6,-1,,,,0.9 !定义混凝土的破坏参数MP,EX,2,32500MP,PRXY,2,0.173MP,DENS,2,2450e-12TB,MISO,2,2,14,0 !多线性各向同性强化模型TBTEMP,0TBPT,,0.0000593,1.9286 !应力应变数据表TBPT,,0.0003,5.5782TBPT,,0.0006,10.5336TBPT,,0.0009,14.8662TBPT,,0.0012,18.576TBPT,,0.0016,22.5536TBPT,,0.002,25.4239TBPT,,0.0033,27.739TBPT,,0.005,22.267TBPT,,0.01,10.7514TBPT,,0.015,7.911TBPT,,0.02,6.658TBPT,,0.025,5.915TBPT,,0.03,5.40255TB,CONC,2,1,9,TBDATA,,0.3,0.9,1.71,-1,,,,0.9!(6)建立几何模型K,1,0.244,0.775,0.244K,2,-0.244,0.775,0.244K,3,0.244,0.775,-0.244K,4,-0.244 ,0.775,-0.244K,5,0.244,-0.775,0.244K,6,-0.244,-0.775,0.244K,7,0.244,-0.775,-0.244K,8,-0.244 ,-0.775,-0.244K,9,0.244,0.8,0.244K,10,-0.244,0.8,0.244K,11,0.244,0.8,-0.244K,12,-0.244 ,0.8,-0.244K,13,0.244,-0.8,0.244K,14,-0.244,-0.8,0.244K,15,0.244,-0.8,-0.244K,16,-0.244 ,-0.8,-0.244V,5,7,8,6,1,3,4,2 !建立混凝土的体1V,1,3,4,2,9,11,12,10 !建立盖板的体2和3V,13,15,16,14,5,7,8,6,VGLUE,ALL !粘结在一起/VIEW,1,1,1,1/ANG,1VSEL,,,,1 !混凝土的体1赋予属性VATT,2,2,2,0VSEL,S,,,2,3,1 !盖板的体2和3赋予属性VATT,1,3,3,0ASEL,S,,,2,5,1 !混凝土的体四个面和盖板的四个面赋予shell单元及其他的属性ASEL,A,,,7,10,1ASEL,A,,,13,16,1AATT,1,4,4,0ALLSEL,ALLLSEL,S,LINE,,5,11,2 !混凝土的竖向线划分为0.05LESIZE,ALL,0.05LSEL,A,LINE,,13,19,2 !两个盖板的竖向线划分为0.025(和其长度相等) LESIZE,ALL,0.025LSEL,A,LINE,,25,28,1LESIZE,ALL,0.025LSEL,INVE !横向截面的线划分为0.048(可整除)LESIZE,ALL,0.0488ASEL,ALLMSHAPE,0,3DMSHKEY,2 !采用映射网格划分AMESH,2,5,1 !对具有壳单元属性的面进行划分AMESH,7,10,1AMESH,13,16,1ALLSEL,ALLVMESH,1,3,1 !划分三个体NUMMRG,NODENUMMRG,ELEMNUMMRG,KPFINISH/SOLANTYPE,0 !静力分析类型NLGEOM,1 !打开大变形NROPT,AUTO,, !自动选择求解方法EQSLV,PCG,,O, !PCG求解器ALLSEL,ALL!(9)施加位移约束NSEL,S,LOC,Y,0.8D,ALL,,,,,,UX,,UZ,,,, !盖板顶加x和z向的水平约束ALLSEL,ALLNSEL,S,LOC,Y,-0.8D,ALL,,,,,,UX,UY,UZ,,, !盖板底加x和y和z向的水平约束ALLSEL,ALL!(10)施加均布荷载(分7个荷载步)ALLSEL,ALLNSEL,S,LOC,Y,0.8SF ,ALL,PRES,-10000E3TIME,1KBC,0NSUBST,10ALLSEL,ALL LSWRITE,1,NSEL,S,LOC,Y,0.8SF ,ALL,PRES,-20000E3 TIME,2KBC,0NSUBST,10 ALLSEL,ALL LSWRITE,2,NSEL,S,LOC,Y,0.8SF ,ALL,PRES,-30000E3 TIME,3KBC,0NSUBST,10 ALLSEL,ALL LSWRITE,3,NSEL,S,LOC,Y,0.8SF ,ALL,PRES,-40000E3 TIME,4KBC,0NSUBST,10 ALLSEL,ALL LSWRITE,4,NSEL,S,LOC,Y,0.8SF ,ALL,PRES,-45000E3 TIME,5KBC,0NSUBST,10 ALLSEL,ALL LSWRITE,5,NSEL,S,LOC,Y,0.8SF ,ALL,PRES,-48000E3 TIME,6KBC,0NSUBST,10ALLSEL,ALLLSWRITE,6,NSEL,S,LOC,Y,0.8SF ,ALL,PRES,-51212.76E3TIME,7KBC,0NSUBST,10ALLSEL,ALLLSWRITE,7,AUTOTS,1 !打开自动时间步长CNVTOL,F,5000,0.05,2,, !定义收敛准则CNVTOL,U,,0.03,0,,NEQIT,50, !平衡迭代次数PRED,ON,,ON !打开预测器OUTRES,ALL,ALL !输出所有子步结果LSSOLVE,1,7,1我还试过钢管也用solid45单元，建出的模型所有钢管与混凝土及盖板之间的接触部分都是共面的，因为我不考虑钢管和混凝土之间的滑移，所以我觉得接触部分是不是应该就没有问题，但是求解还是不收敛的，做了很久了，一直找不到原因，希望各位可以帮帮忙！！！由于你打开了大位移选项，NLGEOM,1 ，你关上后肯定立即可以算下去，你的几何非线性必须考虑么？考虑几何非线时，你的子步太小，NSUBST,10可能导致不收敛，建议试试增大至NSUBST,100 另NROPT,AUTO,,与EQSLV,PCG,,O, 可能是矛盾的，建议你查查手册。
ANSYS 钢筋混凝土建模
ANSYS 钢筋混凝土建模一、简介钢筋混凝土有限元建模的方法与结果评价（前后处理），是对钢筋混凝土结构进行数值模拟的重要步骤，能否把握模型的可行性、合理性，如何从计算结果中寻找规律，是有限元理论应用于实际工程的关键一环。
Blackeage以自己做过的一组钢筋混凝土暗支撑剪力墙的数值模拟为例，从若干方面提出一些经验与建议。
希望大家一起讨论、批评指正(******************.cn)。
程序：ANSYS单元：SOLID65、BEAM188建模方式：分离暗支撑剪力墙结构由北京工业大学曹万林所提出，简言之就是一种在普通钢筋配筋情况下，加配斜向钢筋的剪力墙结构。
二、单元选择以前经常采用的钢筋混凝土建模方法是通过SOLID65模拟混凝土，通过SOLID65的实常数指定钢筋配筋率，后来发现这种整体式的模型并不理想，而且将钢筋周围的SOLID65单元选择出来，再换算一个等效的配筋率，工作量也并不小。
最关键的是采用整体式模型之后，得不出什么有意义的结论，弄一个荷载-位移曲线出来又和实验值差距比较大。
只有计算的开裂荷载与实验还算是比较接近，但这个手算也算得出来的东西费劲去装模作样的建个模型又有什么意义？所以，这次我尝试采用分离式的模型，钢筋与混凝土单元分别建模，采用节点共享的方式。
建模时发现，只要充分、灵活地运用APDL的技巧，处理好钢筋与混凝土单元节点的位置，效率还是很高的。
暗支撑剪力墙数值模型看过很多的资料，分离式模型是用LINK8与SOLID65的组合方式，这样做到是非常直观，因为LINK8是spar类型的单元，每个节点有3个自由度，这与SOLID65单元单节点自由度数量是一致的。
但是问题也就由此产生，当周围的混凝土开裂或是压碎时，SOLID65将不能对LINK8的节点提供足够地约束（如下图箭头方向），从而导致总刚矩阵小主元地出现影响计算精度，或者干脆形成瞬变体系导致计算提前发散。
LINK8+SOLID65的问题如果采用梁单元模拟暗钢筋，就算包裹钢筋的混凝土破坏了，钢筋单元本身仍可对连接点提供一定的侧向刚度（其实钢筋本身就是有一定抗弯刚度的），保证计算进行下去。
基于ANSYS的钢筋混凝土结构三维实体建模技术
基于ANSYS的钢筋混凝土结构三维实体建模技术信息发布： ZXL 发布时间： 2006-4-6 16:59:21 被阅： 666 次来源：易路网文字〖大中小〗自动滚屏（右键暂停）1 引言ANSYS有限元分析软件具有功能极为强大的前后处理及计算分析能力，能够模拟结构等多种物理场间的耦合效应，可以实现对结构的“全过程仿真”和“全过程分析”。
但ANSYS中提供的材料模型大多数为基于经典材料力学理论的均质材料模型，与混凝土的本构关系有很大区别。
借助ANSYS对钢筋混凝土这种复杂材料的结构进行分析时，需要针对分析对象的结构层次、分析类型、荷载水平等合理采取单元类型和材料模型，才能够取得满意的分析结果。
结构分析的主要步骤为：①建立结构模型。
定义单元类型、赋予单元实常数构成结构几何体；定义几何体材料属性从而构成结构模型；对结构模型进行划分网格。
②施加载荷和边界条件。
根据工程具体分析要求，等效简化边界约束，布置承受荷载。
③校验、求解。
运行过程中通过参数效验、检查模型建立的准确性，及时进行修正。
④评价和分析结果。
运行结果可以通过控制输出模式，根据需要得到各种有关的力学性能参数，以表格、图形、动画等各种形式表现输出，对计算结果进行分析评价。
2 钢筋混凝土结构有限元分析单元类型通常钢筋混凝土结构有限元分析单元分为两个层次：杆系单元和实体单元。
前者着重分析单元力（包括力和弯矩）与位移（包括位移和转角）之间的关系，而后者着重分析单元的应力—应变关系。
单元类型的选取应兼顾计算规模、材料模型的精度等多方面的因素。
对于全结构规模较大，可将结构离散成杆系单元进行分析。
对于复杂区域（梁柱节点）或重要的构件等可将杆系结构体系计算的力和位移施加到实体单元模型上，分析局部应力和应变。
在结构分析中应尽可能多地采用三维实体单元模型，力求最大程度的真实模拟实际结构构件。
3 钢筋混凝土结构有限元分析中的模型钢筋混凝土结构不同于一般均质材料，它是由钢筋和混凝土两种材料构成的，一般钢筋是被包围在混凝土之中，而且相对体积较少，因此建立结构有限元模型需考虑这些特性。
ANSYS钢筋混凝土分离式建模
!跨中施加110KN的集中力FINISH $/CLEAR $/PREP7!AS0=380.1 $AS1=50.3 $A=30 $B=150!H=300 $L=2650 $L0=125!ET,1,SOLID65!KEYOPT,1,1,0!KEYOPT,1,5,1!KEYOPT,1,6,3!KEYOPT,1,7,1!ET,2,LINK180!ET,3,SOLID185,,3R,1,AS0 $R,2,AS1 $R,3MP,EX, 1,2.4E4 $MP,PRXY,1,0.2 $FC=25!TB,CONCR,1,1,9!TBDATA,,0.35,0.75,3.1125,-1!TB,MISO,1,,15!TBPT,,0.0002,4.8 $TBPT,,0.0004,9.375 $TBPT,,0.0006,13.51! TBPT,,0.0008,17.02 $TBPT,,0.001,19.83, $TBPT,,0.0012,21.95! TBPT,,0.0014,23.43 $TBPT,,0.0016,24.365 $TBPT,,0.0018,24.856! TBPT,,0.002,FC $TBPT,,0.0038,FC !TBPLOT!MP,EX,2,2E5 $MP,PRXY,2,0.25!TB,BKIN,2 $TBDATA,,360!MP,EX,3,2E5 $MP,PRXY,3,0.25!TB,BKIN,3 $TBDATA,,210!N,1,,B $N,9 $FILL,1,9!NGEN,11,9,1,9,1,,,A!NGEN,2,1000,1,99,1,75!NGEN,3,1000,1001,1099,1,50!NGEN,7,1000,3001,3099,1,75!NGEN,4,1000,9001,9099,1,200/3!NGEN,7,1000,12001,12099,1,75!NGEN,2,1000,18001,18099,1,50!/VIEW,1,-1,-1,1!TYPE,2 $REAL,2 $MAT,3!*DO,II,11,16,1 $E,II,II+1 $*ENDDO!*DO,II,83,88,1 $E,II,II+1 $*ENDDO!*DO,II,11,74,9 $E,II,II+9 $*ENDDO!*DO,II,17,80,9 $E,II,II+9 $*ENDDO!EGEN,20,1000,1,28,1!*DO,II,83,18083,1000 $E,II,II+1000 $*ENDDO!*DO,II,89,18089,1000 $E,II,II+1000 $*ENDDO!TYPE,2 $REAL,1 $MAT,2*DO,II,11,18011,1000 $E,II,II+1000 $*ENDDO!*DO,II,17,18017,1000 $E,II,II+1000 $*ENDDO!/ESHAPE,1 $EP!BLC4,,,L/2,B,H $BLC4,75,,100,B,-40 $WPOFFS,,,H $BLC4,625,,200,B,40 $WPCSYS,-1! WPOFFS,75 $WPROTA,,,90 $VSBW,ALL!WPOFFS,,,100 $VSBW,ALL $WPOFFS,,,450 $VSBW,ALL!WPOFFS,,,200 $VSBW,ALL $WPOFFS,,,450 $VSBW,ALL!WPCSYS,-1 $ALLSEL!LSEL,S,LOC,Y,0 $LSEL,A,LOC,Y,150 $LSEL,R,LOC,X,0 $LESIZE,ALL,,,10!LSEL,S,LOC,Z,0 $LSEL,A,LOC,Z,300 $LSEL,R,LOC,X,0 $LESIZE,ALL,,,8!LSEL,S,LOC,Z,300 $LSEL,R,LOC,Y,0 $LSEL,R,LOC,X,0,75 $LESIZE,ALL,75!LSEL,S,LOC,Z,300 $LSEL,R,LOC,Y,0 $LSEL,R,LOC,X,75,175 $LESIZE,ALL,25!LSEL,S,LOC,Z,300 $LSEL,R,LOC,Y,0 $LSEL,R,LOC,X,175,625 $LESIZE,ALL,75/2!LSEL,S,LOC,Z,300 $LSEL,R,LOC,Y,0 $LSEL,R,LOC,X,625,825 $LESIZE,ALL,100/3! LSEL,S,LOC,Z,300 $LSEL,R,LOC,Y,0 $LSEL,R,LOC,X,825,1275 $LESIZE,ALL,75/2! LSEL,S,LOC,Z,300 $LSEL,R,LOC,Y,0 $LSEL,R,LOC,X,1275,1325 $LESIZE,ALL,50!LSEL,S,LOC,Z,340 $LSEL,R,LOC,X,625 $LESIZE,ALL,,,8!LSEL,S,LOC,Y,0 $LSEL,R,LOC,X,625 $LSEL,R,LOC,Z,300,340 $LESIZE,ALL,,,1!LSEL,S,LOC,Z,-40 $LSEL,R,LOC,X,75 $LESIZE,ALL,,,8!LSEL,S,LOC,Y,0 $LSEL,R,LOC,Z,-40 $LESIZE,ALL,,,4!LSEL,S,LOC,Y,0 $LSEL,R,LOC,X,75 $LSEL,R,LOC,Z,0,-40 $LESIZE,ALL,,,1!VSEL,S,LOC,Z,0,H $VATT,1,3,1MSHAPE,0,3D $MSHKEY,1VMESH,ALL $ALLSELVSEL,S,LOC,Z,-40,0 $VSEL,A,LOC,Z,H,H+40 $VATT,2,3,3MSHAPE,0,3D $MSHKEY,1 $VMESH,ALL/VIEW,1,-0.2,-1,1 $EPLOT $ALLSELNUMMRG,ALL $NUMCMP,ALL $EPLOT/SOLU!NSEL,S,LOC,Z,-40 $NSEL,R,LOC,X,L0 $D,ALL,UY,,,,,UZ $ALLSEL!ASEL,S,LOC,X,L/2 $DA,ALL,SYMM $ALLSEL!NSEL,S,LOC,Z,H+40 $NSEL,R,LOC,X,725!*GET,NODE1,NODE,,COUNT $F,ALL,FZ,-110000/NODE1 $ALLSEL!ANTYPE,STATIC $NLGEOM,ON $NSUBST,80 $OUTRES,ALL,ALL $AUTOTS,1 $LNSRCH,1! CNVTOL,F,,0.05,2,0.5 $ALLSEL!SOLVE $FINISH/POST1 $SET,LAST $PRRSOL,FZ!SET,LAST $PLDISP,1!ESEL,S,TYPE,,2 $ETABLE,SAXL,LS,1 $PLLS,SAXL,SAXL!ESEL,S,TYPE,,1 $/DEVICE,VECTOR,ON $PLCRACK!/POST26 $NSOL,2,NODE(L/2,B/2,0),U,Z!PROD,3,2,,,,,,-1 $PROD,4,1,,,LOAD,,,110!/AXLAB,X,MID-UZ(MM) $/AXLAB,Y,P(KN)!XVAR,3 $PLVAR,4。
- 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,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,,uylsel,alllsel,s,loc,z,3000 lsel,r,loc,y,0dl,all,,uylsel,allksel,s,loc,x,0 ksel,r,loc,y,0 dk,all,uxallsel,all!施加荷载qmz=0.3asel,s,loc,y,200 sfa,all,1,pres,qmz allsel,all nsubst,40 outres,all,all time,qmz*10 neqit,40solvefinish/post1etable,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,-1 tb,miso,1,,14tbpt,,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.00tbpt,,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,3d mshkey,1 vmesh,all allsel,allfinish/solu!施加约束lsel,s,loc,z,0 lsel,r,loc,y,0 dl,all,,uydl,all,,uzlsel,alllsel,s,loc,z,3000 lsel,r,loc,y,0 dl,all,,uylsel,allksel,s,loc,x,0 ksel,r,loc,y,0 dk,all,ux allsel,all!施加荷载qmz=0.3asel,s,loc,y,200 sfa,all,1,pres,qmz allsel,all outres,all,all time,qmz*10 nsubst,40 neqit,40solvefinish/post1pldisp,1 etable,zxyl,ls,1 plls,zxyl,zxyl,1 finish/post26nsol,2,33,u,y prod,3,1,,,,,,1/100 prod,4,2,,,,,,-1 xvar,4plvar,3Welcome To Download !!!欢迎您的下载，资料仅供参考！。