def _create_kei_func_grid(shape, spacing, alpha):
from scipy.special import kei
dx, dy = np.meshgrid(np.arange(shape[1]) * spacing[1],
np.arange(shape[0]) * spacing[0])
return kei(np.sqrt(dx ** 2 + dy ** 2) / alpha)
def _set_kei_func_grid(self):
from scipy.special import kei
dx, dy = np.meshgrid(
np.arange(self._grid.number_of_node_columns) * self._grid.dx,
np.arange(self._grid.number_of_node_rows) * self._grid.dy)
return kei(np.sqrt(dx ** 2 + dy ** 2) / self.alpha)
def _create_kei_func_grid(shape, xy_spacing, alpha):
from scipy.special import kei
dx, dy = np.meshgrid(
np.arange(shape[1]) * xy_spacing[1], np.arange(shape[0]) * xy_spacing[0]
)
return kei(np.sqrt(dx ** 2 + dy ** 2) / alpha)
def _set_kei_func_grid(self):
from scipy.special import kei
alpha = get_flexure_parameter(self._eet, self._youngs, 2)
dx, dy = np.meshgrid(
np.arange(self._grid.number_of_node_columns) * self._grid.dx,
np.arange(self._grid.number_of_node_rows) * self._grid.dy)
return kei(np.sqrt(dx ** 2 + dy ** 2) / alpha)
def uniform_loading(r, x, a, rou_p, l):
if x <= r:
kerpa = sp.kerp(a)
berx = sp.ber(x / l)
keipa = sp.keip(a)
beix = sp.beip(x / l)
w = rou_basalt * uniform_h * (a * kerpa * berx - a * keipa * beix +
1) / rou_p
else:
berpa = sp.berp(a)
kerx = sp.ker(x / l)
beipa = sp.beip(a)
keix = sp.kei(x / l)
w = rou_basalt * uniform_h * (a * berpa * kerx -
a * beipa * keix) / rou_p
return w
def get_fcw(_fcw):
ucw_star_mag = sqrt(1 / 2 * _fcw * self.V1) * ub_mag #eqn 50
l = kappa * ucw_star_mag / omega #eqn 29
zeta_0 = kb / 30 / l #eqn 31
ztemp = 2 * sqrt(zeta_0) #eqn 55
K = 1 / ztemp * 1 / sqrt(
square(ker(ztemp)) + square(kei(ztemp))) #eqn 55
temp_54_a = 0.0971 * sqrt(kb / self.Ab) * K / power(
_fcw, 3.0 / 4.0) #eqn 54
# temp_54_b = self.V2/2.0/power(self.alpha,1.0/4.0) #eqn 54
temp_54_b = self.V2 * power(self.V1, 1.0 / 4.0) / 2.0 / sqrt(
self.alpha) #eqn 54 THESIS!
rhs_54 = square(temp_54_a) + 2.0 * temp_54_a * temp_54_b * cos(
self.phi_bar) #eqn 54
# lhs_54 = power(self.alpha, 3.0/4.0)/4.0 - square(temp_54_b) #eqn 54
lhs_54 = self.alpha * sqrt(self.V1) / 4.0 - square(
temp_54_b) #eqn 54
self.uw_star_mag = sqrt(kappa * ucw_star_mag * ub_mag *
sqrt(zeta_0) * K)
return lhs_54 - rhs_54
