Optimizing design of lattice materials based on finite element simulation

The optimized design of simple cross-truss and column lattice structures was carried out by the SolidWorks simulation module.The effective density of the structure was calculated according to the weight reduction requirements proposed by the project.Then,the vari-ation curve between the maximum bearing stress of the unit structure and the structural variables was obtained by simulation.Meanwhile,the mathematical equation between the maximum bearing stress and the structural variables could be obtained through MATLAB fitting.The results indicated that with the decrease in the number of cells,the compressive strength of the prepared column lattice increased(400 to 4 cells,compressive strength 29 MPa to 160 MPa).However,the yield strength increased with the number of cells.The compression strength of the simple cross-truss lattice samples indicated an increase trend with the decrease of the pillar size(an increase of the number of units),reaching 91 MPa(pillar diameter 0.52 mm,number of units 25).While the yield strength increased with the increasing of the number of units.In addition,the additive manufacturing processes of simple cubic lattice and simple cross-pillar lattice were investigated using selective laser melting.The compression performance obtained from the experiment is compared with the simulation results,which are in good agreement.The results of this paper can provide an important reference for optimizing design of lattice materials.

Study on hot deformation behavior of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy using a combination of strain-compensated Arrhenius constitutive model and finite element simulation method

Isothermal hot compression experiments were conducted on homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy to investigate hot deformation behavior at the temperature range of 673-773 K and the strain rate range of 0.001-1 s^(-1)by using a Gleeble-1500D thermo mechanical simulator.Metallographic characterization on samples deformed to true strain of 0.70 illustrates the occurrence of flow localization and/or microcrack at deformation conditions of 673 K/0.01 s^(-1),673 K/1 s^(-1)and 698 K/1 s^(-1),indicating that these three deformation conditions should be excluded during hot working of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy.Based on the measured true stress-strain data,the strain-compensated Arrhenius constitutive model was constructed and then incorporated into UHARD subroutine of ABAQUS software to study hot deformation process of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy.By comparison with measured force-displacement curves,the predicted results can describe well the rheological behavior of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy,verifying the validity of finite element simulation of hot compression process with this complicated constitutive model.Numerical results demonstrate that the distribution of values of material parameters(α,n,Q and ln A)within deformed sample is inhomogeneous.This issue is directly correlated to the uneven distribution of equivalent plastic strain due to the friction effect.Moreover,at a given temperature the increase of strain rate would result in the decrease of equivalent plastic strain within the central region of deformed sample,which hinders the occurrence of dynamic recrystallization(DRX).

Plastic characterization of metals by combining nanoindentation test and finite element simulation

Materials with the same elastic modulus E and representative stress and strain （σr,εr） present similar indentation-loading curves, whatever the value of strain hardening exponent n. Based on this definition, a good approach was proposed to extract the plastic properties or constitutive equations of metals from nanoindentation test combining finite element simulation. Firstly, without consideration of strain hardening, the representative stress was determined by varying assumed representative stress over a wide range until a good agreement was reached between the computed and experimental loading curves. Similarly, the corresponding representative strain was determined with different hypothetical values of strain hardening exponent in the range of 0-0.6. Through modulating assumed strain hardening exponent values to make the computed unloading curve coincide with that of the experiment, the real strain hardening exponent was acquired. Once the strain hardening exponent was determined, the initial yield stress ay of metals could be obtained by the power law constitution. The validity of the proposed methodology was verified by three real metals： AISI 304 steel, Fe andA1 alloy.

Prediction of hot tearing susceptibility of direct chill casting of AA6111 alloys via finite element simulations

To predict hot tearing susceptibility(HTS)during solidification and improve the quality of Al alloy castings,constitutive equations for AA6111 alloys were developed using a direct finite element(FE)method.A hot tearing model was established for direct chill(DC)casting of industrial AA6111 alloys via coupling FE model and hot tearing criterion.By applying this model to real manufacture processes,the effects of casting speed,bottom cooling,secondary cooling,and geometric variations on the HTS were revealed.The results show that the HTS of the billet increases as the speed and billet radius increase,while it reduces as the interfacial heat transfer coefficient at the bottom or secondary water-cooling rate increases.This model shows the capabilities of incorporating maximum pore fraction in simulating hot tearing initiation,which will have a significant impact on optimizing casting conditions and chemistry for minimizing HTS and thus controlling the casting quality.

Finite element simulation of inertia friction welding of superalloy bars

