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)