马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
小弟最近搞论文学习ansys,想做一个锚杆挡墙的模型,但总是运行后出现警告说某节点位移过大,调整了很久的边界约束也没效果,如下的命令流大家能帮我看一下吗?自己的水平实在很水,希望大家别见笑,望高手指点!
finish
/CLEAR,start
/fILEname,slope
/units,si
/pnum,kp,1
/pnum,line,1
/number,0
/PLOPTS,DATE,OFF
/PREP7
*afun,deg
!定义单元类型
et,1,solid65
et,2,link8
et,3,solid45
!keyopt,1,3,2
!岩土体的材料属性
mp,ex,3,3e8
mp,nuxy,3,0.35
mp,dens,3,2200
tb,dp,3
tbdata,3,3e4,25,0
!肋柱、挡土板的属性
mp,ex,1,2.8e10
mp,nuxy,1,0.25
mp,dens,1,2500
tb,CONC,1,,,
tbdata,,0.3,1,1.96e6,22e6 !张开剪切传递系数为0.3,抗拉强度fc=1.96e6pa,抗压强度fs=22e6pa,
!锚杆属性
mp,ex,2,2E11
mp,nuxy,2,0.2
mp,dens,2,7800
!锚杆实常数
R,1,3.14*0.036*0.036/4,,
!创建边坡模型线
k,1,0,0,-0.2
k,2,0,2.5,-0.2
k,3,0,2.5,0
k,4,0,0,0
k,5,-9,0,-0.2
k,6,-9,2.5,-0.2
k,7,-9,2.5,0
k,8,-9,0,0
k,9,0,1.1,-0.6
k,10,0,1.4,-0.6
k,11,0,1.4,-0.2
k,12,0,1.1,-0.2
k,13,-9,1.1,-0.6
k,14,-9,1.4,-0.6
k,15,-9,1.4,-0.2
k,16,-9,1.1,-0.2
k,17,-5.81,1.25,6.16
k,18,-8.31,1.25,6.16
k,19,0,0,20
k,20,0,2.5,20
k,21,-9,2.5,20
k,22,-9,0,20
k,23,-4,1.25,-0.6
k,24,-6.5,1.25,-0.6
!生成挡土板体
v,1,2,3,4,5,6,7,8
!生成立柱
v,9,10,11,12,13,14,15,16
!生成土体
a,4,3,7,8,
a,19,20,21,22
a,4,3,20,19
a,8,7,21,22
a,4,8,22,19
a,3,7,21,20
va,13,14,15,16,17,18
!创建锚杆模型线
l,23,17
l,24,18
hptcreat,line,33,25,ratio,0.0001
hptcreat,line,33,26,ratio,0.9999
hptcreat,line,34,27,ratio,0.0001
hptcreat,line,34,28,ratio,0.9999
!自由网格划分单元尺寸控制
!smrtsize,5,,0.5,2
!划分立柱体
mat,1
type,1
!设定映射网格划分
mshkey,1
!设定六面体网格划分
mshape,0,3d
vmesh,2
!mopt,trans,1.5
!划分挡板体
mat,1
type,1
!设定映射网格划分
mshkey,1
!设定六面体网格划分
mshape,0,3d
vmesh,1
!mopt,trans,1.5
!划分土体
mat,3
type,3
!设定映射网格划分
mshkey,1
!设定六面体网格划分
mshape,0,3d
vmesh,3
!mopt,tran,1.5
!划分锚杆线
lsel,s,,,33,34,1
lesize,all,,,20
mat,2
type,2
real,1
lmesh,33,34,1
allsel
eplot
!锚固点附近加密划分
!krefine,17,18,1,3
!krefine,23,24,1,3
!给挡土板和土体侧面(y方向)施加约束
nsel,s,loc,y,0
nsel,a,loc,y,2.5
nplot
d,all,uy
allsel,all
!给土体(Z=20)面施加约束
nsel,s,loc,z,20
nplot
d,all,all
allsel,all
!给土体底面施加约束
asel,s,area,,16
nsla,s,1
nplot
d,all,all
allsel,all
!挡土板底面施加约束
asel,s,area,,6
nsla,s,1
nplot
d,all,ux
d,all,uy
d,all,uz
allsel,all
!给立柱底面施加约束
asel,s,area,,12
nsla,s,1
nplot
d,all,all
allsel,all
!施加重力
allsel,all
acel,9.8,0,0
allsel
进入求解器
/solu
antype,static
nsubst,100
pred,on,
nropt,full
nlgeom,on
lnsrch,on
outres,all,last
!solcontrol,1,1
cnvtol,f,0.01 !力收敛准则设定
cnvtol,u,0.05 !位移收敛准则设定
allsel,all
solve
save
finish |