Hello everybody,
i am trying to change my timestep after few simulations steps but i am
getting a kink into my solution once i reduce my timestep.
I did this actually simple by adjusting my timestep once a certain time is
reached when i apply the dirichlet boundary conditions:
if(time.current() > trans*first_loading)
{
*time.adjust_timestep_size(0);*
const double delta_length = parameters.stretch1;
const unsigned int n_stretch_steps = static_cast<unsigned int>(
parameters.load_time_first/time.get_delta_t());
const double delta_u_x = delta_length/n_stretch_steps;
VectorTools::interpolate_boundary_values(dof_handler_ref,
boundary_id,
ConstantFunction<dim>(-
delta_u_x,n_components),
constraints,
fe.component_mask(
x_displacement));
// Fix z displacement to avoid rotation
VectorTools::interpolate_boundary_values(dof_handler_ref,
boundary_id,
ZeroFunction<dim>(
n_components),
constraints,
fe.component_mask(
z_displacement));
}
The adjust_timestep_size() void just does change delta_t to a new delta_t
but i cannot find my error.
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/193e058e-4f71-4b33-8749-a1633ce40f98n%40googlegroups.com.
