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.
