def IFsteel_Hill48R_20160608(): """ IF steel parameters tuned for Hill48R Hardening parameters were tuned for bulge test data """ ## hardening model f_hrd = return_swift(n=0.255392940,m=5e-2,ks= 6.00162331e+02, e0=4.23405410e-04, qq=1e3) f_yld = wrapHill48R([2.092113652699876, 1.8999691982523326, 2.8779441147053473]) return Constitutive(f_yld=f_yld,f_hrd=f_hrd)
def IFsteel(): """ Return constitutive description for IF steel parameters """ ## hardening model f_hrd = return_swift(n=0.28985,m=5e-2, ks=518.968, e0=0.0007648, qq=1e3) ## yield function characterized by three r-values f_yld = wrapHill48R([2.2, 2.0, 2.9]) # r0, r45, r90 return Constitutive(f_yld=f_yld, f_hrd=f_hrd)
def IFsteel(): """ Return constitutive description for IF steel parameters """ ## hardening model f_hrd = return_swift(n=0.28985, m=5e-2, ks=518.968, e0=0.0007648, qq=1e3) ## yield function characterized by three r-values f_yld = wrapHill48R([2.2, 2.0, 2.9]) # r0, r45, r90 return Constitutive(f_yld=f_yld, f_hrd=f_hrd)
def main(ax=None): if type(ax).__name__ == 'NoneType': fig = plt.figure() ax = fig.add_subplot(111) import time t0 = time.time() ## two parameters from Jeong et al. Acta Mat 112, 2016 ## for the interstitial-free steel r0 = 2.20 r45 = 2.0 r90 = 2.9 ## the rest of parameters are assumed as 1. rb = 1. y0 = 1 y45 = 1 y90 = 1 yb = 1. h48g = wrapHill48R([r0, r45, r90]) yld2000 = wrapYLD(r=[r0, r45, r90, rb], y=[y0, y45, y90, yb], m=6, k=2) funcs = vm, h48g, yld2000 labs = ['von Mises', 'Hill48R', 'yld2000'] ls = ['-', '--', '-.', ':'] xs = [] ys = [] t_indv = [] for i in xrange(len(funcs)): t_indv0 = time.time() x, y = locus(funcs[i]) t_indv.append(time.time() - t_indv0) xs.append(x) ys.append(y) for i in xrange(len(funcs)): print 'time %s :' % labs[i], t_indv[i] print 'total time:', time.time() - t0 for i in xrange(len(funcs)): ax.plot(xs[i], ys[i], label=labs[i], ls=ls[i]) ax.legend(loc='best') fn = 'vm_check.pdf' print '%s has been saved' % fn ax.set_xlim(-0.25, 1.75) ax.set_ylim(-0.25, 1.75) ax.set_aspect('equal') try: fig.savefig(fn) except: pass
def main(ax=None): if type(ax).__name__=='NoneType': fig = plt.figure() ax = fig.add_subplot(111) import time t0=time.time() ## two parameters from Jeong et al. Acta Mat 112, 2016 ## for the interstitial-free steel r0 = 2.20 r45 = 2.0 r90 = 2.9 ## the rest of parameters are assumed as 1. rb =1. y0 =1 y45 =1 y90 =1 yb =1. h48g = wrapHill48R([r0,r45,r90]) yld2000 = wrapYLD(r=[r0,r45,r90,rb],y=[y0,y45,y90,yb],m=6,k=2) funcs = vm, h48g, yld2000 labs = ['von Mises', 'Hill48R','yld2000'] ls = ['-','--','-.',':'] xs=[];ys=[] t_indv=[] for i in xrange(len(funcs)): t_indv0= time.time() x,y = locus(funcs[i]) t_indv.append(time.time()-t_indv0) xs.append(x) ys.append(y) for i in xrange(len(funcs)): print 'time %s :'%labs[i], t_indv[i] print 'total time:',time.time()-t0 for i in xrange(len(funcs)): ax.plot(xs[i],ys[i],label=labs[i],ls=ls[i]) ax.legend(loc='best') fn='vm_check.pdf' print '%s has been saved'%fn ax.set_xlim(-0.25,1.75) ax.set_ylim(-0.25,1.75) ax.set_aspect('equal') try: fig.savefig(fn) except: pass
def IFsteel_Hill48R_20160608(): """ IF steel parameters tuned for Hill48R Hardening parameters were tuned for bulge test data """ ## hardening model f_hrd = return_swift(n=0.255392940, m=5e-2, ks=6.00162331e+02, e0=4.23405410e-04, qq=1e3) f_yld = wrapHill48R( [2.092113652699876, 1.8999691982523326, 2.8779441147053473]) return Constitutive(f_yld=f_yld, f_hrd=f_hrd)