ansys桩土分析命令流

/config,nres,2000000

/units,si !国际单位制

fk1=1!法向接触刚度

ftoln=0.1!初始渗透因子

/prep7

et,1,solid45

A=3.14*0.3*0.3

cylind,0,0.192,0,11,0,90 !旋喷桩上部内侧v1

cylind,0.192,0. 2,0,11,0,90!钢管v2

cylind,0.2,0.3,0,11,0,90 !旋喷桩上部外侧v3

cylind,0,0.3,11,22,0,90 !旋喷桩下部v4

cylind,0.3,1.5,0,2.6,0,90 !2.6m淤泥质粉质粘土v5

cylind,0.3,1.5,2.6,7.92,0,90 !5.32m粉土v6

cylind,0.3,1.5,7.92,16.3,0,90 !8.38m淤泥v7

cylind,0.3,1.5,16.3,18.67,0,90!2.37m粉质粘土v8

cylind,0.3,1.5, 18.67,20.47 ,0,90!1.8m粉质粘土v9

cylind,0.3,1.5, 20.47,25.3 ,0,90!4.83m粉土v10

cylind,0,0.3,22,42,0,90 !桩底内侧20m粉质粘土v11

cylind,0.3,1.5,25.3,42,0,90 !桩底外侧20m粉质粘土v12

allsel,all

vsel,s,,,2

mp,ex,2,2e11 !钢管

mp,nuxy,2,0.3

mp,dens,2,7850

allsel,all

vsel,s,,,1,4

vsel,u,,,2

mp,ex,3,1.5e8

mp,nuxy,3,0.2

mp,dens,3,2131

!tb,concr,3 !水泥土破坏模型，材料属性

!tbdata,,0.7,0.9,0.3e6,2e6!张开裂缝的剪力传递系数 ,闭合裂缝的剪力传递系数 ,单轴抗拉强度, 单轴抗压强度

allsel,all

vsel,s,,,5

vplot

mp,ex,5,3e6 !淤泥质粉质黏土

mp,nuxy,5,0.35

mp,dens,5,1950

!mp,mu,5,0.633 !设置摩擦系数为0.633

tb,dp,5 !土层5的D-P参数

tbdata,1,13.2e3,10.1,0 !土的粘聚力，摩擦角，膨胀角

allsel,all

vsel,s,,,6

vplot

mp,ex,6,9e6 !粉土

mp,nuxy,6,0.3

mp,dens,6,1970

!mp,mu,6,0.633 !设置摩擦系数为0.633

tb,dp,6 !土层6的D-P参数

tbdata,1,15.2e3,21.2,30 !土的粘聚力，摩擦角，膨胀角

allsel,all

vsel,s,,,7

vplot

mp,ex,7,1.5e6 !淤泥

mp,nuxy,7,0.4

mp,dens,7,1670

!mp,mu,6,0.633 !设置摩擦系数为0.633

tb,dp,7 !土层7的D-P参数

tbdata,1,12.9e3,9,0 !土的粘聚力，摩擦角，膨胀角

allsel,all

vsel,s,,,8

vplot

mp,ex,8,5e6 !粉质粘土

mp,nuxy,8,0.35

mp,dens,8,1880

!mp,mu,8,0.633 !设置摩擦系数为0.633

tb,dp,8 !土层8的D-P参数

tbdata,1,15.8e3,15.2,0 !土的粘聚力，摩擦角，膨胀角

allsel,all

vsel,s,,,9

vplot

mp,ex,9,8e6 !粉质粘土

mp,nuxy,9,0.35

mp,dens,9,2000

!mp,mu,9,0.633 !设置摩擦系数为0.633

tb,dp,9 !土层9的D-P参数

tbdata,1,37.1e3,19.8,0 !土的粘聚力，摩擦角，膨胀角

allsel,all

vsel,s,,,10

vplot

mp,ex,10,10e6 !粉土

mp,nuxy,10,0.3

mp,dens,10,1980

!mp,mu,10,0.633 !设置摩擦系数为0.633

tb,dp,10 !土层10的D-P参数

tbdata,1,14.3e3,20.1,10 !土的粘聚力，摩擦角，膨胀角

allsel,all

vsel,s,,,11,12

vplot

mp,ex,11,6e6 !粉质粘土

mp,nuxy,11,0.35

mp,dens,11,1910

!mp,mu,10,0.633 !设置摩擦系数为0.633

tb,dp,11 !土层11的D-P参数

tbdata,1,23.4e3,15.9,0 !土的粘聚力，摩擦角，膨胀角

!布尔运算粘结体所有的除桩之外的部分

