modeling of sheet pile wall

[sanjianke]

Hi, there,

I am building a soil model with sheet pile wall in a rigid container of shake table test under earthquake engineering. Basically, I am following the example "excavation supported by cantilevered sheet pile wall". I am not sure about the function of "InitialStateWrapper" and its applicability for dynamic analysis. When I tried to run the model, I always got the error message as follows. Please help me out. I will appreciate it very much.

WARNING FE_Element::SetID() - numDOF and number of dof at the DOF_Groups
WARNING Super LU::solve (void) - Error 931 returned in factorization dgstrf ()
WARNING NewtonRaphson::solveCurrentStep() -theLinearSysOfEqn failed in solve ()
DirectIntegrationAnalysis::analyze() -the Algorithm failed at time 1
Opensees > anaylze failed, returned: -3 error flag

I have copied my simplified input tcl file as following.

wipe

##-----------------------------------------------------------------------------------------
## 2. CREATE PORE PRESSURE NODES AND FIXITIES
##-----------------------------------------------------------------------------------------

model BasicBuilder -ndm 2 -ndf 3
set ppNodesInfo [open ppNodesInfo.dat w]
## define pore pressure nodes:
node 1 0 0
node 2 0.5 0
node 3 1 0
node 4 1.5 0
node 5 2 0
node 6 2.5 0
node 7 3 0
node 8 3.5 0
node 9 4 0
node 10 4.5 0
node 11 4.586 0
node 12 5 0
node 13 5.6805 0
node 14 6 0
node 15 6.5 0
node 16 6.9195 0
node 17 7.5 0
node 18 8 0
node 19 8.5 0
node 20 9 0
node 21 9.5 0
node 22 10 0
node 23 10.5 0
node 24 11 0
node 25 11.5 0
node 26 12 0
node 27 12.5 0
node 28 13 0
node 29 13.5 0
node 30 14 0
node 31 14.5 0
node 32 15 0
node 33 15.5 0
node 34 16 0
node 35 0 0.4
node 36 0.5 0.4
node 37 1 0.4
node 38 1.5 0.4
node 39 2 0.4
node 40 2.5 0.4
node 41 3 0.4
node 42 3.5 0.4
node 43 4 0.4
node 44 4.5 0.4
node 45 4.586 0.4
node 46 5 0.4
node 47 5.6805 0.4
node 48 6 0.4
node 49 6.5 0.4
node 50 6.9195 0.4
node 51 7.5 0.4
node 52 8 0.4
node 53 8.5 0.4
node 54 9 0.4
node 55 9.5 0.4
node 56 10 0.4
node 57 10.5 0.4
node 58 11 0.4
node 59 11.5 0.4
node 60 12 0.4
node 61 12.5 0.4
node 62 13 0.4
node 63 13.5 0.4
node 64 14 0.4
node 65 14.5 0.4
node 66 15 0.4
node 67 15.5 0.4
node 68 16 0.4
node 69 0 0.8
node 70 0.5 0.8
node 71 1 0.8
node 72 1.5 0.8
node 73 2 0.8
node 74 2.5 0.8
node 75 3 0.8
node 76 3.5 0.8
node 77 4 0.8
node 78 4.5 0.8
node 79 4.586 0.8
node 80 5 0.8
node 81 5.6805 0.8
node 82 6 0.8
node 83 6.5 0.8
node 84 6.9195 0.8
node 85 7.5 0.8
node 86 8 0.8
node 87 8.5 0.8
node 88 9 0.8
node 89 9.5 0.8
node 90 10 0.8
node 91 10.5 0.8
node 92 11 0.8
node 93 11.5 0.8
node 94 12 0.8
node 95 12.5 0.8
node 96 13 0.8
node 97 13.5 0.8
node 98 14 0.8
node 99 14.5 0.8
node 100 15 0.8
node 101 15.5 0.8
node 102 16 0.8
node 103 0 1.2
node 104 0.5 1.2
node 105 1 1.2
node 106 1.5 1.2
node 107 2 1.2
node 108 2.5 1.2
node 109 3 1.2
node 110 3.5 1.2
node 111 4 1.2
node 112 4.5 1.2
node 113 4.586 1.2
node 114 5 1.2
node 115 5.6805 1.2
node 116 6 1.2
node 117 6.5 1.2
node 118 6.9195 1.2
node 119 7.5 1.2
node 120 8 1.2
node 121 8.5 1.2
node 122 9 1.2
node 123 9.5 1.2
node 124 10 1.2
node 125 10.5 1.2
node 126 11 1.2
node 127 11.5 1.2
node 128 12 1.2
node 129 12.5 1.2
node 130 13 1.2
node 131 13.5 1.2
node 132 14 1.2
node 133 14.5 1.2
node 134 15 1.2
node 135 15.5 1.2
node 136 16 1.2
node 137 0 1.6
node 138 0.5 1.6
node 139 1 1.6
node 140 1.5 1.6
node 141 2 1.6
node 142 2.5 1.6
node 143 3 1.6
node 144 3.5 1.6
node 145 4 1.6
node 146 4.5 1.6
node 147 4.586 1.6
node 148 5 1.6
node 149 5.6805 1.6
node 150 6 1.6
node 151 6.5 1.6
node 152 6.9195 1.6
node 153 7.5 1.6
node 154 8 1.6
node 155 8.5 1.6
node 156 9 1.6
node 157 9.5 1.6
node 158 10 1.6
node 159 10.5 1.6
node 160 11 1.6
node 161 11.5 1.6
node 162 12 1.6
node 163 12.5 1.6
node 164 13 1.6
node 165 13.5 1.6
node 166 14 1.6
node 167 14.5 1.6
node 168 15 1.6
node 169 15.5 1.6
node 170 16 1.6
node 171 0 2
node 172 0.5 2
node 173 1 2
node 174 1.5 2
node 175 2 2
node 176 2.5 2
node 177 3 2
node 178 3.5 2
node 179 4 2
node 180 4.5 2
node 181 4.586 2
node 182 5 2
node 183 5.6805 2
node 184 6 2
node 185 6.5 2
node 186 6.9195 2
node 187 7.5 2
node 188 8 2
node 189 8.5 2
node 190 9 2
node 191 9.5 2
node 192 10 2
node 193 10.5 2
node 194 11 2
node 195 11.5 2
node 196 12 2
node 197 12.5 2
node 198 13 2
node 199 13.5 2
node 200 14 2
node 201 14.5 2
node 202 15 2
node 203 15.5 2
node 204 16 2
node 205 0 2.4
node 206 0.5 2.4
node 207 1 2.4
node 208 1.5 2.4
node 209 2 2.4
node 210 2.5 2.4
node 211 3 2.4
node 212 3.5 2.4
node 213 4 2.4
node 214 4.5 2.4
node 215 4.586 2.4
node 216 5 2.4
node 217 5.6805 2.4
node 218 6 2.4
node 219 6.5 2.4
node 220 6.9195 2.4
node 221 7.5 2.4
node 222 8 2.4
node 223 8.5 2.4
node 224 9 2.4
node 225 9.5 2.4
node 226 10 2.4
node 227 10.5 2.4
node 228 11 2.4
node 229 11.5 2.4
node 230 12 2.4
node 231 12.5 2.4
node 232 13 2.4
node 233 13.5 2.4
node 234 14 2.4
node 235 14.5 2.4
node 236 15 2.4
node 237 15.5 2.4
node 238 16 2.4
node 239 0 2.8
node 240 0.5 2.8
node 241 1 2.8
node 242 1.5 2.8
node 243 2 2.8
node 244 2.5 2.8
node 245 3 2.8
node 246 3.5 2.8
node 247 4 2.8
node 248 4.5 2.8
node 249 4.586 2.8
node 250 5 2.8
node 251 5.6805 2.8
node 252 6 2.8
node 253 6.5 2.8
node 254 6.9195 2.8
node 255 7.5 2.8
node 256 8 2.8
node 257 8.5 2.8
node 258 9 2.8
node 259 9.5 2.8
node 260 10 2.8
node 261 10.5 2.8
node 262 11 2.8
node 263 11.5 2.8
node 264 12 2.8
node 265 12.5 2.8
node 266 13 2.8
node 267 13.5 2.8
node 268 14 2.8
node 269 14.5 2.8
node 270 15 2.8
node 271 15.5 2.8
node 272 16 2.8
node 273 0 3.2
node 274 0.5 3.2
node 275 1 3.2
node 276 1.5 3.2
node 277 2 3.2
node 278 2.5 3.2
node 279 3 3.2
node 280 3.5 3.2
node 281 4 3.2
node 282 4.5 3.2
node 283 4.586 3.2
node 284 5 3.2
node 285 5.6805 3.2
node 286 6 3.2
node 287 6.5 3.2
node 288 6.9195 3.2
node 289 7.5 3.2
node 290 8 3.2
node 291 8.5 3.2
node 292 9 3.2
node 293 9.5 3.2
node 294 10 3.2
node 295 10.5 3.2
node 296 11 3.2
node 297 11.5 3.2
node 298 12 3.2
node 299 12.5 3.2
node 300 13 3.2
node 301 13.5 3.2
node 302 14 3.2
node 303 14.5 3.2
node 304 15 3.2
node 305 15.5 3.2
node 306 16 3.2
node 307 4.586 3.6
node 308 5 3.6
node 309 5.6805 3.6
node 310 6 3.6
node 311 6.5 3.6
node 312 7 3.6
node 313 7.5 3.6
node 314 8 3.6
node 315 8.5 3.6
node 316 9 3.6
node 317 9.5 3.6
node 318 10 3.6
node 319 10.5 3.6
node 320 11 3.6
node 321 11.5 3.6
node 322 12 3.6
node 323 12.5 3.6
node 324 13 3.6
node 325 13.5 3.6
node 326 14 3.6
node 327 14.5 3.6
node 328 15 3.6
node 329 15.5 3.6
node 330 16 3.6
node 331 4.586 4
node 332 5 4
node 333 5.6805 4
node 334 6 4
node 335 6.5 4
node 336 6.9195 4
node 337 7.5 4
node 338 8 4
node 339 8.5 4
node 340 9 4
node 341 9.5 4
node 342 10 4
node 343 10.5 4
node 344 11 4
node 345 11.5 4
node 346 12 4
node 347 12.5 4
node 348 13 4
node 349 13.5 4
node 350 14 4
node 351 14.5 4
node 352 15 4
node 353 15.5 4
node 354 16 4
node 355 4.586 4.25
node 356 5 4.25
node 357 5.6805 4.25
node 358 6 4.25
node 359 6.5 4.25
node 360 6.9195 4.25
node 361 7.5 4.25
node 362 8 4.25
node 363 8.5 4.25
node 364 9 4.25
node 365 9.5 4.25
node 366 10 4.25
node 367 10.5 4.25
node 368 11 4.25
node 369 11.5 4.25
node 370 12 4.25
node 371 12.5 4.25
node 372 13 4.25
node 373 13.5 4.25
node 374 14 4.25
node 375 14.5 4.25
node 376 15 4.25
node 377 15.5 4.25
node 378 16 4.25
node 379 4.586 4.5
node 380 5 4.5
node 381 5.6805 4.5
node 382 6.9195 4.5
node 383 7.5 4.5
node 384 8 4.5
node 385 8.5 4.5
node 386 9 4.5
node 387 9.5 4.5
node 388 10 4.5
node 389 10.5 4.5
node 390 11 4.5
node 391 11.5 4.5
node 392 12 4.5
node 393 12.5 4.5
node 394 13 4.5
node 395 13.5 4.5
node 396 14 4.5
node 397 14.5 4.5
node 398 15 4.5
node 399 15.5 4.5
node 400 16 4.5


