Elastoplastic deformation during projectile–wall collision

The Smoothed Particle Hydrodynamics method for elastic solid deformation is modified to include von Mises plasticity with linear isotropic hardening and is then used to investigate high speed collisions of elastic and elastoplastic bodies. The Lagrangian mesh-free nature of SPH makes is very well suited to these extreme deformation problems eliminating issues relating to poor element quality at high strains that limits finite element usage for these types of problems. It demonstrates excellent numerical stability at very high strains (of more than 200%). SPH can naturally track history dependent material properties such as the cumulative plastic strain and the degree of work hardening produced by its strain history. The high speed collisions modelled here demonstrate that the method can cope easily with collisions of multiple bodies and can also naturally resolve self-collisions of bodies undergoing high levels of plastic strain. The nature and the extent of the elastic and plastic deformation of a rectangular body impacting on an elastic wall and of an elastic projectile impacting on a thin elastic wall are investigated. The final plastically deformed shapes of the projectile and wall are compared for a range of material properties and the evolution of the maximum plastic strain throughout each collision and the coefficient of restitution are used to make quantitative comparisons. Both the elastoplastic projectile–elastic wall and the elastic projectile–elastoplastic wall type collisions have two distinct plastic flow regimes that create complex relationships between the yield stress and the responses of the solid bodies.     

Modelling of metal forging using SPH

In solid metal forming processes, such as forging, large distortions in the material present challenging problems for numerical simulation using grid based methods. Computations invariably fail after some level of mesh distortion is reached unless suitable re-meshing is implemented to cope with the mesh distortion arising from the material deformation. The issue of mesh distortion and the subsequent re-meshing are topics of much research for grid based methods. These problems can be overcome by using a mesh-less numerical framework. In this paper, the application of a mesh-less method called Smoothed Particle Hydrodynamics (SPH) for modelling three-dimensional complex forging processes is demonstrated. It is shown that SPH is a useful simulation method for obtaining insights into the material deformation and flow pattern during forging of realistic industrial components. The effect of process parameters and material properties on the quality of the forged component is evaluated via SPH simulations. This includes the determination of forging force required for adequate die filling which is an important criterion for die designs. Material hardening, controlled by the degree of heat treatment, is found to have a profound effect on the material deformation pattern and the final product. Forging defects such as incomplete die filling, asymmetry in forged components, flashing and lap formation are shown to be predicted by SPH. SPH can thus potentially be used both for assessment of the quality of forged products and evaluation of prototype forging system designs.     

Modelling of the deformation processes and the localization of plastic deformations in the torsion–tension of solids of revolution

Generalized two-dimensional problems of the torsion of elastoplastic solids of revolution of arbitrary shape for large deformations under non-uniform stress-strain conditions are formulated and a method for their numerical solution is proposed. The use of this method to construct strain diagrams of materials based on experiments on the torsion of axisysmmetric samples of variable thickness until fracture occurs is described. Experimental and numerical investigations of processes of elastoplastic deformation, loss of stability and supercritical behaviour of solid cylindrical steel samples of variable thickness under conditions of monotonic kinematic loading with a torque, a tension and a combined load are presented. The mutual influence of torsion and tension on the deformation process and the limit states is estimated, and the universality (the independence of the form of the stress-strain state) of the “stress intensity – Odqvist parameter” diagram for steel for large deformations is proved.     

Modeling Abrasive Wear Caused by Small Solid Particles of Different Sizes

