Convergense problem for kinematic pile-soil interaction

[Roy Michel]

Hi

I'm working with kinematic pile-soil interaction problem on my Master Thesis. I have defined a soil column which is 10 by 10 meters (2D, 2DOF) which is connected to a single pile (2D, 3DOF). The contact formulation is defined by zerolength elements with PySimple1-material.

I want to model the soil is modelled with quad elements of 10 meter width and 0.1m high, but the analyze does not converge. (100 quad elements and 100 elastic beam column elements)

When I run the input file, it runs well for quad elements of 10 meter width and 0.25m high. (40 quad elements and 40 elastic beam column elements)

What is the problem in my model?

I would be grateful if someone could help me!

-Roy


My input-file is:
# Created by Roy Nalbant
#
# One elastic soil column subjected to own weigth (static) and sinusoidal load at bedrock
# Soil column has dimensions 10m x 10m (w x H) and elements are 10m x 0.1m (w x heigth of elements)
# Soil column is connected to a pile with "p-y springs"
#
wipe
#
# Use kN, m, s
#
# Some user defined variables
#
set Es 13000 ;# Young's modulus for the soil
set poisson 0.3 ;#
set rho 2 ;# mass density of the elastic medium

set thick 1 ;# thickness of the elements
set pressure 0 ;# isotropic consolidation pressure on quad elements
set gravX 0 ;# unit weigth in x-direction
set gravY [expr -9.81*$rho] ;# unit weigth in y-direction

set A 1 ;# Cross-section area of the pile
set Ep 9000000 ;# Young's modulus for the pile
set Iz 0.083333333 ;# second moment of area about the local z-axis
set masspile 2500 ;# mass of the pile in kg/m, is later divided by 2000 to get kN pr 0.5m

set pult 225 ;# 9*su, given per meter
set y50 0.05 ;# disp when 50% of pult is mobilized

set omega 10 ;# circular frequency of sine wave
set pi 3.141592654
set tStart 0 ;# start time when the load factor is applied
set tFinish 20 ;# end time when the load factor is applied
set period [expr (2*$pi)/$omega] ;# characteristic period of sine wave
set cFactor 1 ;# load-factor coefficient

set gamma 0.6 ;# Newmark gamma parameter
set beta 0.25 ;# Newmark beta paramter
set alphaM 0 ;# mass-proportional-damping (R-damping)
set betaK 0 ;# currentstate-stiffness-proportional-damping (R-damping)
set betaKinit 0.001 ;# initialstate-stiffness-proportional-damping (R-damping)
set betaKcomm 0 ;# lastcommitstate-stiffness-proportional-damping (R-damping)

set numSteps 20000 ;# number of steps to be analyzed
set deltaT 0.001 ;# time increment



#############################
# Create the model builder for the soil column
model basic -ndm 2 -ndf 2


###
# Define the material for the soil
# nDMaterial ElasticIsotropic $matTag $Es $v

nDMaterial ElasticIsotropic 1 $Es $poisson $rho



###
# Define the nodes
node 1 0 0
node 3 0 0.1
node 5 0 0.2
node 7 0 0.3
node 9 0 0.4
node 11 0 0.5
node 13 0 0.6
node 15 0 0.7
node 17 0 0.8
node 19 0 0.9
node 21 0 1
node 23 0 1.1
node 25 0 1.2
node 27 0 1.3
node 29 0 1.4
node 31 0 1.5
node 33 0 1.6
node 35 0 1.7
node 37 0 1.8
node 39 0 1.9
node 41 0 2
node 43 0 2.1
node 45 0 2.2
node 47 0 2.3
node 49 0 2.4
node 51 0 2.5
node 53 0 2.6
node 55 0 2.7
node 57 0 2.8
node 59 0 2.9
node 61 0 3
node 63 0 3.1
node 65 0 3.2
node 67 0 3.3
node 69 0 3.4
node 71 0 3.5
node 73 0 3.6
node 75 0 3.7
node 77 0 3.8
node 79 0 3.9
node 81 0 4
node 83 0 4.1
node 85 0 4.2
node 87 0 4.3
node 89 0 4.4
node 91 0 4.5
node 93 0 4.6
node 95 0 4.7
node 97 0 4.8
node 99 0 4.9
node 101 0 5
node 103 0 5.1
node 105 0 5.2
node 107 0 5.3
node 109 0 5.4
node 111 0 5.5
node 113 0 5.6
node 115 0 5.7
node 117 0 5.8
node 119 0 5.9
node 121 0 6
node 123 0 6.1
node 125 0 6.2
node 127 0 6.3
node 129 0 6.4
node 131 0 6.5
node 133 0 6.6
node 135 0 6.7
node 137 0 6.8
node 139 0 6.9
node 141 0 7
node 143 0 7.1
node 145 0 7.2
node 147 0 7.3
node 149 0 7.4
node 151 0 7.5
node 153 0 7.6
node 155 0 7.7
node 157 0 7.8
node 159 0 7.9
node 161 0 8
node 163 0 8.1
node 165 0 8.2
node 167 0 8.3
node 169 0 8.4
node 171 0 8.5
node 173 0 8.6
node 175 0 8.7
node 177 0 8.8
node 179 0 8.9
node 181 0 9
node 183 0 9.1
node 185 0 9.2
node 187 0 9.3
node 189 0 9.4
node 191 0 9.5
node 193 0 9.6
node 195 0 9.7
node 197 0 9.8
node 199 0 9.9
node 201 0 10



node 2 10 0
node 4 10 0.1
node 6 10 0.2
node 8 10 0.3
node 10 10 0.4
node 12 10 0.5
node 14 10 0.6
node 16 10 0.7
node 18 10 0.8
node 20 10 0.9
node 22 10 1
node 24 10 1.1
node 26 10 1.2
node 28 10 1.3
node 30 10 1.4
node 32 10 1.5
node 34 10 1.6
node 36 10 1.7
node 38 10 1.8
node 40 10 1.9
node 42 10 2
node 44 10 2.1
node 46 10 2.2
node 48 10 2.3
node 50 10 2.4
node 52 10 2.5
node 54 10 2.6
node 56 10 2.7
node 58 10 2.8
node 60 10 2.9
node 62 10 3
node 64 10 3.1
node 66 10 3.2
node 68 10 3.3
node 70 10 3.4
node 72 10 3.5
node 74 10 3.6
node 76 10 3.7
node 78 10 3.8
node 80 10 3.9
node 82 10 4
node 84 10 4.1
node 86 10 4.2
node 88 10 4.3
node 90 10 4.4
node 92 10 4.5
node 94 10 4.6
node 96 10 4.7
node 98 10 4.8
node 100 10 4.9
node 102 10 5
node 104 10 5.1
node 106 10 5.2
node 108 10 5.3
node 110 10 5.4
node 112 10 5.5
node 114 10 5.6
node 116 10 5.7
node 118 10 5.8
node 120 10 5.9
node 122 10 6
node 124 10 6.1
node 126 10 6.2
node 128 10 6.3
node 130 10 6.4
node 132 10 6.5
node 134 10 6.6
node 136 10 6.7
node 138 10 6.8
node 140 10 6.9
node 142 10 7
node 144 10 7.1
node 146 10 7.2
node 148 10 7.3
node 150 10 7.4
node 152 10 7.5
node 154 10 7.6
node 156 10 7.7
node 158 10 7.8
node 160 10 7.9
node 162 10 8
node 164 10 8.1
node 166 10 8.2
node 168 10 8.3
node 170 10 8.4
node 172 10 8.5
node 174 10 8.6
node 176 10 8.7
node 178 10 8.8
node 180 10 8.9
node 182 10 9
node 184 10 9.1
node 186 10 9.2
node 188 10 9.3
node 190 10 9.4
node 192 10 9.5
node 194 10 9.6
node 196 10 9.7
node 198 10 9.8
node 200 10 9.9
node 202 10 10



