《Abquse模拟库水位升降的例子(共26页).docx》由会员分享，可在线阅读，更多相关《Abquse模拟库水位升降的例子(共26页).docx（26页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上Abaqus 模拟库水位升降的例子含附件Abaqus可以通过子程序或者编辑关键字来实现库水位升降的模拟。如下例子，一均质土石坝受水位变化时浸润线分布，设该土石坝高12m，水位变化如图：数值模拟中实现这个问题，需要对边界条件上加载孔隙水压力，首先要把时间定下来，然后把每个时间对应的水位高度定下来，然后就是写对应的程序了，关键就在于，需要在各个位置的节点处定义不同的幅值曲线。这个问题使用子程序会很方便，也可以不用，只需要定义一堆关键字吧，但是GUI方式是完全没法实现的。下面为实现该模拟的关键字（完整文件见附件）：首先定义幅值曲线，以下为例说明幅值曲线的含义：*Amplit
2、ude 为关键字，即幅值曲线name=Amp-1 幅值曲线的名称定义为Amp-1time=TOTAL TIME 时间有两种，一种是step time，即以该分析步时间为基准，并只在该分析步时间生效；另一种是total time，即作用于所有分析步中，以第一个（initial之后的第一个分析步）分析步为初始时间。0., 0., 180., 110., 1980., 110., 2160., 0.上一行为曲线的定义，即表示了下图的曲线。分别定义上游坝面与水接触的各个节点的幅值曲线：*Amplitude, name=Amp-1, time=TOTAL TIME0., 0., 180., 110.,
3、1980., 110., 2160., 0.*Amplitude, name=Amp-2, time=TOTAL TIME0., -10., 180., 100., 1980., 100., 2160., -10.*Amplitude, name=Amp-3, time=TOTAL TIME0., -20., 180., 90., 1980., 90., 2160., -20.*Amplitude, name=Amp-4, time=TOTAL TIME0., -30., 180., 80., 1980., 80., 2160., -30.*Amplitude, name=Amp-5, tim
4、e=TOTAL TIME0., -40., 180., 70., 1980., 70., 2160., -40.*Amplitude, name=Amp-6, time=TOTAL TIME0., -50., 180., 60., 1980., 60., 2160., -50.*Amplitude, name=Amp-7, time=TOTAL TIME0., -60., 180., 50., 1980., 50., 2160., -60.*Amplitude, name=Amp-8, time=TOTAL TIME0., -70., 180., 40., 1980., 40., 2160.,
5、 -70.*Amplitude, name=Amp-9, time=TOTAL TIME0., -80., 180., 30., 1980., 30., 2160., -80.*Amplitude, name=Amp-10, time=TOTAL TIME0., -90., 180., 20., 1980., 20., 2160., -90.*Amplitude, name=Amp-11, time=TOTAL TIME0., -100., 180., 10., 1980., 10., 2160., -100.*Amplitude, name=Amp-12, time=TOTAL TIME0.
6、, -110., 180., 0., 1980., 0., 2160., -110.*Amplitude, name=Amp-13, time=TOTAL TIME0., -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。分别为上游坝面与水接触的
7、各个节点的添加孔压边界条件：*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,
8、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,
9、 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-dow
10、n* Section: Section-1*Solid Section, elset=_PickedSet2, material=soil,*End Part* * ASSEMBLY*Assembly, name=Assembly* *Instance, name=Part-1-1, part=Part-1*End Instance* *End Assembly*Amplitude, name=Amp-1, time=TOTAL TIME0., 0., 180., 110., 1980., 110., 2160., 0.*Amplitude, name=Amp-2, time=TOTAL TI
11、ME0., -10., 180., 100., 1980., 100., 2160., -10.*Amplitude, name=Amp-3, time=TOTAL TIME0., -20., 180., 90., 1980., 90., 2160., -20.*Amplitude, name=Amp-4, time=TOTAL TIME0., -30., 180., 80., 1980., 80., 2160., -30.*Amplitude, name=Amp-5, time=TOTAL TIME0., -40., 180., 70., 1980., 70., 2160., -40.*Am
12、plitude, name=Amp-6, time=TOTAL TIME0., -50., 180., 60., 1980., 60., 2160., -50.*Amplitude, name=Amp-7, time=TOTAL TIME0., -60., 180., 50., 1980., 50., 2160., -60.*Amplitude, name=Amp-8, time=TOTAL TIME0., -70., 180., 40., 1980., 40., 2160., -70.*Amplitude, name=Amp-9, time=TOTAL TIME0., -80., 180.,
13、 30., 1980., 30., 2160., -80.*Amplitude, name=Amp-10, time=TOTAL TIME0., -90., 180., 20., 1980., 20., 2160., -90.*Amplitude, name=Amp-11, time=TOTAL TIME0., -100., 180., 10., 1980., 10., 2160., -100.*Amplitude, name=Amp-12, time=TOTAL TIME0., -110., 180., 0., 1980., 0., 2160., -110.*Amplitude, name=
14、Amp-13, time=TOTAL TIME0., -120., 180., -10., 1980., -100., 2160., -120.* * MATERIALS*Material, name=soil*Density2.,*Elastic10000., 0.25*Permeability, specific=10. 0.06912,0.*Permeability, Type=SATURATION 0., 0.14297 0., 0.17885 0., 0.20018 0., 0.22573 0., 0.25948 0., 0.30581 0.02478, 0.37139 0.0635
15、5, 0.46358 0., 0.58765 0., 0.73489 0., 0.87212 0.75905, 0.9593 1., 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.994 0., 1.* * BOUNDARY CONDITIONS* * Name: BC-1 Type: Displacement/Rotation* Boundary中涉及PickedSet的要注意，本例约束底面、侧面和顶面所有位移，*Boun
16、dary_PickedSet4, 1, 1_PickedSet4, 2, 2* Name: BC-2 Type: Displacement/Rotation*Boundary_PickedSet5, 1, 1_PickedSet5, 2, 2* Name: BC-3 Type: Displacement/Rotation*Boundary_PickedSet6, 2, 2* * PREDEFINED FIELDS* * Name: Predefined Field-1 Type: Saturation*Initial Conditions, type=SATURATION_PickedSet7
17、, 0.167* Name: Predefined Field-2 Type: Void ratio*Initial Conditions, TYPE=RATIO_PickedSet8, 1.* Name: Predefined Field-3 Type: Pore pressure Using Field: inipre*Initial Conditions, TYPE=PORE PRESSUREPart-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.
18、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, -70.Part-1-1.20, -60.Part-1-1.21, -50.Part-1-1.22, -40.Part-1-1.23, -30.Part-1-1.24,
19、-20.Part-1-1.25, -10.Part-1-1.26, 0.Part-1-1.27, -120.Part-1-1.28, -110.Part-1-1.29, -100.Part-1-1.30, -90.Part-1-1.31, -80.Part-1-1.32, -70.Part-1-1.33, -60.Part-1-1.34, -50.Part-1-1.35, -40.Part-1-1.36, -30.Part-1-1.37, -20.Part-1-1.38, -10.Part-1-1.39, 0.Part-1-1.40, -120.Part-1-1.41, -110.Part-1
20、-1.42, -100.Part-1-1.43, -90.Part-1-1.44, -80.Part-1-1.45, -70.Part-1-1.46, -60.Part-1-1.47, -50.Part-1-1.48, -40.Part-1-1.49, -30.Part-1-1.50, -20.Part-1-1.51, -10.Part-1-1.52, 0.Part-1-1.53, -120.Part-1-1.54, -110.Part-1-1.55, -100.Part-1-1.56, -90.Part-1-1.57, -80.Part-1-1.58, -70.Part-1-1.59, -6
21、0.Part-1-1.60, -50.Part-1-1.61, -40.Part-1-1.62, -30.Part-1-1.63, -20.Part-1-1.64, -10.Part-1-1.65, 0.Part-1-1.66, -120.Part-1-1.67, -110.Part-1-1.68, -100.Part-1-1.69, -90.Part-1-1.70, -80.Part-1-1.71, -70.Part-1-1.72, -60.Part-1-1.73, -50.Part-1-1.74, -40.Part-1-1.75, -30.Part-1-1.76, -20.Part-1-1
22、.77, -10.Part-1-1.78, 0.Part-1-1.79, -120.Part-1-1.80, -110.Part-1-1.81, -100.Part-1-1.82, -90.Part-1-1.83, -80.Part-1-1.84, -70.Part-1-1.85, -60.Part-1-1.86, -50.Part-1-1.87, -40.Part-1-1.88, -30.Part-1-1.89, -20.Part-1-1.90, -10.Part-1-1.91, 0.Part-1-1.92, -120.Part-1-1.93, -110.Part-1-1.94, -100.
23、Part-1-1.95, -90.Part-1-1.96, -80.Part-1-1.97, -70.Part-1-1.98, -60.Part-1-1.99, -50.Part-1-1.100, -40.Part-1-1.101, -30.Part-1-1.102, -20.Part-1-1.103, -10.Part-1-1.104, 0.Part-1-1.105, -120.Part-1-1.106, -110.Part-1-1.107, -100.Part-1-1.108, -90.Part-1-1.109, -80.Part-1-1.110, -70.Part-1-1.111, -6
24、0.Part-1-1.112, -50.Part-1-1.113, -40.Part-1-1.114, -30.Part-1-1.115, -20.Part-1-1.116, -10.Part-1-1.117, 0.Part-1-1.118, -120.Part-1-1.119, -110.Part-1-1.120, -100.Part-1-1.121, -90.Part-1-1.122, -80.Part-1-1.123, -70.Part-1-1.124, -60.Part-1-1.125, -50.Part-1-1.126, -40.Part-1-1.127, -30.Part-1-1.
25、128, -20.Part-1-1.129, -10.Part-1-1.130, 0.Part-1-1.131, -120.Part-1-1.132, -110.Part-1-1.133, -100.Part-1-1.134, -90.Part-1-1.135, -80.Part-1-1.136, -70.Part-1-1.137, -60.Part-1-1.138, -50.Part-1-1.139, -40.Part-1-1.140, -30.Part-1-1.141, -20.Part-1-1.142, -10.Part-1-1.143, 0.Part-1-1.144, -120.Par
26、t-1-1.145, -110.Part-1-1.146, -100.Part-1-1.147, -90.Part-1-1.148, -80.Part-1-1.149, -70.Part-1-1.150, -60.Part-1-1.151, -50.Part-1-1.152, -40.Part-1-1.153, -30.Part-1-1.154, -20.Part-1-1.155, -10.Part-1-1.156, 0.Part-1-1.157, -120.Part-1-1.158, -110.Part-1-1.159, -100.Part-1-1.160, -90.Part-1-1.161
27、, -80.Part-1-1.162, -70.Part-1-1.163, -60.Part-1-1.164, -50.Part-1-1.165, -40.Part-1-1.166, -30.Part-1-1.167, -20.Part-1-1.168, -10.Part-1-1.169, 0.Part-1-1.170, -120.Part-1-1.171, -110.Part-1-1.172, -100.Part-1-1.173, -90.Part-1-1.174, -80.Part-1-1.175, -70.Part-1-1.176, -60.Part-1-1.177, -50.Part-
28、1-1.178, -40.Part-1-1.179, -30.Part-1-1.180, -20.Part-1-1.181, -10.Part-1-1.182, 0.Part-1-1.183, -120.Part-1-1.184, -110.Part-1-1.185, -100.Part-1-1.186, -90.Part-1-1.187, -80.Part-1-1.188, -70.Part-1-1.189, -60.Part-1-1.190, -50.Part-1-1.191, -40.Part-1-1.192, -30.Part-1-1.193, -20.Part-1-1.194, -1
29、0.Part-1-1.195, 0.Part-1-1.196, -120.Part-1-1.197, -110.Part-1-1.198, -100.Part-1-1.199, -90.Part-1-1.200, -80.Part-1-1.201, -70.Part-1-1.202, -60.Part-1-1.203, -50.Part-1-1.204, -40.Part-1-1.205, -30.Part-1-1.206, -20.Part-1-1.207, -10.Part-1-1.208, 0.Part-1-1.209, -120.Part-1-1.210, -110.Part-1-1.
30、211, -100.Part-1-1.212, -90.Part-1-1.213, -80.Part-1-1.214, -70.Part-1-1.215, -60.Part-1-1.216, -50.Part-1-1.217, -40.Part-1-1.218, -30.Part-1-1.219, -20.Part-1-1.220, -10.Part-1-1.221, 0.Part-1-1.222, -120.Part-1-1.223, -110.Part-1-1.224, -100.Part-1-1.225, -90.Part-1-1.226, -80.Part-1-1.227, -70.P
31、art-1-1.228, -60.Part-1-1.229, -50.Part-1-1.230, -40.Part-1-1.231, -30.Part-1-1.232, -20.Part-1-1.233, -10.Part-1-1.234, 0.Part-1-1.235, -120.Part-1-1.236, -110.Part-1-1.237, -100.Part-1-1.238, -90.Part-1-1.239, -80.Part-1-1.240, -70.Part-1-1.241, -60.Part-1-1.242, -50.Part-1-1.243, -40.Part-1-1.244
32、, -30.Part-1-1.245, -20.Part-1-1.246, -10.Part-1-1.247, 0.Part-1-1.248, -120.Part-1-1.249, -110.Part-1-1.250, -100.Part-1-1.251, -90.Part-1-1.252, -80.Part-1-1.253, -70.Part-1-1.254, -60.Part-1-1.255, -50.Part-1-1.256, -40.Part-1-1.257, -30.Part-1-1.258, -20.Part-1-1.259, -10.Part-1-1.260, 0.Part-1-
33、1.261, -120.Part-1-1.262, -110.Part-1-1.263, -100.Part-1-1.264, -90.Part-1-1.265, -80.Part-1-1.266, -70.Part-1-1.267, -60.Part-1-1.268, -50.Part-1-1.269, -40.Part-1-1.270, -30.Part-1-1.271, -20.Part-1-1.272, -10.Part-1-1.273, 0.Part-1-1.274, -120.Part-1-1.275, -110.Part-1-1.276, -100.Part-1-1.277, -
34、90.Part-1-1.278, -80.Part-1-1.279, -70.Part-1-1.280, -60.Part-1-1.281, -50.Part-1-1.282, -40.Part-1-1.283, -30.Part-1-1.284, -20.Part-1-1.285, -10.Part-1-1.286, 0.Part-1-1.287, -120.Part-1-1.288, -110.Part-1-1.289, -100.Part-1-1.290, -90.Part-1-1.291, -80.Part-1-1.292, -70.Part-1-1.293, -60.Part-1-1
35、.294, -50.Part-1-1.295, -40.Part-1-1.296, -30.Part-1-1.297, -20.Part-1-1.298, -10.Part-1-1.299, 0.Part-1-1.300, -120.Part-1-1.301, -110.Part-1-1.302, -100.Part-1-1.303, -90.Part-1-1.304, -80.Part-1-1.305, -70.Part-1-1.306, -60.Part-1-1.307, -50.Part-1-1.308, -40.Part-1-1.309, -30.Part-1-1.310, -20.P
36、art-1-1.311, -10.Part-1-1.312, 0.Part-1-1.313, -120.Part-1-1.314, -110.Part-1-1.315, -100.Part-1-1.316, -90.Part-1-1.317, -80.Part-1-1.318, -70.Part-1-1.319, -60.Part-1-1.320, -50.Part-1-1.321, -40.Part-1-1.322, -30.Part-1-1.323, -20.Part-1-1.324, -10.Part-1-1.325, 0.Part-1-1.326, -120.Part-1-1.327,
37、 -110.Part-1-1.328, -100.Part-1-1.329, -90.Part-1-1.330, -80.Part-1-1.331, -70.Part-1-1.332, -60.Part-1-1.333, -50.Part-1-1.334, -40.Part-1-1.335, -30.Part-1-1.336, -20.Part-1-1.337, -10.Part-1-1.338, 0.Part-1-1.339, -120.Part-1-1.340, -110.Part-1-1.341, -100.Part-1-1.342, -90.Part-1-1.343, -80.Part
38、-1-1.344, -70.Part-1-1.345, -60.Part-1-1.346, -50.Part-1-1.347, -40.Part-1-1.348, -30.Part-1-1.349, -20.Part-1-1.350, -10.Part-1-1.351, 0.Part-1-1.352, -120.Part-1-1.353, -110.Part-1-1.354, -100.Part-1-1.355, -90.Part-1-1.356, -80.Part-1-1.357, -70.Part-1-1.358, -60.Part-1-1.359, -50.Part-1-1.360, -
39、40.Part-1-1.361, -30.Part-1-1.362, -20.Part-1-1.363, -10.Part-1-1.364, 0.* -* * STEP: Step-1*Step, name=Step-1, nlgeom=NO*Geostatic* * LOADS* * Name: Load-1 Type: Gravity*Dload, GRAV, 10., 0., -1.* * OUTPUT REQUESTS*Restart, write, frequency=0* * FIELD OUTPUT: F-Output-1*Output, field, variable=PRES
40、ELECT* * HISTORY OUTPUT: H-Output-1*Output, history, variable=PRESELECT*End Step* -* * STEP: Step-2*Step, name=Step-2, nlgeom=NO, amplitude=RAMP, inc=1000, unsymm=YES*Soils, consolidation, end=PERIOD, utol=100., creep=none1., 180., 0.0018, 180.,*SFLOWPart-1-1.Surf-fot,QD,700Part-1-1.Surf-down,QD,700
41、* * BOUNDARY CONDITIONS* * Name: BC-5 Type: Pore pressure*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, am
42、plitude=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,
43、 8,1.*Boundary, amplitude=Amp-13Part-1-1.1, 8, 8,1.* * OUTPUT REQUESTS*Restart, write, frequency=0* * FIELD OUTPUT: F-Output-1*Output, field, variable=PRESELECT* * HISTORY OUTPUT: H-Output-1*Output, history, variable=PRESELECT*End Step* -* * STEP: Step-3*Step, name=Step-3, nlgeom=NO, amplitude=RAMP, inc=1000, unsymm=YES*Soils, consolidation, end=PERIOD, utol=100., creep=none1., 1800., 0.018, 1800.,*SFLOWPart-1-1.Surf-fot,QD,700Part-1-1.Surf-down,QD,700* * OUTPUT REQUESTS*Restart, write, frequency=0* * FIELD OUTPUT: F-Output-1*Output, field, variable=PRESELECT* *