FLAC3D命令流

命令流

1

; --------------------

计算第一步

p4 15 15 0 p5 0 15 15 p6 15 0 10 p7 20

config

gen zone bri p0 0 0 0 pl 10 0 0 p2 0 10 0 p3 0 0 10 &

;... STATE: STATE1 ....

20 20 &

size 10 10 10 rat 1.0 0.9 1.1 group

brick_1

grid 10,10

model elastic

group 'User:Soil' notn ull

gen zone bri p0 20 0 0 p1 add 10 0 0 p2 add 0 20 0 p3

model elastic notnull group 'User:Soil'

add 0 0 15 &

size 10 10 10 rat 1.0 0.9 1.1 group

brick_2

gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1.0

0.9

1.1 group brick_3

plot sur

2-1

new

gen zon bri size 3 3 3

model elas

prop bulk 3e8 shear 1e8

ini dens 2000

fix z ran z -.1 .1

fix x ran x -.1 .1

fix x ran x 2.9 3.1

fix y ran y -.1 .1

fix y ran y 2.9 3.1

set grav 0 0 -10

solve

app n stress -10e4 ran z 3 x 1 2 y 1 2

hist gp vel 0 0 3

hist gp vel 0 3 3

plo hist 1 red

plo add hist 2 blue

solve

3-1

; ----------------------

工程信息

;Project Record Tree export

;Title:Simple test

prop

den

sity=1500.0

bulk=3E6

shear=1E6

not

null

group

'User:Soil'

fix

x y j 1

fix

x i 1

fix

x i 11

set gravity=9.81

history 999 un bala need

solve

save

; --------------------

计算第二步

;... STATE: STATE2 ....

in itial xdisp 0 ydisp 0

in itial xvel 0 yvel 0

model null i 4 7 j 8 10

group 'n ull' i 4 7 j 8 10

group delete 'n ull'

history 1 xdisp i=4, j=11

solve

save

; -------------------

绘图命令

;*** plot comma nds ****

;plot n ame: syy

plot hold grid syy fill

;plot n ame: Un bala need force

plot hold history 999

;plot n ame: grid

plot hold grid magnify 20.0 lred grid displaceme nt

;plot n ame: Xdis-A

plot hold history 1 line

5-

1

new

< br>;

定义球体半径和半径方向上单元网格数

def parm

rad=10.0 rad_size=5 end parm

建立八分之一球体外接立

方体网格

gen zone pyramid p0 rad 0 0 p1 rad 0 rad p2 rad rad 0

p3 0 0 0 &

p4 rad rad rad size rad_size rad_size rad_size group

1

gen zone

pyramid

p0 0

rad 0

p1 rad rad

0

p2

0 rad rad

p3 0 0 0 &

p4 rad rad rad size rad_size rad_size rad_size group

2

gen zone

pyramid

p0 0

0 rad

p1 0

rad

rad p2

rad 0

rad

p3 0 0 0 &

p4 rad rad rad size rad_size rad_size rad_size group

3

利用

FISH

语言将内部立方体节点调整到球面

J

def make_sphere

p_gp=gp_head

loop while p_gp#null

;

获取节点点坐标值：

P=(px,py,pz) px=gp_xpos(p_gp)

py=gp_ypos(p_gp) pz=gp_zpos(p_gp)

dist=sqrt(px*px+py*py+pz*pz)

if dist>0 then

;

节点位置调整

maxp=max(px,max(py,pz))

k=(maxp/rad)*(rad/dist) gp_xpos(p_gp)=k*px

gp_ypos(p_gp)=k*py gp_zpos(p_gp)=k*pz

end_if

p_gp=gp_next(p_gp)

end_loop

end

make_sphere gen zone ref

gen zone ref dip 90

gen zone ref dip 90 dd 90

J

;

显示球体网格

;

利用镜像生成完整球体网格

plot surf pl set back wh pl bl gr

5-

3

n

gen

zon

bri

size

1

1

2

group

soil

