刚建立的NALC墙板命令流 有错误 请高手指教 非常感谢

wargod

我是新手，刚建立的，里面有错误，请指教。其中一个提醒的是：发现对TB,MISO输入应力应变曲线时，第一点的切线模量应等于你给的材料的弹性模量，以后的各个点之间的斜率应不大于你给的弹性模量
finish$/clear$/prep7
et,1,solid65,,,,,,,1$et,2,link8
mp,ex,1,1780$mp,prxy,1,0.2$mp,dens,1,500
fc=5.0$ft=0.45
tb,concr,1 $tbdata,,0.5,0.95,ft,-1,,  
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.1e5$mp,prxy,2,0.3$mp,dens,1,8900   
tb,biso,2$tbdata,,210,0
pi=acos(-1)
r,1,0.25*pi*6*6
blc4,,,100,3000,600
wpoff,50  
wprota,,,90   
vsbw,all
wpcsys,-1
wpoff,50
wpoff,,,100
vsbw,all
wpoff,,,150
vsbw,all
wpoff,,,100
vsbw,all
wpoff,,,150
vsbw,all
wpcsys,-1     
wprota,,-90
wpoff,,,50
vsbw,all
wpoff,,,100
vsbw,all
wpoff,,,300
vsbw,all
wpoff,,,1050
vsbw,all
wpoff,,,1050
vsbw,all
wpoff,,,300
vsbw,all
wpoff,,,100
vsbw,all
wpcsys,-1
elemsiz=50                                                   
lsel,s,loc,x,50$lsel,r,loc,z,100                             
cm,sgj,line$latt,2,1,2$lesize,all,elemsiz
lsel,s,loc,x,50$lsel,r,loc,z,250
cm,sgjj,line$latt,2,1,2$lesize,all,elemsiz
lsel,s,loc,x,50$lsel,r,loc,z,350
cm,sgjjj,line$latt,2,1,2$lesize,all,elemsiz
lsel,s,loc,x,50$lsel,r,loc,z,500
cm,sgjjjj,line$latt,2,1,2$lesize,all,elemsiz
lsel,s,loc,x,50$lsel,r,loc,y,50                           
cm,pj,line$latt,2,1,2$lesize,all,elemsiz
lsel,s,loc,x,50$lsel,r,loc,y,150
cm,pjj,line$latt,2,1,2$lesize,all,elemsiz
lsel,s,loc,x,50$lsel,r,loc,y,450
cm,pjjj,line$latt,2,1,2$lesize,all,elemsiz
lsel,s,loc,x,50$lsel,r,loc,y,1500
cm,pjjjj,line$latt,2,1,2$lesize,all,elemsiz
lsel,s,loc,x,50$lsel,r,loc,y,2550
cm,pjjjjj,line$latt,2,1,2$lesize,all,elemsiz
lsel,s,loc,x,50$lsel,r,loc,y,2850
cm,pjjjjjj,line$latt,2,1,2$lesize,all,elemsiz
lsel,s,loc,x,50$lsel,r,loc,y,2950
cm,pjjjjjjj,line$latt,2,1,2$lesize,all,elemsiz
lsel,all
cmsel,s,sgj$cmsel,a,sgjj$cmsel,a,sgjjj$cmsel,a,sgjjjj$cmsel,a,pj$cmsel,a,pjj$cmsel,a,pjjj$cmsel,a,pjjjj$cmsel,a,pjjjjj$cmsel,a,pjjjjjj$cmsel,a,pjjjjjjj$cm,gj,line
lmesh,all$lsel,all                                                               
/eshape,1$eplot                                                                             
vatt,1,,1$mshkey,1$esize,elemsiz$vmesh,all$allsel,all
nsel,s,loc,y,0$d,all,all
acel,,9.8
allsel,all  
/solv$antype,0$nsubst,100$outres,all,all$autos,on
neqit,50
cnvtol,u,,0.015
solve
/post1$set,last$pldisp,1
esel,s,type,,2
etable,saxl,ls,1$plls,saxl,saxl
esel,s,type,,1
/device,vector,on$plcrack