边坡计算实例.docx

边坡计算实例.docx

《边坡计算实例.docx》由会员分享，可在线阅读，更多相关《边坡计算实例.docx（9页珍藏版）》请在冰豆网上搜索。

边坡计算实例

边坡计算实例

genzoneradcyl&

p00025p1add8000p2add00-25p3add0800&

dim24.524.51212ratio1111.1

genzoneradcyl&

p0000p1add8000p2add00-15p3add0800&

dim12121212ratio1111.1fillsize571010

genzonebrick&

p00-40-15p1add1200p2add0400p3add0015&

size5127

genzonebrick&

p012-40-15p1add6800p2add0400p3add0015&

size10127ratio1.111

genzonebrick&

p012-400p1add6800p2add0400p3add12.5025&

p4add68400p5add12.54025p6add68025p7add684025&

size101210rat1.111

;installwatertable

inipp0grad00-1e4rangez-15.10

;callwater.fis

;assignMohrCoulombmodelandproperties

modmo

probulk2e8she1e8fric45coh1e6ten1e6

;boundaryconditions

fixxrangex-.1.1

fixxrangex79.980.1

fixyrangey-40.1-39.9

fixyrangey79.980.1

fixxrangez-15.1-14.9

fixyrangez-15.1-14.9

fixzrangez-15.1-14.9

;initializedensity:

saturatedandunsaturated

definidens

pnt=zonehead

loopwhilepnt#null

ifzpp（pnt）#0.0then

zdensity（pnt）=2600

else

zdensity（pnt）=2500

endif

pnt=znext（pnt）

endloop

end

inidens

setg00-10

;initializegravitystresses

inisyy-6.25e5grad002.5e4

inisxx-6.25e5grad002.5e4

iniszz-6.25e5grad002.5e4

;histories

histunbal

histgpxdisp24.5-4025

histgpzdisp24.5-4025

histgpxdisp24.5025

histgpydisp24.5025

histgpzdisp24.5025

histgpydisp024.525

histgpzdisp024.525

;equilibrate

solve

savesleq.sav

;reducestrength

procoh3.7e4

inixd0

iniyd0

inizd0

inixv0

iniyv0

inizv0

step4000

saveslfail.sav

;createplotviewsofresults

;view:

slgeom（modelgrid）

plotcreateslgeom

setoutslgeom.plt

plotsetrot20050

plotsetang20

plotsetcent40205

plotaddsurfacelcyan

plotaddaxesbrown

plothard

;view:

slwater（watertablewithmodelgrid）

plotcreateslwater

setoutslwater.plt

plotsetrot20050

plotsetang20

plotsetcent40205

plotaddwaterlgra

plotaddsketchblack

plothard

;view:

slppres（porepressurecontours）

plotcreateslppres

setoutslppres.plt

plotsetrot20050

plotsetang20

plotsetcent40205

plotaddcontppint2.5e4shon

plothard

;view:

sldisp（displacementcontours）

plotcreatesldisp

setoutsldisp.plt

plotsetrot20050

plotsetang20

plotsetcent40205

plotaddcontdispshon

plothard

;view:

slhist（displacementhistories）

plotcreateslhist

setoutslhist.plt

plotaddhist2378

plothard

ret

;DataFile"SLPS.DAT"

configextra5

tit

PlaneStrainSimulation

grid1517

modmo

monulli15j817

gen0-150012012-15i16j18

genss80080-15rat1.11i616j18

gens24.5258025srat1.11i616j818

fixx

i1

fixxyj1

fixxi16

probul2e8she1e8den2500fric45coh1e6ten1e6

setgrav10

waterde

n1000

table10012018.2512.525173521502474.5258025

watertable1

definidens

loopi（1,izones）

loopj（1,jzones）

ifpp（i,j）#0then

density（i,j）=2600.

endif

endloop

endloop

end

inidens

inisyy-1e6var01e6

inisxx-1e6var01e6

iniszz-1e6var01e6

hisunbal

histxdi6j18

histydi6j18

solve

saveslps1.sav

inixd0yd0

inixv0yv0

procoh3.7e4

step1250

saveslps2.sav

setoutslps1.plt

plotpenh23

setoutslps2.plt

plotpenbouydispzfillabs

tit

flacmodelgrid

setoutslgria.plt

plotpengridwateriwhite

ret

;DataFile"SLAXI.DAT"

configextra5axi

tit

AxisymmetrySimulation

grid1517

modmo

monulli15j817

gen0-150012012-15i16j18

genss80080-15rat1.11i616j18

gens24.5258025srat1.11i616j818

fixxi1

fixxyj1

fixxi16

probul2e8she1e8den2500fric45coh1e6ten1e6

setgrav10

waterden1000

table10012018.2512.525173521502474.5258025

watertable1

definidens

loopi（1,izones）

loopj（1,jzones）

ifpp（i,j）#0then

density（i,j）=2600.

endif

endloop

endloop

end

inidens

inisyy-1e6var01e6

inisxx-1e6var01e6

iniszz-1e6var01e6

hisunbal

histxdi6j18

histydi6j18

solve

saveslax1.sav

inixd0yd0

inixv0yv0

procoh3.7e4

step1250

saveslax2.sav

setoutslax1.plt

plotpenh23min-1.7

ret

;DataFile"WATER.FIS"

defparm

nptab1=7

nprof=10

end

parm

table112018.2512.5251735215024752515025

defaux1

loopn（1,nptab1）

rr=xtable（1,n）

zz=ytable（1,n）

loopk（1,nprof）

caseofk

alfa=0.5*（pi+.1）*float（k-2）/float（nprof-2）

xtable（n+10,k）=rr*cos（alfa）

ytable（n+10,k）=rr*sin（alfa）

case1

xtable（n+10,k）=rr

ytable（n+10,k）=-40.0

case2

xtable（n+10,k）=rr

ytable（n+10,k）=0.0

endcase

endloop

endloop

end

aux1

defaux2

command

setgrav00-10

waterdens1000

endcommand

loopn（1,nptab1-1）

loopk（1,nprof-1）

xx1=xtable（n+10,k）

yy1=ytable（n+10,k）

zz1=ytable（1,n）

xx2=xtable（n+10,k+1）

yy2=ytable（n+10,k+1）

zz2=zz1

xx3=xtable（n+11,k+1）

yy3=ytable（n+11,k+1）

zz3=ytable（1,n+1）

xx4=xtable（n+11,k）

yy4=ytable（n+11,k）

zz4=zz3

command

watertablefacexx1,yy1,zz1xx2,yy2,zz2xx3,yy3,zz3

watertablefacexx1,yy1,zz1xx3,yy3,zz3xx4,yy4,zz4

endcommand

endloop

endloop

end

aux2

ret