report.write('\\begin{center}\n')
report.write('\includegraphics[width=120mm]{' +
report_graphics_outDir + steelGrphFile + '}\n')
report.write('\end{center}\n')
scConcr = sect1.fiberSectionParameters.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 = section_report.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'
nominalCover=0.25 - 0.19)
sccData.positvRebarRows = def_simple_RC_section.LongReinfLayers([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.defRCSection(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()
