How to SEMI-AUTOMATICALLY apply Periodic Boundary Conditions in ABAQUS

Do you know that you can apply Periodic Boundary Conditions (PBCs) on any 2D domain without resorting to tedious manual methods or use of an ABAQUS plugin? This is what this video is about. It introduces a semi-automatic way of applying PBCs on 2D domains. The video opens up with an theoretical introduction to PBCs before a description of the PBC2DGen Semi-automatic script. Finally, two case studies (uniaxial tension along X-axis and pure Shear XY simulations) were undertaken to demonstrate the use of the code in real problems.
#cmvideos #abaqustutorial #periodicboundaryconditions #pbc #pbcgen2d #pbcgenlite
Time stamps
--------------------------------------
You can quickly navigate through the videos using these chapters.
00:00 - Intro
09:08 - Boundary Conitions for two case studies studied.
09:31 - Info: All scripts are available in description section.
09:41 - Details of the Scripting of PBCGen2D_Semi-automatic
23:05 - Case Study I: Uniaxial Tensile XX Simulation setup
23:42 - Case Study II: Pure Shear XY Plane setup
25:17 - Visualization of Simulation results
26:04 - Extracting Stress-strain Data
28:16 - Outro
-------------------------------------
📥* DOWNLOAD PBCGEN2D SEMI-AUTO SCRIPT AND EXCEL FILES
If you want to get hold of the script for imposing PBCs semi-automatically as well as template Microsoft Excel files used for generating stress-strain plots, then visit this link (you will need to sign up first to CMVideos Insider Group) before you are directed to the site where the resources are. Here is the link: cmvig.cmvideos.org/
--------------------------------------
📥* DOWNLOAD LINK TO SOFTWARE
1. PBCGen2D: A software for applying periodic boundary condition
Download softwaree: cmvideos.ck.page/products/pbc...
2. PBCGenlite: A lite version of software for applying periodic boundary condition
Download software: cmvideos.ck.page/products/pbc...
3. MontCarlGen2D: A software for creating random inclusions in 2D RVEs
Download software here: cmvideos.ck.page/products/mon...
📹 RELATED VIDEOS
-----------------------------------
1. How to manually apply Periodic Boundary Conditions in ABAQUS - • How to manually apply ...
2. How to AUTOMATICALLY apply PBCs in ABAQUS using PBCGenLite - • How to AUTOMATICALLY a...
3. Theory of Periodic Boundary Condition implementation in FEM - • Theory of Periodic Bou...
-----------------------------------
🕮 CM VIDEOS E-BOOK SERIES
I am writing a series of E-books to complement the contents on this channel. If you want to join my writing journey and get weekly newsletter about this, please subscribe via this link: cmvig.cmvideos.org/ebook-news...
-----------------------------------
💁‍♂️ TIP JAR: If you want to support the channel, you can leave a tip here. I am saving for new camera and lighting, so all tips will be reinvested into this. Thanks a lot - cmvideos.ck.page/products/tha...
CMVIDEOS WEBSITE
--------------------------------------
*Website under development at: www.cmvideos.org
👨‍👩‍👦 CM Video Insider
--------------------------------------
Become a CM Video Insider (100% free) to access modelling resources, advanced workshops, bonus trainings, behind the scene contents and regular CM Video updates: cmvig.cmvideos.org/
✔️*Link to Michael Okereke's Finite Element Applications textbook
--------------------------------------
Okereke, M. and Keates, S., 2018. Finite element applications. Cham: Springer International Publishing AG. amzn.to/3Bv54Di
---------------------------------------
My KZread Accessories
• Rode NT USB Mini Studio Micro USB - amzn.to/3SXha0y
• LOGITUBO Full HD 1080P Webcam - amzn.to/3MoqHvf
• Microsoft Surface Book 2 - amzn.to/3CwlT2a
• Microsoft Arc Bluetooth Mouse - Black - amzn.to/3EwPbAu
• Microsoft Surface Pen - amzn.to/3CkYsc8
• Iiyama LCD Monitor - amzn.to/3CRfpw7
• External Hard Drive - USB 3.0, Black - amzn.to/3RX5xWf
• Webcam Stand - amzn.to/3SZgwQ0
• Samsung Galaxy S21 Smartphone - amzn.to/3iL4Hgl
• Jabra Evolve 40 MS Stereo Headset - amzn.to/3EC2krP
• Green Screen 5x6.5ft - amzn.to/3RXldZw
• Scripting book: Hardcover Notebook - Black - amzn.to/3TbvDFG
Other Accessories
------------------------------------------
Business Projector - amzn.to/3ThN8o4
Please note the above are affiliate links and I do get commissions when you make purchases from those links.
🔔 Subscribe
--------------------------------------
If this is the kind of content you like, please do subscribe to this channel and click the notification button so when contents like this are made, you will be the first to see it. If you also have ideas of videos you want me to make, do suggest those in the comment section and I will consider them.