###
# Define the elements
#
#element quad $eleTag $iNode $jNode $kNode $lNode $thick $type $matTag $pressure mDensity $gravX $gravY
#
element quad 1 1 2 4 3 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 2 3 4 6 5 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 3 5 6 8 7 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 4 7 8 10 9 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 5 9 10 12 11 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 6 11 12 14 13 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 7 13 14 16 15 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 8 15 16 18 17 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 9 17 18 20 19 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 10 19 20 22 21 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 11 21 22 24 23 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 12 23 24 26 25 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 13 25 26 28 27 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 14 27 28 30 29 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 15 29 30 32 31 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 16 31 32 34 33 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 17 33 34 36 35 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 18 35 36 38 37 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 19 37 38 40 39 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 20 39 40 42 41 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 21 41 42 44 43 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 22 43 44 46 45 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 23 45 46 48 47 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 24 47 48 50 49 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 25 49 50 52 51 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 26 51 52 54 53 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 27 53 54 56 55 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 28 55 56 58 57 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 29 57 58 60 59 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 30 59 60 62 61 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 31 61 62 64 63 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 32 63 64 66 65 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 33 65 66 68 67 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 34 67 68 70 69 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 35 69 70 72 71 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 36 71 72 74 73 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 37 73 74 76 75 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 38 75 76 78 77 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 39 77 78 80 79 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 40 79 80 82 81 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 41 81 82 84 83 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 42 83 84 86 85 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 43 85 86 88 87 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 44 87 88 90 89 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 45 89 90 92 91 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 46 91 92 94 93 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 47 93 94 96 95 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 48 95 96 98 97 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 49 97 98 100 99 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 50 99 100 102 101 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 51 101 102 104 103 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 52 103 104 106 105 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 53 105 106 108 107 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 54 107 108 110 109 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 55 109 110 112 111 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 56 111 112 114 113 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 57 113 114 116 115 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 58 115 116 118 117 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 59 117 118 120 119 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 60 119 120 122 121 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 61 121 122 124 123 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 62 123 124 126 125 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 63 125 126 128 127 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 64 127 128 130 129 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 65 129 130 132 131 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 66 131 132 134 133 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 67 133 134 136 135 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 68 135 136 138 137 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 69 137 138 140 139 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 70 139 140 142 141 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 71 141 142 144 143 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 72 143 144 146 145 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 73 145 146 148 147 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 74 147 148 150 149 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 75 149 150 152 151 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 76 151 152 154 153 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 77 153 154 156 155 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 78 155 156 158 157 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 79 157 158 160 159 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 80 159 160 162 161 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 81 161 162 164 163 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 82 163 164 166 165 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 83 165 166 168 167 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 84 167 168 170 169 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 85 169 170 172 171 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 86 171 172 174 173 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 87 173 174 176 175 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 88 175 176 178 177 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 89 177 178 180 179 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 90 179 180 182 181 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 91 181 182 184 183 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 92 183 184 186 185 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 93 185 186 188 187 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 94 187 188 190 189 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 95 189 190 192 191 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 96 191 192 194 193 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 97 193 194 196 195 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 98 195 196 198 197 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 99 197 198 200 199 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY
element quad 100 199 200 202 201 $thick "PlaneStrain" 1 $pressure 0 $gravX $gravY



###
# dummy-nodes
node 301 10 0
node 302 10 0.1
node 303 10 0.2
node 304 10 0.3
node 305 10 0.4
node 306 10 0.5
node 307 10 0.6
node 308 10 0.7
node 309 10 0.8
node 310 10 0.9
node 311 10 1
node 312 10 1.1
node 313 10 1.2
node 314 10 1.3
node 315 10 1.4
node 316 10 1.5
node 317 10 1.6
node 318 10 1.7
node 319 10 1.8
node 320 10 1.9
node 321 10 2
node 322 10 2.1
node 323 10 2.2
node 324 10 2.3
node 325 10 2.4
node 326 10 2.5
node 327 10 2.6
node 328 10 2.7
node 329 10 2.8
node 330 10 2.9
node 331 10 3
node 332 10 3.1
node 333 10 3.2
node 334 10 3.3
node 335 10 3.4
node 336 10 3.5
node 337 10 3.6
node 338 10 3.7
node 339 10 3.8
node 340 10 3.9
node 341 10 4
node 342 10 4.1
node 343 10 4.2
node 344 10 4.3
node 345 10 4.4
node 346 10 4.5
node 347 10 4.6
node 348 10 4.7
node 349 10 4.8
node 350 10 4.9
node 351 10 5
node 352 10 5.1
node 353 10 5.2
node 354 10 5.3
node 355 10 5.4
node 356 10 5.5
node 357 10 5.6
node 358 10 5.7
node 359 10 5.8
node 360 10 5.9
node 361 10 6
node 362 10 6.1
node 363 10 6.2
node 364 10 6.3
node 365 10 6.4
node 366 10 6.5
node 367 10 6.6
node 368 10 6.7
node 369 10 6.8
node 370 10 6.9
node 371 10 7
node 372 10 7.1
node 373 10 7.2
node 374 10 7.3
node 375 10 7.4
node 376 10 7.5
node 377 10 7.6
node 378 10 7.7
node 379 10 7.8
node 380 10 7.9
node 381 10 8
node 382 10 8.1
node 383 10 8.2
node 384 10 8.3
node 385 10 8.4
node 386 10 8.5
node 387 10 8.6
node 388 10 8.7
node 389 10 8.8
node 390 10 8.9
node 391 10 9
node 392 10 9.1
node 393 10 9.2
node 394 10 9.3
node 395 10 9.4
node 396 10 9.5
node 397 10 9.6
node 398 10 9.7
node 399 10 9.8
node 400 10 9.9
node 401 10 10




