混凝土非线性分析应力应变曲线输出解读Word格式文档下载.docx

1、graph variables. 在对话框上the first variable对应的交互框中输入2。点击确定键，则预想的应力-应变曲线就显示在屏幕上。ok!试试看！欢迎大家继续批评指正！/prep7 et,1,solid65 et,2,link8 et,3,solid45r,1r,2,236e-6r,3,75.4e-6 r,4,151.2e-6 !混凝土本构关系mp,ex,1,2.522e10mp,prxy,1,0.2 mp,dens,1,2500tb,concr,1,1 tbdata,0.5,0.95,5e6,-1tb,miso,1,1,7tbtemp,0 tbpt,0.00068,17.

2、15e6tbpt,0.00136,28.03e6tbpt,0.00204,33.62e6tbpt,0.00272,35.96e6tbpt,0.0034,36.53e6tbpt,0.00408,36.53e6纵筋的本构关系mp,ex,2,2.0e11mp,prxy,2,0.28 tb,bkin,2 tbdata,362e6,2e9！箍筋的本构关系mp,EX,3,2.1E11mp,PRXY,3,.27tb,bkin,3tbtemp,0tbdata,272.5e6,2.1e9 建立模型local,11,1,0,0,0,0,-90,0 ！建立局部坐标系，采用柱坐标，绕y轴顺时针转90度csys,11*

3、do,j,1,8n,j,0.185,j*360/8,-0.6*enddo ngen,3,100,1,8,1,0,0,0.3 ！复制节点，3次，编号增量，源节点18，0,0,0.3ngen,4,100,201,208,1,0,0,0.15ngen,7,100,501,508,1,0,0,0.3type,2 ！指定单元属性，建立基础箍筋单元real,3mat,3*do,i,1,6*do,j,1,7e,（i-1）*100+j,（i-1）*100+j+1*enddoe,（i-1）*100+1,（i-1）*100+8建立柱上箍筋单元real,4*do,i,1,5e,（i+5）*100+j,（i+5）*1

4、00+j+1e,（i+5）*100+1,（i+5）*100+8建立纵向钢筋real,2mat,2*do,i,1,11e,（i-1）*100+j,i*100+jallsel,all ！合并节点，压缩节点nummrg,allnumcmp,allwpcsys,1,11cylind,0,0.185,-0.6,0,0,360 ！建立圆柱，内外径i-or，高度z1-z2，起终角度cylind,0,0.185,0,0.45,0,360cylind,0.185,0.3,-0.6,0,0,360cylind,0,0.185,0.45,2.25,0,360cylind,0,0.185,2.25,2.55,0,36

5、0cylind,0.185,0.20,2.25,2.55,0,360cylind,0,0.185,2.55,2.565,0,360wprota,90 ！绕y轴逆时针转90度vsel,s,1,7,1 ！选择体，选择类型，item，comb，最小编号，最大编号，增量vsbw,all ！体相减numcmp,all wprota,90 ！绕z轴逆时针转90度vsel,s,1,14,1 vsbw,all vsel,s,1,28,1 vglue,all wprota,-90 wprota,-90 allsel,allvsel,s,1,16,1vsel,a,21,24,1vatt,1,1,1 vsel,s,

6、17,20,1vsel,a,25,28,1vatt,2,1,3 allsel,all lsel,s,1,8,1 ！选择线lsel,a,11,22,1 ！补选（a）线lsel,a,25,40,1 lsel,a,43,46,1 lesize,all,2 ！线尺寸为2lsel,s,47,50,1 lsel,a,53,56,1lsel,a,58,73,1lesize,all,2lsel,s,98,107,1 lsel,a,110,112,2lesize,all,1lsel,s,23,24,1 lsel,a,9,10,1lsel,s,78,82,1 lesize,all,3lsel,s,83,87,1

7、lesize,all,6lsel,s,41,42,1 lsel,a,88,92,1lsel,s,93,97,1 mshkey,1 mshape,0,3d vsel,s,1,28vmesh,all wprota,90 csys,0/solu ！进入加载、求解nsel,s,loc,x,0.3 ！选择半径在x=3处的节点，并施加约束d,all,allOUTRES,ALL,2 ！输出获得的解pred,on ！打开预应力选项cnvtol,f,0.05,2 !定义收敛条件,使用缺省的VALUE autots,1 !打开自动时间步控制 lnsrch,1 !打开线性搜索 ncnv,2 !如果不收敛时结束而不退

8、出 neqit,50 !每一子步中方程的迭代次数限值 savetime,1 ！定义第1载荷步nsubst,10asel,s,29,32,1sfa,all,pres,5468699 ！在已选择的面上施加压力荷载acel,9.8 ！重力加速度y=9.8kbc,1 ！阶跃加载方式lswrite,1 time,10 ！定义第2荷载步cp,1,ux,298,300,303 !将节点298,300,303ux向位移耦合d,298,ux,0.015 ！再借点298处施加ux位移0.015KBC,0NSUBST,20,1000000,10 lswrite,2time,20 ！定义第3载荷步d,298,ux,-

