砌体结构Ansys命令流

!monobrickwall!solid65-整体式建模!*Finish/clear/units,si/com,Structural/config,nres,400000/prep7/title,mono-brick-wall!*depth_sup=0.24!弹性支垫240mm!单元属性!forbrick,mortarandlevelingcoatet,1,solid65!*KEYOPT,1,1,0KEYOPT,1,5,0KEYOPT,1,6,0KEYOPT,1,7,1!forelasticsupportet,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.600mp,prxy,2,0.30!realconstantsr,1!*--------------------------------modeling--------------------------------------!geometrymodelblock,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,!brickwallVsel,s,volu,,2$vatt,2,1,2,!brickwallvsel,all!必须再次选择所有单元VSWEEP,ALLsavenumcmp,allallsel!*/device,vector,1/eshape,1!*-----------------------------------耦合---------------------------------------allsel,all!elasticsupportnodessetesel,s,type,,2$nsle,s,all$nsel,r,loc,z,3.12nsel,u,loc,x,-0.24$cm,slavenode,node!othernodesetesel,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_nodeslave_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_nodemaster_node(i)=ndnext(master_node(i-1))*enddo!********将与从属节点耦合的节点数组初始化*****************del,cp_node*dim,cp_node,array,count1_node*do,i,1,count1_nodecp_node(i)=0*enddo!*********开始选择程序****************allselcmsel,s,masternode*do,i,1,count1_nodekk=1k=1*dowhile,kkk=nnear(slave_node(i))nsel,s,cp,,allkk=nsel(k)+0.001allselcmsel,s,masternodensel,u,node,,kcm,masternode,node*enddocp_node(i)=k*enddo!*******选择完毕****************!*******开始耦合****************allsel,all/prep7*do,i,1,count1_nodecp,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法：brickwallNRsolution!displacementload!用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!*------------------------------------loadsteps--------------------------------!*------------------------------------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!*--------------------------------2ndstep--------------------------------------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!*--------------------------------3rdstep--------------------------------------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!*--------------------------------4thstep--------------------------------------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!*allsellssolve,1,4,1FINISH也可以用弧长法，不过没有一定经验这比较难控制：!monobrickwall-1arclengthsolution!displacementload!*---------------------------------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!*--------------loadsteps----------------!*--------------firststep-----------------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!*-------------secondstep----------------esel,s,type,,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uz!*nsel,s,loc,x,-0