【求助】接触问题不能收敛

yangheng8855

【求助】一个支挡结构的结构与土的接触问题，试过了本论给的几乎全部参考意见，如调子步数，收敛准则，非线性搜索，迭代次数等都不行。本人已经改无可改，黔驴技穷了，真心求教各位版主，高手。
命令如下：

!准备工作，定义单元类型，材料属性

fk1=10    !法向接触刚度

fk1=1     !法向接触刚度

ftoln=1    !初始渗透因子

finish

/clear

/filname,ZJTLDTQ,1

/prep7

et,1,solid45

r,1

et,2,170

et,3,174

keyopt,3,9,0

keyopt,3,12,2

r,2,,,fk1,ftoln, , ,

rmore,,,,,,5

r,3,,,fk2,ftoln, , ,

rmore,,,,,,5

mp,ex,1,1e9

mp,prxy,1,0.3

mp,dens,1,2400

mp,ex,2,1e8

mp,prxy,2,0.3

mp,dens,2,2200

mp,ex,3,3e7

mp,prxy,3,0.36

mp,dens,3,2100

mp,ex,4,3e7

mp,prxy,4,0.36

mp,dens,4,2100

mp,ex,5,2.5e10

mp,prxy,5,0.2

mp,dens,5,2500

mp,ex,6,3e9

mp,prxy,6,0.23

mp,dens,6,2300

mp,mu,7,0.4

tb,dp,1

tbdata,,1e6,35,0

tb,dp,2

tbdata,,1e8,40,0

tb,dp,3

tbdata,,2e4,30,28

tb,dp,4

tbdata,,2e4,30,27

!定义桩长

Hpile=14

!建模，基岩

k,,0,0,0

k,,0,20

k,,20,26.5

k,,28.6,32.1

k,,50.6,52.8

k,,56.2,59.3

k,,62.0,59.5

k,,64.2,59.9

k,,66.0,64.0

k,,69.8,65.9

k,,73.4,66.0

k,,77.5,64.3

k,,97.5,71.9

k,,97.5,0,0

*do,i,1,13

l,i,i+1

*enddo

l,1,14

al,all

vext,1,,,,,-5

!创建承台

block,41.2,45.2,48.2,49.8,0,-5  

numcmp,all

!创建第二层土体

k,,0,24.3

k,,20,30.8

k,,31,39.3

k,,46,52.5

k,,51.7,57.5

k,,56.2,61.7

k,,62.8,61.7

k,,64.2,65.9

k,,66.8,66.3

k,,69.7,68.5

k,,73.6,68.4

k,,76.7,66.5

k,,97.5,73

Vplot

L,2,37

L,37,38

L,38,39

L,39,33

L,35,40

*do,i,40,48

l,i,i+1

*enddo

L,49,13

Lsel,s,,,2,12

Lsel,a,,,55,69

Lsel,a,,,47,48

al,all

vext,23,,,,,-5

numcmp,all

!将实心土变成空心，创建桩

block,41.7,44.7,48.2-hpile,48.2,-2,-4     

vsbv,1,4

block,41.7,44.7,48.2-hpile,48.2,-2,-4

vsbv,3,1

block,41.7,44.7,48.2-hpile,48.2,-2,-4   

vplot

numcmp,all

!创建挡土墙

allsel

Vplot

k,,41.2,50.3,0

k,,41.5,50.3,0

k,,41.5,50.8,0

k,,41.8,50.8,0

k,,42.3,61.8,0

k,,42.9,61.8,0

k,,45.1,57,0

k,,47.1,57,0

a,36,102,103,104,105,106,107,108,109,40,35

vext,68,,,,,-5

numcmp,all   

!创建第三层土体

k,,52.7,61.3

k,,50.6,60.2

a,40,41,42,121,122,109

vext,81,,,,,-5

k,,0,28.5

k,,27.6,40

a,129,130,36,33,39,38,37

vext,89,,,,,-5

numcmp,all

!创建挡土墙墙背土体

a,107,108,109,122,121,42

vext,98,,,,,-5

numcmp,all  

!粘结为一体

vglue,1,2

vglue,3,6

vglue,3,7

vglue,2,8

Allsel

numcmp,all

!网格划分准备工作，切分

!切割基岩

wpoff,35

wprota,,,90

vsbw,all

wpoff,,,21.2

vsbw,all

wpoff,,48.2-hpile

wprota,,90

vsbw,4

wpcsys,-1

wpoff,41.7,48.2-hpile

wprota,,,90

vsbw,6

vsbw,14

