ABAQUS-二次开发-边界条件

图 1 水位涨落周期变化曲线

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

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

Part-1-1.nod27,8,,30

*BOUNDARY,AMPLITUDE=down5

Part-1-1.nod26,8,,40

*BOUNDARY,AMPLITUDE=down6

Part-1-1.nod25,8,,50

*BOUNDARY,AMPLITUDE=down7

Part-1-1.nod24,8,,60

*BOUNDARY,AMPLITUDE=down8

Part-1-1.nod23,8,,70

*BOUNDARY,AMPLITUDE=down9

Part-1-1.nod22,8,,80

*BOUNDARY,AMPLITUDE=down10

Part-1-1.nod21,8,,90

*BOUNDARY,AMPLITUDE=down11

Part-1-1.nod20,8,,100

*BOUNDARY,AMPLITUDE=down12

Part-1-1.nod19,8,,110

*BOUNDARY,AMPLITUDE=down13

Part-1-1.nod18,8,,120

*BOUNDARY,AMPLITUDE=down14

Part-1-1.nod17,8,,130

*BOUNDARY,AMPLITUDE=down15

Part-1-1.nod16,8,,140

*BOUNDARY,AMPLITUDE=down16

Part-1-1.nod15,8,,150

*BOUNDARY,AMPLITUDE=down17

Part-1-1.nod14,8,,160

*BOUNDARY,AMPLITUDE=down18

Part-1-1.nod13,8,,170

*BOUNDARY,AMPLITUDE=down19

Part-1-1.nod12,8,,180

*BOUNDARY,AMPLITUDE=down20

Part-1-1.nod11,8,,190

*BOUNDARY,AMPLITUDE=down21

Part-1-1.nod10,8,,200

*BOUNDARY,AMPLITUDE=down22

Part-1-1.nod9,8,,210

*BOUNDARY,AMPLITUDE=down23

Part-1-1.nod8,8,,220

*BOUNDARY,AMPLITUDE=down24

Part-1-1.nod7,8,,230

*BOUNDARY,AMPLITUDE=down25

Part-1-1.nod6,8,,240

*BOUNDARY,AMPLITUDE=down26

Part-1-1.nod5,8,,250

*BOUNDARY,AMPLITUDE=down27

Part-1-1.nod4,8,,260

*BOUNDARY,AMPLITUDE=down28

Part-1-1.nod3,8,,270

*BOUNDARY,AMPLITUDE=down29

Part-1-1.nod2,8,,280

*BOUNDARY,AMPLITUDE=down30

Part-1-1.nod1,8,,290

*BOUNDARY,AMPLITUDE=down31

Part-1-1.nod0,8,,300

以上是一部分关键字，其实根本问题就是找到对应节点，对各个节点附上与时间相关的边界。这个问题，其实只能算是半成品式的二次开发吧。用过LS-DYNA可能更清楚我写的这个东西是怎么一回事。