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
