ansys在设置桩土接触时接触面怎么选取呀 [复制链接]

桩的接触面好选择，图我应该选择哪几个接触面呢

你这是万水教材上的例子吧？书上说的比较清楚了

是的，我照着操作做不出来，我用的书是《土木工程有限元分析》

最近做了一个桩土分析，做不出来，求大神们指导一下
inish
/clear
/prep7
et,1,plane42
et,2,solid45
mp,ex,1,3.8e10  
mp,prxy,1,0.2
mp,dens,1,2500

mp,ex,2,4.9e6
mp,prxy,2,0.3
mp,dens,2,1780
tb,dp,2
tbdata,1,20e3,7.8,0  

mp,ex,3,1.9e6
mp,prxy,3,0.3
mp,dens,3,1670
tb,dp,3
tbdata,1,9e3,3.7,0  

mp,ex,4,3.8e6
mp,prxy,4,0.3
mp,dens,4,1790
tb,dp,4
tbdata,1,2000,21.8,0

mp,ex,5,9.4e6
mp,prxy,5,0.3
mp,dens,5,1810
tb,dp,5
tbdata,1,41.4e3,16.5,0

mp,ex,6,4.9e6
mp,prxy,6,0.3
mp,dens,6,1800
tb,dp,6
tbdata,1,2000,23.8,0

mp,ex,7,11.5e6
mp,prxy,7,0.3
mp,dens,7,1820
tb,dp,7
tbdata,1,62.9e3,18.2,0

mp,ex,8,5.8e6
mp,prxy,8,0.3
mp,dens,8,1820
tb,dp,8
tbdata,1,2000,42,0

mp,ex,9,12.7e6
mp,prxy,9,0.3
mp,dens,9,1840
tb,dp,9
tbdata,1,77.8e3,19.3,0

rectng,0,0.5,0,-52
rectng,0.5,10,0,-2
rectng,0.5,10,-2,-12.5
rectng,0.5,10,-12.5,-14.6
rectng,0.5,10,-14.6,-21.7
rectng,0.5,10,-21.7,-24.6
rectng,0.5,10,-24.6,-37.7
rectng,0.5,10,-37.7,-42.5
rectng,0.5,10,-42.5,-52
rectng,0,0.5,-52,-78
rectng,0.5,10,-52,-78

type,2
mat,1
vrotat,1,,,,,,4,1,90,1
mat,2
vrotat,2,,,,,,4,1,90,1
mat,3
vrotat,3,,,,,,4,1,90,1
mat,4
vrotat,4,,,,,,4,1,90,1
mat,5
vrotat,5,,,,,,4,1,90,1
mat,6
vrotat,6,,,,,,4,1,90,1
mat,7
vrotat,7,,,,,,4,1,90,1
mat,8
vrotat,8,,,,,,4,1,90,1
mat,9
vrotat,9,,,,,,4,1,90,1
mat,9
vrotat,10,,,,,,4,1,90,1
mat,9
vrotat,11,,,,,,4,1,90,1
模型建立

vsel,s,,,2      
type,2
mat,2
vplot,2
lsel,s,,,54,57,1
lesize,all,,,16
lsel,s,,,5,7,2
lsel,a,,,50,52,2
lesize,all,,,10
lsel,s,,,6,8,2
lsel,a,,,51,53,2
lesize,all,0.5
vsweep,2

vsel,s,,,3      
type,2
mat,3
vplot,3
lsel,s,,,62,65,1
lesize,all,,,16
lsel,s,,,9,11,2
lsel,a,,,58,60,2
lesize,all,,,10
lsel,s,,,10,12,2
lsel,a,,,59,61,2
lesize,all,0.5
vsweep,3

vsel,s,,,4
type,2
mat,4
vplot,4
lsel,s,,,70,73,1
lesize,all,,,16
lsel,s,,,13,15,2
lsel,a,,,66,68,2
lesize,all,,,10
lsel,s,,,14,16,2
lsel,a,,,67,69,2
lesize,all,0.5
vsweep,4

vsel,s,,,5
type,2
mat,5
vplot,5
lsel,s,,,78,81,1
lesize,all,,,16
lsel,s,,,17,19,2
lsel,a,,,74,76,2
lesize,all,,,10
lsel,s,,,18,20,2
lsel,a,,,75,77,2
lesize,all,0.5
vsweep,5

vsel,s,,,6
type,2
mat,6
vplot,6
lsel,s,,,86,89,1
lesize,all,,,16
lsel,s,,,21,23,2
lsel,a,,,82,84,2
lesize,all,,,10
lsel,s,,,22,24,2
lsel,a,,,83,85,2
lesize,all,0.5
vsweep,6