## output nodal information to data file
puts $ppNodesInfo "1 0 0"
puts $ppNodesInfo "2 0.5 0"
puts $ppNodesInfo "3 1 0"
puts $ppNodesInfo "4 1.5 0"
puts $ppNodesInfo "5 2 0"
puts $ppNodesInfo "6 2.5 0"
puts $ppNodesInfo "7 3 0"
puts $ppNodesInfo "8 3.5 0"
puts $ppNodesInfo "9 4 0"
puts $ppNodesInfo "10 4.5 0"
puts $ppNodesInfo "11 4.586 0"
puts $ppNodesInfo "12 5 0"
puts $ppNodesInfo "13 5.5 0"
puts $ppNodesInfo "14 6 0"
puts $ppNodesInfo "15 6.5 0"
puts $ppNodesInfo "16 7 0"
puts $ppNodesInfo "17 7.5 0"
puts $ppNodesInfo "18 8 0"
puts $ppNodesInfo "19 8.5 0"
puts $ppNodesInfo "20 9 0"
puts $ppNodesInfo "21 9.5 0"
puts $ppNodesInfo "22 10 0"
puts $ppNodesInfo "23 10.5 0"
puts $ppNodesInfo "24 11 0"
puts $ppNodesInfo "25 11.5 0"
puts $ppNodesInfo "26 12 0"
puts $ppNodesInfo "27 12.5 0"
puts $ppNodesInfo "28 13 0"
puts $ppNodesInfo "29 13.5 0"
puts $ppNodesInfo "30 14 0"
puts $ppNodesInfo "31 14.5 0"
puts $ppNodesInfo "32 15 0"
puts $ppNodesInfo "33 15.5 0"
puts $ppNodesInfo "34 16 0"
puts $ppNodesInfo "35 0 0.4"
puts $ppNodesInfo "36 0.5 0.4"
puts $ppNodesInfo "37 1 0.4"
puts $ppNodesInfo "38 1.5 0.4"
puts $ppNodesInfo "39 2 0.4"
puts $ppNodesInfo "40 2.5 0.4"
puts $ppNodesInfo "41 3 0.4"
puts $ppNodesInfo "42 3.5 0.4"
puts $ppNodesInfo "43 4 0.4"
puts $ppNodesInfo "44 4.5 0.4"
puts $ppNodesInfo "45 4.586 0.4"
puts $ppNodesInfo "46 5 0.4"
puts $ppNodesInfo "47 5.6805 0.4"
puts $ppNodesInfo "48 6 0.4"
puts $ppNodesInfo "49 6.5 0.4"
puts $ppNodesInfo "50 6.9195 0.4"
puts $ppNodesInfo "51 7.5 0.4"
puts $ppNodesInfo "52 8 0.4"
puts $ppNodesInfo "53 8.5 0.4"
puts $ppNodesInfo "54 9 0.4"
puts $ppNodesInfo "55 9.5 0.4"
puts $ppNodesInfo "56 10 0.4"
puts $ppNodesInfo "57 10.5 0.4"
puts $ppNodesInfo "58 11 0.4"
puts $ppNodesInfo "59 11.5 0.4"
puts $ppNodesInfo "60 12 0.4"
puts $ppNodesInfo "61 12.5 0.4"
puts $ppNodesInfo "62 13 0.4"
puts $ppNodesInfo "63 13.5 0.4"
puts $ppNodesInfo "64 14 0.4"
puts $ppNodesInfo "65 14.5 0.4"
puts $ppNodesInfo "66 15 0.4"
puts $ppNodesInfo "67 15.5 0.4"
puts $ppNodesInfo "68 16 0.4"
puts $ppNodesInfo "69 0 0.8"
puts $ppNodesInfo "70 0.5 0.8"
puts $ppNodesInfo "71 1 0.8"
puts $ppNodesInfo "72 1.5 0.8"
puts $ppNodesInfo "73 2 0.8"
puts $ppNodesInfo "74 2.5 0.8"
puts $ppNodesInfo "75 3 0.8"
puts $ppNodesInfo "76 3.5 0.8"
puts $ppNodesInfo "77 4 0.8"
puts $ppNodesInfo "78 4.5 0.8"
puts $ppNodesInfo "79 4.586 0.8"
puts $ppNodesInfo "80 5 0.8"
puts $ppNodesInfo "81 5.6805 0.8"
puts $ppNodesInfo "82 6 0.8"
puts $ppNodesInfo "83 6.5 0.8"
puts $ppNodesInfo "84 6.9195 0.8"
puts $ppNodesInfo "85 7.5 0.8"
puts $ppNodesInfo "86 8 0.8"
puts $ppNodesInfo "87 8.5 0.8"
puts $ppNodesInfo "88 9 0.8"
puts $ppNodesInfo "89 9.5 0.8"
puts $ppNodesInfo "90 10 0.8"
puts $ppNodesInfo "91 10.5 0.8"
puts $ppNodesInfo "92 11 0.8"
puts $ppNodesInfo "93 11.5 0.8"
puts $ppNodesInfo "94 12 0.8"
puts $ppNodesInfo "95 12.5 0.8"
puts $ppNodesInfo "96 13 0.8"
puts $ppNodesInfo "97 13.5 0.8"
puts $ppNodesInfo "98 14 0.8"
puts $ppNodesInfo "99 14.5 0.8"
puts $ppNodesInfo "100 15 0.8"
puts $ppNodesInfo "101 15.5 0.8"
puts $ppNodesInfo "102 16 0.8"
puts $ppNodesInfo "103 0 1.2"
puts $ppNodesInfo "104 0.5 1.2"
puts $ppNodesInfo "105 1 1.2"
puts $ppNodesInfo "106 1.5 1.2"
puts $ppNodesInfo "107 2 1.2"
puts $ppNodesInfo "108 2.5 1.2"
puts $ppNodesInfo "109 3 1.2"
puts $ppNodesInfo "110 3.5 1.2"
puts $ppNodesInfo "111 4 1.2"
puts $ppNodesInfo "112 4.5 1.2"
puts $ppNodesInfo "113 4.586 1.2"
puts $ppNodesInfo "114 5 1.2"
puts $ppNodesInfo "115 5.6805 1.2"
puts $ppNodesInfo "116 6 1.2"
puts $ppNodesInfo "117 6.5 1.2"
puts $ppNodesInfo "118 6.9195 1.2"
puts $ppNodesInfo "119 7.5 1.2"
puts $ppNodesInfo "120 8 1.2"
puts $ppNodesInfo "121 8.5 1.2"
puts $ppNodesInfo "122 9 1.2"
puts $ppNodesInfo "123 9.5 1.2"
puts $ppNodesInfo "124 10 1.2"
puts $ppNodesInfo "125 10.5 1.2"
puts $ppNodesInfo "126 11 1.2"
puts $ppNodesInfo "127 11.5 1.2"
puts $ppNodesInfo "128 12 1.2"
puts $ppNodesInfo "129 12.5 1.2"
puts $ppNodesInfo "130 13 1.2"
puts $ppNodesInfo "131 13.5 1.2"
puts $ppNodesInfo "132 14 1.2"
puts $ppNodesInfo "133 14.5 1.2"
puts $ppNodesInfo "134 15 1.2"
puts $ppNodesInfo "135 15.5 1.2"
puts $ppNodesInfo "136 16 1.2"
puts $ppNodesInfo "137 0 1.6"
puts $ppNodesInfo "138 0.5 1.6"
puts $ppNodesInfo "139 1 1.6"
puts $ppNodesInfo "140 1.5 1.6"
puts $ppNodesInfo "141 2 1.6"
puts $ppNodesInfo "142 2.5 1.6"
puts $ppNodesInfo "143 3 1.6"
puts $ppNodesInfo "144 3.5 1.6"
puts $ppNodesInfo "145 4 1.6"
puts $ppNodesInfo "146 4.5 1.6"
puts $ppNodesInfo "147 4.586 1.6"
puts $ppNodesInfo "148 5 1.6"
puts $ppNodesInfo "149 5.6805 1.6"
puts $ppNodesInfo "150 6 1.6"
puts $ppNodesInfo "151 6.5 1.6"
puts $ppNodesInfo "152 6.9195 1.6"
puts $ppNodesInfo "153 7.5 1.6"
puts $ppNodesInfo "154 8 1.6"
puts $ppNodesInfo "155 8.5 1.6"
puts $ppNodesInfo "156 9 1.6"
puts $ppNodesInfo "157 9.5 1.6"
puts $ppNodesInfo "158 10 1.6"
puts $ppNodesInfo "159 10.5 1.6"
puts $ppNodesInfo "160 11 1.6"
puts $ppNodesInfo "161 11.5 1.6"
puts $ppNodesInfo "162 12 1.6"
puts $ppNodesInfo "163 12.5 1.6"
puts $ppNodesInfo "164 13 1.6"
puts $ppNodesInfo "165 13.5 1.6"
puts $ppNodesInfo "166 14 1.6"
puts $ppNodesInfo "167 14.5 1.6"
puts $ppNodesInfo "168 15 1.6"
puts $ppNodesInfo "169 15.5 1.6"
puts $ppNodesInfo "170 16 1.6"
puts $ppNodesInfo "171 0 2"
puts $ppNodesInfo "172 0.5 2"
puts $ppNodesInfo "173 1 2"
puts $ppNodesInfo "174 1.5 2"
puts $ppNodesInfo "175 2 2"
puts $ppNodesInfo "176 2.5 2"
puts $ppNodesInfo "177 3 2"
puts $ppNodesInfo "178 3.5 2"
puts $ppNodesInfo "179 4 2"
puts $ppNodesInfo "180 4.5 2"
puts $ppNodesInfo "181 4.586 2"
puts $ppNodesInfo "182 5 2"
puts $ppNodesInfo "183 5.6805 2"
puts $ppNodesInfo "184 6 2"
puts $ppNodesInfo "185 6.5 2"
puts $ppNodesInfo "186 6.9195 2"
puts $ppNodesInfo "187 7.5 2"
puts $ppNodesInfo "188 8 2"
puts $ppNodesInfo "189 8.5 2"
puts $ppNodesInfo "190 9 2"
puts $ppNodesInfo "191 9.5 2"
puts $ppNodesInfo "192 10 2"
puts $ppNodesInfo "193 10.5 2"
puts $ppNodesInfo "194 11 2"
puts $ppNodesInfo "195 11.5 2"
puts $ppNodesInfo "196 12 2"
puts $ppNodesInfo "197 12.5 2"
puts $ppNodesInfo "198 13 2"
puts $ppNodesInfo "199 13.5 2"
puts $ppNodesInfo "200 14 2"
puts $ppNodesInfo "201 14.5 2"
puts $ppNodesInfo "202 15 2"
puts $ppNodesInfo "203 15.5 2"
puts $ppNodesInfo "204 16 2"
puts $ppNodesInfo "205 0 2.4"
puts $ppNodesInfo "206 0.5 2.4"
puts $ppNodesInfo "207 1 2.4"
puts $ppNodesInfo "208 1.5 2.4"
puts $ppNodesInfo "209 2 2.4"
puts $ppNodesInfo "210 2.5 2.4"
puts $ppNodesInfo "211 3 2.4"
puts $ppNodesInfo "212 3.5 2.4"
puts $ppNodesInfo "213 4 2.4"
puts $ppNodesInfo "214 4.5 2.4"
puts $ppNodesInfo "215 4.586 2.4"
puts $ppNodesInfo "216 5 2.4"
puts $ppNodesInfo "217 5.6805 2.4"
puts $ppNodesInfo "218 6 2.4"
puts $ppNodesInfo "219 6.5 2.4"
puts $ppNodesInfo "220 6.9195 2.4"
puts $ppNodesInfo "221 7.5 2.4"
puts $ppNodesInfo "222 8 2.4"
puts $ppNodesInfo "223 8.5 2.4"
puts $ppNodesInfo "224 9 2.4"
puts $ppNodesInfo "225 9.5 2.4"
puts $ppNodesInfo "226 10 2.4"
puts $ppNodesInfo "227 10.5 2.4"
puts $ppNodesInfo "228 11 2.4"
puts $ppNodesInfo "229 11.5 2.4"
puts $ppNodesInfo "230 12 2.4"
puts $ppNodesInfo "231 12.5 2.4"
puts $ppNodesInfo "232 13 2.4"
puts $ppNodesInfo "233 13.5 2.4"
puts $ppNodesInfo "234 14 2.4"
puts $ppNodesInfo "235 14.5 2.4"
puts $ppNodesInfo "236 15 2.4"
puts $ppNodesInfo "237 15.5 2.4"
puts $ppNodesInfo "238 16 2.4"
puts $ppNodesInfo "239 0 2.8"
puts $ppNodesInfo "240 0.5 2.8"
puts $ppNodesInfo "241 1 2.8"
puts $ppNodesInfo "242 1.5 2.8"
puts $ppNodesInfo "243 2 2.8"
puts $ppNodesInfo "244 2.5 2.8"
puts $ppNodesInfo "245 3 2.8"
puts $ppNodesInfo "246 3.5 2.8"
puts $ppNodesInfo "247 4 2.8"
puts $ppNodesInfo "248 4.5 2.8"
puts $ppNodesInfo "249 4.586 2.8"
puts $ppNodesInfo "250 5 2.8"
puts $ppNodesInfo "251 5.6805 2.8"
puts $ppNodesInfo "252 6 2.8"
puts $ppNodesInfo "253 6.5 2.8"
puts $ppNodesInfo "254 6.9195 2.8"
puts $ppNodesInfo "255 7.5 2.8"
puts $ppNodesInfo "256 8 2.8"
puts $ppNodesInfo "257 8.5 2.8"
puts $ppNodesInfo "258 9 2.8"
puts $ppNodesInfo "259 9.5 2.8"
puts $ppNodesInfo "260 10 2.8"
puts $ppNodesInfo "261 10.5 2.8"
puts $ppNodesInfo "262 11 2.8"
puts $ppNodesInfo "263 11.5 2.8"
puts $ppNodesInfo "264 12 2.8"
puts $ppNodesInfo "265 12.5 2.8"
puts $ppNodesInfo "266 13 2.8"
puts $ppNodesInfo "267 13.5 2.8"
puts $ppNodesInfo "268 14 2.8"
puts $ppNodesInfo "269 14.5 2.8"
puts $ppNodesInfo "270 15 2.8"
puts $ppNodesInfo "271 15.5 2.8"
puts $ppNodesInfo "272 16 2.8"
puts $ppNodesInfo "273 0 3.2"
puts $ppNodesInfo "274 0.5 3.2"
puts $ppNodesInfo "275 1 3.2"
puts $ppNodesInfo "276 1.5 3.2"
puts $ppNodesInfo "277 2 3.2"
puts $ppNodesInfo "278 2.5 3.2"
puts $ppNodesInfo "279 3 3.2"
puts $ppNodesInfo "280 3.5 3.2"
puts $ppNodesInfo "281 4 3.2"
puts $ppNodesInfo "282 4.5 3.2"
puts $ppNodesInfo "283 4.586 3.2"
puts $ppNodesInfo "284 5 3.2"
puts $ppNodesInfo "285 5.6805 3.2"
puts $ppNodesInfo "286 6 3.2"
puts $ppNodesInfo "287 6.5 3.2"
puts $ppNodesInfo "288 6.9195 3.2"
puts $ppNodesInfo "289 7.5 3.2"
puts $ppNodesInfo "290 8 3.2"
puts $ppNodesInfo "291 8.5 3.2"
puts $ppNodesInfo "292 9 3.2"
puts $ppNodesInfo "293 9.5 3.2"
puts $ppNodesInfo "294 10 3.2"
puts $ppNodesInfo "295 10.5 3.2"
puts $ppNodesInfo "296 11 3.2"
puts $ppNodesInfo "297 11.5 3.2"
puts $ppNodesInfo "298 12 3.2"
puts $ppNodesInfo "299 12.5 3.2"
puts $ppNodesInfo "300 13 3.2"
puts $ppNodesInfo "301 13.5 3.2"
puts $ppNodesInfo "302 14 3.2"
puts $ppNodesInfo "303 14.5 3.2"
puts $ppNodesInfo "304 15 3.2"
puts $ppNodesInfo "305 15.5 3.2"
puts $ppNodesInfo "306 16 3.2"
puts $ppNodesInfo "307 4.586 3.6"
puts $ppNodesInfo "308 5 3.6"
puts $ppNodesInfo "309 5.6805 3.6"
puts $ppNodesInfo "310 6 3.6"
puts $ppNodesInfo "311 6.5 3.6"
puts $ppNodesInfo "312 6.9195 3.6"
puts $ppNodesInfo "313 7.5 3.6"
puts $ppNodesInfo "314 8 3.6"
puts $ppNodesInfo "315 8.5 3.6"
puts $ppNodesInfo "316 9 3.6"
puts $ppNodesInfo "317 9.5 3.6"
puts $ppNodesInfo "318 10 3.6"
puts $ppNodesInfo "319 10.5 3.6"
puts $ppNodesInfo "320 11 3.6"
puts $ppNodesInfo "321 11.5 3.6"
puts $ppNodesInfo "322 12 3.6"
puts $ppNodesInfo "323 12.5 3.6"
puts $ppNodesInfo "324 13 3.6"
puts $ppNodesInfo "325 13.5 3.6"
puts $ppNodesInfo "326 14 3.6"
puts $ppNodesInfo "327 14.5 3.6"
puts $ppNodesInfo "328 15 3.6"
puts $ppNodesInfo "329 15.5 3.6"
puts $ppNodesInfo "330 16 3.6"
puts $ppNodesInfo "331 4.586 4"
puts $ppNodesInfo "332 5 4"
puts $ppNodesInfo "333 5.6805 4"
puts $ppNodesInfo "334 6 4"
puts $ppNodesInfo "335 6.5 4"
puts $ppNodesInfo "336 6.9195 4"
puts $ppNodesInfo "337 7.5 4"
puts $ppNodesInfo "338 8 4"
puts $ppNodesInfo "339 8.5 4"
puts $ppNodesInfo "340 9 4"
puts $ppNodesInfo "341 9.5 4"
puts $ppNodesInfo "342 10 4"
puts $ppNodesInfo "343 10.5 4"
puts $ppNodesInfo "344 11 4"
puts $ppNodesInfo "345 11.5 4"
puts $ppNodesInfo "346 12 4"
puts $ppNodesInfo "347 12.5 4"
puts $ppNodesInfo "348 13 4"
puts $ppNodesInfo "349 13.5 4"
puts $ppNodesInfo "350 14 4"
puts $ppNodesInfo "351 14.5 4"
puts $ppNodesInfo "352 15 4"
puts $ppNodesInfo "353 15.5 4"
puts $ppNodesInfo "354 16 4"
puts $ppNodesInfo "355 4.586 4.25"
puts $ppNodesInfo "356 5 4.25"
puts $ppNodesInfo "357 5.6805 4.25"
puts $ppNodesInfo "358 6 4.25"
puts $ppNodesInfo "359 6.5 4.25"
puts $ppNodesInfo "360 6.9195 4.25"
puts $ppNodesInfo "361 7.5 4.25"
puts $ppNodesInfo "362 8 4.25"
puts $ppNodesInfo "363 8.5 4.25"
puts $ppNodesInfo "364 9 4.25"
puts $ppNodesInfo "365 9.5 4.25"
puts $ppNodesInfo "366 10 4.25"
puts $ppNodesInfo "367 10.5 4.25"
puts $ppNodesInfo "368 11 4.25"
puts $ppNodesInfo "369 11.5 4.25"
puts $ppNodesInfo "370 12 4.25"
puts $ppNodesInfo "371 12.5 4.25"
puts $ppNodesInfo "372 13 4.25"
puts $ppNodesInfo "373 13.5 4.25"
puts $ppNodesInfo "374 14 4.25"
puts $ppNodesInfo "375 14.5 4.25"
puts $ppNodesInfo "376 15 4.25"
puts $ppNodesInfo "377 15.5 4.25"
puts $ppNodesInfo "378 16 4.25"
puts $ppNodesInfo "379 4.586 4.5"
puts $ppNodesInfo "380 5 4.5"
puts $ppNodesInfo "381 5.6805 4.5"
puts $ppNodesInfo "382 6.9195 4.5"
puts $ppNodesInfo "383 7.5 4.5"
puts $ppNodesInfo "384 8 4.5"
puts $ppNodesInfo "385 8.5 4.5"
puts $ppNodesInfo "386 9 4.5"
puts $ppNodesInfo "387 9.5 4.5"
puts $ppNodesInfo "388 10 4.5"
puts $ppNodesInfo "389 10.5 4.5"
puts $ppNodesInfo "390 11 4.5"
puts $ppNodesInfo "391 11.5 4.5"
puts $ppNodesInfo "392 12 4.5"
puts $ppNodesInfo "393 12.5 4.5"
puts $ppNodesInfo "394 13 4.5"
puts $ppNodesInfo "395 13.5 4.5"
puts $ppNodesInfo "396 14 4.5"
puts $ppNodesInfo "397 14.5 4.5"
puts $ppNodesInfo "398 15 4.5"
puts $ppNodesInfo "399 15.5 4.5"
puts $ppNodesInfo "400 16 4.5"

