concrDiagram.commitState() sigmaLst.append([eps,concrDiagram.getStress()*1e-6]) print eps,' ',concrDiagram.getStress()*1e-6 concrDiagram.revertToStart() quit() # # plotting of section geometric and mechanical properties # from materials.sections.fiber_section import plotSectionGeometry as pg # pg.plotSectionGeometry(geomSectFibers,'./secEx74.eps') #sectInf=sectionReport.SectionInfo(preprocessor=preprocessor,section=sctFibers) #Obtains section parameters for report #sectInf.writeReport(archTex='./secEx74.tex', pathFigura='./secEx74.eps') # #plot cross-section strains and stresses from postprocess import utils_display # ''' # fiberSet: set of fibers to be represented # title: general title for the graphic # fileName: name of the graphic file (defaults to None: no file generated) # nContours: number of contours to be generated (defaults to 100) # ''' # #utils_display.plotStressStrainFibSet(fiberSet=fibersSccEl1,title='cross-section concrete fibers',fileName='problem.jpeg') utils_display.plotStressStrainFibSet(fiberSet=setsRCEl1.concrFibers.fSet,title='cross-section concrete fibers',fileName='problem.jpeg') # #report of the section material # sectParam=sectionReport.SectionInfo(preprocessor=preprocessor,sectName='example_7.5_EC2W_0.3mm',sectDescr='Test example 7.5 EC2 Worked examples - $w_k \\approx 0.3 mm$. Section definition',concrete=concrete,rfSteel=rfSteel,concrDiag=concrDiagram,rfStDiag=steelDiagram,geomSection=geomSectFibers,width=width,depth=depth) #Obtains section parameters for report # sectParam.writeReport(archTex='figures/sections/secEx75A.tex', pathFigura='figures/sections/secEx75A.eps')
stressStrainDiag.show(plt) for f in setsRCEl1.concrFibers.fSet: print 'y= ', f.getLocY(), 'z= ',f.getLocZ(), 'stress =', f.getForce()/f.getArea()*1e-6, 'strain= ',f.getStrain() ''' # #plot cross-section strains and stresses from postprocess import utils_display # ''' # fiberSet: set of fibers to be represented # title: general title for the graphic # fileName: name of the graphic file (defaults to None: no file generated) # nContours: number of contours to be generated (defaults to 100). If # nContours=0 or nContours=None, then each fiber is represented with a # colored circle. # # ''' #utils_display.plotStressStrainFibSet(fiberSet=fibersSccEl1,title='cross-section fibers',fileName=None,nContours=None) utils_display.plotStressStrainFibSet(fiberSet=setsRCEl1.concrFibers.fSet, title='cross-section concrete fibers', fileName=None, nContours=None, pointSize=50) utils_display.plotStressStrainFibSet( fiberSet=setsRCEl1.reinfFibers.fSet, title='cross-section reinforcement fibers', fileName=None, nContours=None, pointSize=100)