ABAQUS-二次开发-边界条件

这个例子模拟三峡库区的水位涨落时，涉水土质滑坡的渗流场分布情况的，是以前做的，蛮好玩的。如果大家关心些新闻的话，三峡库区有涨水、蓄水、排水、枯水这个一年一次循环的周期，关键问题就在于怎么在数值模拟中考虑这个

时间单位这里用的是“天”，任何问题，只要把单位统一化，就可以实现自己所定义的单位系统下的问题，这个大家应该是很熟悉的了。

t=2天

t=4天

t=8天

t=16天

t=20天t=26天

图2 浸润线位置随涨水时间t变化图

数值模拟中实现这个问题，需要对边界条件上加载孔隙水压力的方式动手脚了，首先把时间定下来，然后把每个时间对应的水位高度定下来，然后就是写对应的程序了，关键就在于，需要在各个位置的节点处定义不同的pp_t幅值曲线。这个问题使用子程序会很方便，也可以不用，只需要定义一堆关键字吧，但是GUI方式是完全没法实现的。其实软件关键字的背后也就是他的脚本语言，就好像是FLAC里写一句struct cable，这个命令本身只有给你看懂那么一个傻瓜意义意义。看看关键字怎么定义的：

*AMPLITUDE,name=down1,VALUE=ABSOLUTE

0,0,30,-300,180,-300

*AMPLITUDE,name=down2,VALUE=ABSOLUTE

0,10,1,0,30,-290,180,-290

*AMPLITUDE,name=down3,VALUE=ABSOLUTE

0,20,2,0,30,-280,180,-280

*AMPLITUDE,name=down4,VALUE=ABSOLUTE

0,30,3,0,30,-270,180,-270

*AMPLITUDE,name=down5,VALUE=ABSOLUTE

0,40,4,0,30,-260,180,-260

*AMPLITUDE,name=down6,VALUE=ABSOLUTE

0,50,5,0,30,-250,180,-250

*AMPLITUDE,name=down7,VALUE=ABSOLUTE

0,60,6,0,30,-240,180,-240

*AMPLITUDE,name=down8,VALUE=ABSOLUTE

0,70,7,0,30,-230,180,-230

*AMPLITUDE,name=down9,VALUE=ABSOLUTE

0,80,8,0,30,-220,180,-220

*AMPLITUDE,name=down10,VALUE=ABSOLUTE

0,90,9,0,30,-210,180,-210

*AMPLITUDE,name=down11,VALUE=ABSOLUTE

0,100,10,0,30,-200,180,-200

*AMPLITUDE,name=down12,VALUE=ABSOLUTE

0,110,11,0,30,-190,180,-190

*AMPLITUDE,name=down13,VALUE=ABSOLUTE

0,120,12,0,30,-180,180,-180

*AMPLITUDE,name=down14,VALUE=ABSOLUTE

0,130,13,0,30,-170,180,-170

*AMPLITUDE,name=down15,VALUE=ABSOLUTE

0,140,14,0,30,-160,180,-160

*AMPLITUDE,name=down16,VALUE=ABSOLUTE

0,150,15,0,30,-150,180,-150

*AMPLITUDE,name=down17,VALUE=ABSOLUTE 0,160,16,0,30,-140,180,-140

*AMPLITUDE,name=down18,VALUE=ABSOLUTE 0,170,17,0,30,-130,180,-130

*AMPLITUDE,name=down19,VALUE=ABSOLUTE 0,180,18,0,30,-120,180,-120

*AMPLITUDE,name=down20,VALUE=ABSOLUTE 0,190,19,0,30,-110,180,-110

*AMPLITUDE,name=down21,VALUE=ABSOLUTE 0,200,20,0,30,-100,180,-100

*AMPLITUDE,name=down22,VALUE=ABSOLUTE 0,210,21,0,30,-90,180,-90

*AMPLITUDE,name=down23,VALUE=ABSOLUTE 0,220,22,0,30,-80,180,-80

*AMPLITUDE,name=down24,VALUE=ABSOLUTE 0,230,23,0,30,-70,180,-70

*AMPLITUDE,name=down25,VALUE=ABSOLUTE 0,240,24,0,30,-60,180,-60

*AMPLITUDE,name=down26,VALUE=ABSOLUTE 0,250,25,0,30,-50,180,-50

*AMPLITUDE,name=down27,VALUE=ABSOLUTE 0,260,26,0,30,-40,180,-40

*AMPLITUDE,name=down28,VALUE=ABSOLUTE 0,270,27,0,30,-30,180,-30

*AMPLITUDE,name=down29,VALUE=ABSOLUTE 0,280,28,0,30,-20,180,-20

*AMPLITUDE,name=down30,VALUE=ABSOLUTE 0,290,29,0,30,-10,180,-10

*AMPLITUDE,name=down31,VALUE=ABSOLUTE 0,300,29.9,0,180,0

*BOUNDARY,AMPLITUDE=down1

Part-1-1.nod30,8,,0

*BOUNDARY,AMPLITUDE=down2

Part-1-1.nod29,8,,10

*BOUNDARY,AMPLITUDE=down3

Part-1-1.nod28,8,,20

*BOUNDARY,AMPLITUDE=down4