钢筋混凝土梁的弹塑性分析ansys命令流
ansys实例命令流-弹塑性分析命令流
/FILNAME,Elastic-Plasitc,1/TITLE, Elastic-Plasitc Analysis!前处理。
/PREP7!**定义梁单元189。
ET,1,BEAM189 !定义单元。
!**梁截面1。
SECTYPE, 1, BEAM, HREC, , 0 !定义梁截面。
SECOFFSET, CENTSECDATA,50,100,6,6,6,6,0,0,0,0 !定义梁截面完成。
!**定义材料。
MPTEMP,,,,,,,, !定义弹塑性材料模型。
MPTEMP,1,0MPDATA,EX,1,,2.05e5MPDATA,PRXY,1,,0.3TB,BISO,1,1,2,TBTEMP,0TBDATA,,150,18600,,,, !定义弹塑性材料模型。
!**建立几何模型。
K,1, , , ,K,2 ,900,K,3 ,,50LSTR, 1, 2!**网格划分。
FLST,5,1,4,ORDE,1 !定义网格密度。
FITEM,5,1CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,,_YLESIZE,_Y1, , ,50, , , , ,1 !定义网格密度完成。
CM,_Y,LINE !网格划分。
LSEL, , , , 1CM,_Y1,LINECMSEL,S,_YCMSEL,S,_Y1LATT,1, ,1, , 3, ,1CMSEL,S,_YCMDELE,_YCMDELE,_Y1LMESH, 1 !网格划分完成。
!施加载荷及求解。
FINISH/SOL!**施加约束。
FLST,2,1,3,ORDE,1 !施加约束。
FITEM,2,1/GODK,P51X, , , ,0,UX,UY,UZ,ROTX, , ,FLST,2,1,3,ORDE,1FITEM,2,2/GODK,P51X, , , ,0,UY,UZ,ROTX, , , , !施加约束完成。
!**加载。
FLST,2,50,2,ORDE,2FITEM,2,1FITEM,2,-50SFBEAM,P51X,1,PRES,100, , , , , ,LSWRITE,1, !定义载荷步1完成。
钢筋混凝土梁的弹塑性分析ansys命令流
!(1)工作环境设置/FILENAME,RC-BEAM !指定工作文件名/TITLE,ALAL YSIS OF A RC-BEAM !指定图形标题!(2)进行前处理器/prep7!(3)定义单元类型ET,1,LINK8 !定义钢筋单元ET,2,SOLID65 !定义混凝土单元ET,3,MESH200 !用于拉伸成体单元KEYOPT,3,1,6!(4)定义钢筋截面积r,1,28.3r,2,50.3r,3,314.1!(5)为solid65单元定义一个实参数组r,4,!(6)定义混凝土材料MP,EX,2,2.55E10MP,PRXY,2,0.3TB,CONC,2,1,9, !定义混凝土的破坏参数TBDA TA,,0.3,0.55,1.55E6,-1,,TBDA TA,,,,0.6!(7)定义钢筋材料模型及参数mp,ex,2,2e5 !纵向受拉钢筋材料mp,prxy,2,0.3tb,bkin,2,1,2,1tbdata,,350mp,ex,3,2e5 !横向箍筋,架立钢筋材料mp,prxy,3,0.25tb,bkin,3,1,2,1tbdata,,200!(8) 创建以及复制节点/pnum,node,1/pnum,elem,1n,1n,9,200fill,1,9ngen,11,9,1,9,1,,30ngen,11,99,1,99,1,,,-150/view,1,1,1,1!(9)建立箍筋单元type,1real,1mat,3!水平箍筋*do,i,11,16,1e,i,i+1e,i+(83-11),i+(83-11)+1*enddo!坚直箍筋*do,i,11,74,9e,i,i+9e,i+6,i+6+9*enddo!产生整个模型的箍筋egen,11,99,all!(10)建立架立筋以及纵筋单元!创建上部的架立钢筋单元*do,i,83,node(25,270,-1500+150),99 e,i,i+99e,i+6,i+6+99*enddo!纵向受拉钢筋单元的属性type,1real,3mat,2!创建纵筋单元*do,i,11,node(25,30,-1500+150),99 e,i,i+99e,i+3,i+3+99e,i+6,i+6+99*enddo/view,1,1,1,1/pnum,elem,0/pnum,node,0/eshape,1eplot!(11)建立混凝土剖面并划分网格k,1K,2,200,k,3,200,300k,4,,300a,1,2,3,4lsel,s,loc,y,0lsel,a,loc,y,300lesize,all,,,8lsel,alllsel,s,loc,x,0lsel,a,loc,x,200lesize,all,,,10type,3amesh,all!(12)拉伸形成混凝土单元type,2real,3mat,1extopt,esize,20extopt,aclear,1vext,all,,,,,-1500/pnum,mat,1/pnum,node,0/pnum,elem,0eplotallsel!(13)合并压缩节点编号nummrg,allnumcmp,alleplot!(14)施加支座约束nsel,s,loc,y,0nsel,r,loc,z,-1500+75d,all,uyd,all,ux!(15)施加对称面约束asel,s,loc,z,0da,all,symm!(16)施加载荷nsel,allnsel,s,loc,y,300nsel,r,loc,z,-450d,all,uy,-30allsel!(17)退出前处理器fini!(18) 进入求解器/solu!(19)求解器选项设置nlgeom,onnsubst,200outres,all,allneqit,50pred,oncnvtol,f,,0.05,2,0.5allsel!(20)求解并退出求解器solvefini!(21)进入通用处理器并读入最后一个子步的结果/post1set,last!(22)后处理操作plnsol,u,y,0,1/device,vector,1plcrack,0,0!(23)退出通用后处理器finish。
钢管混凝土ansys命令流
碳纤维加固桥梁结构的健康监测摘要:提出了一种针对老化桥梁的健康监测的可靠性方法。
首先, 根据桥梁的类型提出了临界失效标准。
然后基于这些标准定义安全界限。
假定大量的失效标准是正态分布的随机变量。
基于这样的假设,基本可靠性指标和基本失效概率为每个关键失效模式估计。
已经发现基本失效概率值,桥梁系统失效概率值按这段时间考虑计算。
最后,用系统失效概率来获得桥梁系统的可靠性指标。
系统的可靠性指数用来表达的桥梁在使用期间的状况较为片面。
选择斯里兰卡国家铁路桥网络的一座铁路桥梁作为一个案例研究来说明该可靠度程序。
对于这座桥,失效的最重要的形式是疲劳和腐蚀。
应用程序的建议方法表明,桥梁在现有负载的情况下的现状是令人满意的。
对可靠性的预测,不是为了抵御破坏,而是一种对基础设施管理进行有效的积极的健康监测方式。
关键字:碳纤维结构健康监测状态评估维护可靠性指标失效概率1、介绍在当今迅速变化的世界中,人类的生活质量和国家经济发展依赖于基础设施的数量、质量、和有效性]1[。
在世界各国,民用基础设施通常是需要维护的,康复,或更换的]2[。
特别是对基础设施构成一个相当大的投资的高速公路和相关的地面运输系统,直接影响到一个国家的工业产值。
目前大多数桥梁养护策略基于由桥梁检查员在不同时间间隔进行的视觉检查]3[。
由于人类检测取决于个人检查员,其结果存在一定程度的不确定性。
因此,绝对的依赖检验报告是不可靠的,而且可能导致错误的决策和桥梁拥有者的高维护成本。
基于桥梁的剩余使用寿命维护成本的优化。
另一方面,桥梁的检查员发现结构和非结构化缺陷,破坏和恶化过程的能力不应该被忽视。
事实上,人类经验应该被纳入维护决策,但要有足够的主观的知识。
只在那时,任何可接受的战略才能成为可能。
当前桥梁结构健康监测方法存在许多缺点和普遍的不确定性。
对结构健康监测和桥梁维护的知识不足造成由于频繁的维护活动,不可预见的桥梁损害,和偶然的桥梁失效引起的重大的经济损失。
ANSYS工程结构数值分析---命令流实例
!------------------------------------------------------!EX8.26 钢筋混凝土简支梁数值分析!分离式模型,关闭压碎,keyopt(1)=0,keyopt(7)=1 !力加载,位移收敛准则,误差1.5%,1/4模型分析!--------------------------------------------finish/clear/config,nres,2000/prep7!1.定义单元与材料性质--------------------et,1,solid65,,,,,,,1 !K1=0,k7=1et,2,link8mp,ex,1,13585mp,prxy,1,0.2fc=14.3ft=1.43tb,concr,1tbdata,,0.5,0.95,ft,-1tb,miso,1,,11tbpt,,0.0002,fc*0.19tbpt,,0.0004,fc*0.36tbpt,,0.0006,fc*0.51tbpt,,0.0008,fc*0.64tbpt,,0.001,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.0e5mp,prxy,2,0.3tb,biso,2tbdata,,300,0pi=acos(-1)r,1,0.25*pi*22*22r,2,0.25*pi*22*22/2r,3,0.25*pi*10*10r,4,0.25*pi*10*10/2!2.创建几何模型blc4,,,150/2,300,2000/2*do,i,1,9wpoff,,,100vsbw,all*enddo wpcsys,-1 wpoff,,,50 vsbw,all wpcsys,-1 wprota,,-90 wpoff,,,30 vsbw,all wpoff,,,240 vsbw,all wpcsys,-1 wpoff,30 wprota,,,90 vsbw,all wpcsys,-1!3.划分钢筋网格elemsiz=50 lsel,s,loc,x,30 lsel,r,loc,y,30 cm,zj,linelatt,2,1,2 lesize,all,elemsiz lsel,s,loc,x,75 lsel,r,loc,y,30 cm,zjb,line latt,2,2,2 lesize,all,elemsiz lsel,s,loc,x,30 lsel,r,loc,y,270 cm,jlj,linelatt,2,3,2 lesize,all,elemsiz lsel,s,tan1,z lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 lsel,u,loc,z,50 cm,gj,linelatt,2,3,2 lesize,all,elemsiz lsel,s,loc,z,0 lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 cm,gjb,line latt,2,4,2lesize,all,elemsizlsel,allcmsel,s,zjcmsel,a,zjbcmsel,a,jljcmsel,a,gjcmsel,a,gjbcm,gj,linelmesh,alllsel,all!4.划分混凝土网格vatt,1,,1mshkey,1esize,elemsizvmesh,allallsel,all!5.施加荷载和约束lsel,s,loc,y,0lsel,r,loc,z,900dl,all,,uyasel,s,loc,z,0da,all,symmasel,s,loc,x,75da,all,symm!-----------------p0=180000q0=p0/150/100asel,s,loc,z,0,50asel,r,loc,y,300sfa,all,1,pres,q0allsel,all!6.求解控制设置/soluantype,0nsubst,80outres,all,allautos,onneqit,50cnvtol,u,,0.015solve!7.进入POST1查看结果/post1set,lastesel,s,type,,2etable,saxl,ls,1plls,saxl,saxlesel,s,type,,1/device,vector,onplcrackeall!8.进入时程后处理/post26nsol,2,205,u,yprod,3,1,,,,,,p0/1000prod,4,2,,,,,,-1xvar,4plvar,3!EX2.6D 60kg/m钢轨实体建模FINISH/CLEAR/PREP7!1.定义参数H1=30.5H2=48.5H3=48.5B1=45.75B2=29.25B3=16.5B4=36.5B5=10XD1=1/3XD2=1/9XD3=1/3XD4=1/20R0=400R1=2R2=4R3=40R4=20R5=25R7=13R8=80R9=300R10=2FAI=43H=H1+2*H2+H3!2.创建轨底主要组成部分的线K,1K,2,B1+B2K,3,B1+B2,H1-B1*XD1-B2*XD2K,4,B1,H1-B1*XD1K,5,0,H1*DO,I,1,4L,I,I+1*ENDDO!3.计算确定轨头下三个圆弧的位置CTA=A TAN(XD3)BTA=(90-FAI)*ACOS(-1)/180-CTAA=1/TAN(CTA)C1=B3/2+R0-R6*SIN(CTA)-(R5-R6)*COS(BTA)C2=R5*SIN(BTA)+R6*COS(CTA)-R6*SIN(BTA)-H2C=(C1-C2/TAN(CTA))/(R0-R5)C3=A*C-SQRT(A*A-C*C+1)C3=C3/(A*A+1)REFA=ASIN(C3)C1=R6*SIN(CTA)+(R5-R6)*COS(BTA)+(R0-R5)*COS(REFA) !4.创建轨腰和轨头结构线ROUX=B3/2+R0-C1ROUY=ROUX*TAN(CTA)ANGI=ASIN((H1+H2)/R0)K,6,B3/2+R0-R0*COS(ANGI)K,7,B3/2+R0-R0*COS(REFA),H1+H2+R0*SIN(REFA) LARC,6,7,2,R0DX1=R6*(COS(BTA)-SIN(CTA))DY1=R6*(COS(CTA)-SIN(BTA))K,8,ROUX-DX1,H1+2*H2+ROUY-DY1K,9,ROUX,H1+2*H2+ROUYLARC,7,8,2,R5LARC,8,9,2,R6K,10,B4,H1+2*H2+B4*XD3K,11,B4-(H3-B4*XD3)*XD4,HL,9,10L,10,11ANGI=ASIN(B5/R9)K,12,0,HK,13,B5,H-(R9-R9*COS(ANGI)) LARC,12,13,1,R9DX1=B4-B5+R8*SIN(ANGI)ANGI1=ASIN(DX1/R8)Y1=R8*COS(ANGI)-R8*COS(ANGI1) K,14,B4,KY(13)-Y1LARC,13,14,1,R8!5.倒角或弧线连接LFILLT,1,2,R1LFILLT,2,3,R2LFILLT,3,4,R3LPTN,4,5LFILLT,17,18,R4LFILLT,8,9,R10LPTN,9,11LFILLT,20,22,R7LDELE,15,16,1,1LDELE,19,21,2,1LSYMM,X,ALLNUMMRG,ALLNUMCMP,ALLAL,ALLASUMVOFFST,1,300!EX2.9B 斜向圆台建模finish/clear/prep7r0=3r1=2h=4cyl4,,,r0wpoff,r1,,hcyl4,,,r1*do,i,1,4askin,i,i+4*enddova,all!斜向延伸创建圆台vgen,1,1,,,,3*r0,,,,1wpcsyscyl4,,,r0vext,7,,,r1,,h,r1/r0,r1/r0!EX2.9B 已知函数方程时的建模!玫瑰线1==================================== !极坐标方程为ρ=asin(4Φ/3)finish/clear/prep7a=20csys,1*do,i,0,360*3fei=i*acos(-1)/180r0=a*sin(fei*4/3)k,,r0,i*enddo*do,i,1,360*3l,i,i+1*enddo!玫瑰线2=================================== !极坐标方程为ρ=asin(4Φ)finish/clear/prep7a=20csys,1*do,i,0,360fei=i*acos(-1)/180r0=a*sin(fei*4)k,,r0,i*enddo*do,i,1,360l,i,i+1*enddo!长辐圆内旋轮线============================ !参数方程x=(a-b)cost+λcos(a-b)/bt!参数方程y=(a-b)sint-λsin(a-b)/btfinish/clear/prep7a=50b=5lmda=15*afun,deg*do,i,0,360x1=(a-b)*cos(i)+lmda*cos((a-b)/b*i)y1=(a-b)*sin(i)-lmda*sin((a-b)/b*i)k,,x1,y1*enddo*do,i,1,360l,i,i+1*enddo!椭圆抛物面================================ !方程为z=x*x/a/a+y*y/b/b!使用蒙皮方法finish/clear/prep7a=20b=10n=20*do,i,1,niz=i/10lsel,nonewpoff,,,izcyl4,,,1adele,alllsscale,all,,,a*sqrt(iz),b*sqrt(iz),,,,1*enddoallsel*do,i,1,4*(n-1),4askin,i,i+4askin,i+1,i+5askin,i+2,i+6askin,i+3,i+7*enddo!EX2.9C 旋转图素建模!类花瓶建模--------------------finish/clear/prep7K,1K,2,100K,3,300,400K,4,150,550k,5,130,880K,6,300,1150k,7,400,1100k,8,500,1150bsplin,allk,1001,,1000arotat,all,,,,,,1,1001,360,5!类轮建模----------------------finish/clear/prep7!先创建类轮结构的一个断面,然后旋转之blc4,,,4,6cyl4,2,6,1.5asba,1,2k,100,-3k,101,-3,10vrotat,all,,,,,,100,101!EX3.1 自定义箱形截面finish/clear/prep7k,1k,2,2k,3,2,2.2k,4,3,2.3k,5,3,2.5k,6,0,2.5a,1,2,3,4,5,6blc4,,0.2,1.7,2asba,1,2wpoff,1.7wprota,,,90asbw,allwpoff,,,0.3asbw,allwpoff,,0.2wprota,,90asbw,allwpoff,,,-2asbw,allarsym,x,allwpcsysaglue,allet,1,plane82lesize,all,,,1amesh,allsecwrite,myboxfinish/clear/prep7et,1,beam189mp,ex,1,3.0e10mp,prxy,1,0.167sectype,1,beam,mesh secread,mybox,,,meshk,1k,2,10k,100,5,5l,1,2lesize,all,,,20latt,1,,1,,100,,1lmesh,all/eshape,1Eplot!EX3.2 自定义多种材料截面finish/clear/prep7Ro=1.5Ri=1.0csys,1cyl4,,,ricyl4,,,roaptn,allwprota,,90asbwa,allwprota,,,90asbw,allwpcsyset,1,plane82mymat1=4mymat2=7mp,ex,mymat1,1.0mp,ex,mymat2,2.0asel,s,loc,x,0,riaatt,mymat1,,1asel,s,loc,x,ri,roaatt,mymat2,,1allselesize,0.25mshape,0,2dmshkey,1amesh,allsecwrite,mycsolid,sectfinish/clear/prep7et,1,beam189mym1=4mym2=7mp,ex,mym1,3.0e10mp,prxy,mym1,0.167 mp,ex,mym2,2.1e11mp,prxy,mym2,0.3 sectype,1,beam,mesh secread,mycsolid,sect,,mesh k,1k,2,,,10l,1,2lesize,all,,,20latt,,,1,,,,1lmesh,all/eshape,1/pnum,mat,1eplot!EX3.3A 单个方位关键点示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,ISECDATA,100,40,160,10,10,8K,1K,2,,,1000L,1,2K,100,,500,500LATT,1,,1,,100,,1LGEN,4,1,,,500LESIZE,ALL,,,10LMESH,ALL/ESHAPE,1EPLOT!EX3.3B 同时使用KB和KE方位关键点示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,ISECDATA,100,40,160,10,10,8L0=1000DL=500DXC=400K,1K,2,,,L0L,1,2K,100,,DLK,200,DXC,-DLK,300,2*DXC,DLK,301,2*DXC+DLK,400K,500,8*DXCLGEN,5,1,,,DXCLATT,1,,1,,100,,1LSEL,S,,,2LATT,1,,1,,200,,1LSEL,S,,,3LATT,1,,1,,300,301,1LSEL,S,,,4LATT,1,,1,,400,,1LSEL,S,,,5LATT,1,,1,,500,,1LSEL,ALLLESIZE,ALL,,,50LMESH,ALL/ESHAPE,1EPLOT!EX3.3C 双方位关键点翻转示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,I SECDATA,100,40,160,10,10,8L0=1000DL=500DXC=400K,1K,2,,,L0K,3,DXCK,4,DXC,,L0/2K,5,DXC,,L0L,1,2L,3,4L,4,5K,100,,DLk,101,,-DLK,200,DXC,DLk,201,2*DXCK,202,DXC,-DLLATT,1,,1,,100,101,1lesize,all,,,100LSEL,S,,,2LATT,1,,1,,200,201,1lesize,all,,,50LSEL,S,,,3LATT,1,,1,,201,202,1lesize,all,,,50LSEL,ALLLMESH,ALL/ESHAPE,1EPLOT!EX3.4A 两种单元形状和两种网格划分比较finish/clear/prep7et,1,plane82k,1k,2,8k,3,7,6k,4,1,6a,1,2,3,4esize,1mshape,0mshkey,0!mshape,0!mshkey,1!mshape,1!mshkey,0!mshape,1!mshkey,1amesh,all!EX3.4B 中间节点位置控制网格划分比较finish/clear/prep7cyl4,,,4,,8,60lesize,all,,,2mshape,0mshkey,1mshmid,0!mshmid,1!mshmid,2amesh,all!EX3.5 线上单元尺寸设置示例!下边密上边稀finish/clear/prep7et,1,plane82blc4,,,10,10lsel,s,tan1,ylesize,all,,,10lsel,s,loc,x,0lesize,all,,,9,1/8lsel,s,loc,x,10lesize,all,,,9,8lsel,allmshape,0mshkey,1amesh,all!中间密外边稀finish/clear/prep7et,1,plane82blc4,,,10,10lsel,s,tan1,ylesize,all,,,10,-1/5lsel,s,tan1,xlesize,all,,,9,-1/8lsel,allmshape,0mshkey,1amesh,all!EX3.6 合并线和连接线以进行映射网格划分finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,3mshape,0mshkey,1lccat,1,2lcomb,4,5amesh,all!EX3.7 合并线或连接线的网格划分设置finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,,10lesize,6,,,8lesize,1,,,4lesize,2,,,3lesize,4,,,5lesize,5,,,2mshape,0mshkey,1lcomb,1,2!lccat,1,2!如采用lcatt,1,2则采用7个划分数lcomb,4,5amesh,all!EX3.7B 简化映射网格划分finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,,6mshape,0amap,1,2,5,3,4!EX3.8A 过渡四边形映射网格finish/clear/prep7k,1k,2,10,-1k,3,8,6k,4,1,3a,1,2,3,4lesize,1,,,8lesize,3,,,3lesize,4,,,7lesize,2,,,2mshape,0,2dmshkey,1amesh,all!EX3.8B 过渡四边形映射网格finish/clear/prep7et,1,plane42k,1k,2,10,-1k,3,8,6k,4,1,3a,1,2,3,4lesize,1,,,11lesize,3,,,3lesize,4,,,2lesize,2,,,2mshape,0,2dmshkey,1amesh,all!EX3.9A 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,all,,,3lesize,5,,,15allselmshape,0,3dvmesh,all!EX3.9B 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,all,,,4lesize,7,,,12mshape,0,3dmshkey,1vmesh,all!EX3.9C 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lsel,s,,,1,8,7lesize,all,,,6lsel,s,,,2,7,5lesize,all,,,3lsel,s,,,4,5lesize,all,,,7lsel,s,,,3,6,3lesize,all,,,2lsel,s,,,9,12,1lesize,all,,,5allselmshape,0,3dmshkey,1vmesh,all!EX3.9D 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,2,,,2lesize,4,,,2lesize,5,,,2lesize,1,,,3lesize,3,,,3lesize,7,,,6lesize,6,,,7lesize,11,,,7lesize,9,,,3lesize,10,,,3lesize,12,,,3allselmshape,0,3dmshkey,1vmesh,all!EX3.10 带两孔的长方体的扫掠网格划分finish/clear/prep7a=10r=2et,1,mesh200,6et,2,solid45blc4,,,2*a,a,acyl4,a/2,a/2,r,,,,awprota,,90cyl4,1.5*a,a/2,r,,,,-avsbv,1,2vsbv,4,3wprota,,,90wpoff,,,a/2vsbw,allwpoff,,,a/2vsbw,allwpoff,,,a/2vsbw,allwpcsyswpoff,,a/2,a/2vsbw,allwprota,,90vsbw,allwpcsysesize,1amap,105,15,16,26,63amap,107,16,13,60,26amap,108,13,14,28,60amap,103,15,14,28,63lesize,94,,,4lesize,79,,,5vsel,s,loc,x,0,avsweep,allasel,s,loc,y,aasel,r,loc,x,a,2*alccat,2,45lccat,57,71lccat,68,78lccat,65,72mshape,0,2dmshkey,1amesh,allvsel,s,loc,x,a,2*avsweep,allallsel/view,1,1,2,3Eplot!EX3.11 拉伸类命令生成体单元网格finish/clear/prep7et,1,82et,2,95blc4,,,4,4blc4,6,,4,4blc4,12,,4,4blc4,18,,4,4esize,1amesh,allesize,,8vrotat,1,,,,,,1,4,90vext,2,,,,,10,0.5,0.5voffst,3,10vdrag,4,,,,,,35!EX3.12 圆的网格划分finish/prep7et,1,plane82r0=10cyl4,,,r0cyl4,3*r0,,,,r0 wprota,,90asbw,allwprota,,,90asbw,allwpoff,,,3*r0asbw,allwpcsys,-1asel,s,loc,x,-r0,r0lsla,slesize,all,,,8mshape,0,2dmshkey,1amesh,allallselasel,s,loc,x,2*r0,4*r0 lesize,all,,,8lsel,r,length,,r0lesize,all,,,8,0.1,1 amesh,allallsel!EX3.13 圆环的网格划分finish/clear/prep7et,1,plane82r0=10cyl4,,,r0/3,,r0,90cyl4,2*r0,,r0/10,,r0,90 asel,s,loc,x,-r0,r0lsla,slesize,all,,,8lsel,r,length,,r0*2/3 lesize,all,,,3,,1 mshape,0,2dmshkey,1amesh,allasel,s,loc,x,2*r0,4*r0lesize,5,,,12lesize,7,,,6lsel,s,length,,r0*9/10lesize,all,,,7amesh,all!EX3.14 圆柱面和圆柱体的网格划分!圆柱面finish/clear/prep7r0=10h0=50et,1,shell63cyl4,,,r0adele,1cm,l1cm,linek,50k,51,,,h0l,50,51adrag,l1cm,,,,,,5lsel,s,loc,z,0lesize,all,,,6lsel,s,length,,h0lesize,all,,,8mshape,0,2dmshkey,1amesh,all!圆柱体finish/clear/prep7r0=10h0=50et,1,solid95cyl4,,,r0,,,,h0wprota,,90vsbw,allwprota,,,90vsbw,allmshape,0,3dmshkey,1lsel,s,loc,z,0lsel,s,length,,h0lesize,all,,,8vmesh,all!EX3.15 圆锥的网格划分!扫掠网格划分finish/clear/prep7cone,10,,,15,,90et,1,200,7et,2,95lsel,s,,,5,6lesize,all,,,12,0.5lsel,alllesize,3,,,6mshape,0,2dmshkey,1amesh,3vsweep,1,3,4vsymm,x,allvsymm,y,allvglue,all!六面体映射网格划分finish/clear/prep7cone,10,,,15,,90et,2,95lesize,all,,,8mshape,0,3dmshkey,1vmesh,allvsymm,x,allvsymm,y,allvglue,all!EX3.16 正多边形面的通用网格划分finish/clear/prep7ns=11rpr4,ns,,,10kp0=100+nsk,kp0*if,mod(ns,2),eq,0,then*do,i,1,ns/2l,kp0,2*i-1*enddo*else*do,i,1,nsl,kp0,i*enddo*endiflsel,s,,,ns+1,2*nscm,l1cm,lineallselasbl,1,l1cm*if,mod(ns,2),eq,0,then cmsel,s,l1cmlesize,all,,,8lsel,invelesize,all,,,4*elselesize,all,,,8*endifallselmshape,0mshkey,1amesh,all!EX3.17 球体及球面网格划分finish/clear/prep7r0=10sphere,,r0,,90vsbw,allvdele,2,,,1esize,2et,1,solid95mshape,0,3dmshkey,1vsymm,x,all vsymm,y,all vsymm,z,allvglue,allfinish/clear/prep7r0=10sphere,,r0,,90 vsbw,allvdele,2,,,1vdele,allasel,s,loc,x,0asel,a,loc,y,0asel,a,loc,z,0adele,all,,,1allselesize,2et,1,shell63 mshape,0,2d mshkey,1amesh,allarsym,x,allarsym,y,allarsym,z,all nummrg,all!EX3.18 孔板网格划分finish/clear/prep7a0=300b0=800r0=15blc4,,,a0,b0cyl4,a0/4,b0/8,r0 agen,2,2,,,a0/2 agen,2,2,3,1,,b0/8 agen,2,2,5,1,,b0*5/8 asel,s,,,2,9,1cm,a2cm,areaallselwprota,,-90*do,i,1,5wpoff,,,b0/16asbw,all*enddowpoff,,,b0*5/16*do,i,1,5wpoff,,,b0/16asbw,all*enddowprota,,,90*do,i,1,3wpoff,,,a0/4asbw,all*enddowpcsys,-1numcmp,alllsel,s,radius,,r0lesize,all,,,8lsel,invelesize,all,,,4lsel,allet,1,82mshape,0,2dmshkey,1asel,u,loc,y,b0/16,b0*5/16 asel,u,loc,y,b0*11/16,b0*15/16 lsla,slsel,r,tan1,xlesize,all,50,,,,1amesh,allallselamap,21,30,31,54,62 amap,22,31,32,54,64 amap,32,29,30,62,76 amap,33,29,32,64,76amap,42,34,35,76,82 amap,43,35,36,76,84 amap,41,33,34,53,82 amap,44,33,36,53,84amap,19,22,23,49,64 amap,20,23,24,49,60 amap,9,21,22,64,74 amap,31,21,24,60,74amap,38,26,27,74,84 amap,39,27,28,74,80 amap,37,25,26,50,84 amap,40,25,28,50,80amap,17,14,15,43,59 amap,18,15,16,43,63 amap,7,13,14,59,73 amap,8,13,16,63,73amap,34,18,19,73,79 amap,30,17,18,79,44 amap,35,19,20,73,83 amap,36,17,20,44,83amap,15,6,7,39,63amap,16,7,8,39,58amap,4,6,5,63,71amap,5,5,8,58,71amap,27,10,11,71,83 amap,28,11,12,71,78 amap,26,9,10,83,40 amap,29,9,12,40,78!EX3.19 角支架的网格划分finish/clear/prep7blc4,,,150,50blc4,100,,50,-50cyl4,,25,25cyl4,125,-50,25aadd,allnumcmp,allcyl4,,25,10cyl4,125,-50,10 asel,s,,,2,3cm,a1cm,area asel,allasba,1,a1cm lcomb,1,6lfillt,1,2,20 asbl,4,6adele,1,,,1lfillt,3,4,20al,18,19,20 aadd,all numcmp,all!---------------------- wprota,,90 wpoff,,,-25 asbw,all wpoff,,,75 asbw,all wprota,,,90 asbw,all wpoff,,,125 asbw,all wpcsys,-1 wpoff,25 wprota,,,90 asbw,all kwpave,18 asbw,all kwpave,3 wprota,,90 asbw,all kwpave,21 wprota,,-45 asbw,8 wpcsys,-1!-------------------- et,1,plane82 mshape,0,2d mshkey,1 esize,6lesize,33,,,6 lesize,37,,,6 lesize,42,,,6amesh,3,5,2amap,6,9,10,4,23amap,7,9,12,1,23amesh,11,13,2amesh,1,9,8amesh,2,4,2amap,16,3,18,26,28amap,12,14,15,5,28amap,15,13,14,6,28lcomb,19,22lcomb,27,46amesh,10,14,4!========================= !EX3.20 大板小孔的网格划分A finish/clear/prep7a0=100blc4,,,a0,a0cyl4,,,a0/100asba,1,2csys,1k,50,a0/20k,51,a0/20,90l,50,51asbl,all,1wprota,,90,45asbw,allet,1,plane82mshape,0,2dmshkey,1lsel,s,loc,x,0,a0/20lesize,all,,,8lsel,alllesize,7,,,12,0.1lesize,8,,,12,0.1lesize,9,,,12,10amesh,allcsys,0arsym,x,allarsym,y,allnummrg,all!============================= !EX3.20 大板小孔的网格划分Bfinish/clear/prep7a0=100blc4,,,a0,a0cyl4,,,a0/100asba,1,2csys,1k,50,a0/20k,51,a0/20,90l,50,51asbl,all,1et,1,plane82mshape,0,2dmshkey,1lesize,5,,,8lesize,1,,,10lesize,4,,,8,5lesize,6,,,8,5amesh,1lesize,7,,,20,0.1lesize,8,,,20,0.1amap,2,50,51,2,4csys,0arsym,x,allarsym,y,allnummrg,all!EX3.21 Z形面网格划分及局部细分finish/clear/prep7blc4,,,15,10blc4,10,6,14,12aadd,allwprota,,-90wpoff,,,6asbw,allwpoff,,,4wprota,,,90wpoff,,,10asbw,allwpoff,,,5asbw,allwpcsys,-1et,1,82esize,2mshape,0,2dmshkey,1amesh,allkrefine,9,10,1,1,,off !不光滑处理(节点位置不变)aclear,allamesh,allkrefine,9,10,1,1,,smooth !光滑处理,节点位置可改变,但既有单元不删除aclear,allamesh,allkrefine,9,10,1,1 !光滑处理,可删除既有单元重细分,节点位置改变!EX3.22 分布小孔环形面的网格划分finish/clear/prep7r1=15r2=25r3=3n=8cyl4,,,r1,,r2,180/ncyl4,0.5*(r1+r2),,r3asba,1,2kl,2larc,6,8,4,0.5*(r1+r2)asbl,all,4lesize,4,,,4lesize,10,,,4lesize,6,,,8lesize,3,,,6lesize,8,,,4lesize,5,,,8lesize,9,,,4lesize,7,,,4et,1,plane82mshape,0,2damap,2,6,7,4,8amap,1,5,6,8,1!lrefine,5,6,1,1,smootharsym,y,allcsys,1agen,n,all,,,,360/nnummrg,all!EX3.23 同心圆轴体的网格划分finish/clear/prep7cylind,50,0,100,150cylind,40,0,200,250cylind,30,0,100,250cylind,20,0,50,250cylind,10,0,0,250vptn,allwprota,,,90vsbw,allwprota,,90vsbw,allet,1,95csys,1lsel,s,length,,50lesize,all,,,5lsel,invelsel,u,radius,,10,50lesize,all,,,4n1=6n2=6n3=8n4=10n5=12lsel,r,loc,x,0,10lesize,all,,,n1,,1lsel,s,radius,,10lesize,all,,,n1lsel,s,radius,,20lsel,s,radius,,30lesize,all,,,n3lsel,s,radius,,40lesize,all,,,n4lsel,s,radius,,50lesize,all,,,n5allselmshape,0,3dmshkey,1vmesh,all!EX3.24A 圆柱与长方体组合模型的网格划分!采用六面体映射网格划分finish/clear/prep7a0=30h1=15h2=25r0=7blc4,,,a0/2,a0/2,h1cyl4,,,r0,,,90,h1+h2vptn,allaccat,4,6et,1,95esize,3mshape,0,3dmshkey,1vmesh,allasel,s,accaadele,alllsel,s,lccaldele,allallselvsymm,x,allvsymm,y,allnummrg,all!EX3.24B 圆柱与长方体组合模型的网格划分!采用过渡六面体映射网格划分finish/cleara0=30h1=15h2=25r0=7blc4,,,a0/2,a0/2,h1cyl4,,,r0,,,90,h1+h2vptn,alllsel,s,radius,,r0lsel,a,length,,r0lesize,all,,,6lsel,s,loc,z,0lsel,a,loc,z,h1lesize,all,,,7asel,s,loc,x,a0/2asel,a,loc,y,a0/2accat,allallselet,1,95esize,3mshape,0,3dmshkey,1vmesh,allasel,s,accaadele,alllsel,s,lccaldele,allallselvsymm,x,allvsymm,y,allnummrg,all!EX3.25-任意圈数圆柱形螺旋弹簧的扫略网格划分FINISH/CLEAR/PREP7!1.定义弹簧参数(同上)D=4C=8DZ=C*DT=DZ/6N=4.7*IF,T,LT,D,THEN*ENDIFTKPD=90!2.创建全部螺旋线CSYS,1TDEG=N*360TDEG1=MOD(TDEG,TKPD) N0=(TDEG-TDEG1)/TKPD+1 *AFUN,DEG*DO,I,1,N0CTA=(I-1)*TKPDZ=T/360*CTAK,I,DZ/2,CTA,Z*ENDDO*IF,TDEG1,LT,1.0E-2,THEN *ELSEN0=N0+1CTA=CTA+TDEG1Z=T/360*CTAK,N0,DZ/2,CTA,Z*ENDIF*DO,I,1,N0-1L,I,I+1*ENDDOCM,L1,LINE!3.在螺旋线端部创建簧丝截面CSYS,0WPOFF,DZ/2 WPROTA,,90CYL4,,,D/2wprota,,90asbw,allwprota,,,90asbw,allcm,a1cm,areaVDRAG,a1cm,,,,,,L1et,1,mesh200,7et,2,solid95cmsel,s,a1cmlsla,slesize,all,,,6amesh,allallselesize,d/2vsweep,all!EX3.26 具坑缺陷圆柱的网格划分finish/clear/prep7l0=50r0=10r1=5r2=1h0=5et,1,solid45cyl4,,,r0,,,,l0wpave,0,r0,l0/2wprota,,90con4,,,r1,r2,h0vsbv,1,2wprota,,90vsbw,allwpoff,,,r1+3vsbw,allwpoff,,,-2*(r1+3)vsbw,allwprota,,,90vsbw,allwpoff,,,r1+1vsbw,allwpoff,,,-2*(r1+1)vsbw,allwpcsys,-1wprota,,90vsbw,allwpoff,,,-r1vsbw,allwpoff,,,2*r1vsbw,allwpcsys,-1numcmp,allksel,s,loc,z,l0/2ksel,r,loc,x,0ksel,u,loc,y,0*get,kp1,kp,0,num,min kp2=kpnext(kp1) ksel,s,loc,z,l0/2+r2*get,kp3,kp,0,num,min allselk,1000,,-r0,l0/2+5l,kp3,1000asel,none*get,l1,line,0,num,max arotat,l1,,,,,,kp1,kp2 cm,a1cm,areaallselvsba,all,a1cm!vsel,s,loc,z,17,33!vsel,r,loc,x,-6,6!aslv,s!aplotaccat,128,61accat,160,138 accat,224,202 accat,183,35accat,31,134accat,151,164 accat,215,231 accat,82,195accat,130,5accat,162,145 accat,229,209 accat,192,89accat,12,123accat,117,155 accat,109,219 accat,14,178mshape,0,3d mshkey,1esize,1vmesh,all!EX3.27长方体开圆柱槽的网格划分finish/clear/prep7blc5,,,20,10,30wpoff,,5,5cyl4,,,8,,,,10wpoff,,,10cyl4,,,4,,,,12vsel,s,,,2,3cm,v1cm,voluallselwpcsys,-1vsbv,1,v1cmvoffst,15,5voffst,16,3voffst,17,15vptn,allwprota,,,90vsbw,allwpcsys,-1wpoff,,,5vsbw,allwpoff,,,10vsbw,allwpoff,,,12vsbw,allet,1,mesh200,7esize,2mshape,0,2dmshkey,1asel,s,loc,z,30lccat,6,74lccat,8,68amesh,allasel,s,loc,z,0lccat,3,73lccat,1,67amesh,allallselet,2,solid95mshape,0,3dmshkey,1vsweep,all!ex4.1累加节点自由度约束值finish/clear/prep7et,1,beam3k,1k,2,10l,1,2esize,1lmesh,alld,1,ux,1e-2,,,,uyd,1,rotzd,2,alldlistnsel,s,d,u,0,0.1dofsel,s,uydscale,2.5dlistdofsel,alldcum,addd,1,ux,2e-2dlistdcum,ignod,2,uy,1.0dlist!ex4.2 对线施加约束并转换finish/clear/prep7et,1,95blc4,,,10,10,10dl,7,,ux,0.1dl,5,,alldl,11,6,symmdl,10,6,asymdl,6,,symmdllistesize,2vmesh,alldtrandlist!ex4.3 施加集中荷载与节点坐标系finish/clear/prep7et,1,beam4k,1k,2,5k,3,10l,1,2l,2,3local,12,0,,,,90nrotat,alldk,1,allfk,2,fy,-1000esize,1lmesh,allnrotat,alllplotfk,3,fy,1000f,6,fx,-1000sbctraneplot!ex4.4a 3d单元sf加载示例finish/clear/prep7et,1,95blc4,,,10,10,20esize,,4vmesh,allasel,s,loc,y,10sf,all,pres,1000asel,s,loc,z,20nsla,s,1sf,all,pres,1000!ex4.4b 2d单元sf加载示例finish/clear/prep7et,1,82blc4,,,100,200blc4,30,60,40,80asba,1,2wprota,,-90wpoff,,,60asbw,allwpoff,,,80asbw,allwprota,,,90wpoff,,,30asbw,allwpoff,,,40asbw,allwpcsys,-1esize,5amesh,all/psf,pres,norm,2sf,all,pres,100sfdele,all,presnsel,s,loc,x,0sf,all,pres,100nsel,s,loc,x,15,20esln,s,1nsel,r,loc,x,15sf,all,pres,110nsel,s,loc,x,40,60nsel,r,loc,y,10,30esln,s,1sf,all,pres,100lsel,s,loc,x,100nsll,s,1esln,snsel,s,loc,x,95sf,all,pres,-100eplot!ex4.5 节点号及其荷载函数finish/clear/prep7et,1,45blc4,,,10,10,20esize,5vmesh,all*dim,mypres,,100*do,i,1,100mypres(i)=i*10.0*enddosffun,pres,mypres(1)nsel,s,loc,y,10sf,all,pres,10sflist*do,i,1,100mypres(i)=i*50.0*enddonsel,s,loc,z,20sf,all,pres,0allselsflist!ex4.6a 2d平面单元plane82 finish/clear/prep7et,1,82blc4,,,10,50esize,2amesh,allsfe,1,4,pres,,100,50sfe,6,4,pres,,200,40nsel,s,loc,x,10esln,ssfe,all,2,pres,,-100/psf,pres,norm,2,0,1eplot!ex4.6b 3d体单元solid95 finish/clear/prep7et,1,95blc4,,,10,10,30esize,5vmesh,all/psf,pres,norm,2,0,1sfe,22,3,pres,,100sfe,22,6,pres,,100,110,120,130 sfelist!ex4.6c 3d壳单元shell63 finish/clear/prep7et,1,63wprota,,90blc4,,,10,10esize,5amesh,all/psf,pres,norm,2,0,1sfe,1,1,pres,,100sfe,1,3,pres,,100sfe,1,6,pres,,100sfe,4,1,pres,,100,110,120,130!ex4.7 在梁单元上施加荷载finish/clear/prep7et,1,beam3k,1k,2,10l,1,2esize,,10lmesh,all/pnum,elem,1sfbeam,3,1,pres,50,100 sfbeam,5,1,pres,100 sfbeam,7,1,pres,50,100,,,0.2,0.1sfbeam,9,1,pres,100,,,,0.4,-1sfbeam,3,2,pres,50,100,,,0.2,0.1!ex4.8a 在线上施加面荷载finish/clear/prep7et,1,82blc4,,,10,30esize,5amesh,all/psf,pres,norm,2sfl,4,pres,10,60sfl,2,pres,60sftraneplot!ex4.8b 在线上施加面荷载finish/clear/prep7et,1,shell63wprota,,90blc4,,,10,30esize,5amesh,allsfl,3,pres,100sftraneplot!ex4.9 利用荷载梯度在直角坐标系下的施加方法finish/clear/prep7et,1,82blc4,,,10,60esize,2amesh,all/psf,pres,norm,2sfgrad,pres,,y,0,-5nsel,s,loc,x,0nsel,r,loc,y,0,40sf,all,pres,600sfgrad,pres,,y,30,-20nsel,s,loc,x,10sf,all,pres,0allseleplot!ex4.10利用荷载梯度在柱坐标系下的施加方法finish/clear/prep7csys,1k,1,10,-90k,2,10,90k,3,10,90,30l,1,2l,2,3adrag,1,,,,,,2ldele,2,,,1numcmp,allet,1,63esize,2amesh,allsfgrad,pres,1,y,-90,1sfa,all,2,pres,400sftran!可考察下述命令及其施加方法!①!sfgrad,pres,1,y,270,1sfa,all,2,pres,400sftran!②!local,12,1cscir,12,1sfgrad,pres,12,y,270,1sfa,all,2,pres,400sftran!ex4.11 ldread命令的示例finish/clear/filname,ldtest1/prep7et,1,plane82mp,ex,1,2.1e5mp,prxy,1,0.3blc4,,,10,40esize,2amesh,alld,all,allsfl,3,pres,100/solusolve/post1prrsolfinish/filename,ldtest2/solulsclear,allnsel,s,loc,y,0d,all,allallselldread,reac,,,,,ldtest1,rst solve/post1plnsol,s,y!ex4.12 初应力荷载finish/clear/filname,colu1/prep7et,1,plane82mp,ex,1,2e5mp,nuxy,1,0.3blc4,,,1,10esize,2amesh,allnsel,s,loc,y,0d,all,uyd,1,uxnsel,s,loc,y,10sf,all,pres,-10。
ansys混凝土命令流
vsel,s,loc,x,-0.002,0 !选择左边建筑结构胶
vatt,2,2,1
esize,0.01
mshape,0,3d$mshkey,1
vmesh,all
vsel,s,loc,x,0.102,0.112 !选择右边钢板
vatt,3,3,1
wpoff,0.002,, !工作平面平移
blc4,,,0.01,0.05 ,0.16 !钢板
wpcsys,-1
wpoff,-0.002,0.025,0.05 !工作平面平移
blc4,,,0.002,0.05 ,0.15 !建筑结构胶
blc4,,,-0.01,0.05 ,0.16 !钢板
fini
/cle
/FILNAME, hntd, 1
/TITLE, tiegang
/prep7
!----------定义单元及材料等---------------------
Et,1,solid45 !定义混凝土单元类型
Et,2,solid45 !定义建筑结构胶单元类型
mp,prxy,1,0.2 !定义混凝土泊松比
mp,prxy,2,0.317 !定义建筑结构胶泊松比
mp,prxy,3,0.3 !定义钢板泊松比
mp,dens,1,2550 !定义混凝土密度
mp,dens,2,1700 !定义建筑结构胶密度
sfa,all,,pres,-1200e3/(2*6)
save
!-----------求解---------------------------
/solu
/nerr,10000,100000
solve
ANSYS应用实例:钢筋混凝土简支梁数值模拟-推荐下载
对全部高中资料试卷电气设备,在安装过程中以及安装结束后进行高中资料试卷调整试验;通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术通关,1系电过,力管根保线据护敷生高设产中技工资术艺料0不高试仅中卷可资配以料置解试技决卷术吊要是顶求指层,机配对组置电在不气进规设行范备继高进电中行保资空护料载高试与中卷带资问负料题荷试2下卷2,高总而中体且资配可料置保试时障卷,各调需类控要管试在路验最习;大题对限到设度位备内。进来在行确管调保路整机敷使组设其高过在中程正资1常料中工试,况卷要下安加与全强过,看度并22工且22作尽22下可22都能22可地护以缩1关正小于常故管工障路作高高;中中对资资于料料继试试电卷卷保破连护坏接进范管行围口整,处核或理对者高定对中值某资,些料审异试核常卷与高弯校中扁对资度图料固纸试定,卷盒编工位写况置复进.杂行保设自护备动层与处防装理腐置,跨高尤接中其地资要线料避弯试免曲卷错半调误径试高标方中高案资等,料,编试要5写、卷求重电保技要气护术设设装交备备置底4高调、动。中试电作管资高气,线料中课并敷3试资件且、设卷料中拒管技试试调绝路术验卷试动敷中方技作设包案术,技含以来术线及避槽系免、统不管启必架动要等方高多案中项;资方对料式整试,套卷为启突解动然决过停高程机中中。语高因文中此电资,气料电课试力件卷高中电中管气资壁设料薄备试、进卷接行保口调护不试装严工置等作调问并试题且技,进术合行,理过要利关求用运电管行力线高保敷中护设资装技料置术试做。卷到线技准缆术确敷指灵设导活原。。则对对:于于在调差分试动线过保盒程护处中装,高置当中高不资中同料资电试料压卷试回技卷路术调交问试叉题技时,术,作是应为指采调发用试电金人机属员一隔,变板需压进要器行在组隔事在开前发处掌生理握内;图部同纸故一资障线料时槽、,内设需,备要强制进电造行回厂外路家部须出电同具源时高高切中中断资资习料料题试试电卷卷源试切,验除线报从缆告而敷与采设相用完关高毕技中,术资要资料进料试行,卷检并主查且要和了保检解护测现装处场置理设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况,然后根据规范与规程规定,制定设备调试高中资料试卷方案。
混凝土梁的ANSYS命令流
! 水***** * *****定义支座垫板的几何特性及材料特性****** * *
E , O I4 T 3 S LD 5 , M, 33O 1 PE, . 0 X , E
SV AE
!***** 建立混凝土梁单元**** ****** *** *
附录 2
B O K 0 L 0 2 , . LC , 0 . 0 0 2 , , 5 , 0 BO K 0 L0 2 , 50 0 2 LC , , 5 0 3 , . , . . , 0 BO KO L 2 , 5- . 0 LC , 0 50 3 ,0 1, , . . 5 BO K 0 L0 2 , 5 02 , 5 LC , , 50 3 , 0 0 3 , . . . . V E , . S LA I L VL EA I G U, L
SV AE
左翼缘构造筋
L 部构造钢筋划分网格 [E ,, E,6 1, SlS LN, 1, 1 I 6 LE, IE, 71, S LA LN ,1, 1 , 7 LT , , AT 2 3 2 , LS Z , , 5 E IE A L 00 L .
LE IA L M SI L , l lAL 名E , L
K 2 , , 5 0 1, , 3 5 0 0 , 5 4 . . . L T , SR 2 , 1 2 2 L T , SR 2 , 3 2 4
!*****下面纵向钢筋划分网格* ** * ****** * **
LE , , E , 3 1, S L S LN , 1, 2 I 5
R1 ,,
! 定义混凝土的常数
M, , e0弹性模量要根据其他公式求 PE, 3 l! X12 4 .
M , Y 10 2 P P X , . R , 0 M , S l2E 密度 PDN , 53! E ,
钢筋混凝土结构Ansys命令流
钢筋混凝土结构Ansys命令流整体式finish/clear/config,nres,2000/prep7!定义单元与材料性质et,1,solid65,,,,,,,1mp,ex,1,13585mp,prxy,1,0.2mp,ex,2,2.0e5mp,prxy,2,0.3fc=14.3ft=1.43tb,concr,1tbdata,,0.5,0.95,ft,fctb,miso,1,,11tbpt,,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.85tb,biso,2tbdata,,300,0pi=acos(-1)v1=0.25*pi*22*22*1.5/(75*60) v2=0.25*pi*10*10/(75*60)v3=0.25*pi*10*10*9.5/(1000*75) v4=0.25*pi*10*10*9.5/(1000*60) r,1,2,v1,0,90,2,v4,0,0,2,v3,90,0 r,2,2,v2,0,90,2,v4,0,0,2,v3,90,0 r,3,2,v3,90,0!几何建模blc4,,,150/2,300,2000/2 wprota,,-90wpoff,,,60vsbw,allwpoff,,,180vsbw,allwpcsys,-1wpoff,,,50vsbw,allwpoff,,,850vsbw,allwpcsys,-1!定义属性划分单元网格vsel,s,loc,y,0,60vatt,1,1,1vsel,s,loc,y,240,300vatt,1,2,1vsel,s,loc,y,60,240vatt,1,3,1esize,50mshkey,1vmesh,all!device,vector,on/eshape,1 !此两个命令查看钢筋设置!施加荷载和约束lsel,s,loc,y,0lsel,r,loc,z,900dl,all,,uyasel,s,loc,z,0da,all,symmasel,s,loc,x,75da,all,symmp0=180000q0=p0/150/100asel,s,loc,z,0,50asel,r,loc,y,300sfa,all,1,pres,q0allsel,all!求解控制/soluantype,0nsubst,50outres,all,allautots,onneqit,50cnvtol,u,,0.015solve!进入post1查看结果/post1pldisp,1etable,reb1,smisc,2pletab,reb1/device,vector,onplcrack!时程后处理/post26nsol,2,69,u,yprod,3,2,,,,,,-1prod,4,1,,,,,,p0/1000xvar,3plvar,4分离式! 关闭压碎,keyopt(1)=0,keyopt(7)=1 finish/clear/config,nres,2000/prep7! 定义单元与材料性质et,1,solid65,,,,,,,1et,2,link8mp,ex,1,13585mp,prxy,1,0.2fc=14.3ft=1.43tb,concr,1tbdata,,0.5,0.95,ft,fc tb,miso,1,,11 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.0e5mp,prxy,2,0.3tb,biso,2tbdata,,300,0pi=acos(-1)r,1,0.25*pi*22*22r,2,0.25*pi*22*22/2r,3,0.25*pi*10*10r,4,0.25*pi*10*10/2!创建几何模型blc4,,,150/2,300,2000/2 *do,i,1,9 wpoff,,,100vsbw,all*enddowpcsys,-1 wpoff,,,50vsbw,all wpcsys,-1 wprota,,-90 wpoff,,,30 vsbw,all wpoff,,,240 vsbw,all wpcsys,-1 wpoff,30 wprota,,,90 vsbw,all wpcsys,-1!划分钢筋网格elemsiz=50 lsel,s,loc,x,30 lsel,r,loc,y,30 cm,zj,line latt,2,1,2 lesize,all,elemsiz lsel,s,loc,x,75 lsel,r,loc,y,30 cm,zjb,line latt,2,2,2 lesize,all,elemsiz lsel,s,loc,x,30 lsel,r,loc,y,270 cm,jlj,line latt,2,3,2 lesize,all,elemsiz lsel,s,tan1,zlsel,r,loc,y,30,270 lsel,r,loc,x,30,70 lsel,u,loc,z,50 cm,gj,line latt,2,3,2 lesize,all,elemsiz lsel,s,loc,z,0 lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 cm,gjb,line latt,2,4,2 lesize,all,elemsiz lsel,all cmsel,s,zj cmsel,a,zjb cmsel,a,jlj cmsel,a,gj cmsel,a,gjb cm,gj,line lmesh,alllsel,all/eshape,1 eplot!划分混凝土网格vatt,1,,1 mshkey,1 esize,elemsiz vmesh,all allsel,all!施加荷载和约束lsel,s,loc,y,0lsel,r,loc,z,900dl,all,,uyasel,s,loc,z,0da,all,symmasel,s,loc,x,75da,all,symmp0=180000q0=p0/150/100 asel,s,loc,z,0,50 asel,r,loc,y,300 sfa,all,1,pres,q0 allsel,all!求解控制/soluantype,0nsubst,100 outres,all,all autos,onneqit,50cnvtol,u,,0.015 solve!进入post1查看结果/post1set,lastpldisp,1esel,s,type,,2 etable,saxl,ls,1 plls,saxl,saxlesel,s,type,,1/device,vector,on plcrack!时程后处理/post26nsol,2,205,u,y prod,3,2,,,,,,-1 prod,4,1,,,,,,p0/1000 xvar,3plvar,4预应力混凝土结构!实体切分法finish/clear/prep7!定义变量emst=1.95e5as=139tf=180e3denss=7921e-12 emcon=3.45e4 densc=2300e-12r0=9000b=150h=200d0=40dd=200-2*d0 spanl=3000d1=(39-3*sqrt(29))/35*dd-d0 !定义单元与材料性质et,1,solid95et,2,link8mp,ex,1,emconmp,prxy,1,0.2mp,dens,1,denscmp,ex,2,emstmp,prxy,2,0.3mp,dens,2,denssr,1r,2,as,tf/(emst*as)*1.0271!切分体形成力筋线blc4,,,b,h,spanllsel,nonek,9,-10,h-d0k,10,-10,-d1,spanl/2k,11,-10,h-d0,spanll,9,10l,10,11lfillt,13,14,r0lcomb,allk,14,b+10,h-d0l,9,14adrag,13,,,,,,14lsel,allvsba,1,7wpoff,b/2wprota,,,90vsbw,allwpcsys,-1 wpoff,,,spanl/2 vsbw,all!划分单元网格lsel,s,radius,,r0 lsel,r,loc,x,b/2latt,2,2,2lesize,all,40 lmesh,allesize,40vatt,1,1,1 mshape,0 mshkey,1 vmesh,allallsel,all!定义约束和荷载并求解lsel,s,loc,z,0lsel,a,loc,z,spanl lsel,r,loc,y,0dl,all,,uydk,kp(0,0,0),ux,,,,uz dk,kp(0,0,spanl),ux allsel,all/soluacel,,9800solvefinish!查看结果/post1pldisp,1etable,sigs,ls,1 plls,sigs,sigs,1 path,kzhx,2 ppath,1,,0,0,spanl/2 ppath,2,,b,0,spanl/2 pdef,sigc,s,z plpath,sigcpath,kzsx,2ppath,1,,b/2,0,spanl/2 ppath,2,,b/2,h,spanl/2 pdef,sigc,s,z plpath,sigc!节点耦合法finish/clear/prep7!定义变量emst=1.95e5as=139tf=180e3denss=7921e-12 emcon=3.45e4 densc=2300e-12r0=9000b=150h=200d0=40dd=200-2*d0spanl=3000d1=(39-3*sqrt(29))/35*dd-d0 !定义单元与材料性质et,1,solid95et,2,link8mp,ex,1,emconmp,prxy,1,0.2mp,dens,1,denscmp,ex,2,emstmp,prxy,2,0.3mp,dens,2,denssr,1r,2,as,tf/(emst*as)*1.0271!定义单元与材料性质et,1,solid95et,2,link8mp,ex,1,emconmp,prxy,1,0.2mp,dens,1,denscmp,ex,2,emstmp,prxy,2,0.3mp,dens,2,denssr,1r,2,as,tf/(emst*as)*1.01225!分别创建体和力筋线blc4,,,b,h,spanllsel,nonek,9,b/2,h-d0k,10,b/2,-d1,spanl/2k,11,b/2,h-d0,spanll,9,10l,10,11lfillt,13,14,r0lcomb,all wpoff,,,spanl/2lsbw,allcm,sline,linelsel,allvsbw,allwpcsys,-1!分别划分单元网格cmsel,s,slinelatt,2,2,2lesize,all,40lmesh,allesize,40vatt,1,1,1lsel,s,length,,h lesize,all,,,10lsel,allmshape,0mshkey,1vmesh,all!定义节点耦合自由度cmsel,s,slinensll,s,1cm,gsnode,node*get,senum,node,,count *dim,snode,,senum*dim,cnode,,senum*get,nodei,node,,num,min snode(1)=nodei*do,i,2,senumsnode(i)=ndnext(snode(i-1)) *enddoallsel,allcmsel,u,gsnode*do,i,1,senumnodei=snode(i)nodej=nnear(nodei) cnode(i)=nodejnsel,u,,,nodej*enddoallsel,all*do,i,1,senumcp,next,all,snode(i),cnode(i) *enddo!定义约束和荷载并求解lsel,s,loc,z,0lsel,a,loc,z,spanllsel,r,loc,y,0dl,all,,uydk,kp(0,0,0),ux,,,,uzdk,kp(0,0,spanl),ux allsel,all/soluacel,,9800solvefinish!查看结果/post1pldisp,1etable,sigs,ls,1 plls,sigs,sigs,1 path,kzhx,2 ppath,1,,0,0,spanl/2 ppath,2,,b,0,spanl/2 pdef,sigc,s,z plpath,sigcpath,kzsx,2ppath,1,,b/2,0,spanl/2 ppath,2,,b/2,h,spanl/2 pdef,sigc,s,z plpath,sigc!约束方程法finish/clear/prep7!定义变量emst=1.95e5as=139tf=180e3denss=7921e-12 emcon=3.45e4 densc=2300e-12r0=9000b=150h=200d0=40dd=200-2*d0spanl=3000d1=(39-3*sqrt(29))/35*dd-d0 !定义单元与材料性质et,1,solid95et,2,link8mp,ex,1,emconmp,prxy,1,0.2mp,dens,1,denscmp,ex,2,emstmp,prxy,2,0.3mp,dens,2,denssr,1r,2,as,tf/(emst*as)*1.0271!定义单元与材料性质et,1,solid95et,2,link8mp,ex,1,emconmp,prxy,1,0.2mp,dens,1,denscmp,ex,2,emstmp,prxy,2,0.3mp,dens,2,denssr,1r,2,as,tf/(emst*as)*1.01225!分别创建体和力筋线blc4,,,b,h,spanllsel,nonek,9,b/2,h-d0k,10,b/2,-d1,spanl/2k,11,b/2,h-d0,spanl l,9,10l,10,11lfillt,13,14,r0 lcomb,all wpoff,,,spanl/2 lsbw,allcm,sline,linelsel,allvsbw,allwpcsys,-1!分别划分单元网格cmsel,s,slinelatt,2,2,2lesize,all,40 lmesh,allesize,40vatt,1,1,1lsel,s,length,,h lesize,all,,,10lsel,allmshape,0 mshkey,1 vmesh,all!定义约束方程cmsel,s,slinensll,s,1ceintf,,ux,uy,uz allsel,all!定义约束和荷载并求解lsel,s,loc,z,0lsel,a,loc,z,spanl lsel,r,loc,y,0dl,all,,uydk,kp(0,0,0),ux,,,,uz dk,kp(0,0,spanl),ux allsel,all/soluacel,,9800solvefinish!查看结果/post1pldisp,1etable,sigs,ls,1 plls,sigs,sigs,1 path,kzhx,2 ppath,1,,0,0,spanl/2 ppath,2,,b,0,spanl/2 pdef,sigc,s,z plpath,sigcpath,kzsx,2ppath,1,,b/2,0,spanl/2 ppath,2,,b/2,h,spanl/2 pdef,sigc,s,z plpath,sigc。
钢筋混凝土梁ansys分析附命令流
钢筋混凝土非线性分析2015大作业上海交通大学陈明1、参数选择梁的截面宽度为200mm,上部配置2Φ8受压筋,混凝土的净保护层厚度为25 mm(从纵向钢筋外边缘算起),箍筋两端区采用8@100的双肢箍,中间区取8@200 双肢箍1)梁的截面高度选300mm;2)两加载间的距离选1000mm;3)混凝土选C30;4)纵向受拉钢筋配筋选218;2、描述选用的有限元模型及单元的特点采用ansys软件进行模拟计算,钢筋混凝土模型采用分离式模型,不考虑钢筋与混凝土之间的相对滑移。
混凝土采用solid65单元模拟,solid65用于模拟三维有钢筋或无钢筋的混凝土模型。
该单元能够计算拉裂和压碎。
在混凝土应用中,该单元的实体功能可以用于建立混凝土模型,同时,还可用加筋功能建立钢筋混凝土模型。
另外,该单元还可以应用于加强复合物和地质材料。
该单元由八个节点定义,每个节点有三个自由度:节点坐标系的x,y,z方向的平动。
至多可以定义三种不同规格的钢筋。
钢筋单元采用link180单元模拟,link180是一个适用于各类工程应用的三维杆单元。
根据具体情况,该单元可以被看作桁架单元、索单元、链杆单元或弹簧单元等等。
本单元是一个轴向拉伸一压缩单元,每个节点有三个自由度:节点坐标系的x,y,z方向的平动。
本单元是一种顶端铰接结构,不考虑单元弯曲。
本单元具有塑性、蠕变、旋转、大变形和大应变功能。
缺省时,当考虑大变形时任何分析中LINK180单元都包括应力刚化选项。
3、描述选用的混凝土与钢筋粘结滑移本构关系的具体形式、参数等。
钢筋的应力应变关系曲线考虑到极限塑性应变最大值为0.01,钢筋本构模型采用多线性模型kinh,初始弹性模量为Es=200000Mpa,强化系数为0.001。
混凝土的应力应变关系曲线混凝土选用各向同性的miso模型,当计入下降端时,程序报错,所以只取了前面的上升段,用5段折线模拟混凝土应力应变曲线。
不考虑混凝土与钢筋之间的相对滑移4、迭代方法和收敛标准。
ANSYS,ADINA在钢筋混凝土结构非线性分析中的应用与算例分析
表1跨中截面的应变和应力计算 结果对比
结论与展望
结论与展望
通过以上分析和算例对比,我们可以得出以下结论:ANSYS和ADINA在钢筋混 凝土结构非线性分析中均具有较高的准确性和可靠性,二者的计算结果非常接近。 这些软件提供了丰富的建模功能、强大的求解器和灵活的材料本构模型选择,使 得它们能够广泛应用于各种复杂的钢筋混凝土结构非线性分析问题中。
ANSYS和ADINA在钢筋混凝土结构非线性分析中的应用步骤和常见问题
钢筋混凝土结构非线性分析基本原理和方法
ANSYS和ADINA都是广泛用于工程模拟的有限元软件,它们在钢筋混凝土结构 非线性分析中具有相似的应用步骤。具体而言,可以分为以下几个步骤:
钢筋混凝土结构非线性分析基本原理和方法
1、建立模型:根据实际结构,建立相应的计算模型,包括几何形状、材料属 性、边界条件等。
谢谢观看
ANSYS,ADINA在钢筋混凝土结构 非线性分析中的应用与算例分
析
目录
01 引言
03 表1跨中截面的应变 和应力计算结果对比
02 算例分析 04 结论与展望 Nhomakorabea 内容摘要
ANSYS和ADINA在钢筋混凝土结构非线性分析中的应用与算例分析
引言
引言
钢筋混凝土结构非线性分析在工程实践中具有重要意义,它能够帮助工程师 们更准确地预测结构的性能和行为。在众多非线性分析软件中,ANSYS和ADINA是 最常用的两个工具。本次演示将介绍这两个软件在钢筋混凝土结构非线性分析中 的应用原理、方法和算例,并对其优势和不足进行分析,以期为工程师们提供有 益的参考。
ANSYS弹性及塑性分析(详细、全面分析)1
目录什么是塑性 (1)路径相关性 (1)率相关性 (1)工程应力、应变与真实应力、应变 (1)什么是激活塑性 (2)塑性理论介绍 (2)屈服准则 (2)流动准则 (3)强化准则 (3)塑性选项 (5)怎样使用塑性 (6)ANSYS输入 (7)输出量 (7)程序使用中的一些基本原则 (8)加强收敛性的方法 (8)查看结果 (9)塑性分析实例(GUI方法) (9)塑性分析实例(命令流方法) (14)弹塑性分析在这一册中,我们将详细地介绍由于塑性变性引起的非线性问题--弹塑性分析,我们的介绍人为以下几个方面:∙什么是塑性∙塑性理论简介∙ANSYS程序中所用的性选项∙怎样使用塑性∙塑性分析练习题什么是塑性塑性是一种在某种给定载荷下,材料产生永久变形的材料特性,对大多的工程材料来说,当其应力低于比例极限时,应力一应变关系是线性的。
另外,大多数材料在其应力低于屈服点时,表现为弹性行为,也就是说,当移走载荷时,其应变也完全消失。
由于屈服点和比例极限相差很小,因此在ANSYS程序中,假定它们相同。
在应力一应变的曲线中,低于屈服点的叫作弹性部分,超过屈服点的叫作塑性部分,也叫作应变强化部分。
塑性分析中考虑了塑性区域的材料特性。
路径相关性:即然塑性是不可恢复的,那么这种问题的就与加载历史有关,这类非线性问题叫作与路径相关的或非保守的非线性。
路径相关性是指对一种给定的边界条件,可能有多个正确的解—内部的应力,应变分布—存在,为了得到真正正确的结果,我们必须按照系统真正经历的加载过程加载。
率相关性:塑性应变的大小可能是加载速度快慢的函数,如果塑性应变的大小与时间有关,这种塑性叫作率无关性塑性,相反,与应变率有关的性叫作率相关的塑性。
大多的材料都有某种程度上的率相关性,但在大多数静力分析所经历的应变率范围,两者的应力-应变曲线差别不大,所以在一般的分析中,我们变为是与率无关的。
工程应力,应变与真实的应力、应变:塑性材料的数据一般以拉伸的应力—应变曲线形式给出。
用Ansys或Abaqus分析钢管混凝土结构或构件
用Ansys或Abaqus分析钢管混凝土结构或构件用Ansys或Abaqus分析钢管混凝土结构或构件以上两个软件国外都有人用来分析钢管混凝土结构,但建模的方法不尽相同。
关键在于钢管和混凝土本构关系的选取以及两者之间的界面处理方法,各位有没有这方面的经验能向我们大家介绍一下。
==========程序中大概只有Drucker-Prager比较适合描述受约束混凝土的本构关系,因为这个模型可以考虑 hydrostatic stress (流体静应力)的影响。
在程序中,需要输入cohesion, angle of internal friction,(one more for ANSYS is the angle of dilatancy)。
值得注意的是,两个软件确定这几个参数的公式各不相同,很是令人头疼。
其实user manuals不可能给出明确的表达式,因为到目前为止,好像没有研究把钢管的强度,混凝土的强度,含钢率等等因素(i.e. the confinement)全部在Drucker-Prager 中考虑进去。
至于两种材料的界面,日本的 Hanbin Ge曾用link element来模拟,但在他的文章中,没有详细的描述。
轴压状况下,好像可以忽略滑移。
偏压可能情况有所不同。
==========韩教授书上的混凝土应力-应变关系,可以简单理解为单向受力的混凝土本构关系(考虑了钢管的约束),因此不能用于多向应力状态下混凝土的有限元分析。
材料非线性有限元分析,需要定义材料的屈服面,流动准则,强化准则,等等。
对受约束的混凝土,还要考虑体积膨胀,钢管对它的约束等因素。
显然,不是一个简单的应力-应变曲线所能概括的。
==========三向有限元分析,需要定义屈服面、流动准则和强化准则等等,而考虑钢管约束的混凝土本构关系,只是应力-应变关系。
对钢管混凝土的有限元分析,主要困难是如何定义屈服面,和模拟两个材料之间的滑移,我曾经用过接触分析(contact analysis)来求轴压构件的承载力,发现最大承载力能够比较精确地求得,但是精确的荷载-位移曲线很难获得,因为商用软件(Ansys\Marc)里面的D-P模型是塑性模型。
用ANSYS建立钢筋混凝土梁模型 (2)
用ANSYS建立钢筋混凝土梁模型问题描述:钢筋混凝土梁在受到中间位移荷载的条件下的变形以及个组成部分的应力情况。
一、用合并节点的方法模拟钢筋混凝土梁1.用solid65号单元以及beam188单元时材料特性钢材的应力应变关系混凝土的弹性模量采用线弹性B=150mmH=300mmL=2000mmP=5mm位移图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.57e1tbpt,,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,2.secnum,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中的实现周岑孙利民(同济⼤学⼟⽊⼯程防灾国家重点实验室200092)摘要钢筋混凝⼟结构是现代⼟⽊⼯程中最常⽤的结构形式。
本⽂针对运⽤ANSYS进⾏钢筋混凝⼟结构的弹塑性分析,通过与理论解⽐较,依据分析对象的结构层次(结构、构件)、分析类型(静⼒单调加载、反复加载)、荷载⽔平(线弹性、弹塑性),讨论了单元类型、材料模型及模型参数的选取,必要时甚⾄采⽤UPF等⼆次开发⼯具进⾏分析。
分析表明,合理的模型可以得到令⼈满意的结果。
关键词钢筋混凝⼟结构弹塑性 ANSYSRealization of RC Structure Elasto-plastic Analysis withANSYSZhou Cen Sun Limin(State Key Laboratory for Disaster Reduction in Civil Engineering )Abstract:RC structure is the most common structure type in modern civil engineering. In this paper, how to analyze RC structure elasto-plastic analysis with ANSYS is discussed. Compared with theoretical results, it is discussed how to select element type, material model and parameter based on the structurelevel (whole structure or member), analysis type (under static monotone load or cyclic load), load level(linear elastic or elasto-plastic) and UPF if necessary. The analysis shows that satisfactory results maybe obtained from rational models.1 前⾔钢筋混凝⼟材料由于结实且价格低廉,已经成为⼟⽊⼯程结构中采⽤得最多的材料。
钢筋混凝土梁设计的ANSYS软件分析
钢筋混凝土梁设计的ANSY S 软件分析彭兵田,刘立时,肖 洁(广东省航盛工程有限公司,广东 广州 510230)收稿日期:2002-04-30作者简介:彭兵田(1974-),男,湖南双峰人,助理工程师,从事路桥施工管理。
摘 要:本文通过ANSY S 软件对简支梁的受力和应变情况作了全面分析,提出将梁体看作理想弹塑性材料进行计算和设计的优点。
同时也指出了利用ANSY S 软件对简支梁分析时的不足。
为钢筋混凝土梁的设计提供了另一种思路。
关键词:ANSY S;钢筋混凝土梁;理想弹塑性中图分类号:U416.02 文献标识码:B 文章编号:1002-4972(2002)06-0012-02ANSYS Softw are Analysis of R einforced Concrete B eam DesignPE NG Bing -tian ,LI U Li -shi ,XI AO Jie(G uangdong Provincial Hangsheng Engineering C o.,Ltd.,G uangzhou 510230,China )Abstract :Through a com prehensive analysis of the stress and strain of sim ple beam by ANSY S s oftware ,itputs forward the advantages in calculation and design when beams are considered as ideal elastic and plastic materials.Meanwhile ,the inadequateness of ANSY S s oftware for sim ple beam analysis is als o pointed out.Another idea in rein 2forced concrete beam design is provided.K ey w ords :ANSY S ;rein forced concrete beam ;ideal elasticity and plasticity 目前,在我国的混凝土规范中,钢筋混凝土梁的设计分正截面抗弯和斜截面抗剪分别进行的。
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工程结构数值分析(钢筋混凝土简支梁分离式整体模型命令流)------------* FINISH/CLEAR,START/TITLE,discrete model of integrity/PREP7!1***************定义单元与材料性质******************ET,1,SOLID65,,,,,,,1ET,2,LINK8MP,EX,1,13585$MP,PRXY,1,0.2 !混凝土MP,EX,2,2.0E5$MP,PRXY,2,0.3 !钢筋FC=14.3$FT=1.43TB,CONCR,1$TBDATA,,0.5,0.95,FT,-1TB,MISO,1,,11TBPT,,0.0002,FC*0.19$TBPT,,0.0004,FC*0.36TBPT,,0.0006,FC*0.51$TBPT,,0.0008,FC*0.64TBPT,,0.0010,FC*0.75$TBPT,,0.0012,FC*0.84TBPT,,0.0014,FC*0.91$TBPT,,0.0016,FC*0.96TBPT,,0.0018,FC*0.99$TBPT,,0.0020,FC*1.00TBPT,,0.0033,FC*1.00TB,BISO,2$TBDATA,,360,0PI=ACOS(-1)R,1,PI*22*22$R,2,PI*10*10$R,3,PI*10*10R,4,0!2****************创建几何模型*********************BLC4,,,150,300,2000*DO,I,1,19$WPOFF,,,100$VSBW,ALL$*ENDDO !切箍筋位置WPCSYS,-1WPOFF,,,950$VSBW,ALL$WPOFF,,,100$VSBW,ALL !切加载位置WPCSYS,-1WPROTA,,-90$WPOFF,,,30$VSBW,ALL !切出拉区钢筋竖向位置WPOFF,,,240$VSBW,ALL !切出压出钢筋竖向位置WPCSYS,-1WPROTA,,,90WPOFF,,,30$VSBW,ALLWPOFF,,,45$VSBW,ALLWPOFF,,,45$VSBW,ALLWPCSYS,-1!3****************划分钢筋网格*********************ELEMSIZ=50LSEL,S,LOC,X,30$LSEL,R,LOC,Y,30 !拉区左侧钢筋,定义组件,设置属性CM,ZJLZ,LINE$LATT,2,1,2$LESIZE,ALL,ELEMSIZLSEL,S,LOC,X,120$LSEL,R,LOC,Y,30 !拉区右侧钢筋,定义组件,设置属性CM,ZJL Y,LINE$LATT,2,1,2$LESIZE,ALL,ELEMSIZLSEL,S,LOC,X,75$LSEL,R,LOC,Y,30 !拉区中间钢筋,定义组件,设置属性CM,ZJLM,LINE$LATT,2,1,2$LESIZE,ALL,ELEMSIZLSEL,S,LOC,X,30$LSEL,R,LOC,Y,270 !压区左侧钢筋,定义组件,设置属性CM,ZJYL,LINE$LATT,2,3,2$LESIZE,ALL,ELEMSIZLSEL,S,LOC,X,120$LSEL,R,LOC,Y,270 !压区右侧钢筋,定义组件,设置属性CM,ZJYR,LINE$LATT,2,3,2$LESIZE,ALL,ELEMSIZLSEL,S,TAN1,ZLSEL,R,LOC,Y,30,270$LSEL,R,LOC,X,30,120 !箍筋,定义组件,设置属性LSEL,U,LOC,Z,950$LSEL,U,LOC,Z,1050$LSEL,U,LOC,X,75CM,GJ,LINE$LATT,2,3,2$LESIZE,ALL,ELEMSIZLSEL,ALLCMSEL,S,ZJLZ$CMSEL,A,ZJLM$CMSEL,A,ZJL YCMSEL,A,ZJYL$CMSEL,A,ZJYR$CMSEL,A,GJLMESH,ALL$LSEL,ALL!/ESHAPE,1$EPLOT!4****************划分混凝土网格*********************V ATT,1,4,1$MSHKEY,1$ESIZE,ELEMSIZ$VMESH,ALL$ALLSEL,ALL!5****************施加约束*********************LSEL,S,LOC,Y,0$LSEL,R,LOC,Z,100DL,ALL,,ALLLSEL,S,LOC,Y,0$LSEL,R,LOC,Z,1900DL,ALL,,ALL!6****************施加荷载*********************PO=180E3$QO=PO/150/100ASEL,S,LOC,Z,950,1050$ASEL,R,LOC,Y,300SFA,ALL,1,PRES,QOALLSEL,ALL!7******************************求解******************************************* /SOLUANTYPE,0$NSUBST,100OUTRES,ALL,ALL$AUTOS,ONNEQIT,50CNVTOL,U,,0.015SOLVE一、书上P482是1/4模型,上述模型是整体模型数值分析命令流二、改动的地方有①TBPT,,0.0033,FC*1.00,书上为TBPT,,0.0033,FC*0.85。
ansys建模实例-混凝土箱梁(命令流)
/TITLE, BRIDGE XX/PREP7! 定义单元类型ET,1,SHELL63ET,2,SOLID45ET,3,LINK8!定义1,钢,2,混凝土,3,无限刚度MP,EX,1,2.1E11MP,DENS,1,7。
8E3MP,PRXY,1,0.3mp,alpx,1,1.0e—5MP,EX,2,3。
5E10MP,DENS,2,2.6E3MP,PRXY,2,0。
1667MP,EX,3,3。
5E15MP,DENS,3,7.8E3MP,PRXY,3,0。
3*set,nz,154!定义箱梁N1R,1,7,7,7,7!定义箱梁N2R,2,8,8,8,8!定义箱梁N3R,3,4,4,4,4!定义混凝土solid65的实常数R, 4!定义钢绞线的实常数R, 5, 140!降温=预加力/线膨胀系数*钢筋截面积*弹性模量tem1=—1190.4e6/(1。
0e-5*2。
11e11)tem2=-1339。
2e6/(1。
0e—5*2.11e11)!定义纵向梯度参数表*dim, z, array, nz*SET,z(1),130,260,371,482,593,702,811,920,1029,1138,1247,1356,1465,1574,1683,1787.19,1891。
38,1995.57*SET,z(19),2099。
76,2203。
95,2308。
14,2412.33,2516.52,2620.71,2724。
9,2829。
09,2933.28,3037。
47,3141.66,3245.85,*SET,z(31),3350.04,3454。
23,3558.42,3662.61,3766。
8,3870。
99*SET,z(37),3975。
18,4079.37,4183.56,4287.75,*SET,z(41),4391.94,4496.13,4600。
32,4704。
51,4808.7,4912。
89,5017.08,5121.27,5225。
ansys分析混凝土的若干问题
ansys分析混凝土的若干问题1. 讨论两种Ansys求极限荷载的方法(1)力加载可以通过对应的方法(比如说特征值屈曲)估计结构的极限荷载的大致范围,然后给结构施加一个稍大的荷载,打开自动荷载步二分法进行非线性静力分析,最后计算会因不收敛终止,则倒数第二个子步对应的就是结构的极限荷载;另外,也可以选择弧长法,采用足够的子步(弧长法可以一直分析到极限承载力之后的过程)同样可以从绘制的荷载位移曲线或计算结果中找出结构的极限荷载。
(2)位移加载给结构施加一个比较大的位移,打开自动荷载步二分法进行非线性分析,保证足够的子步数,这样也可以分析到极限荷载以后,通过绘制荷载位移曲线或查看相应结果文件也可知道结构的极限荷载。
希望众高手讨论一下(1)弧长法求极限荷载的收敛性问题,如何画到荷载位移曲线的下降段?(2)位移法求极限荷载的具体步骤?2. 需要注意的问题1. 由于SOLID 65单元本身是基于弥散裂缝模型和最大拉应力开裂判据,因此在很多情况下会因为应力集中而使混凝土提前破坏,从而和试验结果不相吻合,因此,在实际应用过程中应该对单元分划进行有效控制,根据作者经验,当最小单元尺寸大于5cm 时,就可以有效避免应力集中带来的问题;2. 支座是另一个需要注意的问题。
在有限元分析中,很多时候约束是直接加在混凝土节点上,这样很可能在支座位置产生很大的应力集中,从而使支座附近的混凝土突然破坏,造成求解失败。
因此,在实际应用过程中,应该适当加大支座附近单元的尺寸或者在支座上加一些弹性垫块,避免支座的应力集中;3. 六面体的SOLID 65 单元一般比四面体的单元计算要稳定且收敛性好,因此,只要条件允许,应该尽量使用六面体单元;4. 正确选择收敛标准,一般位移控制加载最好用位移的无穷范数控制收敛,而用力控制加载时可以用残余力的二范数控制收敛。
在裂缝刚刚出现和接近破坏的阶段,可以适当放松收敛标准,保证计算的连续性;3. 关于下降段的问题1)在实际混凝土中都有下降段,但是在计算的时候要特别小心下降段的问题。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
!(1)工作环境设置
/FILENAME,RC-BEAM !指定工作文件名
/TITLE,ALAL YSIS OF A RC-BEAM !指定图形标题
!(2)进行前处理器
/prep7
!(3)定义单元类型
ET,1,LINK8 !定义钢筋单元
ET,2,SOLID65 !定义混凝土单元
ET,3,MESH200 !用于拉伸成体单元KEYOPT,3,1,6
!(4)定义钢筋截面积
r,1,28.3
r,2,50.3
r,3,314.1
!(5)为solid65单元定义一个实参数组
r,4,
!(6)定义混凝土材料
MP,EX,2,2.55E10
MP,PRXY,2,0.3
TB,CONC,2,1,9, !定义混凝土的破坏参数TBDA TA,,0.3,0.55,1.55E6,-1,,
TBDA TA,,,,0.6
!(7)定义钢筋材料模型及参数
mp,ex,2,2e5 !纵向受拉钢筋材料
mp,prxy,2,0.3
tb,bkin,2,1,2,1
tbdata,,350
mp,ex,3,2e5 !横向箍筋,架立钢筋材料mp,prxy,3,0.25
tb,bkin,3,1,2,1
tbdata,,200
!(8) 创建以及复制节点
/pnum,node,1
/pnum,elem,1
n,1
n,9,200
fill,1,9
ngen,11,9,1,9,1,,30
ngen,11,99,1,99,1,,,-150
/view,1,1,1,1
!(9)建立箍筋单元
type,1
real,1
mat,3
!水平箍筋
*do,i,11,16,1
e,i,i+1
e,i+(83-11),i+(83-11)+1
*enddo
!坚直箍筋
*do,i,11,74,9
e,i,i+9
e,i+6,i+6+9
*enddo
!产生整个模型的箍筋
egen,11,99,all
!(10)建立架立筋以及纵筋单元
!创建上部的架立钢筋单元
*do,i,83,node(25,270,-1500+150),99 e,i,i+99
e,i+6,i+6+99
*enddo
!纵向受拉钢筋单元的属性
type,1
real,3
mat,2
!创建纵筋单元
*do,i,11,node(25,30,-1500+150),99 e,i,i+99
e,i+3,i+3+99
e,i+6,i+6+99
*enddo
/view,1,1,1,1
/pnum,elem,0
/pnum,node,0
/eshape,1
eplot
!(11)建立混凝土剖面并划分网格
k,1
K,2,200,
k,3,200,300
k,4,,300
a,1,2,3,4
lsel,s,loc,y,0
lsel,a,loc,y,300
lesize,all,,,8
lsel,all
lsel,s,loc,x,0
lsel,a,loc,x,200
lesize,all,,,10
type,3
amesh,all
!(12)拉伸形成混凝土单元type,2
real,3
mat,1
extopt,esize,20
extopt,aclear,1
vext,all,,,,,-1500
/pnum,mat,1
/pnum,node,0
/pnum,elem,0
eplot
allsel
!(13)合并压缩节点编号nummrg,all
numcmp,all
eplot
!(14)施加支座约束nsel,s,loc,y,0
nsel,r,loc,z,-1500+75
d,all,uy
d,all,ux
!(15)施加对称面约束asel,s,loc,z,0
da,all,symm
!(16)施加载荷
nsel,all
nsel,s,loc,y,300
nsel,r,loc,z,-450
d,all,uy,-30
allsel
!(17)退出前处理器
fini
!(18) 进入求解器
/solu
!(19)求解器选项设置
nlgeom,on
nsubst,200
outres,all,all
neqit,50
pred,on
cnvtol,f,,0.05,2,0.5
allsel
!(20)求解并退出求解器
solve
fini
!(21)进入通用处理器并读入最后一个子步的结果/post1
set,last
!(22)后处理操作
plnsol,u,y,0,1
/device,vector,1
plcrack,0,0
!(23)退出通用后处理器
finish。