马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
小弟用ANSYS建立一个开挖土体的模型,划分网格的时候,拉伸一个体时突然ANSYS自动就关了,想不明白怎么回事,还请大家看看!!!
下面是命令流
/prep7
*afun,deg !定义角度单位为度
!确定基坑参数
ar_h=0.98
ar_d=7.5
h1=-1.5
h2=-4.5
h3=-16.5
sh=h1+h2+h3+ar_d
o=18.4548
l=25.831
h=5.7863
R_i=l/2/sin(o)
a=l/2-h*sin(o)/cos(o)/2+h/2/cos(o)
b=R_i-(l/2-h*sin(o)/cos(o)/2)*cos(o)/sin(o)
x1=l/2-h*sin(o)/(2*cos(o))
hkw=-0.75
!定义材料参数,三层土,一种水泥墙
et,1,solid45
et,2,solid45
et,3,solid45
!第一层
mp,ex,1,10e6
mp,nuxy,1,0.4
mp,dens,1,1957.2
tb,dp,1
tbdata,1,20.1e3,16.8,30
!第二层
mp,ex,2,10e6
mp,nuxy,2,0.4
mp,dens,2,1783.9
tb,dp,2
tbdata,1,13.2e3,11.3,30
!第三层
mp,ex,3,10e6
mp,nuxy,3,0.4
mp,dens,3,1743.12
tb,dp,3
tbdata,1,11.3e3,9.6,30
!拱
et,4,solid45
mp,ex,4,6e8
mp,prxy,4,0.25
mp,dens,4,2160
!面单元
et,5,plane42
!建立关键点
csys,0
k,1,0,0
k,3,l/2,h/2
k,5,0,b
k,6,0,b+ar_h
k,7,a,0
k,8,a+ar_h,0
k,9,0,3*b
k,10,3*a,3*b
k,11,3*a,0
k,13,a+5*ar_h,0
k,14,a+5*ar_h,b+5*ar_h
k,15,0,b+5*ar_h
k,16,3*a,b+5*ar_h
k,17,a+5*ar_h,3*b
LOCAL,11,1,0,-(R_i-b),0,,,,,,
csys,11
k,12,R_i+ar_h,90-o
LOCAL,12,1,x1,0,0,,,,,,
csys,12
l,3,7
l,12,8
csys,11
l,5,3
l,6,12
csys,0
l,1,5
l,1,7
l,7,8
l,8,13
l,13,11
l,5,6
l,6,15
l,15,9
l,13,14
l,15,14
l,11,16
l,16,10
l,9,17
l,16,10
l,14,16
l,14,17
l,17,10
l,3,12
al,1,7,2,21
al,3,21,4,10
al,6,1,3,5
al,13,14,11,4,2,8
al,15,18,13,9
al,16,20,19,18
al,19,17,12,14
!**************开始拉伸成体****************
vext,3,,,,,hkw
allsel
asel,s,loc,z,hkw
vext,all,,,,,hkw
allsel
asel,s,loc,z,2*hkw
vext,all,,,,,hkw
allsel
asel,s,loc,z,3*hkw
vext,all,,,,,hkw
allsel
asel,s,loc,z,4*hkw
vext,all,,,,,hkw
allsel
asel,s,loc,z,5*hkw
vext,all,,,,,hkw
allsel
asel,s,loc,z,6*hkw
vext,all,,,,,hkw
allsel
asel,s,loc,z,7*hkw
vext,all,,,,,hkw
allsel
asel,s,loc,z,8*hkw
vext,all,,,,,hkw
allsel
asel,s,loc,z,9*hkw
vext,all,,,,,hkw
allsel
vext,4,7,,,,h1
allsel
vext,58,,,,,h2
vext,65,,,,,h2
vext,73,,,,,h2
vext,69,,,,,h2
allsel
asel,s,loc,z,h1+h2
vext,all,,,,,h3
allsel
vext,1,2,,,,-ar_d
allsel
asel,s,loc,z,-ar_d
vext,all,,,,,-(-h1-h2-h3-ar_d)
allsel
vplot
!**********下面要对属于同一种材料的体粘贴*************
!粘贴大圆和小圆
vglue,24,25
vglue,11,12,13,14
vglue,15,16,17,18
vglue,20,21,22,23,9,10,26,27,28
vglue,1,2
vglue,3,4,5,6,7,8
vsel,all
numcmp,all !压缩编号,(粘贴操作就有这个毛病,会自动对粘贴后的体重新编号)
vsel,all
allsel
!*****************下面用扫略的方法划分网格**************
!*********内部土体的扫略********
allsel
type,1
mat,1
lesize,1,,,5,,,,,1
lesize,3,,,15,,,,,1
lesize,5,,,15,,,,,1
lesize,6,,,15,,,,,1
extopt,esize,2 !设定扫略路线的分段数
vsweep,1,3,8 !扫略体1,坑外桩设源面为a3,目标面为a8
allsel
type,1
mat,1
lesize,23,,,5,,,,,1
lesize,24,,,15,,,,,1
lesize,25,,,15,,,,,1
lesize,22,,,15,,,,,1
extopt,esize,2 !设定扫略路线的分段数
vsweep,2,8,13 !扫略体2,坑外桩设源面为a8,目标面为a13
allsel
type,2
mat,2
lesize,31,,,5,,,,,1
lesize,32,,,15,,,,,1
lesize,33,,,15,,,,,1
lesize,30,,,15,,,,,1
extopt,esize,2 !设定扫略路线的分段数
vsweep,3,13,18 !扫略体3,坑外桩设源面为a13,目标面为a18
allsel
type,2
mat,2
lesize,39,,,5,,,,,1
lesize,40,,,15,,,,,1
lesize,41,,,15,,,,,1
lesize,38,,,15,,,,,1
extopt,esize,2 !设定扫略路线的分段数
vsweep,4,18,23 !扫略体4,坑外桩设源面为a18,目标面为a23
allsel
type,2
mat,2
lesize,47,,,5,,,,,1
lesize,48,,,15,,,,,1
lesize,49,,,15,,,,,1
lesize,46,,,15,,,,,1
extopt,esize,2 !设定扫略路线的分段数
vsweep,5,23,28 !扫略体5,坑外桩设源面为a23,目标面为a28
allsel
type,2
mat,2
lesize,55,,,5,,,,,1
lesize,56,,,15,,,,,1
lesize,57,,,15,,,,,1
lesize,54,,,15,,,,,1
extopt,esize,2 !设定扫略路线的分段数
vsweep,6,28,33 !扫略体6,坑外桩设源面为a28,目标面为a33
allsel
type,2
mat,2
lesize,63,,,5,,,,,1
lesize,64,,,15,,,,,1
lesize,65,,,15,,,,,1
lesize,62,,,15,,,,,1
extopt,esize,2 !设定扫略路线的分段数
vsweep,7,33,38 !扫略体7,坑外桩设源面为a33,目标面为a38
!*******此处扫略时程序就自动关掉,还请大家看看怎么回事******
!allsel
!type,2
!mat,2
!lesize,71,,,5,,,,,1
!lesize,72,,,15,,,,,1
!lesize,73,,,15,,,,,1
!lesize,70,,,15,,,,,1
!extopt,esize,2 !设定扫略路线的分段数
!vsweep,8,38,43 !扫略体8,坑外桩设源面为a38,目标面为a43 |