close $ppNodesInfo
puts "Finished creating all pore pressure nodes..."

# define fixities for pore pressure nodes at base
fix 1 1 1 0
fix 2 1 1 0
fix 3 1 1 0
fix 4 1 1 0
fix 5 1 1 0
fix 6 1 1 0
fix 7 1 1 0
fix 8 1 1 0
fix 9 1 1 0
fix 10 1 1 0
fix 11 1 1 0
fix 12 1 1 0
fix 13 1 1 0
fix 14 1 1 0
fix 15 1 1 0
fix 16 1 1 0
fix 17 1 1 0
fix 18 1 1 0
fix 19 1 1 0
fix 20 1 1 0
fix 21 1 1 0
fix 22 1 1 0
fix 23 1 1 0
fix 24 1 1 0
fix 25 1 1 0
fix 26 1 1 0
fix 27 1 1 0
fix 28 1 1 0
fix 29 1 1 0
fix 30 1 1 0
fix 31 1 1 0
fix 32 1 1 0
fix 33 1 1 0
fix 34 1 1 0

# define fixities for nodes above water table:
fix 331 0 0 1
fix 332 0 0 1
fix 333 0 0 1
fix 334 0 0 1
fix 345 0 0 1
fix 346 0 0 1
fix 347 0 0 1
fix 348 0 0 1
fix 349 0 0 1
fix 350 0 0 1
fix 351 0 0 1
fix 352 0 0 1
fix 353 0 0 1
fix 354 0 0 1
fix 355 0 0 1
fix 356 0 0 1
fix 357 0 0 1
fix 358 0 0 1
fix 359 0 0 1
fix 360 0 0 1
fix 361 0 0 1
fix 362 0 0 1
fix 363 0 0 1
fix 364 0 0 1
fix 365 0 0 1
fix 366 0 0 1
fix 367 0 0 1
fix 368 0 0 1
fix 369 0 0 1
fix 370 0 0 1
fix 371 0 0 1
fix 372 0 0 1
fix 373 0 0 1
fix 374 0 0 1
fix 375 0 0 1
fix 376 0 0 1
fix 377 0 0 1
fix 378 0 0 1
fix 379 0 0 1
fix 380 0 0 1
fix 381 0 0 1
fix 382 0 0 1
fix 383 0 0 1
fix 384 0 0 1
fix 385 0 0 1
fix 386 0 0 1
fix 387 0 0 1
fix 388 0 0 1
fix 389 0 0 1
fix 390 0 0 1
fix 391 0 0 1
fix 392 0 0 1
fix 393 0 0 1
fix 394 0 0 1
fix 395 0 0 1
fix 396 0 0 1
fix 397 0 0 1
fix 398 0 0 1
fix 399 0 0 1
fix 400 0 0 1

