4-3.txt

是用apdl语言编的

!在四个角点增加钢性垫片，在中心点处直接施加-2mm，压碎关闭计算完毕


/clear,nostart
/config,nres,5000
/prep7
!定义单元类型
et,1,solid65            !混凝土单元
et,2,solid45
keyopt,1,6,3
keyopt,1,5,1
*afun,deg
R,1,2,0.01,90,0,2,0.01,0,0  !整体模型，实常数
!定义材料属性
mp,ex,1,24000		!混凝土材料属性
mp,prxy,1,0.2		
tb,conc,1,1,9
tbdata,,0.35,1,3.1125,-1


mp,ex,2,2e5		!纵向受拉钢筋，钢支座材料属性
mp,prxy,2,0.25
tb,bkin,2,1,2,1
tbdata,,360,20000

blc4,,,1000,1000,100
blc4,,,100,100,-50
wpoff,1000,0,0
blc4,,,-100,100,-50
wpoff,0,1000,0
blc4,,,-100,-100,-50
wpoff,-1000,0,0
blc4,,,100,-100,-50
/view,1,1,1,1
/ANG,1,240,ZS,1 
/replot 
/REP,FAST

lsel,s,loc,z,50
lesize,all,,,2
lsel,all
lsel,s,loc,x,500
lesize,all,,,20
lsel,all
lsel,s,loc,y,500
lesize,all,,,20
vsel,s,loc,z,50
vmesh,all

type,2
mat,2
lsel,s,loc,z,-25
lesize,all,,,1
lsel,all
lsel,s,loc,z,-50
lesize,all,,,2
vsel,s,loc,z,-25
vmesh,all

nummrg,all
numcmp,all
fini

/solu
asel,s,loc,z,-50
da,all,all

nsel,all
nsel,s,loc,z,100
nsel,r,loc,y,500
nsel,r,loc,x,500
d,all,uz,-2


nlgeom,on
NSUBST,100   
OUTRES,ALL,ALL  
NEQIT,40
pred,on
CNVTOL,F, ,0.05,2
allsel
SOLVE  
fini 


/post1

PLNSOL,U,z,0,1
Esel,s,type,,1
PLNSOL,S,1,0,1			!画第一主应力

ETABLE,rebar_1,NMISC,43 	!获取rebar1平均等效应力
PLETAB,REBAR_1,AVG
ETABLE,rebar_2,NMISC,47 	!获取rebar2平均等效应力
PLETAB,REBAR_2,AVG


etable,seqv_i,nmisc,5		!获取solid单元i节点的平均等效应力
PLETAB,seqv_i,AVG

ETABLE,11,NMISC,53	!获取第一个积分点的状态
ETABLE,22,NMISC,60	!获取第二个积分点的状态
ETABLE,33,NMISC,67	!获取第三个积分点的状态
ETABLE,44,NMISC,74	!获取第四个积分点的状态
ETABLE,55,NMISC,81	!获取第五个积分点的状态
ETABLE,66,NMISC,88	!获取第六个积分点的状态
ETABLE,77,NMISC,95	!获取第七个积分点的状态
ETABLE,88,NMISC,102	!获取第八个积分点的状态
PLETAB,11,AVG
PLETAB,22,AVG
PLETAB,33,NOAV
PLETAB,44,AVG
PLETAB,55,NOAV

/DEVICE,VECTOR,1
PLCRACK,0,0			!显示所有开裂位置
PLCRACK,0,1			!显示第一开裂位置