COMPLAS 2021 is the 16th conference of the COMPLAS Series.
The COMPLAS conferences started in 1987 and since then have become established events in the field of computational plasticity and related topics. The first fifteen conferences in the COMPLAS series were all held in the city of Barcelona (Spain) and were very successful from the scientific, engineering and social points of view. We intend to make the 16th edition of the conferenceanother successful edition of the COMPLAS meetings.
The objectives of COMPLAS 2021 are to address both the theoretical bases for the solution of nonlinear solid mechanics problems, involving plasticity and other material nonlinearities, and the numerical algorithms necessary for efficient and robust computer implementation. COMPLAS 2021 aims to act as a forum for practitioners in the nonlinear structural mechanics field to discuss recent advances and identify future research directions.
Scope
COMPLAS 2021 is the 16th conference of the COMPLAS Series.
Cold gas dynamic spraying (CGDS), as a high strain rate shearing and innovative solid-state technique enables to rapidly develop additive manufacturing and coating for metal deposition. This paper investigates the development and evolution of various interfacial bonding characteristics during high strain rate shearing process through Multiphysics numerical simulations of single particle impact. Two different particle-based modeling strategies such as smoothed particle hydrodynamics (SPH), molecular dynamics (MD) are investigated using commercial software ABAQUS/Explicit and LAMMPS, respectively. To separate the difficulties related to complex metallurgy of alloys, our first investigations focus on pure aluminum. The Johnson-Cook (J-C) constitutive model is used to describe the high strain rate self-consolidation process in SPH modeling. Embedded Atom Method (EAM) is used to describe the interactions between Aluminum atoms in MD modeling. The predictions from the different particle-based models are compared with each other and with experimental results. Through the investigations, SPH numerical approach has strong advantage in capturing the phenomena that occur during the cold spray process. It is able to describe the complex features of particle and substrate, especially in the interface vicinity. At the same time, MD numerical approach gives the fundamental understanding of the deposition behavior at the atomistic level. The key finding is the strong relationship between the un-uniform distribution of shear strain and jet formation during high-speed collision. Plastic strain along with an increase of temperature lead to thermal softening of pure Aluminum resulting in metallurgical bonding at the interface.
Abstract Cold gas dynamic spraying (CGDS), as a high strain rate shearing and innovative solid-state technique enables to rapidly develop additive manufacturing and coating for metal [...]
Casting defects can be predicted in advance of practice and countermeasures can be taken to improve casting quality and increase productivity. Applying the casting filters is a method of improvement methods for defects caused by unsuitable molten metal flow. Casting filters have the effect of removing inclusions in molten metal and rectifying the flow. However, specific conditions such as the type, pore size, and setting position of the casting filter are not clear. Casting filter conditions are determined by conventional empirical rules that are not theoretical. In the other view, the use of casting CAE is essential to realize front-loading for the process design process, in which casting defects are predicted in advance of practice and countermeasures are taken. In the previous study, K. Taki et al. performed direct observation of mold filling and flow simulations passing through the casting filter. The particle-based COMINA CAE software was used for the flow simulation of molten metal in complex interior geometries in the filter. The calculations used a model of the filter that was reproduced on an X-ray CT system. To inspect the filter performance it was necessary to make a small and simplified filter model, which is called the 1/4 model of filter. In the present study, the flow dynamics through the filter are investigated using various 1/4 models. The 1/4 model maintains permeability on the surface and porosity of volume while halving the dimensions. As a result, we succeed in reproducing the flow behaviour of molten metal when it passed through the filter by setting the particle size of molten metal to 1/16 of the filter’s pore diameter. Further, we try to evaluate the performance of the filter by extending the calculation target from only the area around the filter to the entire mold. If mold filling behavior for the mold with filter could be simulated, it wouldbe used effectively in casting geometry design and defect countermeasure.
Abstract Casting defects can be predicted in advance of practice and countermeasures can be taken to improve casting quality and increase productivity. Applying the casting filters [...]
F. Itakura, T. Yamada, Y. Maeda, A. Hasuno, Y. Mochida
particles2023.
Abstract
Disturbance of the molten metal flow during ladle pouring before the plunger advancing in the aluminium alloy die casting process can cause entrapment defects of air and oxide film. Slow pouring to control the turbulence of the flow front reduces productivity due to increased cycle time. Further, the risk of cold flake formation increases caused by large temperature drops in accordance with the long cycle time. On the other hand, rapid pouring is desired to improve productivity, but the risk of air entrapment increases. Therefore, quick and quiet pouring is desired in the ladle pouring process. In the present study, we focus on variable tilting speed as a method to achieve good ladle pouring. The effects of variable ladle tilting speed and switching time on the wave behavior of molten metal are investigated in visualization experiments and simulations. The flow behaviours in ladle pouring are simulated using ”COLMINA CAE”, which is the casting analysis software by particle-based SPH method. Furthermore, the plunger advancing process is also examined. From the simulation results, the variable tilting speed from fast to low can suppress the rise of the maximum wave height of molten aluminium alloy. However, the pouring completion time is longer. Further, the falling position of molten metal poured from the ladle varied with changing tilting speed. And then, the wave height is influenced not only by ladle pouring but also by the plunger advancing process. These trends of wave behaviour obtained in the simulation are similar to that of the actual phenomenon. Therefore, the present simulation method can accurately estimate the ladle pouring process and plunger advancing process. So, casting CAE is an effective tool for exploring die casting conditions.
Abstract Disturbance of the molten metal flow during ladle pouring before the plunger advancing in the aluminium alloy die casting process can cause entrapment defects of air and oxide [...]
The process of laser beam welding is simulated using the Weakly-Compressible Smoothed Particle Hydrodynamics (WCSPH) and the Incompressible SPH (ISPH) methods. The presented models consider significant physical effects such as heat conduction, temperature-dependent surface tension with wetting, the phase transitions melting and solidification, and an evaporation-induced recoil pressure. Here, particular emphasis is placed on the modeling differences between the WCSPH and ISPH methods. Then, both methods are evaluated in terms of their accuracy and performance in the simulation of deep penetration laser beam welding with oscillating laser power.
Abstract The process of laser beam welding is simulated using the Weakly-Compressible Smoothed Particle Hydrodynamics (WCSPH) and the Incompressible SPH (ISPH) methods. The presented [...]
N. Sani, J. Quist, K. Jarateg, A. Bilock, L. Cordova, F. Edelvik
particles2023.
Abstract
Additive Manufacturing (AM) has been a subject of significant attention from both industrial manufacturers and research communities. However, several challenges hinder the widespread implementation of this technology in the industry. Powder recoating is a crucial step in powder-bed AM process that involves achieving a uniformly packed bed of powder particles that are later melted by an energy source, such as a laser or electron beam. One of the main challenges is calibrating the contact model parameters accurately to match the flowability and spreadability of specific powder alloys. This paper proposes a Discrete Element Method (DEM) model calibration framework based on surrogate model optimisation. The study utilises a Revolution Powder Analyser (RPA) as the experimental reference system. The proposed method is demonstrated with two AM powder samples, Ti64 and Inconel 718. The results indicate that particle-particle friction, rolling resistance, and van der Waals (vdW) surface energy significantly affect the system responses. Furthermore, the validation results show good correspondence between the simulation with calibrated parameters and experimental data. Overall, proposed calibration framework has the potential to optimise powder recoating and to improve the accuracy and effectiveness of the additive manufacturing.
Abstract Additive Manufacturing (AM) has been a subject of significant attention from both industrial manufacturers and research communities. However, several challenges hinder the [...]
Granular materials play a crucial role in various geotechnical, mining, and bulk handling applications. Understanding their mechanical properties is essential for optimal use in these industries. Traditional experimental methods like Cone Penetration Test (CPT) and open pile testing have limitations on their repeatability and offer little insight into the contact mechanics. The Discrete Element Method (DEM) is a powerful tool for investigating and simulating granular material behaviour at the element scale and provides deeper understanding in geometry-material interactions. However, due to computational costs, spherical particles are often preferred, though they may not always capture realistic particle interactions. In the current study, the packing density and the penetration resistance of particle beds with different particle shapes, including sphere, multi-spheres and polyhedrons, are compared using a plate penetration test modelled in DEM. Sensitivity analyses are performed for sliding friction, consolidation pressure, and Particle Size Distribution (PSD). Results indicate that polyhedral shapes show lower penetration resistance compared to spherical and multi-spherical shapes. Sliding friction has the most significant impact on resistance, while consolidation pressure has minimal effect on porosity. The study highlights the importance of particle shape in granular media modelling and emphasizes the need for further research in this area.
Abstract Granular materials play a crucial role in various geotechnical, mining, and bulk handling applications. Understanding their mechanical properties is essential for optimal [...]
Incompressible fluid analysis using the ISPH or MPS methods requires the solution of the pressure Poisson equation, which takes up most of the overall computation time. In addition, the iteration number for solving pressure Poisson equations may increase as the simulation model scale increases. This is a common problem in particle methods and the other implicit time integration solvers. In different methods, FEM, etc., good quality preconditioning, such as multigrid preconditioning, can significantly improve the convergence of iterative solution methods. There are two types of multigrid preconditioners, algebraic multigrid and geometric multigrid methods, but there are few examples of their application in particle methods. In this study, we attempted to develop a framework for a geometric multigrid preconditioner for solving the pressure Poisson equation in the ISPH. First, we focused on the geometric multigrid preconditioner using background cells, which are used for neighboring particle search, and implemented it on a GPU environment. Through a simple dam-break problem, we compared the computation time between the Conjugate gradient (CG) solver with a traditional preconditioner and the CG solver with a geometric multigrid preconditioner. We confirmed that the background cell-based geometric multigrid preconditioner is practical for the ISPH method.
Abstract Incompressible fluid analysis using the ISPH or MPS methods requires the solution of the pressure Poisson equation, which takes up most of the overall computation time. In [...]
This work extends a multi-phase mixing model framework designed for a Smoothed Particle Hydrodynamics context. Specifically, we propose a higher-order variation using the first-order accurate Generalised Finite Difference differential operators to construct an incompressible scheme for simulating fluid-solid coupled systems resolved via a continuum mixture model. The proposed scheme incorporates inter-phase shear between phases and the viscosity dependency of the solid phase concentration. The scheme is verified by simulating a modified lid-driven cavity case at Re = 1000. In this simulation, our method was capable of treating initially discontinuous concentration fields with a maximum solid volume concentration of 0.5 and a solid-to-fluid density ratio of 4.
Abstract This work extends a multi-phase mixing model framework designed for a Smoothed Particle Hydrodynamics context. Specifically, we propose a higher-order variation using the [...]
I. Bensemmane, W. Leclerc, N. Ferguen, M. Guessama
particles2023.
Abstract
Thermal Barrier Coatings (TBCs) are multilayer systems used in Ni-based superalloy components for gas turbine blades submitted to high temperatures resulting in the development of high thermal stresses. The intricacy of their microstructure coupled to severe environmental conditions lead to their premature failure according to complex mechanisms involving notably creep and Coefficient of Thermal Expansion (CTE) mismatch effects. This work aims to investigate the development of thermal residual stress within TBCs during a heating step using a numerical model based on Discrete Element Method (DEM). The suitability of such an approach is investigated in terms of stress field distribution in comparison to Finite Element Method (FEM). Results reveal the capability of the proposed DEM approach to simulate creep phenomenon in a TBC system under thermal loading and predict accurately thermal stresses leading to failure.
Abstract Thermal Barrier Coatings (TBCs) are multilayer systems used in Ni-based superalloy components for gas turbine blades submitted to high temperatures resulting in the development [...]
Internal erosion, caused by seepage flow inside the soil, accelerates soil failure during a natural disaster. Numerical simulation can be an effective tool to quantitatively evaluate the relationship between internal erosion and the instability of the ground as a whole. Internal erosion and multiphase flow simulation of fluid and granular materials with a particle size distribution require coupling simulations that can represent the interaction between particles and pore water and the movement of particles. There are two main types of coupling models: ”Resolved coupling model,” which can calculate detailed flow and fluid forces, and ”Unresolved coupling model,” which is based on empirical drag and seepage flow models. Previous studies have indicated that both models should be judged appropriately based on the ratio of particle-fluid spatial resolution. However, applying a resolved coupling model to the vast number of soil particles that make up the ground is impractical from a computational cost perspective, and empirical unresolved coupling model has difficulty in representing localized failures such as internal erosion. Therefore, developing a new coupling model that satisfies both computational accuracy and efficiency is desirable. In this study, we applied ISPH (Incompressible Smoothed Particle Hydrodynamics) for fluid analysis and DEM (Discrete Element Method) for soil particles to develop a fluid-soil coupling simulation model that can directly represent the movement of soil particles during the internal erosion process. Through numerical experiments using a particle layer with the vertical upward flow, we understand the limitations of the conventional coupling model and propose a new hybrid type of semi-resolved coupling model that combines these two models appropriately.
Abstract Internal erosion, caused by seepage flow inside the soil, accelerates soil failure during a natural disaster. Numerical simulation can be an effective tool to quantitatively [...]