On 08/17/2017 07:35 PM, Jie Cheng wrote:

Thanks for your answer. Now I understand that the problem was deal.II is expecting an integer as the surface name so that it can be converted to boundary id. So I changed my surface name from "top" to "1" (without quotes). But all the faces that are at the boundary still have the boundary id = 0. The method that I used to check the boundary id is cell->face(face)->boundary_id().

That should be the way to do it. I don't know whether our Abaqus reader actually reads boundary values, though. You may have to take a look at the source code.


I expected some boundary ids to be 1. But deal.II was able to figure out how many faces are at the boundary. What did I miss? Do you have a sample Abaqus file that deal.II can process (read mesh and boundary IDs)?

Take a look at the ones here:
  https://github.com/dealii/dealii/tree/master/tests/grid/grids/abaqus

Best
 W.


--
------------------------------------------------------------------------
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 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.

Reply via email to