砌体结构Ansys命令流.docx

!mono brick wall!solid65-整体式建模!*Finish/clear/units,si/com,Structural/config,nres,400000/prep7/title,mono-brick-wall!*depth_sup=0.24 !弹性支垫240mm!单元属性!for brick, mortar and leveling coatet,1,solid65!*KEYOPT,1,1,0KEYOPT,1,5,0KEYOPT,1,6,0KEYOPT,1,7,1!for elastic supportet,2,solid45!*-砌体材料-!MU10,M5砌体材料属性mp,prxy,1,0.15mp,dens,1,1700mp,ex,1,0.24e10!砌体屈服准则TB,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.02!TBPLOT,MKIN,1!砌体破坏准则!hntrl=0.13E6hntrl=0.21E6tb,concr,1,4,tbdata,0.2,0.9,hntrl,-1 !不考虑压碎!*-弹性支撑-!弹性支撑暂取钢筋材料属性mp,ex,2,2.0e11mp,dens,2,7.600 mp,prxy,2,0.30!real constantsr,1!*-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, !brick wallvsel,all !必须再次选择所有单元VSWEEP,ALL savenumcmp,allallsel!*/device,vector,1/eshape,1!*-耦合-allsel,all!elastic support nodes setesel,s,type,2$nsle,s,all$nsel,r,loc,z,3.12nsel,u,loc,x,-0.24$cm,slavenode,node!other node setesel,s,type,1$nsle,s,all$cm,masternode,node/input,ucouple,macsavefinish耦合文件ucouple,mac：!*-开始-allsel !最好保留这句命令!*将从属节点编号依次存入数组*cmsel,s,slavenode*get,count1_node,node,0,count*del,slave_node*dim,slave_node,array,count1_node*get,slave_node(1),node,0,num,min*do,i,2,count1_node slave_node(i)=ndnext(slave_node(i-1)*enddo !*将主节点编号依次存入数组*allselcmsel,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_node master_node(i)=ndnext(master_node(i-1)*enddo!*将与从属节点耦合的节点数组初始化*del,cp_node*dim,cp_node,array,count1_node*do,i,1,count1_node cp_node(i)=0*enddo!*开始选择程序*allselcmsel,s,masternode*do,i,1,count1_node kk=1 k=1 *dowhile,kk k=nnear(slave_node(i) nsel,s,cp,all kk=nsel(k)+0.001 allsel cmsel,s,masternode nsel,u,node,k cm,masternode,node *enddo cp_node(i)=k *enddo!*选择完毕*!*开始耦合*allsel,all/prep7*do,i,1,count1_node cp,next,UX,slave_node(i),cp_node(i) cp,next,UY,slave_node(i),cp_node(i) cp,next,UZ,slave_node(i),cp_node(i)*enddo!*耦合完毕*求解方法NR法：brick wall NR solution!displacement load!用NR法!*-solution-/soluantype,0!-constraint-nsel,s,loc,z,0d,all,all!*esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uy!*!esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24!d,all,uzlumpm,onsolcontrol,onpred,onlnsrch,onautots,1 !打开自动时间步控制ncnv,2 !如果不收敛时结束而不退出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,0,0.03*compressionallsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,2,200,1neqit,15KBC,0nlgeom,1cnvtol,f,0.01,2,0.05allsellswrite,1!*-2nd step-nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6!d,all,ux,2.79E-3d,all,ux,1.0E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,100neqit,15KBC,0nlgeom,1cnvtol,f,0.05,2,0.05allsellswrite,2!*-3rd step-nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6!d,all,ux,2.79E-3d,all,ux,2.5E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,100neqit,15KBC,0nlgeom,1cnvtol,f,0.09,2,0.25allsellswrite,3!*-4th step-nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6!d,all,ux,2.79E-3d,all,ux,4.0E-3allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,100neqit,15KBC,0nlgeom,1cnvtol,f,0.1,2,0.5allsellswrite,4!*allsel lssolve,1,4,1FINISH也可以用弧长法，不过没有一定经验这比较难控制：!mono brick wall-1 arclength solution!displacement load!*-solution-/soluantype,0!-constraint-nsel,s,loc,z,0d,all,all!*esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uy!*lumpm,onsolcontrol,onncnv,2 !如果不收敛时结束而不退出compression=0.24E7!*-load steps-!*-first step-acel,9.8nsel,s,loc,z,3.36$nsel,r,loc,x,0.001,3.12esln,s,0,allsfe,all,6,pres,0,0.03*compressionallsel,all!*OUTRES,ERASEOUTRES,ALL,ALLarclen,on,1E-5, !使用弧长法NSUBST,200neqit,15KBC,0nlgeom,1cnvtol,f,0.1,2,0.5allsel,alllswrite,1!*-second step-esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uz!*nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6d,all,ux,2.50E-3!d,all,ux,3.51E-3allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLarclen,on,1E-5, !使用弧长法kbc,0NSUBST,250!arctrm,u,0.0028,313,uxneqit,15nlgeom,1!cnvtol,f,0.005,2,0.01!cnvtol,f,0.05,2,0.5cnvtol,f,0.1,2,0.5allsellswrite,2!*-third step-esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uz!*nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6!d,all,ux,2.79E-3!d,all,ux,3.18E-3d,all,ux,3.00E-3allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLarclen,on,1E-7, !使用弧长法kbc,0NSUBST,2000!arctrm,u,0.0028,313,uxneqit,15nlgeom,1!cnvtol,f,0.005,2,0.01!cnvtol,f,0.05,2,0.5cnvtol,f,0.5,2,0.8allsellswrite,3allsel!*-fourth step-esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uz!*nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6d,all,ux,3.51E-3allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLarclen,on,1E-7, !使用弧长法kbc,0NSUBST,300!arctrm,u,0.0028,313,uxneqit,15nlgeom,1!cnvtol,f,0.005,2,0.01!cnvtol,f,0.05,2,0.5cnvtol,f,0.6,2,0.9allsellswrite,4allsel!*lssolve,1,4,1FINISH后处理：/post26esel,s,type,1eplotlines,5000nsol,2,326,u,x,ux!nsol,2,313,u,x,ux*do,i,1,14 jj=node(i-1)*0.24,0.24,0) kk=node(i-1)*0.24,0,0)! ll=node(225.5+(i-1)*225.5,3300,0) rforce,3,jj,f,x,fxadd,6,3,6 rforce,4,kk,f,x,fxadd,6,4,6! rforce,5,ll,f,x,fx!add,6,5,6*enddoprod,7,2,1000prod,8,6,-1/1000/axlab,x,Displacement(mm)/axlab,y,F(kN)/XRANGE,0,4/YRANGE,0,100xvar,7plvar,8!rforce,5,340,f,x,fx!esol,5,87,221,f,x,fx2、循环水平荷载作用下的命令流：出求解文件不同外，其他都相同；!cyclic NR solution!displacement load!*-solution-/soluantype,0!-constraint-nsel,s,loc,z,0d,all,all!*esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uy!*esel,s,type,2$nsle,s,all$nsel,r,loc,x,3.36d,all,uy!*lumpm,onsolcontrol,onpred,onlnsrch,onautots,1 !打开自动时间步控制ncnv,2 !如果不收敛时结束而不退出compression=0.24E7!*-load steps-!*-1st step-acel,9.8nsel,s,loc,z,3.36$nsel,r,loc,x,0.001,3.119esln,s,0,allsfe,all,6,pres,0,0.03*compressionallsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,10neqit,15KBC,0nlgeom,1cnvtol,f,0.1,2,0.5allsel,alllswrite,1!*-初裂荷载20%做1次循环-!*-用NR法-!*-2nd step-esel,s,type,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uz!*esel,s,type,2$nsle,s,all$nsel,r,loc,x,3.36d,all,uz!*nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6d,all,ux,1.48E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,20,2000,10neqit,15NROPT ,FULL,OFFKBC,0nlgeom,1cnvtol,f,0.1,2,0.5allsellswrite,2!*-3rd step-nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,3.36d,all,ux,-1.48E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,40,2000,20neqit,15NROPT ,FULL,OFFKBC,0nlgeom,1cnvtol,f,0.1,2,0.5allsellswrite,3!*-4th step-nsel,s,loc,x,3.36ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,1.48E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,40,2000,20neqit,15NROPT ,FULL,OFFKBC,0nlgeom,1cnvtol,f,0.1,2,0.5allsellswrite,4!*-one cycle-!*-初裂荷载50%做1次循环-!*-用NR法-!*-5th step-nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,2.85E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,100,2000,10neqit,15NROPT ,FULL,OFFKBC,0nlgeom,1cnvtol,f,0.1,2,0.5allsellswrite,5!*-6th step-nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,3.36d,all,ux,-2.85E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,30neqit,15NROPT ,FULL,OFFKBC,0nlgeom,1cnvtol,f,0.1,2,0.5allsellswrite,6!*-7th step-nsel,s,loc,x,3.36ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,2.85E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,30neqit,15NROPT ,FULL,OFFKBC,0nlgeom,1cnvtol,f,0.1,2,0.5allsellswrite,7!*-two cycle-!*-初裂荷载80%做1次循环-!*-用NR法-!*-8th step-nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,3.5E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,1E-7, !使用弧长法kbc,0NSUBST,100,2000,20neqit,15KBC,0nlgeom,1cnvtol,f,0.5,2,0.9allsellswrite,8!*-9th step-nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,3.36d,all,ux,-3.5E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,1E-7, !使用弧长法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,0.5,2,0.9allsellswrite,9!*-10th step-nsel,s,loc,x,3.36ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,3.5E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,1E-7, !使用弧长法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,0.5,2,0.9allsellswrite,10!*-three cycle-!*-初裂荷载100%做1次循环-!*-用NR法-!*-11th step-nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,4.42E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,1E-7, !使用弧长法kbc,0NSUBST,100,2000,20neqit,15KBC,0nlgeom,1cnvtol,f,0.5,2,0.9allsellswrite,11!*-12th step-nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,3.36d,all,ux,-4.42E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,1E-7, !使用弧长法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,0.5,2,0.9allsellswrite,12!*-13th step-nsel,s,loc,x,3.36ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,4.42E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,1E-7, !使用弧长法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,0.5,2,0.9allsellswrite,13!*-four cycle-!*-61.4KN的荷载做1次循环-!*-用NR法-!*-14th step-nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,6.72E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,1E-7, !使用弧长法kbc,0NSUBST,100,2000,20neqit,15KBC,0nlgeom,1cnvtol,f,0.5,2,0.9allsellswrite,14!*-15th step-nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,3.36d,all,ux,-6.72E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,1E-7, !使用弧长法kbc,0NSUBST,200