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()