def _get_aramis_field_data(self):
t_fail = self.time_asc[-1]
ad = AramisFieldData(aramis_info=self.aramis_info)
current_step = (np.abs(ad.step_times - t_fail).argmin())
# print 'ad.step_times - t_fail', ad.step_times - t_fail
ad.current_step = int(current_step * 0.95)
return ad
def _get_aramis_field_data(self):
print '_get_aramis_field_data'
t_fail = self.t_asc[-1]
ad = AramisFieldData(aramis_info=self.aramis_info,
integ_radius=self.integ_radius,
integ_radius_crack=self.integ_radius_crack,
transform_data=True,
scale_data_factor=self.scale_data_factor)
current_step = (np.abs(ad.step_times - t_fail).argmin())
# print 'ad.step_times - t_fail', ad.step_times - t_fail
ad.current_step = current_step
return ad
def _get_aramis_field_data(self):
print '_get_aramis_field_data'
t_fail = self.t_asc[-1]
ad = AramisFieldData(aramis_info=self.aramis_info,
integ_radius=self.integ_radius,
integ_radius_crack=self.integ_radius_crack,
transform_data=True,
scale_data_factor=self.scale_data_factor)
current_step = (np.abs(ad.step_times - t_fail).argmin())
# print 'ad.step_times - t_fail', ad.step_times - t_fail
ad.current_step = current_step
return ad
p.ylabel('N [kN]')
p.legend(loc=2)
p.twinx()
p.plot(e.t_aramis_cut, e.M_t_aramis, color='black', label='M')
p.plot(e.t_aramis_cut, e.MF_t_aramis, color='red', label='M0')
p.plot(e.t_aramis_cut, e.MN_t_aramis, color='aqua', label='MII')
p.ylim(-0.25, 0.5)
# p.xlim(0, 950)
p.ylabel('M [kNm]')
p.legend(ncol=3)
p.subplot(224)
for e in e_list:
aramis_file_path = e.get_cached_aramis_file('Xf15s3-Yf15s3')
AI = AramisInfo(data_dir=aramis_file_path)
ad = AramisFieldData(aramis_info=AI, integ_radius=3)
max_step = e.n_steps
a = e.crack_bridge_strain_all
n_fa = ad.d_ux.shape[0]
h = np.linspace(e.pos_fa[0], e.pos_fa[1], num=n_fa)
# print 'h', h
for step in range(0, max_step, 10):
ad.current_step = step
if a == None:
mid_idx = ad.d_ux.shape[1] / 2
eps_range = 3
eps = np.mean(ad.d_ux[:, mid_idx - eps_range:mid_idx + eps_range], axis=1)
else:
# p.tick_params
p.plot(e.eps_N * 1000, e.N_t_N, color='black', label='cbs (N)')
p.ylim(0, 30)
p.grid()
p.xlabel('crack bridge strain [1*E-3]')
p.ylabel('N [kN]')
p.legend(loc=1)
p.title(test_files)
# eps(N) for aramis steps
p.subplot(232)
for e in e_list:
aramis_file_path = e.get_cached_aramis_file('Xf15s3-Yf15s3')
AI = AramisInfo(data_dir=aramis_file_path)
ad = AramisFieldData(aramis_info=AI,
integ_radius=3)
absa = AramisCDT(aramis_info=AI,
aramis_data=ad)
max_step = e.n_steps
a = e.crack_bridge_strain_all
n_fa = ad.d_ux.shape[0]
h = np.linspace(e.pos_fa[0], e.pos_fa[1], num=n_fa)
t_N = e.t_N_arr
if len(t_N) == 0:
N_end_idx = 0
else:
N_end_idx = np.shape(t_N) [0]
for step in range(0, N_end_idx, 5):