ran

z

1

1

group

rock

ran z 0 1 expgrid 3d

6-

1

gen zon bri size 3 3 3

model mohr

prop bu 3e6 sh 1e6 coh 10e3 fric 15

fix z ran z -.1 .1

fix x ran x -.1 .1

fix x ran x 2.9 3.1

fix y ran y -.1 .1 fix y ran y 2.9 3.1 ini dens 2000

hist unbal set grav 10 solve elastic save

6-

2

rest

ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0

app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2 solve

save

6-

3

rest

ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0

app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2

hist id=2 gp zdis 1 1 3

hist id=3 gp zdis 1 1 2

hist id=4 gp xdis 1 1 3

hist id=5 gp xdis 1 1 3

hist id=6 zone szz 1 1 3

hist id=7 zone szz 1.5 1.5 2.5

hist id=8 zone sxz 1.5 1.5 2.5

solve save

x ran x 1 fix y ran y 0 fix y ran y 1 ini dens

2000

ini szz -40e3 grad 0 0 20e3 ran z 0 2 ini syy

6-

4 rest set log on set logfile

print

zone

stress

print

gp

dis

set

log off

-20e3

grad

0

0

10e3

ran

z

0

2

ini

sxx

-20e3

grad

0 0 10e3 ran z 0 2

6-

5 rest ini xd 0 yd 0 zd 0 xv 0

yv

0

zv

0

app

nstress

-100e3

ran

z

2.9

3.1

x 1 2 y 1 2 plot set rot 20 0 30 plot con

szz ou on magf 10 plot add hist 1

set movie avi step 1 file movie

start solve movie finish

7-

1 n gen zon bri size 1 1 2 model elas

prop bulk 3e7 shear 1e7 fix z ran z 0 fix

x ran x 0 fix x ran x 1 fix y ran y 0 fix

y ran y 1 ini dens 2000 set grav 0 0 -10

solve plo con sz

7-

2 n gen zon bri size 1 1 2

model mohr

prop

bulk

3e7

shear

1e7

c

1e10

f

15

tension

1e10 fix z ran z 0

fix x ran x 0 fix x ran x 1 fix y ran y 0

fix y ran y 1 ini dens 2000 set grav 0 0

-10 solve prop bulk 3e7 shear 1e7 c 10e3

f 15 ten 0 solve plo con sz

7-

3

n

gen

zone

brick

size

1

1

2

model

mohr

prop

bulk

3e7

shear

1e7

coh

10e3

fri

15

ten

0 fix z ran z 0 fix x ran x 0 fix x ran x

1

fix

y

ran

y

0

fix

y

ran

y

1

ini

dens

2000

set grav 0 0 -10 solve elas plo con sz

7-

4 new gen zone brick size 1 1 2 model

mohr prop bulk 3e7 shear 1e7 coh 10e3 fri

15 ten 0 fix z ran z 0 fix x ran x 0 fix

set grav 0 0 -10 solve plo con sz

7-

5

n

gen zon bri size 1 1 2

model m

prop bulk 3e7 shear 1e7 c 10e10 f 15 ten 1e10 fix z ran

z 0

fix x ran x 0

fix x ran x 1

fix y ran y 0

fix y ran y 1

ini dens 2000 ran z 0 1

ini dens 1500 ran z 1 2

ini szz -35e3 grad 0 0 20e3 ran z 0 1 ini syy -17.5e3

grad 0 0 10e3 ran z 0 1 ini sxx -17.5e3 grad 0 0 10e3

ran z 0 1 ini szz -15e3 grad 0 0 15e3 ran z 1 2

ini syy -7.5e3 grad 0 0 7.5e3 ran z 1 2 ini sxx -7.5e3

grad 0 0 7.5e3 ran z 1 2

ini pp 10e3 grad 0 0 -10e3 ran z 0 1 set grav 0 0 -10

solve plo con sz

7-6

n

gen zon bri size 1 1 2

model m

