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)
