def test2(r0=2.20,r45=1.95,r90=2.9):
"""
uniaxial tension tests
Arguments
---------
r0 =2.1
r90 =2.7
"""
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
from mk.library.lib import rot
import numpy as np
from mk.tests.mechtests import inplaneTension
h48 = wrapHill48R([r0,r45,r90])
psis = np.linspace(0,+np.pi/2.,100)
## uniaxial tension stress state in the laboratory axes
sUniaxial=np.zeros(6)
sUniaxial[0]=1.
print '%6s %6s %6s %6s %6s %6s %6s %6s'%(
'th','phi','s11','s22','s12','s1','s2','s3')
psis, rvs, phis = inplaneTension(fYLD=h48)
fig = plt.figure(figsize=(10,3))
ax1 = fig.add_subplot(131)
ax2 = fig.add_subplot(132)
ax3 = fig.add_subplot(133)
import vm_check
ax2.plot(psis*180/np.pi,rvs)
ax3.plot(psis*180/np.pi,phis)
vm_check.main(ax=ax1)
ax2.set_ylabel(r'R value')
ax3.set_ylabel(r'$\Phi$')
fig.tight_layout()
for ax in [ax2,ax3]:
ax.set_xlabel(r'$\theta$ [deg] from RD')
ax.set_xticks(np.arange(0,90.001,30))
fn = 'yf2_test2.pdf'
fig.savefig(fn,bbox_inches='tight')
print '%s has been saved'%fn
def test2(r0=2.20, r45=1.95, r90=2.9):
"""
uniaxial tension tests
Arguments
---------
r0 =2.1
r90 =2.7
"""
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
from mk.library.lib import rot
import numpy as np
from mk.tests.mechtests import inplaneTension
h48 = wrapHill48R([r0, r45, r90])
psis = np.linspace(0, +np.pi / 2., 100)
## uniaxial tension stress state in the laboratory axes
sUniaxial = np.zeros(6)
sUniaxial[0] = 1.
print '%6s %6s %6s %6s %6s %6s %6s %6s' % ('th', 'phi', 's11', 's22',
's12', 's1', 's2', 's3')
psis, rvs, phis = inplaneTension(fYLD=h48)
fig = plt.figure(figsize=(10, 3))
ax1 = fig.add_subplot(131)
ax2 = fig.add_subplot(132)
ax3 = fig.add_subplot(133)
import vm_check
ax2.plot(psis * 180 / np.pi, rvs)
ax3.plot(psis * 180 / np.pi, phis)
vm_check.main(ax=ax1)
ax2.set_ylabel(r'R value')
ax3.set_ylabel(r'$\Phi$')
fig.tight_layout()
for ax in [ax2, ax3]:
ax.set_xlabel(r'$\theta$ [deg] from RD')
ax.set_xticks(np.arange(0, 90.001, 30))
fn = 'yf2_test2.pdf'
fig.savefig(fn, bbox_inches='tight')
print '%s has been saved' % fn
def test3():
"""
"""
import vm_check
vm_check.main()
def test3():
"""
"""
import vm_check
vm_check.main()