def ex(ifig=50,
exp_ref=['exp_dat/Bsteel/EXP_BULGE_JINKIM.txt',
'exp_dat/Bsteel/uni/avgstr_000.txt'],
exp_lab=['Exp bulge','Exp uniaxial'],
mod_ref='STR_STR.OUT'):
"""
Arguments
=========
exp_ref = ['exp_dat'/Bsteel/EXP_BULGE_JINKIM.txt']
exp_lab = ['Exp bulge', 'Exp uniaxial')
mod_ref = 'STR_STR.OUT'
"""
fig=plt.figure(ifig);ax=fig.add_subplot(111)
myrs = ResidualStress()
strains = myrs.straine #(20,6)
# strains = myrs.strainm #(4,6)
strain_eff = strains.T[0]+strains.T[1]
strain_eff = strain_eff.T
stress = []
for i in range(len(strain_eff)):
dum = myrs.analysis(iopt=0,istp=i)
stress.append((dum[0] + dum[1])/2.)
ax.plot(strain_eff,stress,'-x',label='Model SF/IG/ehkl')
stress = []
for i in range(len(strain_eff)):
dum = myrs.analysis(iopt=1,istp=i)
stress.append((dum[0] + dum[1])/2.)
ax.plot(strain_eff,stress,'-x',label='Exp SF/IG/ehkl')
stress = []
for i in range(len(strain_eff)):
dum = myrs.analysis(iopt=2,istp=i)
stress.append((dum[0] + dum[1])/2.)
ax.plot(strain_eff,stress,'-x',label='Model SF/IG + Exp ehkl')
for i in range(len(exp_ref)):
exp_ref = exp_ref[i]
x,y=np.loadtxt(exp_ref).T
ax.plot(x,y,'--',label=exp_lab[i])
dum=np.loadtxt(mod_ref,skiprows=1).T;
x=dum[2]+dum[3];y=(dum[8]+dum[9])/2.;
ax.plot(x,y,'--',label='EVPSC biaxial')
__deco__(ax,iopt=3)
ax.legend(loc='best',fancybox=True).get_frame().set_alpha(0.5)
fig.tight_layout()
plt.show()
return stress
def ex(ifig=50,
exp_ref=['exp_dat/Bsteel/EXP_BULGE_JINKIM.txt',
'exp_dat/Bsteel/uni/avgstr_000.txt'],
exp_lab=['Exp bulge','Exp uniaxial'],
mod_ref='STR_STR.OUT'):
"""
Arguments
=========
exp_ref = ['exp_dat'/Bsteel/EXP_BULGE_JINKIM.txt']
exp_lab = ['Exp bulge', 'Exp uniaxial')
mod_ref = 'STR_STR.OUT'
"""
fig=plt.figure(ifig);ax=fig.add_subplot(111)
myrs = ResidualStress()
strains = myrs.straine #(20,6)
# strains = myrs.strainm #(4,6)
strain_eff = strains.T[0]+strains.T[1]
strain_eff = strain_eff.T
stress = []
for i in xrange(len(strain_eff)):
dum = myrs.analysis(iopt=0,istp=i)
stress.append((dum[0] + dum[1])/2.)
ax.plot(strain_eff,stress,'-x',label='Model SF/IG/ehkl')
stress = []
for i in xrange(len(strain_eff)):
dum = myrs.analysis(iopt=1,istp=i)
stress.append((dum[0] + dum[1])/2.)
ax.plot(strain_eff,stress,'-x',label='Exp SF/IG/ehkl')
stress = []
for i in xrange(len(strain_eff)):
dum = myrs.analysis(iopt=2,istp=i)
stress.append((dum[0] + dum[1])/2.)
ax.plot(strain_eff,stress,'-x',label='Model SF/IG + Exp ehkl')
for i in xrange(len(exp_ref)):
exp_ref = exp_ref[i]
x,y=np.loadtxt(exp_ref).T
ax.plot(x,y,'--',label=exp_lab[i])
dum=np.loadtxt(mod_ref,skiprows=1).T;
x=dum[2]+dum[3];y=(dum[8]+dum[9])/2.;
ax.plot(x,y,'--',label='EVPSC biaxial')
__deco__(ax,iopt=3)
ax.legend(loc='best',fancybox=True).get_frame().set_alpha(0.5)
fig.tight_layout()
plt.show()
return stress