Abrasive wear is one of the mechanisms which cause the decrease of efficiency of hydraulic machines. The working fluid of a hydraulic machine, e.g., a turbine of a hydroelectric power plant, transports small solid particles of different sizes. Those small particles damage the surface of the hydraulic machine when contacting. In contrast to classical approaches in fluid dynamics, here, we present an approach where only mesh-free methods are applied. The Smoothed Particle Hydrodynamics (SPH) method is used for modeling the fluid in this study. The SPH method is a mesh-free method which has its advantages in describing transient fluid flows with free surfaces and large motions. The loading of the fluid consists of small solid particles of different sizes. A coupled approach for describing the loading is used. For the larger abrasive particles the Discrete Element Method and for smaller ones a transport equation is utilized. In doing so it is possible to model a loading of the fluid consisting of small particles of different sizes. The abrasive wear is described with an abrasive wear model. The wear model takes into account different parameters like the size, the velocity of the abrasive particles, and of course material parameters of both the target and the particles. On impact of an abrasive particle, the amount of removed material is stored at the boundary and in doing so the removed material over time is identified. In this work, a representative numerical example is presented. The simulations were performed with the code Pasimodo, developed at the Institute of Engineering and Computational Mechanics. It is the aim of this work to point out that it is possible to model abrasive wear due to abrasive particles with different sizes with a mesh-free approach. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

三维PTT黏弹性液滴撞击固壁面问题的改进SPH模拟

基于光滑粒子动力学(smoothed particle hydrodynamics, SPH)方法,对三维Phan-Thien Tanner(PTT)黏弹性液滴撞击固壁面问题进行了数值模拟.为了有效地防止粒子穿透固壁,且缩减三维数值模拟所消耗的计算时间,提出了一种适合三维数值模拟的改进固壁边界处理方法.为了消除张力不稳定性问题,采用一种简化的人工应力技术.应用改进SPH方法对三维PTT黏弹性液滴撞击固壁面问题进行了数值模拟,精细地捕捉了液滴在不同时刻的自由面,讨论了PTT黏弹性液滴不同于Newton(牛顿)液滴的流动特征,分析了PTT拉伸参数对液滴宽度、高度和弹性收缩比等的影响.模拟结果表明,改进SPH方法能够有效而准确地描述三维PTT黏弹性液滴撞击固壁面问题的复杂流变特性和自由面变化特征.     

Updated Lagrangian Taylor-SPH method for large deformation in dynamic problems

In this paper, the updated Lagrangian Taylor-SPH meshfree method is applied to the numerical analysis of large deformation and failure problems under dynamic conditions. The Taylor-SPH method is a meshfree collocation method developed by the authors over the past years. The governing equations, a set of first-order hyperbolic partial differential equations, are written in mixed form in terms of stress and velocity. This set of equations is first discretized in time by means of a Taylor series expansion in two steps and afterwards in space using a corrected form of the SPH method. Two sets of particles are used for the computation resulting on the elimination of the classical tensile instability. In the paper presented herein the authors propose an updated Lagrangian Taylor-SPH approach to address the large deformations of the solid, and therefore the continuous re-positioning of the particles. In order to illustrate the performance and efficiency of the proposed method, some numerical examples based on elastic and viscoplastic materials involving large deformations under dynamic conditions are solved using the proposed algorithm. Results clearly show that the updated Lagrangian Taylor-SPH method is an accurate tool to model large deformation and failure problems under dynamic loadings.     

Modeling and simulation of injection molding process of polymer melt by a robust SPH method

In the present study, the injection molding process of polymer melt based on the generalized Newtonian fluid model is investigated by a robust smoothed particle hydrodynamics (SPH) method. The numerical method is proposed by introducing a Rusanov flux into the continuity equation to improve the prediction of the pressure distribution and employing a corrected kernel gradient to improve the computational accuracy. In addition, a robust treatment of solid boundary is presented and verified by the spin-down problem. The merits of the robust SPH method are firstly illustrated by 2D dam breaking flow. Then the numerical method is extended to deal with the flow phenomena related to injection molding process of polymer melt. A number of numerical examples including 2D injection moldings of a thin plate mold, a circular disc with core, a ring-shaped channel, and a S-shaped cavity, and 3D injection moldings of a Z-shaped cavity and a four-legged fork are conducted. The numerical results are in agreement with the experiments, which demonstrate that the SPH method proposed here is capable of handling with injection molding process of polymer melt in a robust manner. Moreover, the robust SPH method allows to recover the fluctuations-free pressure and velocity fields which in most cases cannot be easily obtained by the traditional SPH method.     

