[deal.II] Re: geometry

Daniel, Thanks very much for your help. For improving my knowledge and information, It would be very kind of you if you let me know why my code does not produce right geometry and which part is wrong. Sincerely, Benhour On Sunday, July 17, 2016 at 7:10:41 AM UTC-5, Daniel Arndt wrote: > > Benho

[deal.II] Re: geometry

Benhour, you can find a modification of step-3 in which the domain is a quarter circle attached. The four cells from GridGenerator::half_hyper_ball have been replaced by just three. The vertex (0,-radius) has been removed and all the remaining vertices with a negative y-component have been move

[deal.II] Re: geometry

Daniel, In that code, I had a half circle in the right part of the coordinate axis. By writing v(1) <=0, I wanted to choose the part with negative y components from the geometry, however, The attached file has been obtained. I really do not know how I should modify the source you sent to me beca

[deal.II] Re: geometry

Benhour, It would be helpful if you tell us what the problem with the newly created problem is. Looking at the implementation of GridGenerator::half_hyper_ball [1], you can observe that there are cells which have vertices both with positive and negative y-components. This will likely lead to und

[deal.II] Re: geometry

Dear J-P, Thanks for your response. I used the code for creating a quarter of a circle that comes as follow: Triangulation<2> triangulation; const Point<2> center; const double radius = 1.; GridGenerator::half_hyper_ball(triangulation, center, radius); Triangulation<2>::active_cell_iterator cell

[deal.II] Re: geometry

Dear Benhour, You could use one of the GridGenerator options to create a half circle, and then use GridGenerator::create_triangulation_with_removed_cells to remove the excess cells that

[deal.II] Re: geometry

Benhour, Did you try what I suggested in https://groups.google.com/d/msg/dealii/QPLvZIQx7yw/LQ7t0-xWBgAJ ? Best, Daniel -- 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

[deal.II] Re: Geometry and boundary conditions

Hi Benhour, Have you looked at any of the tutorials on geometry creation and solid mechanics? I believe that these specific points are covered there. Regards, J-P On Thursday, July 7, 2016 at 5:21:30 PM UTC+2, benhour.amiria...@gmail.com wrote: > > Dear Daniel, > Thanks very much for your resp

[deal.II] Re: Geometry and boundary conditions

Dear Daniel, Thanks very much for your response. I should model a whole circle, however for simplicity I want to model one quarter of it. In fact I have 2 quarter of circles with different radius and with the same center. I have a boundary load on the perimeter(curve) side of the circle, axial s

[deal.II] Re: Geometry and boundary conditions

Benhour, Have a look at the implementation of GridGenerator::half_hyper_ball [1]. You probably just want to use the first 4 vertices and then p, p+Point<2>(0,-1) *(radius/std::sqrt(2.0)*a), p+Point<2>(0,-1) *radius, additionally. Can you specify in formulas what symmetry boundary conditions mea