###
# fixing the dummy nodes in the y-direction
fix 301 0 1
fix 302 0 1
fix 303 0 1
fix 304 0 1
fix 305 0 1
fix 306 0 1
fix 307 0 1
fix 308 0 1
fix 309 0 1
fix 310 0 1
fix 311 0 1
fix 312 0 1
fix 313 0 1
fix 314 0 1
fix 315 0 1
fix 316 0 1
fix 317 0 1
fix 318 0 1
fix 319 0 1
fix 320 0 1
fix 321 0 1
fix 322 0 1
fix 323 0 1
fix 324 0 1
fix 325 0 1
fix 326 0 1
fix 327 0 1
fix 328 0 1
fix 329 0 1
fix 330 0 1
fix 331 0 1
fix 332 0 1
fix 333 0 1
fix 334 0 1
fix 335 0 1
fix 336 0 1
fix 337 0 1
fix 338 0 1
fix 339 0 1
fix 340 0 1
fix 341 0 1
fix 342 0 1
fix 343 0 1
fix 344 0 1
fix 345 0 1
fix 346 0 1
fix 347 0 1
fix 348 0 1
fix 349 0 1
fix 350 0 1
fix 351 0 1
fix 352 0 1
fix 353 0 1
fix 354 0 1
fix 355 0 1
fix 356 0 1
fix 357 0 1
fix 358 0 1
fix 359 0 1
fix 360 0 1
fix 361 0 1
fix 362 0 1
fix 363 0 1
fix 364 0 1
fix 365 0 1
fix 366 0 1
fix 367 0 1
fix 368 0 1
fix 369 0 1
fix 370 0 1
fix 371 0 1
fix 372 0 1
fix 373 0 1
fix 374 0 1
fix 375 0 1
fix 376 0 1
fix 377 0 1
fix 378 0 1
fix 379 0 1
fix 380 0 1
fix 381 0 1
fix 382 0 1
fix 383 0 1
fix 384 0 1
fix 385 0 1
fix 386 0 1
fix 387 0 1
fix 388 0 1
fix 389 0 1
fix 390 0 1
fix 391 0 1
fix 392 0 1
fix 393 0 1
fix 394 0 1
fix 395 0 1
fix 396 0 1
fix 397 0 1
fix 398 0 1
fix 399 0 1
fix 400 0 1
fix 401 0 1






###
# Define the uniaxial material for the p-y springs
# uniaxialMaterial PySimple1 $matTag $soilType $pult $Y50 $Cd

uniaxialMaterial PySimple1 2 2 [expr $pult/8] $y50 1
uniaxialMaterial PySimple1 3 2 [expr $pult/4] $y50 1


###
# Defining the zerolength elements (springs)
element zeroLength 201 2 301 -mat 2 -dir 1
element zeroLength 202 4 302 -mat 3 -dir 1
element zeroLength 203 6 303 -mat 3 -dir 1
element zeroLength 204 8 304 -mat 3 -dir 1
element zeroLength 205 10 305 -mat 3 -dir 1
element zeroLength 206 12 306 -mat 3 -dir 1
element zeroLength 207 14 307 -mat 3 -dir 1
element zeroLength 208 16 308 -mat 3 -dir 1
element zeroLength 209 18 309 -mat 3 -dir 1
element zeroLength 210 20 310 -mat 3 -dir 1
element zeroLength 211 22 311 -mat 3 -dir 1
element zeroLength 212 24 312 -mat 3 -dir 1
element zeroLength 213 26 313 -mat 3 -dir 1
element zeroLength 214 28 314 -mat 3 -dir 1
element zeroLength 215 30 315 -mat 3 -dir 1
element zeroLength 216 32 316 -mat 3 -dir 1
element zeroLength 217 34 317 -mat 3 -dir 1
element zeroLength 218 36 318 -mat 3 -dir 1
element zeroLength 219 38 319 -mat 3 -dir 1
element zeroLength 220 40 320 -mat 3 -dir 1
element zeroLength 221 42 321 -mat 3 -dir 1
element zeroLength 222 44 322 -mat 3 -dir 1
element zeroLength 223 46 323 -mat 3 -dir 1
element zeroLength 224 48 324 -mat 3 -dir 1
element zeroLength 225 50 325 -mat 3 -dir 1
element zeroLength 226 52 326 -mat 3 -dir 1
element zeroLength 227 54 327 -mat 3 -dir 1
element zeroLength 228 56 328 -mat 3 -dir 1
element zeroLength 229 58 329 -mat 3 -dir 1
element zeroLength 230 60 330 -mat 3 -dir 1
element zeroLength 231 62 331 -mat 3 -dir 1
element zeroLength 232 64 332 -mat 3 -dir 1
element zeroLength 233 66 333 -mat 3 -dir 1
element zeroLength 234 68 334 -mat 3 -dir 1
element zeroLength 235 70 335 -mat 3 -dir 1
element zeroLength 236 72 336 -mat 3 -dir 1
element zeroLength 237 74 337 -mat 3 -dir 1
element zeroLength 238 76 338 -mat 3 -dir 1
element zeroLength 239 78 339 -mat 3 -dir 1
element zeroLength 240 80 340 -mat 3 -dir 1
element zeroLength 241 82 341 -mat 3 -dir 1
element zeroLength 242 84 342 -mat 3 -dir 1
element zeroLength 243 86 343 -mat 3 -dir 1
element zeroLength 244 88 344 -mat 3 -dir 1
element zeroLength 245 90 345 -mat 3 -dir 1
element zeroLength 246 92 346 -mat 3 -dir 1
element zeroLength 247 94 347 -mat 3 -dir 1
element zeroLength 248 96 348 -mat 3 -dir 1
element zeroLength 249 98 349 -mat 3 -dir 1
element zeroLength 250 100 350 -mat 3 -dir 1
element zeroLength 251 102 351 -mat 3 -dir 1
element zeroLength 252 104 352 -mat 3 -dir 1
element zeroLength 253 106 353 -mat 3 -dir 1
element zeroLength 254 108 354 -mat 3 -dir 1
element zeroLength 255 110 355 -mat 3 -dir 1
element zeroLength 256 112 356 -mat 3 -dir 1
element zeroLength 257 114 357 -mat 3 -dir 1
element zeroLength 258 116 358 -mat 3 -dir 1
element zeroLength 259 118 359 -mat 3 -dir 1
element zeroLength 260 120 360 -mat 3 -dir 1
element zeroLength 261 122 361 -mat 3 -dir 1
element zeroLength 262 124 362 -mat 3 -dir 1
element zeroLength 263 126 363 -mat 3 -dir 1
element zeroLength 264 128 364 -mat 3 -dir 1
element zeroLength 265 130 365 -mat 3 -dir 1
element zeroLength 266 132 366 -mat 3 -dir 1
element zeroLength 267 134 367 -mat 3 -dir 1
element zeroLength 268 136 368 -mat 3 -dir 1
element zeroLength 269 138 369 -mat 3 -dir 1
element zeroLength 270 140 370 -mat 3 -dir 1
element zeroLength 271 142 371 -mat 3 -dir 1
element zeroLength 272 144 372 -mat 3 -dir 1
element zeroLength 273 146 373 -mat 3 -dir 1
element zeroLength 274 148 374 -mat 3 -dir 1
element zeroLength 275 150 375 -mat 3 -dir 1
element zeroLength 276 152 376 -mat 3 -dir 1
element zeroLength 277 154 377 -mat 3 -dir 1
element zeroLength 278 156 378 -mat 3 -dir 1
element zeroLength 279 158 379 -mat 3 -dir 1
element zeroLength 280 160 380 -mat 3 -dir 1
element zeroLength 281 162 381 -mat 3 -dir 1
element zeroLength 282 164 382 -mat 3 -dir 1
element zeroLength 283 166 383 -mat 3 -dir 1
element zeroLength 284 168 384 -mat 3 -dir 1
element zeroLength 285 170 385 -mat 3 -dir 1
element zeroLength 286 172 386 -mat 3 -dir 1
element zeroLength 287 174 387 -mat 3 -dir 1
element zeroLength 288 176 388 -mat 3 -dir 1
element zeroLength 289 178 389 -mat 3 -dir 1
element zeroLength 290 180 390 -mat 3 -dir 1
element zeroLength 291 182 391 -mat 3 -dir 1
element zeroLength 292 184 392 -mat 3 -dir 1
element zeroLength 293 186 393 -mat 3 -dir 1
element zeroLength 294 188 394 -mat 3 -dir 1
element zeroLength 295 190 395 -mat 3 -dir 1
element zeroLength 296 192 396 -mat 3 -dir 1
element zeroLength 297 194 397 -mat 3 -dir 1
element zeroLength 298 196 398 -mat 3 -dir 1
element zeroLength 299 198 399 -mat 3 -dir 1
element zeroLength 300 200 400 -mat 3 -dir 1
element zeroLength 301 202 401 -mat 2 -dir 1







