关于不收敛

郎军701

我建立了一个模型，开始画粗点网格，计算可以收敛，变形也跟试验一样，可是我画细点，就出错，不能收敛，这该怎么改啊？求大侠

senja

网格细是什么意思，单方向的细，也就是细长。还是其他的？
收敛性跟网格形状有较大关系，还跟子步数有关。
另外，求解过程出现什么样的提示、警告之类的信息，浙西楼主应该发上来

郎军701

senja 发表于 2012-6-14 10:47
网格细是什么意思，单方向的细，也就是细长。还是其他的？
收敛性跟网格形状有较大关系，还跟子步数有关。
...

同一个模型，我第一个单独施加轴向位移，可以求解，第二个是单独施加切向位移，不收敛，（不知道是不是约束的问题），第三个是同时施加轴向位移和切向位移，可以求解，变形图也与试验一样

我附加了第2个单独加载切向位移的命令流（画的有点粗，你可以直接改粗细），网格画粗的话可以求解，变形也与试验一样，但画细就不收敛，真不明白，求斑竹能抽点时间帮我看看卜。非常感谢

可圈可点

senja 发表于 2012-6-14 10:47
网格细是什么意思，单方向的细，也就是细长。还是其他的？
收敛性跟网格形状有较大关系，还跟子步数有关。
...

你好，我想请教你一些网格划分的问题，就是我网格划分完了之后，计算是不收敛的，不知道怎么修改网格，请帮忙指导一下吧，新手，请多指教。。
finish
/clear
/prep7
Et,1,solid45    !承台-桩
et,2,solid186    !土体

MP,EX,1,3.5e10    !承台
MP,nuxy,1,0.17
Mp,dens,1,2500

MP,EX,2,2.10e11   !桩
MP,nuxy,2,0.3
Mp,dens,2,7850


MP,EX,3,2.6e7       !淤泥质亚粘土
MP,nuxy,3,0.27
Mp,dens,3,790
tb,dp,2
Tbdata,1,11200,21,30    !定义DP材料及材料参数

MP,EX,4,7.4e6       !基岩，参数待定
MP,nuxy,4,0.35
Mp,dens,4,1990

MP,EX,5,7.4e6       !桩基混凝土，参数待定
MP,nuxy,5,0.35
Mp,dens,5,1990


K,1,0,0
k,2,40,0
k,3,40,20
k,4,0,20

*do,i,1,3,1
  l,i,i+1
*enddo
l,4,1

a,1,2,3,4
aplot

!在GUI中拉伸成体
VEXT,1,,,,,100


!建立第一排四根桩
wpoff,11,57.5,25
wprot,0,90,0
cyl4,0,0,1,,,,49.5           !第一排一桩

wpoff,6,0,0
cyl4,0,0,1,,,,49.5       !第一排二桩

wpoff,6,0,0
cyl4,0,0,1,,,,49.5   !第一排三桩

wpoff,6,0,0
cyl4,0,0,1,,,,49.5              !第一排四桩

wpcsys,-1
allsel,all,all
vplot

!复制排桩基
allsel
vgen,6,2,5,1,,,10

!减体
vsel,all
vplot
vsel,s,,,2,25,1
CM,zhuang,VOLU
allsel
VSBV,1,zhuang
vsel,all

!划分土体
allsel
wpcsys,-1
wpoff,0,5,0
wprota,,90
allsel
vsbw,26
allsel
vplot

!建立桩
!建立第一排四根桩
CSYS,0
wpcsys,-1
wpoff,11,57.5,25
wprot,0,90,0
cyl4,0,0,1,,,,49.5           !第一排一桩

wpoff,6,0,0
cyl4,0,0,1,,,,49.5       !第一排二桩

wpoff,6,0,0
cyl4,0,0,1,,,,49.5   !第一排三桩

wpoff,6,0,0
cyl4,0,0,1,,,,49.5              !第一排四桩

wpcsys,-1
allsel,all,all
vplot

!复制排桩基
allsel
vgen,6,3,6,1,,,10
allsel

!建立承台
CSYS,0
wpcsys,-1
wpoff,9,55.5,20
block,0,22,0,3.5,0,60
allsel

!承台与体进行减运算
vsel,all
vsel,s,,,3,26,1
CM,zhuang,VOLU
allsel
VSBV,27,zhuang,,DELETE,KEEP
nummrg,all
numcmp,all
allsel