vsel,s,,,7
type,2
mat,7
vplot,7
lsel,s,,,94,97,1
lesize,all,,,16
lsel,s,,,25,27,2
lsel,a,,,90,92,2
lesize,all,,,10
lsel,s,,,26,28,2
lsel,a,,,91,93,2
lesize,all,0.5
vsweep,7

vsel,s,,,8
type,2
mat,8
vplot,8
lsel,s,,,102,105,1
lesize,all,,,16
lsel,s,,,98,100,2
lsel,a,,,29,31,2
lesize,all,,,10
lsel,s,,,30,32,2
lsel,a,,,99,101,2
lesize,all,0.5
vsweep,8

vsel,s,,,9
type,2
mat,9
vplot,9
lsel,s,,,110,113,1
lesize,all,,,16
lsel,s,,,33,35,2
lsel,a,,,106,108,2
lesize,all,,,10
lsel,s,,,34,36,2
lsel,a,,,107,109,2
lesize,all,0.5
vsweep,9

vsel,s,,,10
type,2
mat,9
vplot,10
lsel,s,,,117,118,1
lesize,all,,,16
lsel,s,,,37,39,2
lsel,a,,,114,116,2
lesize,all,,,6
lsel,s,,,38,40,2
lsel,a,,,115
lesize,all,0.5
vsweep,10

vsel,s,,,11
type,2
mat,9
vplot,11
lsel,s,,,123,126,1
lesize,all,,,16
lsel,s,,,41,43,2
lsel,a,,,119,121,2
lesize,all,,,10
lsel,s,,,42,44,2
lsel,a,,,120,122,2
lesize,all,0.5
vsweep,11
划分网格

接下来merge 除了桩以外所有节点，不会命令流，直接GUI操作的

vsel,s,,,1
type,2
mat,1
vplot,1
lsel,s,,,1,3,2
lsel,a,,,12
lsel,a,,,45
lsel,a,,,47,48,1
lesize,all,,,6
lsel,s,,,2,4,2
lsel,a,,,46,1
lesize,all,0.5  
vsweep,1
划分桩单元
allsel
ET,3,170
ET,4,173
keyopt,4,9,0
keyopt,4,12,2
mp,mu,2,0.2
mp,mu,3,0.2
mp,mu,4,0.2
mp,mu,5,0.2
mp,mu,6,0.2
mp,mu,7,0.2
mp,mu,8,0.2
mp,mu,9,0.2
定义接触单元
asel,s,,,13
aplot
cm,target,area
type,3
real,10
nsla,s,1
esurf,all

asel,s,,,19
aplot
cm,contact,area
type,4
real,10
nsla,s,1
esurf,all

asel,s,,,13
aplot
cm,target,area
type,3
real,11
nsla,s,1
esurf,all

asel,s,,,24
aplot
cm,contact,area
type,4
real,11
nsla,s,1
esurf,all

asel,s,,,13
aplot
cm,target,area
type,3
real,12
nsla,s,1
esurf,all

asel,s,,,29
aplot
cm,contact,area
type,4
real,12
nsla,s,1
esurf,all

asel,s,,,13
aplot
cm,target,area
type,3
real,13
nsla,s,1
esurf,all

asel,s,,,34
aplot
cm,contact,area
type,4
real,13
nsla,s,1
esurf,all

asel,s,,,13
aplot
cm,target,area
type,3
real,14
nsla,s,1
esurf,all

asel,s,,,39
aplot
cm,contact,area
type,4
real,14
nsla,s,1
esurf,all

asel,s,,,13
aplot
cm,target,area
type,3
real,15
nsla,s,1
esurf,all

asel,s,,,44
aplot
cm,contact,area
type,4
real,15
nsla,s,1
esurf,all

asel,s,,,13
aplot
cm,target,area
type,3
real,16
nsla,s,1
esurf,all

asel,s,,,49
aplot
cm,contact,area
type,4
real,16
nsla,s,1
esurf,all

asel,s,,,13
aplot
cm,target,area
type,3
real,17
nsla,s,1
esurf,all

asel,s,,,54
aplot
cm,contact,area
type,4
real,17
nsla,s,1
esurf,all

asel,s,,,12
aplot
cm,target,area
type,3
real,18
nsla,s,1
esurf,all

asel,s,,,58
aplot
cm,contact,area
type,4
real,18
nsla,s,1
esurf,all
定义面面接触

后面加载用的GUI，对面56,60（土底面）,；17,22,27,32,37,42,47,52,61（土侧面无限远处）施加全约束。对64,11,59,10，，5,,9,50,8,45,7,40,6,35,5,30,4,25,3,20,2,1,15施加正对称约束

然后再桩顶（面14）加载2075380Pa，打开大变形选项，子步设为30，用N-R Full求解器，为什么不能算。求大神帮忙看看