###
# tie the lateral nodes
#
#
# equalDOF $rNodeTag $cNodeTag $dof1 $dof2

equalDOF 1 2 1 2
equalDOF 3 4 1 2
equalDOF 5 6 1 2
equalDOF 7 8 1 2
equalDOF 9 10 1 2
equalDOF 11 12 1 2
equalDOF 13 14 1 2
equalDOF 15 16 1 2
equalDOF 17 18 1 2
equalDOF 19 20 1 2
equalDOF 21 22 1 2
equalDOF 23 24 1 2
equalDOF 25 26 1 2
equalDOF 27 28 1 2
equalDOF 29 30 1 2
equalDOF 31 32 1 2
equalDOF 33 34 1 2
equalDOF 35 36 1 2
equalDOF 37 38 1 2
equalDOF 39 40 1 2
equalDOF 41 42 1 2
equalDOF 43 44 1 2
equalDOF 45 46 1 2
equalDOF 47 48 1 2
equalDOF 49 50 1 2
equalDOF 51 52 1 2
equalDOF 53 54 1 2
equalDOF 55 56 1 2
equalDOF 57 58 1 2
equalDOF 59 60 1 2
equalDOF 61 62 1 2
equalDOF 63 64 1 2
equalDOF 65 66 1 2
equalDOF 67 68 1 2
equalDOF 69 70 1 2
equalDOF 71 72 1 2
equalDOF 73 74 1 2
equalDOF 75 76 1 2
equalDOF 77 78 1 2
equalDOF 79 80 1 2
equalDOF 81 82 1 2
equalDOF 83 84 1 2
equalDOF 85 86 1 2
equalDOF 87 88 1 2
equalDOF 89 90 1 2
equalDOF 91 92 1 2
equalDOF 93 94 1 2
equalDOF 95 96 1 2
equalDOF 97 98 1 2
equalDOF 99 100 1 2
equalDOF 101 102 1 2
equalDOF 103 104 1 2
equalDOF 105 106 1 2
equalDOF 107 108 1 2
equalDOF 109 110 1 2
equalDOF 111 112 1 2
equalDOF 113 114 1 2
equalDOF 115 116 1 2
equalDOF 117 118 1 2
equalDOF 119 120 1 2
equalDOF 121 122 1 2
equalDOF 123 124 1 2
equalDOF 125 126 1 2
equalDOF 127 128 1 2
equalDOF 129 130 1 2
equalDOF 131 132 1 2
equalDOF 133 134 1 2
equalDOF 135 136 1 2
equalDOF 137 138 1 2
equalDOF 139 140 1 2
equalDOF 141 142 1 2
equalDOF 143 144 1 2
equalDOF 145 146 1 2
equalDOF 147 148 1 2
equalDOF 149 150 1 2
equalDOF 151 152 1 2
equalDOF 153 154 1 2
equalDOF 155 156 1 2
equalDOF 157 158 1 2
equalDOF 159 160 1 2
equalDOF 161 162 1 2
equalDOF 163 164 1 2
equalDOF 165 166 1 2
equalDOF 167 168 1 2
equalDOF 169 170 1 2
equalDOF 171 172 1 2
equalDOF 173 174 1 2
equalDOF 175 176 1 2
equalDOF 177 178 1 2
equalDOF 179 180 1 2
equalDOF 181 182 1 2
equalDOF 183 184 1 2
equalDOF 185 186 1 2
equalDOF 187 188 1 2
equalDOF 189 190 1 2
equalDOF 191 192 1 2
equalDOF 193 194 1 2
equalDOF 195 196 1 2
equalDOF 197 198 1 2
equalDOF 199 200 1 2
equalDOF 201 202 1 2




###
# fix the base
fix 1 1 1
fix 2 1 1

###
#preventing the soil to move vertically to study the effects in the horizontal direction

