def ESIGME(DP, r, parts, E ,v , El, XY, IJ, T=0): #DP:np.array([[u, v],[],[],[]])#node0/1/2/3 step = 2.0/parts xis = np.arange(-1,1,step) etas = np.arange(-1,1,step) DP = DP.reshape((-1,1)) D = ref.D0(E,v) ES = [[ [ref.xietaToXY(El, XY, IJ, xi, eta), ((D*Bc(r, xi, eta))*DP).A1] for eta in etas] for xi in xis] return ES
def ESIGME(DP, r, scale, E ,v , El, XY, IJ, T): #DP:np.array([[w,thetax,thetay],[],[],[]])#node0/1/2/3 step = 2.0/scale xis = np.arange(-1,1,step) etas = np.arange(-1,1,step) DP = DP.reshape((-1,1)) D = ref.D0(E,v) ES = [[ [ref.xietaToXY(El, XY, IJ, xi, eta), ((D*Bc(r, xi, eta)*DP).A1)*(T/2.0)] for eta in etas] for xi in xis] return ES
def ESIGME(DP, r, parts, E, v, El, XY, IJ, T=0): #DP:np.array([[u, v],[],[],[]])#node0/1/2/3 step = 2.0 / parts xis = np.arange(-1, 1, step) etas = np.arange(-1, 1, step) DP = DP.reshape((-1, 1)) D = ref.D0(E, v) ES = [[[ ref.xietaToXY(El, XY, IJ, xi, eta), ((D * Bc(r, xi, eta)) * DP).A1 ] for eta in etas] for xi in xis] return ES