ALLSEL,ALL

vsel,s,,,5,12,1

vplot

vglue,all

vsel,s,,,1,4

vplot

vglue,all

allsel,all

/pnum,line,1

save

allsel,all

allsel,all

lsel,s,,,3,4

lsel,a,,,13,14

lesize,all,,,6

lsel,s,,,7,8

lsel,a,,,18,19

lesize,all,0.1

lsel,s,,,55,58

lesize,all,,,1

type,1

mat,2 !钢管

vmesh,6

allsel,all

vsel,s,,,1,7,6

vsel,a,,,8

vplot

type,1

mat,3 !水泥土

lsel,s,,,1,2

lesize,all,,,6

lsel,s,,,55,56 lesize,all,,,1

lsel,s,,,22,24,2 lesize,all,,,10,0.8,1 lsel,s,,,3

lesize,all,,,6

lsel,s,,,13

lesize,13,,,6

lsel,s,,,25

lesize,25,,,6

lsel,s,,,7,9

lsel,a,,,18,19

lsel,a,,,30,31

lsel,a,,,67,69 lesize,all,0.1

lsel,s,,,4,6 !体8 lsel,a,,,14,26,12 lesize,all,,,6

lsel,s,,,57,58 lesize,all,,,1

lsel,s,,,65,66 lesize,all,,,10,0.8,1 aadd,35,31,2

aslv,s

lsla,s

FLST,2,3,4,ORDE,3 FITEM,2,6

FITEM,2,57

FITEM,2,65

LCCAT,P51X

FLST,2,3,4,ORDE,3 FITEM,2,5

FITEM,2,58

FITEM,2,66

LCCAT,P51X vsweep,1

vsweep,7

vsweep,8

allsel,all

vsel,u,,,1,6,5 vsel,u,,,7,8

vsel,u,,,18

vplot

lsel,s,,,44,46,2 lesize,all,,,6

lsel,s,,,43,45,2 lesize,all,,,5,0.8,1 lsel,s,,,51,54, lesize,all,0.1

lsel,s,,,47,49,2 lesize,all,,,6

lsel,s,,,48,50,2 lesize,all,,,5,0.8,1 type,1

mat,5

vsweep,5

vsel,s,,,13

lsel,s,,,136,139 lesize,all,0.1

lsel,s,,,59,61,2 lesize,all,,,6

lsel,s,,,60,62,2 lesize,all,,,5,0.8,1 type,1

mat,6

vsweep,13

vsel,s,,,14

lsel,s,,,140,143 lesize,all,0.1

lsel,s,,,71,73,2 lesize,all,,,6

lsel,s,,,72,74,2 lesize,all,,,5,0.8,1 type,1

mat,7

vsweep,14

vsel,s,,,15

lsel,s,,,144,147 lesize,all,0.1

lsel,s,,,83,85,2 lesize,all,,,6

lsel,s,,,84,86,2 lesize,all,,,5,0.8,1 type,1

mat,8

vsweep,15

vsel,s,,,16

lsel,s,,,148,151 lesize,all,0.1

lsel,s,,,95,97,2 lesize,all,,,6

lsel,s,,,96,98,2 lesize,all,,,5,0.8,1 type,1

mat,9

vsweep,16

vsel,s,,,19

lsel,s,,,158,163 lesize,all,0.1

lsel,s,,,107,109,2 lesize,all,,,6

lsel,s,,,108,110,2 lesize,all,,,5,0.8,1 type,1

mat,10

vsweep,19

vsel,s,,,17

lsel,s,,,152,156 lesize,all,0.1

lsel,s,,,118,128,10 lesize,all,,,6

lsel,s,,,155,157,2 lesize,all,,,5,0.8,1 type,1

mat,11

vsweep,17

vsel,s,,,18

lsel,s,,,154,158,2

lsel,a,,,159

lsel,a,,,123

lesize,all,0.1

lsel,s,,,118,120,

lesize,all,,,6

lsel,s,,,115,117

lesize,all,,,6

type,1

mat,11

vsweep,18

!用接触向导建立接触对

et,2,170

et,3,173

keyopt,3,9,0! 包括由几何形状和接触表面偏移所引起的穿透keyopt,3,12,2 !不分开接触，允许滑动，法向不分开

!桩建立目标面

r,20 ,,,fk1,ftoln, , ,

mp,mu,100,0.85

allsel,all

vsel,s,,,1,6,5

vsel,a,,,7,8

aslv,s

aplot

asel,s,,,14

asel,a,,,19,39,20

aplot

cm,target,area

type,2

real,20

mat,100

