ANSYS实例命令流程序_接触分析例题_岩土方面
本帖最后由 Chelsea 于 2011-3-13 07:50 编辑例子1是为了获得初始应力，计算可以完成。命令流为：!以便获得初始应力，收敛！！考虑土与岩石之间的摩擦，点面接触
fini !考虑土体的相互作用。自由划分网格，局部细化,不收敛！！！！！！！！！！！
/cle
/filname,diji11
/title,地基应力、应变的分布规律
/units,si
/PREP7
!ET,1,PLANE2
ET,1,PLANE42
KEYOPT,1,3,2
KEYOPT,1,5,0
KEYOPT,1,6,0
MP,EX,1,8e9 !岩石性质
MP,PRXY,1,0.29
MP,DENS,1,2600
MP,EX,2,5.0e7 !风化槽土体性质?
MP,PRXY,2,0.25
TB,dp,2
TBDATA,1,15e3,20,0
MP,DENS,2,2000
et,2,contac48,,,1
mp,mu,3,2.5
r,1,5.0e7,5.0e6,0.01,1.0,0.001,0.001
r,2,5.0e7,5.0e6,0.01,1.0,0.001,0.001
K,1,0,0,0,
K,2,2.5,-12,0,
K,3,6.5,-12,0,
K,4,13.5,-26.5,0,
K,5,20,-30,0,
K,6,24.5,-27,0,
K,7,31.5,-12,0
K,8,35.5,-12,0,
K,9,38,0,0,
k,10,-40,0,0,
k,11,-40,-60,0,
k,12,80,-60,0,
k,13,80,0,0,
K,14,6.5,-12,0,
K,15,13.5,-26.5,0,
K,16,20,-30,0,
K,17,24.5,-27,0,
K,18,31.5,-12,0
*do,i,1,8
lstr,i,i+1
*enddo
LSTR,1,10
*do,i,10,12
lstr,i,i+1
*enddo
LSTR,13,9
*do,i,14,17
lstr,i,i+1
*enddo
LARC,14,18,16,14.57
FLST,2,13,4
*do,i,1,13
FITEM,2,i
*enddo
AL,P51X
FLST,2,5,4
*do,i,14,18
FITEM,2,i
*enddo
AL,P51X
asel,s,area,,1
aatt,1,,1
asel,s,area,,2
aatt,2,,1
esize,4
asel,all
amesh,all
lsel,s,line,,3,4,1
nsll,s,1
CM,_TARGET1,node
lsel,s,line,,14,15,1
nsll,s,1
CM,_CONTACT1,node
lsel,s,line,,5,6,1
nsll,s,1
CM,_TARGET2,node
lsel,s,line,,16,17,1
nsll,s,1
CM,_CONTACT2,node
TYPE,2
mat,3
real,1
gcgen,_CONTACT1,_TARGET1
real,2
gcgen,_CONTACT2,_TARGET2
fini
FINISH
/SOLU
antype,static
nlgeom,on
nropt,full,,off
neqit,40
pred,on
autots,on
time,1
nsubst,100,1000,10
alls
LSEL,s,line,,10,12,2
CM,_Y,LINE
CMSEL,,_Y
DL,_Y, ,ux,0,0
LSEL,s,line,,11
CM,_Y,LINE
CMSEL,,_Y
DL,_Y, ,uy,0,0
alls
acel,,9.8
iswrite,on
lnsrch,on
outres,all,all
solve
例2：!读入初始应力，收敛！！考虑土和混凝土与岩石之间的摩擦，点面接触
fini !考虑土体和混凝土之间的相互作用。不收敛！！！！！！！！！！！
/cle
/filname,diji
/title,地基应力、应变的分布规律
/units,si
/PREP7
ET,1,PLANE42
KEYOPT,1,3,2
KEYOPT,1,5,0
KEYOPT,1,6,0
MP,EX,1,8e9 !岩石性质
MP,PRXY,1,0.29
MP,DENS,1,2600
MP,EX,2,5.0e7 !风化槽土体性质?
MP,PRXY,2,0.25
TB,dp,2
TBDATA,1,15e3,20,0
MP,DENS,2,2000
MP,EX,3,2.55e10 !混凝土性质
MP,PRXY,3,0.167
MP,DENS,3,2400
et,2,contac48,,,1
mp,mu,4,2.5
mp,mu,5,4.5
mp,mu,6,3.0
r,1,5.0e7,5.0e6,0.01,1.0,0.001,0.001
r,2,5.0e7,5.0e6,0.01,1.0,0.001,0.001
r,3,8.0e9,8.0e8,0.01,1.0,0.001,0.001
r,4,8.0e9,8.0e8,0.01,1.0,0.001,0.001
r,5,5.0e7,5.0e6,0.01,1.0,0.001,0.001
K,1,0,0,0,
K,2,2.5,-12,0,
K,3,6.5,-12,0,
K,4,13.5,-26.5,0,
K,5,20,-30,0,
K,6,24.5,-27,0,
K,7,31.5,-12,0
K,8,35.5,-12,0,
K,9,38,0,0,
k,10,-40,0,0,
k,11,-40,-60,0,
k,12,80,-60,0,
k,13,80,0,0,
K,14,6.5,-12,0,
K,15,13.5,-26.5,0,
K,16,20,-30,0,
K,17,24.5,-27,0,
K,18,31.5,-12,0
K,19,0,0,0,
K,20,2.5,-12,0,
K,21,6.5,-12,0,
K,22,31.5,-12,0
K,23,35.5,-12,0,
K,24,38,0,0,
K,25,20,0,0,
*do,i,1,8
lstr,i,i+1
*enddo
LSTR,1,10
*do,i,10,12
lstr,i,i+1
*enddo
LSTR,13,9
*do,i,14,17
lstr,i,i+1
*enddo
LARC,14,18,16,14.57
*do,i,19,20
lstr,i,i+1
*enddo
LARC,21,22,16,14.57
*do,i,22,24
lstr,i,i+1
*enddo
LSTR,25,19
FLST,2,13,4
*do,i,1,13
FITEM,2,i
*enddo
AL,P51X
FLST,2,5,4
*do,i,14,18
FITEM,2,i
*enddo
AL,P51X
FLST,2,7,4
*do,i,19,25
FITEM,2,i
*enddo
AL,P51X
asel,s,area,,1
aatt,1,,1
asel,s,area,,2
aatt,2,,1
esize,4
asel,s,area,,1,2
amesh,all
lsel,s,line,,3,4,1
nsll,s,1
CM,_TARGET1,node
lsel,s,line,,14,15,1
nsll,s,1
CM,_CONTACT1,node
lsel,s,line,,5,6,1
nsll,s,1
CM,_TARGET2,node
lsel,s,line,,16,17,1
nsll,s,1
CM,_CONTACT2,node
TYPE,2
mat,4
real,1
gcgen,_CONTACT1,_TARGET1
real,2
gcgen,_CONTACT2,_TARGET2
asel,s,area,,3
aatt,3,,1
esize,4
amesh,all
lsel,s,line,,1,8,7
nsll,s,1
CM,_TARGET3,node
lsel,s,line,,19,23,4
nsll,s,1
CM,_CONTACT3,node
lsel,s,line,,2,7,5
nsll,s,1
CM,_TARGET4,node
lsel,s,line,,20,22,2
nsll,s,1
CM,_CONTACT4,node
TYPE,2
mat,5
real,3
gcgen,_CONTACT3,_TARGET3
real,4
gcgen,_CONTACT4,_TARGET4
lsel,s,line,,18
nsll,s,1
CM,_TARGET5,node
lsel,s,line,,21
nsll,s,1
CM,_CONTACT5,node
TYPE,2
mat,6
real,5
gcgen,_CONTACT5,_TARGET5
fini
/SOLU
antype,static
nlgeom,on
nropt,full,,off
neqit,40
pred,on
autots,on
time,1
esel,s,,,1,451
nsle,s
isfile,read,diji11,ist,,1
nsubst,100,1000,10
LSEL,s,line,,10,12,2
CM,_Y,LINE
CMSEL,,_Y
DL,_Y, ,ux,0,0
LSEL,s,line,,11
CM,_Y,LINE
CMSEL,,_Y
DL,_Y, ,uy,0,0
acel,,9.8
alls
LSEL,s,line,,25
CM,_Y1,LINE
CMSEL,,_Y1
SFL,_Y1,PRES,4.5e5,6e5,
LSEL,s,line,,24
CM,_Y1,LINE
CMSEL,,_Y1
SFL,_Y1,PRES,4.5e5,4.5e5
alls
iswrite,on
lnsrch,on
outres,all,all
solve 支持一下，看了你好多帖子，希望以后可以看到你更多的帖子 太感谢了，呵呵 回复 2 # Chelsea 的帖子
楼主你好，我在做一个柔柔接触的时候，分析结果又很大的刚体位移，不知如何解决，能请您帮助看看吧
我固定筒中的内核的一端，然后在内核的另一端施加位移荷载。几个荷载步运行后发现DMX很大，可是我是假的位移荷载只有10个单位，查阅资料说是刚体位移，请问这个问题如何解决啊。
不知道是不是我的命令流有问题啊？
页:
[1]