On Jun 19, 2020, at 12:53 PM, Jacob Faibussowitsch <jacob....@gmail.com<mailto:jacob....@gmail.com>> wrote:
Hello, Thank you both for your comprehensive replies. Matt: Thanks for the rundown, I will take a look at PyLith and its cohesive element impls. Blaise: The project is still in very early stages as you can tell, so everything is up in the air. Currently we are weighing using a cohesive element approach vs an element extinction approach but we can certainly add a variational model into the mix. We had leaned towards the cohesive element approach as we also need to link Quantum Monte Carlo generated data to the model, which seemed logically the easiest to do with cohesive elements. Are there any examples of your codes online? The only reference I could find online was in this link https://www.math.lsu.edu/~bourdin/defectmechanics/ So far, you are listing numerical approach. What is your fracture _model_? Both cohesive elements and element erosion approach require some a priori knowledge of the crack path. Is this acceptable for you? My code is released under a BSD license at https://github.com/bourdin/mef90 Regards, Blaise Best regards, Jacob Faibussowitsch (Jacob Fai - booss - oh - vitch) Cell: (312) 694-3391 On Jun 19, 2020, at 10:15 AM, Blaise A Bourdin <bour...@lsu.edu<mailto:bour...@lsu.edu>> wrote: On Jun 18, 2020, at 5:28 AM, Matthew Knepley <knep...@gmail.com<mailto:knep...@gmail.com>> wrote: On Wed, Jun 17, 2020 at 4:05 PM Jacob Faibussowitsch <jacob....@gmail.com<mailto:jacob....@gmail.com>> wrote: Hello, I am looking to perform large scale fracture and crack propagation simulations and have a few questions regarding PETSc support for this. Specifically I am looking for cohesive surface element support with a few twists: 1. Is there support for zero thickness surface elements? For example modeling virtually flat patches of adhesives holding together two larger structures being pulled apart. This is how PyLith works: https://github.com/geodynamics/pylith<https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fgeodynamics%2Fpylith&data=02%7C01%7Cbourdin%40lsu.edu%7C72bd6fed59354581bd8108d81479be5a%7C2d4dad3f50ae47d983a09ae2b1f466f8%7C0%7C0%7C637281860412268459&sdata=bZO7nz9UDVYMC6ipO0EmVxBXoG72VeQjW2SyeTssh6E%3D&reserved=0> 2. Is there support for “joining” two possibly distinct meshes with cohesive surface elements? For example say I have two distinct cylinders representing fibers which would “touch" to form an X shape. No, it would have to be coded. 3. In a similar vein, is there support for a mesh to fracture entirely along a crack formed through the cohesive elements? Imagine the aforementioned X configuration separating entirely into two separate cylinders again. No, it would have to be coded. 4. Is there a mechanism by which you can classify existing elements as cohesive elements? See 1. 5. Is there an already implemented way of imposing tie-constraints between independent meshes? This would potentially be used to tie high order cohesive cells which would have a non-conforming interface to the “regular” mesh. There is nothing for non-conforming interfaces. From googling I have come across DMPlexCreateHybridMesh(), DMPlexConstructCohesiveCells(), and DMPlexCreateCohesiveSubmesh(). While these do implement cohesive cells these functions don’t at first glance seem to allow one to implement the above. Having worked with cohesive elements for more than a decade, I would be cautious about a new code using them for fracture. To me, it appears that variational fracture codes, like those from Blaise Bourdin and J. J. Marigo's group have much better geometric flexibility, and Maurini's work on the solver clears up the hardest part. I definitely concur with this and would be happy to help… Blaise -- A.K. & Shirley Barton Professor of Mathematics Adjunct Professor of Mechanical Engineering Adjunct of the Center for Computation & Technology Louisiana State University, Lockett Hall Room 344, Baton Rouge, LA 70803, USA Tel. +1 (225) 578 1612, Fax +1 (225) 578 4276 Web http://www.math.lsu.edu/~bourdin