Пікірлер: 20

  • @josepabloquesadamolina8036
    @josepabloquesadamolina8036 Жыл бұрын

    Great content! Thanks Dr. Michael Okereke.

  • @MichaelOkereke

    @MichaelOkereke

    Жыл бұрын

    My pleasure @Jose Pablo!

  • @Valera197615
    @Valera197615 Жыл бұрын

    Very thanks Dr. Michael Okereke.

  • @MichaelOkereke

    @MichaelOkereke

    Жыл бұрын

    You are very welcome @Valerii

  • @kiamashoufi1296
    @kiamashoufi1296 Жыл бұрын

    Dear Dr.Okereke, thank you for your informative tutorials. I'm wondering if there is a semi/automatic method to impose PBC on 3D RVEs which are made from two parts. I'm working on a 3D RVE of a UD composite and I have to model fibers and matrix in two different parts to be able to define cohesive surfaces between them. Moreover, my 3D RVE has fibers on edge. Do you have any idea how should I define PBC in this case ? Thanks a lot.

  • @MichaelOkereke

    @MichaelOkereke

    Жыл бұрын

    This is definitely not a straightforward way of imposing PBCs on 3D RVEs. I published about this a while back so you can read about it in this publication: doi.org/10.1016/j.compstruct.2016.10.114 You should get some pointers in that publication. I do not have a publicly available code for doing such but am working on it. For the situation with cohesive surfaces, yes, you have to allow all components to be independent (unattached) so you can specify them independently. It can prove difficult but it is doing. I do not have a video about this but can explore this in future.

  • @syamilzafriify
    @syamilzafriify Жыл бұрын

    Thanks for this insightful video. I'm wondering weither I can use this PBCGen2D software for a hexagonal shape marix unit cell ?

  • @MichaelOkereke

    @MichaelOkereke

    Жыл бұрын

    Hello @Syamil, I will say no. The PBCGen2D software works on 2D Quadrilateral RVEs so it will unfortunately not work in this case.

  • @martinascordio1648
    @martinascordio1648 Жыл бұрын

    Thanks sir, your videos are very helpful. I was wondering how to do the same procedure on a lattice RVE in which nodes must be paired considering opposite surfaces. It's a huge effort doing that manually expecially when each face has around 200 nodes. What do you think? Thanks a lot

  • @MichaelOkereke

    @MichaelOkereke

    Жыл бұрын

    Hello @Martina, I am glad that the video is helpful. For 2D RVE lattices, I did a video in the past using a the automatic way of implementing PBCs (see here kzread.info/dash/bejne/Zo140K2ol73Ffso.html). For surfaces, as you describe here, it is a different implementation. This manual or semi-automatic approach will not be suitable as there are too many nodes to be able to use them. You will have to use an automatic approach where a script searches through the surfaces and note node numbers, nodes coordinate positions and their collinearity with spatially complementary opposite nodal sets. I plan to bring some theoretical and numerical highlight to these things in future.

  • @martinascordio1648

    @martinascordio1648

    Жыл бұрын

    @@MichaelOkereke thank you very much Dr. Michael Okereke

  • @jiyadferoke2600
    @jiyadferoke2600 Жыл бұрын

    Sir...iam using ansys workbench for my research work..how I can apply PBC in ansys workbench or APDL

  • @MichaelOkereke

    @MichaelOkereke

    Жыл бұрын

    Hello @Jiyad, I am not sure I know how this is done in ANSYS. It is something I am still investigating. Thanks for the question.

  • @RossMcLendon
    @RossMcLendon11 ай бұрын

    Very interesting video. Are you aware of the Micromechanics Plugin for Abaqus/CAE that DS provides? It was released at the end of 2016 and automates the imposition of PBCs (and a few other types of constraints) as well as defining loads and performing homogenization on RVE models. Is there a reason that you only constrain the normal component of displacement across the periodic faces? If the tangential component(s) of displacement aren't constrained, you're not actually going to get a periodic displacement field - the tangential displacements will vary between faces in a non-periodic manner. I also think that this will lead to a zero shear condition along the cut faces of the RVE as the tangential components of stress will necessarily go to zero at the boundary (although the zero shear in the derivative fields may not be immediately apparent due to low mesh refinement) and in general the strain and stress field will be non-periodic. One nice way to check the correctness of PBCs is to create a model with multiple tessellations and compare the results when you model one cell and two cells. The stress/strain solution should be identical at a given point.

  • @MichaelOkereke

    @MichaelOkereke

    11 ай бұрын

    Hello @RossMcLendon, Thanks for your detailed comment and pointing out a few things. On *Micromechanics Plugin for ABAQUS CAE*, I actually did not know about it. I only researched into it to find your article about it here: www.linkedin.com/pulse/micromechanics-plugin-abaqus-ross-mclendon/ I know ANSYS has a similar module as this but not ABAQUS. It will be interesting to have a play with it and more importantly understand the theoritical framework that underpins its development. I believe the ABAQUS license I hold with my organisation might not cover that plugin. On *Constraining the periodic faces of rve*, I am not sure you understand completely the framework I am working on. It is not constrained only in the normal direction, rather it is for both normal and transverse in otherwords x and y-axis. The principle here is to isolate first identify the corner nodes (the vertices of the quadrilateral) and use those to use as anchor points to apply the required loading and boundary conditions on the model. The nodes away from vertices are then implemented with a *EQUATION driven periodic boundary conditions formulation that constraints periodic nodes (directly opposite semi-colinear nodes) in both x- (normal) and y- (transverse) directions. It is actually a well established approach that has been published widely in the community and peer reviewed. A lot of authors have adopted this approach and as a student at Oxford University many years ago, I first read about it by an author called Kouznetsova and implemented this as part of my PhD work. I have since published about it extensively and here is a paper where I laid down the principles in detail doi.org/10.1016/j.commatsci.2012.12.036. I also gave detailed description about this in a 2018 Postgraduate textbook I authored and published by Springer which you can read here: link.springer.com/book/10.1007/978-3-319-67125-3 You may like to look at Chapter 8 about Boundary Conditions. On *Correctness of PBCs*, I have already explored this in the said publications above and shown it to be true. The journal paper has a parametric study where I explored different mesh descritations beyond just unit cells and it showed that the PBC implementation is right. Thanks once again for your detailed comments. I will love to catch up separately to continue this discussion. Cheers!

  • @RossMcLendon

    @RossMcLendon

    11 ай бұрын

    @@MichaelOkereke Okay, watching again I think you are defining them correctly - I was mistakenly thinking that you were describing the equation for N_A and N_B as involving only DOF 1 and N_C and N_D as involving only DOF 2. On review I see that you're describing two equations for N_C and N_D, one for DOF 1 and another for DOF 2. Assuming you also define similar equations for N_A and N_B (with N1 and N2), then what you're describing should be periodic. The Micromechanics Plugin uses a somewhat different approach to define PBCs, but the two results will be in essence equivalent apart from how loads are defined (with the plugin, you directly specify strain components rather than specifying relative displacements between periodic faces).

  • @MichaelOkereke

    @MichaelOkereke

    11 ай бұрын

    Certainly @RossMcLendon, I agree. There is definitely a kinematical constraint equation involving also N1 and N2 to make the model truly periodic. This is why I divided them into internal nodes (NA, NB, NC and ND) and the corner nodes (N1, N2, N3 and N4). By intelligent combination of displacements at corner nodes kinematically linked to the internal nodes, we can then define a canonical equation (which I discussed in the video) which is then subsequently implemented using *EQUATION to enforce the right periodicity argument. I am aware of your strain-tensor based approach for imposing PBCs. I had a PhD student who used that approach as part of his PBC implementation for a textile composite study. In fact, in that method, he was able to deal with non-periodic meshes and you can read up about this here: doi.org/10.1016/j.compstruct.2016.10.114. I intend to make a video on that approach to complement what I have presented here. I first got exposed to this displacement-based approach and love the intuitive way it works, something users who are already comfortable with Dirichlet BC-type approach, can easily understand. Your approach will certainly allow for applying combined loading easily based the strain tensor. I know there are other benefits. I will write you outside this forum to follow up on these interesting conversations. I look forward to playing with the Micromechanics plugin and see how useful and complementary it is to what I present on this channel. Thanks once again for taking time to comment on this.

  • @manojkumarneesu2773
    @manojkumarneesu2773 Жыл бұрын

    sir i have a doubt, I am modelling a composite laminate of 8 layers along with the frame on top and bottom. I am using VUMAt code and applying a predefined velocity in initial step on the ball and applying pressure on top frame in step one. I need all my boundary conditions and pressure on the frame has to be applied before the ball hits plate. so how can i do it. Please help me in this case

  • @MichaelOkereke

    @MichaelOkereke

    Жыл бұрын

    Hello @Manoj, this is too specific a request and I don't usually honour this sort of requests. It seems a project you are working on and I'd leave you to figure it out. I tend to make general videos to give ideas to viewers rather than specifically take a project and offer solutions here. As an academic I want my students to do their projects themselves. Good luck with the work 👍.

  • @manojkumarneesu2773

    @manojkumarneesu2773

    Жыл бұрын

    @@MichaelOkereke oke sir i respect your words....thanks 🙂