| Is coupling PD with FEM the way forward to solve in an efficient way crack propagation problems? Tao Ni1; Greta Ongaro2; Pablo Seleson3; Mirco Zaccariotto4; Ugo Galvanetto2; 1HOHAI UNIVERSITY, Nanjing, China; 2UNIVERSITY OF PADOVA, Padova, Italy; 3OAK RIDGE NATIONAL LABORATORY, Oak Ridge, United States; 4UNIVERSITY OF PADUA, Padova, Italy; PAPER: 217/Geomechanics/Plenary (Oral) SCHEDULED: 15:55/Thu. 24 Oct. 2019/Athena (105/Mezz. F) ABSTRACT: Environmental, economic and safety concerns require more and more precise capabilities to perform the life cycle assessment of engineering structures. Therefore, structural engineers should be capable of describing all stages of the structural life even that involving the propagation of cracks and branching under complex loading conditions. The description of the propagation of cracks in structural materials, however, is still an open problem. The unavoidable presence of discontinuities prevents a direct application of the methods based on Classical Continuum Mechanics (CCM). Recently, Peridynamics (PD) has been proposed [1, 2] as a theory in which cracks are not part of the problem but part of the solution; PD is based on integral equations that do not make strong assumptions on the continuity of the displacement field. The integrals of the peridynamic theory are computed on a neighborhood of each material point, which is affected, in a nonlocal way, even by points that are not in direct contact with it. As a consequence of such a nonlocality, computational methods based on PD are usually more computationally expensive than those based on CCM. Several researchers are trying to couple computational methods based on CCM with those based on PD to obtain a computational tool able to simulate crack propagation in an efficient way [3-5]. Coupling two different continuum theories is not straightforward. In our presentation, coupling is realized at the discrete level between the standard displacement version of the Finite Element Method and a meshless version of the Ordinary State based PD. The domain is divided in two portions, one discretized with FEM and the other with OSBPD. If a perfect bonding between the displacements of the two portions is imposed, some out of balance forces are generated. The paper evaluates the magnitude of the out of balance forces and discusses some ways to reduce them. References: 1) S.A. Silling, Reformulation of elasticity theory for discontinuities and long-range forces, J. Mech. Phys. Solids, 48 Issue: 1, 175-209, (2000). 2) S. A. Silling, M. Epton, O. Weckner, J. Xu, E. Askari, Peridynamic states and constitutive modeling, J. Elasticity 88 (2007): 151-184. 3) G. Lubineau,Y. Azdoud, F. Han, et al., A morphing strategy to couple non-local to local continuum mechanics, J. Mech. Phys. Solids, 60 Issue: 6, 1088-1102, (2012). 4) U. Galvanetto, T. Mudric, A. Shojaei, M. Zaccariotto, An effective way to couple FEM meshes and Peridynamics grids for the solution of static equilibrium problems, Mechanics Research Communications 76, 41-47, (2016). 5) M. Zaccariotto, T. Mudric, D.Tomasi, A. Shojaei, U. Galvanetto, Coupling of FEM meshes with Peridynamic grids, Comp. Meth. Appl. Mech. Eng., Volume: 330, Pages: 471-497, (2018). |
