def STIFF(FEMinput): #生成总体刚度矩阵 IJMH = FEMinput['IJMH'] E = FEMinput['TEVW']['E'] PR = FEMinput['TEVW']['pr'] t = FEMinput['TEVW']['t'] XY = FEMinput['XY'] KS = np.zeros((FEMinput['XY'].shape[0]*2, FEMinput['XY'].shape[0]*2)) XY = ref.arrIndexSort(XY) IJMH = ref.lstIndexSort(IJMH) KEs = map(lambda EN: STE(XY[EN[1][:,1],:], E, PR, t, EN[0], IJMH), IJMH)#单元刚度矩阵List KS = np.matrix(reduce(INCKE, KEs, KS)) return KS
def STIFF(FEMinput): #生成总体刚度矩阵 IJMH = FEMinput['IJMH'] E = FEMinput['TEVW']['E'] PR = FEMinput['TEVW']['pr'] t = FEMinput['TEVW']['t'] XY = FEMinput['XY'] KS = np.zeros((FEMinput['XY'].shape[0] * 2, FEMinput['XY'].shape[0] * 2)) XY = ref.arrIndexSort(XY) IJMH = ref.lstIndexSort(IJMH) KEs = map(lambda EN: STE(XY[EN[1][:, 1], :], E, PR, t, EN[0], IJMH), IJMH) #单元刚度矩阵List KS = np.matrix(reduce(INCKE, KEs, KS)) return KS
def GKS(FEMinput): #Calculate global stiffness matrix IJ = FEMinput['IJ'] E = FEMinput['TEV']['E'] V = FEMinput['TEV']['V'] T = FEMinput['TEV']['T'] XY = FEMinput['XY'] KS = np.zeros((XY.shape[0]*2, XY.shape[0]*2)) XY = ref.arrIndexSort(XY) IJ = ref.lstIndexSort(IJ) KE = KEsquare(E, V, T) KEsl = [KEdcr(KE, EN[0], IJ) for EN in IJ]#List of Element stiffness matrix KS = np.matrix(reduce(insKE, KEsl, KS)) return KS
def GKS(FEMinput): #Calculate global stiffness matrix IJ = FEMinput['IJ'] E = FEMinput['TEV']['E'] V = FEMinput['TEV']['V'] T = FEMinput['TEV']['T'] XY = FEMinput['XY'] KS = np.zeros((XY.shape[0] * 2, XY.shape[0] * 2)) XY = ref.arrIndexSort(XY) IJ = ref.lstIndexSort(IJ) KE = KEsquare(E, V, T) KEsl = [KEdcr(KE, EN[0], IJ) for EN in IJ] #List of Element stiffness matrix KS = np.matrix(reduce(insKE, KEsl, KS)) return KS