图 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可能更清楚我写的这个东西是怎么一回事。下载本文
