Simulation of the Cold Spray Particle Deposition Process
Cold spray is a rapidly developing coating technology for depositing materials in the solid state. The cold spray particle deposition process was simulated by modeling the high velocity impacts of spherical particles onto a flat substrate under various conditions. We, for the first time, proposed the Couple Eulerian Lagrangian (CEL) numerical approach to solve the high strain rate deformation problem. The capability of the CEL numerical approach in modeling the Cold Spray deposition process was verified through a systematic parameter study, including impact velocity, initial particle temperature, friction coefficient and materials combination. The simulation results by using the CEL numerical approach agree with the experimental results published in the literature. Comparing with other numerical approaches, which are Lagrangian, ALE and SPH, the CEL analyses are generally more accurate and more robust in higher deformation regimes. Besides simulating the single particle impact problem, we also extended our study into the simulation of multiple impacts. A FCC like particles arrangement model that inspired by the crystal structure was built to investigate the porosity rate and residual stress of deposited particles under various conditions. We observed not only the 3D profiles of voids, but also their distributions and developments during different procedures. Higher impact velocity and higher initial temperature of particles are both of benefit to produce a denser cold spray coating. The compressive residual stresses existed in the interface between the particle and substrate is mainly caused by the large and fast plastic deformation. Another simplified model for multiple impacts was created for the simulation of surface erosion. A severe surface erosion is the result of a high impact velocity, a high friction coefficient and a low contact angle. Two element failure models suitable for high strain rate dynamic problems were introduced in this study. For a ductile material as Copper, it followed two fracture modes in our study, which are tensile failure mode and shear failure mode. The former one mainly occurred beneath the substrate surface and the periphery of substrate craters, nevertheless the latter one was found predominately at the surface of craters. Four steps were found during the propagation of crack: void formation; crack formation; crack growth; coalescence and failure. A simple criterion equation was derived based on the simulation results for predicting the initiation of damage, which the erosion velocity v_{ero} is a function of contact angle and erosion velocity for normal impact v_{pi/2} The equivalent plastic strain could also be a parameter for identifying the onset of damage, identified as being 1.042 for Copper in our study.