A thermo-mechanical coupling.finite element model was built to investigate the inertia friction welding of GH4169 bars. The remeshing and map solution techniques were adopted. Ttle whole welding process was investigated by adopting an innovative heat generation model and the flywheel rotational speed measured via the experiment. The simulated evolution of axial shortening shows a good agreement with the experiment. In addition, extensive .strain concentration presents in the interface and flash, and the largest ,strain exists near the flash root. Moreover, an intere.sting thermal reflux phenomenon during the cooling stage was found.

Finite element simulation of food transport through the esophageal body

The peristaltic transport of swallowed material in the esophagus is a neuro-muscular function involving the nerve control, bolus-structure interaction, and structure-mechanics relationship of the tissue. In this study, a finite element model (FEM) was developed to simulate food transport through the esophagus. The FEM consists of three components, i.e., tissue, food bolus and peristaltic wave, as well as the interactions between them. The transport process was simulated as three stages, i.e., the filling of fluid, contraction of circular muscle and traveling of peristaltic wave. It was found that the maximal passive intraluminal pressure due to bolus expansion was in the range of 0.8-10 kPa and it increased with bolus volume and fluid viscosity. It was found that the highest normal and shear stresses were at the inner surface of muscle layer. In addition, the peak pressure required for the fluid flow was predicted to be 1-15 kPa at the bolus tail. The diseases of systemic sclerosis or osteogenesis imperfecta, with the remodeled microstructures and mechanical properties, might induce the malfunction of esophageal transport. In conclusion, the current simulation was demonstrated to be able to capture the main characteristics in the intraluminal pressure and bolus geometry as measured experimentally. Therefore, the finite element model established in this study could be used to further explore the mechanism of esophageal transport in various clinical applications.

RESEARCH ON THE SELECTION OF FRICTION MODELS IN THE FINITE ELEMENT SIMULATION OF WARM EXTRUSION

During the process of finite element simulation of precision warm forging, the selection of friction models has a direct effect on the precision accuracy of finite element simulation results. Among all the factors which influence the selection of friction models, the distribution rule of normal stress at the tool-workpiece interface is a key one. To find out the distribution rule of normal stress at the tool-workpiece interface, this paper has made a systematic research on three typical plastic deformation processes: forward extrusion, backward extrusion, and lateral extrusion by a method of finite element simulation. Then on the base of synthesizing and correcting traditional friction models, a new general friction model which is fit for warm extrusion is developed at last.

Finite element simulation of aluminum alloy cross valve forming by multi-way loading

Deformation behavior,temperature evolution and coupled effects have a significant influence on forming process and quality of component formed,which are very complex in forming process of aluminum alloy 7075 cross valve under multi-way loading due to the complexity of loading path and the multiplicity of associated processing parameters.A model of the process was developed under DFEORM-3D environment based on the coupled thermo-mechanical finite element method.The comparison between two process models,the conventional isothermal process model and the non-isothermal process model developed in this study,was carried out,and the results indicate that the thermal events play an important role in the aluminum alloy forming process under multi-way loading.The distributions and evolutions of the temperature field and strain filed are obtained by non-isothermal process simulation.The plastic zone and its extension in forming process of cross valve were analyzed.The results may provide guidelines for the determination of multi-way loading forming scheme and loading conditions of the forming cross valve components.

Finite element simulation on the deep drawing of titanium thin-walled surface part

The deep drawing of titanium thin-walled surface part was simulated based on a self-developed three-dimensional finite element model. After an investigation on forming rules, a virtual orthogonal experimental design was adopted to determine the significance of processing parameters, such as die radius, blank holder force, and friction coefficient, on the forming process. The distributions of thickness and equivalent plastic strain of the drawn part were evaluated. The results show that die radius has a relative major influence on the deep drawing process, followed by friction coefficient and blank holder force.

Finite element simulation and microstructure of two-pass inner spinning process of curved-generatrix cone cylindrical parts with annealing/quenching

A two-pass annealing/quenching internal spinning process with small-end rotations is proposed to form a curved generatrix conical thin-walled shell.That is,annealing at 360°C for 2 h followed by the 1st pass spinning,and finally quenching in ice water after holding for 1 h at 498°C followed by the 2nd pass spinning.ABAQUS finite element software is used to simulate the internal spinning process of the products formed under different forming parameters.The distribution laws of spinning force,the stress and strain under different forming processes were compared and analyzed.The mechanical properties and microstructure of the products are subsequently analyzed.The results show that the strain and the residual stress in the skin area of the formed products under two-pass spinning process more uniform,and the hardness and the mechanical performance are improved.The microstructure of the products formed with the 0.15 mm thickness reduction at the 2nd pass is excellent.And the second phase grain size distributed uniformly in the range of 36μm.Whereas,the second phase particles are broken seriously and the size distribution inhomogeneity is increased when the thickness reduction in the skin area is greater than 0.20 mm at the 2nd pass spinning process.

