Farzin, it appears what you want to do is almost what DoFTools::make_periodicity_constraints( https://www.dealii.org/current/doxygen/deal.II/namespaceDoFTools.html#a929249499b1e5624728d212e90a8e037) gives you. You "only" have to apply the additional displacement +H(\chi_+-\chi_-). You can check for the support points of the constraints using DoFTools::map_dofs_to_support_points( https://www.dealii.org/current/doxygen/deal.II/namespaceDoFTools.html#a5514e4f59ea659f63953d62ca429eaff) assuming that you are using a nodal finite element. With this information, you should be able to calculate the displacement for each constrained dof pair and add it using AffineConstraints::set_inhomogeneity ( https://www.dealii.org/current/doxygen/deal.II/classAffineConstraints.html#a370b2b485a03fe2bb82c1098a6bdfc4c ).
Best, Daniel Am Fr., 12. März 2021 um 12:29 Uhr schrieb Farzin Mozafari < mozafari.far...@gmail.com>: > > > Dear all, > > > > I am a new user of Deal ii. I am going to solve a multiscale problem in > solid mechanics using Deal ii. To this end, I am going to model a simple > cube with a linear elastic material under periodic boundary condition as > the RVE of the material containing its microstructural features. Then by > applying simple shear and uniaxial strain-controlled boundary conditions, I > am going to determine the associated components of the averaged stress and > elasticity tensor. > > > > The way that usually people use is to involve an H matrix (as explained in > the attached file). In the literature, this method is called the > generalized periodic boundary condition. Indeed, this method makes us able > to directly apply shear and uniaxial loading through the H matrix > components (i.e. Displacement gradient matrix). I am not sure how can I > implement this method in deal II. > > > > I would be grateful if anyone gives me a piece of advice to carry it out > properly. > > > > Regards > > Farzin > > > > > > -- > 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/7D60CE96-E78D-45FA-A5C7-C51AAD2E3918%40hxcore.ol > <https://groups.google.com/d/msgid/dealii/7D60CE96-E78D-45FA-A5C7-C51AAD2E3918%40hxcore.ol?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/CAOYDWbLbQmfrgPJxsoYc8k%2B54vF6rugJ%2B4hHuA61mfDv9HKA3A%40mail.gmail.com.
