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).

Seismic Wavelet Analysis Based on Finite Element Numerical Simulation

The practice of exploration and production has proved that explosives are excited in different surrounding rocks and the seismic wavelets collected have different characteristics. In this paper, by establishing a numerical model of the explosion in the well, using finite element analysis technology for numerical simulation, the simulation calculated the stress structure in the near-source area of the earthquake excitation, and extracted the seismic wavelet. The results show that the simulation seismic wavelet characteristics of different thin interbedded sand and mudstone structures have changed significantly. Through excitation simulation, the amplitude and spectrum information of seismic wavelets can be compared and analyzed, and the excitation parameters can be optimized. .

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.

Numerical simulations of earthquake rupture induced by pressure perturbation

The subsurface fluid injection can cause pressure increase within faults,leading to earthquake occurrences.However,the factors controlling earthquake rupture due to pressure perturbation remain poorly understood.To resolve this problem,we simulate the physical processes of earthquake nucleation and rupture on strike-slip faults perturbated by pressure migration based on the slip-weakening law.Multiple kinds of factors,including background stress,fluid injection rates,the area of the pressurized region,fault geometry,and fault friction coefficients,are considered in our simulations.Our simulation results reveal that the ratio of shear stress to normal stress rather than their absolute values controls the rupture behavior.With the large stress ratios,high injection rates,and large pressurized areas,earthquakes are prone to propagate as runaway ruptures.Additionally,faults with large aspect ratios of length to width are also favorable for causing runaway ruptures.In contrast,the factors of fault strike,dip angles and friction coefficients have minor influence on rupture behavior.

Physics-based seismic analysis of ancient wood structure:fault-to-structure simulation

Based on the domain reduction method,this study employs an SEM-FEM hybrid workflow which integrates the advantages of the spectral element method(SEM)for flexible and highly efficient simulation of seismic wave propagation in a three-dimensional(3D)regional-scale geophysics model and the finite element method(FEM)for fine simulation of structural response including soil-structure interaction,and performs a physics-based simulation from initial fault rupture on an ancient wood structure.After verification of the hybrid workflow,a large-scale model of an ancient wood structure in the Beijing area,The Tower of Buddhist Incense,is established and its responses under the 1665 Tongxian earthquake and the 1730 Yiheyuan earthquake are simulated.The results from the simulated ground motion and seismic response of the wood structure under the two earthquakes demonstrate that this hybrid workflow can be employed to efficiently provide insight into the relationships between geophysical parameters and the structural response,and is of great significance toward accurate input for seismic simulation of structures under specific site and fault conditions.

Indentation behavior of a semi-infinite piezoelectric semiconductor under a rigid flat-ended cylindrical indenter

This paper theoretically studies the axisymmetric frictionless indentation of a transversely isotropic piezoelectric semiconductor(PSC)half-space subject to a rigid flatended cylindrical indenter.The contact area and other surface of the PSC half-space are assumed to be electrically insulating.By the Hankel integral transformation,the problem is reduced to the Fredholm integral equation of the second kind.This equation is solved numerically to obtain the indentation behaviors of the PSC half-space,mainly including the indentation force-depth relation and the electric potential-depth relation.The results show that the effect of the semiconductor property on the indentation responses is limited within a certain range of variation of the steady carrier concentration.The dependence of indentation behavior on material properties is also analyzed by two different kinds of PSCs.Finite element simulations are conducted to verify the results calculated by the integral equation technique,and good agreement is demonstrated.

Hermite Finite Element Method for Vibration Problem of Euler-Bernoulli Beam on Viscoelastic Pasternak Foundation

Viscoelastic foundation plays a very important role in civil engineering. It can effectively disperse the structural load into the foundation soil and avoid the damage caused by the concentrated load. The model of Euler-Bernoulli beam on viscoelastic Pasternak foundation can be used to analyze the deformation and response of buildings under complex geological conditions. In this paper, we use Hermite finite element method to get the numerical approximation scheme for the vibration equation of viscoelastic Pasternak foundation beam. Convergence and error estimation are rigourously established. We prove that the fully discrete scheme has convergence order O(τ2+h4), where τis time step size and his space step size. Finally, we give four numerical examples to verify the validity of theoretical analysis.

COUPLED SIMULATION OF 3D ELECTRO-MAGNETO-FLOW FIELD IN HALL-HEROULT CELLS USING FINITE ELEMENT METHOD

