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.
This work presents ongoing research on the influence of fluid-structure interaction (FSI) effects on the ductile crack growth in blast-loaded steel plates. Thin steel plates with Xshaped, pre-formed defects are used to allow for large, inelastic strains and ductile fracture. FSI effects were studied by comparing the numerical predictions of the uncoupled and the coupled FSI approach, where experimental data served as a backdrop to evaluate the accuracy of the numerical simulations. Numerical simulations are conducted in the EUROPLEXUS software. The clear conclusion from this study is that ductile fracture and crack propagation are influenced by FSI effects during the dynamic response of the plate. That is, the crack growth was very sensitive to the actual loading on the plate. Moreover, because the increase in CPU cost may be significant when uniformly refining the mesh, adaptive mesh refinement (AMR) was found very promising in reducing the CPU cost and maintaining the solution’s accuracy. The performance of AMR is an interesting finding in the view of numerical simulations of coarsely meshed (prior to AMR) shell structures exposed to blast loading
Abstract This work presents ongoing research on the influence of fluid-structure interaction (FSI) effects on the ductile crack growth in blast-loaded steel plates. Thin steel plates [...]
An innovative 2D axisymmetric fluid-structure interaction model of wire drawing is developed to numerically investigate the interaction between the thin lubricant film and the plastically deforming steel wire. The deformation of the wire is obtained from the linear momentum balance and the lubricant film has been calculated by the Navier-Stokes equations. Moreover, the coupling between wire and lubricant is performed by the IQN-ILS technique and a no-slip condition is imposed on the sliding fluid-structure interaction interface. In order to reduce the computational cost, a layering technique is implemented in the axially moving structure domain. This results on the one hand in monitoring the stresses and displacements of the structure and on the other hand in an observation of the hydrodynamic pressure build-up and wall shear stresses in the lubricant. Additionally, the evolution of the fluid film thickness is presented.
Abstract An innovative 2D axisymmetric fluid-structure interaction model of wire drawing is developed to numerically investigate the interaction between the thin lubricant film and [...]
G. Buron, F. Thouverez, L. Jézéquel, A. Beley, F. Thévenon
coupled2023.
Abstract
In this work, we propose iterative modal solvers to generate multiphysics finite element reduced order models. We consider the strongly coupled problems defined through differential-algebraic equations with sparse discrete operators. Piezoelectric models are common examples of such problems. The approach we propose is based on the Model Order Reduction (MOR) after Implicit Schur method [1] which is used for the Krylov subspace reduction of piezoelectric devices. While their work uses the knowledge of the loading applied to the model to generate a Krylov subspace reduction basis, we propose to build a reduction basis with a priori unknown loading by modal synthesis. The basis is built from the eigenvectors of the problem after the static condensation by Schur complement of one of the physics. Typically, the Schur complement matrix is computed explicitly and it leads to dense operators [2] which limit the problem scales that can be studied due to large memory requirements and costly computations for the eigensolver used afterward. For Krylov-based eigensolvers, the most computationally difficult step is to obtain a basis spanning the eigenspace of the problem on the considered eigenvalue range. By generalizing the MOR after Implicit Schur method, this basis can be constructed by an iterative procedure using the original sparse operators instead of the dense condensed operators. The original model may be significantly larger compared to the condensed model for typical cases. However, keeping the sparsity is a critical computational advantage for the considered problems. This method is minimally intrusive for the eigensolvers that only require the implementation of a matrix-vector product. Comparing this implicit Schur complement approach to the explicit Schur complement approach shows large computational cost reductions. It also underlines the problem scale limitations of the explicit approach even on high performance computing hardware.
Abstract In this work, we propose iterative modal solvers to generate multiphysics finite element reduced order models. We consider the strongly coupled problems defined through differential-algebraic [...]
We consider waveform iterations for dynamic coupled problems with respect to the role of time window length. We review existing theoretical results about the error of waveform iterations and the role of the time window length. Furthermore, we present numerical results for waveform iterations with both time adaptive sub solvers and with fixed time steps. This way, we are able to give a recommendation on the choice of the time window. The use of time windows can lead to an increase in efficiency. For fixed time grids, we can reliably achieve a small performance increase. For time adaptive solvers, more research is needed.
Abstract We consider waveform iterations for dynamic coupled problems with respect to the role of time window length. We review existing theoretical results about the error of waveform [...]
In a partitioned fluid-structure interaction simulation separate flow and structure solvers, each with their own spatial domain, are coupled by exchanging data on the common interface. Its computational cost is dominated by the execution of these solvers, and the cost associated with the coupling algorithm and communication are often deemed negligible. From this point of view, the computational cost is in literature typically expressed by the number of required coupling iterations per time step or equivalently the number of solver executions. However, this reasoning implicitly assumes a constant solver cost and ignores the varying number of internal subproblem iterations, i.e., solver iterations in the nonlinear solvers. This work addresses this shortcoming and shows that the computational cost of a partitioned fluid-structure interaction simulation is significantly impacted by the number of subproblem iterations performed in each solver call. Specifically, it is demonstrated that performing subproblem iterations until the solver is fully converged in each call does typically minimize the number of coupling iterations, but does not lead to minimal computational time. Instead, under the assumption of constant subproblem iteration cost, the optimum is found by minimizing a weighted sum of both coupling and subproblem iterations. The weighting factors are determined by the problem itself as well as the computer architecture.
Abstract In a partitioned fluid-structure interaction simulation separate flow and structure solvers, each with their own spatial domain, are coupled by exchanging data on the common [...]
Engineers increasingly need tools that help them automate complex simulation workflows. Besides performance, robustness and usability requirements, tools should also be easily accessible. To fulfil these requirements, Dapta Ltd is developing a cloud-based framework, which is designed to be an all-in-one solution to create, visualise, test and automate simulation workflows. Here we demonstrate the use of the dapta platform with open-source software libraries, focusing on an FSI multiphysics example.
Abstract Engineers increasingly need tools that help them automate complex simulation workflows. Besides performance, robustness and usability requirements, tools should also be easily [...]
Many multiphysics problem can be described by the coupling of several models through physical surfaces. Relying on existing model-specific solvers is very desirable, however they must be coupled in a way that ensures an accurate and stable coupled simulation. In this contribution, we present a multistep coupling scheme which relies on the history of the exchanged quantities to enable a high-order accurate coupling with time adaptation. Explicit and implicit variants are discussed in details. Numerical experiments conducted with an opensource demonstrator on a conjugate heat transfer problem show that high-order convergence is attained, and that stability is favourable compared to other classical approaches.
Abstract Many multiphysics problem can be described by the coupling of several models through physical surfaces. Relying on existing model-specific solvers is very desirable, however [...]
The intensity and frequency of natural hazards such as landslides, debris flow, and mud flows have increased significantly over the last years due to climate change and global warming. These catastrophic events are responsible for numerous destructions of infrastructures and landscapes and often even claim human lives. Therefore, in addition to the prediction, the design and installation of protective structures are of tremendous importance. In recent decades, highly flexible protective structures have been favored due to their enormous energy absorption capacity while adapting well to the environment. However, dimensioning such protective structures is a very complex task requiring advanced numerical simulation techniques. To capture the behavior of such natural hazards on the one hand and the highly flexible protection structures, including complex elements such as sliding cables or brakes on the other hand, a partitioned coupling approach is proposed in this work. This way, the most appropriate solvers, treated as black-box solvers, can be selected for each physics involved while the interaction is shifted to the shared interface.
Abstract The intensity and frequency of natural hazards such as landslides, debris flow, and mud flows have increased significantly over the last years due to climate change and global [...]
V. Singer, A. Laresse, A. Börst, R. Wüchner, K. Bletzinger
coupled2023.
Abstract
The intensity and frequency of natural hazards such as landslides, debris flow, and mud flows have increased significantly over the last years due to climate change and global warming. These catastrophic events are responsible for numerous destructions of infrastructures and landscapes and often even claim human lives. Therefore, in addition to the prediction, the design and installation of protective structures are of tremendous importance. In recent decades, highly flexible protective structures have been favored due to their enormous energy absorption capacity while adapting well to the environment. However, dimensioning such protective structures is a very complex task requiring advanced numerical simulation techniques. To capture the behavior of such natural hazards on the one hand and the highly flexible protection structures, including complex elements such as sliding cables or brakes on the other hand, a partitioned coupling approach is proposed in this work. This way, the most appropriate solvers, treated as black-box solvers, can be selected for each physics involved while the interaction is shifted to the shared interface.
Abstract The intensity and frequency of natural hazards such as landslides, debris flow, and mud flows have increased significantly over the last years due to climate change and global [...]
Physics Informed Neural Networks (PINNs) have frequently been used for the numerical approximation of Partial Differential Equations (PDEs). The goal of this paper is to construct PINNs along with a computable upper bound of the error, which is particularly relevant for model reduction of Parameterized PDEs (PPDEs). To this end, we suggest to use a weighted sum of expansion coefficients of the residual in terms of an adaptive wavelet expansion both for the loss function and an error bound. This approach is shown here for elliptic PPDEs using both the standard variational and an optimally stable ultra-weak formulation. Numerical examples show a very good quantitative effectivity of the wavelet-based error bound.
Abstract Physics Informed Neural Networks (PINNs) have frequently been used for the numerical approximation of Partial Differential Equations (PDEs). The goal of this paper is to construct [...]