def tab1_cpb(iopt): """ Mg-Th coefficients k c12 c13 c22 c23 c33 1% 0.3539 0.4802 0.2592 0.9517 0.2071 0.4654 5% 0.2763 0.3750 0.0858 0.9894 0.0659 0.1238 10% 0.0598 0.6336 0.2332 1.4018 0.5614 0.7484 iopt=0: returns 1% iopt=1: returns 5% iopt=2: returns 10% """ k = 0. c11 = 0. c12 = 0. c13 = 0. c22 = 0. c23 = 0. c33 = 0. c44 = 0. c55 = 0. c66 = 0. if iopt==0: k,c12,c13,c22,c23,c33=0.3539,0.4802,0.2592,0.9517,0.2071,0.4654 elif iopt==1: k,c12,c13,c22,c23,c33=0.2763,0.3750,0.0858,0.9894,0.0659,0.1238 elif iopt==2: k,c12,c13,c22,c23,c33=0.0598,0.6336,0.2332,1.4018,0.5614,0.7484 else: raise IOError, 'Unexpected case.' c = cpb_lib.calcC(c12=c12,c13=c13,c22=c22,c23=c23,c33=c33) return c, k
def tab1_cpb(iopt): """ Mg-Th coefficients k c12 c13 c22 c23 c33 1% 0.3539 0.4802 0.2592 0.9517 0.2071 0.4654 5% 0.2763 0.3750 0.0858 0.9894 0.0659 0.1238 10% 0.0598 0.6336 0.2332 1.4018 0.5614 0.7484 iopt=0: returns 1% iopt=1: returns 5% iopt=2: returns 10% """ k = 0. c11 = 0. c12 = 0. c13 = 0. c22 = 0. c23 = 0. c33 = 0. c44 = 0. c55 = 0. c66 = 0. if iopt == 0: k, c12, c13, c22, c23, c33 = 0.3539, 0.4802, 0.2592, 0.9517, 0.2071, 0.4654 elif iopt == 1: k, c12, c13, c22, c23, c33 = 0.2763, 0.3750, 0.0858, 0.9894, 0.0659, 0.1238 elif iopt == 2: k, c12, c13, c22, c23, c33 = 0.0598, 0.6336, 0.2332, 1.4018, 0.5614, 0.7484 else: raise IOError, 'Unexpected case.' c = cpb_lib.calcC(c12=c12, c13=c13, c22=c22, c23=c23, c33=c33) return c, k
def tab2_cpb(iopt): """ Mg-Li coefficients k c12 c13 c22 c23 c33 1% 0.2026 0.5871 0.6975 0.9783 0.2840 0.1497 5% 0.2982 0.6103 0.8056 1.0940 0.5745 0.1764 10% 0.1763 0.5324 0.8602 1.0437 0.8404 0.2946 iopt=0: returns 1% iopt=1: returns 5% iopt=2: returns 10% """ k = 0. c11 = 0. c12 = 0. c13 = 0. c22 = 0. c23 = 0. c33 = 0. c44 = 0. c55 = 0. c66 = 0. if iopt==0: k,c12,c13,c22,c23,c33=0.2026,0.5871,0.6975,0.9783,0.2840,0.1497 elif iopt==1: k,c12,c13,c22,c23,c33=0.2982,0.6103,0.8056,1.0940,0.5745,0.1764 elif iopt==2: k,c12,c13,c22,c23,c33=0.1763,0.5324,0.8602,1.0437,0.8404,0.2946 else: raise IOError, 'Unexpected case.' c = cpb_lib.calcC(c12=c12,c13=c13,c22=c22,c23=c23,c33=c33) return c, k
def tab2_cpb(iopt): """ Mg-Li coefficients k c12 c13 c22 c23 c33 1% 0.2026 0.5871 0.6975 0.9783 0.2840 0.1497 5% 0.2982 0.6103 0.8056 1.0940 0.5745 0.1764 10% 0.1763 0.5324 0.8602 1.0437 0.8404 0.2946 iopt=0: returns 1% iopt=1: returns 5% iopt=2: returns 10% """ k = 0. c11 = 0. c12 = 0. c13 = 0. c22 = 0. c23 = 0. c33 = 0. c44 = 0. c55 = 0. c66 = 0. if iopt == 0: k, c12, c13, c22, c23, c33 = 0.2026, 0.5871, 0.6975, 0.9783, 0.2840, 0.1497 elif iopt == 1: k, c12, c13, c22, c23, c33 = 0.2982, 0.6103, 0.8056, 1.0940, 0.5745, 0.1764 elif iopt == 2: k, c12, c13, c22, c23, c33 = 0.1763, 0.5324, 0.8602, 1.0437, 0.8404, 0.2946 else: raise IOError, 'Unexpected case.' c = cpb_lib.calcC(c12=c12, c13=c13, c22=c22, c23=c23, c33=c33) return c, k