Two full 3D steady mathematical models are developed by finite element method （FEM） to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the acquired electromagnetic force as source body forces in Navier-Stokes equations. Effects caused by the ferromagnetic shell, busbar system around, and open boundary problem as well as inside induced current were considered in terms of the magnetic field. Furthermore, a new modeling method is found to set up solid models and then mesh them entirely with so-called structuralized grids, namely hex-mesh. Examples of 75kA prebaked cell with two kinds of busbar arrangements are presented. Results agree with those disclosed in the literature and confirm that the coupled simulation is valid. It is also concluded that the usage of these models facilitates the consistent analysis of the electric field to magnetic field and then flow motion to the greater extent, local distributions of current density and magnetic flux density are very much dependent on the cell structure, the steel shell is a shield to reduce the magnetic field and flow pattern is two dimensional in the main body of the metal pad.

Applications of finite element simulation in orthopedic and trauma surgery

Research in different areas of orthopedic and trauma surgery requires a methodology that allows both a more economic approach and the ability to reproduce different situations in an easy way. Simulation models have been introduced recently in bioengineering and could become an essential tool in the study of any physiological unity, regardless of its complexity. The main problem in modeling with finite elements simulation is to achieve an accurate reproduction of the anatomy and a perfect correlation of the different structures, in any region of the human body. Authors have developed a mixed technique, joining the use of a three-dimensional laser scanner Roland Picza captured together with computed tomography(CT) and 3D CT images, to achieve a perfect reproduction of the anatomy. Finite element(FE) simulation lets us know the biomechanical changes that take place after hipprostheses or osteosynthesis implantation and biological responses of bone to biomechanical changes. The simulation models are able to predict changes in bone stress distribution around the implant, so allowing preventing future pathologies. The development of a FE model of lumbar spine is another interesting application of the simulation. The model allows research on the lumbar spine, not only in physiological conditions but also simulating different load conditions, to assess the impact on biomechanics. Different degrees of disc degeneration can also be simulated to determine the impact on adjacent anatomical elements. Finally, FE models may be useful to test different fixation systems, i.e., pedicular screws, interbody devices or rigid fixations compared with the dynamic ones. We have also developed models of lumbar spine and hip joint to predict the occurrence of osteoporotic fractures, based on densitometric determinations and specific biomechanical models, including approaches from damage and fracture mechanics. FE simulations also allow us to predict the behavior of orthopedic splints applied to the correction of deformities, providing the recovering force-displacement and angle-moment curves that characterize the mechanical behavior of the splint in the overall range of movement.

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.

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.

DEFORMATION ANALYSIS OF SHEET METAL SINGLE-POINT INCREMENTAL FORMING BY FINITE ELEMENT METHOD SIMULATION

Single-point incremental forming （SPIF） is an innovational sheet metal forming method without dedicated dies, which belongs to rapid prototyping technology. In generalizing the SPIF of sheet metal, the deformation analysis on forming process becomes an important and useful method for the planning of shell products, the choice of material, the design of the forming process and the planning of the forming tool. Using solid brick elements, the finite element method（FEM） model of truncated pyramid was established. Based on the theory of anisotropy and assumed strain formulation, the SPIF processes with different parameters were simulated. The resulted comparison between the simulations and the experiments shows that the FEM model is feasible and effective. Then, according to the simulated forming process, the deformation pattern of SPIF can be summarized as the combination of plane-stretching deformation and bending deformation. And the study about the process parameters＇ impact on deformation shows that the process parameter of interlayer spacing is a dominant factor on the deformation. Decreasing interlayer spacing, the strain of one step decreases and the formability of blank will be improved. With bigger interlayer spacing, the plastic deformation zone increases and the forming force will be bigger.

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 Numerical Simulation of Ground Subsidence in Liangjia Colliery

Being aimed at the ground subsidence due to underground coal mining,a numerical model of rock was established and an appropriate method of numerical simulation was put forward.Using the measured subsidence data on the ground,the equivalent mechanical parameters of the rock stratums can be back-calculated by the properly treatment of coal excavation area,then the ground subsidence of other coal mining area can be predicted by FFM.It provided reference for the treatment of the buildings on the ground of this colliery.

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.

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.

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 of influences of grain interaction on rolling textures of fcc metals

A rate dependent crystal plasticity constitutive model considering self and latent hardening in finite element analysis was developed to simulate rolling textures of pure aluminum. By changing the assignment of orientations to finite elements, i.e. assigning the same set of orientations to all elements or different orientations to different elements, the influences of grain interaction on the formation of rolling textures were numerically simulated with this kind of crystal plasticity finite element model. The simulation results reveal that the grains without considering grain interaction rotate faster than those considering grain interaction, and the rotation of grain boundary is slowed down due to the grain interaction. For a good simulation more elements should be assigned to one grain, in which the effects of both the boundary and interior parts of grain contribute to the formation of rolling textures.

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.