# define equal degrees of freedom for pore pressure nodes
equalDOF 1 34 1 2
equalDOF 35 68 1 2
equalDOF 69 102 1 2
equalDOF 103 136 1 2
equalDOF 137 170 1 2
equalDOF 171 204 1 2
equalDOF 205 238 1 2
equalDOF 239 272 1 2
equalDOF 273 306 1 2

puts "Finished creating equalDOF for pore pressure nodes..."

##-----------------------------------------------------------------------------------------
## 3. CREATE LAGRANGE MULTIPLIER NODES FOR BEAM CONTACT ELEMENTS
##-----------------------------------------------------------------------------------------
#
for {set k 1} {$k <= 9} {incr k 1} {
node [expr 998+$k] 0.00 0.00
}
puts "Finished creating all -ndf 2 nodes..."

##-----------------------------------------------------------------------------------------
### 4. CREATE SOIL MATERIALS
##-----------------------------------------------------------------------------------------

## define grade of slope (%)
set grade 0.0
set slope [expr atan($grade/100.0)]
set g -9.81

nDMaterial PressureDependMultiYield02 3 2 2 1.1e5 2.4e5 35 0.1 \
101.0 0.5 26 0.028 0.05 0.1 \
0.05 20 5.0 3.0 1.0 \
0.0 0.65 0.9 0.02 0.7 101.0