fix 3 0 1
fix 4 0 1
fix 5 0 1
fix 6 0 1
fix 7 0 1
fix 8 0 1
fix 9 0 1
fix 10 0 1
fix 11 0 1
fix 12 0 1
fix 13 0 1
fix 14 0 1
fix 15 0 1
fix 16 0 1
fix 17 0 1
fix 18 0 1
fix 19 0 1
fix 20 0 1
fix 21 0 1
fix 22 0 1
fix 23 0 1
fix 24 0 1
fix 25 0 1
fix 26 0 1
fix 27 0 1
fix 28 0 1
fix 29 0 1
fix 30 0 1
fix 31 0 1
fix 32 0 1
fix 33 0 1
fix 34 0 1
fix 35 0 1
fix 36 0 1
fix 37 0 1
fix 38 0 1
fix 39 0 1
fix 40 0 1
fix 41 0 1
fix 42 0 1
fix 43 0 1
fix 44 0 1
fix 45 0 1
fix 46 0 1
fix 47 0 1
fix 48 0 1
fix 49 0 1
fix 50 0 1
fix 51 0 1
fix 52 0 1
fix 53 0 1
fix 54 0 1
fix 55 0 1
fix 56 0 1
fix 57 0 1
fix 58 0 1
fix 59 0 1
fix 60 0 1
fix 61 0 1
fix 62 0 1
fix 63 0 1
fix 64 0 1
fix 65 0 1
fix 66 0 1
fix 67 0 1
fix 68 0 1
fix 69 0 1
fix 70 0 1
fix 71 0 1
fix 72 0 1
fix 73 0 1
fix 74 0 1
fix 75 0 1
fix 76 0 1
fix 77 0 1
fix 78 0 1
fix 79 0 1
fix 80 0 1
fix 81 0 1
fix 82 0 1
fix 83 0 1
fix 84 0 1
fix 85 0 1
fix 86 0 1
fix 87 0 1
fix 88 0 1
fix 89 0 1
fix 90 0 1
fix 91 0 1
fix 92 0 1
fix 93 0 1
fix 94 0 1
fix 95 0 1
fix 96 0 1
fix 97 0 1
fix 98 0 1
fix 99 0 1
fix 100 0 1
fix 101 0 1
fix 102 0 1
fix 103 0 1
fix 104 0 1
fix 105 0 1
fix 106 0 1
fix 107 0 1
fix 108 0 1
fix 109 0 1
fix 110 0 1
fix 111 0 1
fix 112 0 1
fix 113 0 1
fix 114 0 1
fix 115 0 1
fix 116 0 1
fix 117 0 1
fix 118 0 1
fix 119 0 1
fix 120 0 1
fix 121 0 1
fix 122 0 1
fix 123 0 1
fix 124 0 1
fix 125 0 1
fix 126 0 1
fix 127 0 1
fix 128 0 1
fix 129 0 1
fix 130 0 1
fix 131 0 1
fix 132 0 1
fix 133 0 1
fix 134 0 1
fix 135 0 1
fix 136 0 1
fix 137 0 1
fix 138 0 1
fix 139 0 1
fix 140 0 1
fix 141 0 1
fix 142 0 1
fix 143 0 1
fix 144 0 1
fix 145 0 1
fix 146 0 1
fix 147 0 1
fix 148 0 1
fix 149 0 1
fix 150 0 1
fix 151 0 1
fix 152 0 1
fix 153 0 1
fix 154 0 1
fix 155 0 1
fix 156 0 1
fix 157 0 1
fix 158 0 1
fix 159 0 1
fix 160 0 1
fix 161 0 1
fix 162 0 1
fix 163 0 1
fix 164 0 1
fix 165 0 1
fix 166 0 1
fix 167 0 1
fix 168 0 1
fix 169 0 1
fix 170 0 1
fix 171 0 1
fix 172 0 1
fix 173 0 1
fix 174 0 1
fix 175 0 1
fix 176 0 1
fix 177 0 1
fix 178 0 1
fix 179 0 1
fix 180 0 1
fix 181 0 1
fix 182 0 1
fix 183 0 1
fix 184 0 1
fix 185 0 1
fix 186 0 1
fix 187 0 1
fix 188 0 1
fix 189 0 1
fix 190 0 1
fix 191 0 1
fix 192 0 1
fix 193 0 1
fix 194 0 1
fix 195 0 1
fix 196 0 1
fix 197 0 1
fix 198 0 1
fix 199 0 1
fix 200 0 1
fix 201 0 1
fix 202 0 1








#############################################################
# GRAVITY APPLICATION (elastic)
# create the SOE, ConstraintHandler, Integrator, Algorithm and Numberer
system ProfileSPD
test NormDispIncr 1.0e-5 10 0
constraints Transformation
integrator LoadControl 1 1 1 1
algorithm ModifiedNewton
numberer RCM

# create the Analysis
analysis Static

#analyze
analyze 2

puts "Gravity (soil) application complete"

#######################################################################################


##############################################################
# Create the model builder for the pile
model basic -ndm 2 -ndf 3


###
# defining the nodes for the pile
node 501 10 0
node 502 10 0.1
node 503 10 0.2
node 504 10 0.3
node 505 10 0.4
node 506 10 0.5
node 507 10 0.6
node 508 10 0.7
node 509 10 0.8
node 510 10 0.9
node 511 10 1
node 512 10 1.1
node 513 10 1.2
node 514 10 1.3
node 515 10 1.4
node 516 10 1.5
node 517 10 1.6
node 518 10 1.7
node 519 10 1.8
node 520 10 1.9
node 521 10 2
node 522 10 2.1
node 523 10 2.2
node 524 10 2.3
node 525 10 2.4
node 526 10 2.5
node 527 10 2.6
node 528 10 2.7
node 529 10 2.8
node 530 10 2.9
node 531 10 3
node 532 10 3.1
node 533 10 3.2
node 534 10 3.3
node 535 10 3.4
node 536 10 3.5
node 537 10 3.6
node 538 10 3.7
node 539 10 3.8
node 540 10 3.9
node 541 10 4
node 542 10 4.1
node 543 10 4.2
node 544 10 4.3
node 545 10 4.4
node 546 10 4.5
node 547 10 4.6
node 548 10 4.7
node 549 10 4.8
node 550 10 4.9
node 551 10 5
node 552 10 5.1
node 553 10 5.2
node 554 10 5.3
node 555 10 5.4
node 556 10 5.5
node 557 10 5.6
node 558 10 5.7
node 559 10 5.8
node 560 10 5.9
node 561 10 6
node 562 10 6.1
node 563 10 6.2
node 564 10 6.3
node 565 10 6.4
node 566 10 6.5
node 567 10 6.6
node 568 10 6.7
node 569 10 6.8
node 570 10 6.9
node 571 10 7
node 572 10 7.1
node 573 10 7.2
node 574 10 7.3
node 575 10 7.4
node 576 10 7.5
node 577 10 7.6
node 578 10 7.7
node 579 10 7.8
node 580 10 7.9
node 581 10 8
node 582 10 8.1
node 583 10 8.2
node 584 10 8.3
node 585 10 8.4
node 586 10 8.5
node 587 10 8.6
node 588 10 8.7
node 589 10 8.8
node 590 10 8.9
node 591 10 9
node 592 10 9.1
node 593 10 9.2
node 594 10 9.3
node 595 10 9.4
node 596 10 9.5
node 597 10 9.6
node 598 10 9.7
node 599 10 9.8
node 600 10 9.9
node 601 10 10


###
# Defining the fixity of the pile
fix 501 1 1 1


