Academic literature on the topic 'Coupling SPH-FE'

Purpose – The purpose of this paper is to confirm that the axisymmetric finite element and smoothed particle hydrodynamics (FE-SPH) adaptive coupling method is effective to solve explosion problem in concrete based on the experiments. Design/methodology/approach – Axisymmetric FE-SPH adaptive coupling method is first presented to simulate dynamic deformation process of concrete under internal blast loading. Using calculation codes of FE-SPH coupling method, numerical model of explosion is approximated initially by finite element method (FEM), and distorted finite elements are automatically converted into meshless particles to simulate damage, splash of concrete by SPH method, when equivalent plastic strain of elements reaches a specified value. Findings – In this paper, damage process and pressure curve of concrete around explosive are analyzed and buried depth of explosive in concrete influence on damage effect under internal blast loading are obtained. Numerical analyses show that FE-SPH coupling method integrates high computational efficiency of FEM and advantages of SPH method, such as natural simulation to damage, splash and other characteristics of explosion in concrete. Originality/value – This work shows that FE-SPH coupling method has good performance to solve the explosion problem.

Le problème de l’hydroplannage a fait l’objet de peu de travaux de simulation jusqu’à présent du fait de sa complexité : couplage fluide-structure, complexité de la structure du pneu du fait des matériaux en présence, contact avec l’asphalte, complexité de l’écoulement fluide résultant (interface extrêmement complexe,assèchement de la route, ventilation, développement éventuel de la turbulence et de cavitation). Dans ce contexte, Michelin, Centrale Nantes et NextFlowSoftware ont cherché récemment à évaluer la capacité du solveur SPH développé par ces deux derniers pour classifier des pneumatiques en fonction de la géométrie de leurs structures surfaciques, sans prendre en compte la phase gazeuse. Cela a permis de démontrer la faisabilité de telles simulations par méthode SPH, et même d’obtenir de bons résultats avec pour avantages de s’absoudre des difficultés liées au maillage. L’autre avantage conséquent d’utiliser la méthode SPH pour modéliser le fluide dans ce contexte est apparu dans sa capacité à se coupler relativement aisément à des solveurs classiques de type Eléments Finis pour le problème structurel. L’objectif du doctorat est triple, poursuivre la qualification du couplage SPH–Eléments Finis, en particulier en termes énergétiques, développer des schémas permettant d’assurer un bon compromis stabilité / précision / temps de calcul. Deuxièmement quantifier l’influence des différents phénomènes physiques en jeu pour déterminer lesquels doivent être modélisés. Enfin adapter des modélisations SPH permettant de prendre en compte simultanément les différents phénomènes influant pour réaliser des simulations du problème complet
The aquaplaning problem has been the topic of simulation works emphasizing its complexity: fluid structure interactions, structures modelling, materials involved, contact with asphalt and the complexity of the resulting fluid flow (extremely complex interface, drying up the road, ventilation, possible development of turbulence and cavitation). As additional difficulty, the tire is a highly deformable body and fluid-structure interaction effects should be considered, leading to a challenging problem for the numerical modelling. Then Michelin, Ecole Centrale Nantes and NextFlow Software have recently tested the ability of the SPH solver developed by the two latter to classify tires based on their surface structure geometries, without considering the gas phase. In this context, the interest of SPH for modelling efficiently the aquaplaning flow has been underlined. The meshless and Lagrangian feature of SPH naturally avoid the problem of fluid/solid grid compatibility. The other significant advantage of the SPH method, in this context, appears in its ability to be relatively easily coupled to with conventional Finite Element solvers. The aim of this workis three fold. First, qualify the SPH-FE coupling strategy, especially in terms of energy and then develop schemes to ensure a good compromise among stability, accuracy and computation time. Second, quantify the influence of different involved physical phenomena to determine which should be modelled. Finally, adapt SPH models to simultaneously consider different phenomena and to performe simulations of the complete problem

