Abquse模拟库水位升降的例子(含完整程序)

Abaqus模拟库水位升降的例子—含附件Abaqus可以通过子程序或者编辑关键字来实现库水位升降的模拟。如下例子，一均质土石坝受水位变化时浸润线分布，设该土石坝高12m，水位变化如图：数值模拟中实现这个问题，需要对边界条件上加载孔隙水压力，首先要把时间定下来，然后把每个时间对应的水位高度定下来，然后就是写对应的程序了，关键就在于，需要在各个位置的节点处定义不同的幅值曲线。这个问题使用子程序会很方便，也可以不用，只需要定义一堆关键字吧，但是GUI方式是完全没法实现的。下面为实现该模拟的关键字（完整文件见附件）：首先定义幅值曲线，以下为例说明幅值曲线的含义：*Amplitude为关键字，即幅值曲线name=Amp-1幅值曲线的名称定义为Amp-1time=TOTALTIME时间有两种，一种是steptime，即以该分析步时间为基准，并只在该分析步时间生效；另一种是totaltime，即作用于所有分析步中，以第一个（initial之后的第一个分析步）分析步为初始时间。0.,0.,180.,110.,1980.,110.,2160.,0.上一行为曲线的定义，即表示了下图的曲线。02468101205001000150020002500水位时间库水位随时间的变化分别定义上游坝面与水接触的各个节点的幅值曲线：*Amplitude,name=Amp-1,time=TOTALTIME0.,0.,180.,110.,1980.,110.,2160.,0.*Amplitude,name=Amp-2,time=TOTALTIME0.,-10.,180.,100.,1980.,100.,2160.,-10.*Amplitude,name=Amp-3,time=TOTALTIME0.,-20.,180.,90.,1980.,90.,2160.,-20.*Amplitude,name=Amp-4,time=TOTALTIME0.,-30.,180.,80.,1980.,80.,2160.,-30.*Amplitude,name=Amp-5,time=TOTALTIME0.,-40.,180.,70.,1980.,70.,2160.,-40.*Amplitude,name=Amp-6,time=TOTALTIME0.,-50.,180.,60.,1980.,60.,2160.,-50.*Amplitude,name=Amp-7,time=TOTALTIME0.,-60.,180.,50.,1980.,50.,2160.,-60.*Amplitude,name=Amp-8,time=TOTALTIME0.,-70.,180.,40.,1980.,40.,2160.,-70.*Amplitude,name=Amp-9,time=TOTALTIME0.,-80.,180.,30.,1980.,30.,2160.,-80.*Amplitude,name=Amp-10,time=TOTALTIME0.,-90.,180.,20.,1980.,20.,2160.,-90.*Amplitude,name=Amp-11,time=TOTALTIME02468101205001000150020002500孔压时间孔压随时间的变化0.,-100.,180.,10.,1980.,10.,2160.,-100.*Amplitude,name=Amp-12,time=TOTALTIME0.,-110.,180.,0.,1980.,0.,2160.,-110.*Amplitude,name=Amp-13,time=TOTALTIME0.,-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-1Part-1-1.13,8,8,1.*Boundary,amplitude=Amp-2Part-1-1.12,8,8,1.*Boundary,amplitude=Amp-3Part-1-1.11,8,8,1.*Boundary,amplitude=Amp-4Part-1-1.10,8,8,1.*Boundary,amplitude=Amp-5Part-1-1.9,8,8,1.*Boundary,amplitude=Amp-6Part-1-1.8,8,8,1.*Boundary,amplitude=Amp-7Part-1-1.7,8,8,1.*Boundary,amplitude=Amp-8Part-1-1.6,8,8,1.*Boundary,amplitude=Amp-9Part-1-1.5,8,8,1.*Boundary,amplitude=Amp-10Part-1-1.4,8,8,1.*Boundary,amplitude=Amp-11Part-1-1.3,8,8,1.*Boundary,amplitude=Amp-12Part-1-1.2,8,8,1.*Boundary,amplitude=Amp-13Part-1-1.1,8,8,1.Abquse计算浸润线结果如下，由图可见，库水位在第180天上升至11m时坝体渗流并未稳定，在第240天坝体渗流稳定，库水位下降时，形成向上游防渗，对坝坡不利：T=1T=40.7T=130T=180T=210T=240T=1980T=2000T=2110T=2160附件：****PARTS***Part,name=Part-1*Element,type=CPE4P*Surface,type=ELEMENT,name=Surf-up*Surface,type=ELEMENT,name=Surf-fot*Surface,type=ELEMENT,name=Surf-down**Section:Section-1*SolidSection,elset=_PickedSet2,material=soil,*EndPart******ASSEMBLY***Assembly,name=Assembly***Instance,name=Part-1-1,part=Part-1*EndInstance***EndAssembly*Amplitude,name=Amp-1,time=TOTALTIME0.,0.,180.,110.,1980.,110.,2160.,0.*Amplitude,name=Amp-2,time=TOTALTIME0.,-10.,180.,100.,1980.,100.,2160.,-10.*Amplitude,name=Amp-3,time=TOTALTIME0.,-20.,180.,90.,1980.,90.,2160.,-20.*Amplitude,name=Amp-4,time=TOTALTIME0.,-30.,180.,80.,1980.,80.,2160.,-30.*Amplitude,name=Amp-5,time=TOTALTIME0.,-40.,180.,70.,1980.,70.,2160.,-40.*Amplitude,name=Amp-6,time=TOTALTIME0.,-50.,180.,60.,1980.,60.,2160.,-50.*Amplitude,name=Amp-7,time=TOTALTIME0.,-60.,180.,50.,1980.,50.,2160.,-60.*Amplitude,name=Amp-8,time=TOTALTIME0.,-70.,180.,40.,1980.,40.,2160.,-70.*Amplitude,name=Amp-9,time=TOTALTIME0.,-80.,180.,30.,1980.,30.,2160.,-80.*Amplitude,name=Amp-10,time=TOTALTIME0.,-90.,180.,20.,1980.,20.,2160.,-90.*Amplitude,name=Amp-11,time=TOTALTIME0.,-100.,180.,10.,1980.,10.,2160.,-100.*Amplitude,name=Amp-12,time=TOTALTIME0.,-110.,180.,0.,1980.,0.,2160.,-110.*Amplitude,name=Amp-13,time=TOTALTIME0.,-120.,180.,-10.,1980.,-100.,2160.,-120.****MATERIALS***Material,name=soil*Density2.,*Elastic10000.,0.25*Permeability,specific=10.0.06912,0.*Permeability,Type=SATURATION0.000116,0.142970.000232,0.178850.000561,0.200180.001442,0.225730.003736,0.259480.009699,0.305810.02478,0.371390.06355,0.463580.156538,0.587650.335775,0.734890.575975,0.872120.75905,0.95931.,1.*Sorption-100.,0.167-72.5,0.207-54.622,0.317-43.887,0.466-34.019,0.68525-24.5,0.895-13.639,0.9940.,1.****BOUNDARYCONDITIONS****Name:BC-1Type:Displacement/Rotation**Boundary中涉及PickedSet的要注意，本例约束底面、侧面和顶面所有位移，*Boundary_PickedSet4,1,1_PickedSet4,2,2**Name:BC-2Type:Displacement/Rotation*Boundary_PickedSet5,1,1_PickedSet5,2,2**Name:BC-3Type:Displacement/Rotation*Boundary_PickedSet6,2,2****PREDEFINEDFIELDS****Name:PredefinedField-1Type:Saturation*InitialConditions,type=SATURATION_PickedSet7,0.167**Name:PredefinedField-2Type:Voidratio*InitialConditions,TYPE=RATIO_PickedSet8,1.**Name:PredefinedField-3Type:PorepressureUsingField:inipre*InitialConditions,TYPE=POREPRESSUREPart-1-1.1,-120.Part-1-1.2,-110.Part-1-1.3,-100.Part-1-1.4,-90.Part-1-1.5,-80.Part-1-1.6,-70.Part-1-1.7,-60.Part-1-1.8,-50.Part-1-1.9,-40.Part-1-1.10,-30.Part-1-1.11,-20.Part-1-1.12,-10.Part-1-1.13,0.Part-1-1.14,-120.Part-1-1.15,-110.Part-1-1.16,-100.Part-1-1.17,-90.Part-1-1.18,-80.Part-1-1.19,-7