Rigid body water impact-experimental tests and numerical simulations using the SPH method

Statistics show that water impact of an aircraft in emergency is likely to have tragic consequences and therefore new researches on this topic are recommendable. In 2005, the GARTEUR AG15 was established to improve the SPH method for application to helicopter ditching. As a contribution, water impact drop tests using rigid bodies were performed at the Politecnico di Milano LAST Crash Lab to collect data and validate the numerical models. During the tests, impact decelerations were measured and suitably pressure transducers were developed to measure the impact pressures. Numerical simulations were carried out by adopting the SPH method to model the fluid region. A close experimental-numerical correlation was obtained. Findings are reported and guidelines for further investigations are proposed.     

Comparison study of MPM and SPH in modeling hypervelocity impact problems

Due to the high nonlinearities and extreme large deformation, the hypervelocity impact simulation is a challenging task for numerical methods. Meshfree particle methods, such as the smoothed particle hydrodynamics (SPH) and material point method (MPM), are promising for the simulation of hypervelocity impact problems. In this paper, the material point method is applied to the simulation of hypervelocity impact problems, and a three-dimensional MPM computer code, MPM3D, is developed. The Johnson–Cook material model and Mie–Grüneisen equation of state are implemented. Furthermore, the basic formulations of MPM are compared with SPH, and their performances are compared numerically by using MPM3D and LS-DYNA SPH module.     

Fatigue life prediction of composite laminates by FEA simulation method

This paper is to simulate the fatigue damage evolution in composite laminates and predict fatigue life of the laminates with different lay-up sequences on the basis of the fatigue characteristics of longitudinal, transverse and in-plane shear directions by finite element analysis (FEA) method. In FEA model, considering the scatter of the material’s properties, each element was assigned with different material’s properties. The stress analysis was carried out in MSC Patran/Nastran, and a modified Hashin’s failure criterion was applied to predict the failure of the elements. A new stiffness degradation model was proposed and applied in the simulation and then a strength degradation model was deduced, which is coupled with the presented stiffness degradation model. The reduced or discounted elastic constants were determined based on the failure mechanism of the laminates and the restrictive conditions of orthotropic property. The fatigue behavior and fatigue life of six kinds of E-glass/epoxy composite laminates with different lay-up sequences were experimentally studied, and the S–N curves and stiffness degradation models in longitudinal, transverse and in-plane shear direction were obtained. These fatigue data were adopted in the simulation to simulate fatigue behavior and estimate life of the laminates. The simulation results, including the fatigue life predicted and the residual stiffness, were coincident with the experimental results well except for the quasi-isotropic laminate for the lack of consideration of the out-of-plane fatigue character in the simulation.     

Orbital debris impact simulation using a parallel hybrid particle-element code

Simulations of three dimensional orbital debris impact problems, using a parallel hybrid particle-finite element code, show good agreement with experiment and good speedup in parallel computation. The simulations included single and multi-plate shields as well as aluminum and composite shielding materials, at an impact velocity of eleven kilometers per second.     

Thermal conductivity of aligned CNT/polymer composites using mesoscopic simulation

Thermal conductivity of CNT/polymer composites depends on alignment, dispersion, volume fraction and size of CNTs as well as polymer size. By coupling smoothed particle hydrodynamics and dissipative particle dynamics, thermal conductivities of random and aligned composites along with their meso morphologies are studied in detail. Thermal conductivity along the alignment of CNT can be significantly enhanced to 16 times that of polymer by increasing volume fraction, dispersion degree and length of CNT, meanwhile thermal conductivity perpendicular to the alignment of CNT is affected modestly by these factors. Enhancement of thermal conductivity of random composites could only be efficiently achieved by increasing the volume fraction of CNT. Particularly, thermal conductivity $\kappa$ is proportional to the square of volume fraction of CNT v in well dispersed random and aligned composites, i.e. $\kappa \propto v^2$.     

Meshfree particle simulation of micro channel flows with surface tension

This paper presents a study of micro channel flows using a meshfree particle approach. The approach is based on smoothed particle hydrodynamics (SPH) and its variant, adaptive smoothed particle hydrodynamics (ASPH). The incompressible flow in the micro channels is modeled as an artificially compressible flow. The surface tension is incorporated into the equations of motion. The classic Poiseuille flow and a practical micro channel flow problem of flip-chip underfill encapsulation process are investigated. It is found that the adaptive kernel can well match the computational geometry with long channels and can greatly save computational time. The simulation results are in close agreement with the analytical solutions.     

