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.
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 [...]
The term “4D Printing” (4DP) is defined as the ability for a part produced using an additive manufacture process to change its shape when activated by or exposed to one or more stimuli over time. This emerging technology offers unique advantages over conventional Additive Manufacturing (AM) by extending the three dimensions of space into the fourth dimension of time. 4DP parts can be programmed to actuate passively without the need for an external power source such as an electromechanical or other active system, thereby reducing the probability of failure and the complexity of components. This work attempts to address some of the challenges faced by the design engineer in a project team when producing technical documentation to specify the desired shape transformation of a 4DP part with a structured graphical representation at an appropriate level of abstraction. In this paper the requirements for a shape transforming 4DP part are represented as the allowable variation in dimensional size and tolerance in geometric form of the functionally critical features on the part for each function that the transformed shape serves. In this paper, the authors describe how the proposed standard to specify the desired shape transformations of a 4DP part could use graphical symbols in a structured specification by means of a Transformation Control Frame (TCF) to define the rules of transforming between shapes and a Bill of Transformations (BoT) to enumerate all the Transformation Control Frames (TCF) necessary to describe the intended sequence of shape transformations. To illustrate how the graphical symbols could be applied, a SMA actuated gripper is presented as a use-case.
Abstract The term “4D Printing” (4DP) is defined as the ability for a part produced using an additive manufacture process to change its shape when activated by or exposed to one [...]
In the present contribution, we propose an effective numerical thermal modeling solution for melt pool simulations in Laser-based Powder Bed Fusion of Metals processes. The proposed model employs an anisotropic conductivity to represent melt pool dynamics effectsin a homogeneous material model. The numerical implementation of the proposed physical model is first experimentally calibrated and then validated with respect to a series of melt pool measurements as acquired by using a short-wave infrared (SWIR) camera monitoring system.
Abstract In the present contribution, we propose an effective numerical thermal modeling solution for melt pool simulations in Laser-based Powder Bed Fusion of Metals processes. The [...]
We utilized an Adjustable Ring-Mode (ARM) laser to achieve an almost fully equiaxed microstructure in powder bed Fusion-laser beam Scalmalloy®. ARM laser-built specimens exhibited over 90% fine-grained material, while circular laser-built specimens yielded less than 50% fine-grained material, using the same laser power, speed, and hatch spacing. To gain insights into these interesting results, we employed a Cellular Automata (CA) solidification simulation, incorporating the nucleation role of L12 Al3(ScxZr1-x) precipitates through a particle-based nucleation model. The simulation was coupled with the corresponding temperature field derived from finite difference analyses of the circular and ARM laser beams. The simulation results revealed a significantly thicker precipitation zone (equiaxed grains) under the ARM laser compared to the circular beam, primarily attributed to reduced temperature and cooling rates. The excellent correlation between simulation and experimental results demonstrates promising potential for the predictive application of the developed model. It can be effectively utilized to optimize heat source modulation and process parameters, thereby enabling the adaptation of microstructure and mechanical properties
Abstract We utilized an Adjustable Ring-Mode (ARM) laser to achieve an almost fully equiaxed microstructure in powder bed Fusion-laser beam Scalmalloy®. ARM laser-built specimens [...]
In this paper, the effects of geometrical imperfections observed in a lattice structure fabricated by metal 3D printer on the compressive response were investigated by using FE simulation. Geometrical imperfections which are due to excessive heat transfer and the melting of unnecessary metal powder during the fabrication process was observed using a 3D X-Ray microscope (XRM) machine. Based on the observation, two types of geometrical imperfections (strut diameter deviation and the center-axis offset) were measured, and the quantities of these imperfections on the mechanical properties of lattice block were discussed. By introducing imperfections to the FE model, a likelihood of reduced mechanical properties can be potentially adverted. In addition, by comparing the amount of geometrical imperfections, the initial stiffness and plastic collapse strength in the models based on different strut diameters, we proposed appropriate manufacturing conditions for the lattice blocks that would minimize the reduction of their mechanical properties.
Abstract In this paper, the effects of geometrical imperfections observed in a lattice structure fabricated by metal 3D printer on the compressive response were investigated by using [...]
The dynamic vibration response of sandwich beams with an anti-tetra-chiral lattice as a lightweight sandwiched core have been studied by using a nonlinear finite element analysis (FEA). Since the anti-tetra-chiral structure has a weak shear stiffness, its vibration response is strongly affected by the shear deformation. In our calculation, a 3-point bending flexural test was conducted for calculating the effective shear stiffness as well as the effective Young’s modulus of the chiral core. The natural frequency of the sandwich beam has been calculated by FEA, and predicted by using the Rayleigh-Ritz method, assuming that the sandwich beam is composed of composite continuum materials with equivalent Young’s modulus and shear modulus. Moreover, the natural frequency and damping ration of the sandwich beam produced by a 3D printer bas been measured through a vibration test, and compared with numerical results in order to clarify the effectiveness of the chiral sandwich beam as a mechanical component.
Abstract The dynamic vibration response of sandwich beams with an anti-tetra-chiral lattice as a lightweight sandwiched core have been studied by using a nonlinear finite element analysis [...]