Abquse模拟库水位升降的例子(含完整程序)
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Abaqus 模拟库水位升降的例子—含附件
Abaqus可以通过子程序或者编辑关键字来实现库水位升降的模拟。如下例子，一均质土石坝受水位变化时浸润线分布，设该土石坝高12m，水位变化如图：
数值模拟中实现这个问题，需要对边界条件上加载孔隙水压力，首先要把时间定下来，然后把每个时间对应的水位高度定下来，然后就是写对应的程序了，关键就在于，需要在各个位置的节点处定义不同的幅值曲线。这个问题使用子程序会很方便，也可以不用，只需要定义一堆关键字吧，但是GUI方式是完全没法实现的。下面为实现该模拟的关键字（完整文件见附件）：
首先定义幅值曲线，以下为例说明幅值曲线的含义：
*Amplitude 为关键字，即幅值曲线
name=Amp-1 幅值曲线的名称定义为Amp-1
time=TOTAL TIME 时间有两种，一种是step time，即以该分析步时间为基准，并只在该分析步时间生效；另一种是total time，即作用于所有分析步中，以第一个（initial之后的第一个分析步）分析步为初始时间。
0., 0., 180., 110., 1980., 110., 2160., 0.
上一行为曲线的定义，即表示了下图的曲线。
分别定义上游坝面与水接触的各个节点的幅值曲线：*Amplitude, name=Amp-1, time=TOTAL TIME
0., 0., 180., 110., 1980., 110., 2160., 0.
*Amplitude, name=Amp-2, time=TOTAL TIME
0., -10., 180., 100., 1980., 100., 2160., -10. *Amplitude, name=Amp-3, time=TOTAL TIME
0., -20., 180., 90., 1980., 90., 2160., -20.
*Amplitude, name=Amp-4, time=TOTAL TIME
0., -30., 180., 80., 1980., 80., 2160., -30.
*Amplitude, name=Amp-5, time=TOTAL TIME
0., -40., 180., 70., 1980., 70., 2160., -40.
*Amplitude, name=Amp-6, time=TOTAL TIME
0., -50., 180., 60., 1980., 60., 2160., -50.
*Amplitude, name=Amp-7, time=TOTAL TIME
0., -60., 180., 50., 1980., 50., 2160., -60.
*Amplitude, name=Amp-8, time=TOTAL TIME
0., -70., 180., 40., 1980., 40., 2160., -70.
*Amplitude, name=Amp-9, time=TOTAL TIME
0., -80., 180., 30., 1980., 30., 2160., -80.
*Amplitude, name=Amp-10, time=TOTAL TIME
0., -90., 180., 20., 1980., 20., 2160., -90.
*Amplitude, name=Amp-11, time=TOTAL TIME
0., -100., 180., 10., 1980., 10., 2160., -100.
*Amplitude, name=Amp-12, time=TOTAL TIME
0., -110., 180., 0., 1980., 0., 2160., -110.
*Amplitude, name=Amp-13, time=TOTAL TIME
0., -120., 180., -10., 1980., -100., 2160., -120.
接下来就是定义各个节点的边界条件——孔压：
*Boundary 为关键字，表示边界条件
amplitude=Amp-1 采用名为Amp-1的幅值曲线定义边界条件
Part-1-1.13, 8, 8, 1. Part-1-1.13为Part-1-1部件的编号为13的节点；8, 8表示孔压自由度；1表示大小为1* Amp-1。
分别为上游坝面与水接触的各个节点的添加孔压边界条件：
*Boundary, amplitude=Amp-1
Part-1-1.13, 8, 8, 1.
*Boundary, amplitude=Amp-2
Part-1-1.12, 8, 8,1.
*Boundary, amplitude=Amp-3
Part-1-1.11, 8, 8,1.
*Boundary, amplitude=Amp-4
Part-1-1.10, 8, 8,1.
*Boundary, amplitude=Amp-5
Part-1-1.9, 8, 8,1.
*Boundary, amplitude=Amp-6
Part-1-1.8, 8, 8,1.
*Boundary, amplitude=Amp-7
Part-1-1.7, 8, 8,1.
*Boundary, amplitude=Amp-8
Part-1-1.6, 8, 8,1.
*Boundary, amplitude=Amp-9
Part-1-1.5, 8, 8,1.