Structural damage accumulation and stable postcritical deformation of composite materials

The macroscopic failure of composite materials is preceded by complex multilevel processes accompanied by accumulation and localization of damaged centers and formation of a failure cluster. Therefore, the study of these mechanisms is one of the basic problems for the mechanics of modern composite materials used in aerospace engineering. The formation of a theory of the stable postcritical deformation of the work-softening media is considered. The pseudo-plastic deformation affected by structural damage of granular composites is investigated within the framework of the considered two-level structurally phenomenological model of heterogeneous media. The stable evolution of the interconnected processes is accompanied by stress redistributions, partial or complete unloading, and strain or damage localization that are one of the main causes of implementation of the postcritical deformation stage. The numerical calculation results of inelastic deformation and failure of the periodic unidirectional fiber-reinforced composites are presented under conditions of the displacement-controlled transverse proportional loading mode. The main mechanisms of the work-softening behavior for the indicated type of materials are described in the macro-homogeneous stress-strain states. Macroscopically, the failure of heterogeneous media as a result of postcritical deformation and the loss of stability of damage accumulation depends on the stiffness of the loading system. When a deformable body is fixed on the closed surface with sufficiently but not infinitely large coefficients of stiffness, it is possible to observe the equilibrium development of the localized volumes of work-softening and damage. The constitutive equations for the work-softening isotropic, transverse isotropic, and orthotropic media are presented. The effect of the loading system on the stability of deformation, damage accumulation, and failure under monotone and nonmonotone triaxial loading was studied. The growth of failure strains with increase in stiffness of the loading system and unequal resistance of heterogeneous body are registered and investigated. A preventive unloading method is offered for the mathematical modeling of the damage accumulation during the testing of the materials on the servo-controlled systems. The displacement-controlled mode is simulated by a series of soft loading and unloading cycles. The detected phenomenon of failure where the unloading leads to stress-strain diagrams with a negative slope of the descending branch was not found either in the displacement or stress-controlled monotone loading mode.Submitted to the 10th International Conference on Mechanics of Composite Materials, April 20–23, 1998, Riga, Latvia.Perm' State Technical University, Russia. Translated from Mekhanika Kompozitnykh Materialov, Vol. 34, No. 2, pp. 234–250, March–April, 1998.     

Coupled SPH and Phase Field method for hydraulic fracturing

Simulation of fracture initiation and propagation using classical mesh-based methods involves computationally expensive operations for pre-processing and (adaptive) remeshing of complex geometries. To overcome the difficulties we propose to use a Smoothed Particle Hydrodynamics (SPH) method to model hydraulic fracturing. SPH is a meshless Lagrangian method highly suitable for large deformations [1]. The present contribution discusses a numerical approach to model fractures initiation and propagation, by coupling SPH with the Phase Field method [2]. The proposed hybrid method overcomes the instability problems that can present SPH due to kernel incompleteness. We first validate the proposed model with a stationary elastic fracture and compare the results with the classical SPH to the SPH-Phase Field approach. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

圆柱壳在径向冲击载荷作用下的弹性脉冲屈曲

当圆柱壳承受径向脉冲载荷时，如果其径厚比大于一特定值，圆柱壳将产生弹性动力屈曲．本文根据有关实验结果，假定变形模态，采用Ｌａｇｒａｎｇｅ方法分析了有限长薄圆柱壳（ａ／ｈ＝４８０）在余弦冲击载荷作用下的弹性脉冲动力屈曲．导出了动力屈曲方程组，借助数值方法求解方程，并与有关计算结果进行了比较．     

Extension of SPH to simulate non-isothermal free surface flows during the injection molding process

