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)
