Parallelized Implementation of the Finite Particle Method for Explicit Dynamics in GPU

As a novel kind of particle method for explicit dynamics,the finite particle method(FPM)does not require the formation or solution of global matrices,and the evaluations of the element equivalent forces and particle displacements are decoupled in nature,thus making this method suitable for parallelization.The FPM also requires an acceleration strategy to overcome the heavy computational burden of its explicit framework for time-dependent dynamic analysis.To this end,a GPU-accelerated parallel strategy for the FPM is proposed in this paper.By taking advantage of the independence of each step of the FPM workflow,a generic parallelized computational framework for multiple types of analysis is established.Using the Compute Unified Device Architecture(CUDA),the GPU implementations of the main tasks of the FPM,such as evaluating and assembling the element equivalent forces and solving the kinematic equations for particles,are elaborated through careful thread management and memory optimization.Performance tests show that speedup ratios of 8,25 and 48 are achieved for beams,hexahedral solids and triangular shells,respectively.For examples consisting of explicit dynamic analyses of shells and solids,comparisons with Abaqus using 1 to 8 CPU cores validate the accuracy of the results and demonstrate a maximum speed improvement of a factor of 11.2.

OPTIMIZATION OF AUTOBODY PANEL STAMPING PROCESS BASED ON DYNAMIC EXPLICIT FINITE ELEMENT METHOD

Taking CPU time cost and analysis accuracy into account, dynamic explicit finite ele- ment method is adopted to optimize the forming process of autobody panels that often have large sizes and complex geometry. In this paper, for the sake of illustrating in detail how dynamic explicit finite element method is applied to the numerical simulation of the autobody panel forming process,an example of optimization of stamping process pain meters of an inner door panel is presented. Using dynamic explicit finite element code Ls-DYNA3D, the inner door panel has been optimized by adapting pa- rameters such as the initial blank geometry and position, blank-holder forces and the location of drawbeads, and satisfied results are obtained.

Study on Contact Algorithm of Dynamic Explicit FEM for Sheet Forming Simulation

Based on existing algorithms, a newly developed contact search algorithm is proposed. The new algorithm consists of global search, local searching, local tracking and penetration calculation processes. It requires no iteration steps. It can deal with not only general tool surfaces with vertical walls, but also tool surfaces meshed with elements having very poor aspect ratios. It is demonstrated that the FE code employing this new contact search algorithm becomes more reliable, efficient and accurate for sheet metal forming simulation than conventional ones.

A nonlinear explicit dynamic GBT formulation for modeling impact response of thin-walled steel members

A nonlinear explicit dynamic finite element formulation based on the generalized beam theory(GBT)is proposed and developed to simulate the dynamic responses of prismatic thin-walled steel members under transverse impulsive loads.Considering the rate strengthening and thermal softening effects on member impact behavior,a modified Cowper-Symonds model for constructional steels is utilized.The element displacement field is built upon the superposition of GBT cross-section deformation modes,so arbitrary deformations such as cross-section distortions,local buckling and warping shear can all be involved by the proposed model.The amplitude function of each cross-section deformation mode is approximated by the cubic non-uniform B-spline basis functions.The Kirchhoff s thin-plate assumption is utilized in the construction of the bending related displacements.The Green-Lagrange strain tensor and the second Piola-Kirchhoff(PK2)stress tensor are employed to measure deformations and stresses at any material point,where stresses are assumed to be in plane-stress state.In order to verify the effectiveness of the proposed GBT model,three numerical cases involving impulsive loading of the thin-walled parts are given.The GBT results are compared with those of the Ls-Dyna shell finite element.It is shown that the proposed model and the shell finite element analysis has equivalent accuracy in displacement and stress.Moreover,the proposed model is much more computationally efficient and structurally clearer than the shell finite elements.

Sheet Metal Forming Simulation and Its Application in Stamping Process of Automobile Panels

The paper starts with a brief overview to the necessity of sheet metal forming simulation and the complexity of automobile panel forming, then leads to finite element analysis (FEA) which is a powerful simulation tool for analyzing complex three-dimensional sheet metal forming problems. The theory and features of the dynamic explicit finite element methods are introduced and the available various commercial finite element method codes used for sheet metal forming simulation in the world are discussed,and the civil and international status quo of automobile panel simulation as well. The front door outer panel of one certain new automobile is regarded as one example that the dynamic explicit FEM code Dynaform is used for the simulation of the front door outer panel forming process. Process defects such as ruptures are predicted. The improving methods can be given according to the simulation results. Foreground of sheet metal forming simulation is outlined.