set thick(3) 4
set xWgt(3) 0.0
set yWgt(3) [expr $g*cos($slope)]
set eVoid(3) 0.558
set porosity(3) [expr $eVoid(3)/(1+$eVoid(3))]
set uBulk(3) 2.2e6
set hPerm(3) 1.0e-2
set vPerm(3) 1.0e-2

nDMaterial PressureDependMultiYield02 2 2 2.393 1.1e5 2.4e5 35 0.1 \
101.0 0.5 26 0.028 0.05 0.1 \
0.05 20 5.0 3.0 1.0 \
0.0 0.65 0.9 0.02 0.7 101.0
set thick(2) 4
set xWgt(2) 0.0
set yWgt(2) [expr $g*cos($slope)]
set eVoid(2) 0.558
set porosity(2) [expr $eVoid(2)/(1+$eVoid(2))]
set uBulk(2) 2.2e6
set hPerm(2) 1.0e-5
set vPerm(2) 1.0e-5

nDMaterial PressureDependMultiYield02 1 2 2.393 1.1e5 2.4e5 35 0.1 \
101.0 0.5 26 0.028 0.05 0.1 \
0.05 20 5.0 3.0 1.0 \
0.0 0.65 0.9 0.02 0.7 101.0

set thick(1) 4
set xWgt(1) 0.0
set yWgt(1) [expr $g*cos($slope)]
set eVoid(1) 0.558
set porosity(1) [expr $eVoid(1)/(1+$eVoid(1))]
set uBulk(1) 2.2e6
set hPerm(1) 1.0e-5
set vPerm(1) 1.0e-5

