frontFillDepth= 1.0
deadLoad= loadCaseManager.setCurrentLoadCase('deadLoad')
wall.createDeadLoad(heelFillDepth= wall.stemHeight,toeFillDepth= frontFillDepth,rho= backFillSoilModel.rho, grav= gravity)
# Dead load. Earth pressure.
Ka= backFillSoilModel.Ka()
zGroundBackFill= 0.0 #Back fill
backFillPressureModel= earth_pressure.EarthPressureModel(K= Ka, zGround= zGroundBackFill, gammaSoil= gSoil, zWater= -1e3, gammaWater= 1000*gravity)
wall.createBackFillPressures(backFillPressureModel, Delta= backFillDelta)
zGroundFrontFill= zGroundBackFill-wall.stemHeight+frontFillDepth #Front fill
frontFillPressureModel= earth_pressure.EarthPressureModel(K= Ka, zGround= zGroundFrontFill, gammaSoil= gSoil, zWater= -1e3, gammaWater= 1000*gravity)
wall.createFrontFillPressures(frontFillPressureModel)
#Live load. Crowd loading.
crowdLoad= loadCaseManager.setCurrentLoadCase('crowdLoad')
crowdEarthPressure= earth_pressure.StripLoadOnBackfill(qLoad= 4e3,zLoad= 0.0, distWall= 2.5/2.0, stripWidth= 2.5)
wall.createPressuresFromLoadOnBackFill(crowdEarthPressure, Delta= backFillDelta)
wall.createLoadOnTopOfStem(xc.Vector([-3e3,0.0,3.6e3]))
#Live load. Rail traffic load.
distRailCLWall= 4.12 #Distance from the center line of the rail track to the wall.
railLoad= loadCaseManager.setCurrentLoadCase('railLoad')
railLoadEarthPressure= earth_pressure.StripLoadOnBackfill(qLoad= 62.5e3,zLoad= -0.6, distWall= distRailCLWall, stripWidth= 2.0)
wall.createPressuresFromLoadOnBackFill(railLoadEarthPressure, Delta= backFillDelta)
#Nosing load
fNosingLoad= 60e3
nosingLoadLength= 1.8+(distRailCLWall-1.0)
nosingLoadSurface= nosingLoadLength*2.0
qNosingLoad= fNosingLoad/nosingLoadLength
print 'qNosingLoad= ', qNosingLoad
#CentrifugalLoad.
soilData=soil01,
vDir=xc.Vector([0, 1, 0]))
earthPressLoadColumn = loads.EarthPressLoad(name='earthPressLoadColumn',
xcSet=columnZconcr,
soilData=soil01,
vDir=xc.Vector([0, 1, 0]))
soil02 = ep.EarthPressureModel(zGround=zList[lastZpos],
zBottomSoils=[-10],
KSoils=[0.001],
gammaSoils=[densSoil * grav],
zWater=0.05,
gammaWater=densWater * grav)
stripL01 = ep.StripLoadOnBackfill(qLoad=2e5,
zLoad=zList[lastZpos],
distWall=1.5,
stripWidth=1.2)
earthPColumnStrL = loads.EarthPressLoad(name='earthPColumnStrL',
xcSet=columnZconcr,
soilData=None,
vDir=xc.Vector([0, 1, 0]))
earthPColumnStrL.stripLoads = [stripL01]
lineL01 = ep.LineVerticalLoadOnBackfill(qLoad=1e5,
zLoad=zList[lastZpos],
distWall=1.0)
earthPColumnLinL = loads.EarthPressLoad(name='earthPColumnLinL',
xcSet=columnZconcr,
soilData=None,
vDir=xc.Vector([0, 1, 0]))
earthPColumnLinL.lineLoads = [lineL01]
ep_muri_z4_desbal= loads.EarthPressLoad(name= 'ep_muri_z4', xcSet=muri_z4_set,soilData=soilZ4_desbal, vDir=xc.Vector([0,1,0])) ep_murd_z1_1= loads.EarthPressLoad(name= 'ep_murd_z1_1', xcSet=murd_z1_1_set,soilData=soilZ1_1, vDir=xc.Vector([0,-1,0])) ep_murd_z1_2= loads.EarthPressLoad(name= 'ep_murd_z1_2', xcSet=murd_z1_2_set,soilData=soilZ1_2, vDir=xc.Vector([0,-1,0])) ep_murd_z1_3= loads.EarthPressLoad(name= 'ep_murd_z1_3', xcSet=murd_z1_3_set,soilData=soilZ1_3, vDir=xc.Vector([0,-1,0])) ep_murd_z2= loads.EarthPressLoad(name= 'ep_murd_z2', xcSet=murd_z2_set,soilData=soilZ2, vDir=xc.Vector([0,-1,0])) ep_murd_z3= loads.EarthPressLoad(name= 'ep_murd_z3', xcSet=murd_z3_set,soilData=soilZ3, vDir=xc.Vector([0,-1,0])) ep_murd_z4= loads.EarthPressLoad(name= 'ep_murd_z4', xcSet=murd_z4_set,soilData=soilZ4, vDir=xc.Vector([0,-1,0])) ep_murd_z4_desbal= loads.EarthPressLoad(name= 'ep_murd_z4', xcSet=murd_z4_set,soilData=soilZ4_desbal, vDir=xc.Vector([0,-1,0])) # cargas sobre el terreno qb_Qptren_1_set=gridGeom.getSetSurfMultiRegion(lstIJKRange=[gm.IJKRange((pxcarr1,pyhast1,0),(pxcarr2,pyhast1,lastZpos))],setName='qb_Qptren_1_set') dist_carr=xList[pxcarr2]-xList[pxcarr1] strip_QPtren=ep.StripLoadOnBackfill(qLoad=2*Qptren_carr/L_Qptren_carr/dist_carr, zLoad=zGroundZ3,distWall=0, stripWidth=L_Qptren_carr) qb_Qptren_1= loads.EarthPressLoad(name= 'qb_Qptren_1', xcSet=qb_Qptren_1_set,soilData=None, vDir=xc.Vector([0,1,0])) qb_Qptren_1.stripLoads=[strip_QPtren] qb_Qptren_2_set=gridGeom.getSetSurfMultiRegion(lstIJKRange=[gm.IJKRange((pxcarr3,pyhast1,0),(pxcarr4,pyhast1,lastZpos))],setName='qb_Qptren_2_set') qb_Qptren_2= loads.EarthPressLoad(name= 'qb_Qptren_2', xcSet=qb_Qptren_2_set,soilData=None, vDir=xc.Vector([0,1,0])) qb_Qptren_2.stripLoads=[strip_QPtren] qb_qstren_set=gridGeom.getSetSurfMultiRegion(lstIJKRange=[gm.IJKRange((pxcan1,pyhast1,0),(pxcan2,pyhast1,lastZpos))],setName='qb_qstren_set') strip_qstren=ep.StripLoadOnBackfill(qLoad=scterr, zLoad=zGroundZ3,distWall=0, stripWidth=L_Qptren_carr) qb_qstren=loads.EarthPressLoad(name= 'qb_qstren', xcSet=qb_qstren_set,soilData=None, vDir=xc.Vector([0,1,0])) qb_qstren.stripLoads=[strip_qstren] strip_paseos=ep.StripLoadOnBackfill(qLoad=qnferr, zLoad=zGroundZ2,distWall=0, stripWidth=20) qb_paseos_i_set=gridGeom.getSetSurfMultiRegion(lstIJKRange=[gm.IJKRange((pxmnp1,pyhast1,0),(pxmxp1,pyhast1,lastZpos)),gm.IJKRange((pxmnp2,pyhast1,0),(pxmxp2,pyhast1,lastZpos))],setName='qb_paseos_i_set') qb_paseos_i=loads.EarthPressLoad(name= 'qb_paseos_i', xcSet=qb_paseos_i_set,soilData=None, vDir=xc.Vector([0,1,0]))
#Live load: road traffic load.
cLC = loadCaseManager.setCurrentLoadCase('roadTrafficLoad')
roadTrafficLoadEarthPressure = earth_pressure.UniformLoadOnBackfill(
K=K0, qLoad=11.02e3)
for e in side_b_elements.getElements:
elemCentroid = e.getPosCentroid(True)
localKVector = e.getCoordTransf.getG3Vector
pressure = -roadTrafficLoadEarthPressure.getPressure(
elemCentroid.z) * localKVector
e.vector3dUniformLoadGlobal(pressure) #SIA 261:2014 table 8
#Live load: rail traffic load.
railLoad = loadCaseManager.setCurrentLoadCase('LM1')
distRailCLWall = 4.5 #Distance from the center line of the rail track to the wall
railLoadEarthPressure = earth_pressure.StripLoadOnBackfill(
qLoad=50e3, zLoad=10.23 - 0.7, distWall=distRailCLWall, stripWidth=3.0)
for e in side_a_elements.getElements:
elemCentroid = e.getPosCentroid(True)
localKVector = e.getCoordTransf.getG3Vector
pressure = railLoadEarthPressure.getPressure(elemCentroid.z) * localKVector
e.vector3dUniformLoadGlobal(pressure) #SIA 261:2014 table 8
derailmentLoad = loadCaseManager.setCurrentLoadCase('DLM1')
derailmentLoadStripWidth = 0.45
platformWidth = 2.5
distDerailmentLoadWall = platformWidth + derailmentLoadStripWidth / 2.0 #Distance from the center line of the derailment load to the wall
derailmentLoadEarthPressure = earth_pressure.StripLoadOnBackfill(
qLoad=145e3 / derailmentLoadStripWidth,
zLoad=10.23 - 0.4,
distWall=distDerailmentLoadWall,
loadVector=xc.Vector(
[0, 0, -grav * densrell * (zGround - zZap - cantoZap / 2.)]))
SCep_zap = loads.UniformLoadOnSurfaces(name='rell_zap',
xcSet=zapTrasdos,
loadVector=xc.Vector([0, 0,
-qunifTerr]))
# empuje del terreno
soil = ep.EarthPressureModel(K=K0,
zGround=zGround,
gammaSoil=densrell * grav,
zWater=-10.0,
gammaWater=grav)
qunifTerr = ep.StripLoadOnBackfill(qLoad=qunifTerr,
zLoad=zGround,
distWall=0,
stripWidth=10)
ep_aleti = loads.EarthPressLoad(name='ep_aleti',
xcSet=aleti,
soilData=soil,
vDir=xc.Vector([-1, 0, 0]))
ep_aletd = loads.EarthPressLoad(name='ep_aletd',
xcSet=aletd,
soilData=soil,
vDir=xc.Vector([1, 0, 0]))
ep_murestr = loads.EarthPressLoad(name='ep_murestr',
xcSet=murestr,
soilData=soil,
vDir=xc.Vector([0, -1, 0]))