运行第一载荷步时出现下面问题:
显示:the value of UY at node5 is 3.340491532E+11,It is greater than the current limit of 1000000.This generally indicates rigid body motion as a result of an unconstrained model,Verify that your model is properly contrained.
命令流如下:(到加载与约束部分)
/TITLE,linkongmian
/UNIT,SI
EXX=63e9 !岩体弹性模量
Pxy1=0.22 !岩体泊凇比
COHESN=11.65e6 !岩体粘聚力
ARGLE=38.5 !岩体参数
F_1=0.33 !岩体摩擦系数
X1=-5 !岩体模型X起点坐标
X2=5 !岩体模型X起点坐标
JAO1=74 !临空面角度
KEP1=0.4 !关键点2位置
JAO2=30 !结构面角度
CIRKP=3.033 !圆心关键点Z轴坐标
D1=1.2 !桩径
H_1=-1.1 !嵌岩深度
/prep7
ET,1,SOLID45 !定义单元
ET,2,TARGE170
ET,3,CONTA173
MP,mu,1,F_1 !结构面摩擦系数
MP,EX,1,EXX !弹性模量
MP,nuxy,1,Pxy1 !泊松比
MP,dens,1,2700 !岩体密度
tb,dp,1
tbdata,1,COHESN,ARGLE,0 !DP材料参数
block,X1,X2,0,-10,10,-10 !建块1,正等侧显示
/VIEW,1,1,1,1
vplot
WPSTYLE,,,,,,,,1 !工作平面设置
KWPAVE,5 !工作平面称动到5点
wpro,,JAO1,
wpro,,-JAO1, !旋转工作面
wpro,,-90,
wpro,,JAO1,
VSBW, 1 !工作平面分离体
VDELE, 2, , ,1 !删除体2并删除体2上的关建点
wpro,,-JAO1,
KL,15,KEP1, ,
KWPAVE,7 !工作平面移动到7点
wpro,,JAO2, !旋转
FLST,2,1,6,ORDE,1 !拾取
FITEM,2,3
VSBW,P51X,SEPO,
wpro,,-JAO2,
K, ,0,0,CIRKP,
KWPAVE, 1
CYLIND,D1/2,0,H_1,2,0,360,
VSBV, 1, 3
wpro,,,90
FLST,2,2,6,ORDE,2
FITEM,2,2
FITEM,2,4
VSBW,P51X
VDELE,3, , ,1
VDELE, 6, , ,1
SAVE !模型建立完
SMRT,2 !自由体网格划分
MSHAPE,1,3D
MSHKEY,0
CM,_Y,VOLU
VSEL, , , ,1 !选取体1并自由划分
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
CM,_Y,VOLU
VSEL, , , , 5 !选取体5并自由划分
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
WPSTYLE,,,,,,,,0 !关闭工作平面显示
allsel !定义实常数
keyopt,3,9,0 !初始渗透常数
keyopt,3,12,2 !接触模型常数
r,4 !实常数编号
real,4 !实常数系数
asel,s,,,20 !选择目标面
aplot
cm,target,area !定义目标面
type,2
mat,1
nsla,s,1
esurf,all
esel,s,type,,2
eplot
/psymb,esys,1
allsel
asel,s,,,31 !选择接触面
aplot
cm,contact,area !定义接触面
type,3
mat,1
nsla,s,1
esurf,all
esel,s,type,,3
eplot
/psymb,esys,1
allsel
gplot
save
DA,P51X,UY, !在底面加y方向位移约束
FLST,2,1,4,ORDE,1
FITEM,2,9
DL,P51X, ,UZ, !在线上加Z方向位移约束
FLST,2,2,5,ORDE,2
FITEM,2,21
FITEM,2,33
DA,P51X,SYMM !对称约束
SFGRAD,PRES, ,Y,36000,-8100, !在面上加梯度力
DA,P51X,SYMM !对称约束
SFGRAD,PRES, ,Y,36000,-8100,
FLST,2,2,5,ORDE,2
FITEM,2,2
FITEM,2,17
SFA,P51X,1,PRES,
FLST,2,1,5,ORDE,1
FITEM,2,3
SFA,P51X,1,PRES,
SFGRAD,PRES, ,Y,57675,-9613,
FLST,2,2,5,ORDE,2
FITEM,2,9
FITEM,2,30
SFA,P51X,1,PRES,
SFGRAD,PRES, ,Y,0,0,
FLST,2,2,5,ORDE,2
FITEM,2,11
FITEM,2,32
SFA,P51X,1,PRES,120000
ACEL,0,9.8,0, !第一次加载模拟岩体自重应力
然后把下面这两个功能打开了,其他的是ANSYS10.0默认设置,求解马上出现上面问题:
nlgeom,on
nropt,full
模型图如下: |