This article presents an extension of smoothed particle hydrodynamics (SPH) to non-isothermal free surface flows during the injection molding process. Specifically, we use the method presented by Xu and Yu, Appl. Math. Model. 48 (2017) pp. 384–409, in which the corrected kernel gradient is implemented to increase the computational accuracy and the Rusanov flux is introduced into the continuity equation to alleviate large and random pressure oscillations. To model non-isothermal free surface flows, a working SPH discretization of the temperature equation is derived. An enhanced treatment of the wall boundary is further developed, which can model arbitrary-shaped mold walls. The proposed SPH method is first validated by solving non-isothermal Couette flow and non-isothermal injection molding of a circular disc with a core and comparing the SPH results with those obtained by other numerical methods or experiments. We then extend the numerical method to non-isothermal injection molding of F-shaped and N-shaped cavities. The convergence of the method is examined with several different particle sizes. The effects of the operating conditions (e.g., injection temperature, temperature of the mold wall, and injection velocity) on the flow behavior are analyzed. All the results illustrate that the present SPH method is a powerful computational tool for simulations of non-isothermal free surface flows during the injection molding process.     

A coupled NMM-SPH method for fluid-structure interaction problems

The main challenges in the numerical simulation of fluid–structure interaction (FSI) problems include the solid fracture, the free surface fluid flow, and the interactions between the solid and the fluid. Aiming to improve the treatment of these issues, a new coupled scheme is developed in this paper. For the solid structure, the Numerical Manifold Method (NMM) is adopted, in which the solid is allowed to change from continuum to discontinuum. The Smoothed Particle Hydrodynamics (SPH) method, which is suitable for free interface flow problem, is used to model the motion of fluids. A contact algorithm is then developed to handle the interaction between NMM elements and SPH particles. Three numerical examples are tested to validate the coupled NMM-SPH method, including the hydrostatic pressure test, dam-break simulation and crack propagation of a gravity dam under hydraulic pressure. Numerical modeling results indicate that the coupled NMM-SPH method can not only simulate the interaction of the solid structure and the fluid as in conventional methods, but also can predict the failure of the solid structure.     

A MUSCL smoothed particles hydrodynamics for compressible multi‐material flows

By incorporating the Monotone Upwind Scheme of Conservation Law (MUSCL) scheme into the smoothed particles hydrodynamics (SPH) method and making use of an interparticle contact algorithm, we present a MUSCL–SPH scheme of second order for multifluid computations, which extends the Riemann‐solved‐based SPH method. The numerical tests demonstrate high accuracy and resolution of the scheme for both shocks, contact discontinuities, and rarefaction waves in the one‐dimensional shock tube problem. For the two‐dimensional cylindrical Noh and shock‐bubble interaction problems, the MUSCL–SPH scheme can resolve shocks well. Copyright © 2015 John Wiley & Sons, Ltd.     

A Lagrangian-based SPH-DEM model for fluid–solid interaction with free surface flow in two dimensions

A Lagrangian-based SPH-DEM coupling model is proposed to study fluid–solid interaction (FSI) problems with free-surface flow. In this model, SPH uses an incompressible divergence-free scheme for simulating complex flow problems. Based on the Mohr–Coulomb criterion with tension cut, the DEM describes the characteristics of solid deformation and failure by means of contact models between particles. The coupling mechanism between SPH and DEM is realised by the decoupling of the force field during the process of fluid–solid interaction. That is, the motions of fluid and solid particles are reflected by the Navier–Stokes equations and interactions among solid particles are determined by Newton's second law in the DEM. To demonstrate the applicability of the SPH-DEM model, three case studies are used to verify the different fluid interaction situations with rigid bodies, deformable objects, and granular assemblies, respectively. The results of the proposed model shows good agreement with experimental data and indicates that it is capable of capturing the features of solid movement, deformation and failure under complex flow conditions with convincing accuracy and high efficiency.     

Granular materials: Parameter identification and modelling of localisation problems

