Dear Shawn, No, its fine. I see that the changes weren’t as drastic as I thought they’d be, so I think that the boundary ID assignment will work out of the box (I’d forgotten that someone else had already done much of the work necessary to support codim-1).
I’ve made a patch to add this functionality. See here: https://github.com/dealii/dealii/pull/7217 <https://github.com/dealii/dealii/pull/7217> Thank a lot for the grids and the minimal test case! Best, J-P > On 20 Sep 2018, at 08:11, Yuxiang Wang <[email protected]> wrote: > > Oh I see! Thank you for explaining. My guess would be that the ABAQUS > equivalent is "node set". Please let me know if you'd like me to try that! > > Best, > Shawn > > On Wed, Sep 19, 2018 at 10:00 PM Jean-Paul Pelteret <[email protected] > <mailto:[email protected]>> wrote: >> Hi Shawn, > > >> Thank you for your help! I'd love to provide the mesh. Please let me know >> whether the attached one would work. > > You’re welcome. This second mesh is perfect. I’ll let you know when I’ve > submitted a patch that allows this grid to be read in. > >> For the boundary coloring - could you please elaborate on what that means? I >> wasn't aware that ABAQUS has such a coloring feature but it could be me not >> knowing this though. > > This means assigning some boundary ID’s (in all of the documentation we call > it “colouring”). Cubit allows one to do this, but I don’t know about ABAQUS > itself (I would have thought it would be possible). I’ll try to modify the > input file you supplied, which should be possible since the geometry is > simple enough. > > Best, > Jean-Paul > >> >> On Tue, Sep 18, 2018 at 10:57 PM Jean-Paul Pelteret <[email protected] >> <mailto:[email protected]>> wrote: >> Dear Shawn, >> >> The GridIn::read_abaqus() function is currently not equipped to deal with >> the codimension 1 case. I might be able to find some time to look at this in >> the next few days though. In addition to the mesh file that you’ve already >> given us, may you please construct and give us a fully structured grid in >> codimension 1 (so something like a rectangle warped in the third dimension) >> that also has some boundary colouring? That would help a great deal with >> implementing this feature, and would also ensure that the assignment of >> boundary ID’s is possible. >> >> Best, >> Jean-Paul >> >>> On 19 Sep 2018, at 01:38, Yuxiang Wang <[email protected] >>> <mailto:[email protected]>> wrote: >>> >>> Prof. Bangerth, >>> >>> Thank you for the prompt response! Yes please find attached the minimal >>> program. The main body is pretty short (see below), modified from step-49. >>> Please let me know if you cannot reproduce it! >>> >>> ``` >>> int main () >>> { >>> Triangulation<2, 3> triangulation; >>> >>> GridIn<2, 3> gridin; >>> gridin.attach_triangulation(triangulation); >>> std::ifstream f("mesh.inp"); >>> gridin.read_abaqus(f); >>> >>> GridOut gridout; >>> std::ofstream fout("mesh.vtu"); >>> gridout.write_vtu(triangulation, fout); >>> } >>> >>> ``` >>> >>> >>> Shawn >>> >>> On Tue, Sep 18, 2018 at 4:28 PM Wolfgang Bangerth <[email protected] >>> <mailto:[email protected]>> wrote: >>> On 09/18/2018 03:18 PM, Yuxiang Wang wrote: >>> > >>> > Thank you for the pointer! I changed my code as you suggested and it >>> > finished >>> > running. >>> >>> Great! >>> >>> >>> > However, when I opened the output file (mesh.vtu), the 3D mesh was >>> > projected >>> > to a 2D plane of z=0, and lost its 3D features. Was that expected >>> > behavior? >>> >>> No. Are you saying that if you read a mesh, then immediately output it >>> again, >>> that the z-coordinate is lost? If that's correct, then that's a bug and it >>> would be nice to have a small program and input file that demonstrates this. >>> >>> 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/ >>> <http://www.dealii.org/> >>> For mailing list/forum options, see >>> https://groups.google.com/d/forum/dealii?hl=en >>> <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/Erk7pxwTEj4/unsubscribe >>> <https://groups.google.com/d/topic/dealii/Erk7pxwTEj4/unsubscribe>. >>> To unsubscribe from this group and all its topics, send an email to >>> [email protected] >>> <mailto:dealii%[email protected]>. >>> For more options, visit https://groups.google.com/d/optout >>> <https://groups.google.com/d/optout>. >>> >>> >>> -- >>> Yuxiang "Shawn" Wang, PhD >>> [email protected] <mailto:[email protected]> >>> +1 (434) 284-0836 >>> >>> -- >>> The deal.II project is located at http://www.dealii.org/ >>> <http://www.dealii.org/> >>> For mailing list/forum options, see >>> https://groups.google.com/d/forum/dealii?hl=en >>> <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] >>> <mailto:[email protected]>. >>> For more options, visit https://groups.google.com/d/optout >>> <https://groups.google.com/d/optout>. >>> <step-49.cc <http://step-49.cc/>><mesh.inp><CMakeLists.txt><mesh.vtu> >> >> >> -- >> The deal.II project is located at http://www.dealii.org/ >> <http://www.dealii.org/> >> For mailing list/forum options, see >> https://groups.google.com/d/forum/dealii?hl=en >> <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/Erk7pxwTEj4/unsubscribe >> <https://groups.google.com/d/topic/dealii/Erk7pxwTEj4/unsubscribe>. >> To unsubscribe from this group and all its topics, send an email to >> [email protected] >> <mailto:[email protected]>. >> For more options, visit https://groups.google.com/d/optout >> <https://groups.google.com/d/optout>. >> >> >> -- >> Yuxiang "Shawn" Wang, PhD >> [email protected] <mailto:[email protected]> >> +1 (434) 284-0836 >> >> -- >> The deal.II project is located at http://www.dealii.org/ >> <http://www.dealii.org/> >> For mailing list/forum options, see >> https://groups.google.com/d/forum/dealii?hl=en >> <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] >> <mailto:[email protected]>. >> For more options, visit https://groups.google.com/d/optout >> <https://groups.google.com/d/optout>. >> <Job-1.inp> > > > -- > The deal.II project is located at http://www.dealii.org/ > <http://www.dealii.org/> > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > <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/Erk7pxwTEj4/unsubscribe > <https://groups.google.com/d/topic/dealii/Erk7pxwTEj4/unsubscribe>. > To unsubscribe from this group and all its topics, send an email to > [email protected] > <mailto:[email protected]>. > For more options, visit https://groups.google.com/d/optout > <https://groups.google.com/d/optout>. > > > -- > Yuxiang "Shawn" Wang, PhD > [email protected] <mailto:[email protected]> > +1 (434) 284-0836 > > -- > The deal.II project is located at http://www.dealii.org/ > <http://www.dealii.org/> > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > <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] > <mailto:[email protected]>. > For more options, visit https://groups.google.com/d/optout > <https://groups.google.com/d/optout>. -- 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.
