[deal.II] Re: increment Dirichlet boundary condition

2018-02-23 Thread Bryukhanov Ilya
I tried to use code from step-22 and from the sources you recommended. The problem is that the solution is zero over the body except the boundary where Dirichlet conditions were set. It is just one iteration. I don't know why they are not in use. Briefly, my code consists of the following steps

[deal.II] Re: increment Dirichlet boundary condition

2018-02-23 Thread Bryukhanov Ilya
I tried to use code from step-22 and from the sources you recommended. The problem is that the solution is zero over the body except the boundary where Dirichlet conditions were set. It is just one iteration. I don't know why they are not in use. Briefly, my code consists of the following steps

[deal.II] Re: increment Dirichlet boundary condition

2018-02-22 Thread Denis Davydov
Hi IIiya, On Thursday, February 22, 2018 at 11:55:33 AM UTC+1, Bryukhanov Ilya wrote: > > Denis, thanks a lot for the answer! > > I think that on the first step I can use usual > "interpolate_boundary_values + apply_boundary_values" functions, > the outline I gave above is agnostic to "first

[deal.II] Re: increment Dirichlet boundary condition

2018-02-22 Thread Bryukhanov Ilya
Denis, thanks a lot for the answer! I think that on the first step I can use usual "interpolate_boundary_values + apply_boundary_values" functions, that transform the global matrix: columns of the constrained nodes become zero columns with diagonal non-zero entries. In a time loop I have to

[deal.II] Re: increment Dirichlet boundary condition

2018-02-22 Thread Bryukhanov Ilya
Denis, thanks a lot for the answer! I think that on the first step I can use usual "interpolate_boundary_values + apply_boundary_values" functions, that transform the global matrix: columns of the constrained nodes become zero columns with diagonal non-zero entries. In a time loop I have to

[deal.II] Re: increment Dirichlet boundary condition

2018-02-22 Thread Bryukhanov Ilya
Denis, thanks a lot for the answer! I think that on the first step I can use usual "interpolate_boundary_values + apply_boundary_values" functions, and this procedures transform the global matrix so that the columns of constrained nodes become zero columns with diagonal non-zero entries. In a

[deal.II] Re: increment Dirichlet boundary condition

2018-02-21 Thread Denis Davydov
Hi Iliya, there are numerous things that are wrong On Wednesday, February 21, 2018 at 5:32:55 PM UTC+1, Bryukhanov Ilya wrote: > > Denis, thanks a lot for your answer and links that you provided. > > However, my solution becomes in a way that the displacement values on that > boundary

[deal.II] Re: increment Dirichlet boundary condition

2018-02-21 Thread Bryukhanov Ilya
Denis, thanks a lot for your answer and links that you provided. However, my solution becomes in a way that the displacement values on that boundary abnormally grow even without assigning new boundary condition in a time loop. I attach figures on the zero step (correct) and on the first step

[deal.II] Re: increment Dirichlet boundary condition

2018-02-20 Thread Denis Davydov
Hi Ilya, I am not away of principally simple solution, the RHS simply contains your bilinear form times interpolated values at constrained DoFs. What you can do is to loop over a subset of cells (those which are at the boundary), still assemble the local matrix and do