Example pile
From Nemesis
Contents |
Main file: pile.slv
import string
try:
nsfile = open("pile[ns].dat",'r')
npfile = open("pile[np].dat",'r')
ncfile = open("pile[nc].dat",'r')
esfile = open("pile[es].dat",'r')
epfile = open("pile[ep].dat",'r')
except IOError:
sys.exit(0)
domain.dim(3)
db.SQLite("pile")
zs=[0.,1.,2.,3.,4.,5.,6.,7.5,8.0,8.5,9.0,9.5,10.0]
zp=[0.,1.,2.,3.,4.,5.,6.,7.5,8.0,8.5,9.0,9.5,10.0,10.5,11.0,11.5,12.0]
zc=[0.,1.,2.,3.,4.,5.,6.,7.5,8.0,8.5,9.0,9.5,10.0,10.5,11.0,11.5,12.0]
###############################################################################
# nodes
###############################################################################
nd=[]
# soil nodes
temp=[]
for line in nsfile.readlines():
words=string.split(line)
temp.append([int(words[0]),float(words[1]),float(words[2])])
for i in range(len(zs)):
for j in range(len(temp)):
nd.append([(i+1)*1000+temp[j][0],temp[j][1],temp[j][2],zs[i]])
# pile nodes
temp=[]
for line in npfile.readlines():
words=string.split(line)
temp.append([int(words[0]),float(words[1]),float(words[2])])
for i in range(len(zp)):
for j in range(len(temp)):
nd.append([(i+1)*1000+temp[j][0],temp[j][1],temp[j][2],zp[i]])
# common nodes
temp=[]
for line in ncfile.readlines():
words=string.split(line)
temp.append([int(words[0]),float(words[1]),float(words[2])])
for i in range(len(zc)):
for j in range(len(temp)):
nd.append([(i+1)*1000+temp[j][0],temp[j][1],temp[j][2],zc[i]])
###############################################################################
# elems
###############################################################################
el=[]
# soil elems
temp=[]
for line in esfile.readlines():
words=string.split(line)
temp.append([int(words[0]),int(words[1]),int(words[2]),int(words[3])])
for i in range(len(zs)-1):
for j in range(len(temp)):
el.append([(i+1)*1000+temp[j][0],(i+1)*1000+temp[j][0],(i+1)*1000+temp[j][1],(i+1)*1000+temp[j][2],(i+1)*1000+temp[j][3],(i+2)*1000+temp[j][0],(i+2)*1000+temp[j][1],(i+2)*1000+temp[j][2],(i+2)*1000+temp[j][3],1])
# pile elems
temp=[]
for line in epfile.readlines():
words=string.split(line)
temp.append([int(words[0]),int(words[1]),int(words[2]),int(words[3])])
for i in range(len(zp)-1):
for j in range(len(temp)):
el.append([(i+1)*1000+temp[j][0],(i+1)*1000+temp[j][0],(i+1)*1000+temp[j][1],(i+1)*1000+temp[j][2],(i+1)*1000+temp[j][3],(i+2)*1000+temp[j][0],(i+2)*1000+temp[j][1],(i+2)*1000+temp[j][2],(i+2)*1000+temp[j][3],2])
material.plainStrainElasticIsotropic(1, 1E3, 0., 0.)
material.plainStrainElasticIsotropic(2, 1E6, 0., 0.)
for i in range(len(nd)):
node.add(nd[i][0],nd[i][1],nd[i][2],nd[i][3])
for i in range(len(el)):
element.brick8d(el[i][0],el[i][1],el[i][2],el[i][3],el[i][4],el[i][5],el[i][6],el[i][7],el[i][8],el[i][9])
###############################################################################
# constraints
###############################################################################
for i in range(len(nd)):
x=nd[i][1]
y=nd[i][2]
z=nd[i][3]
if z==0:
node.fix(nd[i][0],1)
node.fix(nd[i][0],2)
node.fix(nd[i][0],3)
elif x==-5.8 or x==5.8:
node.fix(nd[i][0],1)
elif y== 0.0 or y==6.5:
node.fix(nd[i][0],2)
###############################################################################
# loadcase
###############################################################################
lc.define(1,"horizontal")
h=len(zp)
p=10
load.node(h*1000+111,1,p)
load.node(h*1000+121,1,p)
load.node(h*1000+131,1,p)
load.node(h*1000+141,1,p)
load.node(h*1000+151,1,p)
load.node(h*1000+161,1,p)
load.node(h*1000+171,1,p)
load.node(h*1000+181,1,p)
load.node(h*1000+191,1,p)
load.node(h*1000+210,1,p)
load.node(h*1000+211,1,p)
load.node(h*1000+212,1,p)
load.node(h*1000+213,1,p)
load.node(h*1000+290,1,p)
load.node(h*1000+220,1,p)
load.node(h*1000+214,1,p)
load.node(h*1000+215,1,p)
load.node(h*1000+216,1,p)
load.node(h*1000+280,1,p)
load.node(h*1000+230,1,p)
load.node(h*1000+240,1,p)
load.node(h*1000+250,1,p)
load.node(h*1000+260,1,p)
load.node(h*1000+270,1,p)
###############################################################################
# loadcase
###############################################################################
analysis.static()
reorder.rCM()
soe.band()
analysis.run(1,1)
db.store("pile_soln")
db.exportToVtk("pile_soln")
Auxiliary File: pile[ns].dat
112 -0.8 0 113 -1.4 0 114 -2.2 0 115 -3.2 0 116 -4.4 0 117 -5.8 0 122 -0.8 0.4483 123 -1.4 0.7845 124 -2.2 1.2328 125 -3.2 1.7931 126 -4.4 2.4655 127 -5.8 3.25 132 -0.8 0.8966 133 -1.4 1.569 134 -2.2 2.4655 135 -3.2 3.5862 136 -4.4 4.931 137 -5.8 6.5 142 -0.4164 0.8966 143 -0.7504 1.569 144 -1.1958 2.4655 145 -1.7525 3.5862 146 -2.4206 4.931 147 -3.2 6.5 152 0 0.8966 153 0 1.569 154 0 2.4655 155 0 3.5862 156 0 4.931 157 0 6.5 162 0.4164 0.8966 163 0.7504 1.569 164 1.1958 2.4655 165 1.7525 3.5862 166 2.4206 4.931 167 3.2 6.5 172 0.8 0.8966 173 1.4 1.569 174 2.2 2.4655 175 3.2 3.5862 176 4.4 4.931 177 5.8 6.5 182 0.8 0.4483 183 1.4 0.7845 184 2.2 1.2328 185 3.2 1.7931 186 4.4 2.4655 187 5.8 3.25 192 0.8 0 193 1.4 0 194 2.2 0 195 3.2 0 196 4.4 0 197 5.8 0
Auxiliary File: pile[np].dat
210 -0.2 0 211 -0.0843 0 212 0 0 213 0.0843 0 290 0.2 0 220 -0.2 0.1121 214 -0.0843 0.1121 215 0 0.1121 216 0.0843 0.1121 280 0.2 0.1121 230 -0.2 0.2281 240 -0.0843 0.2281 250 0 0.2281 260 0.0843 0.2281 270 0.2 0.2281
Auxiliary File: pile[nc].dat
111 -0.4 0 121 -0.349 0.1955 131 -0.2644 0.3002 141 -0.1543 0.369 151 0 0.4 161 0.1543 0.369 171 0.2644 0.3002 181 0.349 0.1955 191 0.4 0
