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
