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.
The Smoothed Particle Hydrodynamics (SPH) is a numerical scheme in which the domain is discretized into Lagrangian particles in the context of continuum mechanics. It has been widely used for fluid dynamics problems, and recently it has also been applied to solid mechanics with reasonable success. In this work, we present a total Lagrangian SPH (TL-SPH) for the application of solid mechanics and contact problems. Total Langrangian stands for the usage of the reference configuration to calculate the spatial derivatives. As a consequence, all particles maintain a perfect distribution, which, in turns, results in highly accurate calculations. In this way, the method is able to eliminate any problem related to tensile instability, which is one of the main shortcomings of SPH. Here, we introduce a simple, yet robust, way to include finite strain elastoplasticity into the TL-SPH method based on the logarithmic strain. Then, the elastic part can be easily defined with the Hencky elastic model, and the plastic part with any yield criteria such as Drucker-Prager. In addition, we develop a contact algorithm capable of simulating solid-solid contact problems. In this way, the simulation of different objects becomes a mix of continuous and discontinuous problems. Finally, we show the applicability of the method with several simple tests to validate the accuracy of the TL-SPH for simulating the elastoplastic solid material, as well as for contact problems including direct impact and friction effects. REFERENCES [1] Morikawa D.S. Toward robust landslide simulations from initiation to post-failure using the Smoothed Particle Hydrodynamics. Kyushu University PhD thesis (2022).
Abstract The Smoothed Particle Hydrodynamics (SPH) is a numerical scheme in which the domain is discretized into Lagrangian particles in the context of continuum mechanics. It has [...]
The Discrete Element Method (DEM) is a numerical approach that deals with the motion and interactions of individual elements. This method is mainly used in particle mechanics because it is overshadowed by other techniques such as the finite element method (FEM) when dealing with continuous problems. For this reason, it is not commonly known among structural engineers. However, its use can be found in cases that combine particle and continuum mechanics problems, such as crack propagation in reinforced concrete members. Calculating and optimizing these types of problems using FEM is challenging due to frequent mesh changes or the requirement for difficult-to-detect parameters in standard practice. This post discusses the possibility of using the DEM method with an extension of the beam-bound model (BBM). In this method, a beam element is inserted between the bounded discrete elements to transfer all types of forces and moments. The problem of this method is to define the correct cross-sectional and material characteristics of the individual beam members. In study, we focus on the determination of these parameters and the verification of this method on model examples of reinforced concrete elements such as simply supported and fixed beams.
Abstract The Discrete Element Method (DEM) is a numerical approach that deals with the motion and interactions of individual elements. This method is mainly used in particle mechanics [...]
The identification of material parameters occurring in material models is essential for structural health monitoring. Due to chemical and physical processes, building structures and materials age during their service life. This, in turn, leads to a deterioration in both the reliability and quality of the structures. The material parameters indicate possible damage and material degradation, as they directly reflect the resistance of the structure to external impacts. We further developed physics-informed neural networks (PINNs) [1] for the calibration of the linear-elastic material model from full-field displacement data and global force data in a realistic regime [2]. For a realistic data regime, the optimization problem had to be conditioned. The advantage of this method is a straightforward inclusion of observation data. Unlike grid-based methods, such as the least square finite element method approach, no computational grid and no interpolation of the data are required. However, directly solving inverse problems using PINNs is computationally expensive and prone to realistic noise levels in the measurement data. Moreover, the PINN must be trained completely from scratch for each new full-field displacement measurement, even if the geometry and material of the structure remain unchanged. In our ongoing work, we are therefore focusing on learning parameterized solutions of parametric partial differential equations using PINNs, such as [3]. We further investigate the ability of parametric PINNs to act as a surrogate for the identification of material parameters from full-field displacement data. By learning parameterized solutions, the PINN does not need to be completely re-trained for each full-field displacement measurement. The calibration of the material model can thus be drastically accelerated, and information about the material condition can be provided near real-time. Furthermore, we also plan to apply the parametric PINN to more complex material models, such as those for hyper-elastic and elasto-plastic materials.
Abstract The identification of material parameters occurring in material models is essential for structural health monitoring. Due to chemical and physical processes, building structures [...]
The collision of ships remains a significant cause of accidents, resulting in severe environmental consequences such as oil spillage from oil tankers. Enhancing the crashworthiness of ship structural design is crucial, and one approach being explored is the filling of the double hull structure with granular material [1]. Coating these particles with environmentally friendly materials can optimize their ability to absorb kinetic energy and transfer loads from the outer to the inner hull [2]. However, the mechanical behavior of the coated particles depends on the type of coating material, posing challenges in developing numerical simulation models. To address this, an open-source Discrete Element Method (DEM) code called MUSEN [3] is utilized to numerically model the behavior of coated particles. Furthermore, MUSEN can be extended using the Bonded Particle Method (BPM) to simulate particle breakage through solid bridges. However, the inclusion of coating material in the model increases the number of parameters, as well as computational time and cost. This requires a robust methodology to characterize the mechanical behaviour regardless of the type of coating material with reduced computational cost, keeping in mind the actual application inside the ship double hull. Sensitivity analyses and parametric studies are conducted to understand the effect of input parameters and identify the influential ones. Subsequently, algorithms such as the Particle Swarm Algorithm are employed for parameter optimization. Finally, models of different fidelity are used to compare the results from multi particle compression tests. The findings from these simulations will be presented in this contribution
Abstract The collision of ships remains a significant cause of accidents, resulting in severe environmental consequences such as oil spillage from oil tankers. Enhancing the crashworthiness [...]
Several strategies have been proposed for high performance computing with DEM. Yet, quite few works have been published on this topic (e.g., [1-3]), compared to the extensive use of DE simulations as a modelling tool in research and industry. This is however not surprising, since granular media at the grain scale appear as highly evolutive and disordered systems, which makes effective parallelisation a challenging task. We present an approach which exploits an inherent limitation of DE simulations, as the opportunity for an efficient and flexible parallel implementation on distributed memory computers. Namely, the “numerical sound speed” set by the representative particle diameter and the time step defines an upper bound for the celerity at which perturbations can propagate across the discrete medium. We show how this numerical artefact (of negligible importance for most applications) can be turned into an original criterion of spatial domain decomposition, which leads to a DEM-suited parallelisation scheme. We present our approach through its actual implementation into the DEM’ocritus code [4,5]. We analyse its performance through benchmarks and parametric analyses of biaxial tests, on assemblies of up to 15 million circular particles.
Abstract Several strategies have been proposed for high performance computing with DEM. Yet, quite few works have been published on this topic (e.g., [1-3]), compared to the extensive [...]
M. Harutyunyan, S. Emmerich, S. Steidel, M. Burger, K. Jareteg, J. Quist
particles2023.
Abstract
Modeling of soil-tool interaction for industrial applications involves the coupling of soil models, mostly based on the Discrete Element Method (DEM), with multibody systems, representing construction machinery such as excavators or wheel loaders. To obtain accurate predictions of reaction forces on tools like wheel loader buckets, it is essential to have an appropriate parametrization procedure, which makes use of data obtained from laboratory tests such as the triaxial compression test or direct shear test for different types of soil. Simulations with suitable DEM-softwares can then be validated against the experimental data to assess the applicability and performance of the numerical methods [2, 3] The DEM software GRAnular Physics Engine (GRAPE) developed at Fraunhofer ITWM in Germany has been successfully used to simulate compression and shear tests and has been proven to yield good predictions of the soil-tool interaction and draft forces [4] for spherical particles such as coarse sand. A year-long collaboration with the Fraunhofer-Chalmers Centre (FCC) in Sweden has successfully resulted in the development of a soil simulation toolbox for FCC’s general purpose DEM solver Demify® incorporating and enhancing the simulation techniques used in GRAPE [1]. Co-simulation is enabled through an FMI interface coupling Demify® with multibody systems modeled e.g., in Simulink to evaluate relevant variables such as forces at critical linkage joints of the construction machines. To model real-life application scenarios, the entire workflow must be considered, from the soil parametrization process, to setting the particulate soil model and performing numerical simulations for the specific problem, to the final post-processing to analyze the load data. We will characterize and discuss the different steps of the workflow and present simulation results obtained with the toolbox Demify® for Heavy Machinery.
Abstract Modeling of soil-tool interaction for industrial applications involves the coupling of soil models, mostly based on the Discrete Element Method (DEM), with multibody systems, [...]
Applications requiring materials with layer-to-layer strength, from basketry to thermal protection systems, use interlaced, three-dimensional woven materials. NASA is developing and deploying woven material heat shields for missions, including Artemis I (3D-MAT for compression pads) and potentially Mars Sample Return - Earth Entry System. These materials are complex, hierarchical and must protect from extreme environments and phenomena, such as deformation, impact and high-enthalpy heating. Woven material performance depends on microstructure, damage and weave geometry. Therefore, fiber-specific models are needed to simulate fiber contacts within the weave hierarchical geometry (fiber to tow to yarn to weave) and the inherent directionality of fibers. Explicit fiber models can simulate how weave microstructure evolution affects thermal and mechanical properties. We parameterize a discrete element bonded particle models (DEM-BPM) of fibers to capture thermal and mechanical behavior within and between fibers, with bonded and contact forces, respectively. We study the proportion of heat transfer and stress via the contact network, fiber bonds and the weave geometry, for example, with respect to yarn warp-weft identity (whether it interlaces weave layers). Our results demonstrate the importance of explicit fiber modeling for connecting microstructure with thermal and mechanical properties.
Abstract Applications requiring materials with layer-to-layer strength, from basketry to thermal protection systems, use interlaced, three-dimensional woven materials. NASA is developing [...]
Particle-In-Cell (PIC) methods such as the Material Point Method (MPM) can be cast in formulations suitable to the requirements of data locality and fine-grained parallelism of modern hardware accelerators such as Graphics Processing Units (GPUs). While continuum mechanics simulations have already shown the capabilities of MPM on a wide range of phenomena, the use of the method in compressible gas dynamics is less frequent. This contribution aims to show the potential of a GPU-based MPM parallel implementation for compressible fluid dynamics, as well as to assess the reliability of this approach in reproducing supersonic gas flows against solid obstacles. The results in the paper represent a stepping stone towards a highly parallel, Multi-GPU, MPM-base solver for M ach > 1 Fluid-Structure Interaction problems.
Abstract Particle-In-Cell (PIC) methods such as the Material Point Method (MPM) can be cast in formulations suitable to the requirements of data locality and fine-grained parallelism [...]
The Material Point Method (MPM) is widely used for challenging applications in engineering, and animation but lags behind some other methods in terms of error analysis and computable error estimates. The complexity and nonlinearity of the equations solved by the method and its reliance both on a mesh and on moving particles makes error estimation challenging. Some preliminary error analysis of a simple MPM method has shown the global error to be first order in space and time for a widely-used variant of the Material Point Method. The overall time dependent nature of MPM also complicates matters as both space and time errors and their evolution must be considered thus leading to the use of explicit error transport equations. The preliminary use of an error estimator based on this transport approach has yielded promising results in the 1D case. One other source of error in MPM is the grid-crossing error that can be problematic for large deformations leading to large errors that are identified by the error estimator used. The extension of the error estimation approach to two space higher dimensions is considered and together with additional algorithmic and theoretical results, shown to give promising results in preliminary computational experiments.
Abstract The Material Point Method (MPM) is widely used for challenging applications in engineering, and animation but lags behind some other methods in terms of error analysis and [...]