ANSYS土钉墙算例命令流.docx

近期做一复合土钉墙基坑开挖施工模拟，先做一典型土钉支护工程分析，发现变形曲线与实验值差异较大，后根据邓肯张本构，随深度调整弹性模量，结果稍好一些，但最大水平位移还不是在基坑顶部。由于土体模型使用DP本构，其弹性模量的取值虽然考虑了初始深度修正，但不能随着开挖过程中应力水平变化而变化，变形曲线与实际观测不大符合，如何在ansys中更切实的实现土的本构模型，土钉，土体，支护桩之间的接触模拟，以及高效的整体三维非线性分析，还是值得更深入钻研的问题。

请诸位高手热烈讨论，共同进步。

!3.2.1工程概况

!美国加州大学Davis分校的沈智刚(C.K.Shen)等在校园内进行了土钉支护的足尺寸原位测试试验，!测得了较为全面系统的数据资料「月。该工程场地的层主要为砂质粉土和粉质粘土，地下水位位于地表!!以下18.3m，基坑垂直开挖深为9.15m，分5步开挖，侮步开挖1.83m。土钉为钻孔注浆钉，直径为0.1m,沿!深度设置5排，设计长度6.1m，土钉与水平线倾角为200，水平间距与竖向间距均为1.93m.基坑开挖的!过程中，在距面层1.5m,4.6m,9.15m、和13.7m处设置了测斜仪量测土体的水平位移，在边坡中部第!!2,3,4层土钉上设置了应变计来量测土钉的轴力。

!先不考虑土钉滑移和面层作用

??/clear

??/prep7

??et,1,42

??et,2,solid45

??et,3,link8

???et,4,shell63

???r,1,0.01!td

???r,2,0.01

???mp,ex,1,30000!土性质

???mp,nuxy,1,0.3

???mp,dens,1,1.8

???tb,dp,1

???tbdata,1,41,25,25

??

??mp,ex,5,2.6e7!土钉性质

??mp,nuxy,5,0.3

??mp,dens,5,2.2

??mp,ex,6,3e5!面层性质

??mp,nuxy,6,0.3

??mp,dens,6,2.5

??

??tdl=6!土钉长度

??tda=20!土钉倾角

??pi=3.1415926

??sina=sin(tda/180*pi)

??jkds=18!计算总深

??jks=9!开挖深度

??jkw=9!开挖宽度

??zw=jkw+tdl*3!计算总宽

??tdhj=0.46!土钉横距

??

??

??*do,i,1,11

??k,i,(i-1)*0.9

??k,12,jkds

??*enddo!firstlinepoint

??kgen,2,all,,,,jkw!secondline

??k,25,,tdl+jkw!thirdline

*do,i,14,23

k,i+12,kx(i)+tdl*sina,tdl+jkw

*enddo

k,36,kx(24),tdl+jkw

*do,i,25,36

k,i+12,kx(i),27

*enddo

*do,j,1,3

*do,i,1,11

a,i+12*(j-1),i+1+12*(j-1),i+13+12*(j-1),i+12+12*(j-1)

*enddo

*enddo

??

lsel,s,loc,x,1

lesize,all,1

lsel,s,loc,y,

lesize,all,0.9

mshkey,1

amesh,all

allsel

EXTOPT,ESIZE,1,0,

EXTOPT,ACLEAR,0

type,2

mat,1

vext,all,,,,,tdhj

asel,,,,1,33

EXTOPT,ESIZE,1,0,

EXTOPT,ACLEAR,0

type,2

mat,1

vext,all,,,,,-tdhj

*get,knum,kp,,count

*status,knum

!38,40,44,46,50土钉线号码

lsel,s,,,35

lsel,a,,,40

lsel,a,,,44

lsel,a,,,48

lsel,a,,,52

cm,ltd,line

ksll

type,3

mat,5

real,1

lmesh,all

/prep7

asel,s,,,1,33

aclear,all

nummrg,all

/solu

nsel,s,loc,x,jkds

d,all,all

nsel,s,loc,y,

nsel,a,loc,y,zw

dsym,symm,y

nsel,s,loc,z,-tdhj

nsel