## create wrapper material for initial state analysis
nDMaterial InitialStateAnalysisWrapper 8 1 2
nDMaterial InitialStateAnalysisWrapper 9 2 2
nDMaterial InitialStateAnalysisWrapper 10 3 2

puts "Finished creating all soil materials..."

##-----------------------------------------------------------------------------------------
## 5. CREATE SOIL ELEMENTS
##-----------------------------------------------------------------------------------------
set ElemInfo [open ElemInfo.dat w]

# elements at elevation 0m~3.2m, on the left of the pile and under GWT
for {set j 1} {$j <= 8} {incr j 1} {
for {set i 1} {$i <= 9} {incr i 1} {
set nI [expr $i+($j-1)*34]
set nJ [expr $nI + 1]
set nK [expr $nI + 35]
set nL [expr $nI + 34]

element quadUP [expr $i+($j-1)*9] $nI $nJ $nK $nL \
$thick(1) 8 $uBulk(1) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(1) $yWgt(1)
puts $ElemInfo "[expr $i+($j-1)*9] $nI $nJ $nK $nL"

}
}

# elements at elevation 0m~3.2m, on the right of the pile
for {set j 1} {$j <= 8} {incr j 1} {
for {set i 11} {$i <= 33} {incr i 1} {
set nI [expr $i+($j-1)*34]
set nJ [expr $nI + 1]
set nK [expr $nI + 35]
set nL [expr $nI + 34]

element quadUP [expr $i-10+($j-1)*23+72] $nI $nJ $nK $nL \
$thick(1) 8 $uBulk(1) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(1) $yWgt(1)
puts $ElemInfo "[expr $i-10+($j-1)*23+72] $nI $nJ $nK $nL"

}
}

