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'