侵爆作用下混凝土靶破坏效应试验与数值模拟

基于混凝土靶体的中低速侵彻后浅埋爆炸试验,开展装药长径比为5的混凝土预留孔爆炸、侵彻后裸装药爆炸及侵彻后带壳装药爆炸数值模拟研究.对比数值模拟与试验结果,验证了数值模拟模型及参数设置的合理性.结合数值模拟结果及量纲分析方法,确定了影响侵彻后带壳装药爆炸毁伤效果的主要影响因素.对比预留孔爆炸及侵彻后爆炸毁伤效果,分析了侵...     

Smoothed particle hydrodynamics for numerical simulation of underwater explosion

 Underwater explosion arising from high explosive detonation consists of a complicated sequence of energetic processes. It is generally very difficult to simulate underwater explosion phenomena by using traditional grid-based numerical methods due to the inherent features such as large deformations, large inhomogeneities, moving interface and so on. In this paper, a meshless, Lagrangian particle method, smoothed particle hydrodynamics (SPH) is applied to simulate underwater explosion problems. As a free Lagrangian method, SPH can track the moving interface between the gas produced by the explosion and the surrounding water effectively. The meshless nature of SPH overcomes the difficulty resulted from large deformations. Some modifications are made in the SPH code to suit the needs of underwater explosion simulation in evolving the smoothing length, treating solid boundary and material interface. The work is mainly focused on the detonation of the high explosive, the interaction of the explosive gas with the surrounding water, and the propagation of the underwater shock. Comparisons of the numerical results for three examples with those from other sources are quite good. Major features of underwater explosion such as the magnitude and location of the underwater explosion shock can be well captured. Received: 2 April 2002 / Accepted: 20 September 2002     

A saturated discrete particle model and characteristic-based SPH method in granular materials

Based on the discrete particle model for solid-phase deformation of granular materials consisting of dry particulate assemblages, a discrete particle–continuum model for modelling the coupled hydro-mechanical behaviour in saturated granular materials is developed. The motion of the interstitial fluid is described by two parallel continuum schemes governed by the averaged incompressible N–S equations and Darcy's law, respectively, where the latter one can be regarded as a degraded case of the former. Owing to the merits in both Lagrangian and mesh-free characters, the characteristic-based smoothed particle hydrodynamics (SPH) method is proposed in this paper for modelling pore fluid flows relative to the deformed solid phase that is modelled as packed assemblages of interacting discrete particles. It is assumed that the formulation is Lagrangian with the co-ordinate system transferring with the movement of the solid particles. The assumed continuous fluid field is discretized into a finite set of Lagrangian (material) points with their number equal to that of solid particles situated in the computational domain. An explicit meshless scheme for granular materials with interstitial water is formulated. Numerical results illustrate the capability and performance of the present model in modelling the fluid–solid interaction and deformation in granular materials saturated with water. Copyright © 2007 John Wiley & Sons, Ltd.     

Coupled DEM-SPH simulations of wet continuous screening

During screening, a liquid stream, besides the vibration, can be applied for the acceleration of the separation. The discrete element method coupled with the smoothed particle hydrodynamics (DEM-SPH) is used to numerically analyse wet continuous screening here. Within the applied DEM-SPH a new simple model for the representation of the screening surface is suggested in this study. In this model, the influence of the screening surface on the fluid is represented using external forces, which act on the SPH particles in close vicinity of the screen. A required validation of the DEM-SPH method for the analysis of a vibrated particle-laden system is performed by comparing obtained DEM-SPH results with the results derived using the DEM coupled with finite volume method. The performed simulations of dry and wet continuous screening demonstrate that flowing water, in most simulated cases, accelerates the separation of particles. The presented study demonstrates the potential of the coupled DEM-SPH method for the analysis of wet screening processes. To our best knowledge, the simulation of wet screening using a two-way coupled numerical DEM-SPH approach not resolving the flow around individual particles is demonstrated in the scientific literature for the first time.     

A robust weakly compressible SPH method and its comparison with an incompressible SPH