Role of Friction in Cold Ring Rolling

Cold ring rolling is an advanced but complex metal forming process under coupled effects with multi-factors, such as geometry sizes of rolls and ring blank, material, forming process parameters and friction, etc. Among these factors, friction between rolls and ring blank plays an important role in keeping the stable forming of cold ring rolling. An analytical method was firstly presented for proximately determining the critical friction coefficient of stable forming and then a method was proposed to determine the critical friction coefficient by combining analytical method with numerical simulation. And the influence of friction coefficient on the quality of end-plane and side spread of ring,rolling force, rolling moment and metal flow characteristic in the cold ring rolling process have been explored using the three dimensional （3D） numerical simulation based on the elastic-plastic dynamic finite element method （FEM） under the ABAQUS software environment, and the results show that increasing the friction on the contact surfaces between rolls and ring blank is useful not only for improving the stability of cold ring rolling but also for improving the geometry and dimension precision of deformed ring.

Explicit solvent molecular dynamics simulations of chaperonin-assisted rhodanese folding

Chaperonins are known to facilitate the productive folding of numerous misfolded proteins. Despite their established importance, the mechanism of chaperonin-assisted protein folding remains unknown. In the present article, all-atom explicit solvent molecular dynamics （MD） simulations have been performed for the first time on rhodanese folding in a series of cavity-size and cavity-charge chaperonin mutants. A compromise between stability and flexibility of chaperonin structure during the substrate folding has been observed and the key factors affecting this dynamic process are discussed.

Transient dynamic analysis and comparison on wedging processes of overrunning clutches with different contact surfaces

The effects of contact surface on dynamic wedging behavior of the roller and inner-ring of the overrunning clutch in a dual-turbine torque converter were investigated to reveal the friction self- locking mechanism and dynamic process. Planar strain clutch models including roller, inner-ring and outer-ring were built, and transient wedging process was analyzed with an explicit dynamics meth- od. The modeling of stress and strain distribution and variation of two kinds of contact surfaces show that there are three stages named slipping, wedging and binding respectively during whole wed- ging process. Meanwhile the geometric structures of contact surfaces greatly influence the peak stress and strain distribution of the wedging process of the roller and inner-ring. The load bearing performance of contact surfaces with logarithmic spiral curve is better than that with straight line. Our study provides theoretical foundation for design and further optimization of wedging contact surface of an overrunning clutch in a dual-turbine hydrodynamic torque converter.

Some aspects affecting the response of RC compartment structures subjected to internal blast loads

Protective compartments are typically used to protect some specific structures from internal explosions,such as industrial buildings that contain devices that may explode in certain circumstances.This research investigates how the response of reinforced concrete(RC)compartment structures subjected to internal blast loads are affected by the following aspects:introduction of material nonlinearity in the analysis,reinforcement ratio,and aspect ratio of the compartment.To achieve this goal,a calibrated and sophisticated FE numerical model is introduced,and a parametric study for the intended aspects is carried out.A discussion of the results and conclusions are offered,which show the role of each aspect in the dynamic performance of the compartment structures.The main conclusions are as follows:introduction of material nonlinearity in this type of analysis and for these structures is very important and significant in obtaining accurate outputs that are similar to actual behavior;the reinforcement ratio has a significant effect on the response and its effect varies depending on the thickness of the compartment;in general,increasing the reinforcement ratio enhances the behavior and reduces the stresses in the compartment;and the aspect ratio of the compartment does not show a clear pattern on the response of such structures under internal blast loads.

Investigation of Bearing Fault Diagnosis Based on Explicit Dynamics Analysis in ANSYS/LS-DYNA

In this paper, deep grove ball bearing GB6206 has been chosen as research object, and the explicit dynamics analysis method in ANSYS/LS-DYNA has been used to study features of fault bearing which with a tiny pit in the inner ring raceway. In the process of building this bearing FEM, the following parameters have been well considered, such as boundary conditions, friction, contaction, loads and so on. Through simulation, the corresponding equivalent stress nephograms and acceleration of nodes on the inner ring raceway has been obtained. According to features of acceleration which occurs neighbor to fault pit, bearing＇s fault diagnosis has been realized. This paper provides a new way in monitoring bearing status and diagnosing fault of bearing.

Scalar Susceptibility of QCD from Dyson-Schwinger Approach