###
# defining the mass of each node (only applied in the x-direction)
mass 501 [expr $masspile/20000] 0 0
mass 502 [expr $masspile/10000] 0 0
mass 503 [expr $masspile/10000] 0 0
mass 504 [expr $masspile/10000] 0 0
mass 505 [expr $masspile/10000] 0 0
mass 506 [expr $masspile/10000] 0 0
mass 507 [expr $masspile/10000] 0 0
mass 508 [expr $masspile/10000] 0 0
mass 509 [expr $masspile/10000] 0 0
mass 510 [expr $masspile/10000] 0 0
mass 511 [expr $masspile/10000] 0 0
mass 512 [expr $masspile/10000] 0 0
mass 513 [expr $masspile/10000] 0 0
mass 514 [expr $masspile/10000] 0 0
mass 515 [expr $masspile/10000] 0 0
mass 516 [expr $masspile/10000] 0 0
mass 517 [expr $masspile/10000] 0 0
mass 518 [expr $masspile/10000] 0 0
mass 519 [expr $masspile/10000] 0 0
mass 520 [expr $masspile/10000] 0 0
mass 521 [expr $masspile/10000] 0 0
mass 522 [expr $masspile/10000] 0 0
mass 523 [expr $masspile/10000] 0 0
mass 524 [expr $masspile/10000] 0 0
mass 525 [expr $masspile/10000] 0 0
mass 526 [expr $masspile/10000] 0 0
mass 527 [expr $masspile/10000] 0 0
mass 528 [expr $masspile/10000] 0 0
mass 529 [expr $masspile/10000] 0 0
mass 530 [expr $masspile/10000] 0 0
mass 531 [expr $masspile/10000] 0 0
mass 532 [expr $masspile/10000] 0 0
mass 533 [expr $masspile/10000] 0 0
mass 534 [expr $masspile/10000] 0 0
mass 535 [expr $masspile/10000] 0 0
mass 536 [expr $masspile/10000] 0 0
mass 537 [expr $masspile/10000] 0 0
mass 538 [expr $masspile/10000] 0 0
mass 539 [expr $masspile/10000] 0 0
mass 540 [expr $masspile/10000] 0 0
mass 541 [expr $masspile/10000] 0 0
mass 542 [expr $masspile/10000] 0 0
mass 543 [expr $masspile/10000] 0 0
mass 544 [expr $masspile/10000] 0 0
mass 545 [expr $masspile/10000] 0 0
mass 546 [expr $masspile/10000] 0 0
mass 547 [expr $masspile/10000] 0 0
mass 548 [expr $masspile/10000] 0 0
mass 549 [expr $masspile/10000] 0 0
mass 550 [expr $masspile/10000] 0 0
mass 551 [expr $masspile/10000] 0 0
mass 552 [expr $masspile/10000] 0 0
mass 553 [expr $masspile/10000] 0 0
mass 554 [expr $masspile/10000] 0 0
mass 555 [expr $masspile/10000] 0 0
mass 556 [expr $masspile/10000] 0 0
mass 557 [expr $masspile/10000] 0 0
mass 558 [expr $masspile/10000] 0 0
mass 559 [expr $masspile/10000] 0 0
mass 560 [expr $masspile/10000] 0 0
mass 561 [expr $masspile/10000] 0 0
mass 562 [expr $masspile/10000] 0 0
mass 563 [expr $masspile/10000] 0 0
mass 564 [expr $masspile/10000] 0 0
mass 565 [expr $masspile/10000] 0 0
mass 566 [expr $masspile/10000] 0 0
mass 567 [expr $masspile/10000] 0 0
mass 568 [expr $masspile/10000] 0 0
mass 569 [expr $masspile/10000] 0 0
mass 570 [expr $masspile/10000] 0 0
mass 571 [expr $masspile/10000] 0 0
mass 572 [expr $masspile/10000] 0 0
mass 573 [expr $masspile/10000] 0 0
mass 574 [expr $masspile/10000] 0 0
mass 575 [expr $masspile/10000] 0 0
mass 576 [expr $masspile/10000] 0 0
mass 577 [expr $masspile/10000] 0 0
mass 578 [expr $masspile/10000] 0 0
mass 579 [expr $masspile/10000] 0 0
mass 580 [expr $masspile/10000] 0 0
mass 581 [expr $masspile/10000] 0 0
mass 582 [expr $masspile/10000] 0 0
mass 583 [expr $masspile/10000] 0 0
mass 584 [expr $masspile/10000] 0 0
mass 585 [expr $masspile/10000] 0 0
mass 586 [expr $masspile/10000] 0 0
mass 587 [expr $masspile/10000] 0 0
mass 588 [expr $masspile/10000] 0 0
mass 589 [expr $masspile/10000] 0 0
mass 590 [expr $masspile/10000] 0 0
mass 591 [expr $masspile/10000] 0 0
mass 592 [expr $masspile/10000] 0 0
mass 593 [expr $masspile/10000] 0 0
mass 594 [expr $masspile/10000] 0 0
mass 595 [expr $masspile/10000] 0 0
mass 596 [expr $masspile/10000] 0 0
mass 597 [expr $masspile/10000] 0 0
mass 598 [expr $masspile/10000] 0 0
mass 599 [expr $masspile/10000] 0 0
mass 600 [expr $masspile/10000] 0 0
mass 601 [expr $masspile/20000] 0 0








###
# fixing the pile nodes in the vertical direction
fix 502 0 1 0
fix 503 0 1 0
fix 504 0 1 0
fix 505 0 1 0
fix 506 0 1 0
fix 507 0 1 0
fix 508 0 1 0
fix 509 0 1 0
fix 510 0 1 0
fix 511 0 1 0
fix 512 0 1 0
fix 513 0 1 0
fix 514 0 1 0
fix 515 0 1 0
fix 516 0 1 0
fix 517 0 1 0
fix 518 0 1 0
fix 519 0 1 0
fix 520 0 1 0
fix 521 0 1 0
fix 522 0 1 0
fix 523 0 1 0
fix 524 0 1 0
fix 525 0 1 0
fix 526 0 1 0
fix 527 0 1 0
fix 528 0 1 0
fix 529 0 1 0
fix 530 0 1 0
fix 531 0 1 0
fix 532 0 1 0
fix 533 0 1 0
fix 534 0 1 0
fix 535 0 1 0
fix 536 0 1 0
fix 537 0 1 0
fix 538 0 1 0
fix 539 0 1 0
fix 540 0 1 0
fix 541 0 1 0
fix 542 0 1 0
fix 543 0 1 0
fix 544 0 1 0
fix 545 0 1 0
fix 546 0 1 0
fix 547 0 1 0
fix 548 0 1 0
fix 549 0 1 0
fix 550 0 1 0
fix 551 0 1 0
fix 552 0 1 0
fix 553 0 1 0
fix 554 0 1 0
fix 555 0 1 0
fix 556 0 1 0
fix 557 0 1 0
fix 558 0 1 0
fix 559 0 1 0
fix 560 0 1 0
fix 561 0 1 0
fix 562 0 1 0
fix 563 0 1 0
fix 564 0 1 0
fix 565 0 1 0
fix 566 0 1 0
fix 567 0 1 0
fix 568 0 1 0
fix 569 0 1 0
fix 570 0 1 0
fix 571 0 1 0
fix 572 0 1 0
fix 573 0 1 0
fix 574 0 1 0
fix 575 0 1 0
fix 576 0 1 0
fix 577 0 1 0
fix 578 0 1 0
fix 579 0 1 0
fix 580 0 1 0
fix 581 0 1 0
fix 582 0 1 0
fix 583 0 1 0
fix 584 0 1 0
fix 585 0 1 0
fix 586 0 1 0
fix 587 0 1 0
fix 588 0 1 0
fix 589 0 1 0
fix 590 0 1 0
fix 591 0 1 0
fix 592 0 1 0
fix 593 0 1 0
fix 594 0 1 0
fix 595 0 1 0
fix 596 0 1 0
fix 597 0 1 0
fix 598 0 1 0
fix 599 0 1 0
fix 600 0 1 0
fix 601 0 1 0