# elements at elevation 3.2m~4.0m,on the right of the pile
for {set j 1} {$j <= 2} {incr j 1} {
for {set i 1} {$i <= 23} {incr i 1} {
set nI [expr 282+$i+($j-1)*24]
set nJ [expr $nI + 1]
set nK [expr $nI + 25]
set nL [expr $nI + 24]

element quadUP [expr 256+$i+($j-1)*23] $nI $nJ $nK $nL \
$thick(2) 9 $uBulk(2) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(2) $yWgt(2)
puts $ElemInfo "[expr 256+$i+($j-1)*23] $nI $nJ $nK $nL"
}
}

# elements at elevation 4.0m~4.25m, on the right of the sheet pile
for {set i 1} {$i <= 23} {incr i 1} {
set nI [expr $i+330]
set nJ [expr $nI + 1]
set nK [expr $nI + 25]
set nL [expr $nI + 24]

element quadUP [expr $i+302] $nI $nJ $nK $nL \
$thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
puts $ElemInfo "[expr $i+302] $nI $nJ $nK $nL"
}

# elements at elevation 4.25m~4.5m, on the right of the sheet pile
element quadUP 326 355 356 380 379 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 327 356 357 381 380 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 328 360 361 383 382 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 329 361 362 384 383 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 330 362 363 385 384 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 331 363 364 386 385 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 332 364 365 387 386 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 333 365 366 388 387 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 334 366 367 389 388 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 335 367 368 390 389 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 336 368 369 391 390 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 337 369 370 392 391 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 338 370 371 393 392 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 339 371 372 394 393 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 340 372 373 395 394 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 341 373 374 396 395 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 342 374 375 397 396 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 343 375 376 398 397 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 344 376 377 399 398 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)
element quadUP 345 377 378 400 399 $thick(3) 10 $uBulk(3) 1.0 [expr 0.0003/9.81] [expr 0.0003/9.81] $xWgt(3) $yWgt(3)

puts $ElemInfo "326 355 356 380 379"
puts $ElemInfo "327 356 357 381 380"
puts $ElemInfo "328 360 361 383 382"
puts $ElemInfo "329 361 362 384 383"
puts $ElemInfo "330 362 363 385 384"
puts $ElemInfo "331 363 364 386 385"
puts $ElemInfo "332 364 365 387 386"
puts $ElemInfo "333 365 366 388 387"
puts $ElemInfo "334 366 367 389 388"
puts $ElemInfo "335 367 368 390 389"
puts $ElemInfo "336 368 369 391 390"
puts $ElemInfo "337 369 370 392 391"
puts $ElemInfo "338 370 371 393 392"
puts $ElemInfo "339 371 372 394 393"
puts $ElemInfo "340 372 373 395 394"
puts $ElemInfo "341 373 374 396 395"
puts $ElemInfo "342 374 375 397 396"
puts $ElemInfo "343 375 376 398 397"
puts $ElemInfo "344 376 377 399 398"
puts $ElemInfo "345 377 378 400 399"

close $ElemInfo
puts "Finished creating all soil elements..."

##-----------------------------------------------------------------------------------------
## 6. DEFINE NODAL MASSES FOR MODELING WATER
##-----------------------------------------------------------------------------------------
#
## define nodal masses for 3 dof nodes, mass
model BasicBuilder -ndm 2 -ndf 3
set m [expr 4.5*0.8*4/10]
mass 273 0.000 $m 0.000
mass 274 0.000 $m 0.000
mass 275 0.000 $m 0.000
mass 276 0.000 $m 0.000
mass 277 0.000 $m 0.000
mass 278 0.000 $m 0.000
mass 279 0.000 $m 0.000
mass 280 0.000 $m 0.000
mass 281 0.000 $m 0.000
mass 282 0.000 $m 0.000
puts "Finished creating -ndf 3 nodal masses..."

