加筋土 点面接触问题 不收敛 附命令流 恳请大家帮忙

markxifeng

情况说明：一个路堤，50cm厚的砂垫层，在砂垫层中间铺设一层土工格栅，平面应变问题，对称结构，取一半进行分析，通过点面接触进行分析，具体情况可通过命令流查看。
FINI
/CLEAR
/filname,shadianceng
/PREP7
!定义所需单元类型
ET,1,plane42
ET,2,TARGE169          ! 2-D TARGET ELEMENTS
ET,3,CONTA175          ! 2-D CONTACT ELEMENTS
ET,4,LINK1  !2D杆单元
KEYOPT,1,3,2 !平面应变分析
!定义单元材料特性
!路堤土属性
MP,EX,1,3.5E7         
MP,PRXY,1,0.35
MP,DENS,1,1600
TB,DP,1
TBDATA,,15000,30
!软土地基土属性
MP,EX,2,2.6E6         
MP,PRXY,2,0.35
MP,DENS,2,1950
TB,DP,2
TBDATA,,12000,15
!砂垫层属性
MP,EX,3,4.2E7         
MP,PRXY,3,0.27
MP,DENS,3,1800
TB,DP,3
TBDATA,,1500,44.4
!塑料土工格栅单元属性
MP,EX,4,1.74e5
MP,PRXY,4,0.3  
!接触面属性
MP,MU,5,0.5
R,1            !土体
R,2,0.002      ! 筋材
R,3,,,0.1,0.1  !接触面1
R,4,,,0.1,0.1  !接触面2
!建立有限元模型
K,1,0,0
K,2,21.5,0
k,3,20.75,0.5
K,4,12.5,6
K,5,0,6
k,6,0,0.5
K,7,0,-12
K,8,50,-12
K,9,50,0
K,10,21.125,0.25
K,11,0,0.25
K,12,21.125,0.25
K,13,0,0.25
L,1,2
L,2,10
L,10,11 !L3下层砂垫层接触面
L,11,1
L,13,12  !L5上层砂垫层接触面
L,12,3
L,3,6
L,6,13
AL,1,2,3,4
AL,5,6,7,8
A,6,3,4,5
A,7,8,9,1
AGLUE,2,3
!切分面4
asel,s,,,4!选择面4
kwpave,2
wprota,,,90
asbw,all
asel,all
wpcsys,-1
AGLUE,1,5,6
numcmp,all
allsel
save
asel,,area,,1!Selects a subset of areas
aatt,3,1,1,0!设置面的单元属性 依次为mat real type esys
aesize,1,0.25!面内部的单元尺寸定义 aesize,anum,size anum 面的编号 size 单元尺寸值
mshape,0,2d!单元形状控制 用四边形单元划分网格
mshkey,1!网格类型选择 用映射网格
amesh,1

asel,,area,,2!Selects a subset of areas
aatt,3,1,1,0!设置面的单元属性 依次为mat real type esys
aesize,2,0.25!面内部的单元尺寸定义 aesize,anum,size anum 面的编号 size 单元尺寸值
mshape,0,2d!单元形状控制 用四边形单元划分网格
mshkey,1!网格类型选择 用映射网格
amesh,2
asel,,area,,3!Selects a subset of areas
aatt,1,1,1,0!设置面的单元属性 依次为mat real type esys
aesize,3,0.5!面内部的单元尺寸定义 aesize,anum,size anum 面的编号 size 单元尺寸值
mshape,0,2d!单元形状控制 用四边形单元划分网格
mshkey,1!网格类型选择 用映射网格
amesh,3
asel,,area,,4!Selects a subset of areas
aatt,2,1,1,0!设置面的单元属性 依次为mat real type esys
aesize,4,0.5!面内部的单元尺寸定义 aesize,anum,size anum 面的编号 size 单元尺寸值
mshape,0,2d!单元形状控制 用四边形单元划分网格
mshkey,1!网格类型选择 用映射网格
amesh,4
asel,,area,,5!Selects a subset of areas
aatt,2,1,1,0!设置面的单元属性 依次为mat real type esys
aesize,5,0.5!面内部的单元尺寸定义 aesize,anum,size anum 面的编号 size 单元尺寸值
mshape,0,2d!单元形状控制 用四边形单元划分网格
mshkey,1!网格类型选择 用映射网格
amesh,5

!定义土工格栅单元
n,100000,0,0.25
NGEN,84,1,100000,,,0.25,,,1,

TYPE,4   
MAT,4
REAL,2
!直接生成法生成筋材单元
*do,i,0,82,1
e,100000+i,100001+i
*enddo

!定义接触面1
esel,s,type,,4
eplot
nsle,s,1
nplot
REAL,3
TYPE,3
mat,5
ESURF                  ! GENERATE conta175 ELEMENTS
allsel
LSEL,S,LINE,,3         
NSLL,,1                ! SELECT TARGET NODES ON BIGGER CYLINDER
REAL,3
TYPE,2
ESURF                  ! GENERATE TARGE169 ELEMENTS  
allsel
!定义接触面2
esel,s,type,,4
eplot
nsle,s,1  
REAL,4
TYPE,3
mat,5
ESURF                  ! GENERATE conta175 ELEMENTS  
allsel
LSEL,S,LINE,,5         
NSLL,,1                ! SELECT TARGET NODES ON BIGGER CYLINDER
REAL,4
TYPE,2
ESURF                  ! GENERATE TARGE169 ELEMENTS  
allsel
!!!施加边界条件
nsel,s,loc,X,0
D,all,UX
nsel,all
nsel,s,loc,X,50
D,all,UX
nsel,all
nsel,s,loc,y,-12
D,all,all
nsel,all
save
allsel

/solu
allsel
acel,,9.8
autots,on                                   !自动时间步打开此选择可以让程序决定子步间荷载增量的大小及其是增加或是减小收敛速度较快(命令autots,1)打开后似乎定义的子步数不起控制作用了
nsubst,50                                   !指定荷载子步
neqit,50                                    !平衡迭代次数在每一子步中的迭代次数
!nlgeom,on                                   !指定大变形
!nropt,full                                    !设置成newton-lapnace迭代
pred,on                                     !启用非线性求解预测器
lnsrch,on                                    !启用线性搜索
!outres,all,all         
solve

markxifeng

自己顶下，呵呵

markxifeng

难道没人做这一块吗？

ycding

做这个感觉还是flac3d做起来方便，呵呵，没有算不出，只有算出来合理不合理。

51746170

发表一点浅见。
1.接触问题本来就不容易收敛。
2.我看你划分的网格，在3号材料中间部分，中间那条线上网格划分没有对齐，重新划一下，试试看行不。
祝你好运。

eva666

期待高手指点

371103215

期待高手指点

ddbb

楼主问题解决没，求指教，谢谢