Soil-Pile Interaction

[mohfazelimm]

I'm interested in Soil-Pile interaction. I modeled a very simple problem (script is below) which has a quad element (as soil), a zerolength element (as interface) and an elasticBeamColumn (as pile). I use ndf 2 for quad element and ndf 3 for two other elements. The i and j nodes of zerolength have idential coordinate(node 3 & 5) and so OpenSees gives this error: “Nodes 3 and 5 have differing dof at ends for zerolength 5”. I know the reason of the error but don’t know how to correct it. I hope to know the answer and good source about quad element +beam column element.
Thanks.

--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
wipe
model BasicBuilder -ndm 2 -ndf 2

#SOIL#####
################)
node 1 0 0
node 2 1 0
node 3 1 1
node 4 0 1

nDMaterial ElasticIsotropic 1 200e3 .25
set thick .02
element quad 1 1 2 3 4 $thick PlaneStrain 1

#BC#####
################)
fix 1 1 1
fix 2 1 1

#PILE#####
################)

model BasicBuilder -ndm 2 -ndf 3

node 5 1 1
node 6 1 3
geomTransf Linear 1

set A .2
set E 200e6
set Iz .005
element elasticBeamColumn 2 5 6 $A $E $Iz 1

#INTERFACE#####
################)

uniaxialMaterial Elastic 1 200e3
element zeroLength 3 3 5 -mat 1 -dir 1

#RECORD#####
################)
recorder Node -file Node.out -time -node 3 5 6 -dof 1 disp

#LOADING#####
################)

set accl "Series -dt .02 -filePath acc.txt -factor 9.8"; # time series
pattern UniformExcitation 2 1 -accel $accl; # create uniform excitation

#ANALYZE#####
################)
integrator LoadControl 0.01
test EnergyIncr 1.0e-6 100 2
algorithm Newton
numberer RCM
constraints Plain
system ProfileSPD
analysis Static
analyze 100

[Kadenjeld]

Прелогаем отдых для настоящих мужчин -
[url=http://intimland.website.org/]Заказать проститутку в Москве[/url]