#of a beam element
reinfSteel= ACI_materials.A615G60
reinfSteel.gammaS= 1.0/0.75
beamRCsect= element_section_map.RCSlabBeamSection(name='beamRCsect',sectionDescr='beam section',concrType=concr, reinfSteelType=reinfSteel,width= sectionGeometry.b,depth= sectionGeometry.h)
beamRCsect.dir1PositvRebarRows= def_simple_RC_section.LongReinfLayers([reinfLayer])
beamRCsect.dir1NegatvRebarRows= def_simple_RC_section.LongReinfLayers([reinfLayer])
beamRCsect.dir2PositvRebarRows= def_simple_RC_section.LongReinfLayers([reinfLayer])
beamRCsect.dir2NegatvRebarRows= def_simple_RC_section.LongReinfLayers([reinfLayer])
sections.append(beamRCsect)
# Spatial distribution of reinforced concrete
# sections (assign RC sections to elements).
reinfConcreteSectionDistribution.assign(elemSet=totalSet.getElements,setRCSects=beamRCsect)
#Checking shear.
lsd.shearResistance.controller= ACI_limit_state_checking.ShearController(limitStateLabel= lsd.shearResistance.label)
lsd.shearResistance.controller.analysisToPerform= predefined_solutions.plain_newton_raphson
lsd.normalStressesResistance.outputDataBaseFileName= 'resVerif'
outCfg= lsd.VerifOutVars(listFile='N',calcMeanCF='Y')
(FEcheckedModel,meanFCs)= reinfConcreteSectionDistribution.runChecking(lsd.shearResistance, matDiagType="d",threeDim= True,outputCfg=outCfg)
#print("mean FCs: ", meanFCs)
meanFC0Teor= 0.932218036004
ratio1= abs(meanFCs[0]-meanFC0Teor)/meanFC0Teor
meanFC1Teor= 0.932218036004
ratio2= abs(meanFCs[1]-meanFC1Teor)/meanFC1Teor
'''
lSqrtFck = concrete.getLambdaSqrtFck()
lSqrtFckRef = 63.246 * ACI_materials.toPascal
ratio1 = abs(lSqrtFck - lSqrtFckRef) / lSqrtFckRef
bw = 12 * 0.0254 # 12 in
d = 15 * 0.0254 # 15 in
section = defSimpleRCSection.RecordRCSimpleSection(name='test',
concrType=concrete,
reinfSteelType=reinfSteel,
width=bw,
depth=d / 0.9)
section.shReinfY = defSimpleRCSection.RecordShearReinforcement(
nShReinfBranches=1.0,
areaShReinfBranch=0.592 * 0.0254**2,
shReinfSpacing=1.0 * 0.3048)
shearController = lsc.ShearController(lsd.shearResistance.label)
shearController.setSection(section)
Vc = shearController.getVcNoShearRebars(Nd=0.0)
VcRef = 17.076 * 4.4482216e3 / 0.75 # kip->N
ratio2 = abs(Vc - VcRef) / VcRef
Vmax = shearController.getV_max(Nd=0.0)
VmaxRef = 85.381 * 4.4482216e3
ratio3 = abs(Vmax - VmaxRef) / VmaxRef
ratio4 = abs(0.0004774184 - section.shReinfY.getAs() *
shearController.effectiveDepth) / 0.0004774184
Vu = shearController.getVu(Nd=0.0)
VuRef = 50.4 * 4.4482216e3