prop bulk 3e7 shear 1e7 c 10e10 f 15 ten 1e10 fix z ran

z 0

fix x ran x 0

fix x ran x 1

fix y ran y 0

fix y ran y 1

ini dens 2000 ran z 0 2

ini

szz -50e3

grad

0

0

20e3

ran

z

0

1

ini

syy

-30e3

grad

0 0 10e3 ran z 0 1 ini sxx -30e3 grad 0 0 10e3 ran z 0

1

ini

pp

30e3

grad

0

0

-10e3

ran

z 0

2

app

nstress

-10e3

ran z 2

7-7 new gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3

0 0 90 & p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150

set grav 0 0 -10

solve plo con sz

&

size

6

6

10

model

elas

pro

bulk

10e10

she 10e10 ini den 2500

apply

sxx

-1e9

grad

0

0

1.1111111e7

range

x

-.1

.1

apply

sxx -1e9 grad 0 0 6.6666666e6 range x 59.9 60.1 apply

syy

-1e9

grad

0

0

8.3333333e6

range

y -.1

.1

apply

syy

-1e9 grad 0 0 8.3333333e6 range y 59.9 60.1 apply szz

-1e8 grad 0 0 8.3333333e5 ran z 0 120 set grav 0 0 -10

step 30000

ini xdisp 0 ydisp 0 zdisp 0 ini xvel 0 yvel 0 zvel 0

plo cont szz

7-8 new gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3

0 0 90 & p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150

& size 6 6 10 model elas

pro bulk 10e10 she 10e10 ini den 2500

ini sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1 ini

sxx

-1e9

grad

0

0

6.6666666e6

range

x

59.9

60.1

ini

syy

-1e9 grad 0 0 8.3333333e6 range y -.1 .1 ini syy -1e9

grad

0

0

8.3333333e6

range

y

59.9

60.1

ini

szz

-1e8

ran

z -.1 .1 fix x y z ran z -.1 .1 set grav 0 0 -10 solve

ini xdisp 0 ydisp 0 zdisp 0 ini xvel 0 yvel 0 zvel 0

plo cont szz

7-9

new

model elastic ran z 0 1 ;

激活

0 m ~ 1 m

的单元

prop

bulk 7.8e6 shear 3.0e6 ran z 0 1 ini dens 1500 ran z 0

gen zone brick p0 0 0 -50 p1 27.5 0 -50 p2 0 5 -50 p3

1 solve

按软件默认精度求解

save

0 0 -10 size 8 1 10 group clay

gen zone brick p0 27.5 0 -50 p1 100 0 -50 p2 27.5 5 -50

p3 27.5 0 -10 ratio 1.1 1 1 size 12 1 10 group clay gen

zone brick p0 0 0 -10 p1 27.5 0 -10 p2 0 5 -10 p3 0 0

0 ratio 1 1 0.8 size 8 1 4 group soil

gen zone brick p0 27.5 0 -10 p1 100 0 -10 p2 27.5 5 -10

model elastic ran z 1 2

prop bulk 7.8e6 shear 3.0e6 ran z 1 2

ini dens 1500 ran z 1 2

solve

p3 27.5 0 0 ratio 1.1 1 0.8 size 12 1 4 group soil gen

zone brick

p0 0 0 0

p1

27.5

0 0 p2

0

5 0

p3

0 0 5

p4

27.5

5 0 &

p5 0 5 5 p6 20 0 5 p7 20 5 5 size 8 1

5 group dam

fix x y z ran z -49.9 -50.1

fix x ran x -.1 .1

fix x ran x 99.9 100.1

fix y

model mohr ran z -50 0

model null ran z 0 5

prop bulk 7.8e6 shear 3.0e6 coh 10e10 tension 1e10 ran

group soil

ini dens 1500 ran group soil

prop

bulk

3.91e6

shear

1.5e6

coh

10e10

tension

1e10

ran

group clay

ini

dens

1800

ran

group

clay

set

grav

0

0

-9.8

hist

id=1