This paper presents a comparative study for the weakly compressible (WCSPH) and incompressible (ISPH) smoothed particle hydrodynamics methods by providing numerical solutions for fluid flows over an airfoil and a square obstacle. Improved WCSPH and ISPH techniques are used to solve these two bluff body flow problems. It is shown that both approaches can handle complex geometries using the multiple boundary tangents (MBT) method, and eliminate particle clustering‐induced instabilities with the implementation of a particle fracture repair procedure as well as the corrected SPH discretization scheme. WCSPH and ISPH simulation results are compared and validated with those of a finite element method (FEM). The quantitative comparisons of WCSPH, ISPH and FEM results in terms of Strouhal number for the square obstacle test case, and the pressure envelope, surface traction forces, and velocity gradients on the airfoil boundaries as well as the lift and drag values for the airfoil geometry indicate that the WCSPH method with the suggested implementation produces numerical results as accurate and reliable as those of the ISPH and FEM methods. Copyright © 2011 John Wiley & Sons, Ltd.     

FEA of oblique impact tests on a motorcycle helmet

In accidents, motorcycle riders full-face helmets often make oblique impacts with road surfaces. Finite element analysis was used to predict the rotational and linear acceleration of a Hybrid II headform, representing a motorcyclist's head, in such impacts, considering the effects of friction at the head/helmet and helmet/road interfaces. Simulations of the oblique impact test in British Standard BS 6658 were validated by comparison with published data. This showed that COST 327 experimental data was largely determined by the friction coefficient (0.55) between the helmet shell and abrasive paper, and hardly affected by that between the head and helmet. Slip was predicted at the shell/paper interface throughout the impact, due to the high angular inertia of the helmet, and the normal force remaining below 3.5 kN. Simulations of more severe motorcycle helmet impacts explored the effects of impact site and direction, impact velocity components, helmet fit and the scalp. In these impacts, the higher velocity component normal to the road caused high frictional forces on the helmet shell, eventually causing it to roll on the road. The peak headform rotational accelerations, at some impact sites, were potentially injurious. The most effective method of reducing head rotational acceleration could be a reduction in the linear acceleration limit of the helmet standards.     

基于有限元与光滑粒子耦合的弹丸挤进过程分析

基于有限元（FEM）与光滑粒子（SPH）耦合算法,建立弹丸身管耦合系统动力学模型。通过对某大口径火炮弹丸挤进过程进行数值模拟,弥补弹丸挤进过程有限元分析方法无法有效模拟弹带大变形缺陷,成功模拟出弹丸挤进阶段弹带塑性流动过程,在此基础上深入研究弹带的应力应变变化规律。分析弹丸初始装填角、弹炮间隙、弹丸装填不到位等因素对挤进过程中弹丸动力响应影响。结果表明,仿真与实弹射击的弹带变形一致性较好。     

Contact algorithms for the material point method in impact and penetration simulation

The inherent no‐slip contact constraint in the standard material point method (MPM) creates a greater penetration resistance. Therefore, the standard MPM was not able to treat the problems involving impact and penetration very well. To overcome these deficiencies, two contact methods for MPM are presented and implemented in our 3D explicit MPM code, MPM3D. In MPM, the impenetrability condition may not satisfied on the redefined regular grid at the beginning of each time step, even if it has been imposed on the deformed grid at the end of last time step. The impenetrability condition between bodies is only imposed on the deformed grid in the first contact method, while it is imposed both on the deformed grid and redefined regular grid in the second contact method. Furthermore, three methods are proposed for impact and penetration simulation to determine the surface normal vectors that satisfy the collinearity conditions at the contact surface. The contact algorithms are verified by modeling the collision of two elastic rings and sphere rolling problems, and then applied to the simulation of penetration of steel ball and perforation of thick plate with a particle failure model. In the simulation of elastic ring collision, the first contact algorithm introduces significant disturbance into the total energy, but the second contact algorithm can obtain the stable solution by using much larger time step. It seems that both contact algorithms give good results for other problems, such as the sphere rolling and the projectile penetration. Copyright © 2010 John Wiley & Sons, Ltd.     

薄圆柱体高频感应温度场的测量和有限元模拟

本文利用自制的NiCr-NiAl热电偶精确地测量了薄圆柱工件端面高频感应温度场的分布;应用有限元分析软件ANSYS/Thermal对高频感应加热过程的温度场进行了模拟.热电偶测量的数据直观可信,但测取工件温度场分布费时费力;有限元分析则可较为快速获知工件各处的温度,弥补了实验工作的不足.本文研究结果显示,两种方法得出的结果具有可比性,能够相互参照.     

FEA of crack-particle interactions during delamination in interlayer toughened polymer composites