Dans le cadre de ce travail, une technique non-intrusive est proposée pour coupler la méthode Smoothed Particle Hydrodynamics (SPH) à la méthode des Eléments Finis afin de résoudre numériquement des problèmes dynamiques et non-linéaires d’interaction fluide-structure en permettant l’utilisation des pas de temps différents dans les deux domaines de calcul (fluide et solide). Ces développements sont motivés par le besoin de simuler numériquement des phénomènes rapides et très non-linéaires qui prennent en compte des impacts en se servant des intégrateurs temporels explicites dans chaque sous-domaine de calcul (Newmark explicite pour le solide et Runge-Kutta 2 pour le fluide). De ce fait, le pas de temps de stabilité est limité par des caractéristiques intrinsèques au modèle numérique du phénomène étudié et en conséquence, il devient important de pouvoir intégrer chaque sous-domaine numérique avec un pas de temps proche de son pas de temps de stabilité. Pour permettre d’utiliser un pas de temps proche du pas de temps de stabilité pour chaque sous-domaine, des méthodes de décomposition de domaines dual-Schur sont implémentées et validées pour des cas en 1-D, 2-D, et 3-D. Des simulations numériques d’impacts de cailloux sur des aubes des turbines hydrauliques sont aussi effectue´es afin de prédire le dommage que cet évènement peut engendrer
A method to couple smoothed particle hydrodynamics and finite elements methods for nonlinear transient fluid–structure interaction simulations by adopting different time-steps depending on the fluid or solid sub-domains is proposed. These developments were motivated by the need to simulate highly non-linear and sudden phenomena that take into acount solid impacts and hence require the use of explicit time integrators on both sub-domains (explicit Newmark for the solid and Runge–Kutta 2 for the fluid). However, due to critical time-step required for the stability of the explicit time integrators in, it becomes important to be able to integrate each sub-domain with a different time-step while respecting the features that a previously developed mono time-step coupling algorithm offered. For this matter, a dual-Schur decomposition method originally proposed for structural dynamics was considered, allowing to couple time integrators of the Newmark family with different time-steps with the use of Lagrange multipliers

L’adhérence mouillée des pneumatiques est une performance essentielle touchant à la sécurité des passagers. Dans cettesituation, le contact pneu/sol devient plus complexe à comprendre et modéliser, faisant intervenir des mécanismes physiques non triviaux tels que le couplage fluide-structure et les écoulements turbulents. Dans le but d’améliorer la compréhension de l’hydroplanage des pneumatiques, la présente thèse vise à mettre en place un plan de comparaisonentre les simulations numériques couplées SPH- Éléments finis et les résultats de tests r-PIV. En effet, la méthode SPH présente de nombreux avantages de par sa nature Lagrangienne et sans maillage pour modéliser le domaine fluide. De plus, son couplage avec la méthode des Éléments finis est relativement aisé. Par ailleurs, sur le plan expérimental, la r-PIV a été introduite récemment pour étudier le roulage d’un pneumatique sur une flaque d’eau. Cette nouvelle approche constitue un puissant outil pour valider les simulations numériques sur la base de comparaisons locales de la circulation de l’eau pour une sculpture donnée. Enfin, les simulations numériques constituent également un moyen d’évaluation de la r-PIV grâce à une vision 3D du phénomène et à l’accès à des données encore inaccessibles expérimentalement
The wet grip performance of tires is an essential criterion affecting the safety of passengers. In this situation, the tire/ground contact becomes more complex to understand and model, involving non-trivial physical mechanisms such as fluid-structure coupling and turbulent flows. In the vision of improving our understanding of tires’ hydroplaning, this thesis aims to set up a comparison strategy between the SPH-Finite Elements coupled numerical simulations and the r-PIV testresults. Indeed, the SPH method has many advantages due to its Lagrangian and meshless nature to model the fluid part. Moreover, its coupling with the finite element method is relatively easy. In addition, the r-PIV was recently introduced for experimental investigations of a tire rolling over a water puddle. This new approach performed effectively as a powerfultool for validating numerical simulations based on local comparisons of the water circulation for a given tire tread. Finally, numerical simulations also evaluate r-PIV thanks to a 3D vision of the phenomenon and access to data that are still inaccessible experimentally

