print 'x= ', x, ' d= ', d, ' h= ', h As=setsRCEl1.tensionFibers.getArea(1.0) #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' print 'As=', As 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=fiBott,ro_eff=ro_s_eff) 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 taking into account the effects of tension stiffening fReinfMax= setsRCEl1.reinfFibers.getFiberWithMaxStrain() epsSMax= fReinfMax.getMaterial().getStrain() # maximum strain among steel fibers eps_sm=epsSMax #crack withs w_k=srmax*(eps_sm-eps_cm) print 'crack widths: ',w_k*1e3, ' mm' dom=preprocessor.getDomain #step 2: only revert to start and new run
# *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
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 withs: ',w_k*1e3, ' mm')
xComp = -0.15899240709999138
epsCMinComp = -0.000677387534344379
epsSMaxComp = 0.001285519718131271
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' xComp= -0.206563714084 epsCMinComp= -0.000534072629205 sgCMinComp= -13882927.0767 epsSMaxComp= 0.000951911878497 sgSMaxComp= 190382375.699 hceffComp= 0.13
