刚建立的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