Beam elements with lateral soil springs

[ggokelly]

hi, im trying to make beam elements which are connect with lateral soil springs.

The model has not soil body but i added lateral soil stiffness to the beam directly with zerolength element

However If i run this pile, the window shows an error like below.

' OpenSees > source M3_J2.tcl
ProfileSPDLinDirectSolver::solve<> - aii<0<i,aii>: <57,0>
WARNING NewtonRaphson::solveCurrentStep<> -the LinearSysOfEqn failed in solve<>
StaticAnalysis::analyze<> - the Algorithm failed at iteration: 0 with domain at load factor 0.01
OpenSees> analyze failed, returned: -3 error flag'

here's my command, please help me out to fix this error.

-----------------------------------------------------------------------------------------------------------------------------------
# date: 2012. 02. 02. Thur
# Model geometry : Look at hand drawing

wipe all

######################################
# #
# #
### Build Model for Pile Structure ###
# #
# #
######################################

model BasicBuilder -ndm 2 -ndf 3
geomTransf Linear 1

# Diameter of pile
set D_beam 5
# thickness of pile
set t_beam 0.05
# Young's mod of pile in kPa
set E_beam 100000000
set I_beam 2.3817
set A_beam [expr 3.1415*$D_beam*$D_beam/4]

node 1 40.0 40.0
node 2 40.0 41.0
node 3 40.0 42.0
node 4 40.0 43.0
node 5 40.0 44.0
node 6 40.0 45.0
node 7 40.0 46.0
node 8 40.0 47.0
node 9 40.0 48.0
node 10 40.0 49.0
node 11 40.0 50.0
node 12 40.0 51.0
node 13 40.0 52.0
node 14 40.0 53.0
node 15 40.0 54.0
node 16 40.0 55.0
node 17 40.0 56.0
node 18 40.0 57.0
node 19 40.0 58.0
node 20 40.0 59.0
node 21 40.0 60.0
node 22 40.0 61.0
node 23 40.0 62.0
node 24 40.0 63.0
node 25 40.0 64.0
node 26 40.0 65.0
node 27 40.0 66.0
node 28 40.0 67.0
node 29 40.0 68.0
node 30 40.0 69.0
node 31 40.0 70.0
node 32 40.0 71.0
node 33 40.0 72.0
node 34 40.0 73.0
node 35 40.0 74.0
node 36 40.0 75.0
node 37 40.0 76.0
node 38 40.0 77.0
node 39 40.0 78.0
node 40 40.0 79.0
node 41 40.0 80.0
node 42 40.0 81.0
node 43 40.0 82.0
node 44 40.0 83.0
node 45 40.0 84.0
node 46 40.0 85.0
node 47 40.0 86.0
node 48 40.0 87.0
node 49 40.0 88.0
node 50 40.0 89.0
node 51 40.0 90.0
node 52 40.0 91.0
node 53 40.0 92.0
node 54 40.0 93.0
node 55 40.0 94.0
node 56 40.0 95.0
node 57 40.0 96.0
node 58 40.0 97.0
node 59 40.0 98.0
node 60 40.0 99.0
node 61 40.0 100.0

element elasticBeamColumn 1 1 2 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 2 2 3 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 3 3 4 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 4 4 5 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 5 5 6 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 6 6 7 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 7 7 8 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 8 8 9 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 9 9 10 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 10 10 11 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 11 11 12 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 12 12 13 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 13 13 14 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 14 14 15 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 15 15 16 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 16 16 17 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 17 17 18 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 18 18 19 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 19 19 20 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 20 20 21 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 21 21 22 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 22 22 23 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 23 23 24 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 24 24 25 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 25 25 26 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 26 26 27 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 27 27 28 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 28 28 29 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 29 29 30 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 30 30 31 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 31 31 32 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 32 32 33 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 33 33 34 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 34 34 35 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 35 35 36 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 36 36 37 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 37 37 38 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 38 38 39 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 39 39 40 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 40 40 41 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 41 41 42 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 42 42 43 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 43 43 44 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 44 44 45 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 45 45 46 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 46 46 47 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 47 47 48 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 48 48 49 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 49 49 50 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 50 50 51 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 51 51 52 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 52 52 53 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 53 53 54 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 54 54 55 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 55 55 56 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 56 56 57 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 57 57 58 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 58 58 59 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 59 59 60 $A_beam $E_beam $I_beam 1
element elasticBeamColumn 60 60 61 $A_beam $E_beam $I_beam 1

######################################################
### Add Node for zerolength Springs in Soil Domain ###
######################################################