The prediction of landsliding requires an exact knowledge of the mechanical behaviour of granular materials. This kind of materials, e. g., sand, have a very complex deformation behaviour, which depend on the stress state and on the loading history. In this work, the deformation behaviour of the solid skeleton is characterised via homogeneous triaxial tests on dry sand specimens. Additionally, an appropriate elasto-plastic material law to describe the solid skeleton in the frame of Theory of Porous Media (TPM) is used, which is implemented in the FE tool PANDAS. Furthermore, a single-surface yield criterion with isotropic hardening, which limits the elastic domain, and a non-associated plastic flow are employed. The determination of the material parameters of the linear elasticity law as well as the single-surface yield criterion are based on test data of triaxial experiments. The material parameters are identified using a derivative-based optimisation method (donlp2), which is coupled with PANDAS. Finally, a simulation of a benchmark test is presented to show shear band localisation effects, where the material behaviour is described by a triphasic porous media model based on the TPM, where the constituents are a deformable solid skeleton and two pore fluids, water and air. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Contact interaction of composite shells,subjected to follower loads,with a rigid convex foundation

Based on a 7-parameter shell model, a numerical algorithm is developed for solving the contact problem for a multilayered composite shell lying on a rigid convex foundation, which is subjected to a follower pressure and undergoes arbitrarily large rotations. A new geometrically exact solid shell element is formulated, which permits one to solve the nonlinear deformation problem for thin-walled composite structures under unilateral contact constraints by using a small number of load steps. The calculation of a homogeneous ring and an angle-ply toroidal shell interacting with plane and cylindrical foundations is considered.     

Optimal Thickness of a Cylindrical Shell Under a Time-Dependent Force

We discuss thickness optimization problems for cylindrical tubes that are loaded by time-dependent applied force. This is a problem of shape optimization that leads to optimal control in linear elasticity theory. We determine the optimal thickness of a cylindrical tube by minimizing the deformation of the tube under the influence of an external force. The main difficulty is that the state equation is a hyperbolic partial differential equation of the fourth order. The first order necessary conditions for the optimal solution are derived. Based on them, a numerical method is set up and numerical examples are presented.     

Description of plastic deformation of structural materials in triaxial loading

A model of nonassociated plasticity is put forward for initially isotropic materials deforming with residual changes in volume under the action of triaxial normal stresses. The model is based on novel plastic loading and plastic potential functions, which define closed, convex, every where smooth surfaces in the 6D space of symmetric second-rank stress tensors. By way of example, the plastic deformation of a cylindrical concrete specimen wrapped with a CFRP tape and loaded in axial compression is described. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 44, No. 2, pp. 183–194, March–April, 2008.     

Change in several mechanical properties of human blood vessels with aging

A method is suggested for the determination of the initial length of a sample of a blood vessel wall at a longitudinal loading. It was found that at loading levels of up to 0.01 kg/mm², the deformation properties of the wall of a large subcutaneous vein and femoral artery differ inappreciably, but on further loading substantial differences were observed between the deformation and strength properties of the walls of vessels.     

Deformation of composites with arbitrarily oriented orthotropic fibers under matrix microdamages

In the present work, a model of nonlinear deformation of stochastic composites under microdamaging is developed for the case of a composite with orthotropic inclusions, when microdefects are accumulated in the matrix. The composite is treated as an isotropic matrix strengthened by triaxial arbitrarily oriented ellipsoidal inclusions with orthotropic symmetry of the elastic properties. It is assumed that the process of loading leads to accumulation of damage in the matrix. Fractured microvolumes are modeled by a system of randomly distributed quasispherical pores. The porosity balance equation and relations for determining the effective elastic modules in the case of orthotropic components are taken as basic relations. The fracture criterion is specified as the limiting value of the intensity of average shear stresses acting in the intact part of the material. On the basis of the analytic and numerical approach, we propose an algorithm for the determination of nonlinear deformation properties of the investigated material. The nonlinearity of composite deformations is caused by the finiteness of deformations. By using the numerical solution, the nonlinear stress–strain diagrams are predicted and discussed for an orthotropic composite material for various cases of orientation of inclusions in the matrix.