#-----------------------------------------------------------------------------------------
# 7. CREATE sheet pile BEAM NODES AND FIXITIES
#-----------------------------------------------------------------------------------------
model BasicBuilder -ndm 2 -ndf 3

# define beam nodes
node 8998 4.543 -0.2
node 8999 4.543 0.2
node 9000 4.543 0.6
node 9001 4.543 1.0
node 9002 4.543 1.4
node 9003 4.543 1.8
node 9004 4.543 2.2
node 9005 4.543 2.6
node 9006 4.543 3.0
node 9007 4.543 3.4
node 9008 4.543 3.8
node 9009 4.543 4.2
node 9010 4.543 4.6

## fix the base node of the sheetpile in the vertial direction
fix 8998 0 1 0
puts "Finished creating all -ndf 3 beam nodes..."

#-----------------------------------------------------------------------------------------
# 8. CREATE BEAM MATERIALS
#-----------------------------------------------------------------------------------------

# beam properties
set bthick 0.086
set area 7.54e-3
set I 1.07e-6
set beamE 2.06e8
set numIntPts 3
set transTag 1
set secTag 1

# geometric transformation
geomTransf Linear $transTag

# beam section
section Elastic $secTag $beamE $area $I

puts "Finished creating all beam materials..."

#-----------------------------------------------------------------------------------------
# 9. CREATE BEAM ELEMENTS
#-----------------------------------------------------------------------------------------
element dispBeamColumn 8998 8998 8999 $numIntPts $secTag $transTag
element dispBeamColumn 8999 8999 9000 $numIntPts $secTag $transTag
element dispBeamColumn 9000 9000 9001 $numIntPts $secTag $transTag
element dispBeamColumn 9001 9001 9002 $numIntPts $secTag $transTag
element dispBeamColumn 9002 9002 9003 $numIntPts $secTag $transTag
element dispBeamColumn 9003 9003 9004 $numIntPts $secTag $transTag
element dispBeamColumn 9004 9004 9005 $numIntPts $secTag $transTag
element dispBeamColumn 9005 9005 9006 $numIntPts $secTag $transTag
element dispBeamColumn 9006 9006 9007 $numIntPts $secTag $transTag
element dispBeamColumn 9007 9007 9008 $numIntPts $secTag $transTag
element dispBeamColumn 9008 9008 9009 $numIntPts $secTag $transTag
element dispBeamColumn 9009 9009 9010 $numIntPts $secTag $transTag

puts "Finished creating all beam elements..."

##-----------------------------------------------------------------------------------------
## 10. CREATE CONTACT MATERIAL FOR BEAM CONTACT ELEMENTS
##-----------------------------------------------------------------------------------------
#
## two-dimensional contact material
nDMaterial ContactMaterial2D 4 0.1 1000.0 0.0 0.0

puts "Finished creating all contact materials..."

##-----------------------------------------------------------------------------------------
## 11. CREATE BEAM CONTACT ELEMENTS
##-----------------------------------------------------------------------------------------
#
# set gap and force tolerances for beam contact elements
set gapTol 1.0e-10
set forceTol 1.0e-10

## define beam contact elements
element BeamContact2D 998 8998 8999 10 999 4 $bthick $gapTol $forceTol
element BeamContact2D 999 8999 9000 44 1000 4 $bthick $gapTol $forceTol
element BeamContact2D 1000 9000 9001 78 1001 4 $bthick $gapTol $forceTol
element BeamContact2D 1001 9001 9002 112 1002 4 $bthick $gapTol $forceTol
element BeamContact2D 1002 9002 9003 146 1003 4 $bthick $gapTol $forceTol
element BeamContact2D 1003 9003 9004 180 1004 4 $bthick $gapTol $forceTol
element BeamContact2D 1004 9004 9005 214 1005 4 $bthick $gapTol $forceTol
element BeamContact2D 1005 9005 9006 245 1006 4 $bthick $gapTol $forceTol
element BeamContact2D 1006 9006 9007 279 1007 4 $bthick $gapTol $forceTol

puts "Finished creating all beam-contact elements..."

##-----------------------------------------------------------------------------------------
## 15. GRAVITY ANALYSIS
##-----------------------------------------------------------------------------------------
## update materials to ensure elastic behavior
updateMaterialStage -material 8 -stage 0
updateMaterialStage -material 9 -stage 0
updateMaterialStage -material 10 -stage 0

constraints Penalty 1.e14 1.e14
test NormDispIncr 1e-5 30 1
algorithm Newton
numberer Plain
system SparseGeneral
set gamma 1.6
integrator Newmark $gamma [expr pow($gamma+0.5,2)/4]
analysis Transient
# turn on initial state analysis feature
InitialStateAnalysis on

## set contact elements to be frictionless for gravity analysis
setParameter -value 0 -eleRange 998 1006 friction
analyze 10 5.0e2
puts "Finished with elastic gravity analysis..."

## update materials to consider plastic behavior
updateMaterialStage -material 8 -stage 1
updateMaterialStage -material 9 -stage 1
updateMaterialStage -material 10 -stage 1

## plastic gravity loading
analyze 40 5.0e-2
puts "Finished with plastic gravity analysis..."

# designate end of initial state analysis (zeros displacements, keeps state variables)
InitialStateAnalysis off

## turn on frictional behavior for beam contact elements
setParameter -value 1 -eleRange 998 1006 friction

[TJC]

Check your BeamContact2D inputs. Your slave and Lagrange multiplier nodes should be -ndf 2. From the looks of it, they are -ndf 3 nodes.

http://opensees.berkeley.edu/wiki/index ... mContact2D

[sanjianke]

Hi, TJC,

Thanks for your reply. I think you are right. But the problem is I cannot set the dof of slave nodes (i.e. soil nodes) in BeamContact2D elements to be 2 because I would like to capture the pore water pressure change and liquefaction characteristics of soil during shaking. That's why I use quadUP element to model the soil. Do you have any idea to solve this problem? Thank you so much.


TJC wrote:
> Check your BeamContact2D inputs. Your slave and Lagrange multiplier nodes
> should be -ndf 2. From the looks of it, they are -ndf 3 nodes.
>
> http://opensees.berkeley.edu/wiki/index ... mContact2D

[yeyifan]

Hello, I also came across this problem recently, how did you solve it.