scSteel=sect1.reinfSteelType
steelDiag=scSteel.plotDesignStressStrainDiagram(preprocessor,path=report_graphics_outDir)
steelGrphFile=scSteel.materialName+'_design_stress_strain_diagram'
report.write('\\begin{center}\n')
report.write('\includegraphics[width=120mm]{'+report_graphics_outDir+steelGrphFile+'}\n')
report.write('\end{center}\n')
scConcr=sect1.concrType
concrDiag=scConcr.plotDesignStressStrainDiagram(preprocessor,path=report_graphics_outDir)
concrGrphFile=scConcr.materialName+'_design_stress_strain_diagram'
report.write('\\begin{center}\n')
report.write('\includegraphics[width=120mm]{'+report_graphics_outDir+concrGrphFile+'}\n')
report.write('\end{center}\n')
report.write('\\newpage\n\n')
#Section 1
# plotting of section geometric and mechanical properties
sect1inf=sectionReport.SectionInfoHASimple(preprocessor,sect1)
texFileName=report_graphics_outDir+sect1.sectionName+'.tex'
epsFileName=report_graphics_outDir+sect1.sectionName+'.eps'
sect1inf.writeReport(texFileName,epsFileName)
report.write('\input{'+texFileName+'}\n')
# plotting of interaction diagrams
diagNMy= sect1.defInteractionDiagramNMy(preprocessor)
grFileName=report_graphics_outDir+sect1.sectionName+'NMy'
plotIntDiag(diag=diagNMy,title=sect1.sectionName+ ' N-My interaction diagram',xAxLab='My [kNm]',yAxLab='N [kN]',grFileNm=grFileName,reportFile=report)
diagNMz= sect1.defInteractionDiagramNMz(preprocessor)
grFileName=report_graphics_outDir+sect1.sectionName+'NMz'
plotIntDiag(diag=diagNMz,title=sect1.sectionName+ ' N-Mz interaction diagram',xAxLab='Mz [kNm]',yAxLab='N [kN]',grFileNm=grFileName,reportFile=report)
#Section 2
# plotting of section geometric and mechanical properties
sect2inf=sectionReport.SectionInfoHASimple(preprocessor,sect2)
texFileName=report_graphics_outDir+sect2.sectionName+'.tex'
nominalCover=0.25 - 0.19)
sccData.positvRebarRows = [posReb]
zinf = sccData.h / 2.0
zsup = -sccData.h / 2.0
prueba = xc.FEProblem()
#prueba.logFileName= "/tmp/borrar.log" # Don't print warnings.
prueba.errFileName = "/tmp/borrar.err" # Don't print errors.
preprocessor = prueba.getPreprocessor
print "divIJ= ", sccData.nDivIJ, "divJK= ", sccData.nDivJK
#sccData.nDivIJ= 100
#sccData.nDivJK= 100
sccData.defRCSimpleSection(preprocessor, 'd')
si = sr.SectionInfoHASimple(preprocessor, sccData)
si.writeReport('./prueba.tex', './prueba.eps')
Agross = sccData.b * sccData.h
IyGross = 1.0 / 12.0 * sccData.b * sccData.h**3
IzGross = 1.0 / 12.0 * sccData.h * sccData.b**3
n = si.scc.getSteelEquivalenceCoefficient(preprocessor)
AsBarsUp = 6 * areaFi40
AsBarsDown = 6 * areaFi6
AsTeor = AsBarsUp + AsBarsDown
As = si.areaMainReinforcement
barsUpTeorCenterOfMass = sccData.getNegRowsCGcover() - sccData.h / 2.0
barsUpCenterOfMass = sccData.negReinfLayers[0].getCenterOfMass()
#barsUpCenterOfMass= sccData.negReinfLayer.getCenterOfMass()
zBarsUp = barsUpCenterOfMass[1] - si.GH[1]
barsDownTeorCenterOfMass = sccData.h / 2.0 - sccData.getPosRowsCGcover()