wpoff,,,3

vsbw,16

vsbw,17

numcmp,all

aadd,53,95

numcmp,all

!划分网格

!基岩

vsel,s,,,13

vatt,1,1,1

esize,0.5

vsweep,all

vsel,s,,,20

vatt,1,1,1

lesize,254,,,50

lesize,292,,,26

lesize,282,,,16

lesize,5,,,16

vsweep,all

!桩

vsel,s,,,1

vatt,5,1,1

esize,0.5

vsweep,all

!基岩

vsel,s,,,16

vsel,a,,,19

vatt,1,1,1

lesize,138,,,20

lesize,271,,,20

vsweep,all

vsel,s,,,6

vatt,1,1,1

lesize,256,,,8

lesize,271,,,21

vsweep,all

!承台

vsel,s,,,8

vatt,5,1,1

esize,0.5

vsweep,all

!第二层

vsel,s,,,4

vatt,2,1,1

lesize,260,,,1

lesize,259,,,7

lesize,241,,,7

lesize,219,,,19

vsweep,all

vsel,s,,,18

vatt,2,1,1

lesize,284,,,1

lesize,280,,,2

lesize,48,,,3

lesize,59,,,6

lesize,60,,,14

lesize,61,,,12

vsweep,all

vsel,s,,,14

vsel,a,,,17

vatt,2,1,1

vsweep,all

!第三层

vsel,s,,,9

vatt,3,1,1

lesize,220,,,19

vsweep,all

!回填土

vsel,s,,,7

vatt,4,1,1

vsweep,all

!第三层

vsel,s,,,2

vatt,3,1,1

vsweep,all

!挡土墙

vsel,s,,,5

vatt,6,1,1

vsweep,all

!外围

!基岩

vsel,s,,,3

vatt,1,1,1

esize,1

vsweep,all

vsel,s,,,10

vatt,1,1,1

esize,1

vsweep,all

!二层

vsel,s,,,12

vatt,2,1,1

esize,0.5

vsweep,all

vsel,s,,,15

vatt,2,1,1

esize,0.8

vsweep,all

!三层

vsel,s,,,11

vatt,3,1,1

esize,0.5

vsweep,all

!设置接触面

!桩

allsel

asel,s,,,6

asel,a,,,54,56,2

asel,a,,,27

asel,a,,,48

nsla,s,1

type,2

real,2

mat,7

esurf,all

!承台

Allsel

asel,s,,,20,22

asel,a,,,151

nsla,s,1

type,2

real,2

mat,7

esurf,all

!挡墙

Allsel

asel,s,,,76,79

nsla,s,1

type,2

real,2

mat,7

esurf,all

!设置目标面

!基岩

allsel

asel,s,,,57,60

asel,a,,,55

nsla,s,1

type,3

real,3

esurf,all

!第二层

Allsel

asel,s,,,62,65

asel,a,,,46,47

asel,a,,,118

asel,a,,,127

asel,a,,,135

asel,a,,,140

nsla,s,1

type,3

real,3

esurf,all

!第三层

Allsel

asel,s,,,82

asel,s,,,95

nsla,s,1

type,3

real,3

esurf,all

!杂填土

Allsel

asel,s,,,88

asel,a,,,94

nsla,s,1

type,3

real,3

esurf,all

!挡墙

Allsel

asel,s,,,80

nsla,s,1

type,3

real,3

esurf,all

!施加约束

asel,a,loc,z,0

da,all,symm  

asel,s,loc,z,-5

da,all,uz

asel,s,loc,x,0

asel,a,loc,y,0

asel,a,loc,x,97.5

da,all,all   

!施加荷载

nsel,s,loc,x,44,47.5

nsel,r,loc,y,61.8

nsel,r,loc,z,-5,0

sf,all,pres,57600   

allsel

       

!历程后处理

/solu

antype,static

nlgeom,on

autots,on

acel,,9.8

nropt,full

Lnsrch,on

time,1

Eqslv,-1

nsubst,10,1000,10

Neqit,50

outres,all,all

solve

htbbzzg

      命令流开头的三个语句 (定义接触刚度、渗透因子等) 应该放到 /clear 命令的后面，否则在 /clear 之后，开头三个命令流定义的参数都被清除掉了。      最好把开头的三个语句移动到 /prep7 之后。

yangheng8855

htbbzzg 发表于 2015-6-9 16:38
命令流开头的三个语句 (定义接触刚度、渗透因子等) 应该放到 /clear 命令的后面，否则在 /clear 之后 ...

感谢，没注意到这些小细节。