unbal solve

prop bulk 7.8e6 shear 3.0e6 coh 10e3 fric 15 ran group

soil

prop

bulk

3.91e6

shear

1.5e6

coh

20e3

fric 20

ran

group

clay solve

save

将节点位移清零

;

将节点速度清零

;

记录地基顶部中心点的沉降

;

记录路基坡脚处的沉降

ini xdis 0 ydis 0 zdis 0

ini xvel 0 yvel 0 zvel 0

hist id=4 gp xdis 27.5 0 0

;

记

hist id=2 gp zdis 0 0 0

录路基坡脚处的水平

hist id=3 gp zdis 27.5 0 0

位移

save

model elastic ran z 2 3

prop bulk 7.8e6 shear 3.0e6 ran z 2 3

ini dens 1500 ran z 2 3

solve

save

model elastic ran z 3 4

prop bulk 7.8e6 shear 3.0e6 ran z 3 4

ini dens 1500 ran z 3 4

solve

save

model elastic ran z 4 5

prop bulk 7.8e6 shear 3.0e6 ran z 4 5

ini dens 1500 ran z 4 5

solve

save

pau

;plo bl gr

;gen zone brick p0 0 0 0 p1 100 0 0 p2 0 5 0 p3

0 0 5 size

gen zone brick p0 0 0 -50

p1 27.5 0 -50 p2 0 5 -50

p3 0 0 -10

size 8 1 10 group clay

gen zone brick p0 27.5 0 -50 p1 100 0 -50 p2 27.5 5

-50 p3 27.5 0 -10 ratio 1.1 1 1 size 12 1 10 group clay

gen zone brick p0 0 0 -10

p1 27.5 0 -10 p2 0 5 -10

p3 0 0 0

ratio 1 1 0.8 size 8 1 4 group soil

gen zone brick p0 27.5 0 -10 p1 100 0 -10 p2 27.5 5

-10 p3 27.5 0 0 ratio 1.1 1 0.8 size 12 1 4 group soil

gen zone brick p0 0 0 0 p1 27.5 0 0 p2 0 5 0 p3 0 0 5

p4 27.5 5 0 &

p5 0 5 5 p6 20 0 5 p7 20 5 5 size 8 1 5 group dam

7-10

set log on

打开

log

记录

set logfile

;

设置记录文件名为：

restore

;

调用保存的文件

print gp dis range id 517 any id 533 any

输出

两个节点

的变形值

restore

print gp dis range id 517 any id 533 any restore

print gp dis range id 517 any id 533 any restore

print gp dis range id 517 any id 533 any restore

print gp dis range id 517 any id 533 any set log

off

;

关闭

log

记录

8-

1 def abc

abc = 1 + 2 * 3

abcd = 1.0 / 2.0 end abc print fish

8-

2 def abc

if aa < 0 then

abc = 0.0 else

abc = 2.0 * aa

endif end abc

8-

3

def abc

loop aa (1, 2.5) command print aa

8-

4 new gen zon bri size 3 3 3 model elastic prop

bu 3e7 sh 1e7 ini dens 2000 fix x y z ran z -.1 .1

fix x ran x -.1 .1 fix x ran x 2.9 3.1 fix y ran

y -.1 .1 fix y ran y 2.9 3.1 set grav 10 solve ini

xd 0 yd 0 zd 0 xv 0 yv 0 zv 0 save

8-

5 rest def E_modify

p_z = zone_head d_k = 704 d_n = 0.38 d_pa =

endcommand

endloop

end

101325.0 //

标准大气压

loop while p_z # null

sigma_3 = -1.0 * z_sig1(p_z)

E_new = d_k * d_pa * (sigma_3 / d_pa)

A

d_n

z_prop(p_z,'young') = E_new p_z = z_next(p_z)

endloop

end E_modify

8-6 rest table 1 name load_settlement def

add_load

p_gp = gp_near(2,1,3)

loop n (1,5)

app_load = n * (-1000e3) file_name =