!划分网格
!选择承台
vsel,all
allsel
VSEL,S,,,27

type,1
mat,1
esize,1
mshape,1,3d
mshkey,0
vmesh,all

!选择钢管桩

vsel,all
allsel
vsel,s,,,3,26
type,1
mat,2
esize,1
mshape,1,3d
mshkey,0
vmesh,all

!划分地基


allsel
vsel,s,,,1
vplot

type,2
mat,4
esize,1.2
mshape,1,3d
mshkey,0
vmesh,all


allsel
vplot
vsel,s,,,2
vplot

type,2
mat,3
esize,2
mshape,1,3d
mshkey,0
vmesh,all


!建立摩擦单元
/prep7
Et,7,170   !刚性面
Et,8,174     !柔性面
KEYOPT,8,9,0
KEYOPT,8,12,0
r,20        !定义单元实常数
real,20
Mp,mu,2,0.32



!定义刚性接触面

allsel
vplot
VSEL,s,MAT,,2
VSEL,R,LOC,Y,8,57.5   ! 选择体
ASLV,S                     !选中依附体上的面
aplot                      !显示所选中的面
Cm,target,are             ! 由所选择的（体面线点等等）组成一个元件
Type,7                    !给随后生成的单元激活一个单元类型号
mat,2                     !指定材料编号
real,20                    ! 单元实常数属性指示器
Nsla,s,1                   !指定仅面内节点背选中
esurf,all                   !在已存在的选中单元的自由面上覆盖产生单元
ESEL,S,TYPE,,7             !选择单元类型为7的单元
EPLOT                       !显示单元网格
/PSYMB,ESYS,1               ! 显示接触单元

asel,s,loc,y,57.5
nsla,s,1                   !选择面上的节点
esln,s,1                   !选择节点上的单元
eplot                      !显示单元
esel,r,type,,7              !选择单元类型为7  单元
eplot                  !显示单元
EDELE,all                    !删除桩顶面的接触单元

allsel
vplot

!!!!!!!!!!!!!!!!!!!!定义柔性接触面!!!!!!!!!!!!!!!!!!!!

VSEL,S,mat,,3,4
ASLV,S
ASEL,R,LOC,Z,24,76
ASEL,u,LOC,y,0,8
ASEL,u,LOC,x,0,10
ASEL,u,LOC,x,30,40
ASEL,u,LOC,y,20
APLOT

Cm,contact,area
type,8
real,20
nsla,s,1
nplot
esurf,all
ESEL,S,TYPE,,8
EPLOT
/PSYMB,ESYS,1
allsel
vplot

esel,s,type,,7
eplot
allsel
esel,s,type,,8
eplot

allsel
save
finish

!加载约束
/solu
vsel,s,mat,,3,4
vplot
aslv,s,1
aplot
asel,r,loc,z,0,0
aplot
da,all,all

allsel
vsel,s,mat,,3,4
vplot
aslv,s,1
aplot
asel,r,loc,z,100
aplot
da,all,all

vsel,s,mat,,3,4
vplot
aslv,s,1
aplot
asel,r,loc,x,0,0
aplot
da,all,all

allsel
vsel,s,mat,,3,4
vplot
aslv,s,1
aplot
asel,r,loc,x,40
aplot
da,all,all


allsel
vsel,s,mat,,3,4
vplot
aslv,s,1
aplot
asel,r,loc,y,0
aplot
da,all,all
allsel
vsel,all
allsel
vplot

!!!!!!!风暴潮浪力 601kN
f=601000
wpcsys,-1
nsel,s,loc,y,40.5
nplot
*GET,N_NUM,NODE,,COUNT

f,all,fx,f*24/N_NUM

allsel,all,all
vplot


!设置分析类型控制选项求解

antype,static      !静态分析
nlgeom,on         !计入大变形
autot,on           !打开自动时间步
acel,,9.8
time,1
nsub,20
OUTRES,ALL,ALL
solve
finish

!读入最后荷载步结果
/post1
set,last
esel,s,type,1
plnsol,u,x
plnsol,s,x      !查看桩-承台结果项

esel,s,type,2   
plnsol,u,z
plnsol,s,eqv      !查看土体结果项

esel,s,type,8       !查看接触状态
etable,st1,nmisc
pletab,st1
etable,cz1,nmisc,45  !查看接触压力
pletab,cz1