Abquse模擬庫(kù)水位升降的例子(含完整程序)

1、Abaqus模擬庫(kù)水位升降的例子一含附件Abaqus可以通過(guò)子程序或者編輯關(guān)鍵字來(lái)實(shí)現(xiàn)庫(kù)水位升降的模擬。如下例子，一均質(zhì)土石壩受水位變化時(shí)浸潤(rùn)線分布，設(shè)該土石壩高12m,水位變化如圖：庫(kù)水位隨時(shí)間的變化時(shí)間數(shù)值模擬中實(shí)現(xiàn)這個(gè)問(wèn)題，需要對(duì)邊界條件上加載孔隙水壓力，首先要把時(shí)間定下來(lái)，然后把每個(gè)時(shí)間對(duì)應(yīng)的水位高度定下來(lái)，然后就是寫(xiě)對(duì)應(yīng)的程序了， 關(guān)鍵就在于，需要在各個(gè)位置的節(jié)點(diǎn)處定義不同的幅值曲線。這個(gè)問(wèn)題使用子程序會(huì)很方便，也可以不用，只需要定義一堆關(guān)鍵字吧，但是GUI方式是完全沒(méi)法實(shí)現(xiàn)的。下面為實(shí)現(xiàn)該模擬的關(guān)鍵字（完整文 件見(jiàn)附件）首先定義幅值曲線，以下為例說(shuō)明幅值曲線的含義："Am

2、plitude為關(guān)鍵字，即幅值曲線name=Amp-1幅值曲線的名稱(chēng)定義為Amp-1time=TOTAL TIME時(shí)間有兩種，一種是 step time ,即以該分析步時(shí)間為基準(zhǔn)，并只在該 分析步時(shí)間生效；另一種是total time ,即作用于所有分析步中，以第一個(gè)（initial之后的第一個(gè)分析步）分析步為初始時(shí)間。0., 0., 180., 110., 1980., 110., 2160., 0.上一行為曲線的定義，即表示了下圖的曲線?？讐弘S時(shí)間的變化時(shí)間分別定義上游壩面與水接觸的各個(gè)節(jié)點(diǎn)的幅值曲線："Amplitude, name=Amp-1, time=TOTAL TIME

3、0., 0.,180.,110.,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,

4、 name=Amp-5, time=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.

5、,2160.,-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.0., -100., 180., 10., 1980., 10., 2160., -100.*Amplitude, name=Amp-12, time=TOTAL TIME0., -110., 180., 0., 1980., 0., 2160.,

6、 -110.*Amplitude, name=Amp-13, time=TOTAL TIME0., -120., 180., -10., 1980., -100., 2160., -120.接下來(lái)就是定義各個(gè)節(jié)點(diǎn)的邊界條件 孔壓：*Boundary 為關(guān)鍵字，表示邊界條件amplitude=Amp-1 采用名為Amp-1 的幅值曲線定義邊界條件8, 8 表示孔壓自由Part-1-1.13, 8, 8, 1. Part-1-1.13 為 Part-1-1 部件的編號(hào)為13的節(jié)點(diǎn)；度； 1 表示大小為1*Amp-1 。分別為上游壩面與水接觸的各個(gè)節(jié)點(diǎn)的添加孔壓邊界條件：*Boundary, amp

7、litude=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、, 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.A

9、bquse計(jì)算浸潤(rùn)線結(jié)果如下，由圖可見(jiàn)，庫(kù)水位在第180天上升至11m時(shí)壩體滲流并未穩(wěn)定，在第240天壩體滲流穩(wěn)定，庫(kù)水位下降時(shí)，形成向上游防滲，對(duì)壩坡不利:T=140,。0e+0。-l,LB7e-F02T=40.7-ru, vuue-ruu-l.L46e+02T=130T=180T=210T=240T=1980T=2000+U.UDU£+UU-9.002e+01T=2110T=2160附件：* PARTS*Part, name=Part-1*Element, type=CPE4P*Surface, type=ELEMENT, name=Surf-up*Surface, type=E

10、LEMENT, name=Surf-fot*Surface, type=ELEMENT, name=Surf-down* 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.

11、,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, time=TOTAL TIME0., -40.,180.,70.,1

12、980.,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.,-70.*Amplitude, name=Amp-9, time=TOTAL TIME0., -80.,18

13、0.,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=Amp-13,

14、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=SATURATION0.000232, 0.178850.000561, 0.200180.001442, 0.225730.003736, 0.259480.009699, 0.305810.02478, 0.371390.0635

15、5, 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.* * BOUNDARY CONDITIONS* * Name: BC-1 Type: Displacement/Rotation* * Boundary中涉及PickedSet的要注意，本例約束底面、側(cè)面和頂面所有位移,*Bound

16、ary_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_Picked

17、Set7, 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

18、.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.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, -20.Part-1-1.25, -10.Part

19、-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-1.42, -100.Part-1-1.45,

20、-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, -60.Part-1-1.60, -50.Part-1-1.61, -40.Part-1-1.62, -30.Part-1

21、-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.75, -30.Part-1-1.76, -20.Part-1-1.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, -9

22、0.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.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

23、.100, -40.Part-1-1.101, -30.Part-1-1.102, -20.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, -60.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.

24、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.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.135, -80.Part-1-1.136, -70.Part-1-1.1

25、37, -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.Part-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.Par

26、t-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, -80.Part-1-1.162, -70.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,

27、-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-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

28、-1.189, -60.Part-1-1.190, -50.Part-1-1.191, -40.Part-1-1.192, -30.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

29、.Part-1-1.208, 0.Part-1-1.209, -120.Part-1-1.210, -110.Part-1-1.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.225, -90.Part-1-1.22

30、6, -80.Part-1-1.227, -70.Part-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

31、-1-1.243, -40.Part-1-1.244, -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.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-1.261, -12

32、0.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, -90.Part-1-

33、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.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.294, -50.Part-1-1.295, -40.Part-1-1.296, -30.

34、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.Part-1-1.311, -10.Part-1-1.312, 0.Part-1-1.315,

35、 -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, -110.Part-1-1.328, -100.Part-1-1.329, -90.Part-1-1.330, -80.Part-1-1.331, -70.Part-

36、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.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, -1

37、0.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, -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*G

38、eostatic*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=PRESELECT* HISTORY OUTPUT: H-Output-1*Output, history, variable=PRESELECT*End Step* * STEP: Step-2*Step, name=Step-2, nlgeom=NO, amplit

39、ude=RAMP, inc=1000, unsymm=YES*Soils, consolidation, end=PERIOD, utol=100., creep=none1.1, 180., 0.0018, 180.,*SFLOWPart-1-1.Surf-fot,QD,700Part-1-1.Surf-down,QD,700* * 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, amplitude=Amp-6Part-1-1.8, 8