<div dir="ltr"><div dir="ltr">On Wed, Jun 17, 2020 at 4:05 PM Jacob Faibussowitsch <<a href="mailto:jacob.fai@gmail.com">jacob.fai@gmail.com</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;">Hello,<div><br></div><div>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:</div><div><br></div><div>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.</div></div></blockquote><div><br></div><div>This is how PyLith works: <a href="https://github.com/geodynamics/pylith">https://github.com/geodynamics/pylith</a></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;"><div>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.</div></div></blockquote><div><br></div><div>No, it would have to be coded.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;"><div>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.</div></div></blockquote><div><br></div><div>No, it would have to be coded.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;"><div>4. Is there a mechanism by which you can classify existing elements as cohesive elements?</div></div></blockquote><div><br></div><div>See 1.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;"><div>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.</div></div></blockquote><div><br></div><div>There is nothing for non-conforming interfaces.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;"><div>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. </div></div></blockquote><div><br></div><div>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</div><div>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</div><div>the solver clears up the hardest part.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;"><div><div>
<div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><div>Best regards,<br><br>Jacob Faibussowitsch<br>(Jacob Fai - booss - oh - vitch)<br>Cell: (312) 694-3391</div></div>

</div>


<br></div></div></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>