'7-6_add_step' + string(n) + '.sav' command

app nstress app_load ran z 2.9 3.1 x 1 2 y 1 2

interface 2 face range cylinder end1 (0,0,-4.9) end2

solve

save

file_name endcommand

xtable(1,n) = -1.0 *

(0,0,-5.1) radius .31

app_load ytable(1,n) = gp_zdisp(p_gp) endloop

end add_load save

8-7 rest def find_max_disp p_gp = gp_head

J

pause

gen

zone

cyl

p0

(0,0,6)

p1

(.3,0,6)

p2

(0,0,1)

p3

(0,.3,6)

&

maxdisp_value = 0.0 maxdisp_gpid = 0 loop while p_gp #

p4 (.3,0,1) p5 (0,.3,1) &

null

size 3 10 6

disp_gp

=

sqrt(gp_xdisp(p_gp)

2

+

gp_ydisp(p_gp)

gen zone cyl p0 (0,0,6.1) p1 (.3,0,6.1) p2 (0,0,6) p3

A

A

2 + gp_zdisp(p_gp)

A

2)

if

disp_gp

>

maxdisp_value

maxdisp_value

=

disp_gp

maxdisp_gpid = gp_id(p_gp) endif

p_gp = gp_next(p_gp) endloop end find_max_disp

print maxdisp_value maxdisp_gpid

rest config zextra 1 def get_sigma_dif

p_z = zone_head loop while p_z # null

sigma_dif = z_sig3(p_z) - z_sig1(p_z)

z_extra(p_z,1) = sigma_dif p_z = z_next(p_z)

endloop end get_sigma_dif plot con zextra 1

9-

1

移来移去法接触面的建立

n

gen zone radcyl p0 (0,0,0) p1(8,0,0) p2 (0,0,-5) p3

(0,8,0) &

p4 (8,0,-5) p5 (0,8,-5) p6 (8,8,0) p7 (8,8,-5) & p8

(.3,0,0)

p9

(0,.3,0)

p10

(.3,0,-5)

p11

(0,.3,-5)

&

size

3 10 6 15 ratio 1 1 1 1.15 gen zone radcyl p0 (0,0,-5)

p1 (8,0,-5) p2 (0,0,-8) p3 (0,8,-5)

&

p4(8,0,-8) p5 (0,8,-8) p6 (8,8,-5) p7 (8,8,-8) & p8

(.3,0,-5) p9 (0,.3,-5) p10 (.3,0,-8) p11 (0,.3,-8) &

size 3 6 6 15 ratio 1 1 1 1.15 fill gen zone reflect dd

270 dip 90 group clay

J

interface 1 face range cylinder end1 (0,0,0) end2

(0,0,-5.1) radius .31 &

cylinder end1 (0,0,0) end2 (0,0,-5.1) radius .29 not

(0,.3,6.1) &

p4 (.3,0,6) p5 (0,.3,6) &

size 3 1 6

gen zone reflect dd 270 dip 90 range z 1 6.1

group pile range z 1 6.1 pause

ini z add -6.0 range group pile save pile_

9-

2

.

旦丰旦土洼

导来导去法

n

gen zone brick size 3 3 3 group 2 range x 1 2 y 1 2 z

1 2

group 1 range gr 2 not

save

del ran group 2 not

interface 1 face range x 1 y 1 2 z 1 2 interface 1 face

range x 2 y 1 2 z 1 2 interface 1 face range x 1 2 y 1

z 1 2 interface 1 face range x 1 2 y 1 2 z 1 interface

1 face range x 1 2 y 1 2 z 2

rest del ran group 2 expgrid 3d

impgrid 3d

model ela

fix

x

y

z

ran

z

0

ini

den

2000

set

grav

0

0

-10

interface

1 prop kn 20e6 ks 20e6 coh 10e3 fri 15 app nstr -200e3

ran x 0 1 y 1 2 z 3 solve

1 -1 0 below int 1 maxedge 0.55 int 2 maxedge 0.55

