Pile and no result

Forum for OpenSees users to post questions, comments, etc. on the use of the OpenSees interpreter, OpenSees.exe

Moderators: silvia, selimgunay, Moderators

Post Reply
mohsenkh
Posts: 27
Joined: Wed Apr 21, 2010 7:40 am

Pile and no result

Post by mohsenkh »

Dear

I modeled a pile within a soil which consist four layers. When I analysis the model, there is no error. But there is no result. Can you help me?

Best Regards;

code of my model is:

wipe
model BasicBuilder -ndm 3 -ndf 6
#UNITS: N, M, S;

#node $nodeTag $x $y $z
node 1 0 0 0
node 2 0 0 -1
node 3 0 0 -2
node 4 0 0 -3
node 5 0 0 -4
node 6 0 0 -5
node 7 0 0 -6
node 8 0 0 -7
node 9 0 0 -8


#node $nodeTag $x $y $z
node 100 0 0 0
node 200 0 0 -1
node 300 0 0 -2
node 400 0 0 -3
node 500 0 0 -4
node 600 0 0 -5
node 700 0 0 -6
node 800 0 0 -7
node 900 0 0 -8

mass 1 756 756 756 756 756 756
mass 2 756 756 756 756 756 756
mass 3 756 756 756 756 756 756
mass 4 756 756 756 756 756 756
mass 5 756 756 756 756 756 756
mass 6 756 756 756 756 756 756
mass 7 756 756 756 756 756 756
mass 8 756 756 756 756 756 756
mass 9 756 756 756 756 756 756


#fix 100 1 1 1 1 1 1
fix 200 1 1 1 1 1 1
fix 300 1 1 1 1 1 1
fix 400 1 1 1 1 1 1
fix 500 1 1 1 1 1 1
fix 600 1 1 1 1 1 1
fix 700 1 1 1 1 1 1
fix 800 1 1 1 1 1 1
fix 900 1 1 1 1 1 1

uniaxialMaterial Steel02 1000 392.4e6 2.0e11 0.01 15 0.925 0.15 0.0005 0.01 0.0005 0.01

uniaxialMaterial Concrete01 2000 -24.525e6 -0.002 -18.0e6 -0.006

#Z=-1 PILE
uniaxialMaterial PySimple1 3000 2 96132.96 1.22e-3 0.1

#Z=-2
uniaxialMaterial PySimple1 3100 2 305363.52 2.11e-3 0.1

#Z=-3
uniaxialMaterial PySimple1 3200 1 443040 0.03 0.3

#Z=-4
uniaxialMaterial PySimple1 3300 1 513240 0.03 0.3

#Z=-5
uniaxialMaterial PySimple1 3400 1 582660 0.03 0.3

#Z=-6
uniaxialMaterial PySimple1 3500 2 4006482.48 5.44e-3 0.1

#z=-7
uniaxialMaterial PySimple1 3600 1 1350000 0.03 0.3

#Z=-8
uniaxialMaterial PySimple1 3700 1 1350000 0.03 0.3

#FRICTIONAL RESISTANCE
#uniaxialMaterial TzSimple1 $matTag $tzType $tult $z50

#Z=-1 PILE
uniaxialMaterial TzSimple1 4000 2 5948.34 127.0e-5

#Z=-2
uniaxialMaterial TzSimple1 4100 2 118966.92 127.0e-5

#Z=-3
uniaxialMaterial TzSimple1 4200 1 80199.86 1.86e-3

#Z=-4
uniaxialMaterial TzSimple1 4300 1 86143.67 1.86e-3

#Z=-5
uniaxialMaterial TzSimple1 4400 1 91087.88 1.86e-3

#Z=-6
uniaxialMaterial TzSimple1 4500 2 53734.12 127.0e-5

#Z=-7
uniaxialMaterial TzSimple1 4600 1 174473.4 1.86e-3

#Z=-8
uniaxialMaterial TzSimple1 4700 1 180527 1.86e-3

#Z=-8 PILE
uniaxialMaterial QzSimple1 5000 1 636172.51 7.8e-3

#SECTION OF PILES
section Fiber 6000 {
patch circ 2000 70 22 0.0 0.0 0.0 0.3 0.0 360.0
layer circ 1000 16 0.000616 0.0 0.0 0.236 0.0 360.0
}

geomTransf Corotational 1 0 1 0

#ELEMENT
#ELEMENT OF PILE
#element dispBeamColumn $eleTag $iNode $jNode $numIntgrPts $secTag $transfTag

#ELEMENT 1
element dispBeamColumn 10 1 2 5 6000 1 -mass 756

#ELEMENT 2
element dispBeamColumn 20 2 3 5 6000 1 -mass 756

#ELEMENT 3
element dispBeamColumn 30 3 4 5 6000 1 -mass 756

#ELEMENT 4
element dispBeamColumn 40 4 5 5 6000 1 -mass 756

#ELEMENT 5
element dispBeamColumn 50 5 6 5 6000 1 -mass 756

#ELEMENT 6
element dispBeamColumn 60 6 7 5 6000 1 -mass 756

#ELEMENT 7
element dispBeamColumn 70 7 8 5 6000 1 -mass 756

#ELEMENT 8
element dispBeamColumn 80 8 9 5 6000 1 -mass 756

#Z=-1 PILE
element zeroLength 9000 200 2 -mat 3000 3000 4000 -dir 1 2 3

#Z=-2
element zeroLength 9100 300 3 -mat 3100 3100 4100 -dir 1 2 3

#Z=-3
element zeroLength 9200 400 4 -mat 3200 3200 4200 -dir 1 2 3

#Z=-4
element zeroLength 9300 500 5 -mat 3300 3300 4300 -dir 1 2 3

#Z=-5
element zeroLength 9400 600 6 -mat 3400 3400 4400 -dir 1 2 3

#Z=-6
element zeroLength 9500 700 7 -mat 3500 3500 4500 -dir 1 2 3

#Z=-7
element zeroLength 9600 800 8 -mat 3600 3600 4600 -dir 1 2 3

#Z=-8
element zeroLength 9700 900 9 -mat 3700 3700 4700 -dir 1 2 3

pattern Plain 1 Linear {
load 1 5000 -8000 -1000 1300 23450 67540
}

recorder Node -file Node91.out -time -node 2 -dof 1 2 3 disp
recorder Node -file Node9.out -time -node 2 -dof 1 2 3 reaction
recorder Element -file Element 80.out -time ele 80 localForce

system Band General
constraints Transformation
numberer RCM
test NormDispIncr 1.0e-2 100
algorithm ModifiedNewton
integrator LoadControl 0.1
analysis Static
analyze 10
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: Pile and no result

Post by vesna »

You have a typo in definition of system and one of recorders.

recorder Element -file Element80.out -time ele 80 localForce

system BandGeneral

After you correct it you will see get some error messages telling you that your model is unstable. I looked at your model and I saw that you have 3D zero-lenght elements with 6DOFs that are defined by 3DOFs only. You either need to constrain end nodes in DOFs that were not defined or to define other DOFs as well.
Post Reply