###
# Defining the linear coordinate transformation (LinearCrdTransf) object
# geomTransf Linear $transfTag

geomTransf Linear 1

###
# defining the elastic beam column elements
# element elasticBeamColumn $eleTag $iNode $jNode $A $E $Iz $transfTag



element elasticBeamColumn 401 501 502 $A $Ep $Iz 1
element elasticBeamColumn 402 502 503 $A $Ep $Iz 1
element elasticBeamColumn 403 503 504 $A $Ep $Iz 1
element elasticBeamColumn 404 504 505 $A $Ep $Iz 1
element elasticBeamColumn 405 505 506 $A $Ep $Iz 1
element elasticBeamColumn 406 506 507 $A $Ep $Iz 1
element elasticBeamColumn 407 507 508 $A $Ep $Iz 1
element elasticBeamColumn 408 508 509 $A $Ep $Iz 1
element elasticBeamColumn 409 509 510 $A $Ep $Iz 1
element elasticBeamColumn 410 510 511 $A $Ep $Iz 1
element elasticBeamColumn 411 511 512 $A $Ep $Iz 1
element elasticBeamColumn 412 512 513 $A $Ep $Iz 1
element elasticBeamColumn 413 513 514 $A $Ep $Iz 1
element elasticBeamColumn 414 514 515 $A $Ep $Iz 1
element elasticBeamColumn 415 515 516 $A $Ep $Iz 1
element elasticBeamColumn 416 516 517 $A $Ep $Iz 1
element elasticBeamColumn 417 517 518 $A $Ep $Iz 1
element elasticBeamColumn 418 518 519 $A $Ep $Iz 1
element elasticBeamColumn 419 519 520 $A $Ep $Iz 1
element elasticBeamColumn 420 520 521 $A $Ep $Iz 1
element elasticBeamColumn 421 521 522 $A $Ep $Iz 1
element elasticBeamColumn 422 522 523 $A $Ep $Iz 1
element elasticBeamColumn 423 523 524 $A $Ep $Iz 1
element elasticBeamColumn 424 524 525 $A $Ep $Iz 1
element elasticBeamColumn 425 525 526 $A $Ep $Iz 1
element elasticBeamColumn 426 526 527 $A $Ep $Iz 1
element elasticBeamColumn 427 527 528 $A $Ep $Iz 1
element elasticBeamColumn 428 528 529 $A $Ep $Iz 1
element elasticBeamColumn 429 529 530 $A $Ep $Iz 1
element elasticBeamColumn 430 530 531 $A $Ep $Iz 1
element elasticBeamColumn 431 531 532 $A $Ep $Iz 1
element elasticBeamColumn 432 532 533 $A $Ep $Iz 1
element elasticBeamColumn 433 533 534 $A $Ep $Iz 1
element elasticBeamColumn 434 534 535 $A $Ep $Iz 1
element elasticBeamColumn 435 535 536 $A $Ep $Iz 1
element elasticBeamColumn 436 536 537 $A $Ep $Iz 1
element elasticBeamColumn 437 537 538 $A $Ep $Iz 1
element elasticBeamColumn 438 538 539 $A $Ep $Iz 1
element elasticBeamColumn 439 539 540 $A $Ep $Iz 1
element elasticBeamColumn 440 540 541 $A $Ep $Iz 1
element elasticBeamColumn 441 541 542 $A $Ep $Iz 1
element elasticBeamColumn 442 542 543 $A $Ep $Iz 1
element elasticBeamColumn 443 543 544 $A $Ep $Iz 1
element elasticBeamColumn 444 544 545 $A $Ep $Iz 1
element elasticBeamColumn 445 545 546 $A $Ep $Iz 1
element elasticBeamColumn 446 546 547 $A $Ep $Iz 1
element elasticBeamColumn 447 547 548 $A $Ep $Iz 1
element elasticBeamColumn 448 548 549 $A $Ep $Iz 1
element elasticBeamColumn 449 549 550 $A $Ep $Iz 1
element elasticBeamColumn 450 550 551 $A $Ep $Iz 1
element elasticBeamColumn 451 551 552 $A $Ep $Iz 1
element elasticBeamColumn 452 552 553 $A $Ep $Iz 1
element elasticBeamColumn 453 553 554 $A $Ep $Iz 1
element elasticBeamColumn 454 554 555 $A $Ep $Iz 1
element elasticBeamColumn 455 555 556 $A $Ep $Iz 1
element elasticBeamColumn 456 556 557 $A $Ep $Iz 1
element elasticBeamColumn 457 557 558 $A $Ep $Iz 1
element elasticBeamColumn 458 558 559 $A $Ep $Iz 1
element elasticBeamColumn 459 559 560 $A $Ep $Iz 1
element elasticBeamColumn 460 560 561 $A $Ep $Iz 1
element elasticBeamColumn 461 561 562 $A $Ep $Iz 1
element elasticBeamColumn 462 562 563 $A $Ep $Iz 1
element elasticBeamColumn 463 563 564 $A $Ep $Iz 1
element elasticBeamColumn 464 564 565 $A $Ep $Iz 1
element elasticBeamColumn 465 565 566 $A $Ep $Iz 1
element elasticBeamColumn 466 566 567 $A $Ep $Iz 1
element elasticBeamColumn 467 567 568 $A $Ep $Iz 1
element elasticBeamColumn 468 568 569 $A $Ep $Iz 1
element elasticBeamColumn 469 569 570 $A $Ep $Iz 1
element elasticBeamColumn 470 570 571 $A $Ep $Iz 1
element elasticBeamColumn 471 571 572 $A $Ep $Iz 1
element elasticBeamColumn 472 572 573 $A $Ep $Iz 1
element elasticBeamColumn 473 573 574 $A $Ep $Iz 1
element elasticBeamColumn 474 574 575 $A $Ep $Iz 1
element elasticBeamColumn 475 575 576 $A $Ep $Iz 1
element elasticBeamColumn 476 576 577 $A $Ep $Iz 1
element elasticBeamColumn 477 577 578 $A $Ep $Iz 1
element elasticBeamColumn 478 578 579 $A $Ep $Iz 1
element elasticBeamColumn 479 579 580 $A $Ep $Iz 1
element elasticBeamColumn 480 580 581 $A $Ep $Iz 1
element elasticBeamColumn 481 581 582 $A $Ep $Iz 1
element elasticBeamColumn 482 582 583 $A $Ep $Iz 1
element elasticBeamColumn 483 583 584 $A $Ep $Iz 1
element elasticBeamColumn 484 584 585 $A $Ep $Iz 1
element elasticBeamColumn 485 585 586 $A $Ep $Iz 1
element elasticBeamColumn 486 586 587 $A $Ep $Iz 1
element elasticBeamColumn 487 587 588 $A $Ep $Iz 1
element elasticBeamColumn 488 588 589 $A $Ep $Iz 1
element elasticBeamColumn 489 589 590 $A $Ep $Iz 1
element elasticBeamColumn 490 590 591 $A $Ep $Iz 1
element elasticBeamColumn 491 591 592 $A $Ep $Iz 1
element elasticBeamColumn 492 592 593 $A $Ep $Iz 1
element elasticBeamColumn 493 593 594 $A $Ep $Iz 1
element elasticBeamColumn 494 594 595 $A $Ep $Iz 1
element elasticBeamColumn 495 595 596 $A $Ep $Iz 1
element elasticBeamColumn 496 596 597 $A $Ep $Iz 1
element elasticBeamColumn 497 597 598 $A $Ep $Iz 1
element elasticBeamColumn 498 598 599 $A $Ep $Iz 1
element elasticBeamColumn 499 599 600 $A $Ep $Iz 1
element elasticBeamColumn 500 600 601 $A $Ep $Iz 1