Abstract Smoothed Particle Hydrodynamics (SPH), a particle-based, meshless method originally developed for modeling astrophysical problems, is being increasingly used for modeling fluid mechanics and solid mechanics problems. Due to its advantages over grid-based methods in the handling of large deformations and crack formation, the method is increasingly being applied to model material removal processes. However, SPH method is computationally expensive. One way to reduce the computational time is to partition the domain into two parts where, the SPH method is used in one segment undergoing large deformations and material separation and in the second segment, the conventional finite element (FE) mesh is used. In this work, the accuracy of this SPH-FEM approach is investigated in the context of orthogonal cutting. The high deformation zone (where chips form and curl) is meshed with the SPH method, while the rest of the workpiece is modeled using the FE method. At the interface, SPH particles are coupled with FE mesh for smooth transfer of stress and displacement. The boundary conditions are applied to tool and FE zone of the workpiece. For comparison purposes, a fully-SPH model (workpiece fully discretized by SPH) is also developed. This is followed by a comparison of the results from the coupled SPH-FE model with the SPH model. A comparison of the chip profile, the cutting force, the von Mises stress and the damage parameter show that the coupled SPH-FE model reproduces the SPH model results accurately. However, the SPH-FE model takes almost 40% less time to run, a significant gain over the SPH model. Similar reduction in computation time is observed for in a micro-cutting application (depth of cut of 300 nm). Based on these results, it is concluded that coupling SPH with FEM in machining models decreases simulation time significantly while still producing accurate results. This observation suggests that three-dimensional machining problems can be modeled using the combined SPH-FEM approach without sacrificing accuracies.

Wave impact or slamming is a phenomenon characterized by large local pressures (10 bar or more) for very short durations (order of milliseconds). Slamming loads can cause severe damage to the structure [1]. Different numerical approximation methods are available for simulating the fluid structure interaction problems. Traditional mesh techniques use nodes and elements for approximating the continuum equations whereas particle methods like smoothed particle hydrodynamics (SPH) approximates the continuum equations using the kernel approximation technique and hence can be used for a wide range of fluid dynamics problems [2]. Since composite materials are finding increased application in the ship building industry because of their low weight and high strength properties, it is important to understand the effect of slamming loads on composite structures [3]. Normally, composite structures are made quasi-rigid to resist slamming loads, but inducing some deformability helps in reducing the incident pressures and at the same time reduces the overall weight of the structure and in turn the material cost. On the other side, inducing deformability effects the durability of the structure. In this paper, the effect of slamming on two-dimensional cylindrical structures is studied using three solvers i.e., 1) SPH solver, 2) Explicit solver and 3) Implicit solver. In the case of SPH solver, water is modelled using SPH particles and cylinder is modelled using finite elements (FE), in this case shell elements. A coupling between the SPH and FE solvers is made to simulate the fluid-structure interactions. Contact is modelled using the contact algorithms. In the case of the explicit solver, water is modelled using hexahedron or brick elements with one element in the thickness direction since symmetry is applicable along the thickness of the cylinder. Shell elements are used for modelling the cylinder and contact is handled using node to surface contact algorithm. In the case of the implicit solver, water is represented by pure two-dimensional elements. Quadratic elements are used to represent the continuum around the cylinder and triangular elements are used to represent the far off field and also to control the mesh movement. Line elements are used to represent the cylinder in this case. Two models are tested in all the three solvers: 1) rigid cylinder and 2) deformable cylinder. A comparative study of these three solvers shows that the implicit solver needed more calculation time compared to other solvers. The SPH method required less particles than the number of nodes in the other two methods to converge on the peak pressure. All three solvers show reduction of peak pressure in case of the deformable cylinder.

In this study, we present results from a numerical model of a full-scale fracture propagation test where the pipe sections are filled with impure, dense liquid-phase carbon dioxide. All the pipe sections had a 24″ outer diameter and a diameter/thickness ratio of ∼32. A near symmetric telescopic set-up with increasing toughness in the West and East directions was applied. Due to the near symmetric conditions in both set-up and results, only the East direction is modelled in the numerical study. The numerical model is built in the framework of the commercial finite element (FE) software LS-DYNA. The fluid dynamics is solved using an in-house computational fluid dynamics (CFD) solver which is coupled with the FE solver through a user-defined loading subroutine. As part of the coupling scheme, the FE model sends the crack opening profile to the CFD solver which returns the pressure from the fluid. The pipeline is discretized by shell elements, while the backfill is represented by the smoothed-particle hydrodynamics (SPH) method. The steel pipe is described by the J2 constitutive model and an energy-based fracture criterion, while the Mohr-Coulomb material model is applied for the backfill material. The CFD solver applies a one-dimensional homogeneous equilibrium model where the thermodynamic properties of the CO2 are represented by the Peng-Robinson equation-of-state (EOS). The results from the simulations in terms of crack velocity and pressure agree well with the experimental data for the low and medium toughness pipe sections, while a conservative prediction is given for the high-toughness section. Further work for strengthening the reliability of the model to predict the arrest vs. no-arrest boundary of a running ductile fracture is addressed.