def IFsteel_yld2000_case1():
"""
Return constitutive decription for
IF steel parameters based on YLD2000-2D
"""
import mk.yieldFunction.tuneYld2000
f_hrd = return_swift(n=0.28985, m=5e-2, ks=518.968, e0=0.0007648, qq=1e3)
f_yld = mk.yieldFunction.tuneYld2000.H48toYld(rv=[2.2, 2.0, 2.9, 1.0], m=6)
return Constitutive(f_yld=f_yld, f_hrd=f_hrd)
def IFsteel_H48YR90_vpsc_20160722():
f_hrd = return_voce(a=479.00408591,
b0=339.71480479,
c=7.68395984,
b1=70.86783572,
m=5e-2,
qq=1e3)
f_yld = mk.yieldFunction.yf2.wrapHill48Y(r90=rv3_vpsc[0], ys=ys3_vpsc)
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_yld1_20160722():
f_hrd = return_voce(a=479.00408591,
b0=339.71480479,
c=7.68395984,
b1=70.86783572,
m=5e-2,
qq=1e3)
f_yld = mk.yieldFunction.tuneYld2000.H48toYld_withYS(
rv=[r_exp[0], r_exp[1], r_exp[2]], ys=y_exp, m=6)
return Constitutive(f_yld=f_yld, f_hrd=f_hrd)
def IFsteel_H48YR90_20160722():
f_hrd = return_voce(a=479.00408591,
b0=339.71480479,
c=7.68395984,
b1=70.86783572,
m=5e-2,
qq=1e3)
y_yld = mk.yieldFunction.yf2.wrapHill48Y(ys=[y_exp[0], y_exp[1], y_exp[2]],
r90=r_exp[2])
return Constitutive(f_yld=f_yld, f_hrd=f_hrd)
def IFsteel_yld2_20160722():
f_hrd = return_voce(a=479.00408591,
b0=339.71480479,
c=7.68395984,
b1=70.86783572,
m=5e-2,
qq=1e3)
f_yld = mk.yieldFunction.tuneYld2000.wrapYLD(
r=[rv3_vpsc[0], rv3_vpsc[1], rv3_vpsc[2], rb],
y=[ys3_vpsc[0], ys3_vpsc[1], ys3_vpsc[2], yb_vpsc],
m=6)
return Constitutive(f_yld=f_yld, f_hrd=f_hrd)
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 IsoMat():
"""
Return constitutive description for
Isotropic von Mises material
Hardening parameters were tuned for bulge test data (of the IF steel)
"""
## hardening model
f_hrd = return_swift(n=0.255392940,
m=5e-2,
ks=6.00162331e+02,
e0=4.23405410e-04,
qq=1e3)
## yield function
f_yld = VonMises
return Constitutive(f_yld=f_yld, f_hrd=f_hrd)
def IFsteel_yld2000_2d_1_20160608():
"""
IF steel parameters tuned for yld2000_2d
-- yld2000-2d model tuned by experimental r0, r45, r90.
rb, y0, y45, y90, and yb were tuned using Hill48R(r0,r45,r90)
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)
import mk.yieldFunction.tuneYld2000
f_yld = mk.yieldFunction.tuneYld2000.H48toYld(
rv=[2.092113652699876, 1.8999691982523326, 2.8779441147053473], m=6)
return Constitutive(f_yld=f_yld, f_hrd=f_hrd)