9、0.015NSUBST,100,1000000,10 lswrite,3time,30 ！定义第4载荷步d,298,ux,0.015lswrite,4time,40 ！定义第5载荷步d,298,ux,-0.015 lswrite,5time,50 ！定义第6载荷步lswrite,6time,60 ！定义第7载荷步lswrite,7time,70 ！定义第8载荷步lswrite,8time,80 lswrite,9time,90lswrite,10time,100lswrite,11time,110d,298,ux,0.03NSUBST,200,1000000,10 lswrite,12time

10、,120d,298,ux,-0.03lswrite,13time,130lswrite,14time,140lswrite,15time,150lswrite,16time,160lswrite,17time,170lswrite,18time,180lswrite,19time,190lswrite,20time,200lswrite,21time,210d,298,ux,0.045NSUBST,300,1000000,10 lswrite,22lssolve,1,22,1命令流：1、单调水平荷载作用下的命令流：建模：mono brick wallsolid65-整体式建模*Finish/c

11、lear/units,si/com,Structural/config,nres,400000/prep7/title,mono-brick-walldepth_sup=0.24 !弹性支垫240mm单元属性for brick, mortar and leveling coatet,1,solid65KEYOPT,1,1,0KEYOPT,1,5,0KEYOPT,1,6,0KEYOPT,1,7,1for elastic supportet,2,solid45*-砌体材料-MU10,M5砌体材料属性mp,prxy,1,0.15mp,dens,1,1700mp,ex,1,0.24e10砌体屈服准则T

12、B,MKIN,1TBTEMP,strainTBdata,1,0.2E-3,1.0E-3,2.0E-3,3.0E-3,4.8E-3TBTEMP,TBDATA,1,0.48E6,0.82E6,1.32E6,1.5E6,1.32E6TBDATA,1,0.48E6,1.31E6,2.11E6,2.4E6,2.11E6/XRANGE,0,0.02TBPLOT,MKIN,1砌体破坏准则hntrl=0.13E6hntrl=0.21E6tb,concr,1,4,tbdata,0.2,0.9,hntrl,-1 !不考虑压碎*-弹性支撑-弹性支撑暂取钢筋材料属性mp,dens,2,7.600 mp,prxy,2,

13、0.30real constants*-modeling-geometry modelblock,0,3.12,0,0.24,0,3.12block,-0.24,3.12,0,0.24,3.12,3.36Vsel,s,volu,1,2,1$aslv,s$lsla,slesize,all,0.24Vsel,s,volu,1$vatt,1,1,1, !brick wallVsel,s,volu,2$vatt,2,1,2, !vsel,all !必须再次选择所有单元VSWEEP,ALL allsel/device,vector,1/eshape,1*-耦合-elastic support nodes

14、 setesel,s,type,2$nsle,s,all$nsel,r,loc,z,3.12nsel,u,loc,x,-0.24$cm,slavenode,nodeother node setesel,s,type,1$nsle,s,all$cm,masternode,node/input,ucouple,macfinish耦合文件ucouple,mac：*-开始-allsel !最好保留这句命令*将从属节点编号依次存入数组*cmsel,s,slavenode*get,count1_node,node,0,count*del,slave_node*dim,slave_node,array,co

15、unt1_node*get,slave_node（1）,node,0,num,min*do,i,2,count1_nodeslave_node（i）=ndnext（slave_node（i-1）*将主节点编号依次存入数组*cmsel,s,masternode*get,count2_node,node,0,count*del,master_node*dim,master_node,array,count2_node*get,master_node（1）,node,0,num,min*do,i,2,count2_nodemaster_node（i）=ndnext（master_node（i-1）*

16、将与从属节点耦合的节点数组初始化*del,cp_node*dim,cp_node,array,count1_node*do,i,1,count1_nodecp_node（i）=0*开始选择程序*do,i,1,count1_node kk=1 k=1 *dowhile,kk k=nnear（slave_node（i）nsel,s,cp,allkk=nsel（k）+0.001 nsel,u,node,kcm,masternode,node cp_node（i）=k *选择完毕*开始耦合*cp,next,UX,slave_node（i）,cp_node（i）cp,next,UY,slave_node

17、（i）,cp_node（i）cp,next,UZ,slave_node（i）,cp_node（i）*耦合完毕*求解方法NR法：brick wall NR solutiondisplacement load用NR法*-solution-/soluantype,0-constraint-nsel,s,loc,z,0esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uyd,all,uzlumpm,onsolcontrol,onpred,onlnsrch,on打开自动时间步控制如果不收敛时结束而不退出compression=0.24E7*-load steps-*-1st-acel,9.8nsel,s,loc,z,3.36$nsel,r,loc,x,0.001,3.119esln,s,0,allsfe,all,6,pres