Dear Tuanny, There is no way to do this when reading in from a mesh file. You can, however, add this constraint manually. See this post https://groups.google.com/forum/#!searchin/dealii/vertex_dof_index%7Csort:relevance/dealii/8fUwabkznaU/VlTERBo-GQAJ and other post that reference "vertex_dof_index" for details on how to do it.
I hope that this helps you solve your problem! Regards, Jean-Paul On Friday, August 18, 2017 at 10:05:36 AM UTC+2, Tuanny Cajuhi wrote: > > Hi Wolfgang, > > thank you for the answers! > > >>> Currently I am just using >>> >>> set_all_boundary_ids and set_manifold_id to identify the nodes around >>> the hole. Is there a way to specify the equivalent node (the node in the >>> middle related to the constrained nodes) in each hole and apply the >>> boundary conditions directly in this node? >>> >> >> That's the wrong way to see it. Mathematically, you're applying loads to >> faces, not individual nodes. You should set the boundary_id of the faces on >> which you're trying to apply a load to a different value than the other >> faces. You can then use >> >> VectorTools::interpolate_boundary_value with that particular boundary_id >> to evaluate the load you want to apply. >> > > Thank you for the correction, I expressed myself wrong but I am applying > it correctly :) > >> >> >> Is it also possible to read the equivalent node from the *inp file? >>> >> >> I think I don't quite know what you mean by the "equivalent node". Do you >> mean the center point of the hole? That's not a part of the mesh, so it's >> not usually mentioned in the .inp file. >> > > I think the > > expression "equivalent node" could be misleading from my side. In Abaqus, > there is an option to set a reference point (for example the center of the > hole) and couple it to a surface or node region, e.g. surfaces of the hole. > The boundary conditions and load can be applied to this reference point. In > the .inp file, something like: > > *Nset, nset=_PickedSet35, internal > 1, > *Elset, elset=__PickedSurf34_S1, internal, instance=Part-1-1 > 68, 81, 94, 101 > *Elset, elset=__PickedSurf34_S2, internal, instance=Part-1-1 > 77, 80, 92, 93, 110 > *Surface, type=ELEMENT, name=_PickedSurf34, internal > __PickedSurf34_S1, S1 > __PickedSurf34_S2, S2 > ** Constraint: constraint_upper_hole > *Coupling, constraint name=constraint_upper_hole, ref node=_PickedSet35, > surface=_PickedSurf34 > *Kinematic > > is generated. I was wondering if there was something similar in deal.II > where I could select the surfaces/edges of the hole and then > couple/constrain the boundaries to one reference point. From that, I'd > simply call the VectorTools::interpolate_boundary_value and apply the > boundary conditions to that reference. Is that possible, unnecessary or > there is a better way to do that? > > Thank you. > > Best, > Tuanny > > -- 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]. For more options, visit https://groups.google.com/d/optout.
