Hello everyone.
I was able to apply the boundary conditions I mentioned earlier:
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.
I am also applying a displacement in y on the top face.
For this problem, both x and y displacement fields are expected to be
linear.
y displacement filed is as expected.
I expect my x displacement field to be a linear function of x, zero for
the left face.
But, somehow, my nodes on the right edge are also being fixed in x. So, x
displacement is zero everywhere.
I am unable to figure out why this is happening.
I am creating std::map boundary_values object to apply my boundary
conditions.
First, in setup_boundary_values(), I fill it with the boundary information
related to bottom 2 nodes.(Similar to step 79)
Next, I use the same boundary_values object to apply y displacement on the
top face as shown below:
VectorTools::interpolate_boundary_values(dof_handler,
1,
Functions::ConstantFunction<2>(u_x_values,2),
boundary_values,u_x_mask);
//imposing u_x=0 on
top face
VectorTools::interpolate_boundary_values(dof_handler,
1,
Functions::ConstantFunction<2>(u_y_values,2),
boundary_values,u_y_mask
);//imposing
u_y=0.05 on top face
Thank You
On Wednesday, October 12, 2022 at 9:59:35 PM UTC+5:30 Wasim Niyaz Munshi
ce21d400 wrote:
> 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/2fb027e7-801a-4962-b5d4-496c21be1e10n%40googlegroups.com.
