二维RVE等效性能求解

"""
注意数组与矩阵的不同
array（[[1,2,3],[1,1,2],[1,2,2]])
a[1][2]访问的是第二行第三个元素：2
b=mat（a）将数组转换为矩阵
b[1]访问第2个向量,但不能通过b[1][2]访问向量中的元素
"""
# -*- coding:UTF-8 -*-
from numpy import *         #导入numpy的库函数
import numpy as np          #np.
C_matrix = mat([[3179.71, 1686.34, 0], [1686.34, 3153.74,0], [0, 0, 738.74]])
S_matrix = C_matrix.I
S_array = array(S_matrix)
E11 = 1/S_array[0][0]
E22 = 1/S_array[1][1]
G12 = 1/S_array[2][2]
V21 = -E22 * S_array[0][1]
print E11, E22, G12, V21, 1-E22/E11