*NODE, nset=nall 1, -0.05, 0.05, 0.1 2, 0.05, 0.05, 0.1 3, 0.05, 0.05, 0. 4, -0.05, 0.05, 0. 5, -0.05, -0.05, 0.1 6, -0.05, -0.05, 0. 7, 0.05, -0.05, 0.1 8, 0.05, -0.05, 0. 9, -0.05, 0.05, 0.2 10, 0.05, 0.05, 0.2 11, -0.05, -0.05, 0.2 12, 0.05, -0.05, 0.2 13, 0., 0.05, 0.1 14, 0.05, 0.05, 0.05 15, 0., 0.05, 0. 16, -0.05, 0.05, 0.05 17, -0.05, 0., 0.1 18, -0.05, 0., 0. 19, -0.05, -0.05, 0.05 20, 0., -0.05, 0.1 21, 0., -0.05, 0. 22, 0.05, -0.05, 0.05 23, 0.05, 0., 0.1 24, 0.05, 0., 0. 25, -0.05, 0.05, 0.15 26, 0., 0.05, 0.2 27, 0.05, 0.05, 0.15 28, -0.05, -0.05, 0.15 29, -0.05, 0., 0.2 30, 0.05, -0.05, 0.15 31, 0., -0.05, 0.2 32, 0.05, 0., 0.2 33, 0., 0.05, 0.05 34, -0.05, 0., 0.05 35, 0., -0.05, 0.05 36, 0.05, 0., 0.05 37, 0., 0., 0. 38, 0., 0., 0.1 39, 0., 0.05, 0.150000007629395 40, -0.05, 0., 0.15 41, 0., -0.05, 0.150000007629395 42, 0.05, 0., 0.15 43, 0., 0., 0.2 44, 0.020904779303051, -0.0209769512685274, 0.0742083988596052 45, 0.0231573802485947, 0.0230810350915129, 0.0781477452085217 46, -0.0231578170976784, -0.0230808064736723, 0.0781472126222231 47, -0.0209056057503491, 0.020977442309704, 0.0742084834593151 48, -0.0153686161607931, -0.0151072835435189, 0.0402925171861798 49, -0.0201985800089054, 0.021006867269397, 0.0307219242953456 50, 0.0153664493191051, 0.0151094756463619, 0.0402955344946488 51, 0.0201976056018008, -0.0210034692382957, 0.0307233664610405 52, 0., 0., 0.150000007629395 *ELEMENT, TYPE=C3D4 1, 34, 48, 18, 49 2, 16, 18, 4, 49 3, 15, 4, 37, 49 4, 16, 47, 34, 49 5, 18, 6, 21, 48 6, 44, 48, 50, 51 7, 33, 50, 47, 49 8, 19, 46, 35, 48 9, 21, 37, 48, 51 10, 18, 21, 37, 48 11, 37, 48, 50, 49 12, 37, 50, 48, 51 13, 36, 50, 24, 51 14, 38, 47, 45, 50 15, 36, 44, 50, 51 16, 14, 36, 45, 50 17, 19, 21, 6, 48 18, 21, 8, 37, 51 19, 35, 8, 21, 51 20, 8, 24, 37, 51 21, 24, 50, 37, 51 22, 22, 44, 36, 51 23, 16, 1, 34, 47 24, 33, 47, 16, 49 25, 13, 45, 38, 47 26, 34, 17, 46, 47 27, 17, 5, 19, 46 28, 35, 20, 44, 46 29, 35, 46, 44, 48 30, 1, 17, 34, 47 31, 44, 46, 38, 48 32, 14, 3, 36, 50 33, 15, 24, 3, 50 34, 36, 44, 45, 50 35, 33, 45, 47, 50 36, 36, 44, 23, 45 37, 2, 14, 23, 45 38, 14, 36, 23, 45 39, 7, 23, 36, 44 40, 20, 38, 44, 46 41, 15, 37, 50, 49 42, 5, 20, 35, 46 43, 33, 15, 50, 49 44, 15, 37, 24, 50 45, 2, 13, 33, 45 46, 19, 35, 21, 48 47, 3, 24, 36, 50 48, 35, 21, 48, 51 49, 17, 38, 46, 47 50, 34, 47, 46, 48 51, 14, 15, 3, 50 52, 20, 7, 22, 44 53, 19, 5, 35, 46 54, 22, 24, 8, 51 55, 35, 48, 44, 51 56, 34, 6, 18, 48 57, 23, 44, 38, 45 58, 18, 48, 37, 49 59, 33, 45, 13, 47 60, 13, 1, 16, 47 61, 34, 17, 19, 46 62, 33, 16, 4, 49 63, 34, 46, 19, 48 64, 22, 7, 36, 44 65, 38, 46, 47, 48 66, 44, 38, 50, 48 67, 35, 44, 22, 51 68, 33, 15, 14, 50 69, 44, 38, 45, 50 70, 34, 47, 48, 49 71, 35, 20, 22, 44 72, 38, 47, 50, 48 73, 4, 18, 37, 49 74, 2, 33, 14, 45 75, 33, 14, 45, 50 76, 35, 22, 8, 51 77, 33, 13, 16, 47 78, 47, 50, 48, 49 79, 22, 36, 24, 51 80, 33, 4, 15, 49 81, 16, 34, 18, 49 82, 34, 19, 6, 48 *ELEMENT, TYPE=C3D5 83, 23, 7, 20, 38, 44 84, 2, 23, 38, 13, 45 85, 38, 20, 5, 17, 46 86, 13, 38, 17, 1, 47 *ELEMENT, TYPE=C3D8 87, 28, 41, 20, 5, 40, 52, 38, 17 88, 41, 30, 7, 20, 52, 42, 23, 38 89, 41, 31, 12, 30, 52, 43, 32, 42 90, 41, 28, 11, 31, 52, 40, 29, 43 91, 40, 52, 38, 17, 25, 39, 13, 1 92, 52, 42, 23, 38, 39, 27, 2, 13 93, 52, 43, 32, 42, 39, 26, 10, 27 94, 52, 40, 29, 43, 39, 25, 9, 26 *elset, elset=alleles, gen 1, 94, 1 *nset, nset=bot 9, 10, 11, 12, 26, 29, 31, 32, 43 *Solid section, elset=alleles, material=mat1 , *MATERIAL, NAME=mat1 *DENSITY 1000, *ELASTIC 1.e6, 0.2 **-----*POROUS BULK MODULI **----- 1.0e16, 1.0e8 **-----*PERMEABILITY, SPECIFIC=1.0e-6 **----- 5.e-10 *PORE FLUID PRESSURE, FIELD=2 **-----*INITIAL CONDITIONS,TYPE=RATIO **-----NALL,0.6 *INITIAL CONDITIONS, TYPE=FIELD, VARIABLE=1 NALL, 0.111 *INITIAL CONDITIONS, TYPE=FIELD, VARIABLE=2 NALL, 0.222 **----------------------------------------------------------------------------------- *STEP, NAME=Gravity, NLGEOM=YES, UNSYMM=YES, AMPLITUDE=RAMP *STATIC 0.1, 1.0, 1e-5, 1.0 *DLOAD, OP=NEW alleles, GRAV, 9.81, 0, 0, -1 *boundary bot, 3, 3 9, 1, 2 10, 2, 2 11, 1, 1 *FIELD, VARIABLE=2 1, 0.0 2, 0.0 3, 0.0 4, 0.0 5, 0.0 6, 0.0 7, 0.0 8, 0.0 9, 0.0 10,0.0 11,0.0 12,0.0 13,0.0 14,0.0 15,0.0 16,0.0 17,0.0 18,0.0 19,0.0 20,0.0 21,0.0 22,0.0 23,0.0 24,0.0 25,0.0 26,0.0 27,0.0 28,0.0 29,0.0 30,0.0 31,0.0 32,0.0 33,0.0 34,0.0 35,0.0 36,0.0 37,0.0 38,0.0 39,0.0 40,0.0 41,0.0 42,0.0 43,0.0 44,0.0 45,0.0 46,0.0 47,0.0 48,0.0 49,0.0 50,0.0 51,0.0 52,0.0 ** *OUTPUT, FIELD *NODE OUTPUT U, RF, CF *OUTPUT, FIELD *ELEMENT OUTPUT S, FV, POR *END STEP **************************************************************************************** **----------------------------------------------------------------------------------- *STEP, NAME=GravityPlusPorPressure, NLGEOM=YES, UNSYMM=YES, AMPLITUDE=RAMP *STATIC 0.1, 1.0, 1e-5, 1.0 *boundary bot, 3, 3 9, 1, 2 10, 2, 2 11, 1, 1 *FIELD, VARIABLE=2 1, 100.0 2, 100.0 3, 100.0 4, 100.0 5, 100.0 6, 100.0 7, 100.0 8, 100.0 9, 100.0 10,100.0 11,100.0 12,100.0 13,100.0 14,100.0 15,100.0 16,100.0 17,100.0 18,100.0 19,100.0 20,100.0 21,100.0 22,100.0 23,100.0 24,100.0 25,100.0 26,100.0 27,100.0 28,100.0 29,100.0 30,100.0 31,100.0 32,100.0 33,100.0 34,100.0 35,100.0 36,100.0 37,100.0 38,100.0 39,100.0 40,100.0 41,100.0 42,100.0 43,100.0 44,100.0 45,100.0 46,100.0 47,100.0 48,100.0 49,100.0 50,100.0 51,100.0 52,100.0 ** *OUTPUT, FIELD *NODE OUTPUT U, CF, RF *OUTPUT, FIELD *ELEMENT OUTPUT S, FV, POR *END STEP ****************************************************************************************