nsla,s,1

esurf,all

esel,s,type,,2

eplot

r,21,,,fk1,ftoln, , ,

mp,mu,101,0.3

allsel,all

vsel,s,,,1,6,5

vsel,a,,,7,8

aslv,s

aplot

asel,s,,,14

asel,a,,,19,39,20 aplot

cm,target,area

type,2

real,21

mat,101

nsla,s,1

esurf,all

esel,s,type,,2

eplot

r,22,,,fk1,ftoln, , , mp,mu,102,0.15 allsel,all

vsel,s,,,1,6,5

vsel,a,,,7,8

aslv,s

aplot

asel,s,,,14

asel,a,,,19,39,20 aplot

cm,target,area

type,2

real,22

mat,102

nsla,s,1

esurf,all

esel,s,type,,2

eplot

r,23,,,fk1,ftoln, , , mp,mu,103,0.3

allsel,all

vsel,s,,,1,6,5

vsel,a,,,7,8

aslv,s

aplot

asel,s,,,14

asel,a,,,19,39,20 aplot

cm,target,area

type,2

real,23

mat,103

nsla,s,1

esurf,all

esel,s,type,,2

eplot

r,24,,,fk1,ftoln, , , mp,mu,104,0.35 allsel,all

vsel,s,,,1,6,5

vsel,a,,,7,8

aslv,s

aplot

asel,s,,,14

asel,a,,,19,39,20 aplot

cm,target,area

type,2

real,24

mat,104

nsla,s,1

esurf,all

esel,s,type,,2

eplot

r,25,,,fk1,ftoln, , , mp,mu,105,0.3

allsel,all

vsel,s,,,1,6,5

vsel,a,,,7,8

aslv,s

aplot

asel,s,,,14

asel,a,,,19,39,20 aplot

cm,target,area

type,2

real,25

mat,105

nsla,s,1

esurf,all

esel,s,type,,2 eplot

!桩建立接触面allsel,all vsel,u,,,1,6,5 vsel,u,,,7,8 vplot

asel,s,,,26 asel,a,,,72 asel,a,,,76 asel,a,,,80 asel,a,,,84 asel,a,,,95 asel,a,,,59 aplot

type,3

real,20

mat,100

nsla,s,1 esurf,all

allsel,all vsel,u,,,1,6,5 vsel,u,,,7,8 vplot

asel,s,,,72 aplot

type,3

real,21

mp,mu,101,0.3 mat,101

nsla,s,1 esurf,all

allsel,all vsel,u,,,1,6,5 vsel,u,,,7,8 vplot

asel,s,,,76 aplot

type,3

real,22

mp,mu,102,0.15

mat,102

nsla,s,1

esurf,all

allsel,all

vsel,u,,,1,6,5 vsel,u,,,7,8

vplot

asel,s,,,80

aplot

type,3

real,23

mp,mu,103,0.3 mat,103

nsla,s,1

esurf,all

allsel,all

vsel,u,,,1,6,5 vsel,u,,,7,8

vplot

asel,s,,,84

aplot

type,3

real,24

mp,mu,104,0.35 mat,104

nsla,s,1

esurf,all

allsel,all

vsel,u,,,1,6,5 vsel,u,,,7,8

vplot

asel,s,,,59,95,36 aplot

type,3

real,25

mp,mu,105,0.3 mat,105

nsla,s,1

esurf,all

/psymb,esys,1

allsel

gplot

save

finish

!求解过程

/solu

!定义面约束条件

asel,s,loc,x,0 !侧面施加对称约束

da,all,symm

asel,s,loc,y,0

da,all,symm

asel,s,loc,z,42 !底面试加约束

da,all,all

allsel,all !土外层施加约束

asel,s,,,86

asel,a,,,94

asel,a,,,82

asel,a,,,78

asel,a,,,74

asel,a,,,70

asel,a,,,25

aplot

da,all,all

ALLSEL,ALL

antype,static

nlgeom,on

solcontrol,on

autot,on

pred,on

neqit,50!每个子步允许的最大迭代次数outpr,all,all

OUTRES,all,all

FirstF=840e3/A

LastF=840e3/A!考虑注浆会增加30％的部分Increment=60e3/A

*do,Force,FirstF,LastF,Increment allsel,all

vsel,s,,,1

vsel,a,,,6

vsel,a,,,7

aslv,s,

asel,s,,,1,

asel,a,,,29

asel,a,,,34

aplot

sfa,all,,pres,Force

allsel,all

nsubst,10,50,5

solve

*enddo

finish

（注：可编辑下载，若有不当之处，请指正，谢谢!）