In quantum chromodynamics （QCD）, the scalar susceptibility represents the modification of the quark condensate, to a small perturbation of the parameter responsible for the explicit breaking of the symmetry, i.e., the current quark mass. By studying the linear response of the dressed quark propagator to the presence of a nonzero quark mass, we derive a model-independent formula for the scalar susceptibility, which contains the dressed quark propagator G（p） and the dressed scalar vertex F（p, 0）. The numerical values of the scalar susceptibility Xs are calculated within the framework of the rainbow-ladder approximation of the Dyson-Schwinger approach by employing two typical forms of model gluon propagator.

Quark Loop Effects on Dressed Gluon Propagator in Framework of Global Color Symmetry Model

Based on the global color symmetry mode/ （GCM）, a method for obtaining the quark loop effects on the dressed gluon propagator in GCM is developed. In the chiral limit, it is found that the dressed gluon propagator containing the quark loop effects in the Nambu-Goldstone and Wigner phases are quite different. In solving the quark self-energy functions in the two different phases and subsequent study of bag constant one shouM use the above dressed gluon propagator as input. The above approach for obtaining the current quark mass effects on the dressed gluon propagator is quite general and can also be used to calculate the chemical potential dependence of the dressed gluon propagator.

Mechanical Design of Long-Term Body-Adhered Medical Devices to Maximize On-Body Survival

Long-term, body-adhered medical devices rely on an adhesive interface to maintain contact with the patient. The greatest threat to on-body adhesion is mechanical stress imparted on the medical device. Several factors contribute to the ability of the device to withstand such stresses, such as the mechanical design, shape, and size of the device. This analysis investigates the impact that design changes to the device have on the stress and strain experienced by the system when acted on by a stressor. The analysis also identifies the design changes that are most effective at reducing the stress and strain. An explicit dynamic finite element analysis method was used to simulate several design iterations and a regression analysis was performed to quantify the relationship between design and resultant stress and strain. The shape, height, size, and taper of the medical device were modified, and the results indicate that, to reduce stress and strain in the system, the device should resemble a square in shape, be short in height, and small in size with a large taper. The square shape experienced 17.5% less stress compared to the next best performing shape. A 10% reduction in device height resulted in a 21% reduction in stress and 24% reduction in strain. A 20% reduction in device size caused a 7% reduction in stress and 2% reduction in strain. A 20% increase in device taper size led to a negligible reduction in stress and a 6% reduction in strain. The height of the device had the greatest impact on the resultant stress and strain.

Energy Separation and Explicit Dynamic Analysis of Low Temperature Impact Toughness of Transmission Tower Material Q420B

To determine the physical significance of the impact toughness parameters and accurately characterize the low temperature impact toughness of transmission tower material Q420 B,the finite element model of Charpy impact test is established on the basis of experiment.The simulation and test results are verified,and the specimen fracture is analyzed by scanning electron microscope.The formation and growth mechanism of the crack are dynamically analyzed.On this basis,energy separation method is used to investigate the effect of low temperature on impact toughness.The results show that the simulation and test results are in good agreement,and the ductile-brittle transition temperature of Q420 B is about-50 ℃.The breaking process of the specimen is divided into the crack formation and propagation.When temperature drops from 20 to-60 ℃,the crack propagation energy decreases from 51.0 to 11.9 J,the crack formation energy reduces from 39.9 to 15.8 J,and the fracture time of the material drops from 1.8 to 0.6 ms.

Numerical simulation of a UAV impacting engine fan blades

A 3D digital model of a small Unmanned Aerial Vehicle(UAV)is obtained by using the method of scanning reverse modeling and joint mapping.A numerical simulation of a small UAV strikes on rotary engine blades,presented in this paper,was performed with a Transient Nonlinear Finite Element code PAM-CRASH software.A test of motor strike on plate was developed and the dynamic response of the plate were obtained to validate the numerical simulation method of a UAV strike on blades.Based on this,dynamic damage response caused by UAV on the engine blades were studied.It is indicated that the impact process between the UAV and a single blade can be divided into two typical stages:cutting and impact.Cutting mainly leads to the failure of the leading edge material,and impact mainly leads to the plastic deformation of the blade.At the same time,it is compared with the damage impacted by bird with the same mass.For the same mass of bird and UAV,the damage caused by UAV striking fan blade is more serious,and 1.345 kg UAV striking fan blade of typical civil aviation engine is enough to cause damage to flight safety.