基于ANSYS的钢筋混凝土粘结滑移建模命令流

finish $/clear $/prep7
et,1,solid65
$et,2,link8 $et,3,combin39
keyopt,1,7,1 $keyopt,3,3,3
mp,ex,1,3.25e10 $mp,prxy,1,0.2 $mp,dens,1,2600

tb,concr,1 $tbdata,,0.5,0.95,1.71e6,19.1e6
mp,ex,2,2.1e11 $mp,prxy,2,0.3 $mp,dens,2,7850

r,1
r,2,0.25*acos(-1)*0.03*0.03
r,3,0,0,0.5E-5,4497.60,1.0E-5,7126.36
rmore,1.5E-5,8567.28,2.0E-5,9324.77,2.5E-5,9726.64

rmore,3.0E-5,9924.14,3.5E-5,9891.92,4.0E-5,9428.07
rmore,4.5E-5,8154.07,5.0E-5,5514.84

blc4,,,0.2,0.3,4.0
wpoff,,0.03 $wprota,,90
vsbw,all $wpoff,0.03 $wprota,,,90
vsbw,all $wpoff,,,0.14 $vsbw,all
wpcsys,-1
lsel,s,length,,4.0 $lesize,all,,,4.0/0.05
lsel,all $
esize,0.05
$mshkey,1
vmesh,all

lsel,none
$k,50,0.03,0.03
k,51,0.03,0.03,4.0 $
k,52,0.17,0.03 $
k,53,0.17,0.03,4.0

l,50,51 $
l,52,53
lesize,all,,,4.0/0.05
latt,2,2,2

lmesh,all
cpintf,ux !创建节点的耦合自由度UX
*get,cpmax,cp,,max !获取藕合组总数

type,3
$real,3
*do,i,1,cpmax
*get,n1,cp,i,term,1,node !获取第i组的第一个节点号
*get,n2,cp,i,term,2,node !获取第i组的第二个节点号
e,n1,n2
*enddo
cpintf,uy
lsel,s,loc,z,0
$lsel,r,loc,y,0
dl,all,,uy $
lsel,s,loc,z,4.0

lsel,r,loc,y,0
$dl,all,,uy

lsel,all $
dk,5,ux,,,,uz $
dk,2,ux
asel,s,loc,y,0.3 $sfa,all,1,pres,5e4

allsel,all
/solu $
outres,all,all $
nsubst,500,,100

cnvtol,u,,0.05
solve