Hi Tariq, a very basic approach could be like this, give the cells in your triangulation different material_id's, iterate (in the standard dealii-way) over the cells, on each cell loop over the faces and look, if the neighbor (if this face is not on the boundary) has the same material_id(), if not, you have found a face on the interface between the materials and you can assemble whatever you want.
But what problem do you want to solve in mechanics or what material law do you use? for linear elastic materials I use another approach for assembling, but it depends on what you need at the interface. Maybe this helps, Best, Martin ________________________________ Von: "[email protected]" <[email protected]> An: [email protected] Gesendet: Dienstag, den 11. Januar 2011, 14:50:05 Uhr Betreff: [deal.II] internal boundaries Dear all, I am trying to solve a problem from structural mechanics involving two (finite) elastic domains with different elastic constants. In order for this problem to be well posed I have to specify internal boundary conditions, in my case a vanishing "jump" of a function of the solution at the boundary. As a consequence the variational formulation of my problem contains terms that involve objects on different sides of the internal boundary: more precisely boundary integrals of functions assembled in a cell on one side multiplied by test functions on the other side. The best way to solve this problem I can think of is to figure out the test/shape function in the same cell that corresponds to the one I really want - which should be possible since both cells are images of the same unit cell . Has anybody experiences with something like this or maybe an idea that solves the problem more elegantly? Thanks in advance, tariq
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
