*HEADING STRESS FREE INCLUSION OF TUNNEL LINER (USING *MODEL CHANGE) Units : N, m *RESTART,WRITE ** *NODE,NSET=SOIL,INPUT=modelchangedemo_node.inp *NSET,NSET=TUN_IX,UNSORTED 101, 102, 103, 104, 105, 106, 2105, 2104, 2103, 2102, 2101 *NSET,NSET=TUN_IY,UNSORTED 100, 101, 102, 103, 104, 105, 106, 2105, 2104, 2103, 2102, 2101, 2100 *NSET,NSET=TUNNEL,GENERATE 100, 106 2100, 2105 *NODE 4000, 0.00, -16.0 4012, 0.00, -24.0 *NGEN,NSET=LINER,LINE=C 4000, 4012, 1, ,0.0, -20.0, 0.0, 0.0, 0.0, -1.0 *NSET,NSET=LIN_JX,GENERATE 4001, 4011 *NSET,NSET=LIN_JY,GENERATE 4000, 4012 *NSET,NSET=RHS,GENERATE 8, 3108, 10 *NSET,NSET=BOT,GENERATE 3100, 3108 *NSET,NSET=XSYMM,GENERATE 100, 150, 10 200, 1000, 100 2100, 2150, 10 2200, 3200, 100 *NSET,NSET=XSYMM 4000, 4012, 4100, 4112 *NODE,NSET=NINF 209, 60.00, -20.00 309, 60.00, -17.77 409, 60.00, -14.98 509, 60.00, -11.50 609, 60.00, -10.10 709, 60.00, -8.347 809, 60.00, -6.158 909, 60.00, -3.420 1009, 60.00, 0.000 2309, 60.00, -22.03 2409, 60.00, -25.02 2509, 60.00, -28.50 2609, 60.00, -29.90 2709, 60.00, -31.65 2809, 60.00, -33.84 2909, 60.00, -36.58 3009, 60.00, -40.00 3109, 60.00, -76.00 3200, 0.000, -76.00 3201, 2.245, -76.00 3202, 5.052, -76.00 3203, 8.560, -76.00 3204, 11.17, -76.00 3205, 14.44, -76.00 3206, 18.52, -76.00 3207, 23.62, -76.00 ** *ELEMENT,TYPE=CPE4,ELSET=SOIL, INPUT=modelchangedemo_element.inp *ELEMENT,TYPE=B21 4000, 4000, 4001 *ELGEN,ELSET=LINER 4000, 12 *ELSET,ELSET=SURFACE,GENERATE 900 , 907 , 1 *ELEMENT,TYPE=CINPE4,ELSET=ELINF 3100,3101,3100,3200,3201 3101,3102,3101,3201,3202 3102,3103,3102,3202,3203 3103,3104,3103,3203,3204 3104,3105,3104,3204,3205 3105,3106,3105,3205,3206 3106,3107,3106,3206,3207 3107,3108,3107,3207,3109 3108,3008,3108,3109,3009 2908,2908,3008,3009,2909 2808,2808,2908,2909,2809 2708,2708,2808,2809,2709 2608,2608,2708,2709,2609 2508,2508,2608,2609,2509 2408,2408,2508,2509,2409 2308,2308,2408,2409,2309 2208,208,2308,2309,209 208,308,208,209,309 308,408,308,309,409 408,508,408,409,509 508,608,508,509,609 608,708,608,609,709 708,808,708,709,809 808,908,808,809,909 908,1008,908,909,1009 *ELSET,ELSET=EALL SOIL,ELINF ** *BEAM SECTION,MATERIAL=CONCRETE,SECTION=RECT, ELSET=LINER 1.0 , 0.15 0.0 , 0.0 , -1.0 *SOLID SECTION,MATERIAL=CLAY,ELSET=SOIL *SOLID SECTION,MATERIAL=CLAY,ELSET=ELINF *MATERIAL,NAME=CLAY *ELASTIC,TYPE=ISOTROPIC 0.2E9, 0.2 *MATERIAL,NAME=CONCRETE *ELASTIC,TYPE=ISOTROPIC 19.0E9 , 0.2 ** *EQUATION 2, TUN_IX, 1, 1.0, LIN_JX, 1, -1.0 2, TUN_IY, 2, 1.0, LIN_JY, 2, -1.0 *INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC EALL, 0.0, 0.0, -1.52E6, -76.00, 0.5 *AMPLITUDE,NAME=RELAX,TIME=TOTAL TIME 0.0, 1.0, 1.0, 1.0, 2.0, 0.6, 3.0, 0.6, 4.0, 0.0 ** -------------------------------------- *STEP step 1: add initial stress state *STATIC *DLOAD SOIL, BY, -20.0E3 *MODEL CHANGE,REMOVE LINER, *CLOAD,AMPLITUDE=RELAX 100, 1, 5.4086E+04 101, 1, 4.3918E+04 102, 1, 8.6901E+04 103, 1, 1.2732E+05 104, 1, 1.6185E+05 105, 1, 1.8949E+05 106, 1, 2.0701E+05 2100, 1, 8.2710E+04 2101, 1, 6.3287E+04 2102, 1, 1.2031E+05 2103, 1, 1.6549E+05 2104, 1, 1.9676E+05 2105, 1, 2.1052E+05 100, 2, 1.6652E+05 101, 2, 3.2459E+05 102, 2, 2.9838E+05 103, 2, 2.5160E+05 104, 2, 1.8487E+05 105, 2, 9.9587E+04 106, 2, -2142. 2100, 2, -2.4756E+05 2101, 2, -4.7534E+05 2102, 2, -4.1880E+05 2103, 2, -3.3399E+05 2104, 2, -2.2922E+05 2105, 2, -1.1476E+05 *BOUNDARY XSYMM, 1 XSYMM, 6 *NODE PRINT,FREQUENCY=100,NSET=TUNNEL U, RF *EL PRINT,FREQUENCY=100,ELSET=LINER, POSITION=AVERAGED AT NODES S, *END STEP ** -------------------------------------- *STEP step 2: relax tunnel stress 40 % *STATIC *END STEP ** -------------------------------------- *STEP step 3: add liner stress free *STATIC *MODEL CHANGE,ADD LINER, *END STEP ** -------------------------------------- *STEP,IN=100 step 4: relax tunnel stress to zero & apply surface load *STATIC *DLOAD SURFACE, P3, 50.0E3 *EL FILE,FREQUENCY=100,ELSET=LINER, POSITION=AVERAGED AT NODES S, *END STEP