9-

3 ;

切割模型法

n

gen zone brick size 3 3 3 group 1 range x 1 2 y 1 2 z

Assign

properties

prop

shear

1e8

bulk

2e8

fric

30

range

2 3 group 2 range group 1 not gen separate 1 int 1 wrap

;

1 2

int 1 maxedge 0.5

plo int red

Material prop shear 1e8 bulk 2e8

range Bin

ini den 2000

int

1 prop

ks 2e9

kn

2e9 fric 15 int

2 prop

ks 2e9 kn 2e9 fric 15 ; Assign Boundary Conditions fix

x range x -0.1 0.1 any x 5.9 6.1 any fix y range y -0.1

0.1 any y 5.9 6.1 any fix z range z -0.1 0.1 Bin

9-

4

; Create Material Zones gen zone brick size 5 5 5 &

p0 (0,0,0) p1 (3,0,0) p2 (0,3,0) p3 (0,0,5)

&

p4

(3,3,0)

p5

(0,5,5)

p6

(5,0,5)

p7

(5,5,5)

gen

zone

brick

size 5

5

5

p0

(0,0,5)

edge

5.0

group Material Create Bin Zones gen zone

brick size 1 5 5 &

p0 (3,0,0) p1 add (3,0,0) p2 add (0,3,0) &

p3 add (2,0,5) p4 add (3,6,0) p5 add (2,5,5)

&

p6 add (3,0,5) p7 add (3,6,5) gen zone brick

size 1 5 5 &

p0 (5,0,5) p1 add (1,0,0) p2 add (0,5,0) &

p3 add (0,0,5) p4 add (1,6,0) p5 add (0,5,5)

&

p6 add (1,0,5) p7 add (1,6,5)

gen zone brick size 5 1 5 &

p0 (0,3,0) p1 add (3,0,0) p2 add (0,3,0) &

p3 add (0,2,5) p4 add (6,3,0) p5 add (0,3,5)

&

p6 add (5,2,5) p7 add (6,3,5) gen zone brick

size 5 1 5 &

p0 (0,5,5) p1 add (5,0,0) p2 add (0,1,0) &

p3 add (0,0,5) p4 add (6,1,0) p5 add (0,1,5)

& p6 add (5,0,5) p7 add (6,1,5) group Bin range group

Material not Create named range synonyms range

name=Bin

group

Bin

range

name=Material

group

Material

;

Assign

models

to

groups

model

mohr

range

Material

model

elas range Bin

gen separate Material interface 1 wrap Material Bin

range plane ori 0 0 0 normal 1 -1 0 above

interface

2

wrap

Material

Bin

range

plane

ori

0

0

0

normal

; Monitor histories

hist unbal hist gp

zdisp (6,6,10) hist

gp zdisp (0,0,10)

hist gp zdisp

(0,0,0)

; Settings set large

set grav 0,0,-10

; Cycling step 4000

save

9-5

rest pile_

model mohr range group clay

prop bulk 8.333e7 shear 3.846e7 coh 30000 fric 0 range

group clay

model elas

range group pile

prop bulk 8.333e7 shear 3.846e7 range group pile

interface 1 prop kn 1e8 ks 1e8 fric 20 coh 30000

interface 2 prop kn 1e8 ks 1e8 fric 20 coh 30000

J

ini dens 1230 range group clay

ini dens 1230 range group pile

model null range z -0.1 0.15

J

fix z range z -8.1 -7.9

fix x range x -8.1 -7.9

fix x range x 7.9 8.1

fix y range y -.1 .1

fix y range y 7.9 8.1

set grav 0 0 -10

ini szz 0. grad 0 0 12300. range z -5.5 0.

ini szz 17600 grad 0 0 15500 range z -8 -5.5

ini sxx 0. grad 0 0 5271.4 range z -5.5 0.

ini sxx 7542.86 grad 0 0 6642.86 range z -8 -5.5