node 62 40.0 40.0
node 63 40.0 41.0
node 64 40.0 42.0
node 65 40.0 43.0
node 66 40.0 44.0
node 67 40.0 45.0
node 68 40.0 46.0
node 69 40.0 47.0
node 70 40.0 48.0
node 71 40.0 49.0
node 72 40.0 50.0
node 73 40.0 51.0
node 74 40.0 52.0
node 75 40.0 53.0
node 76 40.0 54.0
node 77 40.0 55.0
node 78 40.0 56.0
node 79 40.0 57.0
node 80 40.0 58.0
node 81 40.0 59.0
node 82 40.0 60.0
node 83 40.0 61.0
node 84 40.0 62.0
node 85 40.0 63.0
node 86 40.0 64.0
node 87 40.0 65.0
node 88 40.0 66.0
node 89 40.0 67.0
node 90 40.0 68.0
node 91 40.0 69.0
node 92 40.0 70.0
node 93 40.0 71.0
node 94 40.0 72.0
node 95 40.0 73.0
node 96 40.0 74.0
node 97 40.0 75.0
node 98 40.0 76.0
node 99 40.0 77.0
node 100 40.0 78.0
node 101 40.0 79.0
node 102 40.0 80.0

######################################
### Add P-Y Springs in Soil Domain ###
######################################
uniaxialMaterial Elastic 2 2035.7 0
uniaxialMaterial Elastic 3 2013.1 0
uniaxialMaterial Elastic 4 1935.9 0
uniaxialMaterial Elastic 5 1975.3 0
uniaxialMaterial Elastic 6 1958.5 0
uniaxialMaterial Elastic 7 2007.9 0
uniaxialMaterial Elastic 8 1986.6 0
uniaxialMaterial Elastic 9 2120.9 0
uniaxialMaterial Elastic 10 2154.6 0
uniaxialMaterial Elastic 11 2140.1 0
uniaxialMaterial Elastic 12 2309 0
uniaxialMaterial Elastic 13 2263.3 0
uniaxialMaterial Elastic 14 2461.6 0
uniaxialMaterial Elastic 15 2690.7 0
uniaxialMaterial Elastic 16 2957.9 0
uniaxialMaterial Elastic 17 2888.9 0
uniaxialMaterial Elastic 18 3219.4 0
uniaxialMaterial Elastic 19 3621.2 0
uniaxialMaterial Elastic 20 4120 0
uniaxialMaterial Elastic 21 4755 0
uniaxialMaterial Elastic 22 5590.7 0
uniaxialMaterial Elastic 23 5760.7 0
uniaxialMaterial Elastic 24 7184.3 0
uniaxialMaterial Elastic 25 9695.2 0
uniaxialMaterial Elastic 26 13742 0
uniaxialMaterial Elastic 27 22453 0
uniaxialMaterial Elastic 28 56401 0
uniaxialMaterial Elastic 29 126640 0
uniaxialMaterial Elastic 30 30458 0
uniaxialMaterial Elastic 31 17639 0
uniaxialMaterial Elastic 32 12432 0
uniaxialMaterial Elastic 33 9717.3 0
uniaxialMaterial Elastic 34 8008.2 0
uniaxialMaterial Elastic 35 6826.3 0
uniaxialMaterial Elastic 36 5951.1 0
uniaxialMaterial Elastic 37 5270.9 0
uniaxialMaterial Elastic 38 3706.7 0
uniaxialMaterial Elastic 39 3344.5 0
uniaxialMaterial Elastic 40 3125.7 0
uniaxialMaterial Elastic 41 2856.7 0
uniaxialMaterial Elastic 42 2617.5 0
uniaxialMaterial Elastic 43 0 0

