Hi Matt, It is difficult to give a clear answer without some more information on how you want to compute with fibers and where the material model for the fibers is coming from.
If you have a way to compute the fourth-order tensor directly then an approach like the one in the linked paper works. A good alternative is to compute a vector field which represents the anisotropic property - a good description of this for cardiac simulations is given in Pasquale's paper from early this year: https://arxiv.org/pdf/2101.10960.pdf that process (similar to what we do at UNC) involves computing a fiber field as a finite element field, so you don't necessarily need to store quadrature data - you can just compute values that you need cell-by-cell. Best, David ________________________________ From: dealii@googlegroups.com <dealii@googlegroups.com> on behalf of Matthew Rich <mjric...@gmail.com> Sent: Thursday, September 30, 2021 9:49 PM To: deal.II User Group <dealii@googlegroups.com> Subject: [deal.II] Implementing a fiber reinforcement in a biphasic model Hi all, I am trying to add some anisotropy to my model via fibers. These fibers would only be active in tension and have different orientations depending on spatial location in the model. My plan of attack was to simply store and extra value at the quadrature points that had a direction representing the fiber and update directions after each time step. (is this a good way to go about doing that). I see a lot of papers doing fiber models in ABAQUS or COMSOL, but I never see any details on how to do it so I am guessing it is trivial. Then I stumbled upon https://doi.org/10.1016/j.compstruc.2020.106334 which uses dealii to implement a fiber model that has two grids and some sort of junction between the grid representing the bulk material and one presenting the fiber. The formulation seemed way more involved than what I was going to try. Am I making a gross simplification in my proposed approach. Matt -- 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 dealii+unsubscr...@googlegroups.com<mailto:dealii+unsubscr...@googlegroups.com>. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/16d823d0-da31-40fd-8eeb-8c7162fc42e0n%40googlegroups.com<https://groups.google.com/d/msgid/dealii/16d823d0-da31-40fd-8eeb-8c7162fc42e0n%40googlegroups.com?utm_medium=email&utm_source=footer>. -- 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 dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/BN7PR03MB43569AFAA4FDD8B545EB2F59EDB59%40BN7PR03MB4356.namprd03.prod.outlook.com.
