本帖最后由 tita897 于 2010-12-6 23:57 编辑
finish
/clear
A=3.1415926*0.5*0.5 !桩的横截面积
r1=0.5 !桩半径
r2=10 !土模型半径
hp1=4.8 !土层2素填土+吹填砂
hp2=3.7 !土层3冲积粘土
hp3=3.3 !土层4冲积粉砂
hp4=2.4 !岩层5强风化砂岩
hp5=3.3 !岩层6中风化砂岩
hp6=26.7 !岩层6中风化砂岩
hp7=2.5 !桩出土长度
hp8=17.5 !桩入土深度
/prep7
et,1,solid45
mp,ex,1,3.15e10 !桩体材料属性
mp,prxy,1,0.2
mp,dens,1,2500
mp,mu,1,0.2
mp,ex,2,2.8e7 !土层2材料属性 1
mp,prxy,2,0.3
mp,dens,2,1900
mp,ex,3,2.5e7 !土层3材料属性 2
mp,prxy,3,0.4
mp,dens,3,1860
mp,ex,4,3.0e7 !土层4材料属性 3
mp,prxy,4,0.3
mp,dens,4,1950
mp,ex,5,1.0e9 !岩层5材料属性 4
mp,prxy,5,0.35
mp,dens,5,2200
mp,ex,6,1.2e10 !岩层6材料属性 5
mp,prxy,6,0.25
mp,dens,6,2600
tb,dp,2 !土层2的D-P参数 1
tbdata,1,6000,22,0
tb,dp,3 !土层3的D-P参数 2
tbdata,1,13600,18,0
tb,dp,4 !土层4的D-P参数 3
tbdata,1,8000,22,0
tb,dp,5 !岩层5的D-P参数 4
tbdata,1,80000,28,0
tb,dp,6 !岩层6的D-P参数 5
tbdata,1,200000,32,0
cylind,0,r1,hp7,-hp8,0,90 !创建桩体
cylind,r1,r2,0,-hp1,0,90 !创建土层2
cylind,r1,r2,-hp1,-hp1-hp2,0,90 !创建土层3
cylind,r1,r2,-hp1-hp2,-hp1-hp2-hp3,0,90 !创建土层4
cylind,r1,r2,-hp1-hp2-hp3,-hp1-hp2-hp3-hp4,0,90 !创建岩层5
cylind,r1,r2,-hp1-hp2-hp3-hp4,-hp1-hp2-hp3-hp4-hp5,0,90 !创建桩进入中分化
cylind,0,r1,-hp1-hp2-hp3-hp4-hp5,-hp1-hp2-hp3-hp4-hp5-hp6,0,90 !建立桩尖岩层几何体
cylind,r1,r2,-hp1-hp2-hp3-hp4-hp5,-hp1-hp2-hp3-hp4-hp5-hp6,0,90 !建立岩层几何体
!布尔运算粘结体所有的除桩之外的部分
vsel,s,,,2,8,1
vplot
vglue,2,3,4,5,6,7,8
allsel,all
numcmp,all
/pnum,volu,1
vplot
wpcsys,,0
vsbw,1
kwpave,8
vsbw,9
kwpave,16
vsbw,1
kwpave,20
vsbw,9
kwpave,24
vsbw,1
wpcsys,,0
!布尔运算粘结体所有的除桩的部分
vsel,s,,,9,14,1
vplot
vglue,9,10,11,12,13,14
allsel,all
numcmp,all
/pnum,volu,1
vplot
!!!!!!!!!!!!!!划分网格
!!!!!!!!!!!桩最上层
allse
vsel,s,,,9
vplot
type,1
mat,1
lesize,4,,,6
lesize,5,,,6
lesize,6,,,6
lesize,73,,,8
vmesh,all
!!!!!!!!!!!桩1
allse
vsel,s,,,10
vplot
type,1
mat,1
lesize,81,,,12
vmesh,all
!!!!!!!!!!!桩2
allse
vsel,s,,,11
vplot
type,1
mat,1
lesize,87,,,10
vmesh,all
!!!!!!!!!!!桩3
allse
vsel,s,,,12
vplot
type,1
mat,1
lesize,90,,,6
vmesh,all
!!!!!!!!!!!桩4
allse
vsel,s,,,13
vplot
type,1
mat,1
lesize,96,,,6
vmesh,all
!!!!!!!!!!!桩5
allse
vsel,s,,,8
vplot
type,1
mat,1
lesize,92,,,6
vmesh,all
!!!!!!!!!!!!!!!划分土体1
alls
vsel,s,,,1
vplot
type,1
mat,2
lesize,16,,,6
lesize,19,,,12
lesize,17,,,20,0.4
vsweep,all
!!!!!!!!!!!!!!!划分土体2
alls
vsel,s,,,3
vplot
type,1
mat,3
lesize,47,,,10
vsweep,all
!!!!!!!!!!!!!!!划分土体3
alls
vsel,s,,,4
vplot
type,1
mat,4
lesize,52,,,6
vsweep,all
!!!!!!!!!!!!!!!划分土体4
alls
vsel,s,,,5
vplot
type,1
mat,5
lesize,56,,,6
vsweep,all
!!!!!!!!!!!!!!!划分土体5
alls
vsel,s,,,6
vplot
type,1
mat,6
lesize,60,,,6
vsweep,all
!!!!!!!!!!!!!!!划分土体6
alls
vsel,s,,,7
vplot
type,1
mat,6
lesize,46,,,26
vsweep,all
!!!!!!!!!!!桩低岩层4
allse
vsel,s,,,2
vplot
type,1
mat,6
lesize,43,,,6
lesize,44,,,6
lesize,45,,,26
vsweep,all
!用接触向导建立接触对
et,2,170
et,3,173
keyopt,3,9,0 !消除初始渗透
keyopt,3,12,2 !不分开接触,允许滑动,法向不分开
r,7
mp,mu,7,0.2
!桩侧,建立目标面
allsel,all
vsel,s,,,10,13,1
vsel,a,,,8
vplot
aslv,s
asel,s,,,49,53,4
asel,a,,,56,60,4
asel,a,,,57
asel,a,,,1
/pnum,area,1
aplot
cm,target1,area
type,2
real,7
mat,7
nsla,s,1
esurf,all
esel,s,type,,2
eplot
/psymb,esys,1
allsel
! 桩侧,建立接触面
allsel
vsel,s,,,3,6,1
vsel,a,,,1,2,1
vplot
aslv,s
asel,s,,,8
asel,a,,,16
asel,a,,,21,33,4
aplot
cm,contact1,area
type,3
real,7
mat,7
nsla,s,1
esurf,all
esel,s,type,,3
eplot
/psymb,esys,1
allsel
!求解过程
/solu
antype,static
!定义面约束条件
asel,s,loc,x,0 !侧面施加对称约束
asel,a,loc,y,0
da,all,symm
aplot
allsel
asel,s,loc,z,-hp8-hp6 !底面试加约束
da,all,uz
aplot
allsel,all !土外层施加约束
asel,s,,,7
asel,a,,,23,35,4
asel,a,,,38
aplot
da,all,all
allsel,all
acel,0,0,9.8 !定义重力加速度
!载荷步,加x方向的水平荷载
allsel
nsel,s,loc,z,hp7
f,all,fx,24e3
查接触状态
cncheck,detail
cncheck,summary
nlgeom,on ! 打开大变形选项
autot,on ! 打开自动时间步
time,1
nsubst,100
outres,all,all
allsel
solve
save
finish
请高手给我看看我的命令流。我的模型是大直径桩在水平荷载作用下的受力性能。受力点在桩顶,我在加第一级荷载24KN,计算的桩发生的水平位移非常小,和我实测的数据相差十倍,但是土的位移比较符合我实测的桩体位移。我不知道为什么加这么小的水平力,桩和土的位移怎么相差这么大,我的水平荷载要加七个等级。这才是第一级,桩土不应该相差那么大。我认为是我的接触有问题,还有就是土体的泊松比和弹性模量。
希望得到各位的指教。 |