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