21 Apr
2017
21 Apr
'17
11:10 a.m.
After looking deep into the example I created my version of minimum workng example of the script script. Here it is:
filename_in = "linear_viscoelastic.py"
conf = ProblemConf.from_file(filename_in)
pb = Problem.from_conf(conf)
time_solver = pb.get_time_solver()
time_solver.init_time()
for out in time_solver():
step, time, state = out
print ("##################")
print("step : %s time : %s " % (step, time))
u = state()
x = pb.get_mesh_coors().T[0]
y = pb.get_mesh_coors().T[1]
z = pb.get_mesh_coors().T[2]
Vx = state.reshape(state.shape[0]/3,3).T[0]
Vy = state.reshape(state.shape[0]/3,3).T[1]
Vz = state.reshape(state.shape[0]/3,3).T[2]
Vx Vy Vz are different at each new step (at least at first few steps). But coordinates of the mesh stays the same. How can I deform the mesh according to Vx Vy Vz? Or it is the wrong way of getting coordinates? ~