ini sxx add 31428.6 grad 0 0 5714.3 range z -8 -5.5

ini syy 0. grad 0 0 5271.4 range z -5.5 0.

ini syy 7542.86 grad 0 0 6642.86 range z -8 -5.5

ini syy add 31428.6 grad 0 0 5714.3 range z -8 -5.5

rest

;

调用保存文件

water density 1000

water

table

origin

0,0,-5.5

normal

0

0

-1

ini dens 1550 range z -8 -5.5

hist unbal

solve rat

1.e-6

save

model elas prop

range

group

pile

bulk 13.9e9 ini

shear 10.4e9

range group pile

dens 2500 call

range group pile

find_

solve rat 1.e-6

save

ini state 0 ini xdis 0.0 ydis 0.0 zdis 0.0

位移清零

apply szz -0.4e6 range z 0.05 0.15 group pile

桩顶加

第一级荷载

solve

save

print gp disp range id 1

;

输出第一级荷载下的

桩顶位移，假定桩顶中心的

id

号为

1

apply szz -0.6e6 range z 0.05 0.15 group pile

桩顶加

第二级荷载

solve save print gp disp range id 1

输出第

二级荷载下的桩顶位移

依次加载，直到桩破坏

;

速度加载法

rest ini state 0 ini xdis 0 ydis

0 zdis 0

def zs_top ;

检测桩顶竖向荷载

ad = top_head zftot = 0.0

loop while ad # null gp_pnt = mem(ad+1) zf =

gp_zfunbal(gp_pnt)

zftot = zftot + zf ad = mem(ad) endloop zs_top

= zftot / 0.1414 end

fix z range z 0.05 .15 group pile

;

固定桩顶速

度，用速度来确定位移

def ramp

while_stepping if step < ncut then

udapp = float(step) * udmax / float(ncut) else

udapp = udmax endif ad = top_head loop while ad

# null

gp_pnt = mem(ad+1) gp_zvel(gp_pnt) = udapp ad =

mem(ad)

endloop

end

hist

gp

zdis

0,0,0

hist

gp

zvel

0,0,0

hist

zs_top

hist

zone

szz

0,0,-.1

set

mech

damp

comb

set

udmax = -1e-8 ncut 30000 step 225000 save

25.0, 25.0 ) & size=(1, 2, 2)

sel shell id=5 range cylinder end1=(0.0, 0.0,0.0) &

end2=(0.0,25.0,0.0) radius=24.5 not plot blo gro

;

位移控制法

def solve_steps

loop n (1,21)

save_file = string(n) + '-'

command

step 40000

save save_file

pri zone stress ran id 2381 a id 2361 a id 2341

a

;

输出桩顶网格单元的应力

endcommand

endloop

end

solve_steps

10-

1

plot ad sel geom black black cid on scale=0.03

plot ad ax

pau

delete delete all zones

sel beam beg 0 0 0 end 2 0 0 nseg 2 sel beam beg 2 0 0

end 4 0 -1 nseg 3 sel beam id=2 beg 4 0 -1 end 5 0 -2

nseg 2 plot sel geo id on nod on scale 0.04 plot ad ax

10-

2

sel node id=1 0 0 0

sel node id=2 2 0 0

sel node id=3 4 0 -1

sel node id=4 5 0 -2

sel beamsel id=1 cid=1 node 1 2 sel beamsel id=1 cid=2

node 2 3 sel beamsel id=1 cid=3 node 3 4 plot sel geo

id on nod on scale 0.04 plot ad ax

10-

3

def set_vals

ptA = 25.0 * sin( 40.0*degrad )

ptB = 25.0 * cos( 40.0*degrad )

end set_vals

gen zone cylinder p0=( 0.0, 0.0, 0.0 ) &

p1=( ptA, 0.0, ptB ) & p2=( 0.0,

25.0,

0.0

)

&

p3=(

0.0,

0.0,

25.0

)

&

p4=(

ptA,

25.0,

ptB

)

&

p5=(

0.0,