###
# defining that the dummy-nodes to define the p-y-springs are equal to the pile-nodes
equalDOF 301 501 1
equalDOF 302 502 1
equalDOF 303 503 1
equalDOF 304 504 1
equalDOF 305 505 1
equalDOF 306 506 1
equalDOF 307 507 1
equalDOF 308 508 1
equalDOF 309 509 1
equalDOF 310 510 1
equalDOF 311 511 1
equalDOF 312 512 1
equalDOF 313 513 1
equalDOF 314 514 1
equalDOF 315 515 1
equalDOF 316 516 1
equalDOF 317 517 1
equalDOF 318 518 1
equalDOF 319 519 1
equalDOF 320 520 1
equalDOF 321 521 1
equalDOF 322 522 1
equalDOF 323 523 1
equalDOF 324 524 1
equalDOF 325 525 1
equalDOF 326 526 1
equalDOF 327 527 1
equalDOF 328 528 1
equalDOF 329 529 1
equalDOF 330 530 1
equalDOF 331 531 1
equalDOF 332 532 1
equalDOF 333 533 1
equalDOF 334 534 1
equalDOF 335 535 1
equalDOF 336 536 1
equalDOF 337 537 1
equalDOF 338 538 1
equalDOF 339 539 1
equalDOF 340 540 1
equalDOF 341 541 1
equalDOF 342 542 1
equalDOF 343 543 1
equalDOF 344 544 1
equalDOF 345 545 1
equalDOF 346 546 1
equalDOF 347 547 1
equalDOF 348 548 1
equalDOF 349 549 1
equalDOF 350 550 1
equalDOF 351 551 1
equalDOF 352 552 1
equalDOF 353 553 1
equalDOF 354 554 1
equalDOF 355 555 1
equalDOF 356 556 1
equalDOF 357 557 1
equalDOF 358 558 1
equalDOF 359 559 1
equalDOF 360 560 1
equalDOF 361 561 1
equalDOF 362 562 1
equalDOF 363 563 1
equalDOF 364 564 1
equalDOF 365 565 1
equalDOF 366 566 1
equalDOF 367 567 1
equalDOF 368 568 1
equalDOF 369 569 1
equalDOF 370 570 1
equalDOF 371 571 1
equalDOF 372 572 1
equalDOF 373 573 1
equalDOF 374 574 1
equalDOF 375 575 1
equalDOF 376 576 1
equalDOF 377 577 1
equalDOF 378 578 1
equalDOF 379 579 1
equalDOF 380 580 1
equalDOF 381 581 1
equalDOF 382 582 1
equalDOF 383 583 1
equalDOF 384 584 1
equalDOF 385 585 1
equalDOF 386 586 1
equalDOF 387 587 1
equalDOF 388 588 1
equalDOF 389 589 1
equalDOF 390 590 1
equalDOF 391 591 1
equalDOF 392 592 1
equalDOF 393 593 1
equalDOF 394 594 1
equalDOF 395 595 1
equalDOF 396 596 1
equalDOF 397 597 1
equalDOF 398 598 1
equalDOF 399 599 1
equalDOF 400 600 1
equalDOF 401 601 1










####################################################################################

#############################################################
# NOW APPLY LOADING SEQUENCE AND ANALYZE (elastic)
# rezero time
setTime 0.0
wipeAnalysis



#pattern UniformExcitation $patternTag $dir -accel (Sine $tStart $tFinish $period <-shift $shift> <-factor $cFactor>)
pattern UniformExcitation 1 1 -accel "Sine $tStart $tFinish $period -factor $cFactor"


#create the recorders
recorder Node -file disp.out -time -node 202 -dof 1 -dT $deltaT disp


# create the Analysis
constraints Transformation
test NormDispIncr 1.0e-2 1000 0
algorithm Newton
numberer RCM
system ProfileSPD

integrator Newmark $gamma $beta $alphaM $betaK $betaKinit $betaKcomm

analysis VariableTransient



#analyze
set startT [clock seconds]
analyze $numSteps $deltaT [expr $deltaT/100] $deltaT 10
set endT [clock seconds]
puts "Dynamic analysis complete"
puts "Execution time: [expr $endT-$startT] seconds."

wipe

[silvia]

i don't know much about the soil element.
but you have a simple analyize command.
if it fails to converge, it just means that you should change something in your analysis method at this point (very typical for complicated nonlinear problems -- it is hard to find a universal analysis model).
please look at the examples manual for the pushover example of frame 5, and see what it does in case of convergence failure.

[Roy Michel]

Thank you a lot for answering!

I actually managed to find the problem. The problem was my beta-value of 0.25 in the Newmark integrator. I set the value to 0.3025 instead, and it converged perfectly!



-Roy