print("\nMinimum strain in steel fibers: ",(epsSMin*1E3),"E-3")
print("Minimum stress in steel fibers: ",(sgSMin/1e6),"E6")
print("Y coordinate of the steel fiber with minimum strain: ",(YepsSMin))
print("Z coordinate of the steel fiber with minimum strain: ",(ZepsSMin))
print("\nMaximum strain in steel fibers: ",(epsSMax*1E3),"E-3")
print("Maximum stress in steel fibers: ",(sgSMax/1e6),"E6")
print("Y coordinate of the steel fiber with maximum strain: ",(YepsSMax))
print("Z coordinate of the steel fiber with maximum strain: ",(ZepsSMax))
'''
# *Crack width calculation*
# depth of the effective area:
hceff = EC2_limit_state_checking.h_c_eff(depth_tot=depth,
depht_eff=depth - coverInf -
0.026 / 2.0,
depth_neutral_axis=abs(x))
# print('depth of the effective area: ',hceff,' m')
# effective reinforcement ratio
roseff = EC2_limit_state_checking.ro_eff(A_s=A_s, width=width, h_c_eff=hceff)
# print('effective reinforcement ratio: ',roseff)
# maximum crack spacing
srmax = EC2_limit_state_checking.s_r_max(k1=0.8,
k2=0.5,
k3=3.4,
k4=0.425,
cover=coverInf,
fiReinf=0.024,
ro_eff=roseff)
# print('maximum crack spacing: ',srmax,' m')
# mean strain in the concrete between cracks
# -*- coding: utf-8 -*-
import os
#Project directory structure
execfile("../env_config.py")
from postprocess import limit_state_data as lsd
from postprocess import RC_material_distribution
#from materials.ehe import EHE_limit_state_checking
from materials.ec2 import EC2_limit_state_checking
from solution import predefined_solutions
lsd.LimitStateData.envConfig = cfg
#Reinforced concrete sections on each element.
#reinfConcreteSections=RC_material_distribution.RCMaterialDistribution()
reinfConcreteSections.mapSectionsFileName = './mapSectionsReinforcementTenStiff.pkl'
reinfConcreteSections = RC_material_distribution.loadRCMaterialDistribution()
#Checking material for limit state.
limitStateLabel = lsd.freqLoadsCrackControl.label
lsd.freqLoadsCrackControl.controller = EC2_limit_state_checking.CrackStraightController(
limitStateLabel=lsd.freqLoadsCrackControl.label)
lsd.freqLoadsCrackControl.controller.analysisToPerform = predefined_solutions.plain_static_modified_newton
lsd.freqLoadsCrackControl.check(reinfConcreteSections)
reinfFibers = setsRCEl1.reinfFibers.fSet
concrSetFbEl1 = sccEl1.getFiberSets()["concrSetFbEl1"]
reinfSetFbEl1 = sccEl1.getFiberSets()["reinfSetFbEl1"]
concrFibers = [f for f in concrSetFbEl1]
x = sccEl1.getNeutralAxisDepth()
d = sccEl1.getEffectiveDepth()
h = sccEl1.getLeverArm()
As = setsRCEl1.tensionFibers.getArea(1.0)
print 'As= ', As
#maximum depth of the effective area:
hceff = EC2_limit_state_checking.h_c_eff(depth_tot=h,
depht_eff=abs(d),
depth_neutral_axis=abs(x))
print 'depth of the effective area: ', hceff, ' m'
#Aceff_EHE_gross=sccEl1.getGrossEffectiveConcreteArea(hceff)
Aceff = sccEl1.getNetEffectiveConcreteArea(hceff, 'reinfSetFbEl1', 15.0)
print 'effective concrete tension area: ', Aceff, ' m2'
ro_s_eff = As / Aceff #effective ratio of reinforcement
print 'effective ratio of reinforcement=', ro_s_eff
#maximum crack spacing
srmax = EC2_limit_state_checking.s_r_max(k1=0.8,
k2=0.5,
k3=3.4,
k4=0.425,
cover=cover,
fiReinf=0.024,
print "\nMinimum strain in steel fibers: ",(epsSMin*1E3),"E-3" print "Minimum stress in steel fibers: ",(sgSMin/1e6),"E6" print "Y coordinate of the steel fiber with minimum strain: ",(YepsSMin) print "Z coordinate of the steel fiber with minimum strain: ",(ZepsSMin) print "\nMaximum strain in steel fibers: ",(epsSMax*1E3),"E-3" print "Maximum stress in steel fibers: ",(sgSMax/1e6),"E6" print "Y coordinate of the steel fiber with maximum strain: ",(YepsSMax) print "Z coordinate of the steel fiber with maximum strain: ",(ZepsSMax) ''' # *Crack width calculation* #depth of the effective area: hceff=EC2_limit_state_checking.h_c_eff(depth_tot=depth,depht_eff=depth-cover-0.024/2.0,depth_neutral_axis=abs(x)) # print 'depth of the effective area: ',hceff,' m' #effective reinforcement ratio roseff=EC2_limit_state_checking.ro_eff(A_s=A_s,width=width,h_c_eff=hceff) # print 'effective reinforcement ratio: ',roseff #maximum crack spacing srmax=EC2_limit_state_checking.s_r_max(k1=0.8,k2=0.5,k3=3.4,k4=0.425,cover=cover,fiReinf=0.024,ro_eff=roseff) # print 'maximum crack spacing: ',srmax,' m' #mean strain in the concrete between cracks eps_cm=concrete.fctm()/concrete.E0()/2.0 #mean strain in the reinforcemen takin into account the effects of tension stiffening eps_sm=epsSMax #crack withs w_k=srmax*(eps_sm-eps_cm) # print 'crack widths: ',w_k*1e3, ' mm'