element zeroLength 61 1 62 -mat 42 43 -dir 1 2
element zeroLength 62 2 63 -mat 41 43 -dir 1 2
element zeroLength 63 3 64 -mat 40 43 -dir 1 2
element zeroLength 64 4 65 -mat 39 43 -dir 1 2
element zeroLength 65 5 66 -mat 38 43 -dir 1 2
element zeroLength 66 6 67 -mat 37 43 -dir 1 2
element zeroLength 67 7 68 -mat 36 43 -dir 1 2
element zeroLength 68 8 69 -mat 35 43 -dir 1 2
element zeroLength 69 9 70 -mat 34 43 -dir 1 2
element zeroLength 70 10 71 -mat 33 43 -dir 1 2
element zeroLength 71 11 72 -mat 32 43 -dir 1 2
element zeroLength 72 12 73 -mat 31 43 -dir 1 2
element zeroLength 73 13 74 -mat 30 43 -dir 1 2
element zeroLength 74 14 75 -mat 29 43 -dir 1 2
element zeroLength 75 15 76 -mat 28 43 -dir 1 2
element zeroLength 76 16 77 -mat 27 43 -dir 1 2
element zeroLength 77 17 78 -mat 26 43 -dir 1 2
element zeroLength 78 18 79 -mat 25 43 -dir 1 2
element zeroLength 79 19 80 -mat 24 43 -dir 1 2
element zeroLength 80 20 81 -mat 23 43 -dir 1 2
element zeroLength 81 21 82 -mat 22 43 -dir 1 2
element zeroLength 82 22 83 -mat 21 43 -dir 1 2
element zeroLength 83 23 84 -mat 20 43 -dir 1 2
element zeroLength 84 24 85 -mat 19 43 -dir 1 2
element zeroLength 85 25 86 -mat 18 43 -dir 1 2
element zeroLength 86 26 87 -mat 17 43 -dir 1 2
element zeroLength 87 27 88 -mat 16 43 -dir 1 2
element zeroLength 88 28 89 -mat 15 43 -dir 1 2
element zeroLength 89 29 90 -mat 14 43 -dir 1 2
element zeroLength 90 30 91 -mat 13 43 -dir 1 2
element zeroLength 91 31 92 -mat 12 43 -dir 1 2
element zeroLength 92 32 93 -mat 11 43 -dir 1 2
element zeroLength 93 33 94 -mat 10 43 -dir 1 2
element zeroLength 94 34 95 -mat 9 43 -dir 1 2
element zeroLength 95 35 96 -mat 8 43 -dir 1 2
element zeroLength 96 36 97 -mat 7 43 -dir 1 2
element zeroLength 97 37 98 -mat 6 43 -dir 1 2
element zeroLength 98 38 99 -mat 5 43 -dir 1 2
element zeroLength 99 39 100 -mat 4 43 -dir 1 2
element zeroLength 100 40 101 -mat 3 43 -dir 1 2
element zeroLength 101 41 102 -mat 2 43 -dir 1 2

equalDOF 1 62 1 2
equalDOF 2 63 1 2
equalDOF 3 64 1 2
equalDOF 4 65 1 2
equalDOF 5 66 1 2
equalDOF 6 67 1 2
equalDOF 7 68 1 2
equalDOF 8 69 1 2
equalDOF 9 70 1 2
equalDOF 10 71 1 2
equalDOF 11 72 1 2
equalDOF 12 73 1 2
equalDOF 13 74 1 2
equalDOF 14 75 1 2
equalDOF 15 76 1 2
equalDOF 16 77 1 2
equalDOF 17 78 1 2
equalDOF 18 79 1 2
equalDOF 19 80 1 2
equalDOF 20 81 1 2
equalDOF 21 82 1 2
equalDOF 22 83 1 2
equalDOF 23 84 1 2
equalDOF 24 85 1 2
equalDOF 25 86 1 2
equalDOF 26 87 1 2
equalDOF 27 88 1 2
equalDOF 28 89 1 2
equalDOF 29 90 1 2
equalDOF 30 91 1 2
equalDOF 31 92 1 2
equalDOF 32 93 1 2
equalDOF 33 94 1 2
equalDOF 34 95 1 2
equalDOF 35 96 1 2
equalDOF 36 97 1 2
equalDOF 37 98 1 2
equalDOF 38 99 1 2
equalDOF 39 100 1 2
equalDOF 40 101 1 2
equalDOF 41 102 1 2

setTime 0.0

pattern Plain 10 {Series -time {0 0 20 10000} -values {0 0 1 1} -factor 1} {
load 61 850.0 0.0 0.0
}

recorder Node -file Data/dispDepth.out -time -node 61 60 59 58 57 56 55 54 53\
52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27\
26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1 -dof 1 disp

integrator LoadControl 0.01
numberer RCM
system ProfileSPD
constraints Transformation
test NormDispIncr 1e-5 20 1
algorithm Newton
analysis Static

recorder display "OpenSees Pile Reaction Example" 10 10 1000 700 -wipe
prp 0 0 1000 # eye location in local coord sys defined by viewing system
vrp 0 0 0 # point on the view plane in global coord, center of local viewing system
vup 0 1 0 # dirn defining up direction of view plane
vpn 0 0 1 # direction of outward normal to view plane
viewWindow -80 150 -80 150
projection 1
fill 10
display 1 0 10

analysis Static

set startT [clock seconds]
analyze 100
set endT [clock seconds]
wipe all

[fmk]

typically the solver failing at the start is either incorrect material def, hanging nodes or insuficient # of boundary conditions .. for your model
the spring nodes have 3 dof, yet you only supply materials to act in 1 2 dirns .. i.e. dirn 3 is causing the problems on nodes > 60.

[ggokelly]

Thank you very much,

I fixed the code accroding to your comment.