wall.reinforcement.setArmature(7,A10_15.getCopy('C'))
wall.reinforcement.setArmature(8,D1618_15.getCopy('B'))
wall.reinforcement.setArmature(11,A14_15.getCopy('B'))
wallFEModel= wall.createFEProblem('Retaining wall A')
preprocessor= wallFEModel.getPreprocessor
nodes= preprocessor.getNodeHandler
#Soil
kS= 40e6 #Module de réaction du sol (estimé).
#print 'kS= ', kS/1e6
kX= typical_materials.defElasticMaterial(preprocessor, "kX",kS/10.0)
kY= typical_materials.defElasticMaterial(preprocessor, "kY",kS)
#kY= typical_materials.defElastNoTensMaterial(preprocessor, "kY",kS)
backFillSoilModel= ep.RankineSoil(phi= math.radians(32),rho= 2000) #Characteristic values.
backFillDelta= 0.0#2.0/3.0*backFillSoilModel.phi
hi= [0.65,100.0]
rhoi= [2000,2100]
phii= [math.radians(34),math.radians(28)]
ci= [3e3,7e3]
stratifiedSoil= fcs.StratifiedSoil(hi,rhoi,phii,ci)
foundationSoilModel= stratifiedSoil.getEquivalentSoil(Beff= 2.5,gMPhi= 1.2,gMc= 1.5) #Design values.
#Mesh.
wall.genMesh(nodes,[kX,kY])
#Sets.
totalSet= preprocessor.getSets.getSet("total")
# -*- coding: utf-8 -*-
from __future__ import print_function
from geotechnics import earth_pressure
import math
Kp= earth_pressure.RankineSoil(math.pi/6.0).Kp()
#print('Kp= ', Kp)
ratio1= abs(Kp-3)/3.0
import os
fname= os.path.basename(__file__)
if (ratio1<1e-12):
print("test ",fname,": ok.")
else:
print("test ",fname,": ERROR.")
wall.footingReinforcement.setReinforcement(9,A16_10.getCopy('A')) # ln. sup. zapata
wall.footingReinforcement.setReinforcement(7,A20_10.getCopy('A')) # tr. inf. zapata
wall.footingReinforcement.setReinforcement(8,A16_20.getCopy('A')) # ln. inf. zapata
wall.stemReinforcement.setReinforcement(6,A12_20.getCopy('A')) #coronación
wallFEModel= wall.createFEProblem('Retaining wall '+sectionName)
preprocessor= wallFEModel.getPreprocessor
nodes= preprocessor.getNodeHandler
#Soils
kX= typical_materials.defElasticMaterial(preprocessor, "kX",kS/10.0)
kY= typical_materials.defElasticMaterial(preprocessor, "kY",kS)
#kY= typical_materials.defElastNoTensMaterial(preprocessor, "kY",kS)
#Backfill soil properties
backFillSoilModel= ep.RankineSoil(phi= math.radians(phiS),rho= rhoS) #Characteristic values.
#Foundation stratified soil properties
stratifiedSoil= fcs.StratifiedSoil(hi,rhoi,phii,ci)
foundationSoilModel= stratifiedSoil.getEquivalentSoil(Beff= 5,gMPhi= 1.2,gMc= 1.5) #Design values.
#Mesh.
wall.genMesh(nodes,[kX,kY])
#Sets.
totalSet= preprocessor.getSets.getSet("total")
#Actions.
loadCaseManager= load_cases.LoadCaseManager(preprocessor)
loadCaseNames= ['selfWeight','earthPress','earthPressAcc']
# ^ y
# |
# |
# B
# +------+
# | |
# | |
# | |
# | | H
# | |
# | |
# | |
# +------+ ----> x
#
soil = earth_pressure.RankineSoil(phi=fillPhi, rho=gammaSoil)
k0 = 1 - math.sin(fillPhi) # earth pressure at rest.
K = k0 #soil.Ka()
print 'K= ', K
x = list()
earth_pressure = list()
Fv = 0.0 #Friction.
for i in range(0, 11):
stepHeight = H / 10.0
z = i * stepHeight
sigma_h = K * gammaSoil * z
fv = sigma_h * math.tan(delta) * stepHeight
Fv += fv
x.append(z)
earth_pressure.append(sigma_h)