Beispiel #1
0
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
Beispiel #2
0
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
Beispiel #3
0
def test3():
    """
    """
    import vm_check
    vm_check.main()
Beispiel #4
0
def test3():
    """
    """
    import vm_check
    vm_check.main()