Thank you Prof. Bangerth. I was indeed calling the setup_boundary_values() function before call dof_handler.distribute_dofs(). I made the necessary adjustments. It is working fine now
Thanks Wasim Niyaz Research scholar CE Dept. IITM On Tue, 11 Oct, 2022, 7:03 pm Wolfgang Bangerth, <[email protected]> wrote: > > Wasim, > I imagine you are calling your setup_boundary_values() function before you > call dof_handler.distribute_dofs()? > Best > W. > > On 10/11/22 05:16, Wasim Niyaz Munshi ce21d400 wrote: > > *** Caution: EXTERNAL Sender *** > > > > I was using a constraint object to apply the boundary conditions and > component > > mask to deal with hanging nodes and apply x and y displacements > separately. > > Now, I want to solve an elasticity problem on a unit square wherein I > want to > > apply a displacement in y direction on top face and have a pin and > roller > > support at lower left and right node respectively. > > I tried to do this as follows: > > 1. I used the constraint object to apply displacement in y direction on > top face. > > 2. I created a std::map boundary_values object to enforce a pin and > roller > > support bc at lower left and right node respectively. > > I took a look at step79. They seem to have similar dirichlet bc at the > bottom > > but apply a neumann bc at the top. I followed a similar procedure to > have a > > pin and roller bc at lower left and right node respectively(using the > > setup_boundary_values class as they have done). > > > > > > My code throws the following error: > > An error occurred in line <324> of file > > > </home/wasim/dealii-candi/deal.II-v9.4.0/include/deal.II/dofs/dof_accessor.templates.h> > > > in function > > static void > > > dealii::internal::DoFAccessorImplementation::Implementation::process_dof_index(const > > > dealii::DoFHandler<dim, spacedim>&, unsigned int, unsigned int, unsigned > int, > > unsigned int, const std::integral_constant<int, structdim>&, > GlobalIndexType&, > > const DoFPProcessor&) [with int dim = 2; int spacedim = 2; int structdim > = 0; > > GlobalIndexType = unsigned int; DoFPProcessor = > > > dealii::internal::DoFAccessorImplementation::Implementation::get_dof_index<2, > > > 2, 0>::<lambda(const auto:9&, auto:10&)>] > > The violated condition was: > > > ::dealii::deal_II_exceptions::internals::compare_less_than(obj_level, > > dof_handler.object_dof_indices.size()) > > Additional information: > > Index 0 is not in the half-open range [0,0). In the current case, > this > > half-open range is in fact empty, suggesting that you are accessing > an > > element of an empty collection such as a vector that has not been > set > > to the correct size. > > > > Stacktrace: > > ----------- > > #0 ./step-103: void > > > dealii::internal::DoFAccessorImplementation::Implementation::process_dof_index<2, > > > 2, 0, unsigned int, > > > dealii::internal::DoFAccessorImplementation::Implementation::get_dof_index<2, > > > 2, 0>(dealii::DoFHandler<2, 2> const&, unsigned int, unsigned int, > unsigned > > int, unsigned int, std::integral_constant<int, 0> > const&)::{lambda(auto:1 > > const&, auto:2&)#1}>(dealii::DoFHandler<2, 2> const&, unsigned int, > unsigned > > int, unsigned int, unsigned int, std::integral_constant<int, 0> const&, > > unsigned int&, > > > dealii::internal::DoFAccessorImplementation::Implementation::get_dof_index<2, > > > 2, 0>(dealii::DoFHandler<2, 2> const&, unsigned int, unsigned int, > unsigned > > int, unsigned int, std::integral_constant<int, 0> > const&)::{lambda(auto:1 > > const&, auto:2&)#1} const&) > > #1 ./step-103: unsigned int > > > dealii::internal::DoFAccessorImplementation::Implementation::get_dof_index<2, > > > 2, 0>(dealii::DoFHandler<2, 2> const&, unsigned int, unsigned int, > unsigned > > int, unsigned int, std::integral_constant<int, 0> const&) > > #2 ./step-103: dealii::DoFAccessor<2, 2, 2, > false>::vertex_dof_index(unsigned > > int, unsigned int, unsigned int) const > > #3 ./step-103: step103::ElasticProblem::setup_boundary_values() > > #4 ./step-103: step103::ElasticProblem::run() > > #5 ./step-103: main > > -------------------------------------------------------- > > > > make[3]: *** [CMakeFiles/run.dir/build.make:71: CMakeFiles/run] Aborted > (core > > dumped) > > make[2]: *** [CMakeFiles/Makefile2:116: CMakeFiles/run.dir/all] Error 2 > > make[1]: *** [CMakeFiles/Makefile2:123: CMakeFiles/run.dir/rule] Error 2 > > make: *** [Makefile:137: run] Error 2 > > Is it because I am applying Dirichlet bc using two objects, > constraints and > > boundary_values? > > > > > > Thank you > > > > > > > > On Mon, Oct 10, 2022 at 7:28 PM Wolfgang Bangerth < > [email protected] > > <mailto:[email protected]>> wrote: > > > > On 10/10/22 07:42, Wasim Niyaz Munshi ce21d400 wrote: > > > I am trying to solve the elasticity equation on a square domain > > (similar to > > > step 8). > > > For the patch test I am using just one element. > > > I am having problems with the application of boundary conditions. > > > I want to fix the lower left node in both x and y, and also fix > the lower > > > right node in y. > > > The lower right node is free to move in x. > > > The problem I am facing is that both these nodes are on the same > face, > > and I > > > do not know how to specify different values of x displacement for > the 2 > > nodes > > > of the same face. > > > > > > Is there a way to set boundary id to a node rather than to a face? > > > > No -- principally because from a mathematical perspective, you can > only set > > boundary indicators on parts of the boundary that have (d-1) > dimensional > > measure larger than zero. Vertices do not satisfy this. > > > > But it is a common thing to do anyway. You might want to look at the > step-79 > > that does something similar. > > > > Best > > W> > > > > > > -- > > > ------------------------------------------------------------------------ > > Wolfgang Bangerth email: [email protected] > > <mailto:[email protected]> > > www: > http://www.math.colostate.edu/~bangerth/ > > <http://www.math.colostate.edu/~bangerth/> > > > > -- > > The deal.II project is located at http://www.dealii.org/ > > < > https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.dealii.org%2F&data=05%7C01%7CWolfgang.Bangerth%40colostate.edu%7C861d8e11d32747990c3308daab7a1301%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C638010838048599291%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=GczFY9xgCuNhKCygAdQqZMzE3pQtDVR%2BJunoKQQyndQ%3D&reserved=0 > > > > For mailing list/forum options, see > > https://groups.google.com/d/forum/dealii?hl=en > > < > https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgroups.google.com%2Fd%2Fforum%2Fdealii%3Fhl%3Den&data=05%7C01%7CWolfgang.Bangerth%40colostate.edu%7C861d8e11d32747990c3308daab7a1301%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C638010838048755540%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=Al7MfSTwuuRHdmniebbLxp4ugBaX1ddbBdlldSxuX5o%3D&reserved=0 > > > > --- > > You received this message because you are subscribed to a topic in > the > > Google Groups "deal.II User Group" group. > > To unsubscribe from this topic, visit > > https://groups.google.com/d/topic/dealii/jI31ro8W_Cg/unsubscribe > > < > https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgroups.google.com%2Fd%2Ftopic%2Fdealii%2FjI31ro8W_Cg%2Funsubscribe&data=05%7C01%7CWolfgang.Bangerth%40colostate.edu%7C861d8e11d32747990c3308daab7a1301%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C638010838048755540%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=u7Mxyrfrhc1DwitmIF%2BiRyqv2syo1cseWqCO%2Bsvnoyc%3D&reserved=0 > >. > > To unsubscribe from this group and all its topics, send an email to > > [email protected] > > <mailto:dealii%[email protected]>. > > To view this discussion on the web visit > > > https://groups.google.com/d/msgid/dealii/64772c85-044f-3492-5895-b1d6ac398b0f%40colostate.edu > > < > https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgroups.google.com%2Fd%2Fmsgid%2Fdealii%2F64772c85-044f-3492-5895-b1d6ac398b0f%2540colostate.edu&data=05%7C01%7CWolfgang.Bangerth%40colostate.edu%7C861d8e11d32747990c3308daab7a1301%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C638010838048755540%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=gE%2B7eZCtkI2CnOoTmXDBuSrLav%2FLa1PmFvtYv2EIDW8%3D&reserved=0 > >. > > > > -- > > The deal.II project is located at http://www.dealii.org/ > > < > https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.dealii.org%2F&data=05%7C01%7CWolfgang.Bangerth%40colostate.edu%7C861d8e11d32747990c3308daab7a1301%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C638010838048755540%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=eJu3T7DEwblMZTdI3vLLnxO2RfznkOCeZEj8k2rSKPY%3D&reserved=0 > > > > For mailing list/forum options, see > > https://groups.google.com/d/forum/dealii?hl=en > > < > https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgroups.google.com%2Fd%2Fforum%2Fdealii%3Fhl%3Den&data=05%7C01%7CWolfgang.Bangerth%40colostate.edu%7C861d8e11d32747990c3308daab7a1301%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C638010838048755540%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=Al7MfSTwuuRHdmniebbLxp4ugBaX1ddbBdlldSxuX5o%3D&reserved=0 > > > > --- > > 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] > > <mailto:[email protected]>. > > To view this discussion on the web visit > > > https://groups.google.com/d/msgid/dealii/CAM8ps5DTY_KXs7MZq%2B_zaW0aB2JemWgJ4SLf3dEow%3D_eqDYutA%40mail.gmail.com > > < > https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgroups.google.com%2Fd%2Fmsgid%2Fdealii%2FCAM8ps5DTY_KXs7MZq%252B_zaW0aB2JemWgJ4SLf3dEow%253D_eqDYutA%2540mail.gmail.com%3Futm_medium%3Demail%26utm_source%3Dfooter&data=05%7C01%7CWolfgang.Bangerth%40colostate.edu%7C861d8e11d32747990c3308daab7a1301%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C638010838048755540%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=WOjaRCuN9KlVtE1fGxNdAwQ3U0z4x%2FWG4oQ4kGrjWBQ%3D&reserved=0 > >. > > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: [email protected] > www: http://www.math.colostate.edu/~bangerth/ > > -- > 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 a topic in the > Google Groups "deal.II User Group" group. > To unsubscribe from this topic, visit > https://groups.google.com/d/topic/dealii/jI31ro8W_Cg/unsubscribe. > To unsubscribe from this group and all its topics, send an email to > [email protected]. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/b1a9fa0b-2505-c786-c24b-e98d63ea4aa0%40colostate.edu > . > -- 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/CAM8ps5AsSwvFoghHoy1bNJOgnQ8JqDaH7SVbDYx7c59whQe69w%40mail.gmail.com.