A numerical analysis was conducted, using previously obtained experimental results, to establish basic toughening mechanisms and fracture behaviour of an interlayer-toughened composite material, containing particulate modified interlayers. Aims of the analysis were to examine the influence of the particles on the plastic zone size that develops in front of the crack tip and to investigate interactions between the particles and the crack tip during elastic-plastic crack propagation. Under both loading modes particles did not promote yielding but induced micro-cracking. This was concluded to be the most dominant toughening mechanism within high particle content interlayers.     

Research and Improvement on the Accuracy of Discontinuous Smoothed Particle Hydrodynamics (DSPH) Method

Discontinuous smoothed particle hydrodynamics (DSPH) method based on traditional SPH method, which can be used to simulate discontinuous physics problems near interface or boundary. Previous works showed that DSPH method has a good application prospect [Xu et al, 2013], but further verification and improvement are demanded. In this paper, we investigate the accuracy of DSPH method by some numerical models. Moreover, to improve the accuracy of DSPH method, first order and second order multidimensional RDSPH methods are proposed by following the idea of restoring particle consistency in SPH (RSPH) method which has shown good results in the improvement of particle consistency and accuracy for non-uniform particles. This restoring particle consistency in DSPH (RDSPH) method has the advantages from both RSPH method and DSPH method. In addition, the accuracy of RDSPH methods near the interface, boundary and in non-uniform interior region are tested in one-dimensional and twodimensional spaces.     

Numerical simulation of crack deflection and penetration at an interface in a bi-material under dynamic loading by time-domain boundary element method

The hybrid time-domain boundary element method (BEM), together with the multi-region technique, is applied to simulate the dynamic process of crack deflection/ penetration at an interface in a bi-material. The whole bi-material is divided into two regions along the interface. The traditional displacement boundary integral equations (BIEs) are employed with respect to the exterior boundaries; meanwhile, the non-hypersingular traction BIEs are used with respect to the part of the crack in the matrix. Crack propagation along the interface is numerically modelled by releasing the nodes in the front of the moving crack tip and crack propagation in the matrix is modeled by adding new elements of constant length to the moving crack tip. The dynamic behaviours of the crack deflection/penetration at an interface, propagation in the matrix or along the interface and kinking out off the interface, are controlled by criteria developed from the quasi-static ones. The numerical results of the crack growth trajectory for different inclined interface and bonded strength are computed and compared with the corresponding experimental results. Agreement between numerical and experimental results implies that the present time-domain BEM can provide a simulation for the dynamic propagation and deflection of a crack in a bi-material.     

Modelling of Bird Strike on an Aircraft Wing Leading Edge Made from Fibre Metal Laminates – Part 2: Modelling of Impact with SPH Bird Model

Fibre Metal Laminates with layers of aluminium alloy and high strength glass fibre composite have been reported to possess excellent impact properties and be suitable for aircraft parts likely to be subjected to impacts such as runway debris or bird strikes. In a collaborative research project, aircraft wing leading edge structures with a glass-based FML skin have been designed, built, and subjected to bird strike tests that have been modelled with finite element analysis. In this second part of a two-part paper, a finite element model is developed for simulating the bird strike tests, using Smooth Particle Hydrodynamics (SPH) for modelling the bird and the material model developed in Part 1 of the paper for modelling the leading edge skin. The bird parameters are obtained from a system identification analysis of strikes on flat plates. Pre-test simulations correctly predicted that the bird did not penetrate the leading edge skin, and correctly forecast that one FML lay-up would deform more than the other. Post test simulations included a model of the structure supporting the test article, and the predicted loads transferred to the supporting structure were in good agreement with the experimental values. The SPH bird model showed no signs of instability and correctly modelled the break-up of the bird into particles. The rivets connecting the skin to the ribs were found to have a profound effect on the performance of the structure.     

Simulation studies on deformation behaviour of AA2017 alloy in semi-solid state using FEA

Abstract

Thixo-forming is one of the near net shaped manufacturing processes in which the final product is made between the liquidus and solidus temperatures. Forming the product at this temperature range provides large benefits compared to conventional forming methods. In the present work, deformation mechanism of AA2017 alloy in the semi-solid state has been studied. As the mechanical behaviour and deformation mechanism of semi-solid metal is completely different from that of the solid state forming, it is necessary to investigate its nature at the semi solid range. To analyse the metal flow, the stress behaviour and the corresponding strain induced, extensive finite element analysis (FEA) based simulation studies have been performed. The commercial software, DEFORM 2D was used for the simulation. The flow behaviour was modelled by formulating an Arrhenius type constitutive relation.