THREE-DIMENSIONAL FINITE ELEMENT SIMULATION OF TOTAL KNEE JOINT IN GAIT CYCLE

Based on CT scanning pictures from a volunteer＇s knee joint, a three-dimensional finite element model of the healthy human knee joint is constructed including complete femur, tibia, fibular, patellar and the main cartilage and ligaments. This model was validated using experimental and numerical results obtained from other authors. The pressure distribution of contact surfaces of knee joint are calculated and analyzed under the load action of ‘heel strike＇, ‘single limb stance＇ and ‘toe-off＇. The results of the gait cycle are that the contact areas of medial cartilage are larger than that of lateral cartilage; the contact force and contact areas would grow larger with the load increasing; the pressure of lateral meniscus is steady, relative to the significant variation of peak pressure in medial meniscus; and the peak value of contact pressure on all components are usually found at about 4570 of the gait cycle.

Finite element simulation for mechanical response of surface mounted solder joints under different temperature cycling

Nonlinear finite element simulation for mechanical response of surface mounted solder joint under different temperature cycling was carried out. Seven sets of parameters were used in order to evaluate the influence of temperature cycling profile parameters. The results show that temperature cycling history has significant effect on the stress response of the solder joint. Based on the concept of relative damage stress proposed by the authors, it is found that enough high temperature holding time is necessary for designing the temperature cycling profile in accelerated thermal fatigue test.

Finite element simulation of welded thin-walled stainless steel container based on SYSWELD

The temperature field and stress fields of 18 - 8 stainless steel container structure were computed during and after tungsten inert gas （TIG） arc welding based on the SYSWELD software. The convection, radiation and conduction were all considered during the simulation process as well as temperature-dependent material properties. The results show that the peak temperature occurs on the heat source location. Steep temperature gradients are observed ahead of the heat source. Axial tensile stress and hoop compressive stress are observed in the weld seam between cylinder and head. Axial compressive stress and hoop tensile stress are observed near the weld seam between cylinder and heads. Axial compressive stress and hoop tensile stress are observed in the axial weld seam of cylinder. Axial tensile stress and hoop compressive stress are observed near the axial weld seam of cylinder. The aim of the above research is to provide a basic theory and some calculation methods for the thin-walled container welding technology so that the failures of these structures in service due to residual stresses may be minimized.

Finite Element Simulation of In-Stent Restenosis with Tissue Growth Model

In this study, a finite element simulation of in-stent restenosis (ISR) is conducted to simulate the deployment and expansion of a stent in an occluded artery with a contact model and a mechanics-based growth model. A tissue growth model based on the multiplicative decomposition of deformation is applied to investigate the growth of the plaque and artery wall upon the stent’s implantation. Due to the high stresses at the contact points between the stent struts and the tissue, further tissue injury or restenosis is observed. The simulation results show that after the stent deployment, the von Mises stress is significantly larger in the plaque compared to the artery wall, especially in the region that is in contact with the stent. However, the growth of the plaque and artery tends to even out the stress concentration over time. The tissue growth is found to be more significant near the inner wall than the outer layer. A 0.77 mm restenosis is predicted, which agrees with published clinical observations. The features of the artery growth are carefully analyzed, and the underlying mechanism is discussed. This study is the first attempt to apply finite element analysis to artery restenosis, which establishes a framework for predicting ISR’s occurrence and severity. The results also provide insights into understanding the underlying mechanism of in-stent restenosis.

Finite element analysis and simulation for cold precision forging of a helical gear

To investigate the effects of billet geometry on the cold precision forging process of a helical gear, six different billet geometries were designed utilizing the relief-hole principle. And the influences of the billet geometry on the forming load and the deformation uniformity were analyzed by three-dimensional （3D） finite element method （FEM） under the commercial software DEFORM 3D. The billet geometry was optimized to meet lower forming load and better deformation uniformity requirement. Deformation mechanism was studied through the distribution of flow velocity field and effective strain field. The forging experiments of the helical gear were successfully performed using lead material as a model material under the same process conditions used in the FE simulations. The results show that the forming load decreases as the diameter of relief-hole do increases, but the effect of do on the deformation uniformity is very complicated. The forming load is lower and the deformation is more uniform when do is 10 mm.

Hot deformation behavior and finite element simulation of Mg-8.3Gd-4.4Y-1.5Zn-0.8Mn alloy

To study the hot deformation behavior of Mg-8.3 Gd-4.4 Y-1.5 Zn-0.8 Mn(wt%) alloy,hot compression tests were conducted using a Gleeble-3500 thermal simulator at temperatures ranging from 653 to773 K,true strain rates of 0.001-1 s^(-1),and a deformation degree of 60%.Results of hot compression experiments show that the flow stress of the alloy increases with the strain rate.The true stress-true strain curves are corrected by correcting the effect of temperature rise in the deformation process.Activation energy,Q,equal to 287380 J/mol and material constant,n,equal to 4.59 were calculated by fitting the true stress-true strain curves.Then,the constitutive equation was established and verified via finite element simulation.Results of the hot processing map show that the probability of material instability increases with the degree of deformation,which indicates that the material is not suitable for large deformation in a single pass.On the whole,the alloy is appropriate for multipass processing with small deformation and a suitable processing temperature and strain rate are 733 K and 0.01 s-1,respectively.

Numerical simulation analysis for deformation deviation and experimental verification for an antenna thin-wall parts considering riveting assembly with finite element method

In the process of thin-wall parts assembly for an antenna,the parts assembly deformation deviation is occurring due to the riveting assembly.In view of the riveting assembly deformation problems,it can be analyzed through transient and static simulation.In this work,the theoretical deformation model for riveting assembly is established with round head rivet.The simulation analysis for riveting deformation is carried out with the riveting assembly piece including four rivets,which comparing with the measuring points experiment results of riveting test piece through dealing with the experimental data using the point coordinate transform method and the space line fitting method.Simultaneously,the deformation deviation of the overall thin-wall parts assembly structure is analyzed through finite element simulation;and its results are verified by the measuring experiment for riveting assembly with the deformation deviation of some key points on the thin-wall parts.Through the comparison analysis,it is shown that the simulation results agree well with the experimental results,which proves the correctness and effectiveness of the theoretical analysis,simulation results and the given experiment data processing method.Through the study on the riveting assembly for thin-wall parts,it will provide a theoretical foundation for improving thin-wall parts assembly quality of large antenna in future.

Mathematical modeling and simulations of stress mitigation by coating polycrystalline particles in lithium-ion batteries

A chemo-mechanical model is developed to investigate the effects on the stress development of the coating of polycrystalline Ni-rich LiNixMnyCo_(z)O_(2)(x≥0.8)(NMC)particles with poly(3,4-ethylenedioxythiophene)(PEDOT).The simulation results show that the coating of primary NMC particles significantly reduces the stress generation by efficiently accommodating the volume change associated with the lithium diffusion,and the coating layer plays roles both as a cushion against the volume change and a channel for the lithium transport,promoting the lithium distribution across the secondary particles more homogeneously.Besides,the lower stiffness,higher ionic conductivity,and larger thickness of the coating layer improve the stress mitigation.This paper provides a mathematical framework for calculating the chemo-mechanical responses of anisotropic electrode materials and fundamental insights into how the coating of NMC active particles mitigates stress levels.

Effect of stress profile on microstructure evolution of cold-drawn commercially pure aluminum wire analyzed by finite element simulation

The evolution of microstructure in the drawing process of commercially pure aluminum wire （CPAW） does not only depend on the nature of materials, but also on the stress profile. In this study, the effect of stress profile on the texture evolution of the CPAW was systematically investigated by combining the numerical simulation and the microstructure observation. The results show that the tensile stress at the wire center promotes the formation of 〈111〉 texture, whereas the shear stress nearby the rim makes little contribution to the texture formation. Therefore, the 〈111 〉 texture at the wire center is stronger than that in the surface layer, which also results in a higher microhardness at the center of the CPAW under axial loading.2017 Published by Elsevier Ltd on behalf of The editorial office of Journal of Materials Science ＆ Technology.

FINITE ELEMENT MODEL AND SIMULATION OF ROTARY FORGING OF A DISC

The rotary forging process of a disc is simulated by 3-D finite element method.The motion of the rotary the is described as the combination of a revolution round the machine axis and a spin round the rotary die axis. Therefore, the workpiece can be loaded and unloaded partly and cyclically by the cone surface of the rotary the continuously, according with the practical rotary forging process. From the siumulation rasults, the causes of center-thinning during rotary forging of discs are that the locally loading of rotary die made the workpiece center get high radial and tangential tensile stresses, and then the shortening in axial